Detection of Genomic Structural Variations Associated with Drug Sensitivity and Resistance in Acute Leukemia.

Acute leukemia is a particularly problematic collection of hematological cancers, and, while somewhat rare, the survival rate of patients is typically abysmal without bone marrow transplantation. Furthermore, traditional chemotherapies used as standard-of-care for patients cause significant side effects. Understanding the evolution of leukemia to identify novel targets and, therefore, drug treatment regimens is a significant medical need. Genomic rearrangements and other structural variations (SVs) have long been known to be causative and pathogenic in multiple types of cancer, including leukemia. These SVs may be involved in cancer initiation, progression, clonal evolution, and drug resistance, and a better understanding of SVs from individual patients may help guide therapeutic options. Here, we show the utilization of optical genome mapping (OGM) to detect known and novel SVs in the samples of patients with leukemia. Importantly, this technology provides an unprecedented level of granularity and quantitation unavailable to other current techniques and allows for the unbiased detection of novel SVs, which may be relevant to disease pathogenesis and/or drug resistance. Coupled with the chemosensitivities of these samples to FDA-approved oncology drugs, we show how an impartial integrative analysis of these diverse datasets can be used to associate the detected genomic rearrangements with multiple drug sensitivity profiles. Indeed, an insertion in the gene MUSK is shown to be associated with increased sensitivity to the clinically relevant agent Idarubicin, while partial tandem duplication events in the KMT2A gene are related to the efficacy of another frontline treatment, Cytarabine.

Analytic Validation of Optical Genome Mapping in Hematological Malignancies.

Structural variations (SVs) play a key role in the pathogenicity of hematological malignancies. Standard-of-care (SOC) methods such as karyotyping and fluorescence in situ hybridization (FISH), which have been employed globally for the past three decades, have significant limitations in terms of resolution and the number of recurrent aberrations that can be simultaneously assessed, respectively. Next-generation sequencing (NGS)-based technologies are now widely used to detect clinically significant sequence variants but are limited in their ability to accurately detect SVs. Optical genome mapping (OGM) is an emerging technology enabling the genome-wide detection of all classes of SVs at a significantly higher resolution than karyotyping and FISH. OGM requires neither cultured cells nor amplification of DNA, addressing the limitations of culture and amplification biases. This study reports the clinical validation of OGM as a laboratory-developed test (LDT) according to stringent regulatory (CAP/CLIA) guidelines for genome-wide SV detection in different hematological malignancies. In total, 60 cases with hematological malignancies (of various subtypes), 18 controls, and 2 cancer cell lines were used for this study. Ultra-high-molecular-weight DNA was extracted from the samples, fluorescently labeled, and run on the Bionano Saphyr system. A total of 215 datasets, Inc.luding replicates, were generated, and analyzed successfully. Sample data were then analyzed using either disease-specific or pan-cancer-specific BED files to prioritize calls that are known to be diagnostically or prognostically relevant. Sensitivity, specificity, and reproducibility were 100%, 100%, and 96%, respectively. Following the validation, 14 cases and 10 controls were run and analyzed using OGM at three outside laboratories showing reproducibility of 96.4%. OGM found more clinically relevant SVs compared to SOC testing due to its ability to detect all classes of SVs at higher resolution. The results of this validation study demonstrate the superiority of OGM over traditional SOC methods for the detection of SVs for the accurate diagnosis of various hematological malignancies.

Comparative Benchmarking of Optical Genome Mapping and Chromosomal Microarray Reveals High Technological Concordance in CNV Identification and Additional Structural Variant Refinement.

The recommended practice for individuals suspected of a genetic etiology for disorders including unexplained developmental delay/intellectual disability (DD/ID), autism spectrum disorders (ASD), and multiple congenital anomalies (MCA) involves a genetic testing workflow including chromosomal microarray (CMA), Fragile-X testing, karyotype analysis, and/or sequencing-based gene panels. Since genomic imbalances are often found to be causative, CMA is recommended as first tier testing for many indications. Optical genome mapping (OGM) is an emerging next generation cytogenomic technique that can detect not only copy number variants (CNVs), triploidy and absence of heterozygosity (AOH) like CMA, but can also define the location of duplications, and detect other structural variants (SVs), including balanced rearrangements and repeat expansions/contractions. This study compares OGM to CMA for clinically reported genomic variants, some of these samples also have structural characterization by fluorescence in situ hybridization (FISH). OGM was performed on IRB approved, de-identified specimens from 55 individuals with genomic abnormalities previously identified by CMA (61 clinically reported abnormalities). SVs identified by OGM were filtered by a control database to remove polymorphic variants and against an established gene list to prioritize clinically relevant findings before comparing with CMA and FISH results. OGM results showed 100% concordance with CMA findings for pathogenic variants and 98% concordant for all pathogenic/likely pathogenic/variants of uncertain significance (VUS), while also providing additional insight into the genomic structure of abnormalities that CMA was unable to provide. OGM demonstrates equivalent performance to CMA for CNV and AOH detection, enhanced by its ability to determine the structure of the genome. This work adds to an increasing body of evidence on the analytical validity and ability to detect clinically relevant abnormalities identified by CMA. Moreover, OGM identifies translocations, structures of duplications and complex CNVs intractable by CMA, yielding additional clinical utility.

