Embryo classification beyond pregnancy: early prediction of first trimester miscarriage using machine learning

Purpose First trimester miscarriage is a major concern in IVF-ET treatments, accounting for one out of nine clinical pregnancies and for up to one out of three recognized pregnancies. To develop a machine learning classifier for predicting the risk of cleavage-stage embryos to undergo first trimester miscarriage based on time-lapse images of preimplantation development. Methods Retrospective study of a 4-year multi-center cohort of 391 women undergoing intra-cytoplasmatic sperm injection (ICSI) and fresh single or double embryo transfers. The study included embryos with positive indication of clinical implantation based on gestational sac visualization either with first trimester miscarriage or live-birth outcome. Miscarriage was determined based on negative fetal heartbeat indication during the first trimester. Data were recorded and obtained in hospital setting and research was performed in university setting. Results A minimal subset of six non-redundant morphodynamic features were screened that maintained high prediction capacity. Features that account for the distribution of the nucleolus precursor bodies within the small pronucleus and pronuclei dynamics were highly predictive of miscarriage outcome as evaluated using the SHapley Additive exPlanations (SHAP) methodology. Using this feature subset, XGBoost and random forest models were trained following a 100-fold Monte-Carlo cross validation scheme. Miscarriage was predicted with AUC 0.68 to 0.69. Conclusion We report the development of a decision-support tool for identifying the embryos with high risk of miscarriage. Prioritizing embryos for transfer based on their predicted risk of miscarriage in combination with their predicted implantation potential is expected to improve live-birth rates and shorten time-to-pregnancy. Supplementary Information The online version contains supplementary material available at 10.1007/s10815-022-02619-5.


Introduction
Most miscarriages occur during the first trimester of pregnancy (first 13 weeks), accounting for one out of nine clinical pregnancies and for up to one out of three recognized pregnancies (positive pregnancy test) [1,2]. Unlike second trimester miscarriages and third trimester fetal loses that occur due to multiple reasons (congenital defects, placental problems, cervical insufficiency, and infections) [3,4], aneuploidy is the main cause of first trimester miscarriages [5]. Early identification of chromosomal abnormalities that permit normal preimplantation development is performed today via preimplantation genetic testing for aneuploidy (PGT-A). PGT-A is an invasive procedure that involves piercing of the zona pellucida and removal of several trophectoderm cells. Since PGT-A typically requires over 24 h in clinical settings, blastocysts are often cryopreserved, thus eliminating fresh transfers. In addition, failure of amplification and mosaicism in PGT-A contribute to false negative prediction [6,7]. Since euploid embryos have a higher potential to implant, the implementation of non-invasive methods for early identification of embryos with high risk of miscarriage has the potential to improve clinical performance [8].
Machine learning-based decision support tools are increasingly incorporated into the healthcare system [9][10][11] and specifically in the assessment of the developmental potential of preimplantation embryos [12]. Supervised machine learning is the most widely used approach in which retrospective datasets with known outcome (aka labels) are employed for training a classification model. Prediction of non-labeled embryos is then performed by identifying patterns that are associated with positive or negative outcome using the trained model. The incorporation of time-lapse incubation systems in IVF clinics provided high quality video files that record the course of preimplantation embryo development. These data facilitated the prediction of the potential of embryos to reach blastulation [13,14] and to implant within the uterus [15][16][17]. Retrospectively labeled datasets of significant size can support the training of convolutional neural networks (CNN) that test the relationships between multiple features by convoluting image pixels. CNN's were used for performing automated morphokinetric annotation [18,19] and for predicting implantation [20][21][22]. Indeed, we found that CNN classifiers of embryo blastulation and implantation identify dynamic features beyond the discrete morphological elements and morphokinetic events that are broadly used by non-convolutional algorithms [20].
In IVF treatments, embryos that are selected for transfer according to their predicted implantation potential may possess high risk of miscarriage despite exhibiting normal morphological and morphokinetic profiles. The unmet need that arises is to generate technologies that will allow the assessment of embryo miscarriage potential in parallel to existing technologies that assess the implantation potential. Hence, our overall goal was to develop a decision-support tool that will perform accurate and robust evaluation of miscarriage potential non-invasively and within a reasonable time frame to alleviate the requirement for cryopreservation. In line with existing decision support tools, here we employed machine learning models to evaluate miscarriage potential based on time-lapse recordings of embryos from IVF clinics during preimplantation stages. Due to the inherent limitations that are associated with prospective studies, including ethical challenges and study size, we performed a multicenter retrospective study. Specifically, a dataset of 464 transferred and successfully implanted embryos that were obtained from four medical centers and labeled according to miscarriage (MC) or live-birth (LB) outcomes was assembled. The fact all the embryos in the dataset that were available for training the classification models had been carefully selected for transfer in the clinics based on their time-lapse visualization and were further successfully implanted introduces a significant computational challenge for predicting miscarriage outcome. Akin to the prediction of embryo implantation outcome, accuracy of MC prediction is expected to be limited by the fact that no information is provided on the uterine factors nor on the related variation between recipients [23]. Despite these challenges, here we present a machine learning classification model of MC prediction within 3 days from ICSI fertilization and demonstrate AUC ~ 0.7. In this manner, our decision support tool can support both cleavage-stage and blastocyst transfer strategies. The implementation of our decision support tool in the clinic is expected to improve live-birth rates and shorten time-to-pregnancy by deselecting for transfer the embryos with high risk of MC.

