STUDY QUESTION Can reduced representation genome sequencing offer an alternative to single nucleotide polymorphism (SNP) arrays as a generic and genome-wide approach for comprehensive preimplantation genetic testing for monogenic disorders (PGT-M), aneuploidy (PGT-A) and structural rearrangements (PGT-SR) in human embryo biopsy samples? SUMMARY ANSWER Reduced representation genome sequencing, with OnePGT, offers a generic, next-generation sequencing-based approach for automated haplotyping and copy-number assessment, both combined or independently, in human single blastomere and trophectoderm samples. WHAT IS KNOWN ALREADY Genome-wide haplotyping strategies, such as karyomapping and haplarithmisis, have paved the way for comprehensive PGT, i.e. leveraging PGT-M, PGT-A and PGT-SR in a single workflow. These methods are based upon SNP array technology. STUDY DESIGN, SIZE, DURATION This multi-centre verification study evaluated the concordance of PGT results for a total of 225 embryos, including 189 originally tested for a monogenic disorder and 36 tested for a translocation. Concordance for whole chromosome aneuploidies was also evaluated where whole genome copy-number reference data were available. Data analysts were kept blind to the results from the reference PGT method. PARTICIPANTS/MATERIALS, SETTING, METHODS Leftover blastomere/trophectoderm whole genome amplified (WGA) material was used, or secondary trophectoderm biopsies were WGA. A reduced representation library from WGA DNA together with bulk DNA from phasing references was processed across two study sites with the Agilent OnePGT solution. Libraries were sequenced on an Illumina NextSeq500 system, and data were analysed with Agilent Alissa OnePGT software. The embedded PGT-M pipeline utilises the principles of haplarithmisis to deduce haplotype inheritance whereas both the PGT-A and PGT-SR pipelines are based upon read-count analysis in order to evaluate embryonic ploidy. Concordance analysis was performed for both analysis strategies against the reference PGT method. MAIN RESULTS AND THE ROLE OF CHANCE PGT-M analysis was performed on 189 samples. For nine samples, the data quality was too poor to analyse further, and for 20 samples, no result could be obtained mainly due to biological limitations of the haplotyping approach, such as co-localisation of meiotic crossover events and nullisomy for the chromosome of interest. For the remaining 160 samples, 100% concordance was obtained between OnePGT and the reference PGT-M method. Equally for PGT-SR, 100% concordance for all 36 embryos tested was demonstrated. Moreover, with embryos originally analysed for PGT-M or PGT-SR for which genome-wide copy-number reference data were available, 100% concordance was shown for whole chromosome copy-number calls (PGT-A). LIMITATIONS, REASONS FOR CAUTION Inherent to haplotyping methodologies, processing of additional family members is still required. Biological limitations caused inconclusive results in 10% of cases. WIDER IMPLICATIONS OF THE FINDINGS Employment of OnePGT for PGT-M, PGT-SR, PGT-A or combined as comprehensive PGT offers a scalable platform, which is inherently generic and thereby, eliminates the need for family-specific design and optimisation. It can be considered as both an improvement and complement to the current methodologies for PGT. STUDY FUNDING/COMPETING INTEREST(S) Agilent Technologies, the KU Leuven (C1/018 to J.R.V. and T.V.) and the Horizon 2020 WIDENLIFE (692065 to J.R.V. and T.V). H.M. is supported by the Research Foundation Flanders (FWO, 11A7119N). M.Z.E, J.R.V. and T.V. are co-inventors on patent applications: ZL910050-PCT/EP2011/060211- WO/2011/157846 ‘Methods for haplotyping single cells’ and ZL913096-PCT/EP2014/068315 ‘Haplotyping and copy-number typing using polymorphic variant allelic frequencies’. T.V. and J.R.V. are co-inventors on patent application: ZL912076-PCT/EP2013/070858 ‘High-throughput genotyping by sequencing’. Haplarithmisis (‘Haplotyping and copy-number typing using polymorphic variant allelic frequencies’) has been licensed to Agilent Technologies. The following patents are pending for OnePGT: US2016275239, AU2014345516, CA2928013, CN105874081, EP3066213 and WO2015067796. OnePGT is a registered trademark. D.L., J.T. and R.L.R. report personal fees during the conduct of the study and outside the submitted work from Agilent Technologies. S.H. and K.O.F. report personal fees and other during the conduct of the study and outside the submitted work from Agilent Technologies. J.A. reports personal fees and other during the conduct of the study from Agilent Technologies and personal fees from Agilent Technologies and UZ Leuven outside the submitted work. B.D. reports grants from IWT/VLAIO, personal fees during the conduct of the study from Agilent Technologies and personal fees and other outside the submitted work from Agilent Technologies. In addition, B.D. has a patent 20160275239 - Genetic Analysis Method pending. The remaining authors have no conflicts of interest.

