HCV selection and HVR1 evolution in a chimpanzee chronically infected with HCV-1 over 12 years

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.

Local Immunodeficiency: Combining Cross-Immunoreactivity Networks

This article continues the analysis of the recently observed phenomenon of local immunodeficiency (LI), which arises as a result of antigenic cooperation among intrahost viruses organized into a network of cross-immunoreactivity (CR). We study here what happens as the result of combining (connecting) the simplest CR networks, which have a stable state of LI. It turned out that many possibilities occur, particularly resulting in a change of roles of some viruses in the CR network. Our results also give some indications about a boundary of the set of CR networks with stable state of LI in the entire collection of all possible CR networks. Such borderline CR networks are characterized by only a marginally stable (neutral rather than stable) state of the LI, or by the existence of such subnetworks in a CR network that evolve independently of each other (although being connected).

Dynamic changes of HCV genomes under selective pressure from DAAs therapy in relapsed patients

Currently, direct-acting antiviral drugs (DAAs) are widely used as therapeutic methods for hepatitis C virus (HCV)-positive patients, however, patients may experience treatment failure, and the dynamic changes of HCV genomes in these patients are unknown. In this study, three real-world DAAs cohorts were enrolled to observe clinical efficacy. In addition, serum samples from treatment failure patients at baseline and relapse were used to analyze changes of the HCV genomes at near full-length genome level, including resistance-associated variants (RAVs), viral quasispecies diversity and selection analysis. Next-generation sequencing was used as the detection method. The overall sustained virological response at 12 w after the end of treatment was achieved in 91.9% (57/62) of HCV patients, and 3 paired samples obtained from relapsed patients. All the 3 patients harbored baseline NS5A RAVs, the frequency of NS5A RAVs increased in 2 patients and a new NS5A RAV emerged in 1 patient at relapse, and almost all the viral strains existed with NS5A RAVs at relapse. The results of the viral quasispecies diversity analysis revealed that viral quasispecies diversity decreased at relapse compared to baseline, and the results of selection analysis indicated that the virus population experienced a bottleneck phenomenon, recent selective sweep and population expansion or was under purification selection after DAAs treatment. This study indicated that the clinical efficacy was excellent in real-world DAAs cohorts, and the viral strains existed at relapse were selective by DAAs therapy.

Quantitative differences between intra-host HCV populations from persons with recently established and persistent infections

Detection of incident hepatitis C virus (HCV) infections is crucial for identification of outbreaks and development of public health interventions. However, there is no single diagnostic assay for distinguishing recent and persistent HCV infections. HCV exists in each infected host as a heterogeneous population of genomic variants, whose evolutionary dynamics remain incompletely understood. Genetic analysis of such viral populations can be applied to the detection of incident HCV infections and used to understand intra-host viral evolution. We studied intra-host HCV populations sampled using next-generation sequencing from 98 recently and 256 persistently infected individuals. Genetic structure of the populations was evaluated using 245,878 viral sequences from these individuals and a set of selected features measuring their diversity, topological structure, complexity, strength of selection, epistasis, evolutionary dynamics, and physico-chemical properties. Distributions of the viral population features differ significantly between recent and persistent infections. A general increase in viral genetic diversity from recent to persistent infections is frequently accompanied by decline in genomic complexity and increase in structuredness of the HCV population, likely reflecting a high level of intra-host adaptation at later stages of infection. Using these findings, we developed a machine learning classifier for the infection staging, which yielded a detection accuracy of 95.22 per cent, thus providing a higher accuracy than other genomic-based models. The detection of a strong association between several HCV genetic factors and stages of infection suggests that intra-host HCV population develops in a complex but regular and predictable manner in the course of infection. The proposed models may serve as a foundation of cyber-molecular assays for staging infection, which could potentially complement and/or substitute standard laboratory assays.

Analysis of heterogeneous genomic samples using image normalization and machine learning

BACKGROUND

Analysis of heterogeneous populations such as viral quasispecies is one of the most challenging bioinformatics problems. Although machine learning models are becoming to be widely employed for analysis of sequence data from such populations, their straightforward application is impeded by multiple challenges associated with technological limitations and biases, difficulty of selection of relevant features and need to compare genomic datasets of different sizes and structures.

RESULTS

We propose a novel preprocessing approach to transform irregular genomic data into normalized image data. Such representation allows to restate the problems of classification and comparison of heterogeneous populations as image classification problems which can be solved using variety of available machine learning tools. We then apply the proposed approach to two important problems in molecular epidemiology: inference of viral infection stage and detection of viral transmission clusters using next-generation sequencing data. The infection staging method has been applied to HCV HVR1 samples collected from 108 recently and 257 chronically infected individuals. The SVM-based image classification approach achieved more than 95% accuracy for both recently and chronically HCV-infected individuals. Clustering has been performed on the data collected from 33 epidemiologically curated outbreaks, yielding more than 97% accuracy.

