Insertion and recombination events at hypervariable region 1 over 9.6 years of hepatitis C virus chronic infection

Antigenic cooperation in viral populations: Transformation of functions of intra-host viral variants

In this paper, we study intra-host viral adaptation by antigenic cooperation - a mechanism of immune escape that serves as an alternative to the standard mechanism of escape by continuous genomic diversification and allows to explain a number of experimental observations associated with the establishment of chronic infections by highly mutable viruses. Within this mechanism, the topology of a cross-immunoreactivity network forces intra-host viral variants to specialize for complementary roles and adapt to the host's immune response as a quasi-social ecosystem. Here we study dynamical changes in immune adaptation caused by evolutionary and epidemiological events. First, we show that the emergence of a viral variant with altered antigenic features may result in a rapid re-arrangement of the viral ecosystem and a change in the roles played by existing viral variants. In particular, it may push the population under immune escape by genomic diversification towards the stable state of adaptation by antigenic cooperation. Next, we study the effect of a viral transmission between two chronically infected hosts, which results in the merging of two intra-host viral populations in the state of stable immune-adapted equilibrium. In this case, we also describe how the newly formed viral population adapts to the host's environment by changing the functions of its members. The results are obtained analytically for minimal cross-immunoreactivity networks and numerically for larger populations.

Reverse epitope mapping of the E2 glycoprotein in antibody associated hepatitis C virus

The humoral immune system responds to chronic hepatitis C virus (HCV) infection by producing neutralising antibodies (nAb). In this study we generated three HCV pseudoparticles in which E1E2 glycoprotein sequence was targeted by the host humoral immune system. We used patient derived virus free Fabs (VF-Fabs) obtained from HCV genotype 1a (n = 3), genotype 1b (n = 7) and genotype 3a (n = 1) for neutralisation of HCVpp produced in this study both individually and in combination. Based on the available anti-HCV monoclonal nAb mapping information we selected amino acid region 384-619 for conformational epitope mapping. Amongst our notable findings, we observed significant reduction in HCVpp infectivity (p<0.05) when challenged with a combination of inter genotype and subtype VF-Fabs. We also identified five binding motifs targeted by patient derived VF-Fab upon peptide mapping, of which two shared the residues with previously reported epitopes. One epitope lies within an immunodominant HVR1 and two were novel. In summary, we used a reverse epitope mapping strategy to identify preferred epitopes by the host humoral immune system. Additionally, we have combined different VF-Fabs to further reduce the HCVpp infectivity. Our data indicates that combining the antigen specificity of antibodies may be a useful strategy to reduce (in-vitro) infectivity.

Humoral immune system targets clonotypic antibody-associated hepatitis C virus

Hypervariable region 1 (HVR1) is one of the potential neutralization domains in the E2 glycoprotein of hepatitis C virus (HCV). Point mutations of the HVR1 can lead to humoral immune escape in HCV-infected patients. In this study, we segregated the chronically infected viraemic sera from HCV-infected patients into populations of antibody-free virus and antibody-associated virus (AAV) and mapped potential epitopes within the E1E2 gene junction of AAV sequences (residues 364-430). Furthermore, we generated HCV pseudoparticles (HCVpp) derived from AAV sequences to assess their infectivity. We studied the neutralization potential of virus-free Fab obtained from antibody-virus complexes, in the HCVpp system. We observed selective targeting of clonotypic HCV variants from the quasispecies pool. Moreover, we identified potential neutralizing epitopes within the HVR1 and an additional epitope that overlapped with a broadly neutralizing AP33 epitope (amino acid 412-423 in E2). We observed a marked difference in the infectivity of HCVpp generated using E1E2 sequences isolated from AAV. We document reduction in the infectivity of HCVpp-H77 and HCVpp derived from AAV sequences when challenged with virus-free Fab. Our results provide novel insights into the complexities of engagement between HCV and the humoral immune system.

Synonymous Co-Variation across the E1/E2 Gene Junction of Hepatitis C Virus Defines Virion Fitness

