G C A T T A C G G C A T genes Article Study of Quasispecies Complexity and Liver Damage Progression after Liver Transplantation in Hepatitis C Virus Infected Patients Meritxell Llorens-Revull 1,2,3 , Josep Gregori 1,2,4 , Cristina Dopazo 2,5 , Francisco Rodriguez-Frías 2,3,6 , Damir Garcia-Cehic 1,2 , Maria Eugenia Soria 1 , Qian Chen 1 , Ariadna Rando 6 , Celia Perales 1,2 , Juan Ignacio Esteban 1,2,7 , Josep Quer 1,2,3, * and Itxarone Bilbao 2,5,8 1 Liver Diseases-Viral Hepatitis, Liver Unit, Vall d’Hebron Institut de Recerca (VHIR), Vall d’Hebron Hospital Universitari, Vall d’Hebron Barcelona Hospital Campus, Passeig Vall d’Hebron 119-129, 08035 Barcelona, Spain; meritxell.llorens@vhir.org (M.L.-R.); josep.gregori@gmail.com (J.G.); damir.garcia@vhir.org (D.G.-C.); maria.soria@vhir.org (M.E.S.); qian.chen@vhir.org (Q.C.); cperales@cbm.csic.es (C.P.); jiesteban@vhebron.net (J.I.E.) 2 Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd), Instituto de Salud Carlos III, Av. Monforte de Lemos, 3-5, 28029 Madrid, Spain; cdopazo@vhebron.net (C.D.); frarodri@gmail.com (F.R.-F.); ibilbao@vhebron.net (I.B.) 3 Biochemistry, Molecular Biology, Universitat Autònoma de Barcelona (UAB), Campus de la UAB, Plaça Cívica, 08193 Bellaterra, Spain 4 Roche Diagnostics SL, Avinguda de la Generalitat, 171-173, 08174 Sant Cugat del Vallès, Spain   5 Hepatobiliopancreatic Surgery and Transplant Unit, Vall d’Hebron Institut de Recerca (VHIR), Vall d’Hebron Hospital Universitari, Vall d’Hebron Barcelona Hospital Campus, Passeig Vall d’Hebron 119-129, Citation: Llorens-Revull, M.; 08035 Barcelona, Spain Gregori, J.; Dopazo, C.; 6 Biochemistry and Microbiology Departments, Vall d’Hebron Institut de Recerca (VHIR), Vall d’Hebron Rodriguez-Frías, F.; Garcia-Cehic, D.; Hospital Universitari, Vall d’Hebron Barcelona Hospital Campus, Passeig Vall d’Hebron 119-129, Soria, M.E.; Chen, Q.; Rando, A.; 08035 Barcelona, Spain; ariadna.rando@vhir.org Perales, C.; Esteban, J.I.; et al. Study 7 Medicine, Universitat Autònoma de Barcelona (UAB), Campus de la UAB, Plaça Cívica, of Quasispecies Complexity and 08193 Bellaterra, Spain 8 Surgery, Universitat Autònoma de Barcelona (UAB), Campus de la UAB, Plaça Cívica, 08193 Bellaterra, Spain Liver Damage Progression after Liver Transplantation in Hepatitis C Virus * Correspondence: josep.quer@vhir.org or josep.quer@ciberehd.org Infected Patients. Genes 2021, 12, 1731. https://doi.org/10.3390/ Abstract: Cirrhosis derived from chronic hepatitis C virus (HCV) infection is still a common indication genes12111731 for liver transplantation (LT). Reinfection of the engrafted liver is universal in patients with detectable viral RNA at the time of transplant and causes fast progression of cirrhosis (within 5 years) in around Academic Editor: Salvador F. Aliño one-third of these patients. To prevent damage to the liver graft, effective direct-acting antiviral (DAA) therapy is required as soon as possible. However, because of post-LT clinical instability, it Received: 29 September 2021 is difficult to determine the optimal time to start DAAs with a low risk of complications. Evaluate Accepted: 26 October 2021 changes in quasispecies complexity following LT and seek a predictive index of fast liver damage Published: 28 October 2021 progression to determine the timing of DAA initiation. HCV genomes isolated from pre-LT and 15-day post-LT serum samples of ten patients, who underwent orthotopic LT, were quantified Publisher’s Note: MDPI stays neutral and sequenced using a next-generation sequencing platform. Sequence alignments, phylogenetic with regard to jurisdictional claims in trees, quasispecies complexity measures, biostatistics analyses, adjusted R2 values, and analysis of published maps and institutional affil- iations. variance (ANOVA) were carried out. Three different patterns of reinfection were observed (viral bottlenecking, conserved pre-LT population, and mixed populations), suggesting that bottlenecking or homogenization of the viral population is not a generalized effect after liver graft reinfection. None of the quasispecies complexity measures predicted the future degree of liver damage. Higher and more uniform viral load (VL) values were observed in all pre-LT samples, but values were more Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. dispersed in post-LT samples. However, VL increased significantly from the pre-LT to 15-day post-LT This article is an open access article samples in patients with advanced fibrosis at 1-year post-LT, suggesting that a VL increase on day 15 distributed under the terms and may be a predictor of fast liver fibrosis progression. HCV kinetics after LT differ between patients conditions of the Creative Commons and are not fibrosis-dependent. Higher VL at day 15 post-LT versus pre-LT samples may predict fast Attribution (CC BY) license (https:// liver fibrosis progression. creativecommons.org/licenses/by/ 4.0/). Genes 2021, 12, 1731. https://doi.org/10.3390/genes12111731 https://www.mdpi.com/journal/genes Genes 2021, 12, 1731 2 of 15 Keywords: hepatitis C virus; liver transplantation; fibrosis; variability; complexity measures; viral load 1. Introduction Hepatitis C virus (HCV) infection is a cause of end-stage liver disease and an indication for liver transplantation (LT). The reinfection of the engrafted liver is universal in patients with detectable viral RNA during the first days, and even the first hours, after transplanta- tion [1,2]. Symptomatic HCV hepatitis then develops in 1 to 4 months, although the clinical pattern varies. Within 5 years after LT occurs, there is a fast progression to cirrhosis in 10% to 30% of patients [3–6]. Before the introduction of interferon-free regimens, only 30% of non-transplant patients with HCV-related cirrhosis had liver decompensation at 10 years, whereas more than 40% of graft recipients showed decompensation within 12 months after the diagnosis of recurrent cirrhosis, and up to 60% experienced a decompensation episode 3 years later [4,7–10]. The factors that accelerate the post-LT progression of liver damage in HCV patients are uncertain and seem to depend on the characteristics of the virus and the patient [11,12]. Reports have shown that some amino acid signatures in the NS5B region of HCV are specific to patients developing cholestatic fibrosis hepatitis, which is a severe variant of HCV infection recurrence after liver transplantation [13]. However, there is little information on the diversity index and other viral factors that could help predict the accelerated progression of liver damage. The advent of safe and highly effective direct-acting antiviral agents (DAAs) has profoundly changed the management of patients with advanced liver damage and those undergoing LT. Although DAA treatment before LT would be the best option, drug-to-drug interactions between some DAAs and various immunosuppressive agents may jeopardize this approach. Furthermore, some DAAs should be avoided in patients with severely impaired liver function and renal dysfunction because of complications, the most common being renal failure, which is not unusual after LT. An alternative is treatment soon after the procedure or when the risk of chronic rejection has decreased and immunosuppressive rejection medication is stable [14]. Hence, the question arises as to when would be the best time to start DAA treatment. The ELITA consensus statements [9] summarize the factors that should be considered to determine whether pre- and post-LT DAA therapy is justified in patients listed for decompensated cirrhosis without hepatocellular carcinoma (HCC). In general, patients who will progress faster to cirrhosis (fast progressors) should be treated earlier after LT than those with slower progression (slow progressors), who can be treated later, after the patient’s health status has improved. Therefore, finding a predictive index of fast liver damage progression to determine the timing of DAA initiation is of great interest, and quasispecies composition has been related to viral persistence, disease progression, and response to antiviral agents [15–17]. HCV is an enveloped, positive-sense, single-stranded RNA virus of the genus Hep- acivirus in the Flaviviridae family [18]. The lack of proofreading activity of the nonstruc- tural protein 5B (NS5B) RNA-dependent RNA polymerase leads to substantial sequence variation in the HCV genome, which displays mutation rates in the range of 10-3–10-5 mutations per nucleotide copied [19]. Hence, HCV does not exist as a single genome, even in a single infected individual; instead, it is a complex and dynamic distribution of non- identical but related genomes known as a viral quasispecies that undergoes a continuous process of genetic variation, competition, and selection [20–25]. The viral quasispecies complexity has been defined as the “intrinsic property that quantifies the diversity and frequency of haplotypes, independently of the population size that contains them” [26]. The relevance of the viral complexity level has been evidenced by the lower adaptability of viruses whose polymerase shows a higher or lower copying fidelity than the wild type with a comparable population size in the same biological context [27–30]. Therefore, when studying an infection, it is important to be aware that all members of the quasispecies Genes 2021, 12, 1731 3 of 15 may be relevant for establishing chronic infection, and an approach based on studying the diversity and complexity of the entire viral population should be considered. The aim of this study was to evaluate changes occurring in the complexity of the HCV quasispecies soon after LT, seeking a parameter that would be predictive of fibrosis progression at 1 year following the procedure, which could be useful for guiding the timing of the start of DAA therapy. 2. Materials and Methods This study consisted of ten HCV patients who underwent orthotopic LT, fulfilled the inclusion criteria, and had none of the exclusion criteria, as shown in the Supplementary Material. This study was approved by the local institutional review board for clinical research, and all patients gave written informed consent in accordance with the 1975 Declaration of Helsinki. Ten patients were included in the study, five being infected with HCV genotype 1 subtype a (G1a), four with G1b, and one with G3a [31]. Serum samples were collected 6 weeks before LT at the time they were included in the waiting list and 15 days after LT at the moment of hospital discharge. The histological grade of liver fibrosis was evaluated according to the ISHAK fibrosis score [32] in liver biopsy specimens 1 year after transplantation. Viral load (VL) before LT (6 weeks) and 15 days after LT were measured using the Cobas 6800 system (Roche Applied Science, Basel, Switzerland; lower limit of detection, 10 IU/mL). 2.1. RNA Extraction, RT-PCR, Heminested-PCR Amplification and Quantification A 30 µL amount of total RNA was obtained by extracting between 140 and 280 µL of the 6 weeks pre-LT and 15 days post-LT serum sample (depending on the VL of each) using the QIAamp Viral RNA Mini Kit 250 (Qiagen, Hilden, Germany) and following the manufacturer’s instructions. Samples P03 and P05 with the lowest VL required a double serum volume for RNA extraction to avoid bias occurring in the sequencing studies. Reverse transcription polymerase chain reaction (RT-PCR) was performed as previ- ously described [31]. Briefly, the reverse transcription of the NS5B region (nucleotides (nt) 8254–8707) was carried out using the Transcriptor One-Step reverse RT-PCR kit (Roche Ap- plied Science, Basel, Switzerland). HCV RNA was reverse-transcribed into complementary deoxyribonucleic acid (cDNA) (30 min at 50 ◦ C) and PCR-amplified for 35 cycles (10 s at 94 ◦ C, 30 s at 53 ◦ C, and 30 s at 68 ◦ C) with specific oligonucleotides (5Bu8254: CNTAYGAY- ACCMGNTGYTTTGACTC; 5Bd8707: TTNGADGAGCADGATGTWATBAGCTC). A final product of 454 nt was obtained. Nucleotide positions were marked according to isolate H77 accession number AF009606 [33]. The reaction mixture for RT-PCR was prepared as follows: 28.5 µL H2O-PCR 1×, 10 µL buffer 5×, 2.5 µL DMSO, 1.5 µL upstream primer (20 pmol), 1.5 µL downstream primer (20 pmol), 1 µL polymerase, and 5 µL RNA. Heminested PCR was then performed with 5 µL of DNA from the above PCR using the FastStart High Fidelity PCR System dNTPack kit (Roche Applied Science Basel, Switzer- land) with 20 pmol of the labeled upstream primer (13N5Bo8254: GTTGTAAAACGACG- GCCAGTCNTAYGAYACCMGNTGYTTTGACTC) and 20 pmol of the labeled downstream primer (13N5Bo8641: CACAGGAAACAGCTATGACCGARTAYCTGGTCATAGCNTC- CGTGAA), both of which included a complementary universal M13 primer at the 50 end. HCV DNA was amplified in a 35-cycle PCR (30 s at 94 ◦ C, 30 s at 56 ◦ C, and 30 s at 72 ◦ C) with a reaction mixture comprised of 33 µL H2 O-PCR 1×, 5 µL buffer 10×, 1 µL dNTP, 2.5 µL DMSO, 1.5 µL of 20 pmol upstream primer, 1.5 µL of 20 pmol downstream primer, 0.5 µL polymerase, and 5 µL DNA. A final product of 428 nt (including primers) was obtained. PCR products from different isolates were pooled together before deep sequencing. Each isolate was tagged with a different multiplex identifier (MID) by performing a short (15 cycles) PCR. The final product was 498 nt in length (including primers, MID, and adaptors for GS–Junior sequencing), spanning nucleotides 8254–8641. Genes 2021, 12, 1731 4 of 15 Amplification products were analyzed using 2% agarose gel electrophoresis with the QIAquick gel extraction kit (Qiagen, Valencia, CA, USA), quantified using the PicoGreen assay (Invitrogen, Carlsbad, CA, USA), and quality-analyzed using a BioAnalyzer DNA 1000 LabChip (Agilent, Santa Clara, CA, USA) prior to sequencing. 2.2. Ultradeep Pyrosequencing (UDPS) Purified DNA from each sample was mixed, forming equimolar pools. Each pool was sequenced by UDPS-based next-generation sequencing (NGS) on the GS–Junior platform (454 Life Sciences, Roche, Branford, CT, USA), following the manufacturers’ protocol. 2.3. Nucleotide Haplotypes and Diversity Analyses Sequencing data analysis was performed as previously reported [34,35]. Briefly, sequences were demultiplexed by the MID and specific primer and pairwise-aligned with respect to the dominant haplotype, excluding reads that did not cover the full amplicon or harbored more than 2 indeterminations or 3 gaps. Accepted indeterminations and gaps were repaired as per the dominant haplotype. Reads were collapsed to haplotypes and their corresponding frequencies. Only those with abundances above 0.1% and common to the forward and reverse strands were kept. All computations were made in the R language and platform [36] with in-house-developed scripts and with the help of the Biostrings [37], ape [38], seqinr [39], and ade4 [40] packages. Diversity analyses were carried out on sequences that passed the filters, standard- izing to the lowest coverage with a down-sampling and fringe-trimming approach [35]. The following diversity indices were used to define the viral quasispecies complexity at the molecular level [26]: mutation frequency (Mfmax), Hill numbers (qD) [41,42], and nucleotide diversity (π) [43,44]. 2.4. Phylogenetic Analysis For the phylogenetic analysis, read alignments were further filtered, discarding all haplotypes below 0.5% [35]. Haplotypes from the pre- and post-LT quasispecies were clustered by UPGMA (Unweighted Pair Group Method with Arithmetic mean) on the matrix of Kimura-80 genetic distances [45]. 2.5. Biostatistic Analyses Quasispecies complexity changes before and after LT were studied using the Mann– Whitney U test on each diversity index, with the null hypothesis of equal population means and the alternate hypothesis of greater diversity in pre-LT samples. Diversity indices were computed in the pre-LT and 15 day post-LT sample. A principal component analysis was performed on the matrix of diversity measures, allowing the representation of samples characterized by variables that could be strongly correlated on orthogonal axes. The aim of this multidimensional exploratory analysis was to determine whether the diversity changes showed a general trend from pre-LT to post-LT when all indices were included in a multivariate analysis. Associations between the liver damage level at 1 year post-LT and each diversity index in both the pre-LT and 15-day post-LT samples were studied by fitting a linear model of a single factor with two levels, high fibrosis (F3–F4) and low fibrosis (F0, F1, F2). The adjusted R2 was assessed as a measure of prediction accuracy. Associations between the differences in VLs (log VL at 15 days post-LT minus log VL pre-LT) and liver damage level (high vs. low) at 1 year post-LT was carried out using a single fixed factor analysis of variance (ANOVA) with a significance level of 0.05. The statistical methods used in this study were reviewed by Josep Gregori from Liver Diseases-Viral Hepatitis, Liver Unit, Vall d’Hebron Institut de Recerca (VHIR), Vall d’Hebron Hospital Universitari, Vall d’Hebron Barcelona Hospital Campus, Universitat Autònoma de Barcelona, Passeig Vall d’Hebron 119-129, Barcelona, Spain; Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd), Genes 2021, 12, 1731 5 of 15 Instituto de Salud Carlos III, Madrid, Spain; Roche Diagnostics SL, Sant Cugat del Vallès, Barcelona, Spain. 3. Results Patient data, viral load, and complexity measures are reported in Supplementary Table S1. 3.1. Viral Kinetics and Histological Grade HCV RNA was detected in all patient samples: mean and standard deviation (SD) 1.52 × 106 (SD 703492.881) at pre-LT and 2.91 × 106 (SD 4655095.8) at 15 days post-LT. The following histological diagnoses according to the ISHAK classification were established in liver biopsies acquired at 1 year: F0 in 3 patients, F1 in 1, F2 in 2, F3 in 2, and F4 in 2 others. One of the two patients with F4 fibrosis grade (P10) developed cholestatic fibrosis hepatitis 6 months after liver transplantation. 3.2. Phylogenetic Studies Three different patterns were obtained in the comparison of the phylogenetic analyses of pre- and post-LT samples (Figure 1). In the first pattern, observed in 7 of the 10 patients (70%), the master sequence (sequence with the highest percentage of reads) remained before and after LT. In this group, five of the patients were subtype 1a and two were subtype 1b. Interestingly, the master sequence was present in a higher percentage before than after LT in 4 of the 7 patients, remarking that 3 patients were HCV subtype 1a, while only 1 was 1b. In patient P10 who had a cholestatic event at six months after LT, the master sequence remained dominant after LT. In the analysis of the quasispecies in these patients, we found that most of the individual haplotypes were not common to both (pre- and post-LT) populations; instead, some specifically shared haplotypes had successfully developed new subpopulations (Figure 1A and Supplementary Figures S1–S6). The second pattern, seen in 2 of the 10 patients, showed a clear viral “bottleneck” effect. The pre-LT master sequence had changed after LT, and two clearly separate quasispecies with no or only one shared haplotype were observed (Figure 1B and Supplementary Figures S1–S7). Both patients were HCV subtype 1b. Finally, in the third pattern, observed in 1 of the 10 patients, the master sequence was highly conserved in both the pre- and post-LT quasispecies, with greater representation in the pre-LT samples. The vast majority of haplotypes were also conserved in the two populations and present in similar percentages (Figure 1C). The patient showed that the unique pattern corresponded with genotype 3a. 3.3. Unidimensional Analyses For each diversity measure, a Mann–Whitney U-test was carried out, in which the null hypothesis was that diversity changes were random after transplantation, and, therefore, the average difference was 0. The alternative hypothesis was that diversity was reduced after transplantation due to the homogenization of the quasispecies, and the average difference was greater than 0. The differences (pre/post-LT) for each diversity measure with the resulting p-value from the Mann–Whitney U-test are shown in Figure 2A–F. In all cases, the 0 was well inside the observed distribution of differences and the p-values obtained were well above 0.05. Genes 2021, 12, 1731 6 of 15 Genes 2021, 12, x FOR PEER REVIEW 6 of 15 Figure 1. Cont. Genes 2021, 12, 1731 7 of 15 Genes 2021, 12, x FOR PEER REVIEW 7 of 15 Figure 1. Cont. Genes 2021, 12, 1731 8 of 15 Genes 2021, 12, x FOR PEER REVIEW 8 of 15 Figure 1. Phylogenetic trees Figurerepresenting the liver 1. Phylogenetic treesquasispecies representingpopulation before (pre-LT) the liver quasispecies and after population (post-LT) before liver (pre-LT) andtrans- after plantation. Haplotypes(post-LT) from theliver pre-LT quasispeciesHaplotypes transplantation. are represented in pre-LT from the green, quasispecies and haplotypes from the post-LT are represented in green,qua- and sispecies are in orange, haplotypes as examplesfrom of the three infection patterns found. (A) Patient 05; (B) Patient 08; (C) Patient the post-LT quasispecies are in orange, as examples of the three infection patterns 01. The nomenclature used wasfound.Pxx.y_z_vvvv, where (A) Patient xxPatient 05; (B) is the patient 08; (C)identifier Patient 01.(P01, The P02 etc), y is the nomenclature sample used identifier (0 forwhere was Pxx.y_z_vvvv, pre- LT, 1 for post-LT), z is the number of differences between the haplotype and the master sequence, and vvvv is an identifier xx is the patient identifier (P01, P02 etc), y is the sample identifier (0 for pre-LT, 1 for post-LT), z is for haplotypes with the same number of mutations in the quasispecies. The last percent number represents the frequency the number of differences between the haplotype and the master sequence, and vvvv is an identifier at which each haplotype is represented in the quasispecies. for haplotypes with the same number of mutations in the quasispecies. The last percent number represents the frequency at which each haplotype is represented in the quasispecies. fore, the average difference was 0. The alternative hypothesis was that diversity was re- duced after transplantation due to the homogenization of the quasispecies, and the aver- age difference was greater than 0. The differences (pre/post-LT) for each diversity measure with the resulting p-value Genes 2021, 12, 1731 from the Mann–Whitney U-test are shown in Figure 2A–F. In all cases, the 0 was9 of well 15 inside the observed distribution of differences and the p-values obtained were well above 0.05. Figure 2. Figure 2. Boxplots Boxplotswith withdiversity diversityindices showing indices showing differences between differences betweenthe the pre-LT and and pre-LT post-LT qua- post-LT sispecies. (A) D0: hill number of order 0, (B) D1: hill number of order 1, (C) D2: hill number quasispecies. (A) D0: hill number of order 0, (B) D1: hill number of order 1, (C) D2: hill number of of order 2, (D) Dinf: hill number of order infinity, (E) Mfmax: mutation frequency, (F) π: nucleotide diversity. order 2, (D) Dinf: hill number of order infinity, (E) Mfmax: mutation frequency, (F) π: nucleotide Each patient is represented as *. The p-value resulting from the Mann-Whitney U-test is included. diversity. Each patient is represented as *. The p-value resulting from the Mann-Whitney U-test is included. 3.4. Principal Components Analysis 3.4. Principal Components As in the Analysisanalyses, changes towards lower viral complexity after unidimensional transplantation were not confirmed. As in the unidimensional In Figure analyses, 3A–C, changes a significant towards lowermultivariate association viral complexity after transplantation were not confirmed. In Figure 3A–C, a significant multivariate association between LT and the diversity measures would result in a consistent pattern of arrows pointing in approximately the same direction, which is not the case. 3.5. Associations between Diversity Measures and Fibrosis at 1 Year Post-LT Liver fibrosis as a measure of liver damage was dichotomized into low (F0, F1) and high (F3, F4) values. The distribution of values for the various diversity indices in each liver damage group and in pre-LT and 15 days post-LT samples is shown in Table 1. Here, the adjusted R2 value acts as an assessment of the prediction accuracy. Of note, the adjusted R2 was very low for all of the diversity measures, indicating that they lacked predictive capability in both samples. Genes 2021, 12, x FOR PEER REVIEW 10 of 15 Genes 2021, 12, 1731 10 of 15 between LT and the diversity measures would result in a consistent pattern of arrows pointing in approximately the same direction, which is not the case. Figure3.3.Samples Figure Samplesrepresented representedononthe theplanes planesofofthe thefirst firstthree threeprincipal principalcomponents. components.Each Eachpair pairofof samplesisisrepresented samples representedas asan anarrow, arrow,with withpre-LT pre-LTatatthe thetail tailand andpost-LT post-LTatatthe thehead. head.(A) (A)PC2 PC2vs. vs.PC1; PC1; (B)PC3 (B) PC3vs. vs.PC1; PC1;(C) (C)PC3 PC3vs. vs.PC2. PC2. 3.5. Associations between Diversity Measures and Fibrosis at 1 year Post-LT Table 1. Ability of diversity measures (qD of order 0, 1, 2, infinity, Mfmax, and π) to predict high or lowLiver fibrosis fibrosis degreeas at a1 year measure of liver post-LT usingdamage wasand both pre-LT dichotomized into 15 day post-LT low (F0, samples. F1) and Prediction high (F3, accuracy of F4) eachvalues. diversityThe distribution measure of values is represented for theR2. as adjusted various diversity indices in each liver damage group and in pre-LT and 15 days post-LT samples is shown in Table 1. Here, the adjusted R2 value acts as an assessment Diversity of theMeasures prediction accuracy. Of note, the ad- justed R2 was Adjusted very low forAdjusted all of the diversity Adjustedmeasures, indicating Adjusted that they Adjusted Adjusted lacked pre- R 2 (D0) R 2 (D1) R 2 (D2) R 2 (Dinf) R 2 (Mf) R2 (π) dictive capability in both samples. Pre-LT 0.086 0.064 0.118 0.111 0.021 0.027 Table 1. Ability of0.060 Post-LT diversity measures 0.026 (qD of order 0.0190, 1, 2, infinity, 0.025 Mfmax, 0.063 and π) to predict high or 0.051 low fibrosis degree at 1 year post-LT using both pre-LT and 15 day post-LT samples. Prediction accuracy of each diversity measure is represented as adjusted R2. 3.6. Associations between VL and Fibrosis at 1 Year Post-LT Patients showing a VL increase at 15 Diversity Measures days post-LT relative to the pre-LT value de- Adjusted veloped a higher R Adjusted 2 degree of fibrosisR (F3, 2 Adjusted F4) at 1Ryear 2 Adjusted post-LT; R2 therefore, Adjusted they R2 Adjusted could beR2 considered fast progressors. (D0) In contrast, patients (D1) (D2) showing a VL decrease (Dinf) (Mf)after LT could (π) be considered Pre-LT slow progressors, 0.086 as all 0.064 showed minimal 0.118 liver fibrosis 0.111 changes 0.021(F0, F1) at 1 year 0.027 after transplantation. Post-LT 0.060 The VL values 0.026 were high and 0.019 uniform 0.025 in all pre-LT samples 0.063 but were 0.051 more dispersed in the post-LT samples (Figure 4A). In the ANOVA, a statistically signifi- cant association between VL and the progression of liver damage was found (p = 0.0144) (Figure 4B). oped a higher degree of fibrosis (F3, F4) at 1 year post-LT; therefore, they could be consid- ered fast progressors. In contrast, patients showing a VL decrease after LT could be con- sidered slow progressors, as all showed minimal liver fibrosis changes (F0, F1) at 1 year after transplantation. The VL values were high and uniform in all pre-LT samples but Genes 2021, 12, 1731 were more dispersed in the post-LT samples (Figure 4A). In the ANOVA, a statistically 11 of 15 significant association between VL and the progression of liver damage was found (p = 0.0144) (Figure 4B). Viral load Figure 4. Viral loaddifferences differencesbetween betweenpre preand andpost-LT post-LT samples. samples. (A)(A) Arrows Arrows showing showing thethe evolu- evolution of log VL from pre-LT (arrow tails on the left) to 15-day post-LT (arrow heads on the right side ofside tion of log VL from pre-LT (arrow tails on the left) to 15-day post-LT (arrow heads on the right the of theAt box). box). oneAt year one year post-LT, post-LT, patientspatients with highwith high fibrosis fibrosis (F3–F4) (F3–F4) in orange, in orange, patientspatients with lowwith low fibrosis fibrosis in (F0–F2) (F0–F2) green.in(B) green. (B) Scatterplot Scatterplot with differences with differences in at in log VL log15VL at 15 days days post-LT post-LT minus The minus pre-LT. pre- LT. The ANOVA p-value in the comparison of the two levels of liver damage is included. Each ANOVA p-value in the comparison of the two levels of liver damage is included. Each patient is patient is represented as *. represented as *. 4. Discussion 4. Discussion environment in which Liver replacement causes a drastic change in the environment which the virus virus had adapted to proliferate; therefore, changes would be expected to occur in in the the composition composition of the quasispecies following LT. After After the the primary primary HCV target organ organ is is removed, removed, there there is a substantial reduction in the overall VL, and the viral populations remaining in the bloodstream and extra-hepatic reservoirs initiate infection of the new liver graft [1,46]. The number of viral particles (population size) starting the reinfection, together with differences in the liver graft, state of immunosuppression, and patient characteristics, may influence the evolutionary outcome and affect the quasispecies composition [22,47,48]. In this scenario, during graft reinfection, a genetic phenomenon known as “bottleneck- ing” is expected [49] in which the fittest virions are selected in the absence of an effective immune response, thereby increasing their frequency in serum with still limited variabil- ity [50]. In this sense, our phylogenetic studies have revealed three different patterns of viral behavior after liver graft infection. In the first pattern, observed in most patients Genes 2021, 12, 1731 12 of 15 (70%), the master sequence was conserved, but the mutant spectra differed. Only some sequences were the same as the pre-LT ones, suggesting that the reinfection of the liver was caused by the master and the most prevalent pre-LT sequences. In the two patients with the second pattern, the master sequence differed from the pre-LT one, and the mu- tant spectra were completely different, suggesting that reinfection of the liver graft was caused by a minor mutant from the pre-LT quasispecies. Finally, in the third pattern (one patient), the pre-LT quasispecies was maintained after liver transplantation, suggesting massive reinfection of the liver graft by most of the pre-LT virions. Despite pattern 2 and 3 only include patients with subtypes 1b and 3 respectively, the small sample size (two for pattern 2 and one for pattern 3), do not allow us to extract conclusions on the effect of the subtype in the viral quasispecies behavior after LT. However, pattern 1 has been found in patients with subtype 1a and 1b, suggesting that in most of the cases, and independently of the liver damage progression, HCV reinfection follows this particular pattern master sequence, which remains after LT but with different mutant spectra. Interestingly, in the patient that had a cholestatic fibosis hepatitis (P10) master sequence, it remained dominant after LT as previously reported by Gambato et al. in which 62% of cholestatic patients showed the remaining of master sequence compared with the 11% in the patients with mild recurrence [13]. The main limitation of this study is the small sample size. The main reason for the low number of patients included is that most of the patients are efficiently treated with DAAs before and soon after LT, thus limiting the number of samples fulfilling inclusion criteria for our study in which we require patients with at least one year of untreated follow-up. It is widely recognized that liver fibrosis rapidly progresses in some patients, leading to cirrhosis at 1 year post-LT (fast-progressors), whereas others show minimal changes in the transplanted graft at the same time point (slow-progressors). In this situation, we hypothesized that the viral complexity indices might predict the fast or slow progression of liver damage. However, after conducting an exhaustive analysis of changes in the HCV quasispecies complexity measures (qD, Mfmax, and π) obtained before LT and 15 days following LT, none of the measures studied were significantly associated with progression to more aggressive liver damage at 1 year following the procedure. This may be because the viral quasispecies is a highly variable and dynamic population that fluctuates greatly over time, and the results may vary according to the moment at which the sample is collected and analyzed. Hence, we cannot exclude that the analysis of samples taken at a later point after LT might have led to the identification of a complexity index predictive of fast liver damage progression. Interestingly, we observed a significant increase in VL values from pre-LT to 15 days post-LT in all patients who had an advanced stage of fibrosis 1 year after the procedure, in accordance with previous findings [51]. Our results support the notion that the difference in VL before and after LT may be of value in predicting fibrosis progression. Thus, VL changes may be a useful criterion to determine whether to administer DAA treatment as soon as possible post-LT or delay it until the patient is clinically stable. 5. Conclusions To sum up, none of the viral complexity measures studied at 15 days after liver transplantation were significantly associated with liver damage progression at 1 year following the procedure. Three different patterns of liver graft reinfection were observed based on phylogenetic analyses. VL values were high and uniform in all pre-LT samples but were more dispersed in the post-LT samples. An increase in viral load after liver transplantation was associated with fast progression to liver fibrosis and could be an indicator that effective antiviral treatment should be started as soon as possible in these particular patients. Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/genes12111731/s1, Figures S1–S7: Phylogenetic trees representing the liver quasispecies Genes 2021, 12, 1731 13 of 15 population before (pre-LT) and after (post-LT) liver transplantation. Table S1. Clinical data, viral load, and pre/post-LT viral complexity measures. Inclusion and exclusion criteria. Author Contributions: M.L.-R. designed the study, performed the experiments, acquired, analyzed, and interpreted the data, wrote the manuscript and approved the version of the article to be published; J.G. made essential contributions to the study concept and design, developed software tools and helped perform the analysis and interpretation of the data, reviewed the manuscript and approved the version of the article to be published; C.D. substantially contributed to collecting the patient samples, the study concept and design, writing and reviewing the manuscript, and approving the version of the article to be published; F.R.-F. and C.P. critically reviewed the study for intellectual content and approved the manuscript; D.G.-C., M.E.S., Q.C., and A.R. participated equally in performing the research and contributed to the sample collection and storage; J.I.E. directed the study, corrected the manuscript and approved the final version; J.Q. and I.B. played essential roles in the development of this study. Both authors made large contributions to the concept and design of the study, added intellectual content in the interpretation of the data, and participated in writing the manuscript and approving the final version. All authors have read and agreed to the published version of the manuscript. Funding: This study was supported by grants from Instituto de Salud Carlos III cofinanced by the European Regional Development Fund (ERDF) with grant numbers PI19/00533, PI19/00301, Clinical Trial Gov. Identifier: NCT01707849, and from Centro para el Desarrollo Tecnológico Industrial-CDTI of the Spanish Ministry of Economics and Competitiveness (MINECO) grant number, IDI-20200297. C.P. is supported by the Miguel Servet program of Instituto de Salud Carlos III, grant CP14/00121, cofinanced by the ERDF. Astellas Pharma Inc and Novartis Pharma also provided funding for the study, but these companies had had no role in the study design, data collection or analysis, decision to publish, or preparation of the manuscript. Institutional Review Board Statement: The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Hospital Universitari Vall d’Hebrón, de Barcelona, EVL-VHC-HVH.12 / EudraCTnº: 2012-002105-22. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study. Data Availability Statement: The data presented in this study are available on request from the corresponding author. Acknowledgments: The authors thank the Department of Hepatobiliopancreatic Surgery and Trans- plantation for performing sample collection, Judit Carbonell Segura for her expert technical assistance, and Celine Cavallo for English language support and helpful editing suggestions. Conflicts of Interest: Josep Gregori works in Roche Diagnostics. This fact did not change. References 1. Garcia-Retortillo, M.; Xavier, F.; Feliu, A.; Eduardo, M.; Josep, C.; Miquel, N.; Antoni, R.; Rodes, J. Hepatitis C virus kinetics during and immediately after liver transplantation. Hepatology 2002, 35, 680–687. [CrossRef] 2. Wright, T.L.; Donegan, E.; Hsu, H.H.; Ferrell, L.; Lake, J.R.; Kim, M.; Combs, C.; Fennessy, S.; Roberts, J.P.; Ascher, N.L.; et al. Recurrent and acquired hepatitis C viral infection in liver transplant recipients. Gastroenterology 1992, 103, 317–322. [CrossRef] 3. Gane, E.J.; Portmann, B.C.; Naoumov, N.V.; Smith, H.M.; Underhill, J.A.; Donaldson, P.T.; Maertens, G.; Williams, R. Long-term outcome of hepatitis C infection after liver transplantation. N. Engl. J. Med. 1996, 334, 815–820. [CrossRef] [PubMed] 4. Berenguer, M.; Prieto, M.; Rayón, J.M.; Mora, J.; Pastor, M.; Ortiz, V.; Carrasco, D.; San Juan, F.; Burgueño, M.-J.; Mir, J.; et al. Natural History of Clinically Compensated Hepatitis C Virus–Related Graft Cirrhosis After Liver Transplantation. Hepatology 2000, 32, 852–858. [CrossRef] [PubMed] 5. Prieto, M.; Berenguer, M.; Rayón, J.M.; Córdoba, J.; Argúello, L.; Carrasco, D.; García-Herola, A.; Olaso, V.; De Juan, M.; Gobernado, M.; et al. High incidence of allograft cirrhosis in hepatitis C virus genotype 1b infection following transplantation: Relationship with rejection episodes. Hepatology 1999, 29, 250–256. [CrossRef] [PubMed] 6. Thein, H.H.; Yi, Q.; Dore, G.J.; Krahn, M.D. Estimation of stage-specific fibrosis progression rates in chronic hepatitis C virus infection: A meta-analysis and meta-regression. Hepatology 2008, 48, 418–431. [CrossRef] [PubMed] 7. Samuel, D.; Forns, X.; Berenguer, M.; Trautwein, C.; Burroughs, A.; Rizzetto, M.; Trepo, C. Report of the Monothematic EASL Conference on Liver Transplantation for Viral Hepatitis. (Paris, France, January 12–14, 2006). J. Hepatol. 2006, 45, 127–143. [CrossRef] [PubMed] Genes 2021, 12, 1731 14 of 15 8. Gallegos-Orozco, J.F.; Yosephy, A.; Noble, B.; Aqel, B.A.; Byrne, T.J.; Carey, E.J.; Douglas, D.D.; Mulligan, D.; Moss, A.; de Petris, G.; et al. Natural history of post-liver transplantation hepatitis C: A review of factors that may influence its course. Liver Transplant. 2009, 15, 1872–1881. [CrossRef] [PubMed] 9. Belli, L.S.; Duvoux, C.; Berenguer, M.; Berg, T.; Coilly, A.; Colle, I.; Fagiuoli, S.; Saye, K.; Georges, P.P.; Massimo, P.; et al. ELITA consensus statements on the use of DAAs in liver transplant candidates and recipients. J. Hepatol. 2017, 67, 585–602. [CrossRef] 10. Berenguer, M.; Ferrell, L.; Watson, J.; Prieto, M.; Kim, M.; Rayón, M.; Córdoba, J.; Herola, A.; Ascher, N.; Mir, J.; et al. HCV-related fibrosis progression following liver transplantation: Increase in recent years. J. Hepatol. 2000, 32, 673–684. [CrossRef] 11. Berenguer, M.; Lopez-Labrador, F.X.; Wright, T.L. Hepatitis C and liver transplantation. J. Hepatol. 2001, 35, 666–678. [CrossRef] 12. Berenguer, M.; Schuppan, D. Progression of liver fibrosis in post-transplant hepatitis C: Mechanisms, assessment and treatment. J. Hepatol. 2013, 58, 1028–1041. [CrossRef] [PubMed] 13. Gambato, M.; Gregori, J.; Quer, J.; Koutsoudakis, G.; González, P.; Caro-Pérez, N.; García-Cehic, D.; García-González, N.; González-Candelas, F.; Ignacio, E.J.; et al. Hepatitis C virus intrinsic molecular determinants may contribute to the development of cholestatic hepatitis after liver transplantation. J. Gen. Virol. 2019, 100, 63–68. [CrossRef] 14. Crespo, G.; Trota, N.; Londoño, M.C.; Mauro, E.; Baliellas, C.; Castells, L.; Castellote, J.; Tort, J.; Forns, X.; Navasa, M. The efficacy of direct anti-HCV drugs improves early post-liver transplant survival and induces significant changes in waiting list composition. J. Hepatol. 2018, 69, 11–17. [CrossRef] [PubMed] 15. Pawlotsky, J.M. Hepatitis C virus population dynamics during infection. Curr. Top. Microbiol. Immunol. 2006, 299, 261–284. [CrossRef] 16. Perales, C. Quasispecies dynamics and clinical significance of hepatitis C virus (HCV) antiviral resistance. Int. J. Antimicrob. Agents 2020, 56, 105562. [CrossRef] [PubMed] 17. Farci, P.; Strazzera, R.; Alter, H.J.; Farci, S.; Degioannis, D.; Coiana, A.; Peddis, G.; Usai, F.; Serra, G.; Chessa, L.; et al. Early changes in hepatitis C viral quasispecies during interferon therapy predict the therapeutic outcome. Proc. Natl. Acad. Sci. USA 2002, 99, 3081–3086. [CrossRef] [PubMed] 18. Vainio, E.J.; Chiba, S.; Ghabrial, S.A.; Maiss, E.; Roossinck, M.; Sabanadzovic, S.; Suzuki, N.; Xie, J.; Nibert, M.; Lefkowitz, E.J.; et al. ICTV virus taxonomy profile: Partitiviridae. J. Gen. Virol. 2018, 99, 17–18. [CrossRef] 19. Ogata, N.; Alter, H.J.; Miller, R.H.; Purcell, R.H. Nucleotide sequence and mutation rate of the H strain of hepatitis C virus. Proc. Natl. Acad. Sci. USA 1991, 88, 3392–3396. [CrossRef] [PubMed] 20. Vignuzzi, M.; Stone, J.K.; Arnold, J.J.; Cameron, C.E.; Andino, R. Quasispecies diversity determines pathogenesis through cooperative interactions within a viral population. Nature 2006, 439, 344–348. [CrossRef] 21. Holland, J.; Spindler, K.; Horodyski, F.; Grabau, E.; Nichol, S.; VandePol, S. Rapid evolution of RNA genomes. Science 1982, 215, 1577–1585. [CrossRef] [PubMed] 22. Domingo, E.; Sheldon, J.; Perales, C. Viral quasispecies evolution. Microbiol. Mol. Biol. Rev. 2012, 76, 159–216. [CrossRef] [PubMed] 23. Steinhauer, D.A.; Holland, J.J. Rapid Evolution of RNA Viruses. Annu. Rev. Microbiol. 1987, 41, 409–431. [CrossRef] [PubMed] 24. Duarte, E.A.; Novella, I.S.; Weaver, S.C.; Domingo, E.; Wain-Hobson, S.; Clarke, D.K.; Moya, A.; Elena, S.F.; de la Torre, J.C.; Holland, J.J. RNA virus quasispecies: Significance for viral disease and epidemiology. Infect. Agents Dis. 1994, 3, 201–214. [PubMed] 25. Domingo, E.; Sabo, D.; Taniguchi, T.; Weissmann, C. Nucleotide sequence heterogeneity of an RNA phage population. Cell 1978, 13, 735–744. [CrossRef] 26. Gregori, J.; Perales, C.; Rodriguez-Frias, F.; Esteban, J.I.; Quer, J.; Domingo, E. Viral quasispecies complexity measures. Virology 2016, 493, 227–237. [CrossRef] [PubMed] 27. Martell, M.; Esteban, J.I.; Quer, J.; Genescà, J.; Weiner, A.; Esteban, R.; Guardia, J.; Gómez, J. 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. [CrossRef] 28. Bordería, A.; Stapleford, K.; Vignuzzi, M. RNA virus population diversity: Implications for inter-species transmission. Curr. Opin. Virol. 2011, 1, 643–648. [CrossRef] 29. Campagnola, G.; McDonald, S.; Beaucourt, S.; Vignuzzi, M.; Peersen, O.B. Structure-function relationships underlying the replication fidelity of viral RNA-dependent RNA polymerases. J. Virol. 2015, 89, 275–286. [CrossRef] 30. Pfeiffer, J.K.; Kirkegaard, K. Increased Fidelity Reduces Poliovirus Fitness and Virulence under Selective Pressure in Mice. PLoS Pathog. 2005, 1, e11. [CrossRef] [PubMed] 31. Quer, J.; Gregori, J.; Rodríguez-frias, F.; Buti, M.; Madejon, A.; Perez-del-pulgar, S.; Muñoz, J.M.; Cubero, M.; Caballero, A.; Domingo, E.; 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–226. [CrossRef] [PubMed] 32. Ishak, K.; Baptista, A.; Bianchi, L.; Callea, F.; De Groote, J.; Gudat, F.; Denk, H.; Desmet, V.; Korb, G.; MacSween, R.N.M.; et al. Histological grading and staging of chronic hepatitis. J. Hepatol. 1995, 22, 696–699. [CrossRef] 33. Kuiken, C.; Mizokami, M.; Deleage, G.; Yusim, K.; Penin, F.; Shin-I, T.; Charavay, C.; Tao, N.; Crisan, D.; Grande, D.; et al. Hepatitis C databases, principles and utility to researchers. Hepatology 2006, 43, 1157–1165. [CrossRef] Genes 2021, 12, 1731 15 of 15 34. Ramirez, C.; Gregori, J.; Buti, M.; Tabernero, D.; Camos, S.; Casillas, R.; Quer, J.; Esteban, R.; Homs, M.; Rodriguez-Frias, F. A comparative study of ultra-deep pyrosequencing and cloning to quantitatively analyze the viral quasispecies using hepatitis B virus infection as a model. Antivir. Res. 2013, 98, 273–283. [CrossRef] [PubMed] 35. Gregori, J.; Esteban, J.I.; Cubero, M.; Garcia-Cehic, D.; Perales, C.; Casillas, R.; Alvarez-Tejado, M.; Rodríguez-Frías, F.; Guardia, J.; Domingo, E.; et al. Ultra-deep pyrosequencing (UDPS) data treatment to study amplicon HCV minor variants. PLoS ONE 2013, 8, e83361. [CrossRef] [PubMed] 36. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013; pp. 1–12. 37. Pages, H.; Aboyoun, P.; Gentleman, R.; DebRoy, S. Biostrings: String Objects Representing Biological Sequences, and Matching Algorithms; R Package Version 2.24.1; Fred Hutchinson Cancer Research Center: Seattle, WA, USA, 2012. 38. Paradis, E.; Claude, J.; Strimmer, K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics 2004, 20, 289–290. [CrossRef] [PubMed] 39. Charif, D.; Lobry, J.R. SeqinR 1.0-2: A Contributed Package to the R Project for Statistical Computing Devoted to Biological Sequences Retrieval and Analysis; Springer: Berlin/Heidelberg, Germany, 2007; pp. 207–232. 40. Dray, S.; Dufour, A.B. The ade4 package: Implementing the duality diagram for ecologists. J. Stat. Softw. 2007, 22, 1–20. [CrossRef] 41. Hill, M.O. Diversity and evenness: A unifying notation and its consequences. Ecology 1973, 54, 427–432. [CrossRef] 42. Gotelli, N.J.; Hsieh, T.C.; Sander, E.L.; Colwell, R.K. Rarefaction and Extrapolation with Hill Numbers: A Framework for Sampling and Estimation in Species Diversity Studies; ESA: Washington DC, USA, 2014; Volume 84, pp. 45–67. [CrossRef] 43. Nei, M. Molecular Evolutionary Genetics; Columbia University: New York, NY, USA, 1987. 44. Nei, M.; Kumar, S. Molecular Evolution and Phylogenetics; Oxford University: New York, NY, USA, 2000. 45. Hudelot, C.; Gowri-Shankar, V.; Jow, H.; Rattray, M.; Higgs, P.G. RNA-based phylogenetic methods: Application to mammalian mitochondrial RNA sequences. Mol. Phylogenet. Evol. 2003, 28, 241–252. [CrossRef] 46. Dahari, H.; Feliu, A.; Garcia-Retortillo, M.; Forns, X.; Neumann, A.U. Second hepatitis C replication compartment indicated by viral dynamics during liver transplantation. J. Hepatol. 2005, 42, 491–498. [CrossRef] [PubMed] 47. Novella, I.S.; Elena, S.F.; Moya, A.; Domingo, E.; Holland, J.J. Size of genetic bottlenecks leading to virus fitness loss is determined by mean initial population fitness. J. Virol. 1995, 69, 2869–2872. [CrossRef] [PubMed] 48. Domingo, E. Virus as Populations: Composition, Complexity, Dynamics and Biological Implications; Academic Press: London, UK, 2016. 49. Pérez-Del-Pulgar, S.; Gregori, J.; Rodríguez-Frías, F.; González, P.; García-Cehic, D.; Ramírez, S.; Casillas, R.; Domingo, E.; Esteban, J.I.; Forns, X.; et al. Quasispecies dynamics in hepatitis C liver transplant recipients receiving grafts from hepatitis C virus infected donors. J. Gen. Virol. 2015, 96, 3493–3498. [CrossRef] [PubMed] 50. Farci, P.; Shimoda, A.; Coiana, A.; Diaz, G.; Peddis, G.; Melpolder, J.C.; Strazzera, A.; Chien, D.Y.; Munoz, S.J.; Balestrieri, A.; et al. The Outcome of Acute Hepatitis C Predicted by the Evolution of trhe Viral Quasispecies. Science 2000, 288, 339–344. [CrossRef] [PubMed] 51. Massaguer, A.; Ramírez, S.; Carrión J a González, P.; Sánchez-Tapias, J.M.; Forns, X. Evolution of the NS3 and NS5B regions of the hepatitis C virus during disease recurrence after liver transplantation. Am. J. Transplant. 2007, 7, 2172–2179. [CrossRef] [PubMed]