HCV eradication with IFN-based therapy does not completely restore gene expression in PBMCs from HIV/HCV-coinfected patients

To evaluate the impact of hepatitis C virus (HCV) elimination via interferon (IFN)-based therapy on gene expression profiles related to the immune system in HIV/HCV-coinfected patients. We conducted a prospective study in 28 HIV/HCV-coinfected patients receiving IFN-based therapy at baseline (HIV/HCV-b) and week 24 after sustained virological response (HIV/HCV-f). Twenty-seven HIV-monoinfected patients (HIV-mono) were included as a control. RNA-seq analysis was performed on peripheral blood mononuclear cells (PBMCs). Genes with a fold-change (FC) ≥ 1.5 (in either direction) and false discovery rate (FDR) ≤ 0.05 were identified as significantly differentially expressed (SDE). HIV/HCV-b showed six SDE genes compared to HIV-mono group, but no significantly enriched pathways were observed. For HIV/HCV-f vs. HIV/HCV-b, we found 58 SDE genes, 34 upregulated and 24 downregulated in the HIV/HCV-f group. Of these, the most overexpressed were CXCL2, PDCD6IP, ATP5B, IGSF9, RAB26, and CSRNP1, and the most downregulated were IFI44 and IFI44L. These 58 SDE genes revealed two significantly enriched pathways (FDR < 0.05), one linked to Epstein-Barr virus infection and another related to p53 signaling. For HIV/HCV-f vs. HIV-mono group, we found 44 SDE genes that revealed 31 enriched pathways (FDR < 0.05) related to inflammation, cancer/cell cycle alteration, viral and bacterial infection, and comorbidities associated with HIV/HCV-coinfection. Five genes were overrepresented in most pathways (JUN, NFKBIA, PIK3R2, CDC42, and STAT3). HIV/HCV-coinfected patients who eradicated hepatitis C with IFN-based therapy showed profound gene expression changes after achieving sustained virological response. The altered pathways were related to inflammation and liver-related complications, such as non-alcoholic fatty liver disease and hepatocellular carcinoma, underscoring the need for active surveillance for these patients.

For over a decade, interferon (IFN)-based therapy was the mainstay of antiviral treatment against HCV. However, the appearance of direct-acting antivirals (DAAs) has revolutionized HCV therapy because it has allowed almost all treated patients to achieve sustained virological responses (SVR) [22][23][24][25]. HCV clearance after peg-IFN-α/ribavirin treatment decreases the risk of clinical events and death in HIV/HCV-coinfected patients [26,27]. However, liver fibrosis regression after SVR with peg-IFN-α/ribavirin is slow since some of the cured patients maintain a low rate of liver fibrosis progression [28,29] and develop LREs [30][31][32], making monitoring of the cirrhotic patient necessary after the eradication of HCV. Moreover, peg-IFN-α/ribavirin treatment is associated with numerous side effects, some of which may require long-term follow-up [33]. This concern is even greater for HIV/HCV-coinfected patients since, in this population, peg-IFN-α/ribavirin therapy shows higher rates of adverse side effects.
Transcriptome analysis has emerged as a key tool for profiling the immune response against pathogens. Peripheral blood transcriptome analysis may provide relevant information on the host response in chronic hepatitis C [34] and after overcoming HCV infection via spontaneous HCV clearance [35], peg-IFN-α/ribavirin treatment [36], or DAA therapy [37,38]. There is little information about the long-term impact of successful peg-IFN-α/ribavirin treatment on the peripheral blood transcriptome. The only article found on this topic shows that inflammatory gene expression persists in peripheral blood mononuclear cells (PBMCs) after achieving SVR with peg-IFN treatment [39].

Objective
We aimed to evaluate the impact of HCV elimination via IFN-based therapy on gene expression profiles related to the immune system in HIV/HCV-coinfected patients.