Patients and embryos
Data from 391 women who underwent fresh single or double embryo transfers were collected from IVF clinics of four medical centers as recently reported (Fig. 1a) [20]. All embryos in the dataset were fertilized via intracytoplasmic sperm injection (ICSI). PGT-A-tested embryos were discarded. A total of 464 positively implanted embryos were included in this study as determined on week 5 by ultrasound imaging of gestational sacs. MC labeling of 1 3 96 embryos was determined based on negative ultrasound indication of fetal heartbeat on week 13 or earlier. The remaining 368 embryos were labeled LB in accordance with a confirmed live-birth pregnancy outcome. Importantly, MC or LB labeling was based on known outcome of each individual embryo. For example, in the case of a twin implantation that proceeds to become a vanishing twin pregnancy due to anembryonic implantation or early clinical loss, the embryos remained unlabeled since the identity of the miscarried embryo and the singleton pregnancy embryo cannot be separated.

Labeled video database
Anonymized time-lapse video files were recorded by timelapse incubators (Embryoscope, VitroLife) and imported together with the associated clinical metadata. Data curation was performed using a PostgreSQL database that we assembled and formulated as reported recently [20]. In short, the database consisted of a front-end website that supported display, query, and data annotation (Hebrew University IT) and maintained by CHELEM LTD. Morphokinetic annotations and quality assurance were performed by qualified and experienced embryologists [20].

Extraction of morphological features
Static morphological features are routinely used for evaluating the potential of embryos to implant [24][25][26][27]. The embryos were characterized by their two-cell and four-cell fragmentation percentage and blastomere size symmetry, as annotated in adherence to established criteria [28]. Additional morphological and geometrical features of the zona pellucida, ooplasm, pronuclei (PN'i), and nucleolus precursor bodies (NPB's) were semi-automatically extracted per the third frame prior to tPNf (~ 1 h) using a custom designed code (Python). The orientation angle of the first cleavage plane was evaluated using the seven z-stack images [29]. The annotation of all morphological features is summarized in Table 1.

Extraction of morphokinetic features
Time-lapse incubators facilitate the annotation of the time points of discrete events. Morphokinetic annotations and quality assurance were performed by qualified and experienced embryologists as we recently reported [20]. In short, morphokinetic profiles were obtained via majority voting according to established protocols [30,31] and further validated blindly by an expert embryologist. Morphokinetic annotations included the time of PN appearance (tPNa) and fading (tPNf), and cleavage of N blastomeres (tN, N = 2 to 8) [32,33]. We also introduced all 36 pairwise intervals and pairwise ratios between the events.

