Next Article in Journal
Double-Stranded RNA Enhances Matrix Metalloproteinase-1 and -13 Expressions through TLR3-Dependent Activation of Transglutaminase 2 in Dermal Fibroblasts
Next Article in Special Issue
Scd1 Deficiency in Early Embryos Affects Blastocyst ICM Formation through RPs-Mdm2-p53 Pathway
Previous Article in Journal
A New Perspective on Thyroid Hormones: Crosstalk with Reproductive Hormones in Females
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 5
10.3390/ijms23052706
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

LC-MS Analysis Revealed the Significantly Different Metabolic Profiles in Spent Culture Media of Human Embryos with Distinct Morphology, Karyotype and Implantation Outcomes

by
Chupalav Eldarov
1,
Alina Gamisonia
1,
Vitaliy Chagovets
1,
Luiza Ibragimova
1,
Svetlana Yarigina
1,
Veronika Smolnikova
1,
Elena Kalinina
1,
Nataliya Makarova
1,
Victor Zgoda
2,
Gennady Sukhikh
1,3 and
Mikhail Bobrov
1,*
1
National Medical Research Center for Obstetrics, Gynecology and Perinatology Named after Academician V.I. Kulakov, 4 Oparina Str., 117997 Moscow, Russia
2
Orekhovich Institute of Biomedical Chemistry, 10 Pogodinskaya Str., 119121 Moscow, Russia
3
Department of Obstetrics, Gynecology, Perinatology and Reproductology, Institute of Professional Education, Federal State Autonomous Educational Institution of Higher Education I.M. Sechenov First Moscow State Medical University of the Ministry of Health of the Russian Federation (Sechenov University), 8/2 Trubetskaya Str., 119991 Moscow, Russia
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(5), 2706; https://doi.org/10.3390/ijms23052706
Submission received: 20 January 2022 / Revised: 17 February 2022 / Accepted: 25 February 2022 / Published: 28 February 2022
(This article belongs to the Special Issue Frontiers in the Development of Germ Cells and Embryos)
Browse Figures
Versions Notes

Abstract

:
In this study we evaluated possible differences in metabolomic profiles of spent embryo culture media (SECM) of human embryos with distinct morphology, karyotype, and implantation outcomes. A total of 153 samples from embryos of patients undergoing in vitro fertilization (IVF) programs were collected and analyzed by HPLC-MS. Metabolomic profiling and statistical analysis revealed clear clustering of day five SECM from embryos with different morphological classes and karyotype. Profiling of day five SECM from embryos with different implantation outcomes showed 241 significantly changed molecular ions in SECM of successfully implanted embryos. Separate analysis of paired SECM samples on days three and five revealed 46 and 29 molecular signatures respectively, significantly differing in culture media of embryos with a successful outcome. Pathway enrichment analysis suggests certain amino acids, vitamins, and lipid metabolic pathways to be crucial for embryo implantation. Differences between embryos with distinct implantation potential are detectable on the third and fifth day of cultivation that may allow the application of culture medium analysis in different transfer protocols for both fresh and cryopreserved embryos. A combination of traditional morphological criteria with metabolic profiling of SECM may increase implantation rates in assisted reproductive technology programs as well as improve our knowledge of the human embryo metabolism in the early stages of development.

1. Introduction

The application of the assisted reproductive technologies (ART) has been constantly improving over the last four decades since the first IVF-assisted birth in 1978 [1]. While the number of ART cycles steadily increase worldwide, live birth rates vary significantly over different regions—from 12–51% for fresh embryos to 18–57% for cryopreserved ones [2].
One of the ways to increase ART efficiency is double embryo transfer [3,4,5,6], but this approach has well-known risks of multiple pregnancies, associated with different complications for infant and mother health [7,8], which have tendency to increase with age [9,10]. Over the past two decades, there has been a shift towards the elective single embryo transfer (eSET) protocols [11]. While this approach mitigates risks associated with multi-pregnancy, it requires developing an accurate strategy of embryo quality assessment, selection, and transfer. Initially, the question was raised of when the embryo should be transferred. Centers of Disease Control and Prevention (CDC) in the USA reported that in 2013, most of the embryos were transferred on day 3 at the cleavage stage or day 5 at the blastocyst stage, but transfers on day 1, 2, 4, and 6 have been reported as well [12]. A meta-analysis including 1654 patients has shown significantly higher clinical pregnancy and live birth rates among embryos transferred at the blastocyst stage compared to those transferred at cleavage state. The 5th-day blastocysts have better synchronization with an IVF stimulated endometrium and the high implantation potential [13]. One of the most obvious and widely used methods of embryo quality assessment utilizes this correlation and is based on morphological criteria, which considers blastocoel size, form, and shape of trophectoderm and inner cell mass for the evaluation of embryo development and viability [14]. Currently, the morphology-based grading system of embryo quality assessment is also enhanced by morphokinetics analysis using time-lapse light microscopy [15,16,17] and deep machine learning methods [18]. However, embryo morphological phenotype evaluation alone is insufficient for accurate prediction of successful implantation in the situation of choice between multiple embryos of “excellent” (“A”) grade [19,20]. Moreover, aneuploid embryos can form morphologically high-scoring blastocysts or even reach the term of delivery. Despite the existing correlation between karyotype and morphology [21], classic grading or time-lapse microscopy is not an alternative to genetic analysis for accurate assessment of embryo karyotype [22,23,24,25]. At the same time, genetic abnormalities highly interfere with implantation and physiological development [26]. That is why methods of preimplantation genetic testing for aneuploidies (PGT-A) are generally used to check the karyotype of selected embryos before the transfer [27,28,29]. However, this invasive method can result in decrease of pregnancy rates if biopsies are performed at the cleavage stage [30,31]. No such effects were found after blastocysts biopsies, but the technical complexity of PGT-A using NGS or array hybridization makes it applicable only for frozen-thawed embryos to avoid day 6 fresh transfers, which were reported to be associated with lower success rates [32]. The high cost of PGT-A procedures can also be considered as a limiting factor for broad method application.
Limitations of morphology-based and genetic classification methods stimulated the search for new approaches for the non-invasive biochemical or molecular evaluation of embryo viability and quality. One of the promising non-invasive approaches includes measurements of biochemical intermediates, such as glucose, lactate, pyruvate, amino acids, b-HCG, or other proteins in the spent embryo culture medium [33,34,35,36]. Recently, several spectroscopy methods (Raman spectroscopy, NMR spectroscopy, IR spectroscopy) have been used to study the embryonic metabolome [37,38,39] as well as gas chromatography coupled with mass spectrometry [40]. The embryo metabolome might reflect the dynamic changes of cell activity regulated by the expression of certain sets of genes and their target proteins, affecting the molecular phenotype of the embryo and its implantation potential. The metabolic activity of the embryo can be assessed by examining the exometabolome of the spent culture medium, which is formed by dynamic nutrient intake and secretion of metabolites. Thus, exometabolomics of spent embryo culture medium (SECM) might provide useful insights on the functional features of embryonic development and overall embryo viability to improve the success rates of the eSET protocol.
One of the most advanced and powerful methods of metabolites profiling is the combination of high-performance liquid chromatography followed by mass spectrometry (HPLC-MS). This approach is widely approved in modern omics research, but it is still rarely used for human embryo culture media analysis. In this study, the HPLC-MS method was used to evaluate the differences in the metabolomic profiles of culture media from human embryos of different morphological classes and karyotypes. Metabolite content of paired SECM on the 3rd and 5th day of cultivation was also studied to show dynamic changes and check the reproducibility of results on different LC-MS systems. Finally, changes in the metabolomic profiles of SECM from embryos with different implantation outcomes after transfer were analyzed.

2. Results

2.1. Influence of Morphology and Karyotype on the Composition of Spent Embryo Culture Media

