Disease Markers
Volume 2018 (2018), Article ID 2514901, 9 pages
https://doi.org/10.1155/2018/2514901

Research Article

Pretreatment Hepatitis C Virus NS5A/NS5B Resistance-Associated Substitutions in Genotype 1 Uruguayan Infected Patients

Fabián Aldunate,¹ Natalia Echeverría,¹ Daniela Chiodi,² Pablo López,² Adriana Sánchez-Cicerón,² Alvaro Fajardo,¹ Martín Soñora,¹ Juan Cristina,¹ Nelia Hernández,² and Pilar Moreno¹

¹Laboratorio de Virología Molecular, Centro de Investigaciones Nucleares, Facultad de Ciencias, Universidad de la República, 11400 Montevideo, Uruguay
²Clínica de Gastroenterología, Hospital de Clínicas, Facultad de Medicina, Universidad de la República, 11600 Montevideo, Uruguay

Correspondence should be addressed to Pilar Moreno; pmoreno@cin.edu.uy

Received 16 March 2018; Revised 30 May 2018; Accepted 10 June 2018; Published 14 August 2018

Academic Editor: Nora A. Fierro

Copyright © 2018 Fabián Aldunate et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Hepatitis C Virus (HCV) infection treatment has dramatically changed with the advent of direct-acting antiviral agents (DAAs). However, the efficacy of DAAs can be attenuated by the presence of resistance-associated substitutions (RASs) before and after treatment. Indeed, RASs detected in DAA treatment-naïve HCV-infected patients could be useful for clinical management and outcome prediction. Although the frequency of naturally occurring HCV NS5A and NS5B RASs has been addressed in many countries, there are only a few reports on their prevalence in the South American region. The aim of this study was to investigate the presence of RASs to NS5A and NS5B inhibitors in a DAA treatment naïve cohort of Uruguayan patients infected with chronic hepatitis C and compare them with reports from other South American countries. Here, we found that naturally occurring substitutions conferring resistance to NS5A and NS5B inhibitors were present in 8% and 19.2%, respectively, of treatment-naïve HCV genotype 1 infected patients. Importantly, the baseline substitutions in NS5A and NS5B herein identified differ from the studies previously reported in Brazil. Furthermore, Uruguayan strains subtype 1a clustered within all major world clades, showing that HCV variants currently circulating in this country are characterized by a remarkable genetic diversity.

1. Introduction

Hepatitis C Virus (HCV) infection treatment has dramatically improved thanks to the introduction of direct-acting antiviral agents (DAAs). These antivirals have significantly increased response rates (up to 98%) and greatly reduced treatment duration [1]. Currently available DAAs are classified into four categories given their molecular targets in the HCV replication cycle: (1) NS3/4A protease inhibitors (PIs) bind to the active site of the NS3/4A protease; (2) NS5A inhibitors interact with domain 1 of the NS5A dimer, although the exact mechanism of NS5A inhibition remains to be fully elucidated; (3) nucleos(t)ide analog NS5B polymerase inhibitors are incorporated into the nascent RNA chain resulting in chain termination by compromising the binding of the incoming nucleotide; (4) nonnucleoside NS5B polymerase inhibitors interact with either the thumb 1, thumb 2, palm 1, or palm 2 domain of NS5B and inhibit polymerase activity by allosteric mechanisms [2–4]. However, the extreme mutation and high replication rates of HCV, together with the immune system pressure, lead to a remarkable genetic variability that can compromise the high response rates to DAAs due to the preexistence of resistance-associated substitutions (RASs) [5, 6].

Each drug or class of DAA is characterized by specific resistance profiles. The likelihood that a DAA will select for and allow outgrowth of viral populations carrying RASs depends on the DAA’s genetic barrier to resistance (the number and type of mutations needed to generate an amino acid substitution that confers resistance), the viral fitness (replicative capacity) of the resistant variant, and viral genotypes and subtypes [7, 8].

The prevalence of RASs in treatment-naïve patients has been broadly reported worldwide [9–16]. However, apart from Brazil and Argentina, this issue has not been fully addressed in other South American countries yet [9, 17–19]. The lack of information in relation to preexisting baseline RASs, added to the high cost of these new drugs, are the major limiting factors for the broad implementation of these new therapies in Uruguay as well as in other Latin American countries (low- or lower-middle income) [20].

In this study, we explored the presence of resistance variants to NS5A and NS5B inhibitors in a DAA treatment naïve cohort of Uruguayan patients chronically infected with hepatitis C. Here, we aimed to contribute to the knowledge of the circulation of HCV resistant variants in the South American region.

2. Materials and Methods

2.1. Patients and Clinical Samples

Serum samples were obtained from 31 patients with serological markers for HCV, which were recruited between 2015 and 2017 at the Gastroenterology Clinic from Hospital de Clínicas, Montevideo, Uruguay. HCV infection was confirmed by Abbott realtime HCV (Abbott Molecular Inc., Des Plaines, USA). Patients selected for this study were both chronically infected with HCV genotype 1 and DAA treatment-naïve at the time of blood extraction. Written informed consent was obtained from all patients. The studies have been performed according to the World Medical Association Declaration of Helsinki and approved by the appropriate institutional board (Hospital de Clínicas ethical committee).

