Somatic mutation landscapes at single-molecule resolution

Methods and applications of genome-wide profiling of DNA damage and rare mutations

DNA damage is a threat to genome integrity and can be a cause of many human diseases, owing to either changes in the chemical structure of DNA or conversion of the damage into a mutation, that is, a permanent change in DNA sequence. Determining the exact positions of DNA damage and ensuing mutations in the genome are important for identifying mechanisms of disease aetiology when characteristic mutations are prevalent and probably causative in a particular disease. However, this approach is challenging particularly when levels of DNA damage are low, for example, as a result of chronic exposure to environmental agents or certain endogenous processes, such as the generation of reactive oxygen species. Over the past few years, a comprehensive toolbox of genome-wide methods has been developed for the detection of DNA damage and rare mutations at single-nucleotide resolution in mammalian cells. Here, we review and compare these methods, describe their current applications and discuss future research questions that can now be addressed.

Mutational signature analyses in multi-child families reveal sources of age-related increases in human germline mutations

Whole-genome sequencing studies of parent-offspring trios have provided valuable insights into the potential impact of de novo mutations (DNMs) on human health and disease. However, the molecular mechanisms that drive DNMs are unclear. Studies with multi-child families can provide important insight into the causes of inter-family variability in DNM rates but they are highly limited. We characterized 2479 de novo single nucleotide variants (SNVs) in 13 multi-child families of Mexican-American ethnicity. We observed a strong paternal age effect on validated de novo SNVs with extensive inter-family variability in the yearly rate of increase. Children of older fathers showed more C > T transitions at CpG sites than children from younger fathers. Validated SNVs were examined against one cancer (COSMIC) and two non-cancer (human germline and CRISPR-Cas 9 knockout of human DNA repair genes) mutational signature databases. These analyses suggest that inaccurate DNA mismatch repair during repair initiation and excision processes, along with DNA damage and replication errors, are major sources of human germline de novo SNVs. Our findings provide important information for understanding the potential sources of human germline de novo SNVs and the critical role of DNA mismatch repair in their genesis.

Therapeutic targeting of senescent cells in the CNS

Senescent cells accumulate throughout the body with advanced age, diseases and chronic conditions. They negatively impact health and function of multiple systems, including the central nervous system (CNS). Therapies that target senescent cells, broadly referred to as senotherapeutics, recently emerged as potentially important treatment strategies for the CNS. Promising therapeutic approaches involve clearing senescent cells by disarming their pro-survival pathways with 'senolytics'; or dampening their toxic senescence-associated secretory phenotype (SASP) using 'senomorphics'. Following the pioneering discovery of first-generation senolytics dasatinib and quercetin, dozens of additional therapies have been identified, and several promising targets are under investigation. Although potentially transformative, senotherapies are still in early stages and require thorough testing to ensure reliable target engagement, specificity, safety and efficacy. The limited brain penetrance and potential toxic side effects of CNS-acting senotherapeutics pose challenges for drug development and translation to the clinic. This Review assesses the potential impact of senotherapeutics for neurological conditions by summarizing preclinical evidence, innovative methods for target and biomarker identification, academic and industry drug development pipelines and progress in clinical trials.

Somatic mutations in autoinflammatory and autoimmune disease

Somatic mutations (also known as acquired mutations) are emerging as common, age-related processes that occur in all cells throughout the body. Somatic mutations are canonically linked to malignant processes but over the past decade have been increasingly causally connected to benign diseases including rheumatic conditions. Here we outline the contribution of somatic mutations to complex and monogenic immunological diseases with a detailed review of unique aspects associated with such causes. Somatic mutations can cause early- or late-onset rheumatic monogenic diseases but also contribute to the pathogenesis of complex inflammatory and immune-mediated diseases, affect disease progression and define new clinical subtypes. Although even variants with a low variant allele fraction can be pathogenic, clonal dynamics could lead to changes over time in the proportion of mutant cells, with possible phenotypic consequences for the individual. Thus, somatic mutagenesis and clonal expansion have relevant implications in genetic testing and counselling. On the basis of both increased recognition of somatic diseases in clinical practice and improved technical and bioinformatic processes, we hypothesize that there will be an ever-expanding list of somatic mutations in various genes leading to inflammatory conditions, particularly in late-onset disease.

Human DNA polymerase ε is a source of C>T mutations at CpG dinucleotides

C-to-T transitions in CpG dinucleotides are the most prevalent mutations in human cancers and genetic diseases. These mutations have been attributed to deamination of 5-methylcytosine (5mC), an epigenetic modification found on CpGs. We recently linked CpG>TpG mutations to replication and hypothesized that errors introduced by polymerase ε (Pol ε) may represent an alternative source of mutations. Here we present a new method called polymerase error rate sequencing (PER-seq) to measure the error spectrum of DNA polymerases in isolation. We find that the most common human cancer-associated Pol ε mutant (P286R) produces an excess of CpG>TpG errors, phenocopying the mutation spectrum of tumors carrying this mutation and deficiencies in mismatch repair. Notably, we also discover that wild-type Pol ε has a sevenfold higher error rate when replicating 5mCpG compared to C in other contexts. Together, our results from PER-seq and human cancers demonstrate that replication errors are a major contributor to CpG>TpG mutagenesis in replicating cells, fundamentally changing our understanding of this important disease-causing mutational mechanism.

Somatic mutation rates scale with time not growth rate in long-lived tropical trees

The rates of appearance of new mutations play a central role in evolution. However, mutational processes in natural environments and their relationship with growth rates are largely unknown, particular in tropical ecosystems with high biodiversity. Here, we examined the somatic mutation landscapes of two tropical trees, Shorea laevis (slow-growing) and S. leprosula (fast-growing), in central Borneo, Indonesia. Using newly constructed genomes, we identified a greater number of somatic mutations in tropical trees than in temperate trees. In both species, we observed a linear increase in the number of somatic mutations with physical distance between branches. However, we found that the rate of somatic mutation accumulation per meter of growth was 3.7-fold higher in S. laevis than in S. leprosula. This difference in the somatic mutation rate was scaled with the slower growth rate of S. laevis compared to S. leprosula, resulting in a constant somatic mutation rate per year between the two species. We also found that somatic mutations are neutral within an individual, but those mutations transmitted to the next generation are subject to purifying selection. These findings suggest that somatic mutations accumulate with absolute time and older trees have a greater contribution towards generating genetic variation.

Effect of sequencing platforms on the sensitivity of chemical mutation detection using Hawk-Seq™

BACKGROUND

Error-corrected next-generation sequencing (ecNGS) technologies have enabled the direct evaluation of genome-wide mutations after exposure to mutagens. Previously, we reported an ecNGS methodology, Hawk-Seq™, and demonstrated its utility in evaluating mutagenicity. The evaluation of technical transferability is essential to further evaluate the reliability of ecNGS-based assays. However, cutting-edge sequencing platforms are continually evolving, which can affect the sensitivity of ecNGS. Therefore, the effect of differences in sequencing instruments on mutation data quality should be evaluated.

RESULTS