The composition of the metabolites in the culture media may depend both on the total metabolic activity of all the cells of the embryo, and on their number by a certain day of development. To evaluate the influence of the development stage on the embryo exometabolome, SECM samples (n = 153) obtained on day 5 of continuous incubation were selected without considering the karyotype. We assumed that the embryos that reached the optimal parameters of development, regardless of the karyotype, would have minor differences in metabolic activity. The number of samples in each class and their morphological grades are shown in Table 1. The distribution of embryos by classes was based on morphological criteria and the probability of their transfer. Only blastocysts with expansion score ≥3 according to Gardner classification were considered for the analysis. That is, first-class embryos always have priority for transfer, second-class embryos are transferred when first-class embryos are absent or when they show features of an aneuploid karyotype. Third-class embryos usually show development retardation signs and are generally not recommended for transfer. Finally, the fourth class includes embryos that stopped developing at the stage of one to several cells and can be considered as a negative control. Embryos in morulae stage having smaller number of cells and blastocysts with the 1 and 2 expansion grades were excluded due to low probability of clinical pregnancy and lower frequency of their transfer.
After processing the HPLC-MS data, 1046 ions were detected in the mass range m/z 100–804. Multivariate statistics methods are widely used to calculate and visualize variation of thousands of features found by untargeted metabolomics, resulting in fast sample clusterization and outlier filtering. To visualize the difference in metabolic profiles of the samples by morphological class, partial least-squares discriminant analysis (PLS-DA) was applied to the HPLC-MS data. PLS-DA showed a clear clustering of nutrient media samples from different morphological classes. Class 4, used as control, significantly differed from the other three classes, which may indicate a notable contribution to the changes in the composition of the nutrient media by actively developing embryos. It can also be seen that Class 1 tends to cluster away from the “good” and “fair” embryos from Classes 2 and 3, respectively (Figure 1).
To assess possible differences within the same class, a PLS-DA analysis was performed between embryos with different morphology in Class 1: 6AA, 6BA, 4AA, 4AB, 3AA (Figure 2). Grades 3AA having fewer cell numbers showed minor differences and tended to separate from the rest of the embryo samples, but no clear clustering within the class was observed. Thus, the different morphological classes of embryos showed class-specific features that distinguish the groups from each other, suggesting that the exometabolome profile of the blastocyst correlates with cell number and the stage of development. Although we observed the tendency of Classes 1–3 to separate, for further steps of the study we selected samples of only fully developed blastocysts of 3–6 AA grades from Class 1.
One of the most important criteria for selecting an embryo for transfer is the quality of its karyotype. In our study, a rather high incidence of karyotype abnormalities was detected by PGT-A: 12 of the 35 selected embryos from Class 1 and 21 of the 43 from Class 2 (Table 1). Among aneuploid embryos from both classes, twenty had one abnormality, eight had two, and five had three or more abnormalities (mono- or trisomies). To assess the influence of the altered karyotype on the metabolite profile, samples of nutrient media from aneuploid and euploid embryos of the first morphological class were analyzed. The score graph of the multivariate OPLS-DA revealed a good clustering of samples by karyotype (Figure 3). This fact indicates that the detected aneuploidies may affect the metabolic activity of embryos, despite the “normal” morphological phenotype.
To detect the molecular signatures that contribute the most to the differences between the euploid and aneuploid groups, molecular ions with a fold change of ≥2 were filtered by univariate statistics (p < 0.05) and mapped to Human Metabolome Database (HMDB, http://www.hmdb.ca/. Accessed on 20 January 2022). Among 23 molecular ions selected by such criteria, only two were identified by HMDB (Table 2). Their possible structures refer to: 16-hydroxy-10-oxohexadecanoic acid (m/z 309.2042) and 5,6-dihydroxyprostaglandin F1 (m/z 389.2517). Despite the clustering of the samples by the OPLS-DA method, the univariate statistical analysis with used parameters did not show any significant differences in common, easily identifiable components of the medium, such as amino acids. The content of most of the detected metabolites (22 out of 23) was increased in the culture media of aneuploid embryos, which indicates a possible metabolic shift in this group (Table 2). To further investigate the differences between groups, all molecular ions marked as significant for the OPLS-DA model were also mapped to HMDB, and their pre-identification within 15 ppm (Table S1) was used to enrich the KEGG pathways. The results of the pathway enrichment are presented in Table 3. It was shown that most of the enriched pathways relate to the amino acids, fatty acids, vitamins B and A metabolism. The most pronounced topological impact was shown for vitamin B6 metabolism and pentose phosphate pathway.

2.2. Implantation Potential

The excellent morphology and normal karyotype of the embryo do not guarantee its successful implantation, which further complicates the choice of the embryo for transfer. To study the possible correlation between the implantation potential of the transferred embryos and the composition of their nutrient media, a retrospective analysis of samples from euploid embryos with known transfer outcomes was performed. Based on the data obtained at the previous stage, morphologically similar embryos of the excellent class were selected.
In SECM from embryos continuously cultured for 5 days multivariate statistical analysis (OPLS-DA) revealed significant differences between the groups with the corresponding transfer outcomes, which indicates an essential role of the embryo’s metabolic rate in culture for subsequent successful implantation (Figure 4). Among the detected molecular ions, 241 were found to have significant (p < 0.05) fold change ≥2 in groups with different implantation potential (Table S2). To further search for possible implantation-specific signatures, molecular ions with a statistically significant fold change (≥2) between the groups were selected and mapped to the HMDB for identification. From 241 ions, only 100 had HMDB IDs giving 1214 different molecular identifications, of which only 142 had KEGG IDs indicating their possible involvement in characterized metabolic pathways (Table S2). Several selected ions of interest are shown in Table 4.
The known constituents of embryonic culture medium are amino acids, which may be used in synthetic and energetic metabolism. The pronounced decrease was observed in the relative abundance of molecular ions identified as phenylalanine following tryptophan, valine, and proline. HMDB identification was validated by comparing retention times of found amino acids with 13C15N-labeled amino-acid standards (Table S3, Figure S1). The most abundant molecular ions were fragmented for structural identification in MSn experiment and their initial identifications were confirmed (Figure 5). However, MSn analysis was limited by low absolute abundance of the most of active SECM components due to small initial sample volume. All molecular ions marked as significant for the OPLS-DA model were also mapped to HMDB, and their possible identification within 10 ppm of m/z was used to enrich the KEGG pathways. It was found that the most enriched pathways relate to the biosynthesis and metabolism of amino acids phenylalanine, valine, arginine, and several others (Table 5), pathways involved in lipid and nucleotide metabolism were also revealed (Table 5).
In discontinuous protocol HPLC-MS analysis of SECM samples of the third and fifth days revealed 1616 and 1404 molecular ions in the mass range m/z 86–804, respectively. These paired SECM samples were divided into two groups depending on the outcome of embryo implantation: successful (n = 5) or unsuccessful (n = 7). The data were subjected to multivariate statistical analysis using OPLS-DA to visualize clustering and find the ions that contribute most to the variance between the two groups. Explicit clustering of the two groups corresponding to different implantation outcomes was detected on both the third and fifth day of cultivation (Figure 6A,B). Molecular ions with a statistically significant fold change (≥2) between the groups were selected: 50 (17 decreased, 34 increased) and 36 (21 decreased, 15 increased) molecular ions on the third and fifth day, respectively (Table S4). In addition, the ions significant for the OPLS-DA model were mapped to HMDB, and the potential metabolites were enriched with the KEGG database (Table 6 and Table 7). Several amino acids, lipid and purine metabolism pathways were enriched both at day 3 and day 5 of cultivation. Interestingly, phenylalanine-associated pathways were enriched in day 5 SECM in both protocols but were absent in day 3 SECM of the discontinuous protocol. Retinol metabolism pathway was enriched in day 3 and 5 SECM of the discontinuous protocol, while in continuous protocol this pathway did not contribute to variation between groups of SECM with successful and unsuccessful implantation (Table 6 and Table 7).

2.3. Blank Embryo Culture Media Decay Analysis

To investigate possible changes of culture media composition over time and their effect on results, we cultivated blank culture media for five days and compared it with a fresh one. More than 650 molecular ions were initially detected by LC-MS analysis and groups (fresh and day 5 media) showed excellent clusterization by OPLS-DA, but only 60 molecular ions were found to change with statistical significance at p < 0.05 (with FDR correction). We compared these ions with molecular ions contributing to difference between groups with different implantation outcome to detect any batch effect. Only 9 matching molecular ions were found (Table S5). We used HMDB database to identify these molecular ions and found no interference with our results data. Thus, possible natural embryo media decay over 5 days of cultivation does not have any significant impact on results and is negligible.

3. Discussion

The elective single-embryo transfer has been shown in several studies to be the recommended method compared to multiple embryo transfer, since it demonstrates a similar birth rate and is devoid of disadvantages associated with multiple pregnancies [41,42,43]. However, this method requires a comprehensive assessment of the viability and implantation potential of a single embryo, as there is no room for error. It should be noted that currently there are no reliable criteria for choosing the best candidate from two euploid embryos with excellent morphological qualities. To address these issues, routine morphology analysis, even with time-lapse microscopy, should be supported by biochemical or molecular methods. One of such methods is metabolomic profiling which allows taking of a “snapshot” of current metabolites composition which might correlate with the embryo functional condition and implantation potential. While direct measurement of embryo metabolism is complicated, spent culture media exometabolome may provide crucial information on embryo development based on different nutrient consumption and metabolites excretion to media by embryos of different morphological quality or the same quality but different implantation potential.
In this study, metabolite profiles of spent culture media were obtained in two cultivation protocols: five days of continuous cultivation or five days of discontinuous cultivation with a complete refreshment of the media on day 3. The media samples were grouped into four classes according to embryo morphological grade and their quality for transfer (Table 1). This classification was based on the experience of ART specialists of our clinic and literature data. Gardner et al. showed that embryos with ≥3AA grades had the highest probability to give clinical pregnancy [44]. Additionally, in a recent article Zhao et al. 2018, based on Gardner’s classification, divided the embryos according to morphological characteristics into four classes: excellent (3–6, AA), good (3–6 AB, BA), medium (3–6 BB, AC, CA), poor (1–6, combinations of B and C). It was shown that embryos of excellent and good quality have a close probability of clinical pregnancy (65 and 59%, respectively), while in the middle and poor classes these parameters decreased (50 and 33%, respectively) [45]. In our study, we assigned some of the embryos of good quality (mentioned above) to an excellent class. We assumed that embryos with similar efficiency should have similar metabolic features; however, the number of cells in the embryo may also affect the quantitative parameters of metabolites in culture media. Unfortunately, the exact number of cells in each morphological group of human embryos is unknown (at least we could not find such data). This is apparently due to the lack of non-invasive evaluation methods and the inability to use many human embryos in experiments. Therefore, we conditionally divided the samples available at the time of the study, as indicated in Table 1.
Multivariate PLS-DA statistics were used for intergroup (Figure 1) and intragroup (Figure 2) comparative analysis. Intergroup PLS-DA for the samples obtained in continuous protocol showed clustering of ‘excellent’ (Class 1) and ‘poor’ (Class 4) groups of embryos from ‘good’ (Class 2) and ‘fair’ (Class 3) groups that were mostly overlapped. It can also be seen that Class 1 and Class 3 tend to separate while Class 2 overlaps with them. A possible explanation is that Class 2 includes groups of embryos that have a similar number of cells to both Class 1 and Class 3. Since the morphological evaluation by the embryologist does not allow for an accurate cell count of each embryo, this circumstance may also contribute to the observed intersections of clusters of similar morphological groups. This result might indicate that a correlation exists between the embryo’s developmental stage and its metabolic activity. Class 1 embryo is expected to have normal phenotype and thus have optimal metabolic activity, and metabolic rates of Class 2 and 3 embryos may be variable due to different possible origins of their “fair” phenotype and lower cell number, while Class 4 embryo should have the lowest metabolic rates due to significant retardation or the arrest of the development. It may be proposed that in embryos with the smaller cell numbers there is a shift in gene expression leading to the formation of ATP-dependent ion gradients necessary for blastocele formation that takes a large percentage of the embryo’s energy to create and maintain the blastocoel cavity [46]. In these embryos, this process may follow in expense of the others, such as cell division. Hence the embryo must possess some degree of metabolic competency to arrange its full-fledged development.
Intragroup comparison made for Class 1 embryos showed that metabolite profiles within the class are quite similar for different grades (Figure 2). One possible explanation may be that any significant metabolic profile shift would significantly influence the development with a consequent class change. Thus, in Class 1, minor variability is observed, with a slight clustering of type 3AA embryos, which may be associated with a smaller number of cells in this group. Additionally, in 3AA embryos, the metabolic pathways responsible for expansion are just beginning to activate compared to more developed ones. Detected differences between morphological groups indicate the need for careful selection of samples according to the specified criteria since they can have a significant impact on the results of the metabolomic analysis.
The metabolomic profiling of SECM has been done previously by a range of spectroscopic methods including Fourier Transform infrared spectroscopy (FT-IR), near-infrared spectroscopy (NIR), Raman spectroscopy [47,48] and NMR [39]. These studies showed differences in metabolism ratios which correlating with embryo viability and implantation potential [39,49,50,51,52]. Seli et al. used NIR spectroscopy and found that viability indices calculated on metabolic profiles of day 2–3 SECM are independent of embryo morphology and correlate with reproductive potential of embryos [53]. This work received considerable interest and discussion. Subsequent randomized control studies of spent embryo culture media on day 2 and day 5 by NIR spectroscopy in adjunct to morphology did not confirm any significant improvement of ongoing pregnancy rates compared to classic morphology-based criteria [54]. Another randomized control study of day 3 SECM with NIR spectroscopy also failed to find any improvement in transfer outcome by this combined approach [55]. Evaluation of the technical design of the experiment revealed that such results were obtained due to limitations of the existing equipment. Additionally, the signal-to-noise ratio threshold was not strong enough for correct viability predictions and noise filtration, which impaired method robustness and reproducibility [56,57].
Raman spectroscopy, while being a powerful non-destructive analysis method with high specificity, has relatively low sensitivity, which can limit its usage for analysis of complex biological fluids, such as serum plasma or spent culture media [58]. Nutrient consumption and metabolite excretion rates even between excellent, fully developed blastocyst and poor-quality morula can differ by a narrow margin, which requires highly sensitive analysis methods. Some highly active compounds, such as hormones or mediators, usually have extremely small concentrations, and their detection is also complicated. Mass-spectrometry coupled with gas or liquid chromatography is an effective and versatile method of analysis, that provides an alternative for spectroscopy-based metabolomic profiling. Cortezzi et al. used direct MS analysis, without chromatographic separation, of the day 3 culture medium with known implantation outcome to predict their implantation potential. The spectra were obtained in negative-ion mode and analyzed by PLS-DA to test the performance of the statistical model on implantation potential prediction. The model correctly identified 100% of positive implantations and 70% of negative ones. The culture media from embryos with positive or negative outcome showed specific biochemical signatures which contribute the most to the difference between groups, though no identification was provided for the detected ions [59]. While untargeted metabolomics of spent culture media is good for biomarkers search, there is often a need not to discover new, but thoroughly identify and quantify existing biomarkers found by biochemical or molecular methods. Targeted free fatty acid LC-MS profiling revealed a statistically significant decrease of docosahexaenoic acid and an insignificant decrease of other essential free fatty acids in day 6 spent culture from morphologically “good” embryos compared to “poor” ones. Authors suggested targeted metabolomic profiling as a powerful method for non-invasive assessment of embryo quality [60].
The genetic analysis revealed a substantial amount of Class 1 and 2 embryos with an abnormal karyotype (Table 1), hence we tested metabolomic profiling as a possible non-invasive alternative for PGT-A. Multivariate statistics (OPLS-DA) has divided culture media samples of the aneuploid group from the normal group within Class 1 (Figure 3), demonstrating the potential influence of the karyotype on metabolite content. Univariate statistics revealed 23 molecular ions differing between these groups (the fold change > 2). Interestingly, 22 of 23 had increased abundance, probably due to decreased uptake of medium components. Most of these 23 ions were not identified by HMDB search within 10 ppm mass accuracy. It can be assumed that with significant changes in the metabolism of aneuploid embryos, they would not have reached the appropriate stage of development and would not have shown the best phenotype. This may explain the fact that we were not able to detect pronounced changes among the known components of the nutrient media and their possible metabolites. However, multivariate statistics have revealed sets of ions, the cumulative change in the content of which has led to the clustering of groups depending on the karyotype. This suggests that aneuploid embryos may have specific differences in metabolism, and certain molecular signatures can be used to detect aneuploidies.
For the evaluation of the overall metabolic activity of embryos with different karyotypes, all significantly different molecular ions defined by multivariate statistics were mapped to the HMDB database, and metabolites identified were enriched by KEGG biochemical pathways database. Amino acids and fatty acids metabolism pathways were among the most enriched ones, which is fully in accordance with the known importance of amino acid turnover for normal embryo development [61]. Aneuploidies are known to have a devastating effect on embryonal development, affecting all levels of the organization, and only rare karyotype abnormalities like trisomy of chromosome 21 are tolerable enough to give live birth but associated with multiple organ pathologies and disorders [62,63].
Due to the incident nature of the karyotype abnormalities, we combined all types of aneuploidies including monosomy, trisomy, and their combinations in one group and analyzed them against euploid embryos. However, the multivariate statistical analysis still shows excellent clustering of euploid and aneuploid embryos, though the contribution of single metabolites might have been smoothed by sample heterogeneity (Figure 3). Further studies with samples grouping by individual types of aneuploidies will likely reveal specific features in the composition of metabolites that are characteristic of each karyotype abnormality. Research conducted by Sanchez-Ribas et al. confirms this assumption [36]. It was shown that several potential metabolites can distinguish trisomy 21 embryos from normal ones although authors had to group ‘pure’ trisomy 21 samples with those having trisomy 21 in combination with other aneuploidies due to extremely rare occurrence of ‘pure’ trisomy 21 [36]. The preliminary results from another study suggest that metabolomics markers may be sensitive to trisomy 18 and may be able to distinguish trisomy 18 from trisomy 21 [64]. Raman spectroscopy and multivariate statistics were recently used to detect aneuploid embryos. The method showed results consistent with PGT-A testing performed by next-generation sequencing, achieving 95.9% of correct detections of either euploidy or aneuploidy [65]. Significant differences in Raman spectra associated with small RNA and lipids were also observed. Thus, metabolomic profiling coupled with targeted metabolomics could be a feasible method for the non-invasive detection of aneuploidies.
In our study, we have also demonstrated that metabolomic profiling of spent culture media of day 5 embryos from the same morphological class can be used to assess their implantation potential. After five days of incubation molecular ions corresponding to phenylalanine, valine, proline, and tryptophan were found to be decreased in spent culture media from embryos with successful implantation, changes of other amino acid derivatives were observed as well (Table 4 and Table S3). Proline, valine, and phenylalanine identifications were confirmed using isotope-labeled standards (Table S3, Figure S1), phenylalanine and tryptophan were further confirmed by MSn analysis (Figure 5). This result is in accordance with previous reports of the amino acid turnover rate as a predictor of embryo viability [34]. It was also shown that differences in the composition of the culture media of the implanted embryos (Class 1) can be detected on the 3rd day of cultivation, as well as on the 5th day after media replacement. In spent culture media of days 3 and 5, 46 and 29 statistically differing molecular ions were found by univariate statistics. No amino acids were identified among the changing features except glutamine. That may be explained by the fact that culture media was refreshed on day 3 in this study, which might negate all accumulated composition differences observed after 5 days of cultivation in continuous protocol. It can be concluded that differences between embryos with distinct implantation potential can be seen starting at least from day 3, which may allow the application of culture medium analysis in different transfer protocols for both fresh and cryopreserved embryos.
Pathway enrichment of metabolites marked as important by multivariate statistics was performed for all implantation outcomes and karyotype comparisons. The most enriched pathways for compounds contributing to karyotype differences were related to the pentose phosphate pathway, vitamin B6, and amino acid metabolism pathways. Glucose-dependent nucleotide synthesis by the pentose phosphate pathway was found to control activation pathway for trophectoderm specific gene transcription during mammalian embryogenesis [66]. Amino acids turnover and high levels of vitamins A and B6 were previously reported to be associated with aneuploidies, embryo quality, and morphokinetics [61,67]. The enrichment of amino acid and vitamin B6 metabolism pathways, as well as lipid-associated pathways, were also observed (Table 5, Table 6 and Table 7) in culture media of embryos with successful or unsuccessful implantation. Drabkova et al. showed that there were not any significant changes in media amino acid composition in day 3 SECM. Additionally, positive correlations with glutamine and negative with glutamate concentration were reported [68]. Our findings also show increased glutamine depletion in day 3 medium from embryos successfully implanted later (Table S4). Amino acids are the key components of a variety of pathways, even those not directly involved into amino acid biosynthesis or degradation.
Earlier studies of lipid effects on embryo quality revealed a significant impact of fatty acids levels on embryo development [69,70]. Fatty acid derivatives such as eicosanoids may be produced on the blastocyst–endometrium interface and seem to be one of the key factors of the decidualization and implantation process [71]. Our results indicate that lipid metabolism may play significant role during embryo development and specific lipid compound sets should be investigated as potential biomarkers for karyotype and embryo assessment.
In a recent study, Harden S.L. et al. investigated the exometabolome of decidualizing human stromal cells in culture and concluded that the metabolic footprints of cellular microenvironment generated by different decidual subpopulations encode spatiotemporal information that may be important for their differentiation and optimal embryo implantation. It has been shown that in all investigated timepoints of decidualization the most enriched metabolic pathways of stromal cells exometabolome are phenylalanine, tyrosine and tryptophan biosynthesis, valine, leucine, and isoleucine biosynthesis, glycerophospholipid metabolism, linoleic acid metabolism, purine and pyrimidine metabolism [72]. It is interesting to note, that the similar pathways were enriched after analysis of day 5 human embryo exometabolome in our study. It is intriguing to propose that these pathways and their secreted metabolites are necessary for the embryonal and endometrial cells to form complementary microenvironment that guide the initial steps of the implantation. Together, these findings support the importance of found molecular signatures as possible markers of successful implantation. Thus, metabolomic profiling of spent culture media may be a promising method in combination with morphological analysis as it provides additional information on embryo quality. Even two morphologically identical embryos may have different metabolic activity and thus implantation potential, or class 2 embryo may have worth morphology but far better metabolic rate. This fact may be crucial for successful eSET, especially with increasing age, when the number of available embryos is limited.
It should be noted that not all successful implants develop into a full-fledged pregnancy and result in delivery. This may be due to different factors such as incorrect localization of the implanted embryo, abnormalities of the placenta development, and fetus gestation [73]. Therefore, after embryo transfer long-term monitoring and detailed analysis of pregnancy outcomes are necessary for further method development. Another limitation of the method may be the small volume of the probe, which makes precise identification and quantification of low abundant metabolites a challenging task. The method we used showed good results in aneuploidies detection, but we pulled all aneuploidy types to one group because of the rare occurrence of such cases. A longer collection period and more test samples for each type of aneuploidy are required for method refinement and validation. Additionally, we understand that our data was obtained using one type of the commercial culture media provided on the market for human embryo cultivation, thus any differences in composition and nutrient ratio (which are known to exist and often not evident due to unavailable media formulations) may lead to the specific metabolic shifts that in turn will change features of metabolomic profile. Thus, to obtain the set of “universal” biomarkers of embryo quality it seems reasonable to evaluate the difference between metabolomic profiles according to corresponding culture conditions.
In conclusion, our data suggest that the application of HPLC-MS for metabolomic analysis can be used as a powerful, sensitive, and versatile method that may provide valuable data both for embryo karyotype and implantation potential assessment. Further implementation of targeted metabolomics in combination with time-lapse morphology analysis may significantly improve the positive outcome ratio after embryo transfer.

4. Materials and Methods

4.1. Oocyte Collection, Embryo Cultivation, and Selection

Ovarian hyperstimulation was performed according to the protocol using gonadotropin and gonadotropin-releasing hormone antagonists. The stimulation protocol was the same for all participants. Follicular maturation was monitored by ultrasound examination and was estimated as sufficient if there were more than 5 follicles of ≥18 mm in diameter. Oocytes were collected by ultrasonography-guided transvaginal aspiration of the follicles. Oocyte fertilization was performed by intracytoplasmic sperm injection (ICSI). Human embryos were cultivated in 30 µL droplets of Continuous Single Culture medium (FUJIFILM Irvine Scientific, Santa Ana, CA, USA) for 3–5 days under mineral oil (FUJIFILM Irvine Scientific, Santa Ana, CA, USA) at 37 °C in a gas atmosphere of 89% N2, 5% O2, and 6% CO2. Assisted hatching was performed on the 5th day of cultivation, after the evaluation of embryo morphology by the expert embryologist. Embryos of different morphological quality were marked by alphanumerical indices according to the Gardner grade system [14]. The numbers indicate the expansion grade of a blastocyst (1–4 blastocoel cavity development, 5–6 hatching out of the shell), the first letter indicates intracellular mass quality, and the second letter is for trophectoderm layer quality. The letters have the following meaning: A—many densely packed cells, B—fewer cells, loose packing, C—just a few cells. Embryos were grouped into four arbitrary classes: first “Excellent”, second “Good”, third “Fair”, fourth “Poor”. The distribution of morphological grades between classes is presented in Table 1.
The karyotype evaluation was performed for the selected embryos of the 1st and 2nd classes, which are the most likely to be transferred during the single embryo transfer protocol. Blastocyst trophectoderm biopsy was performed for the preimplantation genetic testing for aneuploidy (PGT-A), and karyotype quality was assessed by array comparative genomic hybridization [74]). After the biopsy, day 5 embryos were cryopreserved and stored until transfer. After two menstrual cycles, according to the results of PGT, cryopreserved/thawed euploid embryos were transferred into the uterine cavity with the preimplantation support using micronized progesterone. Data from only the first transfer were considered for patients undergoing more than one embryo transfer, and only one cycle of oocyte donation was included in case of multiple retrievals. Clinical pregnancy after the transfer was defined by the positive hCG test two weeks after the transfer and ultrasound confirmation by visualization of the gestational sac on third week after embryo transfer and at 5–6 weeks to detect the heartbeat.

