Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Relapse or reinfection after failing hepatitis C direct acting antiviral treatment: Unravelled by phylogenetic analysis

  • Lize Cuypers ,

    Contributed equally to this work with: Lize Cuypers, Ana Belén Pérez

    Roles Conceptualization, Formal analysis, Methodology, Visualization, Writing – original draft

    lize.cuypers@kuleuven.be

    Affiliations KU Leuven–University of Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Clinical and Epidemiological Virology, Leuven, Belgium, Peter Medawar Building for Pathogen Research, Nuffield Department of Medicine, University of Oxford, Oxford, United Kingdom

  • Ana Belén Pérez ,

    Contributed equally to this work with: Lize Cuypers, Ana Belén Pérez

    Roles Data curation, Methodology, Writing – original draft

    Affiliation Clinical Microbiology Department, University Hospital San Cecilio Granada, Instituto de Investigación Ibs. Granada, Spain

  • Natalia Chueca,

    Roles Data curation, Methodology

    Affiliation Clinical Microbiology Department, University Hospital San Cecilio Granada, Instituto de Investigación Ibs. Granada, Spain

  • Teresa Aldamiz-Echevarría,

    Roles Data curation

    Affiliation Infectious Diseases Unit, Hospital Gregorio Marañón, Madrid, Spain

  • Juan Carlos Alados,

    Roles Data curation

    Affiliation Clinical Microbiology, University Hospital Jerez, Cadiz, Spain

  • Ana María Martínez-Sapiña,

    Roles Data curation

    Affiliation Clinical Microbiology, Hospital Miguel Servet, Zaragoza, Spain

  • Dolores Merino,

    Roles Data curation

    Affiliation Clinical Microbiology, Hospital Infanta Elena, Huelva, Spain

  • Juan Antonio Pineda,

    Roles Data curation

    Affiliation Infectious Diseases Unit, University Hospital de Valme, Sevilla, Spain

  • Francisco Téllez,

    Roles Data curation

    Affiliation UGC Enfermedades Infecciosas y Microbiología, Hospital La Línea, AGS Campo de Gibraltar, Cadiz, Spain

  • Nuria Espinosa,

    Roles Data curation

    Affiliation Clinical Microbiology, Hospital Virgen del Rocío, Sevilla, Spain

  • Javier Salméron,

    Roles Data curation

    Affiliation Hepatology Unit, University Hospital San Cecilio Granada, Instituto de Investigación Ibs. CIBERehd, Granada, Spain

  • Antonio Rivero-Juarez,

    Roles Data curation

    Affiliation Infectious Diseases Unit. Instituto Maimonides de Investigación Biomédica de Córdoba (IMIBIC). Hospital Universitario Reina Sofía de Córdoba. Universidad de Córdoba, Córdoba, Spain

  • María Jesús Vivancos,

    Roles Data curation

    Affiliation Clinical Microbiology, Hospital Ramón y Cajal, Madrid, Spain

  • Víctor Hontañón,

    Roles Data curation

    Affiliation Clinical Microbiology, University Hospital La Paz, Madrid, Spain

  • Anne-Mieke Vandamme,

    Roles Methodology, Supervision, Writing – original draft

    Affiliations KU Leuven–University of Leuven, Department of Microbiology and Immunology, Rega Institute for Medical Research, Clinical and Epidemiological Virology, Leuven, Belgium, Center for Global Health and Tropical Medicine, Microbiology Unit, Institute for Hygiene and Tropical Medicine, University Nova de Lisboa, Lisbon, Portugal

  •  [ ... ],
  • Féderico García

    Roles Conceptualization, Methodology, Supervision, Writing – original draft

    Affiliation Clinical Microbiology Department, University Hospital San Cecilio Granada, Instituto de Investigación Ibs. Granada, Spain

  • [ view all ]
  • [ view less ]

Abstract

Despite high response rates associated to hepatitis C virus (HCV) treatment, no protective immunity is acquired, allowing for reinfection and continued infectiousness. Distinguishing between relapse and reinfection is crucial for patient counselling and to choose the most appropriate retreatment. Here, refined phylogenetic analysis using multiple genes served to assess genotype and reinfection for 53 patients for whom the virus was sampled before start of therapy and at time of sustained virological response evaluation at week 12. At baseline, genotypes were determined as HCV1a (41.5%), HCV1b (24.5%), HCV4 (18.9%) and HCV3a (15.1%), while six cases revealed to be discordantly assigned by phylogeny and commercial assays. Overall, 60.4% was co-infected with HIV. The large majority was classified as people who inject drugs (78.6%), often co-infected with HIV. Transmission was sexual in seven cases, of which five in HIV-positive men-who-have-sex-with-men. Overall, relapse was defined for 44 patients, while no conclusion was drawn for four patients. Five patients were reinfected with a different HCV strain, of which three with a different genotype, showing that phylogeny is needed not only to determine the genotype, but also to distinguish between relapse and intra-subtype reinfection. Of note, phylogenies are more reliable when longer fragments of the viral genome are being sequenced.

Introduction

