A Fibrosis-Independent Hepatic Transcriptomic Signature Identifies Drivers of Disease Progression in Primary Sclerosing Cholangitis.

BACKGROUND AND AIMS: Primary sclerosing cholangitis (PSC) is a heterogeneous cholangiopathy characterized by progressive biliary fibrosis. RNA sequencing of liver tissue from patients with PSC (n = 74) enrolled in a 96-week clinical trial was performed to identify associations between biological pathways that were independent of fibrosis and clinical events. APPROACH AND RESULTS: The effect of fibrosis was subtracted from gene expression using a computational approach. The fibrosis-adjusted gene expression patterns were associated with time to first PSC-related clinical event (e.g., cholangitis, hepatic decompensation), and differential expression based on risk groups and Ingenuity Pathway Analysis were performed. Baseline demographic data were representative of PSC: median age 48 years, 71% male, 49% with inflammatory bowel disease, and 44% with bridging fibrosis or cirrhosis. The first principle component (PC1) of RNA-sequencing data accounted for 18% of variance and correlated with fibrosis stage (ρ = -0.80; P < 0.001). After removing the effect of fibrosis-related genes, the first principle component was not associated with fibrosis (ρ = -0.19; P = 0.11), and a semisupervised clustering approach identified two distinct patient clusters with differential risk of time to first PSC-related event (P < 0.0001). The two groups had similar fibrosis stage, hepatic collagen content, and α-smooth muscle actin expression by morphometry, Enhanced Liver Fibrosis score, and serum liver biochemistry, bile acids, and IL-8 (all P > 0.05). The top pathways identified by Ingenuity Pathway Analysis were eukaryotic translation inhibition factor 2 (eIF2) signaling and regulation of eIF4/p70S6K signaling. Genes involved in the unfolded protein response, activating transcription factor 6 (ATF6) and eIF2, were differentially expressed between the PSC clusters (down-regulated in the high-risk group by log-fold changes of -0.18 [P = 0.02] and -0.16 [P = 0.02], respectively). Clinical events were enriched in the high-risk versus low-risk group (38% [12/32] vs. 2.4% [1/42], P < 0.0001). CONCLUSIONS: Removing the contribution of fibrosis-related pathways uncovered alterations in the unfolded protein response, which were associated with liver-related complications in PSC.

An integrative analysis of colon cancer identifies an essential function for PRPF6 in tumor growth.

The spliceosome machinery is composed of multimeric protein complexes that generate a diverse repertoire of mRNA through coordinated splicing of heteronuclear RNAs. While somatic mutations in spliceosome components have been discovered in several cancer types, the molecular bases and consequences of spliceosome aberrations in cancer are poorly understood. Here we report for the first time that PRPF6, a member of the tri-snRNP (small ribonucleoprotein) spliceosome complex, drives cancer proliferation by preferential splicing of genes associated with growth regulation. Inhibition of PRPF6 and other tri-snRNP complex proteins, but not other snRNP spliceosome complexes, selectively abrogated growth in cancer cells with high tri-snRNP levels. High-resolution transcriptome analyses revealed that reduced PRPF6 alters the constitutive and alternative splicing of a discrete number of genes, including an oncogenic isoform of the ZAK kinase. These findings implicate an essential role for PRPF6 in cancer via splicing of distinct growth-related gene products.

DNA Methylation and Immune Cell Markers Demonstrate Evidence of Accelerated Aging in Patients with Chronic Hepatitis B Virus or Hepatitis C Virus, with or without Human Immunodeficienct Virus Co-infection.

BACKGROUND: Several chronic diseases accelerate biological aging. We investigated age acceleration and the association between peripheral blood DNA methylation (DNAm) and immune cell markers in patients chronically infected with the hepatitis B virus (HBV) or the hepatitis C virus (HCV) with and without human immunodeficiency virus (HIV) co-infection. METHODS: Age acceleration was measured as the difference between epigenetic age (Horvath clock) and chronological age. The immune marker model of age acceleration was developed using Elastic Net regression to select both the immune markers and their associated weights in the final linear model. RESULTS: Patients with chronic HBV (n = 51) had a significantly higher median epigenetic age compared to chronological age (age accelerated) (P < .001). In patients with chronic HCV infection (n = 63), age acceleration was associated with liver fibrosis as assessed by histology (P < .05), or presence of HIV co-infection (P < .05), but not HCV mono-infection. Age acceleration defined by immune markers was concordant with age acceleration by DNA methylation (correlation coefficient = .59 in HBV; P = .0025). One-year treatment of HBV patients with nucleoside therapy was associated with a modest reduction in age acceleration, as measured using the immune marker model (-.65 years, P = .018). CONCLUSION: Our findings suggest that patients with chronic viral hepatitis have accelerated epigenetic aging, that immune markers define biological age, and have the potential to assess the effects of therapeutic intervention on age acceleration.