4.2. Metabolomic Profiling

4.2.1. Samples

Spent embryonic culture media samples (n = 153) obtained on day 5 of continuous incubation were selected for the study, of which 35, 43, 20, and 55 belonged to the 1st, 2nd, 3rd, and 4th morphological quality classes, respectively. In discontinuous protocol, SECM samples were collected on day 3, and cultivation was continued in fresh media till day 5. Paired SECM samples (n = 24) were obtained on day 3 and day 5 from 12 “excellent” quality embryos. Equal volumes (20 μL) of spent culture media were taken for the experiment, labeled, and stored at −80 °C. To assess possible changes in media composition during the cultivation period additional 6 control samples of blank embryo culture media were analyzed on day 3 and day 5 of incubation.

4.2.2. Sample Preparation

Methanol, acetonitrile, and formic acid (all are HPLC grade) were purchased from Merck (Darmstadt, Germany). Deionized water was prepared with a Milli-Q Water Purification System (Darmstadt, Germany). Samples for further HPLC-MS analysis were prepared by adding 60 µL of methanol to 20 µL of SECM, followed by mixing and centrifugation at 14,000× g. The resulting supernatants were aliquoted and kept at −80 °C until the analysis. Mix of 17 13C15N-labeled amino acids (Cambridge isotope, Cambridge, UK) was added to supernatants with final concentration of 25 mM as internal standard. The same mix in final concentration of 31.25 mM was used as QC standard after each 10 sample injections.

