INTRODUCTION

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, are hallmarks of aging (1, 2). However, thus far, of all mutation types, only chromosomal alterations could readily be studied directly during in vivo aging using cytogenetic methods (3). Because of their small size, random nature, and low abundance, most somatic mutations are difficult to detect, except in single cells or in clonal lineages (4). In the past, using transgenic reporters, mutations have been found to accumulate with age in a tissue-specific manner (5). However, this approach does not allow a genome-wide, direct analysis of somatic mutations in human primary cells. More recently, using single-cell whole-genome sequencing (WGS), somatic mutations were found to accumulate with age in human neurons (6) and B lymphocytes (7). Others also reported increased somatic mutations in human primary cells isolated from intestine, colon, and liver, albeit in clones propagated from human tissue-specific stem cells (8), which may not be representative of the differentiated cells that ultimately provide tissue function. Nevertheless, together, these studies confirmed that mutations in different somatic cell types of humans accumulate with age.

Here, we present single-cell genome-wide somatic mutation profiles of differentiated human liver hepatocytes as compared with adult liver stem cells (LSCs). Human liver is of particular interest for studying genome instability because of its high metabolic activity and its role in detoxification of xenobiotics, which makes this organ the most important target for genotoxicity in the body. In humans, accumulation of de novo mutations could contribute to the observed age-related loss of liver function, most notably a severe reduction in metabolic capacity, and multiple pathologies, including fatty liver disease, cirrhosis, hepatitis, infections, and cancer (9, 10). Our results indicate high spontaneous mutation frequencies in differentiated hepatocytes that significantly increase with age. By contrast, mutation frequencies in adult LSCs, defined as the cells that give rise to clonal outgrowths, were fairly low. In differentiated hepatocytes, a considerable number of mutations were found in functional parts of the genome. These results indicate that the human liver is subject to a high burden of genotoxicity and that adult stem cells are a critical component in maintaining overall genome integrity within a tissue.

The quantitative detection of de novo somatic mutations in single cells after whole-genome amplification (WGA) and WGS remains a challenge because of the high chance of errors. Here, we used a well-validated, highly accurate method, single-cell multiple displacement amplification (single-cell MDA, or SCMDA) (11), to analyze somatic mutations in single primary hepatocytes from human donors varying in age between 5 months and 77 years. These cells were isolated shortly after death through perfusion of whole livers from healthy human individuals after informed consent by the donors family (Lonza Walkersville Inc.). Cell viability was higher than 80% and, after Hoechst staining, individual, diploid hepatocytes were isolated via fluorescence-activated cell sorting (FACS) into individual polymerase chain reaction (PCR) tubes (fig. S1A). In total, we sequenced four single hepatocytes and bulk genomic liver DNA for each of 12 human donors (table S1). Each cell was subjected to our recently developed procedure for WGA and WGS (11, 12). Somatic single-nucleotide variants (SNVs) in single cells were identified relative to bulk genomic DNAs at a depth of 20 using VarScan2, MuTect2, and HaplotypeCaller with certain modifications (Materials and Methods and table S2). Overlapping mutations from this tricaller procedure were exclusively considered for further analysis. The results were essentially confirmed by using two alternative variant callers: SCcaller (11) and LiRA (Linked Read Analysis) (13).

After adjusting for genomic coverage, the number of SNVs per cell for 48 hepatocytes from 12 donors was found to vary between 357 and 5206 with four extreme outliers of 20,557 to 37,897 SNVs per cell excluded from the statistical model (Fig. 1A and table S2). The number of mutations per cell was found to increase significantly with the age of the donor (P = 1.22 109), with median values of 1222 855 SNVs per cell in the young group (36 years, n = 21 cells), and 4054 1168 SNVs per cell in the aged group (46 years, n = 23 cells), excluding the four outliers (Fig. 1A). The median number of mutations per cell in hepatocytes from the youngest donor was in the same range as what we recently reported for primary human fibroblasts from young donors, i.e., 1027 and 926 SNVs per cell from the 5-month-old and 6-year-old donors, respectively (11, 12). However, during aging, mutation levels increased over the same age range up to 2.5 times higher than in our previously analyzed human B lymphocytes (7) or human neurons analyzed by others (6) (fig. S2A).

(A) SNV levels in individual differentiated hepatocytes. The y axis on the left indicates the number of mutations per cell, and the y axis on the right indicates mutation frequency per base pair. The median values with SDs among four cells of each subject are indicated. Data indicate an exponential increase in mutation frequency with donor age (R = 0.892, P = 1.16 106). bp, base pair. (B) SNV levels in LSC-derived parent clones (red) and their kindred cells (light green) from three young donors. The Venn diagrams indicate the fraction of SNVs detected in the parent clones (collectively for each individual; n = 3) that were also detected in the kindred LSCs. The bars indicate the median mutation frequencies in clones (red) and kindred single cells (light green). (C) Comparison of SNV levels in differentiated hepatocytes (dark green dots; n = 24 from six donors) and LSCs (light green; n = 10 from three donors), all within the young donor group 36 years. Mutation frequencies were corrected for the estimated number of cell divisions. (D) SNV levels in LSCs and differentiated hepatocytes from the same participants, corrected for the estimated number of cell divisions.

At this stage, we were interested in the possible cause of the high mutation frequencies in the four outlier cells. Three of the four outliers with the highest SNV levels revealed multiple mutations in genes involved in DNA repair (table S3) (14), which could conceivably underlie the observed accelerated mutation accumulation in these cells. Of note, individual outlier cells with high mutation levels have been detected in other tissues (6, 7).

Together, these findings indicate that the liver is prone to high levels of de novo somatic mutations, which could possibly be related to its major role in the metabolization and detoxification of xenobiotics.

The mutation frequencies observed in human hepatocytes from older subjects were higher than those previously found in human neurons and B lymphocytes (6, 7). They were also higher than the mutation frequencies reported for stem cellderived liver organoids (fig. S2B) (8). It is critically important to validate the results obtained with single-cell mutation analysis to rule out possible amplification artifacts. In our previous studies on human primary fibroblasts, we validated single-cell data by also analyzing unamplified DNA from clones derived from cells in the same population (11). Here, we generated liver-specific clones from young donors by plating the prepurified hepatocyte cell suspensions in selective medium for LSC expansion (Materials and Methods). Under these conditions, the differentiated hepatocytes died within 5 to 7 days, while the residential LSCs could be propagated without differentiation. The latter was confirmed using biomarker analysis (Materials and Methods and fig. S1B) (15, 16). In addition, we obtained from a commercial source one sample of human postnatal LSCs from a 1-year-old donor at passage 9 (approximately 27 population doublings), which were expanded and also grown into clones in the same way.

LSC clones could be established only from young individuals, i.e., hepatocyte samples from the 1-year-old, 5-month-old, and 18-year-old participants. This is in keeping with observations that resident stem cell properties change with age, with a general reduction in proliferative capacity and increased cellular senescence (17).

Both LSC clones and kindred single cells derived from the young individuals were processed and subjected to WGS, as described above for differentiated hepatocytes. We then tested for the fraction of mutations called in the clones that were also found in the single cells derived from them. As shown in the Venn diagrams (Fig. 1B and fig. S3, A and B), most of these mutations were indeed confirmed in the single cells. This is very similar to what we previously reported for human single fibroblasts and clones derived from the same population of cells (11), which underscores the validity of our single-cell mutation detection method, also in liver cells. Of note, most of the mutations found in the single cells, but not in their parental clones, are likely to be also real. These are likely either mutations missed during variant calling in the clone or de novo mutations arising in the individual cells during clone culture and expansion.

Once we confirmed the validity of our single-cell data, we directly compared mutation frequencies between the single cells defined as LSCs and differentiated hepatocytes, both from the young donor group. Previous studies have provided evidence for lower spontaneous mutation frequencies in stem as compared with differentiated cells (18, 19). For this comparison to be valid, we compared mutation frequencies per cell division in both cell types. This was necessary because the number of cell divisions is a major factor in causing base substitution mutations through replication errors. We first estimated the number of cell divisions that had occurred in human somatic cells of the young age group since the zygote, as described previously (20) (Materials and Methods). We then added, only to the LSCs, the estimated additional numbers of cell divisions during culture (Materials and Methods). The results show that, on a per cell division basis, somatic mutation frequencies were indeed lower in the LSCs than in the differentiated hepatocytes (about twofold), i.e., 11 SNVs versus 21 SNVs per cell per mitosis, respectively (P = 1.26 104, two-tailed Students t test) (Fig. 1C and table S2). A reduced mutation rate in LSCs could explain the fairly modest age-related increase reported previously for stem cellderived organoids (figs. S2B and S3C) (8). The tendency of differentiated hepatocytes to accumulate mutations to a much higher level than stem cells is further confirmed by the significantly higher cell-to-cell variation among the former (P = 1.42 103, Levenes test; Fig. 1, C and D). These observations are in keeping with the idea that stem cells are superior to differentiated cells in preserving their genome integrity, possibly through an enhanced capability to prevent or repair DNA damage (21, 22).

Next, we analyzed the mutational spectra in LSCs and differentiated hepatocytes. In differentiated hepatocytes, the most common mutation types were GC-to-AT transitions and GC-to-TA transversions (Fig. 2A and fig. S4, A and B). These mutations are known to be induced by oxidative damage (23), which itself has often been considered as a main driver of aging and age-related diseases (24). However, the most rapidly increased mutation type with age was the AT-to-GC transition (P = 2.16 1010, two-tailed Students t test; table S4 for Pearsons 2 test). This mutation can be caused by mispairing of hydroxymethyluracil (5-hmU), another common oxidative DNA lesion. Alternatively, AT-to-GC mutations are induced by mutagenic alkyl-DNA adducts formed as a result of thymine residue alkylation (25, 26). Notably, certain minor alkyl-pyrimidine derivatives can escape repair, accumulate during aging, and lead to mutations much later (26, 27).

(A) Relative contribution of the indicated six mutation types to the point mutation spectrum for the five indicated liver sample groups. Data are represented as the mean relative contribution of each mutation type in sample groups of young and aged differentiated hepatocytes (21 cells from six donors 36 years, and 23 cells from six donors 46 years), adult LSC-derived parent clones and their kindred single cells separately, and a group of outlier cells (n = 4). (B) Three mutational signatures (L1, L2, and L3) were de novo identified by non-negative matrix factorization analysis from the somatic mutations in the different groups in (A). (C) Contributions of signatures L1, L2, and L3 to all SNVs in young and aged hepatocytes, young LSCs, and outlier cells.

Mutation spectra of the LSCs and LSC clones revealed a lower fraction of GC-to-AT transitions as compared with differentiated hepatocytes from the young group (Fig. 2A and figs. S3D and S4, A and B). This could be due to the virgin state of these cells, not participating in metabolizing xenobiotics, which is associated with oxidative DNA damage. However, we cannot rule out that, instead, the altered spectrum is related to in vitro culturing, which may alter the ratio of GC-to-AT transitions and GC-to-TA transversions. In the human LSCs derived from clones, the relative frequency of the GC-to-AT transition mutations is slightly, albeit significantly, increased as compared with the parent clones themselves (P = 7.43 104, two-tailed Students t test; table S4 for Pearsons 2 test; Fig. 2A and fig. S4A). Kindred single LSCs, which were derived from parent LSC clones, representing the original LSCs, have undergone multiple rounds of cell division with ample opportunity for replication errors, for example, as a consequence of ambient oxygen to which these cells have been inevitably exposed during subculture. Hence, this would suggest that cell culture has the opposite effect of what we observed from the stem cell versus differentiated cell difference, i.e., increasing rather than decreasing the fraction of GC-to-AT transitions.

To analyze mutation spectra more precisely, we performed non-negative matrix factorization (Materials and Methods) to extract three de novo mutation signatures, signatures L1, L2, and L3, from the mutation spectra of the four groups of human liver cells analyzed, i.e., combined LSCs and clones collectively, differentiated hepatocytes from young participants, differentiated hepatocytes from aged participants, and the four combined outlier cells. We compared these signatures to the COSMIC (Catalogue Of Somatic Mutations in Cancer) signatures described for various human tumors (Fig. 2B and table S5). Signature L1 substantially increased in differentiated hepatocytes from the aged group as compared with hepatocytes and LSCs from young individuals (Fig. 2C). This signature highly correlated with the liver-specific and age-associated mutation signature A dominant in human organoids of liver-specific origin in the aforementioned organoid study (8), as well as with COSMIC signature SBS5, strongly associated with aging (fig. S4C and table S5) (28, 29). Signature L2, with its increased level of oxidative GC > TA transversions, dominated the mutation spectrum of both LSCs and differentiated hepatocytes from young donors (Fig. 2C) and was significantly reduced in cells from the aged donors. Signature L2 highly correlated with COSMIC signatures SBS18 and SBS36, known to be associated not only with oxidative stress (fig. S4C and table S5) but also with proliferation signature C (table S5), found in all in vitro propagated cell types in the aforementioned organoid study (8). Since this signature was dominant in the LSCs, it possibly reflects the stem/progenitor-like origin of hepatocytes and remains dominant in differentiated hepatocytes of the young individuals (Fig. 2C). Signature 3, dominant in the outlier cells, highly correlated with COSMIC signature SBS5, the aging signature, but also correlated with SBS6 and SBS1, signatures associated with DNA mismatch repair deficiency (29).

The above analysis was confirmed when we, instead of extracting de novo signatures from our four groups of liver cell mutation spectra, tested which of the reference COSMIC signatures could be found in these groups (fig. S4C).

Next, we analyzed the distribution of the somatic mutations in human liver cells across the genome. After pooling all mutations of the 21 differentiated cells from the young and the 23 differentiated cells from the old individuals, excluding the four outliers, the large majority of mutations distributed randomly across the genome in both groups (Fig. 3A). We then tested the possibility that during aging, mutations in functionally relevant sequences were selected against, as we previously observed for age-related mutation accumulation in B lymphocytes (7). Here, the functional liver genome was defined as the transcribed liver exome, using available data on gene expression levels in 175 previously described total liver samples [Genotype-Tissue Expression (GTEx) Consortium] (30), and its regulatory regions, identified as promoters of active genes or open chromatin regions, e.g., transcription factor binding regions, identified by ATAC (Assay for Transposase-Accessible Chromatin) sequencing in total liver tissue (ENCODE) (31). Of note, since the databases used were from whole liver, these definitions would not necessarily apply to LSCs or other subpopulations. However, it is reasonable to assume that whole liver is a good surrogate even for those fairly rare liver-specific cells.

(A) Circos diagram of genomic SNV distribution in four groups: pooled LSCs, young and aged hepatocytes, and outlier cells. (B) SNV levels in the functional genome and genome overall in differentiated hepatocytes (left) and in LSCs (right) as a function of age. Each data point represents the ratio of the number of mutations per cell to the median number of mutations of the four cells from the 5-month-old subject. Mutations in the functional genome are shown in red and those in the genome overall in blue. (C) Mutation frequency per base pair in the transcribed part of the liver genome (red) and the nontranscribed part (blue) in differentiated hepatocytes (left) and LSCs (right) as a function of age.

The ratio of total to functional SNVs in differentiated hepatocytes was found to remain about 1 across the different age levels (P = 0.5134, Wilcoxon signed-rank test, two tailed) (Fig. 3B), indicating no selection against deleterious somatic mutations in low-proliferating hepatocyte populations during aging. By contrast, the same ratio in pooled adult LSCs was about 2 and significantly different from that in differentiated hepatocytes (P = 5.34 104, Wilcoxon signed-rank test, two tailed). This suggests selection against deleterious mutations during the cell proliferation cycles that gave rise to these stem cells. It also suggests that LSCs may have an increased capacity to protect their genome simply by remaining quiescent. We also compared mutation frequencies in transcribed versus untranscribed liver cell genes. Transcribed liver genes were defined as genes with expression values 1 transcripts per kilobase per million (TPM), while nontranscribed genome included all sequences with expression values <1 TPM in liver tissue (GTEx) (30). The results indicated a significantly lower number of SNVs affecting transcribed liver genes than nontranscribed genes across all donor ages (P = 7.21 108, Wilcoxon signed-rank test, two tailed) as well as in the LSCs and clones (P = 7.63 106, Wilcoxon signed-rank test, two tailed) (Fig. 3B), suggesting active transcription-coupled repair in normal human liver (32).

Somatic mutations have long been implicated as a cause of aging (33, 34). However, thus far, it has not been possible to test this hypothesis directly because of a lack of advanced methods to analyze random somatic mutagenesis in vivo, which requires high-throughput sequencing of single cells. Using our advanced single-cell sequencing method, we show that the number of somatic base substitution mutations in normal human liver significantly increases with age, reaching as much as 3.3 times more mutations per cell in aged humans than in young individuals. Of note, the numbers of mutations in aged liver are significantly higher than what has previously been reported for aged human liver organoids (fig. S2B) (8) and also higher than recent results reported for aged human neurons (fig. S2A) and B cells (7). Since we essentially ruled out that many of these mutations are artifacts of the amplification system, the most likely cause of this high mutagenic activity in the human liver is the high metabolic and detoxification activity in this organ, which is known to be associated with genotoxicity (35).

Out of 48 hepatocyte cells analyzed, 4 cells revealed extremely elevated mutation loads, over 10 times exceeding SNV levels in age-matched normal hepatocytes even from the same subject. These outliers have also been observed in the only two studies of somatic mutations in human tissues in vivo using a single-cell WGS approach (6, 7). Of the four outliers observed in this present study, multiple de novo SNVs were found to reside in DNA repair genes, strongly suggesting that these mutations were responsible for mutator phenotypes similar to what has been shown for cancers (36). While we cannot know when the mutations that gave rise to rapid mutation accumulation in these cells occurred, this may have been fairly recently, with imminent death of the cells likely. On average, almost 60 nonsynonymous mutations in the functional exome of these cells were found, suggesting a likely functional effect (table S6). However, since we could not longitudinally follow mutation loads in the same single cells, our data do not allow any conclusions on the cause and effect of the observed mutations.

Somatic mutation frequencies in normal differentiated hepatocytes were found to be much higher than in residential LSCs. This means that in vitro clonal surrogates for cells do not always accurately represent the mutation loads of in vivo differentiated cells, which makes predictions of a functional impact of somatic mutations from these clonal data difficult. While we do not know the mechanism(s) of reduced spontaneous mutation loads in stem as compared with differentiated cells, such evidence has also been reported by others (18, 19), and it is possible that stem cells have superior genome maintenance systems as compared with their differentiated counterparts. However, a caveat in this respect is that the LSCs that we enriched for may not in fact be the LSCs giving rise to most of the differentiated hepatocytes. Hence, we cannot be sure that a direct comparison between a stem cell and differentiated cells derived from this stem cell was in fact made.

Another important question is the possible functional impact of random somatic mutagenesis on the aging phenotype. While from our current data we cannot conclude direct cause-and-effect relationships, our observation that the functional part of the genome accumulated numerous mutations suggests that aging-related cellular degeneration and death could at least, in part, be due to somatic mutations. While the occurrence of no more than 11 nonsynonymous mutations in the transcribed exome of liver hepatocytes from humans in their 70s suggests a minor contribution of changes in the protein-coding part of the genome, the well over 100 de novo mutations in gene regulatory sequences may point toward an important role for stochastic gene expression changes in age-related loss of organ function and increased disease incidence. These mutations could possibly increase transcriptional noise, a molecular phenotype that appears characteristic for cells from aged individuals (3739).

Last, while in our current work only base substitution mutations were analyzed, other types of mutations are likely to occur as well. The frequency of most of these mutations, e.g., small insertions and deletions, copy number variation, and genome structural variation, is likely to be much lower than the frequencies of base substitutions observed to rise to thousands of mutations per cell. However, their effects are possibly much larger since they affect a larger part of the genome and, when in exomes, almost always lead to loss of function. It is conceivable that, taken together, de novo mutations could have serious effects on the function of human somatic cells in vivo above and beyond their causal relevance in liver cancer.

Frozen human hepatocyte samples were purchased from Lonza Walkersville Inc. Whole livers for hepatocyte isolation were obtained with the informed consent of families of registered organ donors. The obtained liver organs were rejected for transplant due to either lack of a donor match or morphological alterations (e.g., tearing and hematoma). All 12 selected hepatocyte donors were healthy participants of various age, gender, and ethnicity (table S1) without any liver cancer or other liver pathology history. These cells had been isolated using a gold standard, two-step liver/liver lobe perfusion procedure. Cells were suspended in 2 to 5 ml of media and counted with Trypan blue to estimate viability (higher than 80%), and frozen in dimethyl sulfoxide/liquid nitrogen (www.lonza.com). One specimen of frozen human neonatal LSCs from a 1-year-old donor was purchased from Kerafast Inc. (www.kerafast.com). These cells had been derived by the Sherley laboratory (Boston, MA, USA) and characterized to confirm their stem cell identity (4042).

After thawing, hepatocyte suspensions were used to collect single hepatocytes into individual 0.2-ml PCR tubes with 2.5 l of phosphate-buffered saline (PBS) by means of FACS (FACSAria, Becton Dickenson). Selection of the target hepatocyte population was based on the large cell size of hepatocytes (forward-scatter/side-scatter parameters) along with the additional fluorescence staining for DNA content and cell viability. Briefly, bulk hepatocyte suspension samples were prior stained according to the manufacturers protocol with the viable DNA-binding dye Hoechst 33342 (Life Technologies) to discriminate cells with a standard diploid chromosome set and LIVE/DEAD Cell Vitality Assay Kit C12 Resazurin/SYTOX Green (Thermo Fisher Scientific) to select viable healthy cells. Typical FACS layout is shown in fig. S1A. Upon sorting, tubes with single cells were frozen on dry ice and kept at 80C until use.

Neonatal LSCs of passage 9 (one passage corresponds to approximately three cell population doublings for these cells according to the manufacturers protocol) from the 1-year-old donor were purchased from Kerafast Inc. The commercial LSCs were cultured in polarization media [Dulbeccos modified Eagles medium, 10% dialyzed fetal bovine serum (Invitrogen), 1.5 mM xanthosine (Sigma), 1 penicillin/streptomycin, epidermal growth factor human (20 ng/ml; Invitrogen), transforming growth factor human recombinant (0.5 ng/ml Sigma)] according to the manufacturers protocol (Kerafast Inc.) (4042). These cells served as controls to characterize de novo isolated and polarized LSCs.

Additional LSC cultures were isolated and polarized and characterized from the bulk commercial hepatocyte suspensions (Lonza Walkersville Inc.) from young donors using previously described protocols with specific modifications (15, 16) combined with the aforementioned Kerafast protocol for neonatal LSCs. Briefly, bulk suspension hepatocytes (0.5 106 to 1 106 of cells) were transferred to polarization media as described for the neonatal LSCs and cultured on cell-adhesive 12-well plates for 5 to 7 days. Then, all nonattached hepatocytes were removed, and fresh media were added to the small remaining population of attached progenitor cells. After 1 to 1.5 weeks of culture and media changes, attached cells symmetrically divided, growing to mixed clonal populations of polarized adult LSCs. These cultures were frozen at early passage (p = 3 to 5) until further use. Only LSCs from donors of younger age (22 years) could be isolated in this way.

Phenotypes of the polarized cells were analyzed for the presence of specific surface stem cell and epithelial progenitor cell epitopes, e.g. EpCAM (epithelial cell adhesion molecule), Lgr5, CD90, CD29, CD105, and CD73, upon staining with antibodies by means of multicolor flow cytometry analysis (LSRII, Becton Dickinson) as recommended previously (15, 16, 43, 44). Characteristic FACS profiles and specific phenotypes for commercial LSCs (control) and two manually isolated and polarized LSC lineages are shown in fig. S1B.

Single-cell derived parent clones and their kindred single cells were prepared and collected using CellRaft arrays (Cell Microsystems) as described previously (11). Briefly, an LSC suspension was plated on a CellRaft array consisting of 12,000 individual portable rafts for single cells at the required density of 5000 cells per array. After 4 to 8 hours, individual LSCs were elongated and attached to the array surface locating on individual rafts. After attachment, the medium with floating cells was replaced, and single-cell positions were marked and tracked during the following 7 to 10 days to detect dividing cells and growing individual single-cell derived clones. Once the colony/clone reached confluence on the raft (8 to 10 cells per raft), it was dislocated from the array with a positioned automatic needle and transferred with a magnetic wand to a 96-well plate. Upon reaching confluence, single-cell derived clones were trypsinized and subsequently transferred to 24-well plates, then 12-well plates, 6-well plates, and, lastly, 10-cm plates to reach a total amount of 1.5 106 to 3 106 cells per parent clone. Together, the process of establishing a clone from a single cell took about 25 to 30 days.

Individual single cells from the parent clones were collected, also using CellRafts, and transferred to a 0.2-ml PCR tube containing 2.5 l of PBS. The presence of a single raft was observed under a magnifying glass. Upon single-cell collection, tubes were fast frozen on dry ice and kept on 80C until further use.

Single hepatocytes from each subject were subjected to WGA using our modified procedure of low-temperature cell lysis and DNA denaturation followed by MDA as described (11). As positive and negative controls for WGA, we used 1 ng of human genomic DNA and DNA-free PBS solution, respectively. Resultant MDA products were purified using AMPureXP beads (Beckman Coulter), and the amplified DNA concentration was measured with the Qubit High Sensitivity dsDNA kit (Invitrogen Life Sciences). To verify sufficient and uniformly amplified single-cell MDA products, we performed the eight-target locus-dropout test as described previously (11). Selected confirmed samples (four single-cell MDA products per subject) were further subjected to library preparation and WGS.

Human bulk genomic DNA was collected from total cell suspensions using the DNeasy Blood & Tissue Kit (Qiagen) according to the manufacturers protocol. LSC clonederived DNA was extracted from clones of at least 1.5 106 to 2.5 106 cells in a similar way. DNA concentration was quantified with the Qubit High Sensitivity dsDNA kit (Invitrogen Life Sciences), and DNA quality was evaluated by 1% agarose gel electrophoresis.

The libraries for Illumina next-generation WGS were generated from 0.2 to 0.4 g of genomic DNA, clone-derived bulk DNA, and single-cell MDA DNA human samples using the NEBNext Ultra II FS DNA Library Prep Kit for Illumina (New England BioLabs). The libraries were sequenced with 2 350base pair paired-end reads on an Illumina HiSeq X Ten sequencing platform by Novogene Inc.

Next-generation WGS at a minimal depth of 20X base coverage was performed on four individual mature hepatocytes per human subject (12 human subjects, 48 single cells in total) (table S2). Bulk DNA from two or three LSC-derived clones and MDA products from three to four corresponding kindred single cells per donor (three donors, eight parent clones, and 10 kindred single LSCs) were sequenced similarly.

For all samples, adapter and low-quality reads were trimmed by Trim Galore (version 0.3.7). Quality checks were performed before and after read trimming by FastQC (version 0.11.4). The trimmed reads were aligned to the human reference genome (GRCh37 with decoy) by BWA mem (version 0.7.10) (45). Duplications were removed using samtools (version 0.1.19) (46). The known indels and single-nucleotide polymorphism (SNPs) were collected from the 1000 Genomes Project (phase 1) and Single Nucleotide Polymorphism Database (dbSNP) (build 144). Then, the reads around known indels were locally realigned, and their base quality scores were recalibrated on the basis of known indels and SNVs, both via the Genome Analysis Toolkit (GATK, version 3.5.0) (47).

Somatic mutations between each single cell and the corresponding bulk and between each clone and corresponding bulk were identified using three different variant callers: VarScan2 (48), MuTect2 (49), and HaplotypeCaller (47). To obtain high-quality mutation calls and avoid high false-positive rates in individual callers, we applied a comprehensive procedure in filtering. First, we only considered mutations on autosomes. Then, we considered mutations with a GATK phred-scaled quality score of at least 30 and excluded mutations overlapping with known SNPs from dbSNP. Furthermore, we required a minimum base depth of 20X and filtered mutations with variant-supporting reads in bulk. Moreover, mutations present in at least two cells in each individual were also removed to further exclude potential germline mutations. The mutations present in all three variant callers were considered as true de novo mutations. Last, considering that amplification errors and/or nonuniform coverage could induce false-positive mutations in no more than one-eighth of the reads, we used a binomial distribution to filter these potential false-positive mutations, which excluded most mutations present in 25% of the reads or less. To further check the power of the used pipeline in filtering amplification errors, we also called the somatic mutations using our alternative, the SCcaller tool (11) and the LiRA pipeline (13) (figs. S2A and S3B).

The frequency of somatic SNVs per cell was estimated after normalizing genomic coveragefrequency of somatic SNVs per cell=#somatic SNVssurveyed genometotal size of genome

As the reads were aligned to the haploid reference genome, the frequency of somatic SNVs per base pair was calculated by dividing the frequency of somatic SNVs per cell by genome size and ploidy of the genome (ploidy = 2)frequency of somatic SNVs per base pair=frequency of somatic SNVs per celltotal size of genome*ploidy of genome

The surveyed genome per single cell/clone was calculated as the number of nucleotides with read mapping quality 20 and position coverage 20X.

The outliers of the hepatocytes were defined using Tukeys range test: Four cells were defined as extreme outliers as their frequencies were higher than Q3 + 3 * IQR, where Q3 is the third quartile of the frequencies and IQR is the interquartile range. The outlier cells were excluded from the statistical model.

For the LSC-differentiated hepatocyte comparison, the absolute de novo mutation frequencies were corrected for the number of cell divisions undergone since the zygote (table S2). We used 45.1 as the number of developmental mitoses (20) and assumed a subsequent turnover rate of one cell division per year, based on empirical evidence from rodents (50, 51). In total, 45.5, 46.3, and 61.6 cell divisions were estimated for both LSCs and differentiated hepatocytes from 5-month-old, 1-year-old, and 18-year-old individuals, respectively. For LSCs from 5-month-old, 1-year-old, and 18-year-old individuals, we then added, respectively, an estimated 33, 41.7, and 33 cell divisions during the enrichment process of stem cells, and 21.9, 24.5, and 21.9 cell divisions associated with clonal outgrowth of the single LSCs.

To determine the overlap between SNVs called in the clones and the single cells derived from them, genome coverage in the clone was normalized to that in its kindred single cell. Mutations found in a single cell and appearing in at least 1 read in the parent clone were considered as overlapping. When there were no variant-supporting reads in the clone, the mutation was determined as kindred cell specific. This assignment left some mutations with an unknown status more likely to be de novo mutations arising in the individual cells during clone culture and expansion.

The identified mutations in all individuals were pooled into four groups: LSC cells/clones from young donors, hepatocytes from young and aged donors, and outlier hepatocytes. The integrated spectra of six mutation types in each group were plotted using the R package MutationalPatterns (52). Using non-negative matrix factorization (NMF) decomposition in the same package, we revealed group-specific mutational signatures as well as de novo identified three signatures in normal human liver cells. To identify the potential origin of the mutational spectra, the group mutational signatures and newly revealed signatures to the published signatures associated with liver-specific organoids and various cancer tissues. Three tissue-specific organoid signatures were obtained from a recent study (8); 67 cancer mutation signatures were downloaded from the latest version 3 of the COSMIC database (https://cancer.sanger.ac.uk/cosmic/signatures/SBS/) (28, 29). The cosine similarity between newly identified and published signatures was calculated for comparisons (table S5).

All reported mutations were annotated based on the gene definitions of GRCh37.87. Mutations were further extracted from the functional genome, including transcribed genes, promoters, and open chromatin regions. The nonsynonymous and synonymous mutations were identified by analysis of variance (ANOVA) (53), while damaging and tolerated mutations were checked by SIFT (54) and PROVEAN (55). When damaging (Sorting Intolearnt From Tolerant, SIFT) or deleterious (Protein Variation Effect Analyzer, PROVEAN), the mutation was marked as damaging, and when tolerated (SIFT) and neutral (PROVEAN), a tolerated mutation.

The open chromatin regions were identified by ENCODE transcription factor binding regions in whole genome and ATAC sequencing data in the functional genome in liver tissue samples. Raw ATAC sequencing data were downloaded from ENCODE (experiment name: ENCSR373TDL) (31). The adapter and low-quality ATAC sequencing reads were filtered using Trim Galore (version 0.3.7). Clean reads were aligned to the human reference genome (GRCh37) with Bowtie2 (version 2.2.3; option: -X 2000). Duplicated reads were removed with the Picard tool (version 1.119). Open chromatin regions were determined by MACS2 (version 2.1.1; option: callpeak -g hs --nomodel --shift 100 --extsize 200) (56).

Gene expression levels for total human liver tissue were obtained from GTEx (https://gtexportal.org/) (30). We defined the transcribed genes as those with expression level 1 TPM in all samples. Also, we separated the transcribed and nontranscribed genome by TPM 1 and < 1 in all samples, respectively.

The rest is here:
Single-cell analysis reveals different age-related somatic mutation profiles between stem and differentiated cells in human liver - Science Advances