Single-cell genome sequencing of human neurons identifies somatic point mutation and indel enrichment in regulatory elements

Mitochondrial genetics through the lens of single-cell multi-omics

Mitochondria carry their own genetic information encoding for a subset of protein-coding genes and translational machinery essential for cellular respiration and metabolism. Despite its small size, the mitochondrial genome, its natural genetic variation and molecular phenotypes have been challenging to study using bulk sequencing approaches, due to its variation in cellular copy number, non-Mendelian modes of inheritance and propensity for mutations. Here we highlight emerging strategies designed to capture mitochondrial genetic variation across individual cells for lineage tracing and studying mitochondrial genetics in primary human cells and clinical specimens. We review recent advances surrounding single-cell mitochondrial genome sequencing and its integration with functional genomic readouts, including leveraging somatic mitochondrial DNA mutations as clonal markers that can resolve cellular population dynamics in complex human tissues. Finally, we discuss how single-cell whole mitochondrial genome sequencing approaches can be utilized to investigate mitochondrial genetics and its contribution to cellular heterogeneity and disease.

Epigenetic signatures of trophoblast lineage and their biological functions

Trophoblasts play a crucial role in embryo implantation and in interacting with the maternal uterus. The trophoblast lineage develops into a substantial part of the placenta, a temporary extra-embryonic organ, capable of undergoing distinctive epigenetic events during development. The critical role of trophoblast-specific epigenetic signatures in regulating placental development has become known, significantly advancing our understanding of trophoblast identity and lineage development. Scientific efforts are revealing how trophoblast-specific epigenetic signatures mediate stage-specific gene regulatory programming during the development of the trophoblast lineage. These epigenetic signatures have a significant impact on blastocyst formation, placental development, as well as the growth and survival of embryos and fetuses. In evolution, DNA hypomethylation in the trophoblast lineage is conserved, and there is a significant disparity in the control of epigenetic dynamics and the landscape of genomic imprinting. Scientists have used murine and human multipotent trophoblast cells as in vitro models to recapitulate the essential epigenetic processes of placental development. Here, we review the epigenetic signatures of the trophoblast lineage and their biological functions to enhance our understanding of placental evolution, development, and function.

Identification of a mosaic MTOR variant in purified neuronal DNA in a patient with focal cortical dysplasia using a novel depth electrode harvesting technique

OBJECTIVE

Recent studies have identified brain somatic variants as a cause of focal epilepsy. These studies relied on resected tissue from epilepsy surgery, which is not available in most patients. The use of trace tissue adherent to depth electrodes used for stereo electroencephalography (EEG) has been proposed as an alternative but is hampered by the low cell quality and contamination by nonbrain cells. Here, we use our improved depth electrode harvesting technique that purifies neuronal nuclei to achieve molecular diagnosis in a patient with focal cortical dysplasia (FCD).

METHODS

Depth electrode tips were collected, pooled by brain region and seizure onset zone, and nuclei were isolated and sorted using fluorescence-activated nuclei sorting (FANS). Somatic DNA was amplified from neuronal and astrocyte nuclei using primary template amplification followed by exome sequencing of neuronal DNA from the affected pool, unaffected pool, and saliva. The identified variant was validated using droplet digital polymerase chain reaction (PCR).

RESULTS

An 11-year-old male with drug-resistant genetic-structural epilepsy due to left anterior insula FCD had seizures from age 3 years. Stereo EEG confirmed seizure onset in the left anterior insula. The two anterior insula electrodes were combined as the affected pool and three frontal electrodes as the unaffected pool. FANS isolated 140 neuronal nuclei from the affected and 245 neuronal nuclei from the unaffected pool. A novel somatic missense MTOR variant (p.Leu489Met, CADD score 23.7) was identified in the affected neuronal sample. Droplet digital PCR confirmed a mosaic gradient (variant allele frequency = .78% in affected neuronal sample; variant was absent in all other samples).

SIGNIFICANCE

Our findings confirm that harvesting neuronal DNA from depth electrodes followed by molecular analysis to identify brain somatic variants is feasible. Our novel method represents a significant improvement compared to the previous method by focusing the analysis on high-quality cells of the cell type of interest.

Large-scale neurophysiology and single-cell profiling in human neuroscience

