Previous Article | Next Article ![]()
Journal of Virology, April 2007, p. 3391-3401, Vol. 81, No. 7
0022-538X/07/$08.00+0 doi:10.1128/JVI.02640-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Department of Biology, Indiana University, Bloomington, Indiana 47401,1 Department of Biochemistry and Molecular Biology and Center for Medical Genomics, Indiana University School of Medicine, Indianapolis, Indiana,2 Institute of Evolution, University of Haifa, Haifa, Israel,3 Department of Molecular Microbiology and Immunology, Saint Louis University School of Medicine, St. Louis, Missouri,4 Department of Medicine, University of Maryland, Baltimore, Maryland,5 Graduate School of Public Health, University of Pittsburgh, Pittsburgh, Pennsylvania6
Received 29 November 2006/ Accepted 19 January 2007
|
|
|---|
|
|
|---|
Alpha interferon mediates its antiviral and pharmacological effects by binding to type I interferon receptors on the cell surface membrane, which leads to transcription of up to 1,000 interferon-stimulated genes, presumably via the Janus-activated kinase 1 (JAK1)-STAT (signaling transducers of activation and transcription) signaling pathway (28, 30). A potential explanation for a lack of response to interferon therapy of HCV infection is an underlying deficient cellular response to interferon with a blunted response to interferon signaling, this being more common among AA patients than CA patients. To test this hypothesis, global gene expression in peripheral blood mononuclear cells (PBMC) before and during the first 28 days of therapy with peginterferon and ribavirin was analyzed in a cohort of AA and CA patients with genotype 1 HCV infection. These patients were undergoing therapy in the Study of Viral Resistance to Antiviral Therapy for Hepatitis C (Virahep-C), a large, prospective, multicenter study designed to define the differences in response rates among AA and CA patients and to determine clinical, immunological, host genetic, viral genetic, and interferon cell signaling factors that were associated with lack of response to treatment (5). The current analysis summarizes results of global gene expression in PBMC during the first 28 days of therapy, comparing patients with a marked (decrease of more than 3.5 log10 at 28 days following treatment initiation), intermediate (1.4 log10 to 3.5 log10 decrease), or poor (<1.4 log10 decrease) viral response.
|
|
|---|
From the Virahep-C cohort, 72 patients who did not have dose reductions of either peginterferon or ribavirin in the first 28 days of treatment were selected such that 12 patients of each race (CA and AA) were included by virological-response category. The three categories of response were marked, defined as a decrease in HCV RNA levels of more than 3.5 log10 IU/ml or to an undetectable level on day 28; intermediate, defined as a decrease of 1.4 to 3.5 log10 IU/ml on day 28; and poor, defined as less than a 1.4 log10 IU/ml decline on day 28 relative to baseline. These definitions were made a priori in an attempt to analyze the biological basis for virological responses. Of these, RNA adequate to provide gene expression information was not obtained from three patients.
PBMC preparation. PBMC were collected in sodium-heparin cell preparation tubes at day 0 (before treatment) and days 1, 2, 7, 14, and 28 after initiation of treatment. Whole blood was diluted with an equal volume (8 ml) of phosphate-buffered saline (PBS), carefully layered over a 10-ml Ficoll-Hypaque gradient (Amersham/Pharmacia, Piscataway, NJ), and centrifuged at 800 rpm for 20 min at room temperature. The buffy coat layer was transferred to a 15-ml RNase-free tube, diluted with PBS, and centrifuged at 100 x g for 15 min at room temperature. The supernatants were discarded, and the PBMC were retained.
RNA extraction. Samples were shipped overnight by express courier at 4°C to a central repository, where RNA was isolated on arrival. The PBMC were lysed in 1 ml of TRI reagent (Molecular Research Center Inc., Cincinnati, OH). The PBMC lysate was mixed with 1-bromo-3-chloropropane phase separation agent for 1 min and incubated at room temperature for 15 min. After centrifugation for 15 min at 12,000 rpm and 4°C, RNA was precipitated from the supernatant overnight at 20°C with an equal volume of isopropanol and 1/10 volume of 7.5 M ammonium acetate. The precipitate was washed twice with 75% ethanol and then with 95% ethanol. RNA was briefly air dried and then further purified using RNeasy columns (QIAGEN, Valencia, CA). The amount and quality of RNA were determined by spectrophotometry and by electrophoresis through 1% agarose with ethidium bromide, and RNA quality was analyzed by the Agilent Bioanalyzer according to the manufacturer's instructions. Samples that did not show two clear bands of rRNA were discarded.
RNA labeling and hybridization. Preparation of cDNA and cRNA and labeling were carried out according to the protocols recommended by Affymetrix (Santa Clara, CA) in the GeneChip expression analysis technical manual, as previously described (34).
Array analysis and data processing. The microarrays were scanned using a dedicated model 3000 scanner controlled by GCOS software. The average intensity on each array was normalized by global scaling to a target intensity of 1,000. Data were extracted using the Affymetrix Microarray Suite 5 (MAS5) algorithm and exported into a custom-designed database (MicroArray Data Portal) in the Center for Medical Genomics (Indiana University-Purdue University Indianapolis, Indianapolis). All DNA microarray chips were analyzed for unequal distribution or artifacts as described previously (4). Any chip shown to be defective was corrected or dropped from the analysis.
The MicroArray Data Portal, in addition to its role as a database and analytical tool, is an informatics platform with active links from each sequence to several public databases. Sequence information for each gene on the HG-U133A GeneChip was obtained by parsing the HG-U133 target file obtained from the Affymetrix informatics website, http://www.affymetrix.com/analysis/download_center.affx). A GenBank accession number and a Unigene cluster were used to match sequences to their corresponding LocusLink number, gene symbol, and map position and to link with Gene Ontology (GO) terms and Enzyme Nomenclature (19) EC numbers. EC numbers were then used in conjunction with the Ligand database to link genes to KEGG pathways (www.genome.ad.jp/kegg/).
Statistical analysis.
The MAS5 data were filtered to eliminate any gene that was not called present in at least 50% of the samples in any one group (fraction present
0.5) (24). Changes (n-fold) for each gene were calculated using the ratio of the MAS5 signals of the baseline and the posttreatment time. If the signal for the posttreatment time point was greater than the baseline the change was calculated as +averageposttreatment/averagebaseline; otherwise, the change was calculated as averagebaseline/averageposttreatment. The asymptotic standard errors (ASE) were estimated using the delta method, and 95% confidence intervals were calculated by multiplying the ASE by 1.96, with the product added to and subtracted from the change.
Welch's t test using the MAS5 signals was used to test for differences in gene expression between CA and AA, and one-way analysis of variance was used to test for differences among the three response groups.
For each gene, the expression levels of posttreatment time points were compared to the baseline (pretreatment) expression levels using a paired t test of the MAS5 signals. Genes whose P value
0.001 and for which the absolute value of the change was at least 1.5-fold were selected as significant. Because of the filtering and differences in power, the numbers of genes considered to be significant in mutually exclusive and exhaustive subgroups will not necessarily add up to the number of genes considered to be significant in the entire sample. Genes that are significant in both racial groups contribute twice to the sum of genes but only once to the number of significant genes in the entire sample. On the other hand, gene expression differences that meet the change criterion may not meet the criterion of a P value <0.001 in either racial group but, due to the increased number of observations for the entire sample compared to each racial group, do have a P of <0.001 for the entire sample; such genes do not contribute to the sum of significant genes across racial groups but do contribute to the number of significant genes in the entire sample.
All analyses were performed using the R statistical language and environment.
Microarray data accession number. Microarray data presented in this paper have been deposited with NCBI/GEO under accession no. GSE7123.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Baseline participant characteristics by response group
|
|
View this table: [in a new window] |
TABLE 2. Baseline patient characteristics by race
|
0.01 and 2 at a P of
0.001, whereas 109 and 11, respectively, would be expected by chance at these significance levels.
Global gene expression response is greater in marked responders than in poor responders.
Gene expression in PBMC changed substantially during peginterferon and ribavirin therapy, with major changes being evident by days 1 and 2 after the initial injection of peginterferon and administration of ribavirin. The numbers of genes that were significantly modified (absolute value of change greater than 1.5-fold and P
0.001) at each time point for each response group and for each racial group within the response group are shown in Table 3. Many genes were altered in expression at the early time points in patients in all three response groups. The number of differentially expressed genes dropped between day 2 and day 7 and increased again slightly between days 7 and 28. Postbaseline PBMC samples were generally taken before administration of interferon. For only one subject at two time points, samples were taken 4 h after administration of interferon; this did not appear to affect the results for this patient.
|
View this table: [in a new window] |
TABLE 3. Number of genes (proportion present > 0.5) modified (at least 1.5-fold change; P 0.001) during peginterferon and ribavirin treatment
|
Within marked responders, there were more genes changed in AA than in CA at every time point (Table 3). Among poor responders, the same relationship held except at day 28. The relationship was more mixed in intermediate responders, where more genes changed expression in AA than CA in the first 2 days after treatment but more genes changed expression in CA than in AA after that.
The number of genes that changed in expression was greater in the marked responders than in the intermediate or poor responders at all time points. There was a smaller difference in the numbers of genes that changed in expression between the intermediate and poor responders at most points. Figures 1 and 2 show the numbers of genes whose expressions increased or decreased, respectively, using the change filter of 1.5-fold and a P value of
0.001. The numbers of genes that were increased in gene expression were far higher in the marked responders than in either the intermediate or poor responders, and the intermediate responders had numbers intermediate between the other two groups. Slightly more genes were up- than down-regulated. Although the differences among the three response categories of patients held for both up- and down-regulation of genes, the decline over time in numbers of down-regulated genes was much sharper than the decline in up-regulated genes.
![]() View larger version (26K): [in a new window] |
FIG. 1. Number of genes up-regulated (P < 0.001; 1.5-fold change) at each time point compared to baseline in each response category of patient.
|
![]() View larger version (24K): [in a new window] |
FIG. 2. Number of genes down-regulated (P < 0.001; 1.5-fold change) at each time point compared to baseline in each category of patients.
|
![]() View larger version (47K): [in a new window] |
FIG. 3. Increase in mRNA as detected by microarrays for oligo(A) synthetase 1 and 2 and MX1 and MX2 in all three response categories of patients at days 1, 2, 7, 14, and 28 after initiation of treatment.
|
|
View this table: [in a new window] |
TABLE 4. Change in gene expressiona
|
|
View this table: [in a new window] |
TABLE 5. Changes on day 1 for CA and AA and the combined groupa
|
|
|
|---|
In the Virahep-C study, AA had lower rates of sustained response and higher levels of serum HCV RNA than CA at almost all time points (6). The current analyses show that differences in the number of genes whose expression changed at least 1.5-fold between marked and poor responders in both AA and CA patients occurred within 24 h of starting treatment (Tables 4 and 5; Fig. 1 and 2). The level of gene expression and number of genes induced or down-regulated were considerably higher among marked virological responders than among poor responders (see Tables S1 and S2 in the supplemental material). However, it is puzzling that AA had higher numbers of genes induced and slightly higher levels of gene expression at day 1 than others.
Alpha interferon is known to act through induction of a large number of genes, the exact number and pattern of which have only been partially identified (33, 34). In this study, 801 genes were found to be increased during peginterferon and ribavirin therapy; many, but not all, of these were known interferon-induced genes (see Tables S1 and S2 in the supplemental material). While global interferon-induced gene expression was less among poor responders than marked responders, no specific gene could be linked to the differences in responses or to racial differences. Thus, poor or nonresponse appeared to be a global blunting of interferon cell signaling, rather than the lack of induction or function of a specific antiviral gene product. These data are in agreement with recent findings in which the gene expression of nonresponders was lower than that of responders when PBMC were cultured from such patients (16). Despite differences in the microarray systems used and assessment of in vitro versus in vivo responses, the levels of induction (n-fold) for many genes were remarkably similar. However, He et al. (16) found that the levels of gene induction in white patients was higher than in black patients. However, we could not find any difference between the racial groups in this study in levels of gene expression. The present study demonstrates that, controlling for virological response, gene expression changes were actually more common, within 2 weeks of treatment initiation, in AA than CA patients who had a marked or poor response. The reason for this difference is not known.
Previous studies using cell culture systems have suggested that HCV replication or presence of HCV antigens may interfere with specific interferon-induced gene products, such as the well characterized antiviral enzymes OAS, protein kinase R, and adenosine deaminase (8, 10, 12). The present analysis, in contrast, suggests that lack of response to administered interferon was due to an ongoing physiological defect that causes blunted regulation of interferon responsiveness. The blunted response might be due to a prior inflammatory response, interferon receptor deficiency or dysfunction, or lack of afferent cell signaling through the JAK-STAT pathway. In this regard, several recent studies in vitro and in vivo have suggested that a deficiency in STAT1 activation or DNA binding occurs in patients with chronic HCV infection (22). Such findings are compatible with the findings in this study. In fact genes such as the IRF-7 gene (Tables 4 and 5), a key gene in induction of interferon, was induced compared to baseline at lower levels in poor responders than in marked-response patients, as was the cig 5 (viperin) gene, previously identified as being important in the interferon response to hepatitis C virus (17). Toll-like receptor 7 (TLR7) has been shown to be important in the recognition of single-stranded viral RNA and subsequent signaling of the interferon, I
B kinase
/ß/
, and mitogen-activated protein kinase cascades leading to NF-
B and AP-1 activation and to IRF-7 and interferon production (13). Expression of the IRF-7 gene was increased from baseline to levels almost twice as high in marked-response than in poor-response patients and was thus strongly induced by peginterferon/ribavirin combination therapy.
Several limitations of the present findings deserve mention. First and foremost, the analysis of gene expression was conducted on PBMC and not on hepatocytes that harbor replicating HCV. Analysis of hepatocytes, however, requires liver biopsy, an invasive procedure which cannot be done repeatedly in humans during interferon therapy. Furthermore, analyses on liver tissue include PBMC and other nonparenchymal cells, and changes in expression in liver tissue may not reflect effects on hepatocytes only. Responses in PBMC are more likely to reflect a global response and not be under the local control of replicating virus or disease activity, which may modulate interferon responses. The chimpanzee model of HCV infection offers a potential approach to analyzing intrahepatic gene expression during interferon therapy (3, 21, 32). However, chimpanzees respond minimally to human alpha interferon therapy, and interpretation of results has to take into consideration interspecies differences.
A final limitation to this study was that it was based upon viral kinetic analyses done during the first 28 days of therapy and was not based on results of sustained virological responses. This design was purposeful, in that early virological responses are highly predictive of ultimate responses and are not affected by nonbiologic factors, such as dose modification, compliance, and dropout. Only patients who took the full prescribed dose of peginterferon were selected. Furthermore, the differences between responders and nonresponders in the strength of gene induction were found even at day 1, which occurred after an observed administration of peginterferon and ribavirin at the initiation of treatment. Thus, by using early viral responses, purely biological factors associated with response and nonresponse could be assessed.
In this study, a poor virological response to peginterferon and ribavirin therapy of HCV infection was found to be associated with global, blunted changes in interferon-responsive gene expression. These results indicate that the blunted response is not specific to the liver or to virally infected cells. This hyporesponsiveness may be determined by host genetics, or it may be due to an environmentally induced lesser sensitivity to interferon. It is also possible that PBMC are exposed to viral proteins in circulation or to hepatocyte-associated HCV antigens which might alter the immune response of such cells to interferon treatment.
Members of Virahep-C contributing to the study include, from the Beth Israel Deaconess Medical Center, Boston, MA, Nezam Afdhal (principal investigator) and Tiffany Geahigan (research coordinator); from the New York-Presbyterian Medical Center, New York, NY, Robert S. Brown, Jr. (principal investigator), Lorna Dove (coinvestigator), Shana Stovel (study coordinator), and Maria Martin (study coordinator); from the University of California, San Francisco, San Francisco, Norah Terrault, (principal investigator), Stephanie Straley, Eliana Agudelo, Melissa Hinds (clinical research coordinator), and Jake Heberlein (clinical research coordinator); from Rush University, Chicago, IL, Thelma E. Wiley (principal investigator) and Monique Williams (study coordinator); from the University of Maryland, Baltimore, Charles D. Howell (principal investigator), Kelly Gibson (project coordinator), Karen Callison (study coordinator), and Jane Lewis (study coordinator); from the University of Miami, Miami, FL, Lennox J. Jeffers (principal investigator), Shvawn McPherson Baker (coinvestigator), Maria DeMedina (project manager), and Carol Hermitt (project coordinator); from the University of Michigan, Ann Arbor, Hari S. Conjeevaram (principal investigator), Robert J. Fontana (coinvestigator), and Donna Harsh (study coordinator); from the University of North Carolina, Chapel Hill, Michael W. Fried (principal investigator [K24 DK066144]), Scott R. Smith (coinvestigator), Dickens Theodore (coinvestigator), Steven Zacks (coinvestigator), Roshan Shrestha (coinvestigator), Karen Dougherty (coinvestigator), Paris Davis (study coordinator), and Shirley Brown (study coordinator); from St. Louis University, St. Louis, MO, John E. Tavis (principal investigator), Adrian Di Bisceglie (coinvestigator), Ermei Yao (coinvestigator), Maureen Donlin (coinvestigator), Nathan Cannon (graduate student), and Ping Wang (lab technician); from Cedars-Sinai Medical Center, Los Angeles, CA, Huiying Yang (principal investigator), George Tang (project scientist), and Dai Wang (project scientist); from the University of Colorado Health Sciences Center, Denver, Hugo R. Rosen (principal investigator), James R. Burton (coinvestigator), and Jared Klarquist (lab technician); from Veteran's Administration, Portland, OR, Scott Weston (lab technician); from Indiana University, Bloomington, Milton W. Taylor (principal investigator), Corneliu Sanda (postdoctoral associate), Takuma Tsukahara (statistician), and Mary Ferris (lab assistant); from the Data Coordinating Center, Graduate School of Public Health at the University of Pittsburgh, Pittsburgh, PA, Steven H. Belle (principal investigator), Richard A. Bilonick (statistician), Geoffrey Block (coinvestigator), Jennifer Cline (data manager), Marika Haritos (statistician), KyungAh Im (statistician), Stephanie Kelley (data manager), Sherry Kelsey (coinvestigator), Laurie Koozer (project coordinator), Sharon Lawlor (data coordinator), Stephen B. Thomas (coinvestigator), Abdus Wahed (statistician), Yuling Wei (project coordinator), Leland J. Yee (consultant), and Song Zhang (statistician); from the National Institute of Diabetes and Digestive and Kidney Diseases, Patricia Robuck (project scientist), James Everhart (scientific advisor), Jay H. Hoofnagle (scientific advisor), Edward Doo (scientific advisor), T. Jake Liang (scientific advisor), and Leonard B. Seeff (scientific advisor); and from the National Cancer Institute, David E. Kleiner (central pathologist).
We thank Mary Ferris for the excellent record keeping and entering of data into the portal at the Center for Medical Genetics. We thank Ron Jerome and Chunxiao Zhu for expert assistance with the microarray studies, which were carried out using the facilities of the Center for Medical Genomics at Indiana University School of Medicine. We also thank Song Zhang and Jia Li from the data coordinating center, Pittsburgh, for statistical support and Jay H. Hoofnagle for help in editing the manuscript.
Published ahead of print on 31 January 2007. ![]()
Supplemental material for this article may be found at http://jvi.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»