Somatic mosaicism reveals clonal distributions of neocortical development

A personalized multi-platform assessment of somatic mosaicism in the human frontal cortex

Somatic mutations in individual cells lead to genomic mosaicism, contributing to the intricate regulatory landscape of genetic disorders and cancers. To evaluate and refine the detection of somatic mosaicism across different technologies with personalized donor-specific assembly (DSA), we obtained tissue from the dorsolateral prefrontal cortex (DLPFC) of a post-mortem neurotypical 31-year-old individual. We sequenced bulk DLPFC tissue using Oxford Nanopore Technologies (~60X), NovaSeq (~30X), and linked-read sequencing (~28X). Additionally, we applied Cas9 capture methodology coupled with long-read sequencing (TEnCATS), targeting active transposable elements. We also isolated and amplified DNA from flow-sorted single DLPFC neurons using MALBAC, sequencing 115 of these MALBAC libraries on Nanopore and 94 on NovaSeq. We constructed a haplotype-resolved assembly with a total length of 5.77 Gb and a phase block length of 2.67 Mb (N50) to facilitate cross-platform analysis of somatic genetic variations. We observed an increase in the phasing rate from 11.6% to 38.0% between short-read and long-read technologies. By generating a catalog of phased germline SNVs, CNVs, and TEs from the assembled genome, we applied standard approaches to recall these variants across sequencing technologies. We achieved aggregated recall rates from 97.3% to 99.4% based on long-read bulk tissue data, setting an upper bound for detection limits. Moreover, utilizing haplotype-based analysis from DSA, we achieved a remarkable reduction in false positive somatic calls in bulk tissue, ranging from 14.9% to 72.4%. We developed pipelines leveraging DSA information to enhance somatic large genetic variant calling in long-read single cells. By examining somatic variation using long-reads in 115 individual neurons, we identified 468 candidate somatic heterozygous large deletions (1.5Mb - 20Mb), 137 of which intersected with short-read single-cell data. Additionally, we identified 61 putative somatic TEs (60 Alus, one LINE-1) in the single-cell data. Collectively, our analysis spans personalized assembly to single-cell somatic variant calling, providing a comprehensive ab initio ad finem approach and resource in real human tissue.

Multiplexing cortical brain organoids for the longitudinal dissection of developmental traits at single-cell resolution

Dissecting human neurobiology at high resolution and with mechanistic precision requires a major leap in scalability, given the need for experimental designs that include multiple individuals and, prospectively, population cohorts. To lay the foundation for this, we have developed and benchmarked complementary strategies to multiplex brain organoids by pooling cells from different pluripotent stem cell (PSC) lines either during organoid generation (mosaic models) or before single-cell RNA sequencing (scRNA-seq) library preparation (downstream multiplexing). We have also developed a new computational method, SCanSNP, and a consensus call to deconvolve cell identities, overcoming current criticalities in doublets and low-quality cell identification. We validated both multiplexing methods for charting neurodevelopmental trajectories at high resolution, thus linking specific individuals' trajectories to genetic variation. Finally, we modeled their scalability across different multiplexing combinations and showed that mosaic organoids represent an enabling method for high-throughput settings. Together, this multiplexing suite of experimental and computational methods provides a highly scalable resource for brain disease and neurodiversity modeling.

Human Brain Ancestral Barcodes

Dynamic CpG methylation "barcodes" were read from 15,000 to 21,000 single cells from three human male brains. To overcome sparse sequencing coverage, the barcode had ~31,000 rapidly fluctuating X-chromosome CpG sites (fCpGs), with at least 500 covered sites per cell and at least 30 common sites between cell pairs (average of ~48). Barcodes appear to start methylated and record mitotic ages because excitatory neurons and glial cells that emerge later in development were less methylated. Barcodes are different between most cells, with average pairwise differences (PWDs) of ~0.5 between cells. About 10 cell pairs per million were more closely related with PWDs < 0.05. Barcodes appear to record ancestry and reconstruct trees where more related cells had similar phenotypes, albeit some pairs had phenotypic differences. Inhibitory neurons showed more evidence of tangential migration than excitatory neurons, with related cells in different cortical regions. fCpG barcodes become polymorphic during development and can distinguish between thousands of human cells.

