Hypoxia inducible factors regulate hepatitis B virus replication by activating the basal core promoter.

BACKGROUND & AIMS: Hypoxia inducible factors (HIFs) are a hallmark of inflammation and are key regulators of hepatic immunity and metabolism, yet their role in HBV replication is poorly defined. HBV replicates in hepatocytes within the liver, a naturally hypoxic organ, however most studies of viral replication are performed under conditions of atmospheric oxygen, where HIFs are inactive. We therefore investigated the role of HIFs in regulating HBV replication.
METHODS: Using cell culture, animal models, human tissue and pharmacological agents inhibiting the HIF-prolyl hydroxylases, we investigated the impact of hypoxia on the HBV life cycle.
RESULTS: Culturing liver cell-based model systems under low oxygen uncovered a new role for HIFs in binding HBV DNA and activating the basal core promoter, leading to increased pre-genomic RNA and de novo HBV particle secretion. The presence of hypoxia responsive elements among all primate members of the hepadnaviridae highlights an evolutionary conserved role for HIFs in regulating this virus family.
CONCLUSIONS: Identifying a role for this conserved oxygen sensor in regulating HBV transcription suggests that this virus has evolved to exploit the HIF signaling pathway to persist in the low oxygen environment of the liver. Our studies show the importance of considering oxygen availability when studying HBV-host interactions and provide innovative routes to better understand and target chronic HBV infection. LAY
SUMMARY: Viral replication in host cells is defined by the cellular microenvironment and one key factor is local oxygen tension. Hepatitis B virus (HBV) replicates in the liver, a naturally hypoxic organ. Hypoxia inducible factors (HIFs) are the major sensors of low oxygen; herein, we identify a new role for these factors in regulating HBV replication, revealing new therapeutic targets.

See Editorial, pages 16–18

Introduction

HBV is a global health problem, with more than 250 million people chronically infected and at least 880,000 deaths/year from HBV-related liver diseases such as cirrhosis and hepatocellular carcinoma (HCC) (WHO, Global hepatitis Report 2017). The HBV vaccine has no effect on chronic infection and current treatments suppress viral replication but are not curative. Chronic infection is associated with a blunted host innate immune response, high viral antigen expression and exhausted antiviral T cell responses that fail to control HBV replication. In most cases, treatment may be life-long and patients with a functional cure may still develop liver cancer. Effective antiviral drugs have revolutionized treatments for hepatitis C virus and there is a growing impetus to identify curative therapies for HBV.

HBV is the prototype member of the hepadnaviridae family that replicate via episomal copies of a covalently closed circular DNA (cccDNA) genome. cccDNA is frequently referred to as a viral mini-chromosome, where gene expression is regulated by DNA methylation, host RNA Polymerase II and transcription factors (reviewed in). Viral replication is primarily determined by the size of the cccDNA pool and its transcriptional activity, and host factors such as hepatocyte nuclear factor 4 alpha (HNF4α) have been reported to regulate HBV replication (reviewed in). The viral genome has 4 promoters (the basal core promoter [BCP], Sp1, Sp2 and Xp) that transcribe 6 major viral RNAs of decreasing length with a common 3’ polyadenylation signal. These RNAs include: pre-core (pC) that encodes the e antigen (HBeAg); pre-genomic (pgRNA) that is translated to yield core protein (HBcAg) and polymerase; preS1, preS2 and S RNAs encoding the surface envelope glycoproteins and the X transcript for the multi-functional x protein (HBx). Encapsidated pgRNA is reverse-transcribed by the viral polymerase to generate new DNA genomes that can be reimported to the nucleus to maintain the cccDNA pool or are enveloped and secreted as infectious particles. Defining the role of host factors that regulate HBV pgRNA genesis and half-life will increase our understanding of the HBV life cycle and enable the design of more effective antiviral approaches.

