Genomic and oncogenic preference of HBV integration in hepatocellular carcinoma

Abstract

Hepatitis B virus (HBV) can integrate into the human genome, contributing to genomic instability and hepatocarcinogenesis. Here by conducting high-throughput viral integration detection and RNA sequencing, we identify 4,225 HBV integration events in tumour and adjacent non-tumour samples from 426 patients with HCC. We show that HBV is prone to integrate into rare fragile sites and functional genomic regions including CpG islands. We observe a distinct pattern in the preferential sites of HBV integration between tumour and non-tumour tissues. HBV insertional sites are significantly enriched in the proximity of telomeres in tumours. Recurrent HBV target genes are identified with few that overlap. The overall HBV integration frequency is much higher in tumour genomes of males than in females, with a significant enrichment of integration into chromosome 17. Furthermore, a cirrhosis-dependent HBV integration pattern is observed, affecting distinct targeted genes. Our data suggest that HBV integration has a high potential to drive oncogenic transformation.

ARTICLE

Received 29 Jan 2016 | Accepted 24 Aug 2016 | Published 5 Oct 2016

Ling-Hao Zhao1,2,*, Xiao Liu3,*, He-Xin Yan1,2,*, Wei-Yang Li3,4,*, Xi Zeng3,5,*, Yuan Yang1,*, Jie Zhao1, Shi-Ping Liu3, Xue-Han Zhuang3, Chuan Lin1, Chen-Jie Qin1,2, Yi Zhao4, Ze-Ya Pan1, Gang Huang1, Hui Liu1,Jin Zhang1, Ruo-Yu Wang1, Yun Yang1, Wen Wen1,2, Gui-Shuai Lv1,2, Hui-Lu Zhang1,2, Han Wu1,2, Shuai Huang1, Ming-Da Wang1,2, Liang Tang1,2, Hong-Zhi Cao3, Ling Wang6, T.P. Lee5, Hui Jiang3, Ye-Xiong Tan1,2, Sheng-Xian Yuan1, Guo-Jun Hou1, Qi-Fei Tao1, Qin-Guo Xu1, Xiu-Qing Zhang3, Meng-Chao Wu1, Xun Xu3,Jun Wang3,7, Huan-Ming Yang3, Wei-Ping Zhou1 & Hong-Yang Wang1,2,8

Hepatitis B virus (HBV) can integrate into the human genome, contributing to genomic instability and hepatocarcinogenesis. Here by conducting high-throughput viral integration detection and RNA sequencing, we identify 4,225 HBV integration events in tumour and adjacent non-tumour samples from 426 patients with HCC. We show that HBV is prone to integrate into rare fragile sites and functional genomic regions including CpG islands. We observe a distinct pattern in the preferential sites of HBV integration between tumour and non-tumour tissues. HBV insertional sites are signicantly enriched in the proximity of telomeres in tumours. Recurrent HBV target genes are identied with few that overlap. The overall HBV integration frequency is much higher in tumour genomes of males than in females, with a signicant enrichment of integration into chromosome 17. Furthermore, a cirrhosis-dependent HBV integration pattern is observed, affecting distinct targeted genes. Our data suggest that HBV integration has a high potential to drive oncogenic transformation.

1 Eastern Hepatobiliary Surgery Hospital, Second Military Medical University, Shanghai 200438, China. 2 National Center for Liver Cancer, Shanghai 200438, China. 3 BGI-Shenzhen, Shenzhen 518083, China. 4 School of Bioscience and Bioengineering, South China University of Technology, Guangzhou 510006, China. 5 Faculty of Medicine, School of Biomedical Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong 999077, China. 6 Department of Vascular and Endocrine Surgery, Xijing Hospital, Fourth Military Medical University, Xi An 710032, China. 7 Department of Biology, University of Copenhagen, Ole Maales Vej 5,2200 Copenhagen, Denmark. 8 State Key Laboratory of Oncogenes and Related Genes, Cancer Institute of Renji Hospital, Shanghai Jiaotong University, Shanghai 200127, China. * These authors contributed equally to this work. Correspondence and requests for materials should be addressed to H.-Y.W. (email: mailto:[email protected]

Web End [email protected] ) or to W.-P.Z. (email: mailto:[email protected]

Web End [email protected] ) or to H.-M.Y. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

DOI: 10.1038/ncomms12992 OPEN

Genomic and oncogenic preference of HBV integration in hepatocellular carcinoma

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992

Hepatocellular carcinoma (HCC) ranks fth in global cancer incidence and represents the third leading cause of cancer deaths1. Chronic infection with hepatitis B or C

virus (HBV or HCV) represents the major risk factors for the development of HCC2,3. Unlike HCV, an RNA virus which never integrates into the host genome during its lifecycle, HBVDNA frequently integrates into host genome and progressively contributes to hepatocarcinogenesis4. Integration of HBV DNA into the host genome causes genetic damage and chromosomal instability, which is known to be selectively advantageous for tumour progression5. The deregulated host gene expression or a subset of rearranged integration sites has been shown to play a key role in HCC development6. In addition, expression of viral proteins such as X protein and S antigen as well as their oncogenic mutants may further enable the acquisition of neoplastic characteristics. Previous isolation of HBV integration sites using PCR-based methods7 and more recent use of deep sequencing in a small size cohort of HCC patients suggests that the HBV insertional sites occurred randomly throughout the genome, leading to the presumption that there were no preferential sites of integration. However, recurrent HBV integrations have recently been identied in a large cohort of HCC patients through the application of high-throughput next-generation sequencing811. Although recurrent integrations still represented a minority of the total events characterized, the occurrence of recurrent sites of insertion suggests that HBV may have preferential integration sites associated with distinct biological consequences and clinical outcomes. The most common HBV integration event is located at the telomerase reverse transcriptase gene (TERT), which is thought to confer early clonal advantage during chronic HBV infection. Other frequently targeted genes such as KMT2B and CCNE1 have recently been identied as new classes of genes not previously known to play a causal role in cancer10. However, it is of note that recurrent HBV integration sites identied in previous studies may only represent a minority of recurrent events.

