Next Article in Journal
Bioactive Endodontic Hydrogels: From Parameters to Personalized Medicine
Next Article in Special Issue
Anabolic Steroids Activate the NF-κB Pathway in Porcine Ovarian Putative Stem Cells Independently of the ZIP-9 Receptor
Previous Article in Journal
Antifungal Efficacy of Antimicrobial Peptide Octominin II against Candida albicans
Previous Article in Special Issue
Porcine Follicular Fluid-Derived Exosome: The Pivotal Material for Porcine Oocyte Maturation in Lipid Antioxidant Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Early Pregnancy Markers in the Serum of Ewes Identified via Proteomic and Metabolomic Analyses

1
State Key Laboratory of Farm Animal Biotech Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
2
Frontiers Science Center for Molecular Design Breeding, China Agricultural University, Beijing 100193, China
3
Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education of China, Yangzhou University, Yangzhou 225009, China
4
International Joint Research Laboratory in Universities of Jiangsu Province of China for Domestic Animal Germplasm Resources and Genetic Improvement, Yangzhou University, Yangzhou 225009, China
5
Department of Animal Sciences, Purdue University, West Lafayette, IN 47907, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2023, 24(18), 14054; https://doi.org/10.3390/ijms241814054
Submission received: 16 August 2023 / Revised: 5 September 2023 / Accepted: 6 September 2023 / Published: 13 September 2023

Abstract

:
The diagnosis of ewes’ pregnancy status at an early stage is an efficient way to enhance the reproductive output of sheep and allow producers to optimize production and management. The techniques of proteomics and metabolomics have been widely used to detect regulatory factors in various physiological processes of animals. The aim of this study is to explore the differential metabolites and proteins in the serum of pregnant and non-pregnant ewes by proteomics and metabolomics. The serum of ewes at 21, 28 and 33 days after artificial insemination (AI) were collected. The pregnancy stratus of the ewes was finally determined through ultrasound examination and then the ewes were grouped as Pregnant (n = 21) or N on-pregnant (n = 9). First, the serum samples from pregnant or non-pregnant ewes at 21 days after AI were selected for metabolomic analysis. It was found that the level of nine metabolites were upregulated and 20 metabolites were downregulated in the pregnant animals (p < 0.05). None of these differential metabolomes are suitable as markers of pregnancy due to their small foldchange. Next, the proteomes of serum from pregnant or non-pregnant ewes were evaluated. At 21 days after AI, the presence of 321 proteins were detected, and we found that the level of three proteins were upregulated and 11 proteins were downregulated in the serum of pregnant ewes (p < 0.05). The levels of serum amyloid A (SAA), afamin (AFM), serpin family A member 6 (SERPINA6) and immunoglobulin-like domain-containing protein between pregnant and non-pregnant ewes at 21-, 28- and 33-days post-AI were also analyzed via enzyme-linked immunosorbent assay (ELISA). The levels of SAA and AFM were significantly higher in pregnant ewes than in non-pregnant ewes, and could be used as markers for early pregnancy detection. Overall, our results show that SAA and AFM are potential biomarkers to determine the early pregnancy status of ewes.

1. Introduction

Accurate and timely pregnancy diagnosis has a significant economic impact on sheep farming. Early pregnancy diagnosis is crucial to determine the pregnancy status of ewes, and can be used to evaluate the conception rate, the efficacy of artificial insemination (AI) and the presence of possible diseases. Ewes identified as non-pregnant can be rebred or inseminated, thereby shortening the lambing interval and providing an economic advantage to the producer. Furthermore, reproductive health is sensitive to toxic exposure, specifically to endocrine disruptor pollutants, and have long-term adverse effects.
At present, transabdominal ultrasonography is one of the most common early pregnancy detection methods. When applied to determine the pregnancy status of ewes on day 21, the sensitivity was reported to be as low as 44.4%, and its sensitivity reaches 100% only on day 35 of gestation and afterward [1]. Other available methods such as visual early pregnancy examination, transrectal ultrasonography, near-infrared spectroscopy, hormonal assays or radiography are also not ideal [2,3]. The pregnancy-specific protein B (PSPB) has been identified, and its protein size was determined to be between 47 and 53 kDa [4]. In addition, bovine pregnancy-associated glycoprotein (bPAG) from bovine fetal cotyledon has been purified to homogeneity by HPLC [5]. PAG is synthesized by mononuclear and binuclear cells of the trophectoderm of ruminants. A portion of PAG can be released into the maternal blood circulation and detected by PAG antibodies [6]. The concentration of PAG in cows increased after pregnancy, and there was a significant difference in PAG concentration between pregnant and non-pregnant cows on day 22 following insemination [7]. Test kits have been developed based on PAG and successfully used in the cattle industry. The molecular structure of bovine and sheep PAG are relatively similar. Although PAG-based early pregnancy determination methods for sheep have been tested [8,9,10], these test kits have not been widely used in sheep farms due to high costs or limited efficiency to identify early pregnancy status. Furthermore, the reproductive health of females is sensitive to nutrition, infectious diseases, toxic exposure—specifically to endocrine disruptor pollutants—and have long-term adverse effects [11]. An advanced early pregnancy diagnosis method would help famers to distinguish the ewes with subfertility. There is clearly a need to identify new biomarkers for early pregnancy detection in sheep.
Proteomic and metabolomic analyses are powerful tools to identify the dynamic proteins and endogenous metabolites in body fluids, which could provide a deeper understanding of the biological system at the molecular level [12,13,14]. Proteins and metabolites play a crucial role in regulating the molecular pathways associated with reproductive activities [15]. Proteomics and metabolomics are emerging high-throughput approaches which have enabled researchers to identify hundreds and thousands of molecular using a very low quantity of samples from males and females suffering from infertility [16,17]. Proteomic analysis was used to reveal potential early pregnancy or litter size-related biomarkers in pigs, goats, cattle and jennies [18,19,20,21,22,23,24]. More importantly, these techniques allowed researchers to explore the molecules involved in the medical conditions of pregnant females suffering from ectopic pregnancy, hypertension disorder or diabetes mellitus [25,26,27,28,29,30,31]. Research using metabolomics was also conducted to find biomarkers in serum to foretell litter size in sheep during 7 to 70 days of gestation, or metabolites such as FSH, P4, AMH and amino acids to predict the number of ovulated oocytes following ovarian super stimulation [32,33,34]. Thus, the latest approaches empowered us to identify potential biomarkers in the very early pregnancy stages of ewes.
Early and accurate pregnancy examination to shorten the lambing interval of ewes and improve the reproductive efficiency is of great economic significance for farmers. Therefore, we designed experiments to detect early pregnancy markers in ewes using the latest proteomic and metabolomic analyses. Serum samples from pregnant (21, 28 and 33 days post-artificial insemination) and non-pregnant ewes were collected and evaluated. We found that the levels of 29 metabolites and 14 proteins were significantly different between the experimental groups at 21 days post-AI. In addition, the levels of four proteins were examined via ELISA, and it was found that serum amyloid A (SAA) and afamin (AFM) levels were significantly elevated in pregnant compared to non-pregnant ewes on 21, 28 and 33 days post-AI. Our current study provides evidence that SAA and AFM could serve as a marker to advance the detection time and determine early pregnancy in ewes.

2. Results

2.1. Serum Collection and Pregnancy Diagnosis