Despite occasional claims of integration events [1], the hepatitis C virus (HCV) is considered to have no latent forms and to rely on ongoing replication. Therefore, in contrast to HIV, patients infected with HCV can either spontaneously clear their infection or can be cured thanks to treatment with highly effective antiviral regimens which are nowadays based on molecules that directly target viral proteins, the so called direct-acting antivirals (DAAs). Viral cure rates have dramatically increased to 90–95% following the introduction of interferon-free (IFN-free) DAA combination therapies [2]. Although less crucial than during the IFN era, the HCV genotype (GT) still plays a substantial role in deciding the best treatment regimen of a patient [3], hand in hand with the respective stage of liver fibrosis or cirrhosis and the previous experience of the patient to IFN-based treatment regimens.

In spite of high treatment response rates, no protective immunity is built [4], so patients that are cured at the end of therapy can be infected with a new HCV strain, and might be still at risk to develop liver diseases [5], and to form a source of further transmission events. Patients reinfected with a strain determined to be of a different genotype or subtype than the previous strain they were infected with, can easily be identified using genotyping assays. However, when a reinfection with a similar strain of the same subtype occurs, phylogenetic analysis is required to distinguish a reinfection from a virological relapse [67]. While iatrogenic infected HCV patients are expected to clear their infection and to remain cured, viral eradication is not expected soon at global level, given the high number of individuals that is unaware of their infection status, the low treatment uptake and the continued de novo HCV incidence that in many countries exceeds the rate of curative treatment. Especially in specific populations such as people who inject drugs (PWID) and HIV-positive men who have sex with men (MSM), viral eradication is hampered by the assumed high reinfection and transmission rates after successful DAA treatment. However, wide ranges in the number of reinfection cases have been reported for these risk groups [89], depending on the patient population studied, the treatment era in which the study was planned and if reinfection followed viral clearance due to the immune system or due to a prescribed treatment regimen. In general, reinfection rates are reported to be lower for PWIDs, although depending on their active drug use [1011], compared to HIV-positive MSM [10,1213]. For PWIDs, generally rather low reinfection rates from 0 to 5 cases per 100 person-years were reported [1314], although rates can go up to 10% [9]. In case of patients that are co-infected with HIV, of which the majority are MSM, rates easily increase to 15–25% [10,15]. Moreover, recently also a high prevalence of HCV infection has been reported in HIV-negative MSM that are enrolled for HIV pre-exposure prophylaxis [16]. Distinguishing between virologic relapse and infection with a new viral strain is highly important to determine the true treatment efficacy of the current DAA regimens, to define the most appropriate retreatment for these patients [17], and more importantly for appropriate patient counselling. Treatment scale-up, diagnostic testing approaches and strategies to prevent new infections as well as reinfections, need to occur concomitantly to accomplish worldwide viral eradication of HCV.

Methods

Inclusion of patients from the GEHEP-004 cohort

HCVREsp-GEHEP004 is a prospective multicentre cohort including 7189 HCV infected patients treated with IFN-free DAA regimens, attending 54 different Spanish centres, all part of the Group for the Study of Viral Hepatitis (GEHEP). So far, this cohort enrolled around 450 patients that were unsuccessfully treated with DAAs (update from [18]). From this last group, in 53 patients, samples at two time points were referred for resistance testing, and all were included in this analysis. For all other patients, only one of these two time points was sampled. The patients included in this study were infected with various HCV GTs and were treated between 2014 and 2016 with the standard of care DAA regimens at that time, in Spain. In total, 37.7% of all patients (20/53) failed previous therapy based on sofosbuvir + ledipasvir, with or without ribavirin, followed by a treatment consisting of three DAA classes (3D: paritaprevir + ombitasvir + dasabuvir + ribavirin) for 20.8% of the patients (11/53). Next in line was the regimen based on sofosbuvir + daclatasvir + ribavirin. A lower number of patients had initiated therapy with sofosbuvir + simeprevir + ribavirin, simeprevir + daclatasvir + ribavirin, sofosbuvir + ribavirin or the 2D combination (without dasabuvir) + ribavirin (S1 Table). The relatively high therapy failure rate can be attributed to the use of inferior regimens during the early days of DAA therapy in 2014, and due to the incorrect determination of HCV genotypes using commercial assays, resulting in the choice of a potential suboptimal regimen.

Genetic sequencing

HCV was sampled at two time points for each patient, before start of treatment (at baseline) and at time of sustained virological response evaluation at week 12 (further on referred to as at SVR12 evaluation– 12 weeks after end of treatment). All patients had an undetectable viral load at the end of treatment, hence why the detection of virus 12 weeks after end of treatment either indicates the occurrence of reinfection or relapse, the latter considered as treatment failure [19]. In the context of drug resistance testing, genetic sequencing was performed for three different genetic regions in the HCV genome, more particularly NS3, NS5A and NS5B, routinely being sequenced using in house developed assays, respectively amplifying 513, 258 and 312 nucleotides. Depending on the composition of the administered DAA regimen, all three, or only one or two genes, were sequenced. For the purpose of this study, retrospective Sanger sequencing was performed to recover the genetic sequence of at least two regions for as many patients as possible. For the cases here distinguished as a reinfection, next-generation sequencing (NGS) was performed on the baseline and follow-up sample of the respective patient, starting from the same NS5B PCR product as obtained for Sanger sequencing, to rule out the occurrence of a mixed infection at the start of treatment and/or at time of SVR12 evaluation. Until the end of 2016, the 454 GS Junior pyrosequencing methodology (Branford, Connecticut, USA) was used to perform NGS, however afterwards the Illumina platform (San Diego, CA, USA) was applied, based on a modified protocol described in [20]. The HCV genotype of all individual reads was determined to rule out the occurrence of a mixed infection. All sequences generated within this study, have been submitted to Genbank (accession numbers MG983221-MG983474). The Ethics Committee of the San Cecilio Hospital, Granada, approved the study, sequencing experiments were performed in accordance to good laboratory practices, and no informed consent was required as patient information was anonymized prior to analysis.