DLR
We calculated the distances from the linear regression (DLR) as follows. Pairs of the morphokinetic features of LB embryos appear to be linearly correlated as demonstrated for t2-t3 events (Fig. 2). These morphokinetic relationships are defined by the linear regression of the morphokinetic scatter plot (R 2 = 0.87). To estimate the association with the averaged morphokinetic profiles of high quality positivelyimplanted embryos, we calculated the shortest Euclidean distance of each embryo from these linear Regression (Fig. 2).

PN dynamics
Male and female PN'i were visible for 17 ± 2.5 h (average ± STD) between tPNa and tPNf. The trajectories of the PN'I were extracted using optical flow as follows. Fifteen non-adjacent pixels were selected within each PN using the starting points frame. The time-lapse consecutive coordinates of each episode relative to the embryo frame of

Prediction models
Due to the limited size of the dataset, we employed XGBoost (XG) and random forest (RF) models that minimize the risk of overfitting by integrating ensembles of multiple classification processes [35,36]. In addition, XG and RF are based on clusters of decision trees that are suitable for identifying complex patterns which are derived by nonlinear relationships between the features. The XG model was implemented using XGBoost Python library and the RF model was implemented using the Scikit-learn Python package [37]. We set the number of trees and depth 6 and 1 for the XG model, and 65 and 3 for the RF model, and used the default feature importance scores for both models. We set 0.25 learning rate for the XG model and the RF splitting criterion was defined based on Gini impurity. Parameter search optimization over the above models was performed using the Hyperopt Python framework package [38]. Hyperparameter optimization was

Validation and testing methodologies
The embryos were divided into a train/validation set and an uncontaminated test set with balanced numbers of MC and LB labels (Fig. 1b). To effectively increase the variation between train set embryos, we performed a 100-fold Monte-Carlo cross validation scheme. Specifically, the classifier was trained on randomly selected 90% of the relevant train-validation set embryos and validated on the remaining 10%. In this scheme, each embryo was selected for validation multiple times, thus allowing the classifier to "learn" a given embryo as it is contrasted by different combinations of trainset embryos. The fact that the validation set AUC is in high agreement with the test set AUC for both models increases confidence that the Monte-Carlo cross-validation scheme did not generate overfitting (Fig. 8). MC prediction was performed by averaging the cross-validated model's weights.

Statistical analyses and graphical design
Statistical analyses were performed using NumPy, Scikitlearn, pandas, and SciPy Python packages. Figures were generated using Matplotlib and Seaborn Python packages. Manual features were extracted using the Pillow Python package.

Python packages version stamping
See Supplementary Table S1.

Results
Motivated by the prediction of pregnancy using machine learning, we hypothesized that early pregnancy loss can also be predicted based on time-lapse imaging of preimplantation development. Hence, we assembled a dataset of videos of MC-and LB-labeled embryos that were collected from four IVF clinics that accept patients with a diverse ethnic and racial backgrounds (Eastern and Western European Jews, North African and Middle Eastern Jews, Arabs, and others) [39]. This embryo dataset spans a range of maternal age. Clinics H1 and H3 tend to perform double embryo transfers on day-3 of preimplantation development or earlier whereas single embryo transfers are frequently performed in clinics H3 and H4 in which half of the embryos are transferred on day-5 or later (Fig. 1a). We first compared the morphokinetic profiles of MC and LB embryos, and found indistinguishable profiles, suggesting that the risk of MC cannot be assessed using the standard statistical approaches (Fig. 3a). Similarly, there were no differences between MC and LB embryos in terms of their developmental states and day of transfer (Fig. 3b). We conclude that the ethnic diversity, maternal age distributions, and inherent variation in the clinical protocols between the hospitals contribute to the generality of our study. In addition, the overlapping morphokinetic profiles, day of transfer, and the developmental states of the transferred embryos are consistent with communal criteria for selecting both MC and LB embryos for transfer, thus minimizing the effects of confounding parameters and satisfying a natural experiment.