Optical Genome Mapping for Oncology Applications.

Optical genome mapping (OGM) is a next-generation cytogenomic technology that has the potential to replace standard-of-care technologies used in the genetic workup of various malignancies. The ability to detect various classes of structural variations that include copy number variations, deletions, duplications, balanced and unbalanced events (insertions, inversions, and translocation) and complex genomic rearrangements in a single assay and analysis demonstrates the utility of the technology in tumor research and clinical application. Herein, we provide the methodological details for performing OGM and pre- and post-analytical quality control (QC) checks and describe critical steps that should be performed with caution, probable causes for specific QC failures, and potential method modifications that could be implemented as part of troubleshooting. The protocol description and troubleshooting guide should help new and current users of the technology to improve or troubleshoot the problems (if any) in their workflow. © 2023 Wiley Periodicals LLC. Basic Protocol: Optical genome mapping.

Optical Genome Mapping: Integrating Structural Variations for Precise Homologous Recombination Deficiency Score Calculation.

Homologous recombination deficiency (HRD) is characterized by the inability of a cell to repair the double-stranded breaks using the homologous recombination repair (HRR) pathway. The deficiency of the HRR pathway results in defective DNA repair, leading to genomic instability and tumorigenesis. The presence of HRD has been found to make tumors sensitive to ICL-inducing platinum-based therapies and poly(adenosine diphosphate [ADP]-ribose) polymerase (PARP) inhibitors (PARPi). However, there are no standardized methods to measure and report HRD phenotypes. Herein, we compare optical genome mapping (OGM), chromosomal microarray (CMA), and a 523-gene NGS panel for HRD score calculations. This retrospective study included the analysis of 196 samples, of which 10 were gliomas, 176 were hematological malignancy samples, and 10 were controls. The 10 gliomas were evaluated with both CMA and OGM, and 30 hematological malignancy samples were evaluated with both the NGS panel and OGM. To verify the scores in a larger cohort, 135 cases were evaluated with the NGS panel and 71 cases with OGM. The HRD scores were calculated using a combination of three HRD signatures that included loss of heterozygosity (LOH), telomeric allelic imbalance (TAI), and large-scale transitions (LST). In the ten glioma cases analyzed with OGM and CMA using the same DNA (to remove any tumor percentage bias), the HRD scores (mean ± SEM) were 13.2 (±4.2) with OGM compared to 3.7 (±1.4) with CMA. In the 30 hematological malignancy cases analyzed with OGM and the 523-gene NGS panel, the HRD scores were 7.6 (±2.2) with OGM compared to 2.6 (±0.8) with the 523-gene NGS panel. OGM detected 70.8% and 66.8% of additional variants that are considered HRD signatures in gliomas and hematological malignancies, respectively. The higher sensitivity of OGM to capture HRD signature variants might enable a more accurate and precise correlation with response to PARPi and platinum-based drugs. This study reveals HRD signatures that are cryptic to current standard of care (SOC) methods used for assessing the HRD phenotype and presents OGM as an attractive alternative with higher resolution and sensitivity to accurately assess the HRD phenotype.

Clinical Utility of Optical Genome Mapping and 523-Gene Next Generation Sequencing Panel for Comprehensive Evaluation of Myeloid Cancers.