Study subjects
We carried out a prospective study (before and after design) on 28 HIV/HCV-coinfected patient samples collected between 2012 and February 2016. In this study, patients were a subgroup from the GESIDA 3603b study (see Additional file 7: Appendix S1) previously described [40]. The study was carried out strictly following the Declaration of Helsinki and was approved by the Research Ethics Committee of the Instituto de Salud Carlos III (CEI PI 23_2011). Before enrolment, all the participants signed written consent.
Additionally, 27 HIV-monoinfected patients (HIVmono) were used as a control group, matched by gender and age with the HIV/HCV-coinfected patient group. Individuals in the HIV-mono group were negative for HCV and HBV infections and had undetectable HIV viral load and CD4 + >500 cells/µl (normal standard for HIVinfected patients). We included this HIV-mono group as a reference point for HIV-monoinfected patients with stable HIV disease.

Clinical data and Samples
Clinical data were collected prospectively using an online form. Subsequently, the database was monitored to verify the data collected. Peripheral venous blood samples were collected in EDTA tubes, and PBMCs were isolated by Ficoll-Paque gradient. The viable PBMCs were stored in liquid nitrogen at the Spanish HIV HGM BioBank until use. The clinical data and samples used in this study were from the beginning of HCV treatment (HIV/HCV-b) and 24 weeks after SVR (HIV/HCV-f ).

