Epidemic dynamics of two coexisting hepatitis C virus subtypes

A novel approach to identify candidate prognostic factors for hepatitis C treatment response integrating clinical and viral genetic data

The combined therapy of pegylated interferon (IFN) plus ribavirin (RBV) has been for a long time the standard treatment for patients infected with hepatitis C virus (HCV). In the case of genotype 1, only 38%–48% of patients have a positive response to the combined treatment. In previous studies, viral genetic information has been occasionally included as a predictor. Here, we consider viral genetic variation in addition to 11 clinical and 19 viral populations and evolutionary parameters to identify candidate baseline prognostic factors that could be involved in the treatment outcome. We obtained potential prognostic models for HCV subtypes la and lb in combination as well as separately. We also found that viral genetic information is relevant for the combined treatment assessment of patients, as the potential prognostic model of joint subtypes includes 9 viral-related variables out of 11. Our proposed methodology fully characterizes viral genetic information and finds a combination of positions that modulate inter-patient variability.

Relevance of baseline viral genetic heterogeneity and host factors for treatment outcome prediction in hepatitis C virus 1b-infected patients

BACKGROUND

Only about 50% of patients chronically infected with HCV genotype 1 (HCV-1) respond to treatment with pegylated interferon-alfa and ribavirin (dual therapy), and protease inhibitors have to be administered together with these drugs increasing costs and side-effects. We aimed to develop a predictive model of treatment response based on a combination of baseline clinical and viral parameters.

METHODOLOGY

Seventy-four patients chronically infected with HCV-1b and treated with dual therapy were studied (53 retrospectively -training group-, and 21 prospectively -validation group-). Host and viral-related factors (viral load, and genetic variability in the E1-E2, core and Interferon Sensitivity Determining Region) were assessed. Multivariate discriminant analysis and decision tree analysis were used to develop predictive models on the training group, which were then validated in the validation group.

PRINCIPAL FINDINGS

A multivariate discriminant predictive model was generated including the following variables in decreasing order of significance: the number of viral variants in the E1-E2 region, an amino acid substitution pattern in the viral core region, the IL28B polymorphism, serum GGT and ALT levels, and viral load. Using this model treatment outcome was accurately predicted in the training group (AUROC = 0.9444; 96.3% specificity, 94.7% PPV, 75% sensitivity, 81% NPV), and the accuracy remained high in the validation group (AUROC = 0.8148, 88.9% specificity, 90.0% PPV, 75.0% sensitivity, 72.7% NPV). A second model was obtained by a decision tree analysis and showed a similarly high accuracy in the training group but a worse reproducibility in the validation group (AUROC = 0.9072 vs. 0.7361, respectively).

CONCLUSIONS AND SIGNIFICANCE

The baseline predictive models obtained including both host and viral variables had a high positive predictive value in our population of Spanish HCV-1b treatment naïve patients. Accurately identifying those patients that would respond to the dual therapy could help reducing implementation costs and additional side effects of new treatment regimens.

Recombination in hepatitis C virus

Hepatitis C virus (HCV) is a Flavivirus with a positive-sense, single-stranded RNA genome of about 9,600 nucleotides. It is a major cause of liver disease, infecting almost 200 million people all over the world. Similarly to most RNA viruses, HCV displays very high levels of genetic diversity which have been used to differentiate six major genotypes and about 80 subtypes. Although the different genotypes and subtypes share basic biological and pathogenic features they differ in clinical outcomes, response to treatment and epidemiology. The first HCV recombinant strain, in which different genome segments derived from parentals of different genotypes, was described in St. Petersburg (Russia) in 2002. Since then, there have been only a few more than a dozen reports including descriptions of HCV recombinants at all levels: between genotypes, between subtypes of the same genotype and even between strains of the same subtype. Here, we review the literature considering the reasons underlying the difficulties for unequivocally establishing recombination in this virus along with the analytical methods necessary to do it. Finally, we analyze the potential consequences, especially in clinical practice, of HCV recombination in light of the coming new therapeutic approaches against this virus.