We assessed the performance of four sequencing platforms (HiSeq2500, NovaSeq6000, NextSeq2000, and DNBSEQ-G400) with the Hawk-Seq™ protocol for mutagenicity evaluation using DNA samples from mouse bone marrow exposed to benzo[a]pyrene (BP). The overall mutation (OM) frequencies per 106 bp in vehicle-treated samples were 0.22, 0.36, 0.46, and 0.26 for HiSeq2500, NovaSeq6000, NextSeq2000, and DNBSEQ-G400, respectively. The OM frequency of NextSeq2000 was significantly higher than that of HiSeq2500, suggesting the difference to be based on the platform. The relatively higher value in NextSeq2000 was a consequence of the G:C to C:G mutations in NextSeq2000 data (0.67 per 106 G:C bp), which was higher than the mean of the four platforms by a ca. of 0.25 per 106 G:C bp. A clear dose-dependent increase in G:C to T:A mutation frequencies was observed in all four sequencing platforms after BP exposure. The cosine similarity values of the 96-dimensional trinucleotide mutation patterns between HiSeq and the three other platforms were 0.93, 0.95, and 0.92 for NovaSeq, NextSeq, and DNBSeq, respectively. These results suggest that all platforms can provide equivalent data that reflect the characteristics of the mutagens.

CONCLUSIONS

All platforms sensitively detected mutagen-induced mutations using the Hawk-Seq™ analysis. The substitution types and frequencies of the background errors differed depending on the platform. The effects of sequencing platforms on mutagenicity evaluation should be assessed before experimentation.

Frequency and spectrum of mutations in human sperm measured using duplex sequencing correlate with trio-based de novo mutation analyses

De novo mutations (DNMs) are drivers of genetic disorders. However, the study of DNMs is hampered by technological limitations preventing accurate quantification of ultra-rare mutations. Duplex Sequencing (DS) theoretically has < 1 error/billion base-pairs (bp). To determine the DS utility to quantify and characterize DNMs, we analyzed DNA from blood and spermatozoa from six healthy, 18-year-old Swedish men using the TwinStrand DS mutagenesis panel (48 kb spanning 20 genic and intergenic loci). The mean single nucleotide variant mutation frequency (MF) was 1.2 × 10- 7 per bp in blood and 2.5 × 10- 8 per bp in sperm, with the most common base substitution being C > T. Blood MF and substitution spectrum were similar to those reported in blood cells with an orthogonal method. The sperm MF was in the same order of magnitude and had a strikingly similar spectrum to DNMs from publicly available whole genome sequencing data from human pedigrees (1.2 × 10- 8 per bp). DS revealed much larger numbers of insertions and deletions in sperm over blood, driven by an abundance of putative extra-chromosomal circular DNAs. The study indicates the strong potential of DS to characterize human DNMs to inform factors that contribute to disease susceptibility and heritable genetic risks.

Toward Systemic Lipofuscin Removal

Lipofuscin is indigestible garbage that accumulates in the autophagic vesicles and cytosol of postmitotic cells with age. Drs. Brunk and Terman postulated that lipofuscin accumulation is the main or at least a major driving factor in aging. They even posited that the evolution of memory is the reason why we get lipofuscin at all, as stable synaptic connections must be maintained over time, meaning that the somas of neurons must also remain in the same locale. In other words, they cannot dilute out their garbage over time through cell division. Mechanistically, their position certainly makes sense given that rendering a large percentage of a postmitotic cell's lysosomes useless must almost certainly negatively affect that cell and the surrounding microenvironment. It may be the case that lipofuscin accumulation is the main issue with regard to current age-related disease. Degradation in situ may be an insurmountable task currently. However, a method of systemic lipofuscin removal is discussed herein.

Somatic mutations in aging and disease

Time always leaves its mark, and our genome is no exception. Mutations in the genome of somatic cells were first hypothesized to be the cause of aging in the 1950s, shortly after the molecular structure of DNA had been described. Somatic mutation theories of aging are based on the fact that mutations in DNA as the ultimate template for all cellular functions are irreversible. However, it took until the 1990s to develop the methods to test if DNA mutations accumulate with age in different organs and tissues and estimate the severity of the problem. By now, numerous studies have documented the accumulation of somatic mutations with age in normal cells and tissues of mice, humans, and other animals, showing clock-like mutational signatures that provide information on the underlying causes of the mutations. In this review, we will first briefly discuss the recent advances in next-generation sequencing that now allow quantitative analysis of somatic mutations. Second, we will provide evidence that the mutation rate differs between cell types, with a focus on differences between germline and somatic mutation rate. Third, we will discuss somatic mutational signatures as measures of aging, environmental exposure, and activities of DNA repair processes. Fourth, we will explain the concept of clonally amplified somatic mutations, with a focus on clonal hematopoiesis. Fifth, we will briefly discuss somatic mutations in the transcriptome and in our other genome, i.e., the genome of mitochondria. We will end with a brief discussion of a possible causal contribution of somatic mutations to the aging process.

Protection of the genome and the central exome by peripheral non-coding DNA against DNA damage in health, ageing and age-related diseases

DNA in eukaryotic genomes is under constant assault from both exogenous and endogenous sources, leading to DNA damage, which is considered a major molecular driver of ageing. Fortunately, the genome and the central exome are safeguarded against these attacks by abundant peripheral non-coding DNA. Non-coding DNA codes for small non-coding RNAs that inactivate foreign nucleic acids in the cytoplasm and physically blocks these attacks in the nucleus. Damage to non-coding DNA produced during such blockage is removed in the form of extrachromosomal circular DNA (eccDNA) through nucleic pore complexes. Consequently, non-coding DNA serves as a line of defence for the exome against DNA damage. The total amount of non-coding DNA/heterochromatin declines with age, resulting in a decrease in both physical blockage and eccDNA exclusion, and thus an increase in the accumulation of DNA damage in the nucleus during ageing and in age-related diseases. Here, we summarize recent evidence supporting a protective role of non-coding DNA in healthy and pathological states and argue that DNA damage is the proximate cause of ageing and age-related genetic diseases. Strategies aimed at strengthening the protective role of non-coding DNA/heterochromatin could potentially offer better systematic protection for the dynamic genome and the exome against diverse assaults, reduce the burden of DNA damage to the exome, and thus slow ageing, counteract age-related genetic diseases and promote a healthier life for individuals.

Clonal dynamics and somatic evolution of haematopoiesis in mouse

Haematopoietic stem cells maintain blood production throughout life. While extensively characterised using the laboratory mouse, little is known about how the population is sustained and evolves with age. We isolated stem cells and progenitors from young and old mice, identifying 221,890 somatic mutations genome-wide in 1845 single cell-derived colonies, and used phylogenetic analysis to infer the ontogeny and population dynamics of the stem cell pool. Mouse stem cells and progenitors accrue ~45 somatic mutations per year, a rate only about 2-fold greater than human progenitors despite the vastly different organismal sizes and lifespans. Phylogenetic patterns reveal that stem and multipotent progenitor cell pools are both established during embryogenesis, after which they independently self-renew in parallel over life. The stem cell pool grows steadily over the mouse lifespan to approximately 70,000 cells, self-renewing about every six weeks. Aged mice did not display the profound loss of stem cell clonal diversity characteristic of human haematopoietic ageing. However, targeted sequencing revealed small, expanded clones in the context of murine ageing, which were larger and more numerous following haematological perturbations and exhibited a selection landscape similar to humans. Our data illustrate both conserved features of population dynamics of blood and distinct patterns of age-associated somatic evolution in the short-lived mouse.

Building the brain mosaic: an expanded view