Determination of the HCV genotype and between-subtype or genotype recombination

Genotypes were determined using the commercial assay Versant HCV Genotype 2.0 assay (LiPA) for the majority of the baseline samples (86.8% or 46/53), while the genotype of 9.4% (5/53) and 3.8% (2/53) of the baseline samples was defined by the Abbott Real Time HCV Genotype II assay and Trugene HCV Genotyping Kit, respectively. The HCV genotype and subtype of the samples was also determined by manual phylogenetic analysis and the use of subtyping tools COMET and Oxford [2122]. In case of discordance between the different tools, a bootscan analysis with Simplot was performed to assess potential breakpoints for recombination [23]. Resistance-associated substitutions (RASs) were evaluated for the NS3, NS5A and NS5B sequences at time of SVR12 evaluation, based on variants reported in literature [24].

Dataset construction and assessment of within-subtype recombination

The ten most similar sequences to each of the 53 Spanish taxa, were retrieved by the use of the standalone BLAST tool [25]. The resulting datasets were constituted based on the fragment that was sequenced, more specifically resulting in seven datasets, either consisting of sequences coding for one fragment (NS3, NS5A, NS5B) or concatenated in case more than one region was sequenced (NS3-NS5A-NS5B, NS3-NS5A, NS3-NS5B and NS5A-NS5B). The latter strategy increased the genomic fragment length, and therefore the phylogenetic signal. Sequences were aligned per dataset and manually edited to assure high quality, using Seaview and MEGA 7.0 [2627]. An extensive recombination analysis was done using RDP4 and TreePuzzle [2829]. The likelihood-mapping algorithm in TreePuzzle was used to assess the percentage of conflicting phylogenetic signal, represented by the dots at the sides of the triangle. Potential recombination events indicated by RDP4 were only confirmed if two parental sequences could be identified, and if breakpoints were significant. Since this is depending on the strains included in the analysis, absence of evidence of within-subtype recombination might not be equal to evidence of absence of within-subtype recombination.

In-depth phylogenetic analysis

The resulting codon-correct alignments were inferred by neighbour-joining (NJ) and maximum-likelihood (ML) approaches, applying the GTR gamma model to allow for among-site variation, and evaluating tree robustness with 1000 bootstrap replicates. For the NJ approach, trees were inferred as implemented in MEGA 7.0, while for the (approximate) ML algorithm, trees were constructed using FastTree and RAxML [3031]. Additionally, positions known to be impacted by RASs [24], were removed from the alignment, to evaluate them as potential source of confounding in the inference of tree topologies.

When comparing strains before and after treatment in the same patient, evidence of reinfection was defined as a difference in HCV genotype or subtype, or as a significantly different clustering located in different clades in the phylogenetic tree. Evidence of virological relapse was defined as significant clustering in the same clade, while no conclusion was drawn when clades were supported with a bootstrap <70%. Despite this predefined threshold, the lowest bootstrap value to define a relapse was found to be 83%.

Results

Predominance of HCV1a infection, PWIDs and HIV co-infection

At baseline, 22 out of 53 patients (41.5% or 22/53) were phylogenetically determined to be infected with HCV1a, followed by HCV1b (24.4% or 13/53), HCV4 (18.9% or 10/53: 7.5% 4a – 11.3% 4d), and HCV3a (15.1% or 8/53). A large share of patients was co-infected with HIV (60.4% or 32/53). The transmission route of infection was known for 79% of patients (42/53), with all patients lacking information for risk of transmission being mono-infected with HCV. The large majority of patients with known route of transmission were reported to be (former) PWID (78.6% or 33/42), of whom 81.8% (27/33) was found to be co-infected with HIV and more than half of them with HCV1a (51.5% or 17/33). HCV transmission through sexual contact was inferred for seven patients, of which five were reported as MSM and two as heterosexual, all of them concomitantly infected with HIV. One patient was infected through blood transfusion during childhood and another one during haemodialysis. For none of the patients, a dual route of transmission was reported in the clinical database.

HCV genotype misclassifications

HCV genotype and subtype assignment by phylogenetic analysis and the use of well-known subtyping tools was in agreement with the assignment by commercial assays for 66.0% of all baseline samples (35/53). The majority of misclassifications by commercial assays was due to assignments solely on genotype level, lacking information on subtype level. However, also six baseline samples were wrongly classified on HCV genotype and subtype level (S1 Table). The sequence at baseline for one patient of the cohort was concordantly identified to be HCV1a by all genotype determination methods, while the strain at time of SVR12 evaluation was assigned as HCV1a or HCV1b, dependent on the genetic region analysed and the determination method. In-depth analysis of the concatenated alignment covering regions NS5A and NS5B, could offer a definite classification of this patient’s sequence at time of SVR12 evaluation as HCV1a, excluding the occurrence of an inter-subtype recombination event. Moreover, in none of the 53 patients in the study, evidence of a recombination event was identified using Simplot or RDP4. The ability to detect recombination was evaluated using TreePuzzle, showing only a small percentage (0.3–0.5%) of phylogenies characterized by conflicting signal. Evaluating RASs, either naturally occurring or substituted under drug selective pressure, showed that patients in this study mainly harboured variants in the NS3 and NS5A genes, while only for three patients RASs were detected in NS5B (282T/R). Two of the latter patients also harboured NS5A RASs, while in total 10 patients harboured dual variants in NS3 and NS5A. Single NS3 and NS5A variants were detected in 8 and 14 patients, respectively. Patients defined to have experienced a reinfection were analysed separately, showing that four out of five patients acquired NS5A RASs over the course of treatment, while for one of them already a natural variant was present at baseline. Details on the specific variants can be found in S2 Table.