Feature extraction
To identify visual markers of first trimester miscarriage, we assembled an expansive set of 314 features that were derived from measurable static and dynamic properties of preimplantation embryo development and span a range of length scales from the NPB to the zona pellucida (Table-1 and Methods). The feature set included morphological, morphokinetic, and other dynamic features as defined in the "Methods" section and summarized in Table 1. Morphological features account for specific geometrical properties of the embryo as derived from single frames. Morphokinetic features included the time points of morphokinetic events, their pairwise intervals and pairwise ratios, and the DLR features that quantify an effective distance of an embryo from the average profile of positively implanted embryos. The remaining features are associated with PN dynamics as evaluated using optical flow (the "Methods" section). Optical flow tracking was validated against manual tracking, showing only insignificant disagreement that was upperbounded by one quarter of PN radius (Fig. 4a). Interestingly, distinctive PN trajectories were observed and classified into one of the following categories based on their end-to-end distance and contour curvature (Fig. 4b): (i) curved trajectories; (ii) linear trajectories; (iii) stationary trajectories. By obtaining the temporal coordinates of the PN'i trajectories with respect to the embryo frame of reference, we defined specific features that were related to PN step size and velocity, the distribution of distances between PN'i, and their distance from embryo center of mass (Table 1).

Feature screening
The breadth of the feature set allowed us to examine a variety of morphodynamic elements and processes. However, the size of the labeled embryo dataset limited the number of training features in order to avoid overfitting. Hence, we carried out feature screening following a stepwise scheme based on statistical distance criteria between MC and LB embryo distributions followed by MC prediction performance using random forest (RF) and XGBoost (XG) models. Statistical distance between MC and LB embryos were calculated by the Kolmogorov-Smirnov (KS) test and the Kullback-Leibler (KL) divergence (Fig. 5a). Four features were distinguished by high KL-entropy and low KS p-value (green symbols): (1) t2-t3 distance from linear regression (DLR 2,3 ; the "Methods" section); (2) maximal distance between the NPBs of the small PN; (3) the difference between the distances of the two PN'i from the embryo center of mass (at tPNf); (4) the maximal distance between pairs of NPBs in the small PN normalized by small PN radius. These features were thus expected to provide high predictive strength.
KS-KL statistical distance criteria estimate the predictive potential of each feature individually. Since our goal is to identify the minimal group of features with high predictive strength, we applied a series of KS p-value and KL entropy thresholds, selected the subset of qualified features with lower KS p-value and higher KL entropy, and assessed their predictive potential by training RF and XG models. Six KS p-value thresholds, KS i , and three KL entropy thresholds, KL j , were chosen (Fig. 5a, black frame), thus defining eighteen feature subsets: To better assess classification performance by simulating multiple validation sets, we executed a 100-fold Monte-Carlo cross-validation experiment (Methods). The predictive  (Fig. 5b). Next, we selected the ten feature subsets with the highest RF and XG AUC values (Fig. 5b, red frame). In total, these ten subsets included 64 features; a fraction of which were shared by multiple subsets. The features were then ranked by their RF and XG importance scores, as averaged across all relevant Monte-Carlo crossvalidation experiments (Fig. 5c). In this manner, we identified the eight top ranked features (Fig. 5c, red frame).

