Urinary metabolomics reveals the biological characteristics of early pregnancy in pigs
Porcine Health Management volume 8, Article number: 14 (2022)
Embryo implantation in sows is an important event during pregnancy. During this process, blastocysts undergo dramatic morphologic changes, and the endometrium becomes receptive. Studies have shown that developmental changes associated with the crosstalk between peri-implantation embryos and embryo-uterine are driven by various biomolecules secreted by the endometrium and embryos. In sows, changes in the uterus are also reflected in circulating body fluids and urine. Metabolomics reveals the metabolic state of cells, tissues, and organisms. In this study, we collected urine samples from large white sows during the peri-implantation period. The levels of urinary metabolites at different periods were analyzed using ultra-performance liquid chromatography/tandem mass spectrometry (UPLC–MS/MS) analysis techniques.
A total of 32 samples were collected from 8 sows during the estrus period and at each phase of early pregnancy (9, 12, and 15 days of gestation). A total of 530 metabolites were identified with high confidence in all samples. Compared with samples collected during the estrus phase, 269 differential metabolites were found in samples obtained during early pregnancy.
The identified metabolites included lipids and lipid-like molecules, organic acids and their derivatives, organic oxygen compounds, organoheterocyclic compounds, benzenoids, among others. Metabolites, such as choline and pregnanediol-3-glucuronide, play important roles in pregnancy in sows and other animals. These results reveal the metabolic changes in urine of sows during early pregnancy phase. The differential urinary metabolites can be used for assessing peri-implantation status in sows. Understanding these metabolic changes may promote the management of pregnant sows through various interventions such as provision of proper nutrition.
Early pregnancy in pigs is a complex process. The embryonic implantation stage (the first 15 days of pregnancy), is the phase of early pregnancy in sows. During this phase, the embryo and endometrium undergo strong morphological changes [1, 2]. At 8–10 days of gestation, the embryo freely moves throughout the uterine cavity, and gradually differentiates from a filamentous shape into a fetal shape . At around 12 days of pregnancy, tubular embryos elongate rapidly within 1–2 h and become elongated filamentous embryos [4, 5]. This period is also called the peri-implantation period of pregnant sows . By the 15th day of pregnancy, embryonic implantation begins [7, 8]. During this process, pig embryo moves between the uterine horns to provide enough space for attachment of the embryo . During the embryonic attachment stage, steroid hormones are secreted , as well as many different preparations, such as growth factors, prostaglandins and cytokines , which assist in establishing a strong connection between the endometrium and the embryo. In the free blastocyst phase, embryo-maternal communication is mediated through complex signaling factors. The embryo-maternal communication is also fed back in body fluids, including urine as revealed in metabolomic studies . Therefore, understanding the physiological characteristics of pigs during the early pregnancy period is particularly important to improve fertility.
Metabolomics refers to the study of metabolites and their chemical processes in biological samples. While a number of metabolites can be used as sources of metabolic information, mammalian urine provides some advantages as it contains a large number of metabolites  and can be non-invasively collected, thereby minimizing harmful impacts on animals . In addition, as a humoral circulation metabolite, urine can better reflect all biochemical pathways in the body . In recent years, most clinical assays have applied metabolomics in studies on pregnancy and pregnancy-related diseases like diabetes and ectopic pregnancy [15, 16]. Such studies have involved the use of maternal blood, amniotic fluid, and follicular fluid as metabolites [17, 18]. Elsewhere, mammalian urine has been used as a biomarker for studies on gestational diabetes mellitus and abortion . It has been hypothesized that activities such as pregnancy induce metabolomic signatures that can be observed in mammalian urine . However, there is little knowledge regarding urine metabolomic characteristics associated with the early pregnancy phase of mammals, especially pigs.
Pig embryos are too small to be prominent in the peri-implantation period. It is difficult to diagnose the pregnancy status from the appearance, or with the aid of B-ultrasound. Pig embryo implantation in early pregnancy is a multi-complex process that is affected by genetic, nutritional and environmental aspects. The effects of these processes may eventually be reflected in the final metabolites. Therefore, metabolomics analyses of sow urine in early pregnancy is a promising method for identifying the changes in metabolic biological characteristics in early pregnancy. If non-pregnant phenotypes are detected in early pregnancy phases, then, intervention measures can be taken to better improve reproductive efficiency and economic benefits. We used ultra-performance liquid chromatography–mass spectrometry to evaluate the differential changes in urine metabolites of sows in early pregnancy. The selected differential metabolites can be non-invasively used to assess the status of sows during the peri-implantation period. The differential metabolites revealed in this study can be used to model evidence for maternal-embryo communication as well as regulation and control strategies of maternal nutritional requirements for embryonic implantation during early pregnancy in pigs.
Study animals and sample collection
Eight clinically healthy Large White sows (parity 2) from the Wens Foodstuffs Group Co., Ltd (Yunfu, China) pig farm were selected for estrus synchronization treatment after which artificial insemination was performed. Midstream urine samples were collected from sows by the natural urination method on the 0th (DP0, non-pregnant, calculated on the day before insemination), 9th (DP9), 12th (DP12) and 15th (DP15) day of pregnancy. The collected urine samples were quickly placed in liquid nitrogen for metabolite investigation . Immediate sample freezing was necessary to quench any rapid degradation activities, such as oxidation of labile metabolites as well as various enzymatic reactions [21, 22]. The pregnancy status of the sows were confirmed by B-ultrasound at the 4th and 6th weeks.
Metabolite Extraction and quality control
Briefly, 200 μL urine samples were added into a 1.5 mL centrifuge tube using a pipette. An 800 μL acetonitrile:methanol (1:1, v/v) extraction solution was added to the content in the centrifuge tube, vortexed for 30 s at 5℃ and 40 kHz ultrasound for 30 min. Samples were maintained at − 20 °C for 30 min to precipitate the proteins. The supernatants were removed, centrifuged at 13,000×g for 15 min at 4 °C and dried using nitrogen. Then, they were reconstituted with 120 μL acetonitrile:water (1:1, v/v) reconstitution solution, ultrasonically extracted at 5 °C and 40 kHz for 5 min and centrifuged at 13,000×g at 4 °C for 5 min after which the supernatant was used for LC–MS/MS analysis .
To evaluate the stability of the analysis system during the on-boarding process, 20 μL of the supernatant was added to each sample and mixed as the quality control sample (Quality Control, QC). During instrumental analysis, a QC sample was injected for every 8 analytical samples. In data analysis, repeatability of the QC sample was used to investigate the stability of the instrument, and to find variables with large variations in the analysis system to ensure data reliability .
Chromatographic separation of metabolites was performed on a Thermo UHPLC system equipped with an ACQUITY BEH C18 column (100 mm × 2.1 mm I.D., Waters, Milford, USA). The mobile phase included 0.1% formic acid in water (solvent A) and 0.1% formic acid in acetonitrile: isopropanol (1:1, v/v) (solvent B). To balance the system, the solvent gradient varied according to the following conditions: from 0 to 3 min, from 95% (A): 5% (B) to 80% (A): 20% (B); from 0 to 3 min. At 3 to 9 min, 80% (A): 20% (B) to 5% (A): 95% (B); 9 to 13 min, and 5% (A): 95% (B) to 5% (A): 95% (B). From 13 to 13.1 min, 5% (A): 95% (B) to 95% (A): 5% (B), from 13.1 to 16 min, and 95% (A): 5% (B) to 95% (A): 5% (B). Sample injection volume was 2 μL at a flow rate of 0.40 mL/min. Column temperature was kept at 40 ℃. During analysis, these samples were stored at 4 ℃.
Mass spectra data were collected using the Thermo UHPLC-Q Exactive Mass Spectrometer equipped with electrospray ionization (ESI) sources operating in positive and negative ion modes, respectively. Optimum conditions were: Aus gas heater temperature 400℃, sheath gas flow 40 psi; Aus gas flow 30 psi; ion spray voltage floating (ISVF), negative mode-2800 V and positive mode 3500 V, normalized collision energy, MS/MS rolling 20–40–60 V. Data acquisition was performed via the Data Dependent Acquisition (DDA) mode. Detection was performed at the mass range of 70–1050 m/z.
Metabolomics data processing
After UPLC-TOF/MS analysis, raw data was imported into Progenesis QI 2.3 (Waters, USA) for peak detection and comparisons. The preprocessing results generated a data matrix consisting of retention time (RT), mass-to-charge ratio (m/z) values and peak intensity. Online search was done using biochemical databases (Human Metabolome Database (HMDB) (http://www.hmdb.ca/) and Metlin database (https://metlin.scripps.edu/)), with cumulative mass, MS/MS fragment spectrum and isotope ratio differences. Mass spectrometry was used to identify these metabolic characteristics. For metabolites with MS/MS confirmation, only those metabolites with MS/MS fragment scores higher than 30 were considered credibly identified . In the Progenesis QI workflow, we considered confidence scores higher than 50 sufficient for assignment of constituents in the iboga extract .
Multivariate statistical analyses were performed using ropls (Version1.6.2) R package from Bioconductor on Majorbio Cloud Platform (https://cloud.majorbio.com). Principle component analysis (PCA) using an unsupervised method was performed to obtain an overview of the metabolic data. Partial least squares-discriminant analysis (PLS-DA) and orthogonal partial least square discriminant analysis (OPLS-DA) was used for statistical analysis to determine the overall metabolic changes between comparable groups. Before OPLS-DA, all metabolite variables were scaled. Effectiveness of the model was evaluated from the model parameters, R2 and Q2, which provided information for interpretability and predictability of the model and avoided over-fitting risk. Variable Importance in Projection (VIP) was calculated in the OPLS-DA model. p values were determined using the paired Student’s t-test on single dimensional statistical analysis. Statistically significant groups with VIP values more than 1 and p values less than 0.05 were selected. The first 30 different metabolites were selected in the metabolic concentration of each pregnancy period, relative to the estrus period. Various metabolites between the two periods were summarized and mapped to their biochemical pathways through metabolic enrichment and pathway analysis based on database search (KEGG, http://www.genomic.jp/kegg/). Scipy.stats (Python package), using Fisher's exact test, was used to identify statistically significant enrichment pathways. In addition, the identified metabolite expression data (from DP0 to DP15 stage) were normalized to 0, log2(DP9/DP0), log2(DP12/DP0) and log2(DP15/DP0), and the differential metabolites clustered by STEM .
Untargeted metabolic profiling of urine during estrus and early pregnancy
To investigate the metabolic changes in urine during pregnancy in sows, non-targeted metabolomics analysis was performed. Chromatographic separation spectra showed good overlapping of ion chromatograms for each phase samples and quality control samples, in positive ion mode (Fig. 1A) and negative ion mode (Fig. 1B), after analysis by UPLC/Q-TOF MS. After data pre-processing, 5327 features were extracted from the positive ions. After normalization, 4660 features with relative standard deviation (RSD) of < 30% accounted for 87.48% of all QC samples (Additional file 2: Table S1). A total of 8747 features were extracted from negative ions. After normalization, 7560 features with RSD < 30% accounted for 86.43% of all QC samples (Additional file 3: Table S2). The correlation coefficient among the quality control samples was > 0.95. These results imply that there is no significant batch effect on the preprocessed data. Analysis of over 12,220 mass spectra recorded in the total ionization mode revealed 530 compound identifications with high confidence using Progenesis QI (score > 50, Additional file 4: Table S3).
The taxonomic information of metabolites obtained from the HMDB 4.0 database showed that there were diverse metabolites in urine (Fig. 2). With regards to metabolic concentrations (cationic mode and anionic mode), proportions were: Lipids and lipid-like molecules accounted for 28.57% (132 metabolites), organic acids and derivatives accounted for 18.83% (87 metabolites), organic oxygen compounds accounted for 16.45% (76 metabolites), organoheterocyclic compounds accounted for 15.37% (71 metabolites), benzenoids accounted for 8.87% (41 metabolites), phenylpropanoids and polyketides accounted for 6.49% (30 metabolites), nucleosides, nucleotides, and their analogues accounted for 1.95% (9 metabolites), organic nitrogen compounds accounted for 1.52% (7 metabolites), organooxygen compounds accounted for 1.52% (7 metabolites), hydrocarbons accounted for 0.22% (1 metabolite), while lignans, neolignans and related compounds accounted for 0.22% (1 metabolite).
Principle component analysis (PCA) was performed to observe the consistency of experimental samples and quality control samples. The PCA score curve representing QC samples was tightly clustered, compared to other samples (Fig. 3A, B). The QC results showed that our experiment was stable and the obtained data met the conditions for subsequent statistical analyses.
PCA was used to assess sample separation and aggregation between pregnant and non-pregnant sows. The aggregation points indicated a high similarity of the observed variables, whereas discrete points showed significant differences in observed variables (VIP ≥ 1; ratio ≥ 2 or Ratio ≤ 1/2; q ≤ 0.05). In the positive ion mode, the PCA score showed 54.50% variation at PC1 and 10.70% variation at PC2 (Fig. 3A). In the negative ion mode, the PCA score showed 46.80% variation at PC1 and 14.30% variation at PC2 (Fig. 3B). The results showed that urine in estrus and early pregnancy periods exhibited different metabolic characteristics.
The PCA diagram (Fig. 3) showed that metabolites at 9 days of pregnancy (DP9), 12 days of pregnancy (DP12) and 15 days of pregnancy (DP15) were clustered together, while metabolites at estrus period (DP0) were clustered together. There were little differences in urine metabolism of sows in early pregnancy, however, there were significant differences between pregnancy and estrus. Then, we focused on changes in metabolites during pregnancy and estrus periods.
PCA analysis further showed different metabolic characteristics between the urine of pregnant sows and that of estrous sows, indicating that a significant separation of urine metabolites in sows during estrus and early pregnancy. PCA characterized the separation between observation groups in the experimental models and could not identify specific changes between groups. Therefore, to identify specific differences between groups, PLS-DA, a supervised statistical method of discriminant analysis, was employed. The higher the value of PLS-DA model parameters (R2 and Q2), the higher the reliability of the PLS-DA model. Comparative analysis of 9 days of pregnancy and estrus: The R2 of PLS-DA model in the positive ion mode was 0.819 and Q2 was 0.961 (Additional file 1: Figure S1), while in the negative ion mode, R2 of PLS-DA model was 0.809 and Q2 was 0.978 (Additional file 1: Figure S1). Comparative analysis of 12 days of pregnancy and estrus: The R2 of PLS-DA model in the positive ion mode was 0.821 and Q2 was 0.953 (Additional file 1: Figure S2), while in the negative ion mode, the R2 of PLS-DA model was 0.836 and Q2 was 0.964 (Additional file 1: Figure S2). Comparative analysis of 15 days of pregnancy and estrus: The R2 of PLS-DA model in the positive ion mode was 0.651 and Q2 was 0.968 (Additional file 1: Figure S3), whereas in the negative ion mode, the R2 of PLS-DA model was 0.765 and Q2 was 0.983 (Additional file 1: Figure S3). Generally, findings from the PLS-DA model revealed that, both R2 and Q2 were high in early pregnancy. Based on the parameters of the PLS-DA model, model differences were credible, and the PLS-DA model was used for subsequent analyses.
Furthermore, OPLS-DA was used to analyze metabolites in urine during estrus and pregnancy. Comparative analysis of 9 days of pregnancy and estrus in the positive ion mode revealed OPLS-DA parameters as: R_2X = 0.667, R_2Y = 0.996 and Q_2 = 0.949 (Additional file 1: Figure S4) whereas the negative ion mode had OPLS-DA parameters as: R_2X = 0.486, R_2Y = 0.985 and Q_2 = 0.951 (Additional file 1: Figure S4). Comparative analysis of the 12th day of pregnancy and estrus revealed that in the positive ion mode, OPLS-DA parameters were: R_2X = 0.810, R_2Y = 0.992, Q_2 = 0.946 (Additional file 1: Figure S5), while in the negative ion mode, OPLS-DA parameters were: R_2X = 0.765, R_2Y = 0.995, Q_2 = 0.961 (Additional file 1: Figure S5). Comparative analysis of 15 days of pregnancy and the estrus revealed that in the positive ion mode, OPLS-DA parameters were: R_2X = 0.710, R_2Y = 0.984, Q_2 = 0.948 (Additional file 1: Figure S6). On the other hand, negative ion mode showed OPLS-DA parameters as: R_2X = 0.740, R_2Y = 0.997, Q_2 = 0.965 (Additional file 1: Figure S6). In each comparison, 200 permutation tests were used to validate the OPLS-DA model. The intercept of the OPLS-DA model did not reach the over-fitting threshold (R2Y > 0.4, Q2Y > 0.05). The OPLS-DA score chart showed a clear separation of urine between estrus and pregnancy samples. All Q2 values were greater than 0.4, indicating that the OPLS-DA model used was reliable and achieved consistent modeling as well as prediction results.
Expression characteristics of differential metabolites during the peri-implantation period
Different metabolites were screened by multivariate statistical analyses of VIP scores obtained via the OPLS-DA model. Compared to the estrus period, a total of 269 differential metabolites were obtained in early pregnancy. Differential metabolites were defined as: VIP ≥ 1; Ratio ≥ 2 or Ratio ≤ 1/2; q ≤ 0.05. Compared to the estrus stage, there were 186 differential metabolites (Additional file 5: Table S4) at 9 days of gestation, 177 differential metabolites (Additional file 6: Table S5) at 12 days of gestation and 169 differential metabolites (Additional file 7: Table S6) at 15 days of gestation. Significantly different metabolites were used to construct heat maps for unsupervised clustering. Heat maps were used to define the different levels of metabolites in urine during estrus and pregnancy stages. Consistent with OPLS-DA results, significant clustering (Additional file 1: Figure S7) was shown in both early pregnancy and during estrus stages.
The VIP score, based on the OPLS-DA model represented metabolites that were significantly differentially expressed in the two groups, and those variables with VIP scores greater than 1 were considered important in the classification model. Compared to the estrus stage, the top 30 differentially expressed metabolites in the three phases of early pregnancy were selected for analyses (Fig. 4A–C). Nine metabolites overlapped during early pregnancy and estrus periods (Fig. 4D). Of the nine, one metabolite, Pregnanediol-3-glucuronide (PdG), was up-regulated in pregnancy compared to estrus period. Eight metabolites, including N6-carbamoyl-L-threonyladenosine, salicyluric acid, prostaglandin E3, marmelolactone A, caryophyllene epoxide, 3,4,5-trihydroxy-6-[3-(3-phenylpropanoyl)phenoxy]oxane-2-carboxylic acid, glucosyl 6-hydroxy-2,6-dimethyl-2E,7-octadienoate, and 6-(2,4-dihydroxybenzoyloxy)-3,4,5-trihydroxyoxane-2-carboxylic acid were down-regulated. The significant changes among these metabolites during early pregnancy may be associated with regulation of blastocyst implantation.
Metabolic pathways for the metabolites
Variable metabolites from DP9 and DP0, DP12 and DP0, DP15 and DP0, were subjected to further biochemical pathway analyses. Functional enrichment showed a significant shift (p < 0.05) in expression of molecules associated with amino acid metabolism, lipid metabolism and organic acid production during early pregnancy (Fig. 5). Compared to the estrus stage, at 9 days of pregnancy, the major urinary metabolite enrichment pathways of pregnancy were bile secretion, phenylalanine metabolism, drug metabolism—cytochrome P450, purine metabolism, choline metabolism in cancer, cholinergic synapse, steroid hormone biosynthesis, and sphingolipid signaling pathway (Fig. 5A). Compared to estrus period, at 12 days of pregnancy, the major urinary metabolite enrichment pathways were purine metabolism, alanine, aspartate and glutamate metabolism, apoptosis, glycosaminoglycan biosynthesis—chondroitin sulfate/dermatan sulfate, tyrosine metabolism, tryptophan metabolism, bile secretion, amino sugar as well as nucleotide sugar metabolism, and sphingolipid signaling pathway (Fig. 5B). Relative to the estrus period, at 15 days of pregnancy, the main urinary metabolite enrichment pathways were phenylalanine, tyrosine and tryptophan biosynthesis, glycosaminoglycan biosynthesis—chondroitin sulfate/dermatan sulfate, phenylalanine metabolism, tyrosine metabolism, drug metabolism—cytochrome P450, and purine metabolism (Fig. 5C). In the urinary metabolic pathway, in early pregnancy, a shift in amino acid metabolism was observed in biosynthesis and metabolism of phenylalanine, tyrosine and tryptophan. However, the amino acid metabolic pathway also underwent slight changes at different stages, with phenylalanine predominating on day 9 of pregnancy and tyrosine on days 12 and 15 of pregnancy. In terms of lipid metabolism, in early pregnancy, 17-beta-Estradiol-3-glucuronide, 2-Methoxy-estradiol-17b 3-glucuronide, 4-oxo-Retinoic acid, 9,10,13-TriHOME, cholestane-3,7,12,25-tetrol-3-glucuronide, pregnanediol-3-glucuronide, SM(d18:0/16:1(9Z)) and traumatic acid were enriched.
Furthermore, to characterize the expression profiles of 530 metabolites, expression data (DP0 to DP9, DP9 to DP12, and DP12 to DP15) were normalized to 0, log2 (DP9/DP0), log2 (DP12/DP0), and log2 (DP15/ DP0), as analyzed by the Short Time-series Expression Miner (STEM), the 530 metabolites were clustered into eight profiles (0–7) (Fig. 6, Additional file 8: Table S7). The metabolite, PdG, was clustered in profile 7, and its expression levels continued to increase during the peri-implantation period. The eight metabolites that were significantly down-regulated during the peri-implantation stage were clustered in profile 0 (5 metabolites: salicyluric acid, prostaglandin E3, marmelolactone A, caryophyllene epoxide, and glucosyl 6-hydroxy-2,6-dimethyl-2E,7-octadienoate), profile 1 (2 metabolites: 3,4,5-trihydroxy-6-[3-(3-phenylpropanoyl)phenoxy]oxane-2-carboxylic acid, 6-(2,4-dihydroxybenzoyloxy)-3,4,5-trihydroxyoxane-2-carboxylic acid), and profile 4 (1 metabolite: N6-carbamoyl-L-threonyladenosine) (Additional file 8: Table S7). These metabolic changes in early pregnancy reflect changes in nutritional needs of the embryo and uterus, as well as the metabolic state of the mother.
We used the MS-based approach and reliable metabolite identification to detect urine metabolites through UPLC Q-TOF/MS analysis. Urine was obtained from eight sows in estrus and early pregnancy periods and examined to determine changes in metabolite levels of different pathways associated with early pregnancy. A total of 269 differential metabolites were identified in the early pregnancy period while compared to estrus period, a total of 186 different metabolites were found at 9 days of gestation, 177 different metabolites at 12 days of gestation, and 169 different metabolites at 15 days of gestation. These differential metabolites represent changes in various metabolic processes, such as amino acid and lipid metabolism as well as organic acid contents. One of the metabolites for choline is vitamin B. Choline is an essential nutrient for mammals, an important component of phospholipids in cell membranes and one of the raw materials for the body's synthesis of glycosaminoglycans . The metabolite is involved in lipid metabolism, brain development and fetal development [29, 30]. The choline derivative, phosphatidylcholine, is a major constituent of all cell membranes and is required for the biosynthesis of lipoproteins, including very low-density lipoproteins (VLDLs), which facilitate hepatic exports of lipids . Choline regulates growth and metabolism in animals [32, 33]. During embryogenesis, choline is essential for the growth of cellular structures . Choline kinase is detected at all stages of embryonic development, and its activity increases during development and qualitatively parallels the extent of phosphorylcholine formation in intact embryos . Studies have shown that choline chloride alters development, specific gene expression and DNA methylation in bovine embryos cultured in vitro in a concentration-dependent manner. The role of choline chloride in promoting DNA methylation in the blastocyst is not only as a methyl donor, but may also be due to the role of choline itself or choline metabolites in a variety of functions including cell signalling . Transcriptional coding for enzymes involved in the synthesis of these choline derivatives is present during the peri-implantation phase in mammals [36, 37], and there is evidence that these pathways are important for preimplantation development. In mice, addition of acetylcholine to the culture medium allows embryos to develop from the 2-cell block zone to the blastocyst stage , and knockdown of the choline kinase gene (the first enzyme in the phosphatidylcholine synthesis pathway) results in the death of mouse embryos between 3.5 and 7.5 days of development . In addition, studies have confirmed that before the implantation, compared with the blank group, the implantation rate of the choline-added rat group was higher than that of the blank group . During pregnancy, choline ingested by the mother is transported through the placenta to the fetus . Adequate choline intake during pregnancy is essential for proper fetal development . Supply of maternal choline to the fetus plays an important role in fetal brain development, membrane biosynthesis, and neurotransmission . In a research involving rat models, prenatal choline supplementation preserved a sustained memory state of rats and protected offsprings from developing memory defects . In pigs, dietary choline deficiency can cause significant changes in plasma choline metabolites at the end of lactation. These changes in concentrations of nutrients such as choline may affect the development of early newborns . Although examination was not performed on natural pregnant sows, we postulate that this is one of the explanations for the increase in urine choline levels in sows during pregnancy. Pregnancy and lactation are the periods when maternal reserves of choline are depleted. At the same time, the availability of choline for normal development of the brain is critical because it influences lifelong memory enhancement . The broad role and diverse involvement of choline in cells suggests the importance of its pathway to regulate embryonic development and prepare for implantation.
Of the various differentially expressed metabolites, pregnanediol-3-glucuronide (PdG) was highly expressed in all the three phases of early pregnancy. The PdG metabolite plays an important role in regulation of mammalian pregnancy. Volkery et al. investigated the levels of PdG in plasma, saliva, milk, and urine of alpaca without pregnancy and throughout pregnancy. The concentrations of PdG in plasma, milk, and urine were significantly higher than those in non-pregnant alpaca . The ovaries of mammals are luteinized once they are pregnant, and the corpus luteum dissolves at the end of pregnancy . Ecochard et al. characterized the variabilities in hormonal profiles during the luteal phase in normal cycles. They found that PdG levels increased with early lutealization, high levels of PdG were maintained in the mid-luteal phase which decreased gradually to a low level in the late luteal phase (luteolysis) during the non-pregnant state . The occurrence of normal and low luteal PdG levels are potential markers of luteal phase abnormality . In a study by Rene Leiva et al., assessment of urine PdG levels was used as an auxiliary means to confirm the effectiveness of ovulation . Its concentration was found to be positively correlated with the number of pregnant embryos in feces of female golden monkeys , while its accuracy and sensitivity as a pregnancy test marker were validated in monitoring pregnancies in Mus musculus . In this study, we found that PdG was highly expressed in the urine of sows during early pregnancy (Fig. 6, Profile 7), in tandem with levels in various body fluids during early pregnancy in other mammals. PdG, the main terminal metabolite of progesterone, plays an important role in physiological processes like female menstrual cycle, pregnancy embryogenesis and maternal immune responses in some mammals . In this study, sow urinary PdG levels showed a continuous increasing trend during the peri-implantation period, which may be associated with embryonic implantation and development, as the need for pregnancy maintenance increased.
Pathways associated with metabolites identified at each pregnancy phase were analysed. The results showed that tyrosine metabolism, phenylalanine metabolism, and tryptophan metabolism play an important role in pregnancy. Liu et al. performed a longitudinal analysis of urine from 50 healthy pregnant women and found own-regulation of tyrosine metabolism during pregnancy . During pregnancy, tyrosine metabolism can be formed by phenylalanine hydroxylation  and these two pathways are also reflected in this study (Fig. 5). Maternal tyrosine intake during gestation plays a pre-regulatory role in the development of dopaminergic systems in offsprings after birth . McBride et al. determined the reference ranges for phenylalanine and tyrosine levels in healthy pregnancy. They reported that phenylalanine and tyrosine levels rapidly declined in early pregnancy and remained relatively stable in mid- and late pregnancy . Tryptophan is an indispensable amino acid in mammalian protein synthesis. Tryptophan metabolism within different tissues is associated with numerous physiological functions. The liver regulates tryptophan homeostasis by degrading excess tryptophan. Tryptophan degradation into kynurenine by immune cells plays a crucial role in regulation of immune responses during infections, inflammation, and pregnancy . Tryptophan enhances maternal protein synthesis, fetal growth and development and improves kynurenines inhibition of fetal rejection during gestation in mammals . Tryptophan metabolism during mammalian gestation can be achieved through multiple pathways . For instance, through the kynurenine pathway, where tryptophan is converted to kynurenine by the action of indoleamine 2-3 deoxygenase and tryptophan-2,3-double-oxygenase. Then, kynurenine is converted to kynurenic acid by the actions of kynurenine aminotransferase (KAT). In cases of depression, kynurenic acid may have neuro-protective effects . The tryptophan metabolic pathway also generates 3-hydroxyanthranilic acid (HA), an effective antioxidant, intra-cellularly via kynurenic acid, in the extrahepatic tissue . Tryptophan and its metabolites are effective in scavenging free radicals (including reactive oxygen and reactive chlorine), therefore, they can be targeted in the placenta as potent antioxidants . In this study, compared to the estrus period, various differentially expressed metabolites in early pregnancy were found to be enriched in the amino acid metabolism pathway. This change may be required for embryonic development or immune regulation as in other animals. We did not exhaust differential metabolites and related pathways (including some down-regulated metabolites), therefore, additional studies should be conducted in this area. Furthermore, there are some limitations of this study, one being whether the relevant metabolite levels correlate with the number of successful embryos implanted to screen sows with different fecundity levels, and the other being whether the urinary metabolites of sows without implantation after insemination also have relevant changes is an area that deserves further research. In addition, in commercial farms, the collection of urine samples immediately put into liquid nitrogen preservation may not working in daily practise. It can be taken to add 0.42% sodium azide preservative at the final concentration within 2 h of urine sample collection, and store it at − 80 °C.
Urinary metabolomics of sows during estrus and early pregnancy were determined by UHPLC. A total of 269 urinary metabolites, including lipids and organic acids were identified in early pregnancy. In particular, PdG levels were consistently highly expressed at each pregnancy phase of early pregnancy. As the main terminal metabolite of progesterone, its high expressions are necessary to maintain pregnancy. To the best of our knowledge, this is the first investigation of urinary metabolic changes in sows during early pregnancy. Our findings elucidate on maternal urinary changes caused by intrauterine embryonic implantation during early pregnancy in sows. Moreover, they form a basis for evaluating the molecular mechanisms involved in maternal–fetal communication during embryonic implantation and nutritional requirements of early pregnancy pigs.
Ultra-performance liquid chromatography–tandem mass spectrometry
Data dependent acquisition
Day 0 of pregnancy
Day 9 of pregnancy
Day 12 of pregnancy
Day 15 of pregnancy
Principle component analysis
Partial least squares-discriminant analysis
Orthogonal partial least square-discriminant analysis
Variable importance in projection
Quality control sample
Short time-series expression miner
Keys JL, King GJ. Microscopic examination of porcine conceptus-maternal interface between days 10 and 19 of pregnancy. Am J Anat. 1990;188(3):221–38. https://doi.org/10.1002/aja.1001880302.
Bazer FW, Johnson GA. Pig blastocyst-uterine interactions. Differentiation. 2014;87(1–2):52–65. https://doi.org/10.1016/j.diff.2013.11.005.
Anderson LL. Growth, protein content and distribution of early pig embryos. Anat Rec. 1978;190(1):143–53. https://doi.org/10.1002/ar.1091900112.
Le Blomberg A, Garrett WM, Guillomot M, Miles JR, Sonstegard TS, Van Tassell CP, et al. Transcriptome profiling of the tubular porcine conceptus identifies the differential regulation of growth and developmentally associated genes. Mol Reprod Dev. 2006;73(12):1491–502. https://doi.org/10.1002/mrd.20503.
Waclawik A, Kaczmarek MM, Blitek A, Kaczynski P, Ziecik AJ. Embryo-maternal dialogue during pregnancy establishment and implantation in the pig. Mol Reprod Dev. 2017;84(9):842–55. https://doi.org/10.1002/mrd.22835.
Li W, Xi Y, Xue S, Wang Y, Wu L, Liu H, et al. Sequence analysis of microRNAs during pre-implantation between Meishan and Yorkshire pigs. Gene. 2018;646:20–7. https://doi.org/10.1016/j.gene.2017.12.046.
Franczak A, Wojciechowicz B, Kotwica G. Transcriptomic analysis of the porcine endometrium during early pregnancy and the estrous cycle. Reprod Biol. 2013;13(3):229–37. https://doi.org/10.1016/j.repbio.2013.07.001.
Hong L, Liu R, Qiao X, Wang X, Wang S, Li J, et al. Differential microRNA expression in porcine endometrium involved in remodeling and angiogenesis that contributes to embryonic implantation. Front Genet. 2019;10:661. https://doi.org/10.3389/fgene.2019.00661.
Geisert RD, Schmitt RAM. Early embryonic survival in the pig: Can it be improved? J Ani Sci. 2002;80(E-suppl_1):E54–65. https://doi.org/10.2527/animalsci2002.0021881200800ES10009x.
Spencer TE, Bazer FW. Conceptus signals for establishment and maintenance of pregnancy. Anim Reprod Sci. 2004;82–83:537–50. https://doi.org/10.1016/j.anireprosci.2004.04.014.
Maner DF, Saffan BD, Wiggins RA, Thompson JD, Preedy JRK. Interrelationship of estrogen concentrations in the maternal circulation, fetal circulation and maternal urine in late pregnancy. J Clin Endocrinol Metab. 1963;23(5):445–58. https://doi.org/10.1210/jcem-23-5-445.
Blazenovic I, Kind T, Sa MR, Ji J, Vaniya A, Wancewicz B, et al. Structure annotation of all mass spectra in untargeted metabolomics. Anal Chem. 2019;91(3):2155–62. https://doi.org/10.1021/acs.analchem.8b04698.
Wald C. Diagnostics: a flow of information. Nature. 2017;551(7679):S48–50. https://doi.org/10.1038/551S48a.
Khamis MM, Adamko DJ, El-Aneed A. Mass spectrometric based approaches in urine metabolomics and biomarker discovery. Mass Spectrom Rev. 2017;36(2):115–34. https://doi.org/10.1002/mas.21455.
Bardanzellu F, Fanos V. The choice of amniotic fluid in metabolomics for the monitoring of fetus health—update. Expert Rev Proteom. 2019;16(6):487–99. https://doi.org/10.1080/14789450.2019.1615892.
Liu X, Wang X, Sun H, Guo Z, Liu X, Yuan T, et al. Urinary metabolic variation analysis during pregnancy and application in gestational diabetes mellitus and spontaneous abortion biomarker discovery. Sci Rep. 2019;9(1):2605. https://doi.org/10.1038/s41598-019-39259-2.
Orczyk-Pawilowicz M, Jawien E, Deja S, Hirnle L, Zabek A, Mlynarz P. Metabolomics of human amniotic fluid and maternal plasma during normal pregnancy. PLoS ONE. 2016;11(4): e0152740. https://doi.org/10.1371/journal.pone.0152740.
Karaer A, Tuncay G, Mumcu A, Dogan B. Metabolomics analysis of follicular fluid in women with ovarian endometriosis undergoing in vitro fertilization. Syst Biol Reprod Med. 2019;65(1):39–47. https://doi.org/10.1080/19396368.2018.1478469.
Wu JQ, Gao YH. Physiological conditions can be reflected in human urine proteome and metabolome. Expert Rev Proteom. 2015;12(6):623–36.
Fernández-Peralbo MA, Luque de Castro MD. Preparation of urine samples prior to targeted or untargeted metabolomics mass-spectrometry analysis. TrAC Trends Anal Chem. 2012;41:75–85. https://doi.org/10.1016/j.trac.2012.08.011.
Bernini P, Bertini I, Luchinat C, Nincheri P, Staderini S, Turano P. Standard operating procedures for pre-analytical handling of blood and urine for metabolomic studies and biobanks. J Biomol NMR. 2011;49(3–4):231–43. https://doi.org/10.1007/s10858-011-9489-1.
Ryan D, Robards K, Prenzler PD, Kendall M. Recent and potential developments in the analysis of urine: a review. Anal Chim Acta. 2011;684(1–2):8–20. https://doi.org/10.1016/j.aca.2010.10.035.
Saude EJ, Sykes BD. Urine stability for metabolomic studies: effects of preparation and storage. Metabolomics. 2007;3(1):19–27. https://doi.org/10.1007/s11306-006-0042-2.
Want EJ, Wilson ID, Gika H, Theodoridis G, Plumb RS, Shockcor J, et al. Global metabolic profiling procedures for urine using UPLC–MS. Nat Protoc. 2010;5(6):1005–18. https://doi.org/10.1038/nprot.2010.50.
Li Y, Wang D, Zeng C, Liu Y, Huang G, Mei Z. Salivary metabolomics profile of patients with recurrent aphthous ulcer as revealed by liquid chromatography-tandem mass spectrometry. J Int Med Res. 2018;46(3):1052–62. https://doi.org/10.1177/0300060517745388.
Bading-Taika B, Akinyeke T, Magana AA, Choi J, Ouanesisouk M, Torres ERS, et al. Phytochemical characterization of Tabernanthe iboga root bark and its effects on dysfunctional metabolism and cognitive performance in high-fat-fed C57BL/6J mice. J Food Bioact. 2018;3:111–23. https://doi.org/10.31665/JFB.2018.3154.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006;7:191. https://doi.org/10.1186/1471-2105-7-191.
Blusztajn JK, Slack BE, Mellott TJ. Neuroprotective actions of dietary choline. Nutrients. 2017. https://doi.org/10.3390/nu9080815.
Zeisel SH. Choline: critical role during fetal development and dietary requirements in adults. Annu Rev Nutr. 2006;26:229–50. https://doi.org/10.1146/annurev.nutr.26.061505.111156.
Niculescu MD, Craciunescu CN, Zeisel SH. Dietary choline deficiency alters global and gene-specific DNA methylation in the developing hippocampus of mouse fetal brains. FASEB J. 2006;20(1):43–9. https://doi.org/10.1096/fj.05-4707com.
Korsmo HW, Jiang X, Caudill MA. Choline: exploring the growing science on its benefits for moms and babies. Nutrients. 2019. https://doi.org/10.3390/nu11081823.
de Veth MJ, Artegoitia VM, Campagna SR, Lapierre H, Harte F, Girard CL. Choline absorption and evaluation of bioavailability markers when supplementing choline to lactating dairy cows. J Dairy Sci. 2016;99(12):9732–44. https://doi.org/10.3168/jds.2016-11382.
Tayebati SK, Amenta F. Choline-containing phospholipids: relevance to brain functional pathways. Clin Chem Lab Med. 2013;51(3):513–21. https://doi.org/10.1515/cclm-2012-0559.
Pratt HP. Phospholipid synthesis in the preimplantation mouse embryo. J Reprod Fertil. 1980;58(1):237–48. https://doi.org/10.1530/jrf.0.0580237.
Estrada-Cortes E, Negron-Perez VM, Tribulo P, Zenobi MG, Staples CR, Hansen PJ. Effects of choline on the phenotype of the cultured bovine preimplantation embryo. J Dairy Sci. 2020;103(11):10784–96. https://doi.org/10.3168/jds.2020-18598.
Wang QT, Piotrowska K, Ciemerych MA, Milenkovic L, Scott MP, Davis RW, et al. A genome-wide study of gene activity reveals developmental signaling pathways in the preimplantation mouse embryo. Dev Cell. 2004;6(1):133–44. https://doi.org/10.1016/s1534-5807(03)00404-0.
Ikeda S, Namekawa T, Sugimoto M, Kume S. Expression of methylation pathway enzymes in bovine oocytes and preimplantation embryos. J Exp Zool A Ecol Genet Physiol. 2010;313(3):129–36. https://doi.org/10.1002/jez.581.
Cho SK, Yoon SY, Hur CG, Yang HY, Choe C, Kim EJ, et al. Acetylcholine rescues two-cell block through activation of IP3 receptors and Ca2+/calmodulin-dependent kinase II in an ICR mouse strain. Pflugers Arch. 2009;458(6):1125–36. https://doi.org/10.1007/s00424-009-0686-7.
Wu G, Aoyama C, Young SG, Vance DE. Early embryonic lethality caused by disruption of the gene for choline kinase alpha, the first enzyme in phosphatidylcholine biosynthesis. J Biol Chem. 2008;283(3):1456–62. https://doi.org/10.1074/jbc.M708766200.
Moore NP, Wahl M, Schneider S. Implantation loss induced by ethanolamine in the rat is ameliorated by a choline-supplemented diet. Reprod Toxicol. 2018;78:102–10. https://doi.org/10.1016/j.reprotox.2018.04.005.
Baumgartner HK, Trinder KM, Galimanis CE, Post A, Phang T, Ross RG, et al. Characterization of choline transporters in the human placenta over gestation. Placenta. 2015;36(12):1362–9. https://doi.org/10.1016/j.placenta.2015.10.001.
Drews K. Folate metabolism–epigenetic role of choline and vitamin B12 during pregnancy. Ginekol Pol. 2015;86(12):940–6. https://doi.org/10.17772/gp/60561.
Yang Y, Liu Z, Cermak JM, Tandon P, Sarkisian MR, Stafstrom CE, et al. Protective effects of prenatal choline supplementation on seizure-induced memory impairment. J Neurosci. 2000;20(22):RC109.
Mudd AT, Alexander LS, Johnson SK, Getty CM, Malysheva OV, Caudill MA, et al. Perinatal DIEtary choline deficiency in sows influences concentrations of choline metabolites, fatty acids, and amino acids in milk throughout lactation. J Nutr. 2016;146(11):2216–23. https://doi.org/10.3945/jn.116.238832.
Zeisel SH. Choline: needed for normal development of memory. J Am Coll Nutr. 2000;19(5 Suppl):528S-531S. https://doi.org/10.1080/07315724.2000.10718976.
Volkery J, Gottschalk J, Sobiraj A, Wittek T, Einspanier A. Progesterone, pregnanediol-3-glucuronide, relaxin and oestrone sulphate concentrations in saliva, milk and urine of female alpacas (Vicugna pacos) and their application in pregnancy diagnosis. Vet Rec. 2012;171(8):195. https://doi.org/10.1136/vr.100393.
Oliver R, Pillarisetty LS. Anatomy, abdomen and pelvis, ovary corpus luteum. In: StatPearls. Treasure Island (FL); 2021.
Ecochard R, Bouchard T, Leiva R, Abdulla S, Dupuis O, Duterque O, et al. Characterization of hormonal profiles during the luteal phase in regularly menstruating women. Fertil Steril. 2017;108(1):175-182 e1. https://doi.org/10.1016/j.fertnstert.2017.05.012.
Abdulla SH, Bouchard TP, Leiva RA, Boyle P, Iwaz J, Ecochard R. Hormonal predictors of abnormal luteal phases in normally cycling women. Front Public Health. 2018;6:144. https://doi.org/10.3389/fpubh.2018.00144.
Leiva R, McNamara-Kilian M, Niezgoda H, Ecochard R, Bouchard T. Pilot observational prospective cohort study on the use of a novel home-based urinary pregnanediol 3-glucuronide (PDG) test to confirm ovulation when used as adjunct to fertility awareness methods (FAMs) stage 1. BMJ Open. 2019;9(5): e028496. https://doi.org/10.1136/bmjopen-2018-028496.
Muren, Kusuda S, Doi O, Naito H, Hashikawa H. Puberty, ovarian cycle, pregnancy, and postpartum ovulation in captive Sichuan golden monkeys (Rhinopithecus roxellana) based on changes in urinary and fecal gonadal steroid metabolites. Theriogenology. 2017;87:179–86. https://doi.org/10.1016/j.theriogenology.2016.08.020.
Mithileshwari C, Srivastava T, Kumar V, Kumar A, Umapathy G. Non-invasive assessment of fecal progestagens and pregnancy detection in Himalayan musk deer (Moschus chrysogaster). Theriogenology. 2016;85(2):216–23. https://doi.org/10.1016/j.theriogenology.2015.09.009.
Chen J, Hou H, Chen H, Luo Y, He Y, Wang H, et al. Quantitation of serum pregnanediol-3-glucuronide level in different endocrine statuses by LC-MS/MS. J Pharm Biomed Anal. 2020;184: 113171. https://doi.org/10.1016/j.jpba.2020.113171.
Whittaker PG, Lee CH, Cooper BG, Taylor R. Evaluation of phenylalanine and tyrosine metabolism in late human pregnancy. Metabolism. 1999;48(7):849–52. https://doi.org/10.1016/s0026-0495(99)90217-2.
Santana C, Martin L, Rodriguez Diaz M. Tyrosine ingestion during rat pregnancy alters postnatal development of dopaminergic neurons in the offspring. Brain Res. 1994;635(1–2):96–102. https://doi.org/10.1016/0006-8993(94)91427-3.
McBride KL, Pluciniczak J, Rhyand T, Bartholomew D. Phenylalanine and tyrosine measurements across gestation by tandem mass spectrometer on dried blood spot cards from normal pregnant women. Genet Med. 2019;21(8):1821–6. https://doi.org/10.1038/s41436-018-0407-8.
Le Floc’h N, Otten W, Merlot E. Tryptophan metabolism, from nutrition to potential therapeutic applications. Amino Acids. 2011;41(5):1195–205. https://doi.org/10.1007/s00726-010-0752-7.
Badawy AA. Tryptophan metabolism, disposition and utilization in pregnancy. 2015. Biosci Rep. https://doi.org/10.1042/BSR20150197.
Badawy AA. Kynurenine pathway of tryptophan metabolism: regulatory and functional aspects. Int J Tryptophan Res. 2017;10:1178646917691938. https://doi.org/10.1177/1178646917691938.
Savitz J, Drevets WC, Smith CM, Victor TA, Wurfel BE, Bellgowan PS, et al. Putative neuroprotective and neurotoxic kynurenine pathway metabolites are associated with hippocampal and amygdalar volumes in subjects with major depressive disorder. Neuropsychopharmacology. 2015;40(2):463–71. https://doi.org/10.1038/npp.2014.194.
Pae HO, Oh GS, Lee BS, Rim JS, Kim YM, Chung HT. 3-Hydroxyanthranilic acid, one of L-tryptophan metabolites, inhibits monocyte chemoattractant protein-1 secretion and vascular cell adhesion molecule-1 expression via heme oxygenase-1 induction in human umbilical vein endothelial cells. Atherosclerosis. 2006;187(2):274–84. https://doi.org/10.1016/j.atherosclerosis.2005.09.010.
Xu K, Liu G, Fu C. The tryptophan pathway targeting antioxidant capacity in the placenta. Oxid Med Cell Longev. 2018;2018:1054797. https://doi.org/10.1155/2018/1054797.
Yang QJ, Zhao JR, Hao J, Li B, Huo Y, Han YL, et al. Serum and urine metabolomics study reveals a distinct diagnostic model for cancer cachexia. J Cachexia Sarcopenia Muscle. 2018;9(1):71–85. https://doi.org/10.1002/jcsm.12246.
Haug K, Cochrane K, Nainala VC, Williams M, Chang J, Jayaseelan KV, et al. MetaboLights: a resource evolving in response to the needs of its scientific community. Nucleic Acids Res. 2020;48(D1):D440–4. https://doi.org/10.1093/nar/gkz1019.
We thank the Guangdong Wen’s pig breeding company for their assistance in this research.
This work was supported by the Guangdong Provincial Promotion Project on Livestock Industry Transformation and Utilization of Livestock Waste (2018XMZC), the National Natural Science Foundation (Grant No. 31802033), the Department of Science and Technology of Guangdong Province (Grant No. 2019A1515011835), and the Science & Technology Planning Project of Guangzhou (Grant No. 201904010434).
Ethical approval and consent to participate
All procedures involving animals were approved by the Ethical Committee of the Laboratory Animal Center of South China Agricultural University (Approval No. SYXK-2014-0136).
Consent for publication
The authors declare that this research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. All the positive ion features.
Table S2. All the negative ion features.
Table S3. Metabolites identified at high confidence level.
Table S4. MetabsetVip_DP9_vs_DP0.
Table S5. MetabsetVip_DP12_vs_DP0.
Table S6. MetabsetVip_DP15_vs_DP0.
Table S7. Expression profiles(0-7).
Supplementary Figures (Figure S1 to S7).
About this article
Cite this article
Zhou, C., Cai, G., Meng, F. et al. Urinary metabolomics reveals the biological characteristics of early pregnancy in pigs. Porc Health Manag 8, 14 (2022). https://doi.org/10.1186/s40813-022-00256-z