Relating the liver damage with hepatitis C virus polymorphism in core region and human variables in HIV-1-coinfected patients

Hepatitis C virus (HCV) infection is the most important cause of chronic hepatitis, cirrhosis and end-stage liver disease leading to liver transplantation worldwide. Chronic infection by HCV causes liver fibrosis, which is accelerated by unknown mechanisms in patients with human immunodeficiency virus-1 (HIV-1) coinfection. Although the genetic variability of both HCV and HIV has been extensively studied in the context of monoinfections, more limited data is available regarding HCV-HIV coinfection. HCV disease progression among HIV coinfected patients may be influenced not only by demographic, epidemiological and clinical background variables, but also by genetic differences in infecting viruses. To explore this issue, we carried out a study in coinfected patients trying to associate the degree of liver damage to several demographic, clinical, and epidemiological characteristics of the patients, and also to the genetic variability of HCV between patients. For this purpose, we have applied different statistical techniques including the statistical generalized linear model (GLM) framework. The stage of fibrosis was indirectly measured by noninvasive means using the indexes Forns, APRI and FIB-4. HCV genetic variability between patients was estimated by sequencing the core region and by reconstructions of consensus maximum parsimony phylogenetic trees with 50% and 75% bootstrap majority rules. The results showed a direct correlation of the fibrosis biomarkers with the AST/ALT ratio, MoftIDU and with 3a HCV genotype clades, among others.

Baseline prediction of combination therapy outcome in hepatitis C virus 1b infected patients by discriminant analysis using viral and host factors

Background

Current treatment of chronic hepatitis C virus (HCV) infection has limited efficacy −especially among genotype 1 infected patients−, is costly, and involves severe side effects. Thus, predicting non-response is of major interest for both patient wellbeing and health care expense. At present, treatment cannot be individualized on the basis of any baseline predictor of response. We aimed to identify pre-treatment clinical and virological parameters associated with treatment failure, as well as to assess whether therapy outcome could be predicted at baseline.

Methodology

Forty-three HCV subtype 1b (HCV-1b) chronically infected patients treated with pegylated-interferon alpha plus ribavirin were retrospectively studied (21 responders and 22 non-responders). Host (gender, age, weight, transaminase levels, fibrosis stage, and source of infection) and viral-related factors (viral load, and genetic variability in the E1–E2 and Core regions) were assessed. Logistic regression and discriminant analyses were used to develop predictive models. A “leave-one-out” cross-validation method was used to assess the reliability of the discriminant models.

Principal Findings

Lower alanine transaminase levels (ALT, p = 0.009), a higher number of quasispecies variants in the E1–E2 region (number of haplotypes, nHap_E1–E2) (p = 0.003), and the absence of both amino acid arginine at position 70 and leucine at position 91 in the Core region (p = 0.039) were significantly associated with treatment failure. Therapy outcome was most accurately predicted by discriminant analysis (90.5% sensitivity and 95.5% specificity, 85.7% sensitivity and 81.8% specificity after cross-validation); the most significant variables included in the predictive model were the Core amino acid pattern, the nHap_E1–E2, and gamma-glutamyl transferase and ALT levels.

Conclusions and Significance

Discriminant analysis has been shown as a useful tool to predict treatment outcome using baseline HCV genetic variability and host characteristics. The discriminant models obtained in this study led to accurate predictions in our population of Spanish HCV-1b treatment naïve patients.

Viral phylodynamics and the search for an 'effective number of infections'