Viral gene expression is shaped by the cellular microenvironment and one important factor to consider is local oxygen tension. The liver receives oxygenated blood from the hepatic artery and oxygen-depleted blood via the hepatic portal vein, resulting in an oxygen gradient of 8–4% between the periportal and pericentral areas, respectively. This oxygen gradient associates with liver zonation, a phenomenon where hepatocytes show distinct functional and structural heterogeneity across the liver. Mammalian cells adapt to low oxygen through an orchestrated transcriptional response regulated by hypoxia inducible factor (HIF). This transcription factor, comprising HIF-1β and HIF-1α or HIF-2α subunits, is regulated by oxygen-dependent and independent stress signals such as inflammation and oxidative stress., The heterodimeric HIF complex binds a consensus RCGTG(C) motif or hypoxic responsive element (HRE) in the promoter and enhancer regions of responsive genes. When oxygen is abundant, newly synthesised HIFα subunits are hydroxylated by prolyl-hydroxylase domain (PHD) proteins and rapidly targeted for proteasomal degradation. In contrast, when oxygen is limited, these HIFα subunits translocate to the nucleus, dimerize with HIF-1β and regulate the transcription of host genes involved in cell metabolism and immune regulation. The majority of reports studying HBV replication in vitro are performed at atmospheric oxygen levels (18%) where HIFs are inactive, so their role in viral replication has been overlooked. We evaluated the effect of HIF signaling on HBV replication and uncovered a positive role for HIFs in activating viral transcription that could inform future therapeutic strategies.

Materials and methods

HBV, cells and hypoxic culture

HBV (Genotype D) was purified from HepAD38 cells and the infectious titer measured as previously reported. HepG2 cells expressing sodium taurocholate co-transporting polypeptide (NTCP) (Stefan Urban, Heidelberg University) were cultured in DMEM with 100 U/ml penicillin/0.1 mg/ml streptomycin/10% FCS/glutamax/non-essential amino acids. HepG2-NTCP cells were seeded on collagen coated plasticware and infected with HBV (multiplicity of infection 200–400 genome equivalents/cell) with 4% polyethylene glycol 8000 for 6 h. Viral inoculum was removed, cells washed 3x with PBS and infected cells maintained in DMEM-10% FCS in the presence or absence of 2.5% DMSO. For hypoxic cultures, cells were incubated in a Galaxy 48R incubator (Eppendorf) with variable oxygen tension; unless otherwise stated all incubations were for 72 h. Normoxic cells were cultured at 5% CO2 and 18% oxygen.

For further details regarding the materials and methods used, please refer to the CTAT table and supplementary information.

Results

Hypoxia inducible factors activate the HBV basal core promoter

To investigate a role for HIFs in the HBV life cycle, we first assessed whether the viral genome encodes HREs. Screening >7,000 HBV sequences available in the HBV database identified HREs within Enhancer I (1236–1240) and an antisense motif between Enhancer I and II (EnhI/II) (1604–1599) (Fig. 1A). It is interesting to note that both HREs are present in HBV sequences obtained from Bronze age and medieval samples. Analysis of hepadnaviridae NCBI referent genomes shows that, with the exception of HBV genotype H, all human, higher primate, and woodchuck viruses encode both motifs (Fig. 1A). Monkey, bat, rodent and avian viral genomes lack the motifs, suggesting a conserved evolutionary role for HIFs in the regulation of primate members of this family. To evaluate this hypothesis we stabilized HIF expression using the licensed PHD inhibitor FG-4592 (Roxadustat). Since HIFα isoform expression can vary between cell types we show that FG-4592 induced HIF-1α and HIF-2α expression in HepG2 cells derived from a human HCC (Fig. 1B). Transfecting HepG2 cells with plasmids encoding EnhI/II and the basal core promoter (HBV-Luc) or control HRE-luciferase showed that FG-4592 induced a significant, time-dependent activation of both reporters (Fig. 1B).

Fig. 1

HIFs positively regulate HBV transcription.

(A) Conservation of HREs amongst diverse HBV lineages. Conservation of hypoxia responsive elements at positions 1236–1240 and 1605–1599 were assessed amongst 7,313 HBV sequences and consensus plots generated. NCBI referent full-length hepadnaviridae genome sequences were aligned and a Neighbor Joining tree constructed, where branches highlighted in red denote HRE containing genomes. (B) HIFs induce HBV promoter activity. HepG2-NTCP cells expressing an HRE or HBV EnhI/II and BCP construct driving luciferase (HBV-Luc) were treated with PHD inhibitor FG-4592 (30 μM) and promoter activity assessed at 24 h intervals for 72 h. HIF-1α and HIF-2α expression was confirmed by western blot. Statistical analysis was performed using a 2-Way ANOVA (∗∗p ≤0.01, ∗∗∗∗p ≤0.0001). (C) HIFs promote transcription in HBV de novo infection. HepG2 cells supporting episomal copies of HBV DNA (HepG2-HBV-pEpi) were treated with FG-4592 (30μM) for 72 h and pC/pgRNA, host genes CAIX and NDRG1, were quantified by qPCR. HBV-infected (MOI 200) HepG2-NTCP cells were treated with FG-4592 (30 μM) and pC/pgRNA, CAIX and NDRG1 expression was quantified 72 h post infection. (D) HIF KD studies. siRNAs targeting HIF-1α or HIF-2α were delivered into HepG2-NTCP cells that were infected with HBV (MOI 200) and treated with FG-4592 (30 μM) for 72 h. siRNA knockdown was confirmed by qPCR of host genes, CAIX (HIF-1α dependent) and VEGFA (HIF-1α and -2α co-dependent) and cell lysates probed for HIF-1α, HIF-2α and CAIX expression. HBV pC/pgRNA levels quantified by PCR. Data is presented as fold induction relative to the 18% oxygen control for each condition. Statistical analysis was performed using a one-way ANOVA (∗p ≤0.05, ∗∗p ≤0.01). BCP, basal core promoter; HIFs, hypoxia inducible factors; HRE, hypoxic responsive element; KD, knockdown; MOI, multiplicity of infection. (This figure appears in color on the web.)

