An enrichment protocol and analysis pipeline for long read sequencing of the hepatitis B virus transcriptome

Oxygen-dependent histone lysine demethylase 4 restricts hepatitis B virus replication

Mammalian cells have evolved strategies to regulate gene expression when oxygen is limited. Hypoxia-inducible factors (HIF) are the major transcriptional regulators of host gene expression. We previously reported that HIFs bind and activate hepatitis B virus (HBV) DNA transcription under low oxygen conditions; however, the global cellular response to low oxygen is mediated by a family of oxygenases that work in concert with HIFs. Recent studies have identified a role for chromatin modifiers in sensing cellular oxygen and orchestrating transcriptional responses, but their role in the HBV life cycle is as yet undefined. We demonstrated that histone lysine demethylase 4 (KDM4) can restrict HBV, and pharmacological or oxygen-mediated inhibition of the demethylase increases viral RNAs derived from both episomal and integrated copies of the viral genome. Sequencing studies demonstrated that KDM4 is a major regulator of the hepatic transcriptome, which defines hepatocellular permissivity to HBV infection. We propose a model where HBV exploits cellular oxygen sensors to replicate and persist in the liver. Understanding oxygen-dependent pathways that regulate HBV infection will facilitate the development of physiologically relevant cell-based models that support efficient HBV replication.

CTCF regulates hepatitis B virus cccDNA chromatin topology

Hepatitis B Virus (HBV) is a small DNA virus that replicates via an episomal covalently closed circular DNA (cccDNA) that serves as the transcriptional template for viral mRNAs. The host protein, CCCTC-binding factor (CTCF), is a key regulator of cellular transcription by maintaining epigenetic boundaries, nucleosome phasing, stabilisation of long-range chromatin loops and directing alternative exon splicing. We previously reported that CTCF binds two conserved motifs within Enhancer I of the HBV genome and represses viral transcription, however, the underlying mechanisms were not identified. We show that CTCF depletion in cells harbouring cccDNA-like HBV molecules and in de novo infected cells resulted in an increase in spliced transcripts, which was most notable in the abundant SP1 spliced transcript. In contrast, depletion of CTCF in cell lines with integrated HBV DNA had no effect on the abundance of viral transcripts and in line with this observation there was limited evidence for CTCF binding to viral integrants, suggesting that CTCF-regulation of HBV transcription is specific to episomal cccDNA. Analysis of HBV chromatin topology by Assay for Transposase Accessible Chromatin Sequencing (ATAC-Seq) revealed an accessible region spanning Enhancers I and II and the basal core promoter (BCP). Mutating the CTCF binding sites within Enhancer I resulted in a dramatic rearrangement of chromatin accessibility where the open chromatin region was no longer detected, indicating loss of the phased nucleosome up- and down-stream of the HBV enhancer/BCP. These data demonstrate that CTCF functions to regulate HBV chromatin conformation and nucleosomal positioning in episomal maintained cccDNA, which has important consequences for HBV transcription regulation.

The N6-methyladenosine demethylase ALKBH5 regulates the hypoxic HBV transcriptome

Chronic hepatitis B is a global health problem and current treatments only suppress hepatitis B virus (HBV) infection, highlighting the need for new curative treatments. Oxygen levels influence HBV replication and we previously reported that hypoxia inducible factors (HIFs) activate the basal core promoter to transcribe pre-genomic RNA. Application of a probe-enriched long-read sequencing method to map the HBV transcriptome showed an increased abundance of all viral RNAs under low oxygen or hypoxic conditions. Importantly, the hypoxic-associated increase in HBV transcripts was dependent on N6-methyladenosine (m6A) modifications and an m6A DRACH motif in the 5' stem loop of pre-genomic RNA defined transcript half-life under hypoxic conditions. Given the essential role of m6A modifications in the viral transcriptome we assessed the oxygen-dependent expression of RNA demethylases and bioinformatic analysis of published single cell RNA-seq of murine liver showed an increased expression of the RNA demethylase ALKBH5 in the peri-central low oxygen region. In vitro studies with a human hepatocyte derived HepG2 cell line showed increased ALKBH5 gene expression under hypoxic conditions. Silencing the demethylase reduced the levels of HBV pre-genomic RNA and host gene (CA9, NDRG1, VEGFA, BNIP3, FUT11, GAP and P4HA1) transcripts and this was mediated via reduced HIFα expression. In summary, our study highlights a previously unrecognized role for ALKBH5 in orchestrating viral and cellular transcriptional responses to low oxygen.

Hepatitis B virus serum RNA transcript isoform composition and proportion in chronic hepatitis B patients by nanopore long-read sequencing

Introduction

Serum hepatitis B virus (HBV) RNA is a promising new biomarker to manage and predict clinical outcomes of chronic hepatitis B (CHB) infection. However, the HBV serum transcriptome within encapsidated particles, which is the biomarker analyte measured in serum, remains poorly characterized. This study aimed to evaluate serum HBV RNA transcript composition and proportionality by PCR-cDNA nanopore sequencing of samples from CHB patients having varied HBV genotype (gt, A to F) and HBeAg status.

Methods

Longitudinal specimens from 3 individuals during and following pregnancy (approximately 7 months between time points) were also investigated. HBV RNA extracted from 16 serum samples obtained from 13 patients (73.3% female, 84.6% Asian) was sequenced and serum HBV RNA isoform detection and quantification were performed using three bioinformatic workflows; FLAIR, RATTLE, and a GraphMap-based workflow within the Galaxy application. A spike-in RNA variant (SIRV) control mix was used to assess run quality and coverage. The proportionality of transcript isoforms was based on total HBV reads determined by each workflow.

Results

All chosen isoform detection workflows showed high agreement in transcript proportionality and composition for most samples. HBV pregenomic RNA (pgRNA) was the most frequently observed transcript isoform (93.8% of patient samples), while other detected transcripts included pgRNA spliced variants, 3' truncated variants and HBx mRNA, depending on the isoform detection method. Spliced variants of pgRNA were primarily observed in HBV gtB, C, E, or F-infected patients, with the Sp1 spliced variant detected most frequently. Twelve other pgRNA spliced variant transcripts were identified, including 3 previously unidentified transcripts, although spliced isoform identification was very dependent on the workflow used to analyze sequence data. Longitudinal sampling among pregnant and post-partum antiviral-treated individuals showed increasing proportions of 3' truncated pgRNA variants over time.

Conclusions

This study demonstrated long-read sequencing as a promising tool for the characterization of the serum HBV transcriptome. However, further studies are needed to better understand how serum HBV RNA isoform type and proportion are linked to CHB disease progression and antiviral treatment response.