Abyzov, Alexej

Abstract Somatic mutations that arise post-zygotically create genetic diversity among normal human cells and provide key insights into human development and aging. Fibroblast-derived induced pluripotent stem cells (iPSCs) have proved to be a useful system for disease modelling; however, due to their clonal nature, iPSC lines carry somatic mutations inherited from the founder cells, raising concerns about their genomic integrity. At the same time, this clonality enables single-cell–level discovery of somatic mutations and the reconstruction of developmental lineages. In living individuals, though, this approach requires invasive biopsies and is limited to skin-derived lineages. Here, we generated 33 urine-derived iPSC lines from four males representing two father– son relationships, performed shallow whole-genome sequencing of the lines and analyzed somatic mutations. Derived iPSCs representing single cells from urine carried a few hundred of somatic single-nucleotide variants per genome, dominated by endogenous, clock-like mutational signatures and lacking environmental imprints such as UV-associated mutations. Copy-number analysis identified somatic CNVs in most of the lines and revealed higher CNV burdens in fathers than in sons, consistent with age-related structural mosaicism. Shared mutations across lines enabled reconstruction of cell lineage phylogenetic trees. In summary, urine-derived iPSCs showed genomic alterations comparable to those in fibroblast-derived iPSC lines and represent a valuable non-invasive alternative for disease modeling. Overall, this study provides the first genome-wide characterization of somatic mutations in urine-derived iPSCs and establishes them as a practical and non-invasive platform for charting somatic mutation landscapes and tracing developmental lineages in living humans.

Abstract Background Current literature suggestsisocitrate dehydrogenase (IDH)-mutant astrocytoma contains several molecular subgroups. In this study, we are interested in determining the connection between different molecular subgroups with grade and/or survival. Methods A cohort of 470 Mayo Clinic adult patients (≥18 years, 56.2% male) with primary IDH-mutant astrocytoma diagnosed by World Health Organization (WHO) 2021 criteria were examined. Results were validated in an independent cohort of 614 Mayo Clinic Neuropathology consult patients and 235 The Cancer Genome Atlas (TCGA) patients. Results The Mayo Clinic Practice cohort confirmed the association of CDKN2A/B deletion with overall survival (OS, homozygous vs hemizygous vs intact, 2.7 vs 9.6 vs 17.2 years, P < .001). Phosphatase and tensin homolog (PTEN) deletion was also associated with poor OS (7.3 vs 17.4 years, P < .001). Increased number of copy number alterations was associated with OS (continuous variable, HR = 1.027, P < .001). Carrying one or more copies of the germline risk allele at rs55705857 was associated with earlier age of onset (median age 33 vs 35 years, P = .01), and a shorter OS after adjusting for age, grade, sex and treatment (HR = 1.81, P = .007). The Mayo Clinic Neuropathology Consult cohort and TCGA were utilized to validate age of onset and survival, respectively. Unsupervised clustering of the copy number alterations identified several clinically significant groups that may define pathways to disease progression. Losses of chromosomes 11p, 13q, 1p, and 10q were all associated with reduced overall survival in the Mayo Clinic cohort. Conclusions Patients with hemizygous loss of CDKN2A/B, loss of PTEN, increased number of copy number alterations, specific chromosomal arm losses or rs55705857 germline risk allele have reduced overall survival.

Xeroderma Pigmentosum group C (XP-C) is an autosomal recessive disorder caused by mutations in the XPC gene, leading to defective nucleotide excision repair. This defect leads to genomic instability and a profound cancer predisposition. To model this disease, we generated induced pluripotent stem cells (iPSCs) from an XP-C patient carrying a novel homozygous nonsense mutation in the XPC gene (c.1830C>A). The resulting iPSCs demonstrated typical pluripotent characteristics, including expression of key markers and trilineage differentiation capability. However, genomic assessment revealed progressive karyotypic instability during extended culture. While initial whole-genome sequencing detected no major chromosomal abnormalities, subsequent G-banding analysis identified acquired trisomy 12 in two lines (CL12 and CL27) and a derivative X chromosome in a third line (CL30). These abnormalities were absent in early-passage analyses, indicating that they were acquired and selected for during extended culture. The acquisition of a derivative X chromosome in CL30, alongside recurrent trisomy 12, represents a unique cytogenetic signature likely attributable to the underlying XPC defect. We hypothesize that the loss of GG-NER creates a permissive genomic environment, accelerating the accumulation of DNA damage and chromosomal missegregation under replicative stress. This temporal divergence in genomic integrity highlights how culture pressures drive chromosomal evolution in XP-C iPSCs independently of initial reprogramming. Our findings emphasize that XP-C iPSCs require continuous genomic surveillance and provide a model for investigating how DNA repair deficiencies interact with in vitro culture stress.

Somatic mosaicism is increasingly recognized as a fundamental feature of human biology, yet the detection of somatic mutations remains challenging. The SMaHT Network conducted four large-scale benchmarking experiments involving cell-lines and donor tissues, to evaluate sequencing technologies, experimental approaches, and computational methods for detecting different types of somatic mutations, generating community resource with >1,000X short-read and 100-400X long-read data for each of the nine analyzed samples. We determined effective strategies for utilizing short- and long-reads sequencing for mutation detection and demonstrated that using donor-specific assemblies and human pangenome improved calling, extending mutation catalogs to challenging genomic regions. We benchmarked six duplex technologies and showed that single-cell sequencing resolves cell type-specific mutational patterns and heterogeneity. Our results indicate that bulk, single-cell, and duplex analyses are complementary - and leveraging all three provides comprehensive characterization of mosaicism within tissues. Together, these findings provide a roadmap for accurate, genome-wide somatic mutation discovery and analysis.