CONCLUSIONS

Sequence image normalization method allows for a robust conversion of genomic data into numerical data and overcomes several issues associated with employing machine learning methods to viral populations. Image data also help in the visualization of genomic data. Experimental results demonstrate that the proposed method can be successfully applied to different problems in molecular epidemiology and surveillance of viral diseases. Simple binary classifiers and clustering techniques applied to the image data are equally or more accurate than other models.

A Hepatitis C Virus DNA Vaccine Encoding a Secreted, Oligomerized Form of Envelope Proteins Is Highly Immunogenic and Elicits Neutralizing Antibodies in Vaccinated Mice

Hepatitis C virus (HCV) persistently infects approximately 71 million people globally. To prevent infection a vaccine which elicits neutralizing antibodies against the virus envelope proteins (E1/E2) which are required for entry into host cells is desirable. DNA vaccines are cost-effective to manufacture globally and despite recent landmark studies highlighting the therapeutic efficacy of DNA vaccines in humans against cervical cancer, DNA vaccines encoding E1/E2 developed thus far are poorly immunogenic. We now report a novel and highly immunogenic DNA vaccination strategy that incorporates secreted E1 and E2 (sE1 and sE2) into oligomers by fusion with the oligomerization domain of the C4b-binding protein, IMX313P. The FDA approved plasmid, pVax, was used to encode sE1, sE2, or sE1E2 with or without IMX313P, and intradermal prime-boost vaccination studies in BALB/c mice showed that vaccines encoding IMX313P were the most effective in eliciting humoral and cell-mediated immunity against the envelope proteins. Further boosting with recombinant E1E2 proteins but not DNA nor virus-like particles (VLPs) expressing E1E2 increased the immunogenicity of the DNA prime-boost regimen. Nevertheless, the antibodies generated by the homologous DNA prime-boost vaccinations more effectively inhibited the binding of VLPs to target cells and neutralized transduction with HCV pseudoparticles (HCVpp) derived from different genotypes including genotypes 1, 2, 3, 4, 5, and 6. This report provides the first evidence that IMX313P can be used as an adjuvant for E1/E2-based DNA vaccines and represents a translatable approach for the development of a HCV DNA vaccine.

Tree shrew, a potential animal model for hepatitis C, supports the infection and replication of HCV in vitro and in vivo

The tree shrew (Tupaia belangeri chinensis), a small animal widely distributed in Southeast Asia and southwest China, has the potential to be developed as an animal model for hepatitis C. To determine the susceptibility of the tree shrew to hepatitis C virus (HCV) infection in vitro and in vivo, a well-established HCV, produced from the J6/JFH1-Huh7.5.1 culture system, was used to infect cultured primary tupaia hepatocytes (PTHs) and tree shrews. The in vitro results showed that HCV genomic RNA and HCV-specific nonstructural protein 5A (NS5A) could be detected in the PTH cell culture from days 3-15 post-infection, although the viral load was lower than that observed in Huh7.5.1 cell culture. The occurrence of five sense mutations [S391A, G397A, L402F and M405T in the hypervariable region 1 (HVR1) of envelope glycoprotein 2 and I2750M in NS5B] suggested that HCV undergoes genetic evolution during culture. Fourteen of the 30 experimental tree shrews (46.7 %) were found to be infected, although the HCV viremia was intermittent in vivo. A positive test for HCV RNA in liver tissue provided stronger evidence for HCV infection and replication in tree shrews. The results of an immunohistochemistry assay also demonstrated the presence of four HCV-specific proteins (Core, E2, NS3/4 and NS5A) in the hepatocytes of infected tree shrews. The pathological changes observed in the liver tissue of infected tree shrews could be considered to be representative symptoms of mild hepatitis. These results revealed that the tree shrew can be used as an animal model supporting the infection and replication of HCV in vitro and in vivo.

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.

The evolutionary patterns of hepatitis C virus subtype 2a and 6a isolates linked to an outbreak in China in 2012

An HCV outbreak occurred in 2012 in China, affecting hundreds of patients. We characterized HCV subtype 2a and 6a sequences from 60 and 102 patients, respectively, and co-analyzed them with 82 local controls and 103 calibrating references. The close grouping of the patients׳ sequences contrasted sharply with the diversity of local controls. Scaled by the calibrating references, the emergence of patients׳ isolates was estimated at 2-5 years before sampling. In contrast, the controls intermingled with the calibrating references that were much older. For both subtypes, the major and minor clusters could be defined, with the closeness to indicate linked transmission.

CONCLUSION

