Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 22 December 2021
Sec. Human and Medical Genomics
This article is part of the Research Topic Reproductive Genomics View all 13 articles

Blastocyst Morphology Based on Uniform Time-Point Assessments is Correlated With Mosaic Levels in Embryos

  • 1Division of Infertility, Lee Women’s Hospital, Taichung, Taiwan
  • 2Institute of Medicine, Chung Shan Medical University, Taichung, Taiwan
  • 3Department of Obstetrics and Gynecology, Chung Shan Medical University Hospital, Taichung, Taiwan

Avoiding aneuploid embryo transfers has been shown to improve pregnancy outcomes in patients with implantation failure and pregnancy loss. This retrospective cohort study aims to analyze the correlation of time-lapse (TL)-based variables and numeric blastocyst morphological scores (TLBMSs) with different mosaic levels. In total, 918 biopsied blastocysts with time-lapse assessments at a uniform time-point were subjected to next-generation sequencing–based preimplantation genetic testing for aneuploidy. In consideration of patient- and cycle-related confounding factors, all redefined blastocyst morphology components of low-grade blastocysts, that is, expansion levels (odds ratio [OR] = 0.388, 95% confidence interval [CI] = 0.217–0.695; OR = 0.328, 95% CI = 0.181–0.596; OR = 0.343, 95% CI = 0.179–0.657), inner cell mass grades (OR = 0.563, 95% CI = 0.333–0.962; OR = 0.35, 95% CI = 0.211–0.58; OR = 0.497, 95% CI = 0.274–0.9), and trophectoderm grades (OR = 0.29, 95% CI = 0.178–0.473; OR = 0.242, 95% CI = 0.143–0.411; OR = 0.3, 95% CI = 0.162–0.554), were less correlated with mosaic levels ≤20%, <50%, and ≤80% as compared with those of top-grade blastocysts (p < 0.05). After converting blastocyst morphology grades into scores, high TLBMSs were associated with greater probabilities of mosaic levels ≤20% (OR = 1.326, 95% CI = 1.187–1.481), <50% (OR = 1.425, 95% CI = 1.262–1.608), and ≤80% (OR = 1.351, 95% CI = 1.186–1.539) (p < 0.001). The prediction abilities of TLBMSs were similar for mosaic levels ≤20% (AUC = 0.604, 95% CI = 0.565–0.642), <50% (AUC = 0.634, 95% CI = 0.598–0.671), and ≤80% (AUC = 0.617, 95% CI = 0.576–0.658). In conclusion, detailed evaluation with TL monitoring at the specific time window reveals that redefined blastocyst morphology components and converted numeric TLBMSs are significantly correlated with all of the threshold levels of mosaicism. However, the performance of TLBMSs to differentiate blastocysts with aberrant ploidy risk remains perfectible.

Introduction

Several studies have suggested that time-lapse (TL) monitoring can enhance aneuploid embryo identification, particularly based on blastulation kinetics (Campbell et al., 2013a; Campbell et al., 2013b; Basile et al., 2014; Minasi et al., 2016; Del Carmen Nogales et al., 2017; Mumusoglu et al., 2017; Desai et al., 2018). Campbell et al. revealed that the period from insemination to the earliest signs of compaction, to the earliest signs of cavitation (tSB), and to the blastocyst with a full-filled blastocoel (tB) are prolonged in aneuploid embryos compared with those in euploid embryos. A risk classification model of aneuploidy based on two blastulation kinetics parameters (tSB and tB) has been established, which has the potential to select embryos with increased rates of fetal heart beat and live birth during fresh embryo transfer cycles without PGT-A (Campbell et al., 2013a; Campbell et al., 2013b). Similar results were obtained from a dataset of 928-blastocysts and confirmed that euploid embryos have fast kinetics related to blastocyst development (Minasi et al., 2016). Moreover, Mumusoglu et al. (2017) further considered the confounding factors related to patient and ovarian stimulation and demonstrated that tB and the time from insemination to the expanded blastocyst stage can moderately predict aneuploidy (Mumusoglu et al., 2017). Using TL technology, dynamic dysmorphisms of cleavage-stage embryos, such as multinucleation (MN), reverse cleavage (RC), irregular chaotic division (ICD), or direct cleavage (DC), can be evaluated in more detail, and two or more dysmorphisms are correlated with a reduced euploidy probability (Desai et al., 2018). However, these studies have determined the ploidy status, at least in part if not all, of the biopsied embryos through array comparative genomic hybridization (aCGH), which has insufficient sensitivity for the detection of diploid–aneuploid mosaicism (van Echten-Arends et al., 2011; Munne and Wells 2017; Sachdev et al., 2017).

Studies have demonstrated the potential of the postimplantation development of transferred mosaic embryos (Greco et al., 2015; Munne et al., 2017; Viotti et al., 2021). Importantly, the application of high-resolution next-generation sequencing (hr-NGS) may enhance sensitivity to the level of 10 MB, through which the detailed ploidy characteristics of embryos, for example, chromosomal gains and losses, different mosaic levels, and whole-or segmental-chromosome aneuploidy, are efficiently identified (Kane et al., 2016; Munne and Wells 2017). Clinical results have demonstrated that compared with array CGH counterparts, hr-NGS significantly improves the pregnancy outcomes of single-euploid embryo transfer cycles (Friedenthal et al., 2018), and evidence supports that healthy babies can result from embryos with low or high mosaic levels (Lee et al., 2020b; Lin et al., 2020). However, the clinical outcomes of mosaic embryos appear to be significantly affected by mosaic levels (<50% vs. ≥50%), mosaic types (chromosome numbers), and blastocyst morphology (Viotti et al., 2021). It would be interesting to evaluate the correlation of the mosaic status with embryonic characteristics. To date, whether any embryonic feature has the ability to distinguish embryos with different mosaicism levels remains to be investigated.

Numerous experimental studies have delineated a link between blastocyst morphology and chromosomal abnormalities (Alfarawati et al., 2011; Capalbo et al., 2014; Minasi et al., 2016; Barash et al., 2017; Wang et al., 2018; Zhan et al., 2020). Individual morphological components of blastocyst evaluation, that is, expansion degrees, inner cell mass (ICM) grades, and trophectoderm (TE) grades, may be correlated with the ploidy status independently (Minasi et al., 2016; Zhan et al., 2020). The combination of expansion, ICM, and TE grades could therefore result in complication for ranking the embryos. Thus, several studies have combined morphological components to categorize blastocysts based on quality (Capalbo et al., 2014; Wirleitner et al., 2016) and have even considered the blastocyst biopsy day and morphological features for generating a blastocyst score (Zhan et al., 2020). Although these modified approaches demonstrate that blastocyst morphology is highly associated with its ploidy (Capalbo et al., 2014; Zhan et al., 2020), embryo grading variability may occur due to multiple reasons, such as subjective judgments of individual embryologists, laboratory settings, and observation time windows. One of the advantages of TL monitoring is objectively and accurately defined observation time after insemination or ICSI. The hypothesis was therefore raised that blastocyst morphology observed at a definite time period by TL monitoring might be precisely correlated with the mosaic level in an embryo. If the hypothesis is correct, such blastocyst morphology scores may be used to improve embryo selection without the help of AI or labor-intensive interpretation of the whole TL videos.

In the current study, hr-NGS was applied to classify the mosaic level of a blastocyst. We refined embryo assessments at a uniform time-point by using the definition based on TL monitoring and aimed to continuously improve our understanding regarding their relationship with different mosaic levels. In addition, the effectiveness of comprehensive assessment, which converted blastocyst morphology into a numeric score, in predicting ploidy was evaluated by considering the effects of potential confounding variables.

Materials and Methods

