Spatiotemporal Regulation of The Bovine Oviduct Fluid Proteome Around The Time of Ovulation

doi:10.21203/rs.3.rs-971243/v1

Download PDF

Research Article

Spatiotemporal Regulation of The Bovine Oviduct Fluid Proteome Around The Time of Ovulation

https://doi.org/10.21203/rs.3.rs-971243/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Unraveling the oviduct fluid (OF) and its regulation is crucial to understand the microenvironment in which sperm capacitation, fertilization and early embryo development take place. Therefore, our aim was to determine the spatiotemporal changes in the OF proteome according to the anatomical region of the oviduct (ampulla vs. isthmus), the proximity of the ovulating ovary (ipsilateral vs. contralateral side) and the peri-ovulatory stage (pre-ovulatory or Pre-ov vs. post-ovulatory or Post-ov). Oviducts from adult cyclic cows were collected at a local slaughterhouse and pools of OF were analyzed by nanoLC-MS/MS and label-free protein quantification (n=32 OF pools for all region ´ stage ´ side conditions). A total of 3,760 proteins were identified in the OF, of which 37% were predicted to be secreted. The oviduct region was the major source of variation in protein abundance, followed by the proximity of the ovulating ovary and finally the peri-ovulatory stage. Differentially abundant proteins between regions, stages and sides were involved in a broad variety of biological functions, including protein binding, response to stress, cell-to-cell adhesion, calcium homeostasis and the immune system. This work highlights the dynamic regulation of oviduct secretions and provides new protein candidates for interactions between the maternal environment, the gametes and the early embryo.

Molecular Biology

oviduct

spatiotemporal

embryo

proteins

In mammals, pregnancy establishment requires a sequence of events in the different segments of the oviduct around the time of ovulation. After migration through the uterus, sperm reach the isthmus (i.e., the caudal part of the oviduct) and bind to oviduct epithelial cells (OECs) to form a sperm reservoir until ovulation time when they are gradually released toward the ampulla^1,2. From the isthmus to the ampulla, spermatozoa are progressively capacitated and acquire the ability to fertilize oocytes³. In the immediate post-ovulatory period, the oocyte and surrounding cumulus cells go through the ampulla where fertilization takes place⁴. In the 3-5 days following ovulation, cell divisions of the zygote start in the isthmus and the early embryo finally enters the uterus⁵.

The oviduct fluid (OF) offers an optimal environment for these early reproductive events. The OF is a complex mixture of small metabolites, lipids, glycans and proteins originating from secretions of the luminal epithelium, selective ultrafiltrate from the circulating plasma and putative inputs from the follicular fluid⁶. Previous work from our group showed that the steroidomic⁷, metabolomic⁸ and proteomic⁹ composition of the bovine OF varied throughout the estrous cycle and according to the proximity of the Pre-ov follicle (POF) or corpus luteum (CL). Furthermore, variations in gene expression have been reported between ampulla and isthmus^10–13. However, whether these variations in gene expression are translated into variations in OF protein abundance between ampulla and isthmus is currently not known. Recently, a small set of OF proteins were reported to vary in abundance between cyclic and pregnant cows on day 3 post-estrus, with slight differences between the ampulla and isthmus¹⁴. Region-related differences in oviductal secretions were reported in vitro: effects of OF-derived extracellular vesicles (EVs)¹⁵ and OF proteins¹⁶ on sperm viability and capacitation differed between ampulla and isthmus. Moreover, the quality of embryos developed in vitro with bovine OF-derived EVs¹⁷ or ovine OECs¹⁸ differed depending on the ampullary or isthmic origin of the EVs or OECs.

Which factors contribute the most to the regulation of the OF proteome remains unclear. We hypothesize that the OF proteome is monitored by a complex set of factors: the anatomical region of the oviduct, the proximity of POF/CL and the stage of cycle. A dynamic variation between the ampulla and isthmus in the oviduct ipsilateral to ovulation, but also between the Pre-ov and Post-ov stages are expected, in response to the specific requirements of gametes and embryos. Therefore, the objective of this study was to determine the effect on the OF proteome of (1) the region (isthmus vs. ampulla), (2) the proximity of the ovulating ovary (ipsilateral vs. contralateral side) and (3) the peri-ovulatory stage (Pre-ov vs. Post-ov) on the OF proteome using nanoLC-MS/MS analysis.

Collection of bovine oviductal fluid (OF)

Both oviducts from adult cyclic cows were collected at a local slaughterhouse (40 min from the lab), placed immediately on ice and transported to the laboratory within max 2 h after animal death. The oviducts were classified into pre-ovulatory (Pre-ov; presence of a corpus albicans and one pre-ovulatory follicle of 11-20 mm in diameter, corresponding to approximately days 19-21 of cycle) or post-ovulatory (Post-ov; presence of an early corpus luteum made with red loosely organized tissue less than 1 cm in diameter, corresponding to approximately day 1-4 of cycle), as previously described⁶⁰. Oviducts ipsilateral and contralateral to ovulation side were trimmed free of surrounding tissue and vessels and linearized. After removal of the infundibulum and utero-tubal junction, both oviducts were divided into two parts at the ampullary-isthmic junction. The ampulla and isthmus were individually flushed with 200 µL of cold sterile protein-free phosphate buffered saline (PBS) and kept on ice. Cellular debris and cells were eliminated from the OF by two centrifugations (2,000 g for 10 min then 12,000 g for 10 min) at 4°C. Samples were stored at −80°C until proteomic analysis. A total of 18 Pre-ov and 16 Post-ov cows were collected. In the Pre-ov group, in order to exclude animals with ovarian cyst or atretic follicle, the follicular fluid of the Pre-ov follicle was collected and stored at – 80°C for steroid hormone analysis, as previously described⁶¹. Concentrations of P4 were measured by a competitive enzyme-linked immunosorbent assay according to a method previously described⁶² whereas E2 concentrations were measured using the BioSource E2-EASIA kit (BioSource Europe S.A., Louvain-la-Neuve, Belgium). Assays were performed in duplicate. Two cows were excluded from the Pre-ov group based on intra-follicular P4 concentrations higher than 160 ng/mL, as previously described⁷. In the resulting Pre-ov group (n = 16), the ranges of P4 and E2 intra-follicular concentrations were 19.4-134.2 ng/mL and 103->2000 ng/mL, respectively.

