Somatic mutation landscapes at single-molecule resolution

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.

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.

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 and 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.

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.

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.

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.

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.

Impact of CRISPR/HDR editing versus lentiviral transduction on long-term engraftment and clonal dynamics of HSPCs in rhesus macaques

For precise genome editing via CRISPR/homology-directed repair (HDR), effective and safe editing of long-term engrafting hematopoietic stem cells (LT-HSCs) is required. The impact of HDR on true LT-HSC clonal dynamics in a relevant large animal model has not been studied. To track the output and clonality of HDR-edited cells and to provide a comparison to lentivirally transduced HSCs in vivo, we developed a competitive rhesus macaque (RM) autologous transplantation model, co-infusing HSCs transduced with a barcoded GFP-expressing lentiviral vector (LV) and HDR edited at the CD33 locus. CRISPR/HDR-edited cells showed a two-log decrease by 2 months following transplantation, with little improvement via p53 inhibition, in comparison to minimal loss of LV-transduced cells long term. HDR long-term clonality was oligoclonal in contrast to highly polyclonal LV-transduced HSCs. These results suggest marked clinically relevant differences in the impact of current genetic modification approaches on HSCs.

Somatic mutations in human ageing: New insights from DNA sequencing and inherited mutations

The accumulation of somatic mutations is a driver of cancer and has long been associated with ageing. Due to limitations in quantifying mutation burden with age in non-cancerous tissues, the impact of somatic mutations in other ageing phenotypes is unclear. Recent advances in DNA sequencing technologies have allowed the large-scale quantification of somatic mutations in ageing tissues. These studies have revealed a gradual accumulation of mutations in normal tissues with age as well as a substantial clonal expansion driven mostly by cancer-related mutations. Nevertheless, it is difficult to envision how the burden and stochastic nature of age-related somatic mutations identified so far can explain most ageing phenotypes that develop gradually. Studies across species have also found that longer-lived species have lower somatic mutation rates, though these could be due to selective pressures acting on other phenotypes such as perhaps cancer. Recent studies in patients with higher somatic mutation burden and no signs of accelerated ageing further question the role of somatic mutations in ageing. Overall, with a few exceptions like cancer, recent DNA sequencing studies and inherited mutations do not support the idea that somatic mutations accumulating with age drive ageing phenotypes, and the phenotypic role, if any, of somatic mutations in ageing remains unclear.

Principles of digital sequencing using unique molecular identifiers

Massively parallel sequencing technologies have long been used in both basic research and clinical routine. The recent introduction of digital sequencing has made previously challenging applications possible by significantly improving sensitivity and specificity to now allow detection of rare sequence variants, even at single molecule level. Digital sequencing utilizes unique molecular identifiers (UMIs) to minimize sequencing-induced errors and quantification biases. Here, we discuss the principles of UMIs and how they are used in digital sequencing. We outline the properties of different UMI types and the consequences of various UMI approaches in relation to experimental protocols and bioinformatics. Finally, we describe how digital sequencing can be applied in specific research fields, focusing on cancer management where it can be used in screening of asymptomatic individuals, diagnosis, treatment prediction, prognostication, monitoring treatment efficacy and early detection of treatment resistance as well as relapse.

Subclonal Cancer Driver Mutations Are Prevalent in the Unresected Peritumoral Edema of Adult Diffuse Gliomas

Adult diffuse gliomas commonly recur regardless of therapy. As recurrence typically arises from the peritumoral edema adjacent to the resected bulk tumor, the profiling of somatic mutations from infiltrative malignant cells within this critical, unresected region could provide important insights into residual disease. A key obstacle has been the inability to distinguish between next-generation sequencing (NGS) noise and the true but weak signal from tumor cells hidden among the noncancerous brain tissue of the peritumoral edema. Here, we developed and validated True2 sequencing to reduce NGS-associated errors to <1 false positive/100 kb panel positions while detecting 97.6% of somatic mutations with an allele frequency ≥0.1%. True2 was then used to study the tumor and peritumoral edema of 22 adult diffuse gliomas including glioblastoma, astrocytoma, oligodendroglioma, and NF1-related low-grade neuroglioma. The tumor and peritumoral edema displayed a similar mutation burden, indicating that surgery debulks these cancers physically but not molecularly. Moreover, variants in the peritumoral edema included unique cancer driver mutations absent in the bulk tumor. Finally, analysis of multiple samples from each patient revealed multiple subclones with unique mutations in the same gene in 17 of 22 patients, supporting the occurrence of convergent evolution in response to patient-specific selective pressures in the tumor microenvironment that may form the molecular foundation of recurrent disease. Collectively, True2 enables the detection of ultralow frequency mutations during molecular analyses of adult diffuse gliomas, which is necessary to understand cancer evolution, recurrence, and individual response to therapy.

SIGNIFICANCE

True2 is a next-generation sequencing workflow that facilitates unbiased discovery of somatic mutations across the full range of variant allele frequencies, which could help identify residual disease vulnerabilities for targeted adjuvant therapies.

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.

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.

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.

Tissue mosaicism following stem cell aging: blood as an exemplar