4.2.3. HPLC-MS Analysis

LC-MS experiments were carried out at the Institute of Biomedical Chemistry using equipment of “Human proteome” Core Facility Center (Moscow, Russia). For samples obtained after continuous cultivation, HPLC-MS analysis was performed using liquid chromatography system Ultimate 3000 Nano LC (Thermo Fisher Scientific, Waltham, MA, USA). The HPLC system was coupled to the Q Exactive Hybrid Quadrupole-Orbitrap mass spectrometer with ESI (Thermo Fisher Scientific, Waltham, MA, USA). The HPLC-MS method for the analysis was as follows: Zorbax SB-C18 column (150 × 0.5 mm, 5 µm, Agilent, Santa Clara, CA, USA), injection volume 20 µL in three replicates per sample, and column temperature 40 °C. The mobile phase A consisted of the aqueous solution of 0.1% formic acid, and the eluent B was water-acetonitrile-formic acid (20/80/0.1, v/v/v) mixture. The gradient program was the following: 0 min—5% B, 20 min—35% B, 25 min—99% B, 30 min—99% B, 30.5 min—5% B, 35 min—5% B. The flow rate was 0.2 mL/min. The mass spectrometer was operated in the positive-ion mode under the following conditions: mass range of m/z 100–1500, with a resolution of 60,000.
For samples obtained after discontinuous cultivation, HPLC-MS analyses were performed using liquid chromatography system Ultimate 3000 Nano LC (Thermo Fisher Scientific, Waltham, MA USA) coupled to the hybrid quadrupole time-of-flight Maxis Impact Q-ToF mass spectrometer with ESI (Bruker Daltoniks, Germany). The HPLC-MS method for the analysis was as follows: Zorbax SB-C18 column (150 × 0.5 mm, 5 µm, Agilent, Santa Clara, CA, USA), injection volume 20 µL in three replicates per sample, and column temperature 40 °C. The mobile phase A consisted of an aqueous solution of 0.1% formic acid, and the eluent B was 0.1% formic acid in acetonitrile. The gradient program was the following: 0 min—5% B, 15 min—5% B, 25 min—95% B, 30 min—95% B, 31 min—5% B, 34 min—5% B The flow rate was 0.04 mL/min. The mass spectrometer was operated in the positive-ion mode under the following conditions: mass range of m/z 100–1500, with a resolution of 50,000, capillary voltage 4.1 kV, nebulizer gas 0.4 Bar, dry gas 4.0 L/min, dry gas temperature 180 °C.

4.3. Data Processing and Statistical Analysis

Peak detection, noise signal removal, data processing, and analysis were carried out by XCMS package on Workflow4Metabolomics (www.workflow4metabolomics.org accessed on 8 January 2021) web platform [75] with the following parameters: peak detection algorithm was centwave; max tolerated m/z deviation was 10 ppm for Orbitrap data and 15 ppm for qToF data; minimum and maximum peak width was 5 s and 50 s respectively.
Partial least square discriminant analysis (PLS-DA) [76] and Orthogonal partial least square analysis (OPLS-DA) [77] were used for the multivariate statistics in multigroup and pairwise comparisons, respectively. t-test was used for pairwise univariate statistics. The significance threshold was p-value < 0.05. Data normalization and statistical analysis were performed by the Metaboanalyst web platform [78].

5. Patients

Retrospective study included patients (n = 21 and n = 12 for continuous and discontinuous protocols respectively), undergoing IVF programs from March 2017 to October 2019 in the Department of Assisted Reproductive Technologies of the National Medical Research Center of Obstetrics, Gynecology, and Perinatology named after Academician V.I.Kulakov (Ministry of Health of the Russian Federation). The study was approved by the Ethics Committee of the National Medical Research Center for Obstetrics, Gynecology, and Perinatology named after Academician V.I.Kulakov (Protocol 2 from 9 February 2017). All study participants provided written informed consent.
Inclusion criteria were age 24–37 years old, body mass index 18–29 kg/m2, primary infertile women with major indications for IVF, tuboperitoneal infertility, more than 5 oocytes were collected during IVF, the couples had declared no history of chromosomal abnormalities. Exclusion criteria were low ovarian reserve, uterus malformations, abnormal endometrial development, endometriosis, gynecological tumors, a history of recurrent pregnancy loss, any somatic or mental disorders contradictory for pregnancy and childbirth, any acute inflammatory processes, acute exacerbation of any chronic disease, pathozoospermia. Patients with no follicular growth effect after stimulation or unsuccessful oocyte retrieval were excluded as well.

Supplementary Materials

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

