# IFI35, mir-99a and HCV Genotype to Predict Sustained Virological Response to Pegylated-Interferon Plus Ribavirin in Chronic Hepatitis C

Emilie Estrabaud, Kevin Appourchaux, Ivan Bieche, Fabrice Carrat, Martine Lapalus, Olivier Lada, Michelle Martinot-Peignoux, Nathalie Boyer, Patrick Marcellin, Michel Vidaud, et al.

► To cite this version:

Emilie Estrabaud, Kevin Appourchaux, Ivan Bieche, Fabrice Carrat, Martine Lapalus, et al. IFI35, mir-99a and HCV Genotype to Predict Sustained Virological Response to Pegylated-Interferon Plus Ribavirin in Chronic Hepatitis C. PLoS ONE, 2015, 10 (4), pp.e0121395. 10.1371/journal.pone.0121395. hal-01233435

HAL Id: hal-01233435

https://hal.sorbonne-universite.fr/hal-01233435

Submitted on 25 Nov 2015

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# IFI35, mir-99a and HCV Genotype to Predict Sustained Virological Response to Pegylated-Interferon Plus Ribavirin in Chronic Hepatitis C

Emilie Estrabaud<sup>1,2,3,4*</sup>, Kevin Appourchaux<sup>1,2,3,4</sup>, Ivan Bièche<sup>5</sup>, Fabrice Carrat<sup>6</sup>, Martine Lapalus<sup>1,2,3,4</sup>, Olivier Lada<sup>1,2,3,4</sup>, Michelle Martinot-Peignoux<sup>1,2,3,4</sup>, Nathalie Boyer<sup>1,2,3,4</sup>, Patrick Marcellin<sup>1,2,3,4</sup>, Michel Vidaud<sup>5</sup>, Tarik Asselah<sup>1,2,3,4</sup>

1 INSERM, UMR1149, Team «Viral hepatitis », Centre de Recherche sur l'inflammation, BP 416, Paris, France, 2 Université Denis Diderot Paris 7, site Bichat, BP 416, Paris, France, 3 Service d'Hépatologie, PMAD Hôpital Beaujon, 100 Bd du Général Leclerc, Clichy la Garenne, Clichy, France, 4 Laboratory of Excellence Labex INFLAMEX, PRES Paris Sorbonne Cité, France, 5 UMR745 INSERM, Université Paris Descartes, Sorbonne Paris Cité, Faculté des Sciences Pharmaceutiques et Biologiques, Paris, France, 6 Unité de santé publique, Hôpital Saint-Antoine, Assistance Publique Hôpitaux de Paris and UMR-S 707, UPMC Université Paris 06 & INSERM, Paris, France

* emilie.estrabaud@inserm.fr

<figure><img src="image_3.png" /></figure>

OPEN ACCESS

**Citation:** Estrabaud E, Appourchaux K, Bièche I, Carrat F, Lapalus M, Lada O, et al. (2015) IFI35, mir-99a and HCV Genotype to Predict Sustained Virological Response to Pegylated-Interferon Plus Ribavirin in Chronic Hepatitis C. PLoS ONE 10(4): e0121395. doi:10.1371/journal.pone.0121395

**Academic Editor:** Chen-Hua Liu, National Taiwan University Hospital, TAIWAN

**Received:** October 17, 2014

**Accepted:** January 31, 2015

**Published:** April 6, 2015

