Journal of Translational Medicine The prediction of interferon treatment effects based on time series microarray gene expression profiles

Tao Huang firstname.lastname@example.com
Shanghai Institutes for Biological Sciences
Key Laboratory of Systems Biology
Chinese Academy of Sciences
200031ShanghaiPR China

Shanghai Center for Bioinformation Technology
200235ShanghaiPR China

Kang Tu
Shanghai Institutes for Biological Sciences
Key Laboratory of Systems Biology
Chinese Academy of Sciences
200031ShanghaiPR China

Shanghai Center for Bioinformation Technology
200235ShanghaiPR China

Yu Shyr firstname.lastname@example.com
Cancer Biostatistics Center
Vanderbilt University
37232NashvilleUSA

Chao-Chun Wei firstname.lastname@example.com
Shanghai Center for Bioinformation Technology
200235ShanghaiPR China

Lu Xie email@example.com
Shanghai Center for Bioinformation Technology
200235ShanghaiPR China

Yi-Xue Li firstname.lastname@example.com
Shanghai Institutes for Biological Sciences
Key Laboratory of Systems Biology
Chinese Academy of Sciences
200031ShanghaiPR China

Shanghai Center for Bioinformation Technology
200235ShanghaiPR China

Journal of Translational Medicine The prediction of interferon treatment effects based on time series microarray gene expression profiles
10.1186/1479-5876-6-44BioMed Central Page 1 of 9 (page number not for citation purposes) Methodology * Corresponding authors †Equal contributors
Background:The status of a disease can be reflected by specific transcriptional profiles resulting from the induction or repression activity of a number of genes. Here, we proposed a timedependent diagnostic model to predict the treatment effects of interferon and ribavirin to HCV infected patients by using time series microarray gene expression profiles of a published study.Methods:In the published study, 33 African-American (AA) and 36 Caucasian American (CA) patients with chronic HCV genotype 1 infection received pegylated interferon and ribavirin therapy for 28 days. HG-U133A GeneChip containing 22283 probes was used to analyze the global gene expression in peripheral blood mononuclear cells (PBMC) of all the patients on day 0 (pretreatment), 1, 2, 7, 14, and 28. According to the decrease of HCV RNA levels on day 28, two categories of responses were defined: good and poor. A voting method based on Student's t test, Wilcoxon test, empirical Bayes test and significance analysis of microarray was used to identify differentially expressed genes. A time-dependent diagnostic model based on C4.5 decision tree was constructed to predict the treatment outcome. This model not only utilized the gene expression profiles before the treatment, but also during the treatment. Leave-one-out cross validation was used to evaluate the performance of the model.Results:The model could correctly predict all Caucasian American patients' treatment effects at very early time point. The prediction accuracy of African-American patients achieved 85.7%. In addition, thirty potential biomarkers which may play important roles in response to interferon and ribavirin were identified.Conclusion:Our method provides a way of using time series gene expression profiling to predict the treatment effect of pegylated interferon and ribavirin therapy on HCV infected patients. Similar experimental and bioinformatical strategies may be used to improve treatment decisions for other chronic diseases.

Background

Chronic diseases such as infectious disease, cancer, and diabetes are among the most common and costly health problems. The therapy of chronic diseases often lasts for a long time, while the treatment effect may be questionable and yet the side effects may be serious. Hepatitis C virus (HCV) is one of the major causes of chronic hepatitis, cirrhosis, and hepatocellular carcinoma. The current recommended treatment for chronic HCV infection is the combination of pegylated alpha interferon (peginterferon) and the oral antiviral drug ribavirin given for 24 or 48 weeks, but the chance to induce a sustained response is only 54%-56%[1]. Using interferon and ribavirin for a long time may cause serious side effects, such as fever, chills, body aches, headaches, myeloid disorders [2] and neuropsychiatric symptoms [3]. The patients with poor response should better give up such treatment in the early stage. However the underlying mechanisms for different responses are not fully understood and it is hard to foresee treatment effects by conventional methods.