Advances in large-scale single-unit human neurophysiology, single-cell RNA sequencing, spatial transcriptomics and long-term ex vivo tissue culture of surgically resected human brain tissue have provided an unprecedented opportunity to study human neuroscience. In this Perspective, we describe the development of these paradigms, including Neuropixels and recent brain-cell atlas efforts, and discuss how their convergence will further investigations into the cellular underpinnings of network-level activity in the human brain. Specifically, we introduce a workflow in which functionally mapped samples of human brain tissue resected during awake brain surgery can be cultured ex vivo for multi-modal cellular and functional profiling. We then explore how advances in human neuroscience will affect clinical practice, and conclude by discussing societal and ethical implications to consider. Potential findings from the field of human neuroscience will be vast, ranging from insights into human neurodiversity and evolution to providing cell-type-specific access to study and manipulate diseased circuits in pathology. This Perspective aims to provide a unifying framework for the field of human neuroscience as we welcome an exciting era for understanding the functional cytoarchitecture of the human brain.

Reconstructing phylogenetic trees from genome-wide somatic mutations in clonal samples

Phylogenetic trees are a powerful means to display the evolutionary history of species, pathogens and, more recently, individual cells of the human body. Whole-genome sequencing of laser capture microdissections or expanded stem cells has allowed the discovery of somatic mutations in clones, which can be used as natural barcodes to reconstruct the developmental history of individual cells. Here we describe Sequoia, our pipeline to reconstruct lineage trees from clones of normal cells. Candidate somatic mutations are called against the human reference genome and filtered to exclude germline mutations and artifactual variants. These filtered somatic mutations form the basis for phylogeny reconstruction using a maximum parsimony framework. Lastly, we use a maximum likelihood framework to explicitly map mutations to branches in the phylogenetic tree. The resulting phylogenies can then serve as a basis for many subsequent analyses, including investigating embryonic development, tissue dynamics in health and disease, and mutational signatures. Sequoia can be readily applied to any clonal somatic mutation dataset, including single-cell DNA sequencing datasets, using the commands and scripts provided. Moreover, Sequoia is highly flexible and can be easily customized. Typically, the runtime of the core script ranges from minutes to an hour for datasets with a moderate number (50,000-150,000) of variants. Competent bioinformatic skills, including in-depth knowledge of the R programming language, are required. A high-performance computing cluster (one that is capable of running mutation-calling algorithms and other aspects of the analysis at scale) is also required, especially if handling large datasets.

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.

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.

Accurate and sensitive mutational signature analysis with MuSiCal

Mutational signature analysis is a recent computational approach for interpreting somatic mutations in the genome. Its application to cancer data has enhanced our understanding of mutational forces driving tumorigenesis and demonstrated its potential to inform prognosis and treatment decisions. However, methodological challenges remain for discovering new signatures and assigning proper weights to existing signatures, thereby hindering broader clinical applications. Here we present Mutational Signature Calculator (MuSiCal), a rigorous analytical framework with algorithms that solve major problems in the standard workflow. Our simulation studies demonstrate that MuSiCal outperforms state-of-the-art algorithms for both signature discovery and assignment. By reanalyzing more than 2,700 cancer genomes, we provide an improved catalog of signatures and their assignments, discover nine indel signatures absent in the current catalog, resolve long-standing issues with the ambiguous 'flat' signatures and give insights into signatures with unknown etiologies. We expect MuSiCal and the improved catalog to be a step towards establishing best practices for mutational signature analysis.

Pangenomics: A new era in the field of neurodegenerative diseases

A pangenome is composed of all the genetic variability of a group of individuals, and its application to the study of neurodegenerative diseases may provide valuable insights into the underlying aspects of genetic heterogenetiy for these complex ailments, including gene expression, epigenetics, and translation mechanisms. Furthermore, a reference pangenome allows for the identification of previously undetected structural commonalities and differences among individuals, which may help in the diagnosis of a disease, support the prediction of what will happen over time (prognosis) and aid in developing novel treatments in the perspective of personalized medicine. Therefore, in the present review, the application of the pangenome concept to the study of neurodegenerative diseases will be discussed and analyzed for its potential to enable an improvement in diagnosis and prognosis for these illnesses, leading to the development of tailored treatments for individual patients from the knowledge of the genomic composition of a whole population.

DNA damage and transcription stress