Preparation Of Of Samples For Proteomic Analyses

The OF samples were thawed on ice and pools of 100 µL of OF were made by mixing 25 µL of individual samples. Four pools of 4 cows per side ⋅ region ⋅ stage condition were constituted, leading to 32 OF samples distributed in 8 conditions (Pre-ov ipsilateral ampulla; Pre-ov ipsilateral isthmus; Pre-ov contralateral ampulla; Pre-ov contralateral isthmus; Post-ov ipsilateral ampulla; Post-ov ipsilateral isthmus; Post-ov contralateral ampulla and Post-ov contralateral isthmus). In the following, the term “sample” refers to these pools of OF. Protein concentration of OF samples was assessed using the Uptima BC Assay kit (Interchim, Montluçon, France) according to manufacturer’s instructions and using bovine serum albumin as a standard. The homogeneity of samples was checked after migration of the OF pools (5 µg proteins per lane) on a 4-20% SDS-PAGE (80 V, 10 min and 180 V 30 min) and staining with Coomassie blue (see Supplementary Figure S3).

Nanoliquid Chromatography Coupled With Tandem Mass Spectrometry (Nanolc-ms/ms)

For proteomic analysis, 10 µg of proteins per sample were prepared using the PreOmics iST-BCT kit (PreOmics GmbH, Martinsried, Germany) following the manufacturer’s instructions. Briefly, proteins were lysed, denatured, reduced and alkylated for 10 min at 95° then Trypsin/LysC digested for 60 min at 37°C. Purification and elution of peptides were then carried out for 1 h at room temperature. Samples were then purified from salts, contaminants and other detergents, using reversed phase SDBS spin cartridge (PreOmics GmbH, Martinsried, Germany)). Peptide injection and analysis were performed as previously described⁶³. The resulting peptide mixtures were loaded on a 75 µm × 250 mm IonOpticks Aurora 2 C18 column (Ion Opticks Pty Ltd., Bundoora, Australia). A gradient of basic reversed-phase buffers (Buffer A: 0.1% formic acid, 98% H2O MilliQ, 2% acetonitrile; Buffer B: 0.1% formic acid, 100% acetonitrile) was run on a NanoElute HPLC System (Bruker Daltonik GmbH, Bremen, Germany) at a flow rate of 400 nL/min at 50°C. The liquid chromatography (LC) run lasted for 120 min (2–15% of buffer B during 60 min; up to 25% at 90 min; up to 37% at 100 min; up to 95% at 110 min and finally 95% for 10 min to wash the column). The column was coupled online to a TIMS TOF Pro (Bruker Daltonik GmbH, Bremen, Germany) with a CaptiveSpray ion source (Bruker Daltonik). The temperature of the ion transfer capillary was set at 180°C. Ions were accumulated for 114 ms, and mobility separation was achieved by ramping the entrance potential from −160 V to −20 V within 114 ms. The acquisition of the MS and MS/MS mass spectra was done with average resolutions of 60,000 and 50,000 full width at half maximum (mass range 100–1700 m/z), respectively. To enable the PASEF method, precursor m/z and mobility information was first derived from full scan TIMS-MS experiments (with a mass range of m/z 100–1700). The quadrupole isolation width was set to 2 and 3 Th and, for fragmentation, the collision energies varied between 31 and 52 eV depending on the precursor mass and charge. TIMS, MS operation and PASEF were controlled and synchronized using the control instrument software OtofControl 6.2 (Bruker Daltonik). LC-MS/MS data were acquired using the PASEF method with a total cycle time of 1.31 s, including 1 TIMS MS scan and 10 PASEF MS/MS scans. The 10 PASEF scans (100 ms each) containing, on average, 12 MS/MS scans per PASEF scan. Ion mobility-resolved mass spectra, nested ion mobility vs. m/z distributions, as well as summed fragment ion intensities were extracted from the raw data file with DataAnalysis 5.3 (Bruker Daltonik GmbH, Bremen, Germany). Signal-to-noise (S/N) ratio were increased by summations of individual TIMS scans. Mobility peak positions and peak half-widths were determined based on extracted ion mobilograms (±0.05 Da) using the peak detection algorithm implemented in the DataAnalysis software. Features detection were also performed using DataAnalysis 5.3 software and exported in .mgf format.

Protein Identification And Data Validation

Peptides were identified using the MASCOT software (version 2.5.1; Matrix Science, London, UK) against the Uniprot Bos taurus database (May 2019, 23523 sequences) using its automatic decoy database search to calculate a false discovery rate (FDR). The parameters used for database searches included trypsin as enzyme (one missed cleavage allowed), carbamidomethylcysteine as fixed modification, oxidation of methionine and N-terminal protein acetylation as variable modifications. Monoisotopic mass was considered and mass tolerance was set at 15 ppm for MS ions and 0,05 Da for MS/MS ions. Mascot results from the target and decoy databases were incorporated to Scaffold Q+ software (version 5.0.1, Proteome Software, Portland, USA, www.proteomesoftware.com). Threshold for peptide and protein identification were set to 95.0 % as specified by the Peptide Prophet algorithm⁶⁴ and the Protein Prophet algorithm⁶⁵.

Label-free Protein Quantification And Statistical Analysis

All proteins containing at least two unique peptides (False discovery rate (FDR) < 0.01%) were considered for protein quantification. The abundance of proteins was assessed using the Spectrum Count quantitative method and the “weighted spectra” option of the Scaffold software, as previously described⁶³. The weighted spectra option of Scaffold is based on the assignment of each peptide to a weight according to whether they are shared or not between identified proteins. The normalization of data was achieved by multiplying each spectrum count in each sample by the average count over the sample’s total spectrum count.

