Next Article in Journal
Nisin S, a Novel Nisin Variant Produced by Ligilactobacillus salivarius P1CEA3
Next Article in Special Issue
Genome-Wide Association Study of Age at First Calving in U.S. Holstein Cows
Previous Article in Journal
Extracellular Vesicles and Their Zeta Potential as Future Markers Associated with Nutrition and Molecular Biomarkers in Breast Cancer
Previous Article in Special Issue
The Role of Inositols in the Hyperandrogenic Phenotypes of PCOS: A Re-Reading of Larner’s Results
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 7
10.3390/ijms24076811
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Rhythm of the First Language: Dynamics of Extracellular Vesicle-Based Embryo–Maternal Communication in the Pre-Implantation Microenvironment

by
Kasun Godakumara
1,
Paul R. Heath
2 and
Alireza Fazeli
1,3,4,*
1
Institute of Veterinary Medicine and Animal Sciences, Estonian University of Life Sciences, 50411 Tartu, Estonia
2
Sheffield Institute for Translational Neuroscience (SITRAN), University of Sheffield, 385a Glossop Rd., Broomhall, Sheffield S10 2HQ, UK
3
Department of Pathophysiology, Institute of Biomedicine and Translational Medicine, Faculty of Medicine, University of Tartu, 14B Ravila, 50411 Tartu, Estonia
4
Academic Unit of Reproductive and Developmental Medicine, Department of Oncology and Metabolism, The Medical School, University of Sheffield, Sheffield S10 2RX, UK
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(7), 6811; https://doi.org/10.3390/ijms24076811
Submission received: 19 February 2023 / Revised: 27 March 2023 / Accepted: 3 April 2023 / Published: 6 April 2023
(This article belongs to the Special Issue Advance in Reproductive Biology and Related Diseases)
Browse Figures
Versions Notes

Abstract

:
One of the most critical steps in mammalian reproduction is implantation. Embryos with an impaired capacity for embryo–maternal crosstalk are thought to have a reduced potential for implantation. One agent of embryo–maternal communication is extracellular vesicles (EV). EVs are lipid bilayer-bound biological nanoparticles implicated in intercellular communication between many of the known cell types. In the current study, we isolated EVs from trophoblast analogue JAr spheroids and supplemented the EVs with receptive endometrium analogue RL95-2 cells to simulate pre-implantation embryo–maternal dialogue. The transcriptome of the endometrial cells was examined at 30 min, 4 h and 48 h intervals using Oxford Nanopore® technology. At the time points, 30 min, 4 h and 48 h, the endometrial cells showed a significantly altered transcriptome. It seems trophoblast EVs induce a swift and drastic effect on the endometrial transcriptome. The effect peaks at around 4 h of EV supplementation, indicating a generalized effect on cell physiology. Alterations are especially apparent in biological pathways critical to embryonic implantation, such as extracellular matrix–receptor interactions and cytokine–receptor interactions. These observations can be helpful in elucidating the dynamics of embryo–maternal communication in the pre-implantation period.

1. Introduction

One of the most critical steps in mammalian reproduction is implantation. During the process, the embryo first adheres to the luminal epithelium of the uterus and then invades the underlying stroma. Implantation is also considered the bottleneck step in human-assisted reproduction, where two out of three of the failed attempts are caused by a failure to implant [1]. Despite tremendous advances in the field of assisted reproduction technologies (ART), the success rate of embryo transfer remains less than 50%, mostly due to the embryo’s inability to implant [2]. Deciphering the causes for this unique situation is a subject of vigorous scientific inquiry, given that the morphological quality of the embryo does not correlate completely with the rate of implantation or successful birth, even in cases where the uterine environment is optimal for implantation [3].
The endometrium undergoes multiple physiological and morphological changes to be receptive to incoming implantation during the short period known as the window of implantation (WOI) [4]. These changes extend to a complete re-organisation of the uterine immune micro-environment. Conventional wisdom suggests that the endometrium and embryo are both passive or “quiet” and all the events leading up to and during the implantation and invasion of the endometrium by the embryo are regulated by external factors such as hormones and cytokines. This is, in fact, proven scientific theory and is supported by multitudinous reports. However, in recent decades, evidence has started to point to another component of the regulatory mechanism of embryo implantation, embryo–maternal communication [5].
This hypothesis posits that the embryo and the endometrium do not passively follow the directives of the endocrine system but actively participate in the process by communicating with each other using one or more signalling methods [6]. Although endometrial regulation by endocrine agents such as oestrogen and progesterone is well known, it does not represent the complete picture of endometrial receptivity regulation. Another aspect of the regulation of implantation is the action of cytokines, which function in an autocrine or paracrine manner during the process of implantation [7]. For instance, leukaemia inhibitory factor (LIF) is a mediator of oestrogen activity and a member of the interleukin-6 family [8]. LIF is also reported to stimulate and regulate stromal proliferation via the epidermal growth factor (EGF) signalling pathway [9]. In addition, heparin-binding epidermal growth factor-like growth factor (HB-EGF) is another crucial molecule in the implantation process [10]. Functioning in a juxtracrine manner, HB-EGF is a transmembrane protein that releases a soluble growth factor that reacts with the EGF family receptors on the trophoblast surface to regulate adhesion [11,12,13].
Another less well-studied facet of embryo–maternal communication is the communication that is hypothesized to utilize extracellular vesicles (EV) as the mode of signal transfer. EVs are defined as membrane-bound, nano-sized particles produced by living cells and released into the extracellular space [14,15]. All of the known cell types across the biological kingdoms are reported to produce EVs [16]. Many biological fluids, such as serum, cerebrospinal fluid and seminal plasma, are known to contain EVs [17]. One of the primary functions of EVs is to transport biologically significant cargo between cells. The transfer is known to occur between cells of the same individual, between individuals and even cross-kingdom. The cargo encompassed in EVs is extremely varied and heterogeneous, even among EVs extracted from purified cell cultures [14].
The involvement of EVs in embryo–maternal communication, especially in the crucial stage of embryo implantation, is not a well-studied field. However, multiple studies report on the involvement of EVs in other aspects of mammalian reproduction, from gamete production and maturation to foetal–maternal communication [18,19]. Our own studies have previously shown significant EV-mediated embryo–maternal communication and the significant effects that such communications induce in endometrial cells [20,21,22,23,24,25,26]. We have termed these communications “the first language” since they are the first communication between the progenitor and the progeny. We have shown evidence that EVs take part in the “first language” via trophoblast EV-induced transcriptomic alteration in the endometrium, increasing the receptivity of the endometrium. However, the exact mechanisms that EVs utilize to induce these physiological effects are still unclear. In the current study, we endeavour to investigate the transcriptomic effects induced by trophoblast EVs in the endometrial cells as a function of time using mRNA sequencing.

2. Results

The current study was carried out to investigate the effects induced by trophoblast analogue 3D spheroid-derived EVs on endometrial analogue cell transcriptomes as a function of time. Since we have demonstrated that trophoblast spheroid-derived EVs can alter the endometrial transcriptome with 24 h incubation in our previous communications, we have selected 30 min, 4 h and 48 h incubation times for this study. After each incubation, the endometrial RNA was collected and sequenced to determine the effects of trophoblast spheroid-derived EVs on the endometrial transcriptome.

