Kolmogorov, Mikhail

Summary Cancer genome sequencing is essential for understanding tumor evolution and advancing precision medicine. 1 However, reference gaps and germline variants obscure detection of small and large somatic variants and methylation in repetitive regions. 1–3 It is common for tumor cells to gain or lose chromosome arms due to somatic structural changes that occur inside highly repetitive satellite DNA sequences in the centromeres. 4 To identify the full spectrum of somatic variants, including complex rearrangements, we construct and curate near-complete, haplotype-resolved assemblies of the most recent common ancestor of an early-passage broadly-consented hypodiploid pancreatic cancer cell line and matched normal tissues. The tumor assembly completely recapitulates all 35 tumor chromosomes observed with karyotyping, with multiple translocation-induced hybrid chromosomes. The hybrid chromosomes contain putative functional dicentric and fused centromeres, nested foldback inversions causing 14 breakpoints with a haplotype switch in a single event, and centromeric satellite tandem duplications up to 136 kbp. Direct comparison of tumor and normal assembly haplotypes uncovers >7,000 variants altering >1 Mbp of sequence in repetitive regions that have been hidden by reference gaps and germline variants. 44 % of somatic small variants change representation because they alter germline variants on GRCh38, impacting mutational signatures and kataegis/omikli clusters. Most somatic LINE insertions originate from two hypomethylated non-reference germline LINE insertions, highlighting their impact on insertion mutation burden. These assemblies demonstrate that centromeric, acrocentric, and telomeric regions conventionally excluded from analysis harbor extensive somatic and epigenetic changes. Resolving complete tumor genomes enables a deeper understanding of cancer structural plasticity and the endpoints of breakage-fusion-bridge cycles. These assembled, curated paired normal-tumor benchmarks will serve as a critical foundation for developing future algorithms to characterize the most intractable regions of cancer genomes.

Abstract We present chromosome-level, phased diploid genome assemblies of two widely used human fibroblast cell lines: BJ (46,XY) and IMR-90 (46,XX). Using Oxford Nanopore, PacBio HiFi, and Hi-C sequencing data, we generated assemblies spanning 5.9 and 6.0 Gbp with diploid quality values exceeding QV 60. To validate structural integrity, we developed KaryoScope, an alignment-free tool for generating computational karyograms from k-mer feature databases. We identify >50 000 structural variants relative to T2T-CHM13v2.0, the majority of which are heterozygous and cell-line-specific. Combining reference-based and de novo gene annotation, we uncover a previously unreported 1 Mbp homozygous duplication at the 16p11.2 locus in BJ, demonstrating that even karyotypically normal cell lines can harbor clinically relevant submicroscopic rearrangements. We show that mapping publicly available short-read, RNA-seq, and ChIP-seq data to sample-matched diploid assemblies substantially improves read alignment and enables haplotype phasing of 23%–28% of short reads. The BJ and IMR-90 assemblies and associated variant calls are publicly available as a resource for the research community.

ABSTRACT IDH-mutant astrocytomas maintain telomeres through the alternative lengthening of telomeres (ALT) pathway, producing extreme inter-arm telomere length heterogeneity, yet how this heterogeneity shapes structural genome evolution remains unknown. Using Oxford Nanopore long-read sequencing of 20 IDH-mutant astrocytomas, we profiled structural variants (SVs), copy number variants, extrachromosomal DNA (ecDNA) and measured allele-specific telomere lengths from individual long reads. We identified pervasive complex rearrangements, including chromothripsis and foldback events consistent with breakage-fusion-bridge cycles, and widespread ecDNAs. SV breakpoints were enriched at telomeric and centromeric regions regardless of local telomere length, revealing constitutive structural fragility. Arm-level telomere length analysis uncovered a dual-mode model: arms with short telomeres preferentially harbored breakage-associated events, while arms with long ALT-maintained telomeres were enriched for ecDNA and amplification-associated events. These findings identify chromosome-arm-specific telomere length as a determinant of structural genome evolution in ALT-driven tumors.

Abstract Background: Alternative lengthening of telomeres (ALT) is common in isocitrate dehydrogenase (IDH) mutant astrocytomas, but how telomere length variability contributes to genome instability remains unclear. We hypothesized that chromosome arm–specific telomere heterogeneity underlies distinct patterns of somatic structural variation (SV), copy-number alteration (CNA), and extrachromosomal DNA (ecDNA) formation. Methods: Oxford Nanopore long-read whole-genome sequencing (∼20–40x coverage) was performed in a clinical cohort of n = 20 IDH-mutant astrocytomas. Somatic SVs were called with Severus and filtered to a high-confidence tumor-specific set. CNAs were inferred by read-depth segmentation (DNAcopy) and further refined using a long-read CNA framework (Wakhan). Chromosome arm–specific telomere length was estimated with Telogator, and SV breakpoints were annotated for enrichment near subtelomeric and pericentromeric regions. Associations between arm-level telomere length and SV/CNA features (rate and proportion) were tested using linear mixed-effects models with patient as a random effect. High-level amplification regions were defined directly from CNA segments with high copy number, and ecDNA-like circular amplicons were subsequently reconstructed from amplified intervals using Decoil. Results: IDH-mutant astrocytomas displayed pronounced SV heterogeneity, with complex breakends (BNDs) dominating the somatic SV landscape. SVs were spatially patterned, as deletions and insertions preferentially accumulated near subtelomeric and pericentromeric domains, implicating chromosome ends and repeat-rich regions as recurrent substrates for rearrangement. At the tumor level, markedly shortened telomeres were associated with increased SV burden and CNA burden, consistent with telomere attrition as a correlate of genome-wide structural instability in highly rearranged genomes. At higher resolution, telomere length varied substantially across chromosome arms and classified distinct structural outcomes: arms with shorter telomeres were enriched for complex insertional SVs, whereas arms with relatively elongated telomeres more frequently carried focal high-level amplifications and ecDNA-like circular amplicons. Conclusions: These findings reveal a link between arm-specific telomere heterogeneity and the emergence of structural and amplification architectures in IDH-mutant astrocytomas. Telomere remodeling may influence where rearrangements and ecDNA structures arise, offering potential biomarkers of genome instability and tumor evolution. Generative AI was used to improve language clarity; all content was verified by the authors. Citation Format: Maryam Jehangir, Kristen Drucker, T. Rhyker Ranallo-Benavidez, Thomas M. Kollmeyer, Ayse Keskus, Daniel Lachance, Jeanette E. Eckel-Passow, Ofer Shoshani, Mikhail Kolmogorov, Robert B. Jenkins, Floris P. Barthel. Telomere length heterogeneity shapes structural and amplification landscapes in IDH-mutant astrocytoma [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Brain Cancer; 2026 Mar 23-25; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2026;86(6_Suppl):Abstract nr PR009.