The standard-of-care (SOC) for genomic testing of myeloid cancers primarily relies on karyotyping/fluorescent in situ hybridization (FISH) (cytogenetic analysis) and targeted gene panels (usually ≤54 genes) that harbor hotspot pathogenic variants (molecular genetic analysis). Despite this combinatorial approach, ~50% of myeloid cancer genomes remain cytogenetically normal, and the limited sequencing variant profiles obtained from targeted panels are unable to resolve the molecular etiology of many myeloid tumors. In this study, we evaluated the performance and clinical utility of combinatorial use of optical genome mapping (OGM) and a 523-gene next-generation sequencing (NGS) panel for comprehensive genomic profiling of 30 myeloid tumors and compared it to SOC cytogenetic methods (karyotyping and FISH) and a 54-gene NGS panel. OGM and the 523-gene NGS panel had an analytical concordance of 100% with karyotyping, FISH, and the 54-gene panel, respectively. Importantly, the IPSS-R cytogenetic risk group changed from very good/good to very poor in 22% of MDS (2/9) cases based on comprehensive profiling (karyotyping, FISH, and 54-gene panel vs. OGM and 523-gene panel), while additionally identifying six compound heterozygous events of potential clinical relevance in six cases (6/30, 20%). This cost-effective approach of using OGM and a 523-gene NGS panel for comprehensive genomic profiling of myeloid cancers demonstrated increased yield of actionable targets that can potentially result in improved clinical outcomes.

Chromosome-length genome assembly and karyotype of the endangered black-footed ferret (Mustela nigripes).

The black-footed ferret (Mustela nigripes) narrowly avoided extinction to become an oft-cited example of the benefits of intensive management, research, and collaboration to save a species through ex situ conservation breeding and reintroduction into its former range. However, the species remains at risk due to possible inbreeding, disease susceptibility, and multiple fertility challenges. Here, we report the de novo genome assembly of a male black-footed ferret generated through a combination of linked-read sequencing, optical mapping, and Hi-C proximity ligation. In addition, we report the karyotype for this species, which was used to anchor and assign chromosome numbers to the chromosome-length scaffolds. The draft assembly was ~2.5 Gb in length, with 95.6% of it anchored to 19 chromosome-length scaffolds, corresponding to the 2n = 38 chromosomes revealed by the karyotype. The assembly has contig and scaffold N50 values of 148.8 kbp and 145.4 Mbp, respectively, and is up to 96% complete based on BUSCO analyses. Annotation of the assembly, including evidence from RNA-seq data, identified 21,406 protein-coding genes and a repeat content of 37.35%. Phylogenomic analyses indicated that the black-footed ferret diverged from the European polecat/domestic ferret lineage 1.6 million yr ago. This assembly will enable research on the conservation genomics of black-footed ferrets and thereby aid in the further restoration of this endangered species.

Clinical Validation and Diagnostic Utility of Optical Genome Mapping in Prenatal Diagnostic Testing.

The standard-of-care diagnostic prenatal testing includes a combination of cytogenetic methods, such as karyotyping, fluorescence in situ hybridization (FISH), and chromosomal microarray (CMA), using either direct or cultured amniocytes or chorionic villi sampling. However, each technology has its limitations: karyotyping has a low resolution (>5 Mb), FISH is targeted, and CMA does not detect balanced structural variations (SVs). These limitations necessitate the use of multiple tests, either simultaneously or sequentially, to reach a genetic diagnosis. Optical genome mapping (OGM) is an emerging technology that can detect several classes of SVs in a single assay, but it has not been evaluated in the prenatal setting. This validation study analyzed 114 samples that were received in our laboratory for traditional cytogenetic analysis with karyotyping, FISH, and/or CMA. OGM was 100% concordant in identifying the 101 aberrations that included 29 interstitial/terminal deletions, 28 duplications, 26 aneuploidies, 6 absence of heterozygosity regions, 3 triploid genomes, 4 isochromosomes, and 1 translocation; and the method revealed the identity of 3 marker chromosomes and 1 chromosome with additional material not determined by karyotyping. In addition, OGM detected 64 additional clinically reportable SVs in 43 samples. OGM has a standardized laboratory workflow and reporting solution that can be adopted in routine clinical laboratories and demonstrates the potential to replace the current standard-of-care methods for prenatal diagnostic testing.

Genetic Predisposition to Neurological Complications in Patients with COVID-19.

Several studies have identified rare and common genetic variants associated with severe COVID-19, but no study has reported genetic determinants as predisposition factors for neurological complications. In this report, we identified rare/unique structural variants (SVs) implicated in neurological functions in two individuals with neurological manifestations of COVID-19. This report highlights the possible genetic link to the neurological symptoms with COVID-19 and calls for a collective effort to study these cohorts for a possible genetic linkage.

A contiguous de novo genome assembly of sugar beet EL10 (Beta vulgaris L.).