Study Design and Patient Selection

The data were collected at Lee Women’s Hospital from January 2017 to August 2018. In total, 210 infertile couples undergoing PGT-A and TL cultivation were enrolled into this retrospective cohort study, which was approved by the Institutional Review Board of Chung Sun Medical University, Taichung, Taiwan (approval number CS19039). Patients with surgical sperm retrieval (n = 4) and without qualified blastocysts for PGT-A (n = 28) were excluded. A total of 918 biopsied blastocysts from 178 couples were analyzed in this study.

Controlled Ovarian Stimulation

All the procedures were conducted following the standard protocols stated in our previous report (Lee et al., 2019). Briefly, in this study, controlled ovarian stimulation was applied using the gonadotrophin-releasing hormone (GnRH) agonist long protocol or the GnRH antagonist protocol. Patients who were assigned the GnRH agonist long protocol received leuprolide acetate (Lupron; Takeda Chemical Industries, Osaka, Japan) subcutaneously (0.05–0.1 mg/day) from day 21 of the first cycle to the day of human chorionic gonadotropin (hCG; Ovidrel; Serono, Modugno, Italy) injection (250 μg). Exogenous gonadotropin (Gonal-F; Serono or Menopur; Ferring, São Paulo, Brazil) was then administered from day 3 of the second cycle until the leading follicle was ≥18 mm. Patients who were assigned the GnRH antagonist protocol received the gonadotropin regimen, same as that in the agonist long protocol, and additional cetrorelix acetate (Cetrotide; Merck Serono, Geneva, Switzerland) injections (0.25 mg/day) were given subcutaneously before hCG injection. Oocyte maturation was triggered through hCG administration, and ultrasound-guided ovum retrieval was performed approximately 36 h after hCG injection (Lee et al., 2020a).

Cultivation, Insemination, and Micromanipulation

The TL culture system, EmbryoScope+ (Vitrolife, Kungsbacka, Sweden), was used in this study. The cultivation condition consisting of an atmosphere of 5% O2, 5% CO2, and 90% N2 at 37°C was applied for oocytes and embryos. The media covered with paraffin oil (Vitrolife) were equilibrated in an incubator for at least 2 h before use, including the fertilization medium (SAGE Biopharma, Bedminster, NJ, United States) with 15% serum protein substitute (SPS; SAGE Biopharma) for conventional insemination or intracytoplasmic sperm injection (ICSI), cleavage medium (SAGE Biopharma) with 15% SPS, and blastocyst medium (SAGE Biopharma) with 15% SPS for embryo culture. On the basis of the World Health Organization Laboratory Manual for the Examination and Processing of Human Semen, abnormal sperm quality was defined if at least one of the following properties was found: oligospermia (sperm concentration ≤15 × 106/ml), asthenospermia (total motility ≤40%), and teratospermia (normal sperm form ≤4%). ICSI was used for all of the mature oocytes from couples with oligoasthenoteratozoospermia (OAT) or for approximately half of the mature oocytes from couples with non-OAT.

Time-Lapse Assessment and Morphological Grading

Embryo morphology and developmental dynamics were noninvasively observed through the capture of images with 11 focal planes at 10-min intervals using a TL device. According to guidelines by Vitrolife and published nomenclature (Ciray et al., 2014), the exact timings of developmental events and morphological features after complete pronuclei fading (tPNf) were annotated on a daily basis. Supplementary Table S1 provides the definitions of blastulation-related kinetics, including the time to accomplish compaction (tM), tSB, tB, the time interval between tM and tB (tB-tM), and the time interval between tSB and tB (tB-tSB), as well as cleavage aberrations at early cleavage stages (≤8-cell stage), including delayed division (DD), DC, RC, ICD, MN at the two-cell stage (MN2), and MN at the four-cell stage (MN4). In order to synchronized the observation time window of individual blastocysts, the morphological evaluation was performed at 118 h post insemination (hpi) by all of recorded images of embryonic development using redefined descriptions of embryo expansion, TE quality, and ICM quality through TL monitoring (Supplementary Figure S1 and Supplementary Table S2). A score from 1 to 7 was then assigned to the expansion level (level of <1 = 1, level 1 = 3, level 2 = 5, and level 3 = 7), and a score from 0 to 2 was assigned to ICM and TE groups (grade A = 2, grade B = 1, and grade C or less = 0). The TL-based blastocyst morphological score (TLBMS) was then developed to represent the morphological grading of a blastocyst using the following formula: expansion score + (ICM score × TE score).

Determination of Mosaic Levels

After in vitro cultivation, embryos were assigned into different groups according to morphological assessments at 118 hpi, that is TLBMS, embryo expansion, TE quality, and ICM quality. The expanded blastocysts with at least grade B of ICM or at least grade B of TE on day 5 (around 120 hpi) or day 6 (around 140 hpi) were selected for embryo biopsy. Five to eight TE cells were then carefully separated from a blastocyst through micromanipulation by using inverted microscopy with a laser system (Chen et al., 2017). Biopsied cells were rinsed with phosphate-buffered saline first and then placed into an RNAse–DNAse-free polymerase chain reaction tube. The remaining blastocyst was incubated in a tri-gas incubator for ≥3 h and then subjected to cryopreservation by using the Cryotech vitrification method (Cryotech, Tokyo, Japan). The mosaic levels of biopsied samples were determined according to the manufacturer’s instructions of hr-NGS, including the SurePlex DNA amplification system (Illumina, San Diego, CA, United States) for the extraction and preparation of genomic DNA of TE cells, the VeriSeq PGS workflow (Illumina) for the preparation of genomic DNA libraries, the VeriSeq DNA Library Prep Kit (Illumina) for the normalization of each DNA library’s concentration, and the Miseq System with Miseq Reagent Kit v3 (Illumina) for DNA sequencing of individual libraries. Finally, the generated bioinformatic data were analyzed using BlueFuse Multi software (Illumina), and segmental changes in individual chromosomes were defined as altered segments of ≥10 MB in size. The mosaic level was determined by the percentage of aneuploid cells in a TE biopsy specimen. If the embryo had more than one mosaic chromosomal region, the highest value of mosaic levels was used for analysis (Viotti et al., 2021). The ploidy status of each sample was determined according to the following criteria: 1) mosaic levels ≤20% (euploidy); 2) mosaic levels between >20% and <50% (low-level mosaicism); 3) mosaic levels between ≥50% and ≤80% (high-level mosaicism); and 4) mosaic levels >80% (aneuploidy). The hr-NGS data for all embryos have been deposited to the NCBI SRA database, and the BioProject accession number is PRJNA782420.

Statistical Analysis

SPSS (v 20.0; IBM Corporation, United States) or Prism (version 6.0 h; GraphPad Software, United States) were used for statistical analysis. p < 0.05 was considered significant in all analyses. Generalized estimating equation (GEE) analysis was used to assess the associations between the probability of a specific ploidy status and individual independent variables in unadjusted (univariate) or adjusted (multivariate) models. The ploidy status examined included three thresholds of mosaic levels, that is, ≤20%, <50%, and ≤80%. The receiver operating characteristic (ROC) curve was applied to estimate the prediction performance of TLBMS and determine cutoff values. Finally, the rates of euploidy, low-level mosaicism, high-level mosaicism, and aneuploidy between groups were compared using Fisher’s exact test.

Results