To extend these observations we used 2 model systems to investigate a role for HIFs in regulating viral transcriptional activity: HepG2 cells with an episomal replicating HBV (HepG2-pEpi-HBV) and HepG2-NTCP cells that support de novo HBV infection. The BCP drives transcription of pC and pgRNA from 2 start sites that are only 70 base pairs apart and since our PCR cannot discriminate between these viral-encoded RNAs, we label transcripts as pC/pgRNA to represent the sum of both RNAs. FG-4592 induced a significant increase in pC/pgRNA and HIF-regulated genes carbonic anhydrase IX (CAIX) and N-myc downstream regulated 1 (NDRG1) in both models (Fig. 1C). Since HIF-1α and HIF-2α can show non-overlapping and opposing functions we were interested to investigate their individual roles in regulating viral transcription. We transiently silenced HIF-1α or HIF-2α, individually or together, in HepG2-NTCP cells prior to infecting with HBV in the presence or absence of FG-4592. We demonstrate effective silencing by western blotting for HIF expression and by quantifying CAIX or vascular endothelial growth factor (VEGF) mRNA levels (HIF-1α and HIF-1α/HIF-2α regulated host genes, respectively) (Fig. 1D). The FG-4592-dependent increase in pC/pgRNA in the infected cells was ablated by the combined HIF-1α and HIF-2α silencing, demonstrating a role for both HIFα isoforms in regulating promoter activity.

As our cell-based studies predict an association between HIF-transcriptional activity and viral RNAs in the infected liver, we quantified hypoxic gene transcripts and HBV RNA in liver biopsies from patients diagnosed with chronic hepatitis (n = 12) or chronic infection (n = 7) following the 2017 EASL guidelines (Fig. 2A). We recently reported increased expression of hypoxia regulated gene transcripts in chronic hepatitis B (CHB), enabling us to identify hepatic HIF-target genes for inclusion in a customized nCounter (Nanostring) array together with a probe to measure HBV RNA. The HBV probe detected a 400-fold range in viral RNA counts amongst the CHB samples (Fig. 2B), enabling us to stratify patients into 2 groups (nCounter score low <350 and high >350). We found a significant positive association between HIF-gene expression and HBV RNA (Fig. 2C), consistent with our in vitro studies suggesting a role for HIFs in positively regulating HBV transcription.

Fig. 2

Hypoxic gene expression associates with HBV RNA levels in chronic infection.

(A) Clinical characteristics of CHB cohort. (B) HBV RNA in the CHB biopsies was quantified by nCounter and classified into high (Normalized count >350: n = 9 and red symbols) or low (Normalized count <350: n = 10 and blue symbols) based on their deviation from the median, where the dashed line indicates the assay cut-off defined by analyzing liver RNA samples from uninfected control patients (n = 8). (C) HBV RNA positively associates with the expression of 15 HIF-regulated genes (BNIP3, BNIP3L, CA9, EPO, FAM115C, FGF11, GLUT1, IGLON5, LCN15, LOX, NDRG1, PFKFB4, SPAG4, TNS1, VEGFA). Statistical testing was performed using Wilcoxon match-pairs signed rank test (p = 0.00012). CHB, chronic hepatitis B; HIF, hypoxia inducible factor. (This figure appears in color on the web.)

HIFs regulate viral transcripts in HBV transgenic mice.

