Search

The accumulation of somatic DNA mutations over time is a hallmark of aging in many dividing and nondividing cells but has not been studied in postmitotic human cardiomyocytes. Using single-cell whole-genome sequencing, we identified and characterized the landscape of somatic single-nucleotide variants (sSNVs) in 56 single cardiomyocytes from 12 individuals (aged from 0.4 to 82 years). Cardiomyocyte sSNVs accumulate with age at rates that are faster than in many dividing cell types and nondividing neurons. Cardiomyocyte sSNVs show distinctive mutational signatures that implicate failed nucleotide excision repair and base excision repair of oxidative DNA damage, and defective mismatch repair. Since age-accumulated sSNVs create many damaging mutations that disrupt gene functions, polyploidization in cardiomyocytes may provide a mechanism of genetic compensation to minimize the complete knockout of essential genes during aging. Age-related accumulation of cardiac mutations provides a paradigm to understand the influence of aging on cardiac dysfunction.

Dementia in Alzheimer’s disease progresses alongside neurodegeneration(1–4), but the specific events that cause neuronal dysfunction and death remain poorly understood. During normal ageing, neurons progressively accumulate somatic mutations(5) at rates similar to those of dividing cells(6,7) which suggests that genetic factors, environmental exposures or disease states might influence this accumulation(5). Here we analysed single-cell whole-genome sequencing data from 319 neurons from the prefrontal cortex and hippocampus of individuals with Alzheimer’s disease and neurotypical control individuals. We found that somatic DNA alterations increase in individuals with Alzheimer’s disease, with distinct molecular patterns. Normal neurons accumulate mutations primarily in an age-related pattern (signature A), which closely resembles ‘clock-like’ mutational signatures that have been previously described in healthy and cancerous cells(6–10). In neurons affected by Alzheimer’s disease, additional DNA alterations are driven by distinct processes (signature C) that highlight C>A and other specific nucleotide changes. These changes potentially implicate nucleotide oxidation(4,11), which we show is increased in Alzheimer’s-disease-affected neurons in situ. Expressed genes exhibit signature-specific damage, and mutations show a transcriptional strand bias, which suggests that transcription-coupled nucleotide excision repair has a role in the generation of mutations. The alterations in Alzheimer’s disease affect coding exons and are predicted to create dysfunctional genetic knockout cells and proteostatic stress. Our results suggest that known pathogenic mechanisms in Alzheimer’s disease may lead to genomic damage to neurons that can progressively impair function. The aberrant accumulation of DNA alterations in neurodegeneration provides insight into the cascade of molecular and cellular events that occurs in the development of Alzheimer’s disease.

ImportanceMesial temporal lobe epilepsy (MTLE) is the most common focal epilepsy subtype and is often refractory to anti-seizure medications. While most MTLE patients do not have pathogenic germline genetic variants, the contribution of post-zygotic (i.e., somatic) variants in the brain is unknown.ObjectiveTo test the association between pathogenic somatic variants in the hippocampus and MTLE.DesignThis case-control genetic association study analyzed the DNA derived from hippocampal tissue of neurosurgically-treated patients with MTLE and age- and sex-matched neurotypical controls. Participants were enrolled from 1988 through 2019 and clinical data was collected retrospectively. Whole-exome and gene-panel sequencing (depth>500X) were used to identify candidate pathogenic somatic variants. A subset of novel variants were functionally evaluated using cellular and molecular assays.SettingLevel 4 epilepsy centers, multi-center study.ParticipantsNon-lesional and lesional (mesial temporal sclerosis, focal cortical dysplasia, and low-grade epilepsy-associated tumors) drug-resistant MTLE patients who underwent anterior medial temporal lobectomy. All patients with available frozen tissue and appropriate consents were included. Control brain tissue was obtained from neurotypical donors at brain banks.ExposuresDrug-resistant MTLE.Main Outcomes and MeasuresPresence and abundance of pathogenic somatic variants in the hippocampus versus the unaffected temporal neocortex.ResultsSamples were obtained from 105 MTLE patients (52 male, 53 female; age: MED [IQR], 32 [26-44]) and 30 neurotypical controls (19 male, 11 female; age: MED [IQR], 37 [18-53]). Eleven pathogenic somatic variants, enriched in the hippocampus relative to the unaffected temporal neocortex (MED [IQR], 1.92 [1.5-2.7] vs 0.3 [0-0.9], pPTPN11, SOS1, KRAS, BRAF, andNF1, all predicted to constitutively activate Ras/Raf/MAPK signaling. Immunohistochemical studies of variant-positive hippocampal tissue demonstrated increased Erk1/2 phosphorylation, indicative of Ras/Raf/MAPK activation, predominantly in glial cells. Molecular assays showed abnormal liquid-liquid phase separation for thePTPN11variants as a possible dominant gain-of-function mechanism.Conclusions and RelevanceHippocampal somatic variants, particularly those activating Ras/Raf/MAPK signaling, may contribute to the pathogenesis of sporadic, drug-resistant MTLE. These findings may provide a novel genetic mechanism and highlight new therapeutic targets for this common indication for epilepsy surgery.

Somatic mutations are DNA variants that occur after the fertilization of zygotes and accumulate during the developmental and aging processes in the human lifespan. Somatic mutations have long been known to cause cancer, and more recently have been implicated in a variety of non-cancer diseases. The patterns of somatic mutations, or mutational signatures, also shed light on the underlying mechanisms of the mutational process. Advances in next-generation sequencing over the decades have enabled genome-wide profiling of DNA variants in a high-throughput manner; however, unlike germline mutations, somatic mutations are carried only by a subset of the cell population. Thus, sensitive bioinformatic methods are required to distinguish mutant alleles from sequencing and base calling errors in bulk tissue samples. An alternative way to study somatic mutations, especially those present in an extremely small number of cells or even in a single cell, is to sequence single-cell genomes after whole-genome amplification (WGA); however, it is critical and technically challenging to exclude numerous technical artifacts arising during error-prone and uneven genome amplification in current WGA methods. To address these challenges, multiple bioinformatic tools have been developed. In this review, we summarize the latest progress in methods for identification of somatic mutations and the challenges that remain to be addressed in the future.

Replication errors and various genotoxins cause DNA double-strand breaks (DSBs) where error-prone repair creates genomic mutations, most frequently focal deletions, and defective repair may lead to neurodegeneration. Despite its pathophysiological importance, the extent to which faulty DSB repair alters the genome, and the mechanisms by which mutations arise, have not been systematically examined reflecting ineffective methods. Here, we develop PhaseDel, a computational method to detect focal deletions and characterize underlying mechanisms in single-cell whole genome sequences (scWGS). We analyzed high-coverage scWGS of 107 single neurons from 18 neurotypical individuals of various ages, and found that somatic deletions increased with age and in highly expressed genes in human brain. Our analysis of 50 single neurons from DNA repair-deficient diseases with progressive neurodegeneration (Cockayne syndrome, Xeroderma pigmentosum, and Ataxia telangiectasia) reveals elevated somatic deletions compared to age-matched controls. Distinctive mechanistic signatures and transcriptional associations suggest roles for somatic deletions in neurodegeneration.

The accumulation of somatic DNA mutations is a hallmark of aging in many dividing cells and contributes to carcinogenesis. Here we survey the landscape of somatic single-nucleotide variants (sSNVs) in heart muscle cells (cardiomyocytes) which normally do not proliferate but often become polyploid with age. Using single-cell whole-genome sequencing we analyzed sSNVs from 48 single cardiomyocytes from 10 healthy individuals (ages 0.4 - 82 yrs.). Cardiomyocyte sSNVs increased strikingly with age, at rates faster than reported in many dividing cells, or in non-dividing neurons. Analysis of nucleotide substitution patterns revealed age-related “clock-like” mutational signatures resembling those previously described. However, cardiomyocytes showed distinct mutational signatures that are rare or absent in other cells, implicating failed nucleotide excision repair of oxidative damage and defective mismatch repair (MMR) during aging. A lineage tree of cardiomyocytes, constructed using clonal sSNVs, revealed that some tetraploid (10%) and most cardiomyocytes with higher ploidy (>60%) derive from distinct clonal origins, implicating cell fusion as a mechanism contributing to many polyploid cardiomyocytes. Since age-accumulated sSNVs create dozens of damaging exonic mutations, cell fusion to form multiploid cardiomyocytes may represent an evolutionary mechanism of cellular genetic compensation that minimizes complete knockout of essential genes during aging. The rates and patterns of accumulation of cardiac mutations provide a paradigm to understand the influence of genomic aging on age-related heart disease.

Elucidating the lineage relationships among different cell types is key to understanding human brain development. Here we developedParallelRNA andDNA analysis afterDeep-sequencing (PRDD-seq), which combines RNA analysis of neuronal cell types with analysis of nested spontaneous DNA somatic mutations as cell lineage markers, identified from joint analysis of single cell and bulk DNA sequencing by single-cell MosaicHunter (scMH). PRDD-seq enables the first-ever simultaneous reconstruction of neuronal cell type, cell lineage, and sequential neuronal formation (“birthdate”) in postmortem human cerebral cortex. Analysis of two human brains showed remarkable quantitative details that relate mutation mosaic frequency to clonal patterns, confirming an early divergence of precursors for excitatory and inhibitory neurons, and an “inside-out” layer formation of excitatory neurons as seen in other species. In addition our analysis allows the first estimate of excitatory neuron-restricted precursors (about 10) that generate the excitatory neurons within a cortical column. Inhibitory neurons showed complex, subtype-specific patterns of neurogenesis, including some patterns of development conserved relative to mouse, but also some aspects of primate cortical interneuron development not seen in mouse. PRDD-seq can be broadly applied to characterize cell identity and lineage from diverse archival samples with single-cell resolution and in potentially any developmental or disease condition.Significance StatementStem cells and progenitors undergo a series of cell divisions to generate the neurons of the brain, and understanding this sequence is critical to studying the mechanisms that control cell division and migration in developing brain. Mutations that occur as cells divide are known as the basis of cancer, but have more recently been shown to occur with normal cell divisions, creating a permanent, forensic map of the clonal patterns that define the brain. Here we develop new technology to analyze both DNA mutations and RNA gene expression patterns in single cells from human postmortem brain, allowing us to define clonal patterns among different types of human brain neurons, gaining the first direct insight into how they form.

Mosaic variants resulting from postzygotic mutations are prevalent in the human genome and play important roles in human diseases. However, except for cancer-related variants, there is no collection of postzygotic mosaic variants in noncancer disease-related and healthy individuals. Here, we present MosaicBase, a comprehensive database that includes 6698 mosaic variants related to 266 noncancer diseases and 27,991 mosaic variants identified in 422 healthy individuals. Genomic and phenotypic information of each variant was manually extracted and curated from 383 publications. MosaicBase supports the query of variants with Online Mendelian Inheritance in Man (OMIM) entries, genomic coordinates, gene symbols, or Entrez IDs. We also provide an integrated genome browser for users to easily access mosaic variants and their related annotations for any genomic region. By analyzing the variants collected in MosaicBase, we find that mosaic variants that directly contribute to disease phenotype show features distinct from those of variants in individuals with mild or no phenotypes, in terms of their genomic distribution, mutation signatures, and fraction of mutant cells. MosaicBase will not only assist clinicians in genetic counseling and diagnosis but also provide a useful resource to understand the genomic baseline of postzygotic mutations in the general human population. MosaicBase is publicly available at http://mosaicbase.com/ or http://49.4.21.8:8000.

Mosaic variants resulting from postzygotic mutations are prevalent in the human genome and play important roles in human diseases. However, except for cancer-related variant collections, there are no collections of mosaic variants in noncancer diseases and asymptomatic individuals. Here, we present MosaicBase (http://mosaicbase.cbi.pku.edu.cn/ or http://49.4.21.8:8000/), a comprehensive database that includes 6,698 mosaic variants related to 269 noncancer diseases and 27,991 mosaic variants identified in 422 asymptomatic individuals. The genomic and phenotypic information for each variant was manually extracted and curated from 383 publications. MosaicBase supports the query of variants with Online Mendelian Inheritance in Man (OMIM) entries, genomic coordinates, gene symbols, or Entrez IDs. We also provide an integrated genome browser for users to easily access mosaic variants and their related annotations within any genomic region. By analyzing the variants collected in MosaicBase, we found that mosaic variants that directly contribute to disease phenotype showed features distinct from those of variants in individuals with a mild or no phenotype in terms of their genomic distribution, mutation signatures, and fraction of mutant cells. MosaicBase will not only assist clinicians in genetic counseling and diagnosis but also provide a useful resource to understand the genomic baseline of postzygotic mutations in the general human population.

Congenital lipomatous overgrowth with vascular, epidermal, and skeletal (CLOVES) anomalies and Klippel-Trenaunay (KTS) syndromes are caused by somatic gain-of-function mutations in PIK3CA, encoding a catalytic subunit of phosphoinositide 3-kinase. Affected tissue is needed to find mutations, as mutant alleles are not detectable in blood. Because some patients with CLOVES develop Wilms tumor, we tested urine as a source of DNA for mutation detection. We extracted DNA from the urine of 17 and 24 individuals with CLOVES and KTS, respectively, and screened 5 common PIK3CA mutation hotspots using droplet digital polymerase chain reaction. Six of 17 CLOVES participants (35%) had mutant PIK3CA alleles in urine. Among 8 individuals in whom a mutation had been previously identified in affected tissue, 4 had the same mutant allele in the urine. One study participant with CLOVES had been treated for Wilms tumor. We detected the same PIK3CA mutation in her affected tissue, urine, and tumor, indicating Wilms tumors probably arise from PIK3CA mutant cells in patients with CLOVES. No urine sample from a participant with KTS had detectable PIK3CA mutations. We suggest that urine, which has the advantage of being collected non-invasively, is useful when searching for mutations in individuals with CLOVES syndrome.

Arteriovenous malformation (AVM) is a fast-flow, congenital vascular anomaly that may arise anywhere in the body. AVMs typically progress, causing destruction of surrounding tissue and, sometimes, cardiac overload. AVMs are difficult to control; they often re-expand after embolization or resection, and pharmacologic therapy is unavailable. We studied extracranial AVMs in order to identify their biological basis. We performed whole-exome sequencing (WES) and whole-genome sequencing (WGS) on AVM tissue from affected individuals. Endothelial cells were separated from non-endothelial cells by immune-affinity purification. We used droplet digital PCR (ddPCR) to confirm mutations found by WES and WGS, to determine whether mutant alleles were enriched in endothelial or non-endothelial cells, and to screen additional AVM specimens. In seven of ten specimens, WES and WGS detected and ddPCR confirmed somatic mutations in mitogen activated protein kinase kinase 1 (MAP2K1), the gene that encodes MAP-extracellular signal-regulated kinase 1 (MEK1). Mutant alleles were enriched in endothelial cells and were not present in blood or saliva. 9 of 15 additional AVM specimens contained mutant MAP2K1 alleles. Mutations were missense or small in-frame deletions that affect amino acid residues within or adjacent to the protein’s negative regulatory domain. Several of these mutations have been found in cancers and shown to increase MEK1 activity. In summary, somatic mutations in MAP2K1 are a common cause of extracranial AVM. The likely mechanism is endothelial cell dysfunction due to increased MEK1 activity. MEK1 inhibitors, which are approved to treat several forms of cancer, are potential therapeutic agents for individuals with extracranial AVM.

Mutations that occur in cells of the body, called somatic mutations, cause human diseases including cancer and some neurological disorders1. In a recent study published in Nature, Lee et al.2 (hereafter “the Lee study”) reported somatic copy number gains of the APP gene, a known risk locus of Alzheimer’s disease (AD), in the neurons of AD-patients and controls (69% vs 25% of neurons with at least one APP copy gain on average). The authors argue that the mechanism of these copy number gains was somatic integration of APP mRNA into the genome, creating what they called genomic cDNA (gencDNA). We reanalyzed the data from the Lee study, revealing evidence that APP gencDNA originates mainly from contamination by exogenous APP recombinant vectors, rather from true somatic retrotransposition of endogenous APP. Our reanalysis of two recent whole exome sequencing (WES) datasets—one by the authors of the Lee study3 and the other by Park et al.4—revealed that reads claimed to support APP gencDNA in AD samples resulted from contamination by PCR products and mRNA, respectively. Lastly, we present our own single-cell whole genome sequencing (scWGS) data that show no evidence for somatic APP retrotransposition in AD neurons or in neurons from normal individuals of various ages.

SUMMARY Metastasis is a complex and poorly understood process. In pancreatic cancer, loss of the transforming growth factor (TGF)-β/BMP effector SMAD4 is correlated with changes in altered histopathological transitions, metastatic disease, and poor prognosis. In this study, we use isogenic cancer cell lines to identify SMAD4 regulated genes that contribute to the development of metastatic colonization. We perform an in vivo screen identifying FOSL1 as both a SMAD4 target and sufficient to drive colonization to the lung. The targeting of these genes early in treatment may provide a therapeutic benefit., Graphical abstract, In brief Loss of SMAD4 is associated with poor outcome in pancreatic cancer. Using an in vivo, isogenic metastatic colonization assay, Dai et al. identified SMAD4-regulated genes that affect metastasis but not primary tumor growth. FOSL1 is a SMAD4-regulated gene that is necessary and sufficient to drive metastatic colonization.

Capillary malformation is a cutaneous vascular anomaly that is present at birth, darkens over time, and can cause overgrowth of tissues beneath the stain. The lesion is caused by a somatic activating mutation in GNAQ. In a previous study, we were unable to identify a GNAQ mutation in patients with a capillary malformation involving an overgrown lower extremity. We hypothesized that mutations in GNA11 or GNA14, genes closely related to GNAQ, also may cause capillary malformations. Human capillary malformation tissue obtained from 8 patients that had tested negative for GNAQ mutations were studied. Lesions involved an extremity (n = 7) or trunk (n = 1). Droplet digital PCR (ddPCR) was used to detect GNA11 or GNA14 mutant cells (p.Arg183) in the specimens. Single molecule molecular inversion probe sequencing (smMIP-seq) was performed to search for other mutations in GNA11. Mutations were validated by subcloning and sequencing amplimers. We found a somatic GNA11 missense mutation (c.547C > T; p.Arg183Cys) in 3 patients with a diffuse capillary malformation of an extremity. Mutant allelic frequencies ranged from 0.3 to 5.0%. GNA11 or GNA14 mutations were not found in 5 affected tissues or in unaffected tissues (white blood cell DNA). GNA11 mutations are associated with extremity capillary malformations causing overgrowth. Pharmacotherapy that affects GNA11 signaling may prevent the progression of capillary malformations.

Mitochondrial flashes (mitoflashes) represent stochastic and discrete mitochondrial events that each comprises a burst of superoxide production accompanied by transient depolarization and matrix alkalinization in a respiratory mitochondrion. While mitochondrial Ca(2+) is shown to be an important regulator of mitoflash activity, little is known about its specific mechanism of action. Here we sought to determine possible molecular players that mediate the Ca(2+) regulation of mitoflashes by screening mitochondrial proteins containing the Ca(2+)-binding motifs. In silico analysis and targeted siRNA screening identified four mitoflash activators (MICU1, EFHD1, SLC25A23, SLC25A25) and one mitoflash inhibitor (LETM1) in terms of their ability to modulate mitoflash response to hyperosmotic stress. In particular, overexpression or down-regulation of EFHD1 enhanced or depressed mitoflash activation, respectively, under various conditions of mitochondrial Ca(2+) elevations. Yet, it did not alter mitochondrial Ca(2+) handling, mitochondrial respiration, or ROS-induced mitoflash production. Further, disruption of the two EF-hand motifs of EFHD1 abolished its potentiating effect on the mitoflash responses. These results indicate that EFHD1 functions as a novel mitochondrial Ca(2+) sensor underlying Ca(2+)-dependent activation of mitoflashes.

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with strong genetic contributions. To provide a comprehensive resource for the genetic evidence of ASD, we have updated the Autism KnowledgeBase (AutismKB) to version 2.0. AutismKB 2.0 integrates multiscale genetic data on 1379 genes, 5420 copy number variations and structural variations, 11 669 single-nucleotide variations or small insertions/deletions (SNVs/indels) and 172 linkage regions. In particular, AutismKB 2.0 highlights 5669 de novo SNVs/indels due to their significant contribution to ASD genetics and includes 789 mosaic variants due to their recently discovered contributions to ASD pathogenesis. The genes and variants are annotated extensively with genetic evidence and clinical evidence. To help users fully understand the functional consequences of SNVs and small indels, we provided comprehensive predictions of pathogenicity with iFish, SIFT, Polyphen etc. To improve user experiences, the new version incorporates multiple query methods, including simple query, advanced query and batch query. It also functionally integrates two analytical tools to help users perform downstream analyses, including a gene ranking tool and an enrichment analysis tool, KOBAS. AutismKB 2.0 is freely available and can be a valuable resource for researchers.

Genomic mosaicism in parental gametes and peripheral tissues is an important consideration for genetic counseling. We studied a Chinese cohort affected by a severe epileptic disorder, Dravet syndrome (DS). There were 56 fathers who donated semen and 15 parents who donated multiple peripheral tissue samples. We used an ultra-sensitive quantification method, micro-droplet digital PCR (mDDPCR), to detect parental mosaicism of the proband’s pathogenic mutation in SCN1A, the causal gene of DS in 112 families. Ten of the 56 paternal sperm samples were found to exhibit mosaicism of the proband’s mutations, with mutant allelic fractions (MAFs) ranging from 0.03% to 39.04%. MAFs in the mosaic fathers’ sperm were significantly higher than those in their blood (p = 0.00098), even after conditional probability correction (p’ = 0.033). In three mosaic fathers, ultra-low fractions of mosaicism (MAF p = 3.006e-06). Our study provides new insights for genetic counseling.

Mice are commonly used to study the temporomandibular joint (TMJ) and to model human TMJ disease. However, evaluating TMJ pathology in mice using standard histologic methods is time consuming, labor intensive, and dependent upon investigators’ expertise at consistently orienting and sectioning across tiny specimens. We describe a method that uses confocal microscopy to rapidly and reliably assess indicators of mandibular condyle cartilage pathology in mice. We demonstrate the utility of this method for detecting abnormalities in chondrocyte distribution in mice lacking lubricin (Prg4), the major boundary lubricant of articular cartilage. We further show that the method can provide information about recombination sites and efficiency in mandibular cartilage for Cre-driver strains. Because specimen preparation and data acquisition with confocal microscopy are simple and fast, the method can serve as a primary screening tool for TMJ pathology, before proceeding to complicated, time consuming, secondary analyses.

RNA editing increases transcriptome diversity through post-transcriptional insertion, deletion, or modification of certain bases on RNA molecules (Licatalosi and Darnell 2010). RNA editing has been found to play important roles, such as altering specific codons of neurotransmitters or ion channels (Higuchi et al. 1993; Burns et al. 1997), creating new splicing sites (Rueter et al. 1999), modifying miRNA seed sequences or their targeting sites (Kawahara et al. 2007a,b), and protecting pre-mRNA by competing with the RNAi pathway (Wu et al. 2011; Warf et al. 2012). RNA editing is a widespread phenomenon in all metazoans (Danecek et al. 2012). The vast majority of known RNA editing comprises adenosine (A) to inosine (I) deamination, which is catalyzed by a family of adenosine deaminases that act on RNA (ADARs) (Keegan et al. 2001). ADARs are essential in mammals, but their individual effects in A-to-I editing remain unclear. ADAR1+/− mutant mice died before embryonic day 14 with defects in the hematopoietic system (Wang et al. 2000). ADAR2−/− mice became progressively prone to seizures after P12 and died young (Higuchi et al. 2000). In invertebrates, Drosophila melanogaster has a single ADAR2-like gene called Dmel\Adar (also known as dADAR), the deletion of which causes locomotor incoordination and temperature-sensitive paralysis (Palladino et al. 2000). Caenorhabditis elegans has two ADAR genes, adr-1 and adr-2, but their individual roles in A-to-I editing remain unclear. Worms lacking either or both genes are viable and exhibit only minor chemotaxis defects (Tonkin et al. 2002), which makes C. elegans an excellent model to study the effects of ADARs. Next-generation sequencing and bioinformatic technologies have recently been used to identify millions of RNA editing events at the genome scale in humans (Peng et al. 2012; Ramaswami et al. 2012, 2013; Bazak et al. 2014), mice (Danecek et al. 2012; Gu et al. 2012), and fruit flies (Graveley et al. 2011). However, the RNA editome in C. elegans remains largely unknown. Thus far, only 10 genes in C. elegans had been found to be edited in noncoding regions (Morse and Bass 1999; Morse et al. 2002). It is also unclear whether RNA editing is regulated across different developmental stages. In this study, we profiled the RNA editomes of different developmental stages of the Bristol N2 wild type and three ADAR mutants using strand-specific RNA sequencing (RNA-seq) and whole-genome DNA resequencing. We also developed a new computational pipeline using a bisulfite-seq-mapping-like approach that achieved higher sensitivity and specificity. This addressed the computational challenges in identifying editing sites, including mapping reads, anchoring junction sites, and differentiating sequencing errors from true editing events (Bass et al. 2012; Piskol et al. 2013).

Postzygotic single-nucleotide mutations (pSNMs) have been studied in cancer and a few other overgrowth human disorders at whole-genome scale and found to play critical roles. However, in clinically unremarkable individuals, pSNMs have never been identified at whole-genome scale largely due to technical difficulties and lack of matched control tissue samples, and thus the genome-wide characteristics of pSNMs remain unknown. We developed a new Bayesian-based mosaic genotyper and a series of effective error filters, using which we were able to identify 17 SNM sites from ∼80× whole-genome sequencing of peripheral blood DNAs from three clinically unremarkable adults. The pSNMs were thoroughly validated using pyrosequencing, Sanger sequencing of individual cloned fragments, and multiplex ligation-dependent probe amplification. The mutant allele fraction ranged from 5%-31%. We found that C→T and C→A were the predominant types of postzygotic mutations, similar to the somatic mutation profile in tumor tissues. Simulation data showed that the overall mutation rate was an order of magnitude lower than that in cancer. We detected varied allele fractions of the pSNMs among multiple samples obtained from the same individuals, including blood, saliva, hair follicle, buccal mucosa, urine, and semen samples, indicating that pSNMs could affect multiple sources of somatic cells as well as germ cells. Two of the adults have children who were diagnosed with Dravet syndrome. We identified two non-synonymous pSNMs in SCN1A, a causal gene for Dravet syndrome, from these two unrelated adults and found that the mutant alleles were transmitted to their children, highlighting the clinical importance of detecting pSNMs in genetic counseling.

BackgroundIntracranial arteriovenous malformation (AVM) is a common cause of primary intracerebral hemorrhage in young adults. Lesions typically are sporadic and contain somatic mutations in KRAS or BRAF. The purpose of this study was to identify somatic mutations in a cohort of participants with brain AVM and to determine if any genotype-phenotype associations exist.MethodsHuman brain AVM specimens (n = 16) were collected during a clinically-indicated procedure and subjected to multiplex targeted sequencing using molecular inversion probe (MIP-seq) for mutations in KRAS, BRAF, HRAS, NRAS, and MAP2K1. Endothelial cells (ECs) were separated from non-ECs by immune-affinity purification. Droplet digital PCR (ddPCR) was used to confirm mutations and to screen for mutations that may have been missed by MIP-seq. Patient and AVM characteristics were recorded.ResultsWe detected somatic mutations in 10 of 16 specimens (63%). Eight had KRAS mutations [G12D (n = 5), G12V (n = 3)] and two had BRAF mutations [V600E (n = 1), Q636X (n = 1)]. We found no difference in age, sex, presenting symptom, AVM location, or AVM size between patients with a confirmed mutation and those without. Nor did we observe differences in these features between patients with KRAS or BRAF mutations. However, two patients with BRAF mutations presented at an older age than other study participants.ConclusionsSomatic mutations in KRAS and, less commonly in BRAF, are found in many but not all intracranial AVM samples. Currently, there are no obvious genotype-phenotype correlations that can be used to predict whether a somatic mutation will be detected and, if so, which gene will be mutated.

The majority of children with Dravet syndrome (DS) are caused by de novo SCN1A mutations. To investigate the origin of the mutations, we developed and applied a new method that combined deep amplicon resequencing with a Bayesian model to detect and quantify allelic fractions with improved sensitivity. Of 174 SCN1A mutations in DS probands which were considered “de novo” by Sanger sequencing, we identified 15 cases (8.6%) of parental mosaicism. We identified another five cases of parental mosaicism that were also detectable by Sanger sequencing. Fraction of mutant alleles in the 20 cases of parental mosaicism ranged from 1.1% to 32.6%. Thirteen (65% of 20) mutations originated paternally and seven (35% of 20) maternally. Twelve (60% of 20) mosaic parents did not have any epileptic symptoms. Their mutant allelic fractions were significantly lower than those in mosaic parents with epileptic symptoms (P = 0.016). We identified mosaicism with varied allelic fractions in blood, saliva, urine, hair follicle, oral epithelium, and semen, demonstrating that postzygotic mutations could affect multiple somatic cells as well as germ cells. Our results suggest that more sensitive tools for detecting low‐level mosaicism in parents of families with seemingly “de novo” mutations will allow for better informed genetic counseling.

Genomic mosaicism arising from postzygotic mutations has long been associated with cancer and more recently with non-cancer diseases. It has also been detected in healthy individuals including healthy parents of children affected with genetic disorders, highlighting its critical role in the origin of genetic mutations. However, most existing software for the genome-wide identification of single-nucleotide mosaicisms (SNMs) requires a paired control tissue obtained from the same individual which is often unavailable for non-cancer individuals and sometimes missing in cancer studies. Here, we present MosaicHunter (http://mosaichunter.cbi.pku.edu.cn), a bioinformatics tool that can identify SNMs in whole-genome and whole-exome sequencing data of unpaired samples without matched controls using Bayesian genotypers. We evaluate the accuracy of MosaicHunter on both simulated and real data and demonstrate that it has improved performance compared with other somatic mutation callers. We further demonstrate that incorporating sequencing data of the parents can be an effective approach to significantly improve the accuracy of detecting SNMs in an individual when a matched control sample is unavailable. Finally, MosaicHunter also has a paired mode that can take advantage of matched control samples when available, making it a useful tool for detecting SNMs in both non-cancer and cancer studies.