Search

Mosaic mutations can be used to track cell ancestries and reconstruct high-resolution lineage trees during cancer progression and during development, starting from the first cell divisions of the zygote. However, this approach requires sampling and analyzing the genomes of multiple cells, which can be redundant in lineage representation, limiting the scalability of the approach. We describe a strategy for cost- and time-efficient lineage reconstruction using clonal induced pluripotent stem cell lines from human skin fibroblasts. The approach leverages shallow sequencing coverage to assess the clonality of the lines, clusters redundant lines and sums their coverage to accurately discover mutations in the corresponding lineages. Only a fraction of lines needs to be sequenced to high coverage. We demonstrate the effectiveness of this approach for reconstructing lineage trees during development and in hematologic malignancies. We discuss and propose an optimal experimental design for reconstructing lineage trees.

SummaryThere is no clear genetic etiology or convergent pathophysiology for autism spectrum disorders (ASD). Using cortical organoids and single-cell transcriptomics, we modeled alterations in the formation of the forebrain between sons with idiopathic ASD and their unaffected fathers in thirteen families. Alterations in the transcriptome suggest that ASD pathogenesis in macrocephalic and normocephalic probands involves an opposite disruption of the balance between the excitatory neurons of the dorsal cortical plate and other lineages such as the early-generated neurons from the putative preplate. The imbalance stemmed from a divergent expression of transcription factors driving cell fate during early cortical development. While we did not find probands’ genomic variants explaining the observed transcriptomic alterations, a significant overlap between altered transcripts and reported ASD risk genes affected by rare variants suggests a degree of gene convergence between rare forms of ASD and developmental transcriptome in idiopathic ASD.

Accurate discovery of somatic mutations in a cell is a challenge that partially lays in immaturity of dedicated analytical approaches. Approaches comparing cell’s genome to a control bulk sample miss common mutations, while approaches to find such mutations from bulk suffer from low sensitivity. We developed a tool, All2, which enables accurate filtering of mutations in a cell from exhaustive comparison of cells’ genomes to each other without data for bulk(s). Based on all pair-wise comparisons, every variant call (point mutation, indel, and structural variant) is classified as either a germline variant, mosaic mutation, or false positive. As All2 allows for considering dropped-out regions, it is applicable to whole genome and exome analysis of cloned and amplified cells. By applying the approach to a variety of available data, we showed that its application reduces false positives, enables sensitive discovery of high frequency mutations, and is indispensable for conducting high resolution cell lineage tracing. All2 is freely available at https://github.com/abyzovlab/All2.

Accurate discovery of somatic mutations in a cell is a challenge that partially lays in immaturity of dedicated analytical approaches. Approaches comparing a cell’s genome to a control bulk sample miss common mutations, while approaches to find such mutations from bulk suffer from low sensitivity. We developed a tool, All2, which enables accurate filtering of mutations in a cell without the need for data from bulk(s). It is based on pair-wise comparisons of all cells to each other where every call for base pair substitution and indel is classified as either a germline variant, mosaic mutation, or false positive. As All2 allows for considering dropped-out regions, it is applicable to whole genome and exome analysis of cloned and amplified cells. By applying the approach to a variety of available data, we showed that its application reduces false positives, enables sensitive discovery of high frequency mutations, and is indispensable for conducting high resolution cell lineage tracing.

Somatic mosaicism, manifesting as single nucleotide variants (SNVs), mobile element insertions, and structural changes in the DNA, is a common phenomenon in human brain cells, with potential functional consequences. Using a clonal approach, we previously detected 200–400 mosaic SNVs per cell in three human fetal brains (15–21 wk postconception). However, structural variation in the human fetal brain has not yet been investigated. Here, we discover and validate four mosaic structural variants (SVs) in the same brains and resolve their precise breakpoints. The SVs were of kilobase scale and complex, consisting of deletion(s) and rearranged genomic fragments, which sometimes originated from different chromosomes. Sequences at the breakpoints of these rearrangements had microhomologies, suggesting their origin from replication errors. One SV was found in two clones, and we timed its origin to ∼14 wk postconception. No large scale mosaic copy number variants (CNVs) were detectable in normal fetal human brains, suggesting that previously reported megabase-scale CNVs in neurons arise at later stages of development. By reanalysis of public single nuclei data from adult brain neurons, we detected an extrachromosomal circular DNA event. Our study reveals the existence of mosaic SVs in the developing human brain, likely arising from cell proliferation during mid-neurogenesis. Although relatively rare compared to SNVs and present in ∼10% of neurons, SVs in developing human brain affect a comparable number of bases in the genome (∼6200 vs. ∼4000 bp), implying that they may have similar functional consequences.

Post-zygotic mutations incurred during DNA replication, DNA repair, and other cellular processes lead to somatic mosaicism. Somatic mosaicism is an established cause of various diseases, including cancers. However, detecting mosaic variants in DNA from non-cancerous somatic tissues poses significant challenges, particularly if the variants only are present in a small fraction of cells. Here, the Brain Somatic Mosaicism Network conducted a coordinated, multi-institutional study to: (i) examine the ability of existing methods to detect simulated somatic single nucleotide variants (SNVs) in DNA mixing experiments; (ii) generate multiple replicates of whole genome sequencing data from the dorsolateral prefrontal cortex, other brain regions, dura mater, and dural fibroblasts of a single neurotypical individual; (iii) devise strategies to discover somatic SNVs; and (iv) apply various approaches to validate somatic SNVs. These efforts led to the identification of 43 bona fide somatic SNVs that ranged in variant allele fractions from ~0.005 to ~0.28. Guided by these results, we devised best practices for calling mosaic SNVs from 250X whole genome sequencing data in the accessible portion of the human genome that achieve 90% specificity and sensitivity. Finally, we demonstrated that analysis of multiple bulk DNA samples from a single individual allows the reconstruction of early developmental cell lineage trees. Thus, this study provides a unified set of best practices to detect somatic SNVs in non-cancerous tissues. The data and methods are freely available to the scientific community and should serve as a guide to assess the contributions of somatic SNVs to neuropsychiatric diseases.

Post-zygotic mosaic mutations can be used to track cell lineages in humans. By using cell cloning and induced pluripotent cell lines, we analyzed early cell lineages in two living individuals (a patient and a control), and a postmortem human specimen. Of ten reconstructed post-zygotic divisions, none resulted in balanced contributions of daughter lineages to tissues. In both living individuals one of two lineages from the first cleavage was dominant across tissues, with 90% frequency in blood. We propose that the efficiency of DNA repair contributes to lineage imbalance. Allocation of lineages in postmortem brain correlated with anterior-posterior axis, associating lineage history with cell fate choices in embryos. Recurrence of germline variants as mosaic suggested that certain loci may be particularly susceptible to mutagenesis. We establish a minimally invasive framework for defining cell lineages in any living individual, which paves the way for studying their relevance in health and disease.

SummaryAutism spectrum disorder (ASD) is a disorder of brain development. Most cases lack a clear etiology or genetic basis, and the difficulty of re-enacting human brain development has precluded understanding of ASD pathophysiology. Here we use three-dimensional neural cultures (organoids) derived from induced pluripotent stem cells (iPSCs) to investigate neurodevelopmental alterations in individuals with severe idiopathic ASD. While no known underlying genomic mutation could be identified, transcriptome and gene network analyses revealed upregulation of genes involved in cell proliferation, neuronal differentiation, and synaptic assembly. ASD-derived organoids exhibit an accelerated cell cycle and overproduction of GABAergic inhibitory neurons. Using RNA interference, we show that overexpression of the transcription factor FOXG1 is responsible for the overproduction of GABAergic neurons. Altered expression of gene network modules and FOXG1 are positively correlated with symptom severity. Our data suggest that a shift toward GABAergic neuron fate caused by FOXG1 is a developmental precursor of ASD.

Brain mutations, young and old Most neurons that make up the human brain are postmitotic, living and functioning for a very long time without renewal (see the Perspective by Lee). Bae et al. examined the genomes of single neurons from the prenatal developing human brain. Both the type of mutation and the rates of accumulation changed between gastrulation and neurogenesis. These early mutations could be generating useful neuronal diversity or could predispose individuals to later dysfunction. Lodato et al. also found that neurons take on somatic mutations as they age by sequencing single neurons from subjects aged 4 months to 82 years. Somatic mutations accumulated with increasing age and accumulated faster in individuals affected by inborn errors in DNA repair. Postmitotic mutations might only affect one neuron, but the accumulated divergence of genomes across the brain could affect function. Science , this issue p. 550 , p. 555 ; see also p. 521

Few studies have been conducted to understand post-zygotic accumulation of mutations in cells of the healthy human body. We reprogrammed 32 skin fibroblast cells from families of donors into human induced pluripotent stem cell (hiPSC) lines. The clonal nature of hiPSC lines allows a high-resolution analysis of the genomes of the founder fibroblast cells without being confounded by the artifacts of single-cell whole-genome amplification. We estimate that on average a fibroblast cell in children has 1035 mostly benign mosaic SNVs. On average, 235 SNVs could be directly confirmed in the original fibroblast population by ultradeep sequencing, down to an allele frequency (AF) of 0.1%. More sensitive droplet digital PCR experiments confirmed more SNVs as mosaic with AF as low as 0.01%, suggesting that 1035 mosaic SNVs per fibroblast cell is the true average. Similar analyses in adults revealed no significant increase in the number of SNVs per cell, suggesting that a major fraction of mosaic SNVs in fibroblasts arises during development. Mosaic SNVs were distributed uniformly across the genome and were enriched in a mutational signature previously observed in cancers and in de novo variants and which, we hypothesize, is a hallmark of normal cell proliferation. Finally, AF distribution of mosaic SNVs had distinct narrow peaks, which could be a characteristic of clonal cell selection, clonal expansion, or both. These findings reveal a large degree of somatic mosaicism in healthy human tissues, link de novo and cancer mutations to somatic mosaicism, and couple somatic mosaicism with cell proliferation.

A whole-genome and transcriptome analysis of 20 human induced pluripotent stem-cell lines shows that reprogramming does not necessarily add de novo copy number variants to what is already present in the somatic cells from which they originated. The ability to derive induced pluripotent stem cells (iPSCs) from somatic cells raises exciting possibilities for the study of human development and regenerative medicine. These applications require that the clonal cells maintain the genetic background of the individual from whom they are derived, so reports of chromosomal copy number variations (CNVs) in reprogrammed cells carry serious implications for their translational utility. Flora Vaccarino and colleagues now report a whole-genome and transcriptome analysis of 20 human iPSC lines from seven individuals. They found that reprogramming does not necessarily add de novo CNVs to those already present in the somatic genome. Interestingly, they also found a mosaic CNV pattern within individuals, confirming previous findings from cultured human fibroblasts. This work shows that iPSCs can be used as a discovery tool for the investigation of genomic mosaicism due to low-frequency CNVs in human tissues. Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) has been suspected of causing de novo copy number variation1,2,3,4. To explore this issue, here we perform a whole-genome and transcriptome analysis of 20 human iPSC lines derived from the primary skin fibroblasts of seven individuals using next-generation sequencing. We find that, on average, an iPSC line manifests two copy number variants (CNVs) not apparent in the fibroblasts from which the iPSC was derived. Using PCR and digital droplet PCR, we show that at least 50% of those CNVs are present as low-frequency somatic genomic variants in parental fibroblasts (that is, the fibroblasts from which each corresponding human iPSC line is derived), and are manifested in iPSC lines owing to their clonal origin. Hence, reprogramming does not necessarily lead to de novo CNVs in iPSCs, because most of the line-manifested CNVs reflect somatic mosaicism in the human skin. Moreover, our findings demonstrate that clonal expansion, and iPSC lines in particular, can be used as a discovery tool to reliably detect low-frequency CNVs in the tissue of origin. Overall, we estimate that approximately 30% of the fibroblast cells have somatic CNVs in their genomes, suggesting widespread somatic mosaicism in the human body. Our study paves the way to understanding the fundamental question of the extent to which cells of the human body normally acquire structural alterations in their DNA post-zygotically.

International audience; Mutations of the inositol 5' phosphatase oculocerebrorenal syndrome of Lowe (OCRL) give rise to the congenital X-linked disorders oculocerebrorenal syndrome of Lowe and Dent disease, two conditions giving rise to abnormal kidney proximal tubule reabsorption, and additional nervous system and ocular defects in the case of Lowe syndrome. Here, we identify two closely related endocytic proteins, Ses1 and Ses2, which interact with the ASH-RhoGAP-like (ASPM-SPD-2-Hydin homology and Rho-GTPase Activating Domain-like) domain of OCRL. The interaction is mediated by a short amino acid motif similar to that used by the rab-5 effector APPL1 (Adaptor Protein containing pleckstrin homology [PH] domain, PTB domain and Leucine zipper motif 1) APPL1 for OCRL binding. Ses binding is mutually exclusive with APPL1 binding, and is disrupted by the same missense mutations in the ASH-RhoGAP-like domain that also disrupt APPL1 binding. Like APPL1, Ses1 and -2 are localized on endosomes but reside on different endosomal subpopulations. These findings define a consensus motif (which we have called a phenylalanine and histidine [F&H] motif) for OCRL binding and are consistent with a scenario in which Lowe syndrome and Dent disease result from perturbations at multiple sites within the endocytic pathway.

OCRL, whose mutations are responsible for Lowe syndrome and Dent disease, and INPP5B are two similar proteins comprising a central inositol 5-phosphatase domain followed by an ASH and a RhoGAP-like domain. Their divergent NH2-terminal portions remain uncharacterized. We show that the NH2-terminal region of OCRL, but not of INPP5B, binds clathrin heavy chain. OCRL, which in contrast to INPP5B visits late stage endocytic clathrin-coated pits, was earlier shown to contain another binding site for clathrin in its COOH-terminal region. NMR structure determination further reveals that despite their primary sequence dissimilarity, the NH2-terminal portions of both OCRL and INPP5B contain a PH domain. The novel clathrin-binding site in OCRL maps to an unusual clathrin-box motif located in a loop of the PH domain, whose mutations reduce recruitment efficiency of OCRL to coated pits. These findings suggest an evolutionary pressure for a specialized function of OCRL in bridging phosphoinositide metabolism to clathrin-dependent membrane trafficking.

Mutations in the inositol 5-phosphatase OCRL are responsible for Lowe syndrome, an X-linked disorder characterized by bilateral cataracts, mental retardation, neonatal hypotonia, and renal Fanconi syndrome, and for Dent disease, another X-linked condition characterized by kidney reabsorption defects. We have previously described an interaction of OCRL with the endocytic adaptor APPL1 that links OCRL to protein networks involved in the disease phenotype. Here, we provide new evidence showing that among the interactions which target OCRL to membranes of the endocytic pathway, binding to APPL1 is the only one abolished by all known disease-causing missense mutations in the ASH-RhoGAP domains of the protein. Furthermore, we demonstrate that APPL1 and rab5 independently contribute to recruit OCRL to enlarged endosomes induced by the expression of constitutively active Rab5. Thus, binding to APPL1 helps localize OCRL at specific cellular sites, and disruption of this interaction may play a role in disease.

Idiopathic generalized epilepsy (IGE) is a common, complex disease with an almost exclusively genetic etiology but with variable phenotypes. Clinically, IGE can be divided into different syndromes. Varying lines of evidence point to the involvement of several interacting genes in the etiology of IGE. We performed a genome scan in 91 families ascertained through a proband with adolescent-onset IGE. The IGEs included juvenile myoclonic epilepsy (JME), juvenile absence epilepsy (JAE), and epilepsy with generalized tonic clonic seizures (EGTCS). Our linkage results support an oligogenic model for IGE, with strong evidence for a locus common to most IGEs on chromosome 18 (lod score 4.4/5.2 multipoint/ two-point) and other loci that may influence specific seizure phenotypes for different IGEs: a previously identified locus on chromosome 6 for JME (lod score 2.5/4.2), a locus on chromosome 8 influencing non-JME forms of IGE (lod score 3.8/2.5), and, more tentatively, two newly discovered loci for absence seizures on chromosome 5 (lod scores 3.8/2.8 and 3.4/1.9). Our data also suggest that the genetic classification of different forms of IGE is likely to cut across the clinical classification of these subforms of IGE. We hypothesize that interactions of different combinations of these loci produce the related heterogeneous phenotypes seen in IGE families. Ann Neurol 2001;49:328 ‐335

Evidence for genetic influences in epilepsy is strong, but reports identifying specific chromosomal origins of those influences conflict. One early study reported that human leukocyte antigen (HLA) markers were genetically linked to juvenile myoclonic epilepsy (JME); this was confirmed in a later study. Other reports did not find linkage to HLA markers. One found evidence of linkage to markers on chromosome 15, another to markers on chromosome 6, centromeric to HLA. We identified families through a patient with JME and genotyped markers throughout chromosome 6. Linkage analysis assuming equal male-female recombination probabilities showed evidence for linkage (LOD score 2.5), but at a high recombination fraction (theta), suggesting heterogeneity. When linkage analysis was redone to allow independent male-female thetas, the LOD score was significantly higher (4.2) at a male-female theta of.5,.01. Although the overall pattern of LOD scores with respect to male-female theta could not be explained solely by heterogeneity, the presence of heterogeneity and predominantly maternal inheritance of JME might explain it. By analyzing loci between HLA-DP and HLA-DR and stratifying the families on the basis of evidence for or against linkage, we were able to show evidence of heterogeneity within JME and to propose a marker associated with the linked form. These data also suggest that JME may be predominantly maternally inherited and that the HLA-linked form is more likely to occur in families of European origin.