Statistical analysis was assessed on proteins with on average 2 normalized weighted spectra (NWS) or more in at least one condition. To obtain an overview of proteomic data, principal component analysis (PCA) of all samples was carried out using the FactoMineR and ggplot2 packages of RStudio software (version 1.4.1106). Analysis of variance on biological replicates and a heatmap representation of differentially abundant proteins (DAPs; p-value <0.05) were carried out using the FactoMineR and gplot packages of RStudio. Pair-wise comparison between stage ⋅ region ⋅ side conditions were evaluated by Student t-tests using the Scaffold software. Proteins were considered as differentially abundant with a t-test p-value ≤ 0.05 and a min fold-change ratio of 2. Lists of top-20 differentially abundant proteins are the quantitatively major proteins (based on NWS) at one given stage, side or region.

Prediction Of Secretory Pathways And Functional Enrichment Analysis

Identification of signal peptide-containing proteins and prediction of unconventional protein secretion (UPS) was first carried out using the Outcyte 1.0 online tool (http://www.outcyte.com/). The list of gene names was used as input and the default score threshold of 0.5 was considered for UPS prediction. Proteins having a transmembrane domain and predicted to be intracellular were also predicted by this software. The identification of signal peptide sequences was then checked on all proteins using the SignalP online tool (version 5.0; http://www.cbs.dtu.dk/services/SignalP/). Proteins identified were also compared to those previously reported in OF-derived EVs in cattle¹⁹ and cats²⁰ and to OEC-derived EVs in pigs²¹. The gene enrichment analysis was performed separately for each pairwise comparison except for region and stage effects for which the analysis was restricted to the ipsilateral side of ovulation. The gene lists of DAPs were first imported in the Database for Annotation, Visualization and Integrated Discovery (DAVID version 6.8) for gene ontology (GO) analysis. Overrepresented molecular functions (MF) and biological processes (BP) GO terms with a p-value < 0.05 were considered as significant. The KEGG pathways associated with DAPs were further analyzed using the Proteomaps software (http://bionic-vis.biologie.uni-greifswald.de/). Proteomaps graphics were generated from NWS values of overabundant DAPs in each pair-wise comparison. To evaluate the potential roles of DAPs in oviductal events, the Metascape Membership analysis tool was used (https://metascape.org/)²². Enrichment analysis of DAPs matching the membership terms “sperm”, “oocyte”, “ovulation”, “fertilization” and “embryo” were retained. The Homo sapiens genome was used as background dataset for all analyses as it is better annotated than the Bos taurus genome.

Data Availability

We comply with the publication's requirements for sharing materials of the journal. Upon acceptance, raw data will be available in PRIDE

https://www.ebi.ac.uk/pride/markdownpage/submitdatapage

Proteins identified in the bovine OF and predicted secretory pathways

Thirty pools of OF were made (4 cows per pool, 8 region × stage × side conditions) and analyzed by nanoLC-MS/MS. A total of 3,760 proteins were identified in the bovine OF (Supplementary Data S1). The top 40 most abundant proteins quantified in the OF included endoplasmin (HSP90B1), heat shock protein 90 alpha (HSP90AA1), HSPA8, myosin heavy chain 9 (MHY9), albumin (ALB) and oviductin (OVGP1), among others (see top-40 most abundant proteins in Supplementary Figure S1). To determine if some proteins were specific to one region, stage or side, the overlap between conditions was evaluated. Overall, 77% (2,894/3,760) of proteins were shared between conditions. Of the remaining 23%, only 34 proteins were quantified with more than 5 normalized weighted spectra (NWS, mean values) in one specific region, side or stage (see Supplementary Data S1).

To determine the potential pathways by which proteins were exported into the OF, peptide signal sequences were identified and unconventional secretion pathways were predicted using the Outcyte 1.0 and SignalP 5.0 server on the list of protein gene names. Overall, 37% (1,373/3,738) of the proteins identified were shown as potentially secreted proteins. Of those, 10% (360/3,738) contained a signal peptide and were predicted to be classically secreted while 27% (1,013/3,738) were predicted to be potentially secreted via unconventional pathways (see the prediction scores and cleavage site positions in Supplementary Data S1). In addition, 41% (1,551/3,738) of the proteins identified in the bovine OF were previously reported in EVs derived from oviduct epithelial cells (see the list of proteins and origin in Supplementary Data S1)^19–21. Of the 1,551 proteins potentially originated in oEVs, 498 (13%) were previously predicted as potentially secreted by the Outcyte and SignalP tools.

Spatiotemporal Regulation Of The Of Proteome And Functional Analysis

To obtain an overview of proteomic data, a principal component analysis (PCA) of proteins quantified with min 2 normalized weighted spectra (NWS) in at least one condition was carried out (3009 proteins). PCA showed a clear separation between the ampullary and isthmic samples (Figure 1a). Two samples in the Pre-ov ipsilateral isthmus group had to be discarded from further analysis because they had a low MS signal resulting in 50% less proteins identified than in all other samples.

To evaluate more specifically which factors among the region, the peri-ovulatory stage and the proximity of the ovulating ovary most affected the OF proteome, a hierarchical clustering of differentially abundant proteins (DAPs, ANOVA p-value ≤ 0.05) was performed. The heatmap representation of this analysis confirmed that the main differences in protein abundance were between the anatomical regions (ampulla vs. isthmus), followed by the proximity of the ovulating ovary (ipsilateral vs. contralateral) and finally the stage around ovulation time (Pre-ov vs. Post-ov) (Figure 1b). Based on these findings, statistical analyses comparing ampulla vs. isthmus, ipsilateral vs. contralateral and Pre-ov vs. Post-ov datasets were performed. Overall, pair-wise comparisons (t-test p-value ≤ 0.05) considering a fold-change ratio ≥ 2 evidenced highest numbers of DAPs when comparing ampulla and isthmus (237-354 DAPs, i.e., 8-12% of quantified proteins), followed by those between ipsi- and contralateral sides (83-248 DAPs, i.e. 3-8% of quantified proteins) and finally between the Pre-ov and Post-ov stages (50-67 DAPs, i.e. 2% of quantified proteins; see Figure 2). Interestingly, the post-ovulatory isthmus, i.e., the place and time of early embryo development, displayed the highest number of DAPs when comparing the ipsilateral with contralateral oviducts. Top-20 DAPs between regions, sides of ovulation and stages are shown in Figures 3, 4 and 5 respectively (see Supplementary Data S2-4 for all DAPs with accession numbers, molecular weight, p-value and fold-changes).

To evaluate the functions underlying those variations, an enrichment analysis of genes corresponding to DAPs was undertaken in the Database for Annotation, Visualization and Integrated Discovery (DAVID). This analysis for the effect of the region and stage was performed only on samples from the ipsilateral side. Analysis of the gene ontology (GO) terms associated with DAPs between regions indicated that 67 and 68 biological processes (BP) and molecular functions (MF) were overrepresented (p-value ≤ 0,05) at Pre-ov and Post-ov stages, respectively (Supplementary Data S2). Of those, 21 were shared between stages, including protein binding, cell-cell adhesion, protein folding, calcium ion binding and carbohydrate metabolic process as the most significant at both Pre- and Post-ov (in yellow in Supplementary Data S2). BP/MF overrepresented only at Pre-ov included protein disulfide isomerase activity, response to endoplasmic reticulum stress, redox homeostasis, negative regulation of apoptotic process, and microtubule motor activity (in green in Supplementary Data S3). BP/MF overrepresented only at Post-ov included calcium-dependent protein binding, calmodulin binding and cellular response to catecholamine stimulus (in blue in Supplementary Data S3). To further decipher the role of DAPs, Proteomaps were carried out on overabundant proteins in each region of the ipsilateral oviduct. At Pre-ov, proteins overabundant in the isthmus were mainly involved in protein processing in endoplasmic reticulum and protein export, while proteins overabundant in the ampulla were more specifically involved in the immune and endocrine systems in addition to the isthmic pathways (Figure 6). Proteomaps generated at post-ovulatory stages displayed globally the same pathways (Supplementary Figure S2). To point out proteins potentially involved in the reproductive events taking place in the oviduct, proteins overabundant in the ipsilateral isthmus or ampulla were further analyzed through the Metascape²² membership tool. Proteins associated with the membership terms ‘sperm’, ‘oocyte/ovulation’, ‘fertilization’ and ‘embryo’ are shown in Figure 7 (see the exhaustive list of GO terms in Membership analysis Metascape Supplementary Data S2).

Regarding the regulating role of the proximity of the ovulating ovary, analysis of the GO terms associated with DAPs between ipsilateral and contralateral oviducts evidenced protein binding as the only enriched MF shared between conditions (Supplementary Data S3). In the ampulla, overrepresented BP/MF included protein binding, endosomal transport, protein transport and multivesicular body assembly among the most significant GO terms at both stages. In the isthmus, the highest number (61) of specifically enriched GO terms was evidenced at Post-ov including mucus secretion, epidermal growth factor (EGF) receptor signaling pathway, ERBB2 signaling pathway and 1-phosphatidylinositol (PI)-3-kinase activity as the most significant ones. In addition, a number of DAPs between sides were associated with the membership terms related to early reproductive events (Figure 7B and Supplementary Data S3).

Regarding the regulating role of the stage, analysis of the GO terms associated with DAPs between Pre-ov and Post-ov samples in the ipsilateral oviduct evidenced 6 and 5 overrepresented BP/MF in the ampulla and isthmus, respectively, of which only protein binding was shared (Supplementary Data S4). BP overrepresented exclusively in the ampulla included regulation of cilium assembly, intraciliary transport and negative regulation of cell proliferation while artery smooth muscle contraction and vasodilatation were overrepresented only in the isthmus (see Supplementary Data S4). The DAPs associated with membership terms related to oviduct events are shown in Supplementary Data S4.

The present study provides an extensive analysis of the bovine OF proteome considering three factors of variation, i.e., the anatomical oviduct region, the peri-ovulatory stage and the proximity of the ovulating ovary. The main findings of this study are that i) the OF proteome varies in abundance between ampulla and isthmus, in accordance with our hypothesis that each segment of the oviduct responds to specific requirements of the gametes and embryos; ii) the region of the oviduct is more determining than the proximity of the ovulating ovary and the peri-ovulatory stage in the regulation of the OF proteome; iii) OF proteins varying in abundance are involved in a broad variety of biological functions, including protein binding, response to stress, calcium homeostasis, and lipid and carbohydrate metabolism, which may play a key role in the early gamete/embryo interactions.