A contiguous assembly of the inbred 'EL10' sugar beet (Beta vulgaris ssp. vulgaris) genome was constructed using PacBio long-read sequencing, BioNano optical mapping, Hi-C scaffolding, and Illumina short-read error correction. The EL10.1 assembly was 540 Mb, of which 96.2% was contained in nine chromosome-sized pseudomolecules with lengths from 52 to 65 Mb, and 31 contigs with a median size of 282 kb that remained unassembled. Gene annotation incorporating RNA-seq data and curated sequences via the MAKER annotation pipeline generated 24,255 gene models. Results indicated that the EL10.1 genome assembly is a contiguous genome assembly highly congruent with the published sugar beet reference genome. Gross duplicate gene analyses of EL10.1 revealed little large-scale intra-genome duplication. Reduced gene copy number for well-annotated gene families relative to other core eudicots was observed, especially for transcription factors. Variation in genome size in B. vulgaris was investigated by flow cytometry among 50 individuals producing estimates from 633 to 875 Mb/1C. Read-depth mapping with short-read whole-genome sequences from other sugar beet germplasm suggested that relatively few regions of the sugar beet genome appeared associated with high-copy number variation.

Structural variant analysis of a cancer reference cell line sample using multiple sequencing technologies.

BACKGROUND: The cancer genome is commonly altered with thousands of structural rearrangements including insertions, deletions, translocation, inversions, duplications, and copy number variations. Thus, structural variant (SV) characterization plays a paramount role in cancer target identification, oncology diagnostics, and personalized medicine. As part of the SEQC2 Consortium effort, the present study established and evaluated a consensus SV call set using a breast cancer reference cell line and matched normal control derived from the same donor, which were used in our companion benchmarking studies as reference samples. RESULTS: We systematically investigated somatic SVs in the reference cancer cell line by comparing to a matched normal cell line using multiple NGS platforms including Illumina short-read, 10X Genomics linked reads, PacBio long reads, Oxford Nanopore long reads, and high-throughput chromosome conformation capture (Hi-C). We established a consensus SV call set of a total of 1788 SVs including 717 deletions, 230 duplications, 551 insertions, 133 inversions, 146 translocations, and 11 breakends for the reference cancer cell line. To independently evaluate and cross-validate the accuracy of our consensus SV call set, we used orthogonal methods including PCR-based validation, Affymetrix arrays, Bionano optical mapping, and identification of fusion genes detected from RNA-seq. We evaluated the strengths and weaknesses of each NGS technology for SV determination, and our findings provide an actionable guide to improve cancer genome SV detection sensitivity and accuracy. CONCLUSIONS: A high-confidence consensus SV call set was established for the reference cancer cell line. A large subset of the variants identified was validated by multiple orthogonal methods.

Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Clinical Validation and Diagnostic Utility of Optical Genome Mapping for Enhanced Cytogenomic Analysis of Hematological Neoplasms.

The current standard-of-care cytogenetic techniques for the analysis of hematological malignancies include karyotyping, fluorescence in situ hybridization, and chromosomal microarray, which are labor intensive and time and cost prohibitive, and they often do not reveal the genetic complexity of the tumor, demonstrating the need for alternative technology for better characterization of these tumors. Herein, we report the results from our clinical validation study and demonstrate the utility of optical genome mapping (OGM), evaluated using 92 sample runs (including replicates) that included 69 well-characterized unique samples (59 hematological neoplasms and 10 controls). The technical performance (quality control metrics) resulted in 100% first-pass rate, with analytical performance (concordance) showing a sensitivity of 98.7%, a specificity of 100%, and an accuracy of 99.2%. OGM demonstrated robust technical, analytical performance, and interrun, intrarun, and interinstrument reproducibility. The limit of detection was determined to be at 5% allele fraction for aneuploidy, translocation, interstitial deletion, and duplication. OGM identified several additional structural variations, revealing the genomic architecture in these neoplasms that provides an opportunity for better tumor classification, prognostication, risk stratification, and therapy selection. Overall, OGM has outperformed the standard-of-care tests in this study and demonstrated its potential as a first-tier cytogenomic test for hematologic malignancies.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models.

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

Optical Genome Mapping and Single Nucleotide Polymorphism Microarray: An Integrated Approach for Investigating Products of Conception.