RNA extraction, library preparation, and sequencing
The RNeasy Minikit (Qiagen ™ ) was used for total RNA extraction from PBMCs following the manufacturer's instructions. Then, the RNA quantity and quality were evaluated by Nanodrop 2000 and 2100 Bioanalyzer RNA NANO assay (Agilent). Only those samples with an RNA integrity number higher than 7.5 were selected for sequencing. DNAse treatment was performed with the RNase-free DNase Set (Qiagen). Poly-A RNA libraries and sequencing were performed at the Centre for Genomic Regulation in Barcelona (Spain). Briefly, Illumina's TruSeq Stranded mRNA Sample Prep Kit v2 was used on 500 nanograms of total RNA for library synthesis. In this way, coding and multiple forms of noncoding RNAs are captured because of polyadenylation. Ten libraries were multiplexed and pooled on each line of the Illumina HiSeq2500 sequencer to obtain an average of 25 million reads per sample, with a single read, 50nts (1x50), approach.
We followed a specific bioinformatics pipeline that has been set out in detail in Additional file 1: Table S1.
We used FastQC (v. 0.11.8) for the quality control test. The adapter trimming step was carried out with Trimmomatic (v. 0.33). TopHat (v. 2.0.14) with GRCH38 was used as a reference genome to perform the alignment and HTSeq for count extraction (v. 0.6.1).  showed an average of about 24 million reads, from which 99.07% of them mapped to the genome. Raw sequence data are publicly available at the ArrayExpress repository (EMBL-EBI; https:// www. ebi. ac. uk/ array expre ss/) under the accession number E-MTAB-10232. From this initial gene universe, we proposed to analyze 4723 genes related to the immune system according to InnateDB (Additional file 2: Table S2), a public database of human immune genes (www. innat edb. com). However, we only found 3935 (83%) genes in our initial gene universe out of the total listed in the InnateDB (Fig. 1).

Calculation of the sample size
The sample size for each group was calculated according to the RNASeqPower [41], which established a minimum of 25.6 samples per group to reject the null hypothesis (the two groups are equal). Calculates were performed using the following parameters: 25 million reads per sample, an estimated dispersion of 0.4, an alpha of 0.05, and a minimum fold change of 1.5. Thus, the longitudinal study (28 vs. 28) and the independent comparisons (28 vs. 27) showed a statistical power of 0.92.

Statistical analysis
We used the Statistical Package for the Social Sciences (SPSS) v22.0 software (IBM Corp., Chicago, USA) to perform descriptive statistical analyzes of the study populations. Differences were calculated by Mann-Whitney tests and Pearson's chi-squared test.
We used the R statistical package version v3.4.1 (R Foundation for Statistical Computing, Vienna, Austria) for gene expression analysis. The Trimmed Mean of M (TMM) was used to normalize the data (R-package edgeR v3.20.9) after selecting genes by using the list of immune system-related genes (see the previous section). We carried out the TMM normalization for each of the comparisons (HIV/HCV-b vs. HIV-mono, HIV/HCV-f vs. HIV/HCV-b, and HIV/HCV-f vs. HIV-mono). Next, we used the FilterByExpr function with default parameters to apply a minimum expression criterion and select just those genes well represented among all samples. Finally, we explored the gene expression between groups by a generalized linear mixed model (GLMM) for paired samples and a generalized linear model (GLM) for nonpaired comparison (R-package lme4 V. 1.1-23), both with a negative binomial distribution. Finally, we used NetworkAnalyst v3.0 to find enriched pathways and networks.
We calculated the fold-change (FC) and log2-FC for each gene for measuring the changes in different groups. P-values were two-tailed and corrected for multiple testing using the false discovery rate (FDR) with the Benjamini and Hochberg method (q-values) to reduce the risk of spurious results. Significantly differentially expressed (SDE) genes were those with FDR ≤0.05 and log2-FC ≤ − 0.584 (1.5-fold downregulated) or log2-FC ≥0.584 (1.5-fold upregulated). Table 1 shows the baseline characteristics of the patients included in this study. Overall, HIV/HCV-coinfected patients showed higher alcohol intake values, HIV acquisition by injection drug use, and FIB4. Also, HIV/HCVcoinfected individuals had lower values of CD4+ T-cells than HIV-monoinfected patients.

Differences in gene expression at baseline (HIV/HCV-b vs. HIV-mono)
A total of 521 genes fitted the expression criteria for the comparison between HIV/HCV-b vs. HIV-mono group ( Fig. 1), but only six SDE genes were found (Fig. 2a, full description in Additional file 3: Table S3A): four upregulated in the HIV/HCV-b group (IL23A, FKBP15, CALR, and DDIT3) and two downregulated (TLR5 and STXBP1). The analysis of the enriched pathways did not show statistically significant results.

Changes in gene expression during follow-up (HIV/HCV-f vs. HIV/HCV-b)
A total of 1124 genes fulfilled the expression criteria ( Fig. 1), and 58 genes were SDE (Fig. 2b, full description in Additional file 3: Table S3B). Thirty-four genes were upregulated in the HIV/HCV-f group, but only six genes were higher than 2-fold (CXCL2, PDCD6IP, ATP5B, IGSF9, RAB26, and CSRNP1). Twenty-four genes were downregulated in the HIV/HCV-f group, but only two genes were less than 0.5-fold (IFI44 and IFI44L).
We searched for significantly enriched pathways in the SDE genes list, and two significant pathways were identified (Additional file 4: Table S4; q-value <0.05). One was the Epstein-Barr virus infection pathway (q-value=0.025), represented by four upregulated genes (CDKN1A, NFKB2, RELB, and PDIA3) and two downregulated genes (BID and HLA-A). The other was the p53 signaling pathway (q-value= 0.025), represented by two SDE upregulated genes (CDKN1A and PMAIP1) and two downregulated genes (BID and TP73).

Differences in gene expression at the end of the follow-up to the control group (HIV/HCV-f vs. HIV-mono)
Finally, 504 genes fulfilled the expression criteria when comparing HIV/HCV-f vs. HIV-mono (Fig. 1), and 44 genes were SDE (Fig. 2c, full description in Additional file 3: Table S3C). Twenty-six SDE genes were upregulated in the HIV/HCV-f group, but only five were higher than 2-fold upregulated (KLF6, HSPA5, JUN, PRRC2C, and PPP1R15A). Eighteen SDE genes were downregulated, but none were less than 0.5-fold.
The pathway analysis showed 31 enriched pathways related to four main topics (Fig. 3, full description in Additional file 5: Table S5): (i) inflammation, (ii) cancer/ cell cycle alteration, (iii) viral and bacterial infection, and (iv) comorbidities related to HIV/HCV-coinfection.

Discussion
This study shows that the gene expression profile in PBMCs of HIV/HCV-coinfected patients 24 weeks after achieving SVR (HIV/HCV-f ) was significantly different from either HIV/HCV-coinfected patients before initiation of IFN-based therapy (HIV/HCV-b) or the control HIV group (HIV-mono). The highest differences were identified after HCV eradication compared to baseline (HIV/HCV-f vs. HIV/HCV-b). Unexpectedly some genes also remained overexpressed or suppressed after HCV clearance with respect to the HIV-monoinfected patients (HIV/HCV-f vs. HIV-mono). However, HIV/HCV-b showed only discrete expression differences from the HIV-mono, and no significantly enriched pathways were found.

HIV/HCV-b vs. HIV-mono
Unexpectedly, when HIV/HCV-b and HIV-mono were compared, we found only six SDE genes and no biological pathways overrepresented, indicating that HCV infection had a moderate impact on PBMCs gene expression in HIV/HCV-coinfected patients. Kottilil et al. showed that although HCV mono-infection induces significant PBMC gene expression changes related to immunoregulatory and proinflammatory pathways, the differences between HIV/HCV-coinfected and HIV-monoinfected patients were small [42]. Our results also agree with previous findings showing that HCV infection leads to hyporesponsiveness of PBMCs [43,44]. Furthermore, continuous stimulation of TLR2, TLR4, and TLR5 has been shown to lead to monocyte tolerance [45].
All the upregulated genes in the HIV/HCV-b group (IL23A, FKBP15, CALR, and DDIT3) are associated with HCV infection. IL23A is a crucial cytokine in the immune and inflammatory response that plays an essential role in HCV infection. Its increase may enhance the antiviral activity of IFN-based therapy by modulating the molecules expressed by Th17 cells in patients infected with HCV [46]. HCV promotes significant changes in the endoplasmic reticulum (ER) [47]. CALR is a chaperone that interacts with HCV to correctly fold their glycoproteins, avoiding the appearance of chronic stress in the ER when the folding capacity of cellular proteins is exceeded [48]. FKBP15 is a binding protein directly involved in regulating endosome trafficking, which HCV hijacks to produce viral particles [49]. Therefore, a higher FKBP15 expression may be related to the higher activation of the endosome-based secretory pathway for a higher virus production through exosome pathways. Besides, FKBP15 has been identified as mutated in HCV-infected cirrhotic liver and is associated with HCC [50]. DDIT3 is another gene related to the ER that is activated by the HCV core protein after promoting ER stress [51].
All the downregulated regulated genes (TLR5 and STXBP1) have also been linked to HCV infection. STXBP1 is a protein crucial in vesicle docking and fusion. Its downregulation seems to be related to a lower HCV release [52] and a lower cytolytic T lymphocyte activity [53]. TLR5 recognizes explicitly bacterial flagellin. However, stimulation of TLR2 by bacterial products induces TLR5 downregulation in human monocytes [45]. In this sense, chronic hepatitis C increases bacterial translocation to the blood, and downregulation of TLR5 may be beneficial for the host in the context of persistent immune stimulation by reducing the immune response and the associated immunopathology.
Finally, TLR5 was down-regulated by more than twofold in the HIV/HCV-b group, but its expression was normalized at the end of the study when HCV eradication reduces microbial translocation in HIV/HCV-coinfected patients [54]. The remaining five SDE genes (IL23A, FKBP15, CALR, DDIT3, and STXBP1) stayed deregulated at the end of follow-up, indicating that the ER-related and exosome pathways may remain altered after HCV elimination with IFN-based therapy. Therefore, of the six SDE genes, when HCV infection disappears, only the expression of TLR5 was normalized. However, the biological significance of this result is difficult to ascertain since the genes did not group in any biological pathway. Besides, the expression differences with respect to the HIV-mono group were small (FC<2).

HIV/HCV-f vs. HIV/HCV-b groups
Gene expression in PBMCs from HIV/HCV-coinfected patients showed the deregulation of 58 genes after achieving SVR. These SDE genes represent two biological pathways, Epstein-Barr virus (EBV) infection and p53 signaling.
EBV establishes a latent infection in over 90% of human adults, alternating between episodes of latency and reactivation inside B-lymphocytes, which may lead to the malignant transformation of B cells. Thus, EBV has been linked to a broad spectrum of human malignancies [55]. Additionally, EBV coinfection may alter HIV-infected patients' immunological response, increasing inflammation and immune activation [55]. However, these SDE genes deregulated in the EBV pathway are also involved in many other immune pathways that may have been affected after HCV treatment.
On the one hand, BID and HLA-A were downregulated after HCV eradication. BID encodes a proapoptotic member of the Bcl-2 protein family, and its downregulation may indicate a decrease of apoptosis after HCV clearance with HCV therapy. Liver BID suppression improves inflammation and liver fibrosis in experimental NASH [56] and can protect from the development of HCC [57,58]. HLA-A is a member of the major histocompatibility complex (MHC) from class Ia involved in antigen presentation to CD8 + T cells during an immune response [59]. A decrease of HLA-A expression may be related to a reduction of CD8 + T cell response after HCV clearance with HCV therapy. However, increased HLA-A expression levels are related to poor control of HIV through HLA-E induction and the inhibition of NK cells [60]. Therefore, the downregulation of BID and HLA-A genes could improve liver disease and immune response.
On the other hand, CDKN1A, NFKB2, RELB, and PDIA3 genes were upregulated after HCV eradication. CDKN1A encodes p21, a known inhibitor of the cell cycle, linked to senescence arrest and regulated by p53. In the liver, p21 overexpression is related to the arrest of liver regeneration and the development of NAFLD, liver fibrosis, and cirrhosis [61]. However, p21 may have two contradictory functions related to HCC (tumor suppressor or promoter), depending on its subcellular localization [62]. Thus, its nuclear location is associated with a tumor suppression function. Besides, IFN may keep p21 in the nucleus, where it develops its anti-hepatocarcinogenesis role [62]. Therefore, p21 upregulation after IFN-treatment may be associated with a reduced risk of HCC compared to the baseline. P21 also has an antiviral function, and its upregulation blocks HIV replication in HIV-infected patients with elite control phenotype [63,64]. NFKB2 encodes the nuclear factor kappa-B (NF-κB) p100, a subunit of the NF-κB complex, which may function as a repressor or transcriptional activator, depending on its dimerization partner. As a transcriptional repressor, NF-κB2 may support the establishment of HIV latency, while as a transcriptional activator, it may facilitate HIV transcription [65]. NF-κB2 has an essential role in HCC pathogenesis [66] and NAFLD [67]. RELB also encodes a transcription factor that interacts with NF-κB2, and it may also facilitate HIV transcription [65,68]. RelB expression is increased in chronic hepatitis C and promotes liver fibrosis [69,70]. PDIA3 encodes an isomerase enzyme that participates in the synthesis of glycoproteins such as HLA-I. It is a redox sensor by activating , and it regulates cell growth, death, and signal transduction via STAT3 [71]. PDIA3 deregulation is associated with diverse pathological conditions [71]. PDIA3 overexpression can promote liver fibrosis [72] and HCC [73,74]. Overall, the upregulation of CDKN1A, NFKB2, RELB, and PDIA3 indicates a worse prognosis of liver disease and an impaired immune response against HIV.
We have also identified SDE genes related to the p53 signaling pathway. P53 is a tumor suppressor protein that restricts tumorigenesis by regulating cell cycle arrest, senescence, and apoptosis. The p53 signaling pathway has been previously associated with HCV infection as p53 function is disrupted by direct interaction with viral particles [75] or by the host immune response [76]. Thus, p53 function disruption contributes to hepatocellular carcinogenesis. Two genes of this pathway were downregulated (BID and TP73) and two upregulated (CDKN1A and PMAIP1). BID and CDKN1A are both involved in the pathway of EBV infection that has been previously discussed. TP73 encodes p73, which is considered a p53-related transcription factor involved in cell cycle regulation and apoptosis. TP73 is up-regulated in tumors of HCC patients [77], and lower p73 expression is related to higher survival in HCC patients [78]. PMAIP1 encodes a pro-apoptotic protein of the Bcl-2 family that is regulated by p53. PMAIP1 overexpression is linked to NASH and fibrosis histological criteria [79], which could indicate a worse prognosis of liver disease. However, in vitro studies have shown PMAIP1 plays a crucial role in hepatocarcinogenesis by limiting cancer cell survival [80] and an anti-tumoral function [81]. Therefore, PMAIP1 upregulation after HCV eradication would be related to a reduced risk of HCC in these patients.

HIV/HCV-f vs. HIV-mono
The comparison of HIV/HCV-f vs. HIV-mono also revealed the upregulation of inflammatory and carcinogenic genes after HCV clearance. Specifically, we found 44 SDE genes and 31 significantly enriched pathways related to inflammation, cancer, complications associated with HIV/HCV coinfection, and viral and bacterial infections. Overall, five upregulated SDE genes were overrepresented in most pathways (JUN, NFKBIA, PIK3R2, CDC42, and STAT3). JUN encodes a subunit of the AP-1 transcription factor that promotes inflammation and insulin resistance. JUN overexpression promotes liver fibrosis and correlates with progression from steatosis to NASH [82]. NFKBIA encodes an inhibitor of NF-κB. NFKBIA overexpression is related to steatosis development by induction of hepatocyte apoptosis and secretion of TNF-α and IL-1β by Kupffer cells [83]. Together with bacterial translocation, these mediators can activate hepatic stellate cells and promote liver fibrosis and establish the conditions for the development of HCC [84]. Moreover, NFKBIA overexpression may attenuate NF-κB signaling and enhance HIV-1 latency [85,86]. PIK3R2 encodes a negative regulatory protein of the PI3K pathway. Increased PIK3R2 expression has been found in cancer, and it has been proposed as an oncogene [87]. CDC42 encodes a Rho GTPase that regulates the cell cycle, and its overexpression is related to the development of fibrosis [88], NAFLD [89], and HCC [90]. STAT3 encodes a transcription factor that regulates proliferation, cell survival, and the immune response. STAT3 signaling is related to liver inflammation, injury, steatosis, fibrosis, and HCC [91], and it is overexpressed in patients with NASH [92] and HCC [93].
The HIV/HCV-f group showed five SDE genes that were 2-fold upregulated with respect to the HIV-mono (KLF6, HSPA5, JUN, PRRC2C, and PPP1R15A). A brief discussion of these genes is provided in Additional file 6: Table S6.

Consequences of HCV elimination by IFN-based therapy
Overall, our results revealed that HCV elimination after IFN-based treatment positively and negatively affects HIV/HCV-coinfected patients.
On the one hand, IFN-based therapy eliminates HCV infection, improving some aspects of inflammation, immune activation, and liver disease progression. Thus, we found an improvement at the end of followup in non-invasive markers of liver disease (liver stiffness, FIB-4, AST, ALT, and albumin) and a significant increase in CD4 + T-cells/mm 3 , improving the patient's health status. Besides, we have recently described in this cohort (GESIDA 3603b study) a decrease in some plasma and T-cell biomarkers related to inflammation and immune activation in HIV/HCV-coinfected patients after achieving SVR with IFN-based treatment. Still, most biomarkers in plasma and T-cell did not improve and were far from normalization compared to HIV-monoinfected patients [40]. Moreover, we explored the specific signatures of lymphocytes T CD4 + and CD8 + , based on previous transcriptome profiling on these cells [94,95], but none of these genes were identified as differentially expressed after HCV eradication.
On the other hand, HCV elimination with IFN-based therapy does not lead to the restoration of normal physiological conditions in some cases. As previously discussed, some SDE genes and enriched pathways are related to increased inflammation and the development of NASH, NAFLD, or HCC. This may be due to multiple factors, some of which are the following: Firstly, IFN treatment was usually dispensed to patients in advanced stages of liver disease (F ≥2), where severe liver damage already exists and may not completely reverse after reaching SVR. Persistent histological inflammation, which may contribute to hepatic fibrosis and carcinogenesis, was detected several years after successful IFN-based therapy [96]. Besides, HCV elimination after DAAs therapy does not restore the epigenetic changes related to HCC development induced by HCV infection [97,98], which indicates that the failure to normalize the liver injury after HCV elimination is not unique to IFN-based treatment. In our study, the follow-up time of HIV/HCV-coinfected patients was long enough (24 weeks after SVR) to observe an increase in plasma cholesterol levels at the end of follow-up. This finding is consistent with data recently published by Carrero et al. [99], which reports an increase in plasma cholesterol and cardiovascular risk in HIV/HCV-coinfected patients who achieved SVR with IFN-based therapy. However, we have not observed the development of other adverse outcomes, such as the development of NAFLD or HCC.
Secondly, it may be due to the immunomodulatory effects of IFN-based treatment. In this regard, our findings are in line with previous investigations. For example, IFN treatment induces a robust cytokine response mediated by STAT1, NF-κB and JNK signaling pathways in PBMCs of HCV-infected patients with SVR compared to non-responder patients [100]. An inflammatory response was also observed in PMBCs after successful IFN treatment of HCV infection [39]. Moreover, as previously commented, we have reported in this cohort (GESIDA 3603b study) that a large number of biomarkers did not improve after achieving SVR with IFN-based therapy, staying away from the values of the HIV control group [40].
Thirdly, another possible explanation resides in the elimination of HCV itself. Chronic hepatitis C profoundly affects the host immune response and other physiological processes [101]. Therefore, HCV elimination may induce undesired effects in the already severely affected immune system and liver physiology, particularly in the context of HIV coinfection.

Limitations of the study
Our study has some limitations. Therefore, our conclusions should be interpreted with caution: (i) The number of patients per group was limited, which could prevent the detection of small differences between groups and increase the risk of false positives results. However, the GLMM approach properly accounts for the random effect in our model, limiting the false-positive rate and increasing the statistical power. (ii) Our study had a prospective observational design and may introduce biases. (iii) Another limitation of this study was the lack of other control groups that could give our findings more robustness. So, we lack an HCV monoinfected group of patients with similar liver disease characteristics to the HIV/HCV group. This control group would help to elucidate gene expression profiles related to the immune system after the HCV clearance with IFN treatment without HIV interaction. Nor were available a group of HCV/HIV-coinfected patients who failed the IFN treatment or HCV/HIV-coinfected patients receiving DAA treatment (IFN-free treatment) as controls. These control group options would have helped us to elucidate whether the differences observed in this study were due to the effect of HCV clearance or IFN treatment per se. iv) iv) This work is an association study and lacks mechanistic exploration or functional validation. Functional validation by quantitative PCR is not usually necessary when using a large number of biological replicates. Given the complex virushost interactions, it is difficult to know which genetic changes are truly clinically important and contribute to clinical outcomes.

Conclusions
HIV/HCV-coinfected patients who eradicated chronic hepatitis C with peg-IFN-α/ribavirin therapy showed profound gene expression changes in PBMCs, which seem to be related to inflammation and liver-related complications, such as NASH and HCC. Our findings may be influenced by the long-term impact of IFN therapy, although a rebound effect due to HCV elimination itself could not be ruled out. Prospective studies should be carried out to monitor the long-term normalization of gene expression profiles and clinical evolution in these patients who achieved successful SVR after IFN-based regimens.