Resolving the 22q11.2 deletion using CTLR-Seq reveals chromosomal rearrangement mechanisms and individual variance in breakpoints.

We developed a generally applicable method, CRISPR/Cas9-targeted long-read sequencing (CTLR-Seq), to resolve, haplotype-specifically, the large and complex regions in the human genome that had been previously impenetrable to sequencing analysis, such as large segmental duplications (SegDups) and their associated genome rearrangements. CTLR-Seq combines in vitro Cas9-mediated cutting of the genome and pulse-field gel electrophoresis to isolate intact large (i.e., up to 2,000 kb) genomic regions that encompass previously unresolvable genomic sequences. These targets are then sequenced (amplification-free) at high on-target coverage using long-read sequencing, allowing for their complete sequence assembly. We applied CTLR-Seq to the SegDup-mediated rearrangements that constitute the boundaries of, and give rise to, the 22q11.2 Deletion Syndrome (22q11DS), the most common human microdeletion disorder. We then performed de novo assembly to resolve, at base-pair resolution, the full sequence rearrangements and exact chromosomal breakpoints of 22q11.2DS (including all common subtypes). Across multiple patients, we found a high degree of variability for both the rearranged SegDup sequences and the exact chromosomal breakpoint locations, which coincide with various transposons within the 22q11.2 SegDups, suggesting that 22q11DS can be driven by transposon-mediated genome recombination. Guided by CTLR-Seq results from two 22q11DS patients, we performed three-dimensional chromosomal folding analysis for the 22q11.2 SegDups from patient-derived neurons and astrocytes and found chromosome interactions anchored within the SegDups to be both cell type-specific and patient-specific. Lastly, we demonstrated that CTLR-Seq enables cell-type specific analysis of DNA methylation patterns within the deletion haplotype of 22q11DS.

Structural variation in the sequencing era.

Identifying structural variation (SV) is essential for genome interpretation but has been historically difficult due to limitations inherent to available genome technologies. Detection methods that use ensemble algorithms and emerging sequencing technologies have enabled the discovery of thousands of SVs, uncovering information about their ubiquity, relationship to disease and possible effects on biological mechanisms. Given the variability in SV type and size, along with unique detection biases of emerging genomic platforms, multiplatform discovery is necessary to resolve the full spectrum of variation. Here, we review modern approaches for investigating SVs and proffer that, moving forwards, studies integrating biological information with detection will be necessary to comprehensively understand the impact of SV in the human genome.

Mutations in human DNA methyltransferase DNMT1 induce specific genome-wide epigenomic and transcriptomic changes in neurodevelopment.

DNA methyltransferase type 1 (DNMT1) is a major enzyme involved in maintaining the methylation pattern after DNA replication. Mutations in DNMT1 have been associated with autosomal dominant cerebellar ataxia, deafness and narcolepsy (ADCA-DN). We used fibroblasts, induced pluripotent stem cells (iPSCs) and induced neurons (iNs) generated from patients with ADCA-DN and controls, to explore the epigenomic and transcriptomic effects of mutations in DNMT1. We show cell type-specific changes in gene expression and DNA methylation patterns. DNA methylation and gene expression changes were negatively correlated in iPSCs and iNs. In addition, we identified a group of genes associated with clinical phenotypes of ADCA-DN, including PDGFB and PRDM8 for cerebellar ataxia, psychosis and dementia and NR2F1 for deafness and optic atrophy. Furthermore, ZFP57, which is required to maintain gene imprinting through DNA methylation during early development, was hypomethylated in promoters and exhibited upregulated expression in patients with ADCA-DN in both iPSC and iNs. Our results provide insight into the functions of DNMT1 and the molecular changes associated with ADCA-DN, with potential implications for genes associated with related phenotypes.

Neural signatures of sleep in zebrafish.

Slow-wave sleep and rapid eye movement (or paradoxical) sleep have been found in mammals, birds and lizards, but it is unclear whether these neuronal signatures are found in non-amniotic vertebrates. Here we develop non-invasive fluorescence-based polysomnography for zebrafish, and show-using unbiased, brain-wide activity recording coupled with assessment of eye movement, muscle dynamics and heart rate-that there are at least two major sleep signatures in zebrafish. These signatures, which we term slow bursting sleep and propagating wave sleep, share commonalities with those of slow-wave sleep and paradoxical or rapid eye movement sleep, respectively. Further, we find that melanin-concentrating hormone signalling (which is involved in mammalian sleep) also regulates propagating wave sleep signatures and the overall amount of sleep in zebrafish, probably via activation of ependymal cells. These observations suggest that common neural signatures of sleep may have emerged in the vertebrate brain over 450 million years ago.

Synthesis of Cs3Cu2I5 Nanocrystals in a Continuous Flow System.

Achieving the goal of generating all of the world's energy via renewable sources and significantly reducing the energy usage will require the development of novel, abundant, nontoxic energy conversion materials. Here, a cost-efficient and scalable continuous flow synthesis of Cs3Cu2I5 nanocrystals is developed as a basis for the rapid advancement of novel nanomaterials. Ideal precursor solutions are obtained through a novel batch synthesis, whose product served as a benchmark for the subsequent flow synthesis. Realizing this setup enabled a reproducible fabrication of Cs3Cu2I5 nanocrystals. The effect of volumetric flow rate and temperature on the final product's morphology and optical properties are determined, obtaining 21% quantum yield with the optimal configuration. Consequently, the size and morphology of the nanocrystals can be tuned with far more precision and in a much broader range than previously achievable. The flow setup is readily applicable to other relevant nanomaterials. It should enable a rapid determination of a material's potential and subsequently optimize its desired properties for renewable energy generation or efficient optoelectronics.

Discovering melting temperature prediction models of inorganic solids by combining supervised and unsupervised learning.

The melting temperature is important for materials design because of its relationship with thermal stability, synthesis, and processing conditions. Current empirical and computational melting point estimation techniques are limited in scope, computational feasibility, or interpretability. We report the development of a machine learning methodology for predicting melting temperatures of binary ionic solid materials. We evaluated different machine-learning models trained on a dataset of the melting points of 476 non-metallic crystalline binary compounds using materials embeddings constructed from elemental properties and density-functional theory calculations as model inputs. A direct supervised-learning approach yields a mean absolute error of around 180 K but suffers from low interpretability. We find that the fidelity of predictions can further be improved by introducing an additional unsupervised-learning step that first classifies the materials before the melting-point regression. Not only does this two-step model exhibit improved accuracy, but the approach also provides a level of interpretability with insights into feature importance and different types of melting that depend on the specific atomic bonding inside a material. Motivated by this finding, we used a symbolic learning approach to find interpretable physical models for the melting temperature, which recovered the best-performing features from both prior models and provided additional interpretability.

Cross-platform validation of neurotransmitter release impairments in schizophrenia patient-derived NRXN1-mutant neurons.

Heterozygous NRXN1 deletions constitute the most prevalent currently known single-gene mutation associated with schizophrenia, and additionally predispose to multiple other neurodevelopmental disorders. Engineered heterozygous NRXN1 deletions impaired neurotransmitter release in human neurons, suggesting a synaptic pathophysiological mechanism. Utilizing this observation for drug discovery, however, requires confidence in its robustness and validity. Here, we describe a multicenter effort to test the generality of this pivotal observation, using independent analyses at two laboratories of patient-derived and newly engineered human neurons with heterozygous NRXN1 deletions. Using neurons transdifferentiated from induced pluripotent stem cells that were derived from schizophrenia patients carrying heterozygous NRXN1 deletions, we observed the same synaptic impairment as in engineered NRXN1-deficient neurons. This impairment manifested as a large decrease in spontaneous synaptic events, in evoked synaptic responses, and in synaptic paired-pulse depression. Nrxn1-deficient mouse neurons generated from embryonic stem cells by the same method as human neurons did not exhibit impaired neurotransmitter release, suggesting a human-specific phenotype. Human NRXN1 deletions produced a reproducible increase in the levels of CASK, an intracellular NRXN1-binding protein, and were associated with characteristic gene-expression changes. Thus, heterozygous NRXN1 deletions robustly impair synaptic function in human neurons regardless of genetic background, enabling future drug discovery efforts.

Complex mosaic structural variations in human fetal brains.

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.

Evaluation of C4 Gene Copy Number in Pediatric Acute Neuropsychiatric Syndrome.

Pediatric acute-onset neuropsychiatric syndrome (PANS) is an abrupt-onset neuropsychiatric disorder. PANS patients have an increased prevalence of comorbid autoimmune illness, most commonly arthritis. In addition, an estimated one-third of PANS patients present with low serum C4 protein, suggesting decreased production or increased consumption of C4 protein. To test the possibility that copy number (CN) variation contributes to risk of PANS illness, we compared mean total C4A and total C4B CN in ethnically matched subjects from PANS DNA samples and controls (192 cases and 182 controls). Longitudinal data from the Stanford PANS cohort (n = 121) were used to assess whether the time to juvenile idiopathic arthritis (JIA) or autoimmune disease (AI) onset was a function of total C4A or C4B CN. Lastly, we performed several hypothesis-generating analyses to explore the correlation between individual C4 gene variants, sex, specific genotypes, and age of PANS onset. Although the mean total C4A or C4B CN did not differ in PANS compared to controls, PANS patients with low C4B CN were at increased risk for subsequent JIA diagnosis (hazard ratio = 2.7, p value = 0.004). We also observed a possible increase in risk for AI in PANS patients and a possible correlation between lower C4B and PANS age of onset. An association between rheumatoid arthritis and low C4B CN has been reported previously. However, patients with PANS develop different types of JIA: enthesitis-related arthritis, spondyloarthritis, and psoriatic arthritis. This suggests that C4B plays a role that spans these arthritis types.

Lead-Free, Luminescent Perovskite Nanocrystals Obtained through Ambient Condition Synthesis.

Heterovalently substituting toxic lead is an increasingly popular design strategy to obtain environmentally sustainable variants of the exciting material class of halide perovskites. Perovskite nanocrystals (NCs) obtained through solution-based methods exhibit exceedingly high optical quality. Unfortunately, most of these synthesis routes still require reaction under inert gas and at very high temperatures. Herein a novel synthesis routine for lead-free double perovskite (LFDP) NCs is presented. An approach based upon the hot injection and ligand-assisted reprecipitation (LARP) methods to achieve a low-temperature and ambient atmosphere-based synthesis for manganese-doped Cs2 NaBiCl6 NCs is presented. Mn incorporation is critical for the otherwise non-emissive material, with a 9:1 Bi:Mn precursor ratio maximizing the bright orange photoluminescence (PL) and quantum yield (QY). Higher synthesis temperatures slightly increase the material's performance, yet NCs synthesized at room temperature are still emissive, highlighting the versatility of the synthetic approach. While the material's indirect bandgap limits its appeal for optoelectronics, this feature could benefit photocatalysis due to longer carrier lifetimes. Moreover, the developed synthesis is facile and can rapidly be adapted to other more viable material compositions and up-scaled to realize applications directly.

Maternal attachment insecurity, maltreatment history, and depressive symptoms are associated with broad DNA methylation signatures in infants.

The early environment, including maternal characteristics, provides many cues to young organisms that shape their long-term physical and mental health. Identifying the earliest molecular events that precede observable developmental outcomes could help identify children in need of support prior to the onset of physical and mental health difficulties. In this study, we examined whether mothers' attachment insecurity, maltreatment history, and depressive symptoms were associated with alterations in DNA methylation patterns in their infants, and whether these correlates in the infant epigenome were associated with socioemotional and behavioral functioning in toddlerhood. We recruited 156 women oversampled for histories of depression, who completed psychiatric interviews and depression screening during pregnancy, then provided follow-up behavioral data on their children at 18 months. Buccal cell DNA was obtained from 32 of their infants for a large-scale analysis of methylation patterns across 5 × 106 individual CpG dinucleotides, using clustering-based significance criteria to control for multiple comparisons. We found that tens of thousands of individual infant CpGs were alternatively methylated in association with maternal attachment insecurity, maltreatment in childhood, and antenatal and postpartum depressive symptoms, including genes implicated in developmental patterning, cell-cell communication, hormonal regulation, immune function/inflammatory response, and neurotransmission. Density of DNA methylation at selected genes from the result set was also significantly associated with toddler socioemotional and behavioral problems. This is the first report to identify novel regions of the human infant genome at which DNA methylation patterns are associated longitudinally both with maternal characteristics and with offspring socioemotional and behavioral problems in toddlerhood.

Integrated functional genomic analyses of Klinefelter and Turner syndromes reveal global network effects of altered X chromosome dosage.

In both Turner syndrome (TS) and Klinefelter syndrome (KS) copy number aberrations of the X chromosome lead to various developmental symptoms. We report a comparative analysis of TS vs. KS regarding differences at the genomic network level measured in primary samples by analyzing gene expression, DNA methylation, and chromatin conformation. X-chromosome inactivation (XCI) silences transcription from one X chromosome in female mammals, on which most genes are inactive, and some genes escape from XCI. In TS, almost all differentially expressed escape genes are down-regulated but most differentially expressed inactive genes are up-regulated. In KS, differentially expressed escape genes are up-regulated while the majority of inactive genes appear unchanged. Interestingly, 94 differentially expressed genes (DEGs) overlapped between TS and female and KS and male comparisons; and these almost uniformly display expression changes into opposite directions. DEGs on the X chromosome and the autosomes are coexpressed in both syndromes, indicating that there are molecular ripple effects of the changes in X chromosome dosage. Six potential candidate genes (RPS4X, SEPT6, NKRF, CX0rf57, NAA10, and FLNA) for KS are identified on Xq, as well as candidate central genes on Xp for TS. Only promoters of inactive genes are differentially methylated in both syndromes while escape gene promoters remain unchanged. The intrachromosomal contact map of the X chromosome in TS exhibits the structure of an active X chromosome. The discovery of shared DEGs indicates the existence of common molecular mechanisms for gene regulation in TS and KS that transmit the gene dosage changes to the transcriptome.

Energy Transfer in Stability-Optimized Perovskite Nanocrystals.

Outstanding optoelectronic properties and a facile synthesis render halide perovskite nanocrystals (NCs) a promising material for nanostructure-based devices. However, the commercialization is hindered mainly by the lack of NC stability under ambient conditions and inefficient charge carrier injection. Here, we investigate solutions to both problems, employing methylammonium lead bromide (MAPbBr3) NCs encapsulated in diblock copolymer core-shell micelles of tunable size. We confirm that the shell does not prohibit energy transfer, as FRET efficiencies between these NCs and 2D CsPbBr3 nanoplatelets (NPLs) reach 73.6%. This value strongly correlates to the micelle size, with thicker shells displaying significantly reduced FRET efficiencies. Those high efficiencies come with a price, as the thinnest shells protect the encapsulated NCs less from environmentally induced degradation. Finding the sweet spot between efficiency and protection could lead to the realization of tailored energy funnels with enhanced carrier densities for high-power perovskite NC-based optoelectronic applications.

Thickness-Dependent Dark-Bright Exciton Splitting and Phonon Bottleneck in CsPbBr3-Based Nanoplatelets Revealed via Magneto-Optical Spectroscopy.

The optimized exploitation of perovskite nanocrystals and nanoplatelets as highly efficient light sources requires a detailed understanding of the energy spacing within the exciton manifold. Dark exciton states are particularly relevant because they represent a channel that reduces radiative efficiency. Here, we apply large in-plane magnetic fields to brighten optically inactive states of CsPbBr3-based nanoplatelets for the first time. This approach allows us to access the dark states and directly determine the dark-bright splitting, which reaches 22 meV for the thinnest nanoplatelets. The splitting is significantly less for thicker nanoplatelets due to reduced exciton confinement. Additionally, the form of the magneto-PL spectrum suggests that dark and bright state populations are nonthermalized, which is indicative of a phonon bottleneck in the exciton relaxation process.

Comprehensive, integrated, and phased whole-genome analysis of the primary ENCODE cell line K562.

K562 is widely used in biomedical research. It is one of three tier-one cell lines of ENCODE and also most commonly used for large-scale CRISPR/Cas9 screens. Although its functional genomic and epigenomic characteristics have been extensively studied, its genome sequence and genomic structural features have never been comprehensively analyzed. Such information is essential for the correct interpretation and understanding of the vast troves of existing functional genomics and epigenomics data for K562. We performed and integrated deep-coverage whole-genome (short-insert), mate-pair, and linked-read sequencing as well as karyotyping and array CGH analysis to identify a wide spectrum of genome characteristics in K562: copy numbers (CN) of aneuploid chromosome segments at high-resolution, SNVs and indels (both corrected for CN in aneuploid regions), loss of heterozygosity, megabase-scale phased haplotypes often spanning entire chromosome arms, structural variants (SVs), including small and large-scale complex SVs and nonreference retrotransposon insertions. Many SVs were phased, assembled, and experimentally validated. We identified multiple allele-specific deletions and duplications within the tumor suppressor gene FHIT Taking aneuploidy into account, we reanalyzed K562 RNA-seq and whole-genome bisulfite sequencing data for allele-specific expression and allele-specific DNA methylation. We also show examples of how deeper insights into regulatory complexity are gained by integrating genomic variant information and structural context with functional genomics and epigenomics data. Furthermore, using K562 haplotype information, we produced an allele-specific CRISPR targeting map. This comprehensive whole-genome analysis serves as a resource for future studies that utilize K562 as well as a framework for the analysis of other cancer genomes.

PRG2 and AQPEP are misexpressed in fetal membranes in placenta previa and percreta†.

The obstetrical conditions placenta accreta spectrum (PAS) and placenta previa are a significant source of pregnancy-associated morbidity and mortality, yet the specific molecular and cellular underpinnings of these conditions are not known. In this study, we identified misregulated gene expression patterns in tissues from placenta previa and percreta (the most extreme form of PAS) compared with control cases. By comparing this gene set with existing placental single-cell and bulk RNA-Seq datasets, we show that the upregulated genes predominantly mark extravillous trophoblasts. We performed immunofluorescence on several candidate molecules and found that PRG2 and AQPEP protein levels are upregulated in both the fetal membranes and the placental disk in both conditions. While this increased AQPEP expression remains restricted to trophoblasts, PRG2 is mislocalized and is found throughout the fetal membranes. Using a larger patient cohort with a diverse set of gestationally aged-matched controls, we validated PRG2 as a marker for both previa and PAS and AQPEP as a marker for only previa in the fetal membranes. Our findings suggest that the extraembryonic tissues surrounding the conceptus, including both the fetal membranes and the placental disk, harbor a signature of previa and PAS that is characteristic of EVTs and that may reflect increased trophoblast invasiveness.

Metal Halide Perovskite Nanocrystals: Synthesis, Post-Synthesis Modifications, and Their Optical Properties.

Metal halide perovskites represent a flourishing area of research, which is driven by both their potential application in photovoltaics and optoelectronics and by the fundamental science behind their unique optoelectronic properties. The emergence of new colloidal methods for the synthesis of halide perovskite nanocrystals, as well as the interesting characteristics of this new type of material, has attracted the attention of many researchers. This review aims to provide an up-to-date survey of this fast-moving field and will mainly focus on the different colloidal synthesis approaches that have been developed. We will examine the chemistry and the capability of different colloidal synthetic routes with regard to controlling the shape, size, and optical properties of the resulting nanocrystals. We will also provide an up-to-date overview of their postsynthesis transformations, and summarize the various solution processes that are aimed at fabricating halide perovskite-based nanocomposites. Furthermore, we will review the fundamental optical properties of halide perovskite nanocrystals by focusing on their linear optical properties, on the effects of quantum confinement, and on the current knowledge of their exciton binding energies. We will also discuss the emergence of nonlinear phenomena such as multiphoton absorption, biexcitons, and carrier multiplication. Finally, we will discuss open questions and possible future directions.

Haplotype-resolved and integrated genome analysis of the cancer cell line HepG2.

HepG2 is one of the most widely used human cancer cell lines in biomedical research and one of the main cell lines of ENCODE. Although the functional genomic and epigenomic characteristics of HepG2 are extensively studied, its genome sequence has never been comprehensively analyzed and higher order genomic structural features are largely unknown. The high degree of aneuploidy in HepG2 renders traditional genome variant analysis methods challenging and partially ineffective. Correct and complete interpretation of the extensive functional genomics data from HepG2 requires an understanding of the cell line's genome sequence and genome structure. Using a variety of sequencing and analysis methods, we identified a wide spectrum of genome characteristics in HepG2: copy numbers of chromosomal segments at high resolution, SNVs and Indels (corrected for aneuploidy), regions with loss of heterozygosity, phased haplotypes extending to entire chromosome arms, retrotransposon insertions and structural variants (SVs) including complex and somatic genomic rearrangements. A large number of SVs were phased, sequence assembled and experimentally validated. We re-analyzed published HepG2 datasets for allele-specific expression and DNA methylation and assembled an allele-specific CRISPR/Cas9 targeting map. We demonstrate how deeper insights into genomic regulatory complexity are gained by adopting a genome-integrated framework.

Genome amplification and cellular senescence are hallmarks of human placenta development.

Genome amplification and cellular senescence are commonly associated with pathological processes. While physiological roles for polyploidization and senescence have been described in mouse development, controversy exists over their significance in humans. Here, we describe tetraploidization and senescence as phenomena of normal human placenta development. During pregnancy, placental extravillous trophoblasts (EVTs) invade the pregnant endometrium, termed decidua, to establish an adapted microenvironment required for the developing embryo. This process is critically dependent on continuous cell proliferation and differentiation, which is thought to follow the classical model of cell cycle arrest prior to terminal differentiation. Strikingly, flow cytometry and DNAseq revealed that EVT formation is accompanied with a genome-wide polyploidization, independent of mitotic cycles. DNA replication in these cells was analysed by a fluorescent cell-cycle indicator reporter system, cell cycle marker expression and EdU incorporation. Upon invasion into the decidua, EVTs widely lose their replicative potential and enter a senescent state characterized by high senescence-associated (SA) ß-galactosidase activity, induction of a SA secretory phenotype as well as typical metabolic alterations. Furthermore, we show that the shift from endocycle-dependent genome amplification to growth arrest is disturbed in androgenic complete hydatidiform moles (CHM), a hyperplastic pregnancy disorder associated with increased risk of developing choriocarinoma. Senescence is decreased in CHM-EVTs, accompanied by exacerbated endoreduplication and hyperploidy. We propose induction of cellular senescence as a ploidy-limiting mechanism during normal human placentation and unravel a link between excessive polyploidization and reduced senescence in CHM.

One thousand somatic SNVs per skin fibroblast cell set baseline of mosaic mutational load with patterns that suggest proliferative origin.

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.