Phylogenetic analysis of concatenated alignments

Because of low phylogenetic signal when only one genetic region was used, additional genes were sequenced retrospectively. As a result, more than one genetic region was available for both time points sampled for all patients, except for two patients for which amplification was unsuccessful for regions NS3 and NS5B, respectively (S1 Table: patients 30 and 49). Therefore, while the tree suggested the occurrence of a relapse, this evidence was not conclusive since bootstrap values were <70%. However, for one of the two patients, reinfection is not highly likely since this patient declared to be infected by blood transfusion during childhood and to show no persistent risk behaviour to acquire a new HCV infection.

For all other patients, concatenated alignments were constructed. Phylogenetic inference using the four concatenated alignments (NS3 –NS5A –NS5B, NS3 –NS5A, NS3 –NS5B and NS5A –NS5B) showed that 44 out of 51 patients had experienced a virological relapse, while five (9.4%) were classified as reinfected with a different HCV strain (S1 Table). For three of these five patients (S1 Table: patients 5, 14 and 36), the sequences at baseline and at time of SVR12 evaluation differed in HCV genotype and/or subtype. NGS experiments ruled out the occurrence of a mixed infection at baseline and at time of SVR12 evaluation for all three patients. Two patients (S1 Table: patients 1 and 18) were reinfected with the same HCV subtype, phylogenetically clustering in a different clade for the two sampled time points (Fig 1). Of the five patients that were defined to be reinfected, three were classified as PWID, one identified himself as MSM and for one the potential route of transmission was unknown. Three out of these five patients were co-infected with HIV. For two additional patients, no conclusions could be drawn, either due to close clustering of both strains, however supported by bootstrap replicates lower than 70%, or due to an inconclusive clustering between the clades HCV1a and HCV1b (S1 Table: patients 6 and 7, Figs 1 and 2). Conclusions drawn from the phylogenetic analyses were not influenced by the removal of drug resistance-associated positions in the concatenated alignments.

thumbnail
Fig 1. Maximum-likelihood tree of the concatenated alignment covering genetic regions NS5A and NS5B.

Both the entire phylogenetic tree (left) and the HCV1b clade in detail (right) are visualized. Patient 18, coloured in green, is a clear example of reinfection with a different HCV strain, since although both viruses at baseline and at time of SVR12 evaluation are classified as HCV1b, they cluster in a different clade in the tree, with a bootstrap support of 100% for the segregation of the different clades. For patient 4, coloured in blue, both strains cluster together with a high bootstrap support (100%), suggesting that this patient experienced a virological relapse. Bootstrap replicates are only visualized for patients who experienced a relapse, all indicated in blue. Patients for which only a baseline or SVR12 evaluation sequence is included in the NS5A-NS5B alignment, are coloured in red. The sequence at time of SVR12 evaluation and baseline for patients 7 and 49 respectively, cluster outside the large HCV1b clade, which might be due to a potential event of recombination. However, detailed analyses with TreePuzzle and RDP4 for all 53 patients, could not support a recombination event. The in-depth Simplot analysis for patients 7 and 49 showed that both strains were classified as HCV1a. Despite the absence of evidence of recombination, we did not draw any conclusions concerning the occurrence of reinfection or relapse for patient 7. The bar at the bottom represents the number of nucleotide substitutions per site.

https://doi.org/10.1371/journal.pone.0201268.g001

thumbnail
Fig 2. Maximum-likelihood tree of the concatenated alignment covering genetic regions NS3 and NS5A.

Both the entire phylogenetic tree and the HCV4 (4a and 4d) clade in detail are visualized. Clusters of patients coloured in blue, are suggested to have experienced a virological relapse, supported by bootstrap values >70%. However, for patient 6, indicated in red, only a bootstrap value of 38 was obtained, resulting in lack of evidence to distinguish between a relapse and reinfection. For patient 39, the fragment NS3-NS5A was only sequenced for the sample at time of SVR12 evaluation. The bar at the bottom represents the number of nucleotide substitutions per site.

https://doi.org/10.1371/journal.pone.0201268.g002

Discussion

Reinfection rates in real-life settings may be underestimated, especially when the patient is reinfected with the same subtype. Here we report how detailed phylogenetic analysis is needed to discriminate intra-subtype reinfection from relapse. This study demonstrates the importance of genetic sequencing, not only to define the most appropriate (re)treatment, but also to perform phylogeny needed to determine the correct HCV genotype and subtype that a patient is infected with and to distinguish between relapse and reinfection.

