Activation of Blood Vessel Development in Endometrial Stromal Cells In Vitro Cocultured with Human Peri-Implantation Embryos Revealed by Single-Cell RNA-Seq.

Human embryos and endometrium tissue were donated by couples with no genetic disorders or diseases from standard in vitro fertilization (IVF) treatment at the Reproductive Medicine Center, Tongji Hospital, with written informed consent. The primary endometrium was dissected using scalpels and then enzymatically digested in 1% collagenase I and IV (Gibco, 17100-017 and 17014-019) by repeatedly pipetting up and down for 30 min at 37 °C. After filtration through 40 μm strainers and centrifugation, the pellet at the bottom was washed twice with phosphate-buffered saline (PBS). Then, the pellet was resuspended, plated, and prepared for use. The procedure of culturing embryos was modified and performed according to previously published studies [ 7 , 8 ]. Briefly, pre-implantation human embryos were cultured according to a previous procedure [ 12 ]. After embryos hatched from the zona pellucida (late blastocyst stage, day 6.5), we transferred them to the dishes with plated endometrial cells and then cocultured them for 2.5 days. The embryos and endometrial cells underwent incubation with 0.25% trypsin before single cells were dissociated with repeated aspiration through a mouth-operated, drawn capillary pipette. Individual endometrial cells and embryonic cells were placed into a 200 μL PCR tube, and single-cell RNA sequencing (RNA-Seq) was performed using the Smart-seq2 method [ 13 ]. Single-cell RNA-Seq data were processed as follows: TrimGalore! was used to trim the fastq data [ 14 ]. Then, the trimmed data were aligned to human reference genome hg38 using STAR (v2.6.0) with default parameters [ 15 ], followed by read counting using featureCounts v1.6.2 [ 16 ]. The cocultured embryonic cells and endometrial cells were separated using unbiased clustering, as described in our previous work [ 10 ]. The trophoblast cells were characterized using 300 epiblast, primitive endoderm, and trophectoderm markers [ 17 ] identified in previous work [ 10 ]. The endometrial cells, including those before and after coculturing with embryos, were identified as EmSCs using the stromal marker, vimentin (VIM) ( Figure S1, Supplementary Materials ). The EmSCs before coculturing were defined as “PRE” and the EmSCs after coculturing were defined as “POST”. Principal component analysis (PCA), unsupervised hierarchal clustering, and DEG identification were carried out using the DEseq2 package (v1.24.0) [ 18 ]. Transcript per million (TPM) was used to normalize the raw read count for subsequent WGCNA analysis [ 19 ]. The transcriptomes of endometrium tissue at WOI and before WOI were downloaded from Hu et al. [ 20 ]. The trophoblast cells in the cocultured embryo and the cocultured EmSCs were used for CellPhone DB analysis [ 21 ]. Genes with a p -value < 0.05 and fold change ≥2 were characterized as DEGs. TB GreenTM Premix Ex TaqTM II (Takara, RR820A) and a Roche Light Cycler 96 system (Roche) were used to carry out qRT-PCR. Relative levels of RNAs were calculated using the ddCt method with GAPDH as an endogenous control. The primers are listed in Table S3 (Supplementary Materials) . WGCNA was performed on normalized gene expression data measured using the TPM metric, with 2352 DEGs determined by the DEseq2 package. WGCNA was performed according to a previously published protocol [ 19 ]. Briefly, a topological overlap matrix (TOM) was constructed with softPower set to three. The relationships between different modules were visualized using the module and eigengene relation heatmap. Genes contained in the modules were extracted and analyzed using Metascape for a precise representation of the data [ 22 ]. The gene coexpression network was extracted and further processed using MCODE in Cytoscape software (V3.7.2) [ 23 ]. GO enrichment analysis was performed on a platform combining the R package “clusterProfiler” [ 24 ] and the Database for Annotation, Visualization, and Integrated Discovery (DAVID) online tool ( https://david.ncifcrf.gov/ (accessed on 13 January 2021)) [ 25 ]. Terms with a p -value < 0.05 were counted as statistically significant. GSEA was performed to detect the key gene sets under the condition of embryo coculture in the datasets of whole EmSCs [ 26 ]. Boxplots were utilized to illustrate selected key genes using the ggplot2 package (V3.2.1) [ 27 ]. The TPM values of selected genes were evaluated using Mann–Whitney U tests.

Implantation is a key event for embryos to continue development [ 1 ]. Successful implantation relies on the synchronized development of the endometrium and embryo [ 2 ]. Human embryos penetrate and begin to invade into endometrial cells during the “window of implantation” (WOI) [ 3 ]. This process requires delicate spatiotemporal coordination by feto-maternal crosstalk, whose failure is considered as the major reason for implantation failure [ 4 ]. Furthermore, formation of blood vessels or angiogenesis in the endometrium is generally considered important for nutrient and oxygen transport between the mother and fetus [ 5 ]. Notably, stromal cells undergoing decidualization, which interact with epithelial cells, maternal blood vessels, and embryonic trophoblast cells, play crucial roles in recruiting development factors and supporting embryo implantation. However, the implantation rate is still lower than expected in natural and assisted reproduction events, and little is known about this process in humans. Considering the unethical and impractical reality of manipulating the human embryo in vivo, approaches to study implantation involve either inference from animal embryo implantation data or the culture of human peri-implantation embryos in vitro [ 6 ]. Recently, two studies developed a new in vitro culture system for human post-implantation embryos [ 7 , 8 ]. Moreover, the in vitro coculture of human embryos with endometrial cells offers us a better way of understanding their interaction [ 6 , 9 ]. Our group and Zhou et al. [ 10 , 11 ] recently published studies focusing on the development of the trophectoderm in peri-implantation conceptuses using coculture and single-cell RNA sequencing methods, respectively. In the present study, we compared EmSCs and trophectoderm cells before and after coculture using the coculture model in vitro, whereby the implantation progress mimicked that in vivo. By analyzing the gene expression profiles, we found significant changes in the EmSCs after coculturing with human embryos compared to the non-cocultured EmSCs. Subsequently, weighted gene coexpression network analysis (WGCNA) and gene ontology (GO) analysis were also applied to analyze the differentially expressed genes (DEGs). The underlying interactions between embryos and EmSCs at these stages were also explored to gain insight into the implantation process and identify potential biomarkers for diagnosis, prognosis, and drug targeting for pregnancy-related malfunctions.

To explore the potential factors involved in the interactions between blastocysts and EmSCs, the non-coculture EmSCs and EmSCs cocultured with day 6.5 pre-implantation blastocysts were included in this study ( Figure 1 ). Using the smart-seq2 method, we collected 24 “PRE” EmSCs and 57 “POST” EmSCs for single-cell RNA-Seq ( Figure S1 and Table S1, Supplementary Materials ). PCA and unsupervised hierarchal clustering analyses revealed that EmSCs in the “PRE” group and the “POST” group formed distinct clusters ( Figure 2 A). Furthermore, among the total of 2352 DEGs in EmSCs found between the “PRE” group and “POST” group using DEG analysis (false discovery rate (FDR) <0.05 and fold change ≥2), 1783 genes, e.g., LRP1 , MUC15 , and COL15A1 , were upregulated and 569 genes, e.g., TMSB4X , STMN1 , and GOLIM4 , were downregulated in the “POST” group ( Figure 2 B and Table S1, Supplementary Materials ), followed by downstream single-cell transcriptomic analysis ( Figure 1 ). To further elucidate the mechanism underlying changes in EmSCs during the implantation process, GO analysis of the DEGs was performed to investigate the overrepresented signals. The results showed that the 1783 upregulated genes in the “POST” group were significantly enriched in biological processes such as extracellular matrix organization, cell–substrate adhesion, female pregnancy, and placenta development ( Figure 2 C). Furthermore, the 569 downregulated genes were also subjected to GO analysis. Intriguingly, most GO items were related to RNA catabolic pathways ( Figure 2 D), indicating the reservation of RNA in the EmSCs to further embryo embedding. Taken together, these results suggest the activation of extracellular matrix pathways in EmSCs to support embryo implantation. Since a previous study compared the transcriptome of whole endometrium tissue at the receptive stage to the pre-receptive stage in the normal menstruation cycle [ 20 ], we sought to determine whether there were any common changes in endometrial cells between in vivo results during WOI without embryo invasion and our in vitro results with embryo invasion. The DEGs at the receptive stage [ 20 ] and the DEGs after implantation in our study [ 10 ] were analyzed. Interestingly, we found 84 shared upregulated genes and 62 shared downregulated genes ( Figure 2 E,F). The 84 shared upregulated genes, including MMP7 , BIRC5 , MUC15 , and MEGF10 , were then subjected to GO analysis to explore the pathways involved during implantation. These upregulated genes were preferentially enriched in the regulation of epidermal growth factor-activated receptor activity, epithelial cell differentiation, and epidermis development ( Figure 2 E), while the 62 downregulated genes were enriched in GO terms such as chromosome and hypoxia pathways ( Figure 2 F). It was reported that BIRC5 acts as a target of interleukin 11 (IL11) during decidualization [ 28 ], whereas MUC15 and MMP7 mediate the invasion of the trophoblast into the endometrium [ 29 , 30 ], and MEGF10 plays an adhesion role according to an in vitro experiment [ 31 ]. Overall, these results indicate the several distinct pathways activated only in EmSCs with the event of embryo implantation. PCA showed that EmSCs were grouped by “PRE” and “POST” states ( Figure 2 A). To systematically investigate the genetic program dynamics, we performed WGCNA on the 2352 DEGs in EmSCs between the “PRE” group and “POST” group. WGCNA identified four gene modules, each of which contained a set of genes that tended to be coexpressed with the stage ( Figure 3 A). By studying module–trait relationships, we found that two gene modules, i.e., the turquoise module and the yellow module, were positively correlated with the “PRE” state, while the other two, i.e., the blue module and the brown module, positively correlated with the “POST” state ( Figure 3 B). These modules represent the core genetic programs that operate during EmSC development at the embryo implantation stage. We then performed GO analysis on the genes in each module to investigate their biological function ( Figure 3 C,D). In the turquoise module, genes were mainly enriched in cell junction organization, extracellular matrix organization, and other functions ( Figure 3 C), while blood vessel development, female pregnancy, and regulation of blood pressure were abundant in the blue module ( Figure 3 D). To further validate whether the pathways described above were significantly activated after coculturing with embryos, GSEA was applied, whereby it was found that the item “blood vessel development” from the GO gene set was significantly enriched in the “POST” group ( Figure 3 E; Table S2, Supplementary Materials ), suggesting a significant difference in blood vessel development after coculturing. Moreover, the only genes in the blue module positively related to the EmSCs of the “POST” group were enriched in blood vessel development ( Figure 3 D). Therefore, we further analyzed the genes in the blue module using the Cytoscape MCODE plugin. The hub genes with a high MCODE score were related to many functions during cell growth and receptivity establishment, including histone modification ( KAT6A and KMT2D ) [ 32 , 33 ], inhibition of endometriosis development ( TNFRSF1A ) [ 34 ], inhibition of inflammatory response ( POM121C ) [ 35 ], and regulation of epithelial–mesenchymal transition ( CDH11 ) [ 36 ]. Intriguingly, a subset of genes related to blood vessel development was identified, including EPN1 , PEAK1 , PCBP1 , PLCG1 , and NCOR2 ( Figure 3 F), which play important roles in hematopoietic stem-cell emergence, endothelial cell recruitment, and vascular branching [ 37 , 38 , 39 , 40 , 41 ]. Revealing the dynamics of blood vessel-related genes impacted by embryos, genes such as RAMP1 , LTBP1 , and LRP1 showed significantly higher expression in the “POST” group than in the “PRE” group ( Figure 4 A). This was then verified by quantitative reverse-transcription polymerase chain reaction ( Figure 4 B). Overall, these findings suggest that blood vessel development might be activated in endometrial cells at the stage of embryo implantation. To systematically analyze the interactions between EmSCs and trophoblast cells, we used CellphoneDB [ 21 ] to identify the expression of possible crosstalk. These interactions were calculated at the transcriptomic level but extrapolated to protein activity. Similar to a previous report using the first-trimester placenta and decidua [ 21 ], we found that growth factors, angiogenesis proteins, and adhesion and recruitment ligands and receptors were highly expressed in the EmSCs and the trophoblast cells during implantation, indicating that these events existed from implantation to the early pregnant placenta ( Figure 5 A). Furthermore, integrins are reported as one of the most important proteins regulating blastocyst implantation [ 42 ]; thus, we sought to determine which integrin and its target were expressed in these two cell types. Interestingly, the EmSCs at the implantation stage highly expressed TNC and VWF , while their receptor integrin aVb3 was highly expressed in the trophoblast cells with relatively low expression in the EmSCs themselves ( Figure 5 B). On the other hand, LAMC1 , FBN1 , and THBS1 were expressed at a high level in the trophoblast cells, while their receptors ( aVb3 , a5b1 , and a3b1 ) were also highly expressed in EmSCs ( Figure 5 B). In addition to integrins, we found a high expression of fibroblast growth factors (FGFs) and cluster of differentiation 44 (CD44) in EmSCs together with a high expression of FGF receptors and CD44 ligands ( HGF , LGALS9 , and FGF2 ) in trophoblast cells ( Figure 5 C), further indicating the preparation of blood vessel establishment. Taken together, these patterns of ligand–receptor expression strongly suggest complex interactions between the EmSCs and trophoblast cells, providing additional insight into the implantation process.

In this study, we explored the influence of human embryos on EmSCs by characterizing the DEGs after implantation at the single-cell level. Despite no decidualization-induced treatment being applied to the EmSCs before coculturing in this study, the EmSCs after coculture showed an activated extracellular matrix and cell–substrate adhesion, highlighting the similarity in stromal cell activity during decidualization at the WOI [ 43 ] ( Figure 2 E,F). Of note, Wang et al. published their work using endometrial single cells during the whole menstrual cycle [ 44 ]. We also compared our single-cell data with theirs. Intriguingly, proliferative-phase genes, such as STC1 , NFATC2 , BMP2 , and PMAIP1 , were downregulated in the “POST” EmSCs compared to the “PRE” EmSCs ( Figure S4, upper panel, Supplementary Materials ), while they were also downregulated at the secretory phase [ 44 ]. However, secretory-phase genes, such as LMCD1 , FGF7 , and IL15 , showed no difference in the cocultured EmSCs as compared with non-cocultured cells, whereas another secretory-phase gene, CRYAB , even showed higher expression before coculture ( Figure S4, lower panel, Supplementary Materials ). These in vitro results suggest changes in stromal cells from the proliferative phase to the secretory phase; however, these changes are not completely possible due to the different states of stromal cells cultured in vitro and the relatively short coculture time with embryos. Implantation is one of the most important processes for successful pregnancy [ 1 ]. Synchronization between the acquisition of invasion competency in embryos and receptivity ability in the endometrium is required for implantation. Consistent with previous findings [ 21 ], our study indicates that the endometrium prepared a receptive state for embryo implantation, as upregulated genes in the “POST” group were enriched in biological processes such as cell junction and adhesion, extracellular matrix organization, and female pregnancy ( Figure 2 C,D). WGCNA is a useful tool to detect specific modules [ 19 ]. By combining this tool with GO, blood vessel development was identified from the blue module ( Figure S2A and Table S2, Supplementary Materials ). The hub genes in this module were also found to be positively related to the “POST” EmSCs ( Figure 3 F and Figure 4 ). Both vasculogenesis and angiogenesis are of crucial importance during blood vessel development when the embryo implants [ 5 ]; accordingly, we found that many genes related to hematopoietic stem-cell emergence, cell growth factors, and chemokines, which were highly expressed in the “POST” EmSCs, were all highly correlated with blood vessel establishment ( Figure 3 and Figure 4 ; Figures S2B and S3 and Table S2, Supplementary Materials ). These results are consistent with the function of the endometrium at the peri-implantation stage, when the uterine spiral arteries begin remodeling to achieve circulation characteristic for successful placentation and the fetal–maternal interface [ 1 , 3 ]. The two-way interactions between the endometrium and trophoblast initiate the implantation process after the embryo enters the uterus [ 45 ]. CellphoneDB (created by Vento-Tormo et al.) has shown its power in analyzing cell–cell interactions [ 21 ]. The number of immunomodulation pathways in our study, including the atypical chemokine receptor 3 (ACKR3)/C–X–C motif chemokine ligand 12 (CXCL12) pathway, was lower compared to the first-trimester placenta data, including programmed death ligand 1 (PDL1)/programmed death 1 (PD1) and adenosine A2B receptor (ADORA2B)/CD39 [ 21 ], which was probably due to the limited cell types included in our studies. However, other ligands/receptors that might regulate the implantation event were found to be highly expressed, such as the T-cell immunoreceptor with immunoglobulin (Ig) and immunoreceptor tyrosine-based inhibitory motif (ITIM) domains (TIGIT)/poliovirus receptor (PVR) pathway, which can inhibit human immune cell activity and might benefit implantation [ 46 ]. Moreover, killer cell lectin-like receptor G2 (KLRG2), responsible for regulating decidualization and embryonic development, possibly through the Wnt pathway [ 47 ], was also highly expressed in our study ( Figure 5 C). Overall, our study provides insight into the endometrial stomal cell state transition during embryo implantation in vitro, showing that receptivity is established after embryo coculture. These findings can help us better understand the intercellular events during implantation. However, our study also had some limitations. Since it is unpractical and unethical to acquire human embryos in vivo during implantation, our study mimicked the implantation process by applying EmSCs to the post-implantation culture method developed by Deglincerti et al., and Shahbazi et al. [ 7 , 8 ]. Therefore, the actual implantation process, with a more complex environment and more cell types in vivo, may undergo different cellular pathways, several of which were shown in our study ( Figure 2 E,F; Figure S4, Supplementary Materials ). Moreover, the differential expression of blood vessel-related genes presented in our study was analyzed at the transcriptomic level, and further validations such as quantitative real-time PCR and/or immunostaining of key pathway signals are required to understand the activity of receptivity state transition. In conclusion, our work provides one piece of the jigsaw puzzle for a better understanding of human embryo implantation, using the in vitro coculturing model of human embryos and EmSCs, which highlighted not only synchronization in signaling pathways and regulatory networks happening in the EmSCs but also an interaction between peri-implantation trophoblast cells and EmSCs.