Hepatitis C virus is a positive-sense single-stranded RNA virus. The gene junction partitioning the viral glycoproteins E1 and E2 displays concurrent sequence evolution with the 3'-end of E1 highly conserved and the 5'-end of E2 highly heterogeneous. This gene junction is also believed to contain structured RNA elements, with a growing body of evidence suggesting that such structures can act as an additional level of viral replication and transcriptional control. We have previously used ultradeep pyrosequencing to analyze an amplicon library spanning the E1/E2 gene junction from a treatment naïve patient where samples were collected over 10 years of chronic HCV infection. During this timeframe maintenance of an in-frame insertion, recombination and humoral immune targeting of discrete virus sub-populations was reported. In the current study, we present evidence of epistatic evolution across the E1/E2 gene junction and observe the development of co-varying networks of codons set against a background of a complex virome with periodic shifts in population dominance. Overtime, the number of codons actively mutating decreases for all virus groupings. We identify strong synonymous co-variation between codon sites in a group of sequences harbouring a 3 bp in-frame insertion and propose that synonymous mutation acts to stabilize the RNA structural backbone.

Synonymous and nonsynonymous distances help untangle convergent evolution and recombination

When estimating a phylogeny from a multiple sequence alignment, researchers often assume the absence of recombination. However, if recombination is present, then tree estimation and all downstream analyses will be impacted, because different segments of the sequence alignment support different phylogenies. Similarly, convergent selective pressures at the molecular level can also lead to phylogenetic tree incongruence across the sequence alignment. Current methods for detection of phylogenetic incongruence are not equipped to distinguish between these two different mechanisms and assume that the incongruence is a result of recombination or other horizontal transfer of genetic information. We propose a new recombination detection method that can make this distinction, based on synonymous codon substitution distances. Although some power is lost by discarding the information contained in the nonsynonymous substitutions, our new method has lower false positive probabilities than the comparable recombination detection method when the phylogenetic incongruence signal is due to convergent evolution. We apply our method to three empirical examples, where we analyze: (1) sequences from a transmission network of the human immunodeficiency virus, (2) tlpB gene sequences from a geographically diverse set of 38 Helicobacter pylori strains, and (3) hepatitis C virus sequences sampled longitudinally from one patient.

A single amino acid change in the hypervariable region 1 of hepatitis C virus genotype 4a aids humoral immune escape

Longitudinal analysis of chronic hepatitis C virus (HCV) infection has shown that the virus has several adaptive strategies that maintain persistence and infectivity over time. We examined four serum samples from the same chronically infected HCV genotype 4a patient for the presence of IgG antibody-associated virus. RNA was isolated from antibody-associated and antibody-free virions. Subsequent to sequence analysis, 27 aa hypervariable region 1 (HVR1) peptides were used to test the humoral immune escape. We demonstrated that differential peptide binding of Fab was associated with a single amino acid change. We provide direct evidence of natural humoral immune escape by HCV within HVR1.

Network Analysis of the Chronic Hepatitis C Virome Defines Hypervariable Region 1 Evolutionary Phenotypes in the Context of Humoral Immune Responses

Hypervariable region 1 (HVR1) of hepatitis C virus (HCV) comprises the first 27 N-terminal amino acid residues of E2. It is classically seen as the most heterogeneous region of the HCV genome. In this study, we assessed HVR1 evolution by using ultradeep pyrosequencing for a cohort of treatment-naive, chronically infected patients over a short, 16-week period. Organization of the sequence set into connected components that represented single nucleotide substitution events revealed a network dominated by highly connected, centrally positioned master sequences. HVR1 phenotypes were observed to be under strong purifying (stationary) and strong positive (antigenic drift) selection pressures, which were coincident with advancing patient age and cirrhosis of the liver. It followed that stationary viromes were dominated by a single HVR1 variant surrounded by minor variants comprised from conservative single amino acid substitution events. We present evidence to suggest that neutralization antibody efficacy was diminished for stationary-virome HVR1 variants. Our results identify the HVR1 network structure during chronic infection as the preferential dominance of a single variant within a narrow sequence space.

IMPORTANCE

HCV infection is often asymptomatic, and chronic infection is generally well established in advance of initial diagnosis and subsequent treatment. HVR1 can undergo rapid sequence evolution during acute infection, and the variant pool is typically seen to diverge away from ancestral sequences as infection progresses from the acute to the chronic phase. In this report, we describe HVR1 viromes in chronically infected patients that are defined by a dominant epitope located centrally within a narrow variant pool. Our findings suggest that weakened humoral immune activity, as a consequence of persistent chronic infection, allows for the acquisition and maintenance of host-specific adaptive mutations at HVR1 that reflect virus fitness.

Phylogenetic analysis of a circulating hepatitis C virus recombinant strain 1b/1a isolated in a French hospital centre