Concatenated genes were used to distinct between relapse and reinfection, however for two patients no conclusion could be drawn due to a bootstrap support <70% of the respective clade, or the inconsistent clustering between different HCV subtypes, respectively. Additionally, for two patients not more than one genetic region was successfully sequenced. For the other patients, the majority had a virological relapse (44/49). These patients are eligible for retreatment. Until recently, it was advised to choose the DAA regimen based on the presence of natural occurring or emergent RASs [23,32]. Since the approval of the new combinations consisting of sofosbuvir, velpatasvir and voxilaprevir, or glecaprevir and pibrentasvir, AASLD and EASL guidelines no longer routinely recommend RAS testing at time of failure. Moreover, in the era of true pan-genotypic combinations, the interest of determining the HCV genotype at baseline is diminishing as the prescribed regimen no longer depends on it. We hypothesize that in developed countries where first-line regimens will consist of the newest therapy combinations, baseline genotyping will be soon seen as obsolete. However, in the light of these changes, phylogenetic analysis used to determine viral failures as relapse or reinfection, will gain largely in importance, as even reinfections with a different genotype as the strain of the first infection, would be missed.

Almost 10% of patients (5/53) were found to be reinfected with a different HCV strain, two of which (2/5 or 40% of reinfections, while 2/53 or 3.8% of all patients) would have been missed in absence of phylogenetic analysis since the sequences at baseline and time of SVR12 evaluation were of the same HCV subtype. Phylogenetic analysis confirmed reinfection in three patients with discordant genotype or subtype. Three out of these five reinfected patients reported to be infected by intravenous drug use, while one was an MSM, and for the fifth patient no route of transmission was known. Additionally, three of them were known to be co-infected with HIV. It has been reported that equipment used during, before, or after sexual activity [33], high-risk sexual behaviour in general [3435], as well as the use of drugs [3637], could be associated to a higher risk of HCV transmission, showing a potential role of multifactorial risk behaviour instead of transmission through one dominant route.

The main limitation of this study was that it was not specifically designed to evaluate reinfection rates; as the aim of the GEHEP cohort is to evaluate the emergence of RASs at failure. Moreover, the selection of 53 patients out of a cohort of 450 may have introduced a bias with respect to a population-based analysis, however these were the only patients for which samples at two time points were available. Since this was a random set of 53 patients, the general characteristics did not differ from the other 397 patients. Complete information regarding the risk of transmission is unavailable for the remaining 397 patients, not allowing us to evaluate the extent of this potential bias. Additionally, a large share of the study population (>60%) was observed to be co-infected with HIV, which is much higher than reported for the overall HCV infected population in Spain [38]. Nevertheless, we believe these limitations did not heavily impact our main conclusion, supporting the added value of phylogeny to rule out intra-subtype reinfection.

Various genes of the HCV genome have been subjected to phylogenetic analysis, with the envelope genes E1 and E2 dominating the study landscape [3941], more particularly hypervariable region-1 as this region is characterized by a high evolutionary rate [42]. However, in the context of treatment follow-up, these genes are not routinely sequenced, in contrast to proteins NS3, NS5A and NS5B. Despite the benefit of a higher phylogenetic signal, the added value of using viral genetic sequences obtained in clinical practice outranges the use of a more phylogenetically applicable gene such as E2, at least in the context of our study. Moreover, recent studies tend to use more often the NS3, NS5A and/or NS5B genes for phylogenetic purposes, although associated to a lower degree of genetic variability [17,43].

Phylogeny of the separate genetic regions often resulted in inconsistencies between the different examined genes, due to unresolved phylogenies. Especially when using the single fragment NS5A, phylogenies consisted of branches supported with low bootstrap values, since these respective sequences only covered 76 amino acid positions. It is therefore worrying that the rise in large-scale genetic sequencing in the context of clinical follow-up and drug resistance testing, has been and still is targeting only a limited number of amino acid sites in the HCV genome. Such data cannot be used to investigate reinfection, instead, longer regions such as full-length genomes are needed to obtain more robust phylogenies, as most often this approach highly improves the phylogenetic signal [17,4041]. When concatenating multiple genes, combining regions characterized by a different evolutionary rate has been proposed to increase the capacity to identify transmission chains [44], while others underline that there is no ultimate gene to perform evolutionary analyses of genetically diverse viruses [45]. Fortunately, since HCV rarely recombines [46], at least studied on inter-genotype or inter-subtype level, thorough analyses using multiple genes should still have a genome-wide representativeness. In our dataset, we did not observe evidence of recombination, neither inter- nor intra-subtype or–genotype.

Less than 70% of HCV genotype assignments by commercial assays appeared to be consistent with well-known subtyping tools and phylogenetic analysis. Although the majority of misclassifications by commercial assays was due to an incomplete assignment on subtype level, still six baseline samples (11%) were wrongly classified on HCV genotype and/or subtype level, which may have potentially resulted in the choice of a less effective DAA therapy for these patients, as they were treated between 2014 and 2016. Commercial assays are known to be unreliable in case of mixed infections or recombination events since they only target a limited region of the HCV genome. Inconsistencies between commercial assays and genetic sequencing followed by phylogeny as HCV genotyping tool have also been reported in absence of recombination, ranging around 10% of misclassifications [4749]. This study confirms the added value of genetic sequencing next to (or potentially instead of) the use of commercial assays in the assessment of HCV genotypes, as well as to acquire additional information on the presence of RASs which have proven to reduce susceptibility to certain antiviral drugs [50]. Notably, NGS was performed on the samples of both time points from patients identified as being reinfected, to rule out the presence of a mixed infection, as studies in the pre-DAA era have shown that the emergence of new viral strains following therapy failure is often associated to the emerging dominance of pre-existing minority variants rather than an actual reinfection [51].