Information on the dynamics of the effective population size over time can be obtained from the analysis of phylogenies, through the application of time-varying coalescent models. This approach has been used to study the dynamics of many different viruses, and has demonstrated a wide variety of patterns, which have been interpreted in the context of changes over time in the ‘effective number of infections’, a quantity proportional to the number of infected individuals. However, for infectious diseases, the rate of coalescence is driven primarily by new transmissions i.e. the incidence, and only indirectly by the number of infected individuals through sampling effects. Using commonly used epidemiological models, we show that the coalescence rate may indeed reflect the number of infected individuals during the initial phase of exponential growth when time is scaled by infectivity, but in general, a single change in time scale cannot be used to estimate the number of infected individuals. This has important implications when integrating phylogenetic data in the context of other epidemiological data.

Acute hepatitis C infection with evidence of heterosexual transmission

A 62-year-old woman acquired acute hepatitis C virus (HCV) infection after heterosexual contact with a known HCV positive former injecting drug user. There were no known sexual or other risk factors for HCV acquisition. Phylogenetic analysis confirmed the case and index were infected with identical genotype 3a strains, consistent with heterosexual transmission in the absence of specific risk factors.

Cellular immune response to infection by different genotypes of hepatitis C virus

The present study aimed to estimate the cellular immune response to infection by different HCV genotypes in some Egyptian patients with chronic liver disease. HCV-RNA and HCV genotypes were conducted using type-specific primers to amplify the core region of HCV genome. The levels of IFN-gamma and TNF-alpha, PIIINP, SOD enzyme were also determined. The results showed that genotype 4 was the prevalent one (88.57%).Genotype 3 was found mixed with genotype 2 and 4 (11.42% and 34.28% respectively). Genotype 1 was absent from patients under study. No significant differences were observed in the level of IFN-gamma interferon, SOD enzyme and ALT enzyme between those genotypes. On the other hand the level of TNF-alpha was greatly increased in patients infected with mixed HCV genotypes 3, 4 compared with patients with genotype 4 and those infected with mixed genotypes 2, 3. Otherwise PIIINP was decreased significantly of the same group compared with other groups.

Genetic variability in the 5' UTR and NS5A regions of hepatitis C virus RNA isolated from non-responding and responding patients with chronic HCV genotype 1 infection

Sequence variation among different hepatitis C virus (HCV) isolates has adaptive significance and reflects the modes and intensities of selection mechanisms operating on the virus. In this work, we sought to investigate using classical population genetics parameters, the genetic variability of HCV genotype 1 using the 5' UTR and NS5A regions from treatment non-responding and responding groups of patients. Both regions showed low genetic variability and the 5' UTR showed neutral deviation. No differences were observed in the nonsynonymous/synonymous nucleotide substitution ratio among groups for NS5A. The analysis of molecular variance test of the 5' UTR region showed an 11.94% variation among groups. Phylogenetic analysis showed no correlation between sequence variations and therapeutic responses.

Combined therapy of interferon plus ribavirin promotes multiple adaptive solutions in hepatitis C virus

Hepatitis C virus (HCV) presents several regions involved potentially in evading antiviral treatment and host immune system. Two regions, known as PKR-BD and V3 domains, have been proposed to be involved in resistance to interferon. Additionally, hypervariable regions in the envelope E2 glycoprotein are also good candidates to participate in evasion from the immune system. In this study, we have used a cohort of 22 non-responder patients to combined therapy (interferon alpha-2a plus ribavirin) for which samples obtained just before initiation of therapy and after 6 or/and 12 months of treatment were available. A range of 25-100 clones per patient, genome region and time sample were obtained. The predominant amino acid sequences for each time sample and patient were determined. Next, the sequences of the PKR-BD and V3 domains and the hypervariable regions from different time samples were compared for each patient. The highest levels of variability were detected at the three hypervariable regions of the E2 protein and, to a lower extent, at the V3 domain of the NS5A protein. However, no clear patterns of adaptation to the host immune system or to antiviral treatment were detected. In summary, although high levels of variability are correlated to viral adaptive response, antiviral treatment does not seem to promote convergent adaptive changes. Consequently, other regions must be involved in evasion strategies likely based on a combination of multiple mechanisms, in which pools of changes along the HCV genome could confer viruses the ability to overcome strong selective pressures.