HCC has a male predominance and is closely related to cirrhosis1214. Males with chronic HBV infection are at higher risk of developing HCC compared with females. In fact, the sex difference becomes apparent from the stage of early chronic infection, where the male-to-female ratio increased from 1.2 in asymptomatic carriers to 6.3 in chronic hepatitis and 9.8 in HCC in a Taiwanese study15,16. On the other hand, although the majority of HCC patients with HBV infection have concurrent liver cirrhosis, the integration of the viral genome into infected cells can directly induce a non-cirrhotic liver to develop HCC. Despite extensive research, the precise mechanisms whereby HBV integration contributes to cirrhotic and non-cirrhotic hepatocarcinogenesis as well as gender disparity remain largely unclear due to technical limitations and a lack of clinical annotation.

In this study, we conducted high-throughput viral integration detection (HIVID), a highly sensitive method for assaying viral insertion compared with the whole-genome sequencing, and analysed 426 HBV-associated HCC patients with or without cirrhosis17. Our results revealed an excessive HBV integration across the host genome with particular genomic pattern in a gender or cirrhosis-dependent manner.

ResultsCharacteristics of the HBV integration events in human genome. To search for HBV integration sites across the HCC genome and decode viralhost interactions, we conducted HIVID analysis on tumour and adjacent non-tumour liver genomes extracted in a cohort including 426 clinically and pathologically

well-characterized HBV-associated HCC cases: 298 with cirrhosis and 128 without cirrhosis; 360 males and 66 females (Supplementary Figs 13 and Supplementary Table 1). HBV integrations were effectively detected at a single base pair resolution using HBV capture sequencing. A total of 4,225 HBV integration breakpoints were discovered. To conrm the newly discovered events, we randomly selected 180 putative insertions for PCR analysis and successfully validate 81.1% of these integration sites (Supplementary Data 1). Despite a handful of integration hotspots, most insertion sites are distributed in an apparently random manner throughout the genome (Fig. 1a). The integrations maps generated reveal a heterogeneous, widespread viral integration landscape in tumour as well as in non-tumour liver tissue from HCC patients. However, HCC tumour samples and their adjacent non-tumour liver tissues exhibit strikingly distinct patterns of viral insertion. Following the known association between HBV infection and HCC development, tumour samples harbour a much higher frequency of integrated reads than the non-tumour samples. The total HBV integration breakpoints in tumour and non-tumour samples are 3,486 and 739 events, respectively (Supplementary Data 2 and 3). Of 426 paired samples, 76.9% tumour and 37.6% non-tumour tissues contained HBV integration. The average number of integration sites in tumour tissues and adjacent non-tumour liver tissues are10.6 and 4.6, respectively, (P 3.3 10 6, unpaired Students

t test; Supplementary Fig. 4). Analysis of the distribution of hotspots across the human genome further conrms the prevalence of HBV integrations in tumours. Among the 826 and 303 genes with HBV insertions (intragenic region and upstream 10 kb) detected in tumour and non-tumour tissues, respectively, only 64 genes were shared by the tumour and non-tumour samples, underscoring that integration patterns are distinct in the tumour and normal samples (Fig. 1a and Supplementary Data 4).

Correlation between HBV integration and genome instability. Statistically signicant enrichment of integration was observed within CpG islands in tumours compared with normal samples (Fig. 1b). The integration frequency dramatically decreased in genomic loci moving away from CpG islands. However, no apparent prevalence in integration was observed in transcription factor binding site and transcription start site (Supplementary Fig. 5). HBV was also prone to integrate into rare fragile sites, but only in a tumour-specic context (Fig. 1c). These results suggest preferential HBV integration into chromosomal repetitive or fragile regions could provide a selective advantage during tumorigenesis, rather than representing by product event from random insertion. More importantly, HBV integration sites were signicantly enriched in the proximity of telomere in the tumour samples but not in the non-tumour samples (Fig. 2), suggesting that HBV favours targeting chromosomal elements critical for the maintenance of chromosome stability.

As integration of HBVDNA into the human genome is considered an early event during hepatocarcinogenesis and can induce chromosomal instability, we next analysed the distribution of HBV breakpoints across individual chromosomes. Although relatively random distribution of integration sites was observed in non-tumour samples, statistically signicant enrichment of integration sites on chromosome 5, 16, 17 and 19 was observed in tumour samples (Fig. 2), suggesting that preferential integration exists at the chromosomal level.

Expression of target genes and effect to HBV integration. The comprehensive HIVID analysis in this large-sample cohort has allowed us to investigate the frequency of recurring tumour-associated integrations in genes. Notably, there were 88 and 17

2 NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992 ARTICLE

0.08 CpG

0.06

0.04

a

b

P<0.01

P<0.01

Y

1

X

30

M

22

21

20

19

18

210

240

2

Ratio

Ratio

0.02

0.00

0.6

0.5

0.4

0.3 NS

NS

30

60

30

60

90

3

17

Tumour Normal Expected

Breakpoints of fragile region