Selection of a non-redundant feature subset
So far, features were screened based on their individual MC-LB statistical distances and their weighted contributions to RF and XG prediction models. However, redundancy among the eight selected features cannot be excluded.
To obtain further insight, we calculated the pairwise Pearson Correlation Matrix between the distributions of their feature values (Fig. 5d). Indeed, feature pairs [cf1, cf3] and [cf2, cf8] were highly correlated. Specifically, the difference between the smallest steps of the two PN'i (cf1) and the difference between the slowest velocities of the two PN'i (cf3) were highly correlated. Likewise, the correlation between the non-normalized (cf2) and the normalized (cf8) maximal distance between NPBs was also high (Table-1).
To identify the minimal subset of non-redundant features that maintain high predictive strength, we first performed feature selection via backward elimination. Starting with the selected group of eight, each feature was removed in its turn and RF and XG 100-fold Monte-Carlo cross-validation was performed on the remaining subsets of seven features (the "Methods" section). The single feature that had the smallest negative effect on RF and XG AUC's was removed. Backward feature elimination was performed in a stepwise manner, thus identifying a series of subsets with a decreasing number of features from eight-to-one (Fig. 5e, blue symbols). We found that the removal of the first three features did not decrease both RF and XG prediction strengths. However, subsequent removal of the remaining features decreased both XG and RF prediction performance, indicating that the subset of five features represents the minimal collection of non-redundant features with high MC prediction. Consistent with the pairwise Pearson correlation analysis discussed above, cf8, cf2, and cf3 features were indeed discarded in this order (Fig. 5d).
Backward feature selection did not exclude the possibility that reintroducing additional features from the original pool of features that had been screened out would not improve MC prediction. The goal was to assess the integrated contributions of multiple features rather than testing features one by one. Since testing all combinatorial subsets is computationally not feasible, we employed a two-step forward feature selection scheme.
Starting with the subset of five-features, the addition of all feature pair combinations was tested by training RF and XG models as described above. The feature pair that improved prediction the most was identified and the feature with the lower index was arbitrarily selected. In this manner, we performed a five + two, six + two, and seven + two forward feature selection cycles thus increasing subset size from five to eight features (Fig. 5e, red symbols). We found that introducing the sixth feature indeed improved the predictive strength of both RF and XG models. However, subsequent addition of the seventh and eighth features failed to increase AUC values. We thus conclude that minimal redundancy and high predictive strength of MC outcome is provided by the final subset of five backward selected features plus one forward selected feature as summarized in Table-1 (right  column). The minimal subset of non-redundant and highly predictive features is described below. f 1 : the maximal distance between NPB pairs within the small PN normalized to PN radius (Fig. 6a): f 2 : radius of the ooplasm.
[Both f 1 and f 2 were evaluated using the central focal plane image that was recorded 1 h prior to tPN f ]. f 3 : DLR 2,3 (Fig. 2). f 4 : DLR 5,6 . f 5 : the ratio between the maximal step size of the large PN, r l i , and the small PN, r s i , (Fig. 6b): Maximal steps are evaluated across N steps performed by the PN'i from tPN a and tPN f . f 6 : the absolute difference between the distances of the large ( l l ) and small ( l s ) PN from zygote center of mass at tPN f (Fig. 6b): We note that half of the features are associated either with NPB pairwise dynamics ( f 1 ) or with PN dynamics ( f 5 and f 6 ). The remaining features correspond to ooplasm size ( f 2 ) and to cleavage-stage morphokinetic profiles ( f 3 and f 4 ).

