Reassessing the feasibility of the zygote score for predicting embryo viability in IVF/ICSI using the GnRH antagonist protocol compared to the long protocol.

Clinicians have used the GnRH-ant protocol for assisted reproduction in recent decades as an alternative to the GnRH agonist (GnRH-ag) protocol [ 1 – 3 ]. Recent researches have shown that the use of GnRH antagonist is associated with decreasing the risk of severe ovarian hyperstimulation syndrome without reducing live birth rate. [ 4 – 7 ]. However, the impact of ovarian stimulation on oocyte and embryo quality is still unclear.[ 8 ] The zygote score (Z-score) is an established scoring system based on the pronuclear morphology, nucleolar precursor bodies, and alignment of polar bodies in the human zygote at approximately 16 to 18 h after insemination [ 9 ]. Scoring of zygotes has proven to be useful for selection at the zygote stage for embryo transfer or cryopreservation, especially in countries with legal restrictions [ 10 , 11 ].Although numerous studies have associated positive clinical results with the use of the Z-score for embryo selection [ 9 , 12 – 14 ], other studies did not find the Z-score to be a valuable predictor [ 15 – 17 ]. When reviewing the literature describing the use of the Z-score in morphologic selection [ 18 ], most studies analyzed the Z-score in the GnRH-ag protocol, without considering the feasibility of using the Z-score for the GnRH-ant protocol. The importance of zygote scoring with the gonadotropin-releasing hormone antagonist (GnRH-ant) protocol has not been established. The present retrospective study examined a cohort of embryos derived from patients undergoing in vitro fertilization/intra-cytoplasmic sperm injection (IVF/ICSI). The aim of this study was to determine the relationship between the different controlled ovarian stimulation (COS) protocols and the Z- score and to determine the feasibility of using the Z-score to predict embryo survival in the GnRH-ant protocol.