The complexity of the brain is closely tied to its nature as a genetic mosaic, wherein each cell is distinguished by a unique constellation of somatic variants that contribute to functional and phenotypic diversity. Postzygotic variation arising during neurogenesis is recognized as a key contributor to brain mosaicism; however, recent advances have broadened our understanding to include sources of neural genomic diversity that develop throughout the entire lifespan, from embryogenesis through aging. Moving beyond the traditional confines of neurodevelopment, in this review, we delve into the complex mechanisms that enable various origins of brain mosaicism.

Emerging genetic technologies informing personalized medicine in Shwachman-Diamond syndrome and other inherited BMF disorders

ABSTRACT

Ribosomopathy Shwachman-Diamond syndrome (SDS) is a rare autosomal recessive inherited bone marrow failure syndrome (IBMFS) caused by mutations in the Shwachman-Bodian-Diamond syndrome gene, which is associated with an increased risk of myeloid malignancy. Tracking how hematopoietic stem cell (HSC) clonal dynamics change over time, assessing whether somatic genetic rescue mechanisms affect these dynamics, and mapping out when leukemic driver mutations are acquired is important to understand which individuals with SDS may go on to develop leukemia. In this review, we discuss how new technologies that allow researchers to map mutations at the level of single HSC clones are generating important insights into genetic rescue mechanisms and their relative risk for driving evolution to leukemia, and how these data can inform the future development of personalized medicine approaches in SDS and other IBMFSs.

Genome-scale mutational signature analysis in fixed archived tissues

Mutation spectra and mutational signatures in cancerous and non-cancerous tissues can be identified by various established techniques of massively parallel sequencing (or next-generation sequencing) including whole-exome or whole-genome sequencing, and more recently by error-corrected/duplex sequencing. One rather underexplored area has been the genome-scale analysis of mutational signatures as markers of mutagenic exposures, and their impact on cancer driver events applied to formalin-fixed or alcohol-fixed paraffin embedded archived biospecimens. This review showcases successful applications of the next-generation sequencing methodologies in archived fixed tissues, including the delineation of the specific tissue fixation-related DNA damage manifesting as artifactual signatures, distinguishable from the true signatures that arise from biological mutagenic processes. Overall, we discuss and demonstrate how next-generation sequencing techniques applied to archived fixed biospecimens can enhance our understanding of cancer causes including mutagenic effects of extrinsic cancer risk agents, and the implications for prevention efforts aimed at reducing avoidable cancer-causing exposures.

Advances in Shear Stress Stimulation of Stem Cells: A Review of the Last Three Decades

Stem cells are widely used in scientific research because of their ability to self-renew and differentiate into a variety of specialized cell types needed for body functions. However, the self-renewal and differentiation of stem cells are regulated by various stimuli, with mechanical stimulation being particularly notable due to its ability to mimic the physical environment in the body. This study systematically collected 2638 research papers published between 1994 and 2024, employing tools such as VOSviewer, CiteSpace, and GraphPad Prism to uncover research hotspots, publication trends, and collaboration networks. The results indicate a yearly increase in global research on the shear stress stimulation of stem cells, with significant contributions from the United States and China in terms of research investment and output. Future research directions include a deeper understanding of the mechanisms underlying mechanical stimulation's effects on stem cell differentiation, the development of new materials and scaffold designs to better replicate the natural cellular environment, and advancements in regenerative medicine. Despite considerable progress, challenges remain in translating basic research findings into clinical applications.

Robustness and reliability of single-cell regulatory multi-omics with deep mitochondrial mutation profiling

The detection of mitochondrial DNA (mtDNA) mutations in single cells holds considerable potential to define clonal relationships coupled with information on cell state in humans. Previous methods focused on higher heteroplasmy mutations that are limited in number and can be influenced by functional selection, introducing biases for lineage tracing. Although more challenging to detect, intermediate to low heteroplasmy mtDNA mutations are valuable due to their high diversity, abundance, and lower propensity to selection. To enhance mtDNA mutation detection and facilitate fine-scale lineage tracing, we developed the single-cell Regulatory multi-omics with Deep Mitochondrial mutation profiling (ReDeeM) approach, an integrated experimental and computational framework. Recently, some concerns have been raised about the analytical workflow in the ReDeeM framework. Specifically, it was noted that the mutations detected in a single molecule per cell are enriched on edges of mtDNA molecules, suggesting they resemble artifacts reported in other sequencing approaches. It was then proposed that all mutations found in one molecule per cell should be removed. We detail our error correction method, demonstrating that the observed edge mutations are distinct from previously reported sequencing artifacts. We further show that the proposed removal leads to massive elimination of bona fide and informative mutations. Indeed, mutations accumulating on edges impact a minority of all mutation calls (for example, in hematopoietic stem cells, the excess mutations on the edge account for only 4.3%-7.6% of the total). Recognizing the value of addressing edge mutations even after applying consensus correction, we provide an additional filtering option in the ReDeeM-R package. This approach effectively eliminates the position biases, leads to a mutational signature indistinguishable from bona fide mitochondrial mutations, and removes excess low molecule high connectedness mutations. Importantly, this option preserves the large majority of unique mutations identified by ReDeeM, maintaining the ability of ReDeeM to provide a more than 10-fold increase in variant detection compared to previous methods. Additionally, the cells remain well-connected. While there is room for further refinement in mutation calling strategies, the significant advances and biological insights provided by the ReDeeM framework are unique and remain intact. We hope that this detailed discussion and analysis enables the community to employ this approach and contribute to its further development.

High-fidelity, large-scale targeted profiling of microsatellites

Microsatellites are highly mutable sequences that can serve as markers for relationships among individuals or cells within a population. The accuracy and resolution of reconstructing these relationships depends on the fidelity of microsatellite profiling and the number of microsatellites profiled. However, current methods for targeted profiling of microsatellites incur significant "stutter" artifacts that interfere with accurate genotyping, and sequencing costs preclude whole-genome microsatellite profiling of a large number of samples. We developed a novel method for accurate and cost-effective targeted profiling of a panel of more than 150,000 microsatellites per sample, along with a computational tool for designing large-scale microsatellite panels. Our method addresses the greatest challenge for microsatellite profiling-"stutter" artifacts-with a low-temperature hybridization capture that significantly reduces these artifacts. We also developed a computational tool for accurate genotyping of the resulting microsatellite sequencing data that uses an ensemble approach integrating three microsatellite genotyping tools, which we optimize by analysis of de novo microsatellite mutations in human trios. Altogether, our suite of experimental and computational tools enables high-fidelity, large-scale profiling of microsatellites, which may find utility in diverse applications such as lineage tracing, population genetics, ecology, and forensics.

CLEMENT: genomic decomposition and reconstruction of non-tumor subclones

Genome-level clonal decomposition of a single specimen has been widely studied; however, it is mostly limited to cancer research. In this study, we developed a new algorithm CLEMENT, which conducts accurate decomposition and reconstruction of multiple subclones in genome sequencing of non-tumor (normal) samples. CLEMENT employs the Expectation-Maximization (EM) algorithm with optimization strategies specific to non-tumor subclones, including false variant call identification, non-disparate clone fuzzy clustering, and clonal allele fraction confinement. In the simulation and in vitro cell line mixture data, CLEMENT outperformed current cancer decomposition algorithms in estimating the number of clones (root-mean-square-error = 0.58-0.78 versus 1.43-3.34) and in the variant-clone membership agreement (∼85.5% versus 70.1-76.7%). Additional testing on human multi-clonal normal tissue sequencing confirmed the accurate identification of subclones that originated from different cell types. Clone-level analysis, including mutational burden and signatures, provided a new understanding of normal-tissue composition. We expect that CLEMENT will serve as a crucial tool in the currently emerging field of non-tumor genome analysis.