Loss of stem cell regenerative potential underlies aging of all tissues. Somatic mosaicism, the emergence of cellular patchworks within tissues, increases with age and has been observed in every organ yet examined. In the hematopoietic system, as in most tissues, stem cell aging through a variety of mechanisms occurs in lockstep with the emergence of somatic mosaicism. Here, we draw on insights from aging hematopoiesis to illustrate fundamental principles of stem cell aging and somatic mosaicism. We describe the generalizable changes intrinsic to aged stem cells and their milieu that provide the backdrop for somatic mosaicism to emerge. We discuss genetic and nongenetic mechanisms that can result in tissue somatic mosaicism and existing methodologies to detect such clonal outgrowths. Finally, we propose potential avenues to modify mosaicism during aging, with the ultimate aim of increasing tissue resiliency.

Deciphering cell states and genealogies of human haematopoiesis

The human blood system is maintained through the differentiation and massive amplification of a limited number of long-lived haematopoietic stem cells (HSCs)1. Perturbations to this process underlie diverse diseases, but the clonal contributions to human haematopoiesis and how this changes with age remain incompletely understood. Although recent insights have emerged from barcoding studies in model systems2-5, simultaneous detection of cell states and phylogenies from natural barcodes in humans remains challenging. Here we introduce an improved, single-cell lineage-tracing system based on deep detection of naturally occurring mitochondrial DNA mutations with simultaneous readout of transcriptional states and chromatin accessibility. We use this system to define the clonal architecture of HSCs and map the physiological state and output of clones. We uncover functional heterogeneity in HSC clones, which is stable over months and manifests as both differences in total HSC output and biases towards the production of different mature cell types. We also find that the diversity of HSC clones decreases markedly with age, leading to an oligoclonal structure with multiple distinct clonal expansions. Our study thus provides a clonally resolved and cell-state-aware atlas of human haematopoiesis at single-cell resolution, showing an unappreciated functional diversity of human HSC clones and, more broadly, paving the way for refined studies of clonal dynamics across a range of tissues in human health and disease.

The lifelong natural history of clonal hematopoiesis and its links to myeloid neoplasia

ABSTRACT

The study of somatic mutations and the associated clonal mosaicism across the human body has transformed our understanding of aging and its links to cancer. In proliferative human tissues, stem cells compete for dominance, and those with an advantage expand clonally to outgrow their peers. In the hematopoietic system, such expansion is termed clonal hematopoiesis (CH). The forces driving competition, namely heterogeneity of the hematopoietic stem cell (HSC) pool and attrition of their environment, become increasingly prominent with age. As a result, CH becomes progressively more common through life to the point of becoming essentially ubiquitous. We are beginning to unravel the specific intracellular and extracellular factors underpinning clonal behavior, with somatic mutations in specific driver genes, inflammation, telomere maintenance, extraneous exposures, and inherited genetic variation among the important players. The inevitability of CH with age combined with its unequivocal links to myeloid cancers poses a scientific and clinical challenge. Specifically, we need to decipher the factors determining clonal behavior and develop prognostic tools to identify those at high risk of malignant progression, for whom preventive interventions may be warranted. Here, we discuss how recent advances in our understanding of the natural history of CH have provided important insights into these processes and helped define future avenues of investigation.

MtDNA deletions and aging

Aging is the major risk factor in most of the leading causes of mortality worldwide, yet its fundamental causes mostly remain unclear. One of the clear hallmarks of aging is mitochondrial dysfunction. Mitochondria are best known for their roles in cellular energy generation, but they are also critical biosynthetic and signaling organelles. They also undergo multiple changes with organismal age, including increased genetic errors in their independent, circular genome. A key group of studies looking at mice with increased mtDNA mutations showed that premature aging phenotypes correlated with increased deletions but not point mutations. This generated an interest in mitochondrial deletions as a potential fundamental cause of aging. However, subsequent studies in different models have yielded diverse results. This review summarizes the research on mitochondrial deletions in various organisms to understand their possible roles in causing aging while identifying the key complications in quantifying deletions across all models.

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.

Quantitative and qualitative mutational impact of ionizing radiation on normal cells

The comprehensive genomic impact of ionizing radiation (IR), a carcinogen, on healthy somatic cells remains unclear. Using large-scale whole-genome sequencing (WGS) of clones expanded from irradiated murine and human single cells, we revealed that IR induces a characteristic spectrum of short insertions or deletions (indels) and structural variations (SVs), including balanced inversions, translocations, composite SVs (deletion-insertion, deletion-inversion, and deletion-translocation composites), and complex genomic rearrangements (CGRs), including chromoplexy, chromothripsis, and SV by breakage-fusion-bridge cycles. Our findings suggest that 1 Gy IR exposure causes an average of 2.33 mutational events per Gb genome, comprising 2.15 indels, 0.17 SVs, and 0.01 CGRs, despite a high level of inter-cellular stochasticity. The mutational burden was dependent on total irradiation dose, regardless of dose rate or cell type. The findings were further validated in IR-induced secondary cancers and single cells without clonalization. Overall, our study highlights a comprehensive and clear picture of IR effects on normal mammalian genomes.

Analyzing somatic mutations by single-cell whole-genome sequencing