Genomic mosaicism reveals developmental organization of trunk neural crest-derived ganglia

The neural crest generates numerous cell types, but conflicting results leave developmental origins unresolved. Here using somatic mosaic variants as cellular barcodes, we infer embryonic clonal dynamics of trunk neural crest, focusing on the sensory and sympathetic ganglia. From three independent adult neurotypical human donors, we identified 1,278 mosaic variants using deep whole-genome sequencing, then profiled allelic fractions in 187 anatomically dissected ganglia. We found a massive rostrocaudal spread of progenitor clones specific to sensory or sympathetic ganglia, which unlike in the brain, showed robust bilateral distributions. Computational modeling suggested neural crest progenitor fate specification preceded delamination from neural tube. Single-cell multiomic analysis suggested both neurons and glia contributed to the rostrocaudal clonal organization. CRISPR barcoding in mice and live imaging in quail embryos confirmed these clonal dynamics across multiple somite levels. Our findings reveal an evolutionarily conserved clonal spread of cells populating peripheral neural ganglia.

Threshold of somatic mosaicism leading to brain dysfunction with focal epilepsy

Somatic mosaicism in a fraction of brain cells causes neurodevelopmental disorders, including childhood intractable epilepsy. However, the threshold for somatic mosaicism leading to brain dysfunction is unknown. In this study, we induced various mosaic burdens in focal cortical dysplasia type II (FCD II) mice, featuring mTOR somatic mosaicism and spontaneous behavioural seizures. The mosaic burdens ranged from approximately 1000 to 40 000 neurons expressing the mTOR mutant in the somatosensory or medial prefrontal cortex. Surprisingly, approximately 8000-9000 neurons expressing the MTOR mutant, extrapolated to constitute 0.08%-0.09% of total cells or roughly 0.04% of variant allele frequency in the mouse hemicortex, were sufficient to trigger epileptic seizures. The mutational burden was correlated with seizure frequency and onset, with a higher tendency for electrographic inter-ictal spikes and beta- and gamma-frequency oscillations in FCD II mice exceeding the threshold. Moreover, mutation-negative FCD II patients in deep sequencing of their bulky brain tissues revealed somatic mosaicism of the mTOR pathway genes as low as 0.07% in resected brain tissues through ultra-deep targeted sequencing (up to 20 million reads). Thus, our study suggests that extremely low levels of somatic mosaicism can contribute to brain dysfunction.

Building the brain mosaic: an expanded view

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

Increased somatic mosaicism in autosomal and X chromosomes for suicide death

Mosaic chromosomal alterations (mCAs) are classified as mosaic deletions (loss), copy-neutral loss of heterozygosity (CN-LOH), and duplications (gain), attracting special attention as biological aging-related acquired genetic alterations. While these mCAs have been linked with aging and various diseases, no study has investigated their association with suicide risk which is associated with abnormal biological aging. Here, we examined the association between suicide deaths and mCAs, including mosaic loss of the X (mLOX) and Y chromosomes, by leveraging blood-derived single nucleotide polymorphism-array data. The first (410 suicide decedents and 88,870 controls) and the second (363 suicide decedents and 88,870 controls) cohorts were analyzed and integrated using meta-analyses (773 suicide decedents and 177,740 controls). Total mCAs in autosomal chromosomes were significantly increased in suicide (p = 1.28 × 10-6, odds ratio [OR] = 1.78), mostly driven by loss (p = 4.05 × 10-9, OR = 2.70) and gain (p = 1.08 × 10-3, OR = 2.23). mLOX were significantly increased in female suicide (p = 2.66 × 10-21, OR = 4.00). The directions of effects of all mCAs in autosomal and sex chromosomes on suicide were the same in the first and second sets. Subgroup analyses suggest that our findings were mostly driven by suicide itself, and not confounded by comorbid psychiatric disorders or physical diseases, smoking status, sample location, or postmortem sample status. In conclusion, we provide the first evidence for aberrant mCAs in somatic autosomal and X chromosomes in suicide, which may contribute to an improved understanding of the genomic pathophysiology underlying suicide.

Developmental signals control chromosome segregation fidelity during pluripotency and neurogenesis by modulating replicative stress