During the period, 744 fresh cycles from 643 women were included, which consisted of 352 long protocol cycles and 392 GnRH-ant protocol cycles. Table 1 shows the characteristics of these cycles. A total of 3,826 zygotes with normal fertilization produced from 744 fresh cycles following IVF/ICSI were examined. We performed fertilization check using the Z-score, a system that has been used to evaluate fertilization in our center since 2001. Fig 1 shows the distribution of Z-scores. Most zygotes in the long protocol had scores of Z1 and Z2. (N/ (%); Z1: 765(40.1%), Z2: 777(40.8%), Z3:316(16.6%), Z4:48(2.5%)) Fewer zygotes in the GnRH-ant protocol had scores of Z1, and more had scores of Z2. (N/ (%); Z1: 567(29.5%), Z2:874(45.5%), Z3:400(20.8%), Z4:79(4.1%)). The Z-score distributions in the two groups were significantly different ( p < 0.05) and shifted toward poorer Z-scores in the GnRH-ant protocol. In this retrospective and observational study, we did not perform statistical comparisons between cycles from the different protocols due to the risk of selection bias. Instead, we observed the cohort zygotes from the two protocols and analyzed potential factors that may influence the zygote score. Our previous study indicated good (Z-1) zygote had higher implantation potential.[ 12 ]We further analyzed the clinical and laboratory factors associated with a good Z-score of Z1. Univariate analyses ( Table 2 ) showed that a score of Z1 correlated significantly with the COS protocol used ( p < 0.001), oocyte maturity ( p < 0.001), progesterone on the day of hCG detection ( p < 0.001), E2 concentration per oocyte ( p = 0.028), total FSH dosage ( p = 0.025), and duration of FSH stimulation ( p = 0.028). a A p- value less than 0.05 was considered statistically significant We then used multivariable logistic regression analysis to identify factors associated with a score of Z1 ( Table 3 ). These results indicate that a score of Z1 was significantly and independently associated with oocyte maturity ( p < 0.001), number of cohort 2PN ( p < 0.001), and COS protocol used ( p < 0.001). The observation of significantly fewer Z1 zygotes from the GnRH-ant protocol is compatible with our clinical observations. B: intercept; SEM: standard error of the mean; Wald: Wald statistic; Exp (B): adjusted odds ratio; CI: confidence interval. * A p- value less than 0.05 was considered statistically significant The relationships between the Z-score and good day-3 embryo and blastocyst formation were investigated.( Fig 2A and 2B ) The results indicated that the Z-score had similar correlations with good embryo formation on day-3 for both protocols. Long protocol vs. GnRH-ant: (N/ (%); Z1: 539(70.6%) vs . 370(65.3%), Z2: 174(22.4%) vs . 208(23.7%), Z3: 13(4.1%) vs . 21(5.2%), Z4: 0(0%) vs . 1(1%); p = 0.696) ( Fig 2A ). An analysis of the relationship between Z-score and blastocyst formation also indicated similar correlations for the two protocols. Long protocol vs. GnRH-ant (N/ (%); Z1:288(38%) vs . 198(35%), Z2:123(16%) vs . 194(22%), Z3: 9(3%) vs . 45(11%), Z4:0(0%) vs . 2(3%); p = 0.066) ( Fig 2B ). Furthermore, we tested the feasibility of using the Z-score in the GnRH-ant protocol to predict embryo survival on day 5 ( Fig 3A and 3B ). The overall survival rates for the long protocol and GnRH-ant protocol were similar (59 ± 1% vs . 62 ± 1%, p = 0.370) ( Fig 3A ). As expected, the overall survival rate on day 5 declined as the Z-score increased (Z1: 89 ± 1%, Z2: 60 ± 1%, Z3: 24 ± 2%, Z4:29 ± 4%; p <0.05; Fig 3B ). Relationship between the Z- score with the formation of good quality embryos on day 3 (A) and blastocysts on day 5 or 6 (B) following the long protocol and the GnRH-ant protocol. Overall day-5 arrest‐free survival of zygotes resulting from the long protocol (solid line) and the GnRH-ant protocol (dotted line) (A) and of zygotes with different Z-scores (B). Cumulative survival rates were compared for the long protocol and GnRH-ant protocol for zygotes with different Z-scores ( Fig 4A–4D ). The GnRH-ant protocol led to better survival for Z1 zygotes (85 ± 1% vs . 94 ± 1%, p < 0.001)( Fig 4A ), Z2 zygotes (54 ± 2% vs . 65 ± 2%, p < 0.001) ( Fig 4B ), and Z4 zygotes (14 ± 5% vs . 35 ± 4%, p = 0.022) ( Fig 4D ), but not for Z3 zygotes (22 ± 3% vs . 25 ± 2%, p = 0.615) ( Fig 4C ). Day-5 arrest‐free survival of zygotes resulting from the long protocol (solid line) and the GnRH-ant protocol (dotted line) with scores of Z1 (A), Z2 (B), Z3 (C), and Z4 (D). Thus, our study found that good-quality day 3 embryo/blastocyst formation and the cumulative survival rate were correlated with the Z-score but not the COS protocol. Fig 5A and 5B displays the relationships between the number of good zygotes (Z1) transferred and pregnancy outcomes in the long (A) and GnRH-ant (B) protocols. In our study, a mean of 2.19 ± 0.86 embryos were transferred for each woman; therefore, it is difficult to compare the pregnancy outcomes according to Z-score based on single embryo transfer. Instead of single embryo transfer, we compared the number of good zygotes (Z1) transferred and pregnancy outcomes. For the long protocol, the number of Z1 in the transferred cohort embryos was not significantly correlated with the clinical pregnancy rate ( r = 0.976; p = 0.088) or live birth rate ( r = 0.944; p = 0.056). For the GnRH-ant protocol, the number of Z1 in the transferred cohort embryos was significantly correlated with the clinical pregnancy rate ( r = 0.976; p = 0.024) and live birth rate ( r = 0.971; p = 0.029). The correlation between the number of transferred good zygote (Z1) and pregnancy outcomes in long protocol (A) and GnRH-ant (B).

