- Open Access
The proteome of IVF-induced aberrant embryo-maternal crosstalk by implantation stage in ewes
Journal of Animal Science and Biotechnology volume 11, Article number: 7 (2020)
Implantation failure limits the success of in vitro fertilization and embryo transfer (IVF-ET). Well-organized embryo-maternal crosstalk is essential for successful implantation. Previous studies mainly focused on the aberrant development of in vitro fertilized (IVF) embryos. In contrast, the mechanism of IVF-induced aberrant embryo-maternal crosstalk is not well defined.
In the present study, using ewes as the model, we profiled the proteome that features aberrant IVF embryo-maternal crosstalk following IVF-ET. By comparing in vivo (IVO) and IVF conceptuses, as well as matched endometrial caruncular (C) and intercaruncular (IC) areas, we filtered out 207, 295, and 403 differentially expressed proteins (DEPs) in each comparison. Proteome functional analysis showed that the IVF conceptuses were characterized by the increased abundance of energy metabolism and proliferation-related proteins, and the decreased abundance of methyl metabolism-related proteins. In addition, IVF endometrial C areas showed the decreased abundance of endometrial remodeling and redox homeostasis-related proteins; while IC areas displayed the aberrant abundance of protein homeostasis and extracellular matrix (ECM) interaction-related proteins. Based on these observations, we propose a model depicting the disrupted embryo-maternal crosstalk following IVF-ET: Aberrant energy metabolism and redox homeostasis of IVF embryos, might lead to an aberrant endometrial response to conceptus-derived pregnancy signals, thus impairing maternal receptivity. In turn, the suboptimal uterine environment might stimulate a compensation effect of the IVF conceptuses, which was revealed as enhanced energy metabolism and over-proliferation.
Systematic proteomic profiling provides insights to understand the mechanisms that underlie the aberrant IVF embryo-maternal crosstalk. This might be helpful to develop practical strategies to prevent implantation failure following IVF-ET.
In mammals, well-orchestrated embryo-maternal crosstalk during the implantation stage is of prime importance to establish and maintain pregnancy. Despite the diversity of implantation and placentation strategies, the reciprocal interaction that occurs between the embryos and the maternal uterine endometria is shared among species. Briefly, pregnancy recognition signals from peri-implantation embryos act on the endometrium in a paracrine manner to stimulate uterine receptivity, which supports conceptuses development [1,2,3]. Aberrant crosstalk impairs embryo development and endometrium receptivity, thus leading to implantation failure, which is the most prominent factor for pregnancy loss following both natural conception and when using assistant reproductive technology (ART) [4, 5].
Until now, the success rate following in vitro fertilization and embryo transfer (IVF-ET) has remained disappointingly low among species [6, 7]. Implantation failure remains the biggest barrier that limits the success rate. Aiming to improve the IVF success rate, numerous studies have been performed that focused on the mechanisms responsible for impaired development potential of IVF embryos, including many genetic and cellular changes, such as epigenetic modifications [8, 9], genetic information processing , energy metabolism , and cytoskeleton organization . Based on these observations, strategies have been applied to correct these aberrations, thus enhancing IVF embryo developmental potential [9, 12,13,14,15].
However, endometrial receptivity, an early sensor of embryo implantation signals and a prerequisite of successful pregnancy , has not been fully considered as a target for improving IVF outcomes. In fact, limited prior work demonstrated that the endometrium responds aberrantly to in vitro fertilized or cloned embryos, compared with its response to in vivo fertilized embryos [1, 4, 17]. This suggests that not only impaired embryo quality, but also the aberrant endometrial receptivity induced by IVF embryos, might contribute to implantation failure following IVF-ET. In clinical practice, treatment of IVF/ET patients with N-acetyl-L-cysteine (NAC), melatonin and selenomethionine during the peri-implantation stage has been used to improve IVF outcomes, suggesting that the maternal endometrium can be used as target, and provides alternative strategies for improving the IVF success rate . However, the mechanism, underlying aberrant embryo-maternal crosstalk following IVF-ET is not well defined.
High-throughput methodologies, such as transcriptomic and proteomic analyses, have been frequently applied to profile the cellular responses of embryos and the endometria by implantation stage in various animal models, including sheep [19,20,21,22], cattle [23,24,25], mice [26, 27], pigs , and humans [29, 30]. Among these, ruminants have been used extensively as models to explore embryo-maternal interactions by implantation stage [31, 32]. Distinct from rodents or humans, the embryo-maternal crosstalk in ruminants is characterized by interaction occurring in both the caruncular (C) and intercaruncular (IC) areas. Aglandular C areas serve as the sites of superficial attachment and placentation. Glandular IC areas, which contain large numbers of branched and coiled uterine glands, are mainly responsible for the synthesis and secretion of histotroph [33, 34]. In the present study, using ewes as the model, we profiled the proteome of aberrant embryo-maternal crosstalk following IVF-ET. Compared with the high-throughput analysis focusing on mRNA expression, proteomic analysis provides a more direct and accurate understanding, because proteins are the executors of most biological programs . Therefore, the proteome of the conceptuses produced by IVF and their matched endometria by implantation stage will provide a novel and detailed reference to understand the mechanisms that underlie aberrant IVF embryo-maternal crosstalk, and will provide important clues to improve IVF outcomes from both the embryonic and maternal sides.
Animals and treatment
The experiments were performed in accordance with the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching, and all procedures were approved by the Institutional Animal Care and Use Committee at China Agricultural University (Beijing, China). Chinese Small Tail Han ewes with normal estrous cycles were selected for the present study. The procedures of estrous synchronization, superovulation, artificial insemination (AI), collection, and transfer of IVO blastocysts were performed as described in our previous study .
The methods for IVF were conducted as described by Ptak et al. . Following collection of sheep ovaries at slaughter, oocytes were aspirated using 12 G needles and placed into oocyte wash buffer TCM199-hepes (Sigma, St. Louis, MO, USA), 1 mg/mL polyvinyl alcohol (Sigma), 10–20 μg/mL heparin sodium (Sigma), P/S (100 IU/mL penicillin (Sigma) and 100 IU/mL streptomycin (Sigma). Oocytes, surrounded by integrated granulosa cells and with evenly granulated cytoplasm were selected for in vitro maturation (IVM). Oocytes from donors were incubated in maturation medium TCM199–HCO3 (Sigma) containing 10% FBS (fetal bovine serum; GIBCO, Grand Island, NY, USA), 10 μg/mL FSH (follicle stimulating hormone, Vetrepharm, Concord, Canada), 10 μg/mL LH (luteinizing hormone, Sigma), 1 μg/mL estradiol (Sigma), 10 ng/mL epidermal growth factor (EGF), 0.1 mmol/L cysteamine [36, 37] (Sigma), and P/S covered with mineral oil and incubated in a humidified atmosphere of 5% CO2 at 38.6 °C for 24–26 h. Matured oocytes were gently denuded of granulosa cells by 0.05% hyaluronidase (Sigma) and transferred into 50 μL drops of synthetic oviductal fluid (SOF, Sigma) enriched with 20% (v/v) serum, 2.9 mmol/L of Ca lactate, and 16 mmol/L of isoproterenol. Ram semen was thawed and capacitated, then in vitro fertilization (IVF) was performed in drops using sperm at a final concentration about 2 × 106 spermatozoa/mL and 15–20 oocytes per drop. After 20 h, zygotes were transferred into a four-well plate with 500 μL of in vitro culture medium [SOF enriched with bovine serum albumin (BSA)], covered with mineral oil, and incubated a humidified atmosphere of 5% CO2, 5% O2, and 90% N2 at 38.6 °C. Two well-developed day 6.5 IVF blastocysts were transferred to each synchronized recipient ewe, thus, the day of fertilization was defined as day 0. To minimize the differences caused by the surgical procedure and embryo quality, we selected a skillful technician to perform the surgical procedure, and only good-quality blastocysts (Grade 1) were transferred to synchronized recipient ewes.
We collected good-quality IVO embryos from thirty donors at day 6.5 of pregnancy, and day 6.5 IVF embryos from the IVF process. Then, two well-developed blastocysts were transferred per synchronized recipient ewe (forty-eight synchronized ewes for the IVO group, thirty-eight for the IVF group). Sampling procedures were similar to the methods detailed in our previous study . Briefly, all recipients were slaughtered at day17 of pregnancy, then their uteri were collected and the conceptuses were flushed out with phosphate-buffered saline (PBS). Thirty-seven recipients in the IVO group and twenty recipients in the IVF group had filamentous conceptuses. The endometrial caruncular (C) and intercaruncular (IC) areas were collected and processed as described by Attia et al. . Opening the ipsilateral uterine horn longitudinally by scissors, the C areas were first carefully cut out and collected, and then the IC areas were sampled. The same technician took the samples from the IVO and IVF group, which were stored at liquid nitrogen until further analysis (Additional file 1: Figure S1A).
We equally divided thirty-six IVO samples (or eighteen IVF samples) into three pools, with twelve IVO samples in each pool (six IVF samples per pool). Each pool was ground to powder in liquid nitrogen and stored overnight at − 20 °C after adding a five-fold volume of chilled acetone containing 10% trichloroacetic acid (TCA) and 10 mmol/L dithiothreitol (DTT). The samples were then centrifuged at 4 °C, 16,000×g for 20 min and the supernatant was discarded. The precipitates were mixed with 1 mL chilled acetone containing 10 mmol/L DTT and centrifuged at 4 °C, 20,000×g for 30 min after storing for 30 min at − 20 °C. Centrifugation was repeated several times until the supernatant was colorless. The pellets were air-dried, dissolved in lysis buffer containing 1 mmol/L phenylmethanesulfonyl fluoride (PMSF), 2 mmol/L ethylenediaminetetraacetic acid (EDTA), and 10 mmol/L DTT and sonicated at 200 W for 15 min before being centrifuged at 30,000×g at room temperature for 30 min. The protein concentration in the supernatant was then detected by using the Bradford method.
Proteins (50 μg) were taken from each sample, and isopycnic samples were prepared by adding 8 mol/L urea solution. To reduce disulfide bonds, the samples were incubated with 10 mmol/L DTT at 56 °C for 1 h, and then cysteine bonding was blocked using 55 mmol/L iodoacetamide (IAM) in a dark room for 45 min. Thereafter, each sample was diluted 8-fold with 50 mmol/L ammonium bicarbonate and digested with Trypsin Gold at a protein: trypsin ratio of 20:1 at 37 °C for 16 h. Following desalting using a Strata X C18 column (Phenomenex, Torrance, CA, USA), the samples were vacuum dried. Peptides generated from digestion were directly loaded for liquid chromatography/electrospray ionization tandem mass spectroscopy (LC-ESI-MS/MS) analysis.
LC-ESI-MS/MS analysis with LTQ-orbitrap collision induced dissociation (CID)
Each sample was resuspended in buffer A [2% acetonitrile (ACN), 0.1% formic acid (FA)] and centrifuged at 20,000×g for 10 min. The final peptide concentration for each sample was approximately 0.5 μg/mL. The digested samples were fractionated using a Shimadzu LC-20 AD nano-high performance liquid chromatography (HPLC) system (Shimadzu, Kyoto, Japan). Each sample (10 μL) was loaded by the autosampler onto a 2-cm C18 trap column (200 μm inner diameter), and the peptides were eluted onto a resolving 10 cm analytical C18 column (75 μm inner diameter) prepared in-house. The samples were loaded at a flow rate of 15 μL/min for 4 min, and then a 91-min gradient from 2% to 35% buffer B (98% ACN, 0.1%FA) was run at a flow rate of 400 nL/min, followed by a 5-min linear gradient to 80% buffer B that was maintained for 8 min before finally returning to 2% buffer B within 2 min. The peptides were subjected to nanoelectrospray ionization and then detected by MS/MS in an LTQ Orbitrap Velos (Thermo Fisher Scientific, Bremen, Germany) coupled online to a HPLC system. Intact peptides were detected in the Orbitrap analyzer at a resolution of 60,000 m/z. Peptides were selected for MS/MS using the CID operating mode with a normalized collision energy setting of 35%, and ion fragments were detected in the LTQ. One MS scan followed by ten MS/MS scans was applied for the ten most abundant precursor ions above a threshold ion count of 5000 in the MS survey scan. Dynamic exclusion was used, with the following parameters: Repeat counts = 2; repeat duration = 30 s; and exclusion duration = 120 s. The applied electrospray voltage was 1.5 kV. Automatic gain control (AGC) was used to prevent overfilling of the ion trap; 1 × 104 ions were accumulated in the ion trap to generate CID spectra. For MS scans, the m/z scan range was 350 to 2000 Da.
MaxQuant software (version 18.104.22.168) was used to analyze the mass spectra. Bos taurus is the only well-annotated species with a genomic database with a high degree of homology to sheep. Therefore, we generated one reference protein database by integrating the following databases and sequences of cattle proteins and limited publically available sheep proteins, and removed duplicate proteins: GenBank nr (20110403), Uniprot cow proteins (20110503), sheep proteins (http://www.livestockgenomics.csiro.au/sheep/), and cow proteins (http://genomes.arc.georgetown.edu/drupal/bovine/). The MS/MS data were searched against the reference protein database using the search engine embedded in MaxQuant. Up to two missed cleavages were allowed. The first search was set to 20 ppm, and the MS/MS tolerance for CID was set to 0.5 Da. The false discovery rate (FDR) was set to 0.01 for peptide and protein identifications, which was estimated based on the fraction of reverse protein hits [38, 39]. Proteins were considered identified when at least two peptides were identified, at least one of which was uniquely assignable to the corresponding sequence. In the case of identified peptides that were all shared between two proteins, these were combined and reported as one protein group. To control the false match frequency, the contents of the protein table were filtered to eliminate identifications from the reverse database and common contaminants [40, 41]. The minimum peptide length was set to six amino acids. To perform label-free quantification analysis, the MaxQuant software suite containing an algorithm based on the extracted ion currents (XICs) of the peptides was used. Xcalibur 2.1 (Thermo Scientific) was used as quality control program to check the quality of chromatographs. This specific label-free processing method was performed as described by Waanders et al. .
A human endometrial cancer cell line (Ishikawa, ATCC, USA). Ishikawa cells were grown at 37 °C in DMEM (Hyclone, Logan, UT) supplemented with 10% fetal bovine serum (FBS; Hyclone, Logan, UT) and 1% penicillin/streptomycin (Invitrogen) in a humidified 5% CO2 incubator. Cells were treated with H2O2 (50 μmol/L, 200 μmol/L), NAC (10 μmol/L), and melatonin (10− 7 mol/L).
To detect the abundance and distribution of actin filaments, fluorescein isothiocyanate labeled.
Phalloidin (FITC-phalloidin, Sigma) was used. FITC-phalloidin was prepared in accordance with the manufacturer’s instructions. Briefly, FITC-phalloidin was dissolved as a stock solution (0.1 mg/mL) in dimethyl sulfoxide and stored at − 20 °C. The stock solution was later diluted to the working concentration (5 μg/mL) in PBST (0.2% Triton-X100 in PBS) before use. Human endometrial cancer cells (Ishikawa line) were plated on glass coverslips and fixed with 3.7% formaldehyde for 10 min at room temperature, and then permeabilized in PBST for 20 min at room temperature. After washing three times with PBS at 37 °C for 5 min each, the cells were incubated with FITC- phalloidin overnight at 4 °C. After washing three times with PBS at 37 °C for 5 min, cells were counterstained with 4′,6-diamidino-2-phenylindole (DAPI; Vector Laboratories, Burlingame, CA, USA) for 10 min and mounted on glass-bottomed culture dishes (Wuxi Nest Biotechnology Co, Ltd., Jiangsu, China) with Vectashield mounting medium (Vector Laboratories). Fluorescent signals were acquired on an upright microscope (BX51; Olympus, Tokyo, Japan) using an attached digital microscope camera (DP72; Olympus).
Total protein detection
The total protein concentration of the IVO and IVF conceptuses, C areas, or IC areas was quantified using an enhanced BCA protein assay kit (Beyotime Biotechnology, Jiangsu, China), and normalized by the weight of the tissues, providing the total protein concentration per gram of tissue.
Determination of ATP content
The ATP levels in the IVF and IVO conceptuses were detected using an Enhanced ATP Assay Kit S0027 (Beyotime Biotechnology) according to the protocol provided by the manufacturer. In brief, 20 mg of tissues were lysed in ATP lysis buffer, centrifuged for 5 min at 4 °C and 12,000×g, and the supernatant was collected. A portion of the supernatant was used to detect the ATP concentration, and the other portion of the supernatant was used to detect the total protein concentration. Finally, the total ATP concentration was normalized by the total protein concentration.
To facilitate data analysis, all proteins were mapped to the Ensembl Bos taurus gene ID. The protein quantification values of the IVO/IVF conceptuses, C areas, and IC areas are shown in Additional file 2: Table S1. Student’s t-test was used to detect the significance of the differentially expressed proteins (DEPs), and P < 0.1 was considered significant, which would avoid removing putative candidates [43, 44]. Data are presented as mean values ± SEM. DAVID version 6.7 (http://david.abcc.ncifcrf.gov/) enables the generation of specific functional annotations of biological processes affected by treatment from the target gene lists produced in high-throughput experiments [45, 46]. We used DAVID to carry the gene-annotation enrichment analysis. Then, QuickGO (https://www.ebi.ac.uk/QuickGO/) was used to identify the DEPs involved in some enriched processes in DAVID’s GO annotation . Moreover, gene symbols corresponding to DEPs were also sent to the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING version 9; http://string.embl.de/ ) to build a network using edge information from three separate forms of evidence: Databases, experiments, and text mining. We used 0.4 (medium confidence), the default and recommended option to construct protein-protein interaction networks by the STRING on-line platform, as the value for edge confidence provided by STRING. To assess the similarities of the different replicates, and to obtain a visual understanding of the relationship between the different experimental groups, we used REVIGO (http://revigo.irb.hr/) to summarize long lists of Gene Ontology terms and visualized the remaining terms in interactive graphs . Phenotype annotations of DEPs were analyzed based on the MGI database (Mouse Genome Informatics, http://www.informatics.jax.org/phenotypes.shtml). The CLUSTER 3.0 data analysis tool was used to carry out hierarchical clustering based on the clusters of protein expression profiles of different technical and biological replicates. Unsupervised hierarchical clustering analysis was carried using the “hclust” function in R (R version 3.5.1. https://www.R-project.org/.). The P-value of Student’s t-test was calculated by GraphPad Prism software or R for individual analysis.
Results and discussion
The proteome of conceptuses and endometria from the IVO and IVF groups
The experimental design is presented in Fig. 1 and Additional file 1: Figure S1A. To profile the IVF-associated proteome in conceptuses and their matched endometria at the implantation stage, IVO and IVF blastocysts were transferred to synchronized recipient ewes. On day 17 of pregnancy, which is the end point of the peri-implantation period [50, 51], and is frequently selected to explore the mechanisms of embryo-maternal crosstalk at the implantation stage in ewes [16, 52], filamentous conceptuses and their matched endometrial C areas and IC areas were sampled from each pregnant ewe. In both the IVO and IVF groups, the collected samples were divided into three pools for further proteomic analysis (Additional file 1: Figure S1A). Overall, using LC-ESI-MS/MS, we successfully identified 6374, 7495, 7933, 6162, 7401 and 8456 peptides in IVO-conceptuses, IVO- C areas, IVO-IC areas, IVF-conceptuses, IVF- C areas and IVF-IC areas, respectively. The consistence among biological replicates was evaluated by determining Pearson’s correlation coefficients using the summed peptide intensity values for each protein in the conceptuses, C areas, and IC areas. We found that the Pearson’s correlation coefficient was over 0.8 (Additional file 1: Figure S1B), indicating a general consistency in sample preparation and detection. In addition, as expected, unsupervised hierarchical clustering of the protein intensity profiles also revealed that the endometrial C and IC areas clustered closely together, and the cluster expanded to the conceptuses in IVO or IVF groups (Additional file 1: Figure S1C). In addition, the samples from the IVO and IVF groups clustered separately (Fig. 1b).
Comparative analysis of differentially expressed proteins (DEPs) between the IVF and IVO groups showed that the levels of 207, 295, and 403 proteins were significantly changed in the conceptuses, endometrial C areas, and IC areas, respectively (Additional file 3: Table S2). Notably, we found some proteins were specifically expressed in IVF conceptuses (e.g., CHD8, TSTA3), i.e., they were aberrantly activated in IVF conceptuses; while some proteins were specifically expressed in IVO conceptuses (e.g., NHLRC2, SRI), i.e., they were deficient in IVF conceptuses. (Fig. 1c). Similar results were also observed in the IVF endometrial samples (Additional file 1: Figure S1D). Compared with the IVO conceptuses, 130 DEPs were upregulated and 77 DEPs were downregulated in the IVF conceptuses. In contrast, the comparative analysis of DEPs between the IVO and IVF endometrial samples indicated that a greater proportion of the DEPs were downregulated in the IVF C or IC areas: 133 DEPs were upregulated and 162 DEPs were downregulated in the IVF C areas relative to the IVO C areas; 86 DEPs were upregulated and 137 DEPs were downregulated in the IVF IC areas relative to the IVO IC areas. (Fig. 1d). We further divided the DEPs into subcategories based on their fold changes (FC), and noticed that the proportions of dramatically changed (FC > 4) DEPs in the IVF conceptuses were much higher than those in the IVF endometrial samples, implying that IVF-induced aberrations are more dramatic in embryos and the subsequent changes in the endometria are milder but more diverse (Fig. 1e). Next, the endometrium could be considered as an early sensor of embryo implantation, therefore, we focused on the 72 DEPs common to C and IC areas using Venn diagrams. The dysregulated functions in IVF endometria were related to “cellular response to stress”, “protein metabolism”, “cytoskeleton organization”, “cell proliferation”, “genetic information processing” and “energy metabolism” (Fig. 1f). Furthermore, we found seven common DEPs among the conceptuses, C areas and IC areas, and their functions were involved in “protein metabolism”, “membrane organization”, and “cell redox homeostasis” (Fig. 1f), which indicated that the IVF-ET process might affect these physiological processes in the conceptuses and endometria. Next, functional analysis was performed for the DEPs in the different samples.
Enhanced energy metabolism, over-proliferation, and depressed methyl metabolism in IVF conceptuses
To gain further insights into the IVF-induced complications in conceptuses by implantation stage, Gene Ontology (GO) annotation using the DEPs between the IVO and IVF conceptuses was performed (Fig. 2a). We found that IVF-induced DEPs were significantly enriched in biological processes related to metabolism of carbohydrates, lipids, and amino acids, especially the GO terms of energy metabolism, such as hexose metabolism, glycolytic process and cellular carbohydrate metabolism. GO terms related to nucleic acid metabolism, such as RNA metabolic process and nucleobase, nucleoside and nucleotide metabolic process, were also enriched, implying the possibility that genetic information processing might be disrupted in the IVF conceptuses. Next, we performed REVIGO analysis to visualize the interactive relationship among the enriched terms. The results showed that metabolism of amino acids, carbohydrates, and nucleic acids were more closely related (Additional file 1: Figure S2A).
Based on these suggestions, we extracted DEPs responsible for mitochondrial functions and glycolysis, which account for the main proportion of energy metabolism of embryos by implantation stage [53,54,55]. A heat map showed that these DEPs were more abundant in IVF conceptuses. This is in line with the detection of the normalized total ATP concentrations in the IVF and IVO conceptuses, which showed that the total ATP concentration of IVF group was significantly higher than that of the IVO group, indicating that the IVF conceptuses might enhance the energy metabolism. (Fig. 2b). The enhanced energy metabolism led us speculate that IVF conceptuses might undergo increased proliferation, since metabolism is a critical determinant for proliferation during implantation development [56,57,58]. To test this, we analyzed the DEPs involved in cellular proliferation based on the gene list provided by QuickGO. As expected, IVF conceptuses showed the upregulation in metabolism and cell cycle-related proteins, which might result in enhanced proliferation in IVF conceptuses (Fig. 2c). Phenotypic analysis using the MGI database showed that some of these DEPs are important for early embryonic development (Additional file 1: Figure S2C). These results are partially in agreement with our previous observations of disrupted energy metabolism in IVF embryos . However, this disruption seems to vary among species: In mice, genes involved in mitochondrial energy metabolism were likely to be inhibited in IVF embryos by the implantation stage , which was associated with decreased fetal weight throughout gestation [13, 26, 27]. By contrast, proteins associated with energy metabolism and proliferation were upregulated in sheep IVF conceptuses. This might explain the distinct neonatal phenotypes among species following IVF-ET, i.e., intrauterine growth restriction (IUGR) and low birth weight in humans and rodents [59,60,61,62], and large offspring syndrome (LOS) in ruminants [63, 64].
Next, using DEPs from the conceptuses as seed nodes, we constructed interaction networks (Additional file 1: Figure S2B). In addition to the enriched terms of macromolecular metabolism, we also identified that the abundance of SUGT1, DNASE2, and TXN and other redox homeostasis-related proteins were upregulated in IVF conceptuses (Additional file 1: Figure S2E), which might be an adaptive change in response to the enhanced metabolism. Our recently published studies have demonstrated that energy metabolism, especially that by mitochondria, is the major source of reactive oxygen species (ROS), and IVF blastocysts are characterized by increased oxidative stress [13, 65]. The present results, using conceptuses at the implantation stage, suggested that redox homeostasis might be consistently dysregulated in IVF embryos during early development.
DNA methylation dynamics is a prominent epigenetic hallmark of early development. Our own work , as well as other related studies , suggested that the IVF process disrupted the establishment or maintenance of DNA methylation, caused by the inhibited expression of DNA methyltransferases (DNMTs). In the present study, we screened for enzymes involved in DNA methylation modification. Although no detectable changes were observed in the level of DNMTs, we found that the abundance of methyl metabolism-related proteins, such as SHMT2, MTHFD1, AHCY, and METTL7A, which catalyze key steps of methyl metabolism and transfer, decreased in the IVF conceptuses (Fig. 2d), For example, folic acid is the methyl donor of s-adenosylmethionine (SAM), and SAM is the unique active methyl donor of the DNA methylation process. The process is involved in three metabolic cycles: The folate cycle, the methionine cycle, and the methylation/demethylation cycle. In our data, the abundance of MTHFD1 and SHMT2 was reduced in IVF conceptuses, which might decrease the amount of —CH3 generated from the folate cycle. This is in accord with our previous observation that the one carbon pool by folate pathway was disrupted in IVF mouse embryos . The expression of AHCY, which plays a role in the methionine cycle, also decreased, which might decrease the amount of —CH3 provided by the methionine cycle. Meanwhile, the decreased abundance of METTL7A might lead to a reduced source of 5-methyl-cytosine (Fig. 2e). Among these enzymes, MTHFD1 and AHCY were annotated with phenotypes such as “embryonic lethality prior to/during organogenesis”, “abnormal neural tube closure” and “impaired somite development” (Additional file 1: Figure S2D), implying that the downregulated proteins associated with methyl metabolism might contribute to the impaired development of IVF conceptuses, such as increased early embryonic lethality and impaired fetal neural development, as we reported previously . The present findings, together with previous observations [8, 26], indicated that not only the methyl transfer reactions, but also one-carbon cycle-mediated methyl metabolism, were depressed by IVF processes. These results might explain why embryos or offspring following IVF-ET are associated with global hypo-methylation or loss of imprinting [66,67,68], and provide a potential strategy to rescue the impaired DNA methylation modifications in IVF embryos by supplementing substrate or precursor of methyl metabolism [69, 70].
Impaired endometrial remodeling and dysregulated redox homeostasis in IVF endometrial C areas
Next, we investigated if changed embryonic development following IVF processes induces an aberrant endometrial response. Functional profiling was performed using DEPs between the IVO and IVF C areas (Fig. 3a). GO annotation showed a similar enrichment pattern to the IVF conceptuses: Energy metabolism, amino acid metabolism, and RNA metabolism-related terms, were significantly represented. These observations were further confirmed by protein-protein interaction network constructions, in which DEPs involved in mitochondrial functions and translation processes were tightly clustered (Additional file 1: Figure S3A). In addition, REVIGO analysis further suggested that GO terms of protein metabolism, cellular amine metabolism, RNA metabolism and translation were more closely related (Additional file 1: Figure S3B). Based on these considerations, we screened DEPs involved in mitochondrial functions and translation processes, because energy metabolism and protein synthesis are primary aspects of endometrial remodeling during pregnancy [34, 71, 72]. Detailed analysis showed that a greater proportion of the DEPs associated with mitochondrial functions and translation were downregulated in IVF endometrial C areas (Fig. 3b). Based on the observation of downregulated translation-related proteins, we first compared the total abundance of the 1548 proteins detected by LC-ESI-MS/MS between the IVO and IVF C areas, and found that the total abundance of the 1548 identified proteins was decreased in each replicate of the IVF endometrial C samples (Fig. 3d). Next, we measured the total protein contents of the collected samples, the results of which showed that significantly lower amounts of total protein were detected in the IVF endometrial C areas, compared with that in their IVO counterparts (Fig. 3e). In addition, we observed the lower abundance of inducing proliferation related-proteins in IVF C areas than in IVO C areas, including MDH2, ALDH2, ALDH7A1, PGK1, ALDOA, EEF1A1, EIF4G2, ETF1, and PCNA (Fig. 3c). This implied that the proliferation of IVF endometrial C areas might be inhibited. Endometrial proliferation, the hallmark of the remodeling response to conceptus-derived signals and maternal hormone signals, is essential to establish endometrial receptivity by the implantation stage in different species [34, 73]. The important role of DEPs associated with these terms in impaired endometrial receptivity could be revealed by MGI annotation: ALDH2 and EPRS were annotated with “abnormal embryo size”; EPRS, EIF4G2, ARAS, and KARS were annotated with “embryonic lethality prior to/during organogenesis” (Additional file 1: Figure S3C).
The proposed hypothesis of impaired endometrial response was also supported by the expression patterns of proteins encoded by interferon-induced genes. In the IVF endometrial C areas, a proportion of these proteins were not upregulated (Fig. 3f). In ruminant species, interferon τ (IFNτ), secreted by the trophectoderm, is the primary signal for pregnancy recognition [74, 75]. IFNτ acts on the endometria to reduce the expression of the estrogen receptor and oxytocin receptor, thus suppressing the uterine luteolytic mechanism, and maintaining progesterone secretion to provide a receptive environment for conceptuses development . In the present study, the protein abundance of BOLA-NC1 (non-classical MHC class I antigen) was reduced by 29 fold. In humans, HLAG (non-classical MHC-1 molecule) plays an important role in embryo implantation and acts as critical tolerogenic mediator for embryo-maternal crosstalk . In cattle, BOLA-NC1 might play a role in early embryo survival and embryo-maternal immune tolerance through interacting with natural killer (NK) cells [77, 78]. The protein abundance of STAT3 (signal transducer and activator of transcription 3) was also downregulated. In mouse models with mutant STAT3 , or chemically inhibited STAT3 signaling , endometrial receptivity was significantly impaired, leading to implantation failure or embryonic lethality immediately after implantation. Other proteins encoded by interferon-induced genes that are essential for successful implantation in ewes , such as IFIT1, JSP-1, and MX2, also showed the lower abundance in the IVF C areas.
Among the enriched terms and pathways identified using functional profiling, we also noted cellular homeostasis and cytoskeleton. Cellular homeostasis is very importance to maintain a relative stable intracellular environment, including PH, media composition, and oxygen. Our data suggested that the levels of certain key enzymes of the glutathione/glutathione peroxidase (GSH/GPX) system were reduced in the IVF C areas (Fig. 3g). Previous studies in mice, ruminants, and humans reported an essential role of GSH/GPX in maintaining endometrial redox homeostasis by protecting against oxidative stress . The inhibited enzymes of the GSH/GPX system, together with the aberrant expression of cytoskeleton-related proteins in the IVF C areas, led us to ask if these complications are functionally associated. To test this, we detected the cytoskeleton organization of in vitro cultured human endometrial cancer cells (Ishikawa line) under the chemical-induced oxidative stress. A previous report showed that H2O2 impaired cytoskeleton organization in a dose-dependent manner . Similarly, our data indicated that 50 μmol/L H2O2 exposure significantly disrupted the F-actin organization, revealed by decreased fluorescence intensity of FITC-phalloidin staining. In addition, 200 μmol/L H2O2 exposure lead to a severely degenerated F-actin organization, and necrotic morphology. By supplementation with NAC at 10 μmol/L, the rate-limiting precursor for GSH synthesis, the severe F-actin damage induced by 200 μmol/L H2O2 exposure was rescued. By contrast, the rescue effect was not observed after supplementation with melatonin at 10− 7 mol/L, the physiological dosage for scavenging ROS  (Fig. 3h). These results suggested that GSH might play a unique role in protecting endometrial redox homeostasis that cannot be substituted by other antioxidants, thus the impaired GSH/GPX system might contribute largely to the aberrant remodeling of the IVF endometrial C areas.
Disrupted protein homeostasis and impaired ECM interaction in IVF endometrial IC areas
Successful implantation in ruminants depends on synergistic functions of the endometrium C and IC areas. The C areas of the endometrium are the sites of superficial attachment and placentation in ewes, while the IC areas contain large numbers of branched and coiled uterine glands that synthesize, secrete and transport a variety of molecules essential to the development of the conceptuses [22, 25]. Thus, we next compared the proteomic profiles between the IVO and IVF endometrial IC areas. As expected, we found a series of IVF-associated changes in the IC areas that differed from those observed in the C areas. Based on the functional annotations by GO, REVIGO, and STRING network construction, a cascade of terms involved in protein synthesis/degradation were enriched, such as “translation”, “protein metabolic process”, “cellular macromolecule synthetic process”, “small ribosome subunit”, and “proteasome complex”. (Fig. 4a, Additional file 1: Figure S4A, and S4B). Then we screened the expression patterns of proteins associated with these terms. The heat map indicated that many DEPs functionally associated with translation and degeneration of proteins were dysregulated. Notably, a greater proportion of DEPs associated with proteasome systems were downregulated in the IVF IC areas (Fig. 4b), suggesting disrupted protein homeostasis in IVF endometrial IC areas. This observation is in line with the result of the comparison of the total abundance of 1611 proteins detected by LC-ESI-MS/MS between the IVO and IVF C areas (Fig. 4c), as well as the BCA protein assay of total expressed proteins (Fig. 4d).
Furthermore, an impaired response to IFNτ was observed in the IVF IC areas, similar to that observed in the IVF C areas. Many proteins encoded by interferon-induced genes were not upregulated, or showed a decreasing tendency in the IVF IC areas (Fig. 4e). In sheep, MX2 expression rapidly increased in response to IFNτ-inducing by implantation, and was thought to regulate the immune system [19, 83]. In addition, JSP.1 and IFI35 are involved in immune response. JSP.1 is related to presentation of foreign antigens to the immune system. The lower abundance of MX2, JSP.1, and IFI35 in the IVF C areas might lead to dysregulated endometrial immune remodeling, which is essential to prepare future maternal immune tolerance , establish endometrial receptivity and the growth of the conceptus by the implantation stage [85, 86].
In addition, the impaired response of IC areas to the signals from the conceptuses was also in line with the downregulated proteins related to extracellular matrix (ECM) organization, proliferation, and energy metabolism (Fig. 4f, Additional file 1: Figure S4C). It has been well-documented that in ruminants, endometrial glands in the IC areas undergo extensive hyperplasia and hypertrophy during early pregnancy, presumably to meet the increasing demands of the developing conceptus for uterine histotroph [87, 88]. Proliferation and ECM interaction are functionally associated and support endometria remodeling . In addition, the important role of ECM proteins in successful implantation were also reported in mice  and humans , as revealed by influencing a series of cellular behaviors that are essential for implantation, e.g., cell migration, cell growth, cell survival, cell proliferation, angiogenesis, and invasion . Our previous works also indicated that the lower abundance of ECM proteins is associated with pregnancy loss caused by the poor endometrial receptive state .
We profiled the proteome of the IVF conceptuses and their matched endometria, aiming at understanding the mechanism of IVF-induced aberrant embryo-maternal crosstalk during early pregnancy. By functionally profiling the IVF conceptuses, we found that DEPs related to energy metabolism and proliferation were upregulated in IVF conceptuses, which might indicate the enhanced proliferation in the IVF conceptuses. This might be explained by the compensation effect that occurs in IVF embryos. Indeed, the metabolic compensation was thought to be causatively associated with LOS in IVF ruminant offspring . Interestingly, the proteins related to one-carbon cycle-mediated methyl metabolism, which provide methyl groups for methyl transfer reactions of DNA methylation, were downregulated in IVF embryos, which might indicate impaired methyl metabolism. Considering the important role of DNA methylation modification in supporting embryonic or fetal development, impaired methyl metabolism might compromise IVF embryos.
By contrast, the functional profiling of IVF-matched endometria showed the aberrant expression of proteins related to energy metabolism, proliferation, cytoskeleton organization, protein hemostasis, EMC interaction, and the antioxidant system, all of which are essential to establish endometrial receptivity. More importantly, the IVF-matched endometria displayed the reduced abundance of proteins encoded by interferon-induced genes, which might result in a reduced response to pregnancy recognition signals in both the C and IC areas.
Based on these findings, we proposed a model for the disruption of embryo-maternal crosstalk in IVF-ET (Fig. 5). In ruminants, IVF-ET embryos would undergo metabolic and proliferative compensation by upregulating related protein expression by the implantation stage, which might be causatively associated with the LOS phenotype from mid-gestation to the perinatal stage. However, the abnormal state of IVF embryos might further lead to an impaired endometrial response to pregnancy recognition signals, which is essential for successful implantation. Depressed endometrial receptivity would further impair fetal development following implantation. This concept suggests that in addition to enhancing IVF embryo quality by optimizing culture systems, improving the maternal uterine environment will also benefit the pregnancy outcome following IVF-ET, thereby providing a practical strategy in both animal reproductive management and clinical ART.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
Automatic gain control
Assisted reproductive technology
Bovine serum albumin
- C areas:
Collision induced dissociation
Differentially expressed proteins
Epidermal growth factor
Fetal bovine serum
False discovery rate
Fluorescein isothiocyanate labeled
Follicle stimulating hormone
High performance liquid chromatography
- IC areas:
Intrauterine growth restriction
In vitro fertilization and embryo transfer
In vitro maturation
Liquid chromatography-electrospray ionization-tandem mass spectroscopy
Large offspring syndrome
Penicillin and streptomycin
Reactive oxygen species
Synthetic oviductal fluid
Extracted ion currents
Mansouri-Attia N, Sandra O, Aubert J, Degrelle S, Everts RE, Giraud-Delville C, et al. Endometrium as an early sensor of in vitro embryo manipulation technologies. Proc Natl Acad Sci. 2009;106:5687–92.
Achache H, Revel A. Endometrial receptivity markers, the journey to successful embryo implantation. Hum Reprod Update. 2006;12:731–46.
Diedrich K, Fauser BCJM, Devroey P, Griesinger G. Evian annual reproduction (EVAR) workshop group. The role of the endometrium and embryo in human implantation. Hum Reprod Update. 2007;13:365–77.
Bauersachs S, Ulbrich SE, Zakhartchenko V, Minten M, Reichenbach M, Reichenbach H-D, et al. The endometrium responds differently to cloned versus fertilized embryos. Proc Natl Acad Sci U S A. 2009;106:5681–6.
Emiliani S, Delbaere A, Devreker F, Englert Y. Embryo-maternal interactive factors regulating the implantation process: implications in assisted reproductive. Reprod BioMed Online. 2005;10:527–40.
Holm P, Walker SK, Seamark RF. Embryo viability, duration of gestation and birth weight in sheep after transfer of in vitro matured and in vitro fertilized zygotes cultured in vitro or in vivo. Reproduction. 1996;107:175–81.
Kruip TAM, den Daas JHG. In vitro produced and cloned embryos: effects on pregnancy, parturition and offspring. Theriogenology. 1997;47:43–52.
Ptak GE, D’Agostino A, Toschi P, Fidanza A, Zacchini F, Czernik M, et al. Post-implantation mortality of in vitro produced embryos is associated with DNA methyltransferase 1 dysfunction in sheep placenta. Hum Reprod. 2013;28:298–305.
Tan K, An L, Miao K, Ren L, Hou Z, Tao L, et al. Impaired imprinted X chromosome inactivation is responsible for the skewed sex ratio following in vitro fertilization. Proc Natl Acad Sci. 2016;113:3197–202.
Driver AM, Peñagaricano F, Huang W, Ahmad KR, Hackbart KS, Wiltbank MC, et al. RNA-Seq analysis uncovers transcriptomic variations between morphologically similar in vivo- and in vitro-derived bovine blastocysts. BMC Genomics. 2012;13:118.
Schieve LA, Meikle SF, Ferre C, Peterson HB, Jeng G, Wilcox LS. Low and very low birth weight in infants conceived with use of assisted reproductive technology. N Engl J Med. 2002;346:731–7.
Tan K, An L, Wang SM, Wang XD, Zhang ZN, Miao K, et al. Actin disorganization plays a vital role in impaired embryonic development of in vitro-produced mouse preimplantation embryos. PLoS One. 2015;10:1–15.
Ren L, Wang Z, An L, Zhang Z, Tan K, Miao K, et al. Dynamic comparisons of high-resolution expression profiles highlighting mitochondria-related genes between in vivo and in vitro fertilized early mouse embryos. Hum Reprod. 2015;30(12):2892–911 dev228.
Tian X, Zhang L, Tan D, Wang F, Liu G, Reiter RJ. Melatonin promotes the in vitro development of pronuclear embryos and increases the efficiency of blastocyst implantation in murine. J Pineal Res. 2013;55:267–74.
Ali AA, Bilodeau JF, Sirard MA. Antioxidant requirements for bovine oocytes varies during in vitro maturation, fertilization and development. Theriogenology. 2003;59:939–49.
Bazer FW, Spencer TE, Ott TL. Interferon tau: A novel pregnancy recognition signal. Am J Reprod Immunol. 1997;37:412–20.
Macklon NS, Brosens JJ. The human endometrium as a sensor of embryo Quality1. Biol Reprod. 2014;91:1–8.
Lundeberg T, Parasassi T, Pittaluga E, Brunelli R. N-acetyl-L-cysteine for use in in vitro fertilization; 2017.
Gray CA, Abbey CA, Beremand PD, Choi Y, Farmer JL, Adelson DL, et al. Identification of endometrial genes regulated by early pregnancy, progesterone, and interferon tau in the ovine Uterus1. Biol Reprod. 2006;74:383–94.
Koch JM, Ramadoss J, Magness RR. Proteomic profile of uterine luminal fluid from early pregnant ewes. J Proteome Res. 2010;9:3878–85.
Wang Y, Wang C, Hou Z, Miao K, Zhao H, Wang R, et al. Comparative analysis of proteomic profiles between endometrial caruncular and intercaruncular areas in ewes during the peri-implantation period. J Anim Sci Biotechnol. 2013;4:39.
Zhao H, Sui L, Miao K, An L, Wang D, Hou Z, et al. Comparative analysis between endometrial proteomes of pregnant and non-pregnant ewes during the peri-implantation period. J Anim Sci Biotechnol. 2015;6:1–14.
Mamo S, Mehta JP, Forde N, McGettigan P, Lonergan P. Conceptus-endometrium crosstalk during maternal recognition of pregnancy in cattle. Biol Reprod. 2012;87(6):1–9.
Berendt FJ, Fröhlich T, Schmidt SEM, Reichenbach H-D, Wolf E, Arnold GJ. Holistic differential analysis of embryo-induced alterations in the proteome of bovine endometrium in the preattachment period. Proteomics. 2005;5:2551–60.
Mansouri-Attia N, Aubert J, Reinaud P, Giraud-Delville C, Taghouti G, Galio L, et al. Gene expression profiles of bovine caruncular and intercaruncular endometrium at implantation. Physiol Genomics. 2009;39:14–27.
Sui L, An L, Tan K, Wang Z, Wang S, Miao K, et al. Dynamic proteomic profiles of in vivo- and in vitro-produced mouse Postimplantation Extraembryonic tissues and Placentas1. Biol Reprod. 2014;91:1–16.
Nie J, An L, Miao K, Hou Z, Yu Y, Tan K, et al. Comparative analysis of dynamic proteomic profiles between in vivo and in vitro produced mouse embryos during postimplantation period. J Proteome Res. 2013;12:3843–56.
Chae J-I, Kim J, Lee SG, Jeon Y-J, Kim D-W, Soh Y, et al. Proteomic analysis of pregnancy-related proteins from pig uterus endometrium during pregnancy. Proteome Sci. 2011;9:41.
Rashid NA, Lalitkumar S, Lalitkumar PG, Gemzell-Danielsson K. Endometrial receptivity and human embryo implantation. Am J Reprod Immunol. 2011;66(Suppl 1):23–30.
Salamonsen LA, Edgell T, Rombauts LJF, Stephens AN, Robertson DM, Rainczuk A, et al. Proteomics of the human endometrium and uterine fluid: a pathway to biomarker discovery. Fertil Steril. 2013;99:1086–92.
Lee KY, DeMayo FJ. Animal models of implantation. Reproduction. 2004;128:679–95.
Barry JS, Anthony RV. The pregnant sheep as a model for human pregnancy. Theriogenology. 2008;69:55–67.
Spencer TE, Johnson GA, Bazer FW, Burghardt RC. Implantation mechanisms: insights from the sheep. Reproduction. 2004;128:657–68.
Igwebuike UM. A review of uterine structural modifications that influence conceptus implantation and development in sheep and goats. Anim Reprod Sci. 2009;112:1–7.
Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, et al. Erratum: corrigendum: global quantification of mammalian gene expression control. Nature. 2013;495:126–7.
De Matos DG, Gasparrini B, Pasqualini SR, Thompson JG. Effect of glutathione synthesis stimulation during in vitro maturation of ovine oocytes on embryo development and intracellular peroxide content. Theriogenology. 2002;57:1443–51.
De Matos DG, Furnus CC, Moses DF, Martinez AG, Matkovic M. Stimulation of glutathione synthesis of in vitro matured bovine oocytes and its effect on embryo development and freezability. Mol Reprod Dev. 1996;45:451–7.
Graumann J, Hubner NC, Kim JB, Ko K, Moser M, Kumar C, et al. Stable isotope labeling by amino acids in cell culture (SILAC) and proteome quantitation of mouse embryonic stem cells to a depth of 5,111 proteins. Mol Cell Proteomics. 2008;7(4):672–83.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26(12):1367–72.
Li GZ, JPC V, Silva JC, Golick D, Gorenstein MV, Geromanos SJ. Database searching and accounting of multiplexed precursor and product ion spectra from the data independent analysis of simple and complex peptide mixtures. Proteomics. 2009;9(6):1696–719.
Feng J, Naiman DQ, Cooper B. Probability-based pattern recognition and statistical framework for randomization: modeling tandem mass spectrum/peptide sequence false match frequencies. Bioinformatics. 2007;23:2210–7.
Waanders LF, Chwalek K, Monetti M, Kumar C, Lammert E, Mann M. Quantitative proteomic analysis of single pancreatic islets. Proc Natl Acad Sci U S A. 2009;106:18902–7.
Kim Y, Jeon J, Mejia S, Yao CQ, Ignatchenko V, Nyalwidhe JO, et al. Targeted proteomics identifies liquid-biopsy signatures for extracapsular prostate cancer. Nat Commun. 2016;7:1–10.
Lapek JD, Greninger P, Morris R, Amzallag A, Pruteanu-Malinici I, Benes CH, et al. Detection of dysregulated protein-association networks by high-throughput proteomics predicts cancer vulnerabilities. Nat Biotechnol. 2017;35(10):983–9.
Huang DW, Sherman BT, Tan Q, Kir J, Liu D, Bryant D, et al. DAVID Bioinformatics Resources: Expanded annotation database and novel algorithms to better extract biology from large gene lists. Nucleic Acids Res. 2007;35(Web Server issue):W169–75.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Huntley RP, Sawford T, Mutowo-Meullenet P, Shypitsyna A, Bonilla C, Martin MJ, et al. The GOA database: Gene Ontology annotation updates for 2015. Nucleic Acids Res. 2015;43(Database issue):D1057–63.
Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, et al. STRING v9.1: Protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res. 2013;41(Database issue):D808–15.
Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800.
Wan P-C, Bao Z-J, Wu Y, Yang L, Hao Z-D, Yang Y-L, et al. αvβ3 integrin may participate in Conceptus attachment by regulating morphologic changes in the endometrium during Peri-implantation in ovine. Reprod Domest Anim. 2011;46:840–7.
Spencer T, Burghardt R, Johnson G, Bazer F. Conceptus signals for establishment and maintenance of pregnancy. Anim Reprod Sci. 2004;82–83:537–50.
Song G, Satterfield MC, Kim J, Bazer FW, Spencer TE. Gastrin-releasing peptide (GRP) in the ovine uterus: regulation by interferon tau and Progesterone1. Biol Reprod. 2008;79(2):376–86.
Zhou W, Choi M, Margineantu D, Margaretha L, Hesson J, Cavanaugh C, et al. HIF1α induced switch from bivalent to exclusively glycolytic metabolism during ESC-to-EpiSC/hESC transition. EMBO J. 2012;31:2103–16.
Leese H. Metabolic control during preimplantation mammalian development. Hum Reprod Update. 1995;1:63–72.
Harvey AJ. The role of oxygen in ruminant preimplantation embryo development and metabolism. Anim Reprod Sci. 2007;98:113–28.
Krisher RL, Prather RS. A role for the Warburg effect in preimplantation embryo development: metabolic modification to support rapid cell proliferation. Mol Reprod Dev. 2012;79:311–20.
Folmes CDL, Terzic A. Metabolic determinants of embryonic development and stem cell fate. Reprod Fertil Dev. 2015;27:82–8.
Ramalho-Santos J, Varum S, Amaral S, Mota PC, Sousa AP, Amaral A. Mitochondrial functionality in reproduction: from gonads and gametes to embryos and embryonic stem cells. Hum Reprod Update. 2009;15:553–72.
Moise J, Laor A, Armon Y, Gur I, Gale R. The outcome of twin pregnancies after IVF. Hum Reprod. 1998;13:1702–5.
Pinborg A, Loft A, Schmidt L, Langhoff-Roos J, Andersen AN. Maternal risks and perinatal outcome in a Danish national cohort of 1005 twin pregnancies: the role of in vitro fertilization. Acta Obstet Gynecol Scand. 2004;83:75–84.
Nassar AH, Usta IM, Rechdan JB, Harb TS, Adra AM, Abu-Musa AA. Pregnancy outcome in spontaneous twins versus twins who were conceived through in vitro fertilization. Am J Obstet Gynecol. 2003;189:513–8.
Cai L, Izumi S, Koido S, Uchida N, Suzuki T, Matsubayashi H, et al. Abnormal placental cord insertion may induce intrauterine growth restriction in IVF-twin pregnancies. Hum Reprod. 2006;21:1285–90.
Chen Z, Robbins KM, Wells KD, Rivera RM. Large offspring syndrome. Epigenetics. 2013;8:591–601.
Chen Z, Hagen D E, Elsik C G, et al. Characterization of global loss of imprinting in fetal overgrowth syndrome induced by assisted reproduction[J]. Proceedings of the National Academy of Sciences. 2015;112(15):4618–4623.
REN L, ZHANG C, TAO L, HAO J, TAN K, MIAO K, et al. High-resolution profiles of gene expression and DNA methylation highlight mitochondrial modifications during early embryonic development. J Reprod Dev. 2017;63:247–61.
Mann MRW. Selective loss of imprinting in the placenta following preimplantation development in culture. Development. 2004;131:3727–35.
Wright KP, Brown L, Casson PR, Johnson JV, Brown S. ART conditions are associated with global genomic hypomethylation at the blastocyst stage in mice. Fertil Steril. 2008;90:S345.
Katari S, Turan N, Bibikova M, Erinle O, Chalian R, Foster M, et al. DNA methylation and gene expression differences in children conceived in vitro or in vivo. Hum Mol Genet. 2009;18:3769–78.
Vollset SE, Gjessing HK, Tandberg A, Rønning T, Irgens LM, Baste V, et al. Folate Supplementation and Twin Pregnancies. Epidemiology. 2005;16:201–5.
Laanpere M, Altmäe S, Stavreus-Evers A, Nilsson TK, Yngve A, Salumets A. Folate-mediated one-carbon metabolism and its effect on female fertility and pregnancy viability. Nutr Rev. 2010;68:99–113.
Chen Q, Zhang A, Yu F, Gao J, Liu Y, Yu C, et al. Label-free proteomics uncovers energy metabolism and focal adhesion regulations responsive for endometrium receptivity. J Proteome Res. 2015;14:1831–42.
Włodarek J, Jaśkowski J, Nowak W, MIKUŁA JAN, Olechnowicz J. Development of ovarian follicles, quality of oocytes and fertility of cows in view of a negative energy balance in the transition period. Med Weter. 2011;67:219–23.
Paulson RJ. Hormonal induction of endometrial receptivity. Fertil Steril. 2011;96:530–5.
Ashworth CJ, Bazer FW. Changes in ovine Conceptus and endometrial function following Asynchronus embryo transfer or administration of Progesterone1. Biol Reprod. 1989;40:425–33.
Farin CE, Imakawa K, Roberts RM. In situ localization of mRNA for the interferon, ovine Trophoblast Protein-1, during early embryonic development of the sheep. Mol Endocrinol. 1989;3:1099–107.
Le Bouteiller P. Commentary major breakthrough in the HLA-G debate: occurrence of pregnancy in human depends on the HLA-G status of preimplantation embryos. Eur J Immunol. 2002;32:309–10.
O’Gorman GM. Naib a Al, Ellis SA, Mamo S, O’Doherty AM, Lonergan P, et al. regulation of a bovine nonclassical major histocompatibility complex class I gene Promoter1. Biol Reprod. 2010;83:296–306.
Araibi EH, Marchetti B, Dornan ES, Ashrafi GH, Dobromylskyj M, Ellis SA, et al. The E5 oncoprotein of BPV-4 does not interfere with the biosynthetic pathway of non-classical MHC class I. Virology. 2006;353:174–83.
Takeda K, Noguchi K, Shi W, Tanaka T, Matsumoto M, Yoshida N, et al. Targeted disruption of the mouse Stat3 gene leads to early embryonic lethality. Proc Natl Acad Sci U S A. 1997;94:3801–4.
Catalano RD, Johnson MH, Campbell EA, Charnock-Jones DS, Smith SK, Sharkey AM. Inhibition of Stat3 activation in the endometrium prevents implantation: a nonsteroidal approach to contraception. Proc Natl Acad Sci. 2005;102:8585–90.
Serviddio G, Loverro G, Vicino M, Prigigallo F, Grattagliano I, Altomare E, et al. Modulation of endometrial redox balance during the menstrual cycle: relation with sex hormones. J Clin Endocrinol Metab. 2002;87:2843–8.
Zhao Y, Davis HW. Hydrogen peroxide-induced cytoskeletal rearrangement in cultured pulmonary endothelial cells. J Cell Physiol. 1998;174:370–9.
Ott TL, Yin J, Wiley AA, Kim H-T, Gerami-Naini B, Spencer TE, et al. Effects of the estrous cycle and early pregnancy on uterine expression of mx protein in sheep (Ovis aries)1. Biol Reprod. 1998;59:784–94.
Lédée N, Prat-Ellenberg L, Chevrier L, Balet R, Simon C, Lenoble C, et al. Uterine immune profiling for increasing live birth rate: a one-to-one matched cohort study. J Reprod Immunol. 2017;119:23–30.
Hannan NJ, Evans J, Salamonsen LA. Alternate roles for immune regulators: establishing endometrial receptivity for implantation. Expert Rev Clin Immunol. 2011;7:789–802.
Ott TL, Kamat MM, Vasudevan S, Townson DH, Pate JL. Maternal immune responses to conceptus signals during early pregnancy in ruminants. Anim Reprod. 2014;11(n3):237–45.
Stewart MD, Johnson GA, Gray CA, Burghardt RC, Schuler LA, Joyce MM, et al. Prolactin receptor and uterine Milk protein expression in the ovine endometrium during the estrous cycle and Pregnancy1. Biol Rep. 2000;62(Issue 6):1779–1789.
Bazer FW. Uterine protein secretions: relationship to development of the conceptus. J Anim Sci. 1975;41(5):1376–82.
Gray CA, Bartol FF, Tarleton BJ, Wiley AA, Johnson GA, Bazer FW, et al. Developmental biology of uterine Glands1. Biol Reprod. 2001;65:1311–23.
Woessner JF. Matrix metalloproteinases and their inhibitors in connective tissue remodeling. FASEB J. 1991;5:2145–54.
Fisher SJ, Damsky CH. Human cytotrophoblast invasion. Semin Cell Biol. 1993;4:183–8.
Reddy KVR, Mangale SS. Integrin receptors: the dynamic modulators of endometrial function. Tissue Cell. 2003;35:260–73.
Niemann H, Korsawe K, Lemme E, Lucas-Hahn A, Herrmann D, Wrenzycki C. Gene expression patterns in in vitro-produced and somatic nuclear transfer-derived preimplantation bovine embryos: relationship to the large offspring syndrome? Anim Reprod Sci. 2004;82–83:593–603.
We thank the members of our laboratory for their helpful comments on the manuscript.
This work was supported by the grants from National Key R&D Program (2017YFD0501901 and 2017YFD0501905), National Natural Science Foundation of China (No.3167246 and 31972573), National Support Program for Youth Top-notch Talents, and the Earmarked Fund for the Innovative Teams of Beijing Swine Industrialization Research Program.
Ethics approval and consent to participate
The experiments were performed accordance with the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching, and all procedures were approved by the Institutional Animal Care and Use Committee at China Agricultural University (Beijing, China).
Consent for publication
The authors declare that they have no competing interests.
Additional file 1: Figure S1. (a) Strategy for samples pooling. Conceptuses are presented as an example. In the IVO (or IVF) group, the samples of 36 (or 18) conceptuses were equally and randomly divided into three pools as independent biological replicates, each pooled sample was detected twice as technical replicates. (b) Pearson’s Correlation coefficient between the randomly selected two biological replicates using the summed peptide intensity values for each protein in the conceptus, C area and IC area samples to confirm the reliability of the processes of sample collection, LC-ESI-MS/MS detection and proteomic analysis. (c) Unsupervised clustering of the protein expression patterns of each biological replicate of the conceptus, C area and IC area samples from the IVO and IVF groups. (d) Volcano plot of DEPs in the C areas, IC areas and conceptuses between the IVO and IVF groups. The red and green dots represent upregulated or downregulated DEPs, respectively (−Log10(P-value) > 1; mean fold change > 2 or < 0.5). Figure S2. (a) Graph visualization of biological processes (BPs) enriched in the IVF conceptuses based on REVIGO analysis. The bigger circle represents the more clustered genes. The gradation of color represents the significance level. Functionally associated BPs are linked. (b) Interaction networks of DEPs between the IVO and IVF conceptuses were created by a web-based search of the STRING database. Boxed regions represent tightly interconnected functional clusters. (c-d) The representative MGI phenotypes associated with embryonic development annotated with DEPs related to metabolism and cell cycle (c) and methyl metabolism (d). (e) Heat map of DEPs related to redox homeostasis in the IVO and IVF conceptuses. Normalized protein abundance is represented in red (relatively high) and green (relatively low). The table shows the representative MGI phenotypes of DEPs in redox homeostasis. Figure S3. (a) Interaction networks of DEPs between IVO and IVF conceptuses were created by a web-based search of the STRING database. Boxed regions represent tightly interconnected functional clusters (b). Graphical visualization of biological processes (BPs) enriched in IVF C areas based on REVIGO analysis. The bigger circle represents the more clustered genes. The gradation of the color represents the significance level. Functionally associated BPs are linked. (c) The representative MGI phenotypes associated with embryonic development annotated with DEPs related to metabolism and translation. Figure S4. (a) Graphical visualization of biological processes (BPs) enriched in IVF IC areas based on REVIGO analysis The bigger circle represents the more clustered genes. The gradation of the color represents the significance level. Functionally associated BPs are linked. (b) Interaction networks of DEPs between the IVO and IVF IC areas are created by a web-based search of the STRING database. Boxed regions represent tightly interconnected functional clusters. (c) Heat map of DEPs involved in Glycolysis/Gluconeogenesis in the IC areas. The level of Log2(IVF/IVO) is represented in red (relatively high) and green (relatively low).
About this article
Cite this article
Yang, Q., Fu, W., Wang, Y. et al. The proteome of IVF-induced aberrant embryo-maternal crosstalk by implantation stage in ewes. J Animal Sci Biotechnol 11, 7 (2020). https://doi.org/10.1186/s40104-019-0405-y
- Embryo-maternal crosstalk
- Endometrial C areas
- Endometrial IC areas