Human development relies on the correct replication, maintenance and segregation of our genetic blueprints. How these processes are monitored across embryonic lineages, and why genomic mosaicism varies during development remain unknown. Using pluripotent stem cells, we identify that several patterning signals-including WNT, BMP, and FGF-converge into the modulation of DNA replication stress and damage during S-phase, which in turn controls chromosome segregation fidelity in mitosis. We show that the WNT and BMP signals protect from excessive origin firing, DNA damage and chromosome missegregation derived from stalled forks in pluripotency. Cell signalling control of chromosome segregation declines during lineage specification into the three germ layers, but re-emerges in neural progenitors. In particular, we find that the neurogenic factor FGF2 induces DNA replication stress-mediated chromosome missegregation during the onset of neurogenesis, which could provide a rationale for the elevated chromosomal mosaicism of the developing brain. Our results highlight roles for morphogens and cellular identity in genome maintenance that contribute to somatic mosaicism during mammalian development.

Genome-wide detection of somatic mosaicism at short tandem repeats

MOTIVATION

Somatic mosaicism has been implicated in several developmental disorders, cancers, and other diseases. Short tandem repeats (STRs) consist of repeated sequences of 1-6 bp and comprise >1 million loci in the human genome. Somatic mosaicism at STRs is known to play a key role in the pathogenicity of loci implicated in repeat expansion disorders and is highly prevalent in cancers exhibiting microsatellite instability. While a variety of tools have been developed to genotype germline variation at STRs, a method for systematically identifying mosaic STRs is lacking.

RESULTS

We introduce prancSTR, a novel method for detecting mosaic STRs from individual high-throughput sequencing datasets. prancSTR is designed to detect loci characterized by a single high-frequency mosaic allele, but can also detect loci with multiple mosaic alleles. Unlike many existing mosaicism detection methods for other variant types, prancSTR does not require a matched control sample as input. We show that prancSTR accurately identifies mosaic STRs in simulated data, demonstrate its feasibility by identifying candidate mosaic STRs in Illumina whole genome sequencing data derived from lymphoblastoid cell lines for individuals sequenced by the 1000 Genomes Project, and evaluate the use of prancSTR on Element and PacBio data. In addition to prancSTR, we present simTR, a novel simulation framework which simulates raw sequencing reads with realistic error profiles at STRs.

AVAILABILITY AND IMPLEMENTATION

prancSTR and simTR are freely available at https://github.com/gymrek-lab/trtools. Detailed documentation is available at https://trtools.readthedocs.io/.

Somatic nuclear mitochondrial DNA insertions are prevalent in the human brain and accumulate over time in fibroblasts

The transfer of mitochondrial DNA into the nuclear genomes of eukaryotes (Numts) has been linked to lifespan in nonhuman species and recently demonstrated to occur in rare instances from one human generation to the next. Here, we investigated numtogenesis dynamics in humans in 2 ways. First, we quantified Numts in 1,187 postmortem brain and blood samples from different individuals. Compared to circulating immune cells (n = 389), postmitotic brain tissue (n = 798) contained more Numts, consistent with their potential somatic accumulation. Within brain samples, we observed a 5.5-fold enrichment of somatic Numt insertions in the dorsolateral prefrontal cortex (DLPFC) compared to cerebellum samples, suggesting that brain Numts arose spontaneously during development or across the lifespan. Moreover, an increase in the number of brain Numts was linked to earlier mortality. The brains of individuals with no cognitive impairment (NCI) who died at younger ages carried approximately 2 more Numts per decade of life lost than those who lived longer. Second, we tested the dynamic transfer of Numts using a repeated-measures whole-genome sequencing design in a human fibroblast model that recapitulates several molecular hallmarks of aging. These longitudinal experiments revealed a gradual accumulation of 1 Numt every ~13 days. Numtogenesis was independent of large-scale genomic instability and unlikely driven by cell clonality. Targeted pharmacological perturbations including chronic glucocorticoid signaling or impairing mitochondrial oxidative phosphorylation (OxPhos) only modestly increased the rate of numtogenesis, whereas patient-derived SURF1-mutant cells exhibiting mtDNA instability accumulated Numts 4.7-fold faster than healthy donors. Combined, our data document spontaneous numtogenesis in human cells and demonstrate an association between brain cortical somatic Numts and human lifespan. These findings open the possibility that mito-nuclear horizontal gene transfer among human postmitotic tissues produces functionally relevant human Numts over timescales shorter than previously assumed.

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.

