Embryo tracking system for high-throughput sequencing-based preimplantation genetic testing.

STUDY QUESTION: Can the embryo tracking system (ETS) increase safety, efficacy and scalability of massively parallel sequencing-based preimplantation genetic testing (PGT)? SUMMARY ANSWER: Applying ETS-PGT, the chance of sample switching is decreased, while scalability and efficacy could easily be increased substantially. WHAT IS KNOWN ALREADY: Although state-of-the-art sequencing-based PGT methods made a paradigm shift in PGT, they still require labor intensive library preparation steps that makes PGT cost prohibitive and poses risks of human errors. To increase the quality assurance, efficiency, robustness and throughput of the sequencing-based assays, barcoded DNA fragments have been used in several aspects of next-generation sequencing (NGS) approach. STUDY DESIGN, SIZE, DURATION: We developed an ETS that substantially alleviates the complexity of the current sequencing-based PGT. With (n = 693) and without (n = 192) ETS, the downstream PGT procedure was performed on both bulk DNA samples (n = 563) and whole-genome amplified (WGAed) few-cell DNA samples (n = 322). Subsequently, we compared full genome haplotype landscapes of both WGAed and bulk DNA samples containing ETS or no ETS. PARTICIPANTS/MATERIALS, SETTING, METHODS: We have devised an ETS to track embryos right after whole-genome amplification (WGA) to full genome haplotype profiles. In this study, we recruited 322 WGAed DNA samples derived from IVF embryos as well as 563 bulk DNA isolated from peripheral blood of prospective parents. To determine possible interference of the ETS in the NGS-based PGT workflow, barcoded DNA fragments were added to DNA samples prior to library preparation and compared to samples without ETS. Coverages and variants were determined. MAIN RESULTS AND THE ROLE OF CHANCE: Current PGT protocols are quality sensitive and prone to sample switching. To avoid sample switching and increase throughput of PGT by sequencing-based haplotyping, six control steps should be carried out manually and checked by a second person in a clinical setting. Here, we developed an ETS approach in which one step only in the entire PGT procedure needs the four-eyes principal. We demonstrate that ETS not only precludes error-prone manual checks but also has no effect on the genomic landscape of preimplantation embryos. Importantly, our approach increases efficacy and throughput of the state-of-the-art PGT methods. LIMITATIONS, REASONS FOR CAUTION: Even though the ETS simplified sequencing-based PGT by avoiding potential errors in six steps in the protocol, if the initial assignment is not performed correctly, it could lead to cross-contamination. However, this can be detected in silico following downstream ETS analysis. Although we demonstrated an approach to evaluate purity of the ETS fragment, it is recommended to perform a pre-PGT quality control assay of the ETS amplicons with non-human DNA, such that the purity of each ETS molecule can be determined prior to ETS-PGT. WIDER IMPLICATIONS OF THE FINDINGS: The ETS-PGT approach notably increases efficacy and scalability of PGT. ETS-PGT has broad applicative value, as it can be tailored to any single- and few-cell sequencing approach where the starting specimen is scarce, as opposed to other methods that require a large number of cells as the input. Moreover, ETS-PGT could easily be adapted to any sequencing-based diagnostic method, including PGT for structural rearrangements and aneuploidies by low-pass sequencing as well as non-invasive prenatal testing. STUDY FUNDING/COMPETING INTEREST(S): M.Z.E. is supported by the EVA (Erfelijkheid Voortplanting & Aanleg) specialty program (grant no. KP111513) of Maastricht University Medical Centre (MUMC+), and the Horizon 2020 innovation (ERIN) (grant no. EU952516) of the European Commission. TRIAL REGISTRATION NUMBER: N/A.

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.

Mutation-specific effects in germline transmission of pathogenic mtDNA variants.

