Skip to main content

Genome-wide capture sequencing to detect hepatitis C virus at the end of antiviral therapy

Abstract

Background

Viral relapse is a major concern in hepatitis C virus (HCV) antiviral therapy. Currently, there are no satisfactory methods to predict viral relapse, especially in the era of direct acting antivirals in which the virus often quickly becomes undetectable using PCR-based approaches that focus on a small viral region. Next-generation sequencing (NGS) provides an alternative option for viral detection in a genome-wide manner. However, owing to the overwhelming dominance of human genetic content in clinical specimens, direct detection of HCV by NGS has a low sensitivity and hence viral enrichment is required.

Methods

Based on template-dependent multiple displacement amplification (tdMDA), an improved method for whole genome amplification (Wang et al., 2017. Biotechniques 63, 21–27), we evaluated two strategies to enhance the sensitivity of NGS-based HCV detection: duplex-specific nuclease (DSN)-mediated depletion of human sequences and HCV probe-based capture sequencing.

Results

In DSN-mediated depletion, human sequences were significantly reduced in the two HCV serum samples tested, 65.3% → 55.6% → 33.7% (#4727) and 68.6% → 56% → 21% (#4970), respectively for no normalization, self- and driver-applied normalization. However, this approach was associated with a loss of HCV sequences perhaps due to its micro-homology with the human genome. In capture sequencing, HCV-mapped sequencing reads occupied 96.8% (#4727) and 22.14% (#4970) in NGS data, equivalent to 1936x and 7380x enrichment, respectively. Capture sequencing was then applied to ten serum samples collected at the end of HCV antiviral therapy. Interestingly, the number of HCV-mapped reads was significantly higher in relapsed patients (n = 5) than those from patients with sustained virological response (SVR) (n = 5), 102.4 ± 72.3 vs. 2.6 ± 0.55, p = 0.014.

Conclusions

Our data provides concept evidence for a highly sensitive HCV detection by capture sequencing. The abundance difference of HCV sequencing reads at the end of HCV antiviral therapy could be applied to predict treatment outcomes.

Peer Review reports

Background

The development of direct acting antivirals (DAAs) is a milestone in antiviral therapy for hepatitis C virus (HCV) infection. Given the high rate of sustained virological response (SVR) among patients receiving DAA regimens, HCV infection is now claimed to be a curable disease [1]. In reality, however, there remains a small portion of patients, around 5%, who have the virus relapsed after treatment cessation [2]. HCV is circulating as a heterogenous population, sometime termed viral quasispecies [3]. In such a population, relapse is often accompanied with the enrichment of resistance-associated mutations that may foster the emergence of new drug-resistant viral strains [4]. Re-treatment with DAAs among relapsed patients confers a reduced SVR rate [5]. Viral relapse is thus a major concern in the era of DAAs. As HCV completes its life cycle solely in cytoplasm and has no evidence of cellular integration [6], its survival depends on continuous replication. Consequently, patients at the time of HCV relapse should maintain a low-level viral replication that is beyond the detection limit of current PCR-based methods [7]. In this setting, an alternative option is next-generation sequencing (NGS) that may detect the virus in a genome-wide manner. A drawback of using NGS for viral detection is its low sensitivity owing to the overwhelming dominance of human and commensal microbial genomes in clinical specimens [8]. Using specially designed primers, our lab previously developed template-dependent multiple displacement amplification (tdMDA) that eliminates primer-mediated artifacts [9]. Based on this improved whole genome amplification method, the current study attempted to increase the sensitivity of NGS through two enrichment strategies, depletion of human sequences or target capture by probe-based hybridization. The optimized method was then applied to clinical serum samples to explore the feasibility in detecting extremely low viral titers, thereby providing an on-treatment tool in predicting the outcomes of HCV antiviral therapy.

Methods

Patient samples

A total of 12 patient serum samples were included in the current study. Of the 12, two samples from our previous study, #4727 and #4970 (both genotype 1a), were collected prior to antiviral therapy and had viral loads respectively at 8.72 × 106 and 5.24 × 105 copies/mL as quantified by Roche Amplicor HCV Monitor (v2.0) [10]. These two samples were used for methodological optimization. The other ten samples were collected at the end of treatment (EOT) from the patients receiving either pegIFNα2a/ribavirin (n = 3) or a DAA regimen (boceprevir and telaprevir) (n = 7). All patients were infected with HCV genotype 1a. The samples were respectively collected from a completed clinical study (pegIFNα2a) and Hepatitis C Therapeutic Registry and Research Network (TARGET) (DAA regimen) ClinicalTrials.gov Identifier: NCT01474811 [10, 11]. Of the ten patients, five achieved SVR (2 with IFN and 3 with DAA) and the other five patients had virus relapses (1 with IFN and 4 with DAA). The study protocol was reviewed and approved by the Saint Louis University Institutional Review Board (IRB protocol 10,592).

RNA extraction, reverse transcription (RT), tdMDA, and Illumina sequencing

RT-tdMDA was conducted as previously described [9]. In brief, total RNA was extracted from 140 μL of serum (samples #4727 and #4970) and eluted into 60 μL Tris buffer (pH 8.5) using the QIAamp Viral RNA Mini kit (Qiagen, Valencia, CA). An aliquot of 10.6 μL of extracted RNA was mixed with 9.4 μL RT matrix consisting of 1x SuperScript III buffer, 10 mM DTT, 80 μM of exonuclease-resistant random pentamer primers with the 5′ end blocked by C18 spacer [9], 2 mM dNTPs (Epicentre), 20 U (0.5 μL) of RNaseOUT Recombinant Ribonuclease Inhibitor (Invitrogen), and 200 U of SuperScript III reverse transcriptase (Life technologies). After incubation at 37 °C for 30 min and 50 °C for 30 min, the reaction was inactivated at 70 °C for 15 min. An aliquot of 4 μL of RT was used for tdMDA in 40-μL reaction consisting of 1x phi29 DNA polymerase buffer, 1 mM of dNTPs, 80 μM of random pentamer primers as used in the RT, and 20 units of phi29 DNA polymerase (New England Biolabs, Ipswich, MA). The reaction was incubated at 28 °C for 14 h and then terminated by being heated at 65 °C for 15 min. After purification with the QIAamp DNA mini kit (Qiagen), RT-tdMDA product was subjected to library construction with the Nextera XT DNA Sample Preparation kit (Illumina, San Diego, CA) and sequenced on the Illumina NextSeq 500 platform (1 × 250-nt single reads and mid-output) at MOgene (St. Louis, MO).

Enrichment of HCV detection via the depletion of human sequences

This approach utilized duplex-specific nuclease that showed optimal activity at high temperatures with a strong preference for cleavage of double-stranded DNA in comparison to single-strand DNA [12]. Mixed human genomic DNA (male and female) (Promega, Madison, WI) was sheared into ~ 150 bp fragments using the Covaris microTUBE device, followed by purification with AMPure XP beads (Beckman Coulter, Indianapolis, IN), first at 0.8x for removing large-size fragments and then at 2.5x for the harvest (Supplementary Figure 1). Next, using an optimized protocol (Supplementary Figure 2), 300 ng of purified RT-tdMDA product was fragmented through an 8-min heating at 98 °C and then mixed with 1200 ng (1:2 or 1:4) of sheared human genomic DNA (the driver) in a 20-μL reaction containing 1x hybridization buffer [50 mM HEPES (pH 7.5) and 0.5 M NaCl]. After being heated at 98 °C for 3 min, the reaction was incubated at 68 °C for 4 h. At the end of the incubation, 8 μL of pre-warmed solution, including 2 μL of DSN (1 U/μL) (Evrogen via Axxora, Farmingdale, NY) and 2.8 μL of 10x DSN buffer, was added for an additional 30-min incubation. The reaction was terminated by adding 28 μL (1:1) of 2x DSN stop solution and then subjected to two-step purification, first by the QIAquick PCR Purification Kit (Qiagen) and then AMPure XP beads at 0.55x that was determined by titration for the size exclusion of < 300 bp fragments (Supplementary Figure 3). The entire purified product was used for 4-h tdMDA. These experiments were repeated multiple times and the efficiency of enrichment was initially estimated using an in-house PCR based on HCV 5’UTR (30 cycles of single round) [13]. Selective tdMDA product was then processed into Illumina sequencing as described above.

Enrichment of HCV detection via capture sequencing

This strategy was applied to samples #4727 and #4970. First, 384 full-length HCV genotype 1a genomes were extracted from the Los Alamos HCV database to generate a consensus sequence [14], which was 9508 nt in length without 120-nt poly(U)/C tract (Supplementary File 1). The 2x HCV titling probes (each 120 nt with 60 nt overlap) were designed using an in-house script and synthesized at Integrated DNA Technologies (IDT) (Coralville, IA). After regular library preparation from RT-tdMDA product, the library was hybridized with probes using the xGen Hybridization and Wash Kit (IDT), followed by PCR amplification, quantitation, and finally the sequencing run as described above.

Data analysis

Raw sequence reads in fastq format from samples were filtered in PRINSEQ (v0.20) for quality control, including read length ≥ 70 nt, mean read quality score ≥ 25, low complexity with DUST score ≤ 7, ambiguous bases ≤1%, and all types of duplicates [15]. To reduce non-specific HCV mapping, quality reads were first subtracted by human sequences [The National Center for Biotechnology Information (NCBI) GRCh38 build] [16] using the Bowtie 2 mapper [17]. Remaining reads were then mapped onto the HCV genotype 1a probe sequence. Read alignment was converted, sorted, and indexed into BAM file in SAMtools and was viewed in BamView [18, 19]. When necessary, HCV consensus sequences were called using two-step procedure as described in our previous study [20]. HCV consensus sequences, together with HCV sequences generated through cloning and Sanger sequencing from the same samples [10], were used for the construction of the phylogenetic tree in MEGA (Molecular Evolutionary Genetics Analysis) [21].

Analysis of micro-homology between HCV 1a consensus sequence and the human genome

Micro-homology may affect the efficiency of the enrichment method, which can be either DSN-mediated depletion of human genomic sequences or HCV capture sequencing. To analyze its potential roles, HCV 1a probe sequence was split into 60-nt fragments with 30-bp overlaps using pyfasta [22]. These sequences were searched for micro-homology against the human genome (GRCh38 build) by blastn at a word size of 6. Based on the blastn output, human genomics sequences were extracted using HOMER with additional steps for the removal of duplicates [23]. Non-redundant human genomics sequences were calculated for the melting temperatures (Tm) using the nearest neighbor method [24].

Detection of HCV at the end of antiviral therapy by capture sequencing

After comparing the sensitivity between the two enrichment options, capture sequencing was selected to detect HCV from the 10 serum samples collected at EOT. All serum samples were tested for HCV using in-house RT-PCR protocol (total 70 cycles) based on 5′ untranslated region (5’UTR) as described previously [11]. RT-tdMDA from a healthy donor serum and water were included as the negative control. The experimental procedure for capture sequencing was the same as what was described for samples #4727 and #4970.

Statistical analysis

Student’s t-test. Data were expressed as the mean ± SD (standard deviation), and p < 0.05 was considered statistically significant.

Data availability

Raw sequence data in fastq format from the current study were deposited in the NCBI Sequence Read Archive (SRA) under BioProject ID: PRJNA545854.

Results

Enrichment comparison: DSN-mediated depletion vs. capture sequencing

Following RT-tdMDA, samples #4727 and #4970 were subjected to Illumina sequencing under three options, regular library preparation, DSN-mediate depletion of human genomic sequences prior to library construction, and capture hybridization. Under regular sequencing, only 585 (0.05%) and 40 (0.003%) reads were mapped onto the HCV genotype 1a consensus sequence for samples #4727 and #4970 respectively (Supplementary Table 1). Both samples had reads from the human genome dominate the data, 65.3% for #4727 and 68.6% for #4970. In DSN-based normalization, HCV PCR didn’t indicate an enrichment effect (data not shown). This was confirmed by subsequent Illumina sequencing. While human reads were reduced to 55.6 and 56% in self-normalization and 33.7 and 21% with the use of the driver for samples #4727 and #4970 respectively, self- and driver-based normalization were associated with the loss of HCV-mapped reads, 585 → 19 → 3 for #4727 and 40 → 4 → 1 for #4970 (Fig. 1). In contrast, capture sequencing had 96.8 and 22.14% reads mapped onto the HCV genotype 1a consensus sequence, which was equivalent to approximately 1936x and 7380x enrichment for #4727 and #4970 respectively (Fig. 1).

Fig. 1
figure1

HCV enrichment by DSN-mediated normalization and capture sequencing. Reads for human genome was mapped using Bowtie 2 onto NCBI GRCh38 build and calculated as percentages (left Y-axis) while HCV was mapped using 184 reference sequences from the Los Alamos HCV database (right Y-axis)

Micro-homology may contribute the loss of HCV sequences in DSN-mediated sequencing

Two factors were speculated to be responsible for the loss of HCV-mapped reads in DSN-mediated normalization. First, DSN activity appeared to be harmful for phi29 DNA polymerase. When DSN was inactivated through the procedure recommended by the manufacturer, tdMDA after purification with AMPure XP beads was inefficient (data not shown). For a complete removal of DSN, an additional purification step using Qiagen spin columns was required after DSN digestion. More purification steps may result in a loss of low-abundant template, such as HCV sequences. Another factor was the micro-homology between HCV and the human genome. At ≥95% similarity, the HCV genotype 1a consensus sequence baited 1784 unique regions, ranging from 11 to 25 nt, on the human genome with Tm higher than 60 °C. While these regions were scattered across the entire HCV genome (Fig. 2a), it was interesting to note that there was uneven distribution among chromosomes with 78% of hits contributed by chromosome X (Fig. 2b). In the existence of excessive human DNA, digestion of HCV sequences owing to the micro-homology appears to be a more plausible explanation.

Fig. 2
figure2

Micro-homology between human genome and HCV genotype 1a consensus sequence. Micro-homologous regions were shown along with the entire HCV genome with corresponding Tm values (a). Distribution of micro-homologous regions across human chromosomes with Tm values higher than 60 °C

Authentic recovery of HCV genomes by capture sequencing

Samples #4727 and #4970 were studied for viral quasispecies through cloning and Sanger sequencing of the 1380-bp amplicon across HCV Core, E1 and E2 regions [8]. After the removal of primer sequences, the 1340-bp consensus sequence derived from multiple clones shared 90% (#4727) and 88% (#4970) nucleotide similarity with the HCV genotype 1a consensus sequence that was used for probe design. However, the consensus sequences from cloning and capture sequencing had much higher nucleotide similarity, 99.4 and 99.25% for samples #4727 and #4970 respectively. Accordingly, the consensus sequences called from capture sequencing were phylogenetically well-clustered with cloned HCV sequences in both samples (Fig. 3). In addition, capture sequencing gave a full coverage across the entire HCV genome without any gaps (data not shown). Thus, probe design at the level of HCV subtypes could tolerate sequence divergence among HCV strains to ensure accurate recovery of HCV genomes.

Fig. 3
figure3

The Neighbor-Joining tree of 1340-bp HCV sequences. The tree was constructed using the Tamura 3-parameter method with bootstrap values (1000 replicates) shown on major branches. Arrows indicated the consensus sequences called from capture sequencing. Also indicated was the HCV genotype 1a consensus sequence that was used for probe deign in capture sequencing

Detection of HCV at the end of antiviral therapy

The experimental protocol for HCV capture sequencing was applied to the 10 serum samples collected at the end of antiviral therapy. These samples tested negative using a commercial kit (COBAS AmpliPrep/COBAS TaqMan HCV Test, v2.0) and in-house RT-nested PCR (data not shown). Using the same probes from HCV genotype 1a consensus sequence, capture sequencing detected 3, 3, 3, 2, and 2 reads in five patients with SVR, which was a sharp contrast to the five patients with viral relapse that had 208, 148, 48, 56, and 52 reads mapped (Fig. 4) (Supplementary Table 1). No HCV-mapped reads were found in two negative controls (Fig. 4). Quantitatively, SVR patients had an average 2.6 ± 0.55 reads in comparison to 102.4 ± 72.3 from relapsed patients (p = 0.014 two-tailed Student t test). And the similar result was obtained after the normalization to the mean number of total reads among these patients (111.4 ± 81.5 vs. 2.64 ± 0.27, p = 0.017). Read mapping showed the lack of a consistent pattern with the reads scattered across the HCV genome (Fig. 4). Finally, the average rate of duplicate reads among 10 samples was 29.8 ± 7.8%. Based on this rate, the estimated library size was 1.01 ± 0.2 million of reads under a Poisson distribution [25]. Therefore, capture sequencing in these kinds of samples have small library sizes that could be easily achieved with current NGS service.

Fig. 4
figure4

Alignment of HCV-mapped reads on the HCV genotype 1a consensus sequence. Patients were grouped based on treatment outcomes (SVR or relapse) and indicated with treatment regimens (IFN or DAA)

Discussion

Based on tdMDA, the current study has explored two options to enhance the sensitivity of clinical HCV sequencing from patient serum samples. DSN has a good activity at high temperature ranging from 60 °C to 80 °C [12]. At these temperatures, abundant DNA species after denaturation in a sample have more chances to form dsDNA that are digested while rare species are intact in the form of ssDNA [12]. Consequently, a normalization among DNA species is achieved. This unique characteristic of DSN is widely used for NGS library normalization and for the depletion of individual genes [26,27,28,29]. Our study confirmed a DSN-mediated normalization in RT-tdMDA product with or without the use of the sheared human genomic DNA (the driver). As illustrated by the ratios of human sequences in NGS data (Fig. 1), inclusion of the driver appears to be more effective than self-normalization (no driver). However, in either self- or driver-mediated normalization, HCV sequences were digested simultaneously, potentially owing to the micro-homology between HCV and the human genome (Fig. 2). In contrast, capture sequencing significantly enriched the recovery of HCV sequences and had no notable off-target effect. Sample #4727 has a HCV on-target rate of 96.8%, which is 4.37 times higher than that from sample #4970 (22.14%) (Fig. 1). Because a pivotal step in capture sequencing is probe-based hybridization, its efficiency is affected by the target’s abundance and sequence divergence from the probe. In our study, both samples #4727 and #4970 belongs to HCV genotype 1a and have similar genetic distance to the HCV genotype 1a consensus sequence (probe) (Fig. 3). Therefore, difference in HCV on-target rate could be attributed to viral titer that was 16.6 times higher in sample #4727 than #4970. Indeed, sample #4970 had a greater enrichment efficiency (7380x) in comparison to #4727 (1936x), illustrating the power of capture sequencing.

Our data is consistent with recent reports that use similar strategy for HCV sequencing [30, 31]. The enrichment of our samples is even higher than those reported [30, 31]. Besides the complete elimination of primer-originated artifacts, tdMDA uses phi29 DNA polymerase that favors the amplification of long templates [32]. The utilization of tdMDA for whole transcriptome amplification would intrinsically enrich HCV sequences in comparison to direct library preparation of serum RNA. Finally, by comparing the sequences from amplicon-based cloning/Sanger sequencing, the study also provides evidence for the sequence authenticity using capture probes designed at the level of HCV subtype (1a) consensus sequence (Fig. 3).

After the application of capture sequencing to serum samples collected at EOT, there was significant difference between patients with SVR and viral relapse in terms of the number of HCV reads detected (Fig. 4). The reads showed a scattered pattern without a full coverage across HCV genome, supporting extraordinarily low viral loads in blood. Indeed, these samples were all negative by using RT-PCR based on HCV 5′ UTR. In our experience, a 70-cycle of nested RT-PCR is even more sensitive than Amplicor HCV Monitor (v2.0) that has a detection limit at 50 IU/mL, equivalent to 45 copies/mL [13, 33,34,35]. Using commercial kits, however, HCV is reported to be detectable and even quantifiable at EOT in a small percentage of patients who achieve SVR [36, 37]. Regardless of drug regimens, HCV antiviral therapy requires certain duration in order to eradicate the virus, usually 48 weeks for IFN-based therapy or 8–12 weeks for DAAs. In view of evolutionary biology, such duration represents a bottleneck that fosters the generation and accumulation of deleterious mutations, so-called Miller’s ratchet [38,39,40], which eventually renders the loss of replicative fitness, thereby virus extinction. Our previous studies have documented the link between HCV mutation loads and the outcomes of IFN-based antiviral therapy [20, 41]. When deleterious mutations approach a level that prevents the rescue of replicative fitness after treatment cession, the virus goes extinct in spite of a positive HCV detection at EOT. At this point, it is necessary to discriminate quantitative and qualitative HCV genomic characteristics between SVR and relapse when detectable at EOT. The current study includes a small number of patients (5 SVR and 5 relapse), which represents a major limitation. in the current study. Further study with more patient samples is required for a solid clinical and evolutionary interpretation in terms of HCV detection by capture sequencing at EOT.

Conclusions

In summary, DSN-mediated depletion of human sequences is an inefficient strategy for the enrichment of HCV sequences in NGS. In contrast, the current study provides preliminary evidence that tdMDA, in combination with HCV capture sequencing, could be used as a highly sensitive tool in detecting HCV. As predicted by evolutionary principles, tailoring treatment duration based on HCV detection at EOT could reduce and even eliminate the incidence of viral relapse. Given the high cost of medicines like DAAs, the method described here allows a cost-effective personalized management in HCV clinic.

Availability of data and materials

Raw sequence data from the current study were deposited in the NCBI Sequence Read Archive (SRA) under BioProject ID: PRJNA545854. The HCV genotype 1a consensus sequence was provided in the Supplementary information.

Abbreviations

NGS:

next-generation sequencing

tdMDA:

Template-dependent multiple displacement amplification

HCV:

Hepatitis C virus

DSN:

Duplex-specific nuclease

DAA:

Direct acting antivirals

SVR:

Sustained virological response

EOT:

End of treatment

TARGET:

Hepatitis C Therapeutic Registry and Research Network

References

  1. 1.

    Aghemo A, Buti M. Hepatitis C therapy: game over. Gastroenterology. 2016;151:795–8.

    PubMed  Google Scholar 

  2. 2.

    Vermehren J, Park JS, Jacobson IM, Zeuzem S. Challenges and perspectives of direct antivirals for the treatment of hepatitis C virus infection. J Hepatol. 2018;69:1178–87.

    CAS  PubMed  Google Scholar 

  3. 3.

    Domingo E, Perales C. Viral quasispecies. PLoS Genet. 2019;15:e1008271.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Sagnelli E, Starace M, Minichini C, Pisaturo M, Macera M, Sagnelli C, Coppola N. Resistance detection and re-treatment options in hepatitis C virus-related chronic liver diseases after DAA-treatment failure. Infection. 2018;46:761–83.

    PubMed  Google Scholar 

  5. 5.

    Pawlotsky JM. Hepatitis C virus resistance to direct-acting antiviral drugs in interferon-free regimens. Gastroenterology. 2016;151:70–86.

    CAS  PubMed  Google Scholar 

  6. 6.

    Lindenbach BD, Rice CM. Unravelling hepatitis C virus replication from genome to function. Nature. 2005;436:933–8.

    CAS  PubMed  Google Scholar 

  7. 7.

    Maasoumy B, Vermehren J, Welker MW, Bremer B, Perner D, Höner Zu Siederdissen C, et al. Clinical value of on-treatment HCV RNA levels during different sofosbuvir-based antiviral regimens. J Hepatol. 2016;65:473–82.

    CAS  PubMed  Google Scholar 

  8. 8.

    Houldcroft CJ, Beale MA, Breuer J. Clinical and biological insights from viral genome sequencing. Nat Rev Microbiol. 2017;15:183–92.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Wang W, Ren Y, Lu Y, Xu Y, Crosby SD, Di Bisceglie AM, Fan X. Template-dependent multiple displacement amplification for profiling human circulating RNA. Biotechniques. 2017;63:21–7.

    PubMed  PubMed Central  Google Scholar 

  10. 10.

    Chambers TJ, Fan X, Droll DA, Hembrador E, Slater T, Nickells MW, et al. Quasispecies heterogeneity within the E1/E2 region as a pretreatment variable during pegylated interferon therapy of chronic hepatitis C virus infection. J Virol. 2005;79:3071–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Gordon SC, Muir AJ, Lim JK, Pearlman B, Argo CK, Ramani A, et al. Safety profile of boceprevir and telaprevir in chronic hepatitis C: real world experience from HCV-TARGET. J Hepatol. 2015;62:286–93.

    CAS  PubMed  Google Scholar 

  12. 12.

    Shagina I, Bogdanova E, Mamedov IZ, Lebedev Y, Lukyanov S, Shagin D. Normalization of genomic DNA using duplex-specific nuclease. BioTechniques. 2010;48:455–9.

    CAS  PubMed  Google Scholar 

  13. 13.

    Fan X, Solomon H, Poulos JE, Neuschwander-Tetri BA, Di Bisceglie AM. Comparison of genetic heterogeneity of hepatitis C viral RNA in liver tissue and serum. Am J Gastroenterol. 1999;94:1347–54.

    CAS  PubMed  Google Scholar 

  14. 14.

    Kuiken C, Yusim Y, Boykin L, Richardson R. The Los Alamos hepatitis C sequence database. Bioinformatics. 2005;21:379–84.

    CAS  PubMed  Google Scholar 

  15. 15.

    Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27:863–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Guo Y, Dai Y, Yu H, Zhao S, Samuels DC, Shyr Y. Improvements and impacts of GRCh38 human reference on high throughput sequencing data analysis. Genomics. 2017;109:83–90.

    CAS  PubMed  Google Scholar 

  17. 17.

    Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–9.

    PubMed  PubMed Central  Google Scholar 

  19. 19.

    Carver T, Harris SR, Otto TD, Berriman M, Parkhill J, McQuillan JA. BamView: visualizing and interpretation of next-generation sequencing read alignments. Brief Bioinformatics. 2013;14:203–12.

    CAS  PubMed  Google Scholar 

  20. 20.

    Wang W, Zhang X, Xu Y, Weinstock GM, Di Bisceglie AM, Fan X. High-resolution quantification of hepatitis C virus genome-wide mutation load and its correlation with the outcome of peginterferon-alpha2a and ribavirin combination therapy. PLoS One. 2014;9:e100131.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33:1870–4.

    CAS  PubMed  Google Scholar 

  22. 22.

    Pedersen BS, Quinlan AR. hts-nim: scripting high-performance genomic analyses. Bioinformatics. 2018;34:3387–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.

    CAS  PubMed  PubMed Central  Google Scholar 

  24. 24.

    SantaLucia J Jr. A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proc Natl Acad Sci U S A. 1998;95:1460–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Pan W, TTM N, Camunas-Soler J, Song CX, Kowarsky M, Blumenfeld YJ, et al. Simultaneously monitoring immune response and microbial infections during pregnancy through plasma cfRNA sequencing. Clin Chem. 2017;63:1695–704.

    CAS  PubMed  Google Scholar 

  26. 26.

    Ammar R, Torti D, Tsui K, Gebbia M, Durbic T, Bader GD, Giaever G, Nislow C. Chromatin is an ancient innovation conserved between Archaea and Eukarya. Elife. 2012;1:e00078.

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Archer SK, Shirokikh NE, Preiss T. Selective and flexible depletion of problematic sequences from RNA-seq libraries at the cDNA stage. BMC Genomics. 2014;15:401.

    PubMed  PubMed Central  Google Scholar 

  28. 28.

    Christodoulou DC, Gorham JM, Herman DS, Seidman JG. Construction of normalized RNA-seq libraries for next-generation sequencing using the crab duplex-specific nuclease. Curr Protoc Mol Biol. 2011;Chapter 4:Unit4.12.

    PubMed  Google Scholar 

  29. 29.

    Gagic D, Maclean PH, Li D, Attwood GT, Moon CD. Improving the genetic representation of rare taxa within complex microbial communities using DNA normalization methods. Mol Ecol Resour. 2015;15:464–76.

    CAS  PubMed  Google Scholar 

  30. 30.

    Bonsall D, Ansari MA, Ip CL, Trebes A, Brown A, Klenerman P, et al. ve-SEQ: Robust, unbiased enrichment for streamlined detection and whole-genome sequencing of HCV and other highly diverse pathogens. F1000Res. 2015;4:1062.

    PubMed  PubMed Central  Google Scholar 

  31. 31.

    Thomson E, Ip CL, Badhan A, Christiansen MT, Adamson W, Ansari MA. Comparison of next-generation sequencing technologies for comprehensive assessment of full-length hepatitis C viral genomes. J Clin Microbiol. 2016;54:2470–84.

    CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Nelson JR. Random-primed, Phi29 DNA polymerase-based whole genome amplification. Curr Protoc Mol Biol. 2014;105:Unit 15.13.

    PubMed  Google Scholar 

  33. 33.

    Albadalejo J, Alonso R, Antinozzi R, Bogard M, Bourgault AM, Colucci G, et al. Multicenter evaluation of the COBAS AMPLICOR HCV assay, an integrated PCR system for rapid detection of hepatitis C virus RNA in the diagnostic laboratory. J Clin Microbiol. 1998;36:862–5.

    CAS  PubMed  PubMed Central  Google Scholar 

  34. 34.

    DiDomenico N, Link H, Knobel R, Caratsch T, Weschler W, Loewy ZG, Rosenstraus M. COBAS AMPLICOR: fully automated RNA and DNA amplification and detection system for routine diagnostic PCR. Clin Chem. 1996;42:1915–23.

    CAS  PubMed  Google Scholar 

  35. 35.

    Pawlotsky JM. Use and interpretation of virological tests for hepatitis C. Hepatology. 2002;36(Suppl 1):S65–73.

    PubMed  Google Scholar 

  36. 36.

    Childs-Kean LM, Hong J. Detectable viremia at the end of treatment with direct-acting antivirals can be associated with subsequent clinical cure in patients with chronic hepatitis C: a case series. Gastroenterology. 2017;153:1165–6.

    PubMed  Google Scholar 

  37. 37.

    Maasoumy B, Buggisch P, Mauss S, Boeker KHW, Müller T, Günther R, et al. Clinical significance of detectable and quantifiable HCV RNA at the end of treatment with ledipasvir/sofosbuvir in GT1 patients. Liver Int. 2018;38:1906–10.

    CAS  PubMed  Google Scholar 

  38. 38.

    Domingo E, Sheldon J, Perales C. Viral quasispecies evolution. Microbiol Mol Biol Rev. 2012;76:159–216.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 39.

    Más A, López-Galíndez C, Cacho I, Gómez J, Martínez MA. Unfinished stories on viral quasispecies and Darwinian views of evolution. J Mol Biol. 2010;397:865–77.

    PubMed  Google Scholar 

  40. 40.

    Pybus OG, Rambaut A, Belshaw R, Freckleton RP, Drummond AJ, Holmes EC. Phylogenetic evidence for deleterious mutation load in RNA viruses and its contribution to viral evolution. Mol Biol Evol. 2007;24:845–52.

    CAS  PubMed  Google Scholar 

  41. 41.

    Ren Y, Wang W, Zhang X, Xu Y, Di Bisceglie AM, Fan X. Evidence for deleterious hepatitis C virus quasispecies mutation loads that differentiate the response patterns in IFN-based antiviral therapy. J Gen Virol. 2016;97:334–43.

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by the US National Institutes of Health (NIH) grants AI139835 (X.F.), AI117128 (X.F.), and a seed grant from the Saint Louis University Liver Center (X.F.).

Author information

Affiliations

Authors

Contributions

Conceived and designed the experiments: XF, MWF, AMD; Performed the experiments and analysis: PP, YX, XF; Data interpretation: XF, MWF, AMD; Wrote the paper: XF. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiaofeng Fan.

Ethics declarations

Ethics approval and consent to participate

All serum samples used in the current study were from a completed study at the Saint Louis University (SLU) and Hepatitis C Therapeutic Registry and Research Network (TARGET) (ClinicalTrials.gov Identifier: NCT01474811). SLU samples were previously collected for research purpose under the Saint Louis University Institutional Review Board (IRB) protocol SLU10592. Serum samples from the TARGET were provided in a code manner, and decoded for treatment outcomes at the stage of data analysis. No other patient information was obtained from the TARGET. No active patient recruitment was involved in the current study. The entire study protocol was reviewed and approved under the Saint Louis University IRB protocol SLU15565. We confirm that we have read the Journal’s position on issues involved in ethical publication and affirm that this report is consistent with those guidelines.

Consent for publication

Not applicable.

Competing interests

The authors have no conflict of interest to declare with respect to this manuscript.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1

: Figure S1. Preparation of human genome baits. Figure S2. Fragmentation of RT-tdMDA product by heating. Figure S3. Titration of Ampure XP beads in DNA purification. Table S1. Mapping statistics of 20 samples. Supplementary file 1. HCV genotype 1a consensus sequence: 9508 bp.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Peng, P., Xu, Y., Fried, M.W. et al. Genome-wide capture sequencing to detect hepatitis C virus at the end of antiviral therapy. BMC Infect Dis 20, 632 (2020). https://doi.org/10.1186/s12879-020-05355-2

Download citation

Keywords

  • Hepatitis C virus
  • Multiple displacement amplification
  • Capture sequencing
  • Duplex-specific nuclease
  • Direct acting antivirals