Author Contributions

Conceptualization, M.B.; Data curation, L.I. and S.Y.; Formal analysis, C.E., A.G., L.I. and S.Y.; Funding acquisition, G.S.; Investigation, C.E., L.I., S.Y., V.S. and N.M.; Methodology, A.G., V.C., V.S., E.K., N.M. and V.Z.; Project administration, M.B.; Resources, V.C., V.Z. and G.S.; Software, C.E.; Supervision, G.S.; Visualization, C.E.; Writing—original draft, C.E.; Writing—review & editing, E.K. and M.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was approved by the Ethics Committee of the National Medical Research Center for Obstetrics, Gynecology, and Perinatology named after Academician V.I.Kulakov (Protocol 2 from 9 February 2017).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Data is contained within the article or Supplementary Materials.

Acknowledgments

The authors wish to thank the laboratory of systems biology of the Orekhovich Institute of biomedical chemistry for providing Orbitrap LC-MS equipment for sample analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Steptoe, P.C.; Edwards, R.G. Birth after the reimplantation of a human embryo. Lancet 1978, 312, 366. [Google Scholar] [CrossRef]
  2. Wyns, C.; Bergh, C.; Calhaz-Jorge, C.; De Geyter, C.; Kupka, M.S.; Motrenko, T.; Rugescu, I.; Smeenk, J.; Tandler-Schneider, A.; Vidakovic, S.; et al. ART in Europe, 2016: Results generated from European registries by ESHRE. Hum. Reprod. Open 2020, 3, hoaa032. [Google Scholar] [CrossRef]
  3. Long, X.; Wang, Y.; Wu, F.; Li, R.; Chen, L.; Qian, W.; Qiao, J. Pregnancy Outcomes of Single/Double Blastocysts and Cleavage Embryo Transfers: A Retrospective Cohort Study of 24, 422 Frozen-Thawed Cycles. Reprod. Sci. 2020, 27, 2271–2278. [Google Scholar] [CrossRef]
  4. Park, D.S.; Kim, J.W.; Eum, J.H.; Lee, W.S.; Yoon, T.K.; Lyu, S.W. Clinical and pregnancy outcomes of double and single blastocyst transfers related with morphological grades in vitrified-warmed embryo transfer. Taiwan J. Obstet. Gynecol. 2020, 59, 398–402. [Google Scholar] [CrossRef] [PubMed]
  5. Kissin, D.M.; Kulkarni, A.D.; Kushnir, V.A.; Jamieson, D.J. Number of embryos transferred after in vitro fertilization and good perinatal outcome. Obstet. Gynecol. 2014, 123, 239–247. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. McLernon, D.J.; Harrild, K.; Berg, C.; Davies, M.J.; De Neubourg, D.; Dumoulin, J.C.M.; Gerris, J.; Kremer, J.A.; Martikainen, H.; Mol, B.W.; et al. Clinical effectiveness of elective single versus double embryo transfer: Meta-analysis of individual patient data from randomised trials. BMJ 2011, 342, 34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Saccone, G.; Zullo, F.; Roman, A.; Ward, A.; Maruotti, G.; Martinelli, P.; Berghella, V. Risk of spontaneous preterm birth in IVF-conceived twin pregnancies. J. Matern. Neonatal Med. 2019, 32, 369–376. [Google Scholar] [CrossRef] [Green Version]
  8. Young, B.C.; Wylie, B.J. Effects of Twin Gestation on Maternal Morbidity. Semin. Perinatol. 2012, 36, 162–168. [Google Scholar] [CrossRef]
  9. Xue, Y.; Li, K.; Zhang, S. The effect on twinning rate of transferring double vitrified-warmed embryos in women of advanced reproductive age: A retrospective study. PeerJ 2020, 1, e8308. [Google Scholar] [CrossRef]
  10. Mackay, A.P.; Berg, C.J.; Atrash, H.K. Pregnancy-related mortality from preeclampsia and eclampsia. Obstet. Gynecol. 2001, 97, 533–538. [Google Scholar]
  11. Cutting, R. Single embryo transfer for all. Best Pract. Res. Clin. Obstet. Gynaecol. 2018, 53, 30–37. [Google Scholar] [CrossRef] [PubMed]
  12. Hornstein, M.D. State of the ART: Assisted Reproductive Technologies in the United States. Reprod. Sci. 2016, 23, 1630–1633. [Google Scholar] [CrossRef] [PubMed]
  13. Papanikolaou, E.G.; Kolibianakis, E.M.; Tournaye, H.; Venetis, C.A.; Fatemi, H.; Tarlatzis, B.; Devroey, P. Live birth rates after transfer of equal number of blastocysts or cleavage-stage embryos in IVF. A systematic review and meta-analysis. Hum. Reprod. 2008, 23, 91–99. [Google Scholar] [CrossRef] [PubMed]
  14. Gardner, D.K.; Lane, M.; Stevens, J.; Schlenker, T.; Schoolcraft, W.B. Blastocyst score affects implantation and pregnancy outcome: Towards a single blastocyst transfer. Fertil. Steril. 2000, 73, 1155–1158. [Google Scholar] [CrossRef]
  15. Pribenszky, C.; Mátyás, S.; Kovács, P.; Losonczi, E.; Zádori, J.; Vajta, G. Pregnancy achieved by transfer of a single blastocyst selected by time-lapse monitoring. Reprod. Biomed. Online 2010, 21, 533–536. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wong, C.C.; Loewke, K.E.; Bossert, N.L.; Behr, B.; De Jonge, C.J.; Baer, T.M.; Reijo Pera, R.A. Non-invasive imaging of human embryos before embryonic genome activation predicts development to the blastocyst stage. Nat. Biotechnol. 2010, 28, 1115–1121. [Google Scholar] [CrossRef]
  17. Meseguer, M.; Herrero, J.; Tejera, A.; Hilligsøe, K.M.; Ramsing, N.B.; Remoh, J. The use of morphokinetics as a predictor of embryo implantation. Hum. Reprod. 2011, 26, 2658–2671. [Google Scholar] [CrossRef] [Green Version]
  18. Kragh, M.F.; Rimestad, J.; Berntsen, J.; Karstoft, H. Automatic grading of human blastocysts from time-lapse imaging. Comput. Biol. Med. 2019, 115, 103494. [Google Scholar] [CrossRef]
  19. Wang, Q.; Sun, Q.-Y.Y. Evaluation of oocyte quality: Morphological, cellular and molecular predictors. Reprod. Fertil. Dev. 2007, 19, 1–12. [Google Scholar] [CrossRef]
  20. Morbeck, D.E. Blastocyst culture in the Era of PGS and FreezeAlls: Is a ‘C’ a failing grade? Hum. Reprod. Open 2017, 2017, hox017. [Google Scholar] [CrossRef] [Green Version]
  21. Hardarson, T.; Selleskog, U.; Reismer, E.; Wigander, A.; Wennerström, S.; Westin, C.; Ahlström, A.; Hanson, C. Zygotes cleaving directly into more than two cells after 25–27 h in culture are predominantly chromosomally abnormal. In Human Reproduction; Oxford University Press: Oxford, UK, 2006; p. 21. [Google Scholar]
  22. Figueira, R.C.D.S.; Setti, A.S.; Braga, D.P.A.F.; Iaconelli, A.; Borges, E. Blastocyst morphology holds clues concerning the chromosomal status of the embryo. Int. J. Fertil. Steril. 2015, 9, 215–220. [Google Scholar] [CrossRef] [Green Version]
  23. Gardner, D.K.; Balaban, B. Assessment of human embryo development using morphological criteria in an era of time-lapse, algorithms and “OMICS”: Is looking good still important? Mol. Hum. Reprod. 2016, 22, 704–718. [Google Scholar] [CrossRef]
  24. Zaninovic, N.; Irani, M.; Meseguer, M. Assessment of embryo morphology and developmental dynamics by time-lapse microscopy: Is there a relation to implantation and ploidy? Fertil. Steril. 2017, 106, 722–729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Forman, E.J.; Upham, K.M.; Cheng, M.; Zhao, T.; Hong, K.H.; Treff, N.R.; Scott, R.T., Jr. Comprehensive chromosome screening alters traditional morphology-based embryo selection: A prospective study of 100 consecutive cycles of planned fresh euploid blastocyst transfer. Fertil. Steril. 2013, 100, 718–724. [Google Scholar] [CrossRef] [PubMed]
  26. Alfarawati, S.; Fragouli, E.; Colls, P.; Stevens, J.; Gutiérrez-Mateo, C.; Schoolcraft, W.B.; Katz-Jaffe, M.G.; Wells, D.; Path, F.R.C. The relationship between blastocyst morphology, chromosomal abnormality, and embryo gender. Fertil. Steril. 2011, 95, 520–524. [Google Scholar] [CrossRef]
  27. Yang, Z.; Liu, J.; Collins, G.S.; Salem, S.A.; Liu, X.; Lyle, S.S.; Peck, A.C.; Sills, E.S.; Salem, R.D. Selection of single blastocysts for fresh transfer via standard morphology assessment alone and with array CGH for good prognosis IVF patients: Results from a randomized pilot study. Mol. Cytogenet. 2012, 5, 24. [Google Scholar] [CrossRef] [Green Version]
  28. Forman, E.J.; Hong, K.H.; Ferry, K.M.; Tao, X.; Taylor, D.; Levy, B.; Treff, N.R.; Scott, R.T., Jr. In vitro fertilization with single euploid blastocyst transfer: A randomized controlled trial. Fertil. Steril. 2013, 100, 100–107.e1. [Google Scholar] [CrossRef] [Green Version]
  29. Scott, R.T.; Upham, K.M.; Forman, E.J.; Hong, K.H.; Scott, K.L.; Taylor, D.; Tao, X.; Treff, N.R. Blastocyst biopsy with comprehensive chromosome screening and fresh embryo transfer significantly increases in vitro fertilization implantation and delivery rates: A randomized controlled trial. Fertil. Steril. 2013, 100, 697–703. [Google Scholar] [CrossRef]
  30. Mastenbroek, S.; Twisk, M.; van der Veen, F.; Repping, S. Preimplantation genetic screening: A systematic review and meta-analysis of RCTs. Hum. Reprod. Update 2011, 17, 454–466. [Google Scholar] [CrossRef] [Green Version]
  31. Scott, R.T.; Upham, K.M.; Forman, E.J.; Zhao, T.; Treff, N.R. Cleavage-stage biopsy significantly impairs human embryonic implantation potential while blastocyst biopsy does not: A randomized and paired clinical trial. Fertil. Steril. 2013, 100, 624–630. [Google Scholar] [CrossRef]
  32. Coates, A.; Kung, A.; Mounts, E.; Hesla, J.; Bankowski, B.; Barbieri, E.; Ata, B.; Cohen, J.; Munné, S. Optimal euploid embryo transfer strategy, fresh versus frozen, after preimplantation genetic screening with next generation sequencing: A randomized controlled trial. Fertil. Steril. 2017, 107, 723–730.e3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Gardner, D.K.; Lane, M.; Stevens, J.; Schoolcraft, W.B. Noninvasive assessment of human embryo nutrient consumption as a measure of developmental potential. Fertil. Steril. 2001, 76, 1175–1180. [Google Scholar] [CrossRef]
  34. Houghton, F.D.; Hawkhead, J.A.; Humpherson, P.G.; Hogg, J.E.; Balen, A.H.; Rutherford, A.J.; Leese, H.J. Non-invasive amino acid turnover predicts human embryo developmental capacity. Hum. Reprod. 2002, 17, 999–1005, Erratum in Hum. Reprod. 2003, 18, 1756–1757. [Google Scholar] [CrossRef] [Green Version]
  35. Domínguez, F.; Gadea, B.; Esteban, F.J.; Horcajadas, J.A.; Pellicer, A.; Simón, C. Comparative protein-profile analysis of implanted versus non-implanted human blastocysts. Hum. Reprod. 2008, 23, 1993–2000. [Google Scholar] [CrossRef] [Green Version]
  36. Sánchez-Ribas, I.; Riqueros, M.; Vime, P.; Puchades-Carrasco, L.; Jönsson, T.; Pineda-Lucena, A.; Ballesteros, A.; Dominguez, F.; Simon, C. Differential metabolic profiling of non-pure trisomy 21 human preimplantation embryos. Fertil. Steril. 2012, 98, 1157–1164.e2. [Google Scholar] [CrossRef] [PubMed]
  37. Ahlström, A.; Wikland, M.; Rogberg, L.; Barnett, J.S.; Tucker, M.; Hardarson, T. Cross-validation and predictive value of near-infrared spectroscopy algorithms for day-5 blastocyst transfer. Reprod. Biomed. Online 2011, 22, 477–484. [Google Scholar] [CrossRef] [Green Version]
  38. Huang, Z.; Sun, Y.; Wang, J.; Du, S.; Li, Y.; Lin, J.; Feng, S.; Lei, J.; Lin, H.; Chen, R.; et al. Rapid and nondestructive method for evaluation of embryo culture media using drop coating deposition Raman spectroscopy. J. Biomed. Opt. 2013, 18, 127003. [Google Scholar] [CrossRef] [Green Version]
  39. Wallace, M.; Cottell, E.; Cullinane, J.; McAuliffe, F.M.; Wingfield, M.; Brennan, L. 1H NMR based metabolic profiling of day 2 spent embryo media correlates with implantation potential. Syst. Biol. Reprod. Med. 2014, 60, 58–63. [Google Scholar] [CrossRef] [Green Version]
  40. Krisher, R.L.; Heuberger, A.L.; Paczkowski, M.; Stevens, J.; Pospisil, C.; Prather, R.S.; Sturmey, R.G.; Herrick, J.R.; Schollcraft, W.B. Applying metabolomic analyses to the practice of embryology: Physiology, development and assisted reproductive technology. Reprod. Fertil. Dev. 2015, 27, 602–620. [Google Scholar] [CrossRef]
  41. Lee, A.M.; Connell, M.T.; Csokmay, J.M.; Styer, A.K. Elective single embryo transfer- the power of one. Contracept. Reprod. Med. 2016, 1, 11. [Google Scholar] [CrossRef] [Green Version]
  42. Ubaldi, F.M.; Capalbo, A.; Colamaria, S.; Ferrero, S.; Maggiulli, R.; Vajta, G.; Sapienza, F.; Cimadomo, D.; Giuliani, M.; Gravotta, E.; et al. Reduction of multiple pregnancies in the advanced maternal age population after implementation of an elective single embryo transfer policy coupled with enhanced embryo selection: Pre- and post-intervention study. Hum. Reprod. 2015, 30, 2097–2106. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Klenov, V.E.; Boulet, S.L.; Mejia, R.B.; Kissin, D.M.; Munch, E.; Mancuso, A.; Van Voorhis, B.J. Live birth and multiple birth rates in US in vitro fertilization treatment using donor oocytes: A comparison of single-embryo transfer and double-embryo transfer. J. Assist. Reprod. Genet. 2018, 35, 1657–1664. [Google Scholar] [CrossRef] [PubMed]
  44. Gardner, D.K.; Stevens, J.; Sheehan, C.B.; Schoolcraft, W.B. Analysis of blastocyst morphology. In Human Preimplantation Embryo Selection; Elder, K., Cohen, J., Eds.; Informa Healthcare: London, UK, 2007; pp. 79–87. [Google Scholar]
  45. Zhao, Y.Y.; Yu, Y.; Zhang, X.W. Overall Blastocyst Quality, Trophectoderm Grade, and Inner Cell Mass Grade Predict Pregnancy Outcome in Euploid Blastocyst Transfer Cycles. Chin. Med. J. 2018, 131, 1261–1267. [Google Scholar] [CrossRef] [PubMed]
  46. Biggers, J.D.; Bell, J.E.; Benos, D.J. Mammalian blastocyst: Transport functions in a developing epithelium. Am. J. Physiol. 1988, 255 4 Pt 1, 419–432. [Google Scholar] [CrossRef]
  47. Botros, L.; Sakkas, D.; Seli, E. Metabolomics and its application for non-invasive embryo assessment in IVF. Mol. Hum. Reprod. 2008, 14, 679–690. [Google Scholar] [CrossRef] [Green Version]
  48. Lundin, K.; Ahlström, A. Quality control and standardization of embryo morphology scoring and viability markers. Reprod. Biomed. Online 2015, 31, 459–471. [Google Scholar] [CrossRef] [Green Version]
  49. Seli, E.; Sakkas, D.; Scott, R.; Kwok, S.C.; Rosendahl, S.M.; Burns, D.H. Noninvasive metabolomic profiling of embryo culture media using Raman and near-infrared spectroscopy correlates with reproductive potential of embryos in women undergoing in vitro fertilization. Fertil. Steril. 2007, 88, 1350–1357. [Google Scholar] [CrossRef]
  50. Scott, R.; Seli, E.; Miller, K.; Sakkas, D.; Scott, K.; Burns, D.H. Noninvasive metabolomic profiling of human embryo culture media using Raman spectroscopy predicts embryonic reproductive potential: A prospective blinded pilot study. Fertil. Steril. 2008, 90, 77–83. [Google Scholar] [CrossRef]
  51. Vergouw, C.G.; Botros, L.L.; Roos, P.; Lens, J.W.; Schats, R.; Hompes, P.G.A.; Burns, B.H.; Lambalk, C.B. Metabolomic profiling by near-infrared spectroscopy as a tool to assess embryo viability: A novel, non-invasive method for embryo selection. Hum. Reprod. 2008, 23, 1499–1504. [Google Scholar] [CrossRef] [Green Version]
  52. Parks, J.C.; McCallie, B.R.; Janesch, A.M.; Schoolcraft, W.B.; Katz-Jaffe, M.G. Blastocyst gene expression correlates with implantation potential. Fertil. Steril. 2011, 95, 1367–1372. [Google Scholar] [CrossRef]
  53. Seli, E.; Vergouw, C.G.; Morita, H.; Botros, L.; Roos, P.; Lambalk, C.B.; Yamashita, N.; Kato, O.; Sakkas, D. Noninvasive metabolomic profiling as an adjunct to morphology for noninvasive embryo assessment in women undergoing single embryo transfer. Fertil. Steril. 2010, 94, 535–542. [Google Scholar] [CrossRef] [PubMed]
  54. Hardarson, T.; Ahlstrm, A.; Rogberg, L.; Botros, L.; Hillensj, T.; Westlander, G.; Sakkas, D.; Wikland, M. Non-invasive metabolomic profiling of Day 2 and 5 embryo culture medium: A prospective randomized trial. Hum. Reprod. 2012, 27, 89–96. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Vergouw, C.G.; Hanna Kostelijk, E.; Doejaaren, E.; Hompes, P.G.A.; Lambalk, C.B.; Schats, R. The influence of the type of embryo culture medium on neonatal birthweight after single embryo transfer in IVF. Hum. Reprod. 2012, 27, 2619–2626. [Google Scholar] [CrossRef] [Green Version]
  56. Sakkas, D. Embryo selection using metabolomics. Methods Mol. Biol. 2014, 1154, 533–540. [Google Scholar]
  57. Vergouw, C.G.; Heymans, M.W.; Hardarson, T.; Sfontouris, I.A.; Economou, K.A.; Ahlström, A.; Rogberg, L.; Lainas, T.G.; Sakkas, D.; Kieslinger, D.C.; et al. No evidence that embryo selection by near-infrared spectroscopy in addition to morphology is able to improve live birth rates: Results from an individual patient data meta-analysis. Hum. Reprod. 2014, 29, 455–461. [Google Scholar] [CrossRef] [Green Version]
  58. Ember, K.J.I.; Hoeve, M.A.; McAughtrie, S.L.; Bergholt, M.S.; Dwyer, B.J.; Stevens, M.M.; Faulds, K.; Forbes, S.J.; Campbell, C.J. Raman spectroscopy and regenerative medicine: A review. NPJ Regen. Med. 2017, 2, 15. [Google Scholar] [CrossRef]
  59. Cortezzi, S.S.; Cabral, E.C.; Trevisan, M.G.; Ferreira, C.R.; Setti, A.S.; De Almeida Ferreira Braga, D.P.; Figueira, R.C.; Iaconelli, A., Jr.; Eberlin, M.N.; Borges, E., Jr. Prediction of embryo implantation potential by mass spectrometry fingerprinting of the culture medium. Reproduction 2013, 145, 453–462. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Yagi, A.; Miyanaga, S.; Shrestha, R.; Takeda, S.; Kobayashi, S.; Chiba, H.; Kamiya, H.; Hui, S. A fatty acid profiling method using liquid chromatography-high resolution mass spectrometry for improvement of assisted reproductive technology. Clin. Chim. Acta 2016, 456, 100–106. [Google Scholar] [CrossRef] [PubMed]
  61. Picton, H.M.; Elder, K.; Houghton, F.D.; Hawkhead, J.A.; Rutherford, A.J.; Hogg, J.E.; Leese, H.J.; Harris, S.E. Association between amino acid turnover and chromosome aneuploidy during human preimplantation embryo development in vitro. Mol. Hum. Reprod. 2010, 16, 557. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zhu, J.; Tsai, H.J.; Gordon, M.R.; Li, R. Cellular Stress Associated with Aneuploidy. Dev. Cell 2018, 44, 420–431. [Google Scholar] [CrossRef] [Green Version]
  63. Tšuiko, O.; Jatsenko, T.; Parameswaran Grace, L.K.; Kurg, A.; Vermeesch, J.R.; Lanner, F.; Altmäe, S.; Salumets, A. A speculative outlook on embryonic aneuploidy: Can molecular pathways be involved? Dev. Biol. 2019, 447, 3–13. [Google Scholar] [CrossRef] [PubMed]
  64. Bahado-Singh, R.O.; Akolekar, R.; Chelliah, A.; Mandal, R.; Dong, E.; Kruger, M.; Wishart, D.S.; Nikolaides, K. Metabolomic analysis for first-trimester trisomy 18 detection. Am. J. Obstet. Gynecol. 2013, 209, 65.e1–65.e9. [Google Scholar] [CrossRef]
  65. Liang, B.; Gao, Y.; Xu, J.; Song, Y.; Xuan, L.; Shi, T.; Wang, N.; Hou, Z.; Zhao, Y.; Huang, W.; et al. Raman profiling of embryo culture medium to identify aneuploid and euploid embryos. Fertil. Steril. 2019, 111, 753–762.e1. [Google Scholar] [CrossRef] [PubMed]
  66. Chi, F.; Sharpley, M.S.; Nagaraj, R.; Roy, S.S.; Banerjee, U. Glycolysis-Independent Glucose Metabolism Distinguishes TE from ICM Fate during Mammalian Embryogenesis. Dev. Cell 2020, 53, 9–26.e4. [Google Scholar] [CrossRef] [PubMed]
  67. Gode, F.; Akarsu, S.; Dikmen, Z.G.; Tamer, B.; Isik, A.Z. The Effect Follicular Fluid Vitamin A, E, D and B6 on Embryo Morphokinetics and Pregnancy Rates in Patients Receiving Assisted Reproduction. Gynecol. Obstet. Reprod. Med. 2019, 25, 89–95. [Google Scholar] [CrossRef] [Green Version]
  68. Drábková, P.; Andrlová, L.; Hampl, R.; Kanďár, R. Amino acid metabolism in human embryos. Physiol. Res. 2016, 65, 823–832. [Google Scholar] [CrossRef]
  69. Haggarty, P.; Wood, M.; Ferguson, E.; Hoad, G.; Srikantharajah, A.; Milne, E.; Hamilton, M.; Bhattacharya, S. Fatty acid metabolism in human preimplantation embryos. Hum. Reprod. 2006, 21, 766–773. [Google Scholar] [CrossRef] [Green Version]
  70. Mirabi, P.; Chaichi, M.J.; Esmaeilzadeh, S.; Ali Jorsaraei, S.G.; Bijani, A.; Ehsani, M.; Hashemi Karooee, S.F. The role of fatty acids on ICSI outcomes: A prospective cohort study. Lipids Health Dis. 2017, 16, 18. [Google Scholar] [CrossRef] [Green Version]
  71. Kennedy, T. Interactions of eicosanoids and other factors in blastocyst implantation. In Eicosanoids and Reproduction; MTP Press Limited: Vancouver, BC, Canada, 2012; p. 73. [Google Scholar]
  72. Harden, S.L.; Zhou, J.; Gharanei, S.; Diniz-da-Costa, M.; Lucas, E.S.; Cui, L.; Murakami, K.; Fang, J.; Chen, Q.; Brosens, J.J.; et al. Exometabolomic Analysis of Decidualizing Human Endometrial Stromal and Perivascular Cells. Front. Cell Dev. Biol. 2021, 9, 626619. [Google Scholar] [CrossRef]
  73. Nikitina, T.V.; Sazhenova, E.A.; Tolmacheva, E.N.; Sukhanova, N.N.; Kashevarova, A.A.; Skryabin, N.A.; Vasilyev, S.A.; Nemtseva, T.N.; Yuriev, S.Y.; Lebedev, I.N. Comparative Cytogenetic Analysis of Spontaneous Abortions in Recurrent and Sporadic Pregnancy Losses. Biomed. Hub 2016, 1, 446099. [Google Scholar] [CrossRef]
  74. Aleksandrova, N.V.; Shubina, E.S.; Ekimov, A.N.; Kodyleva, T.A.; Mukosey, I.S.; Makarova, N.P.; Kulakova, E.V.; Levkov, L.A.; Barkov, I.Y.; Trofimov, D.Y.; et al. Comparative results of preimplantation genetic screening by array comparative genomic hybridization and new-generation sequencing. Mol. Biol. 2017, 51, 269–273. [Google Scholar] [CrossRef]
  75. Giacomoni, F.; Le Corguillé, G.; Monsoor, M.; Landi, M.; Pericard, P.; Pétéra, M.; Duperier, C.; Tremblay-Franco, M.; Martin, J.F.; Jacob, D.; et al. Workflow4Metabolomics: A collaborative research infrastructure for computational metabolomics. Bioinformatics 2015, 31, 1493–1495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Wold, S.; Ruhe, A.; Wold, H.; Dunn, W.J., III; The Collinearity Problem in Linear Regression. The Partial Least Squares (PLS) Approach to Generalized Inverses. SIAM J. Sci. Stat. Comput. 1984, 5, 735–743. [Google Scholar] [CrossRef] [Green Version]
  77. Bylesjö, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J.K.; Holmes, E.; Trygg, J. OPLS discriminant analysis: Combining the strengths of PLS-DA and SIMCA classification. J. Chemom. 2006, 20, 341–351. [Google Scholar] [CrossRef]
  78. Chong, J.; Wishart, D.S.; Xia, J. Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinform. 2019, 68, e86. [Google Scholar] [CrossRef]
