Abstract
Background:
Treatment with Peginterferon Alpha-2b plus Ribavirin is the current standard therapy for chronic hepatitis C (CHC). However, many host related and viral parameters are associated with different outcomes of combination therapy.Objectives:
The aim of this study was to develop an artificial neural network (ANN) model to predetermine individual responses to therapy based on patient’s demographics and laboratory data.Patients and Methods:
This case-control study was conducted in Tehran, Iran, on 139 patients divided into sustained virologic response (SVR) (n = 50), relapse (n = 50) and non-response (n = 39) groups according to their response to combination therapy for 48 weeks. The ANN was trained 300 times (epochs) using clinical data. To test the ANN performance, the part of data that was selected randomly and not used in training process was entered to the ANN and the outputs were compared with real data.Results:
Hemoglobin (P < 0.001), cholesterol (P = 0.001) and IL-28b genotype (P = 0.002) values had significant differences between the three groups. Significant predictive factor(s) for each group were hemoglobin for SVR (OR: 1.517; 95% CI: 1.233-1.868; P < 0.001), IL-28b genotype for relapse (OR: 0.577; 95% CI: 0.339-0.981; P = 0.041) and hemoglobin (OR: 0.824; 95% CI: 0.693-0.980; P = 0.017) and IL-28b genotype (OR: 2.584; 95% CI: 1.430-4.668;P = 0.001) for non-response. The accuracy of ANN to predict SVR, relapse and non-response were 93%, 90%, and 90%, respectively.Conclusions:
Using baseline laboratory data and host characteristics, ANN has been shown as an accurate model to predict treatment outcome, which can lead to appropriate decision making and decrease the frequency of ineffective treatment in patients with chronic hepatitis C virus (HCV) infection.Keywords
Cholesterol Hepatitis C Chronic Hemoglobin IL28B Protein Human Peginterferon Alfa-2b Ribavirin
1. Background
Hepatitis C virus (HCV) infection is widespread, with an estimated 170 to 180 million individuals infected worldwide and 3–4 million new HCV infections each year (1, 2). The current standard of care in chronic HCV hepatitis is the combination of pegylated interferon (PEG-IFN) and ribavirin (RBV) (3-5), which around 50-60% of patients are responders to this therapy (6). However, treatment response rates differ significantly among infected patients.
A wide variety of predictors are advocated for pretreatment evaluation of response to PEG-IFN plus RBV therapy (5). For instance, while up to 80% of patients with genotypes 2 and 3 infection can be cured, the response rate is only 40–50% in genotype 1 infection. Moreover, patients aged < 40 years experienced higher rates of sustained virologic response (SVR) than those aged > 40 years. Other factors such as viral load and body mass index (BMI) can also affect the response rate to the standard treatment (4, 7-9).
There is a high risk of disease progression to liver cirrhosis and subsequently to hepatocellular carcinoma in patients with an unfavorable therapeutic response or in non-responders (3). It is beneficial to predict the response of patients to PEG-IFN and RBV combination therapy before starting the treatment because therapy can be long, costly, and with many side effects (10).
Previous investigators have used artificial neural network (ANN), as artificial intelligence paradigms, to provide a reliable outcome for clinical problems (11-14). ANN is a mathematical model,which is inspired by biological nervous system. It is composed of simple elements operating in parallel. As in nature, connections between elements largely determine the network function. ANNs recognize complex patterns between inputs and outputs via the learning process. Once the hidden association between input and output has been learned, an ANN can correctly predict output from a given input. The capability of neural networks is due to their special features including nonlinear, adaptive, and parallel processing.
2. Objectives
The aim of this study was to develop an ANN model based on viral and host factors to predict treatment outcomes with PEG-IFN and RBV for each patient.
3. Patients and Methods
3.2. Predictive Variables
To predict the individual response of each patient to chronic hepatitis C (CHC) drug therapy, some viral and host-related factors were selected which their impression on different responses to treatment in HCV infected patients had been proven in former studies. Demographic information (including age, gender, weight, and height) was collected by means of a questionnaire. BMI was calculated as weight (kilograms) divided by height (meters) squared. Patient’s blood samples were taken at the beginning of treatment. Complete blood count, fasting blood sugar, lipid profile (cholesterol and triglyceride), prothrombin time, aspartate transaminase (AST) and alanine transaminase (ALT) values were measured by valid clinical laboratories. Total biopsy score of liver histopathologic feature (using modified histology activity index (ISHAK) scoring system), quantification and qualification of HCV RNA (by reverse transcriptase polymerase chain reaction (RT-PCR) and Amplicor analysis with limit of detection 50 IU/mL), genotype of HCV (by Trugene HCV SNC genotyping assay) and genotype of IL-28B SNP (by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP)) were determined for each patient.
3.4. Statistical Analysis
SPSS11.0 (SPSS Inc., Chicago, Illinois, USA) software was used for statistical analysis. One-way ANOVA analysis and multivariate linear regression were performed to specify significant predictive variables and the odds ratios were calculated to compare the magnitude of various risk factors for the outcomes. Confusion matrices for combination of three groups of data was generated to calculate specificity, sensitivity, positive predictive value (PPV), negative predictive value (NPV), likelihood ratio positive (LR+), likelihood ratio negative (LR-) and accuracy. Significance was defined at the level of P < 0.05.
4. Results
Our study was performed on 139 patients with chronic HCV infection (109 men and 30 women). The mean age was 39.87±12.53. There were 50 cases (35.97%) of SVR, 50 (35.97%) of relapse and 39 (28.06%) of non-responder. The demographic characteristics of the patients are shown in Table 1. There were significant differences between the three groups regarding hemoglobin (P < 0.001), cholesterol (P = 0.001) and IL-28b genotype (P = 0.002). Table 2 shows the proportion of each output category for different genotypes of IL-28b, which indicates the role of protective C-Allele in favorable response to combination therapy. Patients carrying CC and CT alleles are more likely to have at least an initial response to treatment. Multivariate linear regression analysis was used to detect the significant predictive factors of SVR, relapse and non-response. As shown in Table 3, the predictive factor of SVR included hemoglobin (OR: 1.517; 95% CI: 1.233-1.868; P < 0.001), the predictive factor of relapse included IL-28b genotype (OR: 0.577; 95% CI: 0.339-0.981; P = 0.041), and the predictive factor of non-response included hemoglobin (OR: 0.824; 95% CI: 0.693-0.980; P = 0.017) and IL-28b genotype (OR: 2.584; 95% CI: 1.430-4.668; P = 0.001). There were no independent predictive factors to predetermine different responses to combination therapy.
The ANN was trained 300 times (epochs). The mean standard error was 1.2689e - 09. When the training was completed, the network output was similar to the real output. To test the ANN performance, the part of data that was selected randomly and not used in training process was entered the ANN and its output was compared with the real output. Table 4 shows the performance indices of ANN in predicting treatment outcome with Peg-IFN and RBV for each patient, compared to the real output. The accuracy of ANN model for predicting SVR, relapse and non-response were 93%, 90%, and 90%, respectively.
Proportion of SVR, Relapse and Non-Response Groups in Different Genotypes of IL-28b a
| SVR, % | Relapse, % | Non-Response, % | |
|---|---|---|---|
| CC | 43.18 | 43.18 | 13.63 |
| CT | 35.29 | 36.76 | 27.94 |
| TT | 25.92 | 22.22 | 51.85 |
Multivariate Regression Analysis to Detect the Statistically Significant Factor for SVR, Relapsers and Non-Responders a
| Variables | Type of Response | P Value | Odds ratio b(95% CI) |
|---|---|---|---|
| Hb, g/dL | |||
| SVR | < 0.0001 | 1.517 (1.233-1.868) | |
| Relapser | 0.088 | 0.873 (0.744-1.023) | |
| Non-Responder | 0.017 | 0.824 (0.693-0.980) | |
| Chol, mg/dL | |||
| SVR | 0.261 | 1.007 (0.996-1.018) | |
| Relapser | 0.250 | 0.994 (0.984-1.005) | |
| Non-responder | 0.905 | 1.000 (0.989-1.012) | |
| IL-28b | |||
| SVR | 0.263 | 0.677 (0.385-1.188) | |
| Relapser | 0.041 | 0.577 (0.339-0.981) | |
| Non-Responder | 0.001 | 2.584 (1.430-4.668) |
The Performance of ANN in Predicting Treatment Outcome With Peg-IFN and RBV a
| SVR, % | Relapse, % | Non-Response, % | |
|---|---|---|---|
| Sensitivity | 92 | 88 | 79 |
| Specificity | 94 | 92 | 94 |
| PPV | 90 | 86 | 83 |
| NPV | 95 | 93 | 91 |
| LR+ | 15.33 | 11 | 13.16 |
| LR- | 0.08 | 0.13 | 0.22 |
| Accuracy | 93 | 90 | 90 |
5. Discussion
Prediction of response to PEG-IFN plus RBV treatment based on viral and host factors using ANN model was the aim of this study. Hemoglobin was the predictive factor of SVR, IL-28b genotype was the predictive factor of relapse, and hemoglobin and IL-28b genotype were the predictive factors of non-response outcome. The ANN model was able to predict SVR, relapse and non-response outcomes with good accuracies.
The role of many factors in different responses of patients receiving chronic hepatitis C therapy was proved in previous investigations. Among these factors, younger age, female gender, absence of obesity, favorable genotype (genotypes 2 and 3 as opposed to genotypes 1 and 4), minimal or absence of fibrosis and milder hepatitis in case of liver histopathology, low baseline HCV RNA level (< 600 000 IU/mL) were associated with remarkable better response (7-9, 16-21). In this study, there were significant differences between SVR, relapse and non-response groups with respect to hemoglobin, serum level of cholesterol and IL-28b genotype.
Using multivariate LR analysis, higher levels of hemoglobin were associated with increase in SVR rate, which is in agreement with shirakawa et al. results. They found higher pretreatment hemoglobin levels in SVR group compared to Non-SVR (10). It may be against the theory that anemia induced by CHC drug therapy (exclusively due to RBV) can improve the treatment results and occasionally lead to SVR (22). However, RBV dose reduction as a routine interventional method in such patients has been reported by dramatically lower SVR rates and prescription of erythropoiesis-stimulating agents has been shown to be a better approach to improve the general condition of patients and drug compliance (23); therefore, reduced hemoglobin level state may be only an indicator of patient’s better corporal response to medication which increases the chance of SVR.
The level of total cholesterol in SVR group was higher than other groups, whereas it was not an independent predictive factor of treatment outcome. Harrison et al. reported in his retrospective study that elevated serum cholesterol levels have been associated with higher SVR rates through unknown mechanisms. However, increase in SVR rate can be due to statin use in patients with elevated cholesterol level and it needs further trials assessing potential advantages of statins as adjuvant therapy for CHC (24).
In agreement with previous studies, IL-28b genotype is a strong predictor of treatment outcome in HCV patients. The global difference of alleles frequency can explain the ethnic variations in treatment response among different populations (4, 25, 26). In the case of rs12979860 genotype, McCarthy et al. and sharafi et al. reported that patients carrying protective C-allele, had about 6-fold increase in SVR rate compared to CT and TT genotypes. According to our results, the C/C variant of the rs12979860 polymorphism was associated with an increased likelihood of SVR, whereas patients with TT genotype were more likely to be non-responders (27-29).
In former researches, patients who had undetectable HCV RNA at the end of therapy (48 weeks) considered to have SVR or named responders, and non-responders have been classified as patients whom HCV RNA counting did not suppress to undetectable at the end of treatment (11, 12). Apart from these, in responders group, if HCV RNA becomes detectable again at week 24 after cessation of therapy, patient is considered to have relapsed. It is important to differentiate sustained virologic responders and relapsers, because relapsers may profit from longer courses of treatment or retreatment recommendations. Therefore, dividing the data into three SVR, relapse and non-response categories and using IL-28b SNPs polymorphism in the set of inputs made this study unique and validated the results.
In earlier studies, logistic regression (LR) models were mainly used as a non-invasive, technical method to predict treatment outcomes (30-32). On the other hand, in some articles the performance differences between two LR and ANN models were discussed in which ANN showed a significantly better performance (12, 33). Considering all these cases, an ANN model was designed which is a non-linear statistical data modeling tool. ANN has the benefit of being able to learn non-linear interconnectivity of inputs and correlations between inputs and outputs by using a set of observations and put them into continuous functions to generate an accurate predictive model without the need of understanding the underlying relationships (13, 14, 34).
Results and calculated performance parameters for each output category showed that designed ANN was able to develop an accurate, non-invasive and effective method, which can be applied on computer-based models for clinical purposes, receiving routine and inexpensive pretreatment clinical data of CHC infected patients and estimating the final response to treatment. The small number of entrance data (especially non-responders group) may be responsible for subsided accuracies and modeling could be extended using additional groups of data. This model should be validated in other populations before clinical implementation. By using such pretreatment predictive strategies in health and medical services, we can obviously reduce the number of patients who may undergo a course of treatment with potential side effects from which they would not drive a benefit. In conclusion, planning a predictive model based on simple and routine laboratory data, by utilizing the ANN, could clearly provide an estimation of how patients respond to PEG-IFN plus RBV therapy, which would be expected to be applied in interventional decision-making.
Acknowledgements
References
-
1.
Kieffer TL, Kwong AD, Picchio GR. Viral resistance to specifically targeted antiviral therapies for hepatitis C (STAT-Cs). J Antimicrob Chemother. 2010;65(2):202-12. [PubMed ID: 19903720]. https://doi.org/10.1093/jac/dkp388.
-
2.
Cooper C. Viral response to specifically targeted antiviral therapy for hepatitis C and the implications for treatment success. Can J Gastroenterol. 2010;24(6):385-90. [PubMed ID: 20559582].
-
3.
Sporea I, Sirli R, Curescu M, Gheorghe L, Popescu A, Bota S, et al. Outcome of antiviral treatment in patients with chronic genotype 1 HCV hepatitis. A retrospective study in 507 patients. J Gastrointestin Liver Dis. 2010;19(3):261-4. [PubMed ID: 20922189].
-
4.
Sarrazin C, Susser S, Doehring A, Lange CM, Muller T, Schlecker C, et al. Importance of IL28B gene polymorphisms in hepatitis C virus genotype 2 and 3 infected patients. J Hepatol. 2011;54(3):415-21. [PubMed ID: 21112657]. https://doi.org/10.1016/j.jhep.2010.07.041.
-
5.
European Association for the Study of the L. EASL Clinical Practice Guidelines: management of hepatitis C virus infection. J Hepatol. 2011;55(2):245-64. [PubMed ID: 21371579]. https://doi.org/10.1016/j.jhep.2011.02.023.
-
6.
Alavian SM, Jabbari H, Daryani NE, Torabi Nami M. Hepatitis C Virus: The Rising Concerns and Growing Hopes, Report From the HCV Symposium, Fourth Tehran Hepatitis Congress, november 2011, Tehran, Iran. Hepat Mon. 2012;12(7):423-9. https://doi.org/10.5812/hepatmon.6679.
-
7.
Al Ashgar H, Helmy A, Khan MQ, Al Kahtani K, Al Quaiz M, Rezeig M, et al. Predictors of sustained virological response to a 48-week course of pegylated interferon alfa-2a and ribavirin in patients infected with hepatitis C virus genotype 4. Ann Saudi Med. 2009;29(1):4-14. [PubMed ID: 19139619].
-
8.
Antonucci G, Angeletti C, Vairo F, Longo MA, Girardi E. Age and prediction of sustained virological response to hepatitis C virus (HCV) infection treatment based on 28-day decrease in HCV RNA levels. J Infect Dis. 2009;200(9):1484-5. author reply 1485. [PubMed ID: 19817588]. https://doi.org/10.1086/644507.
-
9.
Snoeck E, Wade JR, Duff F, Lamb M, Jorga K. Predicting sustained virological response and anaemia in chronic hepatitis C patients treated with peginterferon alfa-2a (40KD) plus ribavirin. Br J Clin Pharmacol. 2006;62(6):699-709. [PubMed ID: 17118125]. https://doi.org/10.1111/j.1365-2125.2006.02741.x.
-
10.
Shirakawa H, Matsumoto A, Joshita S, Komatsu M, Tanaka N, Umemura T, et al. Pretreatment prediction of virological response to peginterferon plus ribavirin therapy in chronic hepatitis C patients using viral and host factors. Hepatology. 2008;48(6):1753-60. [PubMed ID: 18925643]. https://doi.org/10.1002/hep.22543.
-
11.
Saludes V, Bracho MA, Valero O, Ardevol M, Planas R, Gonzalez-Candelas F, et al. Baseline prediction of combination therapy outcome in hepatitis C virus 1b infected patients by discriminant analysis using viral and host factors. PLoS One. 2010;5(11). [PubMed ID: 21152430]. https://doi.org/10.1371/journal.pone.0014132.
-
12.
Wang CH, Mo LR, Lin RC, Kuo JJ, Chang KK, Wu JJ. Artificial neural network model is superior to logistic regression model in predicting treatment outcomes of interferon-based combination therapy in patients with chronic hepatitis C. Intervirology. 2008;51(1):14-20. [PubMed ID: 18309244]. https://doi.org/10.1159/000118791.
-
13.
Grossi E, Buscema M. Introduction to artificial neural networks. Eur J Gastroenterol Hepatol. 2007;19(12):1046-54. [PubMed ID: 17998827]. https://doi.org/10.1097/MEG.0b013e3282f198a0.
-
14.
Raoufy MR, Hajizadeh S, Gharibzadeh S, Mani AR, Eftekhari P, Masjedi MR. Nonlinear model for estimating respiratory volume based on thoracoabdominal breathing movements. Respirology. 2013;18(1):108-16. [PubMed ID: 22897148]. https://doi.org/10.1111/j.1440-1843.2012.02251.x.
-
15.
Strader DB, Wright T, Thomas DL, Seeff LB, American Association for the Study of Liver D. Diagnosis, management, and treatment of hepatitis C. Hepatology. 2004;39(4):1147-71. [PubMed ID: 15057920]. https://doi.org/10.1002/hep.20119.
-
16.
Matsumoto A, Tanaka E, Suzuki T, Ogata H, Kiyosawa K. Viral and host factors that contribute to efficacy of interferon-alpha 2a therapy in patients with chronic hepatitis C. Dig Dis Sci. 1994;39(6):1273-80. [PubMed ID: 8200260].
-
17.
Namazee N, Sali S, Asadi S, Shafiei M, Behnava B, Alavian SM. Real response to therapy in chronic hepatitis C virus patients: a study from iran. Hepat Mon. 2012;12(9). [PubMed ID: 23087759]. https://doi.org/10.5812/hepatmon.6151.
-
18.
Koziel MJ, Peters MG. Viral hepatitis in HIV infection. N Engl J Med. 2007;356(14):1445-54. [PubMed ID: 17409326]. https://doi.org/10.1056/NEJMra065142.
-
19.
Foster GR, Fried MW, Hadziyannis SJ, Chaneac M. 189 Treatment of chronic hepatitis C with peginterferon ALFA-2A (40KD) (PEGASYS®) and ribavirin (COPEGUS®): patient age has a marked influence on the individual estimated probability of achieving a sustained virological response. Hepatology. 2003;38, Supplement:246.
-
20.
Taliani G, Gemignani G, Ferrari C, Aceti A, Bartolozzi D, Blanc PL, et al. Pegylated interferon alfa-2b plus ribavirin in the retreatment of interferon-ribavirin nonresponder patients. Gastroenterology. 2006;130(4):1098-106. [PubMed ID: 16618404]. https://doi.org/10.1053/j.gastro.2006.02.016.
-
21.
Alavian SM, Tabatabaei SV, Keshvari M, Behnava B, Miri SM, Elizee PK, et al. Peginterferon alpha-2a and ribavirin treatment of patients with haemophilia and hepatitis C virus infection: a single-centre study of 367 cases. Liver Int. 2010;30(8):1173-80. [PubMed ID: 20629950]. https://doi.org/10.1111/j.1478-3231.2010.02296.x.
-
22.
Krishnan SM, Dixit NM. Ribavirin-induced anemia in hepatitis C virus patients undergoing combination therapy. PLoS Comput Biol. 2011;7(2). [PubMed ID: 21304937]. https://doi.org/10.1371/journal.pcbi.1001072.
-
23.
Jun DW, Tak WY, Bae SH, Lee YJ. Recent trends in the treatment of chronic hepatitis C. Korean J Hepatol. 2012;18(1):22-8. [PubMed ID: 22511899]. https://doi.org/10.3350/kjhep.2012.18.1.22.
-
24.
Harrison SA, Rossaro L, Hu KQ, Patel K, Tillmann H, Dhaliwal S, et al. Serum cholesterol and statin use predict virological response to peginterferon and ribavirin therapy. Hepatology. 2010;52(3):864-74. [PubMed ID: 20568303]. https://doi.org/10.1002/hep.23787.
-
25.
Fonseca-Coronado S, Vaughan G, Cruz-Rivera MY, Carpio-Pedroza JC, Ruiz-Tovar K, Ruiz-Pacheco JA, et al. Interleukin-28B genotyping by melt-mismatch amplification mutation assay PCR analysis using single nucleotide polymorphisms rs12979860 and rs8099917, a useful tool for prediction of therapy response in hepatitis C patients. J Clin Microbiol. 2011;49(7):2706-10. [PubMed ID: 21613433]. https://doi.org/10.1128/JCM.00877-11.
-
26.
Tanaka Y, Nishida N, Sugiyama M, Kurosaki M, Matsuura K, Sakamoto N, et al. Genome-wide association of IL28B with response to pegylated interferon-alpha and ribavirin therapy for chronic hepatitis C. Nat Genet. 2009;41(10):1105-9. [PubMed ID: 19749757]. https://doi.org/10.1038/ng.449.
-
27.
McCarthy JJ, Li JH, Thompson A, Suchindran S, Lao XQ, Patel K, et al. Replicated association between an IL28B gene variant and a sustained response to pegylated interferon and ribavirin. Gastroenterology. 2010;138(7):2307-14. [PubMed ID: 20176026]. https://doi.org/10.1053/j.gastro.2010.02.009.
-
28.
Sharafi H, Alavian SM. IL28B polymorphism, Explanation for Different Responses to Therapy in Hepatitis C Patients. Hepat Mon. 2011;11(12):958-9. [PubMed ID: 22368678]. https://doi.org/10.5812/kowsar.1735143X.794.
-
29.
Sharafi H, Pouryasin A, Alavian SM, Behnava B, Keshvari M, Salimi S, et al. Distribution of IL28B Genotypes in Iranian Patients with Chronic Hepatitis C and Healthy Individuals. Hepat Mon. 2012;12(12). [PubMed ID: 23550102]. https://doi.org/10.5812/hepatmon.8387.
-
30.
Takayama T, Ebinuma H, Tada S, Yamagishi Y, Wakabayashi K, Ojiro K, et al. Prediction of effect of pegylated interferon alpha-2b plus ribavirin combination therapy in patients with chronic hepatitis C infection. PLoS One. 2011;6(12). [PubMed ID: 22164207]. https://doi.org/10.1371/journal.pone.0027223.
-
31.
Bekkering FC, Stalgis C, McHutchison JG, Brouwer JT, Perelson AS. Estimation of early hepatitis C viral clearance in patients receiving daily interferon and ribavirin therapy using a mathematical model. Hepatology. 2001;33(2):419-23. [PubMed ID: 11172344]. https://doi.org/10.1053/jhep.2001.21552.
-
32.
Akbar HO, Al Ghamdi A, Qattan F, Fallatah HI, Al Rumani M. Chronic hepatitis C in saudi arabia: three years local experience in a university hospital. Hepat Mon. 2012;12(9). [PubMed ID: 23087760]. https://doi.org/10.5812/hepatmon.6178.
-
33.
Raoufy MR, Vahdani P, Alavian SM, Fekri S, Eftekhari P, Gharibzadeh S. A novel method for diagnosing cirrhosis in patients with chronic hepatitis B: artificial neural network approach. J Med Syst. 2011;35(1):121-6. [PubMed ID: 20703578]. https://doi.org/10.1007/s10916-009-9348-8.
-
34.
Raoufy MR, Eftekhari P, Gharibzadeh S, Masjedi MR. Predicting arterial blood gas values from venous samples in patients with acute exacerbation chronic obstructive pulmonary disease using artificial neural network. J Med Syst. 2011;35(4):483-8. [PubMed ID: 20703542]. https://doi.org/10.1007/s10916-009-9384-4.