Somatic mutations are the cause of cancer and have been implicated in other, noncancerous diseases and aging. While clonally expanded mutations can be studied by deep sequencing of bulk DNA, very few somatic mutations expand clonally, and most are unique to each cell. We describe a detailed protocol for single-cell whole-genome sequencing to discover and analyze somatic mutations in tissues and organs. The protocol comprises single-cell multiple displacement amplification (SCMDA), which ensures efficiency and high fidelity in amplification, and the SCcaller software tool to call single-nucleotide variations and small insertions and deletions from the sequencing data by filtering out amplification artifacts. With SCMDA and SCcaller at its core, this protocol describes a complete procedure for the comprehensive analysis of somatic mutations in a single cell, covering (1) single-cell or nucleus isolation, (2) single-cell or nucleus whole-genome amplification, (3) library preparation and sequencing, and (4) computational analyses, including alignment, variant calling, and mutation burden estimation. Methods are also provided for mutation annotation, hotspot discovery and signature analysis. The protocol takes 12-15 h from single-cell isolation to library preparation and 3-7 d of data processing. Compared with other single-cell amplification methods or single-molecular sequencing, it provides high genomic coverage, high accuracy in single-nucleotide variation and small insertions and deletion calling from the same single-cell genome, and fewer processing steps. SCMDA and SCcaller require basic experience in molecular biology and bioinformatics. The protocol can be utilized for studying mutagenesis and genome mosaicism in normal and diseased human and animal tissues under various conditions.

Interpreting the effect of mutations to protein binding sites from large-scale genomic screens

Predicting the functionality of missense mutations is extremely difficult. Large-scale genomic screens are commonly performed to identify mutational correlates or drivers of disease and treatment resistance, but interpretation of how these mutations impact protein function is limited. One such consequence of mutations to a protein is to impact its ability to bind and interact with partners or small molecules such as ATP, thereby modulating its function. Multiple methods exist for predicting the impact of a single mutation on protein-protein binding energy, but it is difficult in the context of a genomic screen to understand if these mutations with large impacts on binding are more common than statistically expected. We present a methodology for taking mutational data from large-scale genomic screens and generating functional and statistical insights into their role in the binding of proteins both with each other and their small molecule ligands. This allows a quantitative and statistical analysis to determine whether mutations impacting protein binding or ligand interactions are occurring more or less frequently than expected by chance. We achieve this by calculating the potential impact of any possible mutation and comparing an expected distribution to the observed mutations. This method is applied to examples demonstrating its ability to interpret mutations involved in protein-protein binding, protein-DNA interactions, and the evolution of therapeutic resistance.

Single-cell lineage tracing with endogenous markers

Resolving lineage relationships between cells in an organism provides key insights into the fate of individual cells and drives a fundamental understanding of the process of development and disease. A recent rapid increase in experimental and computational advances for detecting naturally occurring somatic nuclear and mitochondrial mutation at single-cell resolution has expanded lineage tracing from model organisms to humans. This review discusses the advantages and challenges of experimental and computational techniques for cell lineage tracing using somatic mutation as endogenous DNA barcodes to decipher the relationships between cells during development and tumour evolution. We outlook the advantages of spatial clonal evolution analysis and single-cell lineage tracing using endogenous genetic markers.

Impact of Rare Structural Variant Events in Newly Diagnosed Multiple Myeloma

PURPOSE

Whole-genome sequencing (WGS) of patients with newly diagnosed multiple myeloma (NDMM) has shown recurrent structural variant (SV) involvement in distinct regions of the genome (i.e., hotspots) and causing recurrent copy-number alterations. Together with canonical immunoglobulin translocations, these SVs are recognized as "recurrent SVs." More than half of SVs were not involved in recurrent events. The significance of these "rare SVs" has not been previously examined.

EXPERIMENTAL DESIGN

In this study, we utilize 752 WGS and 591 RNA sequencing data from patients with NDMM to determine the role of rare SVs in myeloma pathogenesis.

RESULTS

Ninety-four percent of patients harbored at least one rare SV event. Rare SVs showed an SV class-specific enrichment within genes and superenhancers associated with outlier gene expression. Furthermore, known myeloma driver genes recurrently impacted by point mutations were dysregulated by rare SVs.

CONCLUSIONS

Overall, we demonstrate the association of rare SVs with aberrant gene expression supporting a potential driver role in myeloma pathogenesis.

Genome Instability and DNA Repair in Somatic and Reproductive Aging

Genetic material is constantly subjected to genotoxic insults and is critically dependent on DNA repair. Genome maintenance mechanisms differ in somatic and germ cells as the soma only requires maintenance during an individual's lifespan, while the germline indefinitely perpetuates its genetic information. DNA lesions are recognized and repaired by mechanistically highly diverse repair machineries. The DNA damage response impinges on a vast array of homeostatic processes and can ultimately result in cell fate changes such as apoptosis or cellular senescence. DNA damage causally contributes to the aging process and aging-associated diseases, most prominently cancer. By causing mutations, DNA damage in germ cells can lead to genetic diseases and impact the evolutionary trajectory of a species. The mechanisms ensuring tight control of germline DNA repair could be highly instructive in defining strategies for improved somatic DNA repair. They may provide future interventions to maintain health and prevent disease during aging.

Measures of genetic diversification in somatic tissues at bulk and single-cell resolution