Ijms 23 02706 g001 550
Figure 1. Score plot of intergroup PLS-DA of day 5 culture medium samples from embryos of different morphological classes.
Figure 1. Score plot of intergroup PLS-DA of day 5 culture medium samples from embryos of different morphological classes.
Ijms 23 02706 g001
Ijms 23 02706 g002 550
Figure 2. Score plot of intragroup comparison of metabolite profiles of day 5 culture medium samples from the class 1 embryos by PLS-DA. Samples are grouped according to morphological classification for day 5 embryos.
Figure 2. Score plot of intragroup comparison of metabolite profiles of day 5 culture medium samples from the class 1 embryos by PLS-DA. Samples are grouped according to morphological classification for day 5 embryos.
Ijms 23 02706 g002
Ijms 23 02706 g003 550
Figure 3. Score plot of OPLS-DA comparative analysis of metabolomic profiles of day 5 culture media samples from the 1st class embryos with euploid and aneuploid karyotypes.
Figure 3. Score plot of OPLS-DA comparative analysis of metabolomic profiles of day 5 culture media samples from the 1st class embryos with euploid and aneuploid karyotypes.
Ijms 23 02706 g003
Ijms 23 02706 g004 550
Figure 4. Score plot of OPLS-DA comparative multivariate statistical analysis of metabolomic profiles of day 5 culture media samples from embryos of the first morphological class with successful and unsuccessful implantation.
Figure 4. Score plot of OPLS-DA comparative multivariate statistical analysis of metabolomic profiles of day 5 culture media samples from embryos of the first morphological class with successful and unsuccessful implantation.
Ijms 23 02706 g004
Ijms 23 02706 g005 550
Figure 5. Mass spectra of MSn analysis for L-Phenylalanine and L-tryptophan.
Figure 5. Mass spectra of MSn analysis for L-Phenylalanine and L-tryptophan.
Ijms 23 02706 g005
Ijms 23 02706 g006 550
Figure 6. Score plot of OPLS-DA of metabolomic profiles of day 3 (A) and day 5 (B) culture media samples from euploid embryos of 1st morphological class with successful and unsuccessful implantation.
Figure 6. Score plot of OPLS-DA of metabolomic profiles of day 3 (A) and day 5 (B) culture media samples from euploid embryos of 1st morphological class with successful and unsuccessful implantation.
Ijms 23 02706 g006
Table
Table 1. Distribution of 5-day embryos into arbitrary classes. Alphanumerical morphological grades for blastocysts in Classes 1–3 were evaluated based on Gardner classification. The numbers of samples in each group and the distribution of euploid/aneuploid embryos are shown as well.
Table 1. Distribution of 5-day embryos into arbitrary classes. Alphanumerical morphological grades for blastocysts in Classes 1–3 were evaluated based on Gardner classification. The numbers of samples in each group and the distribution of euploid/aneuploid embryos are shown as well.
GroupsClass 1Class 2Class 3Class 4
“Excellent”“Good”“Fair”“Poor”
Grades6AA, 6BA, 4AA, 4AB, 3AA4BA, 4BB, 3AB, 3BB, 3BA3BC, 3CB, 3CC10b-5b, 10c-3c, 2pb, 2pn, atresia
KaryotypeEuploidAneuploidEuploidAneuploidN/DN/D
Number231222212055
Table
Table 2. Molecular ions which distinguish day 5 culture media of aneuploid and euploid embryos. RT—retention time; FC—Fold Change (significantly different between groups with p < 0.05 based on univariate statistics). Italics indicate molecular ions identified by HMDB: m/z 309.204—16-hydroxy-10-oxohexadecanoic acid; m/z 389.252—5,6-dihydroxyprostaglandin F1.
Table 2. Molecular ions which distinguish day 5 culture media of aneuploid and euploid embryos. RT—retention time; FC—Fold Change (significantly different between groups with p < 0.05 based on univariate statistics). Italics indicate molecular ions identified by HMDB: m/z 309.204—16-hydroxy-10-oxohexadecanoic acid; m/z 389.252—5,6-dihydroxyprostaglandin F1.
m/zRetention TimeFold Changem/zRetention TimeFold Change
314.7800.942.42354.8961.023.28
450.7550.962.29330.9341.052.31
448.7890.972.03309.2041.342.76
436.8710.972.64370.28111.983.11
333.8440.972.97248.10613.62.31
411.8600.974.72354.28614.582.03
234.1710.983.25412.32816.642.06
530.7690.992.64389.25216.653.11
316.870.992.04384.29616.66−2.00
410.8570.992.52368.27316.672.79
402.8251.002.66169.97829.82.02
Table
Table 3. Pathway enrichment of molecular ion identifications contributing to the difference in day 5 SECM of aneuploid and euploid embryos. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Table 3. Pathway enrichment of molecular ion identifications contributing to the difference in day 5 SECM of aneuploid and euploid embryos. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Pathway NameMatchedpImpact
Pentose phosphate pathway5/220.510.48
Valine, leucine, and isoleucine biosynthesis2/80.530
Vitamin B6 metabolism2/90.60.57
Phenylalanine metabolism2/100.660.14
alpha-Linolenic acid metabolism2/130.80
Butanoate metabolism2/150.860.03
Histidine metabolism2/160.880.28
Retinol metabolism 2/170.90.22
Caffeine metabolism1/100.910
Biotin metabolism1/100.910.2
Terpenoid backbone biosynthesis2/180.920.01
Pentose and glucuronate interconversions2/180.920.08
Arginine and proline metabolism5/380.930.26
Pantothenate and CoA biosynthesis 2/190.930.01
Glycine, serine, and threonine metabolism4/330.940.15
Valine, leucine, and isoleucine degradation5/400.950.03
Alanine, aspartate, and glutamate metabolism3/280.950.09
Table
Table 4. Selected molecular ions identified in SECM from implanted embryos after 5 days of continuous cultivation. Identifications were performed by mapping molecular ions m/z to HMDB metabolites database (accuracy ≤ 15 ppm). Fold change indicates relative change of mean integrated area of chromatographic peak for this ion when compared with SECM from embryos with unsuccessful implantation.
Table 4. Selected molecular ions identified in SECM from implanted embryos after 5 days of continuous cultivation. Identifications were performed by mapping molecular ions m/z to HMDB metabolites database (accuracy ≤ 15 ppm). Fold change indicates relative change of mean integrated area of chromatographic peak for this ion when compared with SECM from embryos with unsuccessful implantation.
m/zRetention TimeNameFold Change
110.0581.044-Aminobutyraldehyde−3.03
114.0661.12Creatinine−2.42
116.0711.06L-Proline−2.27
118.0861.06L-Valine−5.8
119.0821.35L-2,4-diaminobutyric acid−2.09
130.0871.00Pipecolic acid−2.9
132.0771.55Creatine−2.91
147.0668.1Adipic acid−7.47
153.0551.03p-Hydroxyphenylacetic acid−2.14
163.0721.00L-4-Hydroxyglutamine−5.04
166.0872.58L-Phenylalanine−88.24
219.1117.83N-Acetylserotonin−15.24
227.087.62L-Tryptophan−2.7
264.1069.68L-erythro-tetrahydrobiopterin−3.32
269.0913.67Inosine−2.01
280.0931.05Glycerophosphocholine−2.28
309.2041.33Hexadecanedioic acid−2.61
469.0511.04CDP-ethanolamine−2.36
Table
Table 5. Pathway enrichment of molecular ions identifications contributing to the difference in SECM from embryos with successful and unsuccessful implantation in continuous protocol. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Table 5. Pathway enrichment of molecular ions identifications contributing to the difference in SECM from embryos with successful and unsuccessful implantation in continuous protocol. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Pathway NameMatched pImpact
Vitamins
Vitamin B6 metabolism2/90.480.13
Lipids
Linoleic acid metabolism3/50.040
Steroid biosynthesis5/420.880.09
Amino acids
Valine, leucine and isoleucine biosynthesis2/80.420
Valine, leucine and isoleucine degradation4/400.940.04
Phenylalanine metabolism2/100.540.36
Histidine metabolism2/160.80.22
Arginine and proline metabolism5/380.820.1
Cysteine and methionine metabolism4/330.860.14
Tyrosine metabolism5/420.880.12
Aminoacyl-tRNA biosynthesis5/480.940
Other
Butanoate metabolism2/150.770.14
Propanoate metabolism2/230.930.04
Caffeine metabolism3/100.250.69
Pentose and glucuronate interconversions4/180.390.2
Nicotinate and nicotinamide metabolism1/150.940
Table
Table 6. Pathway enrichment of molecular ions identifications contributing to the difference in 3rd day SECM from embryos with successful and unsuccessful implantation in discontinuous protocol. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Table 6. Pathway enrichment of molecular ions identifications contributing to the difference in 3rd day SECM from embryos with successful and unsuccessful implantation in discontinuous protocol. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Pathway NameMatchedpImpact
Vitamins
Vitamin B6 metabolism2/90.430.54
Retinol metabolism5/170.110.22
Lipids
alpha-Linolenic acid metabolism3/130.340.33
Glycerophospholipid metabolism3/360.940.12
Biosynthesis of unsaturated fatty acids4/360.840.00
Terpenoid backbone biosynthesis1/180.960.00
Sphingolipid metabolism1/210.970.00
Amino Acids
Histidine metabolism1/160.940.19
Arginine biosynthesis1/140.910.08
Arginine and proline metabolism3/380.950.14
Cysteine and methionine metabolism2/330.980.06
Valine, leucine, and isoleucine degradation3/400.960.01
Valine, leucine, and isoleucine biosynthesis3/80.120.00
Sugars
Pentose and glucuronate interconversions1/180.960.14
Starch and sucrose metabolism1/180.960.00
Other
Caffeine metabolism2/100.480.00
Table
Table 7. Pathway enrichment of molecular ions identifications contributing to the difference in 5th day SECM from embryos with successful and unsuccessful implantation in discontinuous protocol. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Table 7. Pathway enrichment of molecular ions identifications contributing to the difference in 5th day SECM from embryos with successful and unsuccessful implantation in discontinuous protocol. Molecular ions filtered by PLS-DA were mapped to HMDB database, list of identified metabolites was then mapped to KEGG pathways database. ‘p’ refers to the raw p-value determined by Fisher’s exact test; ‘Impact’ refers to topological impact on pathway (routing betweenness centrality score).
Pathway NameMatchedpImpact
Vitamins
Vitamin B6 metabolism2/90.630.54
Retinol metabolism3/170.770.37
Amino acids
Phenylalanine metabolism4/100.170.74
Phenylalanine, tyrosine, and tryptophan biosynthesis2/40.221.00
Valine, leucine, and isoleucine, biosynthesis3/80.260.00
Arginine biosynthesis2/140.860.08
Arginine and proline metabolism6/380.890.11
Lipids
Arachidonic acid metabolism10/360.280.26
Terpenoid backbone biosynthesis3/180.810.19
Sugars
Galactose metabolism5/270.760.12
Other
Ubiquinone and other terpenoid-quinone biosynthesis2/90.630.00
Caffeine metabolism 5/100.051.00
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eldarov, C.; Gamisonia, A.; Chagovets, V.; Ibragimova, L.; Yarigina, S.; Smolnikova, V.; Kalinina, E.; Makarova, N.; Zgoda, V.; Sukhikh, G.; et al. LC-MS Analysis Revealed the Significantly Different Metabolic Profiles in Spent Culture Media of Human Embryos with Distinct Morphology, Karyotype and Implantation Outcomes. Int. J. Mol. Sci. 2022, 23, 2706. https://doi.org/10.3390/ijms23052706