The patient- and cycle-characteristics were described in Table 1. All morphokinetic and morphological features of biopsied blastocysts (n = 918) with known mosaic levels, comprising 320 euploid blastocysts, 242 low-level mosaic blastocysts, 111 high-level mosaic blastocysts, and 245 aneuploid blastocysts, were evaluated using TL images. The potential confounding variables associated with embryo ploidy were analyzed, including factors related to patient (at the start of infertility treatment) or cycle characteristics, that is female age, anti-Müllerian hormone (AMH), body mass index (BMI), male age, oocyte numbers, mature oocyte numbers, ovarian stimulation protocols (agonist long vs. antagonist protocols), oocyte sources (autologous vs. donor oocytes), sperm quality (abnormal vs. normal sperm), and artificial insemination methods (ICSI vs. conventional insemination); factors related to the kinetics of blastocoel formation, that is, tM, tSB, tB, tB-tM, and tB-tSB; factors related to embryo dysmorphisms, that is, DD, DC, RC, ICD, MN2, and MN4; and factors related to blastocyst morphology, that is, expansion levels, ICM quality, and TE quality.

TABLE 1
www.frontiersin.org

TABLE 1. Patient and cycle characteristics.

Association Between Patient- or Cycle-Related Confounding Variables and Mosaic Levels

Firstly, this study aimed to identify significant confounding factors derived from non-embryonic sources. The correlation between individual patient confounding factors and mosaic levels was tested through univariate regression analysis by using the GEE model (Table 2). Data revealed that female age (odds ratio [OR] = 0.956, 95% confidence interval [CI] = 0.931–0.981, p < 0.05) was associated with mosaic levels ≤20%; female age (OR = 0.932, 95% CI = 0.905–0.96, p < 0.001), mature oocyte numbers (OR = 1.017, 95% CI = 1–1.033, p < 0.05), and oocyte sources (autologous vs. donor oocytes, OR = 0.57, 95% CI = 0.332–0.976, p < 0.05) were significantly associated with mosaic levels <50%; and female age (OR = 0.92, 95% CI = 0.884–0.958, p < 0.001), oocyte numbers (OR = 1.021, 95% CI = 1.004–1.038, p < 0.05), mature oocyte numbers (OR = 1.024, 95% CI = 1.004–1.044, p < 0.05), and oocyte sources (autologous vs. donor oocytes, OR = 0.351, 95% CI = 0.157–0.785, p < 0.05) were significantly associated with mosaic levels ≤80%. According to the results, female age, oocyte numbers, mature oocyte numbers, and oocyte sources were considered to be significant confounding factors from non-embryonic sources.

TABLE 2
www.frontiersin.org

TABLE 2. The correlations between patient- or cycle-characteristics and embryo ploidy status in the current dataset.

Correlation Between Adjusted Embryo Variables and Mosaic Levels

This study attempted to evaluate the associations between embryonic variables and ploidy status. Considering the significant patient- and cycle-related confounding factors, that is, female age, mature oocyte numbers, and oocyte sources, the multivariate GEE analysis demonstrated that blastocyst expansion levels (level ≤1 vs. level 3, OR = 0.388, 95% CI = 0.217–0.695, p < 0.005; level 2 vs. level 3, OR = 0.488, 95% CI = 0.347–0.688, p < 0.001), ICM grades (grade C or less vs. grade A, OR = 0.563, 95% CI = 0.333–0.962, p < 0.05), and TE grades (grade C or less vs. grade A, OR = 0.29, 95% CI = 0.178–0.473, p < 0.001; grade B vs. grade A, OR = 0.586, 95% CI = 0.39–0.883, p < 0.05) were associated with mosaic levels ≤20%. Furthermore, tSB (OR = 0.978, 95% CI = 0.961–0.996, p < 0.05), tB (OR = 0.975, 95% CI = 0.96–0.991, p < 0.005), tB-tM (OR = 0.957, 95% CI = 0.934–0.98, p < 0.001), blastocyst expansion levels (level ≤1 vs. level 3, OR = 0.328, 95% CI = 0.181–0.596, p < 0.001; level 2 vs. level 3, OR = 0.48, 95% CI = 0.323–0.712, p < 0.001), ICM grades (grade C or less vs. grade A, OR = 0.35, 95% CI = 0.211–0.58, p < 0.001; grade B vs. grade A, OR = 0.604, 95% CI = 0.413–0.882, p < 0.01), and TE grades (grade C or less vs. grade A, OR = 0.242, 95% CI = 0.143–0.411, p < 0.001; grade B vs. grade A, OR = 0.584, 95% CI = 0.358–0.953, p < 0.05) were associated with mosaic levels <50%. Finally, tSB (OR = 0.977, 95% CI = 0.957–0.997, p < 0.05), tB (OR = 0.973, 95% CI = 0.955–0.991, p < 0.005), tB-tM (OR = 0.942, 95% CI = 0.918–0.966, p < 0.001), tB-tSB (OR = 0.963, 95% CI = 0.928–0.999, p < 0.05), MN4 (MN4 vs. non-MN4, OR = 1.655, 95% CI = 1.004–2.729, p < 0.05), blastocyst expansion levels (level ≤1 vs. level 3, OR = 0.343, 95% CI = 0.179–0.657, p < 0.005; level 2 vs. level 3, OR = 0.523, 95% CI = 0.327–0.835, p < 0.01), ICM grades (grade C or less vs. grade A, OR = 0.497, 95% CI = 0.274–0.9, p < 0.05), and TE grades (grade C or less vs. grade A, OR = 0.3, 95% CI = 0.162–0.554, p < 0.001) were associated with mosaic levels ≤80% (Table 3). The results indicated that the grading of blastocyst morphology components were positively associated with all of threshold levels of mosaicism. However, the variables of blastulation kinetics, such as tB, tSB, tB-tM, and tB-tSB, were negatively associated with the mosaic level <50% or ≤80%, and the variable of embryo dysmorphisms, i.e., MN4, had a positive association with the mosaic level ≤80%. Additionally, when incorporating the blastocyst morphology component (e.g., expansion levels, ICM grades, or TE grades) with tB and MN4 in the multivariate logistic regression models, the associations between tB and embryo ploidy were attenuated to be non-significant and MN4 was merely correlated with the mosaic level ≤80% (OR = 1.703–2.16, p < 0.05) (Supplementary Tables S3–S5). The non-aneuploid (mosaic level ≤80%) rate of the blastocysts with MN4 (81.5%) was higher than the blastocysts without MN4 (72.2%) (p < 0.05), which resulted from the significantly raised rate of high-level mosaicism in the MN4 group (19.4%) (Supplementary Figure S2A). By contrast, the blastocyst morphology components still had the statistically significant correlations with mosaic levels ≤20% (OR = 0.225–0.528), <50% (OR = 0.204–0.614) or ≤80% (OR = 0.261–0.554), indicating this redefined blastocyst morphology could be capable of categorizing the embryos with different mosaic levels (Supplementary Figure S3 and Supplementary Tables S3–S5).

TABLE 3
www.frontiersin.org

TABLE 3. The correlations between embryonic variables and embryo ploidy status in consideration of patient- and cycle-confounding factors.

Effectiveness of Time Lapse-Based Blastocyst Morphological Score as a Predictor for Embryo Ploidy