Genetic recombination is now a well-established feature of the hepatitis C virus (HCV) variability and evolution, with the recent identification of circulating recombinant forms. In Amiens University Hospital Centre (France), a discrepancy of genotyping results was observed for 9 samples, between their 5' untranslated region assigned to genotype 1b and their NS5B region assigned to genotype 1a, suggesting the existence of a recombinant strain. In the present study, clinical and phylogenetic analyses of these isolates were conducted and a putative relationship with previously identified HCV 1b/1a recombinants was investigated. The results revealed that all 9 strains displayed a breakpoint within the beginning of the core protein, were closely related between each other and with the H23 strain identified in Uruguay (Moreno et al., 2009). Then, the clinical characteristics of the 9 unlinked individuals infected with this 1b/1a genotype were analysed. This is the first report on the circulation, in a French population, of a HCV recombinant strain 1b/1a. The identification of this genotype in other patients and in other geographical zones would allow to further investigate its prevalence in the population and to better understand its molecular epidemiology.

Antigenic cooperation among intrahost HCV variants organized into a complex network of cross-immunoreactivity

Hepatitis C virus (HCV) has the propensity to cause chronic infection. Continuous immune escape has been proposed as a mechanism of intrahost viral evolution contributing to HCV persistence. Although the pronounced genetic diversity of intrahost HCV populations supports this hypothesis, recent observations of long-term persistence of individual HCV variants, negative selection increase, and complex dynamics of viral subpopulations during infection as well as broad cross-immunoreactivity (CR) among variants are inconsistent with the immune-escape hypothesis. Here, we present a mathematical model of intrahost viral population dynamics under the condition of a complex CR network (CRN) of viral variants and examine the contribution of CR to establishing persistent HCV infection. The model suggests a mechanism of viral adaptation by antigenic cooperation (AC), with immune responses against one variant protecting other variants. AC reduces the capacity of the host's immune system to neutralize certain viral variants. CRN structure determines specific roles for each viral variant in host adaptation, with variants eliciting broad-CR antibodies facilitating persistence of other variants immunoreacting with these antibodies. The proposed mechanism is supported by empirical observations of intrahost HCV evolution. Interference with AC is a potential strategy for interruption and prevention of chronic HCV infection.

Hepatitis C virus genetic variability and evolution

Hepatitis C virus (HCV) has infected over 170 million people worldwide and creates a huge disease burden due to chronic, progressive liver disease. HCV is a single-stranded, positive sense, RNA virus, member of the Flaviviridae family. The high error rate of RNA-dependent RNA polymerase and the pressure exerted by the host immune system, has driven the evolution of HCV into 7 different genotypes and more than 67 subtypes. HCV evolves by means of different mechanisms of genetic variation. On the one hand, its high mutation rates generate the production of a large number of different but closely related viral variants during infection, usually referred to as a quasispecies. The great quasispecies variability of HCV has also therapeutic implications since the continuous generation and selection of resistant or fitter variants within the quasispecies spectrum might allow viruses to escape control by antiviral drugs. On the other hand HCV exploits recombination to ensure its survival. This enormous viral diversity together with some host factors has made it difficult to control viral dispersal. Current treatment options involve pegylated interferon-α and ribavirin as dual therapy or in combination with a direct-acting antiviral drug, depending on the country. Despite all the efforts put into antiviral therapy studies, eradication of the virus or the development of a preventive vaccine has been unsuccessful so far. This review focuses on current available data reported to date on the genetic mechanisms driving the molecular evolution of HCV populations and its relation with the antiviral therapies designed to control HCV infection.

Intensive temporal mapping of hepatitis C hypervariable region 1 quasispecies provides novel insights into hepatitis C virus evolution in chronic infection