Evidence of recombination in intrapatient populations of hepatitis C virus

Hepatitis C virus (HCV) is a major cause of liver disease worldwide and a potential cause of substantial morbidity and mortality in the future. HCV is characterized by a high level of genetic heterogeneity. Although homologous recombination has been demonstrated in many members of the family Flaviviridae, to which HCV belongs, there are only a few studies reporting recombination on natural populations of HCV, suggesting that these events are rare in vivo. Furthermore, these few studies have focused on recombination between different HCV genotypes/subtypes but there are no reports on the extent of intra-genotype or intra-subtype recombination between viral strains infecting the same patient. Given the important implications of recombination for RNA virus evolution, our aim in this study has been to assess the existence and eventually the frequency of intragenic recombination on HCV. For this, we retrospectively have analyzed two regions of the HCV genome (NS5A and E1-E2) in samples from two different groups: (i) patients infected only with HCV (either treated with interferon plus ribavirin or treatment naïve), and (ii) HCV-HIV co-infected patients (with and without treatment against HIV). The complete data set comprised 17712 sequences from 136 serum samples derived from 111 patients. Recombination analyses were performed using 6 different methods implemented in the program RDP3. Recombination events were considered when detected by at least 3 of the 6 methods used and were identified in 10.7% of the amplified samples, distributed throughout all the groups described and the two genomic regions studied. The resulting recombination events were further verified by detailed phylogenetic analyses. The complete experimental procedure was applied to an artificial mixture of relatively closely viral populations and the ensuing analyses failed to reveal artifactual recombination. From these results we conclude that recombination should be considered as a potentially relevant mechanism generating genetic variation in HCV and with important implications for the treatment of this infection.

Genetic variability of hepatitis C virus before and after combined therapy of interferon plus ribavirin

We present an analysis of the selective forces acting on two hepatitis C virus genome regions previously postulated to be involved in the viral response to combined antiviral therapy. One includes the three hypervariable regions in the envelope E2 glycoprotein, and the other encompasses the PKR binding domain and the V3 domain in the NS5A region. We used a cohort of 22 non-responder patients to combined therapy (interferon alpha-2a plus ribavirin) for which samples were obtained before initiation of therapy and after 6 or/and 12 months of treatment. A range of 25-100 clones per patient, genome region and time sample were sequenced. These were used to detect general patterns of adaptation, to identify particular adaptation mechanisms and to analyze the patterns of evolutionary change in both genome regions. These analyses failed to detect a common adaptive mechanism for the lack of response to antiviral treatment in these patients. On the contrary, a wide range of situations were observed, from patients showing no positively selected sites to others with many, and with completely different topologies in the reconstructed phylogenetic trees. Altogether, these results suggest that viral strategies to evade selection pressure from the immune system and antiviral therapies do not result from a single mechanism and they are likely based on a range of different alternatives, in which several different changes, or their combination, along the HCV genome confer viruses the ability to overcome strong selective pressures.

Refined analysis of genetic variability parameters in hepatitis C virus and the ability to predict antiviral treatment response