More than 3,700 proteins were identified in the bovine OF, which is the most comprehensive OF proteome published so far in mammals^9,23−26. From the total of protein identified in our study, 37% were predicted as potentially secreted and 10% contained a signal peptide. These proportions of secreted proteins are comparable to those previously reported in the OF in cattle^23,24 and sheep²⁷. It is unclear how intracellular proteins could be exported in high proportions to the oviductal lumen. Considering the very low volume of OF in both ampulla (< 50 µL) and isthmus (< 10 µL) and to avoid cell damage, the samples were collected by gentle flushing and centrifugations were carried out to eliminate cell debris. The high rate of renewal of the oviduct epithelium during estrus^28,29 may lead to intracellular release of proteins into the OF. Moreover, and in the light of our results, it is very likely that intracellular proteins could be derived from apocrine or non-canonical secretory pathways, as it has been shown for annexins (ANXAs)^19,30. In accordance, a high proportion of intracellular proteins were previously identified in EVs derived from OF^19–21.

The most abundant proteins quantified in the OF included several HSPs (HSP90B1, HSP90AA1, HSP90AB1, HSPA8, HSPA5 also known as GRP78), CD109 molecule, complement C3, myosins (MYH9, MYH14) and the oviduct-specific glycoprotein oviductin (OVGP1), as previously reported^9,24,25. The majority (77%) of proteins quantified were detected in all conditions. The remaining 23% had on average low levels of detection (NWS < 5), suggesting that these proteins could have been detected with a higher sensitivity. Although, little attention has been paid to this small set of proteins, we should not rule out that they may have an important role in the oviduct. The use of mass spectrometer with higher sensitivity could determine whether these proteins were only present in a specific segment of the oviduct or stage, although in small amounts.