STUDY QUESTION: Does germline selection (besides random genetic drift) play a role during the transmission of heteroplasmic pathogenic mitochondrial DNA (mtDNA) mutations in humans? SUMMARY ANSWER: We conclude that inheritance of mtDNA is mutation-specific and governed by a combination of random genetic drift and negative and/or positive selection. WHAT IS KNOWN ALREADY: mtDNA inherits maternally through a genetic bottleneck, but the underlying mechanisms are largely unknown. Although random genetic drift is recognized as an important mechanism, selection mechanisms are thought to play a role as well. STUDY DESIGN, SIZE, DURATION: We determined the mtDNA mutation loads in 160 available oocytes, zygotes, and blastomeres of five carriers of the m.3243A>G mutation, one carrier of the m.8993T>G mutation, and one carrier of the m.14487T>C mutation. PARTICIPANTS/MATERIALS, SETTING, METHODS: Mutation loads were determined in PGD samples using PCR assays and analysed mathematically to test for random sampling effects. In addition, a meta-analysis has been performed on mutation load transmission data in the literature to confirm the results of the PGD samples. MAIN RESULTS AND THE ROLE OF CHANCE: By applying the Kimura distribution, which assumes random mechanisms, we found that mtDNA segregations patterns could be explained by variable bottleneck sizes among all our carriers (moment estimates ranging from 10 to 145). Marked differences in the bottleneck size would determine the probability that a carrier produces offspring with mutations markedly different than her own. We investigated whether bottleneck sizes might also be influenced by non-random mechanisms. We noted a consistent absence of high mutation loads in all our m.3243A>G carriers, indicating non-random events. To test this, we fitted a standard and a truncated Kimura distribution to the m.3243A>G segregation data. A Kimura distribution truncated at 76.5% heteroplasmy has a significantly better fit (P-value = 0.005) than the standard Kimura distribution. For the m.8993T>G mutation, we suspect a skewed mutation load distribution in the offspring. To test this hypothesis, we performed a meta-analysis on published blood mutation levels of offspring-mother (O-M) transmission for the m.3243A>G and m.8993T>G mutations. This analysis revealed some evidence that the O-M ratios for the m.8993T>G mutation are different from zero (P-value <0.001), while for the m.3243A>G mutation there was little evidence that the O-M ratios are non-zero. Lastly, for the m.14487T>G mutation, where the whole range of mutation loads was represented, we found no indications for selective events during its transmission. LARGE SCALE DATA: All data are included in the Results section of this article. LIMITATIONS, REASON FOR CAUTION: The availability of human material for the mutations is scarce, requiring additional samples to confirm our findings. WIDER IMPLICATIONS OF THE FINDINGS: Our data show that non-random mechanisms are involved during mtDNA segregation. We aimed to provide the mechanisms underlying these selection events. One explanation for selection against high m.3243A>G mutation loads could be, as previously reported, a pronounced oxidative phosphorylation (OXPHOS) deficiency at high mutation loads, which prohibits oogenesis (e.g. progression through meiosis). No maximum mutation loads of the m.8993T>G mutation seem to exist, as the OXPHOS deficiency is less severe, even at levels close to 100%. In contrast, high mutation loads seem to be favoured, probably because they lead to an increased mitochondrial membrane potential (MMP), a hallmark on which healthy mitochondria are being selected. This hypothesis could provide a possible explanation for the skewed segregation pattern observed. Our findings are corroborated by the segregation pattern of the m.14487T>C mutation, which does not affect OXPHOS and MMP significantly, and its transmission is therefore predominantly determined by random genetic drift. Our conclusion is that mutation-specific selection mechanisms occur during mtDNA inheritance, which has implications for PGD and mitochondrial replacement therapy. STUDY FUNDING/COMPETING INTEREST(S): This work has been funded by GROW-School of Oncology and Developmental Biology. The authors declare no competing interests.

PGD for the m.14487 T>C mitochondrial DNA mutation resulted in the birth of a healthy boy.

We report on the first PGD performed for the m.14487 T>C mitochondrial DNA (mtDNA) mutation in the MT-ND6 gene, associated with Leigh syndrome. The female carrier gave birth to a healthy baby boy at age 42. This case adds to the successes of PGD for mtDNA mutations.

Preimplantation genetic diagnosis in mitochondrial DNA disorders: challenge and success.

BACKGROUND: Mitochondrial or oxidative phosphorylation diseases are relatively frequent, multisystem disorders; in about 15% of cases they are caused by maternally inherited mitochondrial DNA (mtDNA) mutations. Because of the possible severity of the phenotype, the lack of effective treatment, and the high recurrence risk for offspring of carrier females, couples wish to prevent the transmission of these mtDNA disorders to their offspring. Prenatal diagnosis is problematic for several reasons, and concern the often poor correlation between mutation percentages and disease severity and the uncertainties about the representativeness of a fetal sample. A new option for preventing transmission of mtDNA disorders is preimplantation genetic diagnosis (PGD), which circumvents these problems by transferring an embryo below the threshold of clinical expression. METHODS: We present the data on nine PGD cycles in four female carriers of mitochondrial diseases: three mitochondrial encephalopathy, lactic acidosis and stroke-like episodes (MELAS) (m.3243A>G), and one Leigh (m.8993T>G). Our threshold for transfer after PGD is 15% for the m.3243A>G mutation and 30% for the m.8993T>G mutation. RESULTS: All four female carriers produced embryos eligible for transfer. The m.8993T>G mutation in oocytes/embryos showed more skewing than the m.3243A>G. In about 80% of the embryos the mutation load in the individual blastomeres was fairly constant (interblastomere differences <10%). However, in around 11% (in embryos with the m.3243A>G mutation only), the mutation load differed substantially (>15%) between blastomeres of a single embryo, mostly as a result of one outlier. The m.8993T>G carrier became pregnant and gave birth to a healthy son. CONCLUSIONS: PGD provides carriers of mtDNA mutations the opportunity to conceive healthy offspring.