Cell-type-resolved mosaicism reveals clonal dynamics of the human forebrain

Debate remains around the anatomical origins of specific brain cell subtypes and lineage relationships within the human forebrain1-7. Thus, direct observation in the mature human brain is critical for a complete understanding of its structural organization and cellular origins. Here we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific mosaic variant barcode analysis. From four hemispheres and two different human neurotypical donors, we identified 287 and 780 mosaic variants, respectively, that were used to deconvolve clonal dynamics. Clonal spread and allele fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted than resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome transcriptome analysis at both a cell-type-specific and a single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of mosaic variants across 17 locations within one parietal lobe reveals that restriction of clonal spread in the anterior-posterior axis precedes restriction in the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus, cell-type-resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Cerebral organoids display dynamic clonal growth and tunable tissue replenishment

During brain development, neural progenitors expand through symmetric divisions before giving rise to differentiating cell types via asymmetric divisions. Transition between those modes varies among individual neural stem cells, resulting in clones of different sizes. Imaging-based lineage tracing allows for lineage analysis at high cellular resolution but systematic approaches to analyse clonal behaviour of entire tissues are currently lacking. Here we implement whole-tissue lineage tracing by genomic DNA barcoding in 3D human cerebral organoids, to show that individual stem cell clones produce progeny on a vastly variable scale. By using stochastic modelling we find that variable lineage sizes arise because a subpopulation of lineages retains symmetrically dividing cells. We show that lineage sizes can adjust to tissue demands after growth perturbation via chemical ablation or genetic restriction of a subset of cells in chimeric organoids. Our data suggest that adaptive plasticity of stem cell populations ensures robustness of development in human brain organoids.

Chromosomal single-strand break repair and neurological disease: Implications on transcription and emerging genomic tools

Cells are constantly exposed to various sources of DNA damage that pose a threat to their genomic integrity. One of the most common types of DNA breaks are single-strand breaks (SSBs). Mutations in the repair proteins that are important for repairing SSBs have been reported in several neurological disorders. While several tools have been utilised to investigate SSBs in cells, it was only through recent advances in genomics that we are now beginning to understand the architecture of the non-random distribution of SSBs and their impact on key cellular processes such as transcription and epigenetic remodelling. Here, we discuss our current understanding of the genome-wide distribution of SSBs, their link to neurological disorders and summarise recent technologies to investigate SSBs at the genomic level.

Brain injury drives optic glioma formation through neuron-glia signaling

Tissue injury and tumorigenesis share many cellular and molecular features, including immune cell (T cells, monocytes) infiltration and inflammatory factor (cytokines, chemokines) elaboration. Their common pathobiology raises the intriguing possibility that brain injury could create a tissue microenvironment permissive for tumor formation. Leveraging several murine models of the Neurofibromatosis type 1 (NF1) cancer predisposition syndrome and two experimental methods of brain injury, we demonstrate that both optic nerve crush and diffuse traumatic brain injury induce optic glioma (OPG) formation in mice harboring Nf1-deficient preneoplastic progenitors. We further elucidate the underlying molecular and cellular mechanisms, whereby glutamate released from damaged neurons stimulates IL-1β release by oligodendrocytes to induce microglia expression of Ccl5, a growth factor critical for Nf1-OPG formation. Interruption of this cellular circuit using glutamate receptor, IL-1β or Ccl5 inhibitors abrogates injury-induced glioma progression, thus establishing a causative relationship between injury and tumorigenesis.

Genome-wide detection of somatic mosaicism at short tandem repeats

Motivation

Somatic mosaicism, in which a mutation occurs post-zygotically, has been implicated in several developmental disorders, cancers, and other diseases. Short tandem repeats (STRs) consist of repeated sequences of 1-6bp and comprise more than 1 million loci in the human genome. Somatic mosaicism at STRs is known to play a key role in the pathogenicity of loci implicated in repeat expansion disorders and is highly prevalent in cancers exhibiting microsatellite instability. While a variety of tools have been developed to genotype germline variation at STRs, a method for systematically identifying mosaic STRs (mSTRs) is lacking.