HBV can only infect humans and primates and no immune competent animal models are available that support natural HBV infection. The HBV transgenic mouse model (HBVtg) of chronic infection transcribes viral RNAs from an integrated 1.3 overlength genome and, despite the lack of cccDNA, has been used to study host pathways regulating viral transcription (reviewed in). Since both HIFα subunits dimerize with HIF-1β to regulate gene transcription we transiently silenced HIF-1β by siRNA injection and quantified viral parameters in the liver and periphery. To assess the efficacy of the silencing we measured mRNA levels of HIF-1β, the HIF-regulated host genes VEGF and PHD2, the key oxygen sensor that regulates HIF expression. We observed a coordinated reduction of all 3 gene transcripts in 6 of 7 male mice receiving siRNAs targeting HIF-1β and a negligible effect of the siRNA-HIF-1β in 5 female mice (Fig. 3A, Fig. S1). We noted reduced pC/pgRNA and total viral RNAs in 4 of 6 siRNA-HIF1β treated male mice compared to siControl treated animals (Fig. 3B), demonstrating a role for HIFs in regulating HBV transcription. Alanine aminotransferase values were comparable in all mice suggesting a negligible cytotoxic effect of the siRNA-HIF-1β treatment. Furthermore, visual inspection of H&E-stained liver biopsies from siCtrl or siRNA-HIF-1β treated animals showed no evidence of hepatocyte proliferation. We observed a reduced frequency of HBcAg-expressing hepatocytes, peripheral HBeAg and HBsAg levels in the siRNA-HIF-1β treated male mice (Fig. 3B,C). Core antigen (HBcAg) expression showed a peri-central staining pattern (Fig. 3C), consistent with this area of the liver representing a low oxygen environment. These data highlight a role for HIFs in regulating HBV transcription in this transgenic mouse model.

Fig. 3

HIFs regulate viral transcription in HBV transgenic mice.

(A) siRNAs targeting HIF-1β were delivered into male and female HBV1.3 transgenic mice and animals were culled after 1 week. RNA was isolated from liver biopsies and transcript levels of HIF-1β, Vegfa and Phd2, together with pC/pgRNA and total viral RNA quantified. The amount of each transcript was normalized (Z-score) and presented as a heat map (suppression = red, enhancement=blue). (B) HBcAg expression in the liver was assessed by measuring the number of antigen expressing cells/mm2. Peripheral HBeAg, HBsAg and ALT levels were quantified in siCTRL and HIF-1β silenced mice and tabulated per individual mouse. (C) HBcAg staining from fixed murine liver sections. 20x images are shown with scale bars representing 100 μm. ALT, alanine aminotransferase; HIFs, hypoxia inducible factors; pC/pgRNA, pre-core/pre-genomic RNA. (This figure appears in color on the web.)

Hypoxia inducible factors bind HBV DNA

To assess whether HIFs bind HBV DNA we isolated chromatin from HepG2-pEpi-HBV cells and quantified HIF-1β binding to episomal genomes by chromatin immunoprecipitation and quantitative PCR (ChIP-qPCR). Under normoxic conditions we failed to precipitate HIF-1β with any target gene (Fig. 4A), whereas HNF4α bound viral DNA and the promoter of Apolipoprotein B (APOB), (Fig. 4A). The retinoic acid receptor-related orphan receptor-α (RORα) is an important regulator of circadian rhythm and hepatic lipid metabolism and, to the best of our knowledge, does not interact with the HBV genome. As a control for these ChIP experiments, we showed that RORα bound its circadian target gene BMAL1 but failed to bind HBV DNA (Fig. 4A). In vitro studies to stabilize HIFs are routinely performed under 1% oxygen and culturing cells under these conditions or treatment with FG-4592 induced HIF-1β binding to HBV DNA and the host target gene CAIX (Fig. 4B). HIF-1β showed no binding to the γ-Globin promoter, a known housekeeping gene lacking an HRE motif (Fig. 4B). To assess whether HIFs bind HBV DNA in vivo we studied infected chimeric human liver FNRG mice. HBV-infected mice were sacrificed at 4 weeks post-infection and viral replication was confirmed by measuring hepatic cccDNA along with peripheral HBV DNA, HBeAg and HBsAg (Fig. 4C). ChIP-qPCR using primers spanning the BCP showed enriched HIF-1β, HIF-1α and HIF-2α binding to viral genomes relative to an irrelevant IgG (Fig. 4C). In agreement with our earlier in vitro experiments HNF4α also bound viral DNA whereas there was no evidence for an interaction with RORα (Fig. 4C). In summary, these experiments provide clear evidence for HIF-1β binding HBV DNA following FG-4592 or low oxygen treatments. Furthermore, HIF-1α and HIF-2α complexed with HBV DNA in infected human hepatocytes isolated from liver chimeric mice, without any specific treatments to stabilize HIFs, showing the presence of these complexes in vivo.