**Copyright:** © 2015 Estrabaud et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** EE fellowship was supported by the french national agency for research on AIDS and viral hepatitis (ANRS). The study was supported by fundings from the french national agency for research on AIDS and viral hepatitis (ANRS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

**Competing Interests:** The authors have declare that no competing interests exist.

## Abstract

Although, the treatment of chronic hepatitis C (CHC) greatly improved with the use of direct antiviral agents, pegylated-interferon (PEG-IFN) plus ribavirin remains an option for many patients, worldwide. The intra-hepatic level of expression of interferon stimulated genes (ISGs) and the rs12979860 CC genotype located within IFNL3 have been associated with sustained virological response (SVR), in patients with CHC. The aim of the study was to identify micro-RNAs associated with SVR and to build an accurate signature to predict SVR. Pre-treatment liver biopsies from 111 patients, treated with PEG-IFN plus ribavirin, were studied. Fifty-seven patients had SVR, 36 non-response (NR) and 18 relapse (RR). The expression of 851 human miRNAs and 30 selected mRNAs, including ISGs, was assessed by RT-qPCR. In the first group of patients (screen), 20 miRNAs out of the 851 studied were deregulated between NRs and SVRs. From the 4 miRNAs validated (mir-23a, mir-181a*, mir-217 and mir-99a), in the second group of patients (validation), 3 (mir-23a, mir-181a* and mir-99a) were down-regulated in NRs as compared to SVRs. The ISGs, studied, were accumulated in SVRs and IFNL3 rs12979860 CT/TT carriers compared respectively to NRs and CC carriers. Combining, clinical data together with the expression of selected genes and micro-RNAs, we identified a signature (IFI35, mir-99a and HCV genotype) to predict SVR (AUC:0.876) with a positive predictive value of 86.54% with high sensibility (80%) and specificity (80.4%). This signature may help to characterize patients with low chance to respond to PEG-IFN/ribavirin and to elucidate mechanisms of NR.
**Introduction**

Hepatitis C virus (HCV) is a major cause of chronic liver disease, with more than 170 million infected individuals, worldwide [1,2].

The treatment of chronic hepatitis C (CHC) has greatly improved, with the development of direct antiviral agents (DAAs), reviewed in [3,4]. Whereas interferon-free combinations are very effective, pegylated-interferon (PEG-IFN) remains an option for a large number of patients worldwide mainly because of cost issues and access to treatment. Therefore, the prediction of SVR to IFN based therapy is still highly important to identify (i) patients with high chance to cure and therefore candidates for this therapy and (ii) those with low chance to respond to PEG-IFN plus ribavirin, candidates for IFN free therapies. Moreover, the prediction of non-response may help to elucidate its molecular mechanisms. The SVR rate is highly associated with viral factors; patients with genotype 1 and 4 show significantly lower response rate than patients with genotype 2 and 3. Different liver gene expression profiles have been associated with non-response, reviewed in [5]. Interferon stimulated genes (ISGs), have been reported over-expressed, at baseline, in patients with non-response (NR) compared to sustained virological responders (SVR) [6–9]. It has been suggested that in patients with high expression of ISGs, at baseline, the addition of exogenous interferon was not able to activate ISGs expression further [7,9].

Single nucleotide polymorphisms (SNPs), in particular rs12979860, located upstream *IFNL3*, were strongly associated with SVR [10,11]. Patients, carrying the favourable genotype rs12979860 CC, have higher chance to develop a SVR than patients with CT and TT genotypes.

Micro-RNAs (miRNAs) are small non-coding RNAs controlling the expression of up to 60% of total mRNAs and therefore are involved in the regulation of major cellular pathways [12,13].

Both *in vitro* and *in vivo* models showed variation of cellular miRNAs expression upon HCV infection [14–16]. Moreover, IFNα/β up-regulates several cellular miRNAs (mir-196, mir-296, mir-351, mir-431, mir-1, mir-30 and mir-128 and mir-448) with putative recognition sites within HCV genome [17]. Mir-196, mir-296, mir-351, mir-431 and mir-448 were indeed able to substantially attenuate viral replication [17]. Mir-122 is highly expressed in the liver; its over-expression stimulates HCV replication, *in vitro* [18]. The use of anti mir-122 (SPC3649, miravirsen) inhibits viral replication both in animal model and in HCV infected patients [19–21]. A reduction of mir-122 expression has been reported, at baseline, in patients with primary non-response as compared to patients with early virological response [22]. In a previous work, we showed that mir-122 expression was decreased in patients with IFNL3 CT/TT genotypes and that this association was stronger than the one between mir-122 and the response to the treatment [23].

We aimed to identify a specific miRNAs expression profile associated with SVRs and NRs, in patients with chronic hepatitis C. We studied the expression of liver miRNAs in two independent groups of patients (screen and validation groups). We investigated the expression of validated miRNAs, in serum samples. We built a predictive signature of SVR taking into account clinical and biological data, *IFNL3* rs12979860 genotype and intra-hepatic mRNAs and miRNAs expression.

**Patients and Methods**

**Characteristics of patients**

One hundred eleven patients with CHC, followed at Beaujon Hospital, were consecutively included in the study. Pre-treatment percutaneous liver biopsies were collected. The study
conformed to the ethical guidelines of the 1975 Declaration of Helsinki. All participants provide their written informed consent to participate in this study.

The patients met the following criteria:

i. an established diagnosis of CHC

ii. the absence of other cause of chronic liver disease and viral co-infection

iii. the patients received a complete treatment of either PEG-IFNα-2b (MSD) at a dose of 1.5 µg/kg/week and weight-based ribavirin at a dose of 800–1200 mg/day, or PEG-IFNα-2a at a dose of 180 µg/week (Roche) and weight-based ribavirin 1000–1200 mg/day. Duration of treatment was 48 weeks for HCV genotype 1 and 4 and 24 weeks for genotypes 2 and 3. For HCV genotype 1 and 4, If HCV-RNA was detectable at 24 weeks of therapy, treatment was stopped (non response).

iv. the patients had an adequate follow-up: detection of serum HCV RNA by RT-PCR was performed at week 12, week 24, at the end of treatment and at 24 weeks after the end of treatment. SVR was defined as undetectable HCV RNA 24 weeks after completion of treatment, while NR was defined as detectable serum HCV RNA at the end of the treatment [24,25]. Responder-relapsers (RR) was defined as the reappearance of detectable serum HCV RNA within 24 weeks, after treatment cessation in patients with no detectable serum HCV RNA at the end of treatment [24,25]. All the patients included in this study had an excellent compliance.

A total of 111 patients were included in two independent groups (Table 1): the screen group (n = 28; 14 NRs and 14 SVRs) and the validation group (n = 83; 42 SVRs, 22 NRs and 18 RRs). The differential expression of the 851 human miRNAs in NRs and SVRs was assessed in the screen group, therefore the same number of patients with SVR and NR were included. Liver biopsies from the validation group (n = 83; 42 SVRs, 22 NRs and 18 RRs) were used in an independent set of experiments, to confirm the results. We enlarged the population with SVRs, NRs and RRs, in accordance with real life data. The expression of the validated mRNAs was assessed in available serum samples from 68 patients of our cohort. The 20 miRNAs deregulated in the screen group (showing at least a 2-fold difference between NRs and SVRs) were studied independently in the validation group.

A selective analysis was performed in the total group of 111 patients. The aim was to study the expression of 6 miRNAs (mir-122, mir-196b, mir-296-5p, mir-448, mir-431 and mir-218) and 30 mRNAs (S1 Table). While the 6 miRNAs were previously associated with HCV infection [17,18], the 30 selected genes were mainly involved in IFN signalling, cell adhesion and cell junction or encoding growth factors (S1 Table) [26,27]. Using liver biopsies, the fibrosis was staged according to the METAVIR scoring system: F0, no fibrosis; F1, portal fibrosis without septa; F2, portal fibrosis with rare septa; F3, numerous septa without cirrhosis; F4, cirrhosis [28].

## Histologically normal controls

Percutaneous liver biopsy specimens were taken from 9 adults with mild elevation of alanine aminotransferase (ALT) with no cause of liver disease (medication, alcohol, chronic viral hepatitis, autoimmune processes and metabolic disease) [29]. All the specimens were histologically normal with absence of inflammation, fibrosis, pathological pattern and less than 5% steatosis.
<table><caption>Table 1. Characteristics of the patients.</caption><thead><tr><th>Variable</th><th></th><th>Screen</th><th>Validation</th><th>Total</th><th>SVR/NR Pvalues<sup>b</sup></th></tr></thead><tbody><tr><td>n</td><td></td><td>28</td><td>83</td><td>111</td><td></td></tr><tr><td>Gender: male (%) / female (%)</td><td></td><td>14(50) / 14 (50)</td><td>56 (67) / 27 (33)</td><td>70 (63) / 41 (37)</td><td>0.1315</td></tr><tr><td>Age [mean ± SD (range)]</td><td></td><td>50.5 ± 7.5 (34–68)</td><td>46.6 ± 9.1 (28–73)</td><td>47.1 ± 8.9 (28–73)</td><td></td></tr><tr><td>BMI, kg.m-2 [mean ± SD (range)]</td><td></td><td>25.08 ± 3.64 (19.37–33.05)</td><td>25.49 ± 4.21 (17.64–37.83)</td><td>25.45 ± 4.19 (17.64–37.83)</td><td>0.0967</td></tr><tr><td>Alanine minotransferase (ALT) IU/L [mean ± SD (range)]</td><td></td><td>103.0 ± 97.9 (20–352)</td><td>107.1 ± 67.1 (18–459)</td><td>105.1 ± 72.5 (18–459)</td><td>0.5508</td></tr><tr><td>Aspartate aminotransferase (AST) IU/L [mean ± SD (range)]</td><td></td><td>65.0 ± 57.6 (25–233)</td><td>63.87 ± 35.49 (21–323)</td><td>64.07 ± 34.80 (21–323)</td><td>0.1565</td></tr><tr><td>Gamma-glutamyltransferase (GGT) IU/L [mean ± SD (range)]</td><td></td><td>64.73 ± 90.78 (3–315)</td><td>90.31 ± 86.27 (7–411)</td><td>84.98 ± 87.20 (4–411)</td><td>0.0040</td></tr><tr><td>Total-cholesterol mmol/L [mean ± SD (range)]</td><td></td><td>4.80 ± 0.77 (3.43–6.49)</td><td>4.70 ± 0.92 (3.11–7.36)</td><td>4.73 ± 0.87 (3.11–7.36)</td><td>0.0829</td></tr><tr><td>Glycemia mmol/L [mean ± SD (range)]</td><td></td><td>5.0 ± 0.63 (3.2–6.3)</td><td>5.40 ± 1.50 (3.14–13.6)</td><td>5.2 ± 1.31 (3.2–13.6)</td><td>0.9353</td></tr><tr><td>Viral loads: copies per mL [median (range)]</td><td></td><td>4.63.106 (1.4.104–2.73.107)</td><td>1.87.106 (3.103–3.06.107)</td><td>2.97.106 (3.103–3.06.107)</td><td>0.0032</td></tr><tr><td rowspan="6">HCV genotypes, n(%)</td><td>1</td><td>19 (68)</td><td>38 (46)</td><td>57 (51.2)</td><td rowspan="6">1 vs. 2,3,5,6 0.0029</td></tr><tr><td>2</td><td>3 (10.5)</td><td>9 (10.5)</td><td>12 (10.8)</td></tr><tr><td>3</td><td>1 (3.5)</td><td>13 (15.5)</td><td>14 (13)</td></tr><tr><td>4</td><td>4 (14.5)</td><td>22 (26.5)</td><td>26 (23.2)</td></tr><tr><td>5</td><td>1 (3.5)</td><td>1 (1.5)</td><td>1 (0.9)</td></tr><tr><td>6</td><td></td><td>1 (1.5)</td><td>1 (0.9)</td></tr><tr><td rowspan="3">IFNL3 rs12979860 polymorphism, n(%)</td><td>C/C</td><td>11 (39.3)</td><td>20 (24.1)</td><td>31 (28)</td><td rowspan="3">CT vs. CC 0.2634 TT vs. CC 0.0062</td></tr><tr><td>C/T</td><td>7 (25)</td><td>33 (39.8)</td><td>40 (36)</td></tr><tr><td>T/T</td><td>8 (28.6)</td><td>19 (22.9)</td><td>27 (24.3)</td></tr><tr><td rowspan="5">Fibrosis Score (Metavir), n(%)</td><td>1</td><td>11 (39.2)</td><td>23 (27.7)</td><td>34 (30.6)</td><td rowspan="5">0.5835</td></tr><tr><td>2</td><td>16 (57.2)</td><td>24 (28.9)</td><td>40 (36.1)</td></tr><tr><td>3</td><td>1 (3.6)</td><td>17 (20.5)</td><td>18 (16.2)</td></tr><tr><td>4</td><td></td><td>17 (20.5)</td><td>17 (15.3)</td></tr><tr><td>unknown</td><td></td><td>2 (2.4)</td><td>2 (1.8)</td></tr><tr><td rowspan="3">Treatment response, n(%)<sup>a</sup></td><td>SVR</td><td>14 (50)</td><td>43 (52)</td><td>57 (51.2)</td><td></td></tr><tr><td>NR</td><td>14 (50)</td><td>22 (26.5)</td><td>36 (32.4)</td><td></td></tr><tr><td>RR</td><td></td><td>18 (21.5)</td><td>18 (16.2)</td><td></td></tr></tbody></table>

doi:10.1371/journal.pone.0121395.t001

## IFNL3 genotyping

The patients were genotyped for rs12979860 using direct sequencing (Applied Biosystems), as previously described [10]. Primers used are available on request. Free circulating DNA was extracted from 500μl serum sample (Qiagen). The PCR products were separated on an ABI3130 sequencer, and analysed with SEQSCAPE 2.6 (Applied Biosystems).

## Statistical analysis

As the miRNA and mRNA levels did not fit a Gaussian distribution, (a) both the miRNA and mRNA levels in each subgroup of samples were characterized by their median values and ranges, rather than their mean values and coefficients of variation, and (b) relationships

<footnote><sup>a</sup>Results are expressed as mean ± SD (range). NRs, non-responders; SVRs, sustained virological responders; RRs, responder-relapsers.</footnote>

<footnote><sup>b</sup> statistical assessment of each variable in NRs and SVRs in the total group.</footnote>
between the molecular markers and histological and viral parameters were tested using the nonparametric Mann-Whitney *U* test (Mann and Whitney, 1947). For miRNAs, differences between two populations were judged significant at confidence levels greater than 95% (*p* < 0.05) and a fold between the two populations more than 2 or a confidence level greater than 99% (*p*<0.01).

The selected miRNAs were assessed in the validation sample (excluding patients who were RR) using the same methodology. Deregulated miRNAs identified at this validation stage were candidate covariates for further analysis.

To identify independent predictors of SVR, we performed logistic regression modelling comparing NRs and SVRs in the total group (screen + validation, excluding RR patients). We performed univariate analyses to test the association between each candidate covariate and the response using the Wald test and we calculated the associated odds-ratios (OR). Covariates with a p-value <0.05 were selected for further multivariate analysis. In multivariate analysis a backward selection strategy was used to select independent predictors. The performance of the predictive selected model was assessed using the area under the receiver operating characteristic (ROC) curve. We used bootstrapping with 1,000 replicates of the total sample for internal validation of the predictive performance of the model as described in a previous study [30].

The material and methods of the RNA extraction and the RT-q-PCR are detailed in the *S1 Methods*.

# Results

## Patients

Baseline characteristics of the 111 patients studied are presented in Table 1. Patients are mainly male (n = 70 in the total group). The mean age at liver biopsy was 47.1 years. Genotypes were represented as follows with 51.2%, 10.8%, 13% and 23.2% of patients with respectively genotypes 1, 2, 3 and 4, in the total group of patients. Genotypes 1 and 4 are the most frequent in patients followed at Beaujon Hospital.

Patients from the screen group (n = 28) had mainly mild (F1) (n = 11) and moderate (F2) fibrosis (n = 16). In the screen group, patients with more advanced fibrosis were excluded to avoid identification of miRNAs that may be rather associated with advanced fibrosis than with response to PEG-IFN plus ribavirin. In the validation group, 23, 24, 17 and 17 patients were respectively F1, F2, F3 and F4. The fibrosis stage was unknown for 2 patients. Out of the 111 patients, 57 were SVR, 36 NR and 18 RR.

In the total group (excluding RR patients), Gamma-glutamyltransferase (per IU/L increase, OR = 0.991; 95% CI: 0.986-0.997; p = 0.0040), HCV viral load (per log10 cp/ml increase, OR = 0.311; 95% CI: 0.145-0.675; p = 0.0032) and HVC genotype 1 (vs. others, OR = 0.093; 95% CI: 0.020-0.444; p = 0.0029) were associated with lower rates of SVR. Moreover, *IFNL3* rs12979860 TT genotype was associated with NR (vs. CC or CT, OR = 0.271; 95%:CI 0.101-0.729; p = 0.0097).

## MiRNAs expression

**Global miRNA approach.** Out of the 851 miRNAs studied, 20 were particularly deregulated (confidence levels greater than 95% (*p* < 0.05) and a fold between the two population more than 2 or a confidence level greater than 99% (*p*<0.01) between NRs and SVRs, in the screen group) (Fig. 1A and Table 2A). Four out of the 20 miRNAs selected showed the same deregulation of expression in the screen group and the validation group (Fig. 1B and Table 2B). Mir-99a, mir-181a\* and mir-23a were down-regulated in NRs whereas mir-217 was accumulated in NRs compared to SVRs.
The OR were calculated for each of the 4 miRNAs, 3 were statistically significant (mir-99a, mir-181a* and mir-23a) (Fig. 2). The AUC of the model consisting in the combination of mir-99a, mir-181a* and mir-23a was 0.7346 (Fig. 2A). The sensibility, specificity, positive predictive value and negative predictive value of each miRNAs and the model combining the 3 miRNAs are represented on Fig. 2B.

**Candidate miRNAs approach.** In our study, mir-122 was neither significantly decreased in NRs compared to SVRs in the screen group (p = 0.077) nor in the total group (p = 0.06) (Fig. 2C). None of the miRNAs regulating HCV replication in vitro (mir-196b p = 0.499, mir-296-5p p = 0.253, mir-431 p = 0.128, mir-448 p = 0.971) showed difference of expression prior to treatment, in NRs and SVRs (Fig. 2C). The expression of mir-296-5p (p = 0.002), mir-122 (p<0.001), mir-448 (p<0.001) and mir-128 (p<0.001) was significantly down-regulated in patients with chronic hepatitis C compared to normal non-infected patients (Fig. 2C).

The expression of mir-23a, mir-99a, mir-181a*, mir-217 and mir-122 was investigated, in 68 serum samples available (NR = 26, RR = 10, RR = 32) (Fig. 3). The level of expression of mir-23a and mir-99a were low compared to mir-122 (Fig. 3A). None of mir-23a, mir-99a and mir-122 was differentially expressed between SVRs and NRs (Fig. 3A). Interestingly, only serum mir-122 and alanine amino transferase were strongly correlated (Fig. 3B), as we previously described [23]. Mir-217 and mir-181a* were not detected efficiently in our serum samples (Fig. 3). Moreover, the expression of mir-23a, mir-99a and mir-122 were not correlated in serum and liver samples (Table 3).

## Modification of mRNAs expression in NRs and SVRs

From the 30 mRNAs tested 25 were deregulated between NRs and SVRs (Fig. 4A and Table 4). Interestingly all the ISGs assessed were accumulated in NRs compared to SVRs (Fig. 4A and Table 4). The ORs were calculated for each mRNAs (Table 4).

ISGs were particularly over-expressed in rs12979860 TT (16.125 to 2.07 fold Fig. 4B and S2A and S2B Tables) patients compared to CC (Fig. 4B). The 3 co-receptors of HCV entry occludin (OCLN) CD81 molecule (CD81) and claudin 1 (CLDN1) were only accumulated in NRs compared to SVRs and not in IFNL3 TT compared to CC patients (Table 4 and S2 Table).

## Design of a signature to predict SVR

The clinical and biological data, IFNL3 polymorphism, mRNAs and miRNAs expression were analyzed by multivariate analysis. IFNL3 polymorphism was not independent of ISGs expression (S3 Table). ISGs were particularly over-expressed in NRs compared to SVRs. It is likely that ISGs expression is a stronger predictor than IFNL3 SNPs, considering the results of multivariate analysis. Mir-99a (per unit of relative mir-99a expression increase OR = 6.053 95% CI: 1.897–19.314 p = 0.0024) IFI35 (per unit of relative IFI35 expression increase OR = 0.069 95% CI:0.013–0.361 p = 0.0015) and HCV genotype 1 (vs. other genotypes OR = 0.093 95% CI: 0.020–0.444 p = 0.0029) were combined in a model predictor of SVR (Fig. 5A). The details about the selection of the variables in the multivariate analysis are presented in S3 Table. Overall the variables included in the multivariate analysis only mir-99a, IFI35 and HCV genotype 1 were independent predictors of SVR (p = 0.0051 p = 0.0011 and p = 0.0059 respectively Wald test) (Fig. 5B). Interestingly the model predicts SVR with a positive predictive value of 86.54% (and a negative predictive value of 71.8%), a sensibility of 80% and a specificity of 80.4%. In patients infected with genotype 1 both IFI35 and mir-99a expression were important to predict SVR. However, in patients infected with either genotype 2, 3, 4 or 5, IFI35 expression alone was a predictor of SVR (Fig. 5C).
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Fig 1. Modification of miRNAs expression in SVRs and NRs.** Total RNAs was extracted and analyzed for miRNAs content by RT q-PCR. The ΔCt for each miRNA was calculated and normalized to the ΔCt value of SNORD44, in each biopsy. The histograms represent the mean expression of miRNAs/ SNORD44 within the groups of patients NRs, SVRs and RRs normalized to the expression of the same miRNA within the group of normal patients. The</figcaption></figure>
association of miRNAs expression with viral response was calculated using the Wald test. 1A- Out of 851 miRNAs, 20 miRNAs were particularly deregulated between NRs and SVRs in the screen group (p<0.05 with at least a 2-fold difference between NRs and SVRs or p<0.01). 1B- The expression of the 20 miRNAs previously selected was assessed by RT q-PCR, in the validation group, independently. N = normal non-infected patients, NR = non-responders, RR = responders-relapsers, SVR = sustained virological responders.

doi:10.1371/journal.pone.0121395.g001

For the identification of mir-99a, in the screen group, we use of a validation procedure on an independent sample prevented inflation of type I error and false discovery rate. A simple calculation under the hypothesis that all 20 miRNAs detected in the screening stage were false positive findings (a conservative scenario) showed that the likelihood of 4 (or more) false positive results on the validation sample (using a type I error at 5%) was 1.6%. In other words, it was unlikely to observe 4 miRNAs detected on the validation stage if all miRNAs detected on the test stage were false positive. Further analysis showed that only mir-99a was an independent predictor of SVR and mir-99a was integrated in the signature to predict SVR, with a p-Value of 0.0051 in multivariate analysis. We are confident that this result is unlikely related to inflation of type I error. However, further analysis will be needed to confirm the role of mir-99a in the prediction of SVR.

# Discussion

## Identification of miRNAs associated with SVR

The major novelty of our work consists in the identification of 4 miRNAs (mir-23a, mir-99a, mir-181a*, and mir-217) differentially expressed between NRs and SVRs.

The expression of mir-23a, mir-99a and mir-181a* was increased in SVRs compared to NRs. However, mir-217 was reduced in SVRs compared to NRs. Previous reports described

<table><caption>Table 2. Differentially expressed miRNAs between NRs, SVRs and RR, prior to treatment.</caption><thead><tr><th rowspan="2">miRNAs</th><th colspan="3">Screen</th><th colspan="4">Validation</th></tr><tr><th>NR</th><th>SVR</th><th>pvalue</th><th>NR</th><th>SVR</th><th>RR</th><th>pvalue</th></tr></thead><tbody><tr><td>mir-455-3p</td><td>0.211</td><td>0.498</td><td>0.002</td><td>0.168</td><td>0.213</td><td>0.171</td><td>0.44</td></tr><tr><td>mir-20a*</td><td>0.202</td><td>0.415</td><td>0.009</td><td>0.343</td><td>0.406</td><td>0.355</td><td>0.81</td></tr><tr><td>mir-624*</td><td>0.330</td><td>0.672</td><td>0.006</td><td>0.357</td><td>0.410</td><td>0.269</td><td>0.16</td></tr><tr><td>mir-616</td><td>0.450</td><td>0.904</td><td>0.044</td><td>0.374</td><td>0.351</td><td>0.246</td><td>0.39</td></tr><tr><td>mir-29b</td><td>0.244</td><td>0.482</td><td>0.011</td><td>0.356</td><td>0.276</td><td>0.211</td><td>0.30</td></tr><tr><td>mir-19a</td><td>0.155</td><td>0.301</td><td>0.025</td><td>0.092</td><td>0.101</td><td>0.057</td><td>0.80</td></tr><tr><td><strong>mir-99a</strong></td><td><strong>0.223</strong></td><td><strong>0.412</strong></td><td><strong>0.002</strong></td><td><strong>0.343</strong></td><td><strong>0.478</strong></td><td><strong>0.295</strong></td><td><strong>0.0382</strong></td></tr><tr><td>mir-148a</td><td>0.238</td><td>0.424</td><td>0.000</td><td>0.258</td><td>0.361</td><td>0.224</td><td>0.19</td></tr><tr><td>mir-28</td><td>0.344</td><td>0.581</td><td>0.004</td><td>0.411</td><td>0.497</td><td>0.389</td><td>0.31</td></tr><tr><td>mir-30a-5p</td><td>0.317</td><td>0.487</td><td>0.005</td><td>0.375</td><td>0.352</td><td>0.259</td><td>0.98</td></tr><tr><td>mir-30b</td><td>0.568</td><td>0.870</td><td>0.011</td><td>0.312</td><td>0.412</td><td>0.244</td><td>0.53</td></tr><tr><td>mir-744</td><td>0.449</td><td>0.676</td><td>0.008</td><td>0.422</td><td>0.534</td><td>0.448</td><td>0.30</td></tr><tr><td>mir-30e</td><td>0.251</td><td>0.371</td><td>0.004</td><td>0.367</td><td>0.374</td><td>0.282</td><td>0.91</td></tr><tr><td><strong>mir-181a*</strong></td><td><strong>0.400</strong></td><td><strong>0.591</strong></td><td><strong>0.005</strong></td><td><strong>0.522</strong></td><td><strong>0.622</strong></td><td><strong>0.486</strong></td><td><strong>0.0423</strong></td></tr><tr><td>mir-576_5p</td><td>0.568</td><td>0.772</td><td>0.008</td><td>0.600</td><td>0.697</td><td>0.510</td><td>0.28</td></tr><tr><td><strong>mir-23a</strong></td><td><strong>0.602</strong></td><td><strong>0.763</strong></td><td><strong>0.008</strong></td><td><strong>0.558</strong></td><td><strong>0.804</strong></td><td><strong>0.689</strong></td><td><strong>0.009</strong></td></tr><tr><td>mir-200b*</td><td>1.491</td><td>0.891</td><td>0.005</td><td>1.385</td><td>1.424</td><td>1.619</td><td>0.51</td></tr><tr><td>mir-200b</td><td>3.985</td><td>1.986</td><td>0.024</td><td>2.471</td><td>2.035</td><td>2.604</td><td>0.23</td></tr><tr><td><strong>mir-217</strong></td><td><strong>2.043</strong></td><td><strong>0.770</strong></td><td><strong>0.021</strong></td><td><strong>5.932</strong></td><td><strong>2.042</strong></td><td><strong>2.585</strong></td><td><strong>0.031</strong></td></tr><tr><td>mir-886-3p</td><td>6.899</td><td>2.474</td><td>0.004</td><td>2.081</td><td>1.025</td><td>0.923</td><td>0.1222</td></tr></tbody></table>

doi:10.1371/journal.pone.0121395.t002
A

<figure><img src="image_6.png" /></figure>

<table><thead><tr><th colspan="4">ROC Curves</th></tr></thead><tbody><tr><td>Model (0.7346)</td><td>mir_99a_1 (0.6895)</td></tr><tr><td>mir_181a_2_1 (0.6793)</td><td>mir_23a_1 (0.7125)</td></tr></tbody></table>

B

<table><thead><tr><th></th><th>Se</th><th>Spe</th><th>PPV</th><th>NPV</th></tr></thead><tbody><tr><td>mir_99a_1</td><td>57.10</td><td>66.70</td><td>72.70</td><td>50.00</td></tr><tr><td>mir_181a_2*_1</td><td>57.10</td><td>58.30</td><td>68.10</td><td>46.70</td></tr><tr><td>mir_23a_1</td><td>64.30</td><td>69.40</td><td>76.60</td><td>55.60</td></tr><tr><td>3 miRNAs signature*</td><td>64.30</td><td>69.40</td><td>76.60</td><td>55.60</td></tr></tbody></table>

<figure><img src="image_7.png" /><figcaption>Fig 2. Identification of a 3 miRNAs signature (mir-99a, mir-23a and mir-181a*) predictor of SVR and selective analysis of miRNAs expression in NRs and SVRs. 2A- The receiver Operator Curves (ROC) were calculated for each miRNAs and for the signature based on the combination of the 3 miRNAs. 2B- The</figcaption></figure>
sensibility (Se), specificity (Spe), positive predictive value (PPV) and negative predictive value (NPV) were
calculated for each miRNA isolated and for the combination of the 3 miRNAs. 2C- The expression of mir-122,
mir-196b, mir-296-5p, mir-448, mir-431 and mir-128 was analyzed by RT q-PCR in the total group of patients
(n = 111). The ΔCp (ΔCpt = 2^ΔCpsample) for each miRNAs was calculated and normalized to the ΔCp value of
SNORD44 in each biopsy. The histograms represent the mean expression of miRNA/SNORD44 within the
groups of patients NRs, SVRs and RRs normalized to the expression of the same miRNA within the group of
normal patients. The association of mRNAs expression with viral response was calculated using the Wald
test. N = normal un-infected patients, NR = non-responders, RR = responders-relapsers, SVR = sustained
virological responders.

doi:10.1371/journal.pone.0121395.g002

hepatic expression of mir-99a, mir-23a and mir-181a suggesting processing of the pre-Mir and co-expression of mir-181a* [31–33]. Interestingly, in malignant cells, whereas the expression of some miRNAs is increased, the widespread under-expression of miRNAs is a more common phenomenon [34]. The modification of miRNAs expression profile might suggest that in NRs, hepatic cells may be more likely to undergo malignant process than in SVRs.

Interestingly, Mir-217 has been described as an oncogene [35] and a reduction of mir-99a expression has been reported in HCC tissues [31]. Moreover, restoration of mir-99a, in vitro, dramatically suppressed HCC cell growth by inducing G1 phase cell cycle arrest [31]. Mir-99a is increased during TGFB-induced epithelial to mesenchymal transition (EMT) [36]. Since EMT is involved during liver fibrogenesis, the expression of mir-99a may increase during fibrosis. Surprisingly, in our cohort, the expression of mir-99a is reduced in F2-F4 vs. F1 (p = 0.028) and in F3-F4 vs. F1F2 (p = 0.035) (S1 Fig). The EMT, in the liver may not induce a reduction of mir-99a expression. The reduction of mir-99a in NRs and at later stages of fibrosis may indicate that these two groups of patients have a higher risk to develop HCC.

The over-expression of mir-23a has been described in mouse models developing non-alcoholic steato-hepatitis, and HCC [32]. Moreover, mir-23a down-regulates interleukin-6 receptor [37] whereas high level of circulating IL-6 has been associated with NR [38]. The modulation of mir23a expression may regulate IL6 signalling in patients with CHC.

Interestingly, hepatic mir-181a* was up regulated in patients with SVR compared to NRs after PEG-IFN/RBV therapy [39]. This result supports our data and strongly suggest an up-regulation of mir-181a* expression in patients with SVRs.

Only few putative mRNAs targets of mir-23a, mir-217, mir-99a and mir-181a* have been described, so far. A deeper analysis of these miRNAs activity may help to elucidate the mechanisms of NR to IFN-based therapy.

A previous study investigating the modification of miRNAs expression, didn't report the same set of miRNAs associated with SVR [16]. However, the authors didn't validate their results using an independent group of patients and included only patients with Japanese ancestry [16]. These discrepancies may explain the identification of different miRNAs associated with SVR.

Mir-122 showed a trend of down-regulation in NRs. In a previous work, we reported a reduction of mir-122 in NRs plus RRs as compared to SVRs patients [23]. In agreement, in the current cohort, when RRs were added to the group of NRs, the difference of mir-122 expression between SVRs and NRs plus RRs was statistically significant (p = 0.02).

MiR-196, miR-296-5p, miR-431 and miR-448 are all induced by IFN α/β and contain putative recognition sites within HCV genome [17]. Whereas ISGs were significantly accumulated in NRs, no difference of expression of IFN-induced miRNAs was observed. Interestingly, this data raises the hypothesis that activation of miRNAs induced by IFN and ISGs, in NRs, may involve different signalling pathways.

In the serum, the level of expression of mir-23a and mir-99a were low compared to the one of mir-122 (Fig. 3A). Interestingly, serum mir-99a was down-regulated in patients with HCV
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>Fig 3. Modification of miRNAs expression in the serum of patients with chronic hepatitis C. 3A- Mir-23a, mir-99a, mir-181a*, mir-217 and mir-122 were detected by RT-q-PCR in 68 serums (NR = 26, RR = 10, RR = 32). The ΔCt for each miRNAs were calculated and normalized to the ΔCt value of *c.elegans*</figcaption></figure>
mir-39. The histograms represent the mean expression of miRNAs/mir-39 within the groups of patients NRs, RRs and SVRs. 3B- Correlation between mir-122 and alanine amino transferase, in the serum. The correlation was assessed using the Spearman correlation test. NR = non-responders, RR = responders-relapsers, SVR = sustained virological responders.

doi:10.1371/journal.pone.0121395.g003

compared to normal controls and patients with chronic hepatitis B [40]. This discrepancy may explain the low level of serum mir-99a detected in our cohort. Circulating mir-23a was suggested to be associated with type 2 diabetes [41]. However, in our cohort, none of the patients for whom serum was available, had diabetes. So we were unable to investigate this association. Interestingly, mir-122 was expressed in the serum and strongly correlates with Alanine amino transferase (Fig. 3B), confirming our previous results [23]. The levels of expression of mir-217 and mir-181a* were too low to be detected efficiently (Fig. 3A). To our knowledge, there are no study that reported efficient detection of serum mir-217 and mir-181a* in human serum samples.

## Analysis of mRNAs associated with response

All the ISGs studied were over-expressed in NRs, prior to treatment (Fig. 4A and Table 4). IFI35 is involved in the antiviral functions of IFN against viral infections. The over-expression of IFI35 suppresses the activation of IFNβ and ISG56 promoters [42]. Interestingly, in patients with NR, the higher expression of IFI35 may limit the efficacy of PEG-IFN to activate the expression of ISGs.

An interesting additional finding was the association of ISGs expression with the unfavourable T allele of rs12979860 IFNL3 (Fig. 4B and S2A and S2B Tables), as it was previously described [43–45]. Several studies reported that intra-hepatic ISGs expression is a better predictor of SVR than IFNL3 polymorphisms [44,45]. The expression of a set of 3 ISGs (ISG15, IFI27 and IFI6) has been associated with rs12979860 SNPs [46]. IFNL3 polymorphisms have been suggested to be associated with a specific cell-type modulation of ISGs expression (MxA, PKR, OAS1 and ISG15) in hepatic cells and PBMCs [43].

## Predictive signature of SVR

We identified a signature combining the hepatic expression 3 miRNAs (mir-23a, mir-99a and mir-181a*) to predict SVR (AUC: 0.7346) (Fig. 2A). In our study, the accumulation of intra-hepatic ISGs expression, prior to treatment was more highly associated with NR than with IFNL3 rs12979860 TT genotype. The multivariate analysis of all clinical and biological data, miRNAs

<table><caption>**Table 3. Correlation between hepatic and serum expression of miRNAs.**</caption><tr><td rowspan="2">mir-23a</td><td>R Pearson</td><td>-0.197</td></tr><tr><td>pvalue</td><td>0.325</td></tr><tr><td rowspan="2">mir-99a</td><td>R Pearson</td><td>-0.370</td></tr><tr><td>pvalue</td><td>0.090</td></tr><tr><td rowspan="2">mir-181a*</td><td>R Pearson</td><td>0.293</td></tr><tr><td>p value</td><td>0.573</td></tr><tr><td rowspan="2">mir-217</td><td>R Pearson</td><td>NA*</td></tr><tr><td>p value</td><td>NA*</td></tr><tr><td rowspan="2">mir-122</td><td>R Pearson</td><td>0.038</td></tr><tr><td>p value</td><td>0.781</td></tr></table>

*NA: not applicable.

doi:10.1371/journal.pone.0121395.t003
<figure><img src="image_10.png" /><figcaption>**Fig 4. Differential mRNAs expression in NRs and SVRs prior to treatment.** 4A- The hepatic expression of 25 mRNAs was assessed in the total group of patients by RT q-PCR. The ΔCp (ΔCp<sub>t</sub> = 2<sup>ΔCpsample</sup>) value was calculated for each mRNA expression and normalized to the ΔCp value of RPLP0 in each liver biopsy. The histograms represent the mean value of the expression of mRNA in NRs, SVRs and RRs. The expression values were compared by non-parametric test. 4B- High expression of ISGs is associated with IFNL3 TT genotype. The intra-hepatic expression of 25 mRNAs was analyzed by RT-qPCR</figcaption></figure>
from total RNAs extracts. IFNL3 rs12979860 polymorphism was determined by direct sequencing. The ΔCp (ΔCp<sub>t</sub> = 2<sup>ΔCp<sub>sample</sub></sup>) value was calculated for each mRNA expression and normalized to the ΔCp value of RPLP0 in each liver biopsy. The histograms represent the mean value of the ratio of each mRNA expression in TT + CT/CC patients. ISGs: Interferon stimulated genes non-ISGs: other type of genes. NR = non-responders, RR = responders-relapsers, SVR = sustained virological responders.

doi:10.1371/journal.pone.0121395.g004

and mRNAs expression and *IFNL3* SNPs, allowed us to build a signature (IFI35, mir-99a and HCV genotype) which is associated with a high predictive positive value (86.54%), sensibility (80%) and specificity (80.4%). None of the miRNAs identified showed a correlation between its hepatic and its serum expression (Table 3). Therefore, the expression of mir-99a alone or both mir-23a, mir-99a and mir-181a* has to be assessed in the liver, to provide enough information to predict SVR.

In a multivariable regression model, studying a HCV genotype 1 cohort, *IFNL3* rs12979860 CC genotype was the strongest pre-treatment predictor of SVR (OR = 5.2; 95% CI: 4.1–6.7)

<table><caption>Table 4. Differentially expressed mRNAs between NRs and SVRs in patients with CHC prior to treatment.</caption><thead><tr><th>mRNAs</th><th>OR<sup>a</sup></th><th>95 % CI<sup>b</sup></th><th>pvalue</th></tr></thead><tbody><tr><td>1-GIP2</td><td>0.818</td><td>0.711-0.940</td><td>0.0047</td></tr><tr><td>2-GIP3</td><td>0.035</td><td>0.003-0.358</td><td>0.0048</td></tr><tr><td>3-HERC5</td><td>0.300</td><td>0.130-0.695</td><td>0.0050</td></tr><tr><td>4-IFI27</td><td>0.612</td><td>0.470-0.798</td><td>0.0003</td></tr><tr><td>5-IFI35</td><td>0.069</td><td>0.013-0.361</td><td>0.0015</td></tr><tr><td>6-IFI44</td><td>0.048</td><td>0.005-0.466</td><td>0.0088</td></tr><tr><td>7-IFIT1</td><td>&lt;0.001</td><td>&lt;0.001-0.004</td><td>0.0028</td></tr><tr><td>8-IFIT4</td><td>0.745</td><td>0.544-1.019</td><td>0.0653</td></tr><tr><td>9-IFITM1</td><td>0.266</td><td>0.100-0.711</td><td>0.0083</td></tr><tr><td>10-LGALS3BP</td><td>0.004</td><td>&lt;0.001-0.172</td><td>0.0039</td></tr><tr><td>11-MIX1</td><td>0.843</td><td>0.742-0.957</td><td>0.0082</td></tr><tr><td>12-OAS1</td><td>0.656</td><td>0.492-0.875</td><td>0.0041</td></tr><tr><td>13-OAS2</td><td>0.591</td><td>0.268-1.306</td><td>0.1934</td></tr><tr><td>14-OAS3</td><td>0.765</td><td>0.626-0.934</td><td>0.0086</td></tr><tr><td>15-PLSCR1</td><td>0.327</td><td>0.147-0.729</td><td>0.0062</td></tr><tr><td>16-PSMB9</td><td>0.046</td><td>0.010-0.213</td><td>&lt;.0001</td></tr><tr><td>17-RSAD2</td><td>0.761</td><td>0.590-0.983</td><td>0.0368</td></tr><tr><td>18-STAT1</td><td>0.139</td><td>0.036-0.537</td><td>0.0042</td></tr><tr><td>19-USP18</td><td>0.581</td><td>0.397-0.849</td><td>0.0051</td></tr><tr><td>20-CCL21</td><td>&lt;0.001</td><td>&lt;0.001 - &lt;0.001</td><td>0.0032</td></tr><tr><td>21-CD81</td><td>1.282</td><td>0.792-2.076</td><td>0.3125</td></tr><tr><td>22-CLDN1</td><td>0.523</td><td>0.360-0.760</td><td>0.0007</td></tr><tr><td>23-CXCL9</td><td>1.133</td><td>0.425-3.017</td><td>0.8027</td></tr><tr><td>24-CXCL10</td><td>0.014</td><td>&lt;0.001-0.324</td><td>0.0077</td></tr><tr><td>25-CXCL11</td><td>0.371</td><td>0.183-0.751</td><td>0.0059</td></tr><tr><td>26-IL8</td><td>0.711</td><td>0.514-0.982</td><td>0.0382</td></tr><tr><td>27-ITGA2</td><td>0.159</td><td>0.045-0.565</td><td>0.0045</td></tr><tr><td>28-MDK</td><td>0.005</td><td>&lt;0.001-0.159</td><td>0.0030</td></tr><tr><td>29-OCLN</td><td>&lt;0.001</td><td>&lt;0.001 - &lt;0.001</td><td>0.0070</td></tr><tr><td>30-STMN2</td><td>0.935</td><td>0.777-1.124</td><td>0.4730</td></tr></tbody></table>

doi:10.1371/journal.pone.0121395.t004

<footnote>a Odds ratio;</footnote>

<footnote>b 95% Confidence Interval</footnote>
<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Fig 5.** Identification of a signature predicting SVR. Clinical, miRNAs and mRNAs expression data were analyzed by multivariate analysis. 5A- ORs are presented for the 3 variables of the signature: mir-99a, IFI35 and HCV genotype 1. 5B- Receiver Operator Curves for prediction of viral response of the signature (AUC 0.876) (genotype 1, IFI35 and mir-99a expression) based on IFI35 expression (AUC 0.777) and mir-99a (AUC 0.688). 5C- Schematic representation of the prediction of SVR. Patients with genotype non-1 had high SVR rates when they express low level of IFI35 (I for low, < mean IFI35 expression). HCV genotype 1 infected patients with low level of IFI35 and high level of mir-99a (h for high, > mean mir-99a expression) showed high SVR rates.</figcaption></figure>

doi:10.1371/journal.pone.0121395.g005

with only a sensitivity of 56%, a specificity of 79%, a positive predictive value of 69% and a negative predictive value of 68% [47]. Accordingly, in our cohort, *IFNL3* alone showed only moderate accuracy to predict SVR. The predictive value of *IFNL3* polymorphism is stronger in HCV genotypes 1 and 4 than in genotypes 2 and 3 [10,48]. Therefore, in our cohort combining all genotypes, the accuracy of *IFNL3* may be reduced. The combination of mir-99a, *IFI35* expression and HCV genotype was a higher predictive factor, at baseline, than *IFNL3* polymorphism (Fig. 5).
In the current cohort (111 patients) both the 3-miRNAs-signature (AUC = 0.74) and the multi-variable signature (AUC = 0.876) showed higher accuracy to predict SVR than the 2-ISGs-signature (IFI27 and CXCL9) previously described by our group (AUC = 0.713) [6].

An independent study reported a 4-ISGs-signature (IFI27, GIP3/ISG15, RSAD2 and HTA-TIP2) as stronger predictor of SVR than IFNL3 SNPs [44]. However, the authors only investigated the expression of 8 selected mRNAs in a single group of patients without any validation set [44].

Whereas major improvements have been achieved in the treatment of CHC with the development of DAAs, PEG-IFN and ribavirin remains an option for many patients worldwide. Therefore the prediction of responsiveness to PEG-IFN plus ribavirin remains an important issue. The assessment of the hepatic expression of mir-99a and IFI35 may help to discriminate patients with low chance to respond to PEG-IFN plus ribavirin. In clinical practice, performing a liver biopsy remains useful to stage fibrosis [49], since urgent initiation of treatment is recommended for patients with advanced fibrosis or compensated cirrhosis. Performing different tests (stage fibrosis and predict treatment response) on the same liver biopsy could optimize the use of this invasive procedure. Patients with low chance to respond to PEG-IFN plus ribavirin may be candidates for IFN-free regimen.

## Supporting Information

**S1 Fig. Down-regulation of hepatic mir-99a in patients with advanced fibrosis.** Total RNAs was extracted and analyzed for mir-99a content by RT q-PCR. The ΔCt for mir-99a was calculated and normalized to the ΔCt value of SNORD44 in each biopsy. The histograms represent the mean expression of mir-99a/SNORD44 within the groups of patients with the different stages of fibrosis according to metavir score (F1 to F4), normalized to mir-99a within the group of normal patients. S1A- The level of expression of mir-99a was compared between patients with mild fibrosis (F1, n = 35) and those with moderate fibrosis to cirrhosis (F2-F4, n = 71). S1B- the level of expression of mir-99a was compared in patients with mild and moderate fibrosis (F1-F2, n = 71) and those with advanced fibrosis and cirrhosis (F3-F4, n = 35). (TIF)

**S1 Methods. RNA extraction and RT-q-PCR detailed methods.** The detailed methods for the RNA extraction and the RT-q-PCR methods are described in the S1 Methods.
(DOC)

**S1 Table. List of the genes selected for the analysis.**
(DOC)

**S2 Table. Fold of expression of selected genes between IFNL3 CC and TT patients.** Interferon stimulated genes (ISGs) in 2A and non-ISGs in 2B.
(DOC)

**S3 Table. Description of the multivariate analysis.** The analyses were run separately for miRNAs, baseline characteristics and mRNAs. Multivariate findings were combined in a final multivariate model with backward selection to predict SVR in 93 patients (57 SVRs and 36 NRs).
(DOCX)

## Author Contributions

Conceived and designed the experiments: EE KA IB MV TA. Performed the experiments: EE ML KA OL MMP. Analyzed the data: EE FC. Contributed reagents/materials/analysis tools: MMP NB TA PM. Wrote the paper: EE TA KA.
References

1. Marcellin P, Asselah T, Boyer N. Fibrosis and disease progression in hepatitis C. Hepatology. 2002; 36: S47–56. PMID: 12407576

2. Seeff LB. Natural history of chronic hepatitis C. Hepatology. 2002; 36: S35–46. PMID: 12407575

3. Asselah T, Marcellin P. Second-wave IFN-based triple therapy for HCV genotype 1 infection: simepre-vir, faldaprevir and sofosbuvir. Liver Int. 2014; 34 Suppl 1: 60–68. doi: 10.1111/liv.12424 PMID: 24373080

4. Schinazi R, Halfon P, Marcellin P, Asselah T. HCV direct-acting antiviral agents: the best interferon-free combinations. Liver Int. 2014; 34 Suppl 1: 69–78. doi: 10.1111/liv.12423 PMID: 24373081

5. Estrabaud E, Vidaud M, Marcellin P, Asselah T. Genomics and HCV infection: progression of fibrosis and treatment response. J Hepatol. 2012;

6. Asselah T, Bieche I, Narguet S, Sabbagh A, Laurendeau I, Ripault MP, et al. Liver gene expression sig-nature to predict response to pegylated interferon plus ribavirin combination therapy in patients with chronic hepatitis C. Gut. 2008; 57: 516–524. PMID: 17895355

7. Chen L, Borozan I, Feld J, Sun J, Tannis LL, Coltescu C, et al. Hepatic gene expression discriminates responders and nonresponders in treatment of chronic hepatitis C viral infection. Gastroenterology. 2005; 128: 1437–1444. PMID: 15887125

8. Helbig KJ, Lau DT, Semendric L, Harley HA, Beard MR. Analysis of ISG expression in chronic hepatitis C identifies viperin as a potential antiviral effector. Hepatology. 2005; 42: 702–710. PMID: 16108059

9. Sarasin-Filipowicz M, Oakeley EJ, Duong FH, Christen V, Terracciano L, Filipowicz W, et al. Interferon signaling and treatment outcome in chronic hepatitis C. Proc Natl Acad Sci U S A. 2008; 105: 7034–7039. doi: 10.1073/pnas.0707882105 PMID: 18467494

10. Asselah T, De Muynck S, Broet P, Masliah-Planchon J, Blanluet M, Bieche I, et al. IL28B polymorphism is associated with treatment response in patients with genotype 4 chronic hepatitis C. J Hepatol. 2011; 56: 527–532. doi: 10.1016/j.jhep.2011.09.008 PMID: 21951981

11. Ge D, Fellay J, Thompson AJ, Simon JS, Shianna KV, Urban TJ, et al. Genetic variation in IL28B predicts hepatitis C treatment-induced viral clearance. Nature. 2009; 461: 399–401. doi: 10.1038/nature08309 PMID: 19684573

12. Meister G, Tuschl T. Mechanisms of gene silencing by double-stranded RNA. Nature. 2004; 431: 343–349. PMID: 15372041

13. Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004; 116: 281–297. PMID: 14744438

14. Peng X, Li Y, Walters KA, Rosenzweig ER, Lederer SL, Aicher LD, et al. Computational identification of hepatitis C virus associated microRNA-mRNA regulatory modules in human livers. BMC Genomics. 2009; 10: 373. doi: 10.1186/1471-2164-10-373 PMID: 19671175

15. Ura S, Honda M, Yamashita T, Ueda T, Takatori H, Nishino R, et al. Differential microRNA expression between hepatitis B and hepatitis C leading disease progression to hepatocellular carcinoma. Hepatology. 2009; 49: 1098–1112. doi: 10.1002/hep.22749 PMID: 19173277

16. Murakami Y, Tanaka M, Toyoda H, Hayashi K, Kuroda M, Tajima A, et al. Hepatic microRNA expression is associated with the response to interferon treatment of chronic hepatitis C. BMC Med Genomics. 2010; 3: 48. doi: 10.1186/1755-8794-3-48 PMID: 20969775

17. Pedersen IM, Cheng G, Wieland S, Volinia S, Croce CM, Chisari FV, et al. Interferon modulation of cellular microRNAs as an antiviral mechanism. Nature. 2007; 449: 919–922. PMID: 17943132

18. Jopling CL, Yi M, Lancaster AM, Lemon SM, Sarnow P. Modulation of hepatitis C virus RNA abundance by a liver-specific MicroRNA. Science. 2005; 309: 1577–1581. PMID: 16141076

19. Elmen J, Lindow M, Silahtaroglu A, Bak M, Christensen M, Lind-Thomsen A, et al. Antagonism of microRNA-122 in mice by systemically administered LNA-antimiR leads to up-regulation of a large set of predicted target mRNAs in the liver. Nucleic Acids Res. 2008; 36: 1153–1162. PMID: 18158304

20. Janssen HL, Reesink HW, Lawitz EJ, Zeuzem S, Rodriguez-Torres M, Patel K, et al. Treatment of HCV infection by targeting microRNA. N Engl J Med. 2013; 368: 1685–1694. doi: 10.1056/NEJMoa1209026 PMID: 23534542

21. Lanford RE, Hildebrandt-Eriksen ES, Petri A, Persson R, Lindow M, Munk ME, et al. Therapeutic silencing of microRNA-122 in primates with chronic hepatitis C virus infection. Science. 2010; 327: 198–201. doi: 10.1126/science.1178178 PMID: 19965718

22. Sarasin-Filipowicz M, Krol J, Markiewicz I, Heim MH, Filipowicz W. Decreased levels of microRNA miR-122 in individuals with hepatitis C responding poorly to interferon therapy. Nat Med. 2009; 15: 31–33. doi: 10.1038/nm.1902 PMID: 19122656
23. Estrabaud E, Lapalus M, Broet P, Appourchaux K, De Muynck S, Lada O, et al. Reduction of microRNA 122 expression in IFNL3 CT/TT carriers and during progression of fibrosis in patients with chronic hepatitis C. J Virol. 2014; 88: 6394-6402. doi: 10.1128/JVI.00016-14 PMID: 24672032

24. Ghany MG, Nelson DR, Strader DB, Thomas DL, Seeff LB. An update on treatment of genotype 1 chronic hepatitis C virus infection: 2011 practice guideline by the American Association for the Study of Liver Diseases. Hepatology. 2011; 54: 1433-1444. doi: 10.1002/hep.24641 PMID: 21898493

25. EASL. EASL Clinical Practice Guidelines: management of hepatitis C virus infection. J Hepatol. 2011; 55: 245-264. doi: 10.1016/j.jhep.2011.02.023 PMID: 21371579

26. Asselah T, Bieche I, Laurendeau I, Paradis V, Vidaud D, Degott C, et al. Liver gene expression signature of mild fibrosis in patients with chronic hepatitis C. Gastroenterology. 2005; 129: 2064-2075. PMID: 16344072

27. Bieche I, Asselah T, Laurendeau I, Vidaud D, Degot C, Paradis V, et al. Molecular profiling of early stage liver fibrosis in patients with chronic hepatitis C virus infection. Virology. 2005; 332: 130-144. PMID: 15661146

28. Bedossa P, Poynard T. An algorithm for the grading of activity in chronic hepatitis C. The METAVIR Cooperative Study Group. Hepatology. 1996; 24: 289-293. PMID: 8690394

29. Asselah T, Bieche I, Laurendeau I, Martinot-Peignoux M, Paradis V, Vidaud D, et al. Significant gene expression differences in histologically "Normal" liver biopsies: Implications for control tissue. Hepatology. 2008; 48: 953-962. doi: 10.1002/hep.22411 PMID: 18726958

30. Steyerberg EW, Harrell FE Jr, Borsboom GJ, Eijkemans MJ, Vergouwe Y, Habbema JD. Internal validation of predictive models: efficiency of some procedures for logistic regression analysis. J Clin Epidemiol. 2001; 54: 774-781. PMID: 11470385

31. Li D, Liu X, Lin L, Hou J, Li N, Wang C, et al. MicroRNA-99a inhibits hepatocellular carcinoma growth and correlates with prognosis of patients with hepatocellular carcinoma. J Biol Chem. 2011; 286: 36677-36685. doi: 10.1074/jbc.M111.270561 PMID: 21878637

32. Wang B, Hsu S, Frankel W, Ghoshal K, Jacob ST. Stat3-mediated activation of miR-23a suppresses gluconeogenesis in hepatocellular carcinoma by downregulating G6PC and PGC-1alpha. Hepatology. 2012;

33. Zhou B, Li C, Qi W, Zhang Y, Zhang F, Wu JX, et al. Downregulation of miR-181a upregulates sirtuin-1 (SIRT1) and improves hepatic insulin sensitivity. Diabetologia. 2012; 55: 2032-2043. doi: 10.1007/s00125-012-2539-8 PMID: 22476949

34. Lujambio A, Lowe SW. The microcosmos of cancer. Nature. 2012; 482: 347-355. doi: 10.1038/nature10888 PMID: 22337054

35. de Yebenes VG, Bartolome-Izquierdo N, Nogales-Cadenas R, Perez-Duran P, Mur SM, Martinez N, et al. miR-217 is an oncogene that enhances the germinal center reaction. Blood. 2014; 124: 229-239. doi: 10.1182/blood-2013-12-543611 PMID: 24850757

36. Turcatel G, Rubin N, El-Hashash A, Warburton D. MIR-99a and MIR-99b modulate TGF-beta induced epithelial to mesenchymal plasticity in normal murine mammary gland cells. PLoS One. 2012; 7: e31032. doi: 10.1371/journal.pone.0031032 PMID: 22299047

37. Zhu LH, Liu T, Tang H, Tian RQ, Su C, Liu M, et al. MicroRNA-23a promotes the growth of gastric adenocarcinoma cell line MGC803 and downregulates interleukin-6 receptor. Febs J. 2010; 277: 3726-3734. doi: 10.1111/j.1742-4658.2010.07773.x PMID: 20698883

38. Nattermann J, Vogel M, Berg T, Danta M, Axel B, Mayr C, et al. Effect of the interleukin-6 C174G gene polymorphism on treatment of acute and chronic hepatitis C in human immunodeficiency virus coinfected patients. Hepatology. 2007; 46: 1016-1025. PMID: 17668881

39. Gelley F, Zadori G, Nemes B, Fassan M, Lendvai G, Sarvary E, et al. MicroRNA profile before and after antiviral therapy in liver transplant recipients for hepatitis C virus cirrhosis. J Gastroenterol Hepatol. 2014; 29: 121-127. doi: 10.1111/jgh.12362 PMID: 24033414

40. Akamatsu S, Hayes CN, Tsuge M, Miki D, Akiyama R, Ochi H, et al. Differences in serum microRNA profiles in hepatitis B and C virus infection. J Infect. 2014;

41. Yang Z, Chen H, Si H, Li X, Ding X, Sheng Q, et al. Serum miR-23a, a potential biomarker for diagnosis of pre-diabetes and type 2 diabetes. Acta Diabetol. 2014; 51: 823-831. doi: 10.1007/s00592-014-0617-8 PMID: 24981880

42. Das A, Dinh PX, Panda D, Pattnaik AK. Interferon-inducible protein IFI35 negatively regulates RIG-I antiviral signaling and supports vesicular stomatitis virus replication. J Virol. 2014; 88: 3103-3113. doi: 10.1128/JVI.03202-13 PMID: 24371060

43. Abe H, Hayes CN, Ochi H, Maekawa T, Tsuge M, Miki D, et al. IL28 variation affects expression of interferon stimulated genes and peg-interferon and ribavirin therapy. J Hepatol. 2011; 54: 1094-1101. doi: 10.1016/j.jhep.2010.09.019 PMID: 21145800
44. Dill MT, Duong FH, Vogt JE, Bibert S, Bochud PY, Terracciano L, et al. Interferon-induced gene expression is a stronger predictor of treatment response than IL28B genotype in patients with hepatitis C. Gastroenterology. 2010; 140: 1021-1031. doi: 10.1053/j.gastro.2010.11.039 PMID: 21111740

45. Honda M, Sakai A, Yamashita T, Nakamoto Y, Mizukoshi E, Sakai Y, et al. Hepatic ISG expression is associated with genetic variation in interleukin 28B and the outcome of IFN therapy for chronic hepatitis C. Gastroenterology. 2010; 139: 499-509. doi: 10.1053/j.gastro.2010.04.049 PMID: 20434452

46. Urban TJ, Thompson AJ, Bradrick SS, Fellay J, Schuppan D, Cronin KD, et al. IL28B genotype is associated with differential expression of intrahepatic interferon-stimulated genes in patients with chronic hepatitis C. Hepatology. 2010; 52: 1888-1896. doi: 10.1002/hep.23912 PMID: 20931559

47. Thompson AJ, Muir AJ, Sulkowski MS, Ge D, Fellay J, Shianna KV, et al. Interleukin-28B polymorphism improves viral kinetics and is the strongest pretreatment predictor of sustained virologic response in genotype 1 hepatitis C virus. Gastroenterology. 2010; 139: 120-129 e118. doi: 10.1053/j.gastro.2010.04.013 PMID: 20399780

48. 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: 415-421. doi: 10.1016/j.jhep.2010.07.041 PMID: 21112657

49. Asselah T, Marcellin P, Bedossa P. Improving performance of liver biopsy in fibrosis assessment. J Hepatol. 2014; 61: 193-195. doi: 10.1016/j.jhep.2014.03.006 PMID: 24650692