Longitudinal analysis of serum microRNAs as predictors of cirrhosis regression during treatment of hepatitis B virus infection.

BACKGROUND AND AIMS: Most patients with cirrhosis induced by chronic HBV infection experience fibrosis regression after long-term antiviral treatment, while some remain cirrhotic. Fibrosis regression is associated with lower odds of developing hepatic decompensation and hepatocellular carcinoma, but mechanisms impacting differential fibrosis regression between individuals are unclear. We asked whether soluble molecules, including serum microRNAs, could serve as biomarkers of fibrosis regression. METHODS: We analysed cryopreserved sera from clinical trials in which cirrhotic HBV-infected patients (baseline Ishak fibrosis score of 5-6) received 240 weeks of nucleotide analogue treatment. Liver biopsies at week 240 in these trials showed 71/96 patients (74%) had fibrosis regression (Ishak ≤ 4) while 25/96 (26%) remained cirrhotic (Ishak 5-6). We quantified inflammatory markers (CXCL10, soluble CD163) and miRNAs (n = 179) from serum at baseline, week 48 and week 240 of treatment in a sub-cohort of patients with (n = 14) or without (n = 14) fibrosis regression. RESULTS: CXCL10, sCD163 and miRNAs previously associated with HBV replication and inflammation decreased during treatment but did not differ based on fibrosis regression. Two miRNAs (miR-421 and miR-454-3p) had lower baseline expression in patients with subsequent fibrosis regression. In all, 27 miRNAs differed at week 240 and had higher expression in patients with fibrosis regression (eg miR-199a-3p, miR-423-3p, miR-142-3p, miR-let-7d-5p). Several miRNAs (miR-141-3p, let-7d-5p) that correlated with regression have previously been implicated in the pathophysiology of non-alcoholic steatohepatitis. CONCLUSIONS: In cirrhotic patients with chronic HBV infection treated with antiviral therapy, serum miRNAs have differential expression based on fibrosis regression, suggesting potential utility as biomarkers.

Endosomes are specialized platforms for bacterial sensing and NOD2 signalling.

The detection of microbial pathogens involves the recognition of conserved microbial components by host cell sensors such as Toll-like receptors (TLRs) and NOD-like receptors (NLRs). TLRs are membrane receptors that survey the extracellular environment for microbial infections, whereas NLRs are cytosolic complexes that detect microbial products that reach the cytosol. Upon detection, both sensor classes trigger innate inflammatory responses and allow the engagement of adaptive immunity. Endo-lysosomes are the entry sites for a variety of pathogens, and therefore the sites at which the immune system first senses their presence. Pathogens internalized by endocytosis are well known to activate TLRs 3 and 7-9 that are localized to endocytic compartments and detect ligands present in the endosomal lumen. Internalized pathogens also activate sensors in the cytosol such as NOD1 and NOD2 (ref. 2), indicating that endosomes also provide for the translocation of bacterial components across the endosomal membrane. Despite the fact that NOD2 is well understood to have a key role in regulating innate immune responses and that mutations at the NOD2 locus are a common risk factor in inflammatory bowel disease and possibly other chronic inflammatory states, little is known about how its ligands escape from endosomes. Here we show that two endo-lysosomal peptide transporters, SLC15A3 and SLC15A4, are preferentially expressed by dendritic cells, especially after TLR stimulation. The transporters mediate the egress of bacterially derived components, such as the NOD2 cognate ligand muramyl dipeptide (MDP), and are selectively required for NOD2 responses to endosomally derived MDP. Enhanced expression of the transporters also generates endosomal membrane tubules characteristic of dendritic cells, which further enhanced the NOD2-dependent response to MDP. Finally, sensing required the recruitment of NOD2 and its effector kinase RIPK2 (refs 8, 9) to the endosomal membrane, possibly by forming a complex with SLC15A3 or SLC15A4. Thus, dendritic cell endosomes are specialized platforms for both the lumenal and cytosolic sensing of pathogens.