AMA Style

Eldarov C, Gamisonia A, Chagovets V, Ibragimova L, Yarigina S, Smolnikova V, Kalinina E, Makarova N, Zgoda V, Sukhikh G, et al. LC-MS Analysis Revealed the Significantly Different Metabolic Profiles in Spent Culture Media of Human Embryos with Distinct Morphology, Karyotype and Implantation Outcomes. International Journal of Molecular Sciences. 2022; 23(5):2706. https://doi.org/10.3390/ijms23052706

Chicago/Turabian Style

Eldarov, Chupalav, Alina Gamisonia, Vitaliy Chagovets, Luiza Ibragimova, Svetlana Yarigina, Veronika Smolnikova, Elena Kalinina, Nataliya Makarova, Victor Zgoda, Gennady Sukhikh, and et al. 2022. "LC-MS Analysis Revealed the Significantly Different Metabolic Profiles in Spent Culture Media of Human Embryos with Distinct Morphology, Karyotype and Implantation Outcomes" International Journal of Molecular Sciences 23, no. 5: 2706. https://doi.org/10.3390/ijms23052706

APA Style

Eldarov, C., Gamisonia, A., Chagovets, V., Ibragimova, L., Yarigina, S., Smolnikova, V., Kalinina, E., Makarova, N., Zgoda, V., Sukhikh, G., & Bobrov, M. (2022). LC-MS Analysis Revealed the Significantly Different Metabolic Profiles in Spent Culture Media of Human Embryos with Distinct Morphology, Karyotype and Implantation Outcomes. International Journal of Molecular Sciences, 23(5), 2706. https://doi.org/10.3390/ijms23052706

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