Fig. 4

HIFs bind HBV DNA.

(A) Limited evidence for HIF-1β binding HBV under normoxic conditions. Chromatin extracts from HepG2-HBV-pEpi cells cultured under 18% O2 were immunoprecipitated with anti-HIF-1β, anti-HNF4α, anti-RORα or an irrelevant IgG control. PCR for HBV cccDNA and host genes CAIX, APOB and BMAL1 was performed and data presented relative to % of input of the IgG. Data are the mean + SD of 2 independent experiments and statistical analysis performed with Mann-Whitney U test (∗p ≤0.05). (B) Hypoxia and FG-4592 potentiates HIF binding to HBV DNA. Chromatin extracts from HepG2-HBV-pEpi cells cultured at 18% or 1% O2 or treated with FG-4592 (30 μM) for 24 h were immunoprecipitated with anti-HIF1β or irrelevant IgG control. PCR for HBV cccDNA and HIF-regulated gene CAIX and the non-HIF housekeeper gene γ-Globin was performed and the data presented relative to untreated (18% O2) samples. Data are representative of mean + SD of 2 independent experiments and statistical analysis was performed with a Mann-Whitney U test (∗p ≤0.05). (C) HIFs bind HBV DNA in infected liver chimeric mice. Chromatin from HBV-infected mouse liver tissue was immunoprecipitated with antibodies specific for HIF-1β, HIF-1α, HIF-2α, HNF4α, RORα along with an irrelevant IgG control. Following immunoprecipitation, the DNA was analyzed by qPCR with primers spanning the binding sites of the respective transcription factors on the HBV genome. % of input was calculated and data presented relative to the IgG control and shown as the mean + SD for 2 mice. cccDNA, covalently closed circular DNA; HIFs, hypoxia inducible factors; qPCR, quantitative PCR.

A low oxygen environment activates HBV BCP and associated transcripts

To investigate the functional consequences of a low oxygen environment on HBV transcription activity we cultured HepG2-NTCP cells under 1% oxygen and measured promoter activity and pC/pgRNA levels in de novo infected cells. Under low oxygen conditions HepG2-NTCP cells express both HIFα isoforms and show a time-dependent increase in promoter activity (Fig. 5A), increased pC/pgRNA levels and HBcAg expression (Fig. 5B). We observed a significant increase in the level of secreted HBV DNA under low oxygen conditions using 2 independent model systems (Ad-HBV transduced HepG2-NTCP and HepG2.2.15 cells) (Fig. S2), consistent with reports showing that a modest increase in HBcAg promotes cytoplasmic HBV replication. Low oxygen had no impact on cccDNA levels, consistent with our interpretation that HIFs positively regulate BCP activity. To further explore this conclusion, we infected hypoxic HepG2-NTCP cells with a transcriptionally active but replication deficient HBV reporter virus (rHBV-Gluc) and showed a low oxygen-dependent increase in pC/pgRNA (Fig. 5C). This virus encodes Gaussia luciferase under the control of an exogenous transthyretin promoter within the S open reading frame. We observed a modest reduction in reporter activity in the hypoxic infected cells, demonstrating the promoter-dependency of low oxygen in regulating HBV transcription. HIFs are rapidly degraded under normoxic conditions to facilitate dynamic cellular transcriptional responses to fluctuating oxygen levels. We were interested to assess the impact of re-oxygenation on HBV transcription. HBV-infected HepG2-NTCP cells were maintained in 1% oxygen for 72 h, moved to 18% oxygen and sampled at 0, 6, 24 and 48 h post-oxygenation. We quantified cccDNA and transcript levels of pC/pgRNA, CAIX and NDRG1. cccDNA levels showed no significant change over the period of re-oxygenation, however, pC/pgRNA levels declined rapidly, with a 50% reduction by 11 h compared to 8 h and 2 h for CAIX and NDRG1 RNAs, respectively (Fig. 5D). These data clearly show that low oxygen regulates cccDNA transcriptional activity, highlighting the dynamic nature of this phenotype.

Fig. 5

Low oxygen promotes HBV transcription.