Genetic variation across and within individuals

Germline variation and somatic mutation are intricately connected and together shape human traits and disease risks. Germline variants are present from conception, but they vary between individuals and accumulate over generations. By contrast, somatic mutations accumulate throughout life in a mosaic manner within an individual due to intrinsic and extrinsic sources of mutations and selection pressures acting on cells. Recent advancements, such as improved detection methods and increased resources for association studies, have drastically expanded our ability to investigate germline and somatic genetic variation and compare underlying mutational processes. A better understanding of the similarities and differences in the types, rates and patterns of germline and somatic variants, as well as their interplay, will help elucidate the mechanisms underlying their distinct yet interlinked roles in human health and biology.

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.

Chasing Sequencing Perfection: Marching Toward Higher Accuracy and Lower Costs

Next-generation sequencing (NGS), represented by Illumina platforms, has been an essential cornerstone of basic and applied research. However, the sequencing error rate of 1 per 1000 bp (10-3) represents a serious hurdle for research areas focusing on rare mutations, such as somatic mosaicism or microbe heterogeneity. By examining the high-fidelity sequencing methods developed in the past decade, we summarized three major factors underlying errors and the corresponding 12 strategies mitigating these errors. We then proposed a novel framework to classify 11 preexisting representative methods according to the corresponding combinatory strategies and identified three trends that emerged during methodological developments. We further extended this analysis to eight long-read sequencing methods, emphasizing error reduction strategies. Finally, we suggest two promising future directions that could achieve comparable or even higher accuracy with lower costs in both NGS and long-read sequencing.

Harnessing Single-Cell Technologies in the Search for New Therapies for Diamond-Blackfan Anemia Syndrome

The emergence of multiomic single-cell technologies over the last decade has led to improved insights into both normal hematopoiesis and its perturbation in a variety of hematological disorders. Diamond-Blackfan anemia (DBA) syndrome is one such disorder where single-cell assays have helped to delineate the cellular and molecular defects underlying the disease. DBA is caused by heterozygous loss-of-function germline variants in genes encoding ribosomal proteins (RPs). Despite the widespread role of ribosomes in hematopoiesis, the most frequent and severe cytopenia in DBA is anemia. In this review we discussed how single-cell studies, including clonogenic cell culture assays, fluorescence-activated cell sorting (FACS) and single-cell RNA sequencing (scRNA-seq), have led to insights into the pathogenesis of DBA. The main therapies are regular blood transfusions, glucocorticoids, or hematopoietic stem cell transplantation (HSCT) but all are associated with significant morbidity and mortality. We will therefore outline how single-cell studies can inform new therapies for DBA. Furthermore, we discussed how DBA serves as a useful model for understanding normal erythropoiesis in terms of its cellular hierarchy, molecular regulation during homeostasis, and response to "stress."

Transmissible cancers, the genomes that do not melt down

Evolutionary theory predicts that the accumulation of deleterious mutations in asexually reproducing organisms should lead to genomic decay. Clonally reproducing cell lines, i.e., transmissible cancers, when cells are transmitted as allografts/xenografts, break these rules and survive for centuries and millennia. The currently known 11 transmissible cancer lineages occur in dogs (canine venereal tumour disease), in Tasmanian devils (devil facial tumor diseases, DFT1 and DFT2), and in bivalves (bivalve transmissible neoplasia). Despite the mutation loads of these cell lines being much higher than observed in human cancers, they have not been eliminated in space and time. Here, we provide potential explanations for how these fascinating cell lines may have overcome the fitness decline due to the progressive accumulation of deleterious mutations and propose that the high mutation load may carry an indirect positive fitness outcome. We offer ideas on how these host-pathogen systems could be used to answer outstanding questions in evolutionary biology. The recent studies on the evolution of these clonal pathogens reveal key mechanistic insight into transmissible cancer genomes, information that is essential for future studies investigating how these contagious cancer cell lines can repeatedly evade immune recognition, evolve, and survive in the landscape of highly diverse hosts.

Clonal expansion in normal tissues

Cancer originates from a single ancestral cell that acquires a driver mutation, which confers a growth or survival advantage, followed by the acquisition of additional driver mutations by descendant cells. Recently, it has become evident that somatic cell mutations accumulate in normal tissues with aging and exposure to environmental factors, such as alcohol, smoking, and UV rays, increases the mutation rate. Clones harboring driver mutations expand with age, leading to tissue remodeling. Lineage analysis of myeloproliferative neoplasms and der(1;16)-positive breast cancer revealed that driver mutations were acquired early in our lives and that the development of cancer takes decades, unveiling the previously unknown early process of cancer development. Evidence that clonal hematopoiesis affects various diseases, including nonneoplastic diseases, highlights the potential role of the identification and functional analysis of mutated clones in unraveling unknown pathologies. In this review, we summarize the recent updates on clonal expansion in normal tissues and the natural history of cancer revealed through lineage analysis of noncancerous and cancerous tissues.

A somatic genetic clock for clonal species

Age and longevity are key parameters for demography and life-history evolution of organisms. In clonal species, a widespread life history among animals, plants, macroalgae and fungi, the sexually produced offspring (genet) grows indeterminately by producing iterative modules, or ramets, and so obscure their age. Here we present a novel molecular clock based on the accumulation of fixed somatic genetic variation that segregates among ramets. Using a stochastic model, we demonstrate that the accumulation of fixed somatic genetic variation will approach linearity after a lag phase, and is determined by the mitotic mutation rate, without direct dependence on asexual generation time. The lag phase decreased with lower stem cell population size, number of founder cells for the formation of new modules, and the ratio of symmetric versus asymmetric cell divisions. We calibrated the somatic genetic clock on cultivated eelgrass Zostera marina genets (4 and 17 years respectively). In a global data set of 20 eelgrass populations, genet ages were up to 1,403 years. The somatic genetic clock is applicable to any multicellular clonal species where the number of founder cells is small, opening novel research avenues to study longevity and, hence, demography and population dynamics of clonal 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.

Beyond CAG Repeats: The Multifaceted Role of Genetics in Huntington Disease

Huntington disease (HD) is a dominantly inherited neurodegenerative disorder caused by a CAG expansion on the huntingtin (HTT) gene and is characterized by progressive motor, cognitive, and neuropsychiatric decline. Recently, new genetic factors besides CAG repeats have been implicated in the disease pathogenesis. Most genetic modifiers are involved in DNA repair pathways and, as the cause of the loss of CAA interruption in the HTT gene, they exert their main influence through somatic expansion. However, this mechanism might not be the only driver of HD pathogenesis, and future studies are warranted in this field. The aim of the present review is to dissect the many faces of genetics in HD pathogenesis, from cis- and trans-acting genetic modifiers to RNA toxicity, mitochondrial DNA mutations, and epigenetics factors. Exploring genetic modifiers of HD onset and progression appears crucial to elucidate not only disease pathogenesis, but also to improve disease prediction and prevention, develop biomarkers of disease progression and response to therapies, and recognize new therapeutic opportunities. Since the same genetic mechanisms are also described in other repeat expansion diseases, their implications might encompass the whole spectrum of these disorders.

