Abstract
During maternal recognition of pregnancy (MRP), a conceptus-derived signal leads to the persistence of the corpus luteum and the maintenance of gestation. In the horse, the nature of this signal remains to be elucidated. Several studies have focused on the changes in gene expression during MRP, but little information exists at the protein level. The aim of this study was to identify the proteins at the embryo-maternal interface around signalling of MRP in the horse (day 13) by means of mass spectrometry. A distinct influence of pregnancy was established, with 119 proteins differentially expressed in the uterine fluid of pregnant mares compared to cyclic mares and with upregulation of several inhibitors of the prostaglandin synthesis during pregnancy. By creating an overview of the proteins at the embryo-maternal interface in the horse, this study provides a solid foundation for further targeted studies of proteins potentially involved in embryo-maternal interactions, MRP and pregnancy loss in the horse.
Similar content being viewed by others
Introduction
Maternal recognition of pregnancy (MRP) covers the series of events leading to the persistence of the corpus luteum and a receptive uterine environment to support the maintenance of gestation1. In the cycling mare, pulsatile release of prostaglandin F2α (PGF2α) causes luteolysis, resulting in a decline in progesterone. This mechanism is inhibited during pregnancy by the presence of the conceptus2. In pigs, the conceptus derived signal which initiates MRP has been identified a long time ago as oestrogen3 and in ruminants as interferon tau4,5. However, the nature of this signal remains to be elucidated in the horse despite several decades of elaborate research on this topic6,7. Initial studies have focused on the identity of specific candidate signalling molecules and while the equine embryo produces substantial quantities of oestrogen as well as prostaglandins (PG) and limited amounts of interferons, no convincing evidence exists for their signalling role in MRP7. Potential embryonic signal targets involved in the luteostatic mechanism in the horse are prostaglandin-endoperoxide synthase 2 (PTGS2), an enzyme in the biosynthesis of PGF2α, and oxytocin, which stimulates endometrial PGF2α secretion through a positive feedback loop8. Both PTGS2 and oxytocin receptor expression (OXTR) are repressed during early pregnancy compared to cycling mares, with downregulation of PTGS2 at the RNA level and of OXTR at the protein level9,10,11,12,13.
During the last years, the topic of MRP in the horse has been broadened to all pathways involved in embryo-maternal communication around the timing of MRP. Signalling of MRP is a continuum of events, estimated to occur between days 12 and 14. Recipient mares can still get pregnant when an embryo is transferred to their uterus at day 12, but not at day 14 after ovulation14, while repression of PTGS2 occurs by day 13 of pregnancy11. By day 16, clear differences between pregnant and cyclic horses are observed. Transcriptomics of the equine endometrium and equine conceptuses have substantially contributed to the knowledge on pathways affected around the timing of MRP in the horse7,15,16,17,18,19. Technological advantages, including sequencing, favoured development of genomics and transcriptomics compared to proteomics20. However, mRNA abundances can only explain 40% of the variation in protein levels and the actual protein profile is influenced by post-transcriptional regulation mechanisms21. This appeals for complementing transcriptomics knowledge on MRP with quantitative proteomics. This can now be achieved through mass spectrometry (MS). Recent improvements in MS technologies, including data-independent-acquisition, allow reproducible label-free quantification of proteins in complex biological samples22.
Mass spectrometry of the embryo-maternal interface around MRP has been performed in several farm animals including pigs23,24, sheep25 and cattle26,27,28. In the horse, specific molecules with a potential role in MRP have been targeted by immunohistochemistry13,29,30,31,32 and global screening of uterine proteins has been performed in the context of endometritis33. However, the effect of pregnancy on the uterine secretome has not been assessed by means of high-throughput proteomics in the horse up to now. In a recent study, equine blastocysts were collected by uterine lavage on day 8 and an MS analysis was performed of the proteins secreted during culture of these embryos for 24 h and 48 h and of proteins present in the blastocoel fluid and the embryo capsule34. The authors detected prostaglandin F2 receptor inhibitor (PTGFRN) and a progesterone potentiating protein, FK506 binding protein 4 (FKBP4), in the blastocoel fluid, but it remained to be determined whether these proteins were actively secreted into the uterine lumen.
The aim of this study was to gain new insights into the embryo-maternal communication around the signalling of MRP in the horse. Since signalling of MRP is estimated to occur between Day 12 and Day 14, sampling was performed at Day 13 (±0.5 day). We hypothesize that high-throughput proteomics can provide complementary information to the transcriptomic reports. To this end, proteomics was performed by high definition data independent mass spectrometry (HDMSE) with ion mobility drift time-specific collision-energy35. In this way, proteins were identified and quantified in uterine fluid of pregnant and cyclic mares as well in the yolk sac fluid of the pregnant mares.
Results
Sampling
Only reproductively sound mares with negative bacteriology and cytology of the uterine fluid were used for the sampling. In two cycles, namely one pregnant (P) and one control cyclic (C) cycle, a double ovulation occurred. Response to hCG resulted in ovulation 24–36 h after administration. In four cycles, ovulation only occurred 3 days after hCG; once in a P cycle, where artificial insemination (AI) was performed at the same time and in this case, the mare was inseminated a second time 48 h after the first time and she ovulated the day after. In all other P cycles, ovulation occurred within 48 h after AI. In one mare, a line of fluid was noticed by ultrasound of the uterus 1 day after AI and she was treated by intramuscular administration of oxytocin.
Identification of proteins
The average protein concentration was similar in the uterine fluid (UF) of P (9.2 g/mL) and C (9.8 g/mL) mares, while the average protein concentration in the yolk sac (YS) was only 78 µg/mL.
For the first time, an overview was created of the proteins present in the UF and the YS at day 13 after ovulation in the horse. In the UF samples, a total of 10489 peptides were identified, accounting for 41% of all peptide like ions. Protein identification resulted in 1153 identifiable proteins (Supplementary file 1). After filtering and normalization, a total of 707 normalized proteins with at least two unique peptides were assessed for differential expression.
Differential expression of proteins was assessed for P versus C mares and pregnancy was associated with upregulation of 62 proteins (Table 1) and downregulation of 57 proteins (Table 2). For all proteins in this comparison, the log fold change, the adjusted p-value and the number of peptides are listed in Supplementary file 3.
In the YS samples, a total of 6500 peptide ions were identified, representing 51% of all peptide like ions and resulting in 903 identifiable proteins (Supplementary file 2). For the YS proteins, the primary goal was identification, rather than quantification, as different nature of the fluids impedes assessment of differential expression of proteins in YS versus UF.
Gene Ontology enrichment and pathway analysis
Categorization in the Gene Ontology (GO) terms ‘molecular function’, ‘biological process’ and ‘cellular component’ is provided for all quantified proteins in the comparison of P versus C in Supplementary file 3.
Figure 1 summarizes the GO categories in which the differentially expressed proteins are involved. The main category to which most proteins are assigned is ‘cellular process (GO:0009987)’ for the biological processes and ‘binding (GO:0005488)’ for the molecular functions. This coincides with the results in porcine uterine fluid24, but these are also the major categories when all proteins are taken into account. Overall, the differences in categorization between the groups are small.
Gene Ontology (GO) enrichment revealed no statistical overrepresentation when a Bonferroni correction for multiple testing was used (FDR < 0.05). No up- or downregulated KEGG pathways were detected either at a Benjamini-Hochberg corrected p-value of 0.05.
Embryo-maternal interaction
Comparison of the proteins identified in the UF of the P mares and in the YS of the corresponding embryo revealed 347 common proteins, 806 proteins which were only detected in the UF and 556 proteins which were only found in the YS. Figure 2 represents an overview of these UF specific proteins, YS specific proteins and common proteins, with specific display of the proteins up- and downregulated during pregnancy and of the proteins categorized in the extracellular space.
A list of the 347 common proteins is provided in Supplementary file 4, including the functions in which these proteins are involved. Figure 3 summarizes the GOs in which these common proteins were involved. Similar to the results for the UF in Fig. 1, the main GO categories in which the common proteins are involved are also ‘cellular process (GO:0009987)’ and ‘binding (GO:0005488)’ and differences in categorization are small. Common proteins in YS and UF which were also upregulated in P versus C, showed a higher representation in the biological processes ‘developmental process (GO:0032502)’ and ‘response to stimulus (GO:0050896)’. Molecular functions in which these proteins were more involved are ‘transporter activity (GO:0005215)’ and ‘transcription factor activity - protein binding (GO:0000988)’, while the common proteins which were downregulated in P versus C were rather represented in ‘structural molecule activity (GO:0005198)’.
The embryo-maternal interaction was further visualized by Cytoscape 3.3.0 in Fig. 4. The most prominent GO terms in this network are ‘embryo development’ (GO:0009790) and ‘embryo morphogenesis’ (GO:0048598), with a main contribution of proteins originating from the yolk sac, and embryo implantation (GO:0007566) with the involvement of both uterine and embryonic proteins. In Fig. 5, the contribution of growth factors and cytokines in equine embryo-maternal signalling is visualized. The most extensive networks with various proteins found in the yolk sac fluid and/or the uterine fluid of pregnant mares are involved in ‘regulation of cytokine production’ (GO:0001817), ‘response to cytokine’ (GO:0034097) and the downstream GO’s ‘cytokine receptor binding‘ (GO:0005126), ‘cytokine mediated signalling pathway’ (GO:0019221) and ‘regulation of response to cytokine stimulus’ (GO:0060759).
Discussion
Maternal recognition of pregnancy is an intriguing subject in the horse and extensive research on the molecular processes involved has been performed in the field of transcriptomics7,18. However, information on the downstream translation to proteins is scarce. In this study, quantitative proteomics of the uterine luminal fluid assessing the effect of pregnancy was performed for the first time in the horse. At the same time, proteins in the embryonic yolk sac fluid were mapped to provide insight into the embryo-maternal interaction.
With 119 proteins differentially expressed in the uterine fluid of P versus C mares, a distinct influence of pregnancy was established. In general, a function of more than 40% of the differentially expressed proteins in the UF was categorized as ‘binding (GO:0005488)’, coinciding with the findings in pigs and cattle, where the majority of proteins were also allocated to molecular binding24,28 (Fig. 1). ‘Binding’ also represents the main category to which the common proteins in UF and YS were allocated, with subtly higher representation of proteins upregulated during pregnancy in categories linked to embryo-maternal interaction, namely ‘developmental process (GO:0032502)’, ‘response to stimulus (GO:0050896)’, ‘transporter activity (GO:0005215)’ and ‘transcription factor activity - protein binding (GO:0000988)’ at the expense of the more general GO term ‘structural molecule activity (GO:0005198)’ (Fig. 3). Cellular component categorization allocated 45% of the identified UF proteins to the extracellular space (Fig. 2). This coincides with the findings of Swegen et al.34, who specifically targeted secreted proteins by analysing embryo-conditioned medium. This supports the fact the proteins detected in our study mainly represent the proteins secreted in the uterine fluid rather than endometrial cells shed in the uterine lumen. This also accounts for the proteins which were found to be differentially expressed during pregnancy. Sixty four % of these proteins were categorized in the extracellular space; the other may have originated from occasional shedding of embryonic and endometrial cells into the uterine lumen. Figure 2 represents an overview of all UF specific, YS specific and common proteins, including their differential expression in P vs C and their allocation to the extracellular space. Interestingly, the majority of proteins commonly found in UF and YS are indeed present in the extracellular space. These represent candidate proteins involved in embryo-maternal interaction and signalling. In general, our results greatly coincided with the findings of Swegen et al.34 who worked with day 8 blastocysts to examine proteins present in and secreted by early equine embryos. Figure 6 shows the number of proteins which were commonly found in the blastocoel fluid and the YS and those found to be secreted in embryo-conditioned medium at 24 h and 48 h and in the UF in our study. More than two third of the proteins reported in the blastocoel fluid were also detected in the YS and more than one third of the proteins found to be secreted after 48 h of embryo culture were also detected in the UF. Overlap of the results validates our independent findings on the one hand and indicates conserved expression of several proteins throughout development on the other hand.
In the context of MRP, prostaglandin synthesis is of special interest. For three proteins involved in this pathway, namely prostaglandin reductase 1 (PTGR1), glutathione transferase 1 (GSTP1) and annexin A1 (ANXA1), significantly higher amounts were detected in the uterine fluid of pregnant mares compared with cyclic mares. Apart from acting on 15-oxo-PGE1, 15-oxo-PGE2 and 15-oxo-PGE2-alpha as 15-oxo-prostaglandin 13-reductase, PTGR1 catalyzes leukotriene B4 into its biologically less active metabolite, being the key step in the metabolic inactivation of leukotriene B4, as depicted in Fig. 7.
While glutathione transferases are also generally involved in the biosynthesis of prostaglandins and leukotrienes, as well as progesterone and testosterone36, a specific anti-inflammatory effect of GSTP1 by reduction of PTGS2, formerly known as cyclooxygenase-2 (COX-2), has been described37. Furthermore, transport of GSTP1 across the plasma membrane was demonstrated37. Based upon these observations with recombinant human GSTP1 in mice and the high homology of equine GSTP1 with other species, equine GSTP1 in uterine fluid might cross the plasma membrane and target intracellular PTGS2. Interestingly, GSTP1 was also detected in the YS of the equine conceptuses. In this regard, the pregnancy associated upregulation of GSTP1 observed in the equine uterine fluid could be involved in the luteostatic mechanism by inhibiting PTGS2. Further research is needed to examine this hypothesis, as this is the first report on the presence of GSTP1 in equine uterine fluid.
Another anti-inflammatory factor with an inhibitory effect on prostaglandin synthesis, more specifically on phospholipase A2, is annexin A1 (ANXA1)38,39. Annexin A1 was upregulated in the uterine fluid of the pregnant mares when compared to the cyclic mares and this association of annexins with pregnancy coincides with literature. An increase in ANXA1 was also reported in the uterine luminal fluid of pregnant ewes from day 10 to day 1225. Several annexins have been linked to embryo-maternal interaction. Annexin 4 (ANXA4) was found to increase over time from day 10 to day 13 in both cyclic and pregnant pigs24, and we previously reported greater quantities of ANXA4 in the oviductal fluid of pregnant mares, when compared to cyclic mares40. In our study, we detected annexin 1, 2, 3, 4, 5, 7 8 and 11 in the UF, while ANXA2 and ANXA5 were also found in the YS. Swegen et al.34 also reported the presence of ANXA2 in both the equine blastocoel fluid and the embryo-conditioned medium after 48 h. The only annexin found to be upregulated during pregnancy was ANXA1. Annexin 1 is an inhibitor of phospholipase A2, a rate-limiting enzyme which liberates arachidonic acid for the synthesis of prostaglandins and for which a lower enzyme activity of phospholipase A2 has been demonstrated in pregnant mares compared to cyclic mares on day 1441. In our study, phospholipase A2 group IIA (PLA2G2A) tended to be downregulated in the uterine fluid of pregnant mares with a logFC of −1.31 compared to the cyclic condition, but it was not significant at a 0.05 FDR. Both PLA2G2A and phospholipase A2 group VII (PLA2G7) were detected in the YS; the latter was also found in the equine embryo-conditioned medium after 48 h34. Overall, our data suggest pregnancy associated interference with the luteolytic eicosanoid pathway with upregulation of inhibitory factors at different levels of the prostaglandin synthesis pathway.
A close interaction between prostaglandins and oxytocin has been described in the context of MRP in the horse with downregulation of the oxytocin receptor protein in the pregnant endometrium on day 1413. In the present study, the presence of oxytocin in the uterine luminal fluid was examined, but it was not detected. However, this does not mean it was not present in the original samples; the collection method might have retained some peptides and missing values are intrinsic to mass spectrometry42,43. Phosphoinositide phospholipase C (PLCD1), involved in the oxytocin receptor signalling pathway, was detected, but not significantly affected by pregnancy44. The reduced expression of OXTR during pregnancy has been hypothesized to be induced by an observed decrease in the gene expression of oestrogen receptor 1 (ESR1) in the pregnant equine endometrium15. Several proteins related to ESR1 were also found to be affected in the uterine fluid. Surprisingly, pregnant mares showed a strong upregulation of deoxynucleotidyltransferase terminal interacting protein 2 (DNTTIP2), previously known as oestrogen receptor binding protein (ERBP). Binding of DNTTIP2 to ESR1 enhances its transcription45. Upregulation of DNTTIP2, which would lead to increased transcription of ESR1 in pregnant mares is contradictory to findings in literature and further targeted research is required to clarify this aspect. Downstream of the ESR1, the influence of oestrogen on the ezrin–radixin–moesin (ERM) family of actin-binding proteins has been studied, mainly in the context of breast cancer46,47. The distribution pattern of ERM-proteins in the blastocyst and the uterus has been linked to the implantation potential in mice. Protein analysis of uterine fluid has demonstrated the presence of ezrin (EZR) and moesin (MSN) in cattle26,27 and pigs23. In the horse, upregulation of EZR and downregulation of MSN was detected, while an inverse association with pregnancy was noted in cattle26,27. We detected both EZR and MSN in the YS and they were also found in the blastocoel fluid34.
Apart from the specific interest in proteins involved in prostaglandin synthesis, we also aimed to create a general overview of the proteins present at the embryo-maternal interface and potentially involved in signalling. Supplementary File 4 presents all proteins which were commonly found in the UF of P mares and in the YS and the functions of each protein are included. To visualize their role in embryo-maternal interaction and signalling, the proteins involved in GO terms including ‘embryo’, ‘maternal’ or ‘uterus’ are depicted in Fig. 4 and those linked to GO terms ‘growth factor’ and ‘cytokine’ in Fig. 5. The origin of the proteins can be distinguished in red (UF), yellow (YS) and pink (UF and YS) and up- and downregulation during pregnancy is represented by enlargement or shrinkage of the protein respectively. Interestingly, most proteins which were found to be upregulated during pregnancy were detected both in UF and in YS, while downregulation during pregnancy generally coincided with absence of these proteins in the YS, indicating a potentially important role of the embryo in the production of these proteins during pregnancy. Several common proteins were found at the embryo-maternal interface during MRP in cattle, including aconitase 1 (ACO1), which was specifically detected in the uterine fluid of pregnant and not in cyclic heifers, as well as glucose-6-phospate isomerase (GPI), which has been detected in the uterine fluid of both pregnant and cyclic heifers and for which an embryonic source has been presumed based on transcriptomics27. In our study, both proteins were found in the YS and the UF of P mares, with significant upregulation of GPI in P versus C. Two other proteins which were commonly found in UF and YS, namely FK506 binding protein 4 (FKBP4) (Fig. 4) and heat shock protein 90 (HSP90AB1) (Fig. 5), have been elaborately discussed by Swegen et al.34 concerning their progesterone supportive role. Co-operation of both factors is necessary for activation of the progesterone receptor48, FKBP4 has shown to be crucial for uterine receptivity and implantation in mice49 and FKBP4 deficit has been associated with pregnancy loss in human50. While FKBP4 was detected in equine blastocoel fluid and speculated to be involved in signalling, it was not detected in the embryo-conditioned medium34. Interestingly, we did find both FKBP4 and HSP90AB1, not only in YS, but also in UF, even though their presence was not affected by pregnancy. While further confirmation of the role of specific proteins is required, the overview created in this study can be used as a basis for further targeted studies in the horse.
In addition to the role in prostaglandin and progesterone metabolism, involvement in proteolysis and lipid metabolism was also prominent in the commonly detected proteins in our study and the one of Swegen et al.34, also coinciding with previous findings on transcriptomics around MRP15. Several cathepsins (G, D, L and S) were detected in the UF with downregulation of cathepsin L (CTSL) and S during pregnancy. Pregnancy associated downregulation of CTSL1 was also found at the transcriptome level15. Considering lipid metabolism, we detected differential expression of lipocalin (P19), apolipoprotein A1 (APOA1) and apolipoprotein D (APOD). These proteins are important transporters of essential lipids to the developing conceptus. Retinol binding protein (RBP), which also belongs to the lipocalin family, and APOA1 have been detected in uterine fluid of pregnant and cyclic pigs, cattle and sheep23,24,25,27,28, with increasing amounts between day 10 and day 13 in both pregnant and cyclic pigs24. In the horse, lipocalin (P19), apolipoprotein A1 (APOA1) and apolipoprotein D (APOD) were all downregulated in the uterine fluid of the pregnant mares. Pregnancy associated upregulation of APOA1 was reported at the transcriptome level15 and presence of P19 and APOA1 in the yolk sac fluid illustrates their role in the embryo-maternal dialogue. Therefore, lower amounts in the uterine fluid during pregnancy rather indicate the transport and binding to the conceptus. Lipocalin P19 or uterocalin is a progesterone induced protein, which is abundantly present in the equine uterine secretions during dioestrus and early pregnancy51,52. While the early developing equine conceptus moves around the uterus, it entirely depends upon the uterine secretions for its nutrition and P19 can function as a carrier for essential lipids and amino acids53. Coinciding with our findings, P19 has been detected in the trophoblast and the yolk sac fluid of the equine embryo51,52 and it is one of the most abundant proteins in the embryonic capsule54,55,56. Therefore, the lower amount of P19 in P versus C is probably due to binding of substantial quantities to the embryo.
While a novel and informative overview is created, it has to be borne in mind that no statistically significant results were obtained at the level of molecular functions, biological processes and pathways. Differential expression of individual proteins was observed between the different UF conditions, and these proteins were categorized in GO terms, but statistical analysis showed no significant overrepresentation of any of the GO terms or KEGG pathways. Furthermore, it should be noted that MS intrinsically suffers from missing values and conclusions based on the absence of proteins cannot be made42,43. However, the field of proteomics has greatly evolved in recent years, providing the possibility for statistically robust quantitative comparison of individual protein levels in complex biological samples, like uterine fluid22. HDMSE specifically has been shown to provide good proteome coverage and reproducibility35. At the same time, however, analysis of GO terms and pathways for proteomics is still in its infancy57,58. As many of the here described bioinformatics approaches for proteomic analysis were originally developed for genomics, a similar but more matured field, their performance can be expected to show a similar growth as that of the genomic approaches. Moreover, the similarity between these fields potentially allows an integrated approach in which results from several omics studies can be combined.
In conclusion, proteins present in the equine uterine and embryonic yolk sac fluid around the signalling of MRP at day 13 were identified and quantified at large scale for the first time in the horse. We detected upregulation of several inhibitors of prostaglandin synthesis, including PTGR1, GSTP1 and ANXA1, in the uterine fluid of pregnant mares. Overall, an overview was created of the proteins playing a role at the embryo-maternal interface in the horse. This study provides a solid foundation for further targeted studies of proteins potentially involved in embryo-maternal interactions, maternal recognition of pregnancy and pregnancy loss in the horse.
Methods
Sampling
All animal handlings were approved by the Ethical Committee of the Faculty of Veterinary Medicine (EC2013/118) of Ghent University. All methods were performed in accordance with the relevant guidelines and regulations. A switch back design was followed with 5 mares undergoing two different types of cycles: a pregnant cycle (P) and a cyclic control cycle (C). In this way, the samples were paired using the same mare as its own control for pregnancy and the experimental unit was the mare. The order of P and C cycles was randomly altered for the different mares. No resting cycles were included. During the breeding season, five reproductively sound Warmblood mares between 4 and 13 years old were monitored by transrectal ultrasound. Reproductive soundness was confirmed by negative cytology and bacteriology. Mares displaying uterine oedema and a follicle exceeding 35 mm received 1500 IU hCG intravenously and were either inseminated the next day with fresh semen of the same stallion (P) or left unbred (C). Ovulation was evaluated twice daily by ultrasound. In both groups, sampling was performed 13 days after detection of ovulation. To recover undiluted uterine fluid in order to avoid negative effects of excessive Ringer’s salts on MS59, intra-uterine application of a tampon (OB Mini; Johnson & Johnson, Beerse, Belgium) was performed based upon the method described by Wolf et al.33. A double gloved technique was used to avoid vaginal contamination. The tampon was left in the uterus during 10 minutes and upon removal it was placed in a Falcon tube at 4 °C until further processing. Subsequently, the mare’s uterus was flushed with sterile Ringer’s solution by means of a modified endotracheal tube to recover the embryo (P).
To process the uterine fluid, 1 mL of sterile water (B60, Biosolve, Valkenswaard, The Netherlands) was infused on top of the tampon and the tampon was attached in the upper part of the Falcon tube by fixing the cord with the cap. Subsequently, the Falcon tube was centrifuged for 20 minutes at 1000 × g at 4 °C. The supernatant was collected and stored in a Protein LoBind Eppendorf tube (Eppendorf AG, Hamburg, Germany) at −80 °C. Meanwhile, the embryo was isolated in a petri dish and the yolk sac fluid was collected by aspiration with a 21 G needle and stored in a Protein LoBind Eppendorf at −80 °C.
A total of 15 samples were collected, consisting of uterine fluid (UF) (n = 10) from five biological replicates coinciding with the five mares (1–5) for the P and C treatment cycles, as well as yolk sac fluid (YS) (n = 5) from the P cycles.
Sample preparation for mass spectrometry analysis
After thawing, protein concentration in each sample was determined using the Coomassie (Bradford) Protein Assay Kit (Thermo Fisher Scientific, San José, CA, USA) according to the manufacturer’s instructions. Further processing was performed for 10 µg protein of each uterine fluid sample and for 500 ng protein of the yolk sac samples. Samples were dissolved in 20 μL 0.5 M triethylammonium bicarbonate (TEABC; Sigma-Aldrich, St. Louis, MO, USA). Two µl of reducing agent (10 µM DTT; Invitrogen, Merelbeke, Belgium) were added followed by incubation for 1 h at 60 °C. Subsequently, 1 µl of alkylizing agent (200 mM methyl methanethiosulfonate (MMTS) in isopropanol; Sigma-Aldrich, St. Louis, MO, USA) was added and samples were incubated for 10 min at room temperature. Digestion was performed overnight at 37 °C with trypsin lys C (1:20, trypsin:protein w/w, Promega, Leiden, The Netherlands) in TEABC buffer with 1 mM CACL(2) and 5% acetonitrile (Biosolve, Valkenswaard, The Netherlands). Samples were vacuum-dried and stored at −20 °C until analysis.
Data acquisition by HDMSE analysis
The peptides were separated using a nanoscale UPLC system (nanoAcquityUPLC, Waters, Milford, USA) coupled to a Synapt G2-Si mass spectrometer (Waters). Peptides were first trapped in 0.1% formic acid on a 180 µm × 20 mm C18 Trap column. Separation was performed on a HSS C18 1.8 m, 100 m × 250 mm analytical column at a flow rate of 300 nL/min and a temperature of 45 °C. As mobile phase A a 0.1% formic acid with 4% DMSO in water solution was used and 80% ACN containing 0.1% formic acid constituted mobile phase B. Peptides were separated for 60 min at 1–40% solvent B and for 1 min 40–85% solvent B. Seven minutes of rinsing (85% solvent B) re-equilibrated the column to the initial conditions. Eluted peptides were analysed in positive mode ESI-MS using High Definition MSE (HDMSE) with a collision energy look up table as described in22. The spectral acquisition time of low and elevated energy scans was 0.6 s over an m/z range of 50–2000. [Glu1]-Fibrinopeptide B was used for post-acquisition lock mass correction. All UF samples were analysed in the same run; three technical replicates (R1–R3) were run for each sample and four quality controls (QC) were included in which all samples were pooled.
Identification and quantification of peptides and proteins
All data were processed in Progenesis QIP (Progenesis QIP 2.0, Nonlinear Dynamics, Waters), including normalization and quality control. A database with UniProt IDs was created by conversion of Ensembl gene identifiers for Equus caballus (n = 22295) to Uniprot IDs using http://www.uniprot.org/uploadlists/ and including common contaminants (http://www.thegpm.org/crap/). As only secreted proteins are expected to be found, it can be argued that this database should be limited to only these secreted proteins. However, there is much debate on the accuracy of FDR calculations with such limited databases60,61,62 and as such a cautious approach was taken in which all proteins were assessed. Using Progenesis QIP, peptides were identified against this database with a FDR of 4%63 and allowing maximum one miscleavage. Protein quantification was based on the Hi-3 method64, which uses the average of the three most intense peptides of each protein for its quantification. Resulting normalized abundances for each protein, as well as unique peptide counts were further used for analysis of differential expression.
Analysis of differential expression
Analysis of differential expression was performed for the UF samples. Only normalized abundancies of proteins with at least two unique peptides (n = 707) were included in the analysis. Pairwise comparisons of differential expression were made for P versus C with the individual horses as a blocking factor, using R Bioconductor limma package65 and a FDR of 0.05.
Gene Ontology enrichment and pathway analysis
Gene Ontology (GO) terms (molecular functions, biological processes and cellular locations) were downloaded for each protein with the PANTHER Classification System66. For the pair-wise comparison of P and C, a statistical overrepresentation test against all quantified proteins was performed for all significantly up- and downregulated (FDR < 0.05) proteins. These tests were done for all primary GO classes: molecular functions, biological processes and cellular components. Pathways were analysed with Bioconductor’s67 GAGE package68. LogFC values of all quantified proteins were used as input against Equus caballus background reference pathways from KEGG.
Proteins involved in the embryo-maternal interaction were visualized using Cytoscape 3.3.0. To visualize embryo-maternal signalling, all GO terms and pathways that include ‘cytokine’ or ‘growth factor’ in their description were selected, together with all identified proteins in either the yolk sac or uterine fluid belonging to these GO terms or pathways. The connection between these GO terms, pathways and proteins was then visualized using Cytoscape 3.3.0. The same methodology was used to create a network based on GO terms and pathways including ‘embryo’, ‘maternal’ or ‘uterus’ in their description.
Data availability
All data are available in the Supplementary files.
References
Short, R. Implantation and the maternal recognition of pregnancy. Foetal Autonomy 2, 31 (1969).
Kindahl, H., Knudsen, O., Madej, A. & Edqvist, L. E. Progesterone, prostaglandin F-2 alpha, PMSG and oestrone sulphate during early pregnancy in the mare. J Reprod Fertil Suppl 32, 353–359 (1982).
Bazer, F. W. & Thatcher, W. Theory of maternal recognition of pregnancy in swine based on estrogen controlled endocrine versus exocrine secretion of prostaglandin F 2α by the uterine endometrium. Prostaglandins 14, 397–401 (1977).
Lamming, G. E. et al. Local action of trophoblast interferons in suppression of the development of oxytocin and oestradiol receptors in ovine endometrium. J Reprod Fertil 105, 165–175 (1995).
Lamming, G. E. & Mann, G. E. Control of endometrial oxytocin receptors and prostaglandin F2 alpha production in cows by progesterone and oestradiol. J Reprod Fertil 103, 69–73 (1995).
Allen, W. R. Fetomaternal interactions and influences during equine pregnancy. Reproduction 121, 513–527 (2001).
Klein, C. & Troedsson, M. H. Maternal recognition of pregnancy in the horse: a mystery still to be solved. Reprod Fertil Dev 23, 952–963, https://doi.org/10.1071/RD10294 (2011).
Vanderwall, D. K., Silvia, W. J. & Fitzgerald, B. P. Concentrations of oxytocin in the intercavernous sinus of mares during luteolysis: temporal relationship with concentrations of 13, 14-dihydro-15-keto-prostaglandin F2 alpha. J Reprod Fertil 112, 337–346 (1998).
Goff, A. K., Pontbriand, D. & Sirois, J. Oxytocin stimulation of plasma 15-keto-13, 14-dihydro prostaglandin F-2 alpha during the oestrous cycle and early pregnancy in the mare. J Reprod Fertil Suppl 35, 253–260 (1987).
Starbuck, G. R., Stout, T. A., Lamming, G. E., Allen, W. R. & Flint, A. P. Endometrial oxytocin receptor and uterine prostaglandin secretion in mares during the oestrous cycle and early pregnancy. J Reprod Fertil 113, 173–179 (1998).
Boerboom, D. et al. Expression of key prostaglandin synthases in equine endometrium during late diestrus and early pregnancy. Biol Reprod 70, 391–399, https://doi.org/10.1095/biolreprod.103.020800 (2004).
Ealy, A. D., Eroh, M. L. & Sharp, D. C. III. Prostaglandin H synthase Type 2 is differentially expressed in endometrium based on pregnancy status in pony mares and responds to oxytocin and conceptus secretions in explant culture. Anim Reprod Sci 117, 99–105, https://doi.org/10.1016/j.anireprosci.2009.03.014 (2010).
de Ruijter-Villani, M., van Tol, H. T. & Stout, T. A. Effect of pregnancy on endometrial expression of luteolytic pathway components in the mare. Reprod Fertil Dev 27, 834–845, https://doi.org/10.1071/RD13381 (2015).
Wilsher, S., Clutton-Brock, A. & Allen, W. R. Successful transfer of day 10 horse embryos: influence of donor-recipient asynchrony on embryo development. Reproduction 139, 575–585, https://doi.org/10.1530/REP-09-0306 (2010).
Klein, C., Scoggin, K. E., Ealy, A. D. & Troedsson, M. H. Transcriptional profiling of equine endometrium during the time of maternal recognition of pregnancy. Biol Reprod 83, 102–113, https://doi.org/10.1095/biolreprod.109.081612 (2010).
Klein, C. & Troedsson, M. H. Transcriptional profiling of equine conceptuses reveals new aspects of embryo-maternal communication in the horse. Biol Reprod 84, 872–885, https://doi.org/10.1095/biolreprod.110.088732 (2011).
Merkl, M. et al. Microarray analysis of equine endometrium at days 8 and 12 of pregnancy. Biol Reprod 83, 874–886, https://doi.org/10.1095/biolreprod.110.085233 (2010).
Bauersachs, S. & Wolf, E. Transcriptome analyses of bovine, porcine and equine endometrium during the pre-implantation phase. Anim Reprod Sci 134, 84–94, https://doi.org/10.1016/j.anireprosci.2012.08.015 (2012).
Klein, C. Novel equine conceptus?endometrial interactions on Day 16 of pregnancy based on RNA sequencing. Reprod Fertil Dev. https://doi.org/10.1071/RD14489 (2015).
Wright, P. C., Noirel, J., Ow, S. Y. & Fazeli, A. A review of current proteomics technologies with a survey on their widespread use in reproductive biology investigations. Theriogenology 77, 738–765 e752, https://doi.org/10.1016/j.theriogenology.2011.11.012 (2012).
Vogel, C. & Marcotte, E. M. Insights into the regulation of protein abundance from proteomic and transcriptomic analyses. Nature reviews. Genetics 13, 227–232, https://doi.org/10.1038/nrg3185 (2012).
Distler, U., Kuharev, J. & Tenzer, S. Biomedical applications of ion mobility-enhanced data-independent acquisition-based label-free quantitative proteomics. Expert review of proteomics 11, 675–684, https://doi.org/10.1586/14789450.2014.971114 (2014).
Kayser, J. P., Kim, J. G., Cerny, R. L. & Vallet, J. L. Global characterization of porcine intrauterine proteins during early pregnancy. Reproduction 131, 379–388, https://doi.org/10.1530/rep.1.00882 (2006).
Jalali, B. M., Bogacki, M., Dietrich, M., Likszo, P. & Wasielak, M. Proteomic analysis of porcine endometrial tissue during peri-implantation period reveals altered protein abundance. Journal of proteomics 125, 76–88, https://doi.org/10.1016/j.jprot.2015.05.003 (2015).
Brooks, K., Burns, G. W., Moraes, J. G. & Spencer, T. E. Analysis of the Uterine Epithelial and Conceptus Transcriptome and Luminal Fluid Proteome During the Peri-Implantation Period of Pregnancy in Sheep. Biol Reprod 95(88), 1–17, https://doi.org/10.1095/biolreprod.116.141945 (2016).
Munoz, M. et al. Proteome of the early embryo-maternal dialogue in the cattle uterus. J Proteome Res 11, 751–766, https://doi.org/10.1021/pr200969a (2012).
Forde, N., Bazer, F. W., Spencer, T. E. & Lonergan, P. ‘Conceptualizing’ the Endometrium: Identification of Conceptus-Derived Proteins During Early Pregnancy in Cattle. Biol Reprod 92(156), 1–13, https://doi.org/10.1095/biolreprod.115.129296 (2015).
Forde, N. et al. Proteomic analysis of uterine fluid during the pre-implantation period of pregnancy in cattle. Reproduction, https://doi.org/10.1530/REP-13-0010(2014).
Tachibana, Y. et al. Expression of endometrial immune-related genes possibly functioning during early pregnancy in the mare. J Reprod Dev 59, 85–91 (2013).
Klein, C., Scoggin, K. E. & Troedsson, M. H. The expression of interferon-stimulated gene 15 in equine endometrium. Reprod Domest Anim 46, 692–698, https://doi.org/10.1111/j.1439-0531.2010.01731.x (2011).
Hartt, L. S. et al. Temporal and spatial associations of oestrogen receptor alpha and progesterone receptor in the endometrium of cyclic and early pregnant mares. Reproduction 130, 241–250, https://doi.org/10.1530/rep.1.00596 (2005).
Watson, E. D., Buckingham, J., Bjorksten, T. & Nikolakopoulos, E. Immunolocalization of oxytocin and neurophysin in the mare uterus. J Reprod Fertil Suppl 289–296 (2000).
Wolf, C. A., Maslchitzky, E., Gregory, R. M., Jobim, M. I. & Mattos, R. C. Effect of corticotherapy on proteomics of endometrial fluid from mares susceptible to persistent postbreeding endometritis. Theriogenology 77, 1351–1359, https://doi.org/10.1016/j.theriogenology.2011.10.042 (2012).
Swegen, A. et al. From Peptide Masses to Pregnancy Maintenance: A Comprehensive Proteomic Analysis of The Early Equine Embryo Secretome, Blastocoel Fluid, and Capsule. Proteomics 17, https://doi.org/10.1002/pmic.201600433 (2017).
Distler, U. et al. Drift time-specific collision energies enable deep-coverage data-independent acquisition proteomics. Nature methods 11, 167–170, https://doi.org/10.1038/nmeth.2767 (2014).
Hayes, J. D., Flanagan, J. U. & Jowsey, I. R. Glutathione transferases. Annual review of pharmacology and toxicology 45, 51–88, https://doi.org/10.1146/annurev.pharmtox.45.120403.095857 (2005).
Luo, L. et al. Recombinant protein glutathione S-transferases P1 attenuates inflammation in mice. Molecular immunology 46, 848–857, https://doi.org/10.1016/j.molimm.2008.09.010 (2009).
Wallner, B. P. et al. Cloning and expression of human lipocortin, a phospholipase A2 inhibitor with potential anti-inflammatory activity. Nature 320, 77–81, https://doi.org/10.1038/320077a0 (1986).
Huang, K. S. et al. Two human 35 kd inhibitors of phospholipase A2 are related to substrates of pp60v-src and of the epidermal growth factor receptor/kinase. Cell 46, 191–199 (1986).
Smits, K. et al. The Equine Embryo Influences Immune-Related Gene Expression in the Oviduct. Biol Reprod 94, 36, https://doi.org/10.1095/biolreprod.115.136432 (2016).
Ababneh, M. M. & Troedsson, M. H. Endometrial phospholipase A2 activity during the oestrous cycle and early pregnancy in mares. Reprod Domest Anim 48, 46–52, https://doi.org/10.1111/j.1439-0531.2012.02023.x (2013).
Albrecht, D., Kniemeyer, O., Brakhage, A. A. & Guthke, R. Missing values in gel-based proteomics. Proteomics 10, 1202–1211, https://doi.org/10.1002/pmic.200800576 (2010).
Webb-Robertson, B. J. et al. Review, evaluation, and discussion of the challenges of missing value imputation for mass spectrometry-based label-free global proteomics. J Proteome Res 14, 1993–2001, https://doi.org/10.1021/pr501138h (2015).
Park, E. S. et al. Phospholipase C-delta1 and oxytocin receptor signalling: evidence of its role as an effector. Biochem J 331(Pt 1), 283–289 (1998).
Bu, H. et al. ERBP, a novel estrogen receptor binding protein enhancing the activity of estrogen receptor. Biochem Biophys Res Commun 317, 54–59, https://doi.org/10.1016/j.bbrc.2004.02.179 (2004).
Montt-Guevara, M. M. et al. Androgens Regulate T47D Cells Motility and Invasion through Actin Cytoskeleton Remodeling. Frontiers in endocrinology 7, 136, https://doi.org/10.3389/fendo.2016.00136 (2016).
Song, J. et al. Expression and clinicopathological significance of oestrogen-responsive ezrin-radixin-moesin-binding phosphoprotein 50 in breast cancer. Histopathology 51, 40–53, https://doi.org/10.1111/j.1365-2559.2007.02730.x (2007).
Tranguch, S., Smith, D. F. & Dey, S. K. Progesterone receptor requires a co-chaperone for signalling in uterine biology and implantation. Reprod Biomed Online 13, 651–660 (2006).
Tranguch, S. et al. Cochaperone immunophilin FKBP52 is critical to uterine receptivity for embryo implantation. Proceedings of the National Academy of Sciences of the United States of America 102, 14326–14331, https://doi.org/10.1073/pnas.0505775102 (2005).
Chen, H. Y. et al. Expression of FK506-binding protein 52 (FKBP52) in chorionic villi with early recurrent spontaneous abortion. J Matern Fetal Neonatal Med 28, 1165–1169, https://doi.org/10.3109/14767058.2014.947572 (2015).
Crossett, B., Allen, W. R. & Stewart, F. A 19 kDa protein secreted by the endometrium of the mare is a novel member of the lipocalin family. Biochem J 320(Pt 1), 137–143 (1996).
Ellenberger, C. et al. Immunolocalisation of the uterine secretory proteins uterocalin, uteroferrin and uteroglobin in the mare’s uterus and placenta throughout pregnancy. Theriogenology 70, 746–757, https://doi.org/10.1016/j.theriogenology.2008.04.050 S0093-691X(08)00290-2 [pii] (2008).
Suire, S., Stewart, F., Beauchamp, J. & Kennedy, M. W. Uterocalin, a lipocalin provisioning the preattachment equine conceptus: fatty acid and retinol binding properties, and structural characterization. Biochem J 356, 369–376 (2001).
Quinn, B. A., Hayes, M. A., Waelchli, R. O., Kennedy, M. W. & Betteridge, K. J. Changes in major proteins in the embryonic capsule during immobilization (fixation) of the conceptus in the third week of pregnancy in the mare. Reproduction 134, 161–170, https://doi.org/10.1530/REP-06-0241 (2007).
Stewart, F., Charleston, B., Crossett, B., Barker, P. J. & Allen, W. R. A novel uterine protein that associates with the embryonic capsule in equids. J Reprod Fertil 105, 65–70 (1995).
Smits, K. et al. Influence of the uterine environment on the development of in vitro-produced equine embryos. Reproduction 143, 173–181, https://doi.org/10.1530/REP-11-0217 REP-11-0217 [pii] (2012).
Schmidt, A., Forne, I. & Imhof, A. Bioinformatic analysis of proteomics data. BMC systems biology 8(Suppl 2), S3, https://doi.org/10.1186/1752-0509-8-S2-S3 (2014).
Vaudel, M., Sickmann, A. & Martens, L. Introduction to opportunities and pitfalls in functional mass spectrometry based proteomics. Biochim Biophys Acta 1844, 12–20, https://doi.org/10.1016/j.bbapap.2013.06.019 (2014).
Piwowar, A. M., Lockyer, N. P. & Vickerman, J. C. Salt effects on ion formation in desorption mass spectrometry: an investigation into the role of alkali chlorides on peak suppression in time-of-flight-secondary ion mass spectrometry. Anal Chem 81, 1040–1048, https://doi.org/10.1021/ac8020888 (2009).
Noble, W. S. Mass spectrometrists should search only for peptides they care about. Nature methods 12, 605–608, https://doi.org/10.1038/nmeth.3450 (2015).
Noble, W. S. & Keich, U. Mass spectrometrists should search for all peptides, but assess only the ones they care about Reply. Nature methods 14, 644–644 (2017).
Sticker, A., Martens, L. & Clement, L. Mass spectrometrists should search for all peptides, but assess only the ones they care about. Nature methods 14, 643–644 (2017).
Li, G. Z. et al. Database searching and accounting of multiplexed precursor and product ion spectra from the data independent analysis of simple and complex peptide mixtures. Proteomics 9, 1696–1719, https://doi.org/10.1002/pmic.200800564 (2009).
Silva, J. C., Gorenstein, M. V., Li, G. Z., Vissers, J. P. & Geromanos, S. J. Absolute quantification of proteins by LCMSE: a virtue of parallel MS acquisition. Mol Cell Proteomics 5, 144–156, https://doi.org/10.1074/mcp.M500230-MCP200 (2006).
Ritchie, M. E. et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic acids research 43, e47, https://doi.org/10.1093/nar/gkv007 (2015).
Mi, H., Muruganujan, A., Casagrande, J. T. & Thomas, P. D. Large-scale gene function analysis with the PANTHER classification system. Nature protocols 8, 1551–1566, https://doi.org/10.1038/nprot.2013.092 (2013).
Huber, W. et al. Orchestrating high-throughput genomic analysis with Bioconductor. Nature methods 12, 115–121, https://doi.org/10.1038/nmeth.3252 (2015).
Luo, W., Friedman, M. S., Shedden, K., Hankenson, K. D. & Woolf, P. J. GAGE: generally applicable gene set enrichment for pathway analysis. BMC bioinformatics 10, 161, https://doi.org/10.1186/1471-2105-10-161 (2009).
Acknowledgements
The authors would like to thank Petra Van Damme and Isabel Lemahieu for their excellent technical assistance. Katrien Smits works as a postdoctoral researcher for the Research Foundation Flanders (grant FWO13/PDO/08). This research was funded by the Research Foundation Flanders (FWO project GO35511N) and by Ghent University (GOA project 01G01112).
Author information
Authors and Affiliations
Contributions
K.S. performed the experiments and wrote the manuscript. S.W. provided the bioinformatics analysis. K.V.S. was responsible for the mass spectrometry. M.V.D.V., V.D.L., C.V., K.R. and J.G. contributed to the sampling of the mares. F.V.N., L.P., D.D. and A.V.S. participated in the study design. All authors reviewed the manuscript.
Corresponding author
Ethics declarations
Competing Interests
The authors declare no competing interests.
Additional information
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
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 license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Smits, K., Willems, S., Van Steendam, K. et al. Proteins involved in embryo-maternal interaction around the signalling of maternal recognition of pregnancy in the horse. Sci Rep 8, 5249 (2018). https://doi.org/10.1038/s41598-018-23537-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-018-23537-6
This article is cited by
-
Uterine extracellular vesicles as multi-signal messengers during maternal recognition of pregnancy in the mare
Scientific Reports (2022)
-
Spatiotemporal endometrial transcriptome analysis revealed the luminal epithelium as key player during initial maternal recognition of pregnancy in the mare
Scientific Reports (2021)
-
The influence of progesterone on bovine uterine fluid energy, nucleotide, vitamin, cofactor, peptide, and xenobiotic composition during the conceptus elongation-initiation window
Scientific Reports (2019)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.