Genome-wide haplotyping embryos developing from 0PN and 1PN zygotes increases transferrable embryos in PGT-M.

STUDY QUESTION Can 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 ANSWER Genome-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 ALREADY Although 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, DURATION This 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, METHODS Blastomere 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 CHANCE A 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 CAUTION The 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 FINDINGS Genome-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

Parental genomes segregate into distinct blastomeres during multipolar zygotic divisions leading to mixoploid and chimeric blastocysts

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 and copy-number profiling of single cells by massive parallel sequencing

status: publishe

Single-cell genome-wide concurrent haplotyping and copy-number profiling through genotyping-by-sequencing

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

Genome-wide haplotyping embryos developing from 0PN and 1PN zygotes increases transferrable embryos in PGT-M

STUDY QUESTION: Can 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 ANSWER: Genome-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 ALREADY: Although 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, DURATION: This 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, METHODS: Blastomere 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 CHANCE: A 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 CAUTION: The 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 FINDINGS: Genome-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.status: publishe

Additional file 2 of Parental genomes segregate into distinct blastomeres during multipolar zygotic divisions leading to mixoploid and chimeric blastocysts

Additional file 2: Figure S2. Analysis of blastomeres following multipolar zygotic division. A) Interpretation of haplarithm plots. Overview of chromosome-wise haplarithm patterns for distinct genomic constitutions (i.e. biparental disomy, paternal monosomy and paternal meiotic/dispermic uniparental heterodisomy). Corresponding whole-genome errors (i.e. biparental diploid or androgenetic) are characterized by the manifestation of those patterns throughout the (majority of the) genome. Defined single-cell BAF values of the segmented P1, P2, M1 and M2, form haplotype blocks, demarcated by pairwise breakpoints, i.e., homologous recombinations. Haplotype blocks, as well as the distance between the P1-P2 or M1-M2 in the paternal and maternal haplarithm, respectively, and the positioning of homologous recombinations, denote the origin and nature of copy number. The normalized LogR- values are integrated with haplarithm profiles for copy number profiling. Principles of interpretation are according to [64]. B) An overview of haplarithm profiles of 82 blastomeres and two fragments (grey squares) is depicted per category of whole-genome segregation profiles, as discussed in the main text. Each embryo is identified by a description at the top left of the embryo ID and cross (EmbryoID_Embryocross). At the top right, three chronological time-lapse images of the cleaving zygote are depicted. From left to right, the pictures show the initiation of the cleavage furrow, the ongoing first division and the embryo immediately after cleavage and before cell isolation (when video available). For each embryo, a schematic representation of likely steps leading to the genomic profile of each blastomere (B1-B4) or fragment (F1) is given. Chromosome-wise interpretation (1 - X) per blastomere is visualized in the bar above the haplarithm plots (see legend). Below each bar, the paternal haplarithm (pat-BAF), the maternal haplarithm (mat-BAF) and the normalized LogR values (LogR) are depicted. Paternal cross-over sites are depicted by the arrows (black, green or blue). A combination of parental cross-over sites in one blastomere or different cross-over sites in blastomeres of the same embryo uncover polyspermic fertilization or a meiotic error. Maternal cross-over sites (red, pink, orange) were only depicted in gynogenetic blastomeres and in case of whole-genome meiotic errors

Haplotyping-based preimplantation genetic testing reveals parent-of-origin specific mechanisms of aneuploidy formation.

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

Parental genomes segregate into distinct blastomeres during multipolar zygotic divisions leading to mixoploid and chimeric blastocysts

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

Multi-centre evaluation of a comprehensive preimplantation genetic test through haplotyping-by-sequencing

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.status: publishe