2.2. RNA Extraction, cDNA Synthesis, and NS5A and NS5B Amplification

Viral RNA was extracted from 140 μl of serum using the QIAamp Viral RNA mini kit (QIAgen, Hilden, Germany) according to the manufacturer’s protocol. The viral RNA was heated at 65°C for 5 min and used as a template for a reverse transcription reaction. The reverse transcription reaction mixture contained 5 μl of the RNA template, 1 μl of random hexamer 100 ng/μl (Invitrogen Life Technologies, Carlsbad, CA, USA), 1 μl of dNTP mix (10 mM each), 4 μl of 5X first-strand buffer, 2 μl of 0.1 M DTT, 1 μl of SuperScript II reverse transcriptase (200 U/μl) (Invitrogen Life Technologies, Carlsbad, CA, USA), and 1 μl (40 U/μl) RNaseOUT (Invitrogen Life Technologies, Carlsbad, CA, USA). The reverse transcription was performed at 42°C for 50 min, and then the reverse transcriptase enzyme was inactivated at 70°C for 15 min. PCR amplification of NS5A and NS5B genome regions was performed using primers and conditions previously described [10]. Amplicons were purified using the Illustra GFX PCR DNA and Gel Band Purification Kit (GE Healthcare Life Science, Buckinghamshire, UK) according to the manufacturer’s protocol.

2.3. NS5A and NS5B Sequencing

