# Metabolomic Profile of Oviductal Extracellular Vesicles across the Estrous Cycle in Cattle

Julie Gatien, Pascal Mermillod, Guillaume Tsikis, Ophélie Bernardi, Sarah

Janati Idrissi, Rustem Uzbekov, Daniel Le Bourhis, Pascal Salvetti, Carmen

Alminana-Brines, Marie Saint-Dizier

► To cite this version:

Julie Gatien, Pascal Mermillod, Guillaume Tsikis, Ophélie Bernardi, Sarah Janati Idrissi, et al. Metabolomic Profile of Oviductal Extracellular Vesicles across the Estrous Cycle in Cattle. International Journal of Molecular Sciences, 2019, 20 (24), 10.3390/ijms20246339. [DOI: 10.3390/ijms20246339](https://doi.org/10.3390/ijms20246339)

HAL Id: hal-02623052

https://hal.inrae.fr/hal-02623052

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Metabolomic Profile of Oviductal Extracellular Vesicles across the Estrous Cycle in Cattle

Julie Gatien<fnref n="1" />, Pascal Mermillod<fnref n="2" /> ID, Guillaume Tsikis<fnref n="2" />, Ophélie Bernardi<fnref n="2" />, Sarah Janati Idrissi<fnref n="1" />, Rustem Uzbekov<fnref n="3" />,<fnref n="4" /> ID, Daniel Le Bourhis<fnref n="1" />, Pascal Salvetti<fnref n="1" />, Carmen Almiñana<fnref n="2" />,<fnref n="5" />,† and Marie Saint-Dizier<fnref n="2" />,<fnref n="6" />,*,† ID

Received: 6 December 2019; Accepted: 15 December 2019; Published: 16 December 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** Oviductal extracellular vesicles (oEVs) have been proposed as key modulators of gamete/embryo maternal interactions. The aim of this study was to examine the metabolite content of oEVs and its regulation across the estrous cycle in cattle. Oviductal EVs were isolated from bovine oviducts ipsilateral and contralateral to ovulation at four stages of the estrous cycle (post-ovulatory stage, early and late luteal phases, and pre-ovulatory stage). The metabolomic profiling of EVs was performed by proton nuclear magnetic resonance spectroscopy (NMR). NMR identified 22 metabolites in oEVs, among which 15 were quantified. Lactate, myoinositol, and glycine were the most abundant metabolites throughout the estrous cycle. The side relative to ovulation had no effect on the oEVs' metabolite concentrations. However, levels of glucose-1-phosphate and maltose were greatly affected by the cycle stage, showing up to 100-fold higher levels at the luteal phase than at the peri-ovulatory phases. In contrast, levels of methionine were significantly higher at peri-ovulatory phases than at the late-luteal phase. Quantitative enrichment analyses of oEV-metabolites across the cycle evidenced several significantly regulated metabolic pathways related to sucrose, glucose, and lactose metabolism. This study provides the first metabolomic characterization of oEVs, increasing our understanding of the potential role of oEVs in promoting fertilization and early embryo development.

**Keywords:** oviduct; fallopian tube; extracellular vesicles; exosomes; metabolomics; energy substrates; amino acids; NMR

## 1. Introduction

Exosomes and microvesicles, collectively known as extracellular vesicles (EVs), are membrane-enclosed particles able to transfer a complex selection of molecular compounds from one cell to another, and in this way play key roles in cell-to-cell communications [1]. EVs have been identified in all biological fluids and can be secreted by most cell types. In the female genital tract, EVs have been isolated from follicular, uterine, and oviductal fluids (OF) [2–5]. EVs from the OF (oEVs), also known as oviductosomes, have gained growing attention in recent years as they may

<footnote n="1">¹ Allice, 37380 Nouzilly, France; julie.gatien@alice.fr (J.G.); sarah.janatiidrissi@alice.fr (S.J.I.); daniel.lebourhis@alice.fr (D.L.B.); pascal.salvetti@alice.fr (P.S.)</footnote>

<footnote n="2">² Institut National de la Recherche Agronomique (INRA), CNRS 7247, University of Tours, IFCE, UMR85 Physiologie de la Reproduction et des Comportements, 37380 Nouzilly, France; pascal.mermillod@inra.fr (P.M.); guillaume.tsikis@inra.fr (G.T.); ophelie.bernardi@inra.fr (O.B.); carmen.alminanabrines@uzh.ch (C.A.)</footnote>

<footnote n="3">³ Faculty of Medicine, University of Tours, 37000 Tours, France; rustem.uzbekov@univ-tours.fr</footnote>

<footnote n="4">⁴ Faculty of Bioengineering and Bioinformatics, Moscow State University, 119991 Moscow, Russia</footnote>

<footnote n="5">⁵ VetSuisse Faculty, University of Zurich, 8057 Zurich, Switzerland</footnote>

<footnote n="6">⁶ Faculty of Sciences and Techniques, University of Tours, 37200 Tours, France</footnote>

<footnote>* Correspondence: marie.saint-dizier@inra.fr; Tel.: +33-247-427-508</footnote>

<footnote>† These authors contributed equally to this work.</footnote>
act as natural cargos, bringing key components from the maternal compartment into gametes and
embryos for the establishment of pregnancy [5–9]. Oviductal EVs are released by oviductal epithelial
cells (OECs), as shown from in vitro studies [5,7], and may originate from the pre-ovulatory follicle at
ovulation. Previous studies in mice showed that incubation of spermatozoa with oEVs activates a Ca2+
efflux pomp, PMCA4, which is important for the acquisition of sperm hyperactivated motility and
fertilizing competence [6]. In the feline species, oEVs were shown to fuse with the sperm membrane, to
improve sperm motility and fertilizing ability, and to prevent premature acrosome reaction in vitro [8].
Recently, we showed that oEVs can reduce polyspermy during in vitro fertilization of pig oocytes
without affecting the global fertilization rate [10]. We also demonstrated that the co-incubation of
in vitro produced embryos with oEVs increased blastocyst rates, extended embryo survival over time,
and improved embryo quality [5]. In mice, oEVs from donor oviduct fluid in the transfer medium
were also shown to improve birth rates after the transfer of in vitro produced embryos to recipient
mothers [11]. All together, these data support the idea that oEVs play important roles in sperm
capacitation, fertilization, and early embryo development.

To further investigate the compounds carried by oEVs potentially playing these roles, we analyzed
the proteomic and transcriptomic profiles of oEVs across the estrus cycle revealing proteins, mRNAs,
and various types of small non-coding RNAs with important roles in early reproductive events [5,9].
To date, only a few studies have reported the molecular content of oEVs (reviewed in [12]). In the
present study, we focused on the metabolomic profile of oEVs, the least studied component of EVs
in all research areas, with still no available data on oEVs. If we consider that metabolites represent
an intermediate or end point of any cellular process, they can reflect the phenotype printout of a cell
at a specific moment. Thus, monitoring metabolic changes in any biological fluid under a variety
of physiological or pathological conditions may be seen as the goldmine of biomarkers of health
condition [13]. Moreover, determining the metabolite composition of EVs can reveal important aspects
in cell-to-cell communication.

On the other hand, it is widely accepted that the ovarian steroid hormones progesterone (P4)
and estradiol-17β (E2) regulate all components of the oviduct secretory activity [14] including oEV
content [9] and also influence the metabolic content of the OF [15]. Moreover, in mono-ovular species
like the bovine and human, the side of ovulation may also affect the molecular composition of the
proximal OF by topical hormonal regulations and putative input of the ovulatory follicle. Thus, the
objectives of the present study were to analyze the metabolomic profile of oEVs and examine whether
the oEV metabolite content is under the hormonal control of the estrous cycle and is dependent on the
side relative to ovulation.

**2. Results**

2.1. *Characterization of Oviductal Extracellular Vesicles Used for Metabolomic Profiling*

The oEVs samples collected from cyclic cows and used for NMR analysis are shown in Table 1.

Transmission electron microscopy (TEM) observation confirmed the presence of oEVs at the four
stages of cycle analyzed (Figure 1). Oviductal EVs samples contained a majority of vesicles resembling
exosomes (size of 30–100 nm; 76.4 ± 0.8%) and a small population of vesicles resembling microvesicles
(>100 up to 500 nm: 23.7 ± 0.8%), with similar proportions among stages, as we previously showed in
cattle [9]. Western blotting confirmed the expression of exosomal cytosolic markers (CD81, HSP70, and
ANXA1) in the oEV samples (Figure 1).
<table><caption>**Table 1.** Samples of bovine oviductal extracellular vesicles (oEVs) used for proton nuclear magnetic resonance spectroscopy (NMR) analysis.</caption><thead><tr><th>Stage of the Cycle</th><th>No. Animals <sup>1</sup></th><th>Side Relative to Ovulation</th><th>Number of OF Pools</th><th>Volume (µL) of OF before oEV Isolation</th><th>Protein Concentration (mg mL<sup>-1</sup>)</th></tr></thead><tbody><tr><td rowspan="2">Post-ovulatory</td><td rowspan="2">54</td><td>Ipsilateral</td><td>4</td><td>320–420</td><td>10.4–13.6</td></tr><tr><td>Contralateral</td><td>4</td><td>288–385</td><td>9.7–17.4</td></tr><tr><td rowspan="2">Mid-luteal</td><td rowspan="2">40</td><td>Ipsilateral</td><td>3</td><td>355–400</td><td>11.4–12.6</td></tr><tr><td>Contralateral</td><td>3</td><td>305–410</td><td>11.1–13.5</td></tr><tr><td rowspan="2">Late-luteal</td><td rowspan="2">34</td><td>Ipsilateral</td><td>3</td><td>315–470</td><td>9.4–10.3</td></tr><tr><td>Contralateral</td><td>3</td><td>350–470</td><td>13.1–15.9</td></tr><tr><td rowspan="2">Pre-ovulatory</td><td rowspan="2">40</td><td>Ipsilateral</td><td>3</td><td>450–495</td><td>10.0–15.2</td></tr><tr><td>Contralateral</td><td>3</td><td>500–580</td><td>12.0–13.4</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Characterization of bovine oviductal extracellular vesicles (oEVs). Representative images of exosomes (30–100 nm) and microvesicles (>100 nm) in oEV preparations observed by transmission electron microscopy (TEM) across the estrus cycle (a) and Western blotting characterization of bovine oEVs for known exosomal protein markers (b). A pool of samples from four different stages was used, showing that oEVs were positive for CD81, HSP70, and ANXA1.</figcaption></figure>

## 2.2. Metabolite Profile of oEVs across the Estrous Cycle

Proton nuclear magnetic resonance spectroscopy identified 22 metabolites in oEVs (Table 2). The overall oEV metabolite composition showed various classes of metabolites including essential and non-essential amino acids and various energy substrates among other categories (Table 2).

The concentrations of the 15 metabolites could be quantified (Table 3). Among the seven metabolites that could not be quantified, five (glucose, glutamate, leucine, proline, and threonine) were due to a large overlap of their signals with other metabolites and two (formate and ethanolamine) due to signals below the detection threshold. From the 15 metabolites quantified, lactate, myo-inositol, and glycine were the most abundant metabolites throughout the estrous cycle, reaching on average 25–51.5 nmoL mg⁻<fnref n="1" /> of oEV protein.

<footnote n="1">¹ Number of animals used to obtain 3–4 pools of oviductal fluid (OF) per stage × side with a minimum volume of 280 µL/pool that was used for oEV isolation.</footnote>
<table><caption>**Table 2.** Metabolites identified in bovine oviductal extracellular vesicles across the estrous cycle, ordered by chemical classes ¹.</caption><thead><tr><th>Metabolite</th><th>Chemical Taxonomy<br>Subclass<br>Description</th><th>Identified<br>(I)/Quantified (Q)</th><th>Function</th></tr></thead><tbody><tr><td>Valine</td><td>Amino acids, peptides, and<br>analogues<br>(Proteinogenic)<br>Essential amino acids</td><td>I/Q</td><td>Involved in stress, energy and<br>muscle metabolism</td></tr><tr><td>Threonine</td><td>Amino acids, peptides, and<br>analogues<br>(Proteinogenic)<br>Essential amino acids</td><td>I</td><td>Involved in biosynthesis of proteins</td></tr><tr><td>Methionine</td><td>Amino acids, peptides, and<br>analogues<br>(Proteinogenic)<br>Essential amino acids</td><td>I/Q</td><td>Required for normal growth<br>and development</td></tr><tr><td>Leucine</td><td>Amino acids, peptides, and<br>analogues<br>(Proteinogenic)<br>Essential amino acids</td><td>I</td><td>Involved in biosynthesis of proteins, stress,<br>energy, and muscle metabolism.<br>Stimulates insulin release</td></tr><tr><td>Isoleucine</td><td>Amino acids, peptides, and<br>analogues<br>(Proteinogenic)<br>Essential amino acids</td><td>I/Q</td><td>Involved in stress, energy and<br>muscle metabolism.</td></tr><tr><td>Proline</td><td>Amino acids, peptides, and<br>analogues<br>Conditionally essential<br>amino acids<br>(Proteinogenic)</td><td>I</td><td>Involved in biosynthesis of proteins</td></tr><tr><td>Glycine</td><td>Amino acids, peptides, and<br>analogues<br>Non-essential amino acids<br>(Proteinogenic)</td><td>I/Q</td><td>Involved in the body's production of DNA,<br>phospholipids, and collagen, and in release<br>of energy.</td></tr><tr><td>Alanine</td><td>Amino acids, peptides, and<br>analogues<br>Non-essential amino acids<br>(Proteinogenic)</td><td>I/Q</td><td>One of the most important amino acids<br>released by muscle, functioning as a major<br>energy source. Regulator of glucose<br>metabolism, lymphocyte reproduction<br>and immunity.</td></tr><tr><td>Glutamate</td><td>Amino acids, peptides, and<br>analogues<br>(Non-essential alpha Amino<br>Acids)</td><td>I</td><td>Involved in biosynthesis of proteins<br>Role as neurotransmitter, a chemical used<br>by nerve cells to send signals to other cells.</td></tr><tr><td>Creatine</td><td>Amino acids, peptides, and<br>analogues<br>(endogenous amino acid:<br>synthesized from arginine,<br>glycine, and methionine)</td><td>I/Q</td><td>Role in energy metabolism. Responsible for<br>the production of ATP in skeletal muscle<br>through the process of oxidative<br>phosphorylation inside the mitochondria.</td></tr><tr><td>Ethanolamine</td><td>Amines</td><td>I</td><td>Widely distributed in biological tissue and<br>is a component of lecithin.</td></tr><tr><td>Carnitine</td><td>Quaternary ammonium salts<br>Non-proteinogenic amino<br>acids</td><td>I/Q</td><td>Important in providing energy to muscles<br>Described as a vitamin, an amino acid, or<br>essential metabolite</td></tr><tr><td>Choline</td><td>Quaternary ammonium salts<br>(Essential Vitamin)</td><td>I/Q</td><td>Considered an essential vitamin.<br>Precursor for the neurotransmitter<br>acetylcholine, which is involved in many<br>functions including memory and<br>muscle control.</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>Metabolite</th><th>Chemical Taxonomy<br>Subclass<br>Description</th><th>Identified<br>(I)/Quantified (Q)</th><th>Function</th></tr></thead><tbody><tr><td>Lactate</td><td>Alpha hydroxy acids and<br>derivatives<br>(Energy Substrate)</td><td>I/Q</td><td>Plays a role in several biochemical<br>processes and is produced in the muscles<br>during intense activity</td></tr><tr><td>Myoinositol</td><td>Alcohols and polyols</td><td>I/Q</td><td>Involved in the Inositol phosphate<br>metabolism and the Phosphatidylinositol<br>signaling system.</td></tr><tr><td>Glucose-1-phosphate</td><td>Carbohydrates and<br>carbohydrate conjugates<br>(Energy Substrate)</td><td>I/Q</td><td>Glycogenolysis produces<br>glucose-1-phosphate and no energy.</td></tr><tr><td>Maltose</td><td>Carbohydrates and<br>carbohydrate conjugates<br>(Energy Substrate)</td><td>I/Q</td><td>Maltose can be broken down into two<br>glucose molecules</td></tr><tr><td>Glucose</td><td>Carbohydrates and<br>carbohydrate conjugates<br>(Energy Substrate)</td><td>I</td><td>Primary source of energy</td></tr><tr><td>Succinate</td><td>Dicarboxylic acids and<br>derivatives</td><td>I/Q</td><td>Component of the citric acid or TCA cycle<br>and is capable of donating electrons to the<br>electron transfer chain.</td></tr><tr><td>Formate</td><td>Carboxylic acids</td><td>I</td><td>Essential intermediary metabolite in<br>folate-mediated one-carbon metabolism<br>Responsible for both metabolic acidosis and<br>disrupting mitochondrial electron transport<br>and energy production.</td></tr><tr><td>Acetate</td><td>Carboxylic acid derivatives</td><td>I/Q</td><td>Acetate in the form of acetyl CoA is used in<br>metabolism to yield chemical energy.</td></tr><tr><td>Acetone</td><td>Carbonyl compounds</td><td>I/Q</td><td>One of the ketone bodies produced during<br>ketoacidosis. Since ketosis develops under<br>serious metabolic circumstances, all the<br>mechanisms that balance or moderate the<br>effects of ketosis enhance the chance<br>for survival.</td></tr></tbody></table>

<table><caption>**Table 3.** Concentrations of metabolites in bovine oviductal extracellular vesicles according to the stage of the estrous cycle and the side relative to ovulation (ipsilateral vs. contralateral). Data are given as means ± SEM of nmoL mg⁻<fnref n="1" /> of oEV proteins.</caption><thead><tr><th rowspan="2"></th><th colspan="2">Pre-Ovulatory</th><th colspan="2">Post-Ovulatory</th><th colspan="2">Mid-Luteal</th><th colspan="2">Late-Luteal</th></tr><tr><th>Ipsi</th><th>Contra</th><th>Ipsi</th><th>Contra</th><th>Ipsi</th><th>Contra</th><th>Ipsi</th><th>Contra</th></tr></thead><tbody><tr><td>Acetate</td><td>2.5 ± 0.6</td><td>2.8 ± 0.1</td><td>2.7 ± 0.2</td><td>2.7 ± 0.2</td><td>2.9 ± 0.3</td><td>2.8 ± 0.0</td><td>2.8 ± 0.2</td><td>2.8 ± 0.2</td></tr><tr><td>Acetone</td><td>0.6 ± 0.2</td><td>0.7 ± 0.0</td><td>0.9 ± 0.0</td><td>0.9 ± 0.0</td><td>1.0 ± 0.1</td><td>0.9 ± 0.0</td><td>1.0 ± 0.1</td><td>0.9 ± 0.0</td></tr><tr><td>Alanine</td><td>16.7 ± 3.7</td><td>18.1 ± 0.9</td><td>15.5 ± 1.6</td><td>14.9 ± 1.4</td><td>15.5 ± 2.3</td><td>13.8 ± 0.3</td><td>13.4 ± 0.9</td><td>13.8 ± 1.6</td></tr><tr><td>Carnitine</td><td>4.8 ± 0.8</td><td>6.1 ± 0.2</td><td>4.3 ± 0.4</td><td>4.9 ± 0.6</td><td>5.5 ± 0.6</td><td>5.7 ± 0.9</td><td>5.5 ± 0.8</td><td>5.6 ± 0.8</td></tr><tr><td>Choline</td><td>8.8 ± 1.6</td><td>9.4 ± 0.5</td><td>10.9 ± 0.9</td><td>9.7 ± 0.8</td><td>10.3 ± 1.7</td><td>10.2 ± 0.4</td><td>8.6 ± 0.4</td><td>9.5 ± 1.2</td></tr><tr><td>Creatine</td><td>7.8 ± 1.8</td><td>8.4 ± 0.4</td><td>7.7 ± 0.6</td><td>7.8 ± 0.6</td><td>7.9 ± 0.9</td><td>7.9 ± 0.5</td><td>7.7 ± 0.6</td><td>7.7 ± 0.9</td></tr><tr><td>Glucose-1-phosphate</td><td>1.1 ± 0.6</td><td>0.7 ± 0.5</td><td>0.8 ± 0.5</td><td>3.1 ± 1.5</td><td>9.8 ± 2.1</td><td>39.5 ± 5.7</td><td>87.2 ± 41.8</td><td>117.1 ± 1.4</td></tr><tr><td>Glycine</td><td>25.2 ± 5.0</td><td>31.4 ± 5.7</td><td>30.8 ± 3.1</td><td>25.0 ± 2.1</td><td>39.3 ± 6.4</td><td>26.5 ± 0.9</td><td>27.5 ± 1.8</td><td>25.8 ± 4.2</td></tr><tr><td>Isoleucine</td><td>2.6 ± 0.5</td><td>3.2 ± 0.2</td><td>2.8 ± 0.3</td><td>2.6 ± 0.3</td><td>2.7 ± 0.3</td><td>2.5 ± 0.0</td><td>2.5 ± 0.4</td><td>2.4 ± 0.4</td></tr><tr><td>Lactate</td><td>33.0 ± 6.4</td><td>37.7 ± 2.4</td><td>36.7 ± 3.1</td><td>36.9 ± 2.5</td><td>43.0 ± 4.1</td><td>41.9 ± 2.7</td><td>43.3 ± 6.9</td><td>43.7 ± 6.1</td></tr><tr><td>Maltose</td><td>13.9 ± 5.0</td><td>19.1 ± 5.0</td><td>14.1 ± 1.7</td><td>29.3 ± 10.6</td><td>53.2 ± 4.3</td><td>110.5 ± 13.1</td><td>203.7 ± 9.7</td><td>225.1 ± 32.9</td></tr><tr><td>Methionine</td><td>10.0 ± 1.9</td><td>9.8 ± 0.6</td><td>8.7 ± 1.0</td><td>9.2 ± 1.0</td><td>7.3 ± 0.9</td><td>7.4 ± 0.6</td><td>6.2 ± 0.8</td><td>6.4 ± 0.8</td></tr><tr><td>Myoinositol</td><td>37.6 ± 7.9</td><td>41.5 ± 2.2</td><td>43.8 ± 4.2</td><td>44.4 ± 4.2</td><td>50.9 ± 5.5</td><td>51.5 ± 3.1</td><td>47.8 ± 5.7</td><td>50.4 ± 6.2</td></tr><tr><td>Succinate</td><td>0.7 ± 0.2</td><td>0.7 ± 0.0</td><td>0.9 ± 0.1</td><td>0.8 ± 0.1</td><td>1.1 ± 0.1</td><td>0.9 ± 0.1</td><td>1.0 ± 0.2</td><td>0.8 ± 0.2</td></tr><tr><td>Valine</td><td>3.2 ± 0.8</td><td>3.8 ± 0.2</td><td>3.3 ± 0.4</td><td>3.3 ± 0.4</td><td>3.2 ± 0.5</td><td>2.9 ± 0.0</td><td>2.7 ± 0.2</td><td>2.7 ± 0.2</td></tr></tbody></table>

<footnote n="1">¹ Information detailed on the table was obtained from "The human metabolome database" (http://www.hmdb.ca/).</footnote>
## 2.3. Effect of the Stage of the Estrous Cycle and Side of Ovulation on Metabolite Concentrations in oEVs

Principal component analysis of the metabolite levels showed a clear separation of metabolite profiles between the peri-ovulatory stages (Post-ov and Pre-ov), which clustered together, and the Mid-lut and Late-lut stages (Figure 2). However, the PCA plot could not discriminate between the two sides relative to ovulation (ipsilateral vs. contralateral).

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparative analysis of the oviduct EV metabolite content across the bovine estrous cycle. Principal component analysis (PCA) of metabolites measured in oEV at four different stages of the estrous cycle. S1, post-ovulation (red); S2, mid-luteal phase (blue); S3, late luteal phase (green) and S4, pre-ovulation (purple) from oviducts ipsilateral (triangles) and contralateral (round spots) to ovulation.</figcaption></figure>

The ANOVA confirmed the results of the multivariate analysis. Except for maltose, which was slightly affected by the side of ovulation (p = 0.04), there was no difference between the ipsilateral and contralateral oviducts or cycle stage × side interaction in the concentrations of any metabolite. The stage of the cycle had a significant effect on the glucose-1-phosphate, maltose, methionine, and acetone concentrations (Figure 3). Glucose-1-phosphate was the metabolite most affected by the cycle stage, showing 100-fold higher levels at the luteal phase than at peri-ovulatory phases (p < 0.0001). Similarly, maltose was on average 8- to 14-fold more concentrated at Late-lut than at Post-ov and Pre-ov in both sides relative to ovulation (p < 0.0001). In contrast, concentrations of methionine were significantly higher at Pre-ov and Post-ov than at Late-lut (p < 0.01). In addition, concentrations of acetone were lower at Pre-ov than at other stages of the cycle (p < 0.01).
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Differential concentrations of specific bovine oEV metabolites across the estrus cycle. Maltose intra-oEV concentrations (nmol.mg⁻¹ of EV protein, (a) was affected by the cycle stage and side of ovulation. Glucose-1-P (b), methionine (c) and acetone (d) intra-oEV concentrations were only influenced by the stage of the estrus cycle. For (b–d), ipsilateral and contralateral concentrations data were pooled as the side of ovulation did not show any significant effect on the metabolite level. (Post-ov: post-ovulatory phase; Mid-lut: mid-luteal phase; Late-lut: late luteal phase; Pre-ov: pre-ovulatory phase).</figcaption></figure>

## 2.4. Pathways Associated with oEVs Metabolites

The over-representation analysis (ORA) of all identified oEVs metabolites against the pathway associated metabolite sets of MetaboAnalyst showed that glycine and serine metabolism was the most significant pathway (*p* < 0.001 and false discovery rate (FDR) < 0.1; Table 4 and Figure 4).

The quantitative enrichment analysis using the metabolite concentrations at Pre-ov vs. Late-lut stages is shown in Figure 5. Starch (polysaccharide) and sucrose metabolism, nucleotide sugars metabolism, glycolysis, lactose synthesis, gluconeogenesis, galactose metabolism, spermidine/spermine biosynthesis and betaine synthesis were the most regulated pathways between Pre-ov and Late-lut (*p* < 0.01 and FDR < 0.1; Table 5). Similar metabolic pathways related to polysaccharide, sucrose, glucose, and lactose metabolism were significantly regulated by comparing Pre-ov vs. Mid-lut, Post-ov vs. Mid-lut, and Post-ov vs. Late-lut whereas the Post-ov vs. Pre-ov comparison retrieved no significant results.
<table><caption>**Table 4.** Over representation analysis of all identified metabolites in oEVs against pathway associated metabolite sets. The associated bar chart visualization of this analysis is shown in Figure 4.</caption><thead><tr><th>Pathway Associated Metabolite Sets <sup>1</sup></th><th>Implicated Metabolites <sup>2</sup></th><th>Total <sup>3</sup></th><th>Hits <sup>4</sup></th><th>Raw p <sup>5</sup></th><th>Holm p <sup>6</sup></th><th>FDR <sup>7</sup></th></tr></thead><tbody><tr><td>Glycine and serine metabolism</td><td>Creatine, Glycine, L-glutamic acid, L-alanine, L-threonine, L-methionine</td><td>59</td><td>6</td><td>7.82 x 10<sup>-4</sup></td><td>0.0766</td><td>0.0766</td></tr><tr><td>Arginine and Proline Metabolism</td><td>Creatine, Glycine, L-glutamic acid, L-proline, Succinic acid</td><td>53</td><td>5</td><td>0.00332</td><td>0.322</td><td>0.137</td></tr><tr><td>Alanine metabolism</td><td>Glycine, L-glutamic acid, L-alanine</td><td>17</td><td>3</td><td>0.00421</td><td>0.404</td><td>0.137</td></tr><tr><td>Valine, Leucine and Isoleucine degradation</td><td>L-glutamic acid, L-isoleucine, Succinic acid, L-leucine, L-valine</td><td>60</td><td>5</td><td>0.00575</td><td>0.546</td><td>0.141</td></tr><tr><td>Glutathione Metabolism</td><td>Glycine, L-glutamic acid, L-alanine</td><td>21</td><td>3</td><td>0.0078</td><td>0.734</td><td>0.146</td></tr><tr><td>Carnitine Synthesis</td><td>L-carnitine, Glycine, Succinic acid</td><td>22</td><td>3</td><td>0.00892</td><td>0.829</td><td>0.146</td></tr><tr><td>Glutamate Metabolism</td><td>Glycine, L-glutamic acid, L-alanine, Succinic acid</td><td>49</td><td>4</td><td>0.0152</td><td>1</td><td>0.212</td></tr><tr><td>Phosphatidylethanolamine biosynthesis</td><td>Choline, Ethanolamine</td><td>12</td><td>2</td><td>0.0234</td><td>1</td><td>0.267</td></tr><tr><td>Ketone Body Metabolism</td><td>Succinic acid, Acetone</td><td>13</td><td>2</td><td>0.0273</td><td>1</td><td>0.267</td></tr><tr><td>Glucose-Alanine Cycle</td><td>L-glutamic acid, L-alanine</td><td>13</td><td>2</td><td>0.0273</td><td>1</td><td>0.267</td></tr><tr><td>Phosphatidylcholine biosynthesis</td><td>Choline, Ethanolamine</td><td>14</td><td>2</td><td>0.0314</td><td>1</td><td>0.28</td></tr></tbody></table>

Enrichment Overview (top 50)

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Graphical overview of the over representation analysis (ORA) for all metabolites identified in bovine oEVs generated by MetaboAnalyst 4.0 web-based software. Pathway associated metabolite sets are sorted based on fold enrichment and *p* value. Further details on *p*-values and metabolites included in each pathway are detailed in Table 4.</figcaption></figure>

<footnote n="1">¹ Only pathway associated metabolite sets with *p* < 0.05 are shown in the table; ² Implicated metabolites from oEVs in the associated metabolite set; ³ Total number of metabolites in the metabolite set; ⁴ Hits: number of metabolites from oEVs involved in the metabolite set; ⁵ Raw *p*: original *p* value calculated from the enrichment analysis; ⁶ Holm *p*: adjusted raw *p* value by the Holm-Bonferroni method; ⁷ FDR: false discovery rate.</footnote>
Metabolite Sets Enrichment Overview

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Graphical overview of the quantitative enrichment analysis for all metabolites quantified in bovine oEVs generated by MetaboAnalyst 4.0 web-based software. Pathway associated metabolite sets are sorted based on fold enrichment and p-value after comparison between stages. This bar chart was obtained by comparing Pre-ov vs. Late-lut stages. For each stage, ipsilateral and contralateral data were pooled. Further details on p-values and metabolites included in each pathway are detailed in Table 5.</figcaption></figure>

<table><caption>**Table 5.** Quantitative enrichment analysis of bovine oEV metabolites against pathway associated metabolite sets. Associated bar chart visualization of this analysis is shown in Figure 5.</caption><thead><tr><th>Pathway Associated Metabolite Sets <sup>1</sup></th><th>Implicated Metabolites <sup>2</sup></th><th>Total <sup>3</sup></th><th>Hits <sup>4</sup></th><th>Raw p <sup>5</sup></th><th>Holm p <sup>6</sup></th><th>FDR <sup>7</sup></th></tr></thead><tbody><tr><td>Starch and sucrose metabolism</td><td>D-maltose, Glucose-1-phosphate</td><td>31</td><td>2</td><td>3.21 x 10<sup>-5</sup></td><td>0.0015</td><td>6.835 x 10<sup>-4</sup></td></tr><tr><td>Nucleotide sugars metabolism</td><td>Glucose-1-phosphate</td><td>20</td><td>1</td><td>8.651 x 10<sup>-5</sup></td><td>0.0039</td><td>6.835 x 10<sup>-4</sup></td></tr><tr><td>Glycolysis</td><td>Glucose-1-phosphate</td><td>25</td><td>1</td><td>8.651 x 10<sup>-5</sup></td><td>0.0039</td><td>6.835 x 10<sup>-4</sup></td></tr><tr><td>Lactose synthesis</td><td>Glucose-1-phosphate</td><td>20</td><td>1</td><td>8.651 x 10<sup>-5</sup></td><td>0.0039</td><td>6.835 x 10<sup>-4</sup></td></tr><tr><td>Gluconeogenesis</td><td>Glucose-1-phosphate, L-lactic acid</td><td>35</td><td>2</td><td>8.662 x 10<sup>-5</sup></td><td>0.0039</td><td>6.835 x 10<sup>-4</sup></td></tr><tr><td>Galactose metabolism</td><td>Myoinositol, Glucose-1-phosphate</td><td>38</td><td>2</td><td>8.726 x 10<sup>-5</sup></td><td>0.1933</td><td>6.835 x 10<sup>-4</sup></td></tr><tr><td>Spermidine and spermine biosynthesis</td><td>L-methionine</td><td>18</td><td>1</td><td>0.0047</td><td>0.4906</td><td>0.032</td></tr><tr><td>Betaine metabolism</td><td>Choline, L-methionine</td><td>21</td><td>2</td><td>0.0123</td><td>1</td><td>0.072</td></tr><tr><td>Methionine metabolism</td><td>Choline, Glycine, L-methionine</td><td>43</td><td>3</td><td>0.0438</td><td>1</td><td>1</td></tr></tbody></table>

<footnote n="1">¹ Only pathway associated metabolite sets with *p* < 0.05 are shown in the table; ² Implicated metabolites from oEV in the associated metabolite set; ³ Total number of metabolites in the metabolite set; ⁴ Hits: number of metabolites from oEV involved in the metabolite set; ⁵ Raw *p*: original *p* value calculated from the enrichment analysis; ⁶ Holm *p*: adjusted raw *p* value by the Holm-Bonferroni method; ⁷ FDR: false discovery rate.</footnote>
### 3. Discussion

Currently, there is an increasing number of studies analyzing the metabolomics of EVs in different biological fluids, since both metabolites and EVs are being considered as valuable biomarkers of health conditions and play a crucial part in cell-to-cell communications. In the context of pregnancy, oEVs have been proposed as key modulators of the gametes/embryo-oviduct crosstalk to achieve successful pregnancy [12]. To the best of our knowledge, this study provides the first metabolic characterization of oEVs. Furthermore, we demonstrated that the oEV-metabolite cargo is regulated by the hormonal environment of the estrus cycle, although to a lesser extent than proteins and mRNA in oEVs [5,9].

In the present study, NMR allowed the identification of 22 metabolites in oEVs obtained from cows at four different stages of the estrus cycle. Among them, various classes of metabolites were found such as amino acids, vitamins, energy substrates, and others presented in Table 2. Although the amount of detected metabolites in oEVs was low compared to other metabolomic studies performed by mass spectrometry methods, NMR was selected because of (1) its advantages in highly selective, quantitative, and non-destructive analysis [16], and (2) it was the method used to analyze the metabolite composition of the OF in a previous study performed in our laboratory [15], allowing us a better comparison between data. In fact, most of the metabolites identified here were previously identified and quantified in the bovine OF by NMR [15]. Moreover, the predominant metabolites in oEVs (myoinositol, lactate, and glycine), were also the most abundant metabolites quantified by NMR in the OF by [15]. However, the less abundant metabolites in the OF (i.e., aspartate, choline, mannose, and pyruvate) were not detected in the present study. We can speculate that solute metabolites present in the OF cross the oEV membranes and reach intra-oEV concentrations that mirror those in the OF. Minor metabolites present in the fluid could remain at undetectable levels after being packed into oEVs. Nevertheless, we cannot discard the idea that the methodology used in the present study or the starting oEV volume used limited the detection of those minor metabolites. On the other hand, we identified metabolites such as maltose and acetone, which had not been previously detected in bovine OF [15]. These results suggest that these metabolites may be specifically transported from the tubal fluid into oEVs, reaching higher concentrations inside oEVs than outside. However, the mechanisms by which these metabolites concentrated into oEVs still remains unknown.

In the present study, besides determining whether the hormonal environment of the estrus cycle controls the metabolic profile of oEVs, we also aimed to determine if the side related to ovulation could affect the oEV-metabolite composition, by comparing profiles in ipsilateral versus contralateral oviducts. We thought it was important to include this factor because (1) in previous work from our laboratory, we used pools of ipsilateral and contralateral oviducts [5,9], while others have focused only on the ipsilateral side [17], ignoring the evaluation of the side effect; (2) significant differences in steroid hormones and protein content in the OF were recorded between the two sides relative to ovulation in cattle [18,19]. Moreover, the cow, as a mono-ovular species, is an ideal model to study the effect of the proximity of the ovulatory follicle and post-ovulatory corpus luteum on oviduct physiology. Nevertheless, our results showed that the abundance of metabolites in oEVs was overall not affected by the side relative to ovulation. Only the concentration of maltose was slightly different between the ipsilateral and contralateral oviducts. In the same line, only five metabolites out of the 26 quantified by NMR in the OF were affected by the side of ovulation, with small differences between ovulation sides [15]. Other studies in the field of metabolomics have failed to evidence significant differences in the concentrations of lipids, amino acids, and energy substrates between ipsilateral and contralateral oviducts [20–24]. The contribution of the pre-ovulatory follicle to the content of the OF after ovulation has been under debate for many years [14]. However, the similar metabolite composition between both ovulation sides at Post-ov in oEVs as well as in OF [15], together with other studies above-mentioned does not support this hypothesis.

As spermatozoa and preimplantation embryos are able to uptake oEVs [5,6,8], oEV-metabolites at the pre- and post-ovulatory stages may contribute to their continuous supply of amino acids and energy substrates for the success of early reproductive events. In this sense, we would like to point out
that the three most abundant metabolites identified in oEVs (i.e., myoinositol, lactate, and glycine, may have a higher impact on gametes and embryos. Myoinositol and lactate were found to be major energy substrates in oEVs throughout the estrous cycle. Myoinositol is synthesized from glucose and is an energy source for both spermatozoa [25] and embryos [26]. Myoinositol improved bovine sperm motility [27] and embryo development up to the blastocyst stage when added in vitro [28]. Lactate, together with pyruvate, have been identified as the main energy supply for sperm capacitation [29]. After fertilization, the early pre-implantation embryo uses oxidative metabolism to acquire energy, mostly from lactate and pyruvate as the main sources of energy [30]. Glycine, a conditionally essential amino acid, was the third most abundant metabolite in our study and was by far the most abundant amino acid in the oEVs. Using the ORA program of MetaboAnalyst, several amino acid-associated pathways, primarily serine and glycine metabolism, were identified as over-represented in oEVs. Data on the metabolic pathways associated with the establishment of gestation are scarce. However, glycine is known to be important for early embryo development as it is an essential precursor for the synthesis of proteins and nucleic acids [31,32]. Moreover, glycine is required for the rapid proliferation of cells [33]. Glycine also regulates intracellular pH and may have a protecting role against osmotic stress for preimplantation embryos [34].

Among the metabolites identified, glucose-1-P, maltose, methionine, and acetone displayed differential concentrations across the estrus cycle, which may reflect the requirements of gamete/embryos in the oviduct environment. Maltose and glucose-1-phospate, energy substrates in oEVs, were by far the metabolites most affected by the cycle stage, with 8- to 100-fold higher levels at the luteal phase than at the pre- and post-ovulatory stages. Most pathways identified in the quantitative enrichment analysis including starch and sucrose metabolism, nucleotide sugars metabolism, glycolysis, lactose synthesis, and neoglucogenesis involved one or both of these carbohydrates. The sugar D-maltose is a glucoside consisting of two glucose monomers connected by an α-1,4-glycosidic bond. Maltose is usually the product of starch digestion, after alpha-amylase enzymatic action [35]. Although starch is not present in the oviduct, a large amount of glycogen was reported in the mammalian oviduct, where it acts as a glucose mobilized reservoir, storing up to 55,000 glucose moieties per molecule [36]. In the bovine oviduct epithelium, glycogen was quantified at around 60 µg/mg protein [37]. Most of the glucose residues in glycogen macromolecules are linked by α-1,4-glycosidic bonds while about one tenth of glucose residues are linked by α-1,6-glycosidic bonds. The enzymes alpha-amylase (AMY2B) and glycogen phosphorylase (PYGL, liver form and PYGB, brain form), all able to cut α-1,4-linkages, are produced in oviduct epithelial cells in mammals [30,36] and were reported at high abundance in the bovine OF [18,38]. Thus, it is likely that maltose and glucose-1-phosphate are products of glycogen catabolism in the oviduct epithelium. The abundance of maltose and glucose-1-phosphate in oEVs was upregulated during the luteal phase and reached its highest levels at the Late-lut stage. Similarly, glucose-1-phosphate in the bovine OF was reported to be 6-fold more abundant during the luteal phase than at the peri-ovulatory stages [15]. It is of note that levels of maltose and glucose-1-phosphate in oEVs paralleled those of P4 and were in opposition with those of E2 as previously reported in the bovine OF [19] and circulating plasma [39]. The hormonal regulation of glycogen metabolism in the oviduct remains unclear, but reproductive cycle-dependent changes in oviductal glycogen have been evidenced [36]. In pigs and rabbits, reserves of glycogen in the oviducts were highest near ovulation and decreased during the following luteal phase [40,41]. Furthermore, the hormonal control of endometrial glycogen has been well studied in rats, rabbits, and mink, where P4 promotes glycogen catabolism while E2 stimulates uterine glycogen synthesis and storage [36,42,43]. Altogether, these data suggest that the products of glycogen catabolism in the oviduct epithelium are abundantly released in oEVs under P4 action during the luteal phase and that this process is inhibited by E2 during estrus. The low levels of maltose and glucose-1-phosphate in oEVs during the peri-ovulatory period may play important roles for the success of reproduction. Interestingly, maltose has been shown to inhibit the binding of bovine and porcine spermatozoa to oviduct epithelial cell carbohydrates, a process allowing sperm storage and the delivery of capacitated spermatozoa to the site of fertilization [44,45]. The low
levels of sources of glucose in oEVs in the post-ovulatory period are also concordant with the fact that exposure of the early bovine embryo to high amounts of glucose in vitro disturb embryo development and metabolism [46,47].

Additionally, we found that methionine concentrations in oEVs were higher at Pre-ov and Post-ov than at Late-lut. Methionine is a precursor for polypeptide synthesis and is an important factor in the antioxidant balance [48]. Furthermore, due to its role as the immediate precursor of S-adenosyl methionine, which is the universal methyl donor for epigenetic methylation of DNA and histones, methionine can influence the epigenetic programming of the embryo [49]. There is evidence that preimplantation embryo development is vulnerable to disruptions in methionine metabolism. For instance, disruption of methionine metabolism by ethionine, an antimetabolite of methionine, impairs the morula-to-blastocyst transition during bovine preimplantation development in vitro [50]. Supplementation of an in vitro embryo culture medium with large doses of S-adenosyl methionine was reported to cause a shift in bovine embryo sex ratio in favor of males and to induce a genome-wide hypermethylation [51].

Finally, we would like to discuss the unexpected presence of acetone, although at very low levels at all stages examined (≤1 nmol/mg of oEV protein), in oEVs. In dairy cows with subclinical or clinical ketosis, acetone concentrations in the circulating blood may reach 300–350 µmol/L [52,53]. Ketosis (or acetonemia) is a common disease in high-producing dairy cows caused by a negative energy balance and occurring mostly within two months after calving. Concentrations of acetone, acetoacetate, and β-hydroxybutyrate, three ketone bodies, increase in blood and milk during ketosis [52,53]. This metabolic picture is typically associated with anestrus and lower fertility in dairy cows [54]. The ketosis status of the cows used for the present study could not been assessed at the time of oviduct collection, but they were mostly Holstein dairy cows, commonly regarded as particularly subject to ketosis. Thus, it is likely that the low concentrations of acetone detected in oEVs reflected low to moderate levels of acetone in the blood of animals at the time of slaughtering. However, why acetone levels in oEVs were lower at Pre-ov than at other stages of the estrous cycle is difficult to explain, with no data available in the current literature to shed some light.

In conclusion, our study provides the first comprehensive metabolomic profile of oEVs. Moreover, our results show that the hormonal environment of the estrus cycle affects the oEV-metabolite composition, with significant differences between the peri-ovulatory stages and late-luteal phase for specific metabolites. Our findings bring new insights into the roles of oEVs in supporting sperm function, fertilization, and early embryonic development.

# 4. Materials and Methods

## 4.1. *Collection of Oviductal Fluids*

Bovine oviducts were collected at a local slaughterhouse, transported on ice, and classified into four stages of the estrous cycle based on the ovarian and corpus luteum morphologies, as previously described [19]: post-ovulatory (Post-ov; Days 1–4 of estrous cycle, coinciding with the time of embryo presence in the oviduct), mid luteal (Mid-lut; Days 5–11), late luteal (Late-lut; Days 11–17), and pre-ovulatory (Pre-ov; Days 18–20, coinciding, if inseminated, with the time of sperm capacitation in the oviduct) (n = 34–54 cows/stage). The OF and epithelial cells were collected from the ipsilateral (to corpus luteum or Pre-ov follicle) and contralateral oviducts by gentle squeezing with a glass slide. Mixtures were kept on ice until OF was separated from cells by centrifugation at 2000× g for 15 min. Then, OF was stored at –80 °C until oEV isolation. To avoid the inclusion of cows with cystic follicles in the Pre-ov group, animals with a Pre-ov follicle larger than 20 mm in diameter were excluded at the time of oviduct collection. Furthermore, follicular fluid was collected from pre-ovulatory follicles and assayed for P4 and E2 concentrations by a competitive enzyme-linked immunosorbent assay [55] and the BioSource E2-EASIA Kit (BioSource, Louvain-la-Neuve, Belgium), respectively. Cows with
intra-follicular concentrations of P4 higher than 160 ng/mL and of E2 lower than 40 ng/mL were also excluded, as previously described [19].

## 4.2. Isolation and Characterization of oEVs

### 4.2.1. Preparation of Pools Oviductal Fluids and Isolation of oEVs

Preliminary experiments showed that 280 µL of OF was the minimal volume required for the detection of metabolites in oEVs by NMR. In order to reach at least three samples per estrous cycle stage and side, and based on the average volume of 30 µL of OF collected per oviduct, a minimum of 34 animals per stage were used (Table 1). To isolate oEVs, OF samples were thawed on ice and 3-4 biological pools were prepared by mixing OF from individual OF samples to reach a final minimum volume of 280 µL/pool (Table 1). Pools of OF were centrifuged at 12,000× g for 15 min to remove cellular debris and then ultracentrifuged twice at 100,000 g for 90 min to pellet the oEVs, as previously described [5]. The oEVs pellets were resuspended in 50 µL PBS, assayed for protein concentration by BCA (Bicinchoninic Acid) Protein Assay (ThermoFisher, Nanterre, France) and stored at -80 °C for further analysis.

### 4.2.2. Oviductal EVs Observation by TEM

Aliquots of oEVs samples (one per stage and per side) were fixed with glutaraldehyde 2% for observation under TEM, as previously described [9]. The micrographs were obtained using TEM HITACHI HT 7700 Elexience at 80 kV (with a charge-coupled device camera AMT) and JEM 1011 (JEOL, Tokyo, Japan) equipped with a Gatan digital camera driven by Digital Micrograph software (Gatan, Pleasanton, CA, USA) at 100 kV. The processing of the photos and vesicle size calculation (min 300 vesicles per stage) were carried out by ImageJ software.

### 4.2.3. Oviductal EVs Characterization by Western Blotting

Proteins (30 µg) were separated by SDS-PAGE and transferred onto nitrocellulose membranes (Whatman™, GE Healthcare Life Sciences, Amersham, UK). The membranes were blocked with Tris-buffered saline (TBS) containing Tween 20 (0.5% (w/v)) and supplemented with lyophilized low-fat milk (5% w/v) for 3 h at room temperature (TBS-T). The membranes were incubated with primary antibodies diluted in TBS-T containing 5% milk at 4 °C with gentle shaking overnight. The primary antibodies used were directed against CD81 (1:200 dilution; sc-166029, Santa Cruz Biotechnologies, Heidelberg, Germany), heat shock protein 70 (HSP70, 1:1000; Enzo, ADI-SPA-810), and annexin A1 (ANXA1, 1:400; Santa Cruz Biotechnologies, sc-114387). Subsequently, the membranes were washed with TBS-T and incubated at room temperature under agitation with secondary antibodies. The secondary antibodies used were Goat anti Mouse-HRP (1:5000; A4416, Merck, Saint-Quentin-Fallavier, France) for CD81 and HSP70 and Goat anti Rabbit-HRP (1:5000; A6154, Merck) for ANXA1. Blots were developed using a mixture of two chemiluminescence substrate developing kit (GE Healthcare AmershamTH ECL SelectTH Western blotting detection Reagent RPN2235 and Supersignal West Pico #34087 Chemiluminescent Substrate Thermo Scientific).

## 4.3. Proton Nuclear Magnetic Resonance Spectroscopy (1H-NMR) Analyses

Samples of oEVs were analyzed by NMR as previously described for bovine OF samples [15]. Briefly, 24 µL of oEV samples were diluted in 172 µL of 0.2 M potassium phosphate buffer in deuterium oxide (D2O) (pH 7.4 ± 0.5) and 4 µL of 3-(trimethylsilyl) propionic-acid solution at 3.2 mM, used as a spectrometer field lock signal. The resulting solution was transferred to conventional 3-mm NMR tubes. Briefly, 1H-NMR spectra were recorded at 298 K on a Bruker Ascend 600 MHz spectrometer (Bruker, Sadis, Wissembourg, France), equipped with a TCI cryoprobe. Standard NMR spectra were obtained using a "cpmgpr1d" pulse sequence with a 90° pulse, a relaxation delay of 25 s, number of scans of 256 on a time domain of 64 K data points. Data were processed with 0.2 Hz of line broadening for the
exponential decay function using TopSpin version 3.2 software (Bruker daltonik, Karlshure, Germany). Metabolite quantifications using the ERETIC peak as a quantitative reference were obtained by the specific subroutine of the Bruker TopSpin 3.2 program. The ERETIC signal was calibrated on a known concentration amino-acid solution ([ERETIC] = 198 µM). Spectral assignments were done using the free version of ChenomX 7.1 software (ChenomX, Edmonton, Canada), in house database, and human metabolome database (HMDB version 4.0; http://www.hmdb.ca, assessed 18 March 2019) [56]. When the assignment was uncertain, spiking experiments were used to validate metabolite assignments. Comparing the NMR spectra obtained from a biological sample spiked-in with a known metabolite with that obtained without spiking made it possible to confirm or discard the presence or absence of this metabolite in the biological sample.

## 4.4. Data Processing and Statistical Analysis of NMR Data

First, the concentration of each metabolite was normalized to 1 mg of oEV protein. Then, to obtain an overview of the oEV metabolite levels across the different stages of the estrous cycle, principal component analysis (PCA) was generated using normalized metabolite data in RStudio (R packages: ggplot2, ggfortify). Statistical analyses were performed using the GraphPad Prism software (version 7, San Diego, CA, USA). Normalized values were compared between stages of the cycle and sides relative to ovulation with a two-way analysis of variance (ANOVA) followed by Tukey post-tests. When no side × stage and side effect (*p* > 0.05) were identified, data for ipsilateral and contralateral oviducts were combined and an average value calculated for each stage. A *p*-value < 0.05 was considered as significant. Data are presented as means ± SEM.

## 4.5. Pathway Enrichment Analyses of NMR Data

To investigate the presence of biologically meaningful patterns among the oEV metabolites, metabolic pathway analysis was performed using the functional enrichment tool of the free web-based software MetaboAnalyst 4.0 (https://www.metaboanalyst.ca/; assessed 12–24 June 2019) [57]. First, enrichment over representation analysis (ORA) was performed using the list of all metabolites identified in oEVs. Then, to obtain further insights on the regulated metabolic pathways across the cycle, a quantitative enrichment analysis by paired comparisons between stages using normalized (log transformed) concentrations of quantified metabolites was performed. The *p*-value and the false discovery rate (FDR) of each regulated pathway were used to account for significance and false discovery in the pathway analysis, respectively. In this regard, the thresholds of *p*-value < 0.05 and FDR < 0.1 were considered as significant.

**Author Contributions:** Conceptualization; M.S.-D., C.A., and J.G.; Methodology; J.G., D.L.B., C.A., S.J.I., O.B., G.T., and R.U.; Software; J.G., O.B., and M.S.-D.; Formal analysis; J.G., O.B., C.A., and M.S.-D.; Writing—original draft preparation; M.S.-D. and C.A.; Writing—review and editing; all authors.; Supervision; M.S.-D. and P.S.; Project administration; M.S.-D. and P.M.; Funding acquisition; M.S.-D. and P.S.

**Funding:** This research was funded by APIS-GENE under the project name "Embryomimétisme".

**Acknowledgments:** The authors are grateful to Marc Chodkiewicz for his careful editing of this paper and would like to thank Olivier Desnoes, Thierry Delpuech, and Albert Arnould for their help in collecting the genital tracts at the slaughterhouse. Our data were obtained with the assistance of the IBISA Electron Microscopy Facility of the University of Tours and the University Hospital of Tours. M.S.D. thanks the CNRS for hosting her on delegation. M.S.D. thanks the Centre National de la Recherche Scientifique for hosting her on delegation.

**Conflicts of Interest:** The authors declare no conflicts of interest.

## Abbreviations

NMR Nuclear magnetic resonance

oEVs Oviductal extracellular vesicles

OF Oviductal fluid
## References

1. Yanez-Mo, M.; Siljander, P.R.; Andreu, Z.; Zavec, A.B.; Borras, F.E.; Buzas, E.I.; Buzas, K.; Casal, E.; Cappello, F.; Carvalho, J.; et al. Biological properties of extracellular vesicles and their physiological functions. *J. Extracell. Vesicles* **2015**, *4*, 27066. [CrossRef] [PubMed]

2. da Silveira, J.C.; Veeramachaneni, D.N.; Winger, Q.A.; Carnevale, E.M.; Bouma, G.J. Cell-secreted vesicles in equine ovarian follicular fluid contain miRNAs and proteins: A possible new form of cell communication within the ovarian follicle. *Biol. Reprod.* **2012**, *86*, 71. [CrossRef] [PubMed]

3. Ng, Y.H.; Rome, S.; Jalabert, A.; Forterre, A.; Singh, H.; Hincks, C.L.; Salamonsen, L.A. Endometrial exosomes/microvesicles in the uterine microenvironment: A new paradigm for embryo-endometrial cross talk at implantation. *PLoS ONE* **2013**, *8*, e58502. [CrossRef] [PubMed]

4. Lopera-Vasquez, R.; Hamdi, M.; Maillo, V.; Gutierrez-Adan, A.; Bermejo-Alvarez, P.; Ramirez, M.A.; Yanez-Mo, M.; Rizos, D. Effect of bovine oviductal extracellular vesicles on embryo development and quality in vitro. *Reproduction* **2017**, *153*, 461–470. [CrossRef] [PubMed]

5. Alminana, C.; Corbin, E.; Tsikis, G.; Alcantara-Neto, A.S.; Labas, V.; Reynaud, K.; Galio, L.; Uzbekov, R.; Garanina, A.S.; Druart, X.; et al. Oviduct extracellular vesicles protein content and their role during oviduct-embryo cross-talk. *Reproduction* **2017**, *154*, 153–168. [CrossRef] [PubMed]

6. Al-Dossary, A.A.; Strehler, E.E.; Martin-Deleon, P.A. Expression and secretion of plasma membrane Ca²⁺-ATPase 4a (PMCA4a) during murine estrus: Association with oviductal exosomes and uptake in sperm. *PLoS ONE* **2013**, *8*, e80181. [CrossRef] [PubMed]

7. Lopera-Vasquez, R.; Hamdi, M.; Fernandez-Fuertes, B.; Maillo, V.; Beltran-Brena, P.; Calle, A.; Redruello, A.; Lopez-Martin, S.; Gutierrez-Adan, A.; Yanez-Mo, M.; et al. Extracellular Vesicles from BOEC in In Vitro Embryo Development and Quality. *PLoS ONE* **2016**, *11*, e0148083. [CrossRef]

8. Ferraz, M.; Carothers, A.; Dahal, R.; Noonan, M.J.; Songsasen, N. 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.* **2019**, *9*, 9484. [CrossRef]

9. Alminana, C.; Tsikis, G.; Labas, V.; Uzbekov, R.; da Silveira, J.C.; Bauersachs, S.; Mermillod, P. Deciphering the oviductal extracellular vesicles content across the estrous cycle: Implications for the gametes-oviduct interactions and the environment of the potential embryo. *BMC Genom.* **2018**, *19*, 622. [CrossRef]

10. Alcântara-Neto, A.; Fernandez-Rufete, M.; Corbin, E.; Tsikis, G.; Uzbekov, R.; Garanina, A.S.; Coy, P.; Alminana, C.; Mermillod, P. Oviduct fluid extracellular vesicles regulate polyspermy during porcine in vitro fertilization. *Reprod. Fertil. Dev.* **2019**. [CrossRef]

11. Qu, P.; Zhao, Y.; Wang, R.; Zhang, Y.; Li, L.; Fan, J.; Liu, E. Extracellular vesicles derived from donor oviduct fluid improved birth rates after embryo transfer in mice. *Reprod. Fertil. Dev.* **2019**, *31*, 324–332. [CrossRef] [PubMed]

12. Alminana, C.; Bauersachs, S. Extracellular Vesicles in the Oviduct: Progress, Challenges and Implications for the Reproductive Success. *Bioeng* **2019**, *6*, 32. [CrossRef]

13. Zebrowska, A.; Skowronek, A.; Wojakowska, A.; Widlak, P.; Pietrowska, M. Metabolome of Exosomes: Focus on Vesicles Released by Cancer Cells and Present in Human Body Fluids. *Int. J. Mol. Sci.* **2019**, *20*, 3461. [CrossRef] [PubMed]

14. Hunter, R.H. Components of oviduct physiology in eutherian mammals. *Biol. Rev. Camb. Philos. Soc.* **2012**, *287*, 244–255. [CrossRef] [PubMed]

15. Lamy, J.; Gatien, J.; Dubuisson, F.; Nadal-Desbarats, L.; Salvetti, P.; Mermillod, P.; Saint-Dizier, M. Metabolomic profiling of bovine oviductal fluid across the oestrous cycle using proton nuclear magnetic resonance spectroscopy. *Reprod. Fertil. Dev.* **2018**, *30*, 1021–1028. [CrossRef]

16. Palomo, L.; Casal, E.; Royo, F.; Cabrera, D.; van-Liempd, S.; Falcon-Perez, J.M. Considerations for applying metabolomics to the analysis of extracellular vesicles. *Front. Immunol.* **2014**, *5*, 651. [CrossRef] [PubMed]

17. Burns, G.; Brooks, K.; Wildung, M.; Navakanitworakul, R.; Christenson, L.K.; Spencer, T.E. Extracellular vesicles in luminal fluid of the ovine uterus. *PLoS ONE* **2014**, *9*, e90913. [CrossRef] [PubMed]

18. Lamy, J.; Labas, V.; Harichaux, G.; Tsikis, G.; Mermillod, P.; Saint-Dizier, M. Regulation of the bovine oviductal fluid proteome. *Reproduction* **2016**, *152*, 629–644. [CrossRef]

19. Lamy, J.; Liere, P.; Pianos, A.; Aprahamian, F.; Mermillod, P.; Saint-Dizier, M. Steroid hormones in bovine oviductal fluid during the estrous cycle. *Theriogenology* **2016**, *86*, 1409–1420. [CrossRef]
20. Hugentobler, S.A.; Diskin, M.G.; Leese, H.J.; Humpherson, P.G.; Watson, T.; Sreenan, J.M.; Morris, D.G. Amino acids in oviduct and uterine fluid and blood plasma during the estrous cycle in the bovine. *Mol. Reprod. Dev.* **2007**, *74*, 445–454. [CrossRef]

21. Hugentobler, S.A.; Humpherson, P.G.; Leese, H.J.; Sreenan, J.M.; Morris, D.G. Energy substrates in bovine oviduct and uterine fluid and blood plasma during the oestrous cycle. *Mol. Reprod. Dev.* **2008**, *75*, 496–503. [CrossRef] [PubMed]

22. Kenny, D.A.; Humpherson, P.G.; Leese, H.J.; Morris, D.G.; Tomos, A.D.; Diskin, M.G.; Sreenan, J.M. Effect of elevated systemic concentrations of ammonia and urea on the metabolite and ionic composition of oviductal fluid in cattle. *Biol. Reprod.* **2002**, *66*, 1797–1804. [CrossRef] [PubMed]

23. Banliat, C.; Tomas, D.; Teixeira-Gomes, A.P.; Uzbekova, S.; Guyonnet, B.; Labas, V.; Saint-Dizier, M. Stage-dependent changes in oviductal phospholipid profiles throughout the estrous cycle in cattle. *Theriogenology* **2019**, *135*, 65–72. [CrossRef] [PubMed]

24. Elhassan, Y.M.; Wu, G.; Leanez, A.C.; Tasca, R.J.; Watson, A.J.; Westhusin, M.E. Amino acid concentrations in fluids from the bovine oviduct and uterus and in KSOM-based culture media. *Theriogenology* **2001**, *55*, 1907–1918. [CrossRef]

25. Condorelli, R.A.; La Vignera, S.; Bellanca, S.; Vicari, E.; Calogero, A.E. Myoinositol: Does it improve sperm mitochondrial function and sperm motility? *Urology* **2012**, *79*, 1290–1295. [CrossRef]

26. Hynes, A.C.; Sreenan, J.M.; Kane, M.T. Uptake and incorporation of myo-inositol by bovine preimplantation embryos from two-cell to early blastocyst stages. *Mol. Reprod. Dev.* **2000**, *55*, 265–269. [CrossRef]

27. Boni, R.; Gallo, A.; Cecchini, S. Kinetic activity, membrane mitochondrial potential, lipid peroxidation, intracellular pH and calcium of frozen/thawed bovine spermatozoa treated with metabolic enhancers. *Andrology* **2017**, *5*, 133–145. [CrossRef]

28. Lim, K.T.; Jang, G.; Ko, K.H.; Lee, W.W.; Park, H.J.; Kim, J.J.; Lee, S.H.; Hwang, W.S.; Lee, B.C.; Kang, S.K. Improved in vitro bovine embryo development and increased efficiency in producing viable calves using defined media. *Theriogenology* **2007**, *67*, 293–302. [CrossRef]

29. Ferramosca, A.; Zara, V. Bioenergetics of mammalian sperm capacitation. *Biomed. Res. Int.* **2014**, *2014*, 902953. [CrossRef]

30. Li, S.; Winuthayanon, W. Oviduct: Roles in fertilization and early embryo development. *J. Endocrinol.* **2017**, *232*, R1–R26. [CrossRef]

31. Redel, B.K.; Spate, L.D.; Lee, K.; Mao, J.; Whitworth, K.M.; Prather, R.S. Glycine supplementation in vitro enhances porcine preimplantation embryo cell number and decreases apoptosis but does not lead to live births. *Mol. Reprod. Dev.* **2016**, *83*, 246–258. [CrossRef] [PubMed]

32. Herrick, J.R.; Lyons, S.M.; Greene, A.F.; Broeckling, C.D.; Schoolcraft, W.B.; Krisher, R.L. Direct and Osmolarity-Dependent Effects of Glycine on Preimplantation Bovine Embryos. *PLoS ONE* **2016**, *11*, e0159581. [CrossRef] [PubMed]

33. Jain, M.; Nilsson, R.; Sharma, S.; Madhusudhan, N.; Kitami, T.; Souza, A.L.; Kafri, R.; Kirschner, M.W.; Clish, C.B.; Mootha, V.K. Metabolite profiling identifies a key role for glycine in rapid cancer cell proliferation. *Science* **2012**, *336*, 1040–1044. [CrossRef] [PubMed]

34. Van Winkle, L.J.; Haghighat, N.; Campione, A.L. Glycine protects preimplantation mouse conceptuses from a detrimental effect on development of the inorganic ions in oviductal fluid. *J. Exp. Zool.* **1990**, *253*, 215–219. [CrossRef] [PubMed]

35. Crow, R.R.; Kumar, S.; Varela, M.F. Maltose Chemistry and Biochemistry. In *Dietery Sugars: Analysis, Function and Effects*; Preedy, V.R., Ed.; Royal Society of Chemistry: London, UK, 2012.

36. Dean, M. Glycogen in the uterus and fallopian tubes is an important source of glucose during early pregnancy. *Biol. Reprod.* **2019**, *101*, 297–305. [CrossRef] [PubMed]

37. Gomes, H.; Dias, A.J.B.; Moraes, J.; Sobrinho Paes de Carvalho, C.; Logullo, C. Glucose-6-phosphate Metabolic Preferential Destinations in Bovine Oviduct Cells. *Acta Sci. Veternariae* **2010**, *38*, 377–383.

38. 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* **2017**, *12*, e0188105. [CrossRef]

39. Glencross, R.G.; Munro, I.B.; Senior, B.E.; Pope, G.S. Concentrations of oestradiol-17 beta, oestrone and progesterone in jugular venous plasma of cows during the oestrous cycle and in early pregnancy. *Acta Endocrinol (Copenh)* **1973**, *73*, 374–384. [CrossRef]
40. Gregoraszczuk, E.L.; Cala, M.; Witkowska, E. Glycogen distribution in porcine fallopian tube epithelium during the estrus cycle. *Folia. Biol.* **2000**, *48*, 85–90.

41. Gregoire, A.T.; Hafs, H.D. Glycogen content of rabbit tubular genitalia from mating through implantation. *Anat. Rec.* **1971**, *169*, 253–258. [CrossRef]

42. Bowman, K.; Rose, J. Estradiol stimulates glycogen synthesis whereas progesterone promotes glycogen catabolism in the uterus of the American mink (Neovison vison). *Anim. Sci. J.* **2017**, *88*, 45–54. [CrossRef] [PubMed]

43. Demers, L.M.; Jacobs, R.D.; Greep, R.O. Comparative effects of ovarian steroids on glycogen metabolism of rat, rabbit and guinea pig uterine tissue. *Proc. Soc. Exp. Biol. Med.* **1973**, *143*, 1158–1163. [CrossRef] [PubMed]

44. Green, C.E.; Bredl, J.; Holt, W.V.; Watson, P.F.; Fazeli, A. Carbohydrate mediation of boar sperm binding to oviductal epithelial cells in vitro. *Reproduction* **2001**, *122*, 305–315. [CrossRef] [PubMed]

45. Sostaric, E.; van de Lest, C.H.; Colenbrander, B.; Gadella, B.M. Dynamics of carbohydrate affinities at the cell surface of capacitating bovine sperm cells. *Biol. Reprod* **2005**, *72*, 346–357. [CrossRef]

46. Thompson, J.G.; Partridge, R.J.; Houghton, F.D.; Cox, C.I.; Leese, H.J. Oxygen uptake and carbohydrate metabolism by in vitro derived bovine embryos. *J. Reprod Fertil.* **1996**, *106*, 299–306. [CrossRef]

47. Uhde, K.; van Tol, H.T.A.; Stout, T.A.E.; Roelen, B.A.J. Exposure to elevated glucose concentrations alters the metabolomic profile of bovine blastocysts. *PLoS ONE* **2018**, *13*, e0199310. [CrossRef]

48. Martinez, Y.; Li, X.; Liu, G.; Bin, P.; Yan, W.; Mas, D.; Valdivie, M.; Hu, C.A.; Ren, W.; Yin, Y. The role of methionine on metabolism, oxidative stress, and diseases. *Amino Acids* **2017**, *49*, 2091–2098. [CrossRef]

49. Zhang, N. Role of methionine on epigenetic modification of DNA methylation and gene expression in animals. *Anim. Nutr.* **2018**, *4*, 11–16. [CrossRef]

50. Ikeda, S.; Sugimoto, M.; Kume, S. Importance of methionine metabolism in morula-to-blastocyst transition in bovine preimplantation embryos. *J. Reprod. Dev.* **2012**, *58*, 91–97. [CrossRef]

51. Shojaei Saadi, H.A.; Gagne, D.; Fournier, E.; Baldoceda Baldeon, L.M.; Sirard, M.A.; Robert, C. Responses of bovine early embryos to S-adenosyl methionine supplementation in culture. *Epigenomics* **2016**, *8*, 1039–1060. [CrossRef]

52. Enjalbert, F.; Nicot, M.C.; Bavourthe, C.; Moncoulon, R. Ketone bodies in milk and blood of dairy cows: Relationship between concentrations and utilization for detection of subclinical ketosis. *J. Dairy Sci.* **2001**, *84*, 583–589. [CrossRef]

53. Reist, M.; Koller, A.; Busato, A.; Kupfer, U.; Blum, J.W. First ovulation and ketone body status in the early postpartum period of dairy cows. *Theriogenology* **2000**, *54*, 685–701. [CrossRef]

54. Roche, J.F. The effect of nutritional management of the dairy cow on reproductive efficiency. *Anim. Reprod. Sci.* **2006**, *96*, 282–296. [CrossRef] [PubMed]

55. Rico, C.; Medigue, C.; Fabre, S.; Jarrier, P.; Bontoux, M.; Clement, F.; Monniaux, D. Regulation of anti-Mullerian hormone production in the cow: A multiscale study at endocrine, ovarian, follicular, and granulosa cell levels. *Biol. Reprod.* **2011**, *84*, 560–571. [CrossRef]

56. Wishart, D.S.; Feunang, Y.D.; Marcu, A.; Guo, A.C.; Liang, K.; Vazquez-Fresno, R.; Sajed, T.; Johnson, D.; Li, C.; Karu, N.; et al. HMDB 4.0: The human metabolome database for 2018. *Nucleic Acids Res.* **2018**, *46*, D608–D617. [CrossRef]

57. Chong, J.; Soufan, O.; Li, C.; Caraus, I.; Li, S.; Bourque, G.; Wishart, D.S.; Xia, J. MetaboAnalyst 4.0: Towards more transparent and integrative metabolomics analysis. *Nucleic Acids Res.* **2018**, *46*, W486–W494. [CrossRef]

<figure><img src="image_9.png" /></figure>