Single-strand mismatch and damage patterns revealed by single-molecule DNA sequencing

Neuropathologically directed profiling of PRNP somatic and germline variants in sporadic human prion disease

Creutzfeldt-Jakob Disease (CJD), the most common human prion disease, is associated with pathologic misfolding of the prion protein (PrP), encoded by the PRNP gene. Of human prion disease cases, < 1% were transmitted by misfolded PrP, ~ 15% are inherited, and ~ 85% are sporadic (sCJD). While familial cases are inherited through germline mutations in PRNP, the cause of sCJD is unknown. Somatic mutations have been hypothesized as a cause of sCJD, and recent studies have revealed that somatic mutations accumulate in neurons during aging. To investigate the hypothesis that somatic mutations in PRNP may underlie sCJD, we performed deep DNA sequencing of PRNP in 205 sCJD cases and 170 age-matched non-disease controls. We included 5 cases of Heidenhain variant sporadic CJD (H-sCJD), where visual symptomatology and neuropathology implicate localized initiation of prion formation, and examined multiple regions across the brain including in the affected occipital cortex. We employed Multiple Independent Primer PCR Sequencing (MIPP-Seq) with a median depth of > 5000× across the PRNP coding region and analyzed for variants using MosaicHunter. An allele mixing experiment showed positive detection of variants in bulk DNA at a variant allele fraction (VAF) as low as 0.2%. We observed multiple polymorphic germline variants among individuals in our cohort. However, we did not identify bona fide somatic variants in sCJD, including across multiple affected regions in H-sCJD, nor in control individuals. Beyond our stringent variant-identification pipeline, we also analyzed VAFs from raw sequencing data, and observed no evidence of prion disease enrichment for the known germline pathogenic variants P102L, D178N, and E200K. The lack of PRNP pathogenic somatic mutations in H-sCJD or the broader cohort of sCJD suggests that clonal somatic mutations may not play a major role in sporadic prion disease. With H-sCJD representing a localized presentation of neurodegeneration, this serves as a test of the potential role of clonal somatic mutations in genes known to cause familial neurodegeneration.

NeMu: a comprehensive pipeline for accurate reconstruction of neutral mutation spectra from evolutionary data

The recognized importance of mutational spectra in molecular evolution is yet to be fully exploited beyond human cancer studies and model organisms. The wealth of intraspecific polymorphism data in the GenBank repository, covering a broad spectrum of genes and species, presents an untapped opportunity for detailed mutational spectrum analysis. Existing methods fall short by ignoring intermediate substitutions on the inner branches of phylogenetic trees and lacking the capability for cross-species mutational comparisons. To address these challenges, we present the NeMu pipeline, available at https://nemu-pipeline.com, a tool grounded in phylogenetic principles designed to provide comprehensive and scalable analysis of mutational spectra. Utilizing extensive sequence data from numerous available genome projects, NeMu rapidly and accurately reconstructs the neutral mutational spectrum. This tool, facilitating the reconstruction of gene- and species-specific mutational spectra, contributes to a deeper understanding of evolutionary mechanisms across the broad spectrum of known species.

Neuropathologically-directed profiling of PRNP somatic and germline variants in sporadic human prion disease

Creutzfeldt-Jakob Disease (CJD), the most common human prion disease, is associated with pathologic misfolding of the prion protein (PrP), encoded by the PRNP gene. Of human prion disease cases, ~1% were transmitted by misfolded PrP, ~15% are inherited, and ~85% are sporadic (sCJD). While familial cases are inherited through germline mutations in PRNP, the cause of sCJD is unknown. Somatic mutations have been hypothesized as a cause of sCJD, and recent studies have revealed that somatic mutations accumulate in neurons during aging. To investigate the hypothesis that somatic mutations in PRNP may underlie sCJD, we performed deep DNA sequencing of PRNP in 205 sCJD cases and 170 age-matched non-disease controls. We included 5 cases of Heidenhain variant sporadic CJD (H-sCJD), where visual symptomatology and neuropathology implicate focal initiation of prion formation, and examined multiple regions across the brain including in the affected occipital cortex. We employed Multiple Independent Primer PCR Sequencing (MIPP-Seq) with a median depth of >5,000X across the PRNP coding region and analyzed for variants using MosaicHunter. An allele mixing experiment showed positive detection of variants in bulk DNA at a variant allele fraction (VAF) as low as 0.2%. We observed multiple polymorphic germline variants among individuals in our cohort. However, we did not identify bona fide somatic variants in sCJD, including across multiple affected regions in H-sCJD, nor in control individuals. Beyond our stringent variant-identification pipeline, we also analyzed VAFs from raw sequencing data, and observed no evidence of prion disease enrichment for the known germline pathogenic variants P102L, D178N, and E200K. The lack of PRNP pathogenic somatic mutations in H-sCJD or the broader cohort of sCJD suggests that clonal somatic mutations may not play a major role in sporadic prion disease. With H-sCJD representing a focal presentation of neurodegeneration, this serves as a test of the potential role of clonal somatic mutations in genes known to cause familial neurodegeneration.