Results

We introduce prancSTR, a novel method for detecting mSTRs from individual high-throughput sequencing datasets. Unlike many existing mosaicism detection methods for other variant types, prancSTR does not require a matched control sample as input. We show that prancSTR accurately identifies mSTRs in simulated data and demonstrate its feasibility by identifying candidate mSTRs in whole genome sequencing (WGS) data derived from lymphoblastoid cell lines for individuals sequenced by the 1000 Genomes Project. Our analysis identified an average of 76 and 577 non-homopolymer and homopolymer mSTRs respectively per cell line as well as multiple cell lines with outlier mSTR counts more than 6 times the population average, suggesting a subset of cell lines have particularly high STR instability rates.

Availability

prancSTR is freely available at https://github.com/gymrek-lab/trtools.

Documentation

Detailed documentation is available at https://trtools.readthedocs.io/.

Genomic data resources of the Brain Somatic Mosaicism Network for neuropsychiatric diseases

Somatic mosaicism is defined as an occurrence of two or more populations of cells having genomic sequences differing at given loci in an individual who is derived from a single zygote. It is a characteristic of multicellular organisms that plays a crucial role in normal development and disease. To study the nature and extent of somatic mosaicism in autism spectrum disorder, bipolar disorder, focal cortical dysplasia, schizophrenia, and Tourette syndrome, a multi-institutional consortium called the Brain Somatic Mosaicism Network (BSMN) was formed through the National Institute of Mental Health (NIMH). In addition to genomic data of affected and neurotypical brains, the BSMN also developed and validated a best practices somatic single nucleotide variant calling workflow through the analysis of reference brain tissue. These resources, which include >400 terabytes of data from 1087 subjects, are now available to the research community via the NIMH Data Archive (NDA) and are described here.

Cell lineage analysis with somatic mutations reveals late divergence of neuronal cell types and cortical areas in human cerebral cortex

The mammalian cerebral cortex shows functional specialization into regions with distinct neuronal compositions, most strikingly in the human brain, but little is known in about how cellular lineages shape cortical regional variation and neuronal cell types during development. Here, we use somatic single nucleotide variants (sSNVs) to map lineages of neuronal sub-types and cortical regions. Early-occurring sSNVs rarely respect Brodmann area (BA) borders, while late-occurring sSNVs mark neuron-generating clones with modest regional restriction, though descendants often dispersed into neighboring BAs. Nevertheless, in visual cortex, BA17 contains 30-70% more sSNVs compared to the neighboring BA18, with clones across the BA17/18 border distributed asymmetrically and thus displaying different cortex-wide dispersion patterns. Moreover, we find that excitatory neuron-generating clones with modest regional restriction consistently share low-mosaic sSNVs with some inhibitory neurons, suggesting significant co-generation of excitatory and some inhibitory neurons in the dorsal cortex. Our analysis reveals human-specific cortical cell lineage patterns, with both regional inhomogeneities in progenitor proliferation and late divergence of excitatory/inhibitory lineages.

Chromosome instability and aneuploidy in the mammalian brain

This review investigates the role of aneuploidy and chromosome instability (CIN) in the aging brain. Aneuploidy refers to an abnormal chromosomal count, deviating from the normal diploid set. It can manifest as either a deficiency or excess of chromosomes. CIN encompasses a broader range of chromosomal alterations, including aneuploidy as well as structural modifications in DNA. We provide an overview of the state-of-the-art methodologies utilized for studying aneuploidy and CIN in non-tumor somatic tissues devoid of clonally expanded populations of aneuploid cells.CIN and aneuploidy, well-established hallmarks of cancer cells, are also associated with the aging process. In non-transformed cells, aneuploidy can contribute to functional impairment and developmental disorders. Despite the importance of understanding the prevalence and specific consequences of aneuploidy and CIN in the aging brain, these aspects remain incompletely understood, emphasizing the need for further scientific investigations.This comprehensive review consolidates the present understanding, addresses discrepancies in the literature, and provides valuable insights for future research efforts.

Cell-type-resolved somatic mosaicism reveals clonal dynamics of the human forebrain

Debate remains around anatomic origins of specific brain cell subtypes and lineage relationships within the human forebrain. Thus, direct observation in the mature human brain is critical for a complete understanding of the structural organization and cellular origins. Here, we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific Mosaic Variant Barcode Analysis. From four hemispheres from two different human neurotypical donors, we identified 287 and 780 mosaic variants (MVs), respectively that were used to deconvolve clonal dynamics. Clonal spread and allelic fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted compared with resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome-transcriptome analysis at both a cell-type-specific and single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of MVs across 17 locations within one parietal lobe reveals restrictions of clonal spread in the anterior-posterior axis precedes that of the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus cell-type resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Control-independent mosaic single nucleotide variant detection with DeepMosaic

Mosaic variants (MVs) reflect mutagenic processes during embryonic development and environmental exposure, accumulate with aging and underlie diseases such as cancer and autism. The detection of noncancer MVs has been computationally challenging due to the sparse representation of nonclonally expanded MVs. Here we present DeepMosaic, combining an image-based visualization module for single nucleotide MVs and a convolutional neural network-based classification module for control-independent MV detection. DeepMosaic was trained on 180,000 simulated or experimentally assessed MVs, and was benchmarked on 619,740 simulated MVs and 530 independent biologically tested MVs from 16 genomes and 181 exomes. DeepMosaic achieved higher accuracy compared with existing methods on biological data, with a sensitivity of 0.78, specificity of 0.83 and positive predictive value of 0.96 on noncancer whole-genome sequencing data, as well as doubling the validation rate over previous best-practice methods on noncancer whole-exome sequencing data (0.43 versus 0.18). DeepMosaic represents an accurate MV classifier for noncancer samples that can be implemented as an alternative or complement to existing methods.

svCapture: efficient and specific detection of very low frequency structural variant junctions by error-minimized capture sequencing

Error-corrected sequencing of genomic targets enriched by probe-based capture has become a standard approach for detecting single-nucleotide variants (SNVs) and small insertion/deletions (indels) present at very low variant allele frequencies. Less attention has been given to comparable strategies for rare structural variant (SV) junctions, where different error mechanisms must be addressed. Working from samples with known SV properties, we demonstrate that duplex sequencing (DuplexSeq), which demands confirmation of variants on both strands of a source DNA molecule, eliminates false SV junctions arising from chimeric PCR. DuplexSeq could not address frequent intermolecular ligation artifacts that arise during Y-adapter addition prior to strand denaturation without requiring multiple source molecules. In contrast, tagmentation libraries coupled with data filtering based on strand family size greatly reduced both artifact classes and enabled efficient and specific detection of single-molecule SV junctions. The throughput of SV capture sequencing (svCapture) and base-level accuracy of DuplexSeq provided detailed views of the microhomology profile and limited occurrence of de novo SNVs near the junctions of hundreds of newly created SVs, suggesting end joining as a possible formation mechanism. The open source svCapture pipeline enables rare SV detection as a routine addition to SNVs/indels in properly prepared capture sequencing libraries.

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.

Somatic nuclear mitochondrial DNA insertions are prevalent in the human brain and accumulate over time in fibroblasts

The transfer of mitochondrial DNA into the nuclear genomes of eukaryotes (Numts) has been linked to lifespan in non-human species 1-3 and recently demonstrated to occur in rare instances from one human generation to the next 4. Here we investigated numtogenesis dynamics in humans in two ways. First, we quantified Numts in 1,187 post-mortem brain and blood samples from different individuals. Compared to circulating immune cells (n=389), post-mitotic brain tissue (n=798) contained more Numts, consistent with their potential somatic accumulation. Within brain samples we observed a 5.5-fold enrichment of somatic Numt insertions in the dorsolateral prefrontal cortex compared to cerebellum samples, suggesting that brain Numts arose spontaneously during development or across the lifespan. Moreover, more brain Numts was linked to earlier mortality. The brains of individuals with no cognitive impairment who died at younger ages carried approximately 2 more Numts per decade of life lost than those who lived longer. Second, we tested the dynamic transfer of Numts using a repeated-measures WGS design in a human fibroblast model that recapitulates several molecular hallmarks of aging 5. These longitudinal experiments revealed a gradual accumulation of one Numt every ~13 days. Numtogenesis was independent of large-scale genomic instability and unlikely driven cell clonality. Targeted pharmacological perturbations including chronic glucocorticoid signaling or impairing mitochondrial oxidative phosphorylation (OxPhos) only modestly increased the rate of numtogenesis, whereas patient-derived SURF1-mutant cells exhibiting mtDNA instability accumulated Numts 4.7-fold faster than healthy donors. Combined, our data document spontaneous numtogenesis in human cells and demonstrate an association between brain cortical somatic Numts and human lifespan. These findings open the possibility that mito-nuclear horizontal gene transfer among human post-mitotic tissues produce functionally-relevant human Numts over timescales shorter than previously assumed.

A Chinese SCA36 pedigree analysis of NOP56 expansion region based on long-read sequencing

Introduction: Spinocerebellar ataxias 36 (SCA36) is the neurodegenerative disease caused by the GGCCTG Hexanucleotide repeat expansions in NOP56, which is too long to sequence using short-read sequencing. Single molecule real time (SMRT) sequencing can sequence across disease-causing repeat expansion. We report the first long-read sequencing data across the expansion region in SCA36.

Methods: We collected and described the clinical manifestations and imaging features of Han Chinese pedigree with three generations of SCA36. Also, we focused on structural variation analysis for intron 1 of the NOP56 gene by SMRT sequencing in the assembled genome.

Results: The main clinical features of this pedigree are late-onset ataxia symptoms, with a presymptomatic presence of affective and sleep disorders. In addition, the results of SMRT sequencing showed the specific repeat expansion region and demonstrated that the region was not composed of single GGCCTG hexanucleotides and there were random interruptions.

Discussion: We extended the phenotypic spectrum of SCA36. We applied SMRT sequencing to reveal the correlation between genotype and phenotype of SCA36. Our findings indicated that long-read sequencing is well suited to characterize known repeat expansion.

Mapping the Complex Genetic Landscape of Human Neurons

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

Comprehensive multi-omic profiling of somatic mutations in malformations of cortical development

Malformations of cortical development (MCD) are neurological conditions involving focal disruptions of cortical architecture and cellular organization that arise during embryogenesis, largely from somatic mosaic mutations, and cause intractable epilepsy. Identifying the genetic causes of MCD has been a challenge, as mutations remain at low allelic fractions in brain tissue resected to treat condition-related epilepsy. Here we report a genetic landscape from 283 brain resections, identifying 69 mutated genes through intensive profiling of somatic mutations, combining whole-exome and targeted-amplicon sequencing with functional validation including in utero electroporation of mice and single-nucleus RNA sequencing. Genotype-phenotype correlation analysis elucidated specific MCD gene sets associated with distinct pathophysiological and clinical phenotypes. The unique single-cell level spatiotemporal expression patterns of mutated genes in control and patient brains indicate critical roles in excitatory neurogenic pools during brain development and in promoting neuronal hyperexcitability after birth.

Grave-to-cradle: human embryonic lineage tracing from the postmortem body

Curiosity concerning the process of human creation has been around for a long time. Relevant questions seemed to be resolved with the knowledge of how cells divide after fertilization obtained through in vitro fertilization experiments. However, we still do not know how human life is created at the cellular level. Recently, the value of cadavers as a resource from which to obtain "normal" cells and tissues has been established, and human research using postmortem bodies has attracted growing scientific attention. As the human genome can be analyzed at the level of nucleotides through whole-genome sequencing, individual cells in a postmortem body can be traced back to determine what developmental processes have transpired from fertilization. These retrospective lineage tracing studies have answered several unsolved questions on how humans are created. This review covers the methodologies utilized in lineage tracing research in a historical context and the conceptual basis for reconstructing the division history of cells in a retrospective manner using postzygotic somatic variants in postmortem tissue. We further highlight answers that postmortem research could potentially address and discuss issues that wait to be solved in the future.

PCDH19-related epilepsy in mosaic males: The phenotypic implication of genotype and variant allele frequency

OBJECTIVE