The region-related changes in the OF proteome are presented for the first time, providing valuable information on the specific requirements of gametes and embryos. Furthermore, the region was the greatest source of variation in the oviduct proteome. These results are in agreement with previous transcriptomic data evidencing differentially expressed genes between the ampulla and isthmus in cattle^10–13. Also in accordance, the abundance of specific proteins was differently regulated in each region after insemination in rabbits³¹ and between pregnant and cyclic cows¹⁴. Differences in the vascularization³² and in the proportion of ciliated cells²⁹ between isthmus and ampulla could explain this variation in OF proteome. It is well known that sex steroid hormones are major regulators of oviduct physiology and secretory activity⁶. To our knowledge, the region-specific concentrations of progesterone and estradiol in the oviduct have not been reported to date. However, the expression of progesterone receptors was reported to vary according to the oviduct region: transcripts for nuclear receptors (PGR) were more abundant in the isthmus while those for progesterone membrane component 2 (PGRMC2) were higher in the ampulla in cattle¹⁰. Furthermore, it was shown in gilts that each oviductal region responded differently to estradiol, leading to different regional patterns of secretion³³. Taken together, the observed region-specific profiling of the OF proteome could result from differences in blood flow, proportion of secretory cells, gene expression and steroid hormone receptivity between ampulla and isthmus.

After mating or insemination, sperm cross the uterus and reach the isthmus where they interact with both OECs and oviduct secretions; these interactions are believed to create an optimal environment for sperm storage until ovulation^2,34. Of the BP/MF overrepresented exclusively at Pre-ov, proteins involved in cell redox homeostasis (including PRDX3, PRDX4, PDIA3, PDIA6 among others), negative regulation of apoptotic process (incl. HSPA5, HSP90B1, RGN) and microtubule motor activity (incl. DNAH1, KLC3, KIF21A) may be considered as possible candidates to provide an optimal environment for sperm survival and motility until fertilization. Of those protein candidates, HSP90B1, HSPA5 and PDIA3 were among the most abundant proteins in the isthmus (NWS ≥ 100) and were also 2 to 3 times more abundant than in the ampulla at Pre-ov. HSP90B1, PDIA3 and HSPA5 were previously identified as sperm-binding proteins in the OF in cattle³⁵. In humans, the inhibition of PDIA3 by specific antibodies reduced sperm ability to bind to the zona pellucida³⁶. In addition, both HSP90B1 and HSPA5 promoted sperm-zona binding in human³⁷ and pigs³⁸. Recombinant HSPA5 enhanced human sperm calcium influx, one step of sperm capacitation³⁹. Calreticulin (CALR) was another protein found 2,5 times more abundant in the isthmus compared to the ampulla at Pre-ov. Calreticulin is a calcium-binding protein required for cumulus oocyte complex development⁴⁰. Moreover, addition of recombinant CALR in the fertilization medium blocked polyspermy in pigs⁴¹.

After ovulation, sperm undergo capacitation and migrate gradually toward the ampulla where fertilization takes place, then embryo development starts in the isthmus. Hence, we focused on identified proteins that could play a key role in supporting fertilization, embryo transport and development. Numerous DAPs between regions (incl. HSPA5, HSP90B1, ANXA3, ANXA6, MYO1D, MYO1B, PHKB) were involved in calcium homeostasis, including calcium ion binding and calcium-dependent protein binding. Calcium is well known to play major roles in sperm functions such as motility, capacitation and acrosome reaction⁴². In addition, several proteins (GNAS, GNG2, GNB1) were involved in cell response to catecholamines. Catecholamines have been implicated in the spontaneous contractions of the bovine oviduct⁴³ and may help embryo transport towards the isthmus. Among the overabundant proteins in the ampulla at Post-ov, oviductin (OVGP1) was among the most abundant proteins (NWS > 150) and 2.7 times more abundant than in the isthmus. Similarly, OVGP1 transcripts were more abundant in the ampulla than in the isthmus at day 4 post-ovulation in cows¹². OVGP1 is an oviduct-specific glycoprotein secreted around the time of ovulation that plays major roles in zona pellucida hardening and control of polyspermy in mammals^44,45. Addition of recombinant OVGP1 in a fertilization medium improved sperm-zona pellucida binding in hamster⁴⁶ and humans⁴⁷. Moreover, purified or recombinant OVGP1 positively influenced sperm viability, motility and capacitation in buffaloes⁴⁸ and humans⁴⁷. Furthermore, numerous proteins overabundant in the ampulla were involved in the immune system (incl. CFH, CFD, PLAT, PLTP). Complement factor H (CFH) acts as a soluble inhibitor of complement and plays as such an essential role in maintaining a well-balanced immune response. CFH was reported in boar seminal plasma and it was shown that sperm possessing CFH in their outer acrosomal region evaded complement attack in the female tract⁴⁹. The mRNAs for phospholipid transfer protein (PLTP) were reported in the oviduct epithelium of mice and increased in abundance in the presence of embryos but not in pseudo-pregnant mice, suggesting important roles during early development⁵⁰.

Beyond the oviduct region, the proximity of the ovulating ovary was another factor regulating OF protein abundance. The regulatory effect of the POF/CL on OF proteome around the time of ovulation was previously described in cattle⁹ and horse⁵¹. Differentially expressed genes between sides were also identified in the ampulla and isthmus on day 3 post-estrus in cattle⁵². The proportion of secretory cells was reported to be equivalent between ipsilateral and contralateral oviducts⁵³. However, highly asymmetrical concentrations of sex steroid hormones were reported in the bovine oviduct, with up to 3,5x more estradiol before ovulation⁵⁴ and up to 16 times more progesterone after ovulation⁷ in the ipsilateral oviduct than in the contralateral one. An asymmetry in the local vascularization and anatomy of the oviduct were also reported, with a thicker, more edematous and more transparent wall in the ipsilateral side⁵⁵. Overall, those differences in monovular species could probably be involved in the asymmetrical regulation of oviduct secretions.

