Skip to main content

Proteomic profiles and the function of RBP4 in endometrium during embryo implantation phases in pigs

Abstract

Background

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.

Results

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.

Conclusion

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.

Peer Review reports

Introduction

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 [9]. 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 [10]. 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 [16]. 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.

Result

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).

Fig. 1
figure 1

The overview of iTRAQ experiment design. A Total proteins were extracted from porcine endometrium on pregnancy of days (D) 9, 10, 11, 12, 13, 14, 15 and 18 (n = 3). B Three independent 8-plex iTRAQ experiments were performed with the internal reference (the endometrium on D9) on each of the experiments. C Venn Diagram indicated that 2480 commonly identified proteins among the three iTRAQ experiments

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).

Fig. 2
figure 2

The differently expressed proteins in different comparison groups. YK: Days of pregnancy for Yorkshire pigs. Red shows up-regulated and blue shows down-regulated

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.

Fig. 3
figure 3

Hierarchical cluster analysis of different abundance expression. YK: Days of pregnancy for Yorkshire pigs. Red shows higher fold change and blue shows lower fold change

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).

Fig. 4
figure 4

STEM clustering based on different abundance protein. A Each square represents a kind of expression profile and colored squares mean statistical significance. The upper-left number in the box indicates the order of profile. The bottom-left number indicates the P-value. Polyline in the box presents the changing pattern of protein expression levels. B C D Detailed descriptions of protein expression patterns for “profile 1”, “profile 4” and “profile 5”, respectively

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).

Table 1 Validation of iTRAQ results by using MRM method

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).

Fig. 5
figure 5

KEGG analysis of different abundance proteins. Color represents significance, and dot size represents the number of genes clustered in this KEGG term. Green shows significant difference. The number in the parentheses means the number of differentially abundance protein in this term. YK: Days of pregnancy for Yorkshire pigs. “number”: the number of different abundance proteins

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.

Fig. 6
figure 6

The role of RBP4 in EECs and ESCs. A Validation of iTRAQ results of RBP4 by using RT-qPCR. Gray lines and blank column representing the values were determined by iTRAQ and RT-qPCR, respectively. The r-value shows Pearson’s correlation between two sets of values determined by above two methods. B The results of Western blot. C Results of a CCK8 assay after treatment with Si-RBP4 during EECs proliferation. D, E The effects of Si-RBP4 on cell cycle and apoptosis. F, G After treatment with Si-RBP4, cell migration was determined using a transwell assay. In H, relative cell number represents the average cell number found. H, I The effect of Si-RBP4 on the mRNA expression levels of genes related to embryo implantation. J, K Comparison of the number of implantation sites between Si-RBP4 (mmu) injected uterine horn (left, L) and NC-injected uterine horn (right, R) (n = 5). L The expression of RBP4 is presented as a value relative to β-actin (n = 3). Results are mean ± s.d. *P < 0.05, **P < 0.01 (Student’s t-test). Scale bar: 200 µm. NC, negative control

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).

Discussion

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 [6]. 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 [27].

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 [12]. 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 [37]. 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 [37].

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 [42].

RBP4, a key factor related to placentation, was abundant in the luminal epithelium [43]. 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 [12]. 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 [44]. 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.

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).

Bioinformatics analysis

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).

MRM

Protein extraction

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.

Protein digestion

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 [53].

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 [55]. 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.

Transfection

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) [54]. 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].

Western blotting

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) [54]. 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 [54].

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.

Statistical analysis

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.

Abbreviations

ANOVA:

Analysis of variance

2DE:

Two-dimensional gel electrophoresis

CCK-8:

Cell Counting Kit-8

cDNA:

Complementary DNA

c-Myc:

Transcriptional regulator Myc-like

DAPs:

Differentially abundance proteins

EECs:

Endometrial epithelial cells

ESCs:

Endometrial stromal cells

GO:

Gene Ontology

HCA:

Hierarchical cluster analysis

iTRAQ:

Isobaric tags for relative and absolutequantification

KEGG:

Kyoto Encyclopedia of Genes and Genomes

LIF:

LIF interleukin 6 family cytokine

LIFR:

LIF receptor subunit alpha

NC:

Negative control

OD:

Optical density

P-PEPCK:

Phosphoenolpyruvate carboxykinase 1

PCNA:

Proliferating cell nuclear antigen

PVDF:

Polyvinylidene fluoride

RBP4:

Retinol binding protein 4

RPS20:

Ribosomal protein S20

RT-qPCR:

Real time quantitative PCR

SD:

Standard deviation

STEM:

Short time-series expression miner

YK:

Yorkshire