Genome damage and transcription are intimately linked. Tens to hundreds of thousands of DNA lesions arise in each cell each day, many of which can directly or indirectly impede transcription. Conversely, the process of gene expression is itself a source of endogenous DNA lesions as a result of the susceptibility of single-stranded DNA to damage, conflicts with the DNA replication machinery, and engagement by cells of topoisomerases and base excision repair enzymes to regulate the initiation and progression of gene transcription. Although such processes are tightly regulated and normally accurate, on occasion, they can become abortive and leave behind DNA breaks that can drive genome rearrangements, instability, or cell death.

Cell type-specific enrichment of somatic aneuploidy in the mammalian brain

Somatic mutations alter the genomes of a subset of an individual's brain cells1-3, impacting gene regulation and contributing to disease processes4,5. Mosaic single nucleotide variants have been characterized with single-cell resolution in the brain2,3, but we have limited information about large-scale structural variation, including whole-chromosome duplication or loss1,6,7. We used a dataset of over 415,000 single-cell DNA methylation and chromatin conformation profiles across the adult mouse brain to identify aneuploid cells comprehensively. Whole-chromosome loss or duplication occurred in <1% of cells, with rates up to 1.8% in non-neuronal cell types, including oligodendrocyte precursors and pericytes. Among all aneuploidies, we observed a strong enrichment of trisomy on chromosome 16, which is syntenic with human chromosome 21 and constitutively trisomic in Down syndrome. Chromosome 16 trisomy occurred in multiple cell types and across brain regions, suggesting that nondisjunction is a recurrent feature of somatic variation in the brain.

Single-cell somatic copy number variants in brain using different amplification methods and reference genomes

The presence of somatic mutations, including copy number variants (CNVs), in the brain is well recognized. Comprehensive study requires single-cell whole genome amplification, with several methods available, prior to sequencing. We compared PicoPLEX with two recent adaptations of multiple displacement amplification (MDA): primary template-directed amplification (PTA) and droplet MDA, across 93 human brain cortical nuclei. We demonstrated different properties for each, with PTA providing the broadest amplification, PicoPLEX the most even, and distinct chimeric profiles. Furthermore, we performed CNV calling on two brains with multiple system atrophy and one control brain using different reference genomes. We found that 38% of brain cells have at least one Mb-scale CNV, with some supported by bulk sequencing or single-cells from other brain regions. Our study highlights the importance of selecting whole genome amplification method and reference genome for CNV calling, while supporting the existence of somatic CNVs in healthy and diseased human brain.

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.

DNA Instability in Neurons: Lifespan Clock and Driver of Evolution

In the last ten years, the discovery of neuronal DNA postmitotic instability has changed the theoretical landscape in neuroscience and, more broadly, biology. In 2003, A. M. Olovnikov suggested that neuronal DNA is the "initial substrate of aging". Recent experimental data have significantly increased the likelihood of this hypothesis. How does neuronal DNA accumulate damage and in what genome regions? What factors contribute to this process and how are they associated with aging and lifespan? These questions will be discussed in the review. In the course of Metazoan evolution, the instability of neuronal DNA has been accompanied by searching for the pathways to reduce the biological cost of brain activity. Various processes and activities, such as sleep, evolutionary increase in the number of neurons in the vertebrate brain, adult neurogenesis, distribution of neuronal activity, somatic polyploidy, and RNA editing in cephalopods, can be reconsidered in the light of the trade-off between neuronal plasticity and DNA instability in neurons. This topic is of considerable importance for both fundamental neuroscience and translational medicine.

A pan-tissue survey of mosaic chromosomal alterations in 948 individuals

Genetic mutations accumulate in an organism's body throughout its lifetime. While somatic single-nucleotide variants have been well characterized in the human body, the patterns and consequences of large chromosomal alterations in normal tissues remain largely unknown. Here, we present a pan-tissue survey of mosaic chromosomal alterations (mCAs) in 948 healthy individuals from the Genotype-Tissue Expression project, augmenting RNA-based allelic imbalance estimation with haplotype phasing. We found that approximately a quarter of the individuals carry a clonally-expanded mCA in at least one tissue, with incidence strongly correlated with age. The prevalence and genome-wide patterns of mCAs vary considerably across tissue types, suggesting tissue-specific mutagenic exposure and selection pressures. The mCA landscapes in normal adrenal and pituitary glands resemble those in tumors arising from these tissues, whereas the same is not true for the esophagus and skin. Together, our findings show a widespread age-dependent emergence of mCAs across normal human tissues with intricate connections to tumorigenesis.