The POF/CL effect was expected to be an important regulatory factor, preparing the suitable environment for early embryo development. In accordance, the Post-ov isthmus was the most affected by the proximity of the POF/CL, with 250 DAPs compared with 83-177 DAPs in other pair-wise comparisons between sides of ovulation. In line with the present data, transcriptomic analysis of bovine OECs identified 13 times more differentially expressed genes between sides of ovulation in the isthmus than in the ampulla on Day 3 post-estrus⁵². In a previous proteomic study in cattle, the proportion of DAPs in the entire oviduct was also higher at Post-ov than at Pre-ov⁹. The signaling pathways most significantly altered by the proximity of ovulating ovary included the EGFR, ERBB2 and PI3-kinase (PI3K) pathways after ovulation. Recent results in dogs show that exosomes secreted by oviduct cells mediate the EGFR/mitogen-activated protein kinase (MAPK) signaling pathway, which is one of the prerequisite pathways for further development, in cumulus cells⁵⁶. The PI3K pathway plays also pivotal roles in numerous reproductive processes⁵⁷. The PI3K-Akt pathway was one of the most affected pathways in the oviduct after mating, as shown in a transcriptomic study in pigs⁵⁸.

Studies on the regulation of the OF proteome close to the ovulation time are still scarce because most studies compared the estrus/follicular phase to the luteal phase, with a time interval of around 5-12 days between conditions^11,27,59. The time relative to ovulation (Pre-ov vs. Post-ov) was the factor with least variation between samples with only 2% of DAPs identified in the ampulla and isthmus. Higher rates of DAPs (13-17%) were previously evidenced in the bovine OF between stages but with no region separation⁹. Recently, differences before and after ovulation were reported in the equine OF proteome²⁶. These authors reported 15 proteins downregulated in post-ovulatory OF compared to the pre-ovulatory counterpart, while 156 upregulated in post-ovulatory OF compared to pre-ovulatory OF. It is to note that the Pre-ov and Post-ov stages covered approximately 3-4 days around ovulation time and represented a limited time window relative to the 21 to 25-day length of the estrous cycle, which may explain the limited number of DAPs. Functions such as regulation of cilium assembly, intraciliary transport and negative regulation of cell proliferation were affected by the peri-ovulatory stage in the ampulla. The proportions of ciliated cells and of proliferative cells in the oviduct epithelium was reported to change during the estrous cycle in cows²⁹ with an increase during the follicular phase and a decrease in the luteal phase, which may account for this stage-specific enrichment.

One limitation of the present study is the lack of information on cows’ breed, age and metabolic status at the time of OF collection since the samples were collected at a slaughterhouse. Although, transport and OF collection were carried out on ice with limited time after animal death, however, we cannot exclude some protein modifications due to post-mortem collection. The OF proteome was obtained from a total of 32 animals, providing data representative for this species, yet, the samples were analyzed in pools of 4 animals, precluding any evaluation of individual variations.

In conclusion, this study provides a comprehensive characterization of the bovine OF proteome, revealing a fine-tune regulation of the OF according to the anatomical region of the oviduct, the proximity to the ovary and the stage of the cycle. The strong modulatory effect of the oviduct regions points out the anatomical region as the major regulator of OF proteome for the first time. Moreover, the in deep detail of the spatiotemporal proteomic changes of the OF around the ovulation window, enhance our understanding of the oviductal milieu supporting the early reproductive events. The proteins identified here paved the way for future functional experiments towards improving sperm and embryo biotechnologies.

Funding

This work was funded by INRAe and Agence Nationale de la Recherche under the grant number ANR-18-CE92-0049. This work was also supported by grants from Biogenouest, Infrastructures en Biologie Santé et Agronomie (IBiSA) and Conseil Régional de Bretagne awarded to Protim proteomics core facility.

Acknowledgements

We are grateful to Marc Chodkiewicz for the critical reading of this paper. The authors thank Thierry Delpuech and Albert Arnoult for collecting the oviducts at the slaughterhouse. The authors thank the laboratory of hormonal assay (PRC, INRAE, Nouzilly) for the assessment of P4 concentration in OF and pre-ovulatory follicular fluid.

Author contributions

MSD, PM and JS designed the study; CM and MSD wrote the manuscript; MSD and GT collected and prepared the OF; YL did the hormonal assays; RL and EC did the MS analyses; CM analyzed the proteomic data; CA made complementary data analysis ; CP supervised the MS analyses; all authors corrected the manuscript; MSD administrated and supervised the study.

Competing interests

The authors declare no conflict of interest.