Germline mutations modulate the risk of developing schizophrenia (SCZ). Much less is known about the role of mosaic somatic mutations in the context of SCZ. Deep (239×) whole-genome sequencing (WGS) of brain neurons from 61 SCZ cases and 25 controls postmortem identified mutations occurring during prenatal neurogenesis. SCZ cases showed increased somatic variants in open chromatin, with increased mosaic CpG transversions (CpG>GpG) and T>G mutations at transcription factor binding sites (TFBSs) overlapping open chromatin, a result not seen in controls. Some of these variants alter gene expression, including SCZ risk genes and genes involved in neurodevelopment. Although these mutational processes can reflect a difference in factors indirectly involved in disease, increased somatic mutations at developmental TFBSs could also potentially contribute to SCZ.

We developed a generally applicable method, CRISPR/Cas9-targeted long-read sequencing (CTLR-Seq), to resolve, haplotype-specifically, the large and complex regions in the human genome that had been previously impenetrable to sequencing analysis, such as large segmental duplications (SegDups) and their associated genome rearrangements. CTLR-Seq combines in vitro Cas9-mediated cutting of the genome and pulse-field gel electrophoresis to isolate intact large (i.e., up to 2,000 kb) genomic regions that encompass previously unresolvable genomic sequences. These targets are then sequenced (amplification-free) at high on-target coverage using long-read sequencing, allowing for their complete sequence assembly. We applied CTLR-Seq to the SegDup-mediated rearrangements that constitute the boundaries of, and give rise to, the 22q11.2 Deletion Syndrome (22q11DS), the most common human microdeletion disorder. We then performed de novo assembly to resolve, at base-pair resolution, the full sequence rearrangements and exact chromosomal breakpoints of 22q11.2DS (including all common subtypes). Across multiple patients, we found a high degree of variability for both the rearranged SegDup sequences and the exact chromosomal breakpoint locations, which coincide with various transposons within the 22q11.2 SegDups, suggesting that 22q11DS can be driven by transposon-mediated genome recombination. Guided by CTLR-Seq results from two 22q11DS patients, we performed three-dimensional chromosomal folding analysis for the 22q11.2 SegDups from patient-derived neurons and astrocytes and found chromosome interactions anchored within the SegDups to be both cell type-specific and patient-specific. Lastly, we demonstrated that CTLR-Seq enables cell-type specific analysis of DNA methylation patterns within the deletion haplotype of 22q11DS.

Abstract Summary Copy number variation (CNV) and alteration (CNA) analysis is a crucial component in many genomic studies and its applications span from basic research to clinic diagnostics and personalized medicine. CNVpytor is a tool featuring a read depth-based caller and combined read depth and B-allele frequency (BAF) based 2D caller to find CNVs and CNAs. The tool stores processed intermediate data and CNV/CNA calls in a compact HDF5 file—pytor file. Here, we describe a new track in igv.js that utilizes pytor and whole genome variant files as input for on-the-fly read depth and BAF visualization, CNV/CNA calling and analysis. Embedding into HTML pages and Jupiter Notebooks enables convenient remote data access and visualization simplifying interpretation and analysis of omics data. Availability and implementation The CNVpytor track is integrated with igv.js and available at https://github.com/igvteam/igv.js. The documentation is available at https://github.com/igvteam/igv.js/wiki/cnvpytor. Usage can be tested in the IGV-Web app at https://igv.org/app and also on https://github.com/abyzovlab/CNVpytor.

The CRISPR-Cas9 system has enabled researchers to precisely modify/edit the sequence of a genome. A typical editing experiment consists of two steps: (1) editing cultured cells; (2) cell cloning and selection of clones with and without intended edit, presumed to be isogenic. The application of CRISPR-Cas9 system may result in off-target edits, whereas cloning will reveal culture-acquired mutations. We analyzed the extent of the former and the latter by whole genome sequencing in three experiments involving separate genomic loci and conducted by three independent laboratories. In all experiments we hardly found any off-target edits, whereas detecting hundreds to thousands of single nucleotide mutations unique to each clone after relatively short culture of 10–20 passages. Notably, clones also differed in copy number alterations (CNAs) that were several kb to several mb in size and represented the largest source of genomic divergence among clones. We suggest that screening of clones for mutations and CNAs acquired in culture is a necessary step to allow correct interpretation of DNA editing experiments. Furthermore, since culture associated mutations are inevitable, we propose that experiments involving derivation of clonal lines should compare a mix of multiple unedited lines and a mix of multiple edited lines.

Abstract Mosaic mutations can be used to track cell ancestries and reconstruct high-resolution lineage trees during cancer progression and during development, starting from the first cell divisions of the zygote. However, this approach requires sampling and analyzing the genomes of multiple cells, which can be redundant in lineage representation, limiting the scalability of the approach. We describe a strategy for cost- and time-efficient lineage reconstruction using clonal induced pluripotent stem cell lines from human skin fibroblasts. The approach leverages shallow sequencing coverage to assess the clonality of the lines, clusters redundant lines and sums their coverage to accurately discover mutations in the corresponding lineages. Only a fraction of lines needs to be sequenced to high coverage. We demonstrate the effectiveness of this approach for reconstructing lineage trees during development and in hematologic malignancies. We discuss and propose an optimal experimental design for reconstructing lineage trees.