16

15

14

60

90

90

60

60

30

c

90

Tumour Normal Expected

30

5

4

90

60

13

P<0.01

30

NS

90

12

60

6

11

150

30

7

60

P<0.01

10

9

0.2

0.1

0.0 Common Rare NFRs

8

NS

Figure 1 | Distribution of HBV integration breakpoints throughout the human genomes in 426 paired samples. (a) Distribution of integration breakpoints across the human genome in 426 paired samples. Each bar represents the sample frequency of HBV integration breakpoints at a particular locus in the human genome (hg19). Tumour (red) and non-tumour (dark blue) samples with HBV integrations are shown on the inner and outer circles, respectively. Histogram axis units represent number of samples. Some loci with a high frequency of integration are marked. GENE. (b) Comparison of the breakpoints in the CpG island region of 426 paired samples. The expected (assuming uniform, random distribution, yellow) and the observed (actual numbers, tumour: blue; normal liver tissue: purple) percentages of HBV integration breakpoints of tumour and non-tumour samples in CpG islands region are shown. P values were calculated by w2 test. (c) Distribution of integration breakpoints in the fragile region (FR) of 426 paired tumor and non-tumor tissues. The expected (assuming uniform, random distribution, yellow) and the observed (actual numbers, tumour: blue; normal liver tissue: purple) ratios of HBV integration breakpoints in common fragile region, rare fragile region and non-fragile region are shown. P values were calculated by w2 test. Common: common FR; rare: rare FR; NFRs: non-fragile regions.

genes recurrently affected (nZ2) by HBV integration in the tumour and non-tumour samples, respectively, which has dramatically expanded our understanding of integration sites of HBV. Apart from the recurrent integrations previously observed in TERT, MLL4 and CCNE1 genes, novel affected genes (PTPRD, UNC5D, NRG3, CTNND2 and AHRR) were identied, and their expressions were altered at both the transcript and protein levels (Supplementary Figs 6 and 7). Interestingly, these genes have been reported to have oncogenic or tumour suppressive function. In keeping with previously published ndings, the expression of some affected genes with oncogenic potential, including TERT and KMT2B, showed increased expression following HBV integration. Our results suggest that these novel affected genes are, at least in part, functionally relevant to HBV insertion.

We noted that individuals with HBV integrations (BK40) in the tumour displayed a signicantly shorter survival than those without HBV integrations (BK 0). Furthermore, HBV

integration led to an evident upregulation of TERT and patients with HBV integration at TERT gene had signicantly poorer survival (Fig. 3). Interestingly, tumour-enriched integration into the promoter regions of genes was noted, as compared with the non-tumour genomes (385/3486 versus 57/739; Po0.01, w2 test).

In contrast, the integration breakpoints were preferentially

located in the intron region in non-tumour samples (Supplementary Fig. 8). This integration bias indicates that HBV is prone to integrate into promoters in HCC tumours, further affecting transcription of particular genes.

Breakpoints in HBV genome. To investigate the integration mechanisms underlying viralhost interactions, we surveyed the breakpoints on the HBV genome. Consistent with previous reports, B40% of breakpoints were observed at nucleotides 1,4001,900 around the 30-end of the HBx and 50-end of the Precore/Core genes in both tumour and non-tumour samples(Fig. 4a). To deconvolute HBV integration at the transcriptional level, we randomly selected 12 tumour samples for RNA sequencing and characterized the transcribed viral elements corresponding to nucleotides 1,7001,900 of HBV genome(Supplementary Table 3). Interestingly, a peak of HBV breakpoints in a region from 300 nt to 500 nt was noted at the transcriptional level, where the S gene was located(Fig. 4b). These viral breakpoints may not only rearrange viral DNA and functional properties but also physically subvert normal control of nearby cellular genes, thus leading to the dysregulation of transcription network in HCC.

NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992

Chromosomal ends

**

0.08

Expected_ratio Observed_ratio_of_normal Observed_ratio_of_tumor

Ratio

0.04

*

*

* *

*

* *

0.00

Chr1 Chr2 Chr3 Chr4 Chr5 Chr6 Chr7 Chr8 Chr9 Chr10 Chr11 Chr12 Chr13 Chr14 Chr15 Chr16 Chr17 Chr18 Chr19 Chr20 Chr21

Chr22

0.00

0.10

Ratio

* * *

0.20

0.30

*

Whole chromosomal

Figure 2 | Chromosome enrichment and chromosomal ends enrichment of HBV integration in human genome in 426 paired samples. Each bar of whole-chromosome represents the expected (assuming uniform, random distribution, green) and the observed (actual numbers, tumour: blue, normal liver tissue: purple) ratio of HBV integration breakpoints at a particular chromosome in human genome. Ratios are numbered. Each bar of chromosomal ends represents the expected (assuming uniform, random distribution, green) and the observed (actual numbers, tumor: blue, normal liver tissue, purple) ratio of HBV integration breakpoints at the 2 M region of chromosomal ends in human genome. Ratios are numbered. Dark red star represents statistically signicant difference between normal liver samples and random distribution. Blue star represents statistically signicant difference between tumour samples and random distribution. P values were calculated by w2 test.

a b c

1.0

1.0

0.8

0.6

0.4

0.2

0.0

BK=0

BK>0

n=97 NO-TERT

n=315

n=101

200

P=0.0009088

Probability of survival

0.8

0.6

0.4

0.2

0.0

Experssion level (T/N)

150

100

0

n=319

Probability of survival

TERT

50

P=3.8E0.2

P=9.4E0.3