Intra-tissue genetic heterogeneity is universal to both healthy and cancerous tissues. It emerges from the stochastic accumulation of somatic mutations throughout development and homeostasis. By combining population genetics theory and genomic information, genetic heterogeneity can be exploited to infer tissue organization and dynamics in vivo. However, many basic quantities, for example the dynamics of tissue-specific stem cells remain difficult to quantify precisely. Here, we show that single-cell and bulk sequencing data inform on different aspects of the underlying stochastic processes. Bulk-derived variant allele frequency spectra (VAF) show transitions from growing to constant stem cell populations with age in samples of healthy esophagus epithelium. Single-cell mutational burden distributions allow a sample size independent measure of mutation and proliferation rates. Mutation rates in adult hematopietic stem cells are higher compared to inferences during development, suggesting additional proliferation-independent effects. Furthermore, single-cell derived VAF spectra contain information on the number of tissue-specific stem cells. In hematopiesis, we find approximately 2 × 105 HSCs, if all stem cells divide symmetrically. However, the single-cell mutational burden distribution is over-dispersed compared to a model of Poisson distributed random mutations. A time-associated model of mutation accumulation with a constant rate alone cannot generate such a pattern. At least one additional source of stochasticity would be needed. Possible candidates for these processes may be occasional bursts of stem cell divisions, potentially in response to injury, or non-constant mutation rates either through environmental exposures or cell-intrinsic variation.

Cancer as a Disease of Development Gone Awry

In the 160 years since Rudolf Virchow first postulated that neoplasia arises by the same law that regulates embryonic development, scientists have come to recognize the striking overlap between the molecular and cellular programs used by cancers and embryos. Advances in cancer biology and molecular techniques have further highlighted the similarities between carcinogenesis and embryogenesis, where cellular growth, differentiation, motility, and intercellular cross talk are mediated by common drivers and regulatory networks. This review highlights the many connections linking cancer biology and developmental biology to provide a deeper understanding of how a tissue's developmental history may both enable and constrain cancer cell evolution.

Genetic profiling of rat gliomas and cardiac schwannomas from life-time radiofrequency radiation exposure study using a targeted next-generation sequencing gene panel

The cancer hazard associated with lifetime exposure to radiofrequency radiation (RFR) was examined in Sprague Dawley (SD) rats at the Ramazzini Institute (RI), Italy. There were increased incidences of gliomas and cardiac schwannomas. The translational relevance of these rare rat tumors for human disease is poorly understood. We examined the genetic alterations in RFR-derived rat tumors through molecular characterization of important cancer genes relevant for human gliomagenesis. A targeted next-generation sequencing (NGS) panel was designed for rats based on the top 23 orthologous human glioma-related genes. Single-nucleotide variants (SNVs) and small insertion and deletions (indels) were characterized in the rat gliomas and cardiac schwannomas. Translational relevance of these genetic alterations in rat tumors to human disease was determined through comparison with the Catalogue of Somatic Mutations in Cancer (COSMIC) database. These data suggest that rat gliomas resulting from life-time exposure to RFR histologically resemble low grade human gliomas but surprisingly no mutations were detected in rat gliomas that had homology to the human IDH1 p.R132 or IDH2 p.R172 suggesting that rat gliomas are primarily wild-type for IDH hotspot mutations implicated in human gliomas. The rat gliomas appear to share some genetic alterations with IDH1 wildtype human gliomas and rat cardiac schwannomas also harbor mutations in some of the queried cancer genes. These data demonstrate that targeted NGS panels based on tumor specific orthologous human cancer driver genes are an important tool to examine the translational relevance of rodent tumors resulting from chronic/life-time rodent bioassays.

Cancer Evolution: A Multifaceted Affair

Cancer cells adapt and survive through the acquisition and selection of molecular modifications. This process defines cancer evolution. Building on a theoretical framework based on heritable genetic changes has provided insights into the mechanisms supporting cancer evolution. However, cancer hallmarks also emerge via heritable nongenetic mechanisms, including epigenetic and chromatin topological changes, and interactions between tumor cells and the tumor microenvironment. Recent findings on tumor evolutionary mechanisms draw a multifaceted picture where heterogeneous forces interact and influence each other while shaping tumor progression. A comprehensive characterization of the cancer evolutionary toolkit is required to improve personalized medicine and biomarker discovery.

SIGNIFICANCE

Tumor evolution is fueled by multiple enabling mechanisms. Importantly, genetic instability, epigenetic reprogramming, and interactions with the tumor microenvironment are neither alternative nor independent evolutionary mechanisms. As demonstrated by findings highlighted in this perspective, experimental and theoretical approaches must account for multiple evolutionary mechanisms and their interactions to ultimately understand, predict, and steer tumor evolution.

Scrutiny of genome-wide somatic mutation profiles in centenarians identifies the key genomic regions for human longevity

Somatic mutations accumulate with age and are associated closely with human health, their characterization in longevity cohorts remains largely unknown. Here, by analyzing whole genome somatic mutation profiles in 73 centenarians and 51 younger controls in China, we found that centenarian genomes are characterized by a markedly skewed distribution of somatic mutations, with many genomic regions being specifically conserved but displaying a high function potential. This, together with the observed more efficient DNA repair ability in the long-lived individuals, supports the existence of key genomic regions for human survival during aging, with their integrity being of essential to human longevity.