Prediction of MC outcome
MC prediction was performed by training an XG model using the selected feature subset f 1 to f 6 (Fig. 7a). The earliest time for MC prediction is set by the time of six blastomere cleavage event as required for feature f 6 . To effectively increase the variation between train set embryos, we performed a 100-fold Monte-Carlo cross validation scheme (the "Methods" section). As expected, the area under the receiver operating characteristic (ROC) curve (AUC) of train set embryos was relatively high, indicative of moderate overfitting, which is consistent with the size of the available dataset. Nevertheless, AUC of 68-69% was robustly obtained both on the validation set and the test set embryos. Importantly, comparable prediction performance was also obtained both by an RF model and by an integrated XG-RF classifier, thus increasing confidence in the utility of this feature subset for MC prediction (Fig. 8). To verify that including H1 clinic in both classification model training, for which there were only LB embryos available, did not introduce confounding effects, we removed all H1 embryos, including the three H1 test set embryos, plotted the ROC curve and calculated the AUC values ( Supplementary Fig. S1). Satisfyingly, we found agreement with the AUC values as evaluated for the embryos from all clinics (Fig. 8a-b).
To demonstrate clinical utility and delineate the tradeoff between precision (positive predictive value; PPV) and sensitivity (recall) in relation with improving livebirth rates and shortening time-to-pregnancy, we illustrate two extreme scenarios (Fig. 7b). We first consider a women patient having only a small number of high quality embryos with a predicted high potential of LB outcome. Due to their small number, embryos with both moderate and high predicted LB score should not be discarded by setting a high classification threshold for MC prediction (equivalent to setting a low threshold for LB prediction; Fig. 7c-i). The tradeoff of setting a high MC threshold is the inclusion of embryos with moderate predicted risk of MC outcome. In the opposite case, we consider a woman patient with a large number of embryos with of a high predicted LB potential. Here, the clinicians have the benefit of discarding all embryos with a predicted moderate-to-high risk of MC outcome while also discarding embryos with moderate predicted LB potential by setting a relatively low threshold for MC prediction (equivalent to setting a high threshold for LB prediction; Fig. 7c-ii). Precision is favored in first scenario, which supports increasing live-birth rates at the expense of not shortening time-to-pregnancy due to MC. In the second scenario, the clinicians can sustain losing high-quality embryos in order to shorten time-to-pregnancy by lowering the risk of MC.

Feature importance
PN and NPB morphologies were previously associated with aneuploidy [40][41][42][43][44][45] Here we find that the morpho-dynamic properties as characterized by the trajectories of the PN'i and the NPB's differentiate between the potential of MC versus LB outcomes. The maximal distance between NPB pairs within the small PN ( f 1 ) is greater in MC embryos than in LB embryos ( Fig. 9a; p-value 0.001). In addition, MC outcome is associated with a larger ratio between the maximal step sizes of the large versus the small PN'i ( f 5 ) and a shorter difference between PN'i locations with respect to zygote center of mass ( f 6 ). Morphokinetically, DLR 2,3 distances ( f 3 ) are longer in MC embryos whereas DLR 5,6 distances ( f 4 ) are shorter. This indicates that the potential of LB is more closely associated with the averaged profile of positively-implanted embryos with respect to t2 − t3 correlations (p-value < 0.01) and that the potential of MC is more closely associated with the averaged profile of positively implanted embryos with respect to t5 − t6 correlation (p-value < 0.04). Lastly, we find that ooplasm radius ( f 2 ) of MC embryos is slightly larger than of LB embryos ( Fig. 9a; p-value < 0.01).
To evaluate the impact that features have on MC prediction by the XG classifier, we employed the SHapley Additive exPlanations (SHAP) methodology [46]. Features that obtained high positive SHAP values contribute to MC prediction of the embryos whereas high negative SHAP values contribute to LB prediction. We find that for all six features, the embryos with high feature values are separated from embryos with low features values along the SHAP axis ( Fig. 9b). Specifically, high f 1 , f 2 , f 3 , and f 5 values, which are obtained by embryos with positive SHAP, contribute to MC prediction. On the other hand, high f 4 and f 6 contribute to LB prediction (Fig. 9b). Unlike SHAP that measures the relationship between feature values and prediction, feature importance is assessed using the adjusted SHAP (aSHAP), which scores accurate prediction by adjusting a negative sign to false prediction [20]. Indeed, we find that all the features of the final subset, with the exception of f 5 , hold comparable importance in directing accurate MC prediction (Fig. 9c).

