- Research article
- Open Access
- Open Peer Review
Genetic diversity of NS5A protein from hepatitis C virus genotype 3a and its relationship to therapy response
BMC Infectious Diseasesvolume 10, Article number: 36 (2010)
The quasispecies nature of HCV may have important implications for viral persistence, pathogenicity and resistance to antiviral agents. The variability of one of the viral proteins, NS5A, is believed to be related to the response to IFN therapy, the standard treatment for infection. In this study we analyzed the quasispecies composition of NS5A protein in patients infected with HCV genotype 3a, before IFN therapy.
Viral RNA was isolated from samples of 12 patients: four sustained virological responders (SVR), four non-responders (NR), and four end-of-treatment responders (ETR). cDNA was synthesized, the NS5A region was amplified and the fragments obtained were cloned. Fifteen clones from each patient were sequenced with eight primers, generating 179 contigs.
Higher values for substitution (either synonymous or non-synonymous) and for distance were found in the SVR group. However, the NR group showed relatively more non-synonymous mutations than the other groups, owing to the higher values of dN/dS in complete NS5A and most specific regions. Overall, NS5A protein is undergoing purifying selection, since all dN/dS ratios values are below 0.5.
Our study provides an overview of the genetic variability of complete NS5A protein in HCV genotype 3a.
The hepatitis C virus (HCV) is among the most successful of all persistent human viruses . It is estimated that the global prevalence of HCV infection is 2.2%, corresponding to about 130 million HCV-positive persons worldwide . The HCV genome consists of a single-strand positive-sense RNA of approximately 9.6 Kb that contains an open reading frame coding for a polyprotein precursor of approximately 3000 residues. This precursor is cleaved by viral and host proteinases into the viral proteins: the structural protein core, E1, E2 and p7, and the nonstructural proteins NS2, NS3, NS4A, NS4B, NS5A and NS5B .
To date, six main HCV genotypes have been identified, which differ by about 30% in their nucleotide and amino acid sequences . Genotypes 1, 2 and 3 and their subtypes have a global distribution; genotype 4 is found in Africa, genotype 5 in South Africa, and genotype 6 mainly in Asia . In Brazil, Campiotto et al. (2005) reported the presence of genotypes 1, 2, 3, 4 and 5 .
The RNA-dependent RNA polymerase encoded by the NS5B gene is error-prone and lacks proofreading. As a result, base changes are introduced randomly into the viral genome . Therefore, HCV replication is associated with a high mutation rate, giving rise to a mixed and changing population of mutants known as quasispecies [4, 8]. The quasispecies nature of HCV may have important implications for viral persistence, pathogenicity and resistance to antiviral agents [8–11]. This is most problematic for the infected patient, because quasispecies variation confers remarkable adaptive potential on HCV and has been implicated in the evasion and control of the host response to infection and in differential sensitivity to IFN therapy. The hostile antiviral host environment may drive the proliferation of HCV "evasion variants" from a pre-existing quasispecies pool or through viral genetic adaptation . However, intra-genotype analysis of this diversity in the viral genome shows different degrees of variation; regions such as the 5'UTR and the core are highly conserved, the non-structural regions 2, 3, 5b and the 3'UTR are relatively variable, while the envelope regions E1 and E2 and the NS4 and the NS5A genes exhibit the highest sequence diversity . Sequence analysis of the HCV NS5A coding region has similarly identified specific domains that exhibit sequence variation associated with the outcome of IFN therapy .
NS5A is the nonstructural HCV protein most frequently reported to be implicated in interferon resistance. It is a pleiotropic protein, involved both in viral replication and in many interactions with cellular signaling pathways, including the interferon anti-viral pathway . This study analyzed the NS5A quasispecies pattern in patients infected with HCV genotype 3 before IFN plus ribavirin therapy, with the aim of elucidating its molecular constitution and relationship to treatment response.
The viral load data show that all patients presented high viral loads, as expected for pre-treatment samples, ranging from log 5.76 to log 7.03. There was no correlation between viral load and treatment response (Table 1).
This study generated 15 contig sequences of the full NS5A region from 11 patients; 14 sequences from patient RF80 could be sequenced, totaling 179 contigs.
Quasispecies analysis revealed that the HCV NS5A region was highly variable, as only three patients showed two identical nucleotide sequences (Figure 1A). The amino acid data showed more identical sequences owing to synonymous substitution, as can be seen when the nucleotide and amino acid patterns of patient RF15 are compared. All nucleotide sequences from this patient were different; however, the amino acid sequences showed one of the highest degrees of conservation (Figure 1B).
Nucleotide substitution analysis showed that the sustained virological response (SVR) group had the highest means of nucleotide and amino acid substitutions, except for the amino acid sequences of regions NLS and V3, where the end-of-treatment response (ETR) group showed the highest values (Table 2). However, none of the differences between groups were statistically significant.
To identify the mutation sites in the sequences used in this study, they were represented graphically according to the reference sequence NZL1 (GeneBank D17763) (Figure 2). This representation indicates that no specific mutation could be associated with any kind of treatment response. Figure 2 also shows that some of the nucleotide mutations resulted in stop codons. Of the nine stop codons found in the 179 sequences generated in this study, the site was the same in two or more clones in seven cases. The same stop codon sites were found in NS5A amino acids 4 (RF31 and RF145), 166 (RF60 - 2 clones and RF145) and 447 (RF07 and RF109). Also, the sites where the translation stop codon was observed showed no other mutation, except for one clone from patient RF145, which showed a mutation in aa 447.
For the within-sample (ws) data, the mean rate of synonymous substitution per synonymous site (dS) was higher in the SVR group (Table 2). The rate of non-synonymous substitution per non-synonymous site (dN) was also higher in the SVR group, except for regions CRS and PKR-Bd. These data showed no statistically significant differences.
Between-sample (bs) analysis also showed a higher dS mean in the SVR group. SVR had the highest mean dN values, except for regions NLS and V3. All the bs data showed statistical significance, except for the dS means in CRS between SVR and ETR and between NR and ETR.
In the ws analysis, the dN/dS ratio was higher in the NR group, except for the NLS region, where SVR had the highest ratio. In this analysis, only the difference between SVR and NR in the CRS region was statistically significant.
The bs data showed higher dN/dS ratios in NR for most of NS5A; regions CRS and NLS showed higher values in SVR and ETR, respectively. Statistical significance was observed in all regions, except for PKR-Bd between SVR and NR and for region V3 between SVR and ETR.
The mean genetic distances for the complete NS5A region and the other regions studied were higher in the SVR group in either ws or bs analysis (Table 2). No statistically significant difference was found for genetic distance.
Phylogenetic analysis was carried out using the 179 1356-bp sequences generated in this study, the reference sequence for genotype 3, NZL1 (GenBank accession number D17763), the six full-length NS5A sequences from genotype 3a, with country information, available in GenBank (Accession numbers: AY956467; DQ430819; DQ430820; DQ437509; X76918; GQ300882.1), and 20 Brazilian NS5A sequences of 1308 bp (Accession numbers: EF207999.1 to EF208018.1). The resulting phylogenetic tree is presented in Figure 3. All sequences from clones from the same patient are grouped in a monophyletic branch, with high bootstrap values (95-100%). There was no clustering in a monophyletic branch of the sequences from patients in the same treatment response group.
RNA viruses have high mutation rates owing to lack of proofreading of the RNA polymerase. On the one hand, these high mutation rates can be deleterious for the virus by creating non-viable strains; on the other, when the viral population is considered, they can also create a pool of potentially good mutations. The capacity of HCV to circulate as a pool of different but closely related strains, commonly called quasispecies, allows it to evolve and adapt to new environments and to overcome the challenges of infection [4, 8, 14–16]. The phylogenetic tree constructed in this work illustrates this characteristic and shows, for each patient, very similar but not identical strains with a monophyletic origin.
Our investigation of the genetic variability of the NS5A region from HCV 3a demonstrated the high variability of strains found in the same patient. The highest values of substitutions, either synonymous or non-synonymous, and the highest virus diversity, evaluated by genetic distance, were found in patients with sustained virological responses (SVR). Studies on genotype 1b have reported that the number of mutations in a specific region of NS5A, ISDR, could predict an SVR [17–21]. Since other studies refute these results, a meta-analysis was performed with American, European and Japanese sequences [22–30] that allowed a positive correlation to be established between mutations in the ISDR region and SVR. The same study also revealed a geographic correlation; mutations in the Japanese sequences had a higher probability of leading to SVR than mutations in European sequences. The authors suggested that geographic variation may be related to host factors. This meta-analysis included sequences of genotype 1, but no South American sequences were considered. Since the European and Japanese strains yielded different results, other populations may also present a different pattern.
Although the SVR group showed higher mutation values in the region NS5A by either bs or ws analysis in our study, the dN/dS ratio was higher in the group of non-responders (NR). Therefore, though SVR showed more mutations, NR strains had relatively more non-synonymous substitutions. Amino acid mutations are most commonly deleterious because of the changes they cause in the protein phenotype. However, some of these substitutions are neutral when the mutation does not affect protein function, and are maintained by genetic drift. The higher rates of non-synonymous mutations detected in NR patients may indicate an advantage in evading both the immune system and the treatment, since they may be able to modify epitopes so that they are no longer recognized by the immune system.
For any given genome, the mutation rate determines the ability of a virus to maintain essential information while coping with environmental changes [16, 31–39]. All values obtained for comparisons between dN and dS using the dN/dS ratio were below 0.5. These findings indicate that the NS5A protein is undergoing purifying selection, which maintains protein function and therefore virus activity.
The occurrence of stop codons in the same site has been described previously in the literature. In a similar study, performed on patients infected with HCV genotype 1, two stop codons were found in the same site of the CRS region in two different patients . Also, studies with the Dengue Virus, another Flavivirus, have found the same stop codons in 23% of samples studied. These strains circulated in the population studied for 18 months. These studies suggest that defective genomes can circulate with the assistance of non-defective strains [41, 42]. A recent study with Hepatitis C Virus demonstrated that as the virus circulates in the blood, defective genomes can be encapsidated as infectious particles by trans complementation (acting on a viral RNA other than the one from which it has been translated) of the structural proteins. NS5A protein was also shown to be the only non-structural protein that acts in trans in HCV as well as other Flavivirus, so it could be acting as a helper to defective genomes [43–45].
Of the two analytical approaches used in this study, within-sample (ws) analysis is relevant when the focus of the study is the patient profile, to study a population, in this case a response group; the statistical power is very low and no or few assumptions can be made. However, between-sample (bs) data represent the response group value better, since it is obtained using all 60 sequences from each group. This approach is more accurate when the aim is to identify a pattern from the response group. The bs analyses normalize characteristics that are exclusive to individual patients, by using a pool of strains from the same response group. In most cases they yielded statistical significance. Thus, even though the population size of our study is relatively low (12 patients), statistical significance was obtained, supporting the data on the response group patterns discussed.
This is the first study on the quasispecies composition of the complete NS5A region of HCV genotype 3a. Although we found differences among the response groups, other studies are necessary for a better understanding of the relationship between the variability of this region and the response to treatment with interferon and ribavirin, since most studies are performed with genotype 1, and the genomic differences among these genotypes are significant.
Population and samples
Plasma samples were collected from 12 HCV genotype 3-positive patients, seen at the Hemocenter of São José do Rio Preto, State of São Paulo, Brazil. All samples were collected before treatment, and a 6-months follow-up provided treatment response data. Four patients were sustained virological responders (SVR), i.e the virus was not detected after treatment or during the 6-months follow-up; four were non-responders (NR), since they showed no virological response; and four were end-of-treatment responders (ETR), i.e. a virological response was detected, but at the 6-months follow-up there was a rebound. All patients were infected with HCV genotype 3a. Treatment consisted of INF-α and ribavirin administration for 24 weeks. Patients with a history of alcoholism or infection with another agent that could cause liver damage were excluded. This study was approved by the ethics committee of the Base Hospital of São José do Rio Preto, and all participants signed an informed consent.
The viral load was quantified by Cobas TaqMan HCV Test.
Extraction of HCV RNA and amplification of the NS5A region
Viral RNA was extracted from blood serum samples using a QIAamp Viral RNA Mini Kit (QIAgen), and cDNA was synthesized using a High-Capacity cDNA Archive Kit (Applied Biosystems) according to the manufacturer's instructions. The NS5A region was amplified using a set of primers specific for genotype 3, designed for this study. For the PCR reaction, two sets of forward and reverse primers were designed. The forward primers were H.NS5AP-F (5' GAGCGGTACAGTGGATGAAC 3' - nucleotide [nt] 6089 to 6108 in genotype 3 reference sequence NZL1 - GenBank D17763G) and H.NS5AP-F2 (5' GGTACAGTGGATGAACAGG 3' - nt 6093 to 6111 in NZL1). The reverse primers were H.NS5AP-R (5' CCTCCTTTAATGCAGTCTTG 3' - nt 7821 to 7840 in NZL1) and H.NS5AP-R2 (5' ACGACGTTGAATAGACTAGG 3' nt 7734 to 7753 in NZL1). Two sets of forward and reverse primers were also designed for the nested-PCR reaction. The forward primers were H.NS5AN-F (5' CGCATTGCTGAGTTCTCTAAC 3' from nt 6192 to 6212 in NZL1) and H.NS5AN-F2 (5' CTCTAACTGTCACAAGTCTGC 3' nt 6206 to 6226 in NZL1). The reverse primers were H.NS5AN-R (5' CAACAAGGAGTTGCTGAGTG 3' nt 7703 to 7722 in NZL1) and H.NS5AN-R2 (5' CAGCACTACATGGTGTTATC 3' nt 7659 to 7678 in NZL1). For the amplification reaction, 1 U of a proofreading polymerase was used (Elongase® Enzyme Mix; Invitrogen) along with 10 μl of buffer B [300 mM Tris-SO4, (pH 9.1 at 25°C), 90 mM (NH4)2SO4 and 10 mM MgSO4], 10 μl of dNTP, 30 pmol of sense and anti-sense primers, 10 μl of cDNA for PCR reaction, and 5 μl of PCR product for the nested-PCR reaction, plus Milli-Q autoclaved water to a final volume of 50 μl. The amplified products were analyzed on a 1% agarose gel.
Cloning and sequencing
Cloning was performed using a TOPO XL Cloning TM Kit (Invitrogen). Fragments from 15 clones from each patient were purified using a PureLink Mini-prep Plasmid Purification Kit (Invitrogen). The entire NS5A region was sequenced using eight primers: the vector primers M13F and M13R, and six inner primers, three sense and three anti-sense, designed for this study. The forward primers used in the sequencing reaction were: H.NS5AI-F1 (5' TGGCTGCGTATCATCTGGGA 3' - nt 6283 to 6302 in NZL1), H.NS5AI-F2 (5' ACCTCGATGTTGAGAGACCC 3' - nt 6871 to 6890 in NZL1) and H.NS5AI-F3 (5' TATCCTCCAGCCCTTCCTAT 3' - nt 7198 to 7217 in NZL1). The reverse primers used in the sequencing reaction were: H.NS5AI-R1 (5' CACGGACACTTGAGCTCATC 3' - nt 6679 to 6698 in NZL1), H.NS5AI-R2 (5' TTCTTGAAACACTCTGCAGC 3' - nt 7168 to 7187 in NZL1) and H.NS5AI-R3 (5' GTGGACCAAGAGTCGCAACT 3' - nt 7573 to 7592 in NZL1). The sequencing reaction was performed with Dyenamic ET Terminator (GE) and the products were sequenced in an ABI Prism 377 sequencer (Applied Biosystems). The reaction mixture consisted of 1 μl of Milli-Q autoclaved water, 1 μl of primer (5 pmol/μl), 2 μl of sequencing reagent mix, plus 2 μl of sample. Occasionally, when a good quality sequence could not be obtained, the sequencing reaction had to be doubled and a "hot start" had to be performed for 10 min at 95°C for better results. Cycling was carried out according to the manufacturer's instructions.
Sequence and phylogenetic analysis
The sequences obtained were subjected to BioMol - Electropherogram quality analysis http://adenina.biomol.unb.br/phph/, a phred phrap [47, 48] analysis site, for quality check and contig construction. The contigs obtained for each clone were aligned, along with the reference sequence NZL1 for genotype 3 (GenBank accession number D17763), using Clustal X 1.81 software . All sequences were edited on Bio Edit 184.108.40.206  to remove the vector fragments, leaving only the complete sequence of the NS5A region.
Quasispecies analysis was carried out using LOCSPEQ 1.0 software , specially designed for our group for this kind of analysis.
The number of mutations and the genetic distances were calculated using MEGA 4.0 software . The rates of synonymous substitution per synonymous site (dS) and non-synonymous substitution per non-synonymous site (dN), as well as the dN/dS ratio, were obtained at SNAP - Synonymous Non-synonymous Analysis Program -http://www.hiv.lanl.gov[53, 54].
Phylogenetic analysis was performed using PAUP* version 4 software . A neighbor-joining phylogenetic tree was constructed with Tamura-Nei's substitution model including invariant sites (I) and Gamma distribution shape (G) parameter (TRN+I+G), determined by hierarchical likelihood ratio test score criteria using Modeltest 3.06 . One thousand replicates were used to test the reliability of the tree topology, and bootstrap values >70 were considered significant .
Statistical analysis was performed by ANOVA. Comparisons between groups were made using Tukey's method for multiple comparisons. Values of P < 0.05 were considered significant. Standard error of the mean (SEM) values were calculated in Minitab 15.
Types of analysis
In this work, we chose to analyze the results obtained by two different approaches:
Within-sample analysis (ws)
This consisted of calculating the means of genetic distances, dS, dN and dN/dS ratios among the clones of one patient, and then calculating the mean of the four values obtained for each patient in the response group to obtain the group value.
Between-sample analysis (bs)
This consisted of calculating the means of genetic distances, dS, dN and dN/dS ratios among all the clones of the response group to obtain the group values.
Nucleotide sequence accession numbers
The nucleotide sequence data reported here have been submitted to the GenBank nucleotide sequence database with accession numbers from EU826174 to EU826352.
Bowen DG, Walker CM: Adaptive immune responses in acute and chronic hepatitis C virus infection. Nature. 2005, 436 (7053): 946-952. 10.1038/nature04079.
Alter MJ: Epidemiology of hepatitis C virus infection. World J Gastroenterol. 2007, 13 (17): 2436-2441.
Choo QL, Kuo G, Weiner AJ, Overby LR, Bradley DW, Houghton M: Isolation of a cDNA clone derived from a blood-borne non-A, non-B viral hepatitis genome. Science. 1989, 244 (4902): 359-362. 10.1126/science.2523562.
Martell M, Esteban JI, Quer J, Genesca J, Weiner A, Esteban R, Guardia J, Gomez 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 (5): 3225-3229.
Silva LK, Parana R, Souza SP, Berby F, Kay A, Trepo C, Santana N, Cotrim H, Lyra LG, Reis MG: Hepatitis C virus genotypes in a northeastern area of Brazil. Am J Trop Med Hyg. 2000, 62 (2): 257-260.
Campiotto S, Pinho JR, Carrilho FJ, Da Silva LC, Souto FJ, Spinelli V, Pereira LM, Coelho HS, Silva AO, Fonseca JC, et al: Geographic distribution of hepatitis C virus genotypes in Brazil. Braz J Med Biol Res. 2005, 38 (1): 41-49. 10.1590/S0100-879X2005000100007.
Bukh J, Miller RH, Purcell RH: Genetic heterogeneity of hepatitis C virus: quasispecies and genotypes. Semin Liver Dis. 1995, 15 (1): 41-63. 10.1055/s-2007-1007262.
Domingo E, Escarmis C, Sevilla N, Moya A, Elena SF, Quer J, Novella IS, Holland JJ: Basic concepts in RNA virus evolution. Faseb J. 1996, 10 (8): 859-864.
Forns X, Bukh J: The molecular biology of hepatitis C virus. Genotypes and quasispecies. Clin Liver Dis. 1999, 3 (4): 693-716. 10.1016/S1089-3261(05)70234-8.
Forns X, Purcell RH, Bukh J: Quasispecies in viral persistence and pathogenesis of hepatitis C virus. Trends Microbiol. 1999, 7 (10): 402-410. 10.1016/S0966-842X(99)01590-5.
Farci P, Strazzera R, Alter HJ, 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 (5): 3081-3086. 10.1073/pnas.052712599.
Le Guillou-Guillemette H, Vallet S, Gaudy-Graffin C, Payan C, Pivert A, Goudeau A, Lunel-Fabiani F: Genetic diversity of the hepatitis C virus: impact and issues in the antiviral therapy. World J Gastroenterol. 2007, 13 (17): 2416-2426.
Gale M, Foy EM: Evasion of intracellular host defence by hepatitis C virus. Nature. 2005, 436 (7053): 939-945. 10.1038/nature04078.
Eigen M, Biebricher C: Sequence space and quasispecies distribution. 1988, Boca Raton, FL: CRC Press, 3:
Eigen M: Viral quasispecies. Sci Am. 1993, 269 (1): 42-49. 10.1038/scientificamerican0793-42.
Coffin JM: HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy. Science. 1995, 267 (5197): 483-489. 10.1126/science.7824947.
Enomoto N, Sakuma I, Asahina Y, Kurosaki M, Murakami T, Yamamoto C, Ogura Y, Izumi N, Marumo F, Sato C: Mutations in the nonstructural protein 5A gene and response to interferon in patients with chronic hepatitis C virus 1b infection. N Engl J Med. 1996, 334 (2): 77-81. 10.1056/NEJM199601113340203.
Chayama K, Tsubota A, Kobayashi M, Okamoto K, Hashimoto M, Miyano Y, Koike H, Kobayashi M, Koida I, Arase Y, et al: Pretreatment virus load and multiple amino acid substitutions in the interferon sensitivity-determining region predict the outcome of interferon treatment in patients with chronic genotype 1b hepatitis C virus infection. Hepatology. 1997, 25 (3): 745-749. 10.1002/hep.510250342.
Arase Y, Ikeda K, Chayama K, Murashima N, Tsubota A, Suzuki Y, Saitoh S, Kobayashi M, Kobayashi M, Kobayashi M, et al: Efficacy and changes of the nonstructural 5A GENE by prolonged interferon therapy for patients with hepatitis C virus genotype 1b and a high level of serum HCV-RNA. Intern Med. 1999, 38 (6): 461-466. 10.2169/internalmedicine.38.461.
Kurosaki M, Enomoto N, Murakami T, Sakuma I, Asahina Y, Yamamoto C, Ikeda T, Tozuka S, Izumi N, Marumo F, et al: Analysis of genotypes and amino acid residues 2209 to 2248 of the NS5A region of hepatitis C virus in relation to the response to interferon-beta therapy. Hepatology. 1997, 25 (3): 750-753. 10.1002/hep.510250343.
Nakano I, Fukuda Y, Katano Y, Nakano S, Kumada T, Hayakawa T: Why is the interferon sensitivity-determining region (ISDR) system useful in Japan?. J Hepatol. 1999, 30 (6): 1014-1022. 10.1016/S0168-8278(99)80254-2.
Berg T, Mas Marques A, Hohne M, Wiedenmann B, Hopf U, Schreier E: Mutations in the E2-PePHD and NS5A region of hepatitis C virus type 1 and the dynamics of hepatitis C viremia decline during interferon alfa treatment. Hepatology. 2000, 32 (6): 1386-1395. 10.1053/jhep.2000.20527.
Duverlie G, Khorsi H, Castelain S, Jaillon O, Izopet J, Lunel F, Eb F, Penin F, Wychowski C: Sequence analysis of the NS5A protein of European hepatitis C virus 1b isolates and relation to interferon sensitivity. J Gen Virol. 1998, 79 (Pt 6): 1373-1381.
Sarrazin C, Berg T, Lee JH, Ruster B, Kronenberger B, Roth WK, Zeuzem S: Mutations in the protein kinase-binding domain of the NS5A protein in patients infected with hepatitis C virus type 1a are associated with treatment response. J Infect Dis. 2000, 181 (2): 432-441. 10.1086/315263.
Squadrito G, Orlando ME, Cacciola I, Rumi MG, Artini M, Picciotto A, Loiacono O, Siciliano R, Levrero M, Raimondo G: Long-term response to interferon alpha is unrelated to "interferon sensitivity determining region" variability in patients with chronic hepatitis C virus-1b infection. J Hepatol. 1999, 30 (6): 1023-1027. 10.1016/S0168-8278(99)80255-4.
Zeuzem S, Lee JH, Roth WK: Mutations in the nonstructural 5A gene of European hepatitis C virus isolates and response to interferon alfa. Hepatology. 1997, 25 (3): 740-744. 10.1002/hep.510250341.
Chung RT, Monto A, Dienstag JL, Kaplan LM: Mutations in the NS5A region do not predict interferon-responsiveness in american patients infected with genotype 1b hepatitis C virus. J Med Virol. 1999, 58 (4): 353-358. 10.1002/(SICI)1096-9071(199908)58:4<353::AID-JMV6>3.0.CO;2-C.
Nousbaum J, Polyak SJ, Ray SC, Sullivan DG, Larson AM, Carithers RL, Gretch DR: Prospective characterization of full-length hepatitis C virus NS5A quasispecies during induction and combination antiviral therapy. J Virol. 2000, 74 (19): 9028-9038. 10.1128/JVI.74.19.9028-9038.2000.
Murphy MD, Rosen HR, Marousek GI, Chou S: Analysis of sequence configurations of the ISDR, PKR-binding domain, and V3 region as predictors of response to induction interferon-alpha and ribavirin therapy in chronic hepatitis C infection. Dig Dis Sci. 2002, 47 (6): 1195-1205. 10.1023/A:1015349924116.
Pascu M, Martus P, Hohne M, Wiedenmann B, Hopf U, Schreier E, Berg T: Sustained virological response in hepatitis C virus type 1b infected patients is predicted by the number of mutations within the NS5A-ISDR: a meta-analysis focused on geographical differences. Gut. 2004, 53 (9): 1345-1351. 10.1136/gut.2003.031336.
Domingo E, Holland JJ: Mutation rates and rapid evolution of RNA viruses. The Evolutionary Biology of Viruses. 1994, Edited by Raven. New York, 161-184.
Domingo E, Holland JJ, Ahlquist P: RNA Genetics. 1988, Boca Raton, FL: CRC Press
Domingo E, Holland JJ, Biebricher C, Eigen M: Quasispecies: the concept and the word. Molecular Evolution of the Viruses. Edited by: Gibbs A, Calisher C, García-Arenal F. 1995, Cambridge: Cambridge Univ. Press, 171-180.
Domingo E, Martinez-Salas E, Sobrino F, de la Torre JC, Portela A, Ortin J, Lopez-Galindez C, Perez-Brena P, Villanueva N, Najera R, et al: The quasispecies (extremely heterogeneous) nature of viral RNA genome populations: biological relevance--a review. Gene. 1985, 40 (1): 1-8. 10.1016/0378-1119(85)90017-4.
Gebauer F, de la Torre JC, Gomes I, Mateu MG, Barahona H, Tiraboschi B, Bergmann I, de Mello PA, Domingo E: Rapid selection of genetic and antigenic variants of foot-and-mouth disease virus during persistence in cattle. J Virol. 1988, 62 (6): 2041-2049.
Holland JJ: Continuum of change in RNA virus genomes. Concepts in Viral Pathogenesis. Edited by: Notkins AL, Oldstone MBA. 1984, New York: Springer-Verlag, 137-143.
Holland J, de la Torre JC, Steinhauer DA: RNA virus populations as quasispecies. Curr Top Microbiol Immunol. 1992, 1-20. 176
Temin HM: Is HIV unique or merely different?. J Acquir Immune Defic Syndr. 1989, 2 (1): 1-9.
Wei X, Ghosh SK, Taylor ME, Johnson VA, Emini EA, Deutsch P, Lifson JD, Bonhoeffer S, Nowak MA, Hahn BH, et al: Viral dynamics in human immunodeficiency virus type 1 infection. Nature. 1995, 373 (6510): 117-122. 10.1038/373117a0.
Jardim AC, Yamasaki LH, de Queiroz AT, Bittar C, Pinho JR, Carareto CM, Rahal P, Mello IM: Quasispecies of hepatitis C virus genotype 1 and treatment outcome with Peginterferon and Ribavirin. Infect Genet Evol. 2009, 9 (4): 689-698. 10.1016/j.meegid.2008.11.001.
Craig S, Thu HM, Lowry K, Wang XF, Holmes EC, Aaskov J: Diverse dengue type 2 virus populations contain recombinant and both parental viruses in a single mosquito host. J Virol. 2003, 77 (7): 4463-4467. 10.1128/JVI.77.7.4463-4467.2003.
Aaskov J, Buzacott K, Thu HM, Lowry K, Holmes EC: Long-term transmission of defective RNA viruses in humans and Aedes mosquitoes. Science. 2006, 311 (5758): 236-238. 10.1126/science.1115030.
Sugiyama K, Suzuki K, Nakazawa T, Funami K, Hishiki T, Ogawa K, Saito S, Shimotohno KW, Suzuki T, Shimizu Y, et al: Genetic analysis of hepatitis C virus with defective genome and its infectivity in vitro. J Virol. 2009, 83 (13): 6922-6928. 10.1128/JVI.02674-08.
Tong X, Malcolm BA: Trans-complementation of HCV replication by non-structural protein 5A. Virus Res. 2006, 115 (2): 122-130. 10.1016/j.virusres.2005.07.012.
Appel N, Herian U, Bartenschlager R: Efficient rescue of hepatitis C virus RNA replication by trans-complementation with nonstructural protein 5A. J Virol. 2005, 79 (2): 896-909. 10.1128/JVI.79.2.896-909.2005.
BioMol - Electropherogram quality analysis. [http://www.biomol.unb.br/phph/]
Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8 (3): 186-194.
Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8 (3): 175-185.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999, 41: 95-98.
Marucci EA, Zafalon GF, Jardim AC, Yamasaki LH, Bittar C, Rahal P, Machado JM: Routine libraries for pattern recognition in quasispecies. Genet Mol Res. 2008, 7 (3): 970-981. 10.4238/vol7-3X-Meeting013.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) Software Version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.
Korber B: HIV Signature and Sequence Variation Analysis. Computational Analysis of HIV Molecular Sequences. 2000, Dordrecht, Netherlands: Kluwer Academic Publishers
HCV sequence database. [http://www.hiv.lanl.gov]
Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony (*and Others Methods). 2002, Sunderland, Massachusetts
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.
McCormack GP, Clewley JP: The application of molecular phylogenetics to the analysis of viral genome diversity and evolution. Rev Med Virol. 2002, 12 (4): 221-238. 10.1002/rmv.355.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/36/prepub
This work was financially supported by FAPESP and CAPES.
The authors declare that they have no competing interests.
CB, ACGJ, IMVGCM, PR: contributed to the study design, carried out the molecular biology experiments, sequence alignment, phylogenetic analysis and manuscript planning; LHTY: carried out the molecular biology experiments; ATLQ: performed statistical analysis; CMAC: contributed with significant evolutionary knowledge on data analysis; JRRP: contributed to the study design. All authors read and approved the final manuscript.