Hepatitis C virus (HCV) infects approximately 3% of the world population. The chronicity of hepatitis C seems to depend on the level of genetic variability. We have recently (Torres-Puente et al., J Viral Hepat, 2008; 15: 188) reported genetic variability estimates from a large-scale sequence analysis of 67 patients infected with HCV subtypes 1a (23 patients) and 1b (44 patients) and related them to response, or lack of, to alpha-interferon plus ribavirin treatment.. Two HCV genome regions were analysed in samples prior to antiviral therapy, one compressing the three hypervariable regions of the E2 glycoprotein and another one including the interferon sensitive determining region and the V3 domain of the NS5A protein. Haplotype and nucleotide diversity measures showed a clear tendency to higher genetic variability levels in nonresponder than in responder patients. Here, we have refined the analysis of genetic variability (haplotype and nucleotide diversity, number of haplotypes and mutations) by considering their distribution in each of the biologically meaningful subregions mentioned above, as well as in their surrounding and intervening regions. Variability levels are very heterogeneous among the different subregions, being higher for nonresponder patients. Interestingly, significant differences were detected in the biologically relevant regions, but also in the surrounding regions, suggesting that the level of variability of the whole HCV genome, rather than exclusively that from the hypervariable regions, is the main indicator of the treatment response. Finally, the number of haplotypes and mutations seem to be better discriminators than haplotype and nucleotide diversity, especially in the NS5A region.

Composition of inflammatory infiltrate and its correlation with HBV/HCV antigen expression

AIM: To study the composition of liver inflammatory infiltrate in biopsy material from patients chronically infected with hepatotropic viruses and to evaluate the correlation of inflammatory infiltrate with hepatitis B virus (HBV) and hepatitis C virus (HCV) viral antigen expression in chronic B and C hepatitis.

METHODS: The phenotype of inflammatory cells was evaluated by the EnVision system, using a panel of monoclonal antibodies. HBV and HCV antigens were detected with the use of monoclonal anti-HBs, polyclonal anti-HBc and anti-HCV antibodies, respectively.

RESULTS: The cellular composition of liver inflammatory infiltrate was similar in the patients with B and C hepatitis: ~50%-60% of cells were T helper lymphocytes. Approximately 25% were T cytotoxic lymphocytes; B lymphocytes comprised 15% of inflammatory infiltrate; other cells, including NK, totalled 10%. Expression of HLA antigens paralleled inflammatory activity. Portal lymphadenoplasia was found more often in hepatitis C (54.5%) than in hepatitis B (30.6%). Expression of HBcAg was found more often in chronic B hepatitis of moderate or severe activity. Overall inflammatory activity in HBV-infected cases did not correlate with the intensity of HBsAg expression in hepatocytes. Inflammatory infiltrates accompanied the focal expression of HCV antigens. A direct correlation between antigen expression and inflammatory reaction in situ was noted more often in hepatitis C than B.

CONCLUSION: Irrespective of the etiology and activity of hepatitis, components of the inflammatory infiltrate in liver were similar. Overall inflammatory activity did not correlate with the expression of HBsAg and HCVAg; HBcAg expression, however, accompanied chronic hepatitis B of moderate and severe activity.

Complete genome of a European hepatitis C virus subtype 1g isolate: phylogenetic and genetic analyses

BACKGROUND

Hepatitis C virus isolates have been classified into six main genotypes and a variable number of subtypes within each genotype, mainly based on phylogenetic analysis. Analyses of the genetic relationship among genotypes and subtypes are more reliable when complete genome sequences (or at least the full coding region) are used; however, so far 31 of 80 confirmed or proposed subtypes have at least one complete genome available. Of these, 20 correspond to confirmed subtypes of epidemic interest.

RESULTS

We present and analyse the first complete genome sequence of a HCV subtype 1g isolate. Phylogenetic and genetic distance analyses reveal that HCV-1g is the most divergent subtype among the HCV-1 confirmed subtypes. Potential genomic recombination events between genotypes or subtype 1 genomes were ruled out. We demonstrate phylogenetic congruence of previously deposited partial sequences of HCV-1g with respect to our sequence.

CONCLUSION

In light of this, we propose changing the current status of its subtype-specific designation from provisional to confirmed.

Genetic variability in hepatitis C virus and its role in antiviral treatment response

