|Home | About | Journals | Submit | Contact Us | Français|
This study aims to determine the impact of intracellular hepatitis B virus (HBV) DNA, covalently closed circular DNA (cccDNA) and viral replicative activity in both tumour and non-neoplastic liver on prognosis and to determine the relationship of viral replicative activity and Ishak fibrosis in predicting outcome following resection.
A total of 99 prospectively enrolled patients treated with primary liver resection for HBV-HCC are included. Intracellular HBV DNA and cccDNA were quantitated by real-time PCR. The RNA-sequencing (RNA-seq) was performed in a subset of 21 patients who had either minimal liver fibrosis (Ishak stages 0–2) or end-stage fibrosis (Ishak stage 6).
Tumour tissue contained a lower cccDNA copy number compared with paired non-neoplastic liver, and larger tumours (>3cm) had less cccDNA compared with small tumours (3cm). High viral replicative activity in non-neoplastic liver was associated with higher HCC recurrence rate independent of Ishak fibrosis stage. Genes correlated with viral replicative activity in non-neoplastic liver (620 genes) were distinct from those associated with end-stage fibrosis (1226 genes). Genes associated with viral replicative activity were preferentially distributed in regions on chr3, chr16 and chr19.
Viral replicative activity in non-neoplastic liver is associated with HCC recurrence through mechanisms that are distinct from and independent of Ishak fibrosis stage.
Hepatitis B virus (HBV) is an enveloped DNA virus that primarily infects hepatocytes. In a mature virion, viral DNA is present in a relaxed circular form (rcDNA) that is partially double stranded. After entering the host cell, viral genome is transported into the nucleus, where rcDNA is converted into covalently closed circular DNA (cccDNA) (Levrero et al, 2009). The content of cccDNA, once formed, is stable in quiescent hepatocytes, and can be reduced by dilution through cell division, killing of infected cells or suppression by antiviral cytokines (Belloni et al, 2012).
Persistent HBV infection is a major risk factor for development of HBV-associated hepatocellular carcinoma (HCC). Viral proteins, viral genome integration into the host genome, inflammation, oxidative stress and prolonged fibrotic responses associated with chronic HBV infection may all play important roles in hepatocarcinogenesis (Xu et al, 2010; Zhou et al, 2010; Arzumanyan et al, 2013). Our previous studies have shown that end-stage fibrosis (Ishak stage 6) is independently associated with poor overall survival and higher HCC recurrence rate (Wang et al, 2013b, 2013c). Moreover, our previous study using paraffin-embedded non-neoplastic liver tissue specimens from 111 HBV-HCC patients treated between 1991 and 2008 showed that viral replicative activity, as indicated by the proportion of cccDNA in total HBV DNA, correlated with necroinflammatory damage and Ishak fibrosis stage and that high viral replicative activity in the non-neoplastic liver was associated with poor outcome for HBV-HCC patients (Wang et al, 2013c).
The goal of this study is to determine whether the impact of HBV viral replicative activity on outcome is mediated through mechanisms independent of or secondary to fibrosis. This report describes results from a prospective cohort of patients, completely distinct from the previous retrospective cohort. Only data generated from fresh tissue specimens, obtained from prospectively enrolled primary liver resection patients between 2008 and 2013, are included in this study. Intracellular HBV DNA and cccDNA and viral replicative activity (cccDNA/HBV DNA ratio) were measured using standardised qualitative PCR (Laras et al, 2006; Volz et al, 2007; Pollicino et al, 2011) and compared in 99 paired tumour and non-neoplastic liver tissue specimens of HBV-HCC patients. The impact of intracellular viral DNA and viral replication on clinical outcome following surgical resection with intent to cure is also determined.
This is a prospective analysis of 99 HBV-HCC patients who underwent primary surgical resection with curative intent at the Mount Sinai Medical Center in New York between 2008 and 2013. Of these patients, 90% had documented therapy with antiviral drugs, and most had been treated for a short course of time preceding their HCC diagnosis. Majority of the patients received one type of antiviral (82%), 17% received two types of antiviral and 1% received three types of antiviral treatment. Most commonly used antivirals were entecavir (39%), tenofovir (34%) and adefovir (23%). The other antivirals such as telbivudine made up 4% of all antiviral treatment types. Preoperative serum HBV DNA was determined from the medical record. In our institutional clinical lab serum HBV DNA levels are quantitated using the COBAS AmpliPrep-COBAS TaqMan HBV test with a reported range of 23−2 × 108IUml−1 (Faria et al, 2008). Median follow-up of the survivors is 30 months.
Liver resection specimens were assessed by a single liver pathologist. The degree of hepatic fibrosis in non-neoplastic liver tissue was scored using the modified Ishak method (stages 0–6) as follows: stage 0 (normal liver); stage 1 (fibrosis expansion of a few portal tracts); stage 2 (fibrosis of all portal tracts); stage 3 (fibrous expansion of most portal areas with occasional portal to portal bridging); stage 4 (fibrous expansion of portal areas with marked bridging); stage 5 (marked bridging with occasional nodules); and stage 6 (established cirrhosis with the tissue entirely composed of nodules) (Goodman, 2007; Wang et al, 2013b). Patients with multiple tumours were not excluded from this study. Histologic activity index (HAI, scored from 0 to 18) representing necroinflammatory grade, tumour differentiation and vascular invasion were assessed by the pathologist.
Detailed methodology of patient diagnosis and follow-up is included in Supplementary Materials and Methods.
The HBV DNA and cccDNA were amplified from genomic DNA using the QIAamp DNA extraction kit as previously described (Qiagen, Valencia, CA, USA) (Kuang et al, 2010; Schrader and Iredale, 2011; Wang et al, 2013c). Detailed methodology is included in Supplementary Materials and Methods.
Gene expression was evaluated in a subset of 21 patients using RNA-seq technology in non-neoplastic liver (Supplementary Table 3). All patients with tumour size smaller than 5cm, either with minimal liver fibrosis (Ishak stages 0–2, n=10) or with probable or definite cirrhosis (Ishak stage 6, n=11), were included. Total RNA (1–3μg per sample) was submitted to the Mount Sinai Genomics Core Facility for quality control analysis. Detailed methodology for RNA-seq is included in Supplementary Materials and Methods.
The expression level of genes was estimated by RPKM (reads per kilobase of exon model per million mapped reads) that is normalised to the length of a gene and total mapped reads in one sample. Genes that correlated with viral replicative activity were identified using Spearman's correlation test. At P<0.01, 620 genes correlated with viral replicative activity, corresponding to the false discovery rate (FDR) of 0.2 based on the BH method (Wang et al, 2011). The functional enrichment of these 620 genes was then analysed by searching the Gene Ontology (GO) database and by comparing with gene sets in the Molecular Signatures Database (Subramanian et al, 2005). The distribution of these genes on chromosomes was determined by dividing the genome into 30Mb bins and counting the number of genes in each bin. The enrichment level was evaluated using Fisher's exact test by comparing the distribution of signature genes and overall genes.
Genes that were differentially expressed in Ishak stages 0–2 vs Ishak stage 6 were similarly identified and compared with those that correlated with viral replicative activity.
All statistical analyses were performed using SPSS 19 (Intel, Santa Clara, CA, USA) and SAS 9.2 (SAS Institute, Cary, NC, USA) software. Detailed methodology is included in Supplementary Materials and Methods.
Institutional review board (IRB) approval of all aspects of the research was obtained before study initiation. Written informed consent was obtained for every patient enrolled in this study. The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki, as reflected by the appropriate institutional review committee.
This study included 99 HBV-HCC patients who were prospectively enrolled (Supplementary Table 1). Mean age was 55 years, and the majority of patients were male (83%) and of Asian race (86%). Mean tumour size was 4.2cm, and median tumour size was 3.0cm. Mean HAI was 3.8, and mean Ishak stage was 4.2.
Intracellular HBV DNA and cccDNA and viral replicative activity as indicated by the cccDNA/HBV DNA ratio were compared in tumour and paired non-neoplastic liver specimens (n=99, Table 1). The HBV DNA correlated with cccDNA copy number in both tumour and non-neoplastic liver (Spearman's correlation coefficient: r=0.63 and r=0.56, respectively, P<0.001). Compared with paired non-neoplastic liver, tumour tissue was more likely to have undetectable HBV DNA. Tumour tissue had lower cccDNA copy number than paired non-neoplastic liver tissue, and was more likely to have undetectable cccDNA (38% in tumour vs 6% in non-neoplastic liver). In the samples that had detectable cccDNA, viral replicative activity in tumour correlated with that in non-neoplastic liver (Spearman's correlation r=0.53, P<0.001).
The PCR-based assay to quantitate HBV viral load amplifies all sequences flanked by the specific primers that include integrated viral sequence and viral genome. A subset of samples were treated with plasmid-safe DNase to assess the impact of viral integration. After DNase treatment, albumin copy number decreased >100-fold, from 6607 × 103 copies to 46 × 103 copies, indicating the effectiveness of DNase treatment. The HBV DNA copy number in the same set of sample decreased from 22 × 103 copies to 3 × 103 copies, reflecting a seven-fold decrease. The copy number of cccDNA decreased from 1 × 103 copies to 0.3 × 103 copies, reflecting a three-fold decrease. Therefore, the PCR-based assay amplified both viral genome and integrated viral DNA.
Early tumours (smaller than median, 3cm) had higher cccDNA copy number per cell within the tumour when compared with larger tumours (>3cm) (median: 0.0004 copy per cell, n=57, vs 0 copy per cell, n=42) (Figure 1A). In contrast, the cccDNA copy number in the non-neoplastic liver of patients with small vs large tumours did not differ. In patients with small tumours, the cccDNA copy number did not differ between tumour and paired non-neoplastic liver, whereas in patients with larger tumours, the cccDNA copy in tumour was significantly lower than the paired non-neoplastic liver. Similar findings were observed for intrahepatic HBV DNA (Figure 1B). Small tumour and large tumours did not differ in viral replicative activity (median cccDNA/HBV DNA ratio: 0.09 in small tumour vs 0.06 in larger tumours).
Copy number of cccDNA within tumour was not different when compared according to tumour differentiation, vascular invasion or α-fetoprotein production (data not shown).
Of the patients in this study, 90% had documented therapy with antiviral drugs. Despite antiviral therapy, 49% of patients had detectable serum HBV DNA at the time of their liver resection, and 96% had detectable HBV DNA in liver tissue. The HBV DNA, cccDNA and cccDNA/HBV DNA ratio in tumour or non-neoplastic liver did not differ between patients with undetectable vs those with detectable serum HBV DNA (Supplementary Table 2).
The HBeAg status immediately before surgery was available for 67 patients (76% were HBeAg negative, whereas 24% were HBeAg positive) (Supplementary Table 2). The cccDNA/HBV DNA ratio in non-neoplastic liver was significantly higher in HBeAg-negative patients, indicating lower replicative activity in the liver tissue, whereas all other viral features were not different. Moreover, the proportion of HBeAg-negative patients was significantly higher in patients with low HBV replicative activity in non-neoplastic liver compared with those with high HBV replicative activity (93% vs 62%, P=0.004) (Supplementary Table 3). Therefore, only HBV replicative activity in non-neoplastic liver correlated with serum HBeAg status.
The 5-year overall survival rate of the 99 patients was 72%, and HCC recurrence rate was 44% by Kaplan–Meier analysis. Most recurrence occurred within 24 months of surgical resection. The impact of total intracellular HBV DNA, cccDNA and replicative activity in both tumour and non-neoplastic liver tissue on outcome was examined by comparing groups according to median cutoff values (Table 2 and Figure 2A and B). None of the viral characteristics in tumour or non-neoplastic liver was associated with overall survival.
Viral replicative activity within the non-neoplastic liver was associated with HCC recurrence rate (Table 2 and Figure 2A and B). High viral replicative activity (lower cccDNA/HBV DNA ratio) in the non-neoplastic liver was associated with higher recurrence rate (Figure 2B). High viral replicative activity (lower cccDNA/HBV DNA ratio) in the non-neoplastic liver was also associated with greater HAI (4.1 vs 3.4, P=0.04) and higher fibrosis stage (4.7 vs 3.9, P=0.01), consistent with our previous studies in the retrospective cohort (Supplementary Table 3). When adding fibrosis stage as a covariant, viral replicative activity in non-neoplastic liver was still significantly associated with HCC recurrence by Cox proportional hazards regression (HR: 2.1 (1.0–4.5), P=0.04), suggesting that viral replicative activity contributes to HCC recurrence independent of fibrosis.
We hypothesise that viral replicative activity and fibrosis in non-neoplastic liver may influence different biological pathways implicated in HCC recurrence. The gene signatures associated with viral replicative activity or liver fibrosis in non-neoplastic liver were then analysed in a subgroup of 21 patients using RNA-seq (Supplementary Table 4). Of these patients, 10 had minimal fibrosis (Ishak stages 0–2), whereas 11 had end-stage fibrosis (Ishak stage 6). All 21 patients had HCC <5cm.
There were 620 of genes whose expression levels correlated with viral replicative activity (Supplementary Table 5), noted as viral replicative activity-associated (VRAA) genes. The VRAA genes were not enriched for any GO terms and only enriched for the signature ‘GINESTIER_BREAST_CANCER_ZNF217_AMPLIFIED_DN' in MSigDB (3.6-fold enrichment, P=2.9 × 10−11), a set of downregulated genes when ZNF217 locus at 20q13 is amplified in breast cancer (Ginestier et al, 2006). The ZNF217 is a zinc-finger protein that is often amplified at 20q13 in breast cancer (Krig et al, 2010). Zinc-finger proteins as a whole were enriched in the VRAA gene set (1.6-fold enrichment, P=0.01).
There were 1226 genes that were differentially expressed in non-neoplastic liver of Ishak fibrosis stages 0–2 vs Ishak stage 6 (Supplementary Table 6). These genes were functionally enriched for three signatures related to liver cancer (Supplementary Table 7). Genes associated with Ishak fibrosis stage showed little overlap with genes in the VRAA gene set, and only 26 common genes were identified (Supplementary Table 8).
Whether these differentially expressed genes are enriched in any chromosomal locations was also determined. The VRAA genes were enriched in regions on chr19, chr16 and chr3 (P<0.01) (Figure 3A). Fibrosis-associated genes were enriched in regions on chr19, chr17 and chr6 (P<0.01) (Figure 3B).
Viral genotype and variants were also determined from RNA-seq analysis in the 21 pairs of tumour and non-neoplastic liver tissue specimens (Supplementary Table 4). The paired tumour and non-neoplastic liver tissue specimens had the same genotype and the same subvariant. Genotypes C and B were most prevalent, 56% and 39%, respectively, and only 5% had genotype A. Most patients with end-stage fibrosis (78%) had viral genotype C, whereas most patients with early-stage fibrosis had viral genotype B (56%).
The assembled HBV viral contigs in each sample were aligned to each other to determine whether there are multiple viral variants in each sample. All samples contained contigs with with less than 1 single-nucleotide polymorphism (SNP) per kb HBV viral sequence, with the exception of contigs in one sample where two HBV viral contigs were different at 17 SNPs per kb (1.7% sequence difference). These results suggest that all samples contain a single HBV variant (sequence difference <4%). Among all assembled viral genomes, only one pair of samples contained HBV genomes that were mapped to drug-resistant strains (https://hbvdb.ibcp.fr/HBVdb/HBVdbResistance).
Although viral replicative activity in non-neoplastic liver correlated with the severity of liver fibrosis, RNA-seq analysis demonstrated that very distinct gene signatures were associated with viral replicative activity or liver cirrhosis. These genes were enriched for different biological functions and distributed on different regions along the human chromosome. These results therefore support the hypothesis that viral replicative activity and fibrosis in non-neoplastic liver may influence different biological pathways implicated in HCC recurrence.
Viral replicative activity correlates with the expression of genes that showed a 3.6-fold enrichment for a set of downregulated genes when ZNF217 locus at 20q13 is amplified in breast cancer ‘GINESTIER_BREAST_CANCER_ZNF217_AMPLIFIED_DN'. The ZNF217 is an oncogene and a transcription factor that is part of the corepressor of REST complex (Cowger et al, 2007). Binding of ZNF217 to the promoter region of one of its targets, p15ink4b (a cyclin-dependent kinase inhibitor and a tumour suppressor gene), leads to aberrant silencing of p15ink4b and cell immortalisation (Thillainadesan et al, 2012). Indeed, ZNF217 is found at the core of the 20q13 amplicon that is frequently amplified in breast cancer, ovarian cancer and prostate carcinoma (Rahman et al, 2012; Szczyrba et al, 2013; Bai et al, 2014). Overexpression of ZNF217 is associated with poor prognosis in breast cancer and ovarian cancer (Rahman et al, 2012; Vendrell et al, 2012). The functional importance of ZNF217 is poorly understood in HCC, although chromosomal aberrations including DNA copy number gains of 20q are often observed in HCC (Patil et al, 2005). The observations that VRAA genes are enriched for copy number variation (CNV)-regulated genes but not for any GO biological processes suggest that VRAA genes may not result from biological perturbations, but rather from variations in genomic structure.
Tumour, when compared with non-neoplastic liver, had lower cccDNA copy per cell. This may occur as a result of dilution by division of tumour cells. In a urokinase-type plasminogen activator /severe combined immunodeficiency (SCID) mouse model, hepatocyte proliferation induces loss of cccDNA (Lutgehetmann et al, 2010). Studies in the duck model showed that anti-viral therapy induces a greater cccDNA reduction in animals with higher hepatocyte proliferation (Addison et al, 2002). Our data also suggest a role for cell division in regulating cccDNA copy numbers, as small tumours were found to contain higher cccDNA copy number than large tumours; and in patients with small tumours, intra-tumoural cccDNA copy number did not differ from that in paired non-neoplastic liver.
Our observation of lower cccDNA levels within tumour is consistent with a recent report of 42 HBV-HCC patients (Wang et al, 2013a), but differs from an earlier report of 25 HBV-HCC patients showing higher cccDNA copy number in tumour than non-neoplastic liver (Wong et al, 2006). The reasons for the differences with the earlier study are unclear, although our study quantified HBV DNA and cccDNA copy number per cell by measuring and normalising to albumin copy numbers in all samples, rather than by calculating cell numbers based on the estimation of 6.667pg genomic DNA per cell (Wong et al, 2006).
Our previous retrospective studies showed that viral replicative activity in the non-neoplastic liver is associated with poor overall survival, but not recurrence-free survival, following surgical resection. This current study differs from that study in several ways. This current study includes surgical resection patients from 2008 to 2013, whereas the previous study included patients from 1991 to 2008. Several differences should also be noted such as smaller tumour size in the current study (4.3±3.6 vs 7.3±5.0cm). In this current study, viral DNA was measured using fresh surgical resection specimens, whereas in the previous study, viral DNA was measured using DNA isolated from archived paraffin blocks. Despite all of these differences, this current study confirms our previous observations that viral replicative activity correlates with HAI and Ishak stage, and that high viral replicative activity is associated with poor outcome. As these patients are continuing to be followed, the impact of viral replicative activity on recurrence-free survival and overall survival will be further evaluated.
In summary, this study presents a direct comparison of intracellular HBV DNA, cccDNA and viral replicative activity within tumour and non-neoplastic liver. The data showed that high replicative activity within the non-neoplastic liver is associated with higher HCC recurrence rate through mechanisms that are distinct from and independent of Ishak fibrosis stage.
This study was funded by NIH R21CA170722 (to JZ).
The authors declare no conflict of interest.
Supplementary Information accompanies this paper on British Journal of Cancer website (http://www.nature.com/bjc)
This work is published under the standard license to publish agreement. After 12 months the work will become freely available and the license terms will switch to a Creative Commons Attribution-NonCommercial-Share Alike 4.0 Unported License.