Direct transposition of native DNA for sensitive multimodal single-molecule sequencing

Concurrent readout of sequence and base modifications from long unamplified DNA templates by Pacific Biosciences of California (PacBio) single-molecule sequencing requires large amounts of input material. Here we adapt Tn5 transposition to introduce hairpin oligonucleotides and fragment (tagment) limiting quantities of DNA for generating PacBio-compatible circular molecules. We developed two methods that implement tagmentation and use 90-99% less input than current protocols: (1) single-molecule real-time sequencing by tagmentation (SMRT-Tag), which allows detection of genetic variation and CpG methylation; and (2) single-molecule adenine-methylated oligonucleosome sequencing assay by tagmentation (SAMOSA-Tag), which uses exogenous adenine methylation to add a third channel for probing chromatin accessibility. SMRT-Tag of 40 ng or more human DNA (approximately 7,000 cell equivalents) yielded data comparable to gold standard whole-genome and bisulfite sequencing. SAMOSA-Tag of 30,000-50,000 nuclei resolved single-fiber chromatin structure, CTCF binding and DNA methylation in patient-derived prostate cancer xenografts and uncovered metastasis-associated global epigenome disorganization. Tagmentation thus promises to enable sensitive, scalable and multimodal single-molecule genomics for diverse basic and clinical applications.

Ultra-sensitive molecular residual disease detection through whole genome sequencing with single-read error correction

While whole genome sequencing (WGS) of cell-free DNA (cfDNA) holds enormous promise for molecular residual disease (MRD) detection, its performance is limited by WGS error rate. Here we introduce AccuScan, an efficient cfDNA WGS technology that enables genome-wide error correction at single read level, achieving an error rate of 4.2×10 -7 , which is about two orders of magnitude lower than a read-centric de-noising method. When applied to MRD detection, AccuScan demonstrated analytical sensitivity down to 10 -6 circulating tumor allele fraction at 99% sample level specificity. In colorectal cancer, AccuScan showed 90% landmark sensitivity for predicting relapse. It also showed robust MRD performance with esophageal cancer using samples collected as early as 1 week after surgery, and predictive value for immunotherapy monitoring with melanoma patients. Overall, AccuScan provides a highly accurate WGS solution for MRD, empowering circulating tumor DNA detection at parts per million range without high sample input nor personalized reagents.

One Sentence Summary

AccuScan showed remarkable ultra-low limit of detection with a short turnaround time, low sample requirement and a simple workflow for MRD detection.

Disentangling sources of clock-like mutations in germline and soma

The rates of mutations vary across cell types. To identify causes of this variation, mutations are often decomposed into a combination of the single base substitution (SBS) "signatures" observed in germline, soma and tumors, with the idea that each signature corresponds to one or a small number of underlying mutagenic processes. Two such signatures turn out to be ubiquitous across cell types: SBS signature 1, which consists primarily of transitions at methylated CpG sites caused by spontaneous deamination, and the more diffuse SBS signature 5, which is of unknown etiology. In cancers, the number of mutations attributed to these two signatures accumulates linearly with age of diagnosis, and thus the signatures have been termed "clock-like." To better understand this clock-like behavior, we develop a mathematical model that includes DNA replication errors, unrepaired damage, and damage repaired incorrectly. We show that mutational signatures can exhibit clock-like behavior because cell divisions occur at a constant rate and/or because damage rates remain constant over time, and that these distinct sources can be teased apart by comparing cell lineages that divide at different rates. With this goal in mind, we analyze the rate of accumulation of mutations in multiple cell types, including soma as well as male and female germline. We find no detectable increase in SBS signature 1 mutations in neurons and only a very weak increase in mutations assigned to the female germline, but a significant increase with time in rapidly-dividing cells, suggesting that SBS signature 1 is driven by rounds of DNA replication occurring at a relatively fixed rate. In contrast, SBS signature 5 increases with time in all cell types, including post-mitotic ones, indicating that it accumulates independently of cell divisions; this observation points to errors in DNA repair as the key underlying mechanism. Thus, the two "clock-like" signatures observed across cell types likely have distinct origins, one set by rates of cell division, the other by damage rates.