2.1. Trophoblast-Derived EVs Altered Endometrial Transcriptome Significantly in Each Tested Time Point

The transcriptomes of endometrial cells treated with trophoblast-derived EVs were significantly altered at all the observed time points. The co-incubation time appeared to have a significant effect on the endometrial cellular transcriptome based on the principal component analysis (PCA) plot (Figure 1A) and the heatmap (Figure 1B), where the variability among the groups (30 min incubation, 4 h incubation and 48 h incubation, and the untreated control) are clearly visible. The untreated control (0 h) was used after observing the stable nature of the untreated RL 95-2 transcriptome as a function of time in cell culture (Supplementary File S1). No significant intra-group variation is apparent. The different altered transcriptomic profiles are presented in the volcano plots: 30 min treatment vs. control (Figure 1C), 4 h treatment vs. control (Figure 1D) and 48 h treatment vs. control (Figure 1E). In the 30 min vs. control comparison, there were 455 upregulated genes and 472 downregulated genes. In the 4 h vs. control comparison, there were 315 upregulated genes and 373 downregulated genes. In the 48 h vs. control comparison, there were 658 upregulated and 625 downregulated genes. In the 4 h vs. 30 min comparison, 272 genes were upregulated and 203 genes were downregulated (Figure 1F). In the 48 h vs. 30 min comparison, 53 genes were upregulated and 60 genes were downregulated (Figure 1G). In the 4 h vs. 48 h comparison, 314 genes were upregulated and 214 genes were downregulated (Figure 1H). A gene was deemed upregulated or downregulated when the log2 fold change (logFC) was more than +1 or less than −1, respectively. The only value of significance was the FDR (false discovery rate) < 0.05. Differentially expressed genes are presented in Supplementary File S2.

2.2. Different Sets of Genes Showed Altered Expression at Each Observed Time Point

Observed transcriptomic alterations in endometrial cells treated with trophoblast spheroid-derived EVs were caused by different sets of genes at each observed time point. In upregulated genes (Figure 2A), only 61 genes (13.5%) were upregulated at all three time points. In downregulated genes (Figure 2B), only 220 genes, or 15.6% of the total genes, were downregulated in all three conditions. When considering 4 h vs. 30 min and 48 h vs. 30 min comparisons, there is a slightly higher number of jointly upregulated (60 genes; 26.7%) and downregulated (145 genes; 25.3%) genes. A gene was deemed upregulated or downregulated when the log2 fold change (logFC) was more than +1 or less than −1, respectively.

2.3. Pathways Crucial in Embryo Implantation Were Enriched in a Specific Pattern along the Three Observed Time Points

Gene set enrichment analysis revealed that a group of pathways associated with the reorganization of the extracellular matrix (Figure 3A,C) and signalling with the epidermal growth factor receptor (Figure 3B,D) were significantly enriched in all three observed time points. The significance of enrichment was assumed when the adjusted p-value < 0.05. Interestingly, the normalized enrichment scores (NES), which denote the level of enrichment in each pathway in terms of the fold changes in the participating genes, followed a specific pattern in both groups of the considered pathways. The highest NES were observed at the 4 h time point, implying a greater amount of specific transcriptomic alterations. At the 48 h time point, the pathways showed a clear return to the base levels. All of the enriched pathways are presented in Supplementary File S3.

2.4. The Expressions of Key Genes in the Selected Pathways Were Significantly Altered Due to the Influence of EVs

The expressions of several genes that are reported to play key roles in either ECM reorganization or EGFR-based signalling were quantified using qPCR. In most cases, significant upregulations were observed in the endometrial cells incubated for 4 h with trophoblast-derived EVs compared to untreated controls. Interestingly, there were no significant upregulations or downregulations in the 30 min incubation group or the 48 h incubation group, either compared to untreated controls or each other (Figure 4). Five housekeeping genes were used to investigate the stability of the endometrial transcriptome during the period of incubation. There were no significant alterations to the expressions of the housekeeping genes. The measurement of significance was p-value < 0.05 from a single factor ANOVA.

3. Discussion