Conventional cytogenetic analysis of products of conception (POC) is of limited utility because of failed cultures, as well as microbial and maternal cell contamination (MCC). Optical genome mapping (OGM) is an emerging technology that has the potential to replace conventional cytogenetic methods. The use of OGM precludes the requirement for culturing (and related microbial contamination). However, a high percentage of MCC impedes a definitive diagnosis, which can be addressed by an additional pre-analytical quality control step that includes histological assessment of H&E stained slides from formalin-fixed paraffin embedded (FFPE) tissue with macro-dissection for chorionic villi to enrich fetal tissue component for single nucleotide polymorphism microarray (SNPM) analysis. To improve the diagnostic yield, an integrated workflow was devised that included MCC characterization of POC tissue, followed by OGM for MCC-negative cases or SNPM with histological assessment for MCC-positive cases. A result was obtained in 93% (29/31) of cases with a diagnostic yield of 45.1% (14/31) with the proposed workflow, compared to 9.6% (3/31) and 6.4% (2/31) with routine workflow, respectively. The integrated workflow with these technologies demonstrates the clinical utility and higher diagnostic yield in evaluating POC specimens.

16p13.11p11.2 triplication syndrome: a new recognizable genomic disorder characterized by optical genome mapping and whole genome sequencing.

Highly identical segmental duplications (SDs) account for over 5% of the human genome and are enriched in the short arm of the chromosome 16. These SDs are susceptibility factors for recurrent chromosomal rearrangements mediated by non-allelic homologous recombination (NAHR). Chromosomal microarray analysis (CMA) has been widely used as the first-tier test for individuals with developmental disabilities and/or congenital anomalies and several genomic disorders involving the 16p-arm have been identified with this technique. However, the resolution of CMA and the limitations of short-reads whole genome sequencing (WGS) technology do not allow the full characterization of the most complex chromosomal rearrangements. Herein, we report on two unrelated patients with a de novo 16p13.11p11.2 triplication associated with a 16p11.2 duplication, detected by CMA. These patients share a similar phenotype including hypotonia, severe neurodevelopmental delay with profound speech impairment, hyperkinetic behavior, conductive hearing loss, and distinctive facial features. Short-reads WGS could not map precisely any of the rearrangement's breakpoints that lie within SDs. We used optical genome mapping (OGM) to determine the relative orientation of the triplicated and duplicated segments as well as the genomic positions of the breakpoints, allowing us to propose a mechanism involving recombination between allelic SDs and a NAHR event. In conclusion, we report a new clinically recognizable genomic disorder. In addition, the mechanism of these complex chromosomal rearrangements involving SDs could be unraveled by OGM.

Optical genome mapping identifies rare structural variations as predisposition factors associated with severe COVID-19.

Impressive global efforts have identified both rare and common gene variants associated with severe COVID-19 using sequencing technologies. However, these studies lack the sensitivity to accurately detect several classes of variants, especially large structural variants (SVs), which account for a substantial proportion of genetic diversity including clinically relevant variation. We performed optical genome mapping on 52 severely ill COVID-19 patients to identify rare/unique SVs as decisive predisposition factors associated with COVID-19. We identified 7 SVs involving genes implicated in two key host-viral interaction pathways: innate immunity and inflammatory response, and viral replication and spread in nine patients, of which SVs in STK26 and DPP4 genes are the most intriguing candidates. This study is the first to systematically assess the potential role of SVs in the pathogenesis of COVID-19 severity and highlights the need to evaluate SVs along with sequencing variants to comprehensively associate genomic information with interindividual variability in COVID-19 phenotypes.

Constitutional chromothripsis of the APC locus as a cause of genetic predisposition to colon cancer.

PURPOSE: Approximately 20% of patients with clinical familial adenomatous polyposis (FAP) remain unsolved after molecular genetic analysis of the APC and other polyposis genes, suggesting additional pathomechanisms. METHODS: We applied multidimensional genomic analysis employing chromosomal microarray profiling, optical mapping, long-read genome and RNA sequencing combined with FISH and standard PCR of genomic and complementary DNA to decode a patient with an attenuated FAP that had remained unsolved by Sanger sequencing and multigene panel next-generation sequencing for years. RESULTS: We identified a complex 3.9 Mb rearrangement involving 14 fragments from chromosome 5q22.1q22.3 of which three were lost, 1 reinserted into chromosome 5 and 10 inserted into chromosome 10q21.3 in a seemingly random order and orientation thus fulfilling the major criteria of chromothripsis. The rearrangement separates APC promoter 1B from the coding ORF (open reading frame) thus leading to allele-specific downregulation of APC mRNA. The rearrangement also involves three additional genes implicated in the APC-Axin-GSK3B-ß-catenin signalling pathway. CONCLUSIONS: Based on comprehensive genomic analysis, we propose that constitutional chromothripsis dampening APC expression, possibly modified by additional APC-Axin-GSK3B-ß-catenin pathway disruptions, underlies the patient's clinical phenotype. The combinatorial approach we deployed provides a powerful tool set for deciphering unsolved familial polyposis and potentially other tumour syndromes and monogenic diseases.