Article Text
Abstract
Objectives A major goal of curative hepatitis B virus (HBV) treatments is the reduction or inactivation of intrahepatic viral covalently closed circular DNA (cccDNA). Hence, precise cccDNA quantification is essential in preclinical and clinical studies. Southern blot (SB) permits cccDNA visualisation but lacks sensitivity and is very laborious. Quantitative PCR (qPCR) has no such limitations but inaccurate quantification due to codetection of viral replicative intermediates (RI) can occur. The use of different samples, preservation conditions, DNA extraction, nuclease digestion methods and qPCR strategies has hindered standardisation. Within the ICE-HBV consortium, available and novel protocols for cccDNA isolation and qPCR quantification in liver tissues and cell cultures were compared in six laboratories to develop evidence-based guidance for best practices.
Design Reference material (HBV-infected humanised mouse livers and HepG2-NTCP cells) was exchanged for cross-validation. Each group compared different DNA extraction methods (Hirt extraction, total DNA extraction with or without proteinase K treatment (+PK/−PK)) and nuclease digestion protocols (plasmid-safe ATP-dependent DNase (PSD), T5 exonuclease, exonucleases I/III). Samples were analysed by qPCR and SB.
Results Hirt and −PK extraction reduced coexisting RI forms. However, both cccDNA and the protein-free relaxed circular HBV DNA (pf-rcDNA) form were detected by qPCR. T5 and Exo I/III nucleases efficiently removed all RI forms. In contrast, PSD did not digest pf-rcDNA, but was less prone to induce cccDNA overdigestion. In stabilised tissues (eg, Allprotect), nucleases had detrimental effects on cccDNA.
Conclusions We present here a comprehensive evidence-based guidance for optimising, controlling and validating cccDNA measurements using available qPCR assays.
- CHRONIC HEPATITIS
- HEPATITIS B
- LIVER
- REAL TIME PCR
- LIVER BIOPSY
Data availability statement
All data relevant to the study are included in the article or uploaded as online supplemental information.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
Persistence of the hepatitis B virus (HBV) genome, the covalently closed circular DNA (cccDNA), in the liver is a major obstacle for developing HBV cure strategies and its accurate quantification is essential in preclinical and clinical studies. However, cccDNA quantification by PCR lacks standardisation and is technically challenging due to the heterogeneous population of HBV DNA molecules present in infected cells.
WHAT THIS STUDY ADDS
Distinct cccDNA extraction and quantification protocols were cross-validated using infected tissues and cell cultures. We provide experimental evidence that coexistence of certain HBV DNA forms and potential cccDNA overdigestion due to sample preservation conditions and nuclease treatment still challenges PCR-based quantification of cccDNA. The results support recommendations for optimising, controlling and validating cccDNA measurements according to the sample type.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
This work presents an evidence-based guidance for best practices to quantify cccDNA by qPCR. Information provided will assist the HBV cure research programmes aiming at assessing the impact of therapies on the HBV reservoir in preclinical studies and clinical trials.
Introduction
The hepatitis B virus (HBV) is the cause of chronic hepatitis B (CHB) and despite the availability of effective prophylactic vaccines and treatments that efficiently suppress viral replication, HBV puts 296 million carriers at risk of developing liver cirrhosis and hepatocellular carcinoma.1 Thus, there is a strong need to develop therapies with the potential to cure CHB. The unique HBV replication strategy, however, makes viral eradication challenging.2 Knowledge and precise monitoring of the different viral components present in infected cells and in the circulation is paramount for the development and assessment of novel therapeutic strategies.
HBV is an enveloped virus harbouring a small (3.2 kb) partially double-stranded (ds) relaxed circular DNA (rcDNA) genome which specifically infects human hepatocytes. On infection, the rcDNA is repaired by cellular enzymes to generate a double-stranded episomal DNA genome, the covalently closed circular DNA (cccDNA)3 4 (online supplemental figure 1). By associating with histone and non-histone proteins, the cccDNA forms a stable minichromosome within the cell nuclei that serves as transcriptional template for all viral transcripts. Nucleoside analogues (NAs) and interferon alpha are the only approved treatments for CHB but are mainly suppressive as they do not clear cccDNA from the liver. Because of the central role of the cccDNA in HBV persistence, the elimination of this molecule constitutes the ultimate goal of curative treatments.5 Such a goal may be very difficult to achieve and current investigational approaches rather aim at achieving a functional cure, which is defined as the loss of HBsAg and undetectable HBV DNA in serum.2 Examples of currently explored curative strategies include, among others, the use of capsid assembly modulators and entry inhibitors, which can also prevent cccDNA formation and intrahepatic amplification of the cccDNA pool,6 degradation of HBV transcripts using siRNA technologies, immune-mediated clearance and silencing of cccDNA molecules, or its inactivation by gene editing.7 Assessment of the impact of novel therapies on cccDNA load and activity is therefore of utmost importance. Moreover, our understanding of cccDNA biology in preclinical models and patient biopsies is still limited. Because of the paucity of biopsy material, non-invasive biomarkers (eg, HBcrAg, HBV RNA), are increasingly studied as surrogate cccDNA markers.8 However, for the primary validation of these biomarkers, proper cccDNA quantification is needed.
Supplemental material
A major limitation in HBV research and its translation into the clinic is the lack of standardised PCR-based methods allowing specific cccDNA quantification in HBV-infected samples. The main challenge is the reliable detection of cccDNA by quantitative PCR in the presence of high excess of rcDNA coming from virus input or coexisting replicative intermediates (RI), which are identical in sequence to the cccDNA (online supplemental figure 1). To gain specificity, PCR approaches use primers spanning the cohesive ends region of the rcDNA to preferentially detect the cccDNA (online supplemental figure 2A). Unfortunately, this approach is not completely specific since it still detects certain amounts of the highly abundant RI leading to an overestimation of cccDNA amounts.9 10 Only when HBV replication is low (eg, during NA treatment), cccDNA quantification by PCR appears precise.11 12 Therefore, cccDNA measurements in samples with high viral productivity or between samples with and without treatment remain challenging and reducing the levels of RI is mandatory.
To reduce the levels of HBV RI prior to PCR, various nucleases and DNA extraction methods have been proposed, although their efficacy and specificity in engaging distinct HBV DNA forms in different experimental settings vary substantially. The conventional Hirt method employs sodium dodecyl sulfate lysis and high salt precipitation of high-molecular-weight cellular chromatin and protein-bound DNA13 thereby enriching the recovery of protein-free (pf) DNA molecules and facilitating the removal of rcDNA, which is covalently linked to the viral polymerase. Southern blot (SB) after Hirt DNA extraction allows visualisation of different HBV DNA molecules including cccDNA.14 However, this methodology necessitates large amounts of nucleic acids and therefore lacks sensitivity; it is not suitable for liver biopsy analysis or high-throughput assays and does not allow precise quantification.
The nucleases commonly used display slightly different substrate specificities but share the common characteristic of sparing closed circular dsDNA from digestion.11 15–18 Apart from employing different nucleases, the use of different DNA extraction procedures, qPCR conditions and normalisation increases variability and hinders reliable comparative analyses among laboratories. We, therefore, attempted to harmonise cccDNA quantification processes by comparing the most used and newer protocols through cross-validation experiments. Based on the results obtained, we propose recommendations of best practice for controlling and validating cccDNA measurements also according to the sample type.
Material and methods
For generation, infection, treatment and viral characterisation of human liver chimeric mouse, as well as infection of HepG2-NTCP cell culture samples, nuclease digestion and qPCR measurements, see online supplemental material.
Preparation of liver tissue and cell samples
Liver specimens removed at the time of sacrifice were snap-frozen in 2-methylbutane and stored at −80°C until further use. Every lab received three frozen liver pieces (approximately 13 mg each) from one highly infected, untreated uPA/SCID/beige/IL2RG-/- (USG) mouse. The tissue pieces were homogenised separately in 300 µL 10 mM Tris-HCl (pH 7.5)/10 mM EDTA (pH 8.0) buffer using disposable homogenisers (Biomasher II, DWK Life Sciences, Wertheim, Germany), then pooled and distributed equally among the three DNA extractions. HBV-infected HepG2-NTCP cells were pooled before being divided and shipped to the four participating labs. Every lab received a frozen cell pellet of approximately 1×107 cells, which on arrival was dissolved in 900 µL 10 mM Tris-HCl/10 mM EDTA buffer and distributed equally for the three DNA extractions. Since all samples were pooled and divided equally across the DNA extractions, the results could be compared with each other directly without the need for additional normalisation to cellular DNA.
DNA extraction procedures
Total cellular DNA was isolated with the MasterPure Complete DNA and RNA Purification Kit (Epicentre, available through Lucigen, Middleton, Wisconsin, USA) as recommended by the manufacturer.12 19 The total DNA extraction with proteinase K treatment (+PK) and the one without PK (−PK) were performed identically, except for the PK digestion step, which was omitted in the −PK extraction. Briefly, the samples were split up in several microtubes (according to the manufacturer′s instructions, 7×105 cells or 4.4 mg tissue per 300 µL lysis buffer) and 10 mM Tris-HCl/10 mM EDTA buffer was added for a final volume of 300 µL each. One volume of double concentrated TCL buffer (Lucigen) was added to lyse the cells. Subsequently, RNA digestion was performed with 2 µL RNase A (Lucigen) for 30 min at 37°C. For the+PK extraction, 2 µL of PK (Lucigen) was added and incubated for 1 hour at 56°C. After the addition of 300 µL of MPC buffer (Lucigen), the samples were mixed and incubated on ice for 5 min, followed by 10 min high speed centrifugation to pellet the protein. DNA was recovered from the supernatant through isopropanol precipitation, dissolved in 10 mM Tris-HCl/1 mM EDTA buffer and the DNA from one extraction and sample was pooled again. The DNA content was determined by fluorometry (Qubit, Invitrogen) or spectrophotometry (NanoDrop, Thermo Fisher Scientific, Waltham, Massachusetts, USA). A detailed version of this protocol can be found on the ICE-HBV protocol database (https://ice-hbv.org/protocol/a-modified-kit-based-hbv-protein-free-dna-extraction-from-liver-tissues-and-cell-cultures-for-hbv-cccdna-southern-blot -and-qpcr/). For a comparison of different lysis buffers, the protocol was repeated also using RLT buffer from the RNeasy Mini kit (Qiagen, Hilden, Germany). Some mouse liver tissues were immersed in Allprotect Tissue Reagent (Qiagen) as recommended by the manufacturer and stored at −20°C for several months.
One-third of the liver homogenate and cell suspension was subjected to DNA extraction according to Hirt.13 Briefly, the samples were further diluted in 10 mM Tris-HCl/10 mM EDTA buffer to a final volume of 3 mL. Cell lysis was achieved by adding SDS (final concentration 0.6%) and 30 min of incubation at room temperature with slow agitation. After adding KCl (final concentration 0.5M), the lysates were incubated at 4°C overnight, then centrifuged for 30 min at 4°C and 14 500 g. DNA from the supernatant was extracted twice with phenol and once with phenol/chloroform/isoamyl alcohol (25:24:1), followed by an overnight ethanol precipitation. The DNA was pelleted by centrifugation, washed in 70% ethanol and dissolved in 10 mM Tris-HCl/1 mM EDTA buffer.
SB analyses were performed as previously described.19 20 A detailed version of this protocol can be found on the ICE-HBV protocol database (https://ice-hbv.org/protocol/a-sensitive-and-rapid-southern-blot -assay-based-on-branched-dna-technology-for-the-detection-of-hbv-dna-in-cell-culture-and-liver-tissue-samples/) and as online supplemental file.
Results
Comparison of cccDNA extraction protocols in HBV-infected liver tissues
Currently available and novel protocols for cccDNA extraction and quantification were collected and cross-validated according to the scheme depicted in figure 1A. Frozen liver tissues from HBV-infected humanised mice were used as reference material (online supplemental figure 3). We compared total DNA extraction methods that either included (+PK) or excluded (−PK) the PK digestion step, as well as the classical Hirt DNA extraction (Hirt) that does not include PK digestion.13 Omission of the PK should facilitate the removal of covalently protein-bound HBV DNA and was investigated here as potential alternative to Hirt extraction. A previously published cccDNA-selective PCR strategy21 based on selective primers and a Taqman probe detecting all HBV genotypes except G (online supplemental table 1), and not cross-reacting with human or murine genomic sequences, or HBV DNA integrations,12 was chosen for cross-validation using improved PCR conditions (for details see online supplemental figure 2A–C). The specificity for cccDNA detection was assessed on serial dilutions of serum-derived rcDNA and on mixtures of an HBV plasmid (as cccDNA surrogate) and increasing amounts of serum-derived rcDNA (online supplemental figure 2D, E). These analyses indicated that serum rcDNA can be amplified in a false-positive manner despite using cccDNA-selective PCR primers and probes when rcDNA levels exceed cccDNA levels by more than 250-fold. Under such circumstances, measured cccDNA levels increase and a precise cccDNA quantification is not possible without including strategies reducing rcDNA levels.
Supplemental material
The effect of different DNA extraction methods on cccDNA quantification by qPCR and SB in HBV-infected USG mouse liver tissue. (A) Schematic presentation of the experimental design used for the cross-validation. (B) SB analysis on non-digested DNA extracts using HBV DNA probes in three of the labs. (C) qPCR measurements of total HBV DNA and cccDNA in the DNA extracts. Bars depict the median and range across all four labs. (D) qPCR measurements of cccDNA shown separately for every lab, relative to the amount in the +PK DNA extractions. Bars depict the mean of duplicate measurements. +PK, total DNA extraction with proteinase K digestion; −PK, total DNA extraction without proteinase K digestion; cccDNA, covalently closed circular DNA; HBV, hepatitis B virus; Hirt, Hirt DNA extraction, which does not include a proteinase K digestion; qPCR, quantitative PCR; SB, Southern blot.
To evaluate the three distinct DNA extraction methods in different labs, we first analysed the liver of a highly viraemic mouse. Hirt and −PK extractions resulted in reduced HBV RI levels on the SB (figure 1B) compared with the +PK extraction. The classical pattern of a ‘slow-migrating band’ around 3.5 kb (known to contain pf rcDNA or deproteinated (dp)-rcDNA) and a ‘fast-migrating band’ at 2.1 kb (supercoiled cccDNA) was observed both after Hirt and −PK extraction procedures. The high amount of RI in the+PK extracts, however, prevented the visualisation of distinct bands. Of note, using the highly sensitive bDNA detection method, a substantial smear was still present in −PK and Hirt extracted DNA, which was absent in samples from an uninfected mouse (see example in online supplemental figure 8C), pointing out that those signals were HBV-specific. In line with the SB, qPCR measurements in all labs showed lower amounts of total HBV DNA (1.6log and 1.3log reduction) and cccDNA (0.5log and 0.4log) in Hirt and −PK extractions, respectively (figure 1C), as compared with +PK samples. The reduction of protein-bound RIs clearly improved the selectivity of cccDNA quantification by PCR (figure 1C) across all labs (figure 1D).
Normalisation of cccDNA counts to cellular levels
For all downstream cross-validation experiments, the same amount of liver tissue or cell number was used to facilitate comparison between labs and settings. In everyday practice, however, viral DNA counts have to be normalised to cell number to account for differences in the amount of input material or human hepatocyte levels when analysing chimeric livers. Accordingly, we determined the recovery of genomic DNA (gDNA) and mitochondrial DNA (mtDNA) (online supplemental figure 4A–C). Although supernatants from a classical Hirt extraction are known to harbour reduced levels of high molecular weight nucleic acids, genomic DNA was consistently extracted in our Hirt preparation, although with a trend towards a slightly reduced recovery of these DNA species in both Hirt and −PK conditions compared with the +PK extraction. However, this trend was similar across analyses, indicating that viral DNA could be normalised to gDNA or mtDNA derived from the same extraction. Moreover, a head-to-head comparison of distinct lysis buffers revealed that the TCL lysis buffer (Lucigen) led to threefold lower gDNA levels compared with RLT lysis buffer (Qiagen), although these buffers did not affect mtDNA recovery (online supplemental figure 4D–F). Such differential enrichment of gDNA will skew the normalisation of viral DNA to target cells and has to be considered in the normalisation strategy. A modified Hirt procedure using alkaline extraction,22 23 which was tested on liver tissue and cell culture samples in two labs, showed low but still quantifiable genomic DNA recovery compared with the original Hirt (online supplemental figure 5). By enriching only supercoiled DNA species, this method also depleted all HBV DNA species except for the ‘fast-migrating band’.
Comparison of different nucleases and digestion conditions in chimeric liver tissues
We next tested different nucleases and conditions to compare their efficacy. Based on existing protocols commonly used in different laboratories, T5 exonuclease was used in two settings, one being stronger (‘T5’: 10U, 45 min), the other one being gentler (‘T5 low’: 5U, 30 min) because overdigestion of cccDNA had sporadically been observed within this consortium and gentler digestion conditions were anticipated to reduce this effect. PSD was employed in a classical setting (‘PSD’: 10U, 6 hours in 20 µL total volume) and in a setting using more units in a higher volume (‘PSD high’: 30U, 2 hours in 200 µL total volume) based on experiences in tissue-derived samples. The Exo I/III combination was tested using one setting (20U ExoI, 25U ExoIII, 2 hours). All nuclease treatments reduced total HBV DNA to varying degrees (figure 2A). However, considerable variations were observed among labs. Using the classical +PK DNA extracts, ‘PSD low’ conditions were unable to effectively digest total HBV DNA. This was most likely caused by an inhibition of the enzyme by contaminants present in the extracts or inappropriate pH values, since ‘PSD high’ or the use of column-purified +PK DNA for the ‘PSD low’ condition remedied this issue (data not shown).
The effect of nuclease digestions on different DNA extracts in HBV-infected USG mouse liver tissue. The bar graphs depict qPCR measurements of total HBV DNA (A) and cccDNA (B) in all DNA extracts and labs without nuclease digestion (hatched bars) or after the indicated nuclease treatments. The bars depict the median and range of the values from three labs. Every lab performed the nuclease digestion in duplicates. The values (copies/PCR) after nuclease treatment were first normalised to the non-digested value in the respective DNA extract, then all values derived from Hirt or −PK DNA extracts were normalised to the +PK extract. (C) qPCR of mtDNA DNA via the mitochondrial gene ND2, depicted as arbitrary units normalised to the non-digested values in the respective DNA extract. (D) SB analysis in one of the labs using a second stably HBV-infected and untreated USG mouse and all three DNA extracts and four nuclease digestion conditions using HBV DNA probes (top panel) and densitometric analysis of the cccDNA band (lower panel). The ‘X’ denotes a sample that could not be quantified because of high background staining. The samples were run on two separate blots but identical amounts of non-digested DNA, with the help of which densitometry was normalised between blots. ND2, mitochondrially encoded NADH:ubiquinone oxidoreductase core subunit 2. +PK, total DNA extraction with proteinase K digestion; −PK, total DNA extraction without proteinase K digestion; cccDNA, covalently closed circular DNA; HBV, hepatitis B virus; qPCR, quantitative PCR; SB, Southern blot.
When using cccDNA selective primers, we observed a considerable reduction of cccDNA counts after nuclease treatments even in Hirt and −PK DNA extracts, indicating that residual HBV RI were still detected following these extractions in this highly infected sample (9×108 HBV DNA GE/mL) (figure 2B). Analysis of the human mitochondrial gene ND2 as a cellular surrogate marker for episomal DNA (figure 2C), revealed that mtDNA levels were also lower after nuclease digestion although to a lesser degree, indicating that the reduced detection of cccDNA was not only due to improved primer specificity but partially due to cccDNA loss. Although considerable variations were observed among labs, median cccDNA levels were 3.3-fold lower after T5 or Exo I/III treatments compared with PSD treatments across all DNA extracts and labs. To shed light on the reasons for these differences, we performed a head-to-head comparison of qPCR and SB analyses using the same undigested and digested samples (figure 2D). Both methods proved that all nuclease digestions reduced HBV DNA copy numbers (qPCR) and the typical smear present on the SB. However, PSD did not digest the ‘slow-migrating’ band where pf-rcDNA is generally found, while T5 and Exo I/III digested all HBV DNA forms except for the cccDNA-containing ‘fast-migrating’ band. Furthermore, cccDNA bands displayed very similar densities after exonuclease-based or PSD-based treatment within the same extraction, indicating that the higher cccDNA PCR counts observed after PSD digestion is due to a comeasurement of molecules present in the ‘slow-migrating’ band.
cccDNA and pf-rcDNA are detected by qPCR with the same efficiency
To assess whether the HBV DNA species present in the ‘slow-migrating’ band could be recognised by the cccDNA PCR with the same efficiency as the cccDNA itself, we performed preparative agarose gel electrophoresis, excised and analysed by PCR the different HBV DNA bands individually. Liver samples were analysed from one untreated HBV-infected mouse and two lamivudine (Lam)-treated mice (6 and 12 weeks, respectively; online supplemental figure 3B), which were expected to have lower levels of interfering RI. DNA extracted with the −PK method was subjected to SB and qPCR after ‘PSD high’ and ‘T5’ digestion. The characteristic smear on SB composed of HBV RI in highly infected samples (figure 3A) and total HBV DNA copies by PCR (figure 3B) were reduced in Lam-treated mice, thus enabling specific cccDNA measurements by qPCR before and after PSD digestion (figure 3B). However, T5 treatment removed the ‘slow-migrating’ band, leading to fourfold lower cccDNA PCR measurements than PSD treatment, while the density of the faster migrating cccDNA bands remained the same regardless of the nuclease digestion (figure 3A). Interestingly, cccDNA levels appeared to decrease in a time-dependent fashion in Lam-treated mice compared with the untreated mouse.
The pf-rcDNA is detected with the same efficiency as cccDNA by cccDNA-selective qPCR but is removed through T5 exonuclease treatment. Parallel SB (A) and qPCR (B) analysis in three stably HBV-infected USG mice (after 6 or 12 weeks of lamivudine treatment or untreated). Liver DNA was extracted with the −PK, digested with nucleases or left undigested before analysis. (A) SB using HBV DNA probes (top panel) and densitometric analysis of the pf-rcDNA (middle panel) and cccDNA band (lower panel). DNA amounts were normalised to ND2 via qPCR with human-specific primers before digestion to ensure loading of equal amounts of DNA derived from human cells. (B) qPCR measurements of total HBV DNA (grey bars) and cccDNA (checked bars) normalised to human hepatocytes via human HBB counts in the non-digested DNA extract. The cccDNA copy numbers/cell are depicted above each bar. (C) Preparative agarose gel using the same samples as in (A) and identical settings for the gel electrophoresis. The upper and lower bands were excised separately, DNA was extracted and used for qPCR analysis (total HBV DNA and cccDNA). The ratio of total HBV DNA copies to cccDNA copies was calculated for all samples derived from the upper SB band (top) and the lower band (bottom), respectively. +PK, total DNA extraction with proteinase K digestion; −PK, total DNA extraction without proteinase K digestion; cccDNA, covalently closed circular DNA; HBV, hepatitis B virus; Lam, lamivudine; NT, non-treated; qPCR, quantitative PCR; SB, Southern blot.
Identical samples were used for the preparative agarose gel, where slow-migrating and fast-migrating bands were excised separately and analysed by qPCR for total HBV DNA and cccDNA (online supplemental figure 6). The cccDNA band was detected by both PCR assays in a very similar fashion and the ratio of total HBV DNA to cccDNA copies was identical across all mice and treatments except for the untreated, undigested sample (figure 3C). This observation confirmed that the ‘fast-migrating’ band, solely composed of cccDNA, was detected by both PCR assays. Interestingly, the same pattern was observed also when the ‘slow-migrating’ band was analysed, indicating that the DNA species present in this band served as similarly efficient templates for the cccDNA PCR as cccDNA itself. The slightly higher ratios in the untreated, undigested sample were likely due to the codetection of HBV RI, presenting a visible smear only in this sample. Taken together, the higher cccDNA PCR counts determined in samples retaining the slow-migrating band strongly indicate that pf-rcDNA is an efficient template for the cccDNA PCR. The precise composition of this band, however, needs further investigations. Of note, a potential intermediate for cccDNA formation, the closed minus-strand rcDNA (CM-rcDNA), may also be present in that ‘slow-migrating’ band.17 Nicking of the cccDNA (occurring during active transcription7 or accidentally during the extraction process) would result in the loss of supercoiling and nicked cccDNA would migrate together with pf-rcDNA (see below and online supplemental figure 8C). Thus, quantification strategies aiming at detecting cccDNA contained in the ‘fast-migrating’ band should include T5 or ExoI/III nucleases, whereas PSD treatment will guarantee inclusion of potentially damaged, nicked cccDNA.
Comparison of cccDNA extraction and quantification protocols in HepG2-NTCP cells
Following the same criteria and protocols used for liver tissues, we comparatively assessed the impact of distinct DNA extraction methods on HBV DNA levels in HBV-infected HepG2-NTCP cells. Total HBV DNA qPCR measurements indicated that Hirt and −PK extractions substantially removed protein-bound HBV DNA forms by 1.7 and 1.5 logs, respectively (figure 4A). Accordingly, cccDNA counts were reduced by median 0.9 and 0.8 logs compared with+PK extractions indicating improved selectivity for cccDNA (figure 4A, B). We next applied different nuclease digestion conditions to all three DNA extraction method (figure 4C–F). In +PK extracts, the removal of RI reduced cccDNA counts by PCR. In both Hirt and −PK DNA, however, no further substantial cccDNA decline was observed after nuclease digestion, suggesting that HBV RI had been sufficiently removed through the extraction process allowing for precise cccDNA quantification. The minor reduction of cccDNA levels after all nuclease treatments was mirrored in ND2 measurements indicating a general small loss of circular DNA molecules (figure 4E). The SB in figure 4F reproduced the characteristic HBV DNA pattern similar to liver tissues although with less RI-derived smear. Nevertheless, in Hirt and −PK DNA extracts, cccDNA levels determined by PCR were 1.5-fold lower after T5 and ExoI/III compared with PSD-based methods. Again, the ‘slow-migrating’ band visualised by SB was degraded only after exonuclease-based treatments.
The effect of different DNA extraction methods and nuclease digestions on cccDNA quantification by qPCR and SB in HBV-infected HepG2-NTCP cells. (A–F) Frozen cell pellets harvested at day nine post infection were shipped to the participating labs. (A) qPCR measurements of total HBV DNA and cccDNA in the undigested DNA extracts. Bars depict the median and range across all four labs. (B) PCR measurements of cccDNA shown separately for every lab and relative to the amount in the+PK DNA extractions. Bars depict the mean of duplicate measurements. (C, D) qPCR measurements of total HBV DNA (C) and cccDNA (D) in all DNA extracts and labs either without nuclease digestion (hatched bars) or after the indicated nuclease treatments (performed in duplicates). The bars depict the median and range of the values from three labs. The values (copies/PCR) after nuclease treatment were first normalised to the non-digested value in the respective DNA extract, then all values derived from Hirt or −PK DNA extracts were normalised to the +PK extract. (E) qPCR of mtDNA DNA via the mitochondrial gene ND2, depicted as arbitrary units normalised to the non-digested values in the respective DNA extract. (F) SB using HBV DNA probes (top panel) and densitometric analysis of the cccDNA band (lower panel). The ‘X’ denotes a sample that could not be quantified because of high background staining. The samples were run on two separate blots with identical amounts of non-digested DNA, with the help of which densitometry was normalised between blots. +PK, total DNA extraction with proteinase K digestion; −PK, total DNA extraction without proteinase K digestion; cccDNA, covalently closed circular DNA; HBV, hepatitis B virus; qPCR, quantitative PCR; SB, Southern blot.
Another factor that can influence cccDNA measurements is the presence of the inoculum, which contains high amounts of HBV particles that can persist in culture. To assess the impact of input virus on cccDNA quantification, we infected HepG2-NTCP cells with a high-titre HepAD38-derived inoculum (MOI 4000) and followed them up for 9 days in the absence or presence of the entry inhibitor bulevirtide (BLV).24 The inoculum contained high levels of pf-rcDNA, which persisted until day nine in BLV-treated cells where the establishment of a productive infection was precluded (online supplemental figure 7A). This contaminant could be removed through ‘T5’ treatment, resulting in approx. twofold lower cccDNA levels detected by qPCR (online supplemental figure 7B). Of note, using this highly sensitive bDNA technology, several T5-resistant HBV DNA species were detected by SB. Interestingly, these forms are not visible under alkaline Hirt extraction conditions (online supplemental figure 5B).
Impact of DNA quality and storage conditions on cccDNA quantification
Nuclease treatment bears some risks and should be employed with caution. First, we observed high variations in terms of digestion efficiency between labs and even between replicates within one lab (data not shown). Second, excessive digestion conditions can lead to overdigestion of cccDNA. Third, the quality of the extracted DNA is of utmost importance since nicking of the cccDNA will make it a template for T5 ExoIII, but not for PSD, while all enzymes (except ExoI) will act on damaged cccDNA containing double strand breaks. Examples and explanations for risks associated with nuclease treatments can be found in online supplemental results and figure 8.
Another factor that can influence the quality of extracted DNA derives from the storage conditions. In this study, we used freshly frozen, cryopreserved samples; however, cryopreservation of liver biopsies is rarely feasible in clinical studies. Instead, tissue preservatives are used. To mimic these conditions, cryopreserved small liver pieces from two HBV-infected untreated humanised mice and two siRNA-treated mice19 were moved to Allprotect, a common tissue preservative form Qiagen, and stored at −20°C for several months before being extracted with −PK protocols and digested with ‘PSD high’ or ‘T5’ (figure 5A). Compared with the matched fresh-frozen samples analysed in parallel, the Allprotect-stored samples showed an overall higher reduction of cccDNA levels after nuclease digestion by qPCR, with PSD reaching approximately 1log and T5 2log reductions. For comparison, the cccDNA quantified in freshly frozen liver tissue from siRNA-treated mice were not reduced by PSD treatment indicating sufficiently low HBV DNA levels for specific cccDNA quantification in this setting. Accordingly, the reduction of mtDNA was even more pronounced than the drop in cccDNA following Allprotect storage, while mtDNA levels decreased only modestly in DNA from freshly frozen tissue (figure 5B). These data hint at a general DNA damage involving both nicking and double strand breaks and warrants extreme caution when using nuclease digestion in Allprotect-stored samples.
Nuclease treatment bears the risk of overdigesting cccDNA in Allprotect-stored samples. Comparison of nuclease digestion in fresh-frozen liver tissue or liver tissue preserved in Allprotect. Liver DNA was extracted using the −PK method from two untreated HBV-infected mice and two mice treated with siRNA targeting all HBV transcripts for 6 weeks. DNA was digested with PSD or T5 exonuclease and subjected to qPCR for cccDNA (A) and ND2 (B). Every dot depicts a single mouse, bars the median. +PK, total DNA extraction with proteinase K digestion; −PK, total DNA extraction without proteinase K digestion; cccDNA, covalently closed circular DNA; HBV, hepatitis B virus; qPCR, quantitative PCR; siRNA, small interfering RNA.
Discussion and recommendations for cccDNA quantification
To improve our knowledge on cccDNA biology, validate serological biomarkers and assess the impact of new therapies on cccDNA levels and activity, reliable and standardised quantitative cccDNA assays are paramount. To achieve this goal in a more reliable manner then previously possible, it is critical to compare protocols and perform cross-validation experiments across laboratories. The results of this concerted effort provide novel information to guide researchers and clinical trial investigators in designing and performing protocols to accurately quantify cccDNA.
The comparative analysis of parallel qPCR and SB results was a strength of our study design and allowed us to understand how each method influences the outcome and to identify the most critical points. In figure 6, we summarised the methodologies that were validated within our consortium and the observations relevant to cccDNA quantification. Because all protocols must be adapted to the sample type, infection level and treatment, there is not one protocol that suits all needs. Instead, we provide guidance and general recommendations, which are summarised in table 1.
Schematic presentation of the main findings of the study. The flow chart summarises the experimental flow for cccDNA quantification and highlights optimal and suboptimal results for both qPCR and SB, including the major caveats. cccDNA, covalently closed circular DNA; qPCR, quantitative PCR; SB, quantitative PCR.
Recommendations for cccDNA quantification
Through the comparison of three DNA extraction methods, we could show that not only the Hirt procedure, but also the omission of the PK digestion step from a commercial kit for fast nucleic acid purification led to efficient reduction of HBV DNA RI, thus improving the selectivity of cccDNA PCR measurements. However, −PK or Hirt-based extractions13 22 23 might not be sufficient for specific quantification in highly replicative settings and nuclease treatment will still be needed. Therefore, we recommend performing both a DNA extraction that reduces RI and virions, a nuclease treatment to further reduce non-cccDNA species and a PCR design spanning the cohesive ends region of the rcDNA. The use of such a cccDNA-selective PCR is further supported by the fact that none of the extraction methods tested (−PK, Hirt, alkaline Hirt) was able to guarantee a complete depletion of genomic DNA, which could represent a source of integrated HBV sequences. All nucleases used here reduced RI amounts to various degrees, even when −PK or Hirt DNA extraction procedures had already lowered HBV RI levels. Nevertheless, all nucleases showed specific strengths and weaknesses. PSD did not perform well in +PK-extracted liver DNA, possibly due to impurities present in these extracts, whereas T5 and Exo I/III exonucleases performed similarly well although high concentrations and incubation times bear the risk of cccDNA overdigestion. Additionally, the efficiency of the tested nucleases varied both among labs and even within a given lab. The reasons for these variations are unclear, but could depend on small handling differences, input DNA variations and presence of sample impurities. Moreover, we showed that tissue preservatives are particularly prone to damaging the cccDNA. Consequently, nucleases induced strong reduction of cccDNA PCR levels in Allprotect-stored samples. Since patient liver biopsies are often stored in Allprotect or similar tissue preservatives, omission of nuclease digestions should be considered. Ideally, liver biopsy samples should be fresh-frozen and cryopreserved. Thus, digestion conditions have to be carefully titrated to avoid unspecific side activities potentially leading to cccDNA digestion or to insufficient digestion of other HBV DNA forms. Since it is laborious to determine whether a nuclease treatment is needed for every sample and condition, we recommend including it for fresh and fresh-frozen samples.
Additionally, we observed that PSD did not digest the ‘slow-migrating’ band, although this treatment removed most RI from −PK and Hirt samples. Furthermore, we demonstrated that without nuclease treatment, the pf-rcDNA is still present even in the setting of low replication levels and that this ‘slow-migrating’ band is detected by the cccDNA-selective PCR. We therefore recommend using exonuclease-based treatment when detection of this band is not wanted. The precise composition of this ‘slow-migrating’ band detected by SB is not entirely known but it appears to contain a mixture of molecules, which underwent removal of the covalently attached HBV polymerase but that still are unrepaired or at different phases of the repair process, such as DP-rcDNA or CM-rcDNA.3 4 Fully repaired HBV DNA molecules that are not yet covalently closed and supercoiled (immature cccDNA) will also appear within this band.3 4 In cell culture experiments, a prominent part of the ‘slow-migrating’ band might consist of input virus. Since it is unclear which factors and circumstances (infection status or treatments) may affect the conversion of rcDNA to cccDNA, the intensity of the ‘slow-migrating’ band may also vary between samples and experimental settings. Additionally, the choice of the extraction method will influence the detected population of HBV DNA molecules. In the side-by-side comparison of the original Hirt with a modified alkaline Hirt known to preferentially extract supercoiled DNA, the ‘slow-migrating’ bands—potentially resulting from circular HBV DNA with different degrees of supercoiling—appeared lost (online supplemental figure 5). Thus, further studies are needed to assess the biological meaning of such differences and cccDNA quantification strategies including the ‘slow-migrating’ band are justified.
For in vitro infection studies, we recommend employing exonuclease-based treatments (T5 or Exo I/III) to remove the non-cccDNA species present in the inoculum. This contaminant was present as a prominent ‘slow-migrating’ band on SB and persisted until day nine post infection of HepG2-NTCP cells despite the use of an entry inhibitor (bulevirtide), known to abrogate the establishment of productive HBV infection.24 Although the example shown here used an extremely high MOI (4000), T5 exonuclease efficiently reduced this contaminant in Hirt-extracted samples allowing specific cccDNA measurements by qPCR. Thus, it is advisable to include an entry inhibitor control in cell culture to rule out interferences due to input HBV DNA.
Because cccDNA that contains nicks or double strand breaks will become a template for nucleases and migrate together with pf-rcDNA on SB, it is of utmost importance to assess the quality of the DNA in every extract. We recommend using mitochondrial DNA—an episomal dsDNA—as a surrogate for cccDNA and to quantify a mitochondrial gene by PCR in every sample before and after nuclease digestion. Based on our experience, a reduction of less than 1log indicates sufficiently intact DNA, while more than 1log reduction indicates damaged DNA. Overdigestion of mtDNA after exonuclease-based but not PSD-based digestions, indicates that the DNA contains single-strand nicks. Consequently, cccDNA will also be nicked and reduced after exonuclease treatment. PCR assessment of mtDNA levels will indicate the quality of the DNA in the sample and after the extraction.
We acknowledge that the present methods to quantify cccDNA are not optimal yet and encourage the optimisation or development of new methods for both DNA extraction and cccDNA quantification.25 Ideally, improved PCR-based methods should have higher specificity for cccDNA in the presence of contaminating HBV DNA species and should be robust, sensitive and suitable for all virus isolates. A milestone will be reached when novel procedures allow to avoid the use of nucleases. The implementation of digital PCR (dPCR) instead of qPCR has allowed a significant increase in sensitivity and brought the major advantage of absolute quantification without PCR standards, thus limiting the variability thereof.12 26 However, the limited specificity for cccDNA in a complex mixture of sequence-identical HBV DNA species remains even when using dPCR. Novel methodologies to visualise cccDNA in infected cells are being developed and can complement quantitative analyses and serve to corroborate PCR and SB results.27 Additionally, more research is needed to ascertain the nature of all HBV DNA forms, their role in the life-cycle of HBV and how they are influenced by various manipulations.
Data availability statement
All data relevant to the study are included in the article or uploaded as online supplemental information.
Ethics statements
Patient consent for publication
Acknowledgments
We are grateful to Tassilo Volz, Elena S. Kim, Soujuan Wang, Romina Bester, Chloé Saint Jean, Laura Dimier, Maud Michelet, Emmanuel Combe and Moritz Grunwald for their great technical support and to Guofeng Cheng, William E. Delaney and Dieter Glebe for intellectual input and constructive discussions. We would like to thank ICE-HBV, ANRS and DZIF for their support of this initiative.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Twitter @BingqianQu, @ProtzerUlrike
Contributors LA, BT, JL, MLu, HG, SU, UP, SPF, MLe, FZ and MD initiated the study, designed main experiments and discussed data. CK, JL, BT and MY generated cell culture samples. LA generated liver tissue samples. LA, BT, MY, CK and BQ performed cross-validation experiments. MY, LA and JL performed Southern blot. LA analysed the data. LA and MD wrote the manuscript supported by BT, MLe and FZ. All authors read and corrected the manuscript. MD acted as guarantor of the study.
Funding The study was in part supported by the German Center for Infection Research (DZIF-BMBF 05.816; 05.820) to MD, UP and SU and a collaborative grant (ANRS-DZIF MoU 2016) to FZ, ML, BT, MD, UP, and SU. BT, MLe, and FZ were supported by a public grant overseen by the French National Research Agency (ANR) as part of the second “Investissements d’Avenir” program (reference: ANR-17-RHUS-0003). UP and SU were also supported by the German Research Foundation (DFG) via SFB-TRR 179 272983813. MD was also supported by the DFG via SFB841 A5. HG received funding from the US National Institutes of Health (R01AI110762 and R01AI123271). All funding sources supporting the work are acknowledged.
Competing interests LA, BT, MY, CK, BQ, MLu, JL and MLe declare no conflict of interest. FZ: Consultancy for Aligos, Antios, Arbutus, Assembly, Blue Jay, Enanta, Enochian, Gilead, GSK, Zhimeng. HG: Consultancy for Aligos and Assembly; shareholder of Arbutus. UP: Consultancy for Aligos, Arbutus, Gilead, GSK, Sanofi, Vaccitech, VirBio. Research collaboration with Abbott and Roche. Share holder and board member of SCG Cell Therapy. MD: research collaboration with Gilead, MYR GmbH/Hepatera and HUMABS BioMed; consultancy for Gilead and Aligos. SPF: employment and shareholder Gilead Sciences.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Provenance and peer review Not commissioned; externally peer reviewed.
© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.