Exploring the potential of Oxford Nanopore Technologies sequencing for Mycobacterium tuberculosis sequencing: An assessment of R10 flowcells and V14 chemistry

Oxford Nanopore Technologies (ONT) sequencing is a promising technology. We assessed the performance of the new ONT R10 flowcells and V14 rapid sequencing chemistry for Mtb whole genome sequencing of Mycobacterium tuberculosis (Mtb) DNA extracted from clinical primary liquid cultures (CPLCs). Using the recommended protocols for MinION Mk1C, R10.4.1 MinION flowcells, and the ONT Rapid Sequencing Kit V14 on six CPLC samples, we obtained a pooled library yield of 10.9 ng/μl, generated 1.94 Gb of sequenced bases and 214k reads after 48h in a first sequencing run. Only half (49%) of all generated reads met the Phred Quality score threshold (>8). To assess if the low data output and sequence quality were due to impurities present in DNA extracted directly from CPLCs, we added a pre-library preparation bead-clean-up step and included purified DNA obtained from an Mtb subculture as a control sample in a second sequencing run. The library yield for DNA extracted from four CPLCs and one Mtb subculture (control) was similar (10.0 ng/μl), 2.38 Gb of bases and 822k reads were produced. The quality was slightly better with 66% of the produced reads having a Phred Quality >8. A third run of DNA from six CPLCs with bead clean-up pre-processing produced a low library yield (±1 Gb of bases, 166k reads) of low quality (51% of reads with a Phred Quality score >8). A median depth of coverage above 10× was only achieved for five of 17 (29%) sequenced libraries. Compared to Illumina WGS of the same samples, accurate lineage predictions and full drug resistance profiles from the generated ONT data could not be determined by TBProfiler. Further optimization of the V14 ONT rapid sequencing chemistry and library preparation protocol is needed for clinical Mtb WGS applications.

Focal cortical dysplasia II caused by brain somatic mutation of IRS-1 is associated with ERK signaling pathway activation

Somatic mutations have been identified in 10% to 63% of focal cortical dysplasia type II samples, primarily linked to the mTOR pathway. When the causative genetic mutations are not identified, this opens the possibility of discovering new pathogenic genes or pathways that could be contributing to the condition. In our previous study, we identified a novel candidate pathogenic somatic variant of IRS-1 c.1791dupG in the brain tissue of a child with focal cortical dysplasia type II. This study further explored the variant's role in causing type II focal cortical dysplasia through in vitro overexpression in 293T and SH-SY5Y cells and in vivo evaluation via in utero electroporation in fetal brains, assessing effects on neuronal migration, morphology, and network integrity. It was found that the mutant IRS-1 variant led to hyperactivity of p-ERK, increased cell volume, and was predominantly associated with the MAPK signaling pathway. In vivo, the IRS-1 c.1791dupG variant induced abnormal neuron migration, cytomegaly, and network hyperexcitability. Notably, the ERK inhibitor GDC-0994, rather than the mTOR inhibitor rapamycin, effectively rescued the neuronal defects. This study directly highlighted the ERK signaling pathway's role in the pathogenesis of focal cortical dysplasia II and provided a new therapeutic target for cases of focal cortical dysplasia II that are not treatable by rapamycin analogs.

Severity of effect considerations regarding the use of mutation as a toxicological endpoint for risk assessment: A report from the 8th International Workshop on Genotoxicity Testing (IWGT)

Exposure levels without appreciable human health risk may be determined by dividing a point of departure on a dose-response curve (e.g., benchmark dose) by a composite adjustment factor (AF). An "effect severity" AF (ESAF) is employed in some regulatory contexts. An ESAF of 10 may be incorporated in the derivation of a health-based guidance value (HBGV) when a "severe" toxicological endpoint, such as teratogenicity, irreversible reproductive effects, neurotoxicity, or cancer was observed in the reference study. Although mutation data have been used historically for hazard identification, this endpoint is suitable for quantitative dose-response modeling and risk assessment. As part of the 8th International Workshops on Genotoxicity Testing, a sub-group of the Quantitative Analysis Work Group (WG) explored how the concept of effect severity could be applied to mutation. To approach this question, the WG reviewed the prevailing regulatory guidance on how an ESAF is incorporated into risk assessments, evaluated current knowledge of associations between germline or somatic mutation and severe disease risk, and mined available data on the fraction of human germline mutations expected to cause severe disease. Based on this review and given that mutations are irreversible and some cause severe human disease, in regulatory settings where an ESAF is used, a majority of the WG recommends applying an ESAF value between 2 and 10 when deriving a HBGV from mutation data. This recommendation may need to be revisited in the future if direct measurement of disease-causing mutations by error-corrected next generation sequencing clarifies selection of ESAF values.

Genomic Mosaicism of the Brain: Origin, Impact, and Utility

Genomic mosaicism describes the phenomenon where some but not all cells within a tissue harbor unique genetic mutations. Traditionally, research focused on the impact of genomic mosaicism on clinical phenotype-motivated by its involvement in cancers and overgrowth syndromes. More recently, we increasingly shifted towards the plethora of neutral mosaic variants that can act as recorders of cellular lineage and environmental exposures. Here, we summarize the current state of the field of genomic mosaicism research with a special emphasis on our current understanding of this phenomenon in brain development and homeostasis. Although the field of genomic mosaicism has a rich history, technological advances in the last decade have changed our approaches and greatly improved our knowledge. We will provide current definitions and an overview of contemporary detection approaches for genomic mosaicism. Finally, we will discuss the impact and utility of genomic mosaicism.

The clock-like accumulation of germline and somatic mutations can arise from the interplay of DNA damage and repair

The rates at which mutations accumulate across human cell types vary. 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 thought to be 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 2 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 postmitotic 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.

Sensitive circulating tumor DNA-based residual disease detection in epithelial ovarian cancer

Epithelial ovarian cancer (EOC) is one of the leading causes of cancer-related death in women worldwide, and is characterized by a high rate of recurrence after surgery and chemotherapy. We sought to implement a circulating tumor DNA (ctDNA)-based blood test for more accurate post-operative surveillance of this disease. We analyzed 264 plasma samples collected between June 2016 and September 2021 from 63 EOC patients using tumor-guided plasma cell-free DNA analysis to detect residual disease after treatment. Assay specificity was verified using cross-patient analysis of 1,195 control samples. ctDNA was detected in 51 of 55 (93%) samples at diagnosis, and 18 of 18 (100%) samples at progression. Positive ctDNA in the last on-treatment sample was associated with rapid progression (median 1.02 versus 3.38 yr, HR = 5.63, P < 0.001) and reduced overall survival (median 2.31 versus NR yr, HR = 8.22, P < 0.001) in patients with high-grade serous cancer. In the case of 12 patients, ctDNA assays detected progression earlier than standard surveillance, with a median lead time of 5.9 mo. To approach the physical limits of ctDNA detection, five patients were analyzed using ultra-sensitive assays interrogating 479-1,856 tumor mutations, capable of tracking ctDNA fractions down to 0.0004%. Our results demonstrate that ctDNA assays achieve high sensitivity and specificity in detecting post-operative residual disease in EOC.

DNA mismatch and damage patterns revealed by single-molecule sequencing