Evolutionary toggling of the MAPT 17q21.31 inversion region.

Using comparative sequencing approaches, we investigated the evolutionary history of the European-enriched 17q21.31 MAPT inversion polymorphism. We present a detailed, BAC-based sequence assembly of the inverted human H2 haplotype and compare it to the sequence structure and genetic variation of the corresponding 1.5-Mb region for the noninverted H1 human haplotype and that of chimpanzee and orangutan. We found that inversion of the MAPT region is similarly polymorphic in other great ape species, and we present evidence that the inversions occurred independently in chimpanzees and humans. In humans, the inversion breakpoints correspond to core duplications with the LRRC37 gene family. Our analysis favors the H2 configuration and sequence haplotype as the likely great ape and human ancestral state, with inversion recurrences during primate evolution. We show that the H2 architecture has evolved more extensive sequence homology, perhaps explaining its tendency to undergo microdeletion associated with mental retardation in European populations.

Tumour-secreted miR-9 promotes endothelial cell migration and angiogenesis by activating the JAK-STAT pathway.

Angiogenesis plays a crucial role during tumorigenesis and much progress has been recently made in elucidating the role of VEGF and other growth factors in the regulation of angiogenesis. Recently, microRNAs (miRNAs) have been shown to modulate a variety of physiogical and pathological processes. We identified a set of differentially expressed miRNAs in microvascular endothelial cells co-cultured with tumour cells. Unexpectedly, most miRNAs were derived from tumour cells, packaged into microvesicles (MVs), and then directly delivered to endothelial cells. Among these miRNAs, we focused on miR-9 due to the strong morphological changes induced in cultured endothelial cells. We found that exogenous miR-9 effectively reduced SOCS5 levels, leading to activated JAK-STAT pathway. This signalling cascade promoted endothelial cell migration and tumour angiogenesis. Remarkably, administration of anti-miR-9 or JAK inhibitors suppressed MV-induced cell migration in vitro and decreased tumour burden in vivo. Collectively, these observations suggest that tumour-secreted miRNAs participate in intercellular communication and function as a novel pro-angiogenic mechanism.

The mutation spectrum revealed by paired genome sequences from a lung cancer patient.

Lung cancer is the leading cause of cancer-related mortality worldwide, with non-small-cell lung carcinomas in smokers being the predominant form of the disease. Although previous studies have identified important common somatic mutations in lung cancers, they have primarily focused on a limited set of genes and have thus provided a constrained view of the mutational spectrum. Recent cancer sequencing efforts have used next-generation sequencing technologies to provide a genome-wide view of mutations in leukaemia, breast cancer and cancer cell lines. Here we present the complete sequences of a primary lung tumour (60x coverage) and adjacent normal tissue (46x). Comparing the two genomes, we identify a wide variety of somatic variations, including >50,000 high-confidence single nucleotide variants. We validated 530 somatic single nucleotide variants in this tumour, including one in the KRAS proto-oncogene and 391 others in coding regions, as well as 43 large-scale structural variations. These constitute a large set of new somatic mutations and yield an estimated 17.7 per megabase genome-wide somatic mutation rate. Notably, we observe a distinct pattern of selection against mutations within expressed genes compared to non-expressed genes and in promoter regions up to 5 kilobases upstream of all protein-coding genes. Furthermore, we observe a higher rate of amino acid-changing mutations in kinase genes. We present a comprehensive view of somatic alterations in a single lung tumour, and provide the first evidence, to our knowledge, of distinct selective pressures present within the tumour environment.

The effects of hepatitis B virus integration into the genomes of hepatocellular carcinoma patients.

Hepatitis B virus (HBV) infection is a leading risk factor for hepatocellular carcinoma (HCC). HBV integration into the host genome has been reported, but its scale, impact and contribution to HCC development is not clear. Here, we sequenced the tumor and nontumor genomes (>80× coverage) and transcriptomes of four HCC patients and identified 255 HBV integration sites. Increased sequencing to 240× coverage revealed a proportionally higher number of integration sites. Clonal expansion of HBV-integrated hepatocytes was found specifically in tumor samples. We observe a diverse collection of genomic perturbations near viral integration sites, including direct gene disruption, viral promoter-driven human transcription, viral-human transcript fusion, and DNA copy number alteration. Thus, we report the most comprehensive characterization of HBV integration in hepatocellular carcinoma patients. Such widespread random viral integration will likely increase carcinogenic opportunities in HBV-infected individuals.