STUDY QUESTIONCan genome-wide haplotyping increase success following preimplantation genetic testing for a monogenic disorder (PGT-M) by including zygotes with absence of pronuclei (0PN) or the presence of only one pronucleus (1PN)?SUMMARY ANSWERGenome-wide haplotyping 0PNs and 1PNs increases the number of PGT-M cycles reaching embryo transfer (ET) by 81% and the live-birth rate by 75%.WHAT IS KNOWN ALREADYAlthough a significant subset of 0PN and 1PN zygotes can develop into balanced, diploid and developmentally competent embryos, they are usually discarded because parental diploidy detection is not part of the routine work-up of PGT-M.STUDY DESIGN, SIZE, DURATIONThis prospective cohort study evaluated the pronuclear number in 2229 zygotes from 2337 injected metaphase II (MII) oocytes in 268 cycles. PGT-M for 0PN and 1PN embryos developing into Day 5/6 blastocysts with adequate quality for vitrification was performed in 42 of the 268 cycles (15.7%). In these 42 cycles, we genome-wide haplotyped 216 good quality embryos corresponding to 49 0PNs, 15 1PNs and 152 2PNs. The reported outcomes include parental contribution to embryonic ploidy, embryonic aneuploidy, genetic diagnosis for the monogenic disorder, cycles reaching ETs, pregnancy and live birth rates (LBR) for unaffected offspring.PARTICIPANTS/MATERIALS, SETTING, METHODSBlastomere DNA was whole-genome amplified and hybridized on the Illumina Human CytoSNP12V2.1.1 BeadChip arrays. Subsequently, genome-wide haplotyping and copy-number profiling was applied to investigate the embryonic genome architecture. Bi-parental, unaffected embryos were transferred regardless of their initial zygotic PN score.MAIN RESULTS AND THE ROLE OF CHANCEA staggering 75.51% of 0PN and 42.86% of 1PN blastocysts are diploid bi-parental allowing accurate genetic diagnosis for the monogenic disorder. In total, 31% (13/42) of the PGT-M cycles reached ET or could repeat ET with an unaffected 0PN or 1PN embryo. The LBR per initiated cycle increased from 9.52 to 16.67%.LIMITATIONS, REASONS FOR CAUTIONThe clinical efficacy of the routine inclusion of 0PN and 1PN zygotes in PGT-M cycles should be confirmed in larger cohorts from multicenter studies.WIDER IMPLICATIONS OF THE FINDINGSGenome-wide haplotyping allows the inclusion of 0PN and 1PN embryos and subsequently increases the cycles reaching ET following PGT-M and potentially PGT for aneuploidy (PGT-A) and chromosomal structural rearrangements (PGT-SR). Establishing measures of clinical efficacy could lead to an update of the ESHRE guidelines which advise against the use of these zygotes.STUDY FUNDING/COMPETING INTEREST(S)SymBioSys (PFV/10/016 and C1/018 to J.R.V. and T.V.), the Horizon 2020 WIDENLIFE: 692065 to J.R.V., T.V., E.D., A.D. and M.Z.E. M.Z.E., T.V. and J.R.V. co-invented haplarithmisis (‘Haplotyping and copy-number typing using polymorphic variant allelic frequencies’), which has been licensed to Agilent Technologies. H.M. is fully supported by the (FWO) (ZKD1543-ASP/16). The authors have no competing interests to declare.

Chromoanagenesis is the process by which a single catastrophic event creates complex rearrangements confined to a single or a few chromosomes. It is usually characterized by the presence of multiple deletions and/or duplications, as well as by copy neutral rearrangements. In contrast, an array CGH screen of patients with developmental anomalies revealed three patients in which a single chromosome carries from 8 to 11 large copy number gains confined to a single chromosome or chromosomal arm, but the absence of deletions. Subsequent fluorescence in situ hybiridization and massive parallel sequencing revealed the duplicons to be clustered together in distinct locations across the altered chromosomes. Breakpoint junction sequences showed both microhomology and non-templated insertions of up to 40 bp. Hence, these patients each demonstrate a single altered chromosome of clustered insertional duplications, no deletions, and breakpoint junction sequences showing microhomology and/or non-templated insertions. These observations are difficult to reconcile with current mechanistic descriptions of chromothripsis and chromoanasynthesis. Therefore, we hypothesize those rearrangements to be of a mechanistically different origin. In addition, we suggest that large untemplated insertional sequences observed at breakpoints are driven by a non-canonical non-homologous end joining mechanism.

Chromosome instability is inherent to human IVF embryos, but the full spectrum and developmental fate of chromosome anomalies remain uncharacterized. Using haplotyping-based preimplantation genetic testing for monogenic diseases (PGT-M), we mapped the parental and mechanistic origin of common and rare genomic abnormalities in 2300 cleavage stage and 361 trophectoderm biopsies. We show that while single whole chromosome aneuploidy arises due to chromosome-specific meiotic errors in the oocyte, segmental imbalances predominantly affect paternal chromosomes, implicating sperm DNA damage in segmental aneuploidy formation. We also show that postzygotic aneuploidy affects multiple chromosomes across the genome and does not discriminate between parental homologs. In addition, 6% of cleavage stage embryos demonstrated signatures of tripolar cell division with excessive chromosome loss, however hypodiploid blastomeres can be excluded from further embryo development. This observation supports the selective-pressure hypothesis in embryos. Finally, considering that ploidy violations may constitute a significant proportion of non-viable embryos, using haplotyping-based approach to map these events might further improve IVF success rate.

