Figures
Abstract
Daclatasvir and asunaprevir dual oral therapy is expected to achieve high sustained virological response (SVR) rates in patients with HCV genotype 1b infection. However, presence of the NS5A-Y93H substitution at baseline has been shown to be an independent predictor of treatment failure for this regimen. By using the Invader assay, we developed a system to rapidly and accurately detect the presence of mutant strains and evaluate the proportion of patients harboring a pre-treatment Y93H mutation. This assay system, consisting of nested PCR followed by Invader reaction with well-designed primers and probes, attained a high overall assay success rate of 98.9% among a total of 702 Japanese HCV genotype 1b patients. Even in serum samples with low HCV titers, more than half of the samples could be successfully assayed. Our assay system showed a better lower detection limit of Y93H proportion than using direct sequencing, and Y93H frequencies obtained by this method correlated well with those of deep-sequencing analysis (r = 0.85, P <0.001). The proportion of the patients with the mutant strain estimated by this assay was 23.6% (164/694). Interestingly, patients with the Y93H mutant strain showed significantly lower ALT levels (p=8.8 x 10-4), higher serum HCV RNA levels (p=4.3 x 10-7), and lower HCC risk (p=6.9 x 10-3) than those with the wild type strain. Because the method is both sensitive and rapid, the NS5A-Y93H mutant strain detection system established in this study may provide important pre-treatment information valuable not only for treatment decisions but also for prediction of disease progression in HCV genotype 1b patients.
Citation: Yoshimi S, Ochi H, Murakami E, Uchida T, Kan H, Akamatsu S, et al. (2015) Rapid, Sensitive, and Accurate Evaluation of Drug Resistant Mutant (NS5A-Y93H) Strain Frequency in Genotype 1b HCV by Invader Assay. PLoS ONE 10(6): e0130022. https://doi.org/10.1371/journal.pone.0130022
Academic Editor: Ranjit Ray, Saint Louis University, UNITED STATES
Received: March 21, 2015; Accepted: May 16, 2015; Published: June 17, 2015
Copyright: © 2015 Yoshimi et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Data Availability: Newly yielded sequence data have been deposited in the NCBI Short Read Archive under BioProject PRJNA275480. All the other relevant data are within the paper.
Funding: This work was partially supported by Grants-in-Aid for scientific research and development from the Ministry of Health, Labor and Welfare (http://www.mhlw.go.jp/seisakunitsuite/bunya/hokabunya/kenkyujigyou/) and Ministry of Education, Culture, Sports, Science, and Technology (http://www.mext.go.jp/a_menu/shinkou/hojyo/main5_a5.htm), Government of Japan. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. No additional external funding was received for this study.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Hepatitis C virus (HCV) is a major cause of chronic liver disease, liver cirrhosis, and hepatocellular carcinoma, affecting up to 180 million people worldwide [1,2]. Dual oral treatment with the NS5A inhibitor daclatasvir (DCV) and the NS3 protease inhibitor asunaprevir (ASV) was one of the first interferon (IFN)-free regimens examined in treatment-experienced patients with genotype 1 HCV infection. HCV often acquires resistance against direct acting antiviral agents (DAAs) [3]. Presence of the Y93H mutation prior to treatment has been reported as an important predictor of virologic failure [4,5,6,7]. The pre-existing Y93H mutation has been estimated by direct sequencing to be present in 8.3%–19% of Japanese patients [8,9].
Direct sequencing is commonly used to detect viral mutations. However, it is only capable of detecting viral subpopulations with frequencies no lower than 10% to 20% [10,11,12,13]. In contrast, next generation sequencing (NGS), which has recently been applied to analyze viral mutations, can detect relatively low frequency variants (≤ 1%) [14,15], but it is still complex to perform and prohibitively expensive for widespread clinical use.
The Invader assay is better suited for high-throughput SNP typing [16]. To take advantage of its specificity and quantitative nature, the Invader assay has also been used for analysis of allele-specific transcription [17], detection of copy number variation [18] and drug-resistant hepatitis B virus variants [19].
In this study, we developed a rapid NS5A-Y93H strain detection system based on the Invader assay to evaluate the proportion of HCV genotype 1b patients with pre-existing Y93H mutations.
Materials and Methods
Study subjects
A total of 702 serum samples of Japanese HCV genotype 1b infected patients were screened in the study. All patients were NS5A inhibitor-treatment-naïve chronic hepatitis C patients with genotype 1b. Serum HCV RNA was measured at a central laboratory using the Roche COBAS TaqMan HCV Auto assay (Roche Diagnostics K.K., Tokyo, Japan). HCV genotype was determined at the central laboratory by polymerase chain reaction (PCR) amplification and direct sequencing. The study was approved a priori by the ethical committee of Hiroshima University and conforms to the ethical guidelines of the 1975 Declaration of Helsinki. All patients provided written informed consent.
HCV RNA extraction and cDNA synthesis
Total RNA was extracted from 150 μL of each serum sample using NucleoSpin RNA virus columns (Macherey-Nagel, Düren, Germany) according to the manufacturer's instructions. cDNA was synthesized using the PrimeScript RT reagent Kit with gDNA Eraser and oligo dT primer (TaKaRa, Otsu, Shiga, Japan).
Nested-PCR
When designing PCR primers and Invader probes, a total of 240 NS5A sequences of HCV genotype 1b from a publicly-available database [20] were utilized as a guide for successful primer and probe design. All sequences were aligned using the CLUSTALW program. The major base frequency in each nucleotide position, and thereafter, their average in each consecutive 21-bp window, was calculated and plotted (Fig 1). A higher mean major base frequency was assumed to represent lower variability at a given position and presumably improve its suitability for inclusion in a PCR primer. We used degenerate primers (Table 1), which contain alternative bases at several polymorphic sites. An example of the designing process of degenerate primer is shown in Fig 2. A fragment of 308 bp within the NS5A region was amplified from cDNA by nested PCR. The thermal profile of the initial PCR was 95°C for 2 min, followed by 35 cycles at 95°C for 15 s, at 60°C for 45 s, and at 72°C for 60 s. An aliquot of the PCR product (5 μl) was used in the nested PCR. The thermal profile for the nested PCR was the same as for the initial one, except the number of cycles was changed to 20. Amplification products were checked by agarose gel electrophoresis.
Using 240 HCV-1b NS5A sequence data from a publicly-available HCV database, mean major base frequency of any 21 bp nucleotide sequences located at every positions are plotted. A higher major base frequency was assumed to represent lower variability. Finally, two degenerate primer sets for nested PCR were designed at positions shown with bold arrows, which contained alternative bases at several polymorphic sites to increase PCR success rate.
The stacked bar represents the proportion of each base at a given nucleotide position in the NS5A region of HCV-KT9 calculated from 240 NS5A sequences of HCV genotype 1b. The solid black line represents mean major base frequency of consecutive 21 nucleotides. The Y93_F primer was placed where overall mean major base frequency can be approximately maximized, and thereafter alternative bases were included at several positions to avoid mismatch with minor base substitutions. A trade-off relationship between specificity and coverage should be considered. Basically, the same holds for the design for Invader probe set, except that there is less positional flexibility from the nature of the principle itself.
Invader assay
Three nucleotide sequences, which also contained several degenerate sites, were designed for the Invader assay (Table 1). For the calibration of the proportion of mutant variants, three types of standard nucleotide were prepared: Std-YY, Std-YH, and Std-HH (Table 1). Each sequence includes binary target sequences, corresponding to wild (Y) or mutant (H) variants, and also includes annealing sites with the internal primers at each end. An additional four nucleotides were inserted between the binary target sequences to avoid interference. Dilution series of these standards were assayed in triplicate with each assay. For Invader assays, 384-well reaction plates were used with the ABI Prism 7900 HT sequence detection system (Life Technologies, Foster City, CA, USA). The plate was incubated isothermally at 63°C and fluorescence intensities were detected every 2 min for 40 min or until non-specific fluorescence [21] was detected. Each sample was tested in triplicate in the same plate. A flow diagram representing a method of nested-PCR followed by Invader assay is depicted in Fig 3.
(A) The initial PCR was performed to amplify a fragment of 386 bp length containing a part of the NS5A region from cDNA reverse-transcribed from HCV-RNA which was extracted from the serum of a patient. (B) An aliquot of the inital PCR product was used for the second PCR to amplify a 308 bp fragment. The second PCR product was diluted with water and subjected to the Invader assay. (C) Invader oligonucleotide and allele-specific probes anneal with target to form one base overlap. When the base is complementary to the opposing base in the allele-specific probe, cleavase recognizes the structure and releases 5’ flap. The released flap anneals to a FRET probe. The second cleavage reaction releases fluorophore resulting in the generation of a fluorescent signal. The Invader assay was done in triplicate. (D) Three types of standard nucleotide (Std-YY, Std-YH, and Std-HH) were prepared, which includes binary target sequences, corresponding to wild (Y) or mutant (H) variants, and annealing sites with the internal primers at each end. Tenfold serial dilusions of each standard were subjected to the second PCR. (E) The Invader assay for each standard was also performed in triplicate.
Deep sequencing
The frequency of NS5A-Y93H was determined by ultra-deep sequencing using the Ion Torrent PGM (Life Technologies) according to previously described protocols [22] with some modifications. In brief, the fragment distributions of all amplicons were analyzed on a Bio-Analyzer 2100 (Agilent Technologies, Palo Alto, CA). The amplified fragments were modified by the Ion Xpress Plus Fragment Library Kit, and sequence analysis was performed using Ion PGM. Read mapping onto the HCV-KT9 reference sequence (GenBank accession no. AB435162) was performed using the CLC Genomics Workbench software (CLC bio, Aarhus, Denmark). This technique revealed an average coverage depth of over 10,000 reads per position in the unique regions of the genome. A wild type hepatitis C virus-expressing plasmid pHCV-KT9 [23] was used to estimate the minimum detection threshold of Y93H frequency. The sequence data have been deposited in the NCBI Short Read Archive under BioProject PRJNA275480.
Results
Estimation of Y93H mutant proportion by standard oligonucleotides
After nested PCR, the products were diluted up to 20-fold and used as templates for the Invader assay. Fluorescence intensities were measured intermittently up to 40 minutes of incubation or until emergence of non-specific fluorescence. A two-dimensional plot was made based on the relative fluorescence signal intensities of the two probes (Yakima-Yellow/Rox and FAM/Rox). As shown in Fig 4A, by using a crossing angle θsample between the Std-YY and given sample, the Y93H mutant proportion can be calculated with Equation 1, which is represented by the nearly linear curve of θsample vs Y93H% in Fig 4B. Notably, the Std-YH did not always head in the intermediate direction between Std-YY and Std-HH (Fig 4C), showing that the same intensities of the two kinds of fluorophore signals do not necessarily correspond to an equivalent amount of Y93 wild and Y93H mutant strains. By using a crossing angle θStd-YH between the Std-YY and the Std-YH, A corrected Y93H mutant proportion, Y93H%corr, can be derived from Y93H% with Equation 2 as shown in Fig 4D. Estimation was carried out using data at an appropriate time point before the emergence of non-specific fluorescence.
(A) Observed fluorescence intensities of samples and standards were plotted on a two-dimensional graph, where the X-axis was Yakima-Yellow intensity and the Y-axis was FAM intensity, corresponding to Y93 wild and Y93H mutant strains, respectively. By using crossing angle θsample, the Y93H mutant proportion of a sample can be calculated with Equation 1. (B) The resulting curve of θsample vs Y93H% is nearly linear. (C) Although the PCR product from the Std-YH theoretically contains equal amounts of Y93 wild type and Y93H mutant type target sequences, the plot of the Std-YH did not always head in the intermediate direction between Std-YY and Std-HH. (D) Therefore, by using a crossing angle θStd-YH, The corrected Y93H mutant proportion Y93H%corr can be calculated from Y93H% with Equation 2. Estimation was carried out using data at an appropriate time point before the emergence of non-specific fluorescence.
HCV RNA level and assay success rate
Next, the assay success rate and the lower detection limit in HCV RNA level were evaluated. Fig 5 shows histograms of successful samples and unsuccessful samples based on their HCV titer assayed by a quantitative commercial RT-PCR assay. The overall success rate was 98.9% (694/702). Even in serum samples with low HCV titers, more than half of samples could be successfully assayed. Thus, the lower detection limit of HCV RNA level could be estimated to be almost 1 log order. Among 694 patients successfully tested, the mutant strain (Y93H) was observed in 23.6% (164/694) patients.
White and black bars represent successful and unsuccessful detection cases, respectively. Closed circles represent the success rate for each log level of viremia from 1 through 7.
Sensitivity of the Invader assay for NS5A-Y93H strain detection
Sensitivity for NS5A-Y93H strain detection was evaluated in comparison with the deep-sequencing data of 55 sera of HCV 1b patients. These were selected to include various proportions of NS5A-Y93H strain by this assay. As shown in Fig 6, among samples containing more than 2.8% mutant, detection of Y93H mutant variant by the Invader assay was consistently successful, and the lower detection limit of the proportion of Y93H strains could be estimated at 1 to 2%.
Deep-sequencing was performed using 55 sera of HCV 1b patients with various Y93H frequencies by the Invader assay. White and black bars represent successful detection cases and detection failure cases, respectively. Closed circles represent the success rate.
Correlation for measurement of Y93H proportion
A significant positive correlation was observed between the proportion of Y93H strains determined by this assay system and the proportion determined by NGS (r = 0.85, P <0.001) (Fig 7). However, in samples in which a G allele was predominant at the preceding nucleotide position adjacent to NS5A amino acid 93 (plotted by open circles), there was a tendency to overestimate the Y93H proportion by Invader assay.
Positive correlation was observed for the Y93H strain frequencies between deep-sequencing and the Invader assay of 55 HCV 1b patients (Spearman correlation r = 0.85, p < 0.001). In samples in which G allele at the preceding nucleotide position before NS5A amino acid 93 was predominant (> 85%) measured by NGS (open circle), there was a tendency to overestimate the Y93H proportion by the Invader assay.
Association between clinical phenotypes and Y93H
We investigated possible relationships between several clinical phenotypes and Y93H strains. As shown in Table 2, patients with the Y93H mutant strain showed significantly lower ALT levels (p = 8.8 x 10–4), higher serum HCV RNA levels (p = 4.3 x 10–7), and lower HCC risk (p = 6.9 x 10–3) than those with the wild type strain.
Discussion
Several mutations in HCV structural and non-structural proteins have been reported to be associated with interferon resistance, including Core [24], NS3/NS4A, NS5A [25,26], and E2 [27]. In all classes of DAAs, rapid selection of resistant variants and viral breakthrough have also been observed [3]. Such resistant variants can arise from pre-existing viral subpopulations in the patient [28]. Several DAA combination therapies, with or without pegylated-IFN plus ribavirin, were reported to minimize the emergence of resistance and improve efficacy. Dual oral therapy with ASN and DCV was one of the first interferon-free regimens to undergo clinical evaluation and has recently been approved for use in Japan. ASN and DCV therapy results in a high rate (65%–78%) of sustained virological response (SVR) in genotype 1b patients who were null responders during prior interferon therapy; however, the Y93 mutation has been reported to be an independent predictor for non-SVR with DCV plus ASV therapy by multivariate analysis [6].
Under such circumstances, we attempted to establish a rapid, sensitive, and cost-effective method to detect the presence of the drug resistant NS5A-Y93H strain to inform treatment decisions in a practical clinical setting. HCV exists as a highly heterogeneous population of so-called quasispecies [29]. To avoid the adverse effects of hypervariability of the HCV genome on the assay, PCR-priming sites were selected based on sequence analysis of 240 HCV genotype 1b NS5A sequences from the database. To compensate for common sequence variants, several degenerate sites were included in the PCR primers and Invader probes to improve the success rate. A high overall success rate of 98.9% (694/702) was attained, which was higher than that of the earlier method for Y93H detection by real-time PCR using cycling probes [30]. In order to increase successful detection in samples with low HCV titer, nested PCR was undertaken to reach saturation in as many samples as possible. This assay protocol was oriented to evaluate the relative amount of Y93H mutant strain compared to wild strain but not to evaluate the quantity of both strains because serum HCV RNA levels can be easily quantified by commercially available assays.
Deep sequencing is increasingly being utilized to detect low frequency drug resistant HCV variants. The reported error rate of NGS platforms ranges from a few tenths of a percent to several percent [31]. We have previously reported that the minimum Y93H variant frequency detection threshold for deep-sequencing was 0.3% and that Y93H variants could be detected in 50% (5/10) patients before DCV treatment [32]. Miura et al reported that Y93H mutations with a frequency of 1% or higher could be detected by deep-sequencing. In larger studies by deep-sequencing, higher prevalence rates (25%–30.9%) of Y93H were reported [10,33], which were higher than that of direct sequencing (8.2%–19.0%) [8,9]. In this study, our assay system achieved a better lower detection limit of Y93H, estimated to be 1 to 2%, than direct sequencing, although it is not as sensitive as deep-sequencing. The prevalence rate (23.6%) of Y93H estimated by our assay is comparable with that assayed by real-time PCR using cycling probes (19.7%)[30], and ranked between those of deep sequencing and direct sequencing reported in Japanese population, presumably reflecting the degree of lower detection limit of Y93H. Natural prevalence of NS5A Y93H in Japanese population seems to be comparable with those in Western and South American countries examined by direct sequencing [34,35,36].
There were eight samples which failed in this assay, including four of HCV RNA levels within detection limit of the assay. Amplified products were not detected in all these samples by gel electrophoresis. Then we performed nested PCR using two sets of primers designed within a highly conserved 5' untranslated region [37]. Successful amplification was achieved in three samples of relatively high viral titer. Redesign of the primers may improve success rate. In addition to HCV RNA levels, the assay success may be possibly affected by degradation of viral RNA during several freeze/thawing cycles during the long-term storage.
A significant positive correlation was observed between the Invader assay and NGS with respect to Y93H proportion (Fig 7). However, a tendency of overestimation was observed in samples in which a G allele was predominant at the preceding nucleotide (nt276) in the codon adjacent to NS5A aa 93.
We also investigated possible relationships between several clinical phenotypes and Y93H and found that the existence of the Y93H mutant strain was significantly associated with ALT level, serum HCV titer, and HCC risk. Our data show similarities (HCV RNA) and differences (age, ALT, and platelet counts) when compared with those reported by Itakura et al [8], although they were analyzed by direct sequencing. Interestingly, association between Y93H and favorable IL28B/IFNL4 genotype has been recently reported by our group and others in Japan [8,33,38]. IL28B/IFNL4 variant influences viral response to interferon-based HCV therapy as well as persistent HCV infection [39,40]. These findings suggest that Y93H mutation not only confer resistance to NS5A inhibitors but also interact with genes involved in innate immunity and pro-inflammatory response such as IL28B [8,33,38]. Further research with larger samples and functional analysis is needed.
In conclusion, we developed a rapid detection system for estimating NS5A-Y93H strain frequency. This system attained high assay success rates, and sensitive detection of Y93H compared to direct sequencing. The information of NS5A-Y93H strain may provide important information not only for treatment decisions but also for prediction of disease progression in HCV genotype 1b patients.
Acknowledgments
The authors thank the subjects who agreed to participate in this study. We also thank the team members at Hiroshima University Hospital and Hiroshima Liver Study Group for clinical sample collection. The authors thank Kana Izumoto and Tomoko Shibata for technical assistance, Junko Sakamiya for technical and clerical assistance.
Author Contributions
Conceived and designed the experiments: SY HO DM MI KC. Performed the experiments: SY HO EM TU HK H.Abe NH. Analyzed the data: SY HO EM TU HK SA NH CNH. Contributed reagents/materials/analysis tools: HO DM CNH MI H. Aikata KC. Wrote the paper: SY HO CNH KC.
References
- 1. Shepard CW, Finelli L, Alter MJ Global epidemiology of hepatitis C virus infection. Lancet Infect Dis 2005; 5: 558–567. pmid:16122679
- 2. Chevaliez S, Pawlotsky JM Hepatitis C virus: virology, diagnosis and management of antiviral therapy. World J Gastroenterol 2007; 13: 2461–2466. pmid:17552030
- 3. Sarrazin C, Zeuzem S Resistance to direct antiviral agents in patients with hepatitis C virus infection. Gastroenterology 2010; 138: 447–462. pmid:20006612
- 4. Karino Y, Toyota J, Ikeda K, Suzuki F, Chayama K, Kawakami Y, et al. Characterization of virologic escape in hepatitis C virus genotype 1b patients treated with the direct-acting antivirals daclatasvir and asunaprevir. J Hepatol 2013; 58: 646–654. pmid:23178977
- 5. Kumada H, Suzuki Y, Ikeda K, Toyota J, Karino Y, Chayama K, et al. Daclatasvir plus asunaprevir for chronic HCV genotype 1b infection. Hepatology 2014; 59: 2083–2091. pmid:24604476
- 6. Manns M, Pol S, Jacobson IM, Marcellin P, Gordon SC, Peng CY, et al. All-oral daclatasvir plus asunaprevir for hepatitis C virus genotype 1b: a multinational, phase 3, multicohort study. Lancet 2014; 384: 1597–1605. pmid:25078304
- 7. Suzuki Y, Ikeda K, Suzuki F, Toyota J, Karino Y, Chayama K, et al. Dual oral therapy with daclatasvir and asunaprevir for patients with HCV genotype 1b infection and limited treatment options. J Hepatol 2013; 58: 655–662. pmid:23183526
- 8. Itakura J, Kurosaki M, Takada H, Nakakuki N, Matsuda S, Gondou K, et al. Naturally occurring, resistance-associated hepatitis C virus NS5A variants are linked to IL28B genotype and are sensitive to interferon-based therapy. Hepatol Res 2015.
- 9. Suzuki F, Sezaki H, Akuta N, Suzuki Y, Seko Y, Kawamura Y, et al. Prevalence of hepatitis C virus variants resistant to NS3 protease inhibitors or the NS5A inhibitor (BMS-790052) in hepatitis patients with genotype 1b. J Clin Virol 2012; 54: 352–354. pmid:22658798
- 10. Hirotsu Y, Kanda T, Matsumura H, Moriyama M, Yokosuka O, Omata M HCV NS5A resistance-associated variants in a group of real-world Japanese patients chronically infected with HCV genotype 1b. Hepatol Int 2015.
- 11. Pallier C, Castera L, Soulier A, Hezode C, Nordmann P, Dhumeaux D, et al. Dynamics of hepatitis B virus resistance to lamivudine. J Virol 2006; 80: 643–653. pmid:16378967
- 12. Palmer S, Kearney M, Maldarelli F, Halvas EK, Bixby CJ, Bazmi H, et al. Multiple, linked human immunodeficiency virus type 1 drug resistance mutations in treatment-experienced patients are missed by standard genotype analysis. J Clin Microbiol 2005; 43: 406–413. pmid:15635002
- 13. Rehermann B, Fowler P, Sidney J, Person J, Redeker A, Brown M, et al. The cytotoxic T lymphocyte response to multiple hepatitis B virus polymerase epitopes during and after acute viral hepatitis. J Exp Med 1995; 181: 1047–1058. pmid:7532675
- 14. Caraballo Cortes K, Zagordi O, Laskus T, Ploski R, Bukowska-Osko I, Pawelczyk A, et al. Ultradeep pyrosequencing of hepatitis C virus hypervariable region 1 in quasispecies analysis. Biomed Res Int 2013; 2013: 626083. pmid:23710455
- 15. Zagordi O, Klein R, Daumer M, Beerenwinkel N Error correction of next-generation sequencing data and reliable estimation of HIV quasispecies. Nucleic Acids Res 2010; 38: 7400–7409. pmid:20671025
- 16. Ohnishi Y, Tanaka T, Ozaki K, Yamada R, Suzuki H, Nakamura Y A high-throughput SNP typing system for genome-wide association studies. J Hum Genet 2001; 46: 471–477. pmid:11501945
- 17. Kamiyama M, Kobayashi M, Araki S, Iida A, Tsunoda T, Kawai K, et al. Polymorphisms in the 3' UTR in the neurocalcin delta gene affect mRNA stability, and confer susceptibility to diabetic nephropathy. Hum Genet 2007; 122: 397–407. pmid:17671797
- 18. Hosono N, Kubo M, Tsuchiya Y, Sato H, Kitamoto T, Saito S, et al. Multiplex PCR-based real-time invader assay (mPCR-RETINA): a novel SNP-based method for detecting allelic asymmetries within copy number variation regions. Hum Mutat 2008; 29: 182–189. pmid:17879341
- 19. Tadokoro K, Suzuki F, Kobayashi M, Yamaguchi T, Nagano M, Egashira T, et al. Rapid detection of drug-resistant mutations in hepatitis B virus by the PCR-Invader assay. J Virol Methods 2011; 171: 67–73. pmid:20950650
- 20. Shin IT, Tanaka Y, Tateno Y, Mizokami M Development and public release of a comprehensive hepatitis virus database. Hepatol Res 2008; 38: 234–243. pmid:17877727
- 21. Nakahara H, Sekiguchi K, Hosono N, Kubo M, Takahashi A, Nakamura Y, et al. Criterion values for multiplex SNP genotyping by the invader assay. Forensic Sci Int Genet 2010; 4: 130–136. pmid:20129472
- 22. Akuta N, Suzuki F, Sezaki H, Suzuki Y, Hosaka T, Kobayashi M, et al. Evolution of simeprevir-resistant variants over time by ultra-deep sequencing in HCV genotype 1b. J Med Virol 2014; 86: 1314–1322. pmid:24797623
- 23. Hiraga N, Imamura M, Abe H, Hayes CN, Kono T, Onishi M, et al. Rapid emergence of telaprevir resistant hepatitis C virus strain from wildtype clone in vivo. Hepatology 2011; 54: 781–788. pmid:21626527
- 24. Akuta N, Suzuki F, Sezaki H, Suzuki Y, Hosaka T, Someya T, et al. Association of amino acid substitution pattern in core protein of hepatitis C virus genotype 1b high viral load and non-virological response to interferon-ribavirin combination therapy. Intervirology 2005; 48: 372–380. pmid:16024941
- 25. El-Shamy A, Sasayama M, Nagano-Fujii M, Sasase N, Imoto S, Kim SR, et al. Prediction of efficient virological response to pegylated interferon/ribavirin combination therapy by NS5A sequences of hepatitis C virus and anti-NS5A antibodies in pre-treatment sera. Microbiol Immunol 2007; 51: 471–482. pmid:17446688
- 26. Enomoto N, Sakuma I, Asahina Y, Kurosaki M, Murakami T, Yamamoto C, et al. Comparison of full-length sequences of interferon-sensitive and resistant hepatitis C virus 1b. Sensitivity to interferon is conferred by amino acid substitutions in the NS5A region. J Clin Invest 1995; 96: 224–230. pmid:7542279
- 27. Sarrazin C, Kornetzky I, Ruster B, Lee JH, Kronenberger B, Bruch K, et al. Mutations within the E2 and NS5A protein in patients infected with hepatitis C virus type 3a and correlation with treatment response. Hepatology 2000; 31: 1360–1370. pmid:10827164
- 28. Susser S, Welsch C, Wang Y, Zettler M, Domingues FS, Karey U, et al. Characterization of resistance to the protease inhibitor boceprevir in hepatitis C virus-infected patients. Hepatology 2009; 50: 1709–1718. pmid:19787809
- 29. Martell M, Esteban JI, Quer J, Genesca J, Weiner A, Esteban R, et al. Hepatitis C virus (HCV) circulates as a population of different but closely related genomes: quasispecies nature of HCV genome distribution. J Virol 1992; 66: 3225–3229. pmid:1313927
- 30. Uchida Y, Kouyama J, Naiki K, Mochida S A Novel Simple Assay System to Quantify the Percent HCV-RNA Levels of NS5A Y93H Mutant Strains and Y93 Wild-Type Strains Relative to the Total HCV-RNA Levels to Determine the Indication for Antiviral Therapy with NS5A Inhibitors. PLoS One 2014; 9: e112647. pmid:25397971
- 31. Barzon L, Lavezzo E, Militello V, Toppo S, Palu G Applications of next-generation sequencing technologies to diagnostic virology. Int J Mol Sci 2011; 12: 7861–7884. pmid:22174638
- 32. Kosaka K, Imamura M, Hayes CN, Abe H, Hiraga N, Yoshimi S, et al. Emergence of resistant variants detected by ultra-deep sequencing after asunaprevir and daclatasvir combination therapy in patients infected with hepatitis C virus genotype 1. J Viral Hepat 2014.
- 33. Miura M, Maekawa S, Sato M, Komatsu N, Tatsumi A, Takano S, et al. Deep sequencing analysis of variants resistant to the non-structural 5A inhibitor daclatasvir in patients with genotype 1b hepatitis C virus infection. Hepatol Res 2014; 44: E360–367. pmid:24612030
- 34. Paolucci S, Fiorina L, Mariani B, Gulminetti R, Novati S, Barbarini G, et al. Naturally occurring resistance mutations to inhibitors of HCV NS5A region and NS5B polymerase in DAA treatment-naive patients. Virol J 2013; 10: 355. pmid:24341898
- 35. Peres-da-Silva A, de Almeida AJ, Lampe E NS5A inhibitor resistance-associated polymorphisms in Brazilian treatment-naive patients infected with genotype 1 hepatitis C virus. J Antimicrob Chemother 2015; 70: 726–730. pmid:25414201
- 36. Plaza Z, Soriano V, Vispo E, del Mar Gonzalez M, Barreiro P, Seclen E, et al. Prevalence of natural polymorphisms at the HCV NS5A gene associated with resistance to daclatasvir, an NS5A inhibitor. Antivir Ther 2012; 17: 921–926. pmid:22436385
- 37. Nagaoki Y, Hyogo H, Aikata H, Tanaka M, Naeshiro N, Nakahara T, et al. Recent trend of clinical features in patients with hepatocellular carcinoma. Hepatol Res 2012; 42: 368–375. pmid:22151896
- 38. Akamatsu S, Nelson Hayes C, Ochi H, Uchida T, Kan H, Murakami E, et al. Association between variants in the interferon lambda 4 locus and substitutions in the hepatitis C virus non-structural protein 5A (130 out of 130 characters). J Hepatol 2015.
- 39. Ge D, Fellay J, Thompson AJ, Simon JS, Shianna KV, Urban TJ, et al. Genetic variation in IL28B predicts hepatitis C treatment-induced viral clearance. Nature 2009; 461: 399–401. pmid:19684573
- 40. Thomas DL, Thio CL, Martin MP, Qi Y, Ge D, O’hUigin C, et al. Genetic variation in IL28B and spontaneous clearance of hepatitis C virus. Nature 2009; 461: 798–801. pmid:19759533