Genome and transcriptome sequencing of lung cancers reveal diverse mutational and splicing events.

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

A burst of segmental duplications in the genome of the African great ape ancestor.

It is generally accepted that the extent of phenotypic change between human and great apes is dissonant with the rate of molecular change. Between these two groups, proteins are virtually identical, cytogenetically there are few rearrangements that distinguish ape-human chromosomes, and rates of single-base-pair change and retrotransposon activity have slowed particularly within hominid lineages when compared to rodents or monkeys. Studies of gene family evolution indicate that gene loss and gain are enriched within the primate lineage. Here, we perform a systematic analysis of duplication content of four primate genomes (macaque, orang-utan, chimpanzee and human) in an effort to understand the pattern and rates of genomic duplication during hominid evolution. We find that the ancestral branch leading to human and African great apes shows the most significant increase in duplication activity both in terms of base pairs and in terms of events. This duplication acceleration within the ancestral species is significant when compared to lineage-specific rate estimates even after accounting for copy-number polymorphism and homoplasy. We discover striking examples of recurrent and independent gene-containing duplications within the gorilla and chimpanzee that are absent in the human lineage. Our results suggest that the evolutionary properties of copy-number mutation differ significantly from other forms of genetic mutation and, in contrast to the hominid slowdown of single-base-pair mutations, there has been a genomic burst of duplication activity at this period during human evolution.

Copy number variation of CCL3-like genes affects rate of progression to simian-AIDS in Rhesus Macaques (Macaca mulatta).

Variation in genes underlying host immunity can lead to marked differences in susceptibility to HIV infection among humans. Despite heavy reliance on non-human primates as models for HIV/AIDS, little is known about which host factors are shared and which are unique to a given primate lineage. Here, we investigate whether copy number variation (CNV) at CCL3-like genes (CCL3L), a key genetic host factor for HIV/AIDS susceptibility and cell-mediated immune response in humans, is also a determinant of time until onset of simian-AIDS in rhesus macaques. Using a retrospective study of 57 rhesus macaques experimentally infected with SIVmac, we find that CCL3L CNV explains approximately 18% of the variance in time to simian-AIDS (p<0.001) with lower CCL3L copy number associating with more rapid disease course. We also find that CCL3L copy number varies significantly (p<10(-6)) among rhesus subpopulations, with Indian-origin macaques having, on average, half as many CCL3L gene copies as Chinese-origin macaques. Lastly, we confirm that CCL3L shows variable copy number in humans and chimpanzees and report on CCL3L CNV within and among three additional primate species. On the basis of our findings we suggest that (1) the difference in population level copy number may explain previously reported observations of longer post-infection survivorship of Chinese-origin rhesus macaques, (2) stratification by CCL3L copy number in rhesus SIV vaccine trials will increase power and reduce noise due to non-vaccine-related differences in survival, and (3) CCL3L CNV is an ancestral component of the primate immune response and, therefore, copy number variation has not been driven by HIV or SIV per se.

Distinct resistance mechanisms arise to allosteric vs. ATP-competitive AKT inhibitors.

The AKT kinases have emerged as promising therapeutic targets in oncology and both allosteric and ATP-competitive AKT inhibitors have entered clinical investigation. However, long-term efficacy of such inhibitors will likely be challenged by the development of resistance. We have established prostate cancer models of acquired resistance to the allosteric inhibitor MK-2206 or the ATP-competitive inhibitor ipatasertib following prolonged exposure. While alterations in AKT are associated with acquired resistance to MK-2206, ipatasertib resistance is driven by rewired compensatory activity of parallel signaling pathways. Importantly, MK-2206 resistance can be overcome by treatment with ipatasertib, while ipatasertib resistance can be reversed by co-treatment with inhibitors of pathways including PIM signaling. These findings demonstrate that distinct resistance mechanisms arise to the two classes of AKT inhibitors and that combination approaches may reverse resistance to ATP-competitive inhibition.

Characterization of six human disease-associated inversion polymorphisms.