0 10 20 30 40

Months after surgery Months after surgery

50 0 10 20 30 40 50

Int

HBV integration (TERT)

Non_int

Figure 3 | Clinical correlation analysis of HBV integration in HCC. (a) KaplanMeier survival curves for individuals with (BK40, n 319) versus without

(BK 0, n 97) HBV integration breakpoints by log-rank test. Those who lacked prognostic information were excluded from the analysis (n 10).

(b) Gene expression levels of TERT that frequently harboured HBV integrations in samples with versus without HBV integration events. Gene expression was normalized by the corresponding adjacent, normal control and is represented as the tumour/normal gene expression level. P values of unpaired Students t test are shown. In the box plots, the median (50th percentile) is the middle line, with the bottom and top of the box representing the 25th and 75th percentiles of the data, respectively. The ends of the whiskers represent the lowest and highest data within the 1.5 interquartile range (IQR). IQR was dened as the distance between the lower and upper quartiles of the data. (c) KaplanMeier survival curves for individuals with (n 101) versus without

(n 315) HBV integration in TERT by log-rank test.

Next, we investigated the prevalence of various HBV genotypes across our cohort. The sequence analysis of HBVDNA from 426 patients with HBV integration revealed the predominant integration of genotype C in 74.9% patients and genotype B in remaining 13.4% patients. The distribution of HBV genotypes in non-HCC liver tissues was similar to tumour tissues. We also noted that the integration rate of B type was signicantly higher than C type (Supplementary Fig. 9). No signicant correlation between the circulating level of HBVDNA and the frequency of HBV integration was noted (Supplementary Fig. 10). However,

we observed that circulating HBe antigen level had a signicant linear correlation with HBV integration in tumour tissues. This suggests that the viral replication capacity could closely inuence the HBV integration in human genome (Supplementary Fig. 11).

In our analysis, we observed that multiple HBV fragments could integrate into a single site within the human genome. Further analysis demonstrated that, when compared with non-tumour samples, the breakpoints within tumour samples displayed more diversity with regard to the HBV integration fragments. Although most of the breakpoints were integrated by

4 NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992 ARTICLE

a b

2,400

2,500

2,600

2,700

2,400

2,500

2,700

2,800

2,900

3,000

2,800

2,900

3,000

3,100

T-bk-num DNA

RNA

HBV

2,600

N-bk-num T-sample-num N-sample-num

2,100

2,200

2,300

2,100

2,200

2,300

2,000

2,000

1,900

1,900

3,100

1,800

1,800

3,200

3,200

1,700

1,700

HBV

1,600

1,600

100

100

1,500

1,500

200

200

1,400

1,400

300

300

1,300

1,300

400

400

1,200

500

1,200

500

600

700

1,100

600

1,100

1,000

900

800

700

1,000

900

800

Figure 4 | Distribution of integration breakpoints in the HBV genome. (a) Distribution of integration breakpoints in the HBV genome in 426 paired samples. Each bar represents the number of HBV integration breakpoints (Tumour: blue, Normal liver tissue: red) or sample frequency (tumour: green, normal liver tissue: yellow) at a particular locus in HBV genome. Histograms were constructed for 100-bp intervals. Histogram axis units represent number of breakpoints, and outer DNA numbering is given in bases. HBV genes with different functions are shown. (b) Distribution of integration breakpoints of both DNA and RNA level of 12 paired samples. Histograms were constructed for 100-bp intervals. The number of HBV integration of RNA (red) and DNA (blue) level are shown on the inner and outer circles. Histogram axis units represent number of breakpoints, and outer DNA numbering is given in bases. HBV genes with different functions are shown.

less than 2 types of HBV fragments, tumour tissues carried more breakpoints integrated by more than 2 types of HBV fragments than non-tumour tissues did (Supplementary Fig. 12), suggesting an increased heterogeneity for HBV integration in tumours.

The mechanism of HBV integration has been thought to involve double-strand break repair, linearized viral DNA invasion/end joining and viral replication, but the details are still unclear18. Microhomology (MH)-mediated DNA repair pathways, including MH-mediated end joining, fork stalling and template switching and MH-mediated break-induced repair have been proposed to induce genomic rearrangements19,20 and other viral insertion21. In light of these hypotheses, we searched for the MH sequences between the cellular and the inserted HBV DNA near the integration sites, and found they were signicantly enriched with the increase of the MH length, suggesting the potential involvement of MH-mediated mechanism (Supplementary Fig. 13), which could be triggered by the genomic instability/fragility near the integration sites (Supplementary Fig. 14).

Association of HBV integration with sex and clinical outcome. As HBV integrations may be specic biomarkers for prediction of clinical outcomes, we analysed the association between viral integration and clinicopathological parameters in HCCs (Fig. 5 and Supplementary Table 4). Consistent with a male predominance in HCCs, males had much higher HBV integrations than females in tumour genomes (Supplementary Table 5). Specically, we observed that substantial genomic regions were preferably integrated by HBV in males in Chromosomes 2 and 17, but not in females (Fig. 6), and HBV integrations in these regions were closely associated with prognosis (Supplementary Fig. 15). Notably, many genes in chromosome 17p are known to play important roles in hepatocarcinogenesis22 including the key tumour suppressor TP53 and master inhibitor of liver regeneration MKK4 (refs 2325). Furthermore, we found integrations in males were prone to locate into the region of

core protein (P 0.01, w2 test) in the HCC genome when

compared with females.

