Previous Article | Next Article ![]()
Antimicrobial Agents and Chemotherapy, February 2009, p. 615-621, Vol. 53, No. 2
0066-4804/09/$08.00+0 doi:10.1128/AAC.00947-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
for the Virahep-C Study
Department of Epidemiology, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania,1 Division of Infectious Diseases, Department of Medicine, School of Medicine, University of Pittsburgh, Pittsburgh, Pennsylvania,2 Department of Biostatistics, Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania,3 Roche Molecular Systems, Inc., Alameda, California,4 Cedars Sinai Medical Center, Los Angeles, California,5 University of Colorado Health Science Center, Denver, Colorado,6 Liver Diseases Branch, National Institute of Diabetes and Digestive and Kidney Diseases, Bethesda, Maryland7
Received 8 July 2008/ Returned for modification 15 August 2008/ Accepted 3 October 2008
|
|
|---|
2a plus ribavirin. Baseline viral levels were similar by race, irrespective of allele carriage. However, the rate of change in the viral decline was associated with both allele and race. Among the four subgroups defined by race and specific allele, the fastest rates of decline were observed (in terms of estimated mean viral declines log10 IU/ml during the first four weeks) in CA noncarriers for A*03 (2.75; P = 0.018), in CA carriers for Cw*03 (2.99; P = 0.046), and in CA noncarriers for DQA1*04 (2.66; P = 0.018) or DQB1*0402 (2.65; P = 0.018). MHC alleles are associated with the viral decline during the first 28 days of peginterferon therapy. |
|
|---|
Early disappearance of viremia during the first 4 weeks after the start of treatment is believed to be highly predictive of the ultimate sustained response to antiviral therapy. In a recent study, 91% of individuals who were PCR negative for HCV-RNA by week 4 achieved sustained virologic response (SVR), whereas only 60% of patients who had a <2-log10 drop in HCV levels by week 4 but were PCR negative for HCV by week 12 achieved SVR (7). The importance of the early viral response is further underscored in studies of the kinetics of viral decline, where multiple and very frequent measurements of HCV levels are obtained during the initial hours after the start of treatment (9, 19). Rapid rates of viral decline are associated with SVR (undetectable levels of HCV 6 months after discontinuation of therapy). Accordingly, viral dynamics may serve as a potential predictor of SVR (18, 19). Racial differences have also been reported with respect to viral dynamics, with AAs exhibiting slower viral dynamic declines than CAs (16, 17). The mechanisms behind the different dynamics of the viral decline are not known, but host genetic variation may contribute.
Spontaneous clearance of HCV is associated with vigorous and robust CD4+ and CD8+ T cells (3), and several studies have indicated an important role for T cells in antiviral therapy-induced clearance (6, 26, 28), although other studies have failed to find an association (14). The genes of the human major histocompatibility complex (MHC) encode the highly polmorphic human leukocyte antigens (HLA), which are crucial in antigen presentation and mediation of the T-lymphocyte response. Studies of the potential role of host genetic diversity and the dynamics of the early viral decline have been hampered by the difficulty in obtaining sequential HCV viral levels in a large cohort of individuals. We took advantage of the availability of repeated viral level data from a clinical study of AA and CA patients undergoing anti-HCV therapy and examined whether carriage of specific MHC alleles is associated with the dynamics of the HCV viral decline during the first 28 days of peginterferon therapy for chronic HCV and whether host genetic factors are associated with differences in viral dynamics by race.
|
|
|---|
MHC genotyping. Genotyping of class I MHC A, B, and C loci and the class II loci, DRB1, DQA1, DQB1, DPB1, and DPA1, were conducted by using medium resolution HLA typing methods with biotin-labeled PCR primers and hybridization of the labeled amplicons to sequence-specific oligonucleotide probes immobilized on nylon membranes developed by Roche Molecular Systems. All typing of MHC genes were conducted by Roche Molecular Systems (Alameda, CA).
HCV viral assessment. Viral levels were assessed for all participants in the study at baseline, as well as day 1, 2, 7, 14, and 28 after the initiation of treatment. Although more sensitive assays have been subsequently developed, the Roche HCV Amplicor version 2 assay, which has a lower detection limit of 600 copies/ml, was used in the present study.
Evaluation of population structure. Using data from 161 ancestry-informative single nucleotide polymorphisms, we derived estimates of individual admixture for participants in the genetics study and utilized the structured association method developed by Pritchard et al. (23, 24) to evaluate the population structure. We observed two distinct ancestral groups that had a strong correlation with self-reported race (30). Consequently, we used self-reported race in these analyses.
Statistical analysis.
Carriage of MHC alleles that occurred with a frequency of
5% in either racial group was assessed for the potential effect on the change in viral decline during at baseline and days 1, 2, 7, 14, and 28 after the start of treatment (4, 27).
All viral level measurements were log10 transformed to approximate a normal distribution. Due to the fact that the Roche Amplicor assay had a lower limit of detection of 600 IU/ml, we imputed the viral level values for those below detectable using a uniform distribution in the log scale. The time variable was expressed as "weeks" (rather than days or hours) in our analyses. A square-root transformation was applied to the time variable to attain linearity.
Longitudinal analysis was conducted using a multilevel model for change (also known as a linear random coefficients model or linear growth model) (4, 27). Detailed mathematical representations of the model are available online in Table S1 in the supplemental material, along with the variance components (see Table S2 in the supplemental material). Briefly, a two-level modeling strategy was utilized; a level 1 submodel describes how each person changes over time, and a level 2 submodel relates "interindividual differences in change" to predictors such as allele carriage and race. A multilevel modeling was chosen because the primary question of interest was to determine whether allele and/or race had an effect on the rate of change in viral levels over time (the dynamics of viral decline). We also adjusted for potential confounding factors such as patient age, sex, fibrosis score (Ishak > 3 [versus
3]), and weight, which we have previously reported to be associated with the interferon response (5). Since these variables did not significantly change our final qualitative conclusions, we omitted them from the final models for the simplicity and focused on the effects of allele and race only. By fitting each "patient" as a random effect, we accounted for the interindividual variability of viral levels at the baseline (i.e., to allow intercepts to vary between patients) "patient-by-time interaction" as a random effect to account for the interindividual variability in rates of change of viral levels (i.e., to allow slopes of the decline to vary between patients) during the course of the first 28 days of therapy. The variables "time," "allele," and "race" were included as fixed effects. To evaluate the differential effects of time on HCV RNA by race and allele (i.e., statistical interactions), we also included "race-by-time," "allele-by-time," "race-by-allele," and "allele-by-race-by-time interaction" as fixed effects.
The SAS v.9.1 system (Proc Mixed; SAS Institute, Cary, NC) was used for data analyses, and the R v.2.3.1 statistical package was used for graphical presentation of data (R Language and Environment for Statistical Computing [www.r-project.org]). Statistical significance was set at
= 0.05.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Baseline characteristics of the 373 participants with MHC genotyping data in the Virahep-C genetics study categorized by race
|
![]() View larger version (11K): [in a new window] |
FIG. 1. Average viral levels (based on imputed and observed viral levels) in the cohort overall before and after transformation. (A) Viral levels before transformation of the time variable (time is expressed as days); (B) viral levels after transformation of time (time is expressed as weeks) utilizing a square-root transformation. Dark circles indicate the mean viral levels (in log10 IU/ml) at each time point. Vertical bars depict the 95% confidence intervals (mean ± 1.96 [standard error of the mean]). Solid lines represent the estimated linear relationship based on general linear models. A slightly better linear relationship can be seen between mean viral level and time in panel B than in panel A.
|
MHC allele carriage and the dynamics of viral decline.
We observed a total of 47 alleles in our cohort with a frequency of
5%. These alleles, along with their corresponding frequencies, are presented in Table 2. Of the alleles tested, four were significantly associated with the dynamics of the viral decline (P < 0.05). The models and results presented in Table 3 evaluate the effects of allele carriage on the initial baseline viral level, as well as on the rates of change in viral levels, over time, controlling for the effects of race.
|
View this table: [in a new window] |
TABLE 2. Distribution of MHC class I and class II alleles in the Virahep-C cohorta
|
|
View this table: [in a new window] |
TABLE 3. Summary of multilevel modelsa
|
However, we did observe differences in the dynamics of the viral decline (Table 3). For DQA1*04 allele carriage, the rate of viral decline was higher among noncarriers than carriers (P < 0.001) and the magnitude of the rate of decline depended on race (Table 3 and Fig. 2A). The estimated mean viral level declines for the first 4 weeks were 2.66 for CA noncarriers, 1.51 for AA noncarriers, 1.19 for AA carriers, and 0.98 for CA carriers (log10 IU/ml P = 0.018, Fig. 2A), suggesting that CA noncarriers exhibited the fastest decline. This trend was similar when we examined the estimated mean viral level decline during the first week (Fig. 2A).
![]() View larger version (22K): [in a new window] |
FIG. 2. The predicted mean decline of HCV viral levels by allele carriage and race. The vertical axis presents log10-transformed predicted mean viral levels from the random coefficients model, while the horizontal axis presents time. The time variable was expressed as days. In addition, predicted baseline viral levels along with the rates of viral change, are provided below each figure by race and the carriage of each of the alleles over time (with a week transformation). (A) DQA1*04; (B) DQB1*0402; (C) A*03; (D) Cw*03.
|
CA noncarriers of A*03 had the highest rates of viral decline compared to CA carriers; the estimated mean viral level declines for the first 4 weeks were 2.10 for CA carriers, 2.75 for CA noncarriers, 1.65 for AA carriers, and 1.39 for AA noncarriers (log10 IU/ml, P = 0.018) (Fig. 2C).
CA carriers of Cw*03 had the largest rates of viral decline compared to CA noncarriers, while AA noncarriers had higher rates of viral decline than did the AA carriers. The estimated mean viral level declines for the first 4 weeks were 2.99 for CA carriers, 2.49 for CA noncarriers, 1.51 for AA noncarriers, and 1.19 for AA carriers (log10 IU/ml, Fig. 2D). For all four of these alleles, the viral subgenotype (1a or 1b) did not have an impact on outcome (P > 0.05).
|
|
|---|
We also observed that the carriage of DQA1*04 or DQB1*0402 alleles was associated with slower rates of viral decline. Interestingly, the frequencies of these alleles are greater among AAs than among CAs (Table 2). In addition, carriage of each of these alleles is associated with a greater difference in the level of viral decline among CAs than among AAs. The frequency of DQA1*04 is 10.1% among AAs and 2.8% among CAs; the frequency of DQB1*0402 is 8.1% among AAs and 2.8% among CAs. Differential distribution of the frequencies of these alleles by race may contribute to the observed population differences in response to peginterferon by race on the population level.
The DQA1*04 and DQB1*04 alleles have been previously implicated as a risk for Crohn's disease (31) and cirrhosis (10) in chronic HCV infections. The DQB1*0401/*0402 alleles are associated with viral persistence in several Japanese studies (1, 10, 15, 23). Previous studies that have examined the structure of DQB1*04 have shown that a leucine substitution exists at position 56 (21), which is located in the antigen-binding groove of the DQ molecule and probably contributes to the ability to present HCV peptides. The leucine substitution is present only in the DQB1*0401 and DQB1*0402 alleles. All other DQB1 alleles have a proline at this position. Future studies are needed to specifically examine HCV antigen presentation in the context of these two class II MHC alleles.
In addition, the DQB1*0402 and DQA1*04 alleles exist as a haplotype and, therefore, our findings of similar magnitude and direction of effects for the individual alleles are not surprising. Because the DQA1*04 and DQB1*0402 alleles often exist on a haplotype with DRB1*03 or DRB1*08, we examined whether carriage of the DRB1*08-DQA1*04-DQB1*0402 or DRB1*03-DQA1*04-DQB1*0402 haplotypes was associated with viral dynamics. However, due to the relatively low frequency of individuals with DRB1*08-DQA1*04-DQB1*0402 (4 AA and 10 CA patients) and DRB1*03-DQA1*04-DQB1*0402 (24 AA and 0 CA patients) in the cohort overall and the extreme distributions within each race, we were not able to obtain meaningful results with our longitudinal analyses (data not shown) and, therefore, are not able to exclude the possibility that another allele in linkage disequilibrium with the DQA1*04 or DQB1*0402 may negatively affect viral decline during therapy. Extended haplotypes involving DQB1*0401/*0402 have also been associated with greater HCV-induced liver injury (15), although the immediate functional relevance with respect to viral dynamics is not clear.
CA noncarriers of the class I allele A*03 had greater rates of viral decline than CA carriers. In contrast, AA noncarriers had slower rates of viral decline than AA carriers (Fig. 2C), unlike the effect within CAs. The effects of allele carriage appear to be greater among CA noncarriers and carriers than among AA noncarriers and carriers. Overall, the frequency of the A*03 allele is slightly lower among AAs than CAs (11.2% among AAs and 14.2% among CAs). It is possible that the main effects of this allele are seen among CAs, and that the smaller differences in viral dynamics among AA noncarriers and carriers is due to other undetermined genetic and/or nongenetic factors. It is also possible that AA and CA haplotype differences may contribute to the observed differences. Alternatively, it is possible that our observations are due to chance. Larger studies are needed to evaluate the potential effects of extended haplotypes within each race.
Carriers of Cw*03 had more pronounced viral dynamic declines than CA noncarriers (Fig. 2D). Viral declines among AA carriers and noncarriers were similar. Class I alleles are important in CD8+ T cells responses, which have been shown to play a role in the clearance of HCV during therapy (28). Future studies are needed to understand the role of Cw*03 in HCV antigen presentation. Overall, the frequency of the Cw*03 allele is similar among AAs and CAs (10.9% among AAs and 10.3% among CAs).
In summary, our observations suggest that carriage of DQA1*04 and DQB1*0402, or another allele in linkage disequilibrium with these alleles, may represent a risk for slower rates of viral decline during peginterferon therapy. With both DQA1*04 and DQB1*0402, CA noncarriers had better rates of viral decline during the course of therapy compared to CA carriers, AA noncarriers, and AA carriers, respectively (Fig. 2C and D).
The 28-day viral response represents a critical biological point for the study of the response to therapy due to its prognostic potential (16, 17, 19). Unlike SVR, early viral responses are less likely to be affected by dose adjustments and patient compliance and may therefore represent a more purely biologic picture of the role of host genetics in the interferon response. It is likely that in a complex process such as the response to peginterferon treatment for chronic HCV infection, multiple genes are involved, with different genes possibly operating at different time points during the course of therapy. Therefore, it is possible that the alleles that we describe here play an important role, either independently, as a haplotype, or in combination with other unidentified genetic markers during the early period after the start of therapy, while other markers may exert a more influential role at subsequent time points during the course of therapy. In the future, larger studies are needed to refine the potential influence of extended MHC haplotypes on HCV viral dynamics during treatment.
Previous mathematical models using differential equations have suggested that the shape of the viral decline with standard alpha interferon or peginterferon is nonlinear. For example, biphasic and triphasic declines have been reported (9, 22). However, these models are based on very frequent sampling of viral levels spaced only hours apart, as opposed to days apart as with our study, meaning that the frequency and spacing of viral level measurements are critical to the modeling of viral kinetics. In the absence of intensive sampling of viral levels, we used a square root transformation of the time scale to linearize the trajectory of the viral decline for our statistical modeling. This approach allowed us to focus our analysis on the linear relationship of MHC allele carriage on the overall dynamics of the viral decline, rather than the smaller fluctuations in the viral decline after each administration of interferon. Consequently, these linear individual change trajectories were modeled, and the effects of allele and race on the rate of change were examined. This analysis allowed a more simplified summary measure of the relationship between MHC allele carriage and the trajectories of viral levels over time.
A key factor in genetic association studies is the potential for confounding from population stratification. We have previously utilized recently developed techniques in mathematical genetics to infer individual admixture and found self-reported race to be highly correlated with individual admixture. Consequently, we have utilized self-reported race in our study (30). In addition, due to the hypothesis-generating nature of the present study, we have deliberately elected to not adjust for the number of comparisons made rather than reject an association solely based on statistical grounds (29). Although it is possible that the associations that we have observed are false positives, the consistent patterns within each race provide additional evidence of potential biological effects (Fig. 2).
Our study offers an important first picture of the role of diversity in the human MHC on the dynamics of viral decline during the first 28 days of therapy. Future studies are needed to confirm our observations, as well as expand our understanding of the biological mechanisms of these alleles in the context of HCV therapy. We have previously reported that the A*02, B*58, and DPB1*1701 alleles are associated with the SVR to alpha peginterferon-plus-ribavirin therapy (25). The response to therapy is likely a complex process that occurs over the course of time. Most likely, different combinations of genes are involved at different time points to clear viremia. For example, studies of Caucasian populations have consistently shown that the class II alleles DRB1*11 and DQB1*0301 are associated with natural clearance of HCV viremia (11, 29). However, a clear understanding of the biological mechanisms of these alleles is still lacking. Additional research, particularly with understanding of the differences in the presentation of HCV peptides by different HLA alleles, is needed to better understand the differences that we have observed between early viral dynamics and SVR. Additional nongenetic factors, such as decreased immune function through coinfections, may also play a role. Further research is needed to help evaluate these factors, as well as the functional characteristics of these alleles on HCV elimination.
The members of Virahep-C included Nezam Afdhal and Tiffany Geahigan from the Beth Israel Deaconess Medical Center, Boston, MA; Robert S. Brown, Jr., Lorna Dove, and Shana Stovel from the New York-Presbyterian Medical Center, New York, NY; Norah Terrault, Stephanie Straley, Eliana Agudelo, and Melissa Hinds from the University of California, San Francisco, CA; Thelma E. Wiley and Monique Williams from Rush University, Chicago, IL; Charles D. Howell and Karen Callison from the University of Maryland, Baltimore, MD; Lennox J. Jeffers, Shvawn McPherson Baker, Maria DeMedina, and Carol Hermitt from the University of Miami, Miami; Hari S. Conjeevaram, Robert J. Fontana, and Donna Harsh from the University of Michigan, Ann Arbor, MI; Michael W. Fried, Scott R. Smith, Dickens Theodore, Steven Zacks, Roshan Shrestha, Karen Dougherty, Paris Davis, and Shirley Brown from the University of North Carolina, Chapel Hill, NC; John E. Tavis, Adrian Di Bisceglie, Ermei Yao, Maureen Donlin, Nathan Cannon, and Ping Wang from St. Louis University, St. Louis, MO; Huiying Yang, George Tang, and Dai Wang from Cedars-Sinai Medical Center, Los Angeles, CA; Hugo R. Rosen, James R. Burton, and Jared Klarquist from the University of Colorado Health Sciences Center, Denver; Scott Weston from the Veteran's Administration, Portland, OR; Milton W. Taylor, Corneliu Sanda, Joel Schaley, and Mary Ferris from Indiana University, BloomingtonN; Steven H. Belle, Richard A. Bilonick, Geoffrey Block, Jennifer Cline, Marika Haritos, KyungAh Im, Stephanie Kelley, Sherry Kelsey, Laurie Koozer, Sharon Lawlor, Stephen B. Thomas, Abdus Wahed, Yuling Wei, and Leland J. Yee from the Data Coordinating Center, Graduate School of Public Health at the University of Pittsburgh, Pittsburgh; Patricia Robuck, James Everhart, Jay H. Hoofnagle, Edward Doo, T. Jake Liang, and Leonard B. Seeff from the National Institute of Diabetes and Digestive and Kidney Diseases; and David E. Kleiner from the National Cancer Institute, Bethesda, Maryland.
Published ahead of print on 13 October 2008. ![]()
Supplemental material for this article may be found at http://aac.asm.org/. ![]()
Present address: Department of Discovery, Medicine, and Epidemiology, AstraZeneca Pharmaceuticals LP, Wilmington, Delaware. ![]()
Contributing members of the Virahep-C Study are listed in the Acknowledgments. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»