The human genome is a highly dynamic structure that shows a wide range of genetic polymorphic variation. Unlike other types of structural variation, little is known about inversion variants within normal individuals because such events are typically balanced and are difficult to detect and analyze by standard molecular approaches. Using sequence-based, cytogenetic and genotyping approaches, we characterized six large inversion polymorphisms that map to regions associated with genomic disorders with complex segmental duplications mapping at the breakpoints. We developed a metaphase FISH-based assay to genotype inversions and analyzed the chromosomes of 27 individuals from three HapMap populations. In this subset, we find that these inversions are less frequent or absent in Asians when compared with European and Yoruban populations. Analyzing multiple individuals from outgroup species of great apes, we show that most of these large inversion polymorphisms are specific to the human lineage with two exceptions, 17q21.31 and 8p23 inversions, which are found to be similarly polymorphic in other great ape species and where the inverted allele represents the ancestral state. Investigating linkage disequilibrium relationships with genotyped SNPs, we provide evidence that most of these inversions appear to have arisen on at least two different haplotype backgrounds. In these cases, discovery and genotyping methods based on SNPs may be confounded and molecular cytogenetics remains the only method to genotype these inversions.

Recurrent rearrangements of chromosome 1q21.1 and variable pediatric phenotypes.

BACKGROUND: Duplications and deletions in the human genome can cause disease or predispose persons to disease. Advances in technologies to detect these changes allow for the routine identification of submicroscopic imbalances in large numbers of patients. METHODS: We tested for the presence of microdeletions and microduplications at a specific region of chromosome 1q21.1 in two groups of patients with unexplained mental retardation, autism, or congenital anomalies and in unaffected persons. RESULTS: We identified 25 persons with a recurrent 1.35-Mb deletion within 1q21.1 from screening 5218 patients. The microdeletions had arisen de novo in eight patients, were inherited from a mildly affected parent in three patients, were inherited from an apparently unaffected parent in six patients, and were of unknown inheritance in eight patients. The deletion was absent in a series of 4737 control persons (P=1.1x10(-7)). We found considerable variability in the level of phenotypic expression of the microdeletion; phenotypes included mild-to-moderate mental retardation, microcephaly, cardiac abnormalities, and cataracts. The reciprocal duplication was enriched in nine children with mental retardation or autism spectrum disorder and other variable features (P=0.02). We identified three deletions and three duplications of the 1q21.1 region in an independent sample of 788 patients with mental retardation and congenital anomalies. CONCLUSIONS: We have identified recurrent molecular lesions that elude syndromic classification and whose disease manifestations must be considered in a broader context of development as opposed to being assigned to a specific disease. Clinical diagnosis in patients with these lesions may be most readily achieved on the basis of genotype rather than phenotype.

Shotgun sequence assembly and recent segmental duplications within the human genome.

Complex eukaryotic genomes are now being sequenced at an accelerated pace primarily using whole-genome shotgun (WGS) sequence assembly approaches. WGS assembly was initially criticized because of its perceived inability to resolve repeat structures within genomes. Here, we quantify the effect of WGS sequence assembly on large, highly similar repeats by comparison of the segmental duplication content of two different human genome assemblies. Our analysis shows that large (> 15 kilobases) and highly identical (> 97%) duplications are not adequately resolved by WGS assembly. This leads to significant reduction in genome length and the loss of genes embedded within duplications. Comparable analyses of mouse genome assemblies confirm that strict WGS sequence assembly will oversimplify our understanding of mammalian genome structure and evolution; a hybrid strategy using a targeted clone-by-clone approach to resolve duplications is proposed.

The structure and evolution of centromeric transition regions within the human genome.

An understanding of how centromeric transition regions are organized is a critical aspect of chromosome structure and function; however, the sequence context of these regions has been difficult to resolve on the basis of the draft genome sequence. We present a detailed analysis of the structure and assembly of all human pericentromeric regions (5 megabases). Most chromosome arms (35 out of 43) show a gradient of dwindling transcriptional diversity accompanied by an increasing number of interchromosomal duplications in proximity to the centromere. At least 30% of the centromeric transition region structure originates from euchromatic gene-containing segments of DNA that were duplicatively transposed towards pericentromeric regions at a rate of six-seven events per million years during primate evolution. This process has led to the formation of a minimum of 28 new transcripts by exon exaptation and exon shuffling, many of which are primarily expressed in the testis. The distribution of these duplicated segments is nonrandom among pericentromeric regions, suggesting that some regions have served as preferential acceptors of euchromatic DNA.

Methylation signatures in peripheral blood are associated with marked age acceleration and disease progression in patients with primary sclerosing cholangitis.