In order to provide a comprehensive inspection of blastocyst morphology, this study tried to convert the grading of morphological components into TLBMSs and investigated its correlations with different mosaic levels. After adjustment of patient- and cycle-related confounding factors (i.e., female age, mature oocyte numbers, and oocyte sources) as well as the embryo features except for blastocyst morphology, i.e., tB and MN4, the results revealed that TLBMS was significantly associated with mosaic levels ≤20% (OR = 1.326, 95% CI = 1.187–1.481, p < 0.001), <50% (OR = 1.425, 95% CI = 1.262–1.608, p < 0.001), and ≤80% (OR = 1.351, 95% CI = 1.186–1.539, p < 0.001) (Table 4). The ROC curve demonstrated that TLBMS was statistically associated with mosaic levels. The values of area under the curve (AUC) were 0.604, 0.634, and 0.617 for mosaic levels ≤20%, <50%, and ≤80%, respectively (Table 4). The cutoff values of TLBMS and female age were 6 and 36 years, respectively. Embryos were classified into four group according to cutoff values. Between the groups with TLBMSs ≥6 and <6, the rates of mosaic levels ≤20% (41.4 vs. 30.3% in the age group of <36 years, p < 0.05; 39 vs. 16.8% in the age group of ≥36 years, p < 0.0001), mosaic levels <50% (72.8 vs. 55.1% in the age group of <36 years, p < 0.001; 63.7 vs. 35.5% in the age group of ≥36 years, p < 0.0001), and mosaic levels ≤80% (83.4 vs. 72.4% in the age group of <36 years, p < 0.01; 72.9 vs. 51% in the age group of ≥36 years, p < 0.0001) appeared to be different (Figure 1). Moreover, this study further combined TLBMSs with the occurrence of MN4 for embryo selection. The results demonstrated that a higher rate of high-level mosaicism was observed in the good morphology blastocysts (TLBMSs ≥6) with MN4 (20%) than the good morphology blastocysts without MN4 (9.3%) (p < 0.05), suggesting deselection of MN4 embryos could be considered to reduce the risk of selecting good morphology blastocysts with high-level mosaicism (Supplementary Figure S2D).

TABLE 4
www.frontiersin.org

TABLE 4. Evaluation of the relationship with embryo ploidy status and prediction abilities of time-lapse blastocyst morphology scores.

FIGURE 1
www.frontiersin.org

FIGURE 1. Distribution of blastocysts with different mosaic levels according to time-lapse blastocyst morphology scores. Fisher exact test was used for statistical analysis and women were stratified into age <36 years and ≥36 years groups. Same superscript letters indicated statistically significant. Abbreviations “n”, “yrs”, and “TLBMS” denoted the number of embryos, years, and time-lapse blastocyst morphology score, respectively.

Discussion

As mentioned in a literature review, the diploid–aneuploid mosaicism in an embryo can originate from paternal, maternal, and external factors (Taylor et al., 2014). Firstly, this retrospective study was proposed to identify the possible confounding factors from the non-embryonic source. In our dataset, the univariate GEE demonstrated that several patient or cycle characteristics could be the confounding factors of mosaic levels, that is, female age, oocyte sources, and numbers of oocytes and mature oocytes, in addition to embryo parameters (Table 2). Several reports have shown that maternal age is considerably associated with oocyte aneuploidy incidence (Jones 2008; Gianaroli et al., 2010; Greaney et al., 2018; Gruhn et al., 2019). Gruhn et al. (2019) demonstrated that the proportion of aneuploid oocytes follows a U-shaped curve according to female age, in which the risk of meiotic segregation errors is decreased in the 20–32 years age group (Gruhn et al., 2019). Similarly, the current study revealed that female age (23–44 years) was the major confounding factor, which was negatively associated with mosaic levels at all three thresholds (≤20%, <50%, or ≤80%). Furthermore, advanced age was one of the reasons for the association of autologous oocytes and reduced oocyte numbers with increased mosaic levels in our dataset. Taking female age into account revealed that oocyte numbers and oocyte sources were no longer significantly associated with ploidy status (Supplementary Table S6). Schaeffer et al. (2020) analyzed aneuploidy rates and IVF outcomes by comparing patient and donor, showing that with increasing maternal age there are more substantial differences in aneuploidy rates between the patient and the donor embryos. Moreover, Irani et al. (2020) reveals that euploidy rates are similar between groups with different numbers of retrieved oocytes after age stratification. Notably, the present study showed that male factors, that is, male age and semen quality, were not significantly related to mosaic levels. These results are supported by the findings of a previous study that used young donor oocytes for analysis, which revealed that advanced paternal age was not associated with higher aneuploidy rates (Dviri et al., 2021). Moreover, as per blastocysts obtained, sperm factors did not exert significant effects on the euploidy rate, although using poor-quality sperm for insemination appears to reduce fertilization and blastocyst formation rates (Mazzilli et al., 2017). Furthermore, univariate GEE analysis demonstrated that other patient or cycle variables evaluated in this study, such as BMI, AMH, ovarian stimulation protocols, and artificial insemination methods, were not significantly correlated with embryo ploidy; these findings were similar to previous reports (Deng et al., 2020; Irani et al., 2020; Cascales et al., 2021; Pipari et al., 2021; Stovezky et al., 2021).

This study was designed to assess the importance of embryonic characteristics in embryo ploidy. Considering confounding variables from nonembryonic factors, that is, maternal age, mature oocyte numbers, and oocyte sources, the present study revealed that delayed blastocyst formation was positively associated with increased mosaic levels (Table 3). Two studies, which applied TL monitoring and hr-NGS simultaneously for PGT-A cycles, observed that aneuploid blastocysts exhibited a significant delay in tSB or tB as compared with euploid blastocysts (Lee et al., 2019; Martin et al., 2021). Multivariate GEE analysis in this study further demonstrated a meaningful correlation of tSB and tB with mosaic levels when the threshold was set at 50%. Correspondingly, prolonged blastulation intervals (tB-tM or tB-tSB) accompanied with mosaic levels ≥50% and >80% were present. However, this correlation was not detected between blastocyst kinetics and euploidy (mosaic levels ≤20%). A consistent observation with our previous study (Lee et al., 2019) revealed that blastulation kinetics of low-level mosaic blastocysts are comparable with those of euploid blastocysts (Supplementary Table S7). When low-level mosaic blastocysts were included in the reference category (mosaic level >20%) of a binary logistic regression, the associations of blastocyst kinetics with mosaic level ≤20% appeared to be non-significant.

On the other hand, considering patient- and cycle-related confounding factors, the results of this study indicated that the correlation of individual developmental dysmorphisms with embryo ploidy was not identified, except for MN4 (Table 3). Desai et al. (2018) reported a similar influence of developmental dysmorphisms and concluded that individual dysmorphisms, including MN, RC, DC, and ICD, are not associated with euploidy (Desai et al., 2018). In their study, MN embryos were examined at two-to five-cell stages, with the pooling of MN2 and MN4 embryos for euploidy analysis by using aCGH or hr-NGS, which might lead to an inconsistent result for MN embryos as compared with our study. A mouse model was used to prove that MN occurrence at the early developmental stage affects blastomere ploidy and compromises blastocyst developmental potential. Notably, after MN-derived blastocysts were transferred to surrogate mothers, the pregnancy loss rate was not significantly increased by the adverse effects of early chromosome segregation errors (Mashiko et al., 2020). A drastic decrease in MN incidence from two-to four-cell stages and altered morphokinetic features, such as delayed timings of early divisions and prolonged periods of two-cell and four-cell stages, support that a repair mechanism exists for self-correction of MN embryos during early cell divisions (Balakier et al., 2016). According to our data, in biopsied blastocysts, the MN incidence was decreased from the two-cell stage (23.2%, 213/918) to the four-cell stage (11.8%, 108/918), suggesting a repair mechanism might exist during early cell divisions (Supplementary Table S7). However, if the early chromosome segregation errors continued at the four-cell stage, self-correction of MN embryos might not be completed in time, resulting in the high potential of high-level mosaicism (Supplementary Figure S2 and Supplementary Table S7).