To our knowledge, the current study is the first to report the effect of different COS protocols on zygote morphology and the feasibility of using the Z-score to predict embryo survival in the GnRH-ant protocol. Zygote formation follows the dramatic reorganization of sperm chromatin by many factors stored within the oocyte [ 25 , 26 ]. The morphological characteristics of the zygote indicate gamete quality and the potential for subsequent embryo implantation [ 27 , 28 ]. Zygotes with unequal numbers or sizes of nucleoli likely display asynchrony between male and female pronucleus development [ 29 , 30 ]. Many factors from the oocyte/sperm or the process of fertilization may affect the Z-score. Our IVF laboratory adopted a modified Scott scoring system beginning in 2001 to monitor fertilization. Our previous report (Lan, 2003) indicated that the Z‐score is an additional criterion that can be used to select embryos for extended culture and is useful when selecting embryos for transfer. However, the distribution of Z‐scores changes significantly when the GnRH-ant protocol is used for COS. In particular, we found that the Z-scores from the GnRH-ant protocol were generally poorer than those from the long protocol. Significantly fewer Z1 zygotes resulted from the GnRH-ant protocol, which is compatible with our clinical observations. Some possible explanation for the differences in the Z-score between the COS protocols involves the different endocrine profiles during follicular growth. First, the long protocol induces profound suppression and leads to simultaneous maturation of antral follicles [ 8 , 31 , 32 ]. Therefore, more mature follicles are recruited, and typically more Z1 zygotes are obtained. In contrast, the GnRH-ant protocol generates a more natural pattern of follicular recruitment, with uneven follicle sizes and fewer mature oocytes. Second, previous studies have suggested that the COS protocols have differing effects on ovarian E2 metabolism [ 33 , 34 ] and may affect follicular growth and/or luteal function. More specifically, Khalaf et al . [ 32 ] found that the protein kinase C (PKC) pathway was desensitized in GnRH-ag-treated granulosa cells, and this increased the levels of FSH and cAMP-mediated steroidogenesis. Third, in addition to the different endocrine profiles between the two protocols, the effect of the GnRH receptor on the oocyte/zygote should also be considered. GnRH receptors are expressed in the follicles at the gonadotropin-sensitive stage and in luteal cells. The correlation between the expression of GnRH receptors and follicular stage suggests that GnRH receptors directly influence folliculogenesis and oocyte development [ 35 , 36 ], but the mechanism by which GnRH affects the ovary is not completely understood. Thus, when GnRH-ant is used in IVF cycles, it has an unclear effect on the morphology and quality of oocytes. Animal studies indicate that GnRH-ant may have an adverse effect on primordial follicle survival in some species [ 37 ], but only a few studies examined the effect of GnRH-ant on human oocytes and embryos, with inconsistent results [ 38 – 40 ]. The results of our study also indirectly support the hypothesis that GnRH receptors play a role in oocyte maturation and affect zygote formation. Although more mature oocytes/Z-1 zygotes were produced from long protocol than GnRH-ant possibly because of better endocrine profiles and no effect of GnRH antagonist; however, the overall Z-1 zygotes from long protocol had poorer survival and not correlated with clinical pregnancy. Our results suggested that long protocol only improved morphological maturity of oocytes and zygotes’ morphology without changing the survival potential of oocyte. Some Z1 zygotes from long protocol may originated from immature oocytes with poorer potential and therefore had poorer pregnancy outcomes than Z1 zygotes from GnRH-ant protocol. The findings may also explain why Z-score could not be consistently verified in previous studies which used long protocol. There has been debate regarding the use of the Z-score for predicting outcomes following assisted reproduction [ 9 , 25 , 41 – 49 ]. A recent systematic review [ 18 ] reported no conclusive data on the clinical efficacy of the Z-score in fresh cycles, even though biological results showed a good relationship with embryo viability and suggest a role in cycles with day-1 transfer/freezing. These inconsistent results may have occurred because almost all of these studies used the long protocol rather than GnRH-ant. With GnRH-ant, the number of Z1 in the transferred cohort of embryos was significantly correlated with the clinical pregnancy rate and live birth rate in our study ( r = 0.971; p = 0.029). However, this correlation was not observed for the long protocol. A possible explanation for the Z-score having better predictive value in the GnRH-ant protocol is that zygotes from the GnRH-ant protocol were taken from a more natural pattern of follicle recruitment, while the long protocol only improved the zygotes’ morphology without changing the survival potential. The present study has some limitations. First, this retrospective study has risk of selection bias. A matched control group might improve the quality of the findings. Second, the relationship between pregnancy outcomes and Z-score for single embryo transfer requires further investigation. A mean of approximately two embryos were transferred to most patients in this study, so single-embryo outcomes could not be compared with the Z-score. Third, the timing of assessment is critical, as pronuclear development is a dynamic process. Thus, determination of the Z-score from a single light microscopy observation should be used with caution and only in conjunction with other methods of evaluation [ 27 ]. In conclusion, the Z-score distribution was poorer for the GnRH-ant protocol than the long protocol. This may be because these protocols have different effects on the endocrine profile or GnRH receptors. Additionally, the Z-score is a more feasible parameter for predicting embryo viability in IVF/ICSI with the GnRH-ant protocol than the long protocol.