Ploidy or genome‐wide chromosomal anomalies such as triploidy, diploid/triploid mixoploidy, chimerism, and genome‐wide uniparental disomy are the cause of molar pregnancies, embryonic lethality, and developmental disorders. While triploidy and genome‐wide uniparental disomy can be ascribed to fertilization or meiotic errors, the mechanisms causing mixoploidy and chimerism remain shrouded in mystery. Different models have been proposed, but all remain hypothetical and controversial, are deduced from the developmental persistent genomic constitutions present in the sample studied and lack direct evidence. New single‐cell genomic methodologies, such as single‐cell genome‐wide haplotyping, provide an extended view of the constitution of normal and abnormal embryos and have further pinpointed the existence of mixoploidy in cleavage‐stage embryos. Based on those recent findings, we suggest that genome‐wide anomalies, which persist in fetuses and patients, can for a large majority be explained by a noncanonical first zygotic cleavage event, during which maternal and paternal genomes in a single zygote, segregate to different blastomeres. This process, termed heterogoneic division, provides an overarching theoretical basis for the different presentations of mixoploidy and chimerism.

Single-cell whole-genome haplotyping allows simultaneous detection of haplotypes associated with monogenic diseases, chromosome copy-numbering and subsequently, has revealed mosaicism in embryos and embryonic stem cells. Methods, such as karyomapping and haplarithmisis, were deployed as a generic and genome-wide approach for preimplantation genetic testing (PGT) and are replacing traditional PGT methods. While current methods primarily rely on single-nucleotide polymorphism (SNP) array, we envision sequencing-based methods to become more accessible and cost-efficient. Here, we developed a novel sequencing-based methodology to haplotype and copy-number profile single cells. Following DNA amplification, genomic size and complexity is reduced through restriction enzyme digestion and DNA is genotyped through sequencing. This single-cell genotyping-by-sequencing (scGBS) is the input for haplarithmisis, an algorithm we previously developed for SNP array-based single-cell haplotyping. We established technical parameters and developed an analysis pipeline enabling accurate concurrent haplotyping and copy-number profiling of single cells. We demonstrate its value in human blastomere and trophectoderm samples as application for PGT for monogenic disorders. Furthermore, we demonstrate the method to work in other species through analyzing blastomeres of bovine embryos. Our scGBS method opens up the path for single-cell haplotyping of any species with diploid genomes and could make its way into the clinic as a PGT application.

Background During normal zygotic division, two haploid parental genomes replicate, unite and segregate into two biparental diploid blastomeres. Results Contrary to this fundamental biological tenet, we demonstrate here that parental genomes can segregate to distinct blastomeres during the zygotic division resulting in haploid or uniparental diploid and polyploid cells, a phenomenon coined heterogoneic division. By mapping the genomic landscape of 82 blastomeres from 25 bovine zygotes, we show that multipolar zygotic division is a tell-tale of whole-genome segregation errors. Based on the haplotypes and live-imaging of zygotic divisions, we demonstrate that various combinations of androgenetic, gynogenetic, diploid, and polyploid blastomeres arise via distinct parental genome segregation errors including the formation of additional paternal, private parental, or tripolar spindles, or by extrusion of paternal genomes. Hence, we provide evidence that private parental spindles, if failing to congress before anaphase, can lead to whole-genome segregation errors. In addition, anuclear blastomeres are common, indicating that cytokinesis can be uncoupled from karyokinesis. Dissociation of blastocyst-stage embryos further demonstrates that whole-genome segregation errors might lead to mixoploid or chimeric development in both human and cow. Yet, following multipolar zygotic division, fewer embryos reach the blastocyst stage and diploidization occurs frequently indicating that alternatively, blastomeres with genome-wide errors resulting from whole-genome segregation errors can be selected against or contribute to embryonic arrest. Conclusions Heterogoneic zygotic division provides an overarching paradigm for the development of mixoploid and chimeric individuals and moles and can be an important cause of embryonic and fetal arrest following natural conception or IVF.

Haplotyping is imperative for comprehensive analysis of genomes, imputation of genetic variants and interpretation of error-prone single-cell genomic data. Here we present a novel sequencing-based approach for whole-genome SNP typing of single cells, and determine genome-wide haplotypes, the copy number of those haplotypes as well as the parental and segregational origin of chromosomal aberrations from sequencing- and array-based SNP landscapes of single cells. The analytical workflow is made available as an interactive web application HiVA (https://hiva.esat.kuleuven.be).