Hepatitis C virus (HCV) is a major health problem worldwide, infecting an estimated 170 million people. The high genetic variability of HCV contributes to the chronicity of hepatitis C. Here, we report results from a large-scale sequence analysis of 67 patients infected with HCV genotype 1, 23 with subtype 1a and 44 with subtype 1b. Two regions of the HCV genome were analysed in samples prior to combined therapy with alpha interferon plus ribavirin, one compressing the hypervariable regions (HVR1, HVR2 and HVR3) of the E2 glycoprotein and another one including the interferon-sensitive determining region (ISDR) and the V3 domain of the NS5A protein. Genetic diversity measures showed a clear tendency to higher genetic variability levels in nonresponder patients to antiviral treatment than in responder patients, although highly disperse values were present within each response group for both subtypes. A more detailed analysis of amino acid composition revealed the presence of several subtype-specific variants in a few positions, but no discriminating positions between responder and nonresponder patients were detected. Our results also revealed that most amino acid positions were highly conserved, especially for subtype 1a. We conclude that the outcome of the antiviral treatment might depend not only on the nature of one or a few independent positions, but more likely on the combination of several positions along the HCV genome. Moreover, the own host's ability to generate an appropriate systemic response, in combination with the action of antivirals, is also likely to be essential for treatment outcome.

Effect of antiviral treatment and host susceptibility on positive selection in hepatitis C virus (HCV)

We have conducted a large sequence study of the E1-E2 and NS5A regions of the HCV, subtypes 1a and b, both in patients previously treated with interferon, and untreated patients, who later responded, or not, to a combination therapy based on interferon plus ribavirin. We have examined the role played by the number of positively selected sites on disease progression and its relationship with several variables such as patients' age, sex and their risk of acquiring the disease. We have detected three groups of patients that respond or not to combination therapy: responders of intermediate age, older non-responders and young non-responders, they possess an increasing average number of positively selected sites in the E1-E2 region, respectively. We conclude that the host's genetic factors play an important role in whether the disease is contained or becomes chronic.

Hepatitis C virus and the controversial role of the interferon sensitivity determining region in the response to interferon treatment

The degree of variability of the interferon sensitivity determining region (ISDR) in the hepatitis C virus (HCV) genome has been postulated to predict the response to interferon therapy, mainly in patients infected with subtype 1b, although this prediction has been the subject of a long controversy. This prediction has been tested by analyzing a cohort of 67 Spanish patients infected with HCV genotype 1, 23 of which were infected with subtype 1a and 44 with subtype 1b. A sample previous to therapy with alpha-interferon plus ribavirin was obtained and several clones (between 25 and 96) including the ISDR were sequenced from each patient. A significant correlation between mutations at the ISDR and response to treatment for subtype 1b patients, but not for those infected with subtype 1a, has been detected. Although the results suggest that the same relationship holds true for subtype 1a, lack of statistical power because of the small sample size of this subtype prevented firmer conclusions. However, identical ISDR sequences were found in responder and non-responder patients, suggesting that the stability of the ISDR sequence can occasionally help HCV to evade interferon therapy, although this is not a sufficient condition. More complex interactions, including the ISDR or not, are likely to exist and govern the HCV response to interferon treatment.

Contribution of insertions and deletions to the variability of hepatitis C virus populations

Little is known about the potential effects of insertions and deletions (indels) on the evolutionary dynamics of hepatitis C virus (HCV). In fact, the consequences of indels on antiviral treatment response are a field of investigation completely unexplored. Here, an extensive sequencing project was undertaken by cloning and sequencing serum samples from 25 patients infected with HCV subtype 1a and 48 patients with subtype 1b. For 23 patients, samples obtained after treatment with alpha interferon plus ribavirin were also available. Two genome fragments containing the hypervariable regions in the envelope 2 glycoprotein and the PKR-BD domain in NS5A were sequenced, yielding almost 16 000 sequences. Our results show that insertions are quite rare, but they are often present in biologically relevant domains of the HCV genome. Moreover, their frequency distributions between different time samples reflect the quasispecies dynamics of HCV populations. Deletions seem to be subject to negative selection.