In our previous communications, using the same in vitro system as an analogue for the embryo–maternal microenvironment during the pre-implantation period, we observed drastic alterations to the endometrial transcriptome induced by the trophoblast-derived EVs. These observations implied that EVs play a major role in embryo–maternal communication. When considering the pathways enriched in endometrial cells treated by EVs, we were able to conclude that EV-mediated embryo–maternal communication increases the endometrial receptivity, priming the endometrium for the incoming embryo [22,23]. However, the mechanism of the communication was not apparent from one snapshot of the endometrial transcriptome taken at 24 h of EV–cellular co-incubation. The current study was carried out to create a more holistic image of the dynamics of EV-induced transcriptomic alterations in endometrial cells.
Significant alterations to the endometrial transcriptomes were observed in all three tested time points. Zero-hour untreated controls were used to compare the transcriptomic alterations at each time point after confirming that the RL95-2 transcriptome remains stable during 48 h of cell culture (Supplementary File S1). This observation is especially interesting since the first time point was a mere 30 min post-EV supplementation. The extremely short amount of time required to induce the transcriptomic alterations hints at a very quick mechanism of action in terms of signal transfer.
There are a few methods thought to be used by EVs in intercellular communications. The most tested hypotheses include EV-based communication via EV cargo transfer. EVs are known to carry a large variety of biologically active molecules such as proteins, lipids and RNA that are thought to be released into the cytoplasm of the target cell and bring about the physiological effect [27]. This premise calls for a method of EV uptake and cargo release. It is hypothesized that EVs either fuse directly to the plasma membrane of the target cell and release the cargo or the target cell actively uptakes the EVs, which are subsequently degraded to release the cargo [28,29,30].
Given the time taken for EV uptake, vesicle degradation and cargo release [30,31,32], plus the time that it takes to induce the transcriptomic alteration in the endometrial cells, it is unlikely that any mechanism that involves cargo release is responsible for the majority of the observed changes in the endometrial transcriptome at the 30 min time point. However, there is a higher possibility that EV-based communications that involve EV uptake and cargo release are responsible for at least a part of the endometrial transcriptome alterations observed at the 4 h and 48 h time points. The amount of commonly expressed transcripts (either upregulated or downregulated) at all time points is relatively low, implying different mechanisms are responsible for the changes in the transcriptome in each occurrence.
In our previous studies, we tried to correlate the embryonic miRNA carried in the EVs to the expression of their targets in the endometrium, working on the premise that EV miRNA is released to the cytosol of the receptor cell and downregulates their target mRNA. However, we were only able to correlate a minority of transcriptomic alterations (5–10%) to the actions of miRNA carried in EVs as cargo [22,23]. This observation and the short time taken to induce the transcriptomic effects, coupled with the reports that claim EVs are a very poor source of functional RNA and the very low percentage (1–2%) [30,33,34] of EVs whose cargo is released to the cytosol, leads us to deduce that cargo carried in EVs would be responsible for a minority of the transcriptomic effects observed in endometrial cells treated with trophoblast-derived EVs.
There are speculations about an intercellular signalling method involving EVs that does not involve EV uptake or cargo delivery. EVs are thought to bind to the plasma membrane using specific receptor–ligand complexes and initiate a signalling cascade that achieves the relevant physiological effect in the target cell [35,36,37,38,39,40]. Considering the rapidity of the trophoblast EVs’ effect on endometrial cells, it can be proposed that the initial transcriptomic alterations are achieved using a membrane-bound signalling system involving EVs. The signalling can also take place through the endosomal membrane once the EVs are taken up by the target cell via any of the endocytosis pathways.
Considering the pathways that are significantly enriched by the altered transcriptomes at each time point, gene set enrichment analysis revealed the over-representation of two groups of pathways: extracellular matrix (ECM) reorganization and pathways associated with epidermal growth factor receptor (EGFR)-related signalling. Both of these groups have been reported to be critically important in the implantation process. ECM remodelling is a critical change that the endometrial cells undergo in the process of decidualization. The receptors for the adhesion molecules are overexpressed on the apical surface and the pinopodes to facilitate trophoblast adhesion [41,42,43]. Major subunits of the ECM pathway, such as laminin interactions (R-HSA-3000157) [44], integrin cell surface interactions (R-HSA-216083) [45,46,47] and non-integrin membrane–ECM interactions (R-HSA-3000171) [48], are also significantly enriched and are all known to be implicated in the endometrial modifications in the WOI.
The epidermal growth factor receptor (EGFR) is a transmembrane protein belonging to the epidermal growth factor family, a subfamily of four closely related tyrosine kinases (EGFR (ErbB-1), HER2/neu (ErbB-2), HER3 (ErbB-3) and HER4 (ErbB-4)) [49]. They are reported to be critical for a successful early pregnancy, especially in regulating decidualization. EGFR is active during the window of implantation, and the impaired activation of such has been linked to many birth defects, such as intrauterine growth restriction and low birth weight [50,51,52,53]. Growth factors such as EGF and transforming growth factor alpha (TGFα) are implicated in mitogenic action, and the ablation of EGFR leads to severe subfertility in mice [49]. The cause of this subfertility is reported to be due to a failure in the maintenance and progression of the decidua [49,54].
Apart from being simply significantly enriched, the pattern exhibited by the level of enrichment, denoted by the normalized enrichment score (NES) in pathways associated with ECM- and EGFR-based signalling, is of special interest. Most pathways were enriched only mildly at the 30 min time point (NES between 0.3 and 0.9), implying a mild physiological effect. At the 4 h time point, a more drastic enrichment was observed, with enrichment scores of more than 1.5 in many cases. The cellular physiology seemed to return to its base physiology (denoted by 0 NES) at the 48 h time point. In our previous experiments, with 24 h co-incubation, we observed the NES to be around 1.3 to 1.6 for ECM-related pathways.
Another facet of the observed significant effect of trophoblast cell-derived EVs on the endometrial transcriptome can be perceived by contrasting the transcriptomes between the time points. The 4 h vs. 30 min and the 48 h vs. 30 min comparisons show that pathways are very similarly enriched on both occasions (Figure 3C,D). This can be interpreted as both the 4 h and 48 h transcriptomes being similarly altered compared to the transcriptome of the 30 min time point. Along with the fact that 4 h vs. 48 h shows very low amounts of pathway enrichment, these observations reaffirm the previous elucidation that there is a significant enrichment of critical pathways in embryo implantation as early as 30 min of trophoblast EV introduction, which increases by 4 h and remains in a slightly elevated status at the 48 h time point. Some pathways are even elevated in 48 h compared to 4 h, indicating lasting effects on endometrial transcriptome well beyond the half-life of EVs in cell culture conditions.
Quantitative real-time PCR results of 25 selected genes representing both pathway clusters showed a similar trend in expression profiles compared to the sequencing data. Most of the genes were very slightly differentially expressed from 0 h to 30 min, there were significant upregulations in most genes at the 4 hr time point and at 48 h the expression returned generally to the levels of 30 min and 0 hr. None of the genes showed any significant differential expression in the 48 h vs. 30 min comparison, indicating both the drastic nature of the transcriptomic changes seen at the 4 h time point and the fact that the endometrial transcriptome returned to a semblance of the base level after 48 hrs. All the selected genes were reported to regulate some critical function in embryo implantation, either as individual genes or as part of a pathway [49,55,56,57,58,59,60,61]. Interestingly, none of the genes that were upregulated or downregulated in all three time points nor the pathways they enrich were mentioned regarding human embryo implantation in previous publications, to the best of our knowledge (Supplementary File S4).
Taken as a whole, these observations suggest an underlying timeline of EV-induced transcriptomic changes, the mechanisms of which can be considered in the following manner. The minor pathway enrichments at the 30 min time point are the results of transcriptomic changes brought about by the swifter membrane-bound signalling accomplished by molecules on the EV surface binding to the receptors on the endometrial cellular plasma membrane. Then, at the 4 h time point, both the surface signalling and the effects of EV uptake and cargo release are seen as a drastic upregulation of implantation-critical pathways priming the endometrium for the incoming embryo. This heightened level of receptivity continues until the implantation allies with a constant source of embryonic EVs, such as an embryo that is endeavouring to implant. However, since we only introduced EVs once to the system, the transcriptomic changes began to plateau around 24 h post-EV treatment, and the endometrium returned to a slightly elevated level in 48 h. These observations are useful not only in elucidating the physiological effects of EVs but also in determining the half-life of EVs once introduced to the target cell culture.
The exact mechanisms/methods utilized by trophoblast-derived EVs to bring about the observed transcriptomic alterations can only be speculated based on the results of this study. Given the rapidity of the change in the endometrial transcriptome, the theory that posits that EVs contain membrane-bound signalling molecules that can initiate a cascade of events that eventually lead to the observed transcriptomic alterations would be the most educated guess at the present. However, there is very minimal evidence of any such occurrence in any system thought to be utilizing EV-based intercellular communications. Since both pathway clusters can be activated by a multitude of activators, it is not possible to accurately predict the correct molecule/s that may be responsible for the observed effects. Studies into the proteomes of trophoblast-derived EVs might be the next logical step in discovering the elusive agent of EV-mediated embryo–maternal communication.
The in vitro model used to mimic the pre-implantation microenvironment was made from cancer-derived cell lines. However, the model has been in use since the 1970s and is widely regarded as one of the best representations of the pre-implantation microenvironment. Given that using embryos and receptive endometrial tissue for experimentation is highly questionable ethically, the in vitro model was seen as an acceptable alternative.
The sequencing platform Oxford Nanopore®, although designed especially for ultra-long length reads, is an admirably user-friendly and low-cost method to investigate the differential expression of cellular mRNA. The obstacles introduced by the shallow depth of reads were mitigated with multiple wash/reload cycles, shorter read lengths and longer terms of sequencing. Even with extremely short read lengths (when considering the possibilities of Oxford Nanopore®, Littlemore, Oxford, UK), the coverage of the transcriptome was superior, in our hands, to Illumina sequencing. The Oxford Nanopore® platform is much more cost-effective and less labour-intensive compared to sequence-by-synthesis methods.

4. Materials and Methods

4.1. Cell Culture and Trophoblast Spheroid Preparation