(A) Low oxygen activates HBV BCP. HepG2-NTCP cells expressing HBV EnhI/II and BCP construct driving luciferase were cultured at 18% or 1% oxygen. Viral promoter activity was assessed by measuring luciferase activity after 24, 48 and 72 h. Statistical analysis was performed using a 2-way ANOVA (∗∗p ≤0.01, ∗∗∗∗p ≤0.0001). (B) Low oxygen promotes de novo HBV infection associated transcription. HepG2-NTCP cells were infected with HBV at an MOI of 200 and incubated at 18% or 1% oxygen for 72 h. cccDNA, pC/pgRNA and CAIX mRNA levels were quantified by qPCR and HBcAg assessed by western blotting. (C) HepG2-NTCP cells were infected with a modified HBV encoding Gaussia luciferase and incubated under 18% or 1% oxygen for 72 h. pC/pgRNA levels were quantified and Gaussia luciferase activity measured. (D) The dynamic nature of low oxygen-regulation of HBV pC/pgRNA. HBV-infected HepG2-NTCP cells cultured at 1% oxygen for 72 h were perfused with 18% oxygen and samples collected at the indicated times. HBV cccDNA and pC/pgRNA along with host genes CAIX and NDRG1 were quantified and used to estimate mRNA half-life. cccDNA transcriptional activity at 1% oxygen and during the period of re-oxygenation is shown. All data presented as the mean ± SD derived from at least 3 independent experiments consisting of 3 replicates per condition. Statistical analysis was performed using a Mann-Whitney U test (∗p ≤0.05, ∗∗p ≤0.01, ∗∗∗∗p ≤0.0001). cccDNA, covalently closed circular DNA; HIFs, hypoxia inducible factors; pC/pgRNA, pre-core/pre-genomic RNA; qPCR, quantitative PCR.

In the adult liver hepatocytes are non-proliferating and the majority of in vitro experiments studying HBV replication use DMSO to arrest cell proliferation. We noted that culturing HepG2-NTCP cells under 1% oxygen arrests cells (Fig. S3), providing a physiological method to limit cell proliferation. DMSO can induce wide-ranging effects on host gene expression and we evaluated the effect of low oxygen on pC/pgRNA levels in DMSO-arrested HepG2-NTCP cells, comparing a naïve de novo infection to a pre-established infection. In both cases we noted a blunted response of low oxygen to regulate pC/pgRNA in the DMSO-arrested cells (Fig. S4A,B). Screening a panel of hypoxia genes by PCR-array showed a reduced activation of both HIF-1α and HIF-2α regulated genes in DMSO-treated HepG2 cells cultured in 1% oxygen (Fig. S5).

Primary human hepatocytes (PHHs) are considered the gold standard for studying HBV replication. However, PHHs can de-differentiate and lose hepatocyte-specific function in vitro. One approach to limit their de-differentiation is to culture them in DMSO-containing media. Since DMSO blunts HepG2 cellular response to hypoxia, we cultured PHHs and HepG2-NTCP cells in 5%, 3% or 1% oxygen for 48 h and showed negligible expression of HIF-1α and greater levels of HIF-2α in PHHs compared to HepG2-NTCP (Fig. S5). We observed a more dominant pattern of HIF-2α regulated genes in PHHs compared to HepG2-NTCP cells (Fig. S5C). To understand whether HIF-2α transcriptional activity dominates in the HBV-infected liver we investigated the relative contribution of HIF-1α or HIF-2α activity in CHB biopsy samples. 38 of the hypoxia genes showing increased expression in CHB were defined as HIF-1α (n = 26), HIF-2α (n = 11) or co-regulated (n = 1) genes, showing clear evidence for both HIF-1α and HIF-2α transcriptional activity in vivo. These data show the complexities of studying viral and host transcriptional responses to low oxygen in PHHs ex vivo.

Oxygen-regulation of episomal and integrated HBV DNA

In the human genome HIFs bind motifs that can regulate distant gene promoters. We previously reported a qPCR technique to quantify the relative abundance of HBV RNAs and used this method to assess the effect of low oxygen on the pattern of viral RNAs in infected HepG2-NTCP cells. pC/pgRNA was the major viral-encoded RNA in infected HepG2-NTCP cells, consistent with a recent report using 5’RACE to study HBV transcripts. Low oxygen increased pC/pgRNA levels compared to preS1/S or HBx RNAs in HepG2-NTCP cells (Fig. 6A). Furthermore, following oxygen reperfusion of hypoxic infected cells, preS1/S RNA levels showed no significant change, suggesting a minimal role for low oxygen in regulating their genesis or estimated half-life (Fig. 6B). Recent reports highlight a role for HBV integrants in driving HBsAg expression., Since integration of the linear viral DNA generated during viral replication may separate the HREs from downstream transcription initiation sites, we hypothesize that integrated copies of HBV will be oxygen-insensitive. We used 3 well-characterized hepatoma lines: PLC/PRF5, Hep3B and Huh-1 that harbor transcriptionally active integrated genomes., Low oxygen had a minimal effect on preS1/S2 RNA levels (Fig. 6C), with a modest reduction in transcripts in PLC/PRF5 cells, and no significant changes in HBsAg expression. We confirmed the integrant lines were responsive to low oxygen, showing an approximate 10-fold induction in CAIX RNA levels. These experiments show the differential effect of HIFs in the regulation of transcription of episomal cccDNA and HBV integrants.

