Transmission of clonal hepatitis C virus genomes reveals the dominant but transitory role of CD8⁺ T cells in early viral evolution

The RNA genome of the hepatitis C virus (HCV) diversifies rapidly during the acute phase of infection, but the selective forces that drive this process remain poorly defined. Here we examined whether Darwinian selection pressure imposed by CD8(+) T cells is a dominant force driving early amino acid replacement in HCV viral populations. This question was addressed in two chimpanzees followed for 8 to 10 years after infection with a well-defined inoculum composed of a clonal genotype 1a (isolate H77C) HCV genome. Detailed characterization of CD8(+) T cell responses combined with sequencing of recovered virus at frequent intervals revealed that most acute-phase nonsynonymous mutations were clustered in class I epitopes and appeared much earlier than those in the remainder of the HCV genome. Moreover, the ratio of nonsynonymous to synonymous mutations, a measure of positive selection pressure, was increased 50-fold in class I epitopes compared with the rest of the HCV genome. Finally, some mutation of the clonal H77C genome toward a genotype 1a consensus sequence considered most fit for replication was observed during the acute phase of infection, but the majority of these amino acid substitutions occurred slowly over several years of chronic infection. Together these observations indicate that during acute hepatitis C, virus evolution was driven primarily by positive selection pressure exerted by CD8(+) T cells. This influence of immune pressure on viral evolution appears to subside as chronic infection is established and genetic drift becomes the dominant evolutionary force.     

Compartmentalization of hepatitis C virus quasispecies in blood mononuclear cells of patients with mixed cryoglobulinemic syndrome