Correlation of HBV integration with Cirrhosis. Clinically, some patients with HBV-related HCC have a lower incidence of cirrhosis, suggesting that critical oncogenes or tumour suppressor genes were involved in hepatocarcinogenesis without inicting chronic inammation. Gene annotation revealed a cirrhosis-dependent HBV integration pattern; that is, tumours arising from non-cirrhotic liver displayed a signicantly enriched viral integration in the vicinity of putative oncogenes or tumour suppressors such as KMT2B, CCNE1 and AHRR compared with those in the cirrhotic liver (Supplementary Data 5 ), suggesting that constitutive activation or inactivation of these genes may contribute to the early onset of HCC without cirrhotic responses.

DiscussionHere we demonstrate a large scale analysis of HBVDNA integration sites in liver cancer using HIVID approach. Compared with whole-genome sequencing, HIVID combines HBV fragment capture and massively parallel sequencing, thus dramatically increases the detection efciency and decreases sequencing complexity and cost. In line with published results on whole-genome sequencing of HBV-associated HCC, HBV integration occurs frequently in both tumour and adjacent nontumour tissues, with signicantly higher insertion rates in tumours811. Although integrations of HBV were once thought to be random, our data highlights recurrent integrated regions and greatly expands the list of affected genes. In addition to frequent integrations at TERT, KMT2B and CCNE1, which were known targets of HBV insertion, we discovered massively new recurrent HBV integrations among 426 HCC patients. A signicant percentage of them showed expression or function abnormality in tumours and have been shown to be related to cancer progression in other tumour types. Our data provide an important amendment to our understanding of HBV impact on

NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992

200

150

Number of integration

100

50

0

A

B

C

BCLC stage

Gender

TERT

Cirrhosis

Tumour_size

Metastasis

Degree

HBV_type

CCNA2

Integration

CCNE1

DDX11L1

KMT2B

Gender Female Male

Noncirrhosis Cirrhosis Size < 5 cm Size 5 cm

Nonmetastasis Metastasis Degree III Degree IIIIV

HBVB HBVC HBVB,HBVC Other

Number=0 0 < number < 5 5 number < 10 number 10

Figure 5 | Clinical annotation of HBV integration sites in 426 HCC samples. All panels are aligned with vertical tracks representing 426 individuals. The data are sorted by BCLC stage, gender, cirrhosis, tumour size, metastasis, Edmonson-Steiner classication and HBV type. The bottom heat map shows the distribution of HBV integrations into the ve recurrent HBV targeted genes in HCC samples.

Cirrhosis

Tumour_size

Metastasis

Degree

HBV_type

hepatocarcinogenesis, supporting the non-random integration and selection process of HBV targeted genes.

As HBV integration is common and is associated with genomic instability, we examined the distribution of insertion sites in distinct genomic elements26,27. Not surprisingly, chromosomal fragile sites, known sites of genomic rearrangement in cancer, represent preferred sites for the integration of HBV. In addition, HBV showed an evident preference for CpG islands, in which DNA methylation often precedes the appearance of tumours and are prone to additional chromosomal aberrations, such as loss of heterozygosity. As virus negative tumours have fewer epigenetic and genetic alterations28, epigenetic instability via CpG islands methylation/demethylation may contribute to chromosomal instability in HBV-related HCCs. More importantly, the enrichment of HBV integration events was also observed in the proximity of telomere, which plays an important role in maintaining genome stability. Regardless of the mechanism,

dysfunction of telomere function can lead to extensive DNA amplication, large terminal deletions and generate many types of rearrangements commonly associated with human cancers29.

Although specic sequence features near the fusion breakpoints were not identied in the host genome, there were a handful of hotspots in the tumour genome where multiple HBV fragments could integrate. This is consistent with a non-random integration model followed by a positive selection during hepatocarcinogenesis. Other than HBx proteins surrounding DR1 or DR2 sites, RNA sequencing revealed an enriched expression and favoured integration of preS and S protein in the tumour samples, in accordance with the role of the truncated preS2/S protein in mediating stimulation of cellular genes, in trans30.

Despite a signicant gender skew of HCC in male, little is known about the difference of HBV integration between genders. Herein, we show that HBV was prone to integrate in male tumour

6 NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992 ARTICLE

Chr1 Chr2 Chr3 Chr4 Chr5 Chr6 Chr7 Chr8 Chr9 Chr10 Chr11 Chr12 Chr13 Chr14 Chr15 50

40

30

Breakpoint number

a

b

P value = 0.002

P value = 0.005

0

20

40

60

80

100

120 120 120 120 120 120 120 120 120 120 120 120

140 140 140 140 140 140 140 140 140

160 160 160 160 160 160

180 180 180 180 180

200

220 220

240 240

200

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

0

20

40

60

80

100

100

Male

Male

Female

Female

c

Chr16 Chr17 Chr18 Chr19 Chr20 Chr21 Chr22

0

20

40

FemaleMale Femalemale N-region

0

20

40

60

80

0

20

40

60

80

0

20

40

0

20

40

0

20

40

0

20

40

Integration ratio

20

10

0

1.0

0.8

0.6

0.4

0.2

0.0

60 60

Figure 6 | HBV integration preference in male versus female HCC samples. (a) Distribution of the enriched HBV integration regions in tumour samples. The red bars indicate the enriched regions found only in female samples, green bars represent enriched regions only found in male samples, and blue bars mean enriched regions found in both. In every 1M window of chromosomal region, the observed number of integration sites was tested against the expected number where all integrations distributed randomly across the genome. Enrichment was dened as P value smaller than 0.05 calculated by w2 test, which was performed separately in male and female samples. (b) Comparison of the numbers of breakpoints in tumour tissues of male and female samples. The HBV breakpoints number of male (left) and female (right) are shown. The box plots show the median (horizontal bar), 25th and 75th percentiles, and the whiskers of the plots show the smallest and largest values. P values was calculated by unpaired Students t test. (c) Comparison of integration ratio of male and female HCC samples. The HBV integration ratio of male (left) and female (right) HCC samples are shown. P values were calculated by w2 test.