The TL incubation system is designed to enable four-dimension-like observation, which offers a comprehensive and consistent inspection of embryo morphology and dynamics. One of the aims of this study was to modify the morphological evaluation method of a blastocyst at a specific time window (118 hpi) based on traditional assessments (Gardner et al., 2000) and TL monitoring requirements (Kragh et al., 2019). In our dataset, the fertilization methods (ICSI vs. conventional insemination) appeared to have non-significant associations with the timing of a blastocyst with a cavity beginning to form (OR = 1.338, 95% CI = 0.362–4.946, p > 0.05) and the timing of a blastocyst with a cavity beginning to expand (OR = 1.46, 95% CI = 0.362–5.891, p > 0.05). The same observation time-point was therefore applied for the blastocysts derived from ICSI or conventional insemination. The annotations of ICM and TE were conducted more thoroughly to generalize the Gardner blastocyst grading system to time-lapse imaging (Kragh et al., 2019). The annotations of expansion levels emphasized the specific developmental features, including an embryo with the cavity beginning to form (level 1), to expand (level 2), and to herniate (level 3), which were difficult to determine with the conventional method (Supplementary Figure S1). This study found that TL-based morphological components are still considerably associated with mosaic levels at all thresholds (Table 3). Unlike in previous studies (Capalbo et al., 2014; Minasi et al., 2016; Wang et al., 2018), the expansion levels of blastocysts in this study were evaluated at a uniform time window (118 hpi), which could be considered an evaluation of the developmental speed; blastocysts with expansion levels ≤1 (OR = 0.388, 0.328, and 0.343) and 2 (OR = 0.488, 0.480, and 0.523) had significantly lower probabilities of mosaic levels ≤20%, <50%, and ≤80% compared with blastocysts with expansion level 3. The expansion results of TL-based assessments not only revealed the correlation of the blastocyst expansion status with embryo ploidy (Minasi et al., 2016; Wang et al., 2018), but also reflected the importance of the developmental speed on blastocyst formation with lower mosaic levels, which had been evidenced by morphokinetic observations in the current and previous studies (Campbell et al., 2013a; Mumusoglu et al., 2017). On the other hand, in accordance with previous studies (Fragouli et al., 2014; Minasi et al., 2016; Wang et al., 2018), this study indicated a greater likelihood of euploidy among blastocysts with good quality ICM or TE. However, in previous studies, under-sensitive PGT-A approaches and unfixed periods of ICM and TE assessments, for example, D5, D6, or even D7, may result in inconsistent findings. Therefore, in our study, multivariate GEE analysis was conducted to evaluate the associations of ICM and TE quality, which was assessed at the specific time window of cultivation, with embryo ploidy using hr-NGS. In this circumstance, as compared with blastocysts with grade A of ICM or TE, the blastocysts with grade C or less of ICM or TE had a higher risk to be non-euploid (>20%), high-level mosaic (≥50%), and aneuploid (>80%) embryos. Similar risks of non-euploidy or aneuploidy were found between blastocysts with grade A and grade B of ICM, but the blastocysts with grade B of ICM were still at an increased risk of mosaic levels ≥50%. Moreover, similar risks were only observed in aneuploidy between blastocysts with grade A and grade B of TE. These results suggest that the morphology of TE and ICM influenced, at least, the discrimination between mosaic levels <50% and ≥50% among blastocysts with moderate to high quality of ICM or TE (Table 3). This study thus demonstrated that all the three refined morphology criteria of a blastocyst were meaningfully correlated with embryo ploidy (Supplementary Figure S3).

In order to make an effort to improve intuitive decision-making in clinical practice, the current study designed a TLBMS formula as a ranking tool of blastocyst morphology. Furthermore, patient- and cycle-related factors, that is, female age, mature oocyte numbers, and oocyte sources, and significant embryo factors, that is, tB and MN4, were considered to be confounding factors and adjusted in GEE analysis, which revealed an apparent correlation of TLBMS with embryo ploidy (Table 4). Previously, at least two functional formulae were available, which used multiplication or addition for the calculation of parameters from expansion levels, ICM grades, and TE grades to convert blastocyst morphology grades into numbers (Rehman et al., 2007; Zhan et al., 2020). Compared with our study, Rehman et al. (2007) used the multiplication formula to generate a wide range of blastocyst quality scores (BQSs = 1–54, a total of 17 scores), leading to the nonlinear correlation with pregnancy outcomes. The blastocyst biopsy day was not considered in the formula for BQS, and the correlation of BQS with embryo ploidy was not investigated (Rehman et al., 2007). Zhan et al. (2020) further improved the formula by using addition instead of multiplication for calculation, which narrows down the score range (score 3–15, a total of 13 scores) and ranks blastocysts considering the biopsy day. Each grade of the expansion status, ICM quality, and TE quality was assigned a score based on the implantation rate, resulting in a linear relationship between blastocyst morphology scores and euploidy rates and acceptable prediction abilities (AUC = 0.683–0.698) (Zhan et al., 2020). However, this updated BMS uses a complicated grading system for expansion (nine levels), ICM (seven grades), and TE (seven grades), which may increase participant bias by embryologists. Inconsistent results of the embryo ploidy status might appear when different PGT-A platforms were used in a study (Zhan et al., 2020). Applying this scoring system to our dataset revealed slightly lower ploidy prediction abilities (AUC = 0.585 for mosaic levels ≤20%, AUC = 0.628 for mosaic levels <50%, and AUC = 0.613 for mosaic levels ≤80%; data not shown). Distinctly, our TLBMS considering two separated parts of blastocyst morphology, that is, the expansion level and ICM/TE quality, used a simplified formula with a combination of addition and multiplication; a score from expansion levels with distinct definition in TL monitoring (four scores) and a score from ICM and TE combinations (four scores) were added. The expansion levels were considered to be the major part of our scoring system. Owing to only three expansion levels in this study, the scores differed by 2 each time to enlarge the gap between two continues levels. The blastocysts with the expansion level <2 would only had the scores from expansion levels and the minimum score of 1 was thus used for the expansion level of <1 to avoid a given total score of TLBMS to be zero. Additionally, this study revealed that the blastocysts with a low-quality ICM or a low-quality TE had reduced probabilities to be euploid, euploid and low-level mosaic, or non-aneuploid blastocysts. A score of zero was therefore assigned to the grade C ICM and the grade C TE. The combination score, which was generated by multiplication of a ICM score by a TE score, was used to represent the ICM/TE quality. Consequently, the blastocysts with the ICM of grade C or the TE of grade C (i.e., AC, BC, CB, and CA) would have a score of zero, the BB quality blastocysts had the score of 1, the AB or BA quality blastocysts had the score of 2, and the AA quality blastocysts had the best score of 4 (Supplementary Table S2). As compared with previous studies, TLBMS (score 1–11, a total of 9 scores) provided comparable or even better differentiation capabilities for different thresholds of mosaic levels (AUC = 0.604–0.634; Table 4). Notably, when incorporating TLBMSs with the occurrence of MN4 for embryo selection, this study revealed the high-level mosaicism rate of good morphology blastocysts (TLBMSs ≥6) with MN4 was higher than those without MN4 (Supplementary Figure S2). Compared with low-level mosaic blastocysts, previous studies have demonstrated high-level mosaic blastocysts have adverse pregnancy outcomes (e.g., decreased rates of implantation, ongoing pregnancy, and birth) (Viotti et al., 2021) and an increased miscarriage rate (Lin et al., 2020). In order to reduce the possibility of selecting blastocysts with high-level mosaicism, this study thus suggested to lower the priority of good morphology blastocysts with MN4 for embryo selection.