The purified product was then sequenced using the same sets of primers used for PCR amplification. Bidirectional Sanger sequencing was performed by Macrogen Korea (http://www.macrogen.com).

2.4. NS5A and NS5B Genotype Determination

HCV NS5A and NS5B consensus sequences obtained from Uruguayan patients were aligned with sequences from HCV representing all genotypes and main subtypes isolated in different geographic regions of the world. These sequences were obtained from Los Alamos HCV sequence database and from the NIAID Virus Pathogen Database and Analysis Resource (ViPR) [21, 22]. For strains included in these studies, see Supplementary Material Table S1. Sequences were aligned using the CLUSTAL W software [23]. Once aligned, the best evolutionary model that described our sequence data was assessed using ModelGenerator program [24]. Using the GTR + G + I model (General time reversible + gamma + invariant sites), maximum likelihood phylogenetic trees were constructed for both NS5A and NS5B using the MEGA 5.0 software [25]. For NS5A, 953 nucleotides (positions 6367 to 7319, relative to HCV 1a reference strain, H77 NC_004102) were included in the phylogenetic analysis, whereas for NS5B, only 361 nucleotides corresponding to the Okamoto region (positions 8265 to 8625, relative to strain H77 NC_004102) were included. As a measure of the robustness of each node, we employed the bootstrapping method (1000 pseudoreplicates).

For NS5A 1a Uruguayan sequences (), a second alignment and maximum likelihood phylogenetic tree was generated in order to analyze HCV evolutionary relationships between Uruguayan, Brazilian, and worldwide strains. For non-Uruguayan strains included in this analysis, see Supplementary Material Table S2.

2.5. NS5A and NS5B Sequence Analysis

In order to properly identify substitution changes in NS5A and NS5B regions from HCV strains circulating in Uruguayan patients, we generated world consensus sequences for 1a and 1b subtypes using a wide range of NS5A and NS5B sequences from HCV strains isolated worldwide. For this purpose, NS5A gene sequences corresponding to subtypes 1a () and 1b () were retrieved from Los Alamos HCV sequence database and from the NIAID ViPR [21, 22]. Likewise, datasets of 150 and 124 NS5B sequences were generated for subtypes 1a and 1b, respectively. Using Seqman program, implemented in DNAStar 5.01 package (DNASTAR, Madison, USA), a world consensus nucleotide sequences were generated for each gene and subtype. Each Uruguayan sequence was subsequently aligned to the corresponding reference sequences, and then in silico translated. The amino acid sequences obtained were compared in order to explore the presence of RASs as well as the presence of polymorphisms at a RAS position (RAPs) in Uruguayan HCV strains. RAPs are defined as any change from reference sequence for a specific genotype at a position associated with NS5A resistance [26].

3. Results

3.1. Genetic Variability of NS5A and NS5B Genes from HCV Strains Circulating in Uruguayan Patients

To study the genetic variability of NS5A and NS5B regions of HCV strains circulating in Uruguayan patients, sequences of these regions (accession numbers MH070029-MH070090) were aligned with corresponding sequences from 59 HCV strains isolated elsewhere, representing all genotypes and main subtypes (for strains included in these analyses, see Supplementary Material Table S1). Therefore, maximum likelihood phylogenetic trees were constructed. The results of these studies are shown in Figure 1 (Figure 1(a): NS5A and Figure 1(b): NS5B).

(a)

(b)

Figure 1: Maximum-likelihood phylogenetic trees analyses of the NS5A and NS5B genes of HCV strains circulating in Uruguay. Reference strains in the tree are shown by genotype_country_accession number. Numbers at the branches indicate bootstrap values. The bar at the bottom of the tree denotes distance. Uruguayan strains genotype 1a are shown in red and genotype 1b in blue. (a) NS5A region (953 nucleotides); (b) NS5B region (361 nucleotides).

All strains in the phylogenies were assigned according to their genotype, and each cluster was supported by very high bootstrap values for both analyzed regions. Strains isolated from Uruguayan patients () were assigned to genotype 1, 20 of which corresponded to subtype 1a and 11 to subtype 1b. The results of NS5A (Figure 1(a)) and NS5B (Figure 1(b)) phylogenetic analyses were concordant for both genomic regions in all 31 sequences, suggesting no recombination events between these regions.

To further analyze the evolutionary relationships between the Uruguayan strains and those circulating in Brazil and elsewhere, a second maximum likelihood phylogenetic tree of HCV-1a sequences of NS5A partial region was built (Figure 2). As was previously described, two distinct 1a clades (clades 1 and 2) were observed. Brazilian sequences clustered in a large group of related sequences inside clade 1 [9]. Whereas NS5A Uruguayan strains (in red) did not cluster in a particular clade, rather, they grouped dispersedly within all major world clades.

Figure 2: Maximum-likelihood phylogenetic tree analysis of the NS5A of HCV-1a strains circulating in Uruguay. The maximum-likelihood phylogenetic tree shows the evolutionary relationships between the Uruguayan, Brazilian, and worldwide HCV-1a isolates sequences. Uruguayan strains are shown in red. The bar at the bottom of the tree denotes distance. Sequence names have been removed for clarity.

3.2. NS5A Substitution Analysis

With the purpose of studying the amino acid (AA) substitutions along the NS5A protein, Uruguayan HCV AA sequences were aligned with NS5A world consensus sequences (residues 23 to 354 relative to NS5A protein sequence). AA substitutions at positions previously found to be potentially associated with resistance to NS5A inhibitors, as well as polymorphisms at a RAS position, were identified. These results are summarized in Table 1.

Table 1: Amino acid substitutions in HCV NS5A protein from DAA treatment-naïve patients.

Subtype
Isolate
RAS
RAP
Resistance to DAA
Reference

1a
003

H58P
—
HARVONI prescription

1a
020

K24Q
—
[5, 49]

1a
032

H58P
—
HARVONI Prescription

1b
016
L31 V

DCV, LDV, EBV
[5, 6]

1b
039
L31 M

DCV, LDV, EBV
[5, 47, 48]

Daclatasvir (DCV), ledipasvir (LDV), elbasvir (EBV). Harvoni prescription, https://www.gilead.com/~/media/Files/pdfs/medicines/liver-disease/harvoni/harvoni_pi.pdf.

RASs to NS5A inhibitors (L31M and L31V) were identified in 2 strains out of 25 (8%) fully sequenced samples. RAPs were found in 3 strains (subtype 1a): 2 exhibited the substitution H58P and 1 the substitution K24Q. Although these substitutions were not reported as resistant, some changes at these positions were previously described as RASs in subtype 1a, namely H58D and K24R [27, 28]. Finally, substitution E62D was found in one subtype 1a strain. This change is considered as a secondary substitution because, although it does not confer resistance by itself, when combined with a known RAS it does. In fact, it confers a higher level of resistance than the one achieved by the RAS alone [26]. In addition, several polymorphisms that have not been previously reported to be associated with a resistant phenotype were also detected (see Supplementary Material Table S3).

3.3. NS5B Substitution Analysis

In order to study substitutions along NS5B protein, Uruguayan HCV AA sequences were aligned to the NS5B world consensus sequences. Almost full-length AA sequences were obtained in 26 out of 31 analyzed strains. 23 sequences span residues 36 to 539 whereas the remaining 3 span residues 36 to 557 of NS5B protein. This issue limited our studies, since many of the described RASs are observed as of residue 553.

Importantly, RASs to NS5B inhibitors (Table 2) were observed in 5 strains out of 26 sequenced samples (19.2%). C451R was found in two isolates while A421V was found in only one. In 2 of the 3 strains for which we were able to obtain longer sequences, RASs S556G (subtype 1a) and Q556R (subtype 1b) were observed.

Table 2: Amino acids substitutions in HCV NS5B protein from DAA treatment-naïve patients.

Subtype
Isolate
RAS
RAP
Resistance to DAA
Reference

1a
003
C451R

DSV
Trial Coral I - cohort 2

1a
020
A421V

BCV
[29]

1a
022
C451R

DSV
Trial Coral I - cohort 2

1a
037
S556G

DSV
[54, 55]

037

A553G
—
[8]

1b
026

A421V
—
[29]

026
Q556R

DSV
[12]

1b
036

A421V
—
[29]

Beclabuvir (BCV), dasabuvir (DSV). Trial Coral I - cohort 2 http://www.hcv-trials.com/showStudy.asp?Study=86 (Mantry PS. AASLD, 2015, Abs. 1084).

Finally, we found two RAPs: A421V (in 2 subtype 1b strains) and A553G (in 1 subtype 1a strain). Although A421V has been associated with resistance to beclabuvir (BCV) in patients infected with HCV subtype 1a, this resistant phenotype has not been proven in strains subtype 1b [29]. In position 553, the substitution reported as resistant was A553T [8].

As was the case for NS5A, different polymorphisms not previously associated with a resistant phenotype were also detected in NS5B (see Supplementary Material Table S4).

4. Discussion

The advent of DAAs therapies constitutes one of the major breakthroughs in HCV infected patients management. However, these new treatment options are far from being universally available, in particular for HCV infected patients relying on Latin American public healthcare systems. The main limiting factors for worldwide access to DAAs in our region concern the high cost, the inadequate management of public healthcare systems, the limited access of low-income or uninsured populations to healthcare providers, and the lack of accurate epidemiological information [20, 30–32]. In Uruguay, these therapies became recently available, and although some have been approved for their use by the public health authorities (Viekira pak and sofosbuvir/ledipasvir therapies), they are not currently financially covered, except in specific cases. Despite the high rates of viral response achieved with DAA-based treatments, still 1 to10% of the patients fails to eliminate infection, and in these cases, baseline and emergent resistance variants turn out to be key factors contributing to treatment failure [5, 17, 33].

Unfortunately, we are currently unable to properly assess the number of HCV infected people in Uruguay and even more to figure out the frequency and type of RASs circulating. These facts could compromise the effectiveness of these new therapies in our country.

We have previously reported that naturally occurring substitutions conferring resistance to NS3 inhibitors exist in a significant proportion of Uruguayan patients infected with HCV genotype 1, and we showed that this frequency seemed to be higher than in other South American countries (Brazil and Argentina) [34]. The present study describes the prevalence of baseline NS5A and NS5B RASs in HCV genotype 1 infected DAA-naïve patients in a Uruguayan cohort.

The presence of substitutions conferring resistance to NS5A inhibitors has been widely reported both in therapy-naïve and in relapser patients from Europe [10, 33, 35–38], USA [37, 39, 40], and Asia [41–43]. However, NS5A sequences from South America are poorly analyzed yet [9, 44]. Recent studies have revealed that the mean prevalence of NS5A genotype 1 baseline RASs to different inhibitors ranges from 6% to 16% using population sequencing or deep sequencing [27, 37, 45, 46]. Importantly, the prevalence and type of baseline NS5A RASs varies slightly by geographic regions. For instance, L31M was found in 2.2% of genotype 1a infected patients in Europe, in 4.1% of those in Oceania, and strikingly in no patient from the USA [27]. For this reason, we believe that there is a need to contribute data from our region, for which we still do not have enough information, apart from Brazil [9, 44]. The results of this study indicate the presence of DAA NS5A RASs in 2 HCV strains (8% of the patients enrolled in this study), with baseline RASs detected at position 31 (see Table 1). L31M substitution confers resistance to daclatasvir (DCV), ledipasvir (LDV), and elbasvir (EBV) in both 1a and 1b subtypes [5, 6, 8, 28, 47, 48], whereas substitution L31V does it to DCV in subtypes 1a and 1b, to LDV in subtype 1b, and to EBV in subtype 1a [5, 6, 28]. Given that both L31V and L31M are clinically relevant RASs, their detection at baseline may influence the choice of first-line treatment regimens [28].

The substitutions H58P and K24Q found in two patients are considered as resistance-associated polymorphisms (RAPs). The RASs characterized at these positions were H58D and K24G/N/R [5, 6, 27, 28, 49, 50]. The substitution H58P was found as a baseline RAP in relapsers to LDV (HARVONI prescription, https://www.gilead.com/-/media/files/pdfs/medicines/liver-disease/harvoni/harvoni_pi.pdf?la=en). However, it is sometimes regarded as a RAS [10, 51], despite conferring only 1.2 fold change in resistance in in vitro studies using the 1a replicon system [39].

We did not find M28T/V, Q30R/H, or Y93H substitutions as there were previously reported in Brazil and worldwide [9, 27, 44]. The amino acid substitution E62H was found in one Uruguayan patient. Although this change does not confer resistance by itself but in combination with Q30R, it generates a high resistance level to DCV [52].

The presence of baseline NS5A RASs impacts treatment outcome in some patient groups by affecting SVR rates. The detection of NS5A preexistent RASs may play a relevant role in the choice of first-line treatment regimens or in the simplification/shortening of recommended regimens, in order to bring SVR rates close to the highest achievable [27, 38, 41, 53], in particular in countries such as Uruguay, where only two different DAA-containing treatment regimens are approved for their use.

Regarding NS5B gene, global analysis (with the exception of South America [17, 19]) revealed that NS5B DAA resistance substitutions are infrequent [14]. Our study showed the presence of NS5B inhibitors RASs in 5 out of 26 analyzed HCV infected Uruguayan patients naïve to treatment (19.2%). Substitutions found in this work were A421V and S556G associated in subtype 1a with resistance to BCV and dasabuvir (DSV), respectively [8, 28, 29, 54, 55], and Q556R associated with resistance to DSV both in genotype 1a and 1b [12, 28]. Substitution C451R, observed in two Uruguayan patients, was reported previously in patients who failed to clear the infection after treatment with OBV/PTV/r + DSV ± RBV. In these cases, it appeared in combination with G558R (Trial Coral I-Cohort 2: http://www.hcv-trials.com/showStudy.asp?Study=86).

RAPs in positions 421 and 553 (A421V in two subtype 1b isolates and A553G in one subtype 1b isolate) were also found. Although A421V has been associated with resistance to BCV in patients with subtype 1a, this phenotype has not been proven in strains of subtype 1b [29]. In position 553, the substitutions reported as resistant are A553T in subtype 1a [8] and A553V in subtype 1b [54], conferring resistance to DSV.

In contrast to our results, Noble and coworkers (2016) reported the presence of V321A, A421G, M414V, Y448H, L159F, and C316N in Brazilian isolates [17], yet none of these mutations were found in this study, probably due to the diversity found between Uruguayan and Brazilian strains (Figure 2). Nevertheless, substitution A421V was found in Brazil [17], Argentina [19], and Uruguay. The RAS S282T was detected neither in Brazilian reports nor in this current work (Uruguay) [17, 18, 56]. Our findings further confirm and complement previous studies which evidenced a low prevalence of this substitution in vivo, probably due to its low replicative fitness [14, 18, 57]. Despite our results, it is worth mentioning that the presence of baseline NS5B RASs conferring resistance to nucleotide or nonnucleoside NS5B inhibitors has not been shown to have any impact on virologic responses thus far [53, 58].

These results show both diversity in the baseline polymorphisms found in different Latin American countries and in the evolutionary relationships of Uruguayan isolates (Figure 2). This fact could be linked not only to the isolates’ geographic region and viral intrinsic characteristics but also to the genetic background of the host. It is worth mentioning that we live in a vast continent inhabited by populations with different genotypic characteristics that might, depending on the situation, require different approaches to treatment. Indeed, we have recently found that allele and genotype frequencies at IL28B locus of Uruguayan individuals closely resemble those of an admixed population rather than a uniformly European-descendant one [59]. Altogether, we believe that it could be important to carry out studies throughout the South American region in order to establish the prevalence of RASs in NS5A and NS5B in different countries. In fact, this will aid in understanding that not every treatment regimen might be adequate for every patient and country. The data we presented here might guide not only physicians in making therapeutic decisions but also public health authorities in approving more diverse treatment combinations. These treatment formulations would cover most of the circulating strains in our region, a region with an extremely diverse genetic background population.

5. Conclusion

To our knowledge, the present study revealed for the first time the presence of RASs in the NS5A and NS5B regions of HCV genotype 1 Uruguayan strains from patients who have not been previously treated with DAAs and is one of the few South American countries to report on this matter. It is currently unclear if preexisting viral variants with reduced susceptibility to DAAs are clinically relevant for the prediction of virologic treatment failure. However, individualized DAA therapy based on baseline resistance analysis may be beneficial for optimizing treatment efficacy in patients with HCV genotype 1 infection and risk factors for treatment failure. Therefore, the potential role of baseline resistance testing remains an area of critical research and clinical questions.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Fabián Aldunate and Natalia Echeverría contributed equally to this work.

Acknowledgments

The authors thank Dr. Gonzalo Moratorio (Viral populations and Pathogenesis Unit, Institut Pasteur, Paris) for assistance during the manuscript preparation. This work was supported by Agencia Nacional de Investigación e Innovación (ANII) (project FMV_1_2014_1_104171), PEDECIBA, Comisión Académica de Posgrados, Universidad de la República Uruguay (UdelaR), and Comisión Sectorial de Investigación Científica (CSIC), UdelaR, Uruguay.

Supplementary Materials

Supplementary Material Table S1: hepatitis C Virus NS5A and NS5B sequences used as representatives of each genotype to perform the phylogenetic analysis. Their corresponding genotype, country of isolation, and GenBank accession number are indicated. Supplementary Material Table S2: hepatitis C Virus NS5A subtype 1a sequences used to reveal evolutionary relationships between Uruguayan strains and others isolated elsewhere. Their corresponding country of isolation and GenBank accession number are indicated. Supplementary Material Table S3: amino acid substitutions in NS5A protein not previously associated with resistance to NS5A inhibitors. Supplementary Material Table S4: amino acid substitutions in NS5B protein not previously associated with resistance to polymerase inhibitors. (Supplementary Materials)

References

T. Asselah, P. Marcellin, and R. F. Schinazi, “Treatment of hepatitis C virus infection with direct-acting antiviral agents: 100% cure?” Liver International, vol. 38, pp. 7–13, 2018.

J. M. Pawlotsky, “NS5A inhibitors in the treatment of hepatitis C,” Journal of Hepatology, vol. 59, no. 2, pp. 375–382, 2013.

L. Gerber, T. M. Welzel, and S. Zeuzem, “New therapeutic strategies in HCV: polymerase inhibitors,” Liver International, vol. 33, Supplement 1, pp. 85–92, 2013.

P. de Leuw and C. Stephan, “Protease inhibitor therapy for hepatitis C virus-infection,” Expert Opinion on Pharmacotherapy, vol. 19, no. 6, pp. 577–587, 2018.

D. L. Wyles, “Resistance to DAAs: when to look and when it matters,” Current HIV/AIDS Reports, vol. 14, no. 6, pp. 229–237, 2017.

I. Esposito, J. Trinks, and V. Soriano, “Hepatitis C virus resistance to the new direct-acting antivirals,” Expert Opinion on Drug Metabolism & Toxicology, vol. 12, no. 10, pp. 1197–1209, 2016.

L. Cuypers, G. Li, P. Libin, S. Piampongsant, A. M. Vandamme, and K. Theys, “Genetic diversity and selective pressure in hepatitis C virus genotypes 1–6: significance for direct-acting antiviral treatment and drug resistance,” Viruses, vol. 7, no. 9, pp. 5018–5039, 2015.

E. Lontok, P. Harrington, A. Howe et al., “Hepatitis C virus drug resistance-associated substitutions: state of the art summary,” Hepatology, vol. 62, no. 5, pp. 1623–1632, 2015.

A. Peres-da-Silva, A. J. De Almeida, and E. Lampe, “NS5A inhibitor resistance-associated polymorphisms in Brazilian treatment-naive patients infected with genotype 1 hepatitis C virus,” Journal of Antimicrobial Chemotherapy, vol. 70, no. 3, pp. 726–730, 2015.

S. Paolucci, L. Fiorina, B. Mariani et al., “Naturally occurring resistance mutations to inhibitors of HCV NS5A region and NS5B polymerase in DAA treatment-naïve patients,” Virology Journal, vol. 10, no. 1, p. 355, 2013.

Y. Zhang, Y. Cao, R. Zhang et al., “Pre-existing HCV variants resistant to DAAs and their sensitivity to PegIFN/RBV in Chinese HCV genotype 1b patients,” PLoS One, vol. 11, no. 11, article e0165658, 2016.

Z. Chen, H. Li, H. Ren, and P. Hu, “Global prevalence of pre-existing HCV variants resistant to direct-acting antiviral agents (DAAs): mining the GenBank HCV genome data,” Scientific Reports, vol. 6, no. 1, p. 20310, 2016.

T. L. Applegate, S. Gaudieri, A. Plauzolles et al., “Naturally occurring dominant drug resistance mutations occur infrequently in the setting of recently acquired hepatitis C,” Antiviral Therapy, vol. 20, no. 2, pp. 199–208, 2015.

T. M. Welzel, N. Bhardwaj, C. Hedskog et al., “Global epidemiology of HCV subtypes and resistance-associated substitutions evaluated by sequencing-based subtype analyses,” Journal of Hepatology, vol. 67, no. 2, pp. 224–236, 2017.

J. Itakura, M. Kurosaki, H. Takada et al., “Naturally occurring, resistance-associated hepatitis C virus NS5A variants are linked to interleukin-28B genotype and are sensitive to interferon-based therapy,” Hepatology Research, vol. 45, no. 10, pp. E115–E121, 2015.

J. A. Larousse, P. Trimoulet, P. R. Pinson et al., “Prevalence of hepatitis C virus (HCV) variants resistant to NS5A inhibitors in naïve patients infected with HCV genotype 1 in Tunisia,” Virology Journal, vol. 12, no. 1, p. 84, 2015.

C. F. Noble, F. Malta, G. Lisboa-Neto et al., “Natural occurrence of NS5B inhibitor resistance-associated variants in Brazilian patients infected with HCV or HCV and HIV,” Archives of Virology, vol. 162, no. 1, pp. 165–169, 2017.

S. Franco, M. Casadellà, M. Noguera-Julian et al., “No detection of the NS5B S282T mutation in treatment-naïve genotype 1 HCV/HIV-1 coinfected patients using deep sequencing,” Journal of Clinical Virology, vol. 58, no. 4, pp. 726–729, 2013.

M. M. Sede, N. L. Laufer, and J. Quarleri, “Previous failure of interferon-based therapy does not alter the frequency of HCV NS3 protease or NS5B polymerase inhibitor resistance-associated variants: longitudinal analysis in HCV/HIV co-infected patients,” International Journal of Antimicrobial Agents, vol. 46, no. 2, pp. 219–224, 2015.

J. Trinks, M. Caputo, M. L. Hulaniuk, D. Corach, and D. Flichman, “Hepatitis C virus pharmacogenomics in latin american populations: implications in the era of direct-acting antivirals,” Pharmacogenomics and Personalized Medicine, vol. 10, pp. 79–91, 2017.

C. Kuiken, K. Yusim, L. Boykin, and R. Richardson, “The Los Alamos hepatitis C sequence database,” Bioinformatics, vol. 21, no. 3, pp. 379–384, 2005.

B. Pickett, D. Greer, Y. Zhang et al., “Virus pathogen database and analysis resource (ViPR): a comprehensive bioinformatics database and analysis resource for the coronavirus research community,” Viruses, vol. 4, no. 11, pp. 3209–3226, 2012.

J. D. Thompson, D. G. Higgins, and T. J. Gibson, “CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice,” Nucleic Acids Research, vol. 22, no. 22, pp. 4673–4680, 1994.

T. M. Keane, C. J. Creevey, M. M. Pentony, T. J. Naughton, and J. O. Mclnerney, “Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified,” BMC Evolutionary Biology, vol. 6, no. 1, p. 29, 2006.

K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.

J.-H. Sun, D. R. O'Boyle II, Y. Zhang et al., “Impact of a baseline polymorphism on the emergence of resistance to the hepatitis C virus nonstructural protein 5a replication complex inhibitor, BMS-790052,” Hepatology, vol. 55, no. 6, pp. 1692–1699, 2012.

S. Zeuzem, M. Mizokami, S. Pianko et al., “NS5A resistance-associated substitutions in patients with genotype 1 hepatitis C virus: prevalence and effect on treatment outcome,” Journal of Hepatology, vol. 66, no. 5, pp. 910–918, 2017.

European Association for the Study of the Liver, “EASL recommendations on treatment of hepatitis C 2016,” Journal of Hepatology, vol. 66, no. 1, pp. 153–194, 2017.

S. Kim, K.-H. Han, and S. H. Ahn, “Hepatitis C virus and antiviral drug resistance,” Gut and Liver, vol. 10, no. 6, pp. 890–895, 2016.

K. V. Silva-Vidal and N. Méndez-Sánchez, “The good, the bad and the ugly of the new treatments for hepatitis C virus,” Annals of Hepatology, vol. 13, no. 5, pp. 574-575, 2014.

D. B. Olmedo, S. A. Cader, and L. C. Porto, “IFN-λ gene polymorphisms as predictive factors in chronic hepatitis C treatment-naive patients without access to protease inhibitors,” Journal of Medical Virology, vol. 87, no. 10, pp. 1702–1715, 2015.

D. B. Rein, J. S. Wittenborn, B. D. Smith, D. K. Liffmann, and J. W. Ward, “The cost-effectiveness, health benefits, and financial costs of new antiviral treatments for hepatitis C virus,” Clinical Infectious Diseases, vol. 61, no. 2, pp. 157–168, 2015.

R. Brandão, R. Marcelino, F. Gonçalves et al., “Characterization of NS5A and NS5B resistance-associated substitutions from genotype 1 hepatitis C virus infected patients in a Portuguese cohort,” Viruses, vol. 10, no. 5, p. 223, 2018.

N. Echeverría, G. Betancour, F. Gámbaro et al., “Naturally occurring NS3 resistance-associated variants in hepatitis C virus genotype 1: their relevance for developing countries,” Virus Research, vol. 223, pp. 140–146, 2016.

A. L. McCormick, L. Moynihan, M. J. Macartney et al., “Baseline drug-resistance mutations are detectable in HCV genes NS3 and NS5A but not NS5B in acute and chronic HIV-coinfected patients,” Antiviral Therapy, vol. 20, no. 3, pp. 361–363, 2015.

Z. Plaza, V. Soriano, E. Vispo et al., “Prevalence of natural polymorphisms at the HCV NS5A gene associated with resistance to daclatasvir, an NS5A inhibitor,” Antiviral Therapy, vol. 17, no. 5, pp. 921–926, 2012.

C. Sarrazin, H. Dvory-Sobol, E. S. Svarovskaia et al., “Prevalence of resistance-associated substitutions in HCV NS5A, NS5B, or NS3 and outcomes of treatment with ledipasvir and sofosbuvir,” Gastroenterology, vol. 151, no. 3, pp. 501–512.e1, 2016.

V. Cento, M. Aragri, E. Teti et al., “Optimal cure rate by personalized HCV regimens in real-life: a proof-of-concept study,” Journal of Antimicrobial Chemotherapy, vol. 72, no. 12, pp. 3420–3424, 2017.

R. A. Fridell, C. Wang, J. H. Sun et al., “Genotypic and phenotypic analysis of variants resistant to hepatitis C virus nonstructural protein 5A replication complex inhibitor BMS-790052 in humans: in vitro and in vivo correlations,” Hepatology, vol. 54, no. 6, pp. 1924–1935, 2011.

K. M. Kitrinos, A. C. Corsa, A. Worth et al., “Nonstructural protein 5A resistance profile in patients with chronic hepatitis C treated with ledipasvir-containing regimens without sofosbuvir,” Journal of Viral Hepatitis, vol. 25, no. 2, pp. 126–133, 2018.

R. Kozuka, H. Hai, H. Motoyama et al., “The presence of multiple NS5A RASs is associated with the outcome of sofosbuvir and ledipasvir therapy in NS5A inhibitor-naïve patients with chronic HCV genotype 1b infection in a real-world cohort,” Journal of Viral Hepatitis, vol. 25, no. 5, pp. 535–542, 2018.

Z. Li, Z. W. Chen, H. Li, H. Ren, and P. Hu, “Prevalence of hepatitis C virus-resistant association substitutions to direct-acting antiviral agents in treatment-naïve hepatitis C genotype 1b-infected patients in western China,” Infection and Drug Resistance, vol. 10, pp. 377–392, 2017.

H. Fujii, A. Umemura, T. Nishikawa et al., “Real-world efficacy of daclatasvir and asunaprevir with respect to resistance-associated substitutions,” World Journal of Hepatology, vol. 9, no. 25, pp. 1064–1072, 2017.

F. Malta, K. V. Gaspareto, G. Lisboa-Neto, F. J. Carrilho, M. C. Mendes-Correa, and J. R. R. Pinho, “Prevalence of naturally occurring NS5A resistance-associated substitutions in patients infected with hepatitis C virus subtype 1a, 1b, and 3a, co-infected or not with HIV in Brazil,” BMC Infectious Diseases, vol. 17, no. 1, p. 716, 2017.

D. J. Bartels, J. C. Sullivan, E. Z. Zhang et al., “Hepatitis C virus variants with decreased sensitivity to direct-acting antivirals (DAAs) were rarely observed in DAA-naive patients prior to treatment,” Journal of Virology, vol. 87, no. 3, pp. 1544–1553, 2013.

F. Suzuki, H. Sezaki, N. Akuta 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,” Journal of Clinical Virology, vol. 54, no. 4, pp. 352–354, 2012.

S. Zeuzem, R. Ghalib, K. R. Reddy et al., “Grazoprevir–elbasvir combination therapy for treatment-naive cirrhotic and noncirrhotic patients with chronic hepatitis C virus genotype 1, 4, or 6 infection: a randomized trial,” Annals of Internal Medicine, vol. 163, no. 1, pp. 1–13, 2015.

X. Forns, S. C. Gordon, E. Zuckerman et al., “Grazoprevir and elbasvir plus ribavirin for chronic HCV genotype-1 infection after failure of combination therapy containing a direct-acting antiviral agent,” Journal of Hepatology, vol. 63, no. 3, pp. 564–572, 2015.

F. C. Lahser, K. Bystol, S. Curry et al., “The combination of grazoprevir, a hepatitis C virus (HCV) NS3/4A protease inhibitor, and elbasvir, an HCV NS5A inhibitor, demonstrates a high genetic barrier to resistance in HCV genotype 1a replicons,” Antimicrobial Agents and Chemotherapy, vol. 60, no. 5, pp. 2954–2964, 2016.

P. Krishnan, R. Tripathi, G. Schnell et al., “Resistance analysis of baseline and treatment-emergent variants in hepatitis C virus genotype 1 in the AVIATOR study with paritaprevir-ritonavir, ombitasvir, and dasabuvir,” Antimicrobial Agents and Chemotherapy, vol. 59, no. 9, pp. 5445–5454, 2015.

K. V. Gaspareto, R. M. Ribeiro, F. de M Malta et al., “Resistance-associated variants in HCV subtypes 1a and 1b detected by ion torrent sequencing platform,” Antiviral Therapy, vol. 21, no. 8, pp. 653–660, 2016.

S. Nakamoto, T. Kanda, S. Wu, H. Shirasawa, and O. Yokosuka, “Hepatitis C virus NS5A inhibitors and drug resistance mutations,” World Journal of Gastroenterology, vol. 20, no. 11, pp. 2902–2912, 2014.

K. H. Peiffer, J. Vermehren, L. Kuhnhenn et al., “Interferon-free treatment choice according to baseline RASs leads to high SVR rates in HCV genotype 1 infected patients,” Journal of Infection and Chemotherapy, vol. 24, no. 7, pp. 524–530, 2018.

W. Kati, G. Koev, M. Irvin et al., “In vitro activity and resistance profile of dasabuvir, a nonnucleoside hepatitis C virus polymerase inhibitor,” Antimicrobial Agents and Chemotherapy, vol. 59, no. 3, pp. 1505–1511, 2015.

S. Zeuzem, I. M. Jacobson, T. Baykal et al., “Retreatment of HCV with ABT-450/r–ombitasvir and dasabuvir with ribavirin,” The New England Journal of Medicine, vol. 370, no. 17, pp. 1604–1614, 2014.

M. C. B. Castilho, A. N. Martins, I. S. Horbach et al., “Association of hepatitis C virus NS5B variants with resistance to new antiviral drugs among untreated patients,” Memórias do Instituto Oswaldo Cruz, vol. 106, no. 8, pp. 968–975, 2011.

H. Dutartre, C. Bussetta, J. Boretto, and B. Canard, “General catalytic deficiency of hepatitis C virus RNA polymerase with an S282T mutation and mutually exclusive resistance towards 2-modified nucleotide analogues,” Antimicrobial Agents and Chemotherapy, vol. 50, no. 12, pp. 4161–4169, 2006.

S. Fourati and J. M. Pawlotsky, “Virologic tools for HCV drug resistance testing,” Viruses, vol. 7, no. 12, pp. 6346–6359, 2015.

N. Echeverría, D. Chiodi, P. López et al., “IL28B gene polymorphism rs12979860, but not rs8099917, contributes to the occurrence of chronic HCV infection in Uruguayan patients,” Virology Journal, vol. 15, no. 1, p. 40, 2018.