genome and, more importantly, had a chromosome preference. HBV integration into chromosomes2 and 17 was favoured in male tumours, the latter of which contains the TP53 gene and is frequently lost in the human cancers. In addition to aberrations of TP53 gene on chromosome 17p13.1, other genes on 17p13.3 may also play a role in hepatocarcinogenesis22. Notably, located on 17p adjacent to the TP53 gene, MKK4 was recently identied as a master inhibitor of liver regeneration23,24. As MKK4 is one of the most consistently mutated genes across tumour types25, concurrent inactivation of TP53 and MKK4 following HBV integration may constitute an enabling mechanism for genome destabilization and hepatocarcinogenesis. Consistently, patients with HBV integration in chromosomes 2 and 17 displayed signicantly poorer survival than in other chromosomes.

Other than inactivating tumour suppressive programs, HBV integrations into different genomic locations can result in a substantial enhancement of distinct oncogenic genes, which may have a profound inuence on phenotypic characteristics. In support of this notion, a signicant proportion of HBV-associated HCCs develop in the absence of cirrhosis. We identied several important oncogenes preferentially affected by HBV integration in the non-cirrhotic HCC samples, suggesting that constitutive activation or inactivation of these genes may contribute to the early onset of HCC without inicting cirrhotic responses. Further work is warranted to delineate the complex network of these HBV targeted genes and its overall impact on cellular phenotype.

Collectively, our work provides a large scale and unbiased HBV integration map in HCC, revealing the preference of integration

occurring within regions of the genome prone to DNA mutations or rearrangements and novel target genes recurrently affected by HBV integration. Together with the inserted HBV fragments, these characteristics may endow HBV integration a greater opportunity to induce crucial oncogenic alterations to host genes in a gender- and cirrhosis-dependent manner, and eventually leading to HCC development in patients with chronic HBV infection.

Methods

Sample and DNA extraction. All samples were obtained from Eastern Hepatobiliary Surgery Hospital, Shanghai, between 2009 and 2010. Tumour tissues and paired adjacent non-tumour tissue were resected from patients undergoing primary hepatectomy by experienced surgeons. Specimens were immediately cryopreserved in 80 C following resection. All operations were carried out carefully to avoid

contamination between samples. DNA was extracted from paired tumour and normal adjacent samples using Genomic DNA Mini Kit (Invitrogen, Life Technologies) according to the manufactorys instruction. All HBVHCCs, and paired adjacent non-tumour tissue are clinically and pathologically characterized by two independent pathologists and radiologists. Demographic and clinical features of the patients were systematically collected and summarized in Supplementary Table 1. This study was approved by the ethical committee of EHBH hospital and informed consent was obtained from each patient. The inclusion criteria for this study included: (i) HBV-positive HCCs and paired adjacent non-tumour tissues (ii) obtained from consenting patients and (iii) all samples are HCV-negative and HIV-negative and (iv) without for autoimmune hepatitis and metabolic and/or genetics disorders such as Wilsons disease, hemochromatosis.

RNA extraction and real-time quantitative PCR. Total RNAs were extracted from tumours or paired adjacent non-tumour samples using MagMAX mirVana

NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992

Total RNA Isolation Kit (Ambion, Life Technologies) according to the manufac-torys instruction. The concentration and quality of RNAs were determined by Nanodrop 2000 (Thermo Fisher Scientic). The quantitative PCRs were performed using Power SYBR-Green PCR Master Mix (Life Technologies, Carlsbad, CA, USA), primer pair sets and Applied Biosystems ViiA7 Real-time System (Life Technologies) under the conditions of 10 min at 95 C followed by 40 cycles at 95 C for 15 s and at 60 C for 60 s. Beta actin was used as a housekeeping control. The sequences of the primers are listed in Supplementary Table 2.

HBV capture experiment. We designed the sequence-capture probes according to eight types of HBV genome (A, B, C, D, E, F, G and H) sequences and these probes were produced by MyGenostics. DNA fragments B150200 bp were obtained after shearing genomic DNA by Covaris E-210 (Covaris, Inc., Woburn, MA). These fragments were puried, end blunted, A tailed and adaptor ligated. The products of ligation were amplied by PCR and constructed into DNA sequencing library. Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA ) was used to quantify the concentration of the DNA library. Next, hybridization was performed according to the instruction of Target Enrichment Protocol (GenCapTM Enrichment, MyGenostics, USA). This process was carried out at 65 C for 24 h and the un-targeted fragments were then removed by washing buffer. These eluted fragments were amplied by PCR and were further processed to Paired end 100-bp read-length for sequencing by the HiSeq 2000 sequencer (Illumina Inc., San Diego, CA).