The aim of this study was to investigate the quasispecies heterogeneity of hepatitis C virus (HCV) in the plasma, cryoprecipitate, and peripheral lymphocytes of chronically infected HCV patients with mixed cryoglobulinemia (MC). We studied 360 clones from 10 HCV-positive patients with MC and 8 age-, gender- and HCV genotype-matched subjects with chronic HCV infection but without MC. A partial nucleotide sequence encompassing the E1/E2 region, including hypervariable region 1 (HVR1), was amplified and cloned from plasma, cryoprecipitates, and peripheral blood mononuclear cells (PBMC), and the genetic diversity and complexity and synonymous and nonsynonymous substitution rates were determined. Heterogeneous selection pressure at codon sites was evaluated. Compartmentalization was estimated by phylogenetic and phenetic (Mantel's test) approaches. The patients with MC had 3.3 times lower nonsynonymous substitution rates (1.7 versus 5.7 substitutions/100 sites). Among the subjects with HCV genotype 1, the MC patients had significantly less complexity than the controls, whereas the diversity and complexity were similar in the genotype 2 patients and controls. Site-specific selection analysis confirmed the low frequency of MC patients showing positive selection. There was a significant correlation between positive selection and the infecting HCV genotype. The quasispecies were less heterogeneous in PBMC than in plasma. Significant compartmentalization of HCV quasispecies was observed in the PBMC of four of nine subjects (three with MC) and seven of nine cryoprecipitates. In one subject with MC, we detected a 5-amino-acid insertion at codons 385 to 389 of HVR1. Our results suggest reduced quasispecies heterogeneity in MC patients that is related to a low selection pressure which is probably due to an impaired immune response, the HCV genotype, and/or the duration of the infection. The frequent HCV quasispecies compartmentalization in patients' PBMC suggests a possible pathogenetic significance.     

Sequence Variation in the Gene Encoding the Nonstructural 3 Protein of Hepatitis C Virus: Evidence for Immune Selection

To determine whether the persistent nature of hepatitis C infection is related to the emergence of antigenic variants driven by immune selection, we examined the sequence heterogeneity in a portion of the hepatitis C virus (HCV) nonstructural 3 (NS3) gene of a patient infected over the course of more than 2 years. By PCR amplification, cloning, and sequencing, we observed several variable and conserved regions in the NS3 segment of the HCV genome. All variable regions had higher ratios of nonsynonymous/synonymous mutations and encompassed immunodominant epitopes, and their locations were not essential to maintain the known function of HCV RNA helicase. In contrast, the regions that are critical for HCV RNA helicase activity were found to be conserved with lower heterogeneity or lower ratios of nonsynonymous/synonymous mutations, and none except one of these regions was encoded within immunodominant epitopes. Our results are consistent with immune selection of viral variants at the epitope and molecular levels that may enable HCV to evade host defenses over time. Plotting the relatedness of sequence variants revealed a star topology suggesting that a wild-type HCV sequence is maintained, unlike HIV. Received: 2 November 2000 / Accepted: 1 October 2001     

Escape mutations alter proteasome processing of major histocompatibility complex class I-restricted epitopes in persistent hepatitis C virus infection

Mutations in hepatitis C virus (HCV) genomes facilitate escape from virus-specific CD8+ T lymphocytes in persistently infected chimpanzees. Our previous studies demonstrated that many of the amino acid substitutions in HCV epitopes prevented T-cell receptor recognition or binding to class I major histocompatibility complex molecules. Here we report that mutations within HCV epitopes also cause their destruction by changing the pattern of proteasome digestion. This mechanism of immune evasion provides further evidence of the potency of CD8+ T-cell selection pressure against HCV and should be considered when evaluating the significance of mutations in viral genomes from persistently infected chimpanzees and humans.     

Viral Persistence, Antibody to E1 and E2, and Hypervariable Region 1 Sequence Stability in Hepatitis C Virus-Inoculated Chimpanzees

The relationship of viral persistence, the immune response to hepatitis C virus (HCV) envelope proteins, and envelope sequence variability was examined in chimpanzees. Antibody reactivity to the HCV envelope proteins E1 or E2 was detected by enzyme-linked immunosorbent assay (ELISA) in more than 90% of a human serum panel. Although the ELISAs appeared to be sensitive indicators of HCV infection in human serum panels, the results of a cross-sectional study revealed that a low percentage of HCV-inoculated chimpanzees had detectable antibody to E1 (22%) and E2 (15%). Viral clearance, which was recognized in 28 (61%) of the chimpanzees, was not associated with an antibody response to E1 or E2. On the contrary, antibody to E2 was observed only in viremic chimpanzees. A longitudinal study of animals that cleared the viral infection or became chronically infected confirmed the low level of antibody to E1, E2, and the HVR-1. In 10 chronically infected animals, the sequence variation in the E2 hypervariable region (HVR-1) was minimal and did not coincide with antibody to E2 or to the HVR-1. In addition, low nucleotide and amino acid sequence variation was observed in the E1 and E2 regions from two chronically infected chimpanzees. These results suggest that mechanisms in addition to the emergence of HVR-1 antibody escape variants are involved in maintaining viral persistence. The significance of antibodies to E1 and E2 in the chimpanzee animal model is discussed.     

Departure from neutrality at the mitochondrial NADH dehydrogenase subunit 2 gene in humans,but not in chimpanzees.

To test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution, nucleotide sequences were determined for the 1041 bp of the NADH dehydrogenase subunit 2 (ND2) gene in 20 geographically diverse humans and 20 common chimpanzees. Contingency tests of neutrality were performed using four mutational categories for the ND2 molecule: synonymous and nonsynonymous mutations in the transmembrane regions, and synonymous and nonsynonymous mutations in the surface regions. The following three topological mutational categories were also used: intraspecific tips, intraspecific interiors, and interspecific fixed differences. The analyses reveal a significantly greater number of nonsynonymous polymorphisms within human transmembrane regions than expected based on interspecific comparisons, and they are inconsistent with a neutral equilibrium model. This pattern of excess nonsynonymous polymorphism is not seen within chimpanzees. Statistical tests of neutrality, such as TAJIMA''s D test, and the D and F tests proposed by FU and LI, indicate an excess of low frequency polymorphisms in the human data, but not in the chimpanzee data. This is consistent with recent directional selection, a population bottleneck or background selection of slightly deleterious mutations in human mtDNA samples. The analyses further support the idea that mitochondrial genome evolution is governed by selective forces that have the potential to affect its use as a "neutral" marker in evolutionary and population genetic studies.     

Sequence evolution of the hypervariable region in the putative envelope region E2/NS1 of hepatitis C virus is correlated with specific humoral immune responses.

Sequence evolution of the hypervariable region 1 (HVR1) in the N terminus of E2/NS1 of hepatitis C virus (HCV) was studied retrospectively in six chimpanzees inoculated with the same genotype 1b strain, containing a unique predominant HVR1 sequence. Immediately after inoculation, all animals contained the same HVR predominant sequence. Two animals developed an acute self-limiting infection. Anti-HVR1 immunoglobulin G (IgG) was produced 40 to 60 days after inoculation and rapidly disappeared after normalization of transaminases. Another chimpanzee, previously infected with human immunodeficiency virus type 1, showed a delayed response to HVR1 epitopes after superinfection with HCV. No sequence variation of HVR1 was observed in these two animals during the transient viremia in the acute phase. Three other chimpanzees developed a chronic HCV infection. During follow up, sequence evolution occurred in two animals and their anti-HVR1 response remained at varying but detectable levels. The first mutations occurred immediately after the production of anti-HVR1 during the acute phase. However, IgM anti-HVR1 was not detectable. Remarkably, HVR1 sequences remained conserved for more than 6 years in another chronically infected animal. This correlated with the complete absence of detectable anti-HVR1 during this period. Seven years after inoculation, anti-HVR1 IgG was produced and coincided with an HVR1 alteration. These results strongly suggest the involvement of neutralizing anti-HVR antibodies in sequence evolution of HVR1 through immune selection.     

Comparison of DNA and protein polymorphisms between humans and chimpanzees

To examine the nucleotide diversity at silent (synonymous + intron + untranslated) and non-silent (nonsynonymous) sites in chimpanzees and humans, genes at six nuclear loci from two chimpanzees were sequenced. The average silent diversity was 0.19%, which was significantly higher than that in humans (0.05%). This observation suggests a significantly larger effective population size and a higher extent of neutral polymorphism in chimpanzees than in humans. On the other hand, the non-silent nucleotide diversity is similar in both species, resulting in a larger fraction of neutral mutations at non-silent sites in humans than in chimpanzees. Other types of polymorphism data were collected from the literature or databases to examine whether or not they are consistent with the nuclear DNA sequence polymorphism observed here. The nucleotide diversity at both silent and non-silent sites in mitochondrial (mt) DNA genes was compatible with that of the nuclear genes. Microsatellite loci showed a similar high extent of heterozygosity in both species, perhaps due to the combined effect of a high mutation rate and a recent population expansion in humans. At protein loci, humans are more heterozygous than chimpanzees, and the estimated fraction of neutral alleles in humans (0.84) is much larger than that in chimpanzees (0.26). These data show that the neutral fraction in non-silent changes is relatively large in the human population. This difference may be due to a relaxation of the functional constraint against proteins in the human lineage. To evaluate this possibility, it will be necessary to examine nucleotide sequences in relation to the physiological or biochemical properties of proteins.     

Human SNPs reveal no evidence of frequent positive selection

We compared the single-nucleotide polymorphisms (SNPs) in humans in 182 housekeeping and 148 tissue-specific genes. SNPs were divided into rare and common polymorphisms based on their frequencies. We found that housekeeping genes tend to be less polymorphic than tissue-specific genes for both rare and common SNPs. Using mouse as a second species for computing sequence divergences, we found no evidence of positive selection: for both housekeeping and tissue-specific genes, the ratio of nonsynonymous to synonymous common SNPs per site showed no significant difference from that of divergence. Similarly, we observed no evidence of positive selection for the 289 and 149 genes that have orthologs available for divergence calculation between humans and chimpanzees and between humans and Old World monkeys, respectively. A comparison with previous SNP studies suggests that approximately 20% of the nonsynonymous SNPs in the human population are nearly neutral and that positive selection in the human genome might not be as frequent as previously thought.     

Contrasting histories of G6PD molecular evolution and malarial resistance in humans and chimpanzees

Although mutations in the glucose-6-phosphate dehydrogenase (G6PD) gene result in several blood-related diseases in humans, they also confer resistance to malarial infection. This association between G6PD and malaria was supported by population genetic analyses of the G6PD locus, which indicated that these mutations may have recently risen in frequency in certain geographic regions as a result of positive selection. Here we characterize nucleotide sequence variation in a 5.2-kb region of the G6PD locus in a population sample of 56 chimpanzees, as well as among 7 other nonhuman primates, to compare with that in humans in determining whether other primates that are impacted by malaria also exhibit patterns of G6PD polymorphism or divergence consistent with positive selection. We find that chimpanzees have several amino acid variants but that the overall pattern at G6PD in chimpanzees, as well as in Old and New World primates in general, can be explained by recent purifying selection as well as strong functional constraint dating back to at least 30-40 MYA. These comparative analyses suggest that the recent signature of positive selection at G6PD in humans is unique.     

Sequence analysis of the hepatitis C virus genome recovered from serum, liver, and peripheral blood mononuclear cells of infected chimpanzees.

Of 13 different strains of hepatitis C virus (HCV) in the inoculum used, only 1 persisted in human lymphocyte cell lines infected in vitro (N. Nakajima, M. Hijikata, H. Yoshikura, and Y. K. Shimizu, J. Virol. 70:3325-3329, 1996). To determine whether that particular strain (designated H1-2) has a tropism for lymphocytes in vivo, we sequenced hypervariable region 1 (HVR1) of the genome of HCV recovered from the sera, livers, and peripheral blood mononuclear cells (PBMC) of chimpanzees infected with plasma H77, the same inoculum used for the in vitro studies. In the PBMC collected from two chimpanzees during the early phase of infection, H1-2 was detected as the only or predominant HVR1 sequence. H1-2 was also detected in PBMC obtained during persistent infection from a chimpanzee that had been treated with immunosuppressants. From the livers of these chimpanzees, two to six different strains were recovered but H1-2 was not detected. Thus, H1-2 appeared to have an affinity for lymphocytes not only in vitro but also in vivo. In samples collected from a chimpanzee after 6 years of infection, however, such tissue compartmentalization of the HCV genome was not observed; a single strain became predominant in the serum, liver, and PBMC. An HCV strain capable of replicating in both the liver and PBMC probably emerged during in vivo replication and persisted.     

Constraints on Viral Evolution during Chronic Hepatitis C Virus Infection Arising from a Common-Source Exposure

Extraordinary viral sequence diversity and rapid viral genetic evolution are hallmarks of hepatitis C virus (HCV) infection. Viral sequence evolution has previously been shown to mediate escape from cytotoxic T-lymphocyte (CTL) and neutralizing antibody responses in acute HCV infection. HCV evolution continues during chronic infection, but the pressures driving these changes are poorly defined. We analyzed plasma virus sequence evolution in 5.2-kb hemigenomes from multiple longitudinal time points isolated from individuals in the Irish anti-D cohort, who were infected with HCV from a common source in 1977 to 1978. We found phylogenetically distinct quasispecies populations at different plasma time points isolated late in chronic infection, suggesting ongoing viral evolution and quasispecies replacement over time. We saw evidence of early pressure driving net evolution away from a computationally reconstructed common ancestor, known as Bole1b, in predicted CTL epitopes and E1E2, with balanced evolution toward and away from the Bole1b amino acid sequence in the remainder of the genome. Late in chronic infection, the rate of evolution toward the Bole1b sequence increased, resulting in net neutral evolution relative to Bole1b across the entire 5.2-kb hemigenome. Surprisingly, even late in chronic infection, net amino acid evolution away from the infecting inoculum sequence still could be observed. These data suggest that, late in chronic infection, ongoing HCV evolution is not random genetic drift but rather the product of strong pressure toward a common ancestor and concurrent net ongoing evolution away from the inoculum virus sequence, likely balancing replicative fitness and ongoing immune escape.     

Effect of bottlenecking on evolution of the nonstructural protein 3 gene of hepatitis C virus during sexually transmitted acute resolving infection

Sexual partners of patients infected with the hepatitis C virus (HCV) often have detectable HCV-specific T-cell responses in the absence of seroconversion, suggesting unapparent, spontaneously resolving infection. To determine whether differences in the evolutionary potential of bottlenecked inoculum may explain the low rate of HCV persistence after sexual exposure, we have investigated changes in the entire HCV nonstructural 3 (NS3) gene over time in a chronic carrier and compared his viral quasispecies with that of the acute-phase isolate of his sexual partner, who developed acute resolving hepatitis C. The overall rate of accumulation of mutations, estimated by regression analysis of six consecutive consensus NS3 sequences over 8 years, was 1.5 x 10(-3) mutations per site per year, with small intersample fluctuations related to changes in environmental conditions. Comparison of quasispecies parameters in one isolate of the chronic carrier with those of the acute-phase isolate of the infected partner revealed a higher heterogeneity and lower proportion of nonsynonymous mutations in the former. All NS3 sequences from the acute-phase isolate clustered with a single sequence from the chronic isolate, despite complete HLA mismatch between the patients, suggesting bottlenecking during transmission. The low risk of viral persistence after sexual exposure to HCV may be related to the selection of a limited number of viral particles carrying a particular combination of mutations which may further limit the potential of a relatively homogeneous quasispecies to rapidly diversify and overcome the immune response of the exposed host.     

High-resolution phylogenetic analysis of hepatitis C virus adaptation and its relationship to disease progression

Hepatitis C virus (HCV) persists in the majority of those infected despite host immune responses. Evidence has accrued that selectively fixed mutations in the envelope genes (E1 and E2) are associated with viral persistence, particularly those that occur within the first hypervariable region of E2 (HVR1). However, the individual amino acid residues under selection have not been identified, nor have their selection pressures been measured, despite the importance of this information for understanding disease pathogenesis and for vaccine design. We performed a high-resolution analysis of published gene sequence data from individuals undergoing acute HCV infection, employing two phylogenetic methods to determine site-specific selection pressures. Strikingly, we found a statistically significant association between the number of sites selected and disease outcome, with the fewest selected sites in fulminant HCV cases and the greatest number of selected sites in rapid progressors, reflecting the duration and intensity of the arms race between host and virus. Moreover, sites outside the HVR1 appear to play a major role in viral evolution and pathogenesis, although there was no association between viral persistence and specific mutations in E1 and E2. Our analysis therefore allows fine dissection of immune selection pressures, which may be more diverse than previously thought. Such analyses could play a similarly informative role in studies of other persistent virus infections, such as human immunodeficiency virus.     

Expression of Noncovalent Hepatitis C Virus Envelope E1-E2 Complexes Is Not Required for the Induction of Antibodies with Neutralizing Properties following DNA Immunization

Interactive glycoproteins present on the surface of viral particles represent the main target of neutralizing antibodies. The ability of DNA vaccination to induce antibodies directed at such structures was investigated by using eight different expression plasmids engineered either to favor or to prevent interaction between the hepatitis C virus (HCV) envelope glycoproteins E1 and E2. Independently of the injection route (intramuscular or intraepidermal), plasmids expressing antigens capable of forming heterodimers presumed to be the prebudding form of the HCV envelope protein complex failed to induce any significant, stable antibodies following injection in mice. In sharp contrast, high titers of antibodies directed at both conformational and linear determinants were induced by using plasmids expressing severely truncated antigens that have lost the ability to form native complexes. In addition, only a truncated form of E2 induced antibodies reacting against the hypervariable region 1 of E2 (specifically with the C-terminal part of it) known to contain a neutralization site. When injected intraepidermally into small primates, the truncated E2-encoding plasmid induced antibodies able to neutralize in vitro the binding of a purified E2 protein onto susceptible cells. Because such antibodies have been associated with viral clearance in both humans and chimpanzees, these findings may have important implications for the development of protective immunity against HCV.     

False-positive results obtained from the branch-site test of positive selection

Natural selection operating at the amino acid sequence level can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitutions, where r(N)/r(S) (omega) > 1 and omega < 1 suggest positive and negative selection, respectively. The branch-site test has been developed for detecting positive selection operating at a group of amino acid sites for a pre-specified (foreground) branch of a phylogenetic tree by taking into account the heterogeneity of omega among sites and branches. Here the performance of the branch-site test was examined by computer simulation, with special reference to the false-positive rate when the divergence of the sequences analyzed was small. The false-positive rate was found to inflate when the assumptions made on the omega values for the foreground and other (background) branches in the branch-site test were violated. In addition, under a similar condition, false-positive results were often obtained even when Bonferroni correction was conducted and the false-discovery rate was controlled in a large-scale analysis. False-positive results were also obtained even when the number of nonsynonymous substitutions for the foreground branch was smaller than the minimum value required for detecting positive selection. The existence of a codon site with a possibility of occurrence of multiple nonsynonymous substitutions for the foreground branch often caused the branch-site test to falsely identify positive selection. In the re-analysis of orthologous trios of protein-coding genes from humans, chimpanzees, and macaques, most of the genes previously identified to be positively selected for the human or chimpanzee branch by the branch-site test contained such a codon site, suggesting a possibility that a significant fraction of these genes are false-positives.     

Human immunodeficiency virus seroconversion and evolution of the hepatitis C virus quasispecies

When chronic hepatitis C virus (HCV) infections are complicated by acquisition of human immunodeficiency virus (HIV), liver disease appears to accelerate and serum levels of HCV RNA may rise. We hypothesized that HIV might affect the HCV quasispecies by decreasing both complexity (if HIV-induced immunosuppression lessens pressure for selecting HCV substitutions) and the ratio of nonsynonymous (d(N)) to synonymous (d(S)) substitutions, because d(N) may be lower (if there is less selective pressure). To test this hypothesis, we studied the evolution of HCV sequences in 10 persons with chronic HCV infection who seroconverted to HIV and, over the next 3 years, had slow or rapid progression of HIV-associated disease. From each subject, four serum specimens were selected with reference to HIV seroconversion: (i) more than 2 years prior, (ii) less than 2 years prior, (iii) less than 2 years after, and (iv) more than 2 years after. The HCV quasispecies in these specimens was characterized by generating clones containing 1 kb of cDNA that spanned the E1 gene and the E2 hypervariable region 1 (HVR1), followed by analysis of clonal frequencies (via electrophoretic migration) and nucleotide sequences. We examined 1,320 cDNA clones (33 per time point) and 287 sequences (median of 7 per time point). We observed a trend toward lower d(N)/d(S) after HIV seroconversion in 7 of 10 subjects and lower d(N)/d(S) in those with rapid HIV disease progression. However, the magnitude of these differences was small. These results are consistent with the hypothesis that HIV infection alters the HCV quasispecies, but the number of subjects and observation time may be too low to characterize the full effect.     

The clearance of hepatitis C virus infection in chimpanzees may not necessarily correlate with the appearance of acquired immunity

Clearance of hepatitis C virus (HCV) infection in humans and chimpanzees is thought to be associated with the induction of strong T-cell responses. We studied four chimpanzees infected with HCV derived from an infectious full-length HCV genotype 1b cDNA. Two of the chimpanzees cleared the infection to undetectable levels for more than 12 months of follow-up; the other two became persistently infected. Detailed analyses of HCV-specific immune responses were performed during the courses of infection in these chimpanzees. Only weak and transient T helper responses were detected during the acute phase in all four chimpanzees. A comparison of the frequency of gamma interferon (IFN-gamma)-producing CD4(+) and CD8(+) T cells in peripheral blood by ELISpot assay did not reveal any correlation between viral clearance and T-cell responses. In addition, analyses of IFN-gamma, IFN-alpha, and interleukin-4 mRNA levels in liver biopsies, presumably indicative of intrahepatic T-cell responses, revealed no distinct pattern in these chimpanzees with respect to infection outcome. The present study suggests that the outcome of HCV infection in chimpanzees is not necessarily attributable to HCV sequence variation and that chimpanzees may recover from HCV infection by mechanisms other than the induction of readily detectable HCV-specific T-cell responses.     

Gene diversity patterns at 10 X-chromosomal loci in humans and chimpanzees

We have investigated the pattern and extent of nucleotide diversityin 10 X-chromosomal genes where mutations are known to causemental retardation in humans. For each gene, we sequenced theentire coding region from cDNA in humans, chimpanzees, and orangutans,as well as about 3 kb of genomic DNA in 20 humans sampled worldwideand in 10 chimpanzees representing two "subspecies." Overallnucleotide diversity in these genes is about twofold lower inhumans than in chimpanzees, and nucleotide diversity withinand between species is low, suggesting that a high level offunctional constraint acts on these genes. Strikingly, we findthat a summary of the allele frequency spectrum is significantlycorrelated in humans and chimpanzees, perhaps reflecting verysimilar levels of constraint at these genes in the two species.A possible exception is FMR2, which shows a higher number ofnonsynonymous than synonymous substitutions on the human lineage,suggesting the action of positive selection.