Discussion
Time-lapse incubation systems provide high-quality visualization of preimplantation embryos that facilitate the development of machine learning-based classifiers. Such classifiers are routinely utilized within IVF clinics worldwide for selecting the embryos for transfer with the highest implantation potential. However, these algorithms are sensitive to various morphological elements and morphokinetic features that mark the potential to implant but fail to detect features that are associated with the risk of MC outcome. Using the same video files of preimplantation development, our algorithm predicts MC potential as early as the time of 6-cell cleavage event based on the following six morpho-dynamic features: (1) the normalized maximal distance between NPB within the smaller PN; (2) radius of the ooplasm; (3,4) the DLR's of morphokinetic event pairs t2-t3 and t5-t6; (5) the ratio between the maximal step sizes of the large and small PN'i; (6) the difference in the distance of the large and  small PN's from zygote center. In this manner, our classifier supports MC prediction within 3 days from fertilization in parallel to prediction of embryo implantation. Training a machine learning-based classifier for assessing MC potential requires positively-implanted embryos with MC or LB indications. We used a large dataset of 70 K embryos that included 5,500 embryos with known implantation outcome, out of which 96 MC-labeled and 368 LBlabeled embryos were available for training. These proportions demonstrate the challenge in allocating the required clinically-labeled data for our retrospective multi-center study. Owing to the fact that the embryos in our dataset had been selected for transfer according to consensual morphodynamic guidelines, their visual diversity was lower than the heterogeneity that is demonstrated among transferred and non-transferred embryos. Importantly, the reduced diversity among the embryos that are available for model training is shared by other classification tasks including the prediction of embryo implantation and pregnancy outcome. Hence, our retrospective study was limited by sample size and by the morpho-dynamic heterogeneity of the embryos.
To overcome the abovementioned challenges, we performed a rigorous multistep feature screening of the embryos that minimizes the KS p-value while maximizing the KL divergence while evaluating their predictive potential using RF and XG models. We identified a subset of six morpho-dynamic features that retain prediction accuracy during early stages of preimplantation development. Interestingly, the selected subset was enriched with features that manifest the differences in the dynamic properties of PN and NPB compartments between MC and LB embryos. Since aneuploidy is tightly associated with first trimester MC [5], the contribution of PN and NPB-based features is consistent with a direct consequential impact that chromosomal abnormalities may have on PN and NPB dynamics. It would be insightful to test the capacity of the MC classifier to predict various forms of embryo aneuploidy; however, this task stems beyond the scope of this report. The fact that DLR 2,3 and DLR 5,6 features are included within the selected subset highlights potential downstream effects that chromosomal abnormalities may have on cleavage regulation in a manner that decouples early cleavage events (t2-t3) and late cleavage events (t5-t6). Feature screening as reported here is not a deterministic process and the selection of a different feature subset that is both predictive and non-redundant cannot be excluded. While the reported features were extracted using semi-automated algorithms, clinical implementation would be further supported by developing fully automated tools for the extraction of the reported features or of an equivalent subset of predictive and non-redundant features [47].
We demonstrate the prediction of embryo MC outcome with AUC 68%. This AUC value is comparable with the prediction accuracy of other reported algorithms, including CNN-based classification models of embryo implantation outcome [20,22]. The fact that the accuracy in the prediction of embryo implantation and MC outcome, as assessed via AUC, appears to have a finite upper limit may be indicative of the missing input uterine parameters and/or maternal factors [23]. In terms of clinical application, our work provides a real-time noninvasive decision-support tool that can improve the current policy in IVF treatments. Embryos will first be selected for transfer based on their implantation potential as predicted using existing classifiers. Then, the embryos with a high risk of MC will be de-selected, thus improving live-birth rates while shortening time-to-pregnancy. In conclusion, we hope that by providing a proof-of-concept of assessing the risk of MC, our retrospective study will motivate future work using larger datasets that will eventually lead to prospective clinical studies. interest. I.H.V. and A.B.M. declare having no conflict of interest during the time of scientific collaboration and data collection that are relevant to this work. Since January 2020 A.B.M. serves as the Chief Medical Officer and since March 2020 I.H.V. serves as the Scientific Director of Fairtility LTD, which is a company that incorporates AI into different stages in fertility treatment. Y.K.T. declares no financial or non-financial competing or other conflict of interest. Since September 2021, Y.K.T. discloses being employed by IBM-Research.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.