BACKGROUND & AIMS: A DNA methylation (DNAm) signature derived from 353 CpG sites (the Horvath clock) has been proposed as an epigenetic measure of chronological and biological age. This epigenetic signature is accelerated in diverse tissue types in various disorders, including non-alcoholic steatohepatitis, and is associated with mortality. Here, we assayed whole blood DNAm to explore age acceleration in patients with primary sclerosing cholangitis (PSC). METHODS: Using the MethylationEPIC BeadChip (850K) array, DNAm signatures in whole blood were analyzed in 36 patients with PSC enrolled in a 96-week trial of simtuzumab (Ishak F0-1, n = 13; F5-6, n = 23). Age acceleration was calculated as the difference between DNAm age and chronological age. Comparisons between patients with high and low age acceleration (≥ vs. < the median) were made and Cox regression evaluated the association between age acceleration and PSC-related clinical events (e.g. decompensation, cholangitis, transplantation). RESULTS: Age acceleration was significantly higher in patients with PSC compared to a healthy reference cohort (median, 11.1 years, p <2.2 × 10-16). In PSC, demographics, presence of inflammatory bowel disease, and ursodeoxycholic acid use were similar between patients with low and high age acceleration. However, patients with high age acceleration had increased serum alkaline phosphatase, gamma glutamyltransferase, alanine aminotransferase, enhanced liver fibrosis test scores, and greater hepatic collagen and α-smooth muscle actin expression on liver biopsy (all p <0.05). Moreover, patients with high age acceleration had an increased prevalence of cirrhosis (89% vs. 39%; p = 0.006) and greater likelihood of PSC-related events (hazard ratio 4.19; 95% CI 1.15-15.24). CONCLUSION: This analysis of blood DNAm profiles suggests that compared with healthy controls, patients with PSC - particularly those with cirrhosis - exhibit significant acceleration of epigenetic age. Future studies are required to evaluate the prognostic implications and effect of therapies on global methylation patterns and age acceleration in PSC. LAY SUMMARY: An epigenetic clock based on DNA methylation has been proposed as a marker of age. In liver diseases such as non-alcoholic steatohepatitis, age acceleration based on this epigenetic clock has been observed. Herein, we show that patients with primary sclerosing cholangitis have marked age acceleration, which is further accentuated by worsening fibrosis. This measure of age acceleration could be a useful marker for prognostication or risk stratification in primary sclerosing cholangitis.

Genomic modeling of hepatitis B virus integration frequency in the human genome.

Hepatitis B infection is a world-wide public health burden causing serious liver complications. Previous studies suggest that hepatitis B integration into the human genome plays a crucial role in triggering oncogenic process and may also constitutively produce viral antigens. Despite the progress in HBV biology and sequencing technology, our fundamental understanding of how many hepatocytes in the liver actually carry viral integrations is still lacking. Herein we provide evidence that the HBV virus integrates with a lower-bound frequency of 0.84 per diploid genome in hepatitis B positive hepatocellular cancer patients. Moreover, we calculate that integrated viral DNA generates ~80% of the HBsAg transcripts in these patients. These results underscore the need to re-evaluate the clinical end-point and treatment strategies for chronic hepatitis B patients.

Large-scale viral genome analysis identifies novel clinical associations between hepatitis B virus and chronically infected patients.

Despite the high global prevalence of chronic hepatitis B (CHB) infection, datasets covering the whole hepatitis B viral genome from large patient cohorts are lacking, greatly limiting our understanding of the viral genetic factors involved in this deadly disease. We performed deep sequencing of viral samples from patients chronically infected with HBV to investigate the association between viral genome variation and patients' clinical characteristics. We discovered novel viral variants strongly associated with viral load and HBeAg status. Patients with viral variants C1817T and A1838G had viral loads nearly three orders of magnitude lower than patients without those variants. These patients consequently experienced earlier viral suppression while on treatment. Furthermore, we identified novel variants that either independently or in combination with precore mutation G1896A were associated with the transition from HBeAg positive to the negative phase of infection. These observations are consistent with the hypothesis that mutation of the HBeAg open reading frame is an important factor driving CHB patient's HBeAg status. This analysis provides a detailed picture of HBV genetic variation in the largest patient cohort to date and highlights the diversity of plausible molecular mechanisms through which viral variation affects clinical phenotype.