Mutations accumulate in the genome of every cell of the body throughout life, causing cancer and other diseases1,2. Most mutations begin as nucleotide mismatches or damage in one of the two strands of the DNA before becoming double-strand mutations if unrepaired or misrepaired3,4. However, current DNA-sequencing technologies cannot accurately resolve these initial single-strand events. Here we develop a single-molecule, long-read sequencing method (Hairpin Duplex Enhanced Fidelity sequencing (HiDEF-seq)) that achieves single-molecule fidelity for base substitutions when present in either one or both DNA strands. HiDEF-seq also detects cytosine deamination-a common type of DNA damage-with single-molecule fidelity. We profiled 134 samples from diverse tissues, including from individuals with cancer predisposition syndromes, and derive from them single-strand mismatch and damage signatures. We find correspondences between these single-strand signatures and known double-strand mutational signatures, which resolves the identity of the initiating lesions. Tumours deficient in both mismatch repair and replicative polymerase proofreading show distinct single-strand mismatch patterns compared to samples that are deficient in only polymerase proofreading. We also define a single-strand damage signature for APOBEC3A. In the mitochondrial genome, our findings support a mutagenic mechanism occurring primarily during replication. As double-strand DNA mutations are only the end point of the mutation process, our approach to detect the initiating single-strand events at single-molecule resolution will enable studies of how mutations arise in a variety of contexts, especially in cancer and ageing.

Mutations causing premature termination codons discriminate and generate cellular and clinical variability in HHT

ABSTRACT

For monogenic diseases caused by pathogenic loss-of-function DNA variants, attention focuses on dysregulated gene-specific pathways, usually considering molecular subtypes together within causal genes. To better understand phenotypic variability in hereditary hemorrhagic telangiectasia (HHT), we subcategorized pathogenic DNA variants in ENG/endoglin, ACVRL1/ALK1, and SMAD4 if they generated premature termination codons (PTCs) subject to nonsense-mediated decay. In 3 patient cohorts, a PTC-based classification system explained some previously puzzling hemorrhage variability. In blood outgrowth endothelial cells (BOECs) derived from patients with ACVRL1+/PTC, ENG+/PTC, and SMAD4+/PTC genotypes, PTC-containing RNA transcripts persisted at low levels (8%-23% expected, varying between replicate cultures); genes differentially expressed to Bonferroni P < .05 in HHT+/PTC BOECs clustered significantly only to generic protein terms (isopeptide-bond/ubiquitin-like conjugation) and pulse-chase experiments detected subtle protein maturation differences but no evidence for PTC-truncated protein. BOECs displaying highest PTC persistence were discriminated in unsupervised hierarchical clustering of near-invariant housekeeper genes, with patterns compatible with higher cellular stress in BOECs with >11% PTC persistence. To test directionality, we used a HeLa reporter system to detect induction of activating transcription factor 4 (ATF4), which controls expression of stress-adaptive genes, and showed that ENG Q436X but not ENG R93X directly induced ATF4. AlphaFold accurately modeled relevant ENG domains, with AlphaMissense suggesting that readthrough substitutions would be benign for ENG R93X and other less rare ENG nonsense variants but more damaging for Q436X. We conclude that PTCs should be distinguished from other loss-of-function variants, PTC transcript levels increase in stressed cells, and readthrough proteins and mechanisms provide promising research avenues.

Clonal evolution of hematopoietic stem cells after cancer chemotherapy

Normal hematopoietic stem and progenitor cells (HSPCs) inherently accumulate somatic mutations and lose clonal diversity with age, processes implicated in the development of myeloid malignancies 1 . The impact of exogenous stressors, such as cancer chemotherapies, on the genomic integrity and clonal dynamics of normal HSPCs is not well defined. We conducted whole-genome sequencing on 1,032 single-cell-derived HSPC colonies from 10 patients with multiple myeloma (MM), who had undergone various chemotherapy regimens. Our findings reveal that melphalan treatment distinctly increases mutational burden with a unique mutation signature, whereas other MM chemotherapies do not significantly affect the normal mutation rate of HSPCs. Among these therapy-induced mutations were several oncogenic drivers such as TET2 and PPM1D . Phylogenetic analysis showed a clonal architecture in post-treatment HSPCs characterized by extensive convergent evolution of mutations in genes such as TP53 and PPM1D . Consequently, the clonal diversity and structure of post-treatment HSPCs mirror those observed in normal elderly individuals, suggesting an accelerated clonal aging due to chemotherapy. Furthermore, analysis of matched therapy-related myeloid neoplasm (t-MN) samples, which occurred 1-8 years later, enabled us to trace the clonal origin of t-MNs to a single HSPC clone among a group of clones with competing malignant potential, indicating the critical role of secondary mutations in dictating clonal dominance and malignant transformation. Our findings suggest that cancer chemotherapy promotes an oligoclonal architecture with multiple HSPC clones possessing competing leukemic potentials, setting the stage for the selective emergence of a singular clone that evolves into t-MNs after acquiring secondary mutations. These results underscore the importance of further systematic research to elucidate the long-term hematological consequences of cancer chemotherapy.

[Sperm Mosaic Variants and Their Influence on the Offspring]

Genomic mosaicism arising from mosaic variants is a phenomenon that describes the presence of a cell or cell populations with different genome compositions from the germline cells of an individual. It comprises all types of genetic variants. A large proportion of childhood genetic disorders are defined as being de novo, meaning that the disease-causing mutations are only detected in the proband, not in any of the parents. Population studies show that 80% of the de novo mutations arise from the paternal haplotype, that is, from paternal sperm mosaicism. This review provides a summary of the types and detection strategies of sperm mosaicism. In addition, it provides discussions on how recent studies demonstrated that genomic mosaic mutations in parents, especially those in the paternal sperms, could be inherited by the offspring and cause childhood disorders. According to the previous findings of the author's research team, sperm mosaicism derived from early embryogenesis and primordial germ cell stages can explain 5% to 20% of the de novo mutations related to clinical phenotypes and can serve as an important predictor of both rare and complex disorders. Sperm mosaicism shows great potential for clinical genetic diagnosis and consultations. Based on the published literature, the author suggests that, large-scale screening for de novo sperm mosaic mutations and population-based genetic screening should be conducted in future studies, which will greatly enhance the risk assessment in the offspring and effectively improve the genetic health at the population level. Implementation of direct sperm detection for de novo mutations will significantly increase the efficiency of the stratification of patient cohorts and improve recurrence risk assessment for future births. Future research in the field should be focused on the impact of environmental and lifestyle factors on the health of the offspring through sperms and their modeling of mutation signatures. In addition, targeted in vitro modeling of sperm mutations will also be a promising direction.

An improved epigenetic counter to track mitotic age in normal and precancerous tissues

The cumulative number of stem cell divisions in a tissue, known as mitotic age, is thought to be a major determinant of cancer-risk. Somatic mutational and DNA methylation (DNAm) clocks are promising tools to molecularly track mitotic age, yet their relationship is underexplored and their potential for cancer risk prediction in normal tissues remains to be demonstrated. Here we build and validate an improved pan-tissue DNAm counter of total mitotic age called stemTOC. We demonstrate that stemTOC's mitotic age proxy increases with the tumor cell-of-origin fraction in each of 15 cancer-types, in precancerous lesions, and in normal tissues exposed to major cancer risk factors. Extensive benchmarking against 6 other mitotic counters shows that stemTOC compares favorably, specially in the preinvasive and normal-tissue contexts. By cross-correlating stemTOC to two clock-like somatic mutational signatures, we confirm the mitotic-like nature of only one of these. Our data points towards DNAm as a promising molecular substrate for detecting mitotic-age increases in normal tissues and precancerous lesions, and hence for developing cancer-risk prediction strategies.