All the immortalized cell lines were purchased from the American Type Culture Collection in Teddington, United Kingdom. All the cell lines were cultured in a 37 °C and 5% CO2 environment. The endometrial analogue cell line was prepared using RL 95-2 cells (human adenosquamous carcinoma, ATCC CRL-1671). The culture medium for the RL95-2 cells consisted of Dulbecco’s Modified Eagle Medium (DMEM 12-604F, Lonza, Verviers, Belgium), 5 μg/mL human recombinant insulin (Gibco, Invitrogen, Denmark), 1% Penicillin-Streptomycin (P/S, Gibco™ 15140122, Bleiswijk, The Netherlands), 1% glutamine (Sigma, 59202C, Saint Louis, MO, USA) and 10% foetal bovine serum (FBS, Gibco™, 10500064). Cells were cultured at 37 °C and in 5% CO2.
The trophoblast analogue, human first-trimester choriocarcinoma cell line JAr (HTB-144 ™, Teddington, UK) cells, were cultured in RPMI 1640 (Gibco, Inchinnan, UK), 1% L-glutamine, 10% FBS, 1% penicillin/streptomycin. When the cells reached 80% confluency, they were collected using the trypsin-EDTA method, re-suspended in 5 mL of complete medium (1 × 106 cells/mL concentration) and shaken overnight using a gyratory shaker (295 rpm, 37 °C, 5% CO2). The resulting spheroids were used as the analogues for the hatched trophoblast.

4.2. Depleting Extracellular Vesicles from FBS Using Ultrafiltration

In 2018, Kornilov et al. described a method of depleting EVs from FBS without compromising the material makeup of the reagent. We used the ultrafiltration method with modifications to deplete EVs from FBS that would later be used in the culture medium. FBS was centrifuged in a 100 kDa cut-off ultrafiltration unit (100 kDa, MERCK KGAA, Darmstadt, Germany) for 55 min at 3000× g. The filtrate was found to be up to 95% EV depleted [62]. Complete medium with EV-depleted FBS was prepared when trophoblast EVs were supplemented to endometrial cells.

4.3. Size Exclusion Chromatography for EV Enrichment

Size exclusion chromatography was used for trophoblast spheroid EV enrichment. Sepharose beads (cross-linked 4% agarose matrix of 90 µm beads) were used (Sepharose 4 fast flow™, GE HealthCare Bio-Sciences AB, Uppsala, Sweden) in a 10 cm gravity flow column. A 500 µL sample was added to the top of the column and fractions 7 to 10 of 500 µL fractions were collected in DPBS. Collected fractions were then further concentrated using 10 kDa ultrafiltration units.

4.4. RNA Extraction and Quality Control

RNA was extracted using the phenol–chloroform method with TRIzol reagent acting as the chaotropic agent (TRIzol® reagent; Invitrogen, Inchinnan, UK). Concentrated EVs or adhered cells were used as the RNA source. In the case of EVs, the TRIzol reagent was added directly to the solution. In the case of adherent cells, the conditioned medium was removed, and then the TRIzol reagent was added directly to the cell layer. The samples were then left at room temperature for 10 min to completely dissociate the nucleoprotein complexes, and 300 µL of chloroform for 1 mL of TRIzol was added. The samples were then vigorously shaken for 15 s and centrifuged at 12,000× g for 15 min at 4 °C. When the phases were separated, the aqueous phase containing the RNA was removed, and the RNA was precipitated using 500 µL of isopropanol at room temperature for 20 min. To increase the efficiency of RNA extraction, 10 µg of glycogen (UltraPure™ Glycogen, 10814-010, Thermo Fisher Scientific, Bleiswijk, The Netherlands) was added to the mixture. Once the precipitation was complete, the samples were centrifuged at 18,000× g for 30 min at RT to pellet the RNA. The RNA pellets were washed three times in 70% ethanol and eluted in 20 µL of nuclease-free water. Quantity and the quality of RNA were measured using the Qubit HS RNA kit and the Agilent Pico 6000 kits (Agilent technologies, Santa Clara, CA, USA).

4.5. Library Preparation and RNAseq

Oxford Nanopore® platform was used for RNA sequencing due to its utility and affordability. We used the MinION device with R9.4.1 (FLO-MIN106D) flow cells and the cDNA PCR Barcoding kit (SQK-PCB-109) with 100 ng of input RNA. The initial pore count of the flow cell was 897. Library was prepared to produce 500 bases per read (487 ± 56). Sequencing time was 72 h (with wash/refill cycles every 24 hrs). Twelve samples were run through the flow cell with indexing. Demultiplexing, quality control and adapter trimming were achieved through the Oxford Nanopore® proprietary software.
The 72 h sequencing run produced 16.32 million reads (average reads per sample: 1.272 M ± 0.183 M). The proprietary file format FAST5 was basecalled to FASTQ files using the guppy algorithm. The reads were aligned to human genome assembly version GRCh38.108 using Homo_sapiens.GRCh38.108.gtf annotations and GraphMap algorithm. Read counts at the gene level were obtained by using featureCounts [63]. Genes with at least 10 counts for all the samples in at least one of the experimental groups were retained in the analysis for subsequent differential expression testing. Alignment data is presented in the Table 1.

4.6. Differential Gene Expression Analysis and Network Analysis

Differential expression (DE) analysis was carried out in R version 3.6.1, using the edgeR package version 3.26.8 [64]. Tagwise dispersion estimates were obtained based on the trended dispersions, and statistical comparisons were performed using a generalized linear model followed by likelihood ratio tests, also accounting for the experiment batch. We considered the differential expression of genes with a false discovery rate (FDR) ≤ 0.05 to be statistically significant.
Principal components were calculated using prcomp function from the Stats package and visualized using the ggplot2 package [65]. The pheatmap package [66] was used for heatmap visualization with hierarchical clustering based on Euclidean distance.
Gene set enrichment analysis (GSEA) and pathway over-representation analysis were conducted using the ReactomePA package [67] and Reactome Pathway database annotations [67]. GSEA was used for full gene lists obtained from DE analysis that were ranked by −log10FDR × log2FC, where FDR denotes adjusted p-values and FC the fold change.

4.7. Quantitative Real-Time PCR with Absolute Quantitation

A mixture of random hexamer and oligo (dT) primers was used for cDNA synthesis from enriched RNA (SuperScript® VILO™ cDNA synthesis kit, 11754050). cDNA products were amplified in EvaGreen assay system (Solis BioDyne, Tartu, Estonia) with the following program: 95 °C for 15 min, followed by 40 cycles at 95 °C for 20 s, 60 °C for 20 s, and 72 °C for 20 s. For melting curve analysis, the fluorescence signals were collected continuously from 65 °C to 95 °C at 0.05 °C per second.
For absolute quantitation, spike-in and normalizing of transcripts, 100 bp region from Isopenicillin N-CoA synthetase gene was used (Biomer.net company, Ulm/Donau, Germany, molecular weight: 32,239 g/mol, 100 pmol/µL). Synthetic RNA was serially diluted 20 times. For the first serial dilution, 1 µL of synthetic RNA was added to 39 µL RNase-free water to final concentration of 2.5 µM. Serial dilutions were prepared with a dilution factor of 4×. Serial dilutions were reverse-transcribed and amplified using real-time PCR, and the cycle threshold (Ct) values of dilutions were plotted against the copy number of the transcripts. Exponential calibration curve was fitted. In parallel, 1 µL of synthetic transcript was added to the sample during TRIzol RNA extraction, and then the Ct of synthetic RNA in this sample was assayed to calculate the RNA extraction efficiency and normalizing factor [22,68].