Clinical-grade whole genome sequencing-based haplarithmisis enables all forms of preimplantation genetic testing.

High-throughput sequencing technologies have increasingly led to discovery of disease-causing genetic variants, primarily in postnatal multi-cell DNA samples. However, applying these technologies to preimplantation genetic testing (PGT) in nuclear or mitochondrial DNA from single or few-cells biopsied from in vitro fertilised (IVF) embryos is challenging. PGT aims to select IVF embryos without genetic abnormalities. Although genotyping-by-sequencing (GBS)-based haplotyping methods enabled PGT for monogenic disorders (PGT-M), structural rearrangements (PGT-SR), and aneuploidies (PGT-A), they are labour intensive, only partially cover the genome and are troublesome for difficult loci and consanguineous couples. Here, we devise a simple, scalable and universal whole genome sequencing haplarithmisis-based approach enabling all forms of PGT in a single assay. In a comparison to state-of-the-art GBS-based PGT for nuclear DNA, shallow sequencing-based PGT, and PCR-based PGT for mitochondrial DNA, our approach alleviates technical limitations by decreasing whole genome amplification artifacts by 68.4%, increasing breadth of coverage by at least 4-fold, and reducing wet-lab turn-around-time by ~2.5-fold. Importantly, this method enables trio-based PGT-A for aneuploidy origin, an approach we coin PGT-AO, detects translocation breakpoints, and nuclear and mitochondrial single nucleotide variants and indels in base-resolution.

PGD for hereditary breast and ovarian cancer: the route to universal tests for BRCA1 and BRCA2 mutation carriers.

Preimplantation Genetic Diagnosis (PGD) is a method of testing in vitro embryos as an alternative to prenatal diagnosis with possible termination of pregnancy in case of an affected child. Recently, PGD for hereditary breast and ovarian cancer caused by BRCA1 and BRCA2 mutations has found its way in specialized labs. We describe the route to universal single-cell PGD tests for carriers of BRCA1/2 mutations. Originally, mutation-specific protocols with one or two markers were set up and changed when new couples were not informative. This route of changing protocols was finalized after 2 years with universal tests for both BRCA1 and BRCA2 mutation carriers based on haplotyping of, respectively, 6 (BRCA1) and 8 (BRCA2) microsatellite markers in a multiplex PCR. Using all protocols, 30 couples had a total of 47 PGD cycles performed. Eight cycles were cancelled upon IVF treatment due to hypostimulation. Of the remaining 39 cycles, a total of 261 embryos were biopsied and a genetic diagnosis was obtained in 244 (93%). In 34 of the 39 cycles (84.6%), an embryo transfer was possible and resulted in 8 pregnancies leading to a fetal heart beat per oocyte retrieval of 20.5% and a fetal heart beat per embryonic transfer of 23.5%. The preparation time and costs for set-up and validation of tests are minimized. The informativity of microsatellite markers used in the universal PGD-PCR tests is based on CEPH and deCODE pedigrees, making the tests applicable in 90% of couples coming from these populations.

SNP array-based copy number and genotype analyses for preimplantation genetic diagnosis of human unbalanced translocations.

Preimplantation genetic diagnosis (PGD) for chromosomal rearrangements (CR) is mainly based on fluorescence in situ hybridisation (FISH). Application of this technique is limited by the number of available fluorochromes, the extensive preclinical work-up and technical and interpretative artefacts. We aimed to develop a universal, off-the-shelf protocol for PGD by combining single-nucleotide polymorphism (SNP) array-derived copy number (CN) determination and genotyping for detection of unbalanced translocations in cleavage-stage embryos. A total of 36 cleavage-stage embryos that were diagnosed as unbalanced by initial PGD FISH analysis were dissociated (n=146) and amplified by multiple displacement amplification (MDA). SNP CNs and genotypes were determined using SNP array. Epstein-Barr Virus-transformed cell lines with known CR were used for optimising the genomic smoothing (GS) length setting to increase signal to noise ratio. SNP CN analysis showed 23 embryos (64%) that were unbalanced in all blastomeres for the chromosomes involved in the translocation, 5 embryos (14%) that were normal or balanced in all blastomeres and 8 embryos (22%) that were mosaic. SNP genotyping, based on analysis of informative SNP loci with opposing homozygous parental genotypes, confirmed partial monosomies associated with inheritance of unbalanced translocation in surplus embryos. We have developed a universal MDA-SNP array technique for chromosome CN analysis in single blastomeres. SNP genotyping could confirm partial monosomies. This combination of techniques showed improved diagnostic specificity compared with FISH and may provide more reliable PGD analysis associated with higher embryo transfer rate.