Hepatitis C virus (HCV) is an RNA virus which exists as swarms of closely related viruses known as quasispecies (QS). A number of studies have demonstrated associations between QS hypervariable region 1 (HVR1) characteristics (diversity and complexity) and treatment success. We investigated HCV QS change in chronic infection over intervals of 2-4 weeks in 23 chronically infected individuals to describe the natural history of virus evolution and establish whether HCV QS characteristics could be used to individualize treatment regimens at a molecular level. HVR1 QS diversity, complexity and divergence continue to change in an unpredictable fashion in chronic infection even where there is little phylogenetic change, which is likely to preclude the use of these features in treatment individualization. Our phylogenetic analysis identified no change in the HVR1 QS in 12 subjects, minor change in four subjects and we describe a time-ordered phylogeny for the first time over a period as short as 16 weeks in seven subjects. We identified the existence of multiple subpopulation infections using partitioned analysis of QS and illustrated how subpopulations were sequentially replaced in a number of subjects. We illustrated marked variation in the nucleotide substitution per codon position between patients with sequence change and those without change in the phylogenetic tree. Analysis of codon-specific selection pressures identified a number of codons under purifying selection, suggesting that these code for structurally conserved amino acids. We also identified sections of the HVR1 under positive selection with marked sequence heterogeneity, suggesting that these may be potential epitope-binding sites.

Computational framework for next-generation sequencing of heterogeneous viral populations using combinatorial pooling

MOTIVATION

Next-generation sequencing (NGS) allows for analyzing a large number of viral sequences from infected patients, providing an opportunity to implement large-scale molecular surveillance of viral diseases. However, despite improvements in technology, traditional protocols for NGS of large numbers of samples are still highly cost and labor intensive. One of the possible cost-effective alternatives is combinatorial pooling. Although a number of pooling strategies for consensus sequencing of DNA samples and detection of SNPs have been proposed, these strategies cannot be applied to sequencing of highly heterogeneous viral populations.

RESULTS

We developed a cost-effective and reliable protocol for sequencing of viral samples, that combines NGS using barcoding and combinatorial pooling and a computational framework including algorithms for optimal virus-specific pools design and deconvolution of individual samples from sequenced pools. Evaluation of the framework on experimental and simulated data for hepatitis C virus showed that it substantially reduces the sequencing costs and allows deconvolution of viral populations with a high accuracy.

AVAILABILITY AND IMPLEMENTATION

The source code and experimental data sets are available at http://alan.cs.gsu.edu/NGS/?q=content/pooling.

Analysis of the evolution and structure of a complex intrahost viral population in chronic hepatitis C virus mapped by ultradeep pyrosequencing

Hepatitis C virus (HCV) causes chronic infection in up to 50% to 80% of infected individuals. Hypervariable region 1 (HVR1) variability is frequently studied to gain an insight into the mechanisms of HCV adaptation during chronic infection, but the changes to and persistence of HCV subpopulations during intrahost evolution are poorly understood. In this study, we used ultradeep pyrosequencing (UDPS) to map the viral heterogeneity of a single patient over 9.6 years of chronic HCV genotype 4a infection. Informed error correction of the raw UDPS data was performed using a temporally matched clonal data set. The resultant data set reported the detection of low-frequency recombinants throughout the study period, implying that recombination is an active mechanism through which HCV can explore novel sequence space. The data indicate that polyvirus infection of hepatocytes has occurred but that the fitness quotients of recombinant daughter virions are too low for the daughter virions to compete against the parental genomes. The subpopulations of parental genomes contributing to the recombination events highlighted a dynamic virome where subpopulations of variants are in competition. In addition, we provide direct evidence that demonstrates the growth of subdominant populations to dominance in the absence of a detectable humoral response.

IMPORTANCE

Analysis of ultradeep pyrosequencing data sets derived from virus amplicons frequently relies on software tools that are not optimized for amplicon analysis, assume random incorporation of sequencing errors, and are focused on achieving higher specificity at the expense of sensitivity. Such analysis is further complicated by the presence of hypervariable regions. In this study, we made use of a temporally matched reference sequence data set to inform error correction algorithms. Using this methodology, we were able to (i) detect multiple instances of hepatitis C virus intrasubtype recombination at the E1/E2 junction (a phenomenon rarely reported in the literature) and (ii) interrogate the longitudinal quasispecies complexity of the virome. Parallel to the UDPS, isolation of IgG-bound virions was found to coincide with the collapse of specific viral subpopulations.

Role of macrophages and monocytes in hepatitis C virus infections