HCV sequences from the study patients grouped into three subtype 2a and two subtype 6a clusters, in addition to three 6a solitary branches, representing descendants of eight earlier strains that were distinct and otherwise sporadic. Due to iatrogenic transmission through reusing needles, five strains were highly selected and preferentially spread.

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.

Fine mapping of murine antibody responses to immunization with a novel soluble form of hepatitis C virus envelope glycoprotein complex

The hepatitis C virus (HCV) envelope glycoprotein E1E2 complex is a candidate vaccine antigen. Previous immunization studies of E1E2 have yielded various results on its ability to induce virus-neutralizing antibodies in animal models and humans. The murine model has become a vital tool for HCV research owing to the development of humanized mice susceptible to HCV infection. In this study, we investigated the antibody responses of mice immunized with E1E2 and a novel soluble form of E1E2 (sE1E2) by a DNA prime and protein boost strategy. The results showed that sE1E2 elicited higher antibody titers and a greater breadth of reactivity than the wild-type cell-associated E1E2. However, immune sera elicited by either immunogen were only weakly neutralizing. In order to understand the contrasting results of binding and serum neutralizing activities, epitopes targeted by the polyclonal antibody responses were mapped and monoclonal antibodies (MAbs) were generated. The results showed that the majority of serum antibodies were directed to the E1 region 211 to 250 and the E2 regions 421 to 469, 512 to 539, 568 to 609, and 638 to 651, instead of the well-known immunodominant E2 hypervariable region 1 (HVR1). Unexpectedly, in MAb analysis, ∼ 12% of MAbs isolated were specific to the conserved E2 antigenic site 412 to 423, and 85% of them cross-neutralized multiple HCV isolates. The epitopes recognized by these MAbs are similar but distinct from the previously reported HCV1 and AP33 broadly neutralizing epitopes. In conclusion, E1E2 can prime B cells specific to conserved neutralizing epitopes, but the levels of serum neutralizing antibodies elicited are insufficient for effective virus neutralization. The sE1E2 constructs described in this study can be a useful template for rational antigen engineering.

IMPORTANCE

Hepatitis C virus infects 2 to 3% of the world's population and is a leading cause of liver failures and the need for liver transplantation. The virus envelope glycoprotein complex E1E2 produced by detergent extraction of cells overexpressing the protein was evaluated in a phase I clinical trial but failed to induce neutralizing antibodies in most subjects. In this study, we designed a novel form of E1E2 which is secreted from cells and is soluble and compared it to wild-type E1E2 by DNA immunization of mice. The results showed that this new E1E2 is more immunogenic than wild-type E1E2. Detailed mapping of the antibody responses revealed that antibodies to the conserved E2 antigenic site 412 to 423 were elicited but the serum concentrations were too low to neutralize the virus effectively. This soluble E1E2 provides a new reagent for studying HCV and for rational vaccine design.

Origin of hepatitis C virus genotype 3 in Africa as estimated through an evolutionary analysis of the full-length genomes of nine subtypes, including the newly sequenced 3d and 3e

We characterized the full-length genomes of nine hepatitis C virus genotype 3 (HCV-3) isolates: QC7, QC8, QC9, QC10, QC34, QC88, NE145, NE274 and 811. To the best of our knowledge, NE274 and NE145 were the first full-length genomes for confirming the provisionally assigned subtypes 3d and 3e, respectively, whereas 811 represented the first HCV-3 isolate that had its extreme 3' UTR terminus sequenced. Based on these full-length genomes, together with 42 references representing eight assigned subtypes and an unclassified variant of HCV-3, and 10 sequences of six other genotypes, a timescaled phylogenetic tree was reconstructed after an evolutionary analysis using a coalescent Bayesian procedure. The results indicated that subtypes 3a, 3d and 3e formed a subset with a common ancestor dated to ~202.89 [95% highest posterior density (HPD): 160.11, 264.6] years ago. The analysis of all of the HCV-3 sequences as a single lineage resulted in the dating of the divergence time to ~457.81 (95% HPD: 350.62, 587.53) years ago, whereas the common ancestor of all of the seven HCV genotypes dated to ~780.86 (95% HPD: 592.15, 1021.34) years ago. As subtype 3h and the unclassified variant were relatives, and represented the oldest HCV-3 lineages with origins in Africa and the Middle East, these findings may indicate the ancestral origin of HCV-3 in Africa. We speculate that the ancestral HCV-3 strains may have been brought to South Asia from Africa by land and/or across the sea to result in its indigenous circulation in that region. The spread was estimated to have occurred in the era after Vasco da Gama had completed his expeditions by sailing along the eastern coast of Africa to India. However, before this era, Arabians had practised slave trading from Africa to the Middle East and South Asia for centuries, which may have mediated the earliest spread of HCV-3.