We analyzed a published time series microarray dataset of a virological research in which the effects of pegylated interferon and ribavirin on 33 African-American (AA) and 36 Caucasian American (CA) patients with chronic HCV infection were studied [4]. We established a diagnostic model to predict the outcome of pegylated interferon and ribavirin therapy using time series microarray gene expression profiles for AA and CA patients separately.

Although the focus here is on how HCV infected patients respond to pegylated interferon treatment, the model described is generally applicable to other chronic diseases undergoing long term treatment.

Methods

Original time-series microarray data applied in our study

The original time-series microarray data used in this work is from a study of Milton W. Taylor which was published on Journal of Virology last year [4], and publicly available at GEO http://www.ncbi.nlm.nih.gov/geo under accession number GSE7123. The initial data set consists of the gene expression profiles of 33 African-American and 36 Caucasian American patients with chronic HCV genotype 1 infection on day 0 (pretreatment), and 1, 2, 7, 14, and 28 of pegylated interferon and ribavirin therapy. HG-U133A GeneChip containing 22283 probes was used to analyze the global gene expression in peripheral blood mononuclear cells (PBMC) of the patients at each time point. For each patient the decrease of HCV RNA level was calculated by subtracting baseline level (before treatment) from the level on day 28. Good response was defined as a decrease of more than 1.4 log 10 IU/ml of HCV RNA level; and poor response was defined as less than 1.4 log 10 IU/ ml decline from the base level. Only patients with all the gene expression data of 6 time points were involved in our analysis, including 30 Caucasian Americans (CA) of whom 17 were good responders and 13 were poor, 28 African-Americans (AA) of whom 19 were good responders and 9 were poor.

Data preprocessing

First, we normalized the data of total 348 microarrays using quantile method and log 2 -trasformed them. Only probes that were present in at least 75% microarrays with log 2 intensities greater than 7 were kept for further analysis. This resulted in a subset of approximately 13620 probesets representing 9100 different genes.

Statistical analysis to identify differentially expressed genes

To avoid bias that may be created by single feature-selecting statistical method, we constructed a voting method based on several methods including Student's t test [5], Wilcoxon test [6], empirical Bayes test (eBayes) [7]and significance analysis of microarray(SAM) [8] to identify differentially expressed genes. Only genes that passed three out of the four methods were regarded as differentially expressed. The selection criterion was set at a defined P value for all four statistical tests.

Time-dependent diagnostic model

The classifier used in our program at each time point was C4.5-a decision tree classification method [9]. With leaveone-out cross-validation, the model was trained and tested at each time point. The framework is illustrated in Figure 1 and detailed as follows:

Train the model Each patient in the training set was regarded as an instance and the class label for him was the outcome of the treatment. At each time point, differentially expressed genes between good and poor response group were identified using the voting method described above as the marker probe sets of this time point. At the first time point (time point 0, before treatment), the features were that day's gene expression values of the marker probe sets at that time point; at the following time point during the treatment, the features were the combination of that day's gene expression values of marker probe sets at that time point and features of previous time points. For example, the features at day 1 are the expression values of differentially expressed genes at day 1 and the expression values of differentially expressed genes at day 0.

Every patient is assumed treatable until predicted as nontreatable with sufficient differentially expressed genes at that day. For each time point, if the number of differentially expressed genes was equal or greater than 5, the C4.5 classifier will be constructed at this time point; otherwise, The framework of time-dependent diagnostic model Figure 1 The framework of time-dependent diagnostic model. "nontreatable" means the patient was predicted to have a poor response and should be eliminated from the treatment, "treatable" means the patient was predicted to have a good response and should keep the treatment.

differentially expressed gene number at this time point will be set as null and no C4.5 classifier will be trained at this time point. This check step helps to avoid false negative decision.

Test the model At each time of leave-one-out cross-validation, we used the data of N-1 patients to build a model and then applied it on the one left patient to predict his treatment outcome. If a patient was predicted as treatable by every time point's classifier, this was a positive prediction. If the final outcome according to the HCV RNA level decline was good response for this patient, this was a true positive prediction; otherwise, it was a false positive prediction.