Adult stem cell activity in naked mole rats for long-term tissue maintenance

The naked mole rat (NMR), Heterocephalus glaber, the longest-living rodent, provides a unique opportunity to explore how evolution has shaped adult stem cell (ASC) activity and tissue function with increasing lifespan. Using cumulative BrdU labelling and a quantitative imaging approach to track intestinal ASCs (Lgr5+) in their native in vivo state, we find an expanded pool of Lgr5+ cells in NMRs, and these cells specifically at the crypt base (Lgr5+CBC) exhibit slower division rates compared to those in short-lived mice but have a similar turnover as human LGR5+CBC cells. Instead of entering quiescence (G0), NMR Lgr5+CBC cells reduce their division rates by prolonging arrest in the G1 and/or G2 phases of the cell cycle. Moreover, we also observe a higher proportion of differentiated cells in NMRs that confer enhanced protection and function to the intestinal mucosa which is able to detect any chemical imbalance in the luminal environment efficiently, triggering a robust pro-apoptotic, anti-proliferative response within the stem/progenitor cell zone.

Impact of CRISPR/HDR-editing versus lentiviral transduction on long-term engraftment and clonal dynamics of HSPCs in rhesus macaques

For precise genome editing via CRISPR/homology-directed repair (HDR), effective and safe editing of long-term engrafting hematopoietic stem cells (LT-HSCs) requires both sufficient HDR efficiency and protection of LT-HSC function and number. The impact of HDR on true LT-HSCs clonal dynamics in a relevant large animal model has not previously been studied. To track the HDR-edited cells, autologous rhesus macaque (RM) CD34 + cells were electroporated with the gRNA/Cas9 ribonucleoprotein (RNP) and HDR cassette barcode library structure and reinfused into RMs following myeloablation. For competitive model animals, fractionated CD34 + cells were transduced with a barcoded GFP-expressing lentiviral vector (LV) and electroporated via HDR machinery, respectively. CD33 knockout (KO) neutrophils were prevalent early following engraftment and then rapidly decreased, resulting in less than 1% total editing efficiency. Interestingly, in competitive animals, a higher concentration of i53 mRNA result in a less steep reduction in CD33 KO cells, presented a modest decrease in HDR rate (0.1-0.2%) and total indels (1.5-6.5%). In contrast, the drop off of LV-transduced GFP + cells stabilized at 20% after 2 months. We next retrieved embedded barcodes and revealed that various clones contributed to early hematopoietic reconstitution, then after dominant clones appeared at steady state throughout the animals. In conclusion, CRISPR/HDR edited cells disappeared rapidly after the autologous transplantation in RM despite substantial gene editing outcome, whereas LV-transduced cells were relatively well maintained. Clonality of HDR-edited cells drastically shrank at early stage and then relied on several dominant clones, which can be mildly mitigated by the introduction of i53 mRNA.

A naturally occurring variant of MBD4 causes maternal germline hypermutation in primates

As part of an ongoing genome sequencing project at the Oregon National Primate Research Center, we identified a rhesus macaque with a rare homozygous frameshift mutation in the gene methyl-CpG binding domain 4, DNA glycosylase (MBD4). MBD4 is responsible for the repair of C > T deamination mutations at CpG dinucleotides and has been linked to somatic hypermutation and cancer predisposition in humans. We show here that MBD4-associated hypermutation also affects the germline: The six offspring of the MBD4-null dam have a fourfold to sixfold increase in de novo mutation burden. This excess burden was predominantly C > T mutations at CpG dinucleotides consistent with MBD4 loss of function in the dam. There was also a significant excess of C > T at CpA sites, indicating an important, unappreciated role for MBD4 to repair deamination in CpA contexts. The MBD4-null dam developed sustained eosinophilia later in life, but we saw no other signs of neoplastic processes associated with MBD4 loss of function in humans nor any obvious disease in the hypermutated offspring. This work provides the first evidence for a genetic factor causing hypermutation in the maternal germline of a mammal and adds to the very small list of naturally occurring variants known to modulate germline mutation rates in mammals.

Detection of APP gene recombinant in human blood plasma

The pathogenesis of Alzheimer's disease (AD) is believed to involve the accumulation of amyloid-β in the brain, which is produced by the sequential cleavage of amyloid precursor protein (APP) by β-secretase and γ-secretase. Recently, analysis of genomic DNA and mRNA from postmortem brain neurons has revealed intra-exonic recombinants of APP (gencDNA), which have been implicated in the accumulation of amyloid-β. In this study, we computationally analyzed publicly available sequence data (SRA) using probe sequences we constructed to screen APP gencDNAs. APP gencDNAs were detected in SRAs constructed from both genomic DNA and RNA obtained from the postmortem brain and in the SRA constructed from plasma cell-free mRNA (cf-mRNA). The SRA constructed from plasma cf-mRNA showed a significant difference in the number of APP gencDNA reads between SAD and NCI: the p-value from the Mann-Whitney U test was 5.14 × 10-6. The transcripts were also found in circulating nucleic acids (CNA) from our plasma samples with NGS analysis. These data indicate that transcripts of APP gencDNA can be detected in blood plasma and suggest the possibility of using them as blood biomarkers for Alzheimer's disease.