Mapping recurrent mosaic copy number variation in human neurons

When somatic cells acquire complex karyotypes, they often are removed by the immune system. Mutant somatic cells that evade immune surveillance can lead to cancer. Neurons with complex karyotypes arise during neurotypical brain development, but neurons are almost never the origin of brain cancers. Instead, somatic mutations in neurons can bring about neurodevelopmental disorders, and contribute to the polygenic landscape of neuropsychiatric and neurodegenerative disease. A subset of human neurons harbors idiosyncratic copy number variants (CNVs, "CNV neurons"), but previous analyses of CNV neurons are limited by relatively small sample sizes. Here, we develop an allele-based validation approach, SCOVAL, to corroborate or reject read-depth based CNV calls in single human neurons. We apply this approach to 2,125 frontal cortical neurons from a neurotypical human brain. SCOVAL identifies 226 CNV neurons, which include a subclass of 65 CNV neurons with highly aberrant karyotypes containing whole or substantial losses on multiple chromosomes. Moreover, we find that CNV location appears to be nonrandom. Recurrent regions of neuronal genome rearrangement contain fewer, but longer, genes.

The broken Alzheimer's disease genome

The complex pathobiology of late-onset Alzheimer's disease (AD) poses significant challenges to therapeutic and preventative interventions. Despite these difficulties, genomics and related disciplines are allowing fundamental mechanistic insights to emerge with clarity, particularly with the introduction of high-resolution sequencing technologies. After all, the disrupted processes at the interface between DNA and gene expression, which we call the broken AD genome, offer detailed quantitative evidence unrestrained by preconceived notions about the disease. In addition to highlighting biological pathways beyond the classical pathology hallmarks, these advances have revitalized drug discovery efforts and are driving improvements in clinical tools. We review genetic, epigenomic, and gene expression findings related to AD pathogenesis and explore how their integration enables a better understanding of the multicellular imbalances contributing to this heterogeneous condition. The frontiers opening on the back of these research milestones promise a future of AD care that is both more personalized and predictive.

Extended Enrichment for Ultrasensitive Detection of Low-Frequency Mutations by Long Blocker Displacement Amplification

Detecting low-frequency DNA mutations hotspots cluster is critical for cancer diagnosis but remains challenging. Quantitative PCR (qPCR) is constrained by sensitivity, and allele-specific PCR is restricted by throughput. Here we develop a long blocker displacement amplification (LBDA) coupled with qPCR for ultrasensitive and multiplexed variants detection. By designing long blocker oligos to perfectly match wildtype sequences while mispairing with mutants, long blockers enable 14-44 nt enrichment regions which is 2-fold longer than normal BDA in the experiments. For wild template with a specific nucleotide, LBDA can detect different mutation types down to 0.5 % variant allele frequency (VAF) in one reaction, with median enrichment fold of 1,000 on 21 mutant DNA templates compared to the wild type. We applied LBDA-qPCR to detect KRAS and NRAS mutation hotspots, utilizing a single plex assay capable of covering 81 mutations and tested in synthetic templates and colorectal cancer tissue samples. Moreover, the mutation types were verified through Sanger sequencing, demonstrating concordance with results obtained from next generation sequencing. Overall, LBDA-qPCR provides a simple yet ultrasensitive approach for multiplexed detection of low VAF mutations hotspots, presenting a powerful tool for cancer diagnosis and monitoring.

Luminal breast epithelial cells from wildtype and BRCA mutation carriers harbor copy number alterations commonly associated with breast cancer

Cancer-associated mutations have been documented in normal tissues, but the prevalence and nature of somatic copy number alterations and their role in tumor initiation and evolution is not well understood. Here, using single cell DNA sequencing, we describe the landscape of CNAs in >42,000 breast epithelial cells from women with normal or high risk of developing breast cancer. Accumulation of individual cells with one or two of a specific subset of CNAs (e.g. 1q gain and 16q, 22q, 7q, and 10q loss) is detectable in almost all breast tissues and, in those from BRCA1 or BRCA2 mutations carriers, occurs prior to loss of heterozygosity (LOH) of the wildtype alleles. These CNAs, which are among the most common associated with ductal carcinoma in situ (DCIS) and malignant breast tumors, are enriched almost exclusively in luminal cells not basal myoepithelial cells. Allele-specific analysis of the enriched CNAs reveals that each allele was independently altered, demonstrating convergent evolution of these CNAs in an individual breast. Tissues from BRCA1 or BRCA2 mutation carriers contain a small percentage of cells with extreme aneuploidy, featuring loss of TP53 , LOH of BRCA1 or BRCA2 , and multiple breast cancer-associated CNAs in addition to one or more of the common CNAs in 1q, 10q or 16q. Notably, cells with intermediate levels of CNAs are not detected, arguing against a stepwise gradual accumulation of CNAs. Overall, our findings demonstrate that chromosomal alterations in normal breast epithelium partially mirror those of established cancer genomes and are chromosome- and cell lineage-specific.

Telomere length and micronuclei trajectories in APP/PS1 mouse model of Alzheimer's disease: Correlating with cognitive impairment and brain amyloidosis in a sexually dimorphic manner

Although studies have demonstrated that genome instability is accumulated in patients with Alzheimer's disease (AD), the specific types of genome instability linked to AD pathogenesis remain poorly understood. Here, we report the first characterization of the age- and sex-related trajectories of telomere length (TL) and micronuclei in APP/PS1 mice model and wild-type (WT) controls (C57BL/6). TL was measured in brain (prefrontal cortex, cerebellum, pituitary gland, and hippocampus), colon and skin, and MN was measured in bone marrow in 6- to 14-month-old mice. Variation in TL was attributable to tissue type, age, genotype and, to a lesser extent, sex. Compared to WT, APP/PS1 had a significantly shorter baseline TL across all examined tissues. TL was inversely associated with age in both genotypes and TL shortening was accelerated in brain of APP/PS1. Age-related increase of micronuclei was observed in both genotypes but was accelerated in APP/PS1. We integrated TL and micronuclei data with data on cognition performance and brain amyloidosis. TL and micronuclei were linearly correlated with cognition performance or Aβ40 and Aβ42 levels in both genotypes but to a greater extent in APP/PS1. These associations in APP/PS1 mice were dominantly driven by females. Together, our findings provide foundational knowledge to infer the TL and micronuclei trajectories in APP/PS1 mice during disease progression, and strongly support that TL attrition and micronucleation are tightly associated with AD pathogenesis in a female-biased manner.

De novo detection of somatic mutations in high-throughput single-cell profiling data sets