If a patient was predicted as nontreatable by one of the six classifiers (day 0, 1, 2, 7, 14 and 28), this was a negative prediction. That means this patient should be eliminated from the treatment and the workflow of this patient will stop at that time point. If the real outcome was poor response for this patient, this was a true negative prediction; otherwise, this was a false negative prediction.

The prediction accuracy Q of leave-one-out cross-validation was calculated as follows:

tp, tn, fp and fn stand for true positive, true negative, false positive and false negative, respectively. Detailed information about this model, including processed microarray data, R code and results, can be found in Additional file 1 and file 2.

Relevance and significance of candidate biomarkers

To assess the biological relevance of the identified candidate biomarkers which were important for CA response prediction, we used PubGene to find relationships between these candidate biomarkers and IFN (Interferon)/HCV (Hepatitis C viruses). PubGene http:// www.pubgene.org/ is a tool to carry out automated extraction of explicit and implicit biomedical knowledge from publicly available gene and text databases to create a geneto-gene co-citation network by automated analysis of titles and abstracts in MEDLINE records [10]. Moreover, GO and KEGG category enrichment analyses were applied to validate the functional and pathway relevance of selected candidate biomarkers.

Results

Time-dependent diagnostic model

The true power of time series microarray analysis does not come from the analysis of single time point, but rather, from the analysis of a series of time points to identify a biomarker chain. The main idea of our model is to fully utilize gene expression profiles before and during treatment to predict the final treatment outcome.

The time-dependent diagnostic results of all patients, AA patients and CA patients are shown in Figure 2. It illustrates that true negative CA patients were all correctly detected on day 1 and true negative AA patients were mostly detected on day 1 and 7. The leave-one-out crossvalidation results of all patients, AA patients and CA patients are given in Table 1. In the preliminary research, we found that the numbers of differentially expressed genes in all patients, AA patients and CA patients under the same cut-off P value were quite different. To balance the numbers of differentially expressed genes in different groups, P values given to identify differentially expressed genes in all patients, AA patients and CA patients were 0.00001, 0.001 and 0.0001, respectively.

Simplified Time-dependent diagnostic model

We have known that if only static gene expression profiles before treatment were used the prediction accuracy was rather low (data not shown). However from the above results, it occurred to us that the seemingly complicated models may actually be simplified to day 1 classifierdepending only on gene expression profiles of very early treatment time point. The leave-one-out cross-validation Q tp tn tp tn fp fn = + + + + × 100 accuracy based on day 1 classifier (including day 0 and day 1 gene expression profiles) of CA patients could achieve 100%, the same as the result using data of all the time points. With AA patients, the accuracy dropped some, but still much better than if only using pre-treatment gene expression profile. The leave-one-out cross-val-idation results of all patients, AA patients and CA patients on day 1 are given in Table 2.

Identification of candidate biomarkers of CA patients

As stated above, CA patients of HCV infection are more sensitive to the therapy of interferon and ribavirin, and

The leave-one-out cross-validation results of time-dependent diagnostic models based on all patients, AA patients and CA patients Figure 2 The leave-one-out cross-validation results of time-dependent diagnostic models based on all patients, AA patients and CA patients. On day one, most nontreatable patients could be detected correctly. "T" means true prediction; "F" means false prediction.

after one day treatment the outcome could be one hundred percent predicted. Using the feature selection methods described in Methods section, we identified 30 differentially expressed genes or probes on day one between 17 good response CA patients and 13 poor response CA patients as the candidate biomarkers relevant to interferon therapy response. The time series expression profiles of two representative genes TNFSF10 (tumor necrosis factor (ligand) superfamily, member 10) and KLRK1 (killer cell lectin-like receptor subfamily K, member 1) are shown in Figure 3. In Figure 3A, at most of the time, the difference of TNFSF10 expression level between Good response CA and Poor response CA is much greater than the difference The time series expression of two representative genes TNFSF10 (tumor necrosis factor (ligand) superfamily, member 10) and KLRK1 (killer cell lectin-like receptor subfamily K, member 1) Figure 3 The time series expression of two representative genes TNFSF10 (tumor necrosis factor (ligand) superfamily, member 10) and KLRK1 (killer cell lectin-like receptor subfamily K, member 1). The horizontal axis depicts days that the treatment has lasted, the vertical axis stands for expression mean values with error bar. Different colour indicates different groups of patients.

