This article has Open Peer Review reports available.
Genetic diversity of NS5A protein from hepatitis C virus genotype 3a and its relationship to therapy response
© Bittar et al; licensee BioMed Central Ltd. 2010
Received: 9 April 2009
Accepted: 23 February 2010
Published: 23 February 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.
Viral load (UI/ml)
Sustained virological responder
Sustained virological responder
Sustained virological responder
Sustained virological responder
End of treatment responder
End of treatment responder
End of treatment responder
End of treatment responder
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.
Number of mutations, genetic distance, dN, dS and dN/dS.
Regions and types of treatment response
Genetic distance ws
Genetic distance bs
0.1797* b, * c
0.0336* b, * c
0.0664* a, * c
0.4130* a, * c
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.
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 18.104.22.168  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.
This work was financially supported by FAPESP and CAPES.
- Bowen DG, Walker CM: Adaptive immune responses in acute and chronic hepatitis C virus infection. Nature. 2005, 436 (7053): 946-952. 10.1038/nature04079.View ArticlePubMedGoogle Scholar
- Alter MJ: Epidemiology of hepatitis C virus infection. World J Gastroenterol. 2007, 13 (17): 2436-2441.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Gale M, Foy EM: Evasion of intracellular host defence by hepatitis C virus. Nature. 2005, 436 (7053): 939-945. 10.1038/nature04078.View ArticlePubMedGoogle Scholar
- Eigen M, Biebricher C: Sequence space and quasispecies distribution. 1988, Boca Raton, FL: CRC Press, 3:Google Scholar
- Eigen M: Viral quasispecies. Sci Am. 1993, 269 (1): 42-49. 10.1038/scientificamerican0793-42.View ArticlePubMedGoogle Scholar
- Coffin JM: HIV population dynamics in vivo: implications for genetic variation, pathogenesis, and therapy. Science. 1995, 267 (5197): 483-489. 10.1126/science.7824947.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- Domingo E, Holland JJ, Ahlquist P: RNA Genetics. 1988, Boca Raton, FL: CRC PressGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- Holland J, de la Torre JC, Steinhauer DA: RNA virus populations as quasispecies. Curr Top Microbiol Immunol. 1992, 1-20. 176Google Scholar
- Temin HM: Is HIV unique or merely different?. J Acquir Immune Defic Syndr. 1989, 2 (1): 1-9.PubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.Google Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Korber B: HIV Signature and Sequence Variation Analysis. Computational Analysis of HIV Molecular Sequences. 2000, Dordrecht, Netherlands: Kluwer Academic PublishersGoogle Scholar
- HCV sequence database. [http://www.hiv.lanl.gov]
- Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony (*and Others Methods). 2002, Sunderland, MassachusettsGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14 (9): 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/36/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.