Characterization of somatic mutations at single-cell resolution is essential to study cancer evolution, clonal mosaicism and cell plasticity. Here, we describe SComatic, an algorithm designed for the detection of somatic mutations in single-cell transcriptomic and ATAC-seq (assay for transposase-accessible chromatin sequence) data sets directly without requiring matched bulk or single-cell DNA sequencing data. SComatic distinguishes somatic mutations from polymorphisms, RNA-editing events and artefacts using filters and statistical tests parameterized on non-neoplastic samples. Using >2.6 million single cells from 688 single-cell RNA-seq (scRNA-seq) and single-cell ATAC-seq (scATAC-seq) data sets spanning cancer and non-neoplastic samples, we show that SComatic detects mutations in single cells accurately, even in differentiated cells from polyclonal tissues that are not amenable to mutation detection using existing methods. Validated against matched genome sequencing and scRNA-seq data, SComatic achieves F1 scores between 0.6 and 0.7 across diverse data sets, in comparison to 0.2-0.4 for the second-best performing method. In summary, SComatic permits de novo mutational signature analysis, and the study of clonal heterogeneity and mutational burdens at single-cell resolution.

Geographic variation of mutagenic exposures in kidney cancer genomes

International differences in the incidence of many cancer types indicate the existence of carcinogen exposures that have not yet been identified by conventional epidemiology make a substantial contribution to cancer burden1. In clear cell renal cell carcinoma, obesity, hypertension and tobacco smoking are risk factors, but they do not explain the geographical variation in its incidence2. Underlying causes can be inferred by sequencing the genomes of cancers from populations with different incidence rates and detecting differences in patterns of somatic mutations. Here we sequenced 962 clear cell renal cell carcinomas from 11 countries with varying incidence. The somatic mutation profiles differed between countries. In Romania, Serbia and Thailand, mutational signatures characteristic of aristolochic acid compounds were present in most cases, but these were rare elsewhere. In Japan, a mutational signature of unknown cause was found in more than 70% of cases but in less than 2% elsewhere. A further mutational signature of unknown cause was ubiquitous but exhibited higher mutation loads in countries with higher incidence rates of kidney cancer. Known signatures of tobacco smoking correlated with tobacco consumption, but no signature was associated with obesity or hypertension, suggesting that non-mutagenic mechanisms of action underlie these risk factors. The results of this study indicate the existence of multiple, geographically variable, mutagenic exposures that potentially affect tens of millions of people and illustrate the opportunities for new insights into cancer causation through large-scale global cancer genomics.

Mitochondrial haplotype and mito-nuclear matching drive somatic mutation and selection throughout ageing

Mitochondrial genomes co-evolve with the nuclear genome over evolutionary timescales and are shaped by selection in the female germline. Here we investigate how mismatching between nuclear and mitochondrial ancestry impacts the somatic evolution of the mitochondrial genome in different tissues throughout ageing. We used ultrasensitive duplex sequencing to profile ~2.5 million mitochondrial genomes across five mitochondrial haplotypes and three tissues in young and aged mice, cataloguing ~1.2 million mitochondrial somatic and ultralow-frequency inherited mutations, of which 81,097 are unique. We identify haplotype-specific mutational patterns and several mutational hotspots, including at the light strand origin of replication, which consistently exhibits the highest mutation frequency. We show that rodents exhibit a distinct mitochondrial somatic mutational spectrum compared with primates with a surfeit of reactive oxygen species-associated G > T/C > A mutations, and that somatic mutations in protein-coding genes exhibit signatures of negative selection. Lastly, we identify an extensive enrichment in somatic reversion mutations that 're-align' mito-nuclear ancestry within an organism's lifespan. Together, our findings demonstrate that mitochondrial genomes are a dynamically evolving subcellular population shaped by somatic mutation and selection throughout organismal lifetimes.

Droplet based whole genome amplification for sequencing minute amounts of purified Mycobacterium tuberculosis DNA

Implementation of whole genome sequencing (WGS) for patient care is hindered by limited Mycobacterium tuberculosis (Mtb) in clinical specimens and slow Mtb growth. We evaluated droplet multiple displacement amplification (dMDA) for amplification of minute amounts of Mtb DNA to enable WGS as an alternative to other Mtb enrichment methods. Purified genomic Mtb-DNA (0.1, 0.5, 1, and 5 pg) was encapsulated and amplified using the Samplix Xdrop-instrument and sequenced alongside a control sample using standard Illumina protocols followed by MAGMA-analysis. The control and 5 pg input dMDA samples underwent nanopore sequencing followed by Nanoseq and TB-profiler analysis. dMDA generated 105-2400 ng DNA from the 0.1-5 pg input DNA, respectively. Followed by Illumina WGS, dMDA raised mean sequencing depth from 7 × for 0.1 pg input DNA to ≥ 60 × for 5 pg input and the control sample. Bioinformatic analysis revealed a high number of false positive and false negative variants when amplifying ≤ 0.5 pg input DNA. Nanopore sequencing of the 5 pg dMDA sample presented excellent coverage depth, breadth, and accurate strain characterization, albeit elevated false positive and false negative variants compared to Illumina-sequenced dMDA sample with identical Mtb DNA input. dMDA coupled with Illumina WGS for samples with ≥ 5 pg purified Mtb DNA, equating to approximately 1000 copies of the Mtb genome, offers precision for drug resistance, phylogeny, and transmission insights.

Contrasting somatic mutation patterns in aging human neurons and oligodendrocytes

Characterizing somatic mutations in the brain is important for disentangling the complex mechanisms of aging, yet little is known about mutational patterns in different brain cell types. Here, we performed whole-genome sequencing (WGS) of 86 single oligodendrocytes, 20 mixed glia, and 56 single neurons from neurotypical individuals spanning 0.4-104 years of age and identified >92,000 somatic single-nucleotide variants (sSNVs) and small insertions/deletions (indels). Although both cell types accumulate somatic mutations linearly with age, oligodendrocytes accumulated sSNVs 81% faster than neurons and indels 28% slower than neurons. Correlation of mutations with single-nucleus RNA profiles and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed across the genome similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These stark differences suggest an assortment of active mutagenic processes in oligodendrocytes and neurons.

DNA Damage and Parkinson's Disease

The etiology underlying most sporadic Parkinson's' disease (PD) cases is unknown. Environmental exposures have been suggested as putative causes of the disease. In cell models and in animal studies, certain chemicals can destroy dopaminergic neurons. However, the mechanisms of how these chemicals cause the death of neurons is not understood. Several of these agents are mitochondrial toxins that inhibit the mitochondrial complex I of the electron transport chain. Familial PD genes also encode proteins with important functions in mitochondria. Mitochondrial dysfunction of the respiratory chain, in combination with the presence of redox active dopamine molecules in these cells, will lead to the accumulation of reactive oxygen species (ROS) in dopaminergic neurons. Here, I propose a mechanism regarding how ROS may lead to cell killing with a specificity for neurons. One rarely considered hypothesis is that ROS produced by defective mitochondria will lead to the formation of oxidative DNA damage in nuclear DNA. Many genes that encode proteins with neuron-specific functions are extraordinary long, ranging in size from several hundred kilobases to well over a megabase. It is predictable that such long genes will contain large numbers of damaged DNA bases, for example in the form of 8-oxoguanine (8-oxoG), which is a major DNA damage type produced by ROS. These DNA lesions will slow down or stall the progression of RNA polymerase II, which is a term referred to as transcription stress. Furthermore, ROS-induced DNA damage may cause mutations, even in postmitotic cells such as neurons. I propose that the impaired transcription and mutagenesis of long, neuron-specific genes will lead to a loss of neuronal integrity, eventually leading to the death of these cells during a human lifetime.