between Good response AA and Poor response AA. Moreover, the difference between Good response CA and Poor response CA increases with time, but the difference between Good response AA and Poor response AA decreases with time. In Figure 3B, the same phenomenon can be observed on day 1 and 7. Many other genes showed similar tendency at earlier time points. Maybe this partly explains why the treatment of CA patients can be more accurately predicted than AA patients. The dynamic expression graphics of all thirty genes are provided in Additional file 4.

To further evaluate whether these expression signatures are associated with therapeutic outcome (good or poor response), we conducted clustering of CA patients using differentially expressed genes on day 1 between CA groups of good and poor outcome ( Figure 4C), and compared it with the clustering result of all patients using differentially expressed genes on day 1 between all patients of good and poor outcome( Figure 4A), of AA patients using differentially expressed genes on day 1 between AA groups of good and poor outcome ( Figure 4B).

It can be seen that Figure 4C best clustered its patients. These thirty genes could classify the CA patients into good responders and poor responders very well. Therefore the simplified day 1 diagnostic model can clearly be applied to CA patients.

Relevance and significance of candidate biomarkers

The relationship between those candidate biomarkers and IFN (Interferon)/HCV (Hepatitis C viruses) were explored by using PubGene. The two literature networks are shown in Additional file 5, indicating six genes that have direct connections with both IFN and HCV. They are LCK (lymphocyte-specific protein tyrosine kinase), NT5C2 (5'nucleotidase, cytosolic II), KLRK1 (killer cell lectin-like receptor subfamily K, member 1), CD86 (CD86 antigen (CD28 antigen ligand 2, B7-2 antigen)), PARP1 (poly (ADP-ribose) polymerase family, member 1) and TNFSF10 (tumor necrosis factor (ligand) superfamily, member 10). It has been reported that IFN-alpha can significantly enhance CD86 expression on dendritic cells from chronic hepatitis C patients [11,12]. The ligation of tumor necrosis factor receptor (TNFR1) can initiate apoptosis or programmed cell death which is part of interferon (IFN)-mediated anti-viral action [13].

GO category enrichment analysis results (see Additional file 6) show that many of these candidate biomarkers are involved in immunity, such as T cell differentiation and positive regulation of T cell activation, which are concerned with antivirus.

The KEGG category enrichment analysis (see Additional file 6) illustrates that three candidate biomarkers (LCK, lymphocyte-specific protein tyrosine kinase; KLRK1, killer cell lectin-like receptor subfamily K, member 1; TNFSF10, tumor necrosis factor (ligand) superfamily, member 10) are components of the Natural killer cell mediated cytotoxicity pathway, which is important in antineoplastic, antivirus and immune regulation. It has been reported that impairment of natural killer cell activity is associated with chronic hepatitis C virus infection [14].

Discussion

It is evident that the prediction accuracy of CA patients is higher than AA patients. Constitutively different responses of black and white hepatitis C patients to pegylated interferon and ribavirin therapy had been reported [15] before our study. It shows in our study as, markers that may be important for predicting response did not change as remarkably in AA patients as in CA patients. Our diagnostic model did not perform well on AA patients who did not demonstrate sufficient differentially expressed genes at early time points.

Our results implied that if there is a sensitive change of gene expression profile at early treatment time points, the diagnostic model will be more sensitive and useful. This was confirmed by the simplified diagnostic model. With day 1 model CA patients with HCV infection who are more responsive to interferon and ribavirin therapy can already be predicted of their treatment outcome after 28 days. With less sensitive response the diagnostic model may need to be stretched to include profiles after longer time treatment, like in AA patients. Even with AA patients, profiles of early time points are sufficient for making reasonable predictions of outcome. This is good news for clinical and experimental workers. It means that the strategy of using early time-treatment gene expression profiling to predict outcome for potential long-term treatment is affordable and applicable.

Microarray gene expression analysis has been proved valuable in numerous applications including disease classification, diagnosis, survival analysis, choice of therapy etc, but rarely for more complex clinical problems such as the dynamic prediction of treatment effects we addressed in this paper. We tried several methods including traditional statistical techniques and the latest computer-intensive techniques to predict the final treatment effects based on the static gene expression profiles before treatment and the prediction results were unacceptable. Dynamic prediction chains using time series gene expression profiles have been proved to make more successful prediction model. There were two outcome prediction studies based on 70gene expression dataset generated by kinetic reverse-transcription PCR from 52 multiple sclerosis patients treated with rIFNβ [16,17]. They obtained good results but were limited to the 70 genes. Our model directly applied on large-scale microarray data, and may have found some novel biomarkers. Our results justify further biological studies to evaluate whether these candidate biomarkers could truly predict the effect of interferon and ribavirian therapy. Further investigations may shed light on the mechanisms of different responses between CA patients and AA patients of HCV infection to this kind of therapy.

Conclusion

Our time-dependent diagnostic model suggests a way of using time series gene expression profiling to predict the treatment effect of pegylated interferon and ribavirin therapy on HCV infected patients. Similar experimental and bioinformatics strategies may be used to improve treatment decisions for other chronic diseases. This may be an important strategy in future personalized medicine.

Unsupervised two-way hierarchical clustering based on expression profiles of differentially expressed genes of all patients on day 1(A), AA patients on day 1(B), and CA patients on day 1(C)Figure 4Unsupervised two-way hierarchical clustering based on expression profiles of differentially expressed genes of all patients on day 1(A), AA patients on day 1(B), and CA patients on day 1(C). P values given to identify differentially expressed genes in all patients, AA patients and CA patients were 0.00001, 0.001 and 0.0001, respectively. The bar indicates the response status of patients.

Table 1 :
1The Leave-one-out cross-validation results of all patients, AA patients and CA patients.All patients
AA patients
CA patients

Predicted
Good

Predicted
poor

Predicted
Good

Predicted
poor

Predicted
Good

Predicted
poor

2 × 2 Table
Actual
Good

28
8
Actual
Good

16
3
Actual
Good

17
0

Actual
poor

8
14
Actual
poor

1
8
Actual
poor

0
1 3

Accuracy
72.4%
85.7%
100%

* Good means good response, poor means poor response.

Table 2 :
2The leave-one-out cross-validation results of all patients, AA patients and CA patients on day 1.All patients
AA patients
CA patients

Predicted
Good

Predicted
poor

Predicted
Good

Predicted
poor

Predicted
Good

Predicted
Poor

2 × 2 Table
Actual
Good

28
8
Actual
Good

16
3
Actual
Good

17
0

Actual
poor

8
14
Actual
poor

4
5
Actual
poor

0
1 3

Accuracy
72.4%
75%
100%

* Good means good response, poor means poor response.

(page number not for citation purposes)
Journal of Translational Medicine 2008, 6:44 http://www.translational-medicine.com/content/6/1/44
Publish with Bio Med Central and every scientist can read your work free of charge"BioMed Central will be the most significant development for disseminating the results of biomedical researc h in our lifetime." Sir Paul Nurse, Cancer Research UK Your research papers will be: available free of charge to the entire biomedical community peer reviewed and published immediately upon acceptance cited in PubMed and archived on PubMed Central yours -you keep the copyright Submit your manuscript here: http://www.biomedcentral.com/info/publishing_adv.aspBioMedcentral
AcknowledgementsThe authors would like to thank Editor Francesco Marincola for his responsible attitude and encouraging comments to our manuscript. We also would like to thank Dr. Hongyue Dai for his insightful comment.Authors' contributionsTH and KT carried out the study. TH and LX wrote the manuscript. LX and YxL supervised the project. YS provided statistical suggestions, CcW contributed discussions. All authors read and approved the final manuscript.Additional materialAdditional file 1
. Translational Medicine. 6of Translational Medicine 2008, 6:44 http://www.translational-medicine.

Spontaneous viral clearance in patients with acute hepatitis C can be predicted by repeated measurements of serum viral load. H Hofer, T Watkins-Riedel, O Janata, E Penner, H Holzmann, P Steindl-Munda, A Gangl, P Ferenci, Hepatology. 37Hofer H, Watkins-Riedel T, Janata O, Penner E, Holzmann H, Steindl- Munda P, Gangl A, Ferenci P: Spontaneous viral clearance in patients with acute hepatitis C can be predicted by repeated measurements of serum viral load. Hepatology 2003, 37:60-64.

Clinical and bone marrow effects of interferon alfa therapy in myelofibrosis with myeloid metaplasia. A Tefferi, M A Elliot, S Y Yoon, C Y Li, R A Mesa, T G Call, A Dispenzieri, Blood. 971896Tefferi A, Elliot MA, Yoon SY, Li CY, Mesa RA, Call TG, Dispenzieri A: Clinical and bone marrow effects of interferon alfa ther- apy in myelofibrosis with myeloid metaplasia. Blood 2001, 97:1896.

Managing the neuropsychiatric side effects of interferon-based therapy for hepatitis C. C C Crone, G M Gabriel, T N Wise, Cleve Clin J Med. 713SupplCrone CC, Gabriel GM, Wise TN: Managing the neuropsychiat- ric side effects of interferon-based therapy for hepatitis C. Cleve Clin J Med 2004, 71(Suppl 3):S27-32.

Changes in gene expression during pegylated interferon and ribavirin therapy of chronic hepatitis C virus distinguish responders from nonresponders to antiviral therapy. M W Taylor, T Tsukahara, L Brodsky, J Schaley, C Sanda, M J Stephens, J N Mcclintick, H J Edenberg, L Li, J E Tavis, C Howell, S H Belle, J Virol. 81Taylor MW, Tsukahara T, Brodsky L, Schaley J, Sanda C, Stephens MJ, McClintick JN, Edenberg HJ, Li L, Tavis JE, Howell C, Belle SH: Changes in gene expression during pegylated interferon and ribavirin therapy of chronic hepatitis C virus distinguish responders from nonresponders to antiviral therapy. J Virol 2007, 81:3391-3401.

M O'mahony, Sensory Evaluation of Food: Statistical Methods and Procedures. CRC Press487O'Mahony M: Sensory Evaluation of Food: Statistical Methods and Procedures. CRC Press; 1986:487.

W J Conover, Practical Nonparametric Statistics. 3rd editionConover WJ: Practical Nonparametric Statistics 3rd edition. 1998.

Linear models and empirical bayes methods for assessing differential expression in microarray experiments. G K Smyth, Stat Appl Genet Mol Biol. 33Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol 2004, 3:Article3.

Significance analysis of microarrays applied to the ionizing radiation response. V G Tusher, R Tibshirani, G Chu, Proc Natl Acad Sci. 98Tusher VG, Tibshirani R, Chu G: Significance analysis of micro- arrays applied to the ionizing radiation response. Proc Natl Acad Sci USA 2001, 98:5116-5121.

J R Quinlan, C4.5: Programs for Machine Learning. Morgan Kaufmann PublishersQuinlan JR: C4.5: Programs for Machine Learning Morgan Kaufmann Publishers; 1993.

A literature network of human genes for high-throughput analysis of gene expression. T K Jenssen, A Laegreid, J Komorowski, E Hovig, Nat Genet. 28Jenssen TK, Laegreid A, Komorowski J, Hovig E: A literature net- work of human genes for high-throughput analysis of gene expression. Nat Genet 2001, 28:21-28.

Impaired ability of interferon-alpha-primed dendritic cells to stimulate Th1-type CD4 T-cell response in chronic hepatitis C virus infection. H Miyatake, T Kanto, M Inoue, M Sakakibara, A Kaimori, T Yakushijin, I Itose, M Miyazaki, N Kuzushita, N Hiramatsu, T Takehara, A Kasahara, N Hayashi, J Viral Hepat. 14Miyatake H, Kanto T, Inoue M, Sakakibara M, Kaimori A, Yakushijin T, Itose I, Miyazaki M, Kuzushita N, Hiramatsu N, Takehara T, Kasahara A, Hayashi N: Impaired ability of interferon-alpha-primed den- dritic cells to stimulate Th1-type CD4 T-cell response in chronic hepatitis C virus infection. J Viral Hepat 2007, 14:404-412.

Clinical responders to antiviral therapy of chronic HCV infection show elevated antiviral CD4+ and CD8+ T-cell responses. V Pillai, W M Lee, D L Thiele, N J Karandikar, J Viral Hepat. 14Pillai V, Lee WM, Thiele DL, Karandikar NJ: Clinical responders to antiviral therapy of chronic HCV infection show elevated antiviral CD4+ and CD8+ T-cell responses. J Viral Hepat 2007, 14:318-329.

Analyzing the mechanisms of interferon-induced apoptosis using CrmA and hepatitis C virus NS5A. H J Ezelle, S Balachandran, F Sicheri, S J Polyak, G N Barber, Virology. 281Ezelle HJ, Balachandran S, Sicheri F, Polyak SJ, Barber GN: Analyzing the mechanisms of interferon-induced apoptosis using CrmA and hepatitis C virus NS5A. Virology 2001, 281:124-137.

Decrease in CD3-negative-CD8dim(+) and Vdelta2/Vgamma9 TcR+ peripheral blood lymphocyte counts, low perforin expression and the impairment of natural killer cell activity is associated with chronic hepatitis C virus infection. G Par, D Rukavina, E R Podack, M Horanyi, J Szekeres-Bartho, G Hegedus, M Paal, L Szereday, G Mozsik, A Par, J Hepatol. 37Par G, Rukavina D, Podack ER, Horanyi M, Szekeres-Bartho J, Hegedus G, Paal M, Szereday L, Mozsik G, Par A: Decrease in CD3- negative-CD8dim(+) and Vdelta2/Vgamma9 TcR+ periph- eral blood lymphocyte counts, low perforin expression and the impairment of natural killer cell activity is associated with chronic hepatitis C virus infection. J Hepatol 2002, 37:514-522.

Interferon-stimulated gene expression in black and white hepatitis C patients during peginterferon alfa-2a combination therapy. S Luo, W Cassidy, L Jeffers, K R Reddy, C Bruno, C D Howell, Clin Gastroenterol Hepatol. 3Luo S, Cassidy W, Jeffers L, Reddy KR, Bruno C, Howell CD: Inter- feron-stimulated gene expression in black and white hepati- tis C patients during peginterferon alfa-2a combination therapy. Clin Gastroenterol Hepatol 2005, 3:499-506.

Transcription-based prediction of response to IFNbeta using supervised computational methods. S E Baranzini, P Mousavi, J Rio, S J Caillier, A Stillman, P Villoslada, M M Wyatt, M Comabella, L D Greller, R Somogyi, X Montalban, J R Oksenberg, PLoS Biol. 32Baranzini SE, Mousavi P, Rio J, Caillier SJ, Stillman A, Villoslada P, Wyatt MM, Comabella M, Greller LD, Somogyi R, Montalban X, Oksenberg JR: Transcription-based prediction of response to IFNbeta using supervised computational methods. PLoS Biol 2005, 3:e2.

Class prediction from time series gene expression profiles using dynamical systems kernels. K M Borgwardt, S V Vishwanathan, H P Kriegel, Pac Symp Biocomput. Borgwardt KM, Vishwanathan SV, Kriegel HP: Class prediction from time series gene expression profiles using dynamical systems kernels. Pac Symp Biocomput 2006:547-558.