This single-center study had several limitations. The major limitation was its retrospective nature, which lacked randomization and might have resulted in the risk of selection bias. Therefore, this study collected biopsied blastocysts from female patients with a wide age range (23–44 years) for the analysis and took the possible confounding factors into account. Although the sample size of this study (918 blastocysts from 178 patients) was relatively small as compared with a previous study (Zhan et al., 2020), the GEE model was successful applied for analysis of repeated measurements. In comparison with a previous investigation with a large dataset (Munne et al., 2017), a higher mosaicism rate was revealed in this study (38.5 versus 21.8%), which could result from different clinical and laboratory settings between infertility centers (Sachdev et al., 2017; Lee et al., 2019). Even though the mosaicism detected by NGS was applied in donor oocyte cycles, the rates between individual infertility centers still range widely from 17 to 47% (Sachdev et al., 2017). Further investigation would be conducted to evaluate the effects of clinical settings and laboratory techniques on embryonic mosaicism. Multilayered TL images were evaluated at 118 hpi to reach a consistent observation regarding blastocyst morphology, but the implementation of TL technology might increase the financial cost and need specific training. Moreover, the current method was established based on our protocols. Further considerations on the clinical settings of individual IVF laboratories were required prior to application of TLBMS. Nevertheless, the present study was designed to determine the association of blastocyst morphology with embryo ploidy. The derived TLBMS as a ranking tool was positively associated with mosaic levels ≤20%, <50%, and ≤80% (p < 0.001) (Table 4), which could potentially select a blastocyst with a better ploidy status for both young (<36 years) and older (≥36 years) women undergoing infertility treatment (Figure 1). Embryonic factors, such as developmental potency and embryo ploidy, have major influences on the following successful pregnancy. Our previous study demonstrated that the improved implementation rate of the euploid embryo is accompanied by an increased score of the KIDScore D5 algorithm, suggesting the potential to discriminate euploid blastocysts with different developmental competence (Lee et al., 2019). The clinical application of TLBMS could be advised as follows: 1) select the embryo (e.g., TLBMS ≥6): with a relatively low risk of aneuploidy for PGT-A; 2) reduce the biopsy number if the patient demanded reduction in the financial burden; 3) directly combine the KIDScore D5 algorithm in non-PGT cycles; and 4) enhance the performance of noninvasive ploidy evaluation, such as the artificial intelligence approach (Lee et al., 2021).

Conclusion

Time-lapse assessments at the uniform time-point demonstrate significant correlations of embryonic variables derived from blastulation kinetics, cleavage dysmorphisms, and blastocyst morphology with ploidy status. However, only the morphological variables of a blastocyst (i.e., expansion levels, ICM grades, and TE grades) are meaningfully associated with all of the threshold levels of mosaicism. The TLBMS that is converted from blastocyst morphological components appears to be negatively correlated with aberrant ploidy status but the predictive performance remains limited.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: http://www.ncbi.nlm.nih.gov/bioproject/ under the accession number PRJNA782420

Ethics Statement

The studies involving human participants were reviewed and approved by Institutional Review Board of Chung Sun Medical University. The ethics committee waived the requirement of written informed consent for participation.

Author Contributions

M-SL and T-HL formulated this study. C-HC, C-IL, C-CH, H-HC, S-TH, E-HC, C-IC and P-YL collected and processed the data. C-HC, C-IL, and H-HC carried out analyses. C-HC and C-IL wrote the manuscript. All authors reviewed the manuscript and provided editorial feedback.

Funding

This study was generously supported by grants (NSC 101-2314-B-040-007 and MOST-110-2314-B-040-005) from the Ministry of Science and Technology, Executive Yuan, Taiwan, Republic of China.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.783826/full#supplementary-material

References

Alfarawati, S., Fragouli, E., Colls, P., Stevens, J., Gutiérrez-Mateo, C., Schoolcraft, W. B., et al. (2011). The Relationship between Blastocyst Morphology, Chromosomal Abnormality, and Embryo Gender. Fertil. Sterility 95, 520–524. doi:10.1016/j.fertnstert.2010.04.003

CrossRef Full Text | Google Scholar

Balakier, H., Sojecki, A., Motamedi, G., and Librach, C. (2016). Impact of Multinucleated Blastomeres on Embryo Developmental Competence, Morphokinetics, and Aneuploidy. Fertil. Sterility 106, 608–614. doi:10.1016/j.fertnstert.2016.04.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Barash, O. O., Ivani, K. A., Willman, S. P., Rosenbluth, E. M., Wachs, D. S., Hinckley, M. D., et al. (2017). Association between Growth Dynamics, Morphological Parameters, the Chromosomal Status of the Blastocysts, and Clinical Outcomes in IVF PGS Cycles with Single Embryo Transfer. J. Assist. Reprod. Genet. 34, 1007–1016. doi:10.1007/s10815-017-0944-0

CrossRef Full Text | Google Scholar

Basile, N., Nogales, M. d. C., Bronet, F., Florensa, M., Riqueiros, M., Rodrigo, L., et al. (2014). Increasing the Probability of Selecting Chromosomally normal Embryos by Time-Lapse Morphokinetics Analysis. Fertil. Sterility 101, 699–704. doi:10.1016/j.fertnstert.2013.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Campbell, A., Fishel, S., Bowman, N., Duffy, S., Sedler, M., and Hickman, C. F. L. (2013a). Modelling a Risk Classification of Aneuploidy in Human Embryos Using Non-invasive Morphokinetics. Reprod. BioMedicine Online 26, 477–485. doi:10.1016/j.rbmo.2013.02.006

CrossRef Full Text | Google Scholar

Campbell, A., Fishel, S., Bowman, N., Duffy, S., Sedler, M., and Thornton, S. (2013b). Retrospective Analysis of Outcomes after IVF Using an Aneuploidy Risk Model Derived from Time-Lapse Imaging without PGS. Reprod. BioMedicine Online 27, 140–146. doi:10.1016/j.rbmo.2013.04.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Capalbo, A., Rienzi, L., Cimadomo, D., Maggiulli, R., Elliott, T., Wright, G., et al. (2014). Correlation between Standard Blastocyst Morphology, Euploidy and Implantation: an Observational Study in Two Centers Involving 956 Screened Blastocysts. Hum. Reprod. 29, 1173–1181. doi:10.1093/humrep/deu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Cascales, A., Lledó, B., Ortiz, J. A., Morales, R., Ten, J., Llácer, J., et al. (2021). Effect of Ovarian Stimulation on Embryo Aneuploidy and Mosaicism Rate. Syst. Biol. Reprod. Med. 67, 42–49. doi:10.1080/19396368.2020.1850908

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H.-H., Huang, C.-C., Cheng, E.-H., Lee, T.-H., Chien, L.-F., and Lee, M.-S. (2017). Optimal Timing of Blastocyst Vitrification after Trophectoderm Biopsy for Preimplantation Genetic Screening. PLoS One 12, e0185747. doi:10.1371/journal.pone.0185747

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciray, H. N., Campbell, A., Agerholm, I. E., Aguilar, J., Chamayou, S., Esbert, M., et al. (2014). Proposed Guidelines on the Nomenclature and Annotation of Dynamic Human Embryo Monitoring by a Time-Lapse User Group. Hum. Reprod. 29, 2650–2660. doi:10.1093/humrep/deu278

PubMed Abstract | CrossRef Full Text | Google Scholar

Del Carmen Nogales, M., Bronet, F., Basile, N., Martínez, E. M., Liñán, A., Rodrigo, L., et al. (2017). Type of Chromosome Abnormality Affects Embryo Morphology Dynamics. Fertil. Sterility 107, 229–235. doi:10.1016/j.fertnstert.2016.09.019

CrossRef Full Text | Google Scholar

Deng, J., Kuyoro, O., Zhao, Q., Behr, B., and Lathi, R. B. (2020). Comparison of Aneuploidy Rates between Conventional In Vitro Fertilization and Intracytoplasmic Sperm Injection in In Vitro Fertilization-Intracytoplasmic Sperm Injection Split Insemination Cycles. F&S Rep. 1, 277–281. doi:10.1016/j.xfre.2020.07.006