To analyze the genotypes and phenotypes of mosaic male patients with PCDH19-related epilepsy (PCDH19-RE) and explore the correlation between genotype, variant allele frequency (VAF), and phenotypic severity.

METHODS

Clinical data and peripheral blood samples of 11 male mosaic patients were collected and analyzed in our study. The VAF of the PCDH19 gene from peripheral blood was quantified using amplicon-based deep sequencing. Additional 20 mosaic male patients with PCDH19-RE were collected from the published literature, with 10 patients whose VAFs of the PCDH19 gene were available for analytic purposes.

RESULTS

In our cohort of 11 patients, 10 variants were identified, and four were novel. The VAF of the PCDH19 gene from peripheral blood ranged from 27 to 90%. The median seizure onset age was 6 months (range: 4-9 months). Clinical manifestations included cluster seizures (100%), fever sensitivity (73%), focal seizures (91%), developmental delay/intellectual disability (DD/ID, 82%), and autistic features (45%). Thirty-one mosaic male patients collected from our cohort and the literature developed seizures mostly (87%) within one year of age. Variant types included missense variants (42%), truncating variants (52%), splice variants (3%), and whole PCDH19 deletion (3%). Among 21 patients with a definite VAF from our cohort and the literature, nine had a low VAF ( ≤ 50%) and 12 had a high VAF (> 50%). Seventy-five percent of variants from the high VAF group were missense, whereas 89% of those from the low VAF group were truncations. The median seizure onset age was 6 months in the low VAF group and 9 months in the high VAF group (p = 0.018). Forty-four percent (4/9) of patients from the low VAF group achieved seizure-free for ≥1 year, whereas none of the 12 patients from the high VAF group did (p = 0.021). DD/ID was present in 83% (10/12) of the high VAF group and 56% (5/9) of the low VAF group (p = 0.331).

CONCLUSION

The predominant variant types were truncating and missense variants. Missense variants tended to have higher VAFs. Patients with a high VAF were more likely to have a more severe epileptic phenotype. Our findings shed light on the phenotypic implications of VAF in mosaic males with PCDH19-RE.

Unbiased mosaic variant assessment in sperm: a cohort study to test predictability of transmission

Background:

De novo mutations underlie individually rare but collectively common pediatric congenital disorders. Some of these mutations can also be detected in tissues and from cells in a parent, where their abundance and tissue distribution can be measured. We previously reported that a subset of these mutations is detectable in sperm from the father, predicted to impact the health of offspring.

Methods:

As a cohort study, in three independent couples undergoing in vitro fertilization, we first identified male gonadal mosaicism through deep whole genome sequencing. We then confirmed variants and assessed their transmission to preimplantation blastocysts (32 total) through targeted ultra-deep genotyping.

Results:

Across 55 gonadal mosaic variants, 15 were transmitted to blastocysts for a total of 19 transmission events. This represented an overall predictable but slight undertransmission based upon the measured mutational abundance in sperm. We replicated this conclusion in an independent, previously published family-based cohort.

Conclusions:

Unbiased preimplantation genetic testing for gonadal mosaicism may represent a feasible approach to reduce the transmission of potentially harmful de novo mutations. This—in turn—could help to reduce their impact on miscarriages and pediatric disease.

Funding:

No external funding was received for this work.

Whole-Genome Amplification—Surveying Yield, Reproducibility, and Heterozygous Balance, Reported by STR-Targeting MIPs

Whole-genome amplification is a crucial first step in nearly all single-cell genomic analyses, with the following steps focused on its products. Bias and variance caused by the whole-genome amplification process add numerous challenges to the world of single-cell genomics. Short tandem repeats are sensitive genomic markers used widely in population genetics, forensics, and retrospective lineage tracing. A previous evaluation of common whole-genome amplification targeting ~1000 non-autosomal short tandem repeat loci is extended here to ~12,000 loci across the entire genome via duplex molecular inversion probes. Other than its improved scale and reduced noise, this system detects an abundance of heterogeneous short tandem repeat loci, allowing the allelic balance to be reported. We show here that while the best overall yield is obtained using RepliG-SC, the maximum uniformity between alleles and reproducibility across cells are maximized by Ampli1, rendering it the best candidate for the comparative heterozygous analysis of single-cell genomes.