- Open Access
- Open Peer Review
Clinical outcome of preimplantation genetic diagnosis and screening using next generation sequencing
GigaSciencevolume 3, Article number: 30 (2014)
Next generation sequencing (NGS) is now being used for detecting chromosomal abnormalities in blastocyst trophectoderm (TE) cells from in vitro fertilized embryos. However, few data are available regarding the clinical outcome, which provides vital reference for further application of the methodology. Here, we present a clinical evaluation of NGS-based preimplantation genetic diagnosis/screening (PGD/PGS) compared with single nucleotide polymorphism (SNP) array-based PGD/PGS as a control.
A total of 395 couples participated. They were carriers of either translocation or inversion mutations, or were patients with recurrent miscarriage and/or advanced maternal age. A total of 1,512 blastocysts were biopsied on D5 after fertilization, with 1,058 blastocysts set aside for SNP array testing and 454 blastocysts for NGS testing. In the NGS cycles group, the implantation, clinical pregnancy and miscarriage rates were 52.6% (60/114), 61.3% (49/80) and 14.3% (7/49), respectively. In the SNP array cycles group, the implantation, clinical pregnancy and miscarriage rates were 47.6% (139/292), 56.7% (115/203) and 14.8% (17/115), respectively. The outcome measures of both the NGS and SNP array cycles were the same with insignificant differences. There were 150 blastocysts that underwent both NGS and SNP array analysis, of which seven blastocysts were found with inconsistent signals. All other signals obtained from NGS analysis were confirmed to be accurate by validation with qPCR. The relative copy number of mitochondrial DNA (mtDNA) for each blastocyst that underwent NGS testing was evaluated, and a significant difference was found between the copy number of mtDNA for the euploid and the chromosomally abnormal blastocysts. So far, out of 42 ongoing pregnancies, 24 babies were born in NGS cycles; all of these babies are healthy and free of any developmental problems.
This study provides the first evaluation of the clinical outcomes of NGS-based pre-implantation genetic diagnosis/screening, and shows the reliability of this method in a clinical and array-based laboratory setting. NGS provides an accurate approach to detect embryonic imbalanced segmental rearrangements, to avoid the potential risks of false signals from SNP array in this study.
Chromosomal abnormalities, including numerical errors and structural anomalies, are widespread in human embryos produced in vitro, and the incidence increases dramatically in embryos with advancing maternal age [2–4]. Chromosomal abnormalities are a main reason for spontaneous abortions [5, 6], and repeated implantation failures after transfer of in vitro produced embryos . Furthermore, embryos from carriers of balanced translocations [8, 9] are at particularly high risk of chromosomal abnormalities. Previous papers have been published about the application of preimplantation genetic screening (PGS) in cases with advanced maternal age [10, 11] or recurrent pregnancy loss , and the application of preimplantation genetic diagnosis (PGD) in carriers of translocations . Thus, PGD/PGS and the selection of chromosomally normal embryos for transfer should be an effective approach to improve live birth rate, as well as reduce spontaneous abortion and birth defects. Studies were performed using blastocyst biopsy  and either fluorescence in situ hybridization (FISH) , comparative genomic hybridization (CGH) [16–18] or array based methods, such as array-CGH [13, 19, 20] and single-nucleotide polymorphism array (SNP array) [21–24] for sample analysis. So far, FISH was the most common approach for sample analysis, although the obtainable information was limited: the FISH method could only detect a few chromosomes, and no detailed information could be obtained about the sub-chromosomal abnormalities. Array based techniques (aCGH and SNP array) have been used to screen all 24 human chromosomes. However, array-based methods are still relatively expensive, restricting their clinical applications.
Next generation sequencing (NGS) has now become cost-effective, and this precise and comprehensive genetic analysis tool is being increasingly used in human medicine including noninvasive prenatal diagnosis , and there is increasing interest for its application in PGD/PGS as well. Previous studies [26, 27] performed on trophectodermal biopsies have proven that aneuploid and unbalanced rearrangements can be detected accurately by NGS. However, no data are available to confirm the efficiency and clinical outcome of NGS-based PGD/PGS.
Here, we present the first clinical evaluation for NGS-based PGD/PGS combined with trophectoderm (TE) biopsy and cryopreserved embryo transfer (CET). We also employed SNP array as control for comparison in this study. The clinical outcomes including implantation, clinical pregnancy and miscarriage were recorded. In addition, the accuracy of NGS and SNP array were also evaluated in this study.
A total of 395 couples were subjected to in vitro fertilization-preimplantation genetic diagnosis (IVF-PGD) treatment, including 129 couples with a NGS-based test and 266 couples with a SNP array-based test, for the detection of embryonic chromosomal abnormalities. The NGS test was performed using low coverage whole-genome sequencing with a HiSeq 2000 platform. The SNP array test was performed using Affymetrix Gene Chip Mapping Nsp I 262 K. The average age of the patients was 32.1 years (with an age range of 20–44 years). In the NGS cycles group, 84 patients were confirmed to have a chromosomally abnormal karyotype; 18 of which had a Robertsonian translocation, 59 had a reciprocal translocation, and 7 patients had a inversion. Another 45 couples were included because of an advanced maternal age (AMA; ≥38 years) and/or recurrent miscarriage (RM; ≥2). In the SNP array cycles group, 213 patients (of which 58 had a Robertsonian translocation, 144 had a reciprocal translocation, 11 had an inversion) and another 57 couples with AMA and/or RM were included. A total of 1,512 blastocysts were obtained, with 454 blastocysts for NGS testing and 1,058 blastocysts for SNP array analysis. In addition, 150 blastocysts from the 454 blastocysts were subjected to both NGS and SNP array tests, and blastocysts with inconsistent results would be further validated by qPCR.
NGS and SNP array testing of blastocysts
Of the 454 blastocysts with NGS testing, an average of 8.2 million reads was obtained for each blastocyst, covering 5.5% ± 1.2% of the whole human genome and 98.7% ± 3.1% of the mitochondrial DNA (mtDNA). The ratio between the mean depth of the mitochondrial DNA and the genome for each blastocyst with NGS testing was calculated, which represented the relative copy number of mtDNA. It is displayed in Additional file 1: Table S3 as the Ratio (ChrM_depth/Mean_depth). The box plot of the data is shown in Additional file 2: Figure S2.
A total of 198 (43.6%) euploid and 256 (56.4%) chromosomally abnormal blastocysts were identified by NGS. Among the abnormalities, 85 blastocysts had numerical chromosome aberrations, 127 contained imbalanced structural aberrations, and 44 blastocysts were detected with both numerical and imbalanced structural aberrations. Among the diagnosed blastocysts, the median number of normal/balanced embryos per couple was 1.53 (range from 0 to 6). The group of reciprocal translocation carriers obtained 43.0% (105/244) embryos with imbalanced aberrations, most of the imbalances were translocation related. Normal/balanced embryos constituted 54.6% (65/119) of the cohort in couples with AMA and/or RM (Table 1).
Of the 1,058 blastocysts with SNP array analysis, 468 (44.2%) euploid and 590 (55.8%) chromosomally abnormal blastocysts were identified, 189 blastocysts had numerical chromosome aberrations, 298 contained imbalanced structural aberrations and 103 blastocysts were detected with both numerical and imbalanced structural aberrations. The median number of normal/balanced embryos per couple was 1.76 (range from 0 to 8). The group of reciprocal translocation carriers obtained 40.8% (247/606) embryos with imbalanced aberrations. Normal/balanced embryos constituted 60.2% (106/176) of the cohort in couples with AMA and/or RM (Table 2).
Among the 129 couples in the NGS cycles group, 33 couples had no euploid embryos suitable for transfer; 75 couples underwent embryo transfer and the remaining 21 couples are currently still waiting for transfer. In the SNP array cycles group, 177 couples underwent embryo transfer, 66 couples had no suitable embryos for transfer, and 23 couples are currently still waiting.
Of the 666 normal/balanced blastocysts, 421 blastocysts were warmed after vitrification, 406 survived (96.4% of survival rate) and were transferred in 283 cycles. The numbers of blastocysts transferred per cycle were 1.425 (114/80) and 1.438 (292/203) for NGS and SNP array, respectively. The proportion of transferred embryos that successfully implanted was evaluated by ultrasound 6–7 weeks after embryo transfer, indicating that 60 and 139 embryos resulted in a fetal sac, giving implantation rates of 52.6% (60/114) and 47.6% (139/292) for NGS and SNP array, respectively. Prenatal diagnosis with karyotyping of amniocentesis fluid samples did not find any fetus with chromosomal abnormalities. A total of 164 pregnancies were detected, with 129 singletons and 35 twins. The clinical pregnancy rate per transfer cycle was 61.3% (49/80) and 56.7% (115/203) for NGS and SNP array, respectively (Table 3). A total of 24 miscarriages were detected, giving rates of 14.3% (7/49) and 14.8% (17/115) in NGS and SNP array cycles, respectively. In NGS cycles group, four miscarriages were spontaneous abortions occurred in early pregnancy, and the other three were artificial abortions with two caused by embryo diapause and another one caused by extrauterine (cervical) implantation. In the SNP array cycles group, 12 spontaneous abortions in early pregnancy and five artificial abortions occurred. Testing of the miscarried tissue by comparative genomic hybridization (CGH) revealed that the miscarried embryos in both NGS and SNP array cycles were chromosomally normal. The ongoing pregnancy rates were 52.5% (42/80) and 48.3% (98/203) in NGS and SNP array cycles, respectively. Out of these pregnancies, 24 babies were delivered in 20 NGS cycles; so far, all the babies are healthy and chromosomally normal according to karyotype analysis. In the SNP array cycles group the outcome of all pregnancies went to full term and 75 healthy babies were delivered (Table 3).
The outcome measures of NGS and SNP array cycles were used for Pearson’s chi-squared test. The p-value of clinical pregnancy, ongoing pregnancy, miscarriage and implantation rates between NGS and SNP array groups were 0.480, 0.522, 0.934 and 0.362, respectively (Table 3). All of these p-values far outweighed 0.05 that does not reach a specified significance level. It can be concluded that the four outcome measures between NGS and SNP array cycles were the same with insignificant difference.
The Ratio (ChrM_depth/Mean_depth) of the euploid and the chromosomally abnormal blastocysts with NGS testing in Additional file 1: Table S3 was used for Mann–Whitney test. The p-value was 0.017, which was less than 0.05. It can be concluded that the copy number of mtDNA for the chromosomally abnormal blastocysts were significantly higher than the cohort for the euploid blastocysts by NGS testing.
Cases with both NGS and SNP array analysis
Among the NGS-tested 460 blastocysts, 150 blastocysts were selected randomly for Affymetrix SNP array analysis as well, including 68 normal/balanced blastocysts and 82 blastocysts with chromosomal abnormalities. All of the normal/balanced blastocysts in the NGS test group were confirmed by SNP array. Among the chromosomally abnormal blastocysts, 7 blastocysts gave inconsistent results, which were further validated by qPCR. The detailed qPCR data for the 7 blastocysts can be seen in Additional file 1: Table S1. The qPCR data proved the accuracy of NGS since all the qPCR signals were in accordance with NGS (Table 4; Additional file 1: Table S1). There were 4 embryos suspected to be with false-negative signals and 2 embryos suspected to be with false-positive signals in the SNP array results, which were all detected accurately by NGS. A case with gain of Y chromosome was detected by NGS, but omitted by array for the lack of a Y chromosome marker in the SNP array.
In this study, we provide the first evaluation of clinical outcome of NGS-based PGD/PGS. According to the NGS results of 454 blastocysts, we transferred 114 normal/balanced embryos into the uterus, with implantation, clinical pregnancy and ongoing pregnancy rates of 52.6% (60/114), 61.3% (49/80) and 52.5% (42/80), respectively. Twenty-four babies were born, so far, all of them healthy, indicating the accuracy of NGS as well as the safety of blastocyst biopsy and vitrification. It illustrates that NGS-PGD is applicable for the genetically high-risk populations, such as carriers of Robertsonian translocations and reciprocal translocations. The feasibility of NGS-PGS application for the patients with AMA or RM was evaluated as well, though the benefits of aneuploidy screening for this population remain hypothetical [28, 29].
Array based methods (including aCGH and SNP arrays) provide an approach to screen all the 24 human chromosomes simultaneously, and are well established when applied to PGD/PGS – SNP array is widely used in clinics. Furthermore, an initial study  indicated that the pregnancy rate for translocation carriers can be increased significantly by SNP array-PGD compared with FISH-PGD.
NGS is a rapidly developing method for whole genome testing. Our previous study  indicated that aneuploidy and imbalanced rearrangements can be detected accurately by NGS in TE cells at a much lower depth (0.07X) of sequencing. There were 150 blastocysts with both NGS and SNP array testing in this study, 7 blastocysts were detected with inconsistent results (Table 4), and all the NGS results were confirmed to be accurate by qPCR validation. NGS provides high accuracy in the detection of some imbalanced segmental rearrangements due to the ability for NGS to correct the bias from WGA .
The outcome measures of NGS cycles were the same as SNP array, that the differences of the outcome measures were insignificant between these two methods. The numbers of normal/balanced embryos per couple and transferred blastocysts per cycle for the two methods were also insignificantly different. From the 150 blastocysts with both NGS and SNP array testing, we can find that the same normal/balanced embryos were detected by NGS and SNP array. Thus, the efficiency and accuracy for embryo selection would be similar between NGS and SNP array cycles, as well as the clinical outcome measures.
There were 6 embryos with signals to be inconsistent with NGS and qPCR results in SNP array data. These signals may due to the fact that only a maternal or a paternal allele is present in SNP array results and so forth, which can easily lead to false-negatives and false-positives. A case of a Y chromosome gain was detected by NGS, but omitted by array for the lack of a Y chromosome marker in SNP array. False-negatives in PGD/PGS clinical practice may lead to a higher rate of spontaneous abortion and/or birth defect . False-positives in the clinical application may cause a lower rate of clinical pregnancy or implantation [31, 32]. The differences of outcome measures between NGS and SNP array cycles were insignificant, which may be due to the fact that most inexact signals of SNP array did not impact the efficiency of embryo selection in this study. Nevertheless, NGS was able to detect some segmental imbalances more precisely. So that the potential risks from false-negative and false-positive results can be avoided by the NGS test.
NGS is with a bright prospect. A case report described the use of NGS for PGD recently . Several comments for the application of NGS/MPS in PGD/PGS were published [34, 35]. The cost and time of sequencing is already competitive with array tests, and the estimated reagent cost of sequencing for the detection of chromosomal abnormalities is currently less than $100. Furthermore, the cost-effectiveness of sequencing is rapidly increasing and so is the accuracy due to continuous technical improvements, including some bench-top sequencing platforms . NGS is able to achieve high coverage of mtDNA in this study. The relative copy number of mtDNA for the blastocysts with NGS testing was evaluated and a significant difference was found between the copy number of mtDNA for the euploid blastocysts and the chromosomally abnormal blastocysts. It indicated the potential of embryonic mitochondria analysis by NGS to evaluate the potential for embryonic development . Single-cell level-based methods using NGS is developing rapidly, such as single-cell exome sequencing [38, 39], single-cell genome analysis , single-cell structural variation analysis , and single-cell RNA sequencing . These methods hold promise for the further application of NGS in embryology and reproductive medicine.
In conclusion, NGS-based PGD/PGS combined with blastocyst biopsy and vitrification can be efficiently applied in clinical practice. In this study, the clinical outcome measures of the NGS cycles were the same as SNP array, and NGS was able to detect some segmental imbalances that may be omitted by SNP array, which may avoid the potential risks of false signals. NGS-based PGD/PGS should be suitable for more extensive applications and might even provide services for the populations with genetic risks.
Patients with various conditions including Robertsonian translocation, reciprocal translocation, inversion, advanced maternal age (AMA) and/or with recurrent miscarriage (RM) were recruited in the PGD/PGS program at Reproductive & Genetic Hospital of CITIC Xiangya, China from August 2011 to the current date. SNP array-based PGD/PGS cycles were performed earlier than NGS-based PGD/PGS cycles in this study. The karyotype of these patients was determined by G-banding by karyotyping peripheral blood. Blastocyst samples from all couples were obtained with Institutional Review Board (IRB) approval from BGI-Shenzhen, China and Reproductive & Genetic Hospital of CITIC Xiangya, China with reference number of LL-SC-RG-2012-001 for the IRB approval, and all couples signed a written informed consent form (Additional file 3).
Controlled ovarian hyperstimulation and fertilization
Pituitary desensitization was performed using either a long luteal Gonadotropin-releasing hormone (GnRH) agonist protocol or an antagonist protocol [43, 44] based on patient situations. Oocyte retrieval (OR) was performed 34–36 hours after human chorionic gonadotropin (hCG) injection under general anaesthesia. All eggs were fertilized by intracytoplasmic sperm injection (ICSI) 4–6 hours after OR, and normal fertilization was identified 16–18 hours after injection by the presence of two pronuclei and two polar bodies.
Embryo culture and biopsy
All embryos were cultured in sequential media (G1 and G2, Vitrolife, Goteborg, Sweden) to the blastocyst stage. On D4 after fertilization, an 18 μm hole was made in the zona pellucida of all embryos. On D5 or D6, blastocysts in which trophectoderm (TE) cells had herniated out of the zona pellucida were chosen for biopsy. Approximately 3–8 TE cells were aspirated using a biopsy pipette with 30 μm internal diameter and dissected with a Zilos TK laser (Hamilton Thorne, MA, USA). Biopsied TE cells were washed in G-MOPS medium (Vitrolife, Goteborg, Sweden). They were either used directly for Whole Genome Amplification (WGA) or stored at -20°C until WGA.
Next generation sequencing and data analysis
The biopsied TE cells were used for WGA with WGA4 Kit (Sigma Aldrich Co., St. Louis, USA) following the manufacturer’s instructions. 1-2ug of the WGA product was used for library construction according to the guidelines of TruSeq DNA Sample Prep Kits (Illumina, San Diego, USA) with a 350 bp insert size. The libraries were processed for single-end 50 bp read length sequencing using the Illumina HiSeq2000 platform.
The raw sequencing data were processed and quality controlled. The initial 20 bp sequence of each read was omitted, and the rest of the read sequence was aligned to a Human reference genome (Hg18; NCBI Build36) using the Short Oligonucleotide Analysis Package 2 (SOAP2) . After alignment, we used an algorithm developed in-house for chromosomal abnormality analysis. The procedure involved GC correction for WGA-induced bias removal, binary segmentation for locating copy number variations (CNV) breakpoints, and dynamic threshold determination for filtering of final signals . The CNVs larger than 1 megabase (>1 M) would be detected theoretically. All the results were visualized by digital karyotyping for better presentation .
The relative copy number of mtDNA for each blastocyst with NGS testing was calculated as the ratio between the mean depth of the mitochondrial DNA and the genome, which was displayed in Additional file 1: Table S3 as the Ratio (ChrM_depth/Mean_depth) and in Additional file 2: Figure S2 as the box plot figure. The Ratio (ChrM_depth/Mean_depth) of the euploid and the chromosomally abnormal blastocysts with NGS testing were used for Mann–Whitney test.
SNP array testing
WGA products were also processed for the SNP array analysis according to previous reports [21, 24]. The amplified DNA of individual embryo was hybridized to the Gene Chip Mapping Nsp I 262 K microarray (Affymetrix Inc.). Approximately 260,000 SNP signals for each sample were used for copy number analysis by Gene Chip Genotyping Analysis Software (GTYPE; Affymetrix Inc) using a smoothing size of 16 Mb to eliminate background signals. As the SNP array is unable to detect the Y chromosome, Y chromosome-specific PCR was performed for sex determination and detection of Y chromosomal abnormalities. The SNP array method had been validated before clinical application and small numbers of cells from hESC lines were used, which was described in the report .
qPCR was used to validate the results of NGS and SNP array. It was performed according to previous protocols . The WGA products of 2 YH (46, XY) single cells were used as normal control and 3 pairs of primers in chromosome 9, 12, 22 respectively were selected as internal control (Additional file 1: Table S1). Two positive samples from 2 single cells of cell lines with known karyotype as 47, XY, +21 and 45,XX,-14 were validated before clinical application (Additional file 1: Table S2). The qPCR was performed with known results of NGS and SNP arrays rather than blinded.
Blastocysts vitrification and warming
Blastocysts were vitrified after the biopsy using Kitazato vitrification solution (Kitazato Biopharma Co. Ltd. Shizuoka. Japan) and closed High Security Vitrification straws (Cryo Bio System, France). The vitrification and warming procedure was carried out according to the protocol recommended in the Kitazato vitrification kit. Each blastocyst was stored in an individual straw. After warming and dilution, blastocysts were cultured in blastocyst medium for 1–2 hours. We selected the chromosomally normal/balanced blastocysts for warming, and the surviving re-expanded blastocysts with high morphology grade were selected for transfer.
Luteal support and blastocyst transfer
Luteal support was applied in cryopreserved embryo transfer (CET) cycles. Warmed blastocysts were transferred either 5 days after ovulation of a natural cycle or 5 days after the initiation by progesterone. Briefly, 6 mg Estradiol Valerate was started from Day 3 for 10–15 days, then luteal support was applied when a satisfactory endometrial development (thickness ≥8 mm) was confirmed with ultrasound. No more than two blastocysts were transferred, and single blastocyst transfer to each patient with well cryopreserved embryos was recommended.
The main outcome measures included clinical pregnancy, implantation, miscarriage and ongoing pregnancy rates. These four outcome measures of NGS and SNP array cycles were used for Pearson’s chi-squared test (-test) to evaluate the difference between the two groups. Clinical pregnancy was confirmed when an intrauterine gestational sac with heart-beat was observed by ultrasound examination 30–40 days after embryo transfer. The amniocentesis fluid samples from fetuses or the peripheral blood samples from delivered babies were used for karyotyping to confirm the PGD/PGS results.
Availability of supporting data
The data sets supporting the results of this article are available in the European Genome-Phenome Archive repository at accession EGAD00001001037. Further details on data access are available from the GigaScience database .
Advanced maternal age
Cryopreserved embryo transfer
Comparative genomic hybridization
Copy number variation
Fluorescence in situ hybridization
Human chorionic gonadotropin
Intracytoplasmic sperm injection
Institutional Review Board
In vitro fertilization
Next generation sequencing
Preimplantation genetic diagnosis
Preimplantation genetic screening
Quantitative polymerase chain reaction
Single nucleotide polymorphism
Whole genome amplification.
Vanneste E, Voet T, Le Caignec C, Ampe M, Konings P, Melotte C, Debrock S, Amyere M, Vikkula M, Schuit F, Fryns JP, Verbeke G, D'Hooghe T, Moreau Y, Vermeesch JR: Chromosome instability is common in human cleavage-stage embryos. Nat Med. 2009, 15 (5): 577-583. 10.1038/nm.1924.
Kuliev A, Cieslak J, Ilkevitch Y, Verlinsky Y: Chromosomal abnormalities in a series of 6,733 human oocytes in preimplantation diagnosis for age-related aneuploidies. Reprod Biomed Online. 2003, 6 (1): 54-59. 10.1016/S1472-6483(10)62055-X.
Pellestor F, Andreo B, Arnal F, Humeau C, Demaille J: Maternal aging and chromosomal abnormalities: new data drawn from in vitro unfertilized human oocytes. Hum Genet. 2003, 112 (2): 195-203.
Fragouli E, Wells D, Whalley KM, Mills JA, Faed MJ, Delhanty JD: Increased susceptibility to maternal aneuploidy demonstrated by comparative genomic hybridization analysis of human MII oocytes and first polar bodies. Cytogenet Genome Res. 2006, 114 (1): 30-38. 10.1159/000091925.
Kim JW, Lee WS, Yoon TK, Seok HH, Cho JH, Kim YS, Lyu SW, Shim SH: Chromosomal abnormalities in spontaneous abortion after assisted reproductive treatment. BMC Med Genet. 2010, 11: 153-10.1186/1471-2350-11-153.
Hodes-Wertz B, Grifo J, Ghadir S, Kaplan B, Laskin CA, Glassner M, Munne S: Idiopathic recurrent miscarriage is caused mostly by aneuploid embryos. Fertil Steril. 2012, 98 (3): 675-680. 10.1016/j.fertnstert.2012.05.025.
Raziel A, Friedler S, Schachter M, Kasterstein E, Strassburger D, Ron-El R: Increased frequency of female partner chromosomal abnormalities in patients with high-order implantation failure after in vitro fertilization. Fertil Steril. 2002, 78 (3): 515-519. 10.1016/S0015-0282(02)03298-3.
Munne S, Sandalinas M, Escudero T, Fung J, Gianaroli L, Cohen J: Outcome of preimplantation genetic diagnosis of translocations. Fertil Steril. 2000, 73 (6): 1209-1218. 10.1016/S0015-0282(00)00495-7.
Verlinsky Y, Tur-Kaspa I, Cieslak J, Bernal A, Morris R, Taranissi M, Kaplan B, Kuliev A: Preimplantation testing for chromosomal disorders improves reproductive outcome of poor-prognosis patients. Reprod Biomed Online. 2005, 11 (2): 219-225. 10.1016/S1472-6483(10)60961-3.
Munne S, Chen S, Fischer J, Colls P, Zheng X, Stevens J, Escudero T, Oter M, Schoolcraft B, Simpson JL, Cohen J: Preimplantation genetic diagnosis reduces pregnancy loss in women aged 35 years and older with a history of recurrent miscarriages. Fertil Steril. 2005, 84 (2): 331-335. 10.1016/j.fertnstert.2005.02.027.
Liang L, Wang CT, Sun X, Liu L, Li M, Witz C, Williams D, Griffith J, Skorupski J, Haddad G, Gill J, Wang WH: Identification of chromosomal errors in human preimplantation embryos with oligonucleotide DNA microarray. PLoS One. 2013, 8 (4): e61838-10.1371/journal.pone.0061838.
Garrisi JG, Colls P, Ferry KM, Zheng X, Garrisi MG, Munne S: Effect of infertility, maternal age, and number of previous miscarriages on the outcome of preimplantation genetic diagnosis for idiopathic recurrent pregnancy loss. Fertil Steril. 2009, 92 (1): 288-295. 10.1016/j.fertnstert.2008.05.056.
Fiorentino F, Spizzichino L, Bono S, Biricik A, Kokkali G, Rienzi L, Ubaldi FM, Iammarrone E, Gordon A, Pantos K: PGD for reciprocal and Robertsonian translocations using array comparative genomic hybridization. Hum Reprod. 2011, 26 (7): 1925-1935. 10.1093/humrep/der082.
Schoolcraft WB, Fragouli E, Stevens J, Munne S, Katz-Jaffe MG, Wells D: Clinical application of comprehensive chromosomal screening at the blastocyst stage. Fertil Steril. 2010, 94 (5): 1700-1706. 10.1016/j.fertnstert.2009.10.015.
Rubio C, Bellver J, Rodrigo L, Bosch E, Mercader A, Vidal C, De los Santos MJ, Giles J, Labarta E, Domingo J, Crespo J, Remohí J, Pellicer A, Simón C: Preimplantation genetic screening using fluorescence in situ hybridization in patients with repetitive implantation failure and advanced maternal age: two randomized trials. Fertil Steril. 2013, 99 (5): 1400-1407. 10.1016/j.fertnstert.2012.11.041.
Wells D, Sherlock JK, Handyside AH, Delhanty JD: Detailed chromosomal and molecular genetic analysis of single cells by whole genome amplification and comparative genomic hybridisation. Nucleic Acids Res. 1999, 27 (4): 1214-1218. 10.1093/nar/27.4.1214.
Elias S: Preimplantation genetic diagnosis by comparative genomic hybridization. N Engl J Med. 2001, 345 (21): 1569-1571. 10.1056/NEJM200111223452112.
Wells D, Escudero T, Levy B, Hirschhorn K, Delhanty JD, Munne S: First clinical application of comparative genomic hybridization and polar body testing for preimplantation genetic diagnosis of aneuploidy. Fertil Steril. 2002, 78 (3): 543-549. 10.1016/S0015-0282(02)03271-5.
Fragouli E, Alfarawati S, Daphnis DD, Goodall NN, Mania A, Griffiths T, Gordon A, Wells D: Cytogenetic analysis of human blastocysts with the use of FISH, CGH and aCGH: scientific data and technical evaluation. Hum Reprod. 2011, 26 (2): 480-490. 10.1093/humrep/deq344.
Christopikou D, Tsorva E, Economou K, Shelley P, Davies S, Mastrominas M, Handyside AH: Polar body analysis by array comparative genomic hybridization accurately predicts aneuploidies of maternal meiotic origin in cleavage stage embryos of women of advanced maternal age. Hum Reprod. 2013, 28 (5): 1426-1434. 10.1093/humrep/det053.
Treff NR, Northrop LE, Kasabwala K, Su J, Levy B, Scott RT: Single nucleotide polymorphism microarray-based concurrent screening of 24-chromosome aneuploidy and unbalanced translocations in preimplantation human embryos. Fertil Steril. 2011, 95 (5): 1606-1612. 10.1016/j.fertnstert.2010.11.004. e1601-1602
Brezina PR, Benner A, Rechitsky S, Kuliev A, Pomerantseva E, Pauling D, Kearns WG: Single-gene testing combined with single nucleotide polymorphism microarray preimplantation genetic diagnosis for aneuploidy: a novel approach in optimizing pregnancy outcome. Fertil Steril. 2011, 95 (5): 1786-e1785-1788
Schoolcraft WB, Treff NR, Stevens JM, Ferry K, Katz-Jaffe M, Scott RT: Live birth outcome with trophectoderm biopsy, blastocyst vitrification, and single-nucleotide polymorphism microarray-based comprehensive chromosome screening in infertile patients. Fertil Steril. 2011, 96 (3): 638-640. 10.1016/j.fertnstert.2011.06.049.
Tan YQ, Tan K, Zhang SP, Gong F, Cheng DH, Xiong B, Lu CF, Tang XC, Luo KL, Lin G, Lu GX: Single-nucleotide polymorphism microarray-based preimplantation genetic diagnosis is likely to improve the clinical outcome for translocation carriers. Hum Reprod. 2013, 28 (9): 2581-2592. 10.1093/humrep/det271.
Jiang F, Ren J, Chen F, Zhou Y, Xie J, Dan S, Su Y, Yin B, Su W, Zhang H, Wang W, Chai X, Lin L, Guo H, Li Q, Li P, Yuan Y, Pan X, Li Y, Liu L, Chen H, Xuan Z, Chen S, Zhang C, Zhang H, Tian Z, Zhang Z, Jiang H, Zhao L, Zheng W, Li S, Li Y, Wang J, Wang J, Zhang X: Noninvasive Fetal Trisomy (NIFTY) test: an advanced noninvasive prenatal diagnosis methodology for fetal autosomal and sex chromosomal aneuploidies. BMC Med Genomics. 2012, 5: 57-10.1186/1755-8794-5-57.
Yin X, Tan K, Vajta G, Jiang H, Tan Y, Zhang C, Chen F, Chen S, Pan X, Gong C, Li X, Lin C, Gao Y, Liang Y, Yi X, Mu F, Zhao L, Peng H, Xiong B, Zhang S, Cheng D, Lu G, Zhang X, Lin G, Wang W: Massively parallel sequencing for chromosomal abnormality testing in trophectoderm cells of human blastocysts. Biol Reprod. 2013, 88 (3): 69-10.1095/biolreprod.112.106211.
Zhang C, Chen S, Yin X, Pan X, Lin G, Tan Y, Tan K, Xu Z, Hu P, Li X, Chen F, Xu X, Li Y, Zhang X, Jiang H, Wang W: A single cell level based method for copy number variation analysis by low coverage massively parallel sequencing. PLoS One. 2013, 8 (1): e54236-10.1371/journal.pone.0054236.
Harper J, Sermon K, Geraedts J, Vesela K, Harton G, Thornhill A, Pehlivan T, Fiorentino F, SenGupta S, de Die-Smulders C, Magli C, Moutou C, Wilton L: What next for preimplantation genetic screening?. Hum Reprod. 2008, 23 (3): 478-480. 10.1093/humrep/dem424.
Harper JC, Sengupta SB: Preimplantation genetic diagnosis: state of the art 2011. Hum Genet. 2012, 131 (2): 175-186. 10.1007/s00439-011-1056-z.
Wilton L, Thornhill A, Traeger-Synodinos J, Sermon KD, Harper JC: The causes of misdiagnosis and adverse outcomes in PGD. Hum Reprod. 2009, 24 (5): 1221-1228. 10.1093/humrep/den488.
Gosden R: Genetic test may lead to waste of healthy embryos. Nature. 2007, 446 (7134): 372-
Mir P, Rodrigo L, Mercader A, Buendia P, Mateu E, Milan-Sanchez M, Peinado V, Pellicer A, Remohi J, Simon C, Rubio C: False positive rate of an arrayCGH platform for single-cell preimplantation genetic screening and subsequent clinical application on day-3. J Assist Reprod Genet. 2013, 30 (1): 143-149. 10.1007/s10815-012-9918-4.
Treff NR, Fedick A, Tao X, Devkota B, Taylor D, Scott RT: Evaluation of targeted next-generation sequencing-based preimplantation genetic diagnosis of monogenic disease. Fertil Steril. 2013, 99 (5): 1377-1384. 10.1016/j.fertnstert.2012.12.018. e1376
Handyside AH: 24-chromosome copy number analysis: a comparison of available technologies. Fertil Steril. 2013, 100 (3): 595-602. 10.1016/j.fertnstert.2013.07.1965.
Martin J, Cervero A, Mir P, Martinez-Conejero JA, Pellicer A, Simon C: The impact of next-generation sequencing technology on preimplantation genetic diagnosis and screening. Fertil Steril. 2013, 99 (4): 1054-1061. 10.1016/j.fertnstert.2013.02.001. e1053
Collins FS, Hamburg MA: First FDA authorization for next-generation sequencer. N Engl J Med. 2013, 369 (25): 2369-2371. 10.1056/NEJMp1314561.
St John JC, Facucho-Oliveira J, Jiang Y, Kelly R, Salah R: Mitochondrial DNA transmission, replication and inheritance: a journey from the gamete through the embryo and into offspring and embryonic stem cells. Hum Reprod Update. 2010, 16 (5): 488-509. 10.1093/humupd/dmq002.
Xu X, Hou Y, Yin X, Bao L, Tang A, Song L, Li F, Tsang S, Wu K, Wu H, He W, Zeng L, Xing M, Wu R, Jiang H, Liu X, Cao D, Guo G, Hu X, Gui Y, Li Z, Xie W, Sun X, Shi M, Cai Z, Wang B, Zhong M, Li J, Lu Z, Gu N, Zhang X, Goodman L, Bolund L, Wang J, Yang H, Kristiansen K, Dean M, Li Y, Wang J: Single-cell exome sequencing reveals single-nucleotide mutation characteristics of a kidney tumor. Cell. 2012, 148 (5): 886-895. 10.1016/j.cell.2012.02.025.
Hou Y, Song L, Zhu P, Zhang B, Tao Y, Xu X, Li F, Wu K, Liang J, Shao D, Wu H, Ye X, Ye C, Wu R, Jian M, Chen Y, Xie W, Zhang R, Chen L, Liu X, Yao X, Zheng H, Yu C, Li Q, Gong Z, Mao M, Yang X, Yang L, Li J, Wang W, Lu Z, Gu N, Laurie G, Bolund L, Kristiansen K, Wang J, Yang H, Li Y, Zhang X, Wang J: Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm. Cell. 2012, 148 (5): 873-885. 10.1016/j.cell.2012.02.028.
Zong C, Lu S, Chapman AR, Xie XS: Genome-wide detection of single-nucleotide and copy-number variations of a single human cell. Science. 2012, 338 (6114): 1622-1626. 10.1126/science.1229164.
Voet T, Kumar P, Van Loo P, Cooke SL, Marshall J, Lin ML, Zamani Esteki M, Van der Aa N, Mateiu L, McBride DJ, Bignell GR, McLaren S, Teague J, Butler A, Raine K, Stebbings LA, Quail MA, D'Hooghe T, Moreau Y, Futreal PA, Stratton MR, Vermeesch JR, Campbell PJ: Single-cell paired-end genome sequencing reveals structural variation per cell cycle. Nucleic Acids Res. 2013, 41 (12): 6119-6138. 10.1093/nar/gkt345.
Xue Z, Huang K, Cai C, Cai L, Jiang CY, Feng Y, Liu Z, Zeng Q, Cheng L, Sun YE, Liu JY, Horvath S, Fan G: Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013, 500 (7464): 593-597. 10.1038/nature12364.
Erb TM, Wakim AN: GnRH agonist long protocol vs. a single 3-mg gnRH antagonist: a comparison of 2 protocols for pituitary down-regulation in oocyte donor-controlled ovarian hyperstimulation cycles. J Reprod Med. 2008, 53 (5): 331-337.
Orvieto R, Homburg R, Meltcer S, Rabinson J, Anteby EY, Scharf S: GnRH agonist versus GnRH antagonist in controlled ovarian hyperstimulation: their role in patients with an unfavorable prognosis a priori. Fertil Steril. 2009, 91 (4 Suppl): 1378-1380.
Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967. 10.1093/bioinformatics/btp336.
Leary RJ, Cummins J, Wang TL, Velculescu VE: Digital karyotyping. Nat Protoc. 2007, 2 (8): 1973-1986. 10.1038/nprot.2007.276.
Schmittgen TD, Livak KJ: Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008, 3 (6): 1101-1108. 10.1038/nprot.2008.73.
Tan Y, Yin X, Zhang S, Jiang H, Tan K, Li J, Xiong B, Gong F, Zhang C, Pan X, Chen F, Chen S, Gong C, Lu C, Luo K, Gu Y, Zhang X, Wang W, Xu X, Vajta G, Bolund L, Yang H, Lu G, Du Y, Lin G: Supporting materials for: “clinical outcome of preimplantation genetic diagnosis and screening using next generation sequencing”. GigaSci Database. 2014,http://dx.doi.org/10.5524/100112,
This study was supported by grant from the Major State Basic Research Program of China (No. 2012CB944901), Shenzhen Birth Defect Screening Project Lab (JZF No.  861), Key Laboratory of Project in Guangdong Province (2011A060906007) and National Natural Science Funds for Excellent Young Scholar (81222007).
The authors declare that they have no competing interests.
YT analyzed the data, signed the testing report of each patient and wrote the paper. XY wrote the paper, analyzed the data, and coordinated the work of this study. HJ coordinated the work of this study and revised the manuscript. SZ performed embryo biopsy. KT and BX performed the SNP array experiments and interpreted the results. JL, FC and XP performed the NGS experiment. XP performed the qPCR experiment and analyzed the data. CZ, SC and CG established the algorithm for CNV analysis by NGS and analyzed the NGS data. FG performed IVF clinic and embryo transfer. CL and YG performed embryo culture and intracytoplasmic sperm injection. KL contributed to patients’ management. GV and LB improved the manuscript and participated the designing. XZ, WW, XX and LZ provided support for this study. GL and HY conceived the study, provided financial support and revised the manuscript. GL and YD designed the study, provided financial support, analyzed the data and revised the paper. All authors approved for this manuscript.
Yueqiu Tan, Xuyang Yin contributed equally to this work.