|Home | About | Journals | Submit | Contact Us | Français|
Little is known about hepatitis B virus (HBV) diversity changes within a host during the immunotolerant phase of chronic HBV infection. Such knowledge, nevertheless, may help in understanding how host immunity and HBV interact at the early stage of infection. In this study, serial serum samples were collected from a long-term (>17 years) follow-up cohort of seven patients, and multiple copies of the full-length viral genome from serially sampled sera were recovered and analyzed. Viral genetic diversity was positively correlated with host immunity, represented by levels of alanine aminotransferase (ALT), but was negatively correlated with the viral copy number. During the immunotolerant phase, when the host immunity was feeble (ALT < 20 U/liter), viral nucleotide diversity decreased while copy numbers increased. Rates of evolutionary change derived for different patients were in a very narrow range (1.6 × 10−5 to 5.4 × 10−5/site/year). As the disease progressed toward the immunoclearance phase (ALT > 20 U/liter), viral diversity increased but copy numbers decreased. Evolutionary rates varied among patients in accordance with their levels of ALT, ranging from 9.6 × 10−6 to 3.2 × 10−4/site/year. More than half (19/32 sites) of positively selected sites resided in immune epitopes, suggesting their possible role in host immunity. Our results demonstrate that host immunity is a dominant factor in HBV evolution. Different selective forces, including immune-mediated positive selection and virus-mediated negative selection, operate in tandem in shaping viral population dynamics within a host.
In the study of chronic viral disease, the viral diversity change is established as an important factor in the pathogenesis of hepatitis C virus (HCV) and human immunodeficiency virus (HIV) (14, 23, 30, 34). Because hepatitis B virus (HBV) uses reverse transcriptase for replication, the error-prone nature of polymerase makes the rate of evolutionary change much closer to those of RNA viruses than to those of DNA viruses (21). From this aspect, it was proposed that the evolutionary dynamics of HBV should be similar to those of HCV and HIV (48). Nevertheless, little is known about HBV dynamics in chronic HBV infection (CHI), especially in the early stage of infection.
While mutations introduced during replication are an inherent and common driver of all viral variations, the viral generation time, defined as the average length of time from the release of a virion until it infects another cell and causes the release of another virion, was proposed as another factor determining the long-term substitution rate (18). The HBV generation time was estimated to range from 10 to 100 days, with an average of 24.8 days (28), which is 6 and 10 times longer than those of HCV (27) and HIV (31), respectively. Therefore, many features of viral dynamics in CHI may differ from those of other viral diseases. In addition, CHI is distinctive from other chronic viral diseases by having a prolonged immune tolerance phase without an acute reaction if an individual is infected perinatally. During the early stage of CHI, HBV quasispecies are not or are only minimally subjected to host immune selection (3, 6). In teenaged patients, CHI enters an immunoclearance phase with an active host immune response to control viral replication (5). Studies of CHI found associations of viral variants, including nucleotide and structural variations such as insertions and deletions, with disease progression (4). However, previous studies focused on the relatively later period of infection, e.g., during active hepatitis close to or after hepatitis Be antigen (HBeAg) seroconversion (usually after the age of 20 years if infected perinatally) or during the development of cirrhosis or hepatocellular carcinoma (HCC) (usually after the age of 40 years) (37). The evolution of HBV in the immunotolerant phase during the first 2 decades within a host is largely unknown.
To study the evolution of viral populations in this stage of CHI, we selected a cohort of carriers who have been monitored regularly for more than 17 years, beginning before the age of 10 years. For each patient, multiple full-length HBV sequences were recovered from three time points, spanning the immunotolerant phase without alanine aminotransferase (ALT) elevation to the early immunoclearance phase with mild ALT elevation. We estimated the viral nucleotide diversity, mutation rate, and effective population size within hosts between the two stages. Based on our results, a scenario of HBV evolution during the early stage of infection is proposed.
Seven patients were selected from the Department of Pediatrics, National Taiwan University Hospital (NTUH). Among them, six patients were presumed to have acquired HBV perinatally, as their mothers were hepatitis B surface antigen (HBsAg) positive at the time of enrollment. Only patients with genotype B were used to ensure that differences found in viral evolution were not due to genotypic differences. These patients have been monitored for an average of 18.3 years, starting before the age of 10 years. Their ALT, HBsAg, HBeAg, anti-HBs, and anti-HBe levels were monitored regularly at 6-month intervals. A guideline was followed for ALT measurement to ensure that the measurements were consistent during the last 2 decades (see the supplemental material for details). These patients were all HBsAg and HBeAg positive, with no HBeAg seroconversion observed during the study. No antiviral therapy was administered throughout the course of this study. All patients were enrolled by community surveys before their becoming aware of CHI. Therefore, this cohort is ideal for studying HBV evolution under natural conditions. We selected the following three time points to survey the evolutionary changes in HBV quasispecies: I, the time of enrollment; III, the most recent clinical visit; and II, one time point between points I and III, before the elevation of ALT (see Results).
For each patient, HBV DNA was extracted, and the virus titer was quantified from 100 μl of serum sample at the core facility of NTUH. We first amplified the approximately 3.2-kb HBV genome. The PCR conditions and primers used are listed in Table S1 in the supplemental material. Amplification was performed with Stratagene Easy-A High Fidelity PCR cloning enzyme (Stratagene) and 5 ng of DNA template in a 50-μl reaction mix. PCR products were purified with a PCR-M purification kit (Viogene, Taipei, Taiwan), cloned into a Topo XL PCR cloning vector (Invitrogen), and transformed into competent cells according to the manufacturer's protocols. Clones were picked for a PCR check, and 12 positive clones per cloning reaction were subjected to sequencing. We designed six internal primers for full-length HBV sequencing (see Table S1 in the supplemental material). The sequencing reactions were performed using a Big Dye Terminator kit and then loaded into a model 3730 xl sequencer (Applied Biosystems) at the core facility of NTUH. The base positions were numbered according to a reference sequence (GenBank accession no. NC_003977).
Sequences were assembled using the Seqman program of the Lasergene package (Lasergene, Madison, WI), with visual inspection, and then aligned using Clustal_W (38). Viral nucleotide diversity (θ) was calculated by custom C++ software based on the Libsequence package (39). θ represents the number of variable nucleotide sites in a sample of sequences divided by sequence length. For nonoverlapping regions of HBV, θ was calculated separately for silent (synonymous; θS) and amino acid-altering (nonsynonymous; θA) sites.
The evolutionary rate (μ), given as the number of changes per site per year, and the viral effective population size (Ne) were estimated using an established Bayesian Markov chain Monte Carlo (MCMC) approach (11) that incorporates the time of viral sampling. This approach estimates the rate of nucleotide change and dynamics of the population genetic diversity from serially sampled sequence data (10). In addition, an estimate of the posterior distribution of genealogies related to the sequence data was concurrently inferred for each alignment. From these distributions, the times to the origin for viruses isolated from each patient at each time point (the time to the most recent common ancestor [tMRCA]) were obtained. The analyses were performed using the GTR+Γ model of nucleotide substitutions with a relaxed molecular clock model supplemented in BEAST (11). We performed two independent runs with 2 × 106 MCMC steps, of which the first 10% were discarded as burn-in. The results were compared to confirm that both sampled the same distribution and then were combined. The rate of HBV evolution (μ) was also estimated by comparing mean pairwise nucleotide differences within and between DNA samples recovered from different time points (16).
To detect positive selection, two approaches were applied. Selection can be detected by comparing numbers of synonymous changes per synonymous site (Ks) and numbers of nonsynonymous changes per nonsynonymous site (Ka). The ratio Ka/Ks, or ω, is used to assess selection. A Ka/Ks ratio of >1 indicates positive selection, whereas a Ka/Ks ratio of <1 indicates negative selection. To identify the putative amino acid residues under positive selection, two pairs of likelihood-based models implemented in the PAML package (46) were used. The first pair includes M1a (nearly neutral) and M2a (positive selection) (45), while the second pair includes M7 (beta) and M8 (beta and ω) (46). The result was given after running different models in codeml. Because some parts of the HBV genome encode more than one protein via different open reading frames (ORFs), these regions are not suitable for the Ka/Ks-based test. Therefore, we applied the codon-based test only to nonoverlapping coding regions of the HBV genome.
Amino acid sites experiencing positive selection in which a mutation to one particular amino acid is advantageous can also be inferred when a mutant accumulates in the population faster than expected by random genetic drift. We carried out a computer simulation to identify these sites and to estimate the selective coefficient. The details of the simulation are given in the supplemental material. In short, we assumed that HBV's behavior could be described by a discrete-generation model. Supposing that from time t1 to t2 the frequency of a mutant changes from f1 to f2, then under conditions of drift, the probability can be formulated as follows: P(time = t2 − t1 | = f1, = f2, n = Ne) (26).
Clinical information for all patients is listed in Table Table11 and includes the age at time of enrollment and HBV sampling, the average ALT levels between different sampling points, and virus titers. The mean HBV DNA loads were all higher than 107 copies/ml, indicating active viral replication. In general, the average ALT level between the second and third points was twofold greater than that between the first and second points. The ALT profiles during the follow-up and the nucleotide diversities (θ) at different sampling points are shown in Fig. Fig.1.1. Although θ varied in different individuals, a general trend emerged when the values were aligned. Specifically, except for patient F, θ values at the second point were lower than those derived from the first (P = 0.008; Wilcoxon signed rank paired test) and third (P = 0.021) points (Table (Table2).2). Namely, HBV diversity was high at a “relative early” stage of infection, decreased during the immune tolerance phase, and rose again toward the immunoclearance phase with elevated ALT levels (Fig. (Fig.1).1). Nucleotide diversities at synonymous (θS) and nonsynonymous (θA) sites in nonoverlapping regions of the HBV genome were estimated separately. While θA was significantly reduced (P = 0.032) from the first to the second points, θS was not (P = 0.344). From the second to the third points, both θA and θS significantly increased, with P values of 0.021 and 0.023, respectively.
Several studies proposed that selective pressures exerted by both host immunity and viral load are important factors influencing viral diversity in chronic viral diseases (14, 25, 36). We next examined the relationships among them. Since cellular immunity assays were not available for this analysis, we used the ALT level as an indicator to reflect host immune activity, as it was shown that ALT level provides some estimate of the strength of an immunological response against viral infection (14, 28). There was only a weak correlation between nucleotide diversities and spot ALT levels at the time the sequence was sampled (P = 0.112; Kendall's rank sum correlation) (see Fig. S1 in the supplemental material). Since ALT can be influenced by many factors and since changes in diversity between different time points may result from the cumulative effect of immune activity, a single point measurement might not represent the long-term immune status of a host. We therefore used average ALT values between different time points to assess the correlation. By using the average, we assumed that the stochastic effect could be minimized, and the value itself should represent the long-term trend of host immunity. Because there was no ALT information before enrollment, we plotted only the average ALTs and HBV nucleotide diversities derived from the second and third points. The correlation was highly significant (P = 0.002 and tau = 0.604) (Fig. (Fig.2A).2A). To further decide whether this correlation was contributed by all sites in the genome or by specific sites or regions, we assessed the correlation for nucleotide diversities at synonymous (θS) and nonsynonymous (θA) sites, separately. The average ALT level was significantly correlated with θA (P = 0.002 and tau = 0.604) (Fig. (Fig.2B),2B), but the correlation between the average ALT level and θS was very weak (P = 0.101) (Fig. (Fig.2C2C).
There was a negative correlation between the viral load (log10 copies/ml) and θ, i.e., the higher the viral load the lower the nucleotide diversity, with P = 0.016 and tau = −0.381 (Fig. (Fig.3).3). This negative relationship was also observed when synonymous and nonsynonymous sites were analyzed separately (see Fig. S2 in the supplemental material).
The HBV evolutionary rates (μ) showed great variations among patients (Table (Table3).3). Evolutionary rates derived by the pairwise difference method (16) were lowest for patient D (1.1 × 10−6/site/year) and highest for patient G (1.2 × 10−4/site/year). Applying a Bayesian MCMC approach, μ ranged from 8.0 × 10−6/site/year for patient B to 2.1 × 10−4/site/year for patient G. Most of the 95% highest posterior densities (HPDs) from different patients did not overlap, indicating that variations in μ among patients were significant. The evolutionary rates yielded by different methods were in good agreement with each other, with a correlation coefficient (R2) of 0.67. We thus consider our estimations to be authentic. The viral effective population size (Ne) was largest in patient B, at 4.7 × 105, and smallest in patient G, at 1.6 × 103 (Table (Table3).3). The times to the origin for viruses isolated from each patient (tMRCAs) ranged from 20 years (95% HPD, 17 to 23 years) in patient E to 235 years (95% HPD, 23 to 2,051 years) in patient F. Most patients, except patients E and G, had tMRCAs longer than their age at the most recent clinical visit, suggesting that part of the viral diversity might have been inherited from the donor (most likely from the mother through vertical transmission).
We next estimated Ne and μ for different intervals, from the first to the second time points and from the second to the third time points. Because in some cases μ values calculated from different time intervals by the pairwise difference method were negative due to small differences between different samplings, we applied only the evolutionary rate derived from the Bayesian approach. μ values from the first and second time points were in a very narrow range, with the lowest in patient B, at 1.6 × 10−5/site/year, and the highest in patient G, at 5.4 × 10−5/site/year, and their 95% HPDs largely overlapped. The effective population sizes ranged from 1.8 × 104 (patient A) to 1.5 × 105 (patient B). Furthermore, although the 95% HPDs covered a wide time interval, the tMRCAs were all greater than the ages of the patients at the time of the second sampling. In contrast, rates of change derived from the second to third time points showed great differences among patients, with the highest in patient G, at 3.2 × 10−4/site/year, and the lowest in patient F, at 9.6 × 10−6/site/year. Effective population sizes exhibited >100-fold differences among patients, from the largest of 2.2 × 105 (patient F) to the smallest of 9.0 × 102 (patient G). In addition, patients A, E, and G had tMRCAs shorter than their ages at the third sampling.
Our computer simulation identified 36 nucleotide changes altering 38 amino acids, as two nucleotides changed the amino acids in two reading frames, but covering only 32 amino acid sites, as some of them arose in different individuals independently or two mutations occurred in the same codon, with significantly higher frequencies of change (P < 0.05) than expected by random genetic drift (see Materials and Methods) (Table (Table4).4). Of these amino acid sites, 19 are in known T-cell epitopes (9), 1 causes the precore stop mutation, and the other 12 are in regions with no known immune epitopes. Since the HBV genome contains four ORFs encoding 1,613 amino acids, with 621 of these in known T-cell epitopes, the proportion of putative selected sites residing in epitopes was significantly higher than that in nonepitope regions (P = 0.014; Fisher's exact test) (Table (Table5).5). In addition, we noted that the nonstructural protein, i.e., polymerase, exhibited smaller numbers of selected sites than structural proteins, i.e., surface and core, with P values of 0.019 and 0.032, respectively. (The X protein was difficult to categorize by this structural and nonstructural dichotomy, as its function is not fully understood .)
Among 36 mutations, the frequencies of 6 increased between the first and second time points, those of 29 increased between the second and third time points, and that of 1 increased in both intervals (Table (Table4).4). Since patient-years from the first to second and from the second to third time points were 49.3 and 78.8 years, respectively, significantly (P = 0.014) greatly accelerated amino acid changes between the second and third time points correlated well with elevated ALT levels during this period. The selection coefficients of different mutants ranged from 0.007 to 0.036, with an average of 0.020. There were five nucleotide changes (altering six amino acids) for which selection coefficients could not be given, because they reached fixation within the host during follow-up. We added a “greater-than” symbol to their coefficients to indicate that their selection coefficients may have been greater than those numbers.
Selective amino acid sites were also assessed by the likelihood approach implemented in PAML (Table (Table6).6). In total, two likelihood ratio tests recognized nine amino acid sites under selection, with ω values of >1, and eight of these sites overlapped in the results of our computer simulation. Thus, there were a total of 33 amino sites putatively under positive selection.
During a chronic viral infection, such as HCV or HIV infection, the quasispecies are usually largely reduced upon transmission, probably due to host immune selection or severe genetic drift, and gradually restored after the emergence of escape mutants (12, 14). Therefore, the reduction in viral diversity after a decade of CHI is very peculiar. Significant diversity reduction at nonsynonymous rather than synonymous sites resembles the outcome of natural selection, as it is the only evolutionary force which targets a specific region or site. Selection pressures for the virus come from either competition between viral variants, i.e., differences in replicative efficiency, host immune activity, or both. These possibilities, unfortunately, are confounded under most circumstances. By taking advantage of the immune tolerance phase of CHI, we may now dissect their relative contributions in more detail.
It is generally assumed that there is no immune selection during the tolerance phase of CHI (6). This notion is partially supported by the observation that very few infiltrative lymphocytes were found in the livers of asymptomatic HBV carrier children infected in the perinatal period (3). In addition, on average, ALT levels in hepatitis-free populations of different age groups were never >20 U/liter in previous surveys (13, 41). Therefore, low levels of ALT (<20 U/liter) before the second sampling suggest inactive or minimal immune activity and, in turn, imply small or no selective pressure from the host. During this phase, viral strains with higher replicative ability increased their copy numbers and outcompeted those carrying amino acid-altering mutations with lower replicative efficiencies. Consequently, nucleotide diversity at nonsynonymous sites (θA) was reduced with an increase in viral load. From the second to third time points, the elevated ALT level (>20 U/liter) indicated increased host immunity. Positive selection would favor different mutations under host immune pressures (44) and would result in an increased θA and accelerated amino acid changes.
In addition to selection, nucleotide diversity (θ) can also be modified in accordance with variations in viral population size and the mutation rate (μ). Since θ and the viral copy number are negatively correlated (Fig. (Fig.3)3) and a reduced effective population size (Ne) during this period was shown (Table (Table3),3), the only factor left is the change in μ, which is measured in generation time. The HBV generation time is largely decided by the half-lives of infected cells, which range from about 10 to 100 days and are strongly influenced by the strength of a host's immunity (28). If we assume that the replication error introduced by viral polymerase is constant over time, an almost 10-fold reduction in generation time is equivalent to a 10-fold increase in μ measured on a yearly basis. Namely, from the second to the third time point, the increased host immunity speeded up the HBV life cycle by shortening its generation time, which in turn accelerated the rate of change and increased viral nucleotide diversity at both nonsynonymous and synonymous (θS) sites.
Although the above interpretation explains the observed viral diversity changes in CHI during the follow-up, one issue which remains to be determined is when and how the viral diversity at the first time point was generated. It is possible that viral diversity was inherited from the donor (most likely the mother) and gradually decreased during the immune tolerance phase because of viral competition. The estimated times of viral origin (tMRCAs) for most patients were longer than the ages at which HBV was sampled (Table (Table3),3), which seems to favor this scenario. However, since the lower bounds of 95% confidence intervals (95% HPDs) of tMRCA for most patients were generally smaller than their ages, the possibility that diversity was generated during HBV colonization cannot be ruled out. Within 3 to 6 months of the initial infection, HBV can infect 40% to almost 100% of ~2 × 1011 hepatocytes and produce virus (3, 20, 28). This enormous population expansion should also be accompanied by an inflation in viral nucleotide diversity. The above two scenarios are, of course, not mutually exclusive, and both suggest that the reduction in nucleotide diversity may have occurred long before the follow-up. Nevertheless, their relative significances have important consequences for the long-term adaptation of HBV. The second hypothesis implies that only a small portion of a diverse population was transmitted and that the majority was lost, i.e., there was a genetic bottleneck upon HBV transmission. This may cause implications for preadapted viral strains during colonization. By lowering the chance of their transmission, genetic drift has the potential to prevent the accumulation of adaptive changes in a population (47), thereby impeding the long-term adaptation of HBV. If the diversity in the population is inherited from the donor, such a reduction in fitness is not necessary, and the adaptive mutants have a greater chance of being transmitted (12). Further studies focusing on direct mother-to-infant transmission may help to resolve this issue.
Contrasting with the previous finding that HBeAg-positive hosts with normal ALT levels have little viral evolution (33, 43), we identified 33 amino acid sites which are potential candidates for immune selection. Although immune-selected mutations typically require the demonstration that a naturally occurring mutated T-cell epitope fails to trigger an immune response, several lines of evidence suggest that they are, one way or another, responsible for virus and host immune interactions. First, significantly more selected sites resided in immune epitopes than in nonepitope regions. Second, compared with the nonstructural (polymerase) protein, more selected sites resided in structural (surface and core) proteins (Table (Table5).5). Finally, the occurrences of accelerated amino acid changes corresponded well with changes in ALT profiles. Since these patients were all HBeAg positive during the time of the study, our analysis provides evidence that immune selection may have occurred long before the observed HBeAg seroconversion. In addition, since the average ALT level for six of seven patients was <40 U/liter but >20 U/liter, these results illustrate the occurrence of host immune selection even with a slight increase in average ALT level. Although the cutoff of 20 U/liter may seem arbitrary, the ALT levels of hepatitis-free populations of different age groups have never exceeded this level (13, 41). To this extent, an ALT level of 20 U/liter may represent the action of host immunity against HBV.
Compared to the case for HIV, the evidence for HBV escape mutants due to host immunity is limited (25, 33, 44). This discrepancy is probably rooted in the weaker immune response and longer generation time in the latter than the former. The estimated selection coefficients were all <0.05, and a sequence usually contained more than one mutation in different epitopes, supporting previous observations that T-cell responses to HBV are weak (7, 43) but probably multiply specific (44). Since the generation time of HBV is, on average, nine times longer than that of HIV, then if the population size and selective coefficient are equal, the rate of frequency change for an advantageous mutation would be nine times lower in the former than in the latter. Therefore, multiple sequences derived from long-term follow-up of individuals are necessary to identify such mutants. In addition, because of the compact nature of its genome, it should be reasonable to challenge HBV's ability to tolerate many mutations simultaneously under multiple selection regimens. In contrast, a small number of changes among several epitopes may help HBV to gain a slight advantage under feeble immune pressure. In order to identify possible escape mutants, whole-genome sequencing is thus necessary. Short-term follow-up which focuses only on the dominant variant or part of the HBV genome without exploring the full spectrum of diversity of viral quasispecies may fail to recognize such mutants.
Mutations at the same amino acid sites emerged independently from different hosts in the late stage of this study (the third time point), implying their advantage under immune pressure. However, these mutations were absent from the early stages of this study (the first and second time points), suggesting that their superiority may have either been compromised (44) in the absence of host immune pressure or been reversed under evolving host immunity (24). For example, an N118T (or T118N) mutation in the X protein was due to a transversion of adenine (A) to cytosine (C) at position 1726, which is within the enhancer II region of precore mRNA. It has been shown that mutations at this position not only alter the immune epitope but also change the viral load (8, 17). Therefore, complex selection forces, including immune-mediated positive selection and virus-mediated negative selection, operate in tandem in shaping HBV evolution.
In CHI, according to different estimations, 5% to almost 100% of ~2 × 1011 hepatocytes can be infected with about 25 to 50 covalently closed circular DNA (cccDNA) per cell, which translates to at least 1011 to 1013 replicating particles per day (3, 20, 28). Thus, the estimated HBV effective population sizes (Ne), ranging from 103 to 105, are 106- to 1010-fold smaller than the actual size. Several factors may cause a reduction in the effective population size. One of the most important reasons is that continual selection reduces viral Ne. Based on our analysis, both positive and negative selection pressures were detected at different phases of the disease. HBV is apparently under tremendous selective pressure during an infection.
Another possible explanation is that HBV populations do not evolve randomly (panmixia) but, rather, with some population structure. There is evidence that HBV populations are spatially structured, because in healthy chronic carriers, HBsAg-positive cells were usually separated from each other, with intervening areas of HBsAg-negative hepatocytes, and rarely occurred in large clumps (2, 22, 32). Therefore, HBV populations might exhibit a metapopulation structure. In this model, HBVs within the host are divided into local populations (hepatocytes here), and extinct populations (via cell death) are recolonized (regenerated cells) by viral particles from other local populations. In CHI, depending on the strength of the host's immune activity, about 0.3% to 3% of hepatocytes are killed and replenished every day (28, 40). It is known that such patterns of colonization and replacement can dramatically reduce Ne (35, 42), because only small fractions of local populations (hepatocytes) contribute to the long-term intrahost population of HBV.
Finally, changes in population size over time are another possible cause of a small Ne. When population size fluctuates, Ne is dominated by the smallest population size in the process, because it is given by the harmonic mean of population sizes at different stages (19). In our study, patient B had an average ALT level that was higher than that of patient E but exhibited a larger Ne. Despite the ALT measurement, however, patient B had a smaller standard deviation, which meant fewer fluctuations in host immune activity, and perhaps in viral population size, and was able to maintain a larger Ne than patients with similar average ALT levels but larger standard deviations, e.g., patient E. It seems likely that a combination of factors acts together to reduce the intrahost effective population size of HBV.
Although evolutionary rates estimated for different patients varied dramatically, they actually captured the entire spectrum of HBV mutation rates estimated from previous studies, from the lowest of 1.5 × 10−5/site/year to the highest of 7.7 × 10−4/site/year (15, 29, 48). The results shown in Tables Tables11 and and33 demonstrate how HBV evolutionary rates varied in accordance with host immune status. Immune pressure may influence HBV evolutionary rates in two ways. First, it can drive advantageous (supposedly escape) mutants to high frequencies or fixation in a short time, which increases the evolutionary rate. Second, an elevated immune reaction accelerates hepatocyte death and replenishment and shortens the HBV generation time, which, in turn, increases the evolutionary rate. Therefore, it is important to have knowledge of the immune status in advance before explaining the estimated rate of change. Moreover, HBV evolutionary rates derived from patients at different disease stages or with different immune statuses should be compared with caution.
In summary, the interplay between viral replication and host immunity explains the pattern of HBV dynamics within the host during a relatively early stage of infection. That is, without immune selection, competition between peers increases the viral load and decreases nucleotide diversity; in contrast, host immunity accelerates viral evolution and decreases copy numbers but increases diversity. Our observation, if true, may have important implications for the study of long-term HBV adaptation and therapeutic design for CHI.
We thank Lu Yun, Eleanore Chen, and Shi-Hui Chen for technical support. We also thank Shiou-Hwei Yeh, Chun-Jen Liu, and Taichung Tseng for constructive comments.
This study was supported by NHRI grants to D.-S.C. (NHRI-EX97-CD9201), M.-H. Chang (NHRI-EX97-9418BI), and H.-Y.W. (NHRI-EX99-9934BI).
Published ahead of print on 20 January 2010.
†Supplemental material for this article may be found at http://jvi.asm.org/.