A number of studies conducted over many years have shown that hepatitis C virus (HCV) can infect a variety of cell types. In vivo infection of monocytes, macrophages, and dendritic cells by HCV has been frequently shown by a number of researchers. These studies have demonstrated replication of HCV by detecting the presence of both negative genomic strands and a variety of non-structural HCV proteins in infected cells. In addition, analyses of genome sequences have also shown that different cell types can harbor different HCV variants. Investigators have also done preliminary studies of which cellular genes are affected by HCV infection, but there have not yet been a sufficient number of these studies to understand the effects of infection on these cells. Analyses of in vitro HCV replication have shown that monocytes, macrophages and dendritic cells can be infected by HCV from patient sera or plasma. These studies suggest that entry and cellular locations may vary between different cell types. Some studies suggest that macrophages may preferentially allow HCV genotype 1 to replicate, but macrophages do not appear to select particular hypervariable regions. Overall, these studies agree with a model where monocytes and macrophages act as an amplification system, in which these cells are infected and show few cytopathic effects, but continuously produce HCV. This allows them to produce virus over an extended time and allows its spread to other cell types.

Dynamic changes in viral population structure and compartmentalization during chronic hepatitis C virus infection in children

Classic phylogenetic and modern population-based clustering methods were used to analyze hepatitis C virus (HCV) evolution in plasma and to assess viral compartmentalization within peripheral blood mononuclear cells (PBMCs) in 6 children during 3.2-9.6yr of follow-up. Population structure analysis of cloned amplicons encompassing hypervariable region 1 led to the distinction of two evolutionary patterns, one highly divergent and another one genetically homogeneous. Viral adaptability was reflected by co-evolution of viral communities switching rapidly from one to another in the context of divergence and stability associated with highly homogeneous communities which were replaced by new ones after long periods. Additionally, viral compartmentalization of HCV in PBMCs was statistically demonstrated, suggesting their role as a pool of genetic variability. Our results support the idea of a community-based structure of HCV viral populations during chronic infection and highlight a role of the PBMC compartment in the persistence of such structure.

Productive homologous and non-homologous recombination of hepatitis C virus in cell culture

Genetic recombination is an important mechanism for increasing diversity of RNA viruses, and constitutes a viral escape mechanism to host immune responses and to treatment with antiviral compounds. Although rare, epidemiologically important hepatitis C virus (HCV) recombinants have been reported. In addition, recombination is an important regulatory mechanism of cytopathogenicity for the related pestiviruses. Here we describe recombination of HCV RNA in cell culture leading to production of infectious virus. Initially, hepatoma cells were co-transfected with a replicating JFH1ΔE1E2 genome (genotype 2a) lacking functional envelope genes and strain J6 (2a), which has functional envelope genes but does not replicate in culture. After an initial decrease in the number of HCV positive cells, infection spread after 13–36 days. Sequencing of recovered viruses revealed non-homologous recombinants with J6 sequence from the 5′ end to the NS2–NS3 region followed by JFH1 sequence from Core to the 3′ end. These recombinants carried duplicated sequence of up to 2400 nucleotides. HCV replication was not required for recombination, as recombinants were observed in most experiments even when two replication incompetent genomes were co-transfected. Reverse genetic studies verified the viability of representative recombinants. After serial passage, subsequent recombination events reducing or eliminating the duplicated region were observed for some but not all recombinants. Furthermore, we found that inter-genotypic recombination could occur, but at a lower frequency than intra-genotypic recombination. Productive recombination of attenuated HCV genomes depended on expression of all HCV proteins and tolerated duplicated sequence. In general, no strong site specificity was observed. Non-homologous recombination was observed in most cases, while few homologous events were identified. A better understanding of HCV recombination could help identification of natural recombinants and thereby lead to improved therapy. Our findings suggest mechanisms for occurrence of recombinants observed in patients.

Genetic recombination is the alternative joining of nucleic acids leading to novel combinations of genetic information. While DNA recombination in cells is of importance for evolution and adaptive immunity, RNA recombination often has only transient effects. However, RNA viruses are rapidly evolving and recombination can be an important evolutionary step in addition to mutations introduced by the viral polymerase. Recombination can allow escape from the host immune system and from antiviral treatment, and recombination of live attenuated viral vaccines has led to re-emergence of disease. Hepatitis C virus (HCV) is an important human pathogen that chronically infects more than 130 million worldwide and leads to serious liver disease. For HCV, naturally occurring recombinants are rare but clinically important. HCV recombination constitutes a challenge to antiviral treatment and can potentially provide an escape mechanism for the virus. In this study, we established an assay for HCV RNA recombination and characterized the emerging homologous and non-homologous recombinant viruses. Interestingly, recombination did not depend on viral replication, occurred most efficiently between isolates of the same genotype and did not occur with strong site-specificity. Better diagnosis of clinically important recombinants and an increased knowledge on viral recombination could strengthen antiviral and vaccine development.