4.8. Experimental Plan

Trophoblast spheroids were prepared using the trophoblast analogue human choriocarcinoma JAr cell line. EVs were isolated from JAr cells cultured in EV-depleted media for 24 h. The EVs were co-incubated for 30 min (n = 3), 4 h (n = 3) and 48 h (n = 3) with a monolayer of RL95-2 cells. The ratio between the EVs and cells was 50:1. After the incubation, the RL95-2 cells were collected, and the cellular RNA was sequenced using Oxford Nanopore® platform. Control samples were prepared using untreated RL 95-2 cells.

5. Conclusions

Taken as a whole, the observations on the dynamics of the transcriptomic changes in endometrial cells induced by trophoblast cells indicate that: (A) the transcriptome is altered very rapidly when in contact with trophoblast EVs, indicating a probable mechanism/signalling method that is implemented rather quicker than the more familiar “EV uptake → cargo delivery → physiological effect achieved by a cargo” pathway. (B) Pathways that are critical for embryo implantation are significantly enriched in all three observed time points, indicating a specific physiological purpose of EV-mediated embryo–maternal communication. (C) The pathways are enriched in a pattern that indicates the influence of EVs as a function of time, allowing the deduction of the half-life of EVs in a cell culture system. The exact EV-related mechanisms or substances that induce these transcriptomic changes are hard to speculate with the current knowledge. Further investigations into EV surface molecules may reveal the responsible underlying mechanism. Understanding the timeline of EV-mediated embryo–maternal communication will be beneficial in developing diagnostic and therapeutic applications for embryo-derived EVs.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24076811/s1.

Author Contributions