CrossRef Full Text | Google Scholar

Desai, N., Goldberg, J. M., Austin, C., and Falcone, T. (2018). Are Cleavage Anomalies, Multinucleation, or Specific Cell Cycle Kinetics Observed with Time-Lapse Imaging Predictive of Embryo Developmental Capacity or Ploidy? Fertil. Sterility 109, 665–674. doi:10.1016/j.fertnstert.2017.12.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Dviri, M., Madjunkova, S., Koziarz, A., Madjunkov, M., Mashiach, J., Nekolaichuk, E., et al. (2021). Is There an Association between Paternal Age and Aneuploidy? Evidence from Young Donor Oocyte-Derived Embryos: a Systematic Review and Individual Patient Data Meta-Analysis. Hum. Reprod. Update 27, 486–500. doi:10.1093/humupd/dmaa052

PubMed Abstract | CrossRef Full Text | Google Scholar

Fragouli, E., Alfarawati, S., Spath, K., and Wells, D. (2014). Morphological and Cytogenetic Assessment of Cleavage and Blastocyst Stage Embryos. Mol. Hum. Reprod. 20, 117–126. doi:10.1093/molehr/gat073

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedenthal, J., Maxwell, S. M., Munné, S., Kramer, Y., Mcculloh, D. H., Mccaffrey, C., et al. (2018). Next Generation Sequencing for Preimplantation Genetic Screening Improves Pregnancy Outcomes Compared with Array Comparative Genomic Hybridization in Single Thawed Euploid Embryo Transfer Cycles. Fertil. Sterility 109, 627–632. doi:10.1016/j.fertnstert.2017.12.017

CrossRef Full Text | Google Scholar

Gardner, D. K., Lane, M., Stevens, J., Schlenker, T., and Schoolcraft, W. B. (2000). Blastocyst Score Affects Implantation and Pregnancy Outcome: towards a Single Blastocyst Transfer. Fertil. Sterility 73, 1155–1158. doi:10.1016/s0015-0282(00)00518-5

CrossRef Full Text | Google Scholar

Gianaroli, L., Magli, M. C., Cavallini, G., Crippa, A., Capoti, A., Resta, S., et al. (2010). Predicting Aneuploidy in Human Oocytes: Key Factors Which Affect the Meiotic Process. Hum. Reprod. 25, 2374–2386. doi:10.1093/humrep/deq123

PubMed Abstract | CrossRef Full Text | Google Scholar

Greaney, J., Wei, Z., and Homer, H. (2018). Regulation of Chromosome Segregation in Oocytes and the Cellular Basis for Female Meiotic Errors. Hum. Reprod. Update 24, 135–161. doi:10.1093/humupd/dmx035

PubMed Abstract | CrossRef Full Text | Google Scholar

Greco, E., Minasi, M. G., and Fiorentino, F. (2015). Healthy Babies after Intrauterine Transfer of Mosaic Aneuploid Blastocysts. N. Engl. J. Med. 373, 2089–2090. doi:10.1056/NEJMc1500421

PubMed Abstract | CrossRef Full Text | Google Scholar

Gruhn, J. R., Zielinska, A. P., Shukla, V., Blanshard, R., Capalbo, A., Cimadomo, D., et al. (2019). Chromosome Errors in Human Eggs Shape Natural Fertility over Reproductive Life Span. Science 365, 1466–1469. doi:10.1126/science.aav7321

PubMed Abstract | CrossRef Full Text | Google Scholar

Irani, M., Canon, C., Robles, A., Maddy, B., Gunnala, V., Qin, X., et al. (2020). No Effect of Ovarian Stimulation and Oocyte Yield on Euploidy and Live Birth Rates: an Analysis of 12 298 Trophectoderm Biopsies. Hum. Reprod. 35, 1082–1089. doi:10.1093/humrep/deaa028

PubMed Abstract | CrossRef Full Text | Google Scholar

Jones, K. T. (2008). Meiosis in Oocytes: Predisposition to Aneuploidy and its Increased Incidence with Age. Hum. Reprod. Update 14, 143–158. doi:10.1093/humupd/dmm043

PubMed Abstract | CrossRef Full Text | Google Scholar

Kane, S. C., Willats, E., Bezerra Maia e Holanda Moura, S., Hyett, J., and Da Silva Costa, F. (2016). Pre-Implantation Genetic Screening Techniques: Implications for Clinical Prenatal Diagnosis. Fetal Diagn. Ther. 40, 241–254. doi:10.1159/000449381

PubMed Abstract | CrossRef Full Text | Google Scholar

Kragh, M. F., Rimestad, J., Berntsen, J., and Karstoft, H. (2019). Automatic Grading of Human Blastocysts from Time-Lapse Imaging. Comput. Biol. Med. 115, 103494. doi:10.1016/j.compbiomed.2019.103494

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C.-I., Chen, C.-H., Huang, C.-C., Cheng, E.-H., Chen, H.-H., Ho, S.-T., et al. (2019). Embryo Morphokinetics Is Potentially Associated with Clinical Outcomes of Single-Embryo Transfers in Preimplantation Genetic Testing for Aneuploidy Cycles. Reprod. BioMedicine Online 39, 569–579. doi:10.1016/j.rbmo.2019.05.020

CrossRef Full Text | Google Scholar

Lee, C.-I., Chen, H.-H., Huang, C.-C., Chen, C.-H., Cheng, E.-H., Huang, J. Y., et al. (2020a). Effect of Interval between Human Chorionic Gonadotropin Priming and Ovum Pick-Up on the Euploid Probabilities of Blastocyst. Jcm 9, 1685. doi:10.3390/jcm9061685

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C.-I., Cheng, E.-H., Lee, M.-S., Lin, P.-Y., Chen, Y.-C., Chen, C.-H., et al. (2020b). Healthy Live Births from Transfer of Low-Mosaicism Embryos after Preimplantation Genetic Testing for Aneuploidy. J. Assist. Reprod. Genet. 37, 2305–2313. doi:10.1007/s10815-020-01876-6

CrossRef Full Text | Google Scholar

Lee, C.-I., Su, Y.-R., Chen, C.-H., Chang, T. A., Kuo, E. E.-S., Zheng, W.-L., et al. (2021). End-to-end Deep Learning for Recognition of Ploidy Status Using Time-Lapse Videos. J. Assist. Reprod. Genet. 38, 1655–1663. doi:10.1007/s10815-021-02228-8

CrossRef Full Text | Google Scholar

Lin, P.-Y., Lee, C.-I., Cheng, E.-H., Huang, C.-C., Lee, T.-H., Shih, H.-H., et al. (2020). Clinical Outcomes of Single Mosaic Embryo Transfer: High-Level or Low-Level Mosaic Embryo, Does it Matter? Jcm 9, 1695. doi:10.3390/jcm9061695

PubMed Abstract | CrossRef Full Text | Google Scholar

Martín, Á., Rodrigo, L., Beltrán, D., Meseguer, M., Rubio, C., Mercader, A., et al. (2021). The Morphokinetic Signature of Mosaic Embryos: Evidence in Support of Their Own Genetic Identity. Fertil. Sterility 116, 165–173. doi:10.1016/j.fertnstert.2020.12.031

CrossRef Full Text | Google Scholar