The observation that a large part of the infections was assigned to be GT1 (66%) is overall in agreement with a large Spanish study on HCV genotype prevalence and distribution [38]. However, in our study a higher proportion of cases was infected with HCV1a and HCV4 compared to the cohort of [38], potentially due to the difference in population characteristics, since our study has a larger share of HIV/HCV co-infected patients (60.4% versus 19.1%). For these patients, both intravenous drug use and sexual transmission have been defined as the major routes of transmission, for which a higher proportion of HCV1a and HCV4 infections was reported compared to the overall population [38]. Especially HCV1a is more commonly observed among HCV infected PWIDs [52]. The majority of patients in this study were reported to be infected through this transmission route, in agreement with a rising HCV epidemic dominated by PWIDs in Spain, as reported previously [53]. All patients for which the route of transmission was unknown, were identified to be only infected with HCV and not with HIV, stressing the added value of capturing information concerning the potential route of transmission, as done in HIV clinical care.

In this cohort of unsuccessfully treated patients, the majority experienced a true DAA failure. Since about half of those reinfected, showed the same subtype as at baseline, phylogenetic analysis is needed, not only to determine the correct HCV genotype, but also to distinguish between relapse and reinfection. Given that 11% of genotypes were misclassified using commercial assays, and that 4% had a same-subtype reinfection only detectable by phylogeny, the potential of misclassifying a reinfection as a failure could be as high as 15%.

Supporting information

S1 Graphical Abstract. Graphical overview of the study cohort, methodology and results.

https://doi.org/10.1371/journal.pone.0201268.s001

(SVG)

S1 Table. Overview of the assignment of virologic relapse or reinfection for the Spanish cohort.

For each of the 53 patients in the cohort, the genetic region(s) that was or were sequenced for both time points is (are) listed, as well as the HCV genotype and subtype determined for the patient. In case of a misclassification of the HCV genotype by a commercial assay (excluding the ones due to commercial assays not classifying down to subtype level), cells are marked in grey. Phylogenetic analysis showed either evidence for virologic relapse or reinfection, although for four patients no conclusion could be drawn due to lack of phylogenetic signal or bootstrap values <70% or inconsistent clustering, also in this case the cell is colored grey.

https://doi.org/10.1371/journal.pone.0201268.s002

(DOCX)

S2 Table. Overview of the frequency and distribution of resistance-associated substitutions.

For all 53 patients, resistance-associated substitutions (RASs) were listed at baseline and at time of SVR12 evaluation, in parallel to HCV genotype information, HIV co-infection status and classification of reinfection or relapse. For the latter, the symbol ‘?’ indicates that phylogenetic analysis could not resolve whether a relapse or reinfection occurred. In case no genetic sequencing was performed for a particular gene, this was indicated by the symbol ‘-‘ in the cell, while ‘None’ states that no RASs were detected.

https://doi.org/10.1371/journal.pone.0201268.s003

(DOCX)

Acknowledgments

The authors wish to thank Lore Vinken for her assistance and advice in the performance of the Simplot analysis. The authors are also grateful for all the participants of the HepCREsp-GEHEP004 cohort.