Breakpoints detection and annotation of HBV integration sites. HIVID pipeline was used for the breakpoints detection17. First, clean reads were obtained through removing low quality and duplicated reads, as well as adaptor-contaminated reads. Burrows-Wheeler Aligner (BWA) was used to align clean reads onto human (NCBI build 37, HG19) and HBV genome (AF090842.1, AB033554.1, AB014381.1, M32138.1, AB032431.1, AB036910.1, AB064310 and AY090454.1). The paired end reads that can be mapped to human or HBV reference genome with both ends were removed. The remained reads were paired-end assembled to reconstruct fragment sequences, which was to locate the position of breakpoints more precisely. Subsequently, reference mapping was performed again to re-map the paired-end-assembled reads onto human and HBV genome using BWA (ref. 31). The position of a breakpoint was dened as the junction of human and HBV sequence in a paired-end-assembled read. To minimize the impact of different sequencing data amount of each sample, the normalized support-reads number (Norm value) was introduced, which was equal to the supported reads number of each HBV breakpoint per million clean read pairs. The breakpoints with Norm value Z2 were retained. The Sanger sequencing validation rate for the selected breakpoints was81.1%. ANNOVAR was used to do the annotation for the integrated breakpoints17,31.

Process of detecting the integration breakpoints by RNA sequencing. We randomly selected 12 tumour samples undergoing HIVID analysis to perform RNA sequencing. Total RNA isolated with TRIzol reagent was treated with RNase-free DNaseI (New England BioLabs) at 37 C for 10 min. The Dynabeads mRNA Purication Kit (Life Technologies) was used to isolate mRNA from the total RNA samples. RNA-seq libraries were sequenced as paired-end 100 bp sequence tags using the standard Solexa pipeline. Integration sites were analysed by using the transcriptome data according to previously described method17. We removed reads that perfectly aligned to human or HBV genome and reserved chimeric paired-end reads. These chimeric reads were composed of both human genome sequence and HBV genome sequence, and were used to identify HBV integration breakpoints in the transcriptome.

PCR and sanger sequencing validation. PCR and Sanger sequencing were used to verify the selected HBV integration breakpoints from HIVID. All the samples used for SANGER validation are left-over samples from previous rounds of the capture approach. PCR primers were designed based on the paired-end assembled fragment, in which one primer located in human genome and the other in HBV genome. PCR were performed by GeneAmp PCR System 9700 thermal cycler and then preceded to Sanger sequencing on Applied Biosystems 3730 DNA analyzer

(Life Technologies, Inc.).

Immunohistochemical staining. Formalin xed and parafn-embedded sections (4 mm) were subjected to immunohistochemical staining. The slides were incubated overnight at 4 C with primary antibodies, including rabbit anti-AHRR (1:100, ab108518, Abcam), rabbit anti-PTPRD(1:100, LS-B9625, LifeSpan Biosciences), rabbit anti-NRG3(1:200, ab83704, Abcam) and mouse anti-UNC5D (1:100, ab58141, Abcam). Anti-rabbit or anti-mouse horseradish peroxidase-conjugated secondary antibodies (Santa Cruz Biotechnology, Santa Cruz, CA) were applied. Finally, diaminobenzidine colorimetric reagent solution from Dako (Carpinteria, CA) was used and followed by hematoxylin counterstaining (Sigma Chemical Co). Tissue slides were scanned with an Aperio ScanScope GL, and the Aperio ImageScope software (Aperio Technologies, Vista,CA) was used to assess the

scanned images based on the percentage of positively stained cells and staining intensity. Expression levels of these proteins in all clinical samples were quantied.

Data availability. DNA sequencing data and RNA-seq data that support this study have been deposited in the Sequence Read Archive (SRA) database under the accession codes SRA335342 and SRA447498, respectively. The authors declare that all other data is present in the Article and its Supplementary Information Files or available from the authors upon request.

References

1. Forner, A., Llovet, J. M. & Bruix, J. Hepatocellular carcinoma. Lancet 379, 12451255 (2012).

2. Siegel, R., Naishadham, D. & Jemal, A. Cancer statistics, 2013. CA CancerJ. Clin. 63, 1130 (2013).3. Chen, J. G. & Zhang, S. W. Liver cancer epidemic in China: past, present and future. Semin. Cancer Biol. 21, 5969 (2011).

4. Shafritz, D. A. et al. Integration of hepatitis B virus DNA into the genome of liver cells in chronic liver disease and hepatocellular carcinoma. Studies in percutaneous liver biopsies and post-mortem tissue specimens. N. Engl. J. Med. 305, 10671073 (1981).

5. Hanahan, D. & Weinberg, R. A. Hallmarks of cancer: the next generation. Cell 144, 646674 (2011).

6. Lau, C. C. et al. Viral-human chimeric transcript predisposes risk to liver cancer development and progression. Cancer Cell 25, 335349 (2014).7. Matsubara, K. & Tokino, T. Integration of hepatitis B virus DNA and its implications for hepatocarcinogenesis. Mol. Biol. Med. 7, 243260 (1990).

8. Sung, W. K. et al. Genome-wide survey of recurrent HBV integration in hepatocellular carcinoma. Nat. Genet. 44, 765769 (2012).

9. Ding, D. et al. Recurrent targeted genes of hepatitis B virus in the liver cancer genomes identied by a next-generation sequencing-based approach. PLoS Genet. 8, e1003065 (2012).

10. Li, X. et al. The function of targeted host genes determines the oncogenicity of HBV integration in hepatocellular carcinoma. J. Hepatol. 60, 975984 (2014).

11. Jiang, Z. et al. The effects of hepatitis B virus integration into the genomes of hepatocellular carcinoma patients. Genome Res. 22, 593601 (2012).

12. Yuen, M. F., Hou, J. L. & Chutaputti, A. Hepatocellular carcinoma in the Asia pacic region. J. Gastroenterol. Hepatol. 24, 346353 (2009).