Mashiko, D., Ikeda, Z., Yao, T., Tokoro, M., Fukunaga, N., Asada, Y., et al. (2020). Chromosome Segregation Error during Early Cleavage in Mouse Pre-implantation Embryo Does Not Necessarily Cause Developmental Failure after Blastocyst Stage. Sci. Rep. 10, 854. doi:10.1038/s41598-020-57817-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mazzilli, R., Cimadomo, D., Vaiarelli, A., Capalbo, A., Dovere, L., Alviggi, E., et al. (2017). Effect of the Male Factor on the Clinical Outcome of Intracytoplasmic Sperm Injection Combined with Preimplantation Aneuploidy Testing: Observational Longitudinal Cohort Study of 1,219 Consecutive Cycles. Fertil. Sterility 108, 961–972. doi:10.1016/j.fertnstert.2017.08.033

CrossRef Full Text | Google Scholar

Minasi, M. G., Colasante, A., Riccio, T., Ruberti, A., Casciani, V., Scarselli, F., et al. (2016). Correlation between Aneuploidy, Standard Morphology Evaluation and Morphokinetic Development in 1730 Biopsied Blastocysts: a Consecutive Case Series Study. Hum. Reprod. 31, 2245–2254. doi:10.1093/humrep/dew183

PubMed Abstract | CrossRef Full Text | Google Scholar

Mumusoglu, S., Yarali, I., Bozdag, G., Ozdemir, P., Polat, M., Sokmensuer, L. K., et al. (2017). Time-lapse Morphokinetic Assessment Has Low to Moderate Ability to Predict Euploidy when Patient- and Ovarian Stimulation-Related Factors Are Taken into Account with the Use of Clustered Data Analysis. Fertil. Sterility 107, 413–421. doi:10.1016/j.fertnstert.2016.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Munné, S., Blazek, J., Large, M., Martinez-Ortiz, P. A., Nisson, H., Liu, E., et al. (2017). Detailed Investigation into the Cytogenetic Constitution and Pregnancy Outcome of Replacing Mosaic Blastocysts Detected with the Use of High-Resolution Next-Generation Sequencing. Fertil. Sterility 108, 62–71. doi:10.1016/j.fertnstert.2017.05.002

CrossRef Full Text | Google Scholar

Munné, S., and Wells, D. (2017). Detection of Mosaicism at Blastocyst Stage with the Use of High-Resolution Next-Generation Sequencing. Fertil. Sterility 107, 1085–1091. doi:10.1016/j.fertnstert.2017.03.024

CrossRef Full Text | Google Scholar

Pipari, A., Guillen, A., Cruz, M., Pacheco, A., and Garcia-Velasco, J. A. (2021). Serum Anti-müllerian Hormone Levels Are Not Associated with Aneuploidy Rates in Human Blastocysts. Reprod. BioMedicine Online 42, 1211–1218. doi:10.1016/j.rbmo.2021.03.006

CrossRef Full Text | Google Scholar

Rehman, K. S., Bukulmez, O., Langley, M., Carr, B. R., Nackley, A. C., Doody, K. M., et al. (2007). Late Stages of Embryo Progression Are a Much Better Predictor of Clinical Pregnancy Than Early Cleavage in Intracytoplasmic Sperm Injection and In Vitro Fertilization Cycles with Blastocyst-Stage Transfer. Fertil. Sterility 87, 1041–1052. doi:10.1016/j.fertnstert.2006.11.014

CrossRef Full Text | Google Scholar

Sachdev, N. M., Maxwell, S. M., Besser, A. G., and Grifo, J. A. (2017). Diagnosis and Clinical Management of Embryonic Mosaicism. Fertil. Sterility 107, 6–11. doi:10.1016/j.fertnstert.2016.10.006

CrossRef Full Text | Google Scholar

Schaeffer, E., Porchia, L. M., Neumann, A., Luna, A., Rojas, T., and López‐Bayghen, E. (2020). Embryos Derived from Donor or Patient Oocytes Are Not Different for In Vitro Fertilization Outcomes when PGT Allows Euploid Embryo Selection: a Retrospective Study. Clin. Translational Med. 9, 14. doi:10.1186/s40169-020-0266-1

CrossRef Full Text | Google Scholar

Stovezky, Y. R., Romanski, P. A., Bortoletto, P., and Spandorfer, S. D. (2021). Body Mass index Is Not Associated with Embryo Ploidy in Patients Undergoing In Vitro Fertilization with Preimplantation Genetic Testing. Fertil. Sterility 116, 388–395. doi:10.1016/j.fertnstert.2021.02.029

CrossRef Full Text | Google Scholar

Taylor, T. H., Gitlin, S. A., Patrick, J. L., Crain, J. L., Wilson, J. M., and Griffin, D. K. (2014). The Origin, Mechanisms, Incidence and Clinical Consequences of Chromosomal Mosaicism in Humans. Hum. Reprod. Update 20, 571–581. doi:10.1093/humupd/dmu016

PubMed Abstract | CrossRef Full Text | Google Scholar

Van Echten-Arends, J., Mastenbroek, S., Sikkema-Raddatz, B., Korevaar, J. C., Heineman, M. J., Van Der Veen, F., et al. (2011). Chromosomal Mosaicism in Human Preimplantation Embryos: a Systematic Review. Hum. Reprod. Update 17, 620–627. doi:10.1093/humupd/dmr014

PubMed Abstract | CrossRef Full Text | Google Scholar

Viotti, M., Victor, A. R., Barnes, F. L., Zouves, C. G., Besser, A. G., Grifo, J. A., et al. (2021). Using Outcome Data from One Thousand Mosaic Embryo Transfers to Formulate an Embryo Ranking System for Clinical Use. Fertil. Sterility 115, 1212–1224. doi:10.1016/j.fertnstert.2020.11.041

CrossRef Full Text | Google Scholar

Wang, A., Kort, J., Behr, B., and Westphal, L. M. (2018). Euploidy in Relation to Blastocyst Sex and Morphology. J. Assist. Reprod. Genet. 35, 1565–1572. doi:10.1007/s10815-018-1262-x

CrossRef Full Text | Google Scholar

Wirleitner, B., Schuff, M., Stecher, A., Murtinger, M., and Vanderzwalmen, P. (2016). Pregnancy and Birth Outcomes Following Fresh or Vitrified Embryo Transfer According to Blastocyst Morphology and Expansion Stage, and Culturing Strategy for Delayed

Google Scholar

Zhan, Q., Sierra, E. T., Malmsten, J., Ye, Z., Rosenwaks, Z., and Zaninovic, N. (2020). Blastocyst Score, a Blastocyst Quality Ranking Tool, Is a Predictor of Blastocyst Ploidy and Implantation Potential. F&S Rep. 1, 133–141. doi:10.1016/j.xfre.2020.05.004

CrossRef Full Text | Google Scholar

Keywords: blastocyst morphology, mosaic levels, preimplantation genetic testing for aneuploidy, time-lapse monitoring, high-resolution next-generation sequencing

Citation: Chen C-H, Lee C-I, Huang C-C, Chen H-H, Ho S-T, Cheng E-H, Lin P, Chen C-I, Lee T-H and Lee M-S (2021) Blastocyst Morphology Based on Uniform Time-Point Assessments is Correlated With Mosaic Levels in Embryos. Front. Genet. 12:783826. doi: 10.3389/fgene.2021.783826

Received: 27 September 2021; Accepted: 06 December 2021;
Published: 22 December 2021.

Edited by:

Yan Yun, University of California, Davis, United States

Reviewed by:

Silvestre Sampino, Polish Academy of Sciences, Poland
Alison Campbell, CARE Fertility Nottingham, United Kingdom

Copyright © 2021 Chen, Lee, Huang, Chen, Ho, Cheng, Lin, Chen, Lee and Lee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Tsung-Hsien Lee, jackth.lee@gmail.com; Maw-Sheng Lee, msleephd@gmail.com

These authors have contributed equally to this work

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.