- Research
- Open access
- Published:
Single‑cell RNA sequencing analysis of human embryos from the late Carnegie to fetal development
Cell & Bioscience volume 14, Article number: 118 (2024)
Abstract
Background
The cell development atlas of transition stage from late Carnegie to fetal development (7–9 weeks) remain unclear. It can be seen that the early period of human embryos (7–9 weeks) is a critical research gap. Therefore, we employed single‑cell RNA sequencing to identify cell types and elucidate differentiation relationships.
Results
The single‑cell RNA sequencing analysis determines eighteen cell clusters in human embryos during the 7–9 weeks period. We uncover two distinct pathways of cellular development and differentiation. Initially, mesenchymal progenitor cells differentiated into osteoblast progenitor cells and neural stem cells, respectively. Neural stem cells further differentiated into neurons. Alternatively, multipotential stem cells differentiated into adipocyte, hematopoietic stem cells and neutrophil, respectively. Additionally, COL1A2-(ITGA1 + ITGB1) mediated the cell communication between mesenchymal progenitor cells and osteoblast progenitor cells. NCAM1-FGFR1 facilitated the cell communication between mesenchymal progenitor cells and neural stem cells. Notably, NCAM1-NCAM1 as a major contributor mediated the cell communication between neural stem cells and neurons. Moreover, CGA-FSHR simultaneously mediated the communication between multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil. Distinct cell clusters activated specific transcription factors such as HIC1, LMX1B, TWIST1, and et al., which were responsible for their specific functions. These coregulators, such as HOXB13, VSX2, PAX5, and et al., may mediate cell development and differentiation in human embryos.
Conclusions
We provide the cell development atlas for human embryos (7–9 weeks). Two distinct cell development and differentiation pathways are revealed.
Background
Researches on embryonic development have predominantly focused model organisms such as mice and zebrafish [1,2,3,4]. Single-cell sequencing has established the cell atlas of embryonic development for zebrafish and frog [2, 4, 5]. In mouse embryos, the genomic expression activity of 116,312 single cells (from 6.5 to 8.5 days) is examined. This analysis has provided the gene activation information at various time points within the 48 h of gastrulation formation [3]. Furthermore, a comprehensive atlas of 2 million cells from 61 mouse embryos (from 9.5 to 13.5 days) has been constructed [1]. In Macaca fascicularis, six Carnegie stage embryos (CS8 to CS11) are collected and 56,636 single cells are analyzed. The primary cell types of gastral motor stages including protonogenesis, somatic cell development, intestinal duct formation, neural tube formation and neural crest differentiation are identified [6].
The initial 60 days (8 weeks) of human embryonic development can be categorized into 23 Carnegie stages [7, 8]. The cell differentiation and morphogenesis at various developmental stages can be characterized by the formation of zygote, pronuclei, morula, blastocyst, free blastocyst and the attaching blastocyst (CS1-CS4, Week 1), the implantation and formation of extraembryonic mesoderm, primitive streak, gastrulation (CS5-CS6, Week 2), the development and formation of gastrulation, notochordal process, primitive pit, notochordal canal, somite number 1–3 neural folds, cardiac primordium, head fold (CS7-CS9, Week 3), the somitogenesis of somite number 4–12 neural fold fuses, somite number 13–20 rostral neuropore closes, somite number 21–29 caudal neuropore closes, somite number 30 leg buds, lens placode, pharyngeal arches (CS10-CS13, Week 4), the development and formation of lens pit, optic cup, lens vesicle, nasal pit, hand plate (CS14-CS15, Week 5), the development and formation of nasal pits moved ventrally, auricular hillocks, foot plate, finger rays (CS16-CS17, Week 6), the ossification commences, straightening of trunk (CS18-CS19, Week 7), the upper limbs longer and bent at elbow, hands and feet turned inward, eyelids, external ears, rounded head, body and limbs (CS20-CS23, Week 8) [9]. Stage 23 is the final stage of the embryo, and subsequent development is called fetal development.
At present, the development mechanisms of human embryos from stage CS1 to CS7 (1 day–2.5 weeks) have been analyzed [10,11,12,13]. Furthermore, human embryos aged 16 to 19 days have demonstrated differentiation into eleven distinct cell types including progenitor cells of the blood system and primordial germ cell [14]. In some studies, the scRNA-seq combined spatial transcriptome technology is employed to construct the cell atlas of early human embryos from stage CS12 to CS16 (4–6 weeks). Meanwhile, the single cell sequencing data of CS12–CS16 (4–6 weeks) is integrated with the fetal organ data aged 10 to 26 weeks. This analysis has characterized the ancestral state of each fetal organ and facilitated the discovery of previously unknown cell types [15]. Moreover, the gene expression atlas of four million single cells from 121 human fetal samples (72 to 129 days) are constructed using sci-RNA-seq3 [16]. The chromatin accessibility of more than 800,000 fetal single cells is detected using sci-ATAC-seq3. More than one million sites of chromatin accessibility have been found and many of which are only present in certain cells [17].
However, the cell development atlas and differentiation of transition stage from late Carnegie to fetal development (7–9 weeks) have not been reported. It can be seen that the early period of human embryos (7–9 weeks) is a critical research gap that is largely unknown. In this study, we accomplished single-cell sequencing and identified the distinct cell types and their differentiation relationships of human embryos (7–9 weeks). These findings have significant implications not only for the investigation of embryonic cell function and differentiation, but the reconstruction of embryonic cell differentiation trajectory in human being.
Methods
Collection of human embryos and preparation of single‑cell suspensions
This study was reviewed and approved by the Institutional Review Boards of the Rugao Third People’s Hospital and Nantong University Hospital (2020-K013). For scRNA-seq, the human embryos (one 7-week-old human embryo, one 8-week-old human embryo and one 9-week-old human embryo) were segmented into pieces (2 × 5 × 5 = 50 mm3), respectively. The dissociation of human embryonic cells was performed according to the previous study [18]. The tissue pieces were digested with 2 mL tissue digestible solution at 37 °C for 15 min. The cell debris and other impurities were removed using a 40-μm sterile cell filter. Red blood cell lysate of 2 mL was used to lyse red blood cells at 25 °C for 10 min. The cells were counted using the TC20 Automatic Cell Counter (Bio-Rad).
Library preparation and data processing
The single-cell suspension density was adjusted to 1 × 105 cells/mL. The single-cell suspension was loaded into a microfluidic plate and segmented into a single pore chip. Next, cell barcode beads were introduced into the microchip and cleaned. A single-cell lysis buffer (100 μL) was added to the chip, leading to cell lysis and mRNA capture. The cell barcode beads, along with the captured mRNAs, were ejected from the chip for subsequent reverse transcription, cDNA amplification and library construction. After fragment size selection and purification, 150-bp paired-ended sequencing of the single-cell library was performed by Illumina HiSeq × 10 machine. The data of single cell sequencing was in FASTQ format. CeleScope software package was used for comparison and quantification and gene expression matrix was generated.
Dimension‑reduction and clustering analysis
The dimension reduction and clustering analysis of the single-cell RNA expression matrix was performed using Seurat package. Each matrix was loaded to create a Seurat object. Low quality cells were filtered using the same criteria for each sample (the number of distinct RNA molecules is between 200 and 6000, the total RNA count is no more than 30,000, and the percentage of mitochondrial should be less than 30%). After quality control, all Seurat objects were merged into a single object. Subsequently, the matrix was normalized using SCTransform algorithm, and 3000 highly variable genes were obtained using FindVariableFeatures function. In addition, the RunPCA function was used for principal component analysis and 50 PCs were retained. The FindClusters and FindNeighbors functions were used for cluster analysis, and the dimension‑reduction of cell clusters was visualized by UMAP method.
Cell type annotation and enrichment analyses
By referring to the curated marker lists on the CellMarker and PanglaoDB database, the cell clusters were annotated. In addition, the marker genes which were reported in published literatures were also used in the annotation of cell types (Supplemental Material Table 1). The FindAllMarkers function (min.pct = 0.25, logfc.threshold = 0.25) was used to screen differentially expressed genes. Subsequently, the cell type-specific marker genes were illustrated. The proportion of cell clusters was also counted. About the differentially expressed genes, the samples were further divided into three groups (Group 1: 8-week-old vs 7-week-old human embryo, Group 2: 9-week-old vs 8-week-old human embryo, Group3: 9-week-old vs 7-week-old human embryo). The differentially expressed genes of cell clusters in different groups were screened. The cell cluster-specific genes (up-regulated or down-regulated in multiple groups or only in a group) were further screened. In regard to investigate the differentially expressed genes during development, we firstly screened the high expression genes of cell clusters in Group 1, Group 2 and Group 3, respectively. We further intersect the genes that are highly expressed in these three groups. The intersect genes are the high expression genes during development. Finally, clusterProfiler v3.6.1 package was used for the enrichment analysis. Biological process (BP) with p_adj value < 0.05 were defined as significantly enriched.
Trajectory analysis and RNA velocity
Monocle3 was used to construct the differentiation trajectory of the cell clusters. Genes which were expressed in at least 1% of cells and whose empirical discrete value was greater than the fitting discrete value were selected for pseudotime analysis. DDRTree method was used for dimension‑reduction and the reconstruction of cell differentiation trajectory. For RNA velocity, BAM files containing individual cell populations were converted to LOOM files. LOOM files were then converted into clipped and unclipped matrix files. The results were visualized by UMAP algorithm.
Cell–cell communication analysis
Cell–cell communication was performed using the CellChat package. The Seurat object was used to create a new CellChat object. The human ligand-receptor database (CellChatDB) provided by CellChat was imported. The over-expressed genes (ligand-receptors) of cell clusters were identified. The over-expressed interactions between over-expressed ligand-receptors were then identified. The probability of intercellular communication was calculated and the communication networks were predicted. The “identifyCommunicationPatterns” function was used to analyze the outgoing communication patterns of secretory cells and the incoming communication patterns of target cells.
Transcription factor analysis
The transcription factor analysis was conducted using pySCENIC package. Before analysis, the filtered expression matrix (Seurat object) was converted to Loom format using the seurat-disk R package. The subsequent pySCENIC processes consists of four major stages: co-expressed genes were used to infer gene regulatory networks, candidate target genes were filtered according to known transcription factor binding sites, AUC (area under the recovery curve) was calculated according to the enriched gene list, and the AUC matrix was binarized. Module analysis was investigated using the Connection Specificity Index (CSI) parameter which were reported in previous studies [19, 20]. Regulatory networks were analyzed using Cytoscape [21].
Real-time quantitative PCR assays
Total RNA was extracted from the human embryos (7-week-old, 8-week-old and 9-week-old human embryo, respectively). HiScript II Q RT SuperMix for qPCR (Vazyme Biotech, Nanjing, China) were used to synthesize cDNA. Real-time PCR was performed using AceQ® Universal SYBR® qPCR Master Mix (Vazyme Biotech, Nanjing, China) on the CFX96™ Real-Time PCR System (Bio-Rad, California, America). Primers for the target genes were synthesized by Sangon Biotech (Sangon, Shanghai, China). The sequences of primers are listed in Supplemental Material Table 2. The relative mRNA abundances were calculated and normalized to the mean of β-actin mRNA. 2−ΔΔCt method were used to calculate the fold change of mRNA expression.
Results
Clustering analysis of human embryonic cells
The comprehensive cell development atlas of human embryos (20,156 cells from 7-week-old human embryo, 19,112 cells from 8-week-old human embryo and 28,376 cells from 9-week-old human embryo) were constructed (Fig. 1A). After quality control of scRNA-Seq data (Fig. S1A), 67,644 human embryonic cells were clustered into 18 cell clusters. Violin plots showed the expression of MT-ND1, MT-ND2, MT-CO1, MT-CO2, MT-ATP8, MT-ATP6, MT-CO3, MT-ND3, MT-ND4L, MT-ND4, MT-ND5, MT-ND6, MT-CYB in 7-week-old, 8-week-old and 9-week-old human embryo using single cell sequencing data, respectively (Fig. S1B). qPCR confirmed the higher expression of mitochondrial genes in 8-week-old human embryo than that in 7-week-old and 9-week-old human embryo, respectively (Fig. S1C). The 18 cell clusters include germ cells, multipotential stem cells, mesenchymal progenitor cells, adipocyte, erythrocytes, neurons, myeloid cells, hematopoietic stem cells, endothelial cells, neural stem cells, macrophages, neutrophil, osteoblast progenitor cells, cardiomyocytes, fibroblast, basal cells, chondrocytes and epithelial cells (Fig. 1B). The marker genes which were used for annotation were listed (Fig. S2A). The proportions or number of cell clusters were shown (Fig. 1C, D and Fig. S2B, S2C). Violin plots showed the expression of MORC4, APLP2, DNMT1 in hematopoietic stem cells of 7-week-old, 8-week-old and 9-week-old human embryo using single cell sequencing data, respectively (Fig. S2D). qPCR confirmed the lower expression of MORC4, APLP2, DNMT1 in hematopoietic stem cells of 9-week-old human embryo than that in 7-week-old and 8-week-old human embryo, respectively (Fig. S2E). The top 5 DEGs of cell clusters were also listed (Fig. 1E). The GO analysis of cell clusters was presented (Fig. 1F).
Gene expression changes in human embryo development
On the basis of clustering, we further investigated the gene expression changes of different cell clusters during development. The results showed that the number of upregulated and downregulated genes in germ cells, multipotential stem cells and mesenchymal progenitor cells were dramatic changes during development (Fig. 2A). We further investigated the biological processes of different expression genes (DEGs) of multipotential stem cells and mesenchymal progenitor cells in human embryo development. The results showed that the DEGs of multipotential stem cells participated in the regulation of hemopoiesis, regulation of lymphocyte proliferation, mesenchyme development, leukocyte migration, and et al. (Fig. 2B and C). The results showed that the DEGs of mesenchymal progenitor cells participated in the taxis, kidney development, renal system development, muscle tissue development, and et al. (Fig. 2D and E). These results suggest that these two cell clusters have multifunctional differentiation potentials in human embryo development.
Cell differentiation relationships analysis using Monocle3 and RNA velocity
The differentiation relationships among human embryonic cell clusters were calculated by Monocle3. The cell differentiation relationships were found in mesenchymal progenitor cells, osteoblast progenitor cells, neural stem cells and neurons. The mesenchymal progenitor cells were designated as the root in the differentiation pathway. Mesenchymal progenitor cells differentiated into osteoblast progenitor cells and neural stem cells, respectively. Neural stem cells further differentiated into neurons (Fig. 3A). The differentiation extent and orientation of mesenchymal progenitor cells, osteoblast progenitor cells, neural stem cells and neurons were investigated using RNA velocity. The multi-direction arrows showed that mesenchymal progenitor cells and neural stem cells were highly heterogeneous and could differentiate into each direction. The arrow directions of the osteoblast progenitor cells and neurons were consistent, indicating osteoblast progenitor cells and neurons had stable states and differentiated from mesenchymal progenitor cells and neural stem cells, respectively (Fig. 3C and D). The expression genes along the pseudotime of mesenchymal progenitor cells, osteoblast progenitor cells, neural stem cells and neurons could be enriched for the biological process of ossification differentiation, axon development, extracellular matrix organization, and et al. (Fig. 3B). The expression of genes including CCN1, IGFBP5, ID3, HGF, IGFBP3, and CEBPB between mesenchymal progenitor cells and osteoblast progenitor cells were investigated (Fig. 3E). The expression of genes including MAP2, UCHL1, TUBB2B, GAP43, CHL1, and DPYSL2 associated with mesenchymal progenitor cells and neural stem cells were also investigated (Fig. S3A).
Another differentiation pathway could be found between multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil. The multipotential stem cells were designated as the root in this differentiation pathway. Multipotential stem cells differentiated into adipocyte, hematopoietic stem cells and neutrophil, respectively (Fig. 4A). The multi-direction arrows showed that multipotential stem cells were highly heterogeneous and could differentiate into each direction. The arrow directions of the adipocyte, hematopoietic stem cells and neutrophil were consistent, indicating that adipocyte, hematopoietic stem cells and neutrophil had stable states and differentiated from multipotential stem cells (Fig. 4C and D). The expression genes along the pseudotime of multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil could be enriched for the biological process of myeloid cell homeostasis, erythrocyte differentiation, neutrophil migration, and et al. (Fig. 4B). The expression of genes including FAM210B, FBN1, GATA2, PARP1, CDK6, and CLU associated with multipotential stem cells and hematopoietic stem cells were investigated (Fig. 4E). The expression of genes including PPIB, C1QBP, MDK, LGALS3, LBR and FCER1G associated with multipotential stem cells and neutrophil were also investigated (Fig. S3B).
Analysis of cell–cell communication using CellChat
The number and strength of cell–cell interactions across 18 cell clusters was identified using CellChat (Figs. 5A and S4A). The strength of cell–cell interactions across 18 cell clusters in different development stages was also identified (Fig. S5A). The number of ligand-receptors of mesenchymal progenitor cells with other cell clusters were shown in Fig. S4B. Notably, the COL1A2-(ITGA1 + ITGB1) signaling pathway mediated the cell–cell communication between mesenchymal progenitor cells and osteoblast progenitor cells (Figs. 5B and S5B). COL1A2 was found to be highly expressed in mesenchymal progenitor cells and ITGA1 + ITGB1 in osteoblast progenitor cells (Fig. S4C). The number of ligand-receptors of multipotential stem cells with other cell clusters was shown in Fig. S4D. NCAM1-FGFR1 was highly enriched in mesenchymal progenitor cells and neural stem cells (Figs. 5C and S5C). NCAM1-NCAM1 mediated the cell communication between neural stem cells and neurons (Figs. 5D and S5D). The signaling pathways among the patterns of outgoing communication in secretory cells were identified. Mesenchymal progenitor cells, neural stem cells and neurons belong to pattern 3. Pattern 3 contained COLLAGEN, MK, PTN, NCAM, PERIOSTIN, THY1, VTN, SEMA3, VCAM, NOTCH, CADM, ALCAM, CD6, NECTIN, NRG, L1CAM, NPR1, CNTN, GDNF, PROS, NGF, ACTIVIN, NT, AGT and NPR2. Among the patterns of incoming communication of target cells, multipotential stem cells and hematopoietic stem cells belong to pattern 1. Pattern 1 contained VEGF, FSH, THBS, EPHA, DESMOSOME, HSPG, MPZ, GRN, ADGRE5, CDH, EGF, CCL, LIFR, IL1, CSF3, HGF, IGF, IL4, NRG, EPO, SEMA6, AGT and OCLN (Fig. 5E). In this study, CGA-FSHR signaling pathway was found to mediate the cell communication between multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil (Figs. S4E and S5E). We also found that NAMPT-(ITGA5 + ITGB1) signaling pathway mediated the cell communication between multipotential stem cells and neutrophil (Figs. S4F and S5F).
Analysis of transcription factors in cell clusters
The activity of binary regulators was represented in Fig. 6A. HIC1, LMX1B, TWIST1, SP7 and PRRX2 were turned on in mesenchymal progenitor cells. LHX3, EVX2, FOXA1, POU4F3 and TCF4 were turned on in neural stem cells. LHX1, TLX2, NHLH1, ONECUT2 and POU3F2 were turned on in neurons. GATA4, SP7, LMX1B, TWIST1 and LHX4 were turned on in osteoblast progenitor cells. Furthermore, 152 regulons were clustered into nine modules. The cell clusters and regulons of each module were shown. Module 1 included mesenchymal progenitor cells, neurons, neural stem cells and osteoblast progenitor cells, containing the regulons of TCF4, HOXB13, VSX2, PAX5 and NFIA. Module 2 was organized by mesenchymal progenitor cells, osteoblast progenitor cells and fibroblast, including DBP, PRRX2, HIC1, SP7, TWIST1 and FOXD1. Module 9 consisted of macrophages and neutrophil, containing the regulons of LTF, IKZF3, RUNX3, CEBPE and SPI1 (Fig. 6C). The heatmap showed that common regulators were turned on among mesenchymal progenitor cells, neurons, neural stem cells and osteoblast progenitor cells (Fig. 6B). Interaction mapping of transcription factor networks showed that SOX8, SOX11, HDAC2, SIX1, ISL1, SOX10, NKX2-1, FOXA2, FOXA1 and HMGA2 were important regulators of mesenchymal cell differentiation and were commonly switched on among four cell clusters. In addition, PAX6, POU3F2, ISL2, LHX3, NKX6-2, GSX1, LBX1, POU4F1 and NKX2-2 were associated with neuron fate commitment and also commonly switched on. Finally, other common regulators including SOX4, VSX2, POU4F3, NHLH1, LHX4, TCF4, MYF5, ONECUT2, and et al. were involved in the multi-tissue development including endocrine system development, gland development and neural crest cell development, and et al. (Fig. 6D).
Discussion
The human embryonic development stage (7–9 weeks) is the important transition stage which is from the late Carnegie to fetal development. However, the cell atlas and differentiation of embryonic development stage (7–9 weeks) have not been reported. In this study, we have constructed the cell development atlas of human embryos (7–9 weeks) using single-cell RNA sequencing. This study fills the research gap from the late Carnegie to fetal development, which has important significance for investigating the development and differentiation relationship of human embryonic cells during the transition stage.
In this work, we identified multipotential stem cells and mesenchymal progenitor cells in human embryos. These two cell clusters possessed multiple differentiation potentials in the development of human embryos. Based on this prediction, the differentiation relationship of multipotential stem cells, mesenchymal progenitor cells and other cell clusters were performed using Monocle 3. The differentiation pathways involved in mesenchymal progenitor cells were identified by trajectory and RNA velocity analysis. Mesenchymal progenitor cells differentiated into osteoblast progenitor cells and neural stem cells, respectively. Neural stem cells could further differentiate into neurons (Fig. 3A). Many investigations show that the osteoblast progenitor cells differentiate from mesenchymal progenitor cells [22,23,24,25]. Human primitive mesenchymal stem cells could differentiate into neural stem cells [26]. The trans-differentiation of MSCs into neural lineages could be performed by adding UDP-4 which increases the expression of neural genes and protein [27]. MSCs containing NSC-like cells are reported to differentiate into NSC-like cells [28]. Neural stem cells are responsible for restoring the central nervous system and differentiate into neurons [28]. These results show that the potential differentiation relationships of mesenchymal progenitor cells, neurons, neural stem cells and osteoblast progenitor cells in previous studies. Similar differentiation relationships are also observed in this study. In addition, we also found that the multipotential stem cells could further differentiate into adipocyte, hematopoietic stem cells and neutrophil, respectively (Fig. 4A). In this study, although the differentiation relationships of cell clusters are proposed, the important roles and differentiation of cell clusters need to be further verified. But, on the whole, we have identified high expression genes and differentiation relationships of multipotential stem cells and mesenchymal progenitor cells in human embryos. These results lay a foundation for cell sorting and enrichment of multipotential stem cells and mesenchymal progenitor cells. These results also provide a theoretical basis for in vitro induction of osteoblast progenitor cells, neurons, adipocyte, hematopoietic stem cells and neutrophil. It also provides scientific basis for cell replacement therapy of related cell degenerative diseases.
Cell–cell communication among different cell clusters was performed using CellChat. COL1A2 is the osteoblast differentiation factor and participates in osteoblast maturation [29]. The reduction expression of COL1A2 is found to impact the differentiation of mesenchymal cells into osteoblasts [30]. Integrin α1 (ITGA1) and integrin β1 (ITGB1) play important roles in tumorigenesis, progression and metastasis [31, 32]. In this study, the high expression of COL1A2 in mesenchymal progenitor cells and the high expression of ITGA1 and ITGB1 in osteoblast progenitor cells were observed (Fig. S4C). COL1A2-(ITGA1 + ITGB1) signaling pathway were found to mediate the cell–cell communication between mesenchymal progenitor cells and osteoblast progenitor cells (Fig. 5B). These results show that the important roles of COL1A2-(ITGA1 + ITGB1) signaling pathway in the differentiation process from mesenchymal progenitor cells to osteoblast progenitor cells. NCAM1 (Neural cell adhesion molecule 1) is expressed in neural or mesenchymal stem cells [33, 34]. NCAM1 is reported to play important roles in neuronal migration, axonal branching and synaptogenesis [35]. NCAM1 associates with fibroblast growth factor receptor-1 (FGFR1) [36]. NCAM and FGFR are involved in epithelial-mesenchymal transformation [37]. In this study, NCAM1-FGFR1 signaling pathway was found to mediate the cell–cell communication between mesenchymal progenitor cells and neural stem cells and NCAM1-NCAM1 signaling pathway mediated the cell–cell communication between neural stem cells and neurons (Fig. 5C and D). These results show that the potential important roles of NCAM1 signaling pathway in the differentiation process from mesenchymal progenitor cells to neurons. CGA-FSHR signaling pathway are reported to play important roles in the male reproduction [38], oogenesis [39] and the maintain of secondary sexual characteristics [40]. In this study, CGA-FSHR signaling pathway were involved in the communication between multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil (Fig. S4E). These results show that the multipotential roles of CGA-FSHR signaling pathway in the differentiation process from multipotential stem cells to other cells. NAMPT signaling pathway is involved in the maturation and survival of neutrophil [41]. In this study, we found that NAMPT-(ITGA5 + ITGB1) signaling pathway mediated the communication between multipotential stem cells and neutrophil (Fig. S4F). These results show that the potential important roles of NAMPT-(ITGA5 + ITGB1) signaling pathway in the differentiation process from multipotential stem cells to neutrophil in human embryo development. We also found other signaling pathways among mesenchymal progenitor cells, neural stem cells and neurons, including COLLAGEN, MK [42], PTN [43], PERIOSTIN [44], THY1, VTN, SEMA3 [45], VCAM, NOTCH [46], CADM, ALCAM, CD6, NECTIN [47], NRG [48], L1CAM [49], NPR1 [50], CNTN, GDNF [51], NGF [52], ACTIVIN [53], NT [54], AGT, NPR2 (Fig. 5E). These reported signaling pathways may be involved in the differentiation of mesenchymal progenitor cells, neural stem cells and neurons in human embryos. The signaling pathways which mediated cell communication in multipotential stem cells and hematopoietic stem cells, including VEGF [55], FSH [56], EPHA [57], DESMOSOME, HSPG, MPZ, GRN, ADGRE5, CDH [58], EGF [59], CCL, LIFR, IL1, CSF3 [60], HGF [61], IGF [62], IL4, NRG, EPO [63], SEMA6, AGT [64] and OCLN (Fig. 5E) were also screened in this study. The roles of these signaling pathways need to be further explored. Therefore, we have identified many known and unknown signaling pathways and molecules that regulate cell differentiation and development in human embryos, which provide targeted pathways or molecules for human embryo cell engineering and in vitro tissue engineering. In subsequent studies, we can achieve cell induction differentiation through the regulation of relevant signaling pathways or molecules.
The specific transcription factors which activated in distinct cell clusters were revealed in this study (Fig. 6A). For example, HIC1, LMX1B, TWIST1, SP7 and PRRX2 were turned on in mesenchymal progenitor cells in this study. HIC1 [65] transforms mesenchymal stem cells (MSCs) into tumor stem cells. LMX1B [66] is involved in the migration and specification of neural crest-derived mesenchymal progenitor cells. TWIST1 [67] and PRRX2 [68] are turned on in mesenchymal progenitor cells. LHX3, EVX2, FOXA1, POU4F3 and TCF4 were turned on in neural stem cells in this study. LHX3 [69] directs the specification of spinal motor neurons (MNs) in embryos. EVX2 [70] specifies the neurotransmitter fates of particular neurons. FOXA1 [71] promotes the differentiation of neural stem-like cells. TCF4 [72] participates in the nervous system development and neural precursor proliferation. LHX1, TLX2, NHLH1, ONECUT2 and POU3F2 were turned on in neurons in this study. LHX1 [73] is required for the maintenance and regeneration of serotonergic neurons in planarians. TLX2 [74] is identified as an important regulator in cortical interneuron development. NHLH1 [75] is involved in neuronal differentiation. ONECUT2 [76] plays an important role in maintaining neuronal characteristics. POU3F2 [77] participates in neuronal differentiation. GATA4, SP7, LMX1B, TWIST1 and LHX4 were turned on in osteoblast progenitor cells in this study. GATA4 [78] regulates osteoblast differentiation and mineralization. SP7 [79] is involved in osteogenic and chondrogenic differentiation from MSCs. TWIST1 [80] is involved in the osteogenic differentiation of MSCs. LHX4 was turned on in osteoblast progenitor cells in this study. These findings suggest that the transcription factors which specifically turned on in distinct cell clusters are responsible for the specific function of each cell cluster. Additionally, the common regulons were turned on in different cell clusters using module analysis (Fig. 6C). For example, module 1 consisted of mesenchymal progenitor cells, neurons, neural stem cells and osteoblast progenitor cells, containing TCF4 [72] which participated in the nervous system development and neural precursor proliferation, HOXB13 [81] which was involved in stem cell differentiation, VSX2 [82] which was a determinant for the fate specification of spinal V2a interneurons, PAX5 and NFIA. Module 2 was organized by mesenchymal progenitor cells, osteoblast progenitor cells and fibroblast, including DBP, PRRX2 [83] which participated in mesenchymal cell differentiation, HIC1 [65] which transformed mesenchymal stem cells into tumor stem cells, SP7 [79] which was involved in osteogenic and chondrogenic differentiation from MSCs, TWIST1 and FOXD1 [84] which was a regulator of osteoblast differentiation. Finally, the interaction networks of co-expression transcription factors showed the complex and robust interaction relationships in mesenchymal progenitor cells, neurons, neural stem cells and osteoblast progenitor cells. The enrichment results of common transcription factors showed that the co-expression transcription factors participated in mesenchymal cell differentiation, neuron fate commitment and multi-tissue development including endocrine system development, gland development and neural crest cell development (Fig. 6D). These co-regulators which turned on between mesenchymal progenitor cells, osteoblast progenitor cells, neural stem cells and neurons may mediate the differentiation of mesenchymal progenitor cells to osteoblast progenitor cells and neural stem cells which can further differentiate into neurons. The samples we collected were in the transitional stage of Carnegie stage to fetal development. Therefore, those transcription factors which were specifically turned on in distinct cell clusters play important roles in the stem maintenance and function of related cell clusters. However, the co-expressed transcription factors are vital for the differentiation among correlated cell clusters.
Conclusions
In this study, we construct the cell development atlas of human embryos (7–9 weeks). It reveals eighteen distinct cell clusters and two pathways for cell development and differentiation. The potential differentiation relationships and mechanisms between mesenchymal progenitor cells, osteoblast progenitor cells, neural stem cells and neurons are determined. The differentiation relationships and mechanisms of multipotential stem cells, adipocyte, hematopoietic stem cells and neutrophil are also identified (Fig. 7). Although these differentiation relationships and mechanisms need to be further investigating, this study fills the research gap of the cell development atlas and differentiation relationships from the late Carnegie to fetal development. In the future study, the comprehensive cell atlas of the human embryonic cell during embryogenesis could be constructed based on this work and available public data. Our findings provide indispensable resources for investigating the development and differentiation of human embryos. This research also has important significance for the mechanism investigation of aborted embryos, which are caused by abnormal development of embryonic cells.
Availability of data and materials
References
Cao J, Spielmann M, Qiu X, Huang X, Ibrahim DM, Hill AJ, et al. The single-cell transcriptional landscape of mammalian organogenesis. Nature. 2019;566(7745):496–502.
Farrell JA, Wang YQ, Riesenfeld SJ, Shekhar K, Regev A, Schier AF. Single-cell reconstruction of developmental trajectories during zebrafish embryogenesis. Science. 2018. https://doi.org/10.1126/science.aar3131.
Pijuan-Sala B, Griffiths JA, Guibentif C, Hiscock TW, Jawaid W, Calero-Nieto FJ, et al. A single-cell molecular map of mouse gastrulation and early organogenesis. Nature. 2019;566(7745):490–5.
Wagner DE, Weinreb C, Collins ZM, Briggs JA, Megason SG, Klein AM. Single-cell mapping of gene expression landscapes and lineage in the zebrafish embryo. Science. 2018;360(6392):981–7.
Briggs JA, Weinreb C, Wagner DE, Megason S, Peshkin L, Kirschner MW, et al. The dynamics of gene expression in vertebrate embryogenesis at single-cell resolution. Science. 2018. https://doi.org/10.1126/science.aar5780.
Zhai JL, Guo J, Wan HF, Qi LQ, Liu LZ, Xiao ZY, et al. Primate gastrulation and early organogenesis at single-cell resolution. Nature. 2022;612(7941):732–8.
Ishiyama H, Ishikawa A, Imai H, Matsuda T, Yoneyama A, Yamada S, et al. Spatial relationship between the metanephros and adjacent organs according to the Carnegie stage of development. Anat Rec (Hoboken). 2019;302(11):1901–15.
O’Rahilly R, Muller F. Developmental stages in human embryos: revised and new measurements. Cells Tissues Organs. 2010;192(2):73–84.
Flierman S, Tijsterman M, Rousian M, de Bakker BS. Discrepancies in Embryonic Staging: Towards a Gold Standard. Life (Basel). 2023;13(5):1084.
Deglincerti A, Croft GF, Pietila LN, Zernicka-Goetz M, Siggia ED, Brivanlou AH. Self-organization of the in vitro attached human embryo. Nature. 2016;533(7602):251–4.
Shahbazi MN, Jedrusik A, Vuoristo S, Recher G, Hupalowska A, Bolton V, et al. Self-organization of the human embryo in the absence of maternal tissues. Nat Cell Biol. 2016;18(6):700–8.
Xiang LF, Yin Y, Zheng Y, Ma YP, Li YG, Zhao ZG, et al. A developmental landscape of 3D-cultured human pre-gastrulation embryos. Nature. 2020;577(7791):537–42.
Zhou F, Wang R, Yuan P, Ren YX, Mao YN, Li R, et al. Reconstituting the transcriptome and DNA methylome landscapes of human implantation. Nature. 2019;572(7771):660–4.
Tyser RCV, Mahammadov E, Nakanoh S, Vallier L, Scialdone A, Srinivas S. Single-cell transcriptomic characterization of a gastrulating human embryo. Nature. 2021;600(7888):285–9.
Xu YC, Zhang TJ, Zhou Q, Hu MZ, Qi Y, Xue YF, et al. A single-cell transcriptome atlas profiles early organogenesis in human embryos. Nat Cell Biol. 2023;25(4):604–15.
Cao J, O’Day DR, Pliner HA, Kingsley PD, Deng M, Daza RM, et al. A human cell atlas of fetal gene expression. Science. 2020;370(6518):eaba7721.
Domcke S, Hill AJ, Daza RM, Cao JY, O’Day DR, Pliner HA, et al. A human cell atlas of fetal chromatin accessibility. Science. 2020;370(6518):eaba7612.
Shi JW, Fok KL, Dai PY, Qiao F, Zhang MY, Liu HG, et al. Spatio-temporal landscape of mouse epididymal cells and specific mitochondria-rich segments defined by large-scale single-cell RNA-seq. Cell Discov. 2021;7(1):34.
Fuxman Bass JI, Diallo A, Nelson J, Soto JM, Myers CL, Walhout AJM. Using networks to measure similarity between genes: association index selection. Nat Methods. 2013;10(12):1169–76.
Zhu QQ, Tan MY, Wang CF, Chen YF, Wang CF, Zhang JQ, et al. Single-cell RNA sequencing analysis of the temporomandibular joint condyle in 3 and 4-month-old human embryos. Cell Biosci. 2023;13(1):130.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Chen Q, Shou P, Zheng C, Jiang M, Cao G, Yang Q, et al. Fate decision of mesenchymal stem cells: adipocytes or osteoblasts? Cell Death Differ. 2016;23(7):1128–39.
Hu LF, Yin C, Zhao F, Ali A, Ma JH, Qian AR. Mesenchymal stem cells: cell fate decision to osteoblast or adipocyte and application in osteoporosis treatment. Int J Mol Sci. 2018;19(2):360.
Sacramento CM, Assis RIF, Saito MT, Coletta RD, da Rocha DM, Sallum EA, et al. BMP-2 and asporin expression regulate 5-aza-dC-mediated osteoblast/cementoblast differentiation of periodontal dental ligament mesenchymal progenitor cells. Differentiation. 2022;124:17–27.
Vidoni C, Ferraresi A, Secomandi E, Vallino L, Gardin C, Zavan B, et al. Autophagy drives osteogenic differentiation of human gingival mesenchymal stem cells. Cell Commun Signal. 2019;17(1):98.
Brown C, McKee C, Halassy S, Kojan S, Feinstein DL, Chaudhry GR. Neural stem cells derived from primitive mesenchymal stem cells reversed disease symptoms and promoted neurogenesis in an experimental autoimmune encephalomyelitis mouse model of multiple sclerosis. Stem Cell Res Ther. 2021;12(1):499.
Koutsoumparis AE, Patsiarika A, Tsingotjidou A, Pappas I, Tsiftsoglou AS. Neural differentiation of human dental mesenchymal stem cells induced by ATRA and UDP-4: a comparative study. Biomolecules. 2022;12(2):218.
Kaminska A, Radoszkiewicz K, Rybkowska P, Wedzinska A, Sarnowska A. Interaction of neural stem cells (NSCs) and mesenchymal stem cells (MSCs) as a promising approach in brain study and nerve regeneration. Cells. 2022;11(9):1464.
Komori T. Whole aspect of Runx2 functions in skeletal development. Int J Mol Sci. 2022;23(10):5776.
Hill CR, Jacobs BH, Brown CB, Barnett JV, Goudy SL. Type III transforming growth factor beta receptor regulates vascular and osteoblast development during palatogenesis. Dev Dyn. 2015;244(2):122–33.
Boudjadi S, Bernatchez G, Senicourt B, Beausejour M, Vachon PH, Carrier JC, et al. Involvement of the integrin alpha1beta1 in the progression of colorectal cancer. Cancers (Basel). 2017;9(8):96.
Nurzat Y, Su WJ, Min PR, Li K, Xu H, Zhang YX. Identification of therapeutic targets and prognostic biomarkers among integrin subunits in the skin cutaneous melanoma microenvironment. Front Oncol. 2021;11: 751875.
Bonfanti L. PSA-NCAM in mammalian structural plasticity and neurogenesis. Prog Neurobiol. 2006;80(3):129–64.
Skog MS, Nystedt J, Korhonen M, Anderson H, Lehti TA, Pajunen MI, et al. Expression of neural cell adhesion molecule and polysialic acid in human bone marrow-derived mesenchymal stromal cells. Stem Cell Res Ther. 2016;7(1):113.
Vukojevic V, Mastrandreas P, Arnold A, Peter F, Kolassa IT, Wilker S, et al. Evolutionary conserved role of neural cell adhesion molecule-1 in memory. Transl Psychiatry. 2020;10(1):217.
Francavilla C, Cattaneo P, Berezin V, Bock E, Ami D, de Marco A, et al. The binding of NCAM to FGFR1 induces a specific cellular response mediated by receptor trafficking. J Cell Biol. 2009;187(7):1101–16.
Cirovic S, Vjestica J, Mueller CA, Tatic S, Vasiljevic J, Milenkovic S, et al. NCAM and FGFR1 coexpression and colocalization in renal tumors. Int J Clin Exp Pathol. 2014;7(4):1402–14.
Choi D. Expressions of gonadotropin subunit genes in the reproductively inactive golden hamsters. Dev Reprod. 2022;26(2):37–47.
Zhou N, Wang N, Qin X, Liu Q, Wang JH, Dong H, et al. Expression of follicle-stimulating hormone receptor (FSHR), protein kinase B-2 (AKT2) and adapter protein with PH domain, PTB domain, and leucine zipper (APPL1) in pig ovaries. Pol J Vet Sci. 2017;20(4):661–7.
Kleppe L, Andersson E, Skaftnesmo KO, Edvardsen RB, Fjelldal PG, Norberg B, et al. Sex steroid production associated with puberty is absent in germ cell-free salmon. Sci Rep. 2017;7(1):12584.
Siakaeva E, Pylaeva E, Spyra I, Bordbari S, Hoing B, Kurten C, et al. Neutrophil maturation and survival is controlled by IFN-dependent regulation of NAMPT signaling. Int J Mol Sci. 2019;20(22):5584.
Kadomatsu K, Muramatsu T. Midkine and pleiotrophin in neural development and cancer. Cancer Lett. 2004;204(2):127–43.
Tang CY, Wang M, Wang PJ, Wang L, Wu QF, Guo WX. Neural stem cells behave as a functional niche for the maturation of newborn neurons through the secretion of PTN. Neuron. 2019;101(1):32-44 e36.
Liu YF, Luan Y, Guo ZY, Liu Y, Liu CX. Periostin attenuates oxygen and glucose deprivation-induced death of mouse neural stem cells via inhibition of p38 MAPK activation. Neurosci Lett. 2022;774: 136526.
Sawada M, Ohno N, Kawaguchi M, Huang SH, Hikita T, Sakurai Y, et al. PlexinD1 signaling controls morphological changes and migration termination in newborn neurons. EMBO J. 2018;37(4): e97404.
Zhang YN, Wu XY, Liu QY, Yan XN, Liu HZ, Feng DR. Analysis of COL1A1 and COL1A2 gene variants in two fetuses with osteogenesis imperfecta. Zhonghua Yi Xue Yi Chuan Xue Za Zhi. 2023;40(7):821–7.
Lim ST, Lim KC, Giuliano RE, Federoff HJ. Temporal and spatial localization of nectin-1 and l-afadin during synaptogenesis in hippocampal neurons. J Comp Neurol. 2008;507(2):1228–44.
Kataria H, Alizadeh A, Karimi-Abdolrezaee S. Neuregulin-1/ErbB network: an emerging modulator of nervous system injury and repair. Prog Neurobiol. 2019;180: 101643.
Kiefel H, Pfeifer M, Bondong S, Hazin J, Altevogt P. Linking L1CAM-mediated signaling to NF-kappaB activation. Trends Mol Med. 2011;17(4):178–87.
Aballay A. Neural regulation of immunity: role of NPR-1 in pathogen avoidance and regulation of innate immunity. Cell Cycle. 2009;8(7):966–9.
Wei RF, Qiu X, Wang SX, Li Y, Wang YW, Lu KF, et al. NEDL2 is an essential regulator of enteric neural development and GDNF/Ret signaling. Cell Signal. 2015;27(3):578–86.
Guo T, Mandai K, Condie BG, Wickramasinghe SR, Capecchi MR, Ginty DD. An evolving NGF-Hoxd1 signaling pathway mediates development of divergent neural circuits in vertebrates. Nat Neurosci. 2011;14(1):31–6.
Matulka K, Lin HH, Hribkova H, Uwanogho D, Dvorak P, Sun YM. PTP1B is an effector of activin signaling and regulates neural specification of embryonic stem cells. Cell Stem Cell. 2013;13(6):706–19.
Chang HM, Wu HC, Sun ZG, Lian F, Leung PCK. Neurotrophins and glial cell line-derived neurotrophic factor in the ovary: physiological and pathophysiological implications. Hum Reprod Update. 2019;25(2):224–42.
Fragoso R, Igreja C, Clode N, Henriques A, Appleton C, Zhu ZP, et al. VEGF signaling on hematopoietic precursors restricts B-lymphoid commitment in vitro and in vivo. Exp Hematol. 2008;36(10):1329–36.
Mierzejewska K, Borkowska S, Suszynska E, Suszynska M, Poniewierska-Baran A, Maj M, et al. Hematopoietic stem/progenitor cells express several functional sex hormone receptors-novel evidence for a potential developmental link between hematopoiesis and primordial germ cells. Stem Cells Dev. 2015;24(8):927–37.
Nguyen TM, Arthur A, Zannettino ACW, Gronthos S. EphA5 and EphA7 forward signaling enhances human hematopoietic stem and progenitor cell maintenance, migration, and adhesion via Rac1 activation. Exp Hematol. 2017;48:72–8.
Yang YR, Kabir MH, Park JH, Park JI, Kang JS, Ju S, et al. Plasma proteomic profiling of young and old mice reveals cadherin-13 prevents age-related bone loss. Aging (Albany NY). 2020;12(9):8652–68.
Piryani SO, Kam AYF, Kliassov EG, Chen BJ, Spector NL, Chute JP, et al. Epidermal growth factor and granulocyte colony stimulating factor signaling are synergistic for hematopoietic regeneration. Stem Cells. 2018;36(2):252–64.
Zhang HJ, Coblentz C, Watanabe-Smith K, Means S, Means J, Maxson JE, et al. Gain-of-function mutations in granulocyte colony-stimulating factor receptor (CSF3R) reveal distinct mechanisms of CSF3R activation. J Biol Chem. 2018;293(19):7387–96.
Ilangumaran S, Villalobos-Hernandez A, Bobbala D, Ramanathan S. The hepatocyte growth factor (HGF)-MET receptor tyrosine kinase signaling pathway: Diverse roles in modulating immune cell functions. Cytokine. 2016;82:125–39.
Cheng CW, Adams GB, Perin L, Wei M, Zhou XY, Lam BS, et al. Prolonged fasting reduces IGF-1/PKA to promote hematopoietic-stem-cell-based regeneration and reverse immunosuppression. Cell Stem Cell. 2014;14(6):810–23.
Annese T, Tamma R, Ruggieri S, Ribatti D. Erythropoietin in tumor angiogenesis. Exp Cell Res. 2019;374(2):266–73.
Zhang YC, Cheng T, Yuan WP. Recent studies on PI3K/AKT/mTOR signaling pathway in hematopoietic stem cells. Zhongguo Shi Yan Xue Ye Xue Za Zhi. 2013;21(1):245–9.
Leu YW, Huang TH, Hsiao SH. Epigenetic reprogramming of mesenchymal stem cells. Adv Exp Med Biol. 2013;754:195–211.
Akula M, Park JW, West-Mays JA. Relationship between neural crest cell specification and rare ocular diseases. J Neurosci Res. 2019;97(1):7–15.
Zhu QQ, Ma CH, Wang Q, Song Y, Lv TF. The role of TWIST1 in epithelial–mesenchymal transition and cancers. Tumour Biol. 2016;37(1):185–97.
Higuchi M, Kato T, Yoshida S, Ueharu H, Nishimura N, Kato Y. PRRX1- and PRRX2-positive mesenchymal stem/progenitor cells are involved in vasculogenesis during rat embryonic pituitary development. Cell Tissue Res. 2015;361(2):557–65.
Lee S, Cuvillier JM, Lee B, Shen RK, Lee JW, Lee SK. Fusion protein Isl1-Lhx3 specifies motor neuron fate by inducing motor neuron genes and concomitantly suppressing the interneuron programs. Proc Natl Acad Sci USA. 2012;109(9):3383–8.
Juarez-Morales JL, Schulte CJ, Pezoa SA, Vallejo GK, Hilinski WC, England SJ, et al. Evx1 and Evx2 specify excitatory neurotransmitter fates and suppress inhibitory fates through a Pax2-independent mechanism. Neural Dev. 2016;11:5.
Dong DF, Meng L, Yu QQ, Tan GX, Ding M, Tan YJ. Stable expression of FoxA1 promotes pluripotent P19 embryonal carcinoma cells to be neural stem-like cells. Gene Expr. 2012;15(4):153–62.
Braun K, Haberle BM, Wittmann MT, Lie DC. Enriched environment ameliorates adult hippocampal neurogenesis deficits in Tcf4 haploinsufficient mice. BMC Neurosci. 2020;21(1):50.
Currie KW, Pearson BJ. Transcription factors lhx1/5-1 and pitx are required for the maintenance and regeneration of serotonergic neurons in planarians. Development. 2013;140(17):3577–88.
Lim L, Mi D, Llorca A, Marin O. Development and functional diversification of cortical interneurons. Neuron. 2018;100(2):294–313.
Cogliati T, Good DJ, Haigney M, Delgado-Romero P, Eckhaus MA, Koch WJ, et al. Predisposition to arrhythmia and autonomic dysfunction in Nhlh1-deficient mice. Mol Cell Biol. 2002;22(14):4977–83.
van der Raadt J, van Gestel SHC, Nadif Kasri N, Albers CA. ONECUT transcription factors induce neuronal characteristics and remodel chromatin accessibility. Nucleic Acids Res. 2019;47(11):5587–602.
Lin YJ, Hsin IL, Sun HS, Lin S, Lai YL, Chen HY, et al. NTF3 is a novel target gene of the transcription factor POU3F2 and is required for neuronal differentiation. Mol Neurobiol. 2018;55(11):8403–13.
Miranda-Carboni GA, Guemes M, Bailey S, Anaya E, Corselli M, Peault B, et al. GATA4 regulates estrogen receptor-alpha-mediated osteoblast transcription. Mol Endocrinol. 2011;25(7):1126–36.
Chan WCW, Tan Z, To MKT, Chan D. Regulation and role of transcription factors in osteogenesis. Int J Mol Sci. 2021;22(11):5445.
Miraoui H, Severe N, Vaudin P, Pages JC, Marie PJ. Molecular silencing of Twist1 enhances osteogenic differentiation of murine mesenchymal stem cells: implication of FGFR2 signaling. J Cell Biochem. 2010;110(5):1147–54.
Chung N, Jee BK, Chae SW, Jeon YW, Lee KH, Rha HK. HOX gene analysis of endothelial cell differentiation in human bone marrow-derived mesenchymal stem cells. Mol Biol Rep. 2009;36(2):227–35.
Clovis YM, Seo SY, Kwon JS, Rhee JC, Yeo S, Lee JW, et al. Chx10 consolidates V2a interneuron identity through two distinct gene repression modes. Cell Rep. 2016;16(6):1642–52.
Ueharu H, Higuchi M, Nishimura N, Yoshida S, Shibuya S, Sensui K, et al. Expression of Kruppel-like factor 6, KLF6, in rat pituitary stem/progenitor cells and its regulation of the PRRX2 gene. J Reprod Dev. 2014;60(4):304–11.
Teixeira CC, Liu YX, Thant LM, Pang J, Palmer G, Alikhani M. Foxo1, a novel regulator of osteoblast differentiation and skeletogenesis. J Biol Chem. 2010;285(40):31055–65.
Chen TT, Chen X, Zhang SS, Zhu JW, Tang BX, Wang AK, et al. The genome sequence archive family: toward explosive data growth and diverse data types. Genomics Proteomics Bioinformatics. 2021;19(4):578–83.
Members C-N, Partners. Database resources of the national genomics data center, china national center for bioinformation in 2021. Nucleic Acids Res. 2021;49(D1):D18–28.
Acknowledgements
Not applicable.
Funding
This work was supported by the National Natural Science Foundation of China (82001606 to C. N. W. and 81901528 to X. R. W.), Jiangsu Health Innovation Team Program (2020), "333 Project" Scientific Research Project of Jiangsu Province (BRA2020204 to D. W. X.), Jiangsu “Six Talent Peak” High-level Talent Project (WSW-277 to D. W. X.), Lo Kwee Seong Start Up Fund and startup R&D funding from Nantong University (TDYX2022001 and TDYX2021021 to J.W.S.), Natural Science Foundation of Nantong (JC2023103 to J.W.S.), Large Instruments Open Foundation of Nantong University (KFJN2343 to J.W.S.), the Jiangsu Province College Students' innovation and entrepreneurship training program (202310304123Y to H. Q. C.), Project of Nantong Municipal Health Commission (MS 2023042 to Y.W).
Author information
Authors and Affiliations
Contributions
Conceived and designed the experiments: Jianwu Shi, Dawei Xu, Xingmei Feng, Chengniu Wang, Xiaorong Wang and Hongping Ding; Collected human embryo samples: Chenliang Sun, Wenran Wang, Chunhong Li, JiaYing Li, Yong Wang; Analyzed the data: Jianwu Shi, Yufei Chen, Xiaoran Li, Taowen Ye, Jin Dong, Hanqing Chen, Weizhen Wang; Wrote the paper: Jianwu Shi, Chengniu Wang, Dawei Xu, Xingmei Feng, Xiaorong Wang.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
This study was reviewed and approved by the Institutional Review Boards of the Rugao Third People’s Hospital and Nantong University Hospital (2020-K013). All participants have provided informed consent.
Consent for publication
All authors have read the manuscript and agree to publish.
Competing interests
The authors have declared no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13578_2024_1302_MOESM1_ESM.pdf
Additional file 1: Fig. S1. Quality control and expression of mitochondrial genes. A The number of genes, counts and percentage of mitochondria. B Violin plots showing the expression of mitochondrial genes including MT-ND1, MT-ND2, MT-CO1, MT-CO2, MT-ATP8, MT-ATP6, MT-CO3, MT-ND3, MT-ND4L, MT-ND4, MT-ND5, MT-ND6, MT-CYB in the different development stage of human embryos using single cell sequencing data. C qPCR validated the expression of mitochondrial genes in different development stages. Fig. S2. Marker genes and cell proportion analysis. A The marker genes which were screened from the published literatures. B The cell proportions of mesenchymal progenitor cells, osteoblast progenitor cells, neural stem cells and neurons. C The cell proportions of multipotential stem cells, adipocyte, neutrophil and the number of hematopoietic stem cells. D Violin plots showing the expression levels of MORC4, APLP2, DNMT1 in hematopoietic stem cells of different development stage of human embryos using single cell sequencing data. E qPCR validated the expression of MORC4, APLP2, DNMT1 in the different development stage of human embryos. Fig. S3. The expression of key genes along the pseudotime. A The expression of key genes including MAP2, UCHL1, TUBB2B, GAP43, CHL1 and DPYSL2 between mesenchymal progenitor cells and neural stem cells. B The expression of PPIB, C1QBP, MDK, LGALS3, LBR and FCER1G between multipotential stem cells and neutrophil. Fig. S4. Cell–Cell communication analysis of each cell cluster. A The interaction strength of ligand-receptors of each cell cluster. B The number of ligand-receptors of mesenchymal progenitor cells with other cell clusters. C The expression of COL1A2, ITGA1 and ITGB1 among 18 cell clusters. D The number of ligand-receptors of multipotential stem cells with other cell clusters. E CGA-FSHR signaling pathway in each cell cluster. F NAMPT-signaling pathway in each cell cluster. Fig. S5. The strength and signaling pathways of cell clusters in the different development stage of human embryos using single cell sequencing data. A The interaction strength of cell clusters in different development stages. B COL1A2-signaling pathway of cell clusters in different development stages. C NCAM1-FGFR1 signaling pathway of cell clusters in different development stages. D NCAM1-NCAM1 signaling pathway of cell clusters in different development stages. E CGA-FSHR signaling pathway of cell clusters in different development stages. F NAMPT-signaling pathway of cell clusters in different development stages.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wang, C., Wang, X., Wang, W. et al. Single‑cell RNA sequencing analysis of human embryos from the late Carnegie to fetal development. Cell Biosci 14, 118 (2024). https://doi.org/10.1186/s13578-024-01302-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13578-024-01302-9