13. Chu, C. M. et al. Basal core promoter mutation is associated with progression to cirrhosis rather than hepatocellular carcinoma in chronic hepatitis B virus infection. Br. J. Cancer 107, 20102015 (2012).

14. Perz, J. F. et al. The contributions of hepatitis B virus and hepatitis C virus infections to cirrhosis and primary liver cancer worldwide. J. Hepatol. 45, 529538 (2006).

15. Chu, C. M. et al. Sex difference in chronic hepatitis B virus infection: an appraisal based on the status of hepatitis B e antigen and antibody. Hepatology 3, 947950 (1983).

16. Huang, Y. T. et al. Lifetime risk and sex difference of hepatocellular carcinoma among patients with chronic hepatitis B and C. J. Clin. Oncol. 29, 36433650 (2011).

17. Li, W. et al. HIVID: an efcient method to detect HBV integration using low coverage sequencing. Genomics 102, 338344 (2013).

18. Hino, O., Ohtake, K. & Rogler, C. E. Features of two hepatitis B virus (HBV) DNA integrations suggest mechanisms of HBV integration. J. Virol. 63, 26382643 (1989).

19. Verdin, H. et al. Microhomology-mediated mechanisms underlie non-recurrent disease-causing microdeletions of the FOXL2 gene or its regulatory domain. PLoS Genet. 9, e1003358 (2013).

20. Liu, P. et al. Mechanisms for recurrent and complex human genomic rearrangements. Curr. Opin. Genet. Dev. 22, 211220 (2012).

21. Hu, Z. et al. Genome-wide proling of HPV integration in cervical cancer identies clustered genomic hot spots and a potential microhomology-mediated integration mechanism. Nat. Genet. 47, 158163 (2015).

22. Zhao, X. et al. The minimum LOH region dened on chromosome 17p13.3 in human hepatocellular carcinoma with gene content analysis. Cancer Lett. 190, 221232 (2003).

23. Yeasmin, S. et al. MKK4 acts as a potential tumor suppressor in ovarian cancer. Tumour Biol. 32, 661670 (2011).

24. Wuestefeld, T. et al. A Direct in vivo RNAi screen identies MKK4 as a key regulator of liver regeneration. Cell 153, 389401 (2013).

25. Cunningham, S. C. et al. Targeted deletion of MKK4 in cancer cells: a detrimental phenotype manifests as decreased experimental metastasis and suggests a counterweight to the evolution of tumor-suppressor loss. Cancer Res. 66, 55605564 (2006).

26. Feitelson, M. A. & Lee, J. Hepatitis B virus integration, fragile sites, and hepatocarcinogenesis. Cancer Lett. 252, 157170 (2007).

8 NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms12992 ARTICLE

27. Matsuzaki, Y. et al. HBV genome integration and genetic instability in HBsAg-negative and anti-HCV-positive hepatocellular carcinoma in Japan. Cancer Lett. 119, 5361 (1997).

28. Katoh, H. et al. Epigenetic instability and chromosomal instability in hepatocellular carcinoma. Am. J. Pathol. 168, 13751384 (2006).

29. Bailey, S. M. & Murnane, J. P. Telomeres, chromosome instability and cancer. Nucleic Acids Res. 34, 24082417 (2006).

30. Kekule, A. S. et al. The preS2/S region of integrated hepatitis B virus DNA encodes a transcriptional transactivator. Nature 343, 457461 (1990).

31. Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 17541760 (2009).

Acknowledgements

Supported in part by grants: State key infection disease project of China (2012ZX10002); National Key Basic Research Program of China (2014CB542102); NSFC(81221061, 91329000, 31371440, 30672227, 81372207 , 81201940 and 81272899); the National High Technology Research and Development Program of China 2013AA032202 and 2012AA02A201; Shanghai Eastern Scholarship (201232) and New Excellent Youth and Academic leader Program (16XD1403300 and XYQ2013074); Guangdong Innovative Research Team Program (No. 2009010016); Guangdong Enterprise Key Laboratory of Human Disease Genomics; ShenZhen Engineering Laboratory for Clinical Molecular Diagnostic; China National GeneBank-Shenzhen.

Author contributions

H.Y.W., W.-P.Z., H.M.Y. and J.W. conceived the project, established and coordinated the consortium of studies. L.-H.Z., X.L., H.-X.Y. and W.Y.L. designed and executed the analysis, interpreted the ndings and wrote the rst drafts of the manuscript. X.Z. and Y.Y. contributed to analysis design and execution. J.Z., S.P.L., C.L., Z.Y.P., G.H., H.L. and J.Z. contributed to data collection and preparation. Y.Z., H.L.Z., C.J.Q. and

X.H.Z. contributed to interpretation of ndings, and manuscript preparation. M.C.W. contributed to study design and gave important intellectual content of the study. G.S.L., H.Z.C., L.W., T.L.L. and H.J. contributed to data preparation and analysis. R.Y.W., Y.Y., W.W., H.W., S.H., M.D.W. and L.T. contributed to data analysis and manuscript preparation. Y.X.T., S.X.Y., Q.F.T., G.J.H. and Q.G.X. contributed to data collection and preparation. X.Q.Z. and X.X. contributed to interpretation of ndings.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Zhao, L-H. et al. Genomic and oncogenic preference of HBV integration in hepatocellular carcinoma. Nat. Commun. 7, 12992doi: 10.1038/ncomms12992 (2016).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/by/4.0/

r The Author(s) 2016

NATURE COMMUNICATIONS | 7:12992 | DOI: 10.1038/ncomms12992 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9