At 21-, 28- and 33-days post-AI, blood samples of 30 ewes were collected and the sera were isolated for future examination. At 33 days post-AI, all 30 ewes were examined using a real-time B-mode ultrasound scanner. As presented in Figure 1A–D, 21 ewes were pregnant and nine ewes were not pregnant. The serum samples were identified as belonging either to the pregnant or non-pregnant group based on the ultrasound results and then were used for proteomic and metabolomic analyses.

2.2. Metabolomic Analysis of Serum

To evaluate the metabolomic changes in the sera of ewes on day 21 post-AI, six serum samples from pregnant (P) and non-pregnant (N) ewes, respectively, were examined using Liquid Chromatography Mass Spectrometry (LC-MS). The LC/MS raw data were processed using the Progenesis QI (Waters Corporation, Milford, CT, USA) software (version 2.2). The sample information, metabolite name and mass spectral response intensity were then exported and the metabolites were identified according to databases like HMDB (http://www.hmdb.ca/), accessed on 19 September 2022. There were 1355 shared metabolites in the sera of pregnant and non-pregnant ewes (Table S1). For the positive ion, 803 metabolites were identified, and 772 of them were found in library and 446 metabolites in KEGG. Meanwhile, for negative ion, 549 metabolites were identified, 528 of which were found in library and 302 in KEGG. The correlation of the samples was then analyzed according to PLS-DA (Partial Least Squares Discriminant Analysis) using R Packages (Version 1.6.2). As presented in Figure 2A, the PLS-DA score chart showed that the degree of separation between groups P and N is significant. There were 11 unique metabolites in the pregnant group and 19 unique metabolites in the non-pregnant group (Figure 2B,C).
All the metabolites were checked with the KEGG and HMDB databases. As presented in Figure 2D, the most abundant compounds are phospholipids and lipids according to KEGG Compound Classification. Meanwhile, the lipid and amino acid metabolism also ranked among the top pathways based on KEGG Pathway analyses (Figure 2E). Similarly, when the names of the metabolites were checked in the HMDB database (Figure 2F), 492 compounds were identified as lipids and lipid-like molecules, which are compounds of the most abundance.

2.3. Differential Metabolites Detected between Pregnant and Non-Pregnant Ewes

The differential metabolites found between pregnant and non-pregnant ewes were further analyzed. As shown in Figure 3A, the levels of 20 metabolites in the sera of pregnant ewes were elevated and nine metabolites were reduced compared to that of non-pregnant ewes (listed in Table 1). The differential metabolites were then clustered according to biochemical pathways through metabolic enrichment and pathway analysis based on the KEGG database (Figure 3B,C). It was found that amino acid metabolism is ranked among the top pathways.
The expression profile of the 29 differential metabolites were then used to calculate variable importance in projection (VIP). As presented in Figure 3D, we found that indolophenanthridine,2-mercaptobenzothiazole and N-4acetylcytidine were the top three metabolites. Together with three additional ones, the level of six metabolites listed as their VIP values are presented in Figure 3E–J. It can be seen that the serum metabolism of sheep at different pregnancy statuses involves amino acid as well as lipid metabolism. Also, the fold changes of all differential metabolites are large enough for visual diagnosis kit development.

2.4. Proteomic Analysis of Serum

The proteins in the sera from three pregnant (P) and three non-pregnant (N) ewes 21 days post-AI were examined using data-independent acquisition (DIA) proteomics. Total protein annotation was then performed according to the uniport database and the functional information of proteins was comprehensively obtained. There were 321 proteins with functional annotation (Table S2) and two unique proteins were detected in the serum of pregnant or non-pregnant ewes respectively (Figure 4A). Then, the KEGG pathway analysis was performed, which revealed that these proteins were enriched in 21 KEGG signaling pathways (Figure 4B). The subcellular localization of the proteins was also determined and 147 extracellular proteins and 84 cytoplasmic proteins were found (Figure 4C).
The differences in the level of proteins in the sera of pregnant and non-pregnant ewes were also analyzed. As shown in Figure 4D, three proteins were upregulated and 11 proteins were downregulated in the sera from pregnant ewes compared to those from the non-pregnant ones. The detailed information of the 14 differential proteins can be seen in Figure 4E and Table 2. As shown in Figure 4F, the results of the KEGG pathway enrichment analysis identified an increase in the activity of the signaling pathways related to maternal pregnancy, embryo implantation and embryonic development. Among the 14 differential proteins, the levels of serum amyloid A protein (SAA), afamin (AFM), serpin family A member 6 (SERPINA6), Ig-like domain-containing protein (ILDP), albumin, beta-A globin chain β, SMB domain-containing protein and vitronectin are presented in Figure 5A–H.

2.5. Six Levels of Serum Proteins in Pregnant and Non-Pregnant Ewes at Different Stages of Pregnancy as Determined using ELISA

To confirm the changes of proteins detected by proteomics, the levels of SAA, AFM, SPERIN6 and ILDP in serum samples from pregnant and non-pregnant ewes at 21-, 28- and 33-days post-AI were examined via an ELISA kit. The serum SAA in the pregnant ewes was significantly higher than that of the non-pregnant ewes at 21d (3.37 ± 0.22 vs. 2.70 ± 0.16 ng/mL), 28d (2.63 ± 0.43 vs. 1.99 ± 0.27 ng/mL) and 33d (3.55 ± 0.34 vs.1.96 ± 0.22 ng/mL) post-AI (Figure 6A,B). Similarly, the serum AFM in the pregnant ewes was significantly higher than that in the non-pregnant ewes at 21d (91.73 ± 4.91 vs. 79.72 ± 3.50), 28d (75.93 ± 5.09 vs. 66.22 ± 4.68) and 33d (94.26 ± 4.57 vs. 65.86 ± 3.61) post-AI (Figure 6A, B). The levels of SPERIN6 in the sera of pregnant ewes were significantly lower than that of non-pregnant ewes, only at 33 days post-AI. There was no difference in the serum ILDP levels between pregnant and non-pregnant ewes. Combined with the pathways enriched by the above differential proteins, SAA and AFM can be used as markers for the early pregnancy diagnosis of ewes.

3. Discussion

This study aimed to identify biomarkers for the early pregnancy diagnosis of sheep. The analysis of pregnancy-associated molecules has been the subject of many studies; using different species as the identification of such proteins offers advantages that include not only early pregnancy diagnosis, but also the prediction of litter size and the determination of the presence of reproductive diseases. Accurate and timely pregnancy diagnosis has a significant economic impact on sheep farming. In this set of experiments, the sera of ewes at 21-, 28- and 33-days post-AI were collected. The pregnancy status of ewes was confirmed via ultrasound 33 days post-AI. The serum samples from pregnant and non-pregnant ewes on 21 days post-AI were then used for metabolomic and proteomic examination. We identified 29 metabolites and 14 proteins differentially expressed between the two groups. However, none of the metabolites were suitable as markers for pregnancy due to the limited fold changes and no further investigation was conducted on these molecules. Among the 14 differentially expressed proteins, four proteins, including SAA, AFM, SPERIN6 and ILDP, were selected for ELISA detection. The results showed that the levels of two proteins (SAA and AFM) in the sera were significantly different between pregnant and non-pregnant ewes at 21, 28 and 33 days of gestation; we propose that these proteins can be used as protein markers to detect early pregnancy in ewes.
Knowing whether or not an ewe is pregnant soon after mating or artificial insemination is critical for sheep reproductive management. Numerous clinical and immunologic methods have been developed for sheep pregnancy diagnosis. Previously, methods based on ELISA had been used to assay hormonal, early pregnancy factor (EPF) and pregnancy-associated antigens in order to detect the early pregnancy of ewes [35,36,37,38]. However, the pregnancy diagnosis of ewes in practice still depends primarily on the identification of animals that do not return to estrus or the use of B-mode ultrasonography; these methods provide a pregnancy status from 35 days of pregnancy onwards. To shorten the interval between two successive inseminations and minimize the economic losses due to non-pregnant animals, new methods like the visual ELISA-PAG test, originally developed for cattle, have been tested on sheep; these tests proved highly accurate and efficient [8,35,36]. However, these immunological methods have not been widely used on sheep farms due to their high costs and complicated operation processes. Hence, new biomarkers and easy examination procedures for early pregnancy diagnosis are urgently required in the sheep industry. In this set of experiments, high-throughput metabolomic and proteomic methods were applied to detected unique molecules in the sera of pregnant and non-pregnant ewes. Metabolomics and proteomics are sensitive enough to detect thousands of biological molecules in one sample at the same time; they can be used to screen the relevant markers to detect early pregnancy. With this approach, we identified a total of 321 serum proteins in the two groups, among which 14 proteins were significantly different between pregnant and non-pregnant ewes. Four proteins showing obvious fold changes (SAA, AFM, ILDP and SPERIN6) were further analyzed via an ELISA kit. It was demonstrated that the level of SAA and AFM in the sera of pregnant ewes were significantly elevated since day 21 post-AI.
During early pregnancy, the female animal needs to adjust its endocrine system, immune system, metabolites and hormones to allow for the fetus to establish a connection with the maternal body [39]. The expression of many proteins’ changes in the female during the early stages of pregnancy. In our experiments, SAA and AFM proved to be suitable indicators of early pregnancy. When the body is injured and inflammation occurs, the SAA level increases and regulates the adhesion, migration, proliferation and aggregation of cells [40]. Interestingly, others found that elevated SAA levels may be related to abnormalities in the decidualization process, which may lead to maternal infertility or spontaneous abortion [41]. Therefore, further functional experiments are required to elucidate the role of SAA in the early pregnancy of ewes. The other protein, afamin (AFM), is a glycoprotein with vitamin E-binding properties and a putative function in fertility [42]. In previous clinical studies, afamin levels increased in the maternal serum during pregnancy and was proposed to be a potential predictor of preeclampsia [43,44]. Others also conducted a large number of studies on the relationship between serum AFM levels and inflammatory bowel diseases, tumor diseases, coronary heart disease and ovarian cancer [45,46]. AFM is a functional protein in humans, but its functions in ewes is still unclear.
In the current research, serpin family A member 6 (SERPINA6) was identified as another protein present at elevated concentrations in the sera of pregnant ewes. Currently, investigations regarding SERPINA6 focus on its function as a steroid-binding protein in the blood of mammals, reptiles, amphibians and birds [47,48]. In humans, the levels of cortisol in the plasma is associated with cardio-metabolic, inflammatory and neuro-cognitive traits or even diseases. It has been reported that a genetic variant of the SERPINA6 gene (rs7161521) is associated with diurnal and stress-induced hypothalamic–pituitary–adrenal (HPA) axis activity in children [48]. Meanwhile, other studies have found that a SERPINA6 gene variant encodes corticosteroid-binding globulin (CBG), a protein with corticosteroid-binding properties in blood [49,50]. However, the role of SERPINA6 in animals during pregnancy has not been studied. Thus, more research is needed to elucidate whether SERPINA6 is involved in the regulation of the pregnancy of ewes or its presence in serum is simply due to stress. Meanwhile, the level of an Ig-like domain-containing protein (accession number: A0A452FW61) was significantly reduced in pregnant ewes; the exact identity of this protein could not be determined due to technical limitations.
Taken together, the SAA and AFM proteins are differentially expressed between pregnant and non-pregnant ewes and can be used as biomarkers for early pregnancy detection in the sheep.

4. Materials and Methods

4.1. Animals and Experimental Design

The experiment was conducted using adult female Hu sheep (Ovis aries) at the experimental facilities of the Inner Mongolia Golden Grassland Ecological Technology Group Co., Ltd., Bayannur, China (latitude 40°13′ N, longitude 105°12′ E) in July 2022. All experimental protocols concerning the handling of animals were performed in accordance with the requirements of the Institutional Animal Care and Use Committee at the China Agricultural University.

4.2. Collection of Blood Samples from Ewes

Blood was collected from ewes at 21-, 28- and 33-days post-AI. After all the blood samples were placed in a coagulation-promoting tube for 2 h, they were centrifuged at 4000 rpm/s for 10 min. After centrifugation, the supernatant was collected and stored at a −80 °C freezer for later use.

4.3. Serum Metabolomic Analysis

4.3.1. Serum Samples Selection and Quality Control

The sera with uniform color from 6 pregnant sheep (P1, P2, P3, P4, P5 and P6) and 6 non-pregnant sheep (N1, N2, N3, N4, N5 and N6) on 21 days post-AI were selected for metabolome analysis. A pooled quality control (QC) sample was prepared by mixing equal volumes of all samples. The QC samples were disposed and tested in the same manner as the analytic samples to represent the whole sample set, which was used to monitor the stability of the analysis.

4.3.2. LC-MS/MS Analysis of the Serum Samples

The metabolomics analysis of the ewes’ serum samples was performed as described previously [51,52]. Briefly, 100 μL serum sample was added to a 1.5 mL centrifuge tube with 400 μL solution (acetonitrile/methanol = 1:1(v:v)) containing 0.02 mg/mL internal standard (L-2-chlorophenylalanine) to extract metabolites. The samples were mixed by vortexing for 30 s and then low-temperature sonicated for 30 min (5 °C, 40 KHz). The samples were placed at −20 °C for 30 min to precipitate the proteins. The samples were then centrifuged for 15 min (4 °C, 13,000× g). The supernatant was removed and dried under nitrogen and then re-solubilized with 100 μL solution and extracted via ultrasonication for 5 min (5 °C, 40 KHz), followed by centrifugation at 13,000× g and 4 °C for 10 min. The supernatant was transferred to sample vials for LC-MS/MS analysis. The LC-MS/MS analysis of the samples was conducted on a SCIEX UPLC-Triple TOF 5600 system equipped with an ACQUITY HSS T3 column (Waters, Milford, MA, USA) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Data acquisition was performed with the Data Dependent Acquisition (DDA) mode.

4.3.3. Data Preprocessing and Annotation

The raw data of LC/MS were preprocessed using Progenesis QI (Waters Corporation, Milford, CT, USA) software (version 2.2), and a three-dimensional data matrix in CSV format was exported. Internal standard peaks, as well as any known false positive peaks were removed from the data matrix and peak-pooled. At the same time, the metabolites were searched and identified based on database HMDB (http://www.hmdb.ca/), Metlin (https://metlin.scripps.edu/) and Majorbio database, accessed on 19 September 2022. In order to reduce the errors caused by sample preparation and instrument instability, variables with relative standard deviation (RSD) > 30% of QC samples were removed, and normalized as log10 logarithmization to obtain the data matrix for subsequent analysis.

4.3.4. Differential Metabolites Analysis

The variance analysis on the matrix file were performed after data preprocessing. The R package (Version 1.6.2) performed orthogonal least partial squares discriminant analysis (OPLS-DA), and Student’s t-test and fold difference analysis were performed. The variable importance in projection (VIP) was obtained via the OPLS-DA model and the p-value was determined using the Student’s t test; the metabolites with VIP > 1, p < 0.05 were regarded as significantly different metabolites. All the metabolites were screened and differential metabolites among two groups were summarized and mapped into their biochemical pathways through metabolic enrichment and pathway analysis based on database search (KEGG, http://www.genome.jp/kegg/), accessed on 1 October 2022. The data were analyzed through the free online platform of Majorbio Cloud platform (cloud.majorbio.com), accessed on 1 October 2022.

4.4. Serum Proteomic Analysis

4.4.1. Serum Samples Selection and Quality Control

The serum samples from 3 pregnant and 3 non-pregnant ewes at 21 days post-AI were selected. The serum proteins were extracted in fresh lysis buffer and 0.5% sodium deoxycholate with 1 × phosphatase inhibitor cocktail (PhosSTOP, Sigma-Aldrich, St. Louis, MO, USA). The protein concentration was measured via the BCA assay (Thermo Fisher Scientific, Waltham, MA, USA) and confirmed using Coomassie-stained short SDS gel as previously described [53,54]. Protein quantification was performed according to the kit protocol. After protein quantification, SDS-PAGE electrophoresis was performed.

4.4.2. Proteolytic Hydrolysis and Peptide Quantification

The protein samples were taken and added with lysate, and then triethylammonium bicarbonate buffer (TEAB) was added at a final concentration of 100 mM. The iodoacetamide at a final concentration of 40 mM was added and reacted at room temperature for 40 min. Precooled acetone (acetone/sample v:v = 6:1) was added to each tube and kept at −20 °C for 4 h. The samples were then centrifuged at 10,000× g for 20 min and the pellets were collected. After trypsin digestion, peptide from the samples was drained via rotation vacuum concentration (Christ RVC 2-25, Christ, Erfurt, Germany). The peptide was extracted and re-dissolved with 0.1% trifluoroacetic acid (TFA). The supernatant was desalted with Sep-Pak C18 filter cartridge (Waters, Milford, MA, USA) and dried. Peptide quantification was performed according to the peptide quantification kit (Thermo Fisher Scientific, Waltham, MA, USA).

4.4.3. DIA Mass Detection of a Single Sample

Based on peptide quantification results, the peptide samples were redissolved in spectrometry loading buffer (2% ACN with 0.1% formic acid), including appropriate iRT peptide, which was used to calibrate retention time, and were analyzed using an EASY-nLC system (Thermo Fisher Scientific, Waltham, MA, USA) with a timsTOF Pro2 mass spectrometer (Bruker, Bremen, Germany) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Then, the data-independent acquisition (DIA) data were acquired using a timsTOF Pro2 mass spectrometer operated in DIA-PASEF mode. MS data were collected over an m/z range of 400 to 1200 and an ion mobility range of 0.57 to 1.47 Vs·cm−2. Both accumulation time and ramp time were set to 100 ms. During MS/MS data collection, each cycle contained one MS and ten PASEF MS/MS scans. Exclusion was active after 0.4 min. A total of 64 DIA-PASEF windows were used (25 Th isolation windows).

4.4.4. Protein Identification

Spectronaut software (Version 14) was used to analyze the DIA-PASEF raw data based on the spectra library generated by DDA-PASEF data. Retention times were corrected by iRT and 6 peptides per protein and 3 ions per peptide were selected for quantitative analysis. The parameters are as following: Protein FDR ≤ 0.01, Peptide FDR ≤ 0.01, Peptide Confidence ≥ 99%, XIC width ≤ 75 ppm. The shared and modified peptides were excluded, and the peak areas were calculated and summed to generate the quantitative results.

4.4.5. Statistical and Bioinformatics Analysis

The similarity and difference of proteins between different serum samples were compared via sample correlation and principal component analysis using R Packages software (Version 1.6.2) on Majorbio Cloud platform (https://cloud.majorbio.com), accessed on 1 October 2022. p-values were corrected for multiple testing using the Student’s t-test. The proteins with a p-value < 0.05 and |log2FoldChange| > 1 was considered as significantly differentially expressed proteins. GO enrichment and KEGG pathway enrichment analyses of these predicted target genes were performed using R Packages.
The differentially expressed proteins (DEPs) were identified as thresholds of fold change (>1.2 or <0.83) and p-value < 0.05. Functional annotation of all identified proteins was performed using GO (http://geneontology.org/) and KEGG pathway (http://www.genome.jp/kegg/), accessed on 1 October 2022. DEPs were further used for GO and KEGG enrichment analyses. Protein–protein interaction analysis was performed using the String v11.5.

4.5. Evaluation of Protein Levels in Sera of Ewes at Different Time Points

Blood samples of ewes at 21-, 28- and 33-days post-AI were collected and the serum from each sample was selected for analysis. The levels of SAA, AFM, SPERIN6 and ILDP in the samples were detected via ELISA kits purchased from Beijing Beijianxinchuangyuan Biotechnology Co., Ltd., Beijing, China. The analyses were performed according to the manufacture’s protocol. Briefly, 10 μL of the serum samples and 40 μL of diluent were added to the wells, and then the HRP-conjugate reagent was added to each well and incubated for 60 min at 37 °C. After that, 400 μL wash solution filled each well and this was repeated 5 times. Then, each well was filled with 50 μL chromogen solution A and 50 μL solution B, which were mixed gently and incubated for 15 min at 37 °C. Finally, 50 μL stop solution was added to each well. The Optical Density (O.D.) value of each sample was detected at 450 nm using a microtiter plate reader and the data were calculated based on the standard curve using the O.D. value for each well.

4.6. Statistical Analyses

The statistical analyses for proteomic and metabolomic analyses are described in each section. The levels of differentially expressed metabolites and proteins were evaluated using one-way ANOVA. Statistical significance at p < 0.05 is marked as * and at p < 0.01 is marked as **.

5. Conclusions

Early pregnancy diagnosis allows for optimizing the production and timely management of decisions, offering a greater reproductive output in sheep farming. In this manuscript, the researchers proposed to determine the unique molecules in the sera of pregnant ewes at 21 days post-AI. Altogether, 29 differentially expressed metabolites and 14 differentially expressed proteins were detected in the sera of pregnant ewes, compared with non-pregnant ones. The levels of four differentially expressed proteins were re-examined via ELISA, and it was found that the levels of SAA and AFM were significantly increased in the sera of ewes at day 21, day 28 and day 33 of gestation. Therefore, SAA and AFM can be used as markers of early pregnancy in ewes. Finally, precisely unraveling proteomic and metabolomic biomarkers and hallmarks related to the establishment of early pregnancy might be helpful for recognizing and improving the efficiency of post-implantation interplay between in utero-endometrial and placental anatomo-histological compartments and fetuses, following surgical or transvaginal/transcervical procedures used for the transfer of ex vivo-produced ovine and other mammalian embryos propagated by a variety of modern assisted reproductive technologies, such as somatic cell nuclear transfer and in vitro fertilization.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241814054/s1.

Author Contributions

Conceptualization, L.Z.; methodology, Y.Z. and F.X.; investigation, F.X. and L.S.; data curation, F.X., P.J. and W.M.; writing—original draft preparation, L.Z. and G.L.; writing—review and editing, S.G., X.L., W.S. and Z.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (32002181), National Key R&D Program of China (2022ZD0401405), Chinese Universities Scientific Fund(2022TC143), the Open Project Program of International Joint Research Laboratory in Universities of Jiangsu Province of China for Domestic Animal Germplasm Resources and Genetic Improvement (2022) and the 2115 Talent Development Program of China Agricultural University.

Institutional Review Board Statement

All experimental protocols concerning the handling of animals were performed in accordance with the requirements of the Institutional Animal Care and Use Committee at the China Agricultural University.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AFM, afamin; AI, artificial insemination; bPAG, bovine pregnancy-associated glycoprotein; CBG, corticosteroid-binding globulin; DDA, Data Dependent Acquisition; DEPs, differentially expressed proteins; DIA, data-independent acquisition; ELISA, enzyme-linked immunosorbent assay; EPF, early pregnancy factor; HPA, hypothalamic–pituitary–adrenal; ILDP, Ig-like domain-containing protein; LC-MS, Liquid Chromatography Mass Spectrometry; N, non-pregnant; OPLS-DA, orthogonal least partial squares discriminant analysis; P, pregnant; PLS-DA, Partial Least Squares Discriminant Analysis; PSPB, pregnancy-specific protein B; QC, quality control; RSD, relative standard deviation; SAA, serum amyloid A; SERPINA6, serpin family A member 6; TEAB, triethylammonium bicarbonate buffer; TFA, trifluoroacetic acid; VIP, variable importance in projection.

References

  1. Jones, A.K.; Gately, R.E.; McFadden, K.K.; Zinn, S.A.; Govoni, K.E.; Reed, S.A. Transabdominal ultrasound for detection of pregnancy, fetal and placental landmarks, and fetal age before Day 45 of gestation in the sheep. Theriogenology 2016, 85, 939−945.e1. [Google Scholar] [CrossRef] [PubMed]
  2. Andueza, D.; Alabart, J.L.; Lahoz, B.; Muñoz, F.; Folch, J. Early pregnancy diagnosis in sheep using near-infrared spectroscopy on blood plasma. Theriogenology 2014, 81, 509–513. [Google Scholar] [CrossRef] [PubMed]
  3. Petrujkić, B.T.; Cojkić, A.; Petrujkić, K.; Jeremić, I.; Mašulović, D.; Dimitrijević, V.; Savić, M.; Pešić, M.; Beier, R.C. Transabdominal and transrectal ultrasonography of fetuses in Württemberg ewes: Correlation with gestational age. Anim. Sci. J. 2016, 87, 197–201. [Google Scholar] [CrossRef] [PubMed]
  4. Butler, J.E.; Hamilton, W.C.; Sasser, R.G.; Ruder, C.A.; Hass, G.M.; Williams, R.J. Detection and partial characterization of two bovine pregnancy-specific proteins. Biol. Reprod. 1982, 26, 925–933. [Google Scholar] [CrossRef] [PubMed]
  5. Zoli, A.P.; Beckers, J.F.; Wouters-Ballman, P.; Closset, J.; Falmagne, P.; Ectors, F. Purification and characterization of a bovine pregnancy-associated glycoprotein. Biol. Reprod. 1991, 45, 1–10. [Google Scholar] [CrossRef] [PubMed]
  6. Sousa, N.M.; Ayad, A.; Beckers, J.F.; Gajewski, Z. Pregnancy-associated glycoproteins (PAG) as pregnancy markers in the ruminants. J. Physiol. Pharmacol. 2006, 57 (Suppl. S8), 153–171. [Google Scholar] [PubMed]
  7. Friedrich, M.; Holtz, W. Establishment of an ELISA for measuring bovine pregnancy-associated glycoprotein in serum or milk and its application for early pregnancy detection. Reprod. Domest. Anim. 2010, 45, 142–146. [Google Scholar] [CrossRef] [PubMed]
  8. Steckeler, P.; Weber, F.; Zerbe, H.; Rieger, A.; Voigt, K. Evaluation of a bovine visual pregnancy test for the detection of pregnancy-associated glycoproteins in sheep. Reprod. Domest. Anim. 2019, 54, 280–288. [Google Scholar] [CrossRef]
  9. Karen, A.; Amiri, B.E.; Beckers, J.F.; Sulon, J.; Taverne, M.A.; Szenci, O. Comparison of accuracy of transabdominal ultrasonography, progesterone and pregnancy-associated glycoproteins tests for discrimination between single and multiple pregnancy in sheep. Theriogenology 2006, 66, 314–322. [Google Scholar] [CrossRef]
  10. Barbato, O.; Sousa, N.M.; Debenedetti, A.; Canali, C.; Todini, L.; Beckers, J.F. Validation of a new pregnancy-associated glycoprotein radioimmunoassay method for the detection of early pregnancy in ewes. Theriogenology 2009, 72, 993–1000. [Google Scholar] [CrossRef]
  11. Bianchi, S.; Nottola, S.A.; Torge, D.; Palmerini, M.G.; Necozione, S.; Macchiarelli, G. Association between Female Reproductive Health and Mancozeb: Systematic Review of Experimental Models. Int. J. Environ. Res. Public Health 2020, 17, 2580. [Google Scholar] [CrossRef] [PubMed]
  12. Domon, B.; Aebersold, R. Mass spectrometry and protein analysis. Science 2006, 312, 212–217. [Google Scholar] [CrossRef] [PubMed]
  13. Lippolis, J.D.; Reinhardt, T.A. Centennial paper: Proteomics in animal science. J. Anim. Sci. 2008, 86, 2430–2441. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, Y.M.; Heyman, H.M. Mass Spectrometry-Based Metabolomics. Methods Mol. Biol. 2018, 1775, 107–118. [Google Scholar] [PubMed]
  15. Horgan, R.P.; Clancy, O.H.; Myers, J.E.; Baker, P.N. An overview of proteomic and metabolomic technologies and their application to pregnancy research. BJOG 2009, 116, 173–181. [Google Scholar] [CrossRef] [PubMed]
  16. Panner Selvam, M.K.; Finelli, R.; Agarwal, A.; Henkel, R. Proteomics and metabolomics—Current and future perspectives in clinical andrology. Andrologia 2021, 53, e13711. [Google Scholar] [CrossRef] [PubMed]
  17. Dong, R.; Ye, N.; Zhao, S.; Wang, G.; Zhang, Y.; Wang, T.; Zou, P.; Wang, J.; Yao, T.; Chen, M.; et al. Studies on Novel Diagnostic and Predictive Biomarkers of Intrahepatic Cholestasis of Pregnancy through Metabolomics and Proteomics. Front. Immunol. 2021, 12, 733225. [Google Scholar] [CrossRef]
  18. Xin, H.; Meng, F.; Li, B.; Wang, S.; Hu, B.; Li, J. Global proteomic analysis of serum during early pregnancy in the pig using data-independent acquisition mass spectrometry with verification by parallel reaction monitoring. Reprod. Fertil. Dev. 2022, 34, 1115–1127. [Google Scholar] [CrossRef]
  19. De, A.; Ali, M.A.; Chutia, T.; Onteru, S.K.; Behera, P.; Kalita, G.; Kumar, S.; Gali, J.M. Comparative serum proteome analysis reveals potential early pregnancy-specific protein biomarkers in pigs. Reprod. Fertil. Dev. 2019, 31, 613–631. [Google Scholar] [CrossRef]
  20. Yang, P.; Chen, X.; Tian, X.; Zhou, Z.; Zhang, Y.; Tang, W.; Fu, K.; Zhao, J.; Ruan, Y. A Proteomic Study of the Effect of N-acetylcysteine on the Regulation of Early Pregnancy in Goats. Animals 2022, 12, 2439. [Google Scholar] [CrossRef]
  21. Lee, J.E.; Lee, J.Y.; Kim, H.R.; Shin, H.Y.; Lin, T.; Jin, D.I. Proteomic Analysis of Bovine Pregnancy-specific Serum Proteins by 2D Fluorescence Difference Gel Electrophoresis. Asian-Australas. J. Anim. Sci. 2015, 28, 788–795. [Google Scholar] [CrossRef] [PubMed]
  22. Panda, B.S.K.; Mohapatra, S.K.; Chaudhary, D.; Alhussien, M.N.; Kapila, R.; Dang, A.K. Proteomics and transcriptomics study reveals the utility of ISGs as novel molecules for early pregnancy diagnosis in dairy cows. J. Reprod. Immunol. 2020, 140, 103148. [Google Scholar] [CrossRef] [PubMed]
  23. Ceciliani, F.; Lecchi, C.; Urh, C.; Sauerwein, H. Proteomics and metabolomics characterizing the pathophysiology of adaptive reactions to the metabolic challenges during the transition from late pregnancy to early lactation in dairy cows. J. Proteom. 2018, 178, 92–106. [Google Scholar] [CrossRef] [PubMed]
  24. Deng, L.; Han, Y.; Tang, C.; Liao, Q.; Li, Z. Label-Free Mass Spectrometry-Based Quantitative Proteomics Analysis of Serum Proteins during Early Pregnancy in Jennies (Equus asinus). Front. Vet. Sci. 2020, 7, 569587. [Google Scholar] [CrossRef] [PubMed]
  25. Gerton, G.L.; Fan, X.J.; Chittams, J.; Sammel, M.; Hummel, A.; Strauss, J.F.; Barnhart, K. A serum proteomics approach to the diagnosis of ectopic pregnancy. Ann. N. Y. Acad. Sci. 2004, 1022, 306–316. [Google Scholar] [CrossRef]
  26. Nie, L.; Xin, S.; Zheng, J.; Luo, Y.; Zou, Y.; Liu, X.; Chen, H.; Lei, X.; Zeng, X.; Lai, H. DIA-based proteomics analysis of serum-derived exosomal proteins as potential candidate biomarkers for intrahepatic cholestasis in pregnancy. Arch. Gynecol. Obstet. 2023, 308, 79–89. [Google Scholar] [CrossRef] [PubMed]
  27. Jansen, B.C.; Bondt, A.; Reiding, K.R.; Lonardi, E.; de Jong, C.J.; Falck, D.; Kammeijer, G.S.; Dolhain, R.J.; Rombouts, Y.; Wuhrer, M. Pregnancy-associated serum N-glycome changes studied by high-throughput MALDI-TOF-MS. Sci. Rep. 2016, 6, 23296. [Google Scholar] [CrossRef] [PubMed]
  28. Zeng, S.; Han, M.; Jiang, M.; Liu, F.; Hu, Y.; Long, Y.; Zhu, C.; Zeng, F.; Gan, Q.; Ye, W.; et al. Serum complement proteomics reveal biomarkers for hypertension disorder of pregnancy and the potential role of Clustering. Reprod. Biol. Endocrinol. 2021, 19, 56. [Google Scholar] [CrossRef]
  29. Vora, N.; Kalagiri, R.; Mallett, L.H.; Oh, J.H.; Wajid, U.; Munir, S.; Colon, N.; Raju, V.N.; Beeram, M.R.; Uddin, M.N. Proteomics and Metabolomics in Pregnancy-An Overview. Obstet. Gynecol. Surv. 2019, 74, 111–125. [Google Scholar] [CrossRef]
  30. Chen, M.; Guan, Y.; Huang, R.; Duan, J.; Zhou, J.; Chen, T.; Wang, X.; Xia, Y.; London, S.J. Associations between the Maternal Exposome and Metabolome during Pregnancy. Environ. Health Perspect. 2022, 130, 37003. [Google Scholar] [CrossRef]
  31. Alesi, S.; Ghelani, D.; Rassie, K.; Mousa, A. Metabolomic Biomarkers in Gestational Diabetes Mellitus: A Review of the Evidence. Int. J. Mol. Sci. 2021, 22, 5512. [Google Scholar] [CrossRef] [PubMed]
  32. Goldansaz, S.A.; Markus, S.; Plastow, G.; Wishart, D.S. Predictive blood biomarkers of sheep pregnancy and litter size. Sci. Rep. 2022, 12, 10307. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, Q.; Wang, C.; Wang, L.; Feng, R.; Guo, Y.; Feng, S.; Zhang, L.; Zheng, Z.; Su, G.; Fan, L.; et al. Correlation analysis of serum reproductive hormones and metabolites during multiple ovulation in sheep. BMC Vet. Res. 2022, 18, 290. [Google Scholar] [CrossRef]
  34. Chiaradia, E.; Avellini, L.; Tartaglia, M.; Gaiti, A.; Just, I.; Scoppetta, F.; Czentnar, Z.; Pich, A. Proteomic evaluation of sheep serum proteins. BMC Vet. Res. 2012, 8, 66. [Google Scholar] [CrossRef] [PubMed]
  35. de Miranda E Silva Chaves, C.; Dias da Costa, R.L.; Roncato Duarte, K.M.; Machado, D.C.; Paro de Paz, C.C.; Beltrame, R.T. Visual ELISA for detection of pregnancy-associated glycoproteins (PAGs) in ewe serum. Theriogenology 2017, 97, 78–82. [Google Scholar] [CrossRef] [PubMed]
  36. de Miranda E Silva Chaves, C.; da Costa, R.L.D.; Duarte, K.M.R.; Beltrame, R.T.; Quirino, C.R. Evaluation of a cattle rapid test for early pregnancy diagnosis in sheep. Trop. Anim. Health Prod. 2020, 52, 1345–1349. [Google Scholar] [CrossRef] [PubMed]
  37. Mauffré, V.; Grimard, B.; Eozenou, C.; Inghels, S.; Silva, L.; Giraud-Delville, C.; Capo, D.; Sandra, O.; Constant, F. Interferon stimulated genes as peripheral diagnostic markers of early pregnancy in sheep: A critical assessment. Animal 2016, 66, 1856–1863. [Google Scholar] [CrossRef] [PubMed]
  38. Milisits-Németh, T.; Balogh, O.G.; Egerszegi, I.; Kern, L.; Sasser, R.G.; Gábor, G. Detection of pregnancy in sheep using an ELISA for pregnancy-specific protein B. Acta Vet. Hung. 2018, 66, 329–336. [Google Scholar] [CrossRef]
  39. Entrican, G.; Wheelhouse, N.M. Immunity in the female sheep reproductive tract. Vet. Res. 2006, 37, 295–309. [Google Scholar] [CrossRef]
  40. Canisso, I.F.; Ball, B.A.; Cray, C.; Williams, N.M.; Scoggin, K.E.; Davolli, G.M.; Squires, E.L.; Troedsson, M.H. Serum amyloid A and haptoglobin concentrations are increased in plasma of mares with ascending placentitis in the absence of changes in peripheral leukocyte counts or fibrinogen concentration. Am. J. Reprod. Immunol. 2014, 72, 376–385. [Google Scholar] [CrossRef]
  41. Chiriac, I.E.; Vilceanu, N.; Maghiar, A.; Andrei, C.; Hanganu, B.; Daina, L.G.; Dunarintu, S.; Buhas, L.C. Maternal Serum Amyloid A as a Marker of Preterm Birth/PROM: A Systematic Review and Meta-Analysis. Medicina 2023, 59, 1025. [Google Scholar] [CrossRef] [PubMed]
  42. Hubalek, M.; Buchner, H.; Mörtl, M.G.; Schlembach, D.; Huppertz, B.; Firulovic, B.; Köhler, W.; Hafner, E.; Dieplinger, B.; Wildt, L.; et al. The vitamin E-binding protein afamin increases in maternal serum during pregnancy. Clin. Chim. Acta 2014, 434, 41–47. [Google Scholar] [CrossRef] [PubMed]
  43. Köninger, A.; Enekwe, A.; Mach, P.; Andrikos, D.; Schmidt, B.; Frank, M.; Birdir, C.; Kimmig, R.; Gellhaus, A.; Dieplinger, H. Afamin: An early predictor of preeclampsia. Arch. Gynecol. Obstet. 2018, 298, 1009–1016. [Google Scholar] [CrossRef] [PubMed]
  44. Lichenstein, H.S.; Lyons, D.E.; Wurfel, M.M.; Johnson, D.A.; McGinley, M.D.; Leidli, J.C.; Trollinger, D.B.; Mayer, J.P.; Wright, S.D.; Zukowski, M.M. Afamin is a new member of the albumin, alpha-fetoprotein, and vitamin D-binding protein gene family. J. Biol. Chem. 1994, 269, 18149–18154. [Google Scholar] [CrossRef] [PubMed]
  45. Tramontana, A.; Dieplinger, B.; Stangl, G.; Hafner, E.; Dieplinger, H. First trimester serum afamin concentrations are associated with the development of pre-eclampsia and gestational diabetes mellitus in pregnant women. Clin. Chim. Acta 2018, 476, 160–166. [Google Scholar] [CrossRef] [PubMed]
  46. Dieplinger, H.; Dieplinger, B. Afamin—A pleiotropic glycoprotein involved in various disease states. Clin. Chim. Acta 2015, 446, 105–110. [Google Scholar] [CrossRef] [PubMed]
  47. Utge, S.; Räikkönen, K.; Kajantie, E.; Lipsanen, J.; Andersson, S.; Strandberg, T.; Reynolds, R.M.; Eriksson, J.G.; Lahti, J. Polygenic risk score of SERPINA6/SERPINA1 associates with diurnal and stress-induced HPA axis activity in children. Psychoneuroendocrinology 2018, 93, 1–7. [Google Scholar] [CrossRef] [PubMed]
  48. Vashchenko, G.; Das, S.; Moon, K.M.; Rogalski, J.C.; Taves, M.D.; Soma, K.K.; Van Petegem, F.; Foster, L.J.; Hammond, G.L. Identification of Avian Corticosteroid-binding Globulin (SerpinA6) Reveals the Molecular Basis of Evolutionary Adaptations in SerpinA6 Structure and Function as a Steroid-binding Protein. J. Biol. Chem. 2016, 291, 11300–11312. [Google Scholar] [CrossRef]
  49. Anderson, L.N.; Briollais, L.; Atkinson, H.C.; Marsh, J.A.; Xu, J.; Connor, K.L.; Matthews, S.G.; Pennell, C.E.; Lye, S.J. Investigation of genetic variants, birthweight and hypothalamic-pituitary-adrenal axis function suggests a genetic variant in the SERPINA6 gene is associated with corticosteroid binding globulin in the western Australia pregnancy cohort (Raine) study. PLoS ONE 2014, 9, e92957. [Google Scholar] [CrossRef]
  50. Lin, H.Y.; Underhill, C.; Lei, J.H.; Helander-Claesson, A.; Lee, H.Y.; Gardill, B.R.; Muller, Y.A.; Wang, H.; Hammond, G.L. High frequency of SERPINA6 polymorphisms that reduce plasma corticosteroid-binding globulin activity in Chinese subjects. J. Clin. Endocrinol. Metab. 2012, 97, 678–686. [Google Scholar] [CrossRef]
  51. Kong, X.; Guo, Z.; Yao, Y.; Xia, L.; Liu, R.; Song, H.; Zhang, S. Acetic acid alters rhizosphere microbes and metabolic composition to improve willows drought resistance. Sci. Total Environ. 2022, 844, 157132. [Google Scholar] [CrossRef]
  52. Li, C.; Al-Dalali, S.; Zhou, H.; Xu, B. Influence of curing on the metabolite profile of water-boiled salted duck. Food. Chem. 2022, 397, 133752. [Google Scholar] [CrossRef]
  53. Graves, P.R.; Haystead, T.A. Molecular biologist’s guide to proteomics. Microbiol. Mol. Biol. Rev. 2002, 66, 39–63. [Google Scholar] [CrossRef]
  54. Mun, D.G.; Vanderboom, P.M.; Madugundu, A.K.; Garapati, K.; Chavan, S.; Peterson, J.A.; Saraswat, M.; Pandey, A. DIA-Based Proteome Profiling of Nasopharyngeal Swabs from COVID-19 Patients. J. Proteome Res. 2021, 20, 4165–4175. [Google Scholar] [CrossRef]
Figure 1. Representative images of the ultrasound examination: (A,B) ultrasound images of a pregnant ewe on day 33 post-AI; (C,D) ultrasound images of a non-pregnant ewe on day 33 post-AI.
Figure 1. Representative images of the ultrasound examination: (A,B) ultrasound images of a pregnant ewe on day 33 post-AI; (C,D) ultrasound images of a non-pregnant ewe on day 33 post-AI.
Ijms 24 14054 g001
Figure 2. Analysis of the metabolites detected in the serum samples. (A) PLS-DA analysis score chart; (B,C) Venn diagrams for sample comparison. (D) Classification of metabolites based on KEGG Compounds analyses. (E) Classification of pathways based on KEGG. (F) Pie chart of compounds classification based on the HMDB database.
Figure 2. Analysis of the metabolites detected in the serum samples. (A) PLS-DA analysis score chart; (B,C) Venn diagrams for sample comparison. (D) Classification of metabolites based on KEGG Compounds analyses. (E) Classification of pathways based on KEGG. (F) Pie chart of compounds classification based on the HMDB database.
Ijms 24 14054 g002
Figure 3. Analysis of the differential metabolites between pregnant and non-pregnant ewes on 21 days post-AI. (A) Volcano plot for differential metabolites. The red dots represent the upregulated metabolites and the blue dots represent the downregulated metabolites. (B) KEGG classification of the 29 differential metabolites. (C) KEGG pathway enrichment analyses of the 29 differential metabolites. (D) Expression profile of the 29 differential metabolites and their VIP values. (EJ) The differential metabolites ranked as the top 6 based on their VIP values. Statistical significance at p < 0.05 is marked as *, at p < 0.01 is marked as **, and at p < 0.005 is marked as ***.
Figure 3. Analysis of the differential metabolites between pregnant and non-pregnant ewes on 21 days post-AI. (A) Volcano plot for differential metabolites. The red dots represent the upregulated metabolites and the blue dots represent the downregulated metabolites. (B) KEGG classification of the 29 differential metabolites. (C) KEGG pathway enrichment analyses of the 29 differential metabolites. (D) Expression profile of the 29 differential metabolites and their VIP values. (EJ) The differential metabolites ranked as the top 6 based on their VIP values. Statistical significance at p < 0.05 is marked as *, at p < 0.01 is marked as **, and at p < 0.005 is marked as ***.
Ijms 24 14054 g003
Figure 4. Proteomic analysis of sera from ewes. (A) Venn diagrams for the comparison of samples (P represents the pregnant ewes and N represents the non-pregnant ewes). (B) Total protein function annotation based on KEGG pathway enrichment analysis. (C) Subcellular localization analysis of the proteins. (D) Volcano plot of proteins detected in ewes. The red dots represent the upregulated proteins and the blue dots represent the downregulated proteins. (E) Heatmap of all the 14 differential proteins between the two groups. (F) KEGG pathway enrichment analysis for all the 14 differential proteins between the two groups.
Figure 4. Proteomic analysis of sera from ewes. (A) Venn diagrams for the comparison of samples (P represents the pregnant ewes and N represents the non-pregnant ewes). (B) Total protein function annotation based on KEGG pathway enrichment analysis. (C) Subcellular localization analysis of the proteins. (D) Volcano plot of proteins detected in ewes. The red dots represent the upregulated proteins and the blue dots represent the downregulated proteins. (E) Heatmap of all the 14 differential proteins between the two groups. (F) KEGG pathway enrichment analysis for all the 14 differential proteins between the two groups.
Ijms 24 14054 g004
Figure 5. Levels of differentially expressed proteins in the sera of pregnant and non-pregnant ewes on 21 days post-AI. (A) Serpin family A member 6 (SERPINA6). (B) Serum amyloid A protein (SAA). (C) Albumin. (D) Beta-A globin chain β. (E) Afamin (AFM). (F) Ig-like domain-containing protein (ILDP). (G) SMB domain-containing protein. (H) Vitronectin. Statistical significance at p < 0.05 is marked as *.
Figure 5. Levels of differentially expressed proteins in the sera of pregnant and non-pregnant ewes on 21 days post-AI. (A) Serpin family A member 6 (SERPINA6). (B) Serum amyloid A protein (SAA). (C) Albumin. (D) Beta-A globin chain β. (E) Afamin (AFM). (F) Ig-like domain-containing protein (ILDP). (G) SMB domain-containing protein. (H) Vitronectin. Statistical significance at p < 0.05 is marked as *.
Ijms 24 14054 g005
Figure 6. The level of proteins in sera from ewes at 21-, 28- and 33-days post-AI. (AD): The levels of SAA, AFM, SPERIN6 and ILDP were compared at 3 different time points. Statistical significance at p < 0.01 is marked as **, at p < 0.005 is marked as ***, and at p < 0.001 is marked as ****.
Figure 6. The level of proteins in sera from ewes at 21-, 28- and 33-days post-AI. (AD): The levels of SAA, AFM, SPERIN6 and ILDP were compared at 3 different time points. Statistical significance at p < 0.01 is marked as **, at p < 0.005 is marked as ***, and at p < 0.001 is marked as ****.
Ijms 24 14054 g006
Table 1. Differential metabolites between pregnant and non-pregnant ewes on day 21 post-AI.
Table 1. Differential metabolites between pregnant and non-pregnant ewes on day 21 post-AI.
MetaboliteRegulatep-ValueFC (P/N)
Tiglic acidup0.0085031.0083
4-Acetamido-2-aminobutanoic acidup0.00037521.0192
(3R,4R)-3-Amino-1-hydroxy-4-methylpyrrolidin-2-oneup0.030141.0213
Glycylglycylglycineup0.025391.0668
N4-Acetylcytidineup0.024761.0778
3,4-Dimethyl-5-pentyl-2-furanpropanoic acidup0.045411.0456
DL-2-Aminooctanoic acidup0.015691.0753
Capsidiolup0.011631.0407
Polidocanoldown0.021610.9714
4,5-Dihydrovomifoliolup0.049241.0791
PC(18:4(6Z,9Z,12Z,15Z)/20:2(11Z,14Z))up0.016591.01
1-Octen-3-yl glucosidedown0.015610.9823
Indolophenanthridinedown0.043590.8634
D-Cathineup0.0060311.0305
Imexonup0.031421.0665
2-Propylglutaric acidup0.015131.015
1-Oleoylglycerophosphoserinedown0.014330.9535
Thromboxanedown0.015140.9073
(2R,3S)-2-Octyl-3-pentyloxanedown0.03550.9217
PE-NMe(22:2(13Z,16Z)/16:1(9Z))up0.037171.0413
Dodecyl hydrogen sulfatedown0.046240.9879
PS(15:0/24:1(15Z))up0.016131.0313
(+)-gamma-Hydroxy-L-homoarginineup0.02391.0124
2-Mercaptobenzothiazoleup0.044651.1744
Ecgonineup0.034091.1418
MJDBISSN00000004up0.047051.0085
Cyclopentenyl cytosineup0.023621.0617
3-Methyl-2-oxovaleric aciddown0.036530.9819
D-Galactosedown0.026080.9901
Table 2. Differential proteins between pregnant and non-pregnant ewes on day 21 post-AI.
Table 2. Differential proteins between pregnant and non-pregnant ewes on day 21 post-AI.
Accession
Number
DescriptionRelative
Abundance (P)
Relative
Abundance (N)
FC (P/N)p-Value
(P/N)
A0A835ZZL2SERPIN domain-containing protein2825.41793124.84570.90420.04336
A0A452ELD4SMB domain-containing protein 679.40332186.58920.31070.01607
A0A452DX18Complement C3864.38771408.25090.61380.03305
Q1KYZ7Beta-A globin chain β 404.16861532.06830.26380.04558
W5Q9D5Vitronectin 79.6122377.77810.21070.03671
A0A6P3YFZ8Uncharacterized protein 132.6988176.59940.75140.004184
W5NYF4Peptidoglycan recognition protein 2 56.353989.81420.62740.01077
A0A452FXE6Insulin-like growth factor-binding protein 3 51.489666.57620.77340.01531
A0A452F0G8Leukotriene A (4) hydrolase 9.459316.16750.58510.04163
A0A452FW61Ig-like domain-containing protein 0239.55221.00 × 10−50
W5PXI3Afamin069.31521.00 × 10−50
W5PWE9Albumin 19,073.465512,902.10351.4780.03812
A0A835ZMI7Serum amyloid A protein 614.52010424.20
A0A8C2R1B8Serpin family A member 66 5.39730424.20
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhai, Y.; Xia, F.; Shi, L.; Ma, W.; Lv, X.; Sun, W.; Ji, P.; Gao, S.; Machaty, Z.; Liu, G.; et al. Early Pregnancy Markers in the Serum of Ewes Identified via Proteomic and Metabolomic Analyses. Int. J. Mol. Sci. 2023, 24, 14054. https://doi.org/10.3390/ijms241814054

AMA Style

Zhai Y, Xia F, Shi L, Ma W, Lv X, Sun W, Ji P, Gao S, Machaty Z, Liu G, et al. Early Pregnancy Markers in the Serum of Ewes Identified via Proteomic and Metabolomic Analyses. International Journal of Molecular Sciences. 2023; 24(18):14054. https://doi.org/10.3390/ijms241814054

Chicago/Turabian Style

Zhai, Yaying, Fan Xia, Luting Shi, Wenkui Ma, Xiaoyang Lv, Wei Sun, Pengyun Ji, Shuai Gao, Zoltan Machaty, Guoshi Liu, and et al. 2023. "Early Pregnancy Markers in the Serum of Ewes Identified via Proteomic and Metabolomic Analyses" International Journal of Molecular Sciences 24, no. 18: 14054. https://doi.org/10.3390/ijms241814054

APA Style

Zhai, Y., Xia, F., Shi, L., Ma, W., Lv, X., Sun, W., Ji, P., Gao, S., Machaty, Z., Liu, G., & Zhang, L. (2023). Early Pregnancy Markers in the Serum of Ewes Identified via Proteomic and Metabolomic Analyses. International Journal of Molecular Sciences, 24(18), 14054. https://doi.org/10.3390/ijms241814054

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop