Accumulating somatic mutations have been implicated in age-related cellular degeneration and death. Because of their random nature and low abundance, somatic mutations are difficult to detect except in single cells or clonal lineages. Here we show that in single hepatocytes from human liver, an organ normally exposed to high levels of genotoxic stress, somatic mutation frequencies are high and increase substantially with age. Significantly lower mutation frequencies were observed in liver stem cells and organoids derived from them. These results could explain the increased age-related incidence of liver disease in humans and stress the importance of stem cells in maintaining genome integrity.Genome integrity is critically important for cellular function. Evidence has accumulated that loss of genome integrity and the increasingly frequent appearance of various forms of genome instability, from chromosomal aneuploidy to base substitution mutations, is a hallmark of aging (1, 2). However, somatic mutations occur de novo across cells and tissues and are of low abundance. Hence, thus far, of all mutation types only chromosomal alterations could readily be studied directly during in vivo aging using cytogenetic methods (3). More recently it has become possible to also study low-abundant somatic mutations by whole genome sequencing (WGS) of single cells or clones derived from single cells (4-7). Such studies confirmed that like chromosomal alterations also base substitution mutations accumulate with age in human brain (8,9), intestine and liver (4).Due to its high metabolic and detoxification activity liver is more likely to be adversely affected by various damaging agents than most other organs. In humans this may cause agerelated loss of function, most notably a severe reduction in metabolic capacity and multiple pathologies, including fatty liver disease, cirrhosis, hepatitis, infections and cancer (10,11). To test if somatic mutations occur with a high frequency and accumulate rapidly with age we performed WGS on single primary hepatocytes from human donors varying in age between 5 months and 77 years of age. These cells were isolated from healthy human individuals shortly after death through perfusion of whole donor livers (Lonza Walkersville Inc.). Cell viability was higher than 80% and, after staining with Hoechst, hepatocytes were isolated via FACS into individual PCR tubes as diploid single cells (Fig. S1). In total we sequenced 4 single hepatocytes and bulk genomic liver DNA for each of eight human donors (Table S1). Each cell was subjected to our recently developed procedure for whole genome amplification (WGA) and sequencing (5,6). Somatic single nucleotide variants (SNVs) in single cells were identified relative to bulk genomic DNAs at a depth of ≥ 20X using a combination of VarScan2, MuTect2 and Haplotypecaller software with certain modifications (Methods , Table S2). Overlapping mutations, revealed by all three analytical tools, were exclusively considered for further analysis.After adjusting for genomic ...