This retrospective study included 3,826 zygotes with normal fertilization from 643 women and followed 744 IVF/ICSI cycles at the Chang-Gung Memorial Hospital Kaohsiung Medical Center between January 1, 2010 and April 30, 2014. All of the women had cycles from COS with a luteal phase down-regulated protocol (long protocol) or a GnRH-ant protocol, and all couples had at least one zygote with two pronuclear (PN) formations after insemination/ICSI at 16 to 18 h. All cases with unsuccessful oocyte retrieval or polyspermic zygotes were excluded. Women were excluded if they underwent COS with natural cycles or clomiphene using an ultra-long protocol or a short protocol. All charts were retrospectively reviewed by one physician. The study was approved by an appropriately constituted the Institutional Review Board of the Ethics Committee of Chang Gung Memorial Hospital (Institutional Review Board Number: 201600478B0) and Institutional Review Board waived the need for consent. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. COS, oocyte retrieval, embryo culture, and embryo transfer were performed according to our previously described protocols [ 12 , 19 ]. The long protocol and GnRH-ant protocol for COS were individualized and depended on ovarian reserve, age, baseline serum follicle stimulating hormone (FSH) concentration, and previous response to COS. Briefly, women received the long protocol with pituitary down-regulation using leuprolide acetate (Lupron®; Takeda, Tokyo, Japan) with the initial dose of gonadotropin, either human menopausal gonadotropin (hMG) or FSH (purified or recombinant), individualized for each patient (range: 150 to 300 IU). Further dose adjustments were made according to the individual’s ovarian response, based on the serum estradiol (E2) level and ultrasonographic monitoring of follicular growth. When the lead follicle was 16–18 mm in diameter, leuprolide acetate and FSH were discontinued, and human chorionic gonadotropin (hCG) (Ovidrel®; Serono, Modugno, Italy) was administered. Oocyte retrieval was performed by transvaginal ultrasound-guided follicle aspiration at 36–38 h after hCG administration. In the GnRH-ant protocol, patients were given gonadotropin stimulation and then suppression using a flexible GnRH-ant protocol (Ganirelix acetate: 0.25 mg, MSD; Cetrorelix acetate: 0.25 mg; Serono) when the leading follicle was 14 mm. When 2 more follicles reached diameters of 17 mm, a 6,500 IU dose of hCG (Ovidrel®; Serono, Modugno, Italy) was administered. Oocyte retrieval was performed 36–38 h later by transvaginal aspiration under ultrasound guidance. Oocytes were graded for maturity on the basis of the morphological characteristics of the cumulus mass, corona radiata, ooplasm, and detached membrane granulose cells. [ 20 ] Standard IVF/ICSI procedures were used for oocyte fertilization, as previously described [ 12 ]. Fertilization was evaluated 16–18 h after IVF or ICSI. A zygote with two pronuclei (2PN) was defined as normal fertilization. All zygotes were scored according to the modified Scott scoring system [ 21 ]. A single team of embryologists coordinated all procedures of IVF laboratory to ensure that the culture processes and the embryo assessments were standardized. For PN scoring, zygotes were divided into 4 categories (Z1-Z4) based on 3 major features: size and location of the nuclei, appearance of the cytoplasm, and numbers, sizes, and distribution patterns of nucleolar precursor bodies within the nuclei. ( S1 Fig )) Z1 zygotes have an equal number of nucleoli aligned at the PN junctions. Z2 zygotes have an equal number and size of nucleoli (3 to 7) that are equally scattered in the 2 PNs, but the nucleoli have not yet aligned at the PN junction. Z3 zygotes are characterized by inequality of the nuclei (unequal size, unequal numbers, or unequal alignment at the PN junction). Z4 zygotes have PN that are separated or different in size and small nucleoli that are partially aligned or scattered[ 21 ]. Embryos were cultured on days 1 to 3 in G1 TM medium (Scandinavian IVF Science) and on days 3 to 5 in G2 TM medium (Scandinavian IVF Science). Veeck’s morphological grading system [ 22 ] was adopted for day-3 embryo scoring. A “good” 3-day embryo was defined as one that had a Veeck’s grade of 1, with 8 cells, blastomeres of equal size, and no cytoplasmic fragments. Embryos were assessed by survival, morphology, and rate of cleavage. Embryos that had the same number of blastomeres at two sequential observations and zygotes that remained blocked at the pronuclear stage were classified as “developmentally arrested”. Embryos were transferred on day 2, 3, 4, 5 or 6 after oocyte retrieval, as appropriate for each individual. The risk of embryo arrest was defined as the time to a first event. The start point for determining the duration of embryo survival was the date when PN appeared after 16–18 h of incubation, and the end point was after 5 days of extended culture. The study endpoints were the duration of overall survival and event-free survival ( i . e . embryo arrest). In the analysis of embryo arrest, data from patients who received a day-2, day-3, day-4, day-5 or day-6 embryo transfer were censored after the time of transfer (loss due to transfer) or cryopreservation. Luteal phase support continued until the day pregnancy was confirmed by detection of hCG in the urine. If conception occurred, micronized progesterone supplementation was provided for an additional 4 weeks. Clinical pregnancy was defined as the presence of an intrauterine gestational sac with positive cardiac movement on ultrasound [ 23 ]. The live birth rate per transfer was defined as the proportion of IVF cycles reaching embryo transfer that resulted in the birth of at least one live-born child. Continuous data are given as the mean ± standard deviation (SD). The Mann-Whitney rank sum test was used to compare continuous data, and the χ 2 test or Fisher’s exact test was used to compare categorical variables. A Pearson correlation coefficient ( r ) was calculated for correlations. All statistical analyses were performed with SPSS, version 17.0 (SPSS, Inc., Chicago, IL, USA). All p values were two-sided, and a p value less than 0.05 was considered statistically significant. Logistic regression analysis was used to determine the significance of factors associated with a Z-score of 1 (Z1). Multivariate logistic regression analysis was then performed to further examine factors associated with Z1, considering the following variables: female age, male age, ICSI or IVF, etiology of infertility, body mass index, COS protocol, days of stimulation, FSH dosage, E2 level on the day of hCG detection, E2 level per oocyte, P4 level, maturity of oocyte, and number of same-cohort oocytes. Cumulative survival rates were calculated by life-table analysis using the Kaplan–Meier product limit procedure at each day. The differences between groups with different Z-scores were assessed using the log-rank test. A p -value less than 0.05 was considered significant [ 24 ].

Zygote scoring system of Scott et al.[ 21 ] Z1 zygotes have an equal number of nucleoli aligned at the PN junctions (A). Z2 zygotes have an equal number and size of nucleoli (3 to 7) that are equally scattered in the 2 PNs(B) Z3 zygotes are characterized by inequality of the nuclei (unequal size, unequal numbers, or unequal alignment at the PN junction)(C and D). Z4 zygotes have PN that are separated or different in size and small nucleoli that are partially aligned or scattered.(E and F) (TIF) Click here for additional data file.