Proteomic profiles and the function of RBP4 in endometrium during embryo implantation phases in pigs
BMC Genomics volume 24, Article number: 200 (2023)
Endometrial receptivity plays a vital role in the success of embryo implantation. However, the temporal proteomic profile of porcine endometrium during embryo implantation is still unclear.
In this study, the abundance of proteins in endometrium on days 9, 10, 11, 12, 13, 14, 15 and 18 of pregnancy (D9, 10, 11, 12, 13, 14, 15 and 18) was profiled via iTRAQ technology. The results showed that 25, 55, 103, 91, 100, 120, 149 proteins were up-regulated, and 24, 70, 169, 159, 164, 161, 198 proteins were down-regulated in porcine endometrium on D10, 11, 12, 13, 14, 15 and 18 compared with that on D9, respectively. Among these differentially abundance proteins (DAPs), Multiple Reaction Monitoring (MRM) results indicated that S100A9, S100A12, HRG and IFI6 were differentially abundance in endometrial during embryo implantation period. Bioinformatics analysis showed that the proteins differentially expressed in the 7 comparisons were involved in important processes and pathways related to immunization, endometrial remodeling, which have a vital effect on embryonic implantation.
Our results reveal that retinol binding protein 4 (RBP4) could regulate the cell proliferation, migration and apoptosis of endometrial epithelial cells and endometrial stromal cells to affect embryo implantation. This research also provides resources for studies of proteins in endometrium during early pregnancy.
The embryo implantation period has been recognized as a vital stage that influences the litter size of sows [1, 2]. Endometrial receptivity, the development of embryos and synchronous communication between the endometrium and embryos have an important impact on the process of embryo implantation [3, 4]. From Day (D) 9 of pregnancy to D12 of pregnancy, the blastocyst changes rapidly from spherical (1–2 mm) to ovoid (4–5 mm) and tubular (> 10 mm) and finally elongates into a filamentous form within 2–3 h (> 100 mm) [5, 6]. Embryos move in the uterine cavity and determine their attachment position in the uterus on D12 and then reach the surface of the endometrium on D13 [7, 8]. The apical surface of the endometrial epithelium forms pinopods that become enclosed by caps of trophectoderm cells through D13 and D14 . The function of the maternal–fetal interface changes from histotrophic to histotrophic and hemotrophic nutrient transport after D15, and implantation is usually completed between D18 and D24 of pregnancy . Despite our increased knowledge on the role of the endometrium in embryo implantation, only limited information is available about the molecular regulation in the endometrium during the embryo implantation period.
Previous studies have focused on the gene expression of the endometrium by transcriptome sequencing during the embryo implantation period [11,12,13,14,15]. Due to posttranscriptional and posttranslational regulation, the transcription results cannot fully represent the changes in the protein during porcine embryo implantation . Proteomic approaches have also been used to study differentially expressed proteins in the endometrium based on the two-dimensional gel electrophoresis (2DE) method during early pregnancy in pigs [17, 18]. However, isobaric tags for relative and absolute quantification (iTRAQ), enabling the identification and quantification of lower abundance proteins, has not yet been employed to analyze the proteomic profile of the endometrium during porcine embryo implantation. 2DE is a top-down approach that quantifies the differentially abundant proteins at the protein level before identifying the protein by LC–MS/MS, while the iTRAQ method is a bottom-up approach in which the whole proteome is first digested with trypsin, and the generated peptides are separated by chromatography and identified and measured by mass spectrometry.
In this study, the abundance of proteins in the endometrium during embryo implantation was comprehensively profiled using iTRAQ technology in Yorkshire (YK) pigs. Moreover, we identified a candidate protein that may play a vital role in embryo implantation. The principal goal of the research was to provide a resource for studies on the temporal abundance profile of proteins during embryo implantation.
iTRAQ analysis reveals changes in protein abundance in porcine endometrium on different gestation days (D9-15 and 18)
An iTRAQ-based quantitative proteomic method was used to obtain the proteomic profile of the porcine endometrium on the following gestation days: D9-15 and 18 (Fig. 1A). The porcine endometrium on D9 was used as the internal reference for three independent iTRAQ experiments (Fig. 1B). A total of 3,774, 3,814 and 3,779 proteins were detected in each of the three iTRAQ experiments. Among these proteins, 2,480 proteins were detected in all three iTRAQ experiments. Then, we used these proteins to characterize the proteomic profiles of the porcine endometrium on different gestation days and carried out different statistical comparisons (Fig. 1C).
The proteins with |Fold change|> 1.2 and Q-value < 0.05 in each comparison group were defined as differentially abundance proteins (DAPs). The results showed that 25, 55, 103, 91, 100, 120, and 149 proteins were upregulated, and 24, 70, 169, 159, 164, 161, and 198 proteins were downregulated in porcine endometrium on D10, 11, 12, 13, 14, 15 and 18 compared with those on D9, respectively (Fig. 2).
To explore the differences in the protein abundance pattern of the porcine endometrium during different phases of embryo implantation, hierarchical cluster analysis (HCA) was performed to cluster all proteins and their abundance profiles in different groups (Fig. 3) (Additional file 1: Table S1). The results indicated that the abundance profiles of proteins in the YK10 vs. YK9 and YK11 vs. YK9 comparison groups were directly clustered together, and the YK12 vs. YK9 and YK13 vs. YK9 comparison groups were also directly clustered together. Additionally, the abundance profiles of the YK14 vs. YK9 and YK15 vs. YK9 comparison groups were similar.
Time-series expression analysis
The short time-series expression miner (STEM) was used to investigate the pattern of differentially abundant protein (DAP) expression during the 8 stages (Additional file 2: Table S2). The typical abundance patterns of proteins (P < 0.05) were profile 1, profile 4 and profile 5 (Fig. 4A). The expression patterns of proteins in profile 1 and profile 4 were characterized by decreased levels from D9 to D18 (Fig. 4B and Fig. 4C). The abundance levels of proteins in profile 5 reached a minimum on D12 and D13 and significantly increased thereafter, but there was no significant change between the expression levels of proteins in profile 5 on D9 and D18 (Fig. 4D).
MRM method validation
The multiple reaction monitoring (MRM) method was carried out on 8 chosen differentially expressed proteins to verify the reliability of the iTRAQ results. All Pearson correlation coefficients between the MRM data and iTRAQ data were > 0.65, indicating that the expression patterns of protein determined by MRM were generally consistent with those obtained by iTRAQ (Table 1).
Functional annotation of DAPs
To analyze the functions of DAPs in porcine endometrium during preimplantation, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out on the DAPs in 7 comparisons (YK10 vs. YK9, YK11 vs. YK9, YK12 vs. YK9, YK13 vs. YK9, YK14 vs. YK9, YK15 vs. YK9 and YK18 vs. YK9). Significantly enriched GO terms mainly included “serine − type endopeptidase inhibitor activity”, “extracellular space”, “extracellular region”, and “negative regulation of endopeptidase activity” in all comparisons (Additional file 3: Fig. S1). The KEGG analysis showed that the top 20 significantly enriched pathways in 7 comparisons were related to metabolism, including “Carbon metabolism”, “Cysteine and methionine metabolism” and “Citrate cycle (TCA cycle)” (Fig. 5 and Additional file 4: Table S3). Notably, “oxidative phosphorylation” was a significantly enriched pathway in the YK11 vs. YK9 YK12 vs. YK9, YK13 vs. YK9, YK14 vs. YK9 and YK15 vs. YK9 comparison groups (Fig. 5).
The role of retinol binding protein 4 (RBP4) in endometrial epithelial cells (EECs) and endometrial stromal cells (ESCs)
In the endometrial tissues of pigs, the mRNA and protein expression levels of RBP4 were significantly upregulated on D12, D13, D14, D15 and D18 compared with D9 (Fig. 6A). The Pearson correlation coefficient between the real-time quantitative PCR (RT‒qPCR) value and iTRAQ value of RBP4 was 0.88 (Fig. 6A). The Western blot of RBP4 is shown in Fig. 6B (Additional file 5: Fig. S2). Previous studies have also indicated that retinol binding protein acts as a chaperone to transport vitamin A to the conceptus and that RBP4 plays a vital role in cell function [6, 19], so we investigated the role of RBP4 in EECs and ESCs. Si-RBP4 (which acts as an inhibitor of RBP4) was transfected into EECs and ESCs. The mRNA expression level of RBP4 was decreased after the transfection of Si-RBP4 in EECs and ESCs (P < 0.05) (Additional file 6: Fig. S3). The Cell Counting Kit-8 (CCK-8) results showed that knockdown of RBP4 significantly inhibited the proliferation of EECs and ESCs (P < 0.05) (Fig. 6C). The cell cycle analysis and apoptosis analysis showed that Si-RBP4 induced G0/G1 phase block and apoptosis of EECs and ESCs (P < 0.01) (Fig. 6D, Fig. 6E and Additional file 7: Fig. S4). Furthermore, Si-RBP4 remarkably hindered the migration of EECs and ESCs (P < 0.01) (Fig. 6F and Fig. 6G). Altogether, knockdown of RBP4 inhibited the proliferation and migration of EECs and induced their G0/G1 phase block and apoptosis.
To analyze the regulatory mechanism of RBP4 in EECs and ESCs, we detected changes in the mRNA expression levels of genes related to embryo implantation after transfection with Si-RBP4. We observed the downregulation of proliferating cell nuclear antigen (PCNA) and transcriptional regulator Myc-like (c-Myc) expression at the mRNA level in EECs and ESCs (P < 0.05) and the downregulation of the mRNA expression levels of LIF receptor subunit alpha (LIFR) (P < 0.01), LIF interleukin 6 family cytokine (LIF) and phosphoenolpyruvate carboxykinase 1 P-PEPCK in ESCs (P < 0.05) (Fig. 6H and Fig. 6I). Moreover, loss of function of RBP4 in vivo resulted in a decreased number of implanted mouse embryos (Fig. 6J and Fig. 6K). The expression of RBP4 in the uterine horn injected with siRBP4 was lower than that in the control uterine horn (P < 0.05) (Fig. 6L).
The present study aimed to explore the temporal protein abundance patterns in porcine endometrium during early pregnancy and identify vital proteins that may regulate embryo implantation. Comprehensive profiling of the expression of protein in endometrial tissues from pregnant gilts (D9, D10, D11, D12, D13, D14, D15 and D18) was performed via iTRAQ technology. 2DE is a top-down approach that quantifies the differentially abundant proteins at the protein level before identifying the protein by LC–MS/MS, while the iTRAQ method is a bottom-up approach in which the whole proteome is first digested with trypsin; then, the generated peptides are separated by chromatography and identified and measured by mass spectrometry [20, 21].
We sequenced endometrial tissue from these eight periods. Endometrial tissue was obtained from the distal, proximal, and middle parts of the uterine horn according to previously reported methods. Some studies also chose the site of hyperemic zones because they selected the 14th day of gestation, which is easy to observe (observation is not easy earlier) [22, 23]. To ensure the consistency of sample sites, we used the same site of sample tissue. However, the endometrial samples were consistent with previous studies, which included the lamina epithelialis, lamina propria, and tela submucosa but not the tunica muscularis. However, recent research has shown that the published laser microdissection technique can separate layers of endometrial cell types [24,25,26]. On the other hand, we used whole endometrial fragments, which may have resulted in proteome profiling.
We also identified candidate DAPs that may have a crucial effect on embryo implantation. Some studies were published that focused on the investigation of the gene expression profiles in the porcine endometrium during the early gestation phases [12,13,14,15]. However, investigation of the protein abundance in the endometrium of pigs during the whole process from the beginning to the completion of embryo implantation is important for understanding the mechanism of uterine function. Among DAPs, we found that the number of downregulated proteins in all comparison groups was greater than that of upregulated proteins except YK10 vs. YK9. To further explore the abundance profile of DAPs in regulating embryo implantation, STEM analysis was performed for all DAPs. The results showed that two of three typical expression patterns of proteins were characterized by decreased levels from D9 to D18. Twelve and 15 days were the two key time points during the preimplantation phases.
Embryos begin to attach to the surface of the endometrium on D13 . In our study, the abundance profiles of protein in the YK10 vs. YK9 and YK11 vs. YK9 comparison groups were not directly clustered with the other groups. We speculated that the protein expression profile of the endometrium changes after the embryo reaches the endometrium. The results also indicated that protein abundance in the endometrium changed significantly from Day 12 of pregnancy, which was consistent with previous studies. Another study focused on uterine luminal fluid proteins on D9 of the oestrus cycle and D9, 12 and 15 of pregnancy and observed significant changes between D9 and 12 of pregnancy .
Porcine endometrium undergoes functional changes for the success of embryo implantation [1, 28]. The GO term “binding” contributed to embryo implantation and was related to upregulated genes in Meishan pigs compared to Yorkshire pigs . These GO terms (e.g., “unfolded protein binding”, “thyroid hormone binding” and “laminin binding”) involved in the remodeling of endometrium were enriched in all comparison groups (YK11 vs. YK9, YK12 vs. YK9, YK13 vs. YK9, YK14 vs. YK9, YK15 vs. YK9 and YK18 vs. YK9) [29,30,31]. Many previous studies have demonstrated that changes in maternal immunity are beneficial for the establishment of pregnancy during the embryo implantation period, which participates in the inflammatory response [32,33,34,35]. We further classified the significantly enriched GO terms according to molecular function, cellular component and biological process. “Serine-type endopeptidase inhibitor activity” was the only molecular function GO term that was enriched in all comparison groups. “Extracellular region” and “extracellular space” were enriched cellular component GO terms in all comparison groups, and other studies have found that this signaling pathway was also enriched at D14 of pregnancy [22, 36]. The genes coding for extracellular proteins have cell adhesion and immune functions. Extracellular matrix remodeling is necessary for ectopic endometrium implantation. Additionally, “negative regulation of endopeptidase activity” was enriched in all comparison groups and was related to biological processes.
Here, many pathways (e.g., “Carbon metabolism” and “Cysteine and methionine metabolism”) related to metabolism were enriched from D11 to 18. Moreover, “oxidative phosphorylation” and “citrate cycle (TCA cycle)” were enriched in five comparison groups (YK11 vs. YK9, YK12 vs. YK9, YK13 vs. YK9, YK14 vs. YK9, YK15 vs. YK9), which focused on the day of implantation. We hypothesized that the energy supply was necessary after the initiation of implantation. In addition to the change in the endometrial signaling pathway, there was also an influence of the number of conceptuses . Variations within the number of conceptuses result in variable levels of their signals in the uterine lumen; therefore, the maternal response could be different, which is also related to the number of corpora lutea. Variation in the number of corpora lutea may result in different levels of P4 affecting the endometrium .
According to previous reports on differentially abundant proteins, we found that many proteins regulating embryo implantation were differentially abundant in our results. β-catenin and DPP4 are associated with embryonic development [38, 39]. Compared with previous transcriptome sequencing studies, we found that the protein levels of some key genes (STAT1, ISG15, AKR1B1 and ANXA8) regulating embryo attachment were also differentially abundant [40, 41]. In addition, some hormone-regulating proteins are differentially abundant, such as ER and GPX3 .
RBP4, a key factor related to placentation, was abundant in the luminal epithelium . A previous study also showed that the mRNA expression level of RBP4 was upregulated on D12 or D15 compared with D9 via RNA-seq and RT‒qPCR . Here, we found that the expression level of RBP4 in pigs from D11 to D18 was significantly higher than that on D9. A previous study indicated that endometrial remodeling, including cell proliferation, apoptosis, cell migration and cytoskeleton organization, imposes strong effects on embryo implantation . We identified that knockdown of RBP4 inhibited proliferation and migration but induced apoptosis of EECs and ESCs. The knockdown of RBP4 suppressed the mRNA expression level of c-Myc involved in cell proliferation, migration and apoptosis [45,46,47]. Furthermore, this study verified that loss in function of RBP4 in vivo could decrease the number of implanted embryos. The above data indicate that RBP4 could regulate the function of EECs and ESCs to influence endometrial receptivity.
In conclusion, the temporal expression of proteins in the endometrium during embryo implantation was profiled. Our research provides a comprehensive overview of the establishment process of endometrial receptivity during early pregnancy in pigs. Furthermore, we found that RBP4 may exert a significant impact on the process of embryo implantation by regulating the proliferation, apoptosis and migration of EECs and ESCs.
Materials and methods
Animals and tissues
In this study, twenty-four YK gilts, which were similar in age and genetic background, were randomly assigned to eight different pregnancy groups. Furthermore, all the gilts mated with the same boar in the third oestrus cycle with the artificial insemination method. Uteri were collected from the gilts that were subsequently slaughtered on D9 (n = 3), 10 (n = 3), 11 (n = 3), 12 (n = 3), 13 (n = 3), 14 (n = 3), 15 (n = 3) and 18 (n = 3) of gestation. Then, the uteri were flushed with PBS quickly and subsequently opened longitudinally along the anti-mesometrial side. The endometrial tissues of eight periods were collected from three locations of each uterine horn: proximal, medial, and distal. The tissue samples were frozen in liquid nitrogen rapidly and stored at -80 °C for protein extraction.
For iTRAQ analysis, the twenty-four samples were homogenized (0.1 g) using a ceramic mortar and pestle and precooled with liquid nitrogen. Then, tissues were disrupted in lysis buffer (with enzyme inhibitors) by a TissueLyser machine. The mixtures were centrifuged at 25,000 × g for fifteen minutes. The supernatant was removed and mixed with 1 volume of cold acetone and stored at -20 °C for two hours. Subsequently, the samples were centrifuged at 25,000 × g for fifteen minutes, and the supernatant was removed again. Lysis buffer was added to dissolve the pellets. To reduce the disulfide bond of peptides, 10 mM DTT was added to the solution and maintained at 56 °C for one hour. A total of 55 mM IAM was added to the solution, which was maintained in a dark room for forty-five minutes. One volume of chilled acetone was added to the solution, stored at -20 °C for two hours and then centrifuged at 25,000 × g for fifteen minutes. Then, the solution was removed. Finally, lysis buffer was added to dissolve the pellet to obtain the protein solution.
Protein digestion and iTRAQ labeling
One hundred micrograms of protein from each sample was transferred to a 1.5 ml new tube for digestion. Then, each sample was mixed with 5 μg of Trypsin Gold for four hours at 37 °C. Subsequently, 5 μg of Trypsin Gold was added for digestion at 37 °C for eight hours. After trypsin digestion, the peptides were dried by vacuum centrifugation and redissolved in 0.5 M TEAB. Finally, the iTRAQ labeling of prepared peptide samples was performed according to the iTRAQ Reagent 8-plex Kit manufacturer’s protocol. The samples were labeled YK9-113, YK10-114, YK11-115, YK12-116, YK13-117, YK14-118, YK15-119, and YK18-121.
Peptide fractionation and LC‒ESI‒MS/MS analysis
The iTRAQ-labeled peptides were fractionated by RP chromatography using the Shimadzu LC-20AB HPLC Pump system. Each fraction was resuspended in buffer A (5% CAN, 0.1% FA). Centrifuging the mixture at 20,000 × g for ten minutes, we obtained the supernatant, of which the peptide final concentration was approximately 0.5 g/l. Finally, according to the experimental protocol, we can acquire the results. The data were acquired with a 2.5 kV ion spray voltage, 30 psi curtain gas, 15 psi nebulizer gas and interface heater temperature of 150 °C. In addition, the data were obtained with a Triple TOF5600 System (AB SCIEX, Concord, ON) fitted with a Nanospray III source (AB SCIEX, Concord, ON) and a pulled quartz tip as the emitter (New Objectives, Woburn, MA) and controlled with Analyst 1.6 software (AB SCIEX, Concord, ON).
The proteins were identified using the UniProt Pig database (http://www.uniprot.org). The raw MS/MS data were converted into MGF format by the ProteoWizard tool msConvert, and the exported MGF files were searched using Mascot version 2.3.02 in this project against the selected database. The proteins were quantified using automated software called IQuant. In this study, functional annotation of differentially abundant proteins was conducted on the basis of DAVID (https://david.ncifcrf.gov/). The GO terms were obtained by ticking the “GOTERM_BP_DIRECT”, “GOTERM_CC_DIRECT” and “GOTERM_MF_DIRECT” of “Gene Ontology”. The KEGG pathways were obtained by ticking the “KEGG_PATHWAY” of “Pathways” (www.kegg.jp/kegg/kegg1.html) [48,49,50]. These pathways and GO terms with P value < 0.05 were considered significantly enriched pathways and GO terms.
The profiles of differentially abundant proteins over time were visualized by Short Time-series Expression Miner v 1.3.8 (STEM, http://www.cs.cmu.edu/~jernst/stem/). The STEM clustering method was selected as the clustering method, and other parameters were set as default. The protein abundance profiles were ordered by significance.
The fold change values of comparisons were normalized by using the log2 (value). HCA was conducted using imageGP (http://www.ehbio.com/ImageGP/index.php/Home/Index/PHeatmap.html).
Samples were ground into powder in liquid nitrogen and extracted with lysis buffer (7 M urea, 2 M thiourea, 4% CHAPS, 40 mM Tris–HCl, pH 8.5) containing 1 mM PMSF and 2 mM EDTA (final concentration). After 5 min, 10 mM DTT (final concentration) was added to the samples. The suspension was sonicated at 200 W for 15 min and then centrifuged at 4 °C and 30 000 × g for 15 min. The supernatant was mixed well with a 5 × volume of chilled acetone containing 10% (v/v) TCA and incubated at -20 °C overnight. After centrifugation at 4 °C and 30 000 g, the supernatant was discarded. The precipitate was washed with chilled acetone three times. The pellet was air-dried and dissolved in lysis buffer (7 M urea, 2 M thiourea, 4% NP40, 20 mM Tris–HCl, pH 8.0–8.5) [51, 52].
The suspension was sonicated at 200 W for 15 min and centrifuged at 4 °C and 300,000 × g for 15 min. The supernatant was transferred to another tube. To reduce disulfide bonds in proteins of the supernatant, 10 mM DTT (final concentration) was added and incubated at 56 °C for 1 h. Subsequently, 55 mM IAM (final concentration) was added to block the cysteines and incubated for 1 h in the darkroom. The supernatant was mixed well with a 5 × volume of chilled acetone for 2 h at -20 °C to precipitate proteins. After centrifugation at 4 °C and 30 000 g, the supernatant was discarded, and the pellet was air-dried for 5 min, dissolved in 500 μL 0.5 M TEAB (Applied Biosystems, Milan, Italy) and sonicated at 200 W for 15 min. Finally, samples were centrifuged at 4 °C and 30 000 × g for 15 min [51, 52]. The supernatant was transferred to a new tube and quantified using a Bradford kit (Bio-Rad). The proteins in the supernatant were stored at -80 °C for further analysis.
Total protein (100 μg) was removed from each sample solution, and then the protein was digested with Trypsin Gold (Promega, Madison, WI, USA) at a protein:trypsin ratio of 30:1 at 37 °C for 16 h. After trypsin digestion, peptides were dried by vacuum centrifugation. Peptides were reconstituted in 0.5 M TEAB.
LC- MRM- MS
Samples were digested as described and spiked with 50 fmol of β-galactosidase for data normalization. MRM analyses were performed on a QTRAP 6500 mass spectrometer (SCIEX, Framingham, MA, USA) equipped with an LC-20AD nanoHPLC system (Shimadzu, Kyoto, Japan). The mobile phase consisted of solvent A, 0.1% aqueous formic acid, and solvent B, 98% acetonitrile with 0.1% formic acid. Peptides were separated on a C18 column (0.075 × 150 mm column, 3.6 μm) at 300 nL/min and eluted with a gradient of 5%-30% solvent B for 38 min, 30%-80% solvent B for 4 min, and maintenance at 80% for 8 min. For the QTRAP 6500 mass spectrometer, a spray voltage of 2400 V, nebulizer gas of 23 p.s.i., and a dwell time of 10 ms were used. Multiple MRM transitions were monitored using unit resolution in both Q1 and Q3 quadrupoles to maximize specificity .
Cell preparation and culture conditions
Isolation and culture of porcine primary endometrial epithelial cells and endometrial stromal cells were performed according to previous research [12, 54]. The endometrium was separated and shredded with sterile scissors. After washing twice with PBS, the tissue pieces were incubated at 37 °C for 2.5 h with collagenase I (Gibco, NY, USA) and shaken vigorously every half hour. Undigested tissue pieces were removed by screen filtration. Then, the filtrate was centrifuged at 500 × g for 10 min to remove the supernatant (fraction rich in endometrial stromal cells). The pellet (epithelial‐rich fraction) was resuspended twice in PBS and recentrifuged (500 g, 10 min) twice. The supernatant (fraction rich in endometrial stromal cells) was recentrifuged (1000 g, 10 min) twice to obtain the pellet (stromal-rich fraction). The pellets of EECs and ESCs were suspended in Dulbecco’s modified Eagle’s medium/F-12 (DMEM/F12; 1:1) medium (Gibco) supplemented with 10% fetal bovine serum (Gibco) and 1% penicillin‒streptomycin (Gibco) and cultured in a 37 °C and 5% CO2 incubator . The epithelial and stromal cells were then trypsinized with 0.25% trypsin–EDTA and placed in a cell culture flask (Corning, NY, USA) for subsequent experiments. Ninety percent purity of populations of endometrial cells was isolated.
All RNA oligonucleotides were designed and synthesized by GenePharm (Shanghai, China). EECs and ESCs were transfected with Si-RBP4 (sense: CCGGCUGAUCACUCACAAUTT, antisense: CCGGCUGAUCACUCACAAUTT) and Si-NC (sense: UUCUUCGAACGUGUCACGUTT, antisense: ACGUGACACGUUCGGAGAATT) using Lipofectamine® 2000 (Invitrogen Life Technology, Shanghai, China) according to the manufacturer’s instructions. The transfection efficiency was more than 90%.
RNA extraction and real-time quantitative PCR (RT‒qPCR)
RNA was extracted from cells using TRIzol (Invitrogen) as recommended by the manufacturer, and the concentration and quality were measured by a NanoDrop 2000 (ThermoFisher, Waltham, USA) . Complementary DNA (cDNA) was synthesized with a reverse transcription kit (Takara, Tokyo, Japan). Then, Mix (Toyobo, Japan) and specific primers for every gene were used to perform RT‒qPCR on a Real-time System (Roche, Basel, Switzerland) (Additional file 8: Table S4). The expression levels of genes were normalized to ribosomal protein S20 (RPS20) to obtain the relative expression using the 2–ΔΔCt method according to a previous study [12, 56].
Total protein was isolated from the endometrium using SDS Lysis Buffer (P0013B, Beyotime Ltd. China). Protein content was measured using an enhanced BCA protein assay kit (P0010, Beyotime, Ltd. China). The sample was denatured by heating, separated by SDS‒PAGE and transferred to a PVDF (polyvinylidene fluoride) membrane. Next, the membranes were blocked with 5% skimmed milk powder and separately probed with rabbit anti-RBP4 (catalog number: ab233138, Abcam) overnight at 4 °C with a final dilution of 1:1000 (v/v) in 5% milk. After three washes, the membranes were incubated with goat anti-rabbit secondary antibodies (GB23303, Sevicebio, Wuhan, China) at a 1:2000 dilution (v/v) at 37 °C for 1.5 h. Images of membranes treated with ECL (enhanced chemiluminescence) were captured by a Western blotting Detection System (Tiangen, Beijing, China).
Cell proliferation assay
CCK-8 (Dojindo, Shanghai, China) was used to measure cell proliferation 48 h after transfection with the interference fragment of RBP4 following the manufacturer’s instructions [57, 58]. The optical density (OD) at 450 nm of each well plate was determined using a microplate reader (Bio-Rad, CA, USA).
Cell apoptosis and cell cycle
Cell apoptosis analysis was carried out using the AnnexinV-FITC/PI apoptosis kit (Dojindo, Shanghai, China) according to the manufacturer’s protocol. Analyses were performed using a flow cytometer (Beckman, CA, USA).
For cell cycle analysis, the cells were separated by trypsin digestion. The precipitate was collected by centrifugation, washed twice, and then fixed. The supernatant was removed by centrifugation. Each sample was stained with propidium iodide (PI; KeyGEN BioTECH, Jiangsu) for 30 min in the dark at room temperature. Analyses were performed using a flow cytometer (Beckman, CA, USA).
Cell migration assay
Cell migration was performed based on our previous study and assessed by a transwell assay that had 12 mm polycarbonate membranes with an 8.0 μm pore size (Corning, NY, USA) . The EECs were resuspended in serum-free medium as a single-cell solution at 4 h posttransfection. Approximately 2 × 105 EECs or ESCs were seeded on the upper chambers, and complete medium with 10% FBS was added to the lower chamber as a chemoattractant. After 24 h of incubation at 37 °C with 5% CO2, cells that migrated to the lower chamber were fixed with 4% paraformaldehyde for 5 min, stained with 0.1% crystal violet for 5 min, rinsed three times in PBS and subjected to an OLYMPUS DP80 microscope (Tokyo, Japan). The cell numbers of migration were obtained by counting five fields per membrane and represent the average of three independent experiments .
Intrauterine injection of mice
It is difficult to perform in vivo experiments in pigs, so we used a mouse implantation model based on other studies and focused on pig or human implantation [12, 54, 59]. The date of finding the vaginal suppository after mating was designated as the first day. Intrauterine injection surgery under general anesthesia was performed on eight mice on Day 3 of pregnancy in the evening according to previous studies (Sun et al., 2014; Wang et al., 2019). Ten microliters of 10 μM mol/L Si-RBP4 (mmu) (sense: GGCCCUCCAGAAACAGUUUTT, antisense: AAACUGUUUCUGGAGGGCCTT) and inhibitor negative control (NC) were injected into the left and right uterine horns, respectively. Then, the wounds were sutured, and the mice were placed under a 37 °C warmer until the anesthesia diminished. To reduce the pain that the mice endured, temgesic was injected into the mice at 12, 24, and 48 h after surgery. On D5, the endometrium of three mice was isolated to examine the protein expression level of RBP4. On Day 7, five mice were killed, and their uteri were isolated to record the number of implanted embryos.
Data were expressed as means ± standard deviation (SD) derived from at least three independent experiments. Student’s t-test was used for comparisons between two groups. One-way analysis of variance (ANOVA) test for multiple comparisons. P value < 0.05 was considered as statistically significant. * P < 0.05; ** P < 0.01.
Availability of data and materials
All new sequencing data generated in this study have been deposited to ProteomeXchange Consortium via the PRIDE under accession number PXD030681.
Analysis of variance
Two-dimensional gel electrophoresis
Cell Counting Kit-8
Transcriptional regulator Myc-like
Differentially abundance proteins
Endometrial epithelial cells
Endometrial stromal cells
Hierarchical cluster analysis
Isobaric tags for relative and absolutequantification
Kyoto Encyclopedia of Genes and Genomes
LIF interleukin 6 family cytokine
LIF receptor subunit alpha
Phosphoenolpyruvate carboxykinase 1
Proliferating cell nuclear antigen
Retinol binding protein 4
Ribosomal protein S20
Real time quantitative PCR
Short time-series expression miner
Bazer FW, Johnson GA. Pig blastocyst-uterine interactions. Differentiation. 2014;87(1–2):52–65.
Scofield AM, Clegg FG, Lamming GE. Embryonic mortality and uterine infection in the pig. J Reprod Fertil. 1974;36(2):353–61.
Simon C, Martin JC, Pellicer A. Paracrine regulators of implantation. Baillieres Best Pract Res Clin Obstet Gynaecol. 2000;14(5):815–26.
Vigano P, Mangioni S, Pompei F, Chiodo I. Maternal-conceptus cross talk–a review. Placenta. 2003;24 Suppl B:S56-61.
Bazer FW, Spencer TE, Johnson GA, Burghardt RC. Uterine receptivity to implantation of blastocysts in mammals. Front Biosci (Schol Ed). 2011;3:745–67.
Wathes C, Whittemore C. Whittemore’s Science and Practice of Pig Production. 3rd ed. 2007.
Geisert RD, Lucy MC, Whyte JJ, Ross JW, Mathew DJ. Cytokines from the pig conceptus: roles in conceptus development in pigs. J Anim Sci Biotechnol. 2014;5(1):51.
Keys JL, King GJ. Microscopic examination of porcine conceptus-maternal interface between days 10 and 19 of pregnancy. Am J Anat. 1990;188(3):221–38.
Achache H, Revel A. Endometrial receptivity markers, the journey to successful embryo implantation. Hum Reprod Update. 2006;12(6):731–46.
Geisert RD, Johnson GA, Burghardt RC. Implantation and Establishment of Pregnancy in the Pig. Adv Anat Embryol Cell Biol. 2015;216:137–63.
Wang Y, Xue S, Liu X, Liu H, Hu T, Qiu X, Zhang J, Lei M. Analyses of Long Non-Coding RNA and mRNA profiling using RNA sequencing during the pre-implantation phases in pig endometrium. Sci Rep. 2016;6:20238.
Wang Y, Hua R, Xue S, Li W, Wu L, Kang T, Lei M. mRNA/lncRNA expression patterns and the function of fibrinogen-like protein 2 in Meishan pig endometrium during the preimplantation phases. Mol Reprod Dev. 2019;86(4):354–69.
Li W, Xi Y, Xue S, Wang Y, Wu L, Liu H, Lei M. Sequence analysis of microRNAs during pre-implantation between Meishan and Yorkshire pigs. Gene. 2018;646:20–7.
Krawczynski K, Bauersachs S, Reliszko ZP, Graf A, Kaczmarek MM. Expression of microRNAs and isomiRs in the porcine endometrium: implications for gene regulation at the maternal-conceptus interface. BMC Genomics. 2015;16:906.
Hong L, Liu R, Qiao X, Wang X, Wang S, Li J, Wu Z, Zhang H. Differential microRNA Expression in Porcine Endometrium Involved in Remodeling and Angiogenesis That Contributes to Embryonic Implantation. Front Genet. 2019;10:661.
Perez-Deben S, Bellver J, Alama P, Salsano S, Quinonero A, Sebastian-Leon P, Diaz-Gimeno P, Dominguez F. iTRAQ comparison of proteomic profiles of endometrial receptivity. J Proteomics. 2019;203:103381.
Jalali BM, Bogacki M, Dietrich M, Likszo P, Wasielak M. Proteomic analysis of porcine endometrial tissue during peri-implantation period reveals altered protein abundance. J Proteomics. 2015;125:76–88.
Kolakowska J, Souchelnytskyi S, Saini RKR, Franczak A. Proteomic analysis of the endometrium during early pregnancy in the domestic pig. Reprod Fertil Dev. 2017;29(11):2255–68.
Li H, Cao G, Zhang N, Lou T, Wang Q, Zhang Z, Liu C. RBP4 regulates trophoblastic cell proliferation and invasion via the PI3K/AKT signaling pathway. Mol Med Rep. 2018;18(3):2873–9.
Zhan X, Li B, Zhan X, Schluter H, Jungblut PR, Coorssen JR. Innovating the Concept and Practice of Two-Dimensional Gel Electrophoresis in the Analysis of Proteomes at the Proteoform Level. Proteomes. 2019;7(4):36.
Guo J, Ren Y, Hou G, Wen B, Xian F, Chen Z, Cui P, Xie Y, Zi J, Lin L, et al. A Comprehensive Investigation toward the Indicative Proteins of Bladder Cancer in Urine: From Surveying Cell Secretomes to Verifying Urine Proteins. J Proteome Res. 2016;15(7):2164–77.
Samborski A, Graf A, Krebs S, Kessler B, Bauersachs S. Deep sequencing of the porcine endometrial transcriptome on day 14 of pregnancy. Biol Reprod. 2013;88(4):84.
Zeng S, Ulbrich SE, Bauersachs S. Spatial organization of endometrial gene expression at the onset of embryo attachment in pigs. BMC Genomics. 2019;20(1):895.
Xu QX, Zhang WQ, Lu L, Wang KZ, Su RW. Distinguish Characters of Luminal and Glandular Epithelium from Mouse Uterus Using a Novel Enzyme-Based Separation Method. Reprod Sci. 2022.
Aikawa S, Hirota Y, Fukui Y, Ishizawa C, IIda R, Kaku T, Hirata T, Akaeda S, Hiraoka T, Matsuo M, et al. A gene network of uterine luminal epithelium organizes mouse blastocyst implantation. Reprod Med Biol. 2022;21(1):e12435.
Cummings M, Mappa G, Orsi NM. Laser Capture Microdissection and Isolation of High-Quality RNA from Frozen Endometrial Tissue. Methods Mol Biol. 2018;1723:155–66.
Orzechowska K, Kopij G, Paukszto L, Dobrzyn K, Kiezun M, Jastrzebski J, Kaminski T, Smolinska N. Chemerin effect on transcriptome of the porcine endometrium during implantation determined by RNA-sequencingdagger. Biol Reprod. 2022;107(2):557–73.
Burghardt RC, Burghardt JR, Taylor JD 2nd, Reeder AT, Nguen BT, Spencer TE, Bayless KJ, Johnson GA. Enhanced focal adhesion assembly reflects increased mechanosensation and mechanotransduction at maternal-conceptus interface and uterine wall during ovine pregnancy. Reproduction. 2009;137(3):567–82.
Bazer FW, Spencer TE, Johnson GA, Burghardt RC, Wu G. Comparative aspects of implantation. Reproduction. 2009;138(2):195–209.
Waclawik A, Kaczmarek MM, Blitek A, Kaczynski P, Ziecik AJ. Embryo-maternal dialogue during pregnancy establishment and implantation in the pig. Mol Reprod Dev. 2017;84(9):842–55.
Bauersachs S, Wolf E. Uterine responses to the preattachment embryo in domestic ungulates: recognition of pregnancy and preparation for implantation. Annu Rev Anim Biosci. 2015;3:489–511.
Amjadi F, Zandieh Z, Mehdizadeh M, Aghajanpour S, Raoufi E, Aghamajidi A, Aflatoonian R. The uterine immunological changes may be responsible for repeated implantation failure. J Reprod Immunol. 2020;138:103080.
Bischof RJ, Lee R, Lee CS, Meeusen E. Dynamic changes in the lymphocyte subpopulations of pig uterine lymph nodes. Vet Immunol Immunopathol. 1996;51(3–4):315–24.
Chen X, Li A, Chen W, Wei J, Fu J, Wang A. Differential gene expression in uterine endometrium during implantation in pigs. Biol Reprod. 2015;92(2):52.
Kim M, Seo H, Choi Y, Shim J, Kim H, Lee CK, Ka H. Microarray Analysis of Gene Expression in the Uterine Endometrium during the Implantation Period in Pigs. Asian-Australas J Anim Sci. 2012;25(8):1102–16.
Zhang L, Xiong W, Xiong Y, Liu H, Li N, Du Y, Liu Y. Intracellular Wnt/Beta-Catenin Signaling Underlying 17beta-Estradiol-Induced Matrix Metalloproteinase 9 Expression in Human Endometriosis. Biol Reprod. 2016;94(3):70.
Orzechowska K, Dobrzyn K, Kiezun M, Malinowska A, Swiderska B, Kaminski T, Smolinska N. Chemerin Effect on the Endometrial Proteome of the Domestic Pig during Implantation Obtained by LC-MS/MS Analysis. Cells. 2022;11(7):1161.
Yakovleva AA, Zolotov NN, Sokolov OY, Kost NV, Kolyasnikova KN, Micheeva IG. Dipeptidylpeptidase 4 (DPP4, CD26) activity in the blood serum of term and preterm neonates with cerebral ischemia. Neuropeptides. 2015;52:113–7.
Kiewisz J, Krawczynski K, Lisowski P, Blitek A, Zwierzchowski L, Ziecik AJ, Kaczmarek MM. Global gene expression profiling of porcine endometria on Days 12 and 16 of the estrous cycle and pregnancy. Theriogenology. 2014;82(6):897–909.
Gou R, Zhu L, Zheng M, Guo Q, Hu Y, Li X, Liu J, Lin B. Annexin A8 can serve as potential prognostic biomarker and therapeutic target for ovarian cancer: based on the comprehensive analysis of Annexins. J Transl Med. 2019;17(1):275.
Zeng S, Bick J, Ulbrich SE, Bauersachs S. Cell type-specific analysis of transcriptome changes in the porcine endometrium on Day 12 of pregnancy. BMC Genomics. 2018;19(1):459.
Riesewijk A. Gene expression profiling of human endometrial receptivity on days LH+2 versus LH+7 by microarray technology. Mol Hum Reprod. 2003;9(5):253–64.
van der Weijden VA, Puntar B, Rudolf Vegas A, Milojevic V, Schanzenbach CI, Kowalewski MP, Drews B, Ulbrich SE. Endometrial luminal epithelial cells sense embryo elongation in the roe deer independent of interferon-tau†. Biol Reprod. 2019;101(5):882–92.
Kaczmarek MM, Najmula J, Guzewska MM, Przygrodzka E. MiRNAs in the Peri-Implantation Period: Contribution to Embryo-Maternal Communication in Pigs. Int J Mol Sci. 2020;21(6):2229.
Farrell AS, Sears RC. MYC degradation. Cold Spring Harb Perspect Med. 2014;4(3):a014365.
Sun J, Zhang H, Tao D, Xie F, Liu F, Gu C, Wang M, Wang L, Jiang G, Wang Z, et al. CircCDYL inhibits the expression of C-MYC to suppress cell growth and migration in bladder cancer. Artif Cells Nanomed Biotechnol. 2019;47(1):1349–56.
Thompson EB. The many roles of c-Myc in apoptosis. Annu Rev Physiol. 1998;60:575–600.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–92.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Bronstrup M. Absolute quantification strategies in proteomics based on mass spectrometry. Expert Rev Proteomics. 2004;1(4):503–12.
Nikolov M, Schmidt C, Urlaub H. Quantitative mass spectrometry-based proteomics: an overview. Methods Mol Biol. 2012;893:85–100.
James A, Jorgensen C. Basic design of MRM assays for peptide quantification. Methods Mol Biol. 2010;658:167–85.
Hua R, Zhang X, Li W, Lian W, Liu Q, Gao D, Wang Y, Lei M. Ssc-miR-21–5p regulates endometrial epithelial cell proliferation, apoptosis and migration via the PDCD4/AKT pathway. Journal of cell science. 2020;133(23):jcs248898.
Ka H, Jaeger LA, Johnson GA, Spencer TE, Bazer FW. Keratinocyte growth factor is up-regulated by estrogen in the porcine uterine endometrium and functions in trophectoderm cell proliferation and differentiation. Endocrinology. 2001;142(6):2303–10.
Wang Y, Hu T, Wu L, Liu X, Xue S, Lei M. Identification of non-coding and coding RNAs in porcine endometrium. Genomics. 2017;109(1):43–50.
Fu Y, Zhang Y, Lei Z, Liu T, Cai T, Wang A, Du W, Zeng Y, Zhu J, Liu Z, et al. Abnormally activated OPN/integrin alphaVbeta3/FAK signalling is responsible for EGFR-TKI resistance in EGFR mutant non-small-cell lung cancer. J Hematol Oncol. 2020;13(1):169.
Lu Q, Liu T, Feng H, Yang R, Zhao X, Chen W, Jiang B, Qin H, Guo X, Liu M, et al. Circular RNA circSLC8A1 acts as a sponge of miR-130b/miR-494 in suppressing bladder cancer progression via regulating PTEN. Mol Cancer. 2019;18(1):111.
Sun X, Ruan YC, Guo J, Chen H, Tsang LL, Zhang X, Jiang X, Chan HC. Regulation of miR-101/miR-199a-3p by the epithelial sodium channel during embryo implantation: involvement of CREB phosphorylation. Reproduction. 2014;148(6):559–68.
This work was supported by grants: National Natural Foundation of China (82101756), China Agriculture Research System (CARS-36) and Jining key research and development plan (2021YXNS046).
Ethics approval and consent to participate
All animal experiments were conducted with the regulation (No.5 proclaim of the Standing Committee of Hubei People’s Congress). This study was approved by the institutional review board of Standing Committee of Hubei People’s Congress, P.R. China. All methods were carried out in accordance with relevant guidelines and regulations. This study was carried out in compliance with the ARRIVE guidelines (Animal Research: Reporting of In Vivo Experiments).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Wang, Y., Xue, S., Liu, Q. et al. Proteomic profiles and the function of RBP4 in endometrium during embryo implantation phases in pigs. BMC Genomics 24, 200 (2023). https://doi.org/10.1186/s12864-023-09278-5
- Proteomic profiles