Suarez, S. S. Mammalian sperm interactions with the female reproductive tract. Cell Tissue Res, 363, 185–194 (2016).
Mahé, C. et al. Sperm migration, selection, survival, and fertilizing ability in the mammalian oviduct. Biology of Reproduction, 105, 317–331 (2021).
Ded, L., Hwang, J. Y., Miki, K., Shi, H. F. & Chung, J. J. 3D in situ imaging of the female reproductive tract reveals molecular signatures of fertilizing spermatozoa in mice. eLife, 9, e62043 (2020).
Coy, P., García-Vázquez, F. A., Visconti, P. E. & Avilés, M. Roles of the oviduct in mammalian fertilization., 144, 649–660 (2012).
Li, S. & Winuthayanon, W. Oviduct: roles in fertilization and early embryo development. Journal of Endocrinology, 232, R1–R26 (2017).
Saint-Dizier, M., Schoen, J., Chen, S., Banliat, C. & Mermillod, P. Composing the Early Embryonic Microenvironment: Physiology and Regulation of Oviductal Secretions. International Journal of Molecular Sciences, 21, 223 (2020).
Lamy, J. et al. Steroid hormones in bovine oviductal fluid during the estrous cycle., 86, 1409–1420 (2016).
Lamy, J. et al. Metabolomic profiling of bovine oviductal fluid across the oestrous cycle using proton nuclear magnetic resonance spectroscopy. Reprod. Fertil. Dev, 30, 1021 (2018).
Lamy, J. et al. Regulation of the bovine oviductal fluid proteome., 152, 629–644 (2016).
Saint-Dizier, M., Sandra, O., Ployart, S., Chebrout, M. & Constant, F. Expression of nuclear progesterone receptor and progesterone receptor membrane components 1 and 2 in the oviduct of cyclic and pregnant cows during the post-ovulation period. Reprod Biol Endocrinol, 10, 76 (2012).
Cerny, K. L., Garrett, E., Walton, A. J., Anderson, L. H. & Bridges, P. J. A transcriptomal analysis of bovine oviductal epithelial cells collected during the follicular phase versus the luteal phase of the estrous cycle. Reprod Biol Endocrinol, 13, 84 (2015).
Gonella-Diaza, A. M. et al. Size of the Ovulatory Follicle Dictates Spatial Differences in the Oviductal Transcriptome in Cattle. PLoS ONE, 10, e0145321 (2015).
Maillo, V. et al. Spatial differences in gene expression in the bovine oviduct., 152, 37–46 (2016).
Rodríguez-Alonso, B. et al. Spatial and Pregnancy-Related Changes in the Protein, Amino Acid, and Carbohydrate Composition of Bovine Oviduct Fluid. IJMS 21, 1681(2020).
Franchi, A., Moreno-Irusta, A., Domínguez, E. M., Adre, A. J. & Giojalas, L. C. Extracellular vesicles from oviductal isthmus and ampulla stimulate the induced acrosome reaction and signaling events associated with capacitation in bovine spermatozoa. J Cell Biochem, 121, 2877–2888 (2020).
Kumaresan, A., Ansari, M. R. & Garg, A. Modulation of post-thaw sperm functions with oviductal proteins in buffaloes. Anim. Reprod. Sci, 90, 73–84 (2005).
Lopera-Vasquez, R. et al. Effect of bovine oviductal fluid on development and quality of bovine embryos produced in vitro. Reprod. Fertil. Dev, 29, 621 (2017).
Dashti, S., Zare Shahneh, A., Kohram, H. & Zhandi, M. & Dadashpour Davachi, N. Differential influence of ovine oviduct ampullary and isthmic derived epithelial cells on in vitro early embryo development and kinetic. Small Ruminant Research, 136, 197–201 (2016).
Almiñana, C. et al. Oviduct extracellular vesicles protein content and their role during oviduct–embryo cross-talk., 154, 253–268 (2017).
Ferraz, M., de Carothers, A. M. M., Dahal, A., Noonan, R., Songsasen, N. & M. J. & Oviductal extracellular vesicles interact with the spermatozoon’s head and mid-piece and improves its motility and fertilizing ability in the domestic cat. Sci Rep, 9, 9484 (2019).
Laezer, I. et al. Dynamic profile of EVs in porcine oviductal fluid during the periovulatory period., 159, 371–382 (2020).
Zhou, Y. et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun, 10, 1523 (2019).
Pillai, V. V., Weber, D. M., Phinney, B. S. & Selvaraj, V. Profiling of proteins secreted in the bovine oviduct reveals diverse functions of this luminal microenvironment. PLoS ONE, 12, e0188105 (2017).
Gegenfurtner, K. et al. Influence of metabolic status and genetic merit for fertility on proteomic composition of bovine oviduct fluid. Biology of Reproduction, 101, 893–905 (2019).
Papp, S. M. et al. A novel approach to study the bovine oviductal fluid proteome using transvaginal endoscopy., 132, 53–61 (2019).
Fernández-Hernández, P. et al. The Proteome of Equine Oviductal Fluid Varies Before and After Ovulation: A Comparative Study. Front. Vet. Sci, 8, 694247 (2021).
Soleilhavoup, C. et al. Proteomes of the Female Genital Tract During the Oestrous Cycle. Mol. Cell Proteomics, 15, 93–108 (2016).
Ghosh, A., Syed, S. M. & Tanwar, P. S. In vivo genetic cell lineage tracing reveals that oviductal secretory cells self-renew and give rise to ciliated cells. Development, 144, 3031–3041 (2017).
Ito, S., Kobayashi, Y., Yamamoto, Y., Kimura, K. & Okuda, K. Remodeling of bovine oviductal epithelium by mitosis of secretory cells. Cell Tissue Res, 366, 403–410 (2016).
Popa, S. J., Stewart, S. E. & Moreau, K. Unconventional secretion of annexins and galectins. Seminars in Cell & Developmental Biology, 83, 42–50 (2018).
Yu, H. et al. Identification of Rabbit Oviductal Fluid Proteins Involved in Pre-Fertilization Processes by Quantitative Proteomics., 19, 1800319 (2019).
Gutiérrez, H. A., Latorre, R., MartÍnez BagÁn, E. & SÁnchez Margallo, F. M. & López Albors, O. Quantitative Evaluation of the Vasculature Supplying the Oviduct in Pre-Pubertal and Sexually Mature Sows: Quantitative Evaluation of the Vasculature. Anat. Rec, 298, 1978–1983 (2015).
Buhi, W. C., Ashworth, C. J., Bazer, F. W. & Alvarez, I. M. In vitro synthesis of oviductal secretory proteins by estrogen-treated ovariectomized gilts. J. Exp. Zool, 262, 426–435 (1992).
Holt, W. V. & Fazeli, A. Sperm Storage in the Female Reproductive Tract. Annu Rev Anim Biosci, 4, 291–310 (2016).
Lamy, J. et al. Identification by proteomics of oviductal sperm-interacting proteins., 155, 457–466 (2018).
Wong, C. W. et al. The roles of protein disulphide isomerase family A, member 3 (ERp57) and surface thiol/disulphide exchange in human spermatozoa–zona pellucida binding. Hum. Reprod, 32, 733–742 (2017).
Marín-Briggiler, C. I. et al. Glucose-regulated protein 78 (Grp78/BiP) is secreted by human oviduct epithelial cells and the recombinant protein modulates sperm–zona pellucida binding. Fertility and Sterility, 93, 1574–1584 (2010).
Mondéjar, I., Martínez-Martínez, I., Avilés, M. & Coy, P. Identification of Potential Oviductal Factors Responsible for Zona Pellucida Hardening and Monospermy During Fertilization in Mammals.Biology of Reproduction89, (2013).
Lachance, C., Bailey, J. L. & Leclerc, P. Expression of Hsp60 and Grp78 in the human endometrium and oviduct, and their effect on sperm functions. Hum. Reprod, 22, 2606–2614 (2007).
Tokuhiro, K. et al. Calreticulin is required for development of the cumulus oocyte complex and female fertility. Sci Rep, 5, 14254 (2015).
Saavedra, M. D. et al. Calreticulin from suboolemmal vesicles affects membrane regulation of polyspermy., 147, 369–378 (2014).
Bernecic, N. C., Gadella, B. M., Leahy, T. & de Graaf S. P. Novel methods to detect capacitation-related changes in spermatozoa., 137, 56–66 (2019).
Kotwica, G. et al. The concentrations of catecholamines and oxytocin receptors in the oviduct and its contractile activity in cows during the estrous cycle., 60, 953–964 (2003).
Coy, P. et al. Hardening of the zona pellucida of unfertilized eggs can reduce polyspermic fertilization in the pig and cow., 135, 19–27 (2008).
Algarra, B. et al. The C-terminal region of OVGP1 remodels the zona pellucida and modifies fertility parameters.Scientific Reports6, (2016).
Yang, X., Zhao, Y., Yang, X. & Kan, F. W. K. Recombinant Hamster Oviductin Is Biologically Active and Exerts Positive Effects on Sperm Functions and Sperm-Oocyte Binding. PLOS ONE, 10, e0123003 (2015).
Zhao, Y. & Kan, F. W. K. Human OVGP1 enhances tyrosine phosphorylation of proteins in the fibrous sheath involving AKAP3 and increases sperm-zona binding. Journal of Assisted Reproduction and Genetics, 36, 1363–1377 (2019).
Choudhary, S. et al. Effect of recombinant and native buffalo OVGP1 on sperm functions and in vitro embryo development: a comparative study.Journal of Animal Science and Biotechnology8, (2017).
Sakaue, T. et al. Factor H in Porcine Seminal Plasma Protects Sperm against Complement Attack in Genital Tracts. Journal of Biological Chemistry, 285, 2184–2192 (2010).
Lee, K. F. et al. Phospholipid transfer protein (PLTP) mRNA expression is stimulated by developing embryos in the oviduct. J. Cell. Biochem, 95, 740–749 (2005).
Smits, K. et al. Proteome of equine oviducal fluid: effects of ovulation and pregnancy. Reprod. Fertil. Dev, 29, 1085 (2017).
Locatelli, Y. et al. Relative effects of location relative to the corpus luteum and lactation on the transcriptome of the bovine oviduct epithelium. BMC Genomics, 20, 233 (2019).
Sostaric, E. et al. Sperm binding properties and secretory activity of the bovine oviduct immediately before and after ovulation. Mol. Reprod. Dev, 75, 60–74 (2008).
Wijayagunawardane, M. P. et al. Local distributions of oviductal estradiol, progesterone, prostaglandins, oxytocin and endothelin-1 in the cyclic cow., 49, 607–618 (1998).
Kölle, S. et al. Ciliary Transport, Gamete Interaction, and Effects of the Early Embryo in the Oviduct: Ex Vivo Analyses Using a New Digital Videomicroscopic System in the Cow. Biology of Reproduction, 81, 267–274 (2009).
Lee, S. H., Oh, H. J., Kim, M. J. & Lee, B. C. Exosomes derived from oviduct cells mediate the EGFR/MAPK signaling pathway in cumulus cells. J Cell Physiol, 235, 1386–1404 (2020).
Acosta-Martínez, M. PI3K: An Attractive Candidate for the Central Integration of Metabolism and Reproduction.Front. Endocrin.2, (2012).
Álvarez-Rodríguez, M., Martinez, C. A., Wright, D. & Rodríguez-Martinez, H. The role of semen and seminal plasma in inducing large-scale genomic changes in the female porcine peri-ovulatory tract. Sci Rep, 10, 5061 (2020).
Bauersachs, S. et al. Monitoring gene expression changes in bovine oviduct epithelial cells during the oestrous cycle. Journal of Molecular Endocrinology, 32, 449–466 (2004).
Ireland, J. J., Murphee, R. L. & Coulson, P. B. Accuracy of Predicting Stages of Bovine Estrous Cycle by Gross Appearance of the Corpus Luteum. Journal of Dairy Science, 63, 155–160 (1980).
Adib, A. et al. Progesterone improves the maturation of male-induced preovulatory follicles in anoestrous ewes., 148, 403–416 (2014).
Rico, C. et al. Regulation of Anti-Müllerian Hormone Production in the Cow: A Multiscale Study at Endocrine, Ovarian, Follicular, and Granulosa Cell Levels1. Biology of Reproduction, 84, 560–571 (2011).
Banliat, C. et al. Identification of 56 Proteins Involved in Embryo–Maternal Interactions in the Bovine Oviduct. International Journal of Molecular Sciences, 21, 466 (2020).
Keller, A., Nesvizhskii, A. I., Kolker, E. & Aebersold, R. Empirical Statistical Model To Estimate the Accuracy of Peptide Identifications Made by MS/MS and Database Search. Anal. Chem, 74, 5383–5392 (2002).
Nesvizhskii, A. I., Keller, A., Kolker, E. & Aebersold, R. A Statistical Model for Identifying Proteins by Tandem Mass Spectrometry. Anal. Chem, 75, 4646–4658 (2003).