Fig. 6

Oxygen-dependent differential regulation of episomal or integrated HBV DNA.

(A) Quantification of the major HBV transcripts under low oxygen. Infected HepG2-NTCP were cultured at 18% or 1% oxygen for 72 h. pC/pgRNA, PreS1/S and HBx RNA transcripts were quantified by PCR using specific primers with copy number inferred from a standard curve, adjusted to the host housekeeper gene, β2M. (B) Effect of oxygen reperfusion on HBV preS1/S RNA levels. HBV-infected HepG2-NTCP cells cultured under 1% oxygen for 72 h were perfused with 18% oxygen and samples collected at the indicated times. HBV pgRNA and preS1/S RNA levels were quantified by qPCR. Data presented are the mean ± SD derived from 3 independent experiments consisting of 3 replicates per condition. (C) Hepatoma lines bearing viral integrants were treated as in (A) and preS1/S RNAs quantified and data presented are the mean ± SD derived from 3 independent experiments consisting of 3 replicates per condition. Statistical analysis was performed using a Mann-Whitney U test (∗∗∗p ≤0.001). pC/pgRNA, pre-core/pre-genomic RNA; qPCR, quantitative PCR.

Discussion

Oxygen sensing is a fundamental cellular process that shapes the liver transcriptome. We identify a role for HIFs in activating the HBV core promoter and increasing pC/pgRNA, HBcAg expression and secreted viral DNA. Our ChIP experiments show HIF binding to HBV DNA in hypoxic infected cells and from HBV-infected human liver chimeric mice in vivo, supporting a direct role for HIFs in the regulation of HBV transcription. Furthermore, we translate these observations into CHB and show a positive association between hypoxic gene expression and viral RNAs, consistent with a role for HIFs in regulating HBV in vivo. HBV encodes 2 HREs in its compact genome and we show their evolutionary conservation among primate members of the hepadnaviridae.

Silencing HIF-1β in male HBVtg mice reduced pC/pgRNA and total RNA levels, the frequency of HBcAg-expressing hepatocytes, and peripheral HBeAg and HBsAg levels, revealing a role for HIFs in regulating viral transcription in this model. Analyzing published single cell-RNA sequences of mouse liver for transcriptional activators and repressors reported to regulate HBV replication (reviewed in) showed no zonal patterns of expression, consistent with a role for HIFs in regulating viral transcription in this mouse model. We noticed a discrepancy between the male and female mice in our study, where HIF-1β silencing was ineffective in females, consistent with reports showing a greater response in male rats to chronic intermittent hypoxia compared to females., Furthermore, ovariectomy rendered the 2 sexes equally responsive to hypoxia. Estrogen signaling was reported to inhibit both HIF-1α transcriptional activity and HBV transcription,, providing a potential explanation for the sex-dependent differences noted in our model. A recent in-depth proteomic study investigating the molecular basis of sex difference in zebrafish responses to hypoxia preconditioning identified a complex network of signaling pathways. Hypoxic gene expression and HBV RNA levels in our CHB samples were independent of sex, however, this may reflect the post-menopausal age of females in our small cohort and is worthy of further investigation.

Hypoxia is known to induce a broad range of cellular changes beyond those regulated by HIFs. Importantly, we observed higher pC/pgRNA in de novo HBV-infected cells cultured under low oxygen compared to those treated with the PHD inhibitor FG-4592 (Roxadustat), suggesting a potential role for other members of the oxygenase family in the regulation of HBV replication. To investigate this possibility, we completed a transcriptional and proteomic analysis of hypoxic human hepatoma cells (Fig. S6) and observed negligible changes in the expression of host activators or repressors previously reported to regulate HBV replication. A recent study using chemical mimetics (CoCl2 and dimethyloxalylglycine) to stabilize HIFs reported increased DNASE1 levels and a reduction in HBV genome copies, however, our hypoxic studies showed no impact on DNASE1 expression that may reflect differences in the model systems used.

