Ester Kalef-Ezra

Somatic mutations arise throughout life, with functional consequences tied to the cell populations in which they occur. Genome-wide studies measure somatic variations in bulk tissue, whereas single-cell approaches resolve cell identity but provide limited sensitivity for complex alleles. Here we developed SniffCell, which uses DNA methylation carried on native long reads to assign somatic variant-supporting molecules to methylation-resolvable cell types. SniffCell builds cell-type-discriminatory methylation signatures across eight tissues, assigns long reads to cell types, and provides cell-type-specific variant calling. Across peripheral blood mononuclear cells and brain benchmarks, SniffCell recovered sorted cell identities and validated cell-type-specific variant assignments using purified immune-cell, neuronal, and oligodendrocyte fractions. In blood, SniffCell recovered lineage-restricted antigen receptor rearrangements and localized a somatic tandem-repeat expansion to T cells. In the frontal cortex, SniffCell identified recurrent neuron-specific tandem-repeat expansions in genes including FGF14, LRRC7 and SH3RF3 . Across three brain cohorts comprising 172 donors, recurrent neuron-associated expansions were enriched for GAA-rich motifs. In donors with matched blood, and diverged more strongly from the inherited repeat length, whereas oligodendrocyte-associated alleles more often tracked it. SniffCell transforms native bulk long-read genomes into a cell-type-aware resource for somatic variant discovery and reveals recurrent somatic instability in human tissues at cell-type resolution.

Detecting low variant allele fraction (VAF) mosaic variants without matching controls remains a major challenge in genomics, limited by technical noise, lack of benchmarks, and computational scalability. We present the DRAGEN mosaic caller, a hardware-accelerated approach identifying variants down to ∼1-2% VAF with low false-positive rates and hour-scale runtimes for mosaic SNV/indel detection from bulk sequencing. To support evaluation, we introduce a genome-wide low-VAF benchmark for variants between 1-10% VAF. Application to blood, sperm, and brain tissues revealed patterns, including mosaic hotspots and mutational signatures. The first analysis of HG002 blood showed that many “mosaic” variants defined from HG002 cell lines are likely culture-derived and not in vivo mutations. Importantly, DRAGEN also enables personalized assembly pangenome references to improve alignment and mosaic variant detection in complex regions. Together, this development makes routine low-VAF discovery feasible, opening new opportunities to study mosaic mutations in healthy and disease individuals.

The advent of single cell DNA sequencing revealed astonishing dynamics of genomic variability, but failed at characterizing smaller to mid size variants that on the germline level have a profound impact. In this work we discover previously uncharacterized genomic dynamics in 18 cells from three human brains utilizing single cell long-read whole genome sequencing. This provides key insights into the dynamic of the genomes of individual cells and further highlights brain specific activity of transposable elements, but requires validation in larger studies.

The advent of single cell DNA sequencing revealed astonishing dynamics of genomic variability, but failed at characterizing smaller to mid size variants that on the germline level have a profound impact. In this work we discover novel dynamics in three brains utilizing single cell long-read sequencing. This provides key insights into the dynamic of the genomes of individual cells and further highlights brain specific activity of transposable elements.

Using Oxford Nanopore MinION and the Rapid Barcoding library preparation kit, we perform long read sequencing on droplet multiple displacement amplification (dMDA) products generated from single nuclei. Single nuclei contain around 6.6pg of DNA, to sequence this miniscule amount of DNA using long reads a whole genome amplification (WGA) method must be used to provide enough starting material for library preparation. MDA is currently the only WGA method which produces reads >10kb in length. Sequencing DNA of this length using long reads can allow for direct resolution and identification of most structural variants, including large LINE-1 insertions and deletions. dMDA is an improvement of the MDA technique, partitioning DNA template molecules into droplets prevents the generation of intermolecular chimeras, and a finite amount of dMDA reagents being packed into each droplet, prevents overamplification of certain reads. The expected result of this protocol is to be able to produce and sequence libraries generated from dMDA products.

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5–50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2 , including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

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.

Glutamate dehydrogenase (GDH) is a hexameric enzyme that catalyzes the reversible conversion of glutamate to α-ketoglutarate and ammonia while reducing NAD(P)⁺ to NAD(P)H. It is found in all living organisms serving both catabolic and anabolic reactions. In mammalian tissues, oxidative deamination of glutamate via GDH generates α-ketoglutarate, which is metabolized by the Krebs cycle, leading to the synthesis of ATP. In addition, the GDH pathway is linked to diverse cellular processes, including ammonia metabolism, acid-base equilibrium, redox homeostasis (via formation of fumarate), lipid biosynthesis (via oxidative generation of citrate), and lactate production. While most mammals possess a single GDH1 protein (hGDH1 in the human) that is highly expressed in the liver, humans and other primates have acquired, via duplication, an hGDH2 isoenzyme with distinct functional properties and tissue expression profile. The novel hGDH2 underwent rapid evolutionary adaptation, acquiring unique properties that enable enhanced enzyme function under conditions inhibitory to its ancestor hGDH1. These are thought to provide a biological advantage to humans with hGDH2 evolution occurring concomitantly with human brain development. hGDH2 is co-expressed with hGDH1 in human brain, kidney, testis and steroidogenic organs, but not in the liver. In human cerebral cortex, hGDH1 and hGDH2 are expressed in astrocytes, the cells responsible for removing and metabolizing transmitter glutamate, and for supplying neurons with glutamine and lactate. In human testis, hGDH2 (but not hGDH1) is densely expressed in the Sertoli cells, known to provide the spermatids with lactate and other nutrients. In steroid producing cells, hGDH1/2 is thought to generate reducing equivalents (NADPH) in the mitochondria for the biosynthesis of steroidal hormones. Lastly, up-regulation of hGDH1/2 expression occurs in cancer, permitting neoplastic cells to utilize glutamine/glutamate for their growth. In addition, deregulation of hGDH1/2 is implicated in the pathogenesis of several human disorders.

Mammalian glutamate dehydrogenase (GDH), a nuclear-encoded enzyme central to cellular metabolism, is among the most abundant mitochondrial proteins (constituting up to 10% of matrix proteins). To attain such high levels, GDH depends on very efficient mitochondrial targeting that, for human isoenzymes hGDH1 and hGDH2, is mediated by an unusually long cleavable presequence (N53). Here, we studied the mitochondrial transport of these proteins using isolated yeast mitochondria and human cell lines. We found that both hGDHs were very rapidly imported and processed in isolated mitochondria, with their presequences (N53) alone being capable of directing non-mitochondrial proteins into mitochondria. These presequences were predicted to form two α helices (α1: N 1–10; α2: N 16–32) separated by loops. Selective deletion of the α1 helix abolished the mitochondrial import of hGDHs. While the α1 helix alone had a very weak hGDH mitochondrial import capacity, it could direct efficiently non-mitochondrial proteins into mitochondria. In contrast, the α2 helix had no autonomous mitochondrial-targeting capacity. A peptide consisting of α1 and α2 helices without intervening sequences had GDH transport efficiency comparable with that of N53. Mutagenesis of the cleavage site blocked the intra-mitochondrial processing of hGDHs, but did not affect their mitochondrial import. Replacement of all three positively charged N-terminal residues (Arg3, Lys7 and Arg13) by Ala abolished import. We conclude that the synergistic interaction of helices α1 and α2 is crucial for the highly efficient import of hGDHs into mitochondria.