Conceptualization, K.G., P.R.H. and A.F.; Methodology, K.G., P.R.H. and A.F.; Formal Analysis, K.G.; Investigation, K.G.; Resources, A.F.; Data Curation, K.G.; Writing—Original Draft Preparation, K.G.; Writing—Review and Editing, K.G., P.R.H. and A.F.; Visualization, K.G.; Supervision, A.F.; Project Administration, A.F.; Funding Acquisition, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by European Union’s Horizon 2020 research and innovation program under grant agreement No 857418 (COMBIVET) and agreement No: 101079349 (OH-Boost).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in Sequence Read Archive (SRA). Reference number PRJNA948686.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jarvis, G.E. Early embryo mortality in natural human reproduction: What the data say. F1000Research 2016, 5, 2765. [Google Scholar] [CrossRef] [PubMed]
  2. Apter, S.; Balaban, B.; Campbell, A.; Catt, J.; Coticchio, G.; de los Santos, M.J.; Debrock, S.; Ebner, T. The Vienna consensus: Report of an expert meeting on the development of ART laboratory performance indicators. Reprod. Biomed. Online 2017, 35, 494–510. [Google Scholar] [CrossRef] [Green Version]
  3. Wilcox, A.J.; Baird, D.D.; Weinberg, C.R. Time of Implantation of the Conceptus and Loss of Pregnancy. N. Engl. J. Med. 1999, 340, 1796–1799. [Google Scholar] [CrossRef] [PubMed]
  4. Ashary, N.; Tiwari, A.; Modi, D. Embryo Implantation: War in times of love. Endocrinology 2018, 159, 1188–1198. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Modi, D.; Bhartiya, P. Physiology of embryo-endometrial cross talk. Biomed. Res. J. 2015, 2, 83–104. [Google Scholar] [CrossRef]
  6. Bhagwat, S.R.; Chandrashekar, D.S.; Kakar, R.; Davuluri, S.; Bajpai, A.K.; Nayak, S.; Bhutada, S.; Acharya, K.; Sachdeva, G. Endometrial Receptivity: A Revisit to Functional Genomics Studies on Human Endometrium and Creation of HGEx-ERdb. PLoS ONE 2013, 8, e58419. [Google Scholar] [CrossRef]
  7. McEwan, M.; Lins, R.J.; Munro, S.K.; Vincent, Z.L.; Ponnampalam, A.P.; Mitchell, M.D. Cytokine regulation during the formation of the fetal-maternal interface: Focus on cell-cell adhesion and remodelling of the extra-cellular matrix. Cytokine Growth Factor Rev. 2009, 20, 241–249. [Google Scholar] [CrossRef]
  8. Massimiani, M.; Lacconi, V.; La Civita, F.; Ticconi, C.; Rago, R.; Campagnolo, L. Molecular signaling regulating endometrium–blastocyst crosstalk. Int. J. Mol. Sci. 2020, 21, 23. [Google Scholar] [CrossRef] [Green Version]
  9. Hantak, A.M.; Bagchi, I.C.; Bagchi, M.K. Role of uterine stromal-epithelial crosstalk in embryo implantation. Int. J. Dev. Biol. 2014, 58, 139–146. [Google Scholar] [CrossRef] [Green Version]
  10. Cha, J.; Sun, X.; Dey, S.K. Mechanisms of implantation: Strategies for successful pregnancy. Nat. Med. 2012, 18, 1754–1767. [Google Scholar] [CrossRef]
  11. Raab, G.; Kover, K.; Paria, B.C.; Dey, S.K.; Ezzell, R.M.; Klagsbrun, M. Mouse preimplantation blastocysts adhere to cells expressing the transmembrane form of heparin-binding EGF-like growth factor. Development 1996, 122, 637–645. [Google Scholar] [CrossRef] [PubMed]
  12. Stavreus-Evers, A.; Aghajanova, L.; Brismar, H.; Eriksson, H.; Landgren, B.M.; Hovatta, O. Co-existence of heparin-binding epidermal growth factor-like growth factor and pinopodes in human endometrium at the time of implantation. Mol. Hum. Reprod. 2002, 8, 765–769. [Google Scholar] [CrossRef] [Green Version]
  13. Paria, B.C.; Elenius, K.; Klagsbrun, M.; Dey, S.K. Heparin-binding EGF-like growth factor interacts with mouse blastocysts independently of ErbB1: A possible role for heparan sulfate proteoglycans and ErbB4 in blastocyst implantation. Development 1999, 126, 1997–2005. [Google Scholar] [CrossRef] [PubMed]
  14. Raposo, G.; Stoorvogel, W. Extracellular vesicles: Exosomes, microvesicles, and friends. J. Cell Biol. 2013, 200, 373–383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Février, B.; Raposo, G. Exosomes: Endosomal-derived vesicles shipping extracellular messages. Curr. Opin. Cell Biol. 2004, 16, 415–421. [Google Scholar] [CrossRef] [PubMed]
  16. Ng, Y.H.; Rome, S.; Jalabert, A.; Forterre, A.; Singh, H.; Hincks, C.L.; Salamonsen, L.A. Endometrial exosomes/microvesicles in the uterine microenvironment: A new paradigm for embryo-endometrial cross talk at implantation. PLoS ONE 2013, 8, e58502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Smijs, T.G.; Pavel, S.; O’Farrell, N.; Houlton, A.; Horrocks, B.R.; Yasui, T.; Yanagida, T.; Ito, S.; Konakade, Y.; Takeshita, D.; et al. Unveiling massive numbers of cancer-related urinary-microRNA candidates via nanowires. Sci. Adv. 2017, 1, 95–112. [Google Scholar] [CrossRef]
  18. Aleksejeva, E.; Zarovni, N.; Dissanayake, K.; Godakumara, K.; Vigano, P.; Fazeli, A.; Jaakma, Ü.; Salumets, A. Extracellular vesicle research in reproductive science— Paving the way for clinical achievements. Biol. Reprod. 2022, 106, 408–424. [Google Scholar] [CrossRef]
  19. Hart, A.R.; Khan, N.L.A.; Godakumara, K.; Dissanayake, K.; Piibor, J.; Muhandiram, S.; Eapen, S.; Heath, P.R.; Fazeli, A. The role of extracellular vesicles in endometrial receptivity and their potential in reproductive therapeutics and diagnosis. Reprod. Biol. 2022, 22, 100645. [Google Scholar] [CrossRef]
  20. Godakumara, K.; Dissanayake, K.; Hasan, M.M.; Kodithuwakku, S.P.; Fazeli, A. Role of extracellular vesicles in intercellular communication during reproduction. Reprod. Domest. Anim. 2022, 57, 14–21. [Google Scholar] [CrossRef]
  21. Hasan, M.M.; Reshi, Q.U.A.; Lättekivi, F.; Viil, J.; Godakumara, K.; Dissanayake, K.; Andronowska, A.; Jaakma, Ü.; Fazeli, A. Bovine Follicular Fluid Derived Extracellular Vesicles Modulate the Viability, Capacitation and Acrosome Reaction of Bull Spermatozoa. Biology 2021, 10, 1154. [Google Scholar] [CrossRef] [PubMed]
  22. Es-Haghi, M.; Godakumara, K.; Häling, A.; Lättekivi, F.; Lavrits, A.; Viil, J.; Andronowska, A.; Nafee, T.; James, V.; Jaakma, Ü.; et al. Specific trophoblast transcripts transferred by extracellular vesicles affect gene expression in endometrial epithelial cells and may have a role in embryo-maternal crosstalk. Cell Commun. Signal. 2019, 17, 1–18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Godakumara, K.; Ord, J.; Lättekivi, F.; Dissanayake, K.; Viil, J.; Boggavarapu, N.R.; Faridani, O.R.; Jääger, K.; Velthut-Meikas, A.; Jaakmaet, Ü.; et al. Trophoblast derived extracellular vesicles specifically alter the transcriptome of endometrial cells and may constitute a critical component of embryo-maternal communication. Reprod. Biol. Endocrinol. 2021, 19, 1–14. [Google Scholar] [CrossRef] [PubMed]
  24. Reshi, Q.U.A.; Viil, J.; Ord, J.; Lättekivi, F.; Godakumara, K.; Hasan, M.M.; Nõmm, M.; Jääger, K.; Velthut-Meikas, A.; Jaakmaet, Ü.; et al. Spermatozoa induce transcriptomic alterations in bovine oviductal epithelial cells prior to initial contact. J. Cell Commun. Signal 2020, 14, 439–451. [Google Scholar] [CrossRef]
  25. Dissanayake, K.; Nõmm, M.; Lättekivi, F.; Ord, J.; Ressaissi, Y.; Godakumara, K.; Reshi, Q.U.A.; Viil, J.; Jääger, K.; Velthut-Meikas, A.; et al. Oviduct as a sensor of embryo quality: Deciphering the extracellular vesicle (EV)-mediated embryo-maternal dialogue. J. Mol. Med. 2021, 99, 685–697. [Google Scholar] [CrossRef] [PubMed]
  26. Dissanayake, K.; Nõmm, M.; Lättekivi, F.; Ressaissi, Y.; Godakumara, K.; Lavrits, A.; Midekessa, G.; Vill, J.; Bæk, R.; Jørgensen, M.M.; et al. Individually cultured bovine embryos produce extracellular vesicles that have the potential to be used as non-invasive embryo quality markers. Theriogenology 2020, 149, 104–116. [Google Scholar] [CrossRef]
  27. Gangoda, L.; Boukouris, S.; Liem, M.; Kalra, H.; Mathivanan, S. Extracellular vesicles including exosomes are mediators of signal transduction: Are they protective or pathogenic? Proteomics 2015, 15, 260. [Google Scholar] [CrossRef] [Green Version]
  28. Théry, C.; Zitvogel, L.; Amigorena, S. Exosomes: Composition, biogenesis and function. Nat. Rev. Immunol. 2002, 2, 569–579. [Google Scholar] [CrossRef]
  29. Mathieu, M.; Martin-Jaular, L.; Lavieu, G.; Théry, C. Specificities of secretion and uptake of exosomes and other extracellular vesicles for cell-to-cell communication. Nat. Cell Biol. 2019, 21, 9–17. [Google Scholar] [CrossRef]
  30. Bonsergent, E.; Grisard, E.; Buchrieser, J.; Schwartz, O.; Théry, C.; Lavieu, G. Quantitative characterization of extracellular vesicle uptake and content delivery within mammalian cells. Nat. Commun. 2021, 12, 1864. [Google Scholar] [CrossRef] [PubMed]
  31. Feng, D.; Zhao, W.L.; Ye, Y.Y.; Bai, X.C.; Liu, R.Q.; Chang, L.F.; Zhou, Q.; Sui, S.-F. Cellular internalization of exosomes occurs through phagocytosis. Traffic 2010, 11, 675–687. [Google Scholar] [CrossRef] [PubMed]
  32. Fabbri, M.; Paone, A.; Calore, F.; Galli, R.; Gaudio, E.; Santhanam, R.; Lovat, F.; Fadda, P.; Mao, C.; Nuovo, G.J.; et al. MicroRNAs bind to Toll-like receptors to induce prometastatic inflammatory response. Proc. Natl. Acad. Sci. USA 2012, 109, E2110–E2116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Bonsergent, E.; Lavieu, G. Content release of extracellular vesicles in a cell-free extract. FEBS Lett. 2019, 593, 1983–1992. [Google Scholar] [CrossRef] [PubMed]
  34. Chevillet, J.R.; Kang, Q.; Ruf, I.K.; Briggs, H.A.; Vojtech, L.N.; Hughes, S.M.; Arroyo, J.D.; Cheng, H.H.; Meredith, E.K.; Gallichotte, E.N.; et al. Quantitative and stoichiometric analysis of the microRNA content of exosomes. Proc. Natl. Acad. Sci. USA 2014, 111, 14888–14893. [Google Scholar] [CrossRef] [Green Version]
  35. Mulcahy, L.A.; Pink, R.C.; Carter, D.R.F. Routes and mechanisms of extracellular vesicle uptake. J. Extracell. Vesicles 2014, 3, 24641. [Google Scholar] [CrossRef] [Green Version]
  36. Andreu, Z.; Yáñez-Mó, M. Tetraspanins in extracellular vesicle formation and function. Front. Immunol. 2014, 5, 442. [Google Scholar] [CrossRef] [Green Version]
  37. Simeone, P.; Bologna, G.; Lanuti, P.; Pierdomenico, L.; Guagnano, M.T.; Pieragostino, D.; Del Boccio, P.; Vergara, D.; Marchisio, M.; Miscia, S.; et al. Extracellular Vesicles as Signaling Mediators and Disease Biomarkers across Biological Barriers. Int. J. Mol. Sci. 2020, 21, 2514. [Google Scholar] [CrossRef] [Green Version]
  38. Zhdanov, V.P. Ligand-receptor-mediated attachment of lipid vesicles to a supported lipid bilayer. Eur. Biophys. J. 2020, 49, 395. [Google Scholar] [CrossRef]
  39. Staufer, O.; Hernandez Bücher, J.E.; Fichtler, J.; Schröter, M.; Platzman, I.; Spatz, J.P. Vesicle Induced Receptor Sequestration: Mechanisms behind Extracellular Vesicle-Based Protein Signaling. Adv. Sci. 2022, 9, 2200201. [Google Scholar] [CrossRef]
  40. Torres, V.I.; Barrera, D.P.; Varas-Godoy, M.; Arancibia, D.; Inestrosa, N.C. Selective Surface and Intraluminal Localization of Wnt Ligands on Small Extracellular Vesicles Released by HT-22 Hippocampal Neurons. Front. Cell Dev. Biol. 2021, 9, 2887. [Google Scholar] [CrossRef]
  41. Kaloglu, C.; Onarlioglu, B. Extracellular matrix remodelling in rat endometrium during early pregnancy: The role of fibronectin and laminin. Tissue Cell 2010, 42, 301–306. [Google Scholar] [CrossRef]
  42. Brown, J.J.G.; Papaioannou, V.E. Distribution of hyaluronan in the mouse endometrium during the periimplantation period of pregnancy. Differentiation 1992, 52, 61–68. [Google Scholar] [CrossRef] [PubMed]
  43. Johnson, G.A.; Burghardt, R.C.; Bazer, F.W.; Spencer, T.E. Osteopontin: Roles in Implantation and Placentation1. Biol. Reprod. 2003, 69, 1458–1471. [Google Scholar] [CrossRef]
  44. Miner, J.H.; Li, C.; Mudd, J.L.; Go, G.; Sutherland, A.E. Compositional and structural requirements for laminin and basement membranes during mouse embryo implantation and gastrulation. Development 2004, 131, 2247–2256. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Al-Dossary, A.A.; Bathala, P.; Caplan, J.L.; Martin-DeLeon, P.A. Oviductosome-Sperm Membrane Interaction in Cargo Delivery: Detection of fusion and underlying molecular players using three-dimensional super-resolution structured illumination microscopy (SR-SIM). J. Biol. Chem. 2015, 290, 17710–17723. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Reddy, K.V.R.; Mangale, S.S. Integrin receptors: The dynamic modulators of endometrial function. Tissue Cell 2003, 35, 260–273. [Google Scholar] [CrossRef] [PubMed]
  47. Aplin, J.D. Adhesion Molecules and Implantation. In Implantation and Early Development; Cambridge University Press: Cambridge, UK, 2005; pp. 49–60. [Google Scholar] [CrossRef]
  48. Kirn-Safran, C.B.; D’Souza, S.S.; Carson, D.D. Heparan sulfate proteoglycans and their binding proteins in embryo implantation and placentation. Semin. Cell Dev. Biol. 2008, 19, 187–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Large, M.J.; Wetendorf, M.; Lanz, R.B.; Hartig, S.M.; Creighton, C.J.; Mancini, M.A.; Kovanci, E.; Lee, K.-F.; Threadgill, D.W.; Lydon, J.P.; et al. The Epidermal Growth Factor Receptor Critically Regulates Endometrial Function during Early Pregnancy. PLoS Genet. 2014, 10, 1004451. [Google Scholar] [CrossRef] [Green Version]
  50. Gabriel, R.; Alsat, E.; Evian-Brion, D. Alteration of epidermal growth factor receptor in placental membranes of smokers: Relationship with intrauterine growth retardation. Am. J. Obstet. Gynecol. 1994, 170, 1238–1243. [Google Scholar] [CrossRef]
  51. Fujita, Y.; Kurachi, H.; Morishige, K.; Amemiya, K.; Terakawa, N.; Miyake, A.; Tanizawa, O. Decrease in epidermal growth factor receptor and its messenger ribonucleic acid levels in intrauterine growth-retarded and diabetes mellitus-complicated pregnancies. J. Clin. Endocrinol. Metab. 1991, 72, 1340–1345. [Google Scholar] [CrossRef]
  52. Dissanayake, V.H.W.; Tower, C.; Broderick, A.; Stocker, L.J.; Seneviratne, H.R.; Jayasekara, R.W.; Kalsheker, N.; Pipkin, F.B.; Morgan, L. Polymorphism in the epidermal growth factor gene is associated with birthweight in Sinhalese and white Western Europeans. Mol. Hum. Reprod. 2007, 13, 425–429. [Google Scholar] [CrossRef] [PubMed]
  53. Fondacci, C.; Alsat, E.; Gabriel, R.; Blot, P.; Nessmann, C.; Evain-Brion, D. Alterations of human placental epidermal growth factor receptor in intrauterine growth retardation. J. Clin. Investig. 1994, 93, 1149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Tan, Y.; Li, M.; Cox, S.; Davis, M.K.; Tawfik, O.; Paria, B.C.; Das, S.K. HB-EGF directs stromal cell polyploidy and decidualization via cyclin D3 during implantation. Dev. Biol. 2004, 265, 181–195. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Yoshihara, M.; Mizutani, S.; Kato, Y.; Matsumoto, K.; Mizutani, E.; Mizutani, H.; Fujimoto, H.; Osuka, S.; Kajiyama, H. Recent Insights into Human Endometrial Peptidases in Blastocyst Implantation via Shedding of Microvesicles. Int. J. Mol. Sci. 2021, 22, 13479. [Google Scholar] [CrossRef]
  56. Shi, J.W.; Lai, Z.Z.; Yang, H.L.; Yang, S.L.; Wang, C.J.; Ao, D.; Ruan, L.Y.; Shen, H.H.; Zhou, W.J.; Mei, J. Collagen at the maternal-fetal interface in human pregnancy. Int. J. Biol. Sci. 2020, 16, 2220. [Google Scholar] [CrossRef]
  57. Li, T.; Greenblatt, E.M.; Shin, M.E.J.; Brown, T.J.; Chan, C. Endometrial laminin subunit beta-3 expression associates with reproductive outcome in patients with repeated implantation failure. J. Assist. Reprod. Genet. 2021, 38, 1835–1842. [Google Scholar] [CrossRef]
  58. Guo, X.; Yi, H.; Li, T.C.; Wang, Y.; Wang, H.; Chen, X. Role of Vascular Endothelial Growth Factor (VEGF) in Human Embryo Implantation: Clinical Implications. Biomolecules 2021, 11, 253. [Google Scholar] [CrossRef]
  59. Klonisch, T.; Wolf, P.; Hombach-Klonisch, S.; Vogt, S.; Kuechenhoff, A.; Tetens, F.; Fischer, B. Epidermal growth factor-like ligands and erbB genes in the peri-implantation rabbit uterus and blastocyst. Biol. Reprod. 2001, 64, 1835–1844. [Google Scholar] [CrossRef] [Green Version]
  60. Gazor, R.; Eskandari, M.; Sharafshah, A.; Bahadori, M.H.; Golmohammadi, M.G.; Keshavarz, P. Assessment of EGFR Gene Expression Following Vitrification of 2-cell and Blastocyst Mouse Embryos. Avicenna J. Med. Biotechnol 2018, 10, 120. [Google Scholar]
  61. Wang, W.; Taylor, R.N.; Bagchi, I.C.; Bagchi, M.K. Regulation of human endometrial stromal proliferation and differentiation by C/EBPβ involves cyclin E-cdk2 and STAT3. Mol. Endocrinol. 2012, 26, 2016–2030. [Google Scholar] [CrossRef] [Green Version]
  62. Kornilov, R.; Puhka, M.; Mannerström, B.; Hiidenmaa, H.; Peltoniemi, H.; Siljander, P.; Seppänen-Kaijansinkko, R.; Kaur, S. Efficient ultrafiltration-based protocol to deplete extracellular vesicles from fetal bovine serum. J. Extracell. Vesicles 2018, 7, 1422674. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Liao, Y.; Smyth, G.K.; Shi, W. Feature Counts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014, 30, 923–930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Wickham, H. Ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  66. Kolde, R. Pheatmap: Pretty heatmaps. R Package Version 2012, 1, 726. [Google Scholar]
  67. Yu, G.; He, Q.Y. ReactomePA: An R/Bioconductor package for reactome pathway analysis and visualization. Mol. Biosyst. 2016, 12, 477–479. [Google Scholar] [CrossRef]
  68. Wang, J.H.; Jiang, D.; Rao, H.Y.; Zhao, J.M.; Wang, Y.; Wei, L. Absolute quantification of serum microRNA-122 and its correlation with liver inflammation grade and serum alanine aminotransferase in chronic hepatitis C patients. Int. J. Infect. Dis. 2015, 30, e52–e56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 24 06811 g001 550