References

  1. 1. Moustafa A, Xie C, Kirkness E, Biggs W, Wong E, Turpaz Y, et al. The blood DNA virome in 8,000 humans. PLoS Pathog. 2017;13.
  2. 2. European Association for the Study of the Liver. EASL recommendations on treatment of hepatitis C 2018. J Hepatol. 2018.
  3. 3. Cuypers L, Ceccherini-Silberstein F, Van Laethem K, Li G, Vandamme A-M, Rockstroh JK. Impact of HCV genotype on treatment regimens and drug resistance: a snapshot in time. Rev Med Virol. 2016;26: 408–34. pmid:27401933
  4. 4. Abdel-Hakeem MS, Shoukry NH. Protective immunity against hepatitis C: many shades of gray. Front Immunol. 2014;5: 274. pmid:24982656
  5. 5. van der Meer AJ, Berenguer M. Reversion of disease manifestations after HCV eradication. J Hepatol. 2016;65: S95–108. pmid:27641991
  6. 6. Arrais TC, Van Dooren S, Vandamme AM, Brechot C, Rimlinger F, Silva AE, et al. Change in hepatitis C virus genotype in hemodialysis patients after end-of-treatment response to interferon monotherapy -relapse or re-infection? J Med Virol. 2008; 80: 80–6. pmid:18041025
  7. 7. van de Laar TJ, Molenkamp R, van den Berg C, Schinkel J, Beld MG, Prins M, et al. Frequent HCV reinfection and superinfection in a cohort of injecting drug users in Amsterdam. J Hepatol. 2009;51: 667–74. pmid:19646773
  8. 8. Martin TC, Martin NK, Hickman M, Vickerman P, Page EE, Everett R, et al. Hepatitis C virus reinfection incidence and treatment outcome among HIV-positive MSM. AIDS. 2013;27: 2551–7. pmid:23736152
  9. 9. Simmons B, Saleem J, Hill A, Riley RD, Cooke GS. Risk of late relapse or reinfection with hepatitis C virus after achieving a sustained virological response: a systematic review and meta-analysis. Clin Infect Dis. 2016;62: 683–94. pmid:26787172
  10. 10. Cunningham EB, Applegate TL, Lloyd AR, Dore GJ, Grebely J. Mixed HCV infection and reinfection in people who inject drugs—impact on therapy. Nat Rev Gastroenterol Hepatol. 2015;12: 218–30. pmid:25782091
  11. 11. Marco A, Esteban JI, Solé C, da Silva A, Ortiz J, Roget M, et al. Hepatitis C virus reinfection among prisoners with sustained virological response after treatment for chronic hepatitis C. J Hepatol. 2013;59: 45–51. pmid:23523577
  12. 12. Islam N, Krajden M, Shoveller J, Gustafson P, Gilbert M, Buxton JA, et al. Incidence, risk factors, and prevention of hepatitis C reinfection: a population-based cohort study. Lancet Gastroenterol Hepatol. 2017;2: 200–210. pmid:28404135
  13. 13. Pradat P, Pugliese P, Poizot-Martin I, Valantin MA, Cuzin L, Reynes J, et al. Direct-acting antiviral treatment against hepatitis C virus infection in HIV-infected patients–“En route for eradication”? J Infect. 2017;75: 234–241. pmid:28579302
  14. 14. Grebely J, Hajarizadeh B, Dore GJ. Direct-acting antiviral agents for HCV infection affecting people who inject drugs. Nat Rev Gastroenterol Hepatol. 2017;14: 641–51. pmid:28831184
  15. 15. Lambers FA, Prins M, Thomas X, Molenkamp R, Kwa D, Brinkman K, et al. Alarming incidence of hepatitis C virus re-infection after treatment of sexually acquired acute hepatitis C virus infection in HIV-infected MSM. AIDS. 2011;25: F21–7. pmid:21857492
  16. 16. Hoornenborg E, Achterbergh RCA, Schim van der Loeff MF, Davidovich U, Hogewoning A, de Vries HJC, et al. MSM starting preexposure prohylaxis are at risk of hepatitis C virus infection. AIDS. 2017;31: 1603–10. pmid:28657964
  17. 17. Sarrazin C, Isakov V, Svarovskaia ES, Hedskog C, Martin R, Chodavarapu K, et al. Late relapse versus hepatitis C virus reinfection in patients with sustained virologic response after sofosbuvir-based therapies. Clin Infect Dis. 2017;64: 44–52. pmid:27737953
  18. 18. Pérez AB, Chueca N, Fernández-Caballero JA, Toro MG-D, Martínez-Sapiña A, Navarro D, et al. Real-life retreatment of patients failing DAAs in Spain: data from the HCVREsp GEHEP-004 cohort. J Hepatol. 2017;66: S313.
  19. 19. Poordad FF, Flamm SL. Virological relapse in chronic hepatitis C. Antivir Ther. 2009;14: 303–13. pmid:19474464
  20. 20. Quer J, Gregori J, Rodríguez-Frias F, Buti M, Madejon A, Perez-del-Pulgar S, et al. High-resolution hepatitis C virus subtyping using NS5B deep sequencing and phylogeny, an alternative to current methods. J Clin Microbiol. 2015;53: 219–26. pmid:25378574
  21. 21. Struck D, Lawyer G, Ternes AM, Schmit JC, Bercoff DP. COMET: adaptive context-based modeling for ultrafast HIV-1 subtype identification. Nucleic Acids Res. 2014;42.
  22. 22. de Oliveira T, Deforche K, Cassol S, Salminen M, Paraskevis D, Seebregts C, et al. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics. 2005;21: 3797–800. pmid:16076886
  23. 23. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73: 152–60. pmid:9847317
  24. 24. Sorbo MC, Cento V, Di Maio VC, Howe AYM, Garcia F, Perno CF, et al. Hepatitis C virus drug resistance associated substitutions and their clinical relevance: Update 2018. Drug Resist Updat. 2018;37: 17–39. pmid:29525636
  25. 25. Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10: 421. pmid:20003500
  26. 26. Gouy M, Guindon S, Gascuel O. Seaview version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010;27: 221–4. pmid:19854763
  27. 27. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for bigger datasets. Mol Biol Evol. 2016;33: 1870–4. pmid:27004904
  28. 28. Martin DP, Murrell B, Khoosal A, Muhire B. Detecting and analyzing genetic recombination using RDP4. Methods Mol Biol. 2017;1525: 433–460. pmid:27896731
  29. 29. Schmidt HA, Strimmer K, Vingron M, von Haeseler A. TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics. 2002;18: 502–504. pmid:11934758
  30. 30. Price MN, Dehal PS, Arkin AP. FastTree 2 –approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5.
  31. 31. Stamatakis A. Using RAxML to infer phylogenies. Curr Protoc Bioinformatics. 2015;51: 6–14. pmid:26334924
  32. 32. AASLD 2017. Recommendations for testing, managing, and treating hepatitis C. Update September 2017. http://www.hcvguidelines.org/full-report-view
  33. 33. Apers L, Vanden Berghe W, De Wit S, Kabeya K, Callens S, Buyze J, et al. Risk factors for HCV acquisition among HIV-positive MSM in Belgium. J Acquir Immune Defic Syndr. 2015;68: 585–93. pmid:25763786
  34. 34. Hagan H, Jordan AE, Neurer J, Cleland CM. Incidence of sexually transmitted hepatitis C virus infection in HIV-positive men who have sex with men. AIDS. 2015;29: 2335–45. pmid:26258525
  35. 35. Urbanus AT, Van De Laar TJ, Geskus R, Vanhommerig JW, Van Rooijen MS, Schinkel J, et al. Trends in hepatitis C virus infections among MSM attending a sexually transmitted infection clinic; 1995–2010. AIDS. 2014;28: 781–90. pmid:24832014
  36. 36. Bollepalli S, Mathieson K, Jasiurkowski B, Hillier A, Post J, Bhanu S, et al. A comparison of risk factors for HCV-mono-infection, HIV-mono-infection, and HCV/HIV-co-infection in a community setting. Dig Dis Sci. 2008;53: 517–21. pmid:17562178
  37. 37. Lea T, Lee E, Mao L, de Wit J, Holt M. HIV and hepatitis C virus co-infection among men who have sex with men in Sydney, and associations with sexual and drug use practices. Sex Health. 2013;10: 448–51. pmid:23849088
  38. 38. Aguilera A, Navarro D, Rodríguez-Frias F, Viciana I, Martínez-Sapiña AM, Rodríguez MJ, et al. Prevalence and distribution of hepatitis C virus genotypes in Spain during the 2000–2015 period (the GEHEP 005 study). J Viral Hepat. 2017;24: 725–32. pmid:28248445
  39. 39. Al-Qahtani AA, Baele G, Khalaf N, Suchard MA, Al-Anazi MR, Abdo AA, et al. The epidemic dynamics of hepatitis c virus subtypes 4a and 4d in Saudi Arabia. Sci Rep. 2017;7: 44947. pmid:28322313
  40. 40. Jacka B, Applegate T, Krajden M, Olmstead A, Harrigan PR, Marshall B, et al. Phylogenetic clustering of hepatitis C virus among people who inject drugs in Vancouver, Canada. Hepatology. 2014;60: 1571–80. pmid:25042607
  41. 41. Jacka B, Applegate T, Poon AF, Raghwani J, Harrigan PR, DeBeck K, et al. Transmission of hepatitis C virus infection among younger and older people who inject drugs in Vancouver, Canada. J Hepatol. 2016;64: 1247–55. pmid:26924451
  42. 42. Gray RR, Parker J, Lemey P, Salemi M, Katzourakis A, Pybus OG. The mode and tempo of hepatitis C virus evolution within and among hosts. BMC Evol Biol. 2011; 11: 131. pmid:21595904
  43. 43. Paraschiv S, Banica L, Nicolae I, Niculescu I, Abagiu A, Jipa R, et al. Epidemic dispersion of HIV and HCV in a population of co-infected Romanian injecting drug users. PLoS One. 2017;12.
  44. 44. Olmstead AD, Joy JB, Montoya V, Luo I, Poon AF, Jacka B, et al. A molecular phylogenetics-based approach for identifying recent hepatitis C virus transmission events. Infect Genet Evol. 2015;33: 101–9. pmid:25917496
  45. 45. Lemey P, Vandamme AM. Exploring full-genome sequences for phylogenetic support of HIV-1 transmission events. AIDS. 2005;19: 1551–2. pmid:16135914
  46. 46. González-Candelas F, López-Labrador FX, Bracho MA. Recombination in hepatitis C virus. Viruses. 2011;3: 2006–24. pmid:22069526
  47. 47. Ceccherini Silberstein F, Di Maio VC, Aragri M, Ciotti M, Cento V, Perno CF. Hepatitis C virus gene sequencing as a tool for precise genotyping in the era of new direct acting antivirals. Hepatology. 2016;63: 1058–9. pmid:25974082
  48. 48. Chueca N, Rivadulla I, Lovatti R, Reina G, Blanco A, Fernandez-Caballero JA, et al. Using NS5B sequencing for hepatitis C virus genotyping reveals discordances with commercial platforms. PLoS One. 2016;11.
  49. 49. Fernández-Caballero JA, Alvarez M, Chueca N, Pérez AB, García F. The cobas® HCV GT is a new tool that accurately identifies hepatitis C virus genotypes for clinical practice. PLoS One. 2017;12.
  50. 50. Schneider MD, Sarrazin C. Antiviral therapy of hepatitis C in 2014: do we need resistance testing? Antiviral Res. 2014;105: 64–71. pmid:24583028
  51. 51. Abdelrahman T, Hughes J, Main J, McLauchlan J, Thursz M, Thomoson E. Next-generation sequencing sheds light on the natural history of hepatitis C infection in patients who fail treatment. Hepatology. 2015;61: 88–97. pmid:24797101
  52. 52. Robaeys G, Bielen R, Azar DG, Razavi H, Nevens F. Global genotype distribution of hepatitis C viral infection among people who inject drugs. J Hepatol. 2016;65: 1094–1103. pmid:27520879
  53. 53. Esteban JI, Sauleda S, Quer J. The changing epidemiology of hepatitis C virus infection in Europe. J Hepatol. 2008;48: 148–62. pmid:18022726