References

  1. Bazer FW, Johnson GA. Pig blastocyst-uterine interactions. Differentiation. 2014;87(1–2):52–65.

    Article  CAS  PubMed  Google Scholar 

  2. Scofield AM, Clegg FG, Lamming GE. Embryonic mortality and uterine infection in the pig. J Reprod Fertil. 1974;36(2):353–61.

    Article  CAS  PubMed  Google Scholar 

  3. Simon C, Martin JC, Pellicer A. Paracrine regulators of implantation. Baillieres Best Pract Res Clin Obstet Gynaecol. 2000;14(5):815–26.

    Article  CAS  PubMed  Google Scholar 

  4. Vigano P, Mangioni S, Pompei F, Chiodo I. Maternal-conceptus cross talk–a review. Placenta. 2003;24 Suppl B:S56-61.

    Article  CAS  PubMed  Google Scholar 

  5. Bazer FW, Spencer TE, Johnson GA, Burghardt RC. Uterine receptivity to implantation of blastocysts in mammals. Front Biosci (Schol Ed). 2011;3:745–67.

    Article  PubMed  Google Scholar 

  6. Wathes C, Whittemore C. Whittemore’s Science and Practice of Pig Production. 3rd ed. 2007.

    Google Scholar 

  7. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  8. 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.

    Article  CAS  PubMed  Google Scholar 

  9. Achache H, Revel A. Endometrial receptivity markers, the journey to successful embryo implantation. Hum Reprod Update. 2006;12(6):731–46.

    Article  PubMed  Google Scholar 

  10. Geisert RD, Johnson GA, Burghardt RC. Implantation and Establishment of Pregnancy in the Pig. Adv Anat Embryol Cell Biol. 2015;216:137–63.

    Article  PubMed  Google Scholar 

  11. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. 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.

    Article  CAS  PubMed  Google Scholar 

  13. 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.

    Article  CAS  PubMed  Google Scholar 

  14. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. 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.

    Article  CAS  PubMed  Google Scholar 

  17. 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.

    Article  CAS  PubMed  Google Scholar 

  18. 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.

    Article  CAS  PubMed  Google Scholar 

  19. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. 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.

    Article  CAS  PubMed  Google Scholar 

  22. 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.

    Article  PubMed  Google Scholar 

  23. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. 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.

  25. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. 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.

    Article  CAS  PubMed  Google Scholar 

  27. 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.

    Article  PubMed  Google Scholar 

  28. 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.

    Article  CAS  PubMed  Google Scholar 

  29. Bazer FW, Spencer TE, Johnson GA, Burghardt RC, Wu G. Comparative aspects of implantation. Reproduction. 2009;138(2):195–209.

    Article  CAS  PubMed  Google Scholar 

  30. 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.

    Article  CAS  PubMed  Google Scholar 

  31. 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.

    Article  PubMed  Google Scholar 

  32. 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.

    Article  CAS  PubMed  Google Scholar 

  33. 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.

    Article  CAS  PubMed  Google Scholar 

  34. 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.

    Article  PubMed  Google Scholar 

  35. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. 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.

    Article  PubMed  Google Scholar 

  37. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. 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.

    Article  CAS  PubMed  Google Scholar 

  39. 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.

    Article  CAS  PubMed  Google Scholar 

  40. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  41. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  42. 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.

    Article  CAS  PubMed  Google Scholar 

  43. 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.

    PubMed  Google Scholar 

  44. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Farrell AS, Sears RC. MYC degradation. Cold Spring Harb Perspect Med. 2014;4(3):a014365.

    Article  PubMed  PubMed Central  Google Scholar 

  46. 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.

    Article  CAS  PubMed  Google Scholar 

  47. Thompson EB. The many roles of c-Myc in apoptosis. Annu Rev Physiol. 1998;60:575–600.

    Article  CAS  PubMed  Google Scholar 

  48. Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein Sci. 2019;28(11):1947–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  49. 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.

  50. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Bronstrup M. Absolute quantification strategies in proteomics based on mass spectrometry. Expert Rev Proteomics. 2004;1(4):503–12.

    Article  PubMed  Google Scholar 

  52. Nikolov M, Schmidt C, Urlaub H. Quantitative mass spectrometry-based proteomics: an overview. Methods Mol Biol. 2012;893:85–100.

    Article  CAS  PubMed  Google Scholar 

  53. James A, Jorgensen C. Basic design of MRM assays for peptide quantification. Methods Mol Biol. 2010;658:167–85.

    Article  CAS  PubMed  Google Scholar 

  54. 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.

    Article  CAS  PubMed  Google Scholar 

  55. 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.

    Article  CAS  PubMed  Google Scholar 

  56. 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.

    Article  CAS  PubMed  Google Scholar 

  57. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  58. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  59. 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.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

No applicable.

Funding

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).

Author information

Authors and Affiliations

Authors

Contributions

YW and RH performed the study, analysis the data and drafted the manuscript. SX, QL and DG collected samples and analyzed the data. ML contributed to the design of the study, acquiring the funding, and modifying the manuscript. YW contributed to providing resource and contributed to acquiring the funding. All authors submitted comments on drafted and read and approved the final manuscript.

Corresponding authors

Correspondence to Renwu Hua or Minggang Lei.

Ethics declarations

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

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-023-09278-5

Keywords