Detection of low-frequency mutations in clinical samples by increasing mutation abundance via the excision of wild-type sequences

The efficiency of DNA-enrichment techniques is often insufficient to detect mutations that occur at low frequencies. Here we report a DNA-excision method for the detection of low-frequency mutations in genomic DNA and in circulating cell-free DNA at single-nucleotide resolution. The method is based on a competitive DNA-binding-and-digestion mechanism, effected by deoxyribonuclease I (DNase) guided by single-stranded phosphorothioated DNA (sgDNase), for the removal of wild-type DNA strands. The sgDNase can be designed against any wild-type DNA sequences, allowing for the uniform enrichment of all the mutations within the target-binding region of single-stranded phosphorothioated DNA at mild-temperature conditions. Pretreatment with sgDNase enriches all mutant strands with initial frequencies down to 0.01% and leads to high discrimination factors for all types of single-nucleotide mismatch in multiple sequence contexts, as we show for the identification of low-abundance mutations in samples of blood or tissue from patients with cancer. The method can be coupled with next-generation sequencing, droplet digital polymerase chain reaction, Sanger sequencing, fluorescent-probe-based assays and other mutation-detection methods.

Somatic mutation burden in relation to aging and functional life span: implications for cellular reprogramming and rejuvenation

The accrual of somatic mutations has been implicated as causal factors in aging since the 1950s. However, the quantitative analysis of somatic mutations has posed a major challenge due to the random nature of de novo mutations in normal tissues, which has limited analysis to tumors and other clonal lineages. Advances in single-cell and single-molecule next-generation sequencing now allow to obtain, for the first time, detailed insights into the landscape of somatic mutations in different human tissues and cell types as a function of age under various conditions. Here, we will briefly recapitulate progress in somatic mutation analysis and discuss the possible relationship between somatic mutation burden with functional life span, with a focus on differences between germ cells, stem cells, and differentiated cells.

Genome-wide direct quantification of in vivo mutagenesis using high-accuracy paired-end and complementary consensus sequencing

Error-corrected next-generation sequencing (ecNGS) is an emerging technology for accurately measuring somatic mutations. Here, we report paired-end and complementary consensus sequencing (PECC-Seq), a high-accuracy ecNGS approach for genome-wide somatic mutation detection. We characterize a novel 2-aminoimidazolone lesion besides 7,8-dihydro-8-oxoguanine and the resulting end-repair artifacts originating from NGS library preparation that obscure the sequencing accuracy of NGS. We modify library preparation protocol for the enzymatic removal of end-repair artifacts and improve the accuracy of our previously developed duplex consensus sequencing method. Optimized PECC-Seq shows an error rate of <5 × 10-8 with consensus bases compressed from approximately 25 Gb of raw sequencing data, enabling the accurate detection of low-abundance somatic mutations. We apply PECC-Seq to the quantification of in vivo mutagenesis. Compared with the classic gpt gene mutation assay using gpt delta transgenic mice, PECC-Seq exhibits high sensitivity in quantitatively measuring dose-dependent mutagenesis induced by Aristolochic acid I (AAI). Moreover, PECC-Seq specifically characterizes the distinct genome-wide mutational signatures of AAI, Benzo[a]pyrene, N-Nitroso-N-ethylurea and N-nitrosodiethylamine and reveals the mutational signature of Quinoline in common mouse models. Overall, our findings demonstrate that high-accuracy PECC-Seq is a promising tool for genome-wide somatic mutagenesis quantification and for in vivo mutagenicity testing.

Detection of in vivo mutagenicity in rat liver samples using error-corrected sequencing techniques

BACKGROUND

Mutagenicity, the ability of chemical agents to cause mutations and potentially lead to cancer, is a critical aspect of substance safety assessment for protecting human health and the environment. Metabolic enzymes activate multiple mutagens in living organisms, thus in vivo animal models provide highly important information for evaluating mutagenicity in human. Rats are considered suitable models as they share a similar metabolic pathway with humans for processing toxic chemical and exhibit higher responsiveness to chemical carcinogens than mice. To assess mutagenicity in rats, transgenic rodents (TGRs) are widely used for in vivo gene mutation assays. However, such assays are labor-intensive and could only detect transgene mutations inserted into the genome. Therefore, introducing a technology to directly detect in vivo mutagenicity in rats would be necessary. The next-generation sequencing (NGS) based error-corrected sequencing technique is a promising approach for such purposes.

RESULTS

We investigated the applicability of paired-end and complementary consensus sequencing (PECC-Seq), an error-corrected sequencing technique, for detecting in vivo mutagenicity in the rat liver samples. PECC-Seq allows for the direct detection of ultra-rare somatic mutations in the genomic DNA without being constrained by the genomic locus, tissue, or organism. We tested PECC-Seq feasibility in rats treated with diethylnitrosamine (DEN), a mutagenic compound. Interestingly, the mutation and mutant frequencies between PECC-Seq and the TGR assay displayed a promising correlation. Our results also demonstrated that PECC-Seq could successfully detect the A:T > T:A mutation in rat liver samples, consistent with the TGR assay. Furthermore, we calculated the trinucleotide mutation frequency and proved that PECC-Seq accurately identified the DEN treatment-induced mutational signatures.

CONCLUSIONS

Our study provides the first evidence of using PECC-Seq for in vivo mutagenicity detection in rat liver samples. This approach could provide a valuable alternative to conventional TGR assays as it is labor- and time-efficient and eliminates the need for transgenic rodents. Error-corrected sequencing techniques, such as PECC-Seq, represent promising approaches for enhancing mutagenicity assessment and advancing regulatory science.

The 'faulty male' hypothesis for sex-biased mutation and disease

Biological differences between males and females lead to many differences in physiology, disease, and overall health. One of the most prominent disparities is in the number of germline mutations passed to offspring: human males transmit three times as many mutations as do females. While the classic explanation for this pattern invokes differences in post-puberty germline replication between the sexes, recent whole-genome evidence in humans and other mammals has cast doubt on this mechanism. Here, we review recent work that is inconsistent with a replication-driven model of male-biased mutation, and propose an alternative, 'faulty male' hypothesis. This model proposes that males are less able to repair and/or protect DNA from damage compared to females. Importantly, we suggest that this new model for male-biased mutation may also help to explain several pronounced differences between the sexes in cancer, aging, and DNA repair. Although the detailed contributions of genetic, epigenetic, and hormonal influences of biological sex on mutation remain to be fully understood, a reconsideration of the mechanisms underlying these differences will lead to a deeper understanding of evolution and disease.

Single-cell transcriptomic and genomic changes in the aging human brain

Aging brings dysregulation of various processes across organs and tissues, often stemming from stochastic damage to individual cells over time. Here, we used a combination of single-nucleus RNA-sequencing and single-cell whole-genome sequencing to identify transcriptomic and genomic changes in the prefrontal cortex of the human brain across life span, from infancy to centenarian. We identified infant-specific cell clusters enriched for the expression of neurodevelopmental genes, and a common down-regulation of cell-essential homeostatic genes that function in ribosomes, transport, and metabolism during aging across cell types. Conversely, expression of neuron-specific genes generally remains stable throughout life. We observed a decrease in specific DNA repair genes in aging, including genes implicated in generating brain somatic mutations as indicated by mutation signature analysis. Furthermore, we detected gene-length-specific somatic mutation rates that shape the transcriptomic landscape of the aged human brain. These findings elucidate critical aspects of human brain aging, shedding light on transcriptomic and genomics dynamics.

Stem cell mutations, associated cancer risk, and consequences for regenerative medicine

Mutation accumulation in stem cells has been associated with cancer risk. However, the presence of numerous mutant clones in healthy tissues has raised the question of what limits cancer initiation. Here, we review recent developments in characterizing mutation accumulation in healthy tissues and compare mutation rates in stem cells during development and adult life with corresponding cancer risk. A certain level of mutagenesis within the stem cell pool might be beneficial to limit the size of malignant clones through competition. This knowledge impacts our understanding of carcinogenesis with potential consequences for the use of stem cells in regenerative medicine.

Mutation divergence over space in tumour expansion

Mutation accumulation in tumour evolution is one major cause of intra-tumour heterogeneity (ITH), which often leads to drug resistance during treatment. Previous studies with multi-region sequencing have shown that mutation divergence among samples within the patient is common, and the importance of spatial sampling to obtain a complete picture in tumour measurements. However, quantitative comparisons of the relationship between mutation heterogeneity and tumour expansion modes, sampling distances as well as the sampling methods are still few. Here, we investigate how mutations diverge over space by varying the sampling distance and tumour expansion modes using individual-based simulations. We measure ITH by the Jaccard index between samples and quantify how ITH increases with sampling distance, the pattern of which holds in various sampling methods and sizes. We also compare the inferred mutation rates based on the distributions of variant allele frequencies under different tumour expansion modes and sampling sizes. In exponentially fast expanding tumours, a mutation rate can always be inferred for any sampling size. However, the accuracy compared with the true value decreases when the sampling size decreases, where small sampling sizes result in a high estimate of the mutation rate. In addition, such an inference becomes unreliable when the tumour expansion is slow, such as in surface growth.

Oxidative DNA damage and repair at non-coding regulatory regions

DNA breaks at protein-coding sequences are well-established threats to tissue homeostasis and maintenance. They arise from the exposure to intracellular and environmental genotoxins, causing damage in one or two strands of the DNA. DNA breaks have been also reported in non-coding regulatory regions such as enhancers and promoters. They arise from essential cellular processes required for gene transcription, cell identity and function. One such process that has attracted recent attention is the oxidative demethylation of DNA and histones, which generates abasic sites and DNA single-strand breaks. Here, we discuss how oxidative DNA breaks at non-coding regulatory regions are generated and the recently reported role of NuMA (nuclear mitotic apparatus) protein in promoting transcription and repair at these regions.

The evolution of aging and lifespan

Aging is a nearly inescapable trait among organisms yet lifespan varies tremendously across different species and spans several orders of magnitude in vertebrates alone. This vast phenotypic diversity is driven by distinct evolutionary trajectories and tradeoffs that are reflected in patterns of diversification and constraint in organismal genomes. Age-specific impacts of selection also shape allele frequencies in populations, thus impacting disease susceptibility and environment-specific mortality risk. Further, the mutational processes that spawn this genetic diversity in both germline and somatic cells are strongly influenced by age and life history. We discuss recent advances in our understanding of the evolution of aging and lifespan at organismal, population, and cellular scales, and highlight outstanding questions that remain unanswered.

Effects of psoriasis and psoralen exposure on the somatic mutation landscape of the skin

Somatic mutations are hypothesized to play a role in many non-neoplastic diseases. We performed whole-exome sequencing of 1,182 microbiopsies dissected from lesional and nonlesional epidermis from 111 patients with psoriasis to search for evidence that somatic mutations in keratinocytes may influence the disease process. Lesional skin remained highly polyclonal, showing no evidence of large-scale spread of clones carrying potentially pathogenic mutations. The mutation rate of keratinocytes was similarly only modestly affected by the disease. We found evidence of positive selection in previously reported driver genes NOTCH1, NOTCH2, TP53, FAT1 and PPM1D and also identified mutations in four genes (GXYLT1, CHEK2, ZFP36L2 and EEF1A1) that we hypothesize are selected for in squamous epithelium irrespective of disease status. Finally, we describe a mutational signature of psoralens-a class of chemicals previously found in some sunscreens and which are used as part of PUVA (psoralens and ultraviolet-A) photochemotherapy treatment for psoriasis.

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.

Characterization of Cancer Evolution Landscape Based on Accurate Detection of Somatic Mutations in Single Tumor Cells

Accurate detection of somatic mutations in single tumor cells is greatly desired as it allows us to quantify the single-cell mutation burden and construct the mutation-based phylogenetic tree. Here we developed scNanoSeq chemistry and profiled 842 single cells from 21 human breast cancer samples. The majority of the mutation-based phylogenetic trees comprise a characteristic stem evolution followed by the clonal sweep. We observed the subtype-dependent lengths in the stem evolution. To explain this phenomenon, we propose that the differences are related to different reprogramming required for different subtypes of breast cancer. Furthermore, we reason that the time that the tumor-initiating cell took to acquire the critical clonal-sweep-initiating mutation by random chance set the time limit for the reprogramming process. We refer to this model as a reprogramming and critical mutation co-timing (RCMC) subtype model. Next, in the sweeping clone, we observed that tumor cells undergo a branched evolution with rapidly decreasing selection. In the most recent clades, effectively neutral evolution has been reached, resulting in a substantially large number of mutational heterogeneities. Integrative analysis with 522-713X ultra-deep bulk whole genome sequencing (WGS) further validated this evolution mode. Mutation-based phylogenetic trees also allow us to identify the early branched cells in a few samples, whose phylogenetic trees support the gradual evolution of copy number variations (CNVs). Overall, the development of scNanoSeq allows us to unveil novel insights into breast cancer evolution.

Whole-genome high-fidelity sequencing: A novel approach to detecting and characterization of mutagenicity in vivo

Direct DNA sequencing can be used for characterizing mutagenicity in simple and complex biological models. Recently we described a method of whole-genome sequencing for detecting mutations in simple models of cultured bacteria, mammalian cells, and nematode. In the current proof-of-concept study, we expand and improve our method for evaluating a more complex mammalian biological model in outbred mice. We detail the method by applying it to a small set of animals treated with a mutagen with known mutagenicity profiles, N-ethyl-N-nitrosourea (ENU), for consistency with the known data. Whole-genome high-fidelity sequencing (HiFi Sequencing) showed frequencies and spectra of background mutations in tissues of untreated mice that were consistent with normal ageing and characterized by spontaneous or enzymatic deamination of 5-methylcytosine. In mice treated with a single 40 mg/kg dose of ENU, the frequency of mutations in the genomic DNA of solid tissues increased up to 7-fold, with the greatest increase observed in the spleen and the smallest increase in the liver. The most common mutations detected in ENU-treated mice were T > A transitions and T > C transversions, consistent with the types of mutations caused by alkylating agents. The data suggest that HiFi Sequencing may be useful for characterizing mutagenicity of novel compounds in various biological models.

Neuronal Senescence in the Aged Brain

Cellular senescence is a highly complicated cellular state that occurs throughout the lifespan of an organism. It has been well-defined in mitotic cells by various senescent features. Neurons are long-lived post-mitotic cells with special structures and functions. With age, neurons display morphological and functional changes, accompanying alterations in proteostasis, redox balance, and Ca2+ dynamics; however, it is ambiguous whether these neuronal changes belong to the features of neuronal senescence. In this review, we strive to identify and classify changes that are relatively specific to neurons in the aging brain and define them as features of neuronal senescence through comparisons with common senescent features. We also associate them with the functional decline of multiple cellular homeostasis systems, proposing the possibility that these systems are the main drivers of neuronal senescence. We hope this summary will serve as a steppingstone for further inputs on a comprehensive but relatively specific list of phenotypes for neuronal senescence and in particular their underlying molecular events during aging. This will in turn shine light on the association between neuronal senescence and neurodegeneration and lead to the development of strategies to perturb the processes.