Integration of the HBV genome is associated with HCC development, and the introduction of HREs could regulate downstream targets such as oncogenes that could be activated in the hypoxic tumor environment. HBV integrants are thought to be the major source of HBsAg expression in chronic disease that have been associated with exhausted antiviral T cell responses., Our observation that integrant encoded preS1/S transcripts are oxygen-insensitive is relevant if one considers the hypoxic nature of the HCC environment, and how this could influence HBsAg expression, a biomarker for monitoring patients’ response to new therapies.

In summary, we demonstrate a new role for HIFs to positively regulate HBV transcription. A hypoxic environment was also reported to potentiate hepatitis C virus replication via HIF-dependent, and independent pathways. Hypoxia can have variable effects on virus replication, most likely reflecting the differing oxygen tensions at the sites of replication. HBV, along with other members of the hepadnaviridae family, may have evolved to replicate and exploit the HIF-signaling pathway to persist in the low oxygen environment of the liver.

Our observations raise questions as to how cellular hypoxia translates to humans, both in terms of HBV replication, hepatic immunity and response to new therapies. CHB reflects a dynamic interaction between the virus and host immune system cells, where active hepatitis associates with increased cccDNA transcription. Our data showing a role for HIFs in promoting HBV transcription provide an explanation for how inflammatory responses may potentiate HBV replication. The recent licensing of HIF-prolyl hydroxylase inhibitors as erythropoiesis-stimulating agents for the treatment of anemic patients with chronic kidney disease in China and Japan could impact HBV replication and is worthy of further study. Finally, pharmacological agents that inhibit HIF signaling-pathways may be considered for further exploration in clinical trials. Understanding the biology of HBV under more physiological conditions will reveal novel therapeutic targets and provide an environment for innovative antiviral drug screening.

Abbreviations

BCP, basal core promoter; CAIX, carbonic anhydrase IX; cccDNA, covalently closed circular DNA; CHB, chronic hepatitis B; ChIP, chromatin immunoprecipitation; HCC, hepatocellular carcinoma; HIFs, hypoxia inducible factors; HRE, hypoxic responsive element; NDRG1, N-myc downstream regulated 1; NTCP, sodium taurocholate co-transporting polypeptide; pC, pre-core; pgRNA, PHD, prolyl-hydroxylase domain; qPCR, quantitative PCR; RORα, retinoic acid receptor-related orphan receptor-α; VEGF, vascular endothelial growth factor.

Financial support

JAM laboratory is funded by IA 200838/Z/16/Z, project grant MR/R022011/1 and Chinese Academy of Medical Sciences (CAMS) Innovation Fund for Medical Science (CIFMS), China (grant number: 2018-I2M-2-002). UP laboratory is funded by the German research foundation via the collaborative research center TRR179 (project TP14; UP). TM received a Gerok stipend by TRR179. UP and JAM collaboration was supported by funding by the Institute for Advanced Study with the support of the Technical University of Munich via the German Excellence Initiative and EU 7th Framework Program under grant agreement number 291,763. JW received a stipend by the . MW lab is funded by a Wellcome Trust Senior Clinical Research Fellowship 108070/Z/15/Z. MP was (partially) funded by the Università del Piemonte Orientale, FAR-2017. TFB was supported by ARC, Paris (TheraHCC2 IHUARC201901299), , LabEx HepSys and Inserm Plan Cancer. DRM laboratory is funded by the (NIHR-RP-2016-06-004), the Deanship of Scientific Research, King Abdulaziz University, for Saudi Arabia. JAM, UP, TFB and MH were supported by the (EU H2020-667273-HEPCAR).

Authors’ contributions

PW designed experiments, analyzed and curated data and co-wrote the manuscript; PJL designed experiments, analyzed and curated data and edited the manuscript; JMH designed experiments, analyzed and curated data and edited the manuscript; AM conducted experiments, analyzed and curated data; TM conducted experiments; XZ conducted experiments and analyzed data; HB conducted experiments and analyzed data; RM provided clinical samples; NRF conducted experiments; JW provided reagents; LM conducted experiments; VD conducted experiments; TR conducted experiments; LN analyzed data; MPW analyzed data; MP provided clinical material and advice; MH provided reagents and advice; TFB provided reagents and advice; EH provided advice and edited manuscript; DRM provided advice and reagents; UP provided advice and reagents; PB analyzed data and co-wrote the manuscript and JAM designed the study and wrote the manuscript.

Data availability statement

The authors declare that all data supporting the findings of this study are available in the article along with supplementary Information file.

Conflict of interest

None of the authors have any conflict of interest.

Please refer to the accompanying ICMJE disclosure forms for further details.