No competing interests reported.

Download PDF

Editorial decision: Major revision
21 Dec, 2021
Reviews received at journal
02 Dec, 2021
Reviewers agreed at journal
29 Nov, 2021
Reviewers agreed at journal
17 Nov, 2021
Reviewers invited by journal
13 Nov, 2021
Editor assigned by journal
27 Oct, 2021
Editor invited by journal
27 Oct, 2021
Submission checks completed at journal
27 Oct, 2021
First submitted to journal
15 Oct, 2021

You are reading this latest preprint version

Spatiotemporal Regulation of The Bovine Oviduct Fluid Proteome Around The Time of Ovulation

Status:

Version 1

Abstract

Figures

Introduction

Methods

Collection of bovine oviductal fluid (OF)

Preparation Of Of Samples For Proteomic Analyses

Nanoliquid Chromatography Coupled With Tandem Mass Spectrometry (Nanolc-ms/ms)

Protein Identification And Data Validation

Label-free Protein Quantification And Statistical Analysis

Prediction Of Secretory Pathways And Functional Enrichment Analysis

Results

Proteins identified in the bovine OF and predicted secretory pathways

Spatiotemporal Regulation Of The Of Proteome And Functional Analysis

Discussion

Declarations

References

Additional Declarations

Supplementary Files

Status:

Version 1