Figure 1. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. (A) The PCA plot shows the variance between the groups (30 min; green, 4 h; blue, 48 h; orange and control; yellow). (B) Heatmap shows the differential expression of the genes (30 min; green, 4 h; blue, 48 h; purple and control; black). Volcano plots show the expression of individual genes in each comparison: 30 min treatment vs. control (C), 4 h treatment vs. control (D), 48 h treatment vs. control (E) 4 h vs. 30 min (F), 48 h vs. 30 min (G), 4 h vs. 48 h (H). Downregulation (FDR < 0.05, logFC < −1) is denoted in blue and upregulation (FDR < 0.05, logFC > 1) in orange.
Figure 1. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. (A) The PCA plot shows the variance between the groups (30 min; green, 4 h; blue, 48 h; orange and control; yellow). (B) Heatmap shows the differential expression of the genes (30 min; green, 4 h; blue, 48 h; purple and control; black). Volcano plots show the expression of individual genes in each comparison: 30 min treatment vs. control (C), 4 h treatment vs. control (D), 48 h treatment vs. control (E) 4 h vs. 30 min (F), 48 h vs. 30 min (G), 4 h vs. 48 h (H). Downregulation (FDR < 0.05, logFC < −1) is denoted in blue and upregulation (FDR < 0.05, logFC > 1) in orange.
Ijms 24 06811 g001
Ijms 24 06811 g002 550
Figure 2. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. (A,C) The Venn diagram depicts the relationship between the genes upregulated (FDR < 0.05, logFC > 1). (B,D) The Venn diagram depicts the relationship between the genes downregulated (FDR < 0.05, logFC < −1).
Figure 2. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. (A,C) The Venn diagram depicts the relationship between the genes upregulated (FDR < 0.05, logFC > 1). (B,D) The Venn diagram depicts the relationship between the genes downregulated (FDR < 0.05, logFC < −1).
Ijms 24 06811 g002
Ijms 24 06811 g003 550
Figure 3. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. Resulting transcriptomes were subjected to gene set enrichment analysis to identify enriched Reactome pathways. (A,C) The relative enrichment of pathways belonging to the general theme of extracellular matrix reorganization at each time point. (B,D) The relative enrichment of pathways belonging to the general theme of EGFR-based signalling at each time point.
Figure 3. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. Resulting transcriptomes were subjected to gene set enrichment analysis to identify enriched Reactome pathways. (A,C) The relative enrichment of pathways belonging to the general theme of extracellular matrix reorganization at each time point. (B,D) The relative enrichment of pathways belonging to the general theme of EGFR-based signalling at each time point.
Ijms 24 06811 g003
Ijms 24 06811 g004 550
Figure 4. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. RNA isolated from RL 95-2 cells was subjected to quantitative real-time PCR for key genes that are known to play critical roles in either extracellular matrix reorganization or signalling by EGFR. Expressions of five housekeeping genes were also measured to ascertain the stability of the RL 95-2 transcriptome during the period of incubation. Results are presented as absolute copy number vs. the group (30 min incubation, 0.5 H; 4 h incubation, H_4; 48 h incubation, H_48; Con, untreated control). *, p < 0.05 vs. untreated control. #, p < 0.05 vs. 4 h incubation.
Figure 4. Receptive endometrial analogue RL95-2 cells were treated with trophoblast analogue Jar cell-derived EVs for 30 min, 4 h and 48 h to investigate the effect of trophoblastic EVs on endometrial transcriptome. RNA isolated from RL 95-2 cells was subjected to quantitative real-time PCR for key genes that are known to play critical roles in either extracellular matrix reorganization or signalling by EGFR. Expressions of five housekeeping genes were also measured to ascertain the stability of the RL 95-2 transcriptome during the period of incubation. Results are presented as absolute copy number vs. the group (30 min incubation, 0.5 H; 4 h incubation, H_4; 48 h incubation, H_48; Con, untreated control). *, p < 0.05 vs. untreated control. #, p < 0.05 vs. 4 h incubation.
Ijms 24 06811 g004
Table
Table 1. Alignment details of the libraries used in the study.
Table 1. Alignment details of the libraries used in the study.
Group Total Alignments Successfully Assigned Alignments Percentage
30A737,886197,43526.75
30B724,495248,66134.32
30C832,026221,03526.56
4A803,439214,82526.73
4B635,143202,75431.92
4C790,801242,17530.62
48A628,615187,98329.9
48B747,331242,57832.45
48C758,514195,38225.75
ControlA809,428224,11727.68
ControlB664,114176,64026.59
ControlC641,522184,34328.73
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Godakumara, K.; Heath, P.R.; Fazeli, A. Rhythm of the First Language: Dynamics of Extracellular Vesicle-Based Embryo–Maternal Communication in the Pre-Implantation Microenvironment. Int. J. Mol. Sci. 2023, 24, 6811. https://doi.org/10.3390/ijms24076811

AMA Style

Godakumara K, Heath PR, Fazeli A. Rhythm of the First Language: Dynamics of Extracellular Vesicle-Based Embryo–Maternal Communication in the Pre-Implantation Microenvironment. International Journal of Molecular Sciences. 2023; 24(7):6811. https://doi.org/10.3390/ijms24076811

Chicago/Turabian Style

Godakumara, Kasun, Paul R. Heath, and Alireza Fazeli. 2023. "Rhythm of the First Language: Dynamics of Extracellular Vesicle-Based Embryo–Maternal Communication in the Pre-Implantation Microenvironment" International Journal of Molecular Sciences 24, no. 7: 6811. https://doi.org/10.3390/ijms24076811

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop