Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Precancerous liver diseases do not cause increased mutagenesis in liver stem cells

## Abstract

Inflammatory liver disease increases the risk of developing primary liver cancer. The mechanism through which liver disease induces tumorigenesis remains unclear, but is thought to occur via increased mutagenesis. Here, we performed whole-genome sequencing on clonally expanded single liver stem cells cultured as intrahepatic cholangiocyte organoids (ICOs) from patients with alcoholic cirrhosis, non-alcoholic steatohepatitis (NASH), and primary sclerosing cholangitis (PSC). Surprisingly, we find that these precancerous liver disease conditions do not result in a detectable increased accumulation of mutations, nor altered mutation types in individual liver stem cells. This finding contrasts with the mutational load and typical mutational signatures reported for liver tumors, and argues against the hypothesis that liver disease drives tumorigenesis via a direct mechanism of induced mutagenesis. Disease conditions in the liver may thus act through indirect mechanisms to drive the transition from healthy to cancerous cells, such as changes to the microenvironment that favor the outgrowth of precancerous cells.

## Introduction

Liver cancer is the fifth most common cancer worldwide, causing around 720,000 deaths each year1. Different subtypes of primary liver cancer can be recognized, of which hepatocellular carcinoma (HCC; originating from hepatocytes) and intrahepatic cholangiocarcinoma (CCA; originating from cholangiocytes) form the largest groups, together constituting over 85% of all primary liver cancers2. Several factors have been linked to increased HCC risk, including chronic alcohol consumption3, as well as metabolic associated fatty liver disease (MAFLD), and its more progressive form nonalcoholic steatohepatitis (NASH), which can be caused by obesity, diabetes, drugs/medication and metabolic conditions4. These factors have also been linked to an increased risk for intrahepatic CCA5. Though less common, chronic inflammation and fibrosis of the biliary tracts, known as primary sclerosing cholangitis (PSC), also confers increased risk of developing both HCC and CCA6. Our knowledge on how these environmental conditions drive liver cancer is still incomplete7.

Chronic alcohol consumption is thought to enhance the mutational load through the metabolite acetaldehyde, which has been reported to be directly mutagenic8 and indirectly through the formation of reactive oxygen species9,10,11,12,13,14. Increased burden of somatic mutations has also been observed in non-alcoholic liver disease15. NASH and PSC are both characterized by chronic inflammation16,17, which may cause the production of reactive oxygen and nitrogen species that subsequently induce DNA damage18.

However, accurate measurements of in vivo induced mutations are required to confirm that accelerated mutagenesis underlies liver tumorigenesis. We previously established a sensitive method to accurately determine all somatic mutations that have been acquired throughout life in individual human adult stem cells of the liver and gastrointestinal tract19. We used these catalogs of somatically acquired mutations to perform mutational signature analysis, a powerful computational method for identifying mutational processes that have been active in the life history of cells20. Currently, 60 single base substitution (SBS) signatures, 11 double base substitution (DBS) signatures, 18 indel signatures, and 16 structural variation (SV) signatures have been identified20,21. These signatures are a result of endogenous mutational sources (such as apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like (APOBEC) activity22 or homologous recombination deficiency23,24), but also microbial impact25, oxidative stress26, or anti-cancer therapies27,28. Mutational signature analyses on healthy human stem cells revealed that mutational processes are tissue-specific and continuously active throughout life resulting in a linear accumulation of mutations with age19, at least under “normal” conditions. Because of the link between liver disease and liver cancer, we hypothesized that the precancerous state of liver diseases would be reflected by increased mutation rates and accumulation of mutational patterns that are characteristic to the type of DNA damage inflicted during liver tumorigenesis.

In this study, we aim to identify the mutational processes that contribute to the precancerous state in common human liver diseases. To achieve this goal, we have studied the accumulation of mutations in individual stem cells derived from livers of patients with alcoholic cirrhosis, NASH and PSC who received a liver transplantation. Surprisingly, we find that individual stem cells from liver patients did not show increased mutational burden overall, within liver cancer associated genes, nor to specific mutational signatures when compared to liver stem cells from healthy donors. Our findings suggest that environmental conditions drive liver tumorigenesis through means other than by increasing mutagenesis.

## Results

### Mutation rates do not increase in diseased livers

Both main liver cell types, hepatocytes and cholangiocytes, can act as liver stem cells depending on the type of tissue damage that was inflicted29. Cholangiocytes show a high degree of cellular plasticity during regeneration and disease and act as facultative liver stem cells during impaired hepatocyte regeneration29,30. Cholangiocytes can be grown as intrahepatic cholangiocyte organoids (ICO)31 that show long-term self-renewal, differentiation, and engraftment in mouse and rat models of liver failure32,33. In contrast, there are no suitable protocols for the long-term culture of human hepatocytes. We have previously exploited the proliferative capacity of individual cholangiocytes to determine mutation rates in individual liver stem cells of the healthy liver19. We reasoned that cholangiocytes are also suitable for the study of somatic mutation accumulation as a result of the diseased liver environment, because cholangiocytes are exposed to the same environmental conditions as the other liver cell types34.

To determine whether a diseased liver environment leads to increased somatic mutation accumulation, we performed whole-genome sequencing (WGS) on clonal ICOs derived from liver stem cells from patients with diseased livers (Fig. 1a, b and Supplementary Data 1). Our study included: (i) 14 clones from five patients with cirrhosis as a result of chronic alcohol consumption; (ii) 13 clones from five NASH patients; and (iii) eight clones from three PSC patients. Organoid establishment success rates were lower for diseased liver as compared to healthy livers and was most challenging for material obtained from PSC patients. For each patient, a reference blood or multilineage bulk tissue sample was also sequenced to distinguish germline variants from somatic variants. The somatic mutation catalogs in diseased livers were compared to the mutation catalogs from 14 clones from liver adult stem cells derived from seven healthy donors.

In total, we identified 172,650 small mutations (single/double base substitutions or indels) in healthy (35,006), post-alcoholic (49,681), NASH (43,997), and PSC (19,521) livers (as well as 24,445 from a patient with HCC). Consistent with previous observations19, somatic mutations accumulated linearly with age in healthy liver cells, at a rate of approximately 46 SBSs and nine indels per year (Fig. 2). The rate of SBS accumulation in alcoholic, NASH, and PSC ICOs showed no significant differences to that of healthy ICOs (Z-test, p ≥ 0.41), and similarly, the rate of indel accumulation in disease ICOs was also comparable to that of healthy ICOs (Z-test, p ≥ 0.62). We observed a slight increase in variance in SBS accumulation in alcoholic cirrhosis versus healthy ICOs which may suggest that alcohol consumption leads to mutagenesis in some but not all patients, though the increase in variance was weak (F-test, p = 0.032). Likewise, we found increased variance in SBS accumulation in PSC ICOs (F-test, p = 0.001), though this variance is likely due to having many more ICOs originating from one patient PSC2. We also compared the accumulation of DBSs and SVs. While the number of DBSs and SVs was too low to be conclusive (≤50 DBSs and ≤25 SVs; Supplementary Fig. 2), the rate of mutation accumulation overall did not increase in diseased versus healthy ICOs. Taken together, these results suggest that chronic alcohol consumption or an inflamed liver environment does not lead to increased SBS, indel, DBS or SV accumulation in liver cells.

### The mutation profile of diseased livers is similar to that of healthy livers

The presence of genome-wide patterns of mutations (also known as mutational signatures) reflects past activity of mutational processes in cells. Previously, the mutational signatures SBS12 and SBS16 have been associated with HCC20,35, with SBS16 also being associated with alcohol consumption36. Additionally, SBS2 and SBS13 (APOBEC activity) were found to be active in numerous cancer types including CCA20. We expected that the mutational processes in diseased liver would be similar to those in cancerous liver. We thus examined whether liver disease results in increased presence of one of the above signatures or in other signatures related to liver cancer. Since the main liver cancer types are HCC and intrahepatic CCA, we selected the signatures that could be present in liver and biliary cancer based on the signature catalog from the PCAWG (Pan-Cancer Analysis of Whole Genomes) consortium20 (see “Methods” section for further details). We ultimately quantified the presence of ten SBS and seven indel signatures in our ICOs as well as in the PCAWG HCC and CCA samples (Fig. 3). Too few DBSs and SVs were present in the diseased liver samples to perform signature analysis for these variant types (Supplementary Fig. 2).

We observed similar signature profiles in diseased and healthy ICOs, with the most predominant signatures being age-related (Fig. 3; SBS1, SBS5, SBS40, ID1, ID5, and ID8), which are present in normal cells20. We also found minor contributions of signatures SBS4 and ID3. These signatures are also present at a baseline level in HCC and CCA20. In the HCC and CCA samples from PCAWG, age-related signatures were predominant, with HCC samples also showing increased contribution of the known HCC associated signatures SBS12 and SBS16 compared to the healthy and diseased ICOs, while CCA samples showed increased contribution of APOBEC activity (SBS2, SBS13). Because these signatures were not observed in the disease ICOs, these findings argue against environmentally induced mutational processes as a force driving the transition of healthy to precancerous liver cells. It is likely that other events are required to initiate the HCC/CCA related mutational processes.

Given that certain recurrent chromosome arm gains and losses have been reported in HCC37,38, we also examined whether liver disease leads to similar copy number alterations (CNA). We find that aside from two samples with polyploidization (a known phenomenon in normal liver cells39), the genomes of the disease ICOs were relatively stable, while the genomes of HCC and CCA were clearly more unstable with known recurrent arm gains (e.g., 1q, 8q, 17q) and losses (e.g., 8p and 17p) being observed37,38 (Supplementary Fig. 3). These data could suggest that CNA accumulation due to liver disease does not contribute to the healthy to precancerous liver transition. However, as CNAs are rare in non-cancerous cells40, more data would be required to validate this hypothesis.

### Absence of driver gene mutations in ICOs

Certain mutations may confer liver stem cells a growth advantage, and in diseased livers, these cells may be able to proliferate more to regenerate lost tissue. We thus examined whether the liver disease conditions resulted in positive selection of cells with non-synonymous mutations in specific genes using the dndscv algorithm (see details in “Methods” section). However, amongst all of the liver disease groups, no genes were found to be enriched in non-synonymous mutations (q < 0.01, Supplementary Data 2). In line with this result, we did not observe any coding, promoter or 5′/3′ untranslated region (UTR) mutations in driver genes of HCC and CCA (obtained from Intogen; see “Methods” section), except for one missense variant in an alcoholic cirrhosis ICO sample (TERT c.1588C>G) (Fig. 4). This could potentially be explained by the cells from which our disease ICOs were derived not being actual cancer precursors but only harboring passenger mutations. We acknowledge that mutations could occur in other non-coding elements but have not examined these as their impact is currently difficult to assess41,42. Furthermore, we also acknowledge that our small sample sizes limit our ability to find enriched driver gene mutations in the diseased liver ICOs. In contrast, we found enrichment of non-synonymous mutations in TP53 in the PCAWG CCA samples, and in TP53 and CTNNB1 as together with 13 other genes in the PCAWG HCC samples (q < 0.01, Supplementary Data 2). While dndscv does not consider non-coding variants, we also observed TERT promoter mutations in the PCAWG HCC samples (“upstream_gene_variant”, Fig. 4). These genes have been reported as known cancer driver genes in the respective cancer types43. No mutations in these driver genes were found in 34% (84/284) of PCAWG HCC samples indicating that mutations in these genes are not necessarily a requirement for HCC development.

## Discussion

Despite the association between liver disease and primary liver cancer (which includes HCC and CCA), the underlying mechanisms of tumorigenesis remain debated. The prevalent view is that tumorigenesis results from an increased mutational burden8,9,10,11,12,13,14,15,16,17,18. In line with this view, Brunner et al.15 showed that cirrhotic liver cells from NASH patients exhibited an increase in mutational load, even though this increase was small and variance between patients was high. In contrast with these observations, we did not observe an increase in mutation rate in individual stem cells of precancerous livers. Additionally, while previous studies have associated specific mutational signatures and gene mutations in HCC to alcohol consumption44,45, we did not observe an altered mutational landscape in our alcoholic cirrhosis ICOs.

There are several possible explanations for the unchanged mutational landscape in our liver disease stem cells. Firstly, it could be that the cholangiocytes that give rise to the ICO cultures are not the precursors to (pre-)cancerous liver cells. Hepatocellular carcinoma is derived from hepatocytes and not from cholangiocytes and the cholangiocytes that give rise to ICOs may not be representative for the cholangiocytes that have the potential to develop into cholangiocarcinoma. Secondly, there could be selection for the most stable ASCs in culture, which may not necessarily be the precancerous cells. Thirdly, it may be possible that increased mutagenesis affected only a small proportion of the liver cells and that we have sequenced too few cells to identify the hypermutated cells. Lastly, it may be that our patient cohorts were too small to detect changes in mutational landscape, which may indeed be the case for detecting enrichment of driver gene mutations. However, we could determine via a power analysis (Supplementary Note 1) that the sample sizes in our study were sufficient to detect changes in mutational load and mutational signature contribution similar to other studies which also used tissue derived organoids to investigate mutation accumulation26,46,47. It may be however possible to pick up more subtle mutational impacts by increasing sample sizes.

On the other hand, the stem cells that we studied have been exposed to the same environmental stressors as the cells that grew out as tumors (since both have been in contact with the blood stream), so when there would have been a direct mutational impact, we would expect this to be detectable in the cells that we studied. In line with our findings, Brunner et al.15 also found that alcoholic cirrhosis and NASH livers did not exhibit different mutational signatures in comparison to healthy livers. While the mutational impact of NASH and PSC have not been investigated besides in the Brunner et al.15 study, alcohol exposure has been shown to lead to DNA damage in vitro10,13,48 and in blood cells in mice49, but these studies do not accurately reflect the in vivo liver environmental conditions. It is possible that the rate of cellular endocytosis and/or diffusion is slower in the liver, resulting in less exposure to alcohol than would occur in vitro. Additionally, the aforementioned studies were performed on non-quiescent cells which may rely more on replicative repair, whereas liver cells are generally quiescent and likely rely on non-replicative repair which is faster at repairing alcohol induced DNA damage than replicative repair8. Alternatively, liver cells that acquire alcohol induced DNA damage may undergo apoptosis and be replaced by new cells as a result of liver regeneration50, and these cells in turn lack the mutation footprint caused by the alcohol. Nevertheless, the absence of increased mutational burden in our disease ICOs may suggest that increased mutagenesis is not the primary contributing factor towards tumorigenesis.

Opposed to the view that tumorigenesis arises from mutagenesis, an alternative hypothesis proposes that chronic liver inflammation and cirrhosis (which commonly precedes primary liver cancer51) leads to cell death in the liver, requiring normally quiescent liver adult stem cells to proliferate at a much higher rate to regenerate the damaged liver. As a consequence, cells would accumulate more mutations, especially those caused by background mutational mechanisms related to cell proliferation (e.g., aging-associated mutational signatures). Inflammatory disease conditions would thus provide a “fertile ground” for cells with random and potentially pre-existing (oncogenic) mutations that confer a selective growth advantage to clonally expand52,53,54. Such a phenomenon has been described in mouse models, whereby pancreatic cells within mice with both a pathogenic Kras mutation and pancreatitis transitioned into an epigenetic state similar to pancreatic ductal adenocarcinoma, while pancreatic cells in mice with only one or the other retained their original epigenetic state55. Additionally, Hepatitis C Virus (HCV)-induced cirrhotic livers showed an increase in the number and size of clonal patches with mutations in genes that are frequently mutated in HCC50.

Taken together, our findings suggest that mechanisms other than direct mutagenesis drive the transition from healthy to precancerous liver, and highlights the need to explore other potential hypotheses of liver tumorigenesis, including but not limited to the ‘fertile ground’ hypothesis.

## Methods

### Human tissue material

All human tissue biopsies were obtained in the Erasmus MC—University Medical Center Rotterdam. Liver biopsies from healthy liver donors and patients with alcoholic cirrhosis, nonalcoholic steatohepatitis (NASH) or primary sclerosing cholangitis (PSC) were obtained during liver transplantation procedures. All patients were negative for viral infection and metabolic diseases. The biopsies were collected in cold organ preservation fluid (Belzer UW Cold Storage Solution, Bridge to Life, London, UK) and transported and stored at 4 °C until use. The acquisition of these liver biopsies for research purposes was approved by the Medical Ethical Committee of the Erasmus Medical Center (MEC-2014-060 and MEC-2013-143). Informed consent was provided by all patients involved.

### Generating clonal intrahepatic cholangiocyte organoid cultures from human liver biopsies

Healthy and diseased liver tissue biopsies were washed in cold DMEM (ThermoFisher) supplemented with 1% fetal calf serum (FCS) and 1% penicillin-streptomycin (wash solution). Subsequently, the tissue was transferred to a petri dish and thoroughly minced with scalpel blades. The minced tissue was transferred to 4 ml digestion solution consisting of EBSS with Ca2+/Mg2+ (ThermoFisher) with 1 mg/ml Collagenase type IA (Sigma, C9891) and 0.1 mg/ml DNAse I (Sigma DN25). The tissue was incubated for 30 min at 37 °C with regular shaking. Next, the suspension was passed through a pipet to further break up the tissue and passed through a 70 µm Nylon cell strainer. The cells were washed once with wash solution, followed by two washes in Advanced DMEM F12 supplemented with 1% penicillin–streptomycin, 10 mM HEPES, and 1× Glutamax (all from ThermoFisher). After the final wash, the cell pellet was resuspended in Matrigel (Corning) and plated in 40 µl droplets per well in prewarmed non-adhesive 24-well plates. The plates were placed at 37 °C in a humidified atmosphere and 5% CO2. After Matrigel had solidified, 500 µl liver organoid establishment medium was added to the wells. Establishment medium consisted of Advanced DMEM F12 supplemented with 1% penicillin–streptomycin, 10 mM HEPES, 1× Glutamax, 10% R-Spondin conditioned medium (produced in house), B27 supplement without Vitamin A (ThermoFisher), N2 supplement (ThermoFisher), 10 mM Nicotinamide (Sigma Aldrich), 1.25 mM N-acetylcysteine (Sigma Aldrich), Primocin, 5 µM A83-01 (Tocris Bioscience), 10 µM Forskolin (Tocris Bioscience), 100 ng/ml FGF-10 (Peprotech), recombinant human Noggin (Peprotech), 10 μM Rho kinase inhibitor (Abmole), hES cell cloning & recovery supplement (Stemgent), 25 ng/ml HGF (Peprotech), 10 nM Gastrin (Tocris), 50 ng/ml human EGF (Peprotech), and 0.3 nM Wnt-surrogate Fc protein (U-protein Express BV). After 2–3 days after isolation, the first intrahepatic cholangiocyte organoids started to appear and establishment medium was switched to maintenance medium consisting of Advanced DMEM F12 supplemented with 1% penicillin-streptomycin, 10 mM HEPES, 1× Glutamax, 10% R-Spondin conditioned medium, B27 supplement without Vitamin A, N2 supplement 10 mM Nicotinamide, 1.25 mM N-acetylcysteine, Primocin, 5 µM A83-01, 10 µM Forskolin, 100 ng/ml FGF-10, 25 ng/ml HGF (Peprotech), 10 nM Gastrin (Tocris), and 50 ng/ml human EGF (Peprotech). The cultures were maintained for 10–14 days after isolation, to enrich for adult stem cells. Subsequently, clonal organoid cultures were generated from these organoid cultures by FACS or by manual selection and expansion of individual organoids56. The organoid cultures were further expanded until there was enough material for DNA isolation. DNA was isolated from all organoid cultures, blood samples, and tissue biopsies using the Qiasymphony (Qiagen). Whole-genome sequencing libraries were generated from 200 ng of genomic DNA according to standard protocols (Illumina). The organoid cultures and control samples were sequenced paired-end (2 × 100 bp) to a depth of at least 30× coverage on the Illumina HiSeq Xten. The hepatocellular carcinoma biopsies were sequenced paired-end (2 × 100 bp) to a depth of at least 60× coverage on the Illumina HiSeq Xten. Whole-genome sequencing was performed at the Hartwig Medical Foundation in Amsterdam, the Netherlands.

### Variant calling

Germline and somatic variant calling for all samples was performed using the HMF pipeline (https://github.com/hartwigmedical/pipeline; v4.8)57. Briefly, reads were mapped to GRCh37 using BWA-MEM v0.7.5a with duplicates being marked for filtering. Indels were realigned using GATK v3.4.46 IndelRealigner. GATK Haplotype Caller v3.4.46 was used for calling germline variants in the reference sample. For somatic SNV and indel variant calling, GATK BQSR3 was first used to recalibrate base qualities, followed by Strelka v1.0.14 for the variant calling itself. Somatic structural variant calling was performed using GRIDSS v1.8.0. Copy-number calling was performed using PURity & PLoidy Estimator (PURPLE), that combines B-allele frequency (BAF), read depth, and structural variants to estimate the purity and copy number profile of a tumor sample. For SNVs and indels, downstream analyses were performed only on variants marked as “PASS”.

### Mutation context analysis

The counts of single base substitution (SBS), double base substitution (DBS), indel and structural (SV) variant contexts were determined from somatic VCF files using the R package mutSigExtractor (https://github.com/UMCUGenetics/mutSigExtractor; v1.23). The mutation contexts of all mutation types are described in COSMIC (https://cancer.sanger.ac.uk/cosmic/signatures), except for the SV contexts. The 16 SV contexts were composed of the SV type (deletion, duplication, inversion, translocation) and the SV length (1–10 kb, 10–100 kb, 100 kb–1 Mb, 1–10 Mb, and >10 Mb). Note that SV length is not applicable for translocations. The mutation context spectra for each sample group are shown in Supplementary Figs. 47.

To perform mutational signature analysis, we selected the SBS and DBS signatures that were present in at least 10% of liver cancer (Liver-HCC) or biliary cancer (Biliary-AdenoCA) PCAWG samples (https://dcc.icgc.org/releases/PCAWG/mutational_signatures/Signatures_in_Samples/SP_Signatures_in_Samples)20. We then fitted the SBS and DBS mutation contexts to these selected signatures using the fitToSignatures() function from mutSigExtractor (which employs the non-negative least-squares method) to obtain absolute signature contributions. Relative signature contribution per sample was calculated by dividing the absolute contributions by the total signature contribution.

Mutation context and signature absolute contributions per sample can be found in Supplementary Data 3.

### Selection of liver and biliary cancer driver genes

A catalog of driver genes by cancer type was downloaded from Intogen (https://www.intogen.org/download; release 2020.02.01). From the Compendium_Cancer_Genes.tsv file, we selected genes where CANCER_TYPE was “HC” or “CH” (hepatocellular carcinoma and cholangiocarcinoma, respectively), and CGC_CANCER_GENE was “TRUE”. Additionally, TERT has been reported by the PCAWG consortium43 as a known HCC driver gene and was thus also included.

### Identifying non-synonymous mutations

The mutation type of each somatic SNV/indel was determined by SnpEff (http://snpeff.sourceforge.net/; v4.3t). The following variant types were considered non-synonymous mutations: out-of-frame frameshifts, disruptive inframe frameshifts, nonsense, missense, splice variants. We also considered a mutation as non-synonymous if it was annotated as VUS, likely pathogenic, or pathogenic by ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/; GRCh37, database date 2020-02-24), or if it was a hotspot mutation. The underlying code for annotating non-synonymous mutations can be found at https://github.com/UMCUGenetics/geneDriverAnnotator (v1.0).

The dndscv R package52 was used to identify genes that were enriched for non-synonymous mutations. Briefly, this package computes the (local) background mutation rates and sequence composition of genes to calculate the background mutation rate for each gene. A likelihood ratio test is subsequently performed to identify genes that are significantly hit by nonsynonymous mutations. dndscv was run separately for each disease status group (i.e., separately for healthy ICOs, separately for NASH ICOs, etc) using all the somatic mutations from the respective group.

### Statistics and reproducibility

All statistical analyses were performed in R (v4.0.3). To correlate the number of mutations with the age of each patient from which each biopsy was derived (as shown in Fig. 2 and Supplementary Fig. 1), we first assessed normality of the mutational load per mutation type per disease status group was using the Shapiro test (shapiro.test() function) (Supplementary Table 1). This confirmed that the mutational load was normally distributed (p > 0.05), with near normality for SBS load in PSC samples (p = 0.03) and indel load in NASH samples (p = 0.05). Then, the lme() function from the nlme (v3.1) package was used to fit a linear mixed effects regression, with 95% confidence intervals being calculated using the intervals() function from the nlme package. Here, “patient” was modelled as a random effect to account for having different numbers of organoids per patient. Additionally, the intercept was fixed to zero as it was assumed that a patient has no somatic mutations at birth. A two-sided Z-test was used to calculate the difference between two regressions. The Z-statistic was first calculated using the slope (m) and standard errors (SE) of the two regressions (Eq. 1), which was then used to calculate a p-value using the pnorm() function (Eq. 2). A one-sided F-test was performed to calculate whether the variance of the regression of diseased ICOs was greater than that of the healthy ICOs. The F-statistic was first calculated by dividing the variance of the two regressions (extracted from the output of the lme() function) (Eq. 3), which was then used to calculate a p-value using the pf() function (Eq. 4).

$$Z=\frac{{m}_{1}-{m}_{2}}{\sqrt{{{{{{\mathrm{SE}}}}}}_{1}^{\,2}-{{{{{\mathrm{SE}}}}}}_{2}^{\,2}}}$$
(1)
$$p=2\times {{{{{\mathrm{pnorm}}}}}}(-|Z|)$$
(2)
$$F={{{{{\mathrm{var}}}}}}_{{{{{{\mathrm{disease}}}}}}}/{{{{{\mathrm{var}}}}}}_{{{{{{\mathrm{healthy}}}}}}}$$
(3)
$$p=1-pf(F)$$
(4)

To determine whether there was a significant increase in mutation context load in disease versus healthy ICOs (as shown in Supplementary Figs. 47), Wilcoxon rank sum tests (using the wilcox.test() function) were performed per mutation context. Bonferroni multiple testing correction was then applied to the resulting p-values (using the p.adjust() function).

### Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

The BAM files from the whole-genome sequencing data generated in the current study are available at EGA (https://www.ebi.ac.uk/ega/home) under accession numbers EGAS00001002983 and EGAS00001005384. BAM files from hepatocellular carcinoma and cholangiocarcinoma patients from the Pan-Cancer Analysis Whole Genomes (PCAWG) consortium were obtained under request number DACO-5333. For access to the PCAWG BAM files, researchers will need to request access via the ICGC Data Access Compliance Office (DACO; https://daco.icgc.org/). The VCF and tabular files produced from somatic variant calling are available at https://zenodo.org/record/556238158.

## Code availability

The code for the Hartwig Medical Foundation (HMF) germline and somatic variant calling pipeline is available at https://github.com/hartwigmedical/pipeline. The code used for data processing and generating the figures is available at https://github.com/UMCUGenetics/Diseased_livers59.

## References

1. 1.

Jemal, A. et al. Global cancer statistics. Cancer J. Clin. 61, 69–90 (2011).

2. 2.

Altekruse, S. F., Devesa, S. S., Dickie, L. A., McGlynn, K. A. & Kleiner, D. E. Histological classification of liver and intrahepatic bile duct cancers in SEER registries. J. Registry Manage. 38, 201–205 (2011).

3. 3.

Testino, G., Leone, S. & Borro, P. Alcohol and hepatocellular carcinoma: a review and a point of view. World J. Gastroenterol. 20, 15943–15954 (2014).

4. 4.

Loomba, R., Friedman, S. L. & Shulman, G. I. Mechanisms and disease consequences of nonalcoholic fatty liver disease. Cell 184, 2537–2564 (2021).

5. 5.

Wongjarupong, N. et al. Non-alcoholic fatty liver disease as a risk factor for cholangiocarcinoma: a systematic review and meta-analysis. BMC Gastroenterol. 17, 149 (2017).

6. 6.

Razumilava, N., Gores, G. J. & Lindor, K. D. Cancer surveillance in patients with primary sclerosing cholangitis. Hepatology 54, 1842–1852 (2011).

7. 7.

Duan, X.-Y., Zhang, L., Fan, J.-G. & Qiao, L. NAFLD leads to liver cancer: do we have sufficient evidence? Cancer Lett. 345, 230–234 (2014).

8. 8.

Hodskinson, M. R. et al. Alcohol-derived DNA crosslinks are repaired by two distinct mechanisms. Nature 579, 603–608 (2020).

9. 9.

Tamura, M., Ito, H., Matsui, H. & Hyodo, I. Acetaldehyde is an oxidative stressor for gastric epithelial cells. J. Clin. Biochem. Nutr. 55, 26–31 (2014).

10. 10.

Novitskiy, G., Traore, K., Wang, L., Trush, M. A. & Mezey, E. Effects of ethanol and acetaldehyde on reactive oxygen species production in rat hepatic stellate cells. Alcohol. Clin. Exp. Res. 30, 1429–1435 (2006).

11. 11.

van Loon, B., Markkanen, E. & Hübscher, U. Oxygen as a friend and enemy: How to combat the mutational potential of 8-oxo-guanine. DNA Repair 9, 604–616 (2010).

12. 12.

Obe, G. & Ristow, H. Mutagenic, cancerogenic and teratogenic effects of alcohol. Mutat. Res. 65, 229–259 (1979).

13. 13.

Helander, A. & Lindahl-Kiessling, K. Increased frequency of acetaldehyde-induced sister-chromatic exchanges in human lymphocytes treated with an aldehyde dehydrogenase inhibitor. Mutat. Res. Lett. 264, 103–107 (1991).

14. 14.

Matsuda, T., Kawanishi, M., Matsui, S., Yagi, T. & Takebe, H. Specific tandem GG to TT base substitutions induced by acetaldehyde are due to intra-strand crosslinks between adjacent guanine bases. Nucleic Acids Res. 26, 1769–1774 (1998).

15. 15.

Brunner, S. F. et al. Somatic mutations and clonal dynamics in healthy and cirrhotic human liver. Nature 574, 538–542 (2019).

16. 16.

Michelotti, G. A., Machado, M. V. & Diehl, A. M. NAFLD, NASH and liver cancer. Nat. Rev. Gastroenterol. Hepatol. 10, 656–665 (2013).

17. 17.

Dyson, J. K., Beuers, U., Jones, D. E. J., Lohse, A. W. & Hudson, M. Primary sclerosing cholangitis. Lancet 391, 2547–2559 (2018).

18. 18.

Kawanishi, S., Ohnishi, S., Ma, N., Hiraku, Y. & Murata, M. Crosstalk between DNA damage and inflammation in the multiple steps of carcinogenesis. Int. J. Mol. Sci. 18, 1808 (2017).

19. 19.

Blokzijl, F. et al. Tissue-specific mutation accumulation in human adult stem cells during life. Nature 538, 260–264 (2016).

20. 20.

Alexandrov, L. B. et al. The repertoire of mutational signatures in human cancer. Nature 578, 94–101 (2020).

21. 21.

Li, Y. et al. Patterns of somatic structural variation in human cancer genomes. Nature 578, 112–121 (2020).

22. 22.

Chan, K. et al. An APOBEC3A hypermutation signature is distinguishable from the signature of background mutagenesis by APOBEC3B in human cancers. Nat. Genet. 47, 1067–1072 (2015).

23. 23.

Nguyen, L., Martens, J. W. M., Van Hoeck, A. & Cuppen, E. Pan-cancer landscape of homologous recombination deficiency. Nat. Commun. 11, 5584 (2020).

24. 24.

Nik-Zainal, S. et al. Landscape of somatic mutations in 560 breast cancer whole-genome sequences. Nature 534, 47–54 (2016).

25. 25.

Pleguezuelos-Manzano, C. et al. Mutational signature in colorectal cancer caused by genotoxic pks E. coli. Nature 580, 269–273 (2020).

26. 26.

Kuijk, E. et al. The mutational impact of culturing human pluripotent and adult stem cells. Nat. Commun. 11, 2493 (2020).

27. 27.

Kucab, J. E. et al. A compendium of mutational signatures of environmental agents. Cell 177, 821–836.e16 (2019).

28. 28.

Christensen, S. et al. 5-Fluorouracil treatment induces characteristic T>G mutations in human cancer. Nat. Commun. 10, 4571 (2019).

29. 29.

Post, Y. & Clevers, H. Defining adult stem cell function at its simplest: the ability to replace lost cells through mitosis. Cell Stem Cell 25, 174–183 (2019).

30. 30.

Raven, A. et al. Cholangiocytes act as facultative liver stem cells during impaired hepatocyte regeneration. Nature 547, 350–354 (2017).

31. 31.

Marsee, A. et al. Building consensus on definition and nomenclature of hepatic, pancreatic, and biliary organoids. Cell Stem Cell 28, 816–832 (2021).

32. 32.

Huch, M. et al. Long-term culture of genome-stable bipotent stem cells from adult human liver. Cell 160, 299–312 (2015).

33. 33.

Kuijk, E. W. et al. Generation and characterization of rat liver stem cell lines and their engraftment in a rat model of liver failure. Sci. Rep. 6, 22154 (2016).

34. 34.

Sia, D., Villanueva, A., Friedman, S. L. & Llovet, J. M. Liver cancer cell of origin, molecular class, and effects on patient prognosis. Gastroenterology 152, 745–761 (2017).

35. 35.

Alexandrov, L. B. et al. Signatures of mutational processes in human cancer. Nature 500, 415–421 (2013).

36. 36.

Li, X. C. et al. A mutational signature associated with alcohol consumption and prognostically significantly mutated driver genes in esophageal squamous cell carcinoma. Ann. Oncol. 29, 938–944 (2018).

37. 37.

Kim, T.-M. et al. Clinical implication of recurrent copy number alterations in hepatocellular carcinoma and putative oncogenes in recurrent gains on 1q. Int. J. Cancer 123, 2808–2815 (2008).

38. 38.

Court, C. M. et al. Somatic copy number profiling from hepatocellular carcinoma circulating tumor cells. NPJ Precis. Oncol. 4, 16 (2020).

39. 39.

Donne, R., Saroul-Aïnama, M., Cordier, P., Celton-Morizur, S. & Desdouets, C. Polyploidy in liver development, homeostasis and disease. Nat. Rev. Gastroenterol. Hepatol. 17, 391–405 (2020).

40. 40.

Beroukhim, R. et al. The landscape of somatic copy-number alteration across human cancers. Nature 463, 899–905 (2010).

41. 41.

Gasperini, M., Tome, J. M. & Shendure, J. Towards a comprehensive catalogue of validated and target-linked human enhancers. Nat. Rev. Genet. 21, 292–310 (2020).

42. 42.

Doni Jayavelu, N., Jajodia, A., Mishra, A. & Hawkins, R. D. Candidate silencer elements for the human and mouse genomes. Nat. Commun. 11, 1061 (2020).

43. 43.

ICGC/TCGA Pan-Cancer Analysis of Whole Genomes Consortium. Pan-cancer analysis of whole genomes. Nature 578, 82–93 (2020).

44. 44.

Letouzé, E. et al. Mutational signatures reveal the dynamic interplay of risk factors and cellular processes during liver tumorigenesis. Nat. Commun. 8, 1315 (2017).

45. 45.

Fujimoto, A. et al. Whole-genome mutational landscape and characterization of noncoding and structural mutations in liver cancer. Nat. Genet. 48, 500–509 (2016).

46. 46.

Drost, J. et al. Use of CRISPR-modified human stem cell organoids to study the origin of mutational signatures in cancer. Science 358, 234–238 (2017).

47. 47.

Jager, M. et al. Deficiency of nucleotide excision repair is associated with mutational signature observed in cancer. Genome Res. 29, 1067–1077 (2019).

48. 48.

Voordeckers, K. et al. Ethanol exposure increases mutation rate through error-prone polymerases. Nat. Commun. 11, 3664 (2020).

49. 49.

Garaycoechea, J. I. et al. Alcohol and endogenous aldehydes damage chromosomes and mutate stem cells. Nature 553, 171–177 (2018).

50. 50.

Zhu, M. et al. Somatic mutations increase hepatic clonal fitness and regeneration in chronic liver disease. Cell https://doi.org/10.1016/j.cell.2019.03.026 (2019).

51. 51.

Massarweh, N. N. & El-Serag, H. B. Epidemiology of hepatocellular carcinoma and intrahepatic cholangiocarcinoma. Cancer Control 24, 1073274817729245 (2017).

52. 52.

Martincorena, I. et al. Universal patterns of selection in cancer and somatic tissues. Cell 171, 1029–1041.e21 (2017).

53. 53.

Hernandez–Gea, V., Toffanin, S., Friedman, S. L. & Llovet, J. M. Role of the microenvironment in the pathogenesis and treatment of hepatocellular carcinoma. Gastroenterology 144, 512–527 (2013).

54. 54.

Zhu, L. et al. Multi-organ mapping of cancer risk. Cell 166, 1132–1146.e7 (2016).

55. 55.

Alonso-Curbelo, D. et al. A gene-environment-induced epigenetic program initiates tumorigenesis. Nature 590, 642–648 (2021).

56. 56.

Jager, M. et al. Measuring mutation accumulation in single human adult stem cells by whole-genome sequencing of organoid cultures. Nat. Protoc. 13, 59–78 (2018).

57. 57.

Priestley, P. et al. Pan-cancer whole-genome analyses of metastatic solid tumours. Nature 575, 210–216 (2019).

58. 58.

Nguyen, L. Precancerous liver diseases do not cause increased mutagenesis in liver stem cells. https://doi.org/10.5281/zenodo.5562381. (2021).

59. 59.

Nguyen, L. UMCUGenetics/Diseased_livers: natcomm_submission. https://doi.org/10.5281/zenodo.5564971. (2021).

## Acknowledgements

This study was financially supported by the research program InnoSysTox (project number 114027003), by the Netherlands Organisation for Health Research and Development (ZonMw), by the Dutch Cancer Society (project number 10496) and is part of the Oncode Institute, which is partly financed by the Dutch Cancer Society and was funded by the gravitation program CancerGenomiCs.nl from the Netherlands Organisation for Scientific Research (NWO). The authors would like to thank Hartwig Medical Foundation and the Utrecht Sequencing Facility for performing the whole genome sequencing. The Utrecht Sequencing Facility is subsidized by the University Medical Center Utrecht, Hubrecht Institute, and Utrecht University.

## Author information

Authors

### Contributions

R.L., J.J., J.I., M.D., and M.V. collected liver biopsies. M.J., E.K., and N.B. performed organoid culturing. L.N., M.J., M.L., B.R., R.J., and S.B. performed bioinformatic analyses. M.J., E.K., M.V., R.B., L.L., and E.C. were involved in the conceptual design of this study. L.N., E.K., M.J., and E.C. wrote the manuscript. All authors provided textual comments and have approved the manuscript. E.K., R.B., L.L., and E.C. supervised this study.

### Corresponding authors

Correspondence to Edwin Cuppen or Ewart Kuijk.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

Peer review information Communications Biology thanks Yun-Wen Zheng and the other, anonymous, reviewers for their contribution to the peer review of this work. Primary Handling Editor: Eve Rogers. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Nguyen, L., Jager, M., Lieshout, R. et al. Precancerous liver diseases do not cause increased mutagenesis in liver stem cells. Commun Biol 4, 1301 (2021). https://doi.org/10.1038/s42003-021-02839-y

• Accepted:

• Published:

• DOI: https://doi.org/10.1038/s42003-021-02839-y