An organism-wide ATAC-seq peak catalog for the bovine and its use to identify regulatory variants

We report the generation of an organism-wide catalog of 976,813 cis-acting regulatory elements for the bovine detected by the assay for transposase accessible chromatin using sequencing (ATAC-seq). We regroup these regulatory elements in 16 components by nonnegative matrix factorization. Correlation between the genome-wide density of peaks and transcription start sites, correlation between peak accessibility and expression of neighboring genes, and enrichment in transcription factor binding motifs support their regulatory potential. Using a previously established catalog of 12,736,643 variants, we show that the proportion of single-nucleotide polymorphisms mapping to ATAC-seq peaks is higher than expected and that this is owing to an approximately 1.3-fold higher mutation rate within peaks. Their site frequency spectrum indicates that variants in ATAC-seq peaks are subject to purifying selection. We generate eQTL data sets for liver and blood and show that variants that drive eQTL fall into liver- and blood-specific ATAC-seq peaks more often than expected by chance. We combine ATAC-seq and eQTL data to estimate that the proportion of regulatory variants mapping to ATAC-seq peaks is approximately one in three and that the proportion of variants mapping to ATAC-seq peaks that are regulatory is approximately one in 25. We discuss the implication of these findings on the utility of ATAC-seq information to improve the accuracy of genomic selection.

Toward the functional interpretation of somatic structural variations: bulk- and single-cell approaches

Structural variants (SVs) are genomic rearrangements that can take many different forms such as copy number alterations, inversions and translocations. During cell development and aging, somatic SVs accumulate in the genome with potentially neutral, deleterious or pathological effects. Generation of somatic SVs is a key mutational process in cancer development and progression. Despite their importance, the detection of somatic SVs is challenging, making them less studied than somatic single-nucleotide variants. In this review, we summarize recent advances in whole-genome sequencing (WGS)-based approaches for detecting somatic SVs at the tissue and single-cell levels and discuss their advantages and limitations. First, we describe the state-of-the-art computational algorithms for somatic SV calling using bulk WGS data and compare the performance of somatic SV detectors in the presence or absence of a matched-normal control. We then discuss the unique features of cutting-edge single-cell-based techniques for analyzing somatic SVs. The advantages and disadvantages of bulk and single-cell approaches are highlighted, along with a discussion of their sensitivity to copy-neutral SVs, usefulness for functional inferences and experimental and computational costs. Finally, computational approaches for linking somatic SVs to their functional readouts, such as those obtained from single-cell transcriptome and epigenome analyses, are illustrated, with a discussion of the promise of these approaches in health and diseases.

Comprehensive single-cell genome analysis at nucleotide resolution using the PTA Analysis Toolbox

Detection of somatic mutations in single cells has been severely hampered by technical limitations of whole-genome amplification. Novel technologies including primary template-directed amplification (PTA) significantly improved the accuracy of single-cell whole-genome sequencing (WGS) but still generate hundreds of artifacts per amplification reaction. We developed a comprehensive bioinformatic workflow, called the PTA Analysis Toolbox (PTATO), to accurately detect single base substitutions, insertions-deletions (indels), and structural variants in PTA-based WGS data. PTATO includes a machine learning approach and filtering based on recurrence to distinguish PTA artifacts from true mutations with high sensitivity (up to 90%), outperforming existing bioinformatic approaches. Using PTATO, we demonstrate that hematopoietic stem cells of patients with Fanconi anemia, which cannot be analyzed using regular WGS, have normal somatic single base substitution burdens but increased numbers of deletions. Our results show that PTATO enables studying somatic mutagenesis in the genomes of single cells with unprecedented sensitivity and accuracy.

Cannabis- and Substance-Related Carcinogenesis in Europe: A Lagged Causal Inferential Panel Regression Study

Recent European data facilitate an epidemiological investigation of the controversial cannabis-cancer relationship. Of particular concern were prior findings associating high-dose cannabis use with reproductive problems and potential genetic impacts. Cancer incidence data age-standardised to the world population was obtained from the European Cancer Information System 2000-2020 and many European national cancer registries. Drug use data were obtained from the European Monitoring Centre for Drugs and Drug Addiction. Alcohol and tobacco consumption was sourced from the WHO. Median household income was taken from the World bank. Cancer rates in high-cannabis-use countries were significantly higher than elsewhere (β-estimate = 0.4165, p = 3.54 × 10^-115). Eighteen of forty-one cancers (42,675 individual rates) were significantly associated with cannabis exposure at bivariate analysis. Twenty-five cancers were linked in inverse-probability-weighted multivariate models. Temporal lagging in panel models intensified these effects. In multivariable models, cannabis was a more powerful correlate of cancer incidence than tobacco or alcohol. Reproductive toxicity was evidenced by the involvement of testis, ovary, prostate and breast cancers and because some of the myeloid and lymphoid leukaemias implicated occur in childhood, indicating inherited intergenerational genotoxicity. Cannabis is a more important carcinogen than tobacco and alcohol and fulfills epidemiological qualitative and quantitative criteria for causality for 25/41 cancers. Reproductive and transgenerational effects are prominent. These findings confirm the clinical and epidemiological salience of cannabis as a major multigenerational community carcinogen.

High-quality nuclei isolation from postmortem human heart muscle tissues for single-cell studies

Single-cell approaches have become an increasingly popular way of understanding the genetic factors behind disease. Isolation of DNA and RNA from human tissues is necessary to analyze multi-omic data sets, providing information on the single-cell genome, transcriptome, and epigenome. Here, we isolated high-quality single-nuclei from postmortem human heart tissues for DNA and RNA analysis. Postmortem human tissues were obtained from 106 individuals, 33 with a history of myocardial disease, diabetes, or smoking, and 73 controls without heart disease. We demonstrated that the Qiagen EZ1 instrument and kit consistently isolated genomic DNA of high yield, which can be used for checking DNA quality before conducting single-cell experiments. Here, we provide a method for single-nuclei isolation from cardiac tissue, otherwise known as the SoNIC method, which allows for the isolation of single cardiomyocyte nuclei from postmortem tissue by nuclear ploidy status. We also provide a detailed quality control measure for single-nuclei whole genome amplification and a pre-amplification method for confirming genomic integrity.

In vivo Protein Footprinting Reveals the Dynamic Conformational Changes of Proteome of Multiple Tissues in Progressing Alzheimer's Disease

Numerous studies have investigated changes in protein expression at the system level using proteomic mass spectrometry, but only recently have studies explored the structure of proteins at the proteome level. We developed covalent protein painting (CPP), a protein footprinting method that quantitatively labels exposed lysine, and have now extended the method to whole intact animals to measure surface accessibility as a surrogate of in vivo protein conformations. We investigated how protein structure and protein expression change as Alzheimer's disease (AD) progresses by conducting in vivo whole animal labeling of AD mice. This allowed us to analyze broadly protein accessibility in various organs over the course of AD. We observed that structural changes of proteins related to 'energy generation,' 'carbon metabolism,' and 'metal ion homeostasis' preceded expression changes in the brain. We found that proteins in certain pathways undergoing structural changes were significantly co-regulated in the brain, kidney, muscle, and spleen.

The human brain through the lens of somatic mosaicism

Every cell in the human brain possesses a unique genome that is the product of the accumulation of somatic mutations starting from the first postzygotic cell division and continuing throughout life. Somatic mosaicism in the human brain has been the focus of several recent efforts that took advantage of key technological innovations to start elucidating brain development, aging and disease directly in human tissue. On one side, somatic mutation occurring in progenitor cells has been used as a natural barcoding system to address cell phylogenies of clone formation and cell segregation in the brain lineage. On the other side, analyses of mutation rates and patterns in the genome of brain cells have revealed mechanisms of brain aging and disorder predisposition. In addition to the study of somatic mosaicism in the normal human brain, the contribution of somatic mutation has been investigated in both developmental neuropsychiatric and neurodegenerative disorders. This review starts with a methodological perspective on the study of somatic mosaicism to then cover the most recent findings in brain development and aging, and ends with the role of somatic mutations in brain disease. Thus, this review underlies what we have learned and what is still possible to discover by looking at somatic mosaicism in the brain genome.

Single-neuron whole genome sequencing identifies increased somatic mutation burden in Alzheimer's disease related genes

Accumulation of somatic mutations in human neurons is associated with aging and neurodegeneration. To shed light on the somatic mutational burden in Alzheimer's disease (AD) neurons and get more insight into the role of somatic mutations in AD pathogenesis, we performed single-neuron whole genome sequencing to detect genome-wide somatic mutations (single nucleotide variants (SNVs) and Indels) in 96 single prefrontal cortex neurons from 8 AD patients and 8 elderly controls. We found that the mutational burden is ∼3000 somatic mutations per neuron genome in elderly subjects. AD patients have increased somatic mutation burden in AD-related annotation categories, including AD risk genes and differentially expressed genes in AD neurons. Mutational signature analysis showed somatic SNVs (sSNVs) primarily caused by aging and oxidative DNA damage processes but no significant difference was detected between AD and controls. Additionally, functional somatic mutations identified in AD patients showed significant enrichment in several AD-related pathways, including AD pathway, Notch-signaling pathway and Calcium-signaling pathway. These findings provide genetic insights into how somatic mutations may alter the function of single neurons and exert their potential roles in the pathogenesis of AD.

High-Quality Nuclei Isolation from Postmortem Human Heart Muscle Tissues for Single-Cell Studies

Single-cell approaches have become an increasingly popular way of understanding the genetic factors behind disease. Isolation of DNA and RNA from human tissues is necessary to analyze multi-omic data sets, providing information on the single-cell genome, transcriptome, and epigenome. Here, we isolated high-quality single-nuclei from postmortem human heart tissues for DNA and RNA analysis. Postmortem human tissues were obtained from 106 individuals, 33 with a history of myocardial disease, diabetes, or smoking, and 73 controls without heart disease. We demonstrated that the Qiagen EZ1 instrument and kit consistently isolated genomic DNA of high yield, which can be used for checking DNA quality before conducting single-cell experiments. Here, we provide a method for single-nuclei isolation from cardiac tissue, otherwise known as the SoNIC method, which allows for the isolation of single cardiomyocyte nuclei from postmortem tissue by nuclear ploidy status. We also provide a detailed quality control measure for single-nuclei whole genome amplification and a pre-amplification method for confirming genomic integrity.

Contrasting patterns of somatic mutations in neurons and glia reveal differential predisposition to disease in the aging human brain

Characterizing the mechanisms of somatic mutations in the brain is important for understanding aging and disease, but little is known about the mutational patterns of different cell types. We performed whole-genome sequencing of 71 oligodendrocytes and 51 neurons from neurotypical individuals (0.4 to 104 years old) and identified >67,000 somatic single nucleotide variants (sSNVs) and small insertions and deletions (indels). While both cell types accumulate mutations with age, oligodendrocytes accumulate sSNVs 69% faster than neurons (27/year versus 16/year) whereas indels accumulate 42% slower (1.8/year versus 3.1/year). Correlation with single-cell RNA and chromatin accessibility from the same brains revealed that oligodendrocyte mutations are enriched in inactive genomic regions and are distributed similarly to mutations in brain cancers. In contrast, neuronal mutations are enriched in open, transcriptionally active chromatin. These patterns highlight differences in the mutagenic processes in glia and neurons and suggest cell type-specific, age-related contributions to neurodegeneration and oncogenesis.

Active DNA demethylation promotes cell fate specification and the DNA damage response

Neurons harbor high levels of single-strand DNA breaks (SSBs) that are targeted to neuronal enhancers, but the source of this endogenous damage remains unclear. Using two systems of postmitotic lineage specification-induced pluripotent stem cell-derived neurons and transdifferentiated macrophages-we show that thymidine DNA glycosylase (TDG)-driven excision of methylcytosines oxidized with ten-eleven translocation enzymes (TET) is a source of SSBs. Although macrophage differentiation favors short-patch base excision repair to fill in single-nucleotide gaps, neurons also frequently use the long-patch subpathway. Disrupting this gap-filling process using anti-neoplastic cytosine analogs triggers a DNA damage response and neuronal cell death, which is dependent on TDG. Thus, TET-mediated active DNA demethylation promotes endogenous DNA damage, a process that normally safeguards cell identity but can also provoke neurotoxicity after anticancer treatments.

Active DNA demethylation damages DNA

Active DNA demethylation maintains enhancer activity in nonproliferating cells but can damage DNA.