- Research
- Open Access
- Published:
HCV eradication with IFN-based therapy does not completely restore gene expression in PBMCs from HIV/HCV-coinfected patients
Journal of Biomedical Science volume 28, Article number: 23 (2021)
Abstract
Objective
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.
Methods
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).
Results
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).
Conclusion
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.
Introduction
Hepatitis C virus (HCV) infection is a leading cause of chronic liver disease worldwide [1]. HCV and human immunodeficiency virus (HIV) share transmission routes and many HIV-infected individuals are coinfected with HCV [2]. HIV/HCV-coinfected patients develop chronic hepatitis C over decades of infection, but the cirrhosis progression is faster than in HCV-monoinfected patients [3], leading to higher rates of liver-related events (LREs) such as liver decompensation, end-stage liver disease, hepatocellular carcinoma (HCC), and liver-related death [4,5,6,7]. Chronic HCV infection promotes inflammation, immune activation, and dysregulation of immune function [8, 9], which boosts the development of LRE and other comorbidities in infected people [7,8,9]. HIV infection also aggravates the dysregulation of the immune system [7, 8], since many disturbances remain in patients infected with HIV after suppressive antiretroviral therapy (ART), such as deficits in T-cell functions [10,11,12,13,14,15], immune activation [13], inflammation [16], and dysbiosis [17,18,19]. All these immunological alterations increase the risk of AIDS, non-AIDS-related events, and death [20, 21].
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.
Patients and methods
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.
The selection criteria of HIV/HCV-coinfected patients were: 1) chronic HCV and HIV infection; 2) starting IFN-based therapy (peg-IFN-α/ribavirin or peg-IFN-α/ribavirin/DAAs) and achieving SVR; 3) CD4+ T cell ≥200 cells/µL; 4) stable ART ≥6 months; 4) frozen PBMCs samples available. The exclusion criteria were: 1) hepatitis B virus coinfection; 2) acute hepatitis C; 3) hepatic decompensation; 4) HCC.
Additionally, 27 HIV-monoinfected patients (HIV-mono) 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 HIV-infected 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 non-coding 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).
Genes selected for analysis
Figure 1 shows the flowchart of the gene selection process. The bioinformatics analysis identified a gene universe composed of 27,173 different genes. Samples 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/arrayexpress/) under the accession number E-MTAB-10232.
Flowchart of the study design. A total of 27,173 genes were identified by RNA sequencing. Of these, 4723 genes from the immune system were selected from InnateDB, and 3935 genes were selected for differential expression analysis. The FilterByExpr function was applied for each comparison: a HIV/HCV-b vs. HIV-mono, 521 genes were selected, six were SDE genes not involved in any biological pathway; b HIV/HCV-f vs. HIV/HCV-b, 1124 genes were selected, 58 were SDE genes involved in two biological pathways; c HIV/HCV-f vs. HIV-mono, 504 genes were selected, 44 were SDE genes involved in 31 biological pathways. HIV, human immunodeficiency virus; HCV, hepatitis C virus; HIV/HCV-b, HIV/HCV-coinfected patients at baseline; HIV/HCV-f, HIV/HCV-coinfected patients 24 weeks after SVR; HIV-mono, HIV-monoinfected patients; SDE, significantly differentially expressed
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.innatedb.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 non-paired 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).
Results
Patient characteristics at baseline
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/HCV-coinfected 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.
Volcano plots for differentially expressed genes between groups: a HIV/HCV-b versus HIV-mono, b HIV/HCV-f versus HIV/HCV-b, c HIV/HCV-f versus HIV-mono. Volcano-plot discriminates by FDR and log2(FC). The red vertical lines represent the cut-off of FC 1.5 and the blue horizontal line indicates the cut-off of FDR = 0.05. Green dots represent genes significantly altered with FDR ≤ 0.05 and absolute FC ≥ 1.5. FDR, false discovery rate for multiple comparisons using Benjamini and Hochberg procedure; FC, fold-change; HIV, human immunodeficiency virus; HCV, hepatitis C virus; HIV/HCV-b, HIV/HCV-coinfected patients at baseline; HIV/HCV-f, HIV/HCV-coinfected patients 24 weeks after SVR; HIV-mono, HIV-monoinfected patients
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. Overall, five genes were overrepresented among the majority of the altered pathways (JUN, NFKBIA, PIK3R2, CDC42, and STAT3).
The chord-diagram represents connections between SDE genes (left-side) and pathways (right-side) in comparing HIV/HCV-f versus HIV-mono. Red corresponds to inflammation, green corresponds to cell cycle alteration and cancer, blue corresponds to viral and bacterial infection, and black corresponds to complications related to HIV/HCV-coinfection. HIV, human immunodeficiency virus; HCV, hepatitis C virus; HIV/HCV-b, HIV/HCV-coinfected patients at baseline; HIV/HCV-f, HIV/HCV-coinfected patients 24 weeks after SVR; HIV-mono, HIV-monoinfected patients; SDE, significantly differentially expressed
Finally, the HIV/HCV-b vs. HIV-mono and the HIV/HCV-f vs. HIV-mono comparisons had four upregulated genes in common (IL23A, FKBP15, CALR, and DDIT3) and one downregulated gene (STXBP1) (Fig. 4, Additional file 3: Table S3). Therefore, these five SDE genes remained deregulated at 24-week after SVR. Additionally, three SDE genes were upregulated (HSPA5, KLF6, and PRRC2C), and four downregulated (MKL1, AMIGO3, LHX4, and ITGB7) both in the HIV/HCV-f vs. HIV/HCV-b and the HIV/HCV-f vs. HIV-mono comparisons (Fig. 4, Additional file 3:Table S3).
Venn diagram of the SDE genes in the three comparisons analyzed in the study: a HIV/HCV-b versus HIV-mono, b HIV/HCV-f versus HIV/HCV-b, c HIV/HCV-f versus HIV-mono. The number in each circle corresponds to the number of SDE genes in the different comparisons. The overlapping numbers are the number of SDE genes shared between the different comparisons. The non-overlapping numbers are the unique SDE genes in each comparison. HIV, human immunodeficiency virus; HCV, hepatitis C virus; HIV/HCV-b, HIV/HCV-coinfected patients at baseline; HIV/HCV-f, HIV/HCV-coinfected patients 24 weeks after SVR; HIV-mono, HIV-monoinfected patients
Changes in metabolic, HIV-related disease, and liver disease profiles (HIV/HCV-f vs. HIV/HCV-b)
At the end of follow-up, we found significant increases in levels of cholesterol (p=0.022), CD4+ T-cells/mm3 (p=0.026), and albumin (p=0.045); while values of LSM (p=0.002), FIB4 (p<0.001), AST (p<0.001), and ALT (p<0.001) significantly decreased (Table 2).
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 mTORC1, 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.
Finally, we found several 2-fold deregulated genes. HIV/HCV-f showed two highly downregulated genes (IFI44 and IFI44L) and six highly upregulated genes (CXCL2, PDCD6IP, ATP5B, IGSF9, RAB26, and CSRNP1) compared to HIV/HCV-b. For the sake of brevity, we have discussed these genes in Additional file 6:Table S6.
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 follow-up in non-invasive markers of liver disease (liver stiffness, FIB-4, AST, ALT, and albumin) and a significant increase in CD4+ T-cells/mm3, 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 virus-host 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.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding authors upon reasoned request. HIV/HCV-coinfected patients' raw sequences data are publicly available at the ArrayExpress repository (EMBL-EBI; https://www.ebi.ac.uk/) under the accession number E-MTAB-8275. The raw RNA sequences of HIV-infected patients can also be accessed under the accession number E-MTAB-8249.
References
The Polaris Observatory HCV Collaborators. Global prevalence and genotype distribution of hepatitis C virus infection in 2015: a modelling study. Lancet Gastroenterol Hepatol. 2017;2(3):161–76.
Platt L, Easterbrook P, Gower E, McDonald B, Sabin K, McGowan C, Yanny I, Razavi H, Vickerman P. Prevalence and burden of HCV co-infection in people living with HIV: a global systematic review and meta-analysis. Lancet Infect Dis. 2016;16(7):797–808.
Vallet-Pichard A, Pol S. Natural history and predictors of severity of chronic hepatitis C virus (HCV) and human immunodeficiency virus (HIV) co-infection. J Hepatol. 2006;44(1 Suppl):S28-34.
Lo-Re V, Kallan MJ, Tate JP, Localio AR, Lim JK, Goetz MB, Klein MB, Rimland D, Rodriguez-Barradas MC, Butt AA, et al. Hepatic decompensation in antiretroviral-treated patients co-infected with HIV and hepatitis C virus compared with hepatitis C virus-monoinfected patients: a cohort study. Ann Intern Med. 2014;160(6):369–79.
López-Diéguez M, Montes ML, Pascual-Pareja JF, Quereda C, Von Wichmann MA, Berenguer J, Tural C, Hernando A, González-García J, Serrano L, et al. The natural history of liver cirrhosis in HIV-hepatitis C virus-coinfected patients. AIDS. 2011;25(7):899–904.
Macias J, Berenguer J, Japon MA, Giron JA, Rivero A, Lopez-Cortes LF, Moreno A, Gonzalez-Serrano M, Iribarren JA, Ortega E, et al. Fast fibrosis progression between repeated liver biopsies in patients coinfected with human immunodeficiency virus/hepatitis C virus. Hepatology. 2009;50(4):1056–63.
Ingiliz P, Rockstroh JK. Natural history of liver disease and effect of hepatitis C virus on HIV disease progression. Curr Opin HIV AIDS. 2015;10(5):303–8.
Naggie S. Hepatitis C virus, inflammation, and cellular aging: turning back time. Top Antivir Med. 2017;25(1):3–6.
Shin EC, Sung PS, Park SH. Immune responses and immunopathology in acute and chronic viral hepatitis. Nat Rev Immunol. 2016;16(8):509–23.
Pandiyan P, Younes SA, Ribeiro SP, Talla A, McDonald D, Bhaskaran N, Levine AD, Weinberg A, Sekaly RP. Mucosal Regulatory T Cells and T Helper 17 Cells in HIV-Associated Immune Activation. Front Immunol. 2016;7:228.
Miles B, Miller SM, Connick E. CD4 t follicular helper and regulatory cell dynamics and function in HIV infection. Front Immunol. 2016;7:659.
De Biasi S, Bianchini E, Nasi M, Digaetano M, Gibellini L, Carnevale G, Borghi V, Guaraldi G, Pinti M, Mussini C, et al. Th1 and Th17 proinflammatory profile characterizes invariant natural killer T cells in virologically suppressed HIV+ patients with low CD4+/CD8+ ratio. AIDS. 2016;30(17):2599–610.
Fernandes JR, Berthoud TK, Kumar A, Angel JB. IL-23 signaling in Th17 cells is inhibited by HIV infection and is not restored by HAART: Implications for persistent immune activation. PLoS ONE. 2017;12(11):e0186823.
DaFonseca S, Niessl J, Pouvreau S, Wacleche VS, Gosselin A, Cleret-Buhot A, Bernard N, Tremblay C, Jenabian MA, Routy JP, et al. Impaired Th17 polarization of phenotypically naive CD4(+) T-cells during chronic HIV-1 infection and potential restoration with early ART. Retrovirology. 2015;12:38.
Mahnke YD, Fletez-Brant K, Sereti I, Roederer M. Reconstitution of Peripheral T Cells by Tissue-Derived CCR4+ central memory cells following HIV-1 antiretroviral therapy. Pathog Immun. 2016;1(2):260–90.
de Paula HHS, Ferreira ACG, Caetano DG, Delatorre E, Teixeira SLM, Coelho LE, Joao EG, de Andrade MM, Cardoso SW, Grinsztejn B, et al. Reduction of inflammation and T cell activation after 6 months of cART initiation during acute, but not in early chronic HIV-1 infection. Retrovirology. 2018;15(1):76.
Merlini E, Tincati C, Biasin M, Saulle I, Cazzaniga FA, d’Arminio Monforte A, Cappione AJ 3rd, Snyder-Cappione J, Clerici M, Marchetti GC. Stimulation of PBMC and monocyte-derived macrophages via toll-like receptor activates innate immune pathways in HIV-infected patients on virally suppressive combination antiretroviral therapy. Front Immunol. 2016;7:614.
Mudd JC, Brenchley JM. Gut mucosal barrier dysfunction, microbial dysbiosis, and their role in HIV-1 disease progression. J Infect Dis. 2016;214(Suppl 2):S58-66.
Tincati C, Merlini E, Braidotti P, Ancona G, Savi F, Tosi D, Borghi E, Callegari ML, Mangiavillano B, Barassi A, et al. Impaired gut junctional complexes feature late-treated individuals with suboptimal CD4+ T-cell recovery upon virologically suppressive combination antiretroviral therapy. AIDS. 2016;30(7):991–1003.
Hunt PW, Lee SA, Siedner MJ. Immunologic Biomarkers, Morbidity, and Mortality in Treated HIV Infection. J Infect Dis. 2016;214(Suppl 2):S44-50.
Hart BB, Nordell AD, Okulicz JF, Palfreeman A, Horban A, Kedem E, Neuhaus J, Jacobs DR Jr, Duprez DA, Neaton JD, et al. Inflammation-related morbidity and mortality among hiv-positive adults: how extensive is it? J Acquir Immune Defic Syndr. 2018;77(1):1–7.
Liang TJ, Ghany MG. Therapy of hepatitis C–back to the future. N Engl J Med. 2014;370(21):2043–7.
Panel A-IHG: Hepatitis C Guidance 2018 Update: AASLD-IDSA recommendations for testing, managing, and treating hepatitis C Virus Infection. Clin Infect Dis 2018, 67(10):1477–1492.
European Association for the Study of the Liver. Electronic address eee, European Association for the Study of the L: EASL Recommendations on Treatment of Hepatitis C 2018. J Hepatol. 2018;69(2):461–511.
Berenguer J, Gil-Martin A, Jarrin I, Moreno A, Dominguez L, Montes M, Aldamiz-Echevarria T, Tellez MJ, Santos I, Benitez L, et al. All-oral direct-acting antiviral therapy against hepatitis C virus (HCV) in human immunodeficiency virus/HCV-coinfected subjects in real-world practice: Madrid coinfection registry findings. Hepatology. 2018;68(1):32–47.
Berenguer J, Rodriguez E, Miralles P, Von Wichmann MA, Lopez-Aldeguer J, Mallolas J, Galindo MJ, Van Den Eynde E, Tellez MJ, Quereda C, et al. Sustained virological response to interferon plus ribavirin reduces non-liver-related mortality in patients coinfected with HIV and Hepatitis C virus. Clin Infect Dis. 2012;55(5):728–36.
Berenguer J, Rodriguez-Castellano E, Carrero A, Von Wichmann MA, Montero M, Galindo MJ, Mallolas J, Crespo M, Tellez MJ, Quereda C, et al. Eradication of hepatitis C virus and non-liver-related non-acquired immune deficiency syndrome-related events in human immunodeficiency virus/hepatitis C virus coinfection. Hepatology. 2017;66(2):344–56.
Liu Z, Wei X, Chen T, Huang C, Liu H, Wang Y. Characterization of fibrosis changes in chronic hepatitis C patients after virological cure: a systematic review with meta-analysis. J Gastroenterol Hepatol. 2017;32(3):548–57.
Labarga P, Fernandez-Montero JV, de Mendoza C, Barreiro P, Pinilla J, Soriano V. Liver fibrosis progression despite HCV cure with antiviral therapy in HIV-HCV-coinfected patients. Antivir Ther. 2015;20(3):329–34.
Mira JA, Neukam K, Lopez-Cortes LF, Rivero-Juarez A, Tellez F, Giron-Gonzalez JA, de los Santos-Gil I, Ojeda-Burgos G, Merino D, Rios-Villegas MJ et al: Efficacy of and risk of bleeding during pegylated interferon plus ribavirin treatment in HIV/HCV-coinfected patients with pretreatment thrombocytopenia. Eur J Clin Microbiol Infect Dis 2015; 34(9):1879–1884.
Salmon-Ceron D, Nahon P, Layese R, Bourcier V, Sogni P, Bani-Sadr F, Audureau E, Merchadou L, Dabis F, Wittkop L et al: Human Immunodeficiency Virus/Hepatitis C Virus (HCV) Co-infected Patients With Cirrhosis Are No Longer at Higher Risk for Hepatocellular Carcinoma or End-Stage Liver Disease as Compared to HCV Mono-infected Patients. Hepatology 2018.
Allaire M, Nahon P, Layese R, Bourcier V, Cagnot C, Marcellin P, Guyader D, Pol S, Larrey D, De Ledinghen V, et al. Extrahepatic cancers are the leading cause of death in patients achieving hepatitis B virus control or hepatitis C virus eradication. Hepatology. 2018;68(4):1245–59.
European Association for Study of L. EASL Clinical Practice Guidelines: management of hepatitis C virus infection. J Hepatol. 2014;60(2):392–420.
Rosenberg BR, Freije CA, Imanaka N, Chen ST, Eitson JL, Caron R, Uhl SA, Zeremski M, Talal A, Jacobson IM, et al. Genetic Variation at IFNL4 influences extrahepatic interferon-stimulated gene expression in chronic HCV Patients. J Infect Dis. 2018;217(4):650–5.
Rosenberg BR, Depla M, Freije CA, Gaucher D, Mazouz S, Boisvert M, Bedard N, Bruneau J, Rice CM, Shoukry NH. Longitudinal transcriptomic characterization of the immune response to acute hepatitis C virus infection in patients with spontaneous viral clearance. PLoS Pathog. 2018;14(9):e1007290.
Lu MY, Huang CI, Hsieh MY, Hsieh TJ, Hsi E, Tsai PC, Tsai YS, Lin CC, Hsieh MH, Liang PC, et al. Dynamics of PBMC gene expression in hepatitis C virus genotype 1-infected patients during combined peginterferon/ribavirin therapy. Oncotarget. 2016;7(38):61325–35.
Carlton-Smith C, Holmes JA, Naggie S, Lidofsky A, Lauer GM, Kim AY, Chung RT. of the AAsg: IFN-free therapy is associated with restoration of type I IFN response in HIV-1 patients with acute HCV infection who achieve SVR. J Viral Hepat. 2018;25(5):465–72.
Alao H, Cam M, Keembiyehetty C, Zhang F, Serti E, Suarez D, Park H, Fourie NH, Wright EC, Henderson WA, et al. Baseline intrahepatic and peripheral innate immunity are associated with hepatitis C virus clearance during direct-acting antiviral therapy. Hepatology. 2018;68(6):2078–88.
Waldron PR, Holodniy M. Peripheral blood mononuclear cell gene expression remains broadly altered years after successful interferon-based hepatitis C virus treatment. J Immunol Res. 2015;2015:958231.
Garcia-Broncano P, Medrano LM, Berenguer J, Brochado-Kith O, Gonzalez-Garcia J, Jimenez-Sousa MA, Quereda C, Sanz J, Tellez MJ, Diaz L, et al. Mild profile improvement of immune biomarkers in HIV/HCV-coinfected patients who removed hepatitis C after HCV treatment: a prospective study. J Infect. 2020;80(1):99–110.
Hart SN, Therneau TM, Zhang Y, Poland GA, Kocher JP. Calculating sample size estimates for RNA sequencing data. J Comput Biol. 2013;20(12):970–8.
Kottilil S, Yan MY, Reitano KN, Zhang X, Lempicki R, Roby G, Daucher M, Yang J, Cortez KJ, Ghany M, et al. Human immunodeficiency virus and hepatitis C infections induce distinct immunologic imprints in peripheral mononuclear cells. Hepatology. 2009;50(1):34–45.
Roque-Cuellar MC, Sanchez B, Garcia-Lozano JR, Praena-Fernandez JM, Marquez-Galan JL, Nunez-Roldan A, Aguilar-Reina J. Hepatitis C virus-specific cellular immune responses in sustained virological responders with viral persistence in peripheral blood mononuclear cells. Liver Int. 2014;34(6):e80-88.
Chung H, Watanabe T, Kudo M, Chiba T. Correlation between hyporesponsiveness to Toll-like receptor ligands and liver dysfunction in patients with chronic hepatitis C virus infection. J Viral Hepat. 2011;18(10):e561-567.
Cabral ES, Gelderblom H, Hornung RL, Munson PJ, Martin R, Marques AR. Borrelia burgdorferi lipoprotein-mediated TLR2 stimulation causes the down-regulation of TLR5 in human monocytes. J Infect Dis. 2006;193(6):849–59.
Meng P, Zhao S, Niu X, Fu N, Su S, Wang R, Zhang Y, Qiao L, Nan Y. Involvement of the Interleukin-23/Interleukin-17 axis in chronic hepatitis C virus infection and its treatment responses. Int J Mol Sci. 2016;17:7.
Ravindran MS, Bagchi P, Cunningham CN, Tsai B. Opportunistic intruders: how viruses orchestrate ER functions to infect cells. Nat Rev Microbiol. 2016;14(7):407–20.
Choukhi A, Ung S, Wychowski C, Dubuisson J. Involvement of endoplasmic reticulum chaperones in the folding of hepatitis C virus glycoproteins. J Virol. 1998;72(5):3851–8.
Lai CK, Jeng KS, Machida K, Lai MM. Hepatitis C virus egress and release depend on endosomal trafficking of core protein. J Virol. 2010;84(21):11590–8.
Ikeda A, Shimizu T, Matsumoto Y, Fujii Y, Eso Y, Inuzuka T, Mizuguchi A, Shimizu K, Hatano E, Uemoto S, et al. Leptin receptor somatic mutations are frequent in HCV-infected cirrhotic liver and associated with hepatocellular carcinoma. Gastroenterology. 2014;146(1):222–32.
Wang J, Kang R, Huang H, Xi X, Wang B, Wang J, Zhao Z. Hepatitis C virus core protein activates autophagy through EIF2AK3 and ATF6 UPR pathway-mediated MAP1LC3B and ATG12 expression. Autophagy. 2014;10(5):766–84.
Ren H, Elgner F, Himmelsbach K, Akhras S, Jiang B, Medvedev R, Ploen D, Hildt E. Identification of syntaxin 4 as an essential factor for the hepatitis C virus life cycle. Eur J Cell Biol. 2017;96(6):542–52.
Phatarpekar PV, Billadeau DD. Molecular regulation of the plasma membrane-proximal cellular steps involved in NK cell cytolytic function. J Cell Sci. 2020;133:5.
Lopez-Cortes LF, Trujillo-Rodriguez M, Baez-Palomo A, Benmarzouk-Hidalgo OJ, Dominguez-Molina B, Milanes-Guisado Y, Espinosa N, Viciana P, Gutierrez-Valencia A. Eradication of Hepatitis C Virus (HCV) Reduces Immune Activation, Microbial Translocation, and the HIV DNA Level in HIV/HCV-Coinfected Patients. J Infect Dis. 2018;218(4):624–32.
Shindiapina P, Ahmed EH, Mozhenkova A, Abebe T, Baiocchi RA. Immunology of EBV-Related Lymphoproliferative Disease in HIV-Positive Individuals. Front Oncol. 2020;10:1723.
Eguchi A. De Mollerat Du Jeu X, Johnson CD, Nektaria A, Feldstein AE: Liver Bid suppression for treatment of fibrosis associated with non-alcoholic steatohepatitis. J Hepatol. 2016;64(3):699–707.
Bao L, Zhang M, Han S, Zhan Y, Guo W, Teng F, Liu F, Guo M, Zhang L, Ding G, et al. MicroRNA-500a promotes the progression of hepatocellular carcinoma by post-transcriptionally targeting BID. Cell Physiol Biochem. 2018;47(5):2046–55.
Wree A, Johnson CD, Font-Burgada J, Eguchi A, Povero D, Karin M, Feldstein AE. Hepatocyte-specific Bid depletion reduces tumor development by suppressing inflammation-related compensatory proliferation. Cell Death Differ. 2015;22(12):1985–94.
Crux NB, Elahi S. Human Leukocyte Antigen (HLA) and immune regulation: how do classical and non-classical hla alleles modulate immune response to human immunodeficiency virus and hepatitis C virus infections? Front Immunol. 2017;8:832.
Ramsuran V, Naranbhai V, Horowitz A, Qi Y, Martin MP, Yuki Y, Gao X, Walker-Sperling V, Del Prete GQ, Schneider DK, et al. Elevated HLA-A expression impairs HIV control through inhibition of NKG2A-expressing cells. Science. 2018;359(6371):86–90.
Huda N, Liu G, Hong H, Yan S, Khambu B, Yin XM. Hepatic senescence, the good and the bad. World J Gastroenterol. 2019;25(34):5069–81.
Ohkoshi S, Yano M, Matsuda Y. Oncogenic role of p21 in hepatocarcinogenesis suggests a new treatment strategy. World J Gastroenterol. 2015;21(42):12150–6.
Diez-Fuertes F, De La Torre-Tarazona HE, Calonge E, Pernas M, Alonso-Socas MDM, Capa L, Garcia-Perez J, Sakuntabhai A, Alcami J. Transcriptome sequencing of peripheral blood mononuclear cells from elite controller-long term non progressors. Sci Rep. 2019;9(1):14265.
Leng J, Ho HP, Buzon MJ, Pereyra F, Walker BD, Yu XG, Chang EJ, Lichterfeld M. A cell-intrinsic inhibitor of HIV-1 reverse transcription in CD4(+) T cells from elite controllers. Cell Host Microbe. 2014;15(6):717–28.
Jiang G, Dandekar S. Targeting NF-kappaB signaling with protein kinase C agonists as an emerging strategy for combating HIV latency. AIDS Res Hum Retroviruses. 2015;31(1):4–12.
O’Neil BH, Buzkova P, Farrah H, Kashatus D, Sanoff H, Goldberg RM, Baldwin AS, Funkhouser WK. Expression of nuclear factor-kappaB family proteins in hepatocellular carcinomas. Oncology. 2007;72(1–2):97–104.
Sookoian S, Pirola CJ. Systems biology elucidates common pathogenic mechanisms between nonalcoholic and alcoholic-fatty liver disease. PLoS ONE. 2013;8(3):e58895.
Wang M, Yang W, Chen Y, Wang J, Tan J, Qiao W. Cellular RelB interacts with the transactivator Tat and enhance HIV-1 expression. Retrovirology. 2018;15(1):65.
Zhou D, Huang W, Wei J, Zhang J, Liu Z, Ji R, Ge S, Xiao M, Fan Y, Lu C. RelB promotes liver fibrosis via inducing the release of injury-associated inflammatory cytokines. J Cell Mol Med. 2020;24(11):6008–14.
Elssner C, Goeppert B, Longerich T, Scherr AL, Stindt J, Nanduri LK, Rupp C, Kather JN, Schmitt N, Kautz N, et al. Nuclear translocation of RELB is increased in diseased human liver and promotes ductular reaction and biliary fibrosis in mice. Gastroenterology. 2019;156(4):1190–205.
Hettinghouse A, Liu R, Liu CJ. Multifunctional molecule ERp57: From cancer to neurodegenerative diseases. Pharmacol Ther. 2018;181:34–48.
Zhao Q, Feng Y, Jia X, Yin L, Zheng Y, Ouyang D, Zhou H, Zhang L. Proteome analysis of hepatic non-parenchymal cells of immune liver fibrosis rats. Sci China Life Sci. 2014;57(3):303–14.
Kondo R, Ishino K, Wada R, Takata H, Peng WX, Kudo M, Kure S, Kaneya Y, Taniai N, Yoshida H, et al. Downregulation of protein disulfideisomerase A3 expression inhibits cell proliferation and induces apoptosis through STAT3 signaling in hepatocellular carcinoma. Int J Oncol. 2019;54(4):1409–21.
Ko E, Kim JS, Ju S, Seo HW, Chang Y, Kang JA, Park SG, Jung G. Oxidatively modified protein-disulfide isomerase-associated 3 promotes dyskerin pseudouridine synthase 1-mediated malignancy and survival of hepatocellular carcinoma cells. Hepatology. 2018;68(5):1851–64.
Cheng D, Zhao L, Zhang L, Jiang Y, Tian Y, Xiao X, Gong G. p53 controls hepatitis C virus non-structural protein 5A-mediated downregulation of GADD45α expression via the NF-κB and PI3K-Akt pathways. J Gen Virol. 2013;94(Pt 2):326–35.
Mitchell JK, Midkiff BR, Israelow B, Evans MJ, Lanford RE, Walker CM, Lemon SM, McGivern DR. Hepatitis C Virus Indirectly Disrupts DNA Damage-Induced p53 Responses by Activating Protein Kinase R. mBio. 2017;8:2.
Jin Y, Lee WY, Toh ST, Tennakoon C, Toh HC, Chow PK, Chung AY, Chong SS, Ooi LL, Sung WK, et al. Comprehensive analysis of transcriptome profiles in hepatocellular carcinoma. J Transl Med. 2019;17(1):273.
Tannapfel A, Wasner M, Krause K, Geissler F, Katalinic A, Hauss J, Mossner J, Engeland K, Wittekind C. Expression of p73 and its relation to histopathology and prognosis in hepatocellular carcinoma. J Natl Cancer Inst. 1999;91(13):1154–8.
Bechmann LP, Gieseler RK, Sowa JP, Kahraman A, Erhard J, Wedemeyer I, Emons B, Jochum C, Feldkamp T, Gerken G, et al. Apoptosis is associated with CD36/fatty acid translocase upregulation in non-alcoholic steatohepatitis. Liver Int. 2010;30(6):850–9.
Chang SC, Choo WQ, Toh HC, Ding JL. SAG-UPS attenuates proapoptotic SARM and Noxa to confer survival advantage to early hepatocellular carcinoma. Cell Death Discov. 2015;1:15032.
Zhu H, Yang W, He LJ, Ding WJ, Zheng L, Liao SD, Huang P, Lu W, He QJ, Yang B. Upregulating Noxa by ER stress, celastrol exerts synergistic anti-cancer activity in combination with ABT-737 in human hepatocellular carcinoma cells. PLoS ONE. 2012;7(12):e52333.
Schulien I, Hockenjos B, Schmitt-Graeff A, Perdekamp MG, Follo M, Thimme R, Hasselblatt P. The transcription factor c-Jun/AP-1 promotes liver fibrosis during non-alcoholic steatohepatitis by regulating Osteopontin expression. Cell Death Differ. 2019;26(9):1688–99.
DeAngelis RA, Markiewski MM, Taub R, Lambris JD. A high-fat diet impairs liver regeneration in C57BL/6 mice through overexpression of the NF-kappaB inhibitor IkappaBalpha. Hepatology. 2005;42(5):1148–57.
Luedde T, Schwabe RF. NF-kappaB in the liver–linking injury, fibrosis and hepatocellular carcinoma. Nat Rev Gastroenterol Hepatol. 2011;8(2):108–18.
Taura M, Kudo E, Kariya R, Goto H, Matsuda K, Hattori S, Vaeteewoottacharn K, McDonald F, Suico MA, Shuto T, et al. COMMD1/Murr1 reinforces HIV-1 latent infection through IkappaB-alpha stabilization. J Virol. 2015;89(5):2643–58.
Kogan M, Deshmane S, Sawaya BE, Gracely EJ, Khalili K, Rappaport J. Inhibition of NF-kappaB activity by HIV-1 Vpr is dependent on Vpr binding protein. J Cell Physiol. 2013;228(4):781–90.
Vallejo-Diaz J, Chagoyen M, Olazabal-Moran M, Gonzalez-Garcia A, Carrera AC. The opposing roles of PIK3R1/p85alpha and PIK3R2/p85beta in cancer. Trends Cancer. 2019;5(4):233–44.
Muhanna N, Doron S, Wald O, Horani A, Eid A, Pappo O, Friedman SL, Safadi R. Activation of hepatic stellate cells after phagocytosis of lymphocytes: a novel pathway of fibrogenesis. Hepatology. 2008;48(3):963–77.
Wang R, Wang X, Zhuang L. Gene expression profiling reveals key genes and pathways related to the development of non-alcoholic fatty liver disease. Ann Hepatol. 2016;15(2):190–9.
Grise F, Bidaud A, Moreau V. Rho GTPases in hepatocellular carcinoma. Biochim Biophys Acta. 2009;1795(2):137–51.
Wang H, Lafdil F, Kong X, Gao B. Signal transducer and activator of transcription 3 in liver diseases: a novel therapeutic target. Int J Biol Sci. 2011;7(5):536–50.
Handa P, Vemulakonda AL, Maliken BD, Morgan-Stevenson V, Nelson JE, Dhillon BK, Hennessey KA, Gupta R, Yeh MM, Kowdley KV. Differences in hepatic expression of iron, inflammation and stress-related genes in patients with nonalcoholic steatohepatitis. Ann Hepatol. 2017;16(1):77–85.
Lee C, Cheung ST. STAT3: An Emerging Therapeutic Target for Hepatocellular Carcinoma. Cancers (Basel). 2019;11:11.
Lee MS, Hanspers K, Barker CS, Korn AP, McCune JM. Gene expression profiles during human CD4+ T cell differentiation. Int Immunol. 2004;16(8):1109–24.
Helgeland H, Gabrielsen I, Akselsen H, Sundaram AYM, Flam ST, Lie BA. Transcriptome profiling of human thymic CD4+ and CD8+ T cells compared to primary peripheral T cells. BMC Genomics. 2020;21(1):350.
Nirei K, Kanda T, Nakamura H, Matsuoka S, Takayama T, Sugitani M, Moriyama M. Persistent hepatic inflammation plays a role in hepatocellular carcinoma after sustained virological response in patients with HCV Infection. Int J Med Sci. 2018;15(5):466–74.
Perez S, Kaspi A, Domovitz T, Davidovich A, Lavi-Itzkovitz A, Meirson T, Alison Holmes J, Dai CY, Huang CF, Chung RT, et al. Hepatitis C virus leaves an epigenetic signature post cure of infection by direct-acting antivirals. PLoS Genet. 2019;15(6):e1008181.
Hamdane N, Juhling F, Crouchet E, El Saghire H, Thumann C, Oudot MA, Bandiera S, Saviano A, Ponsolles C, Roca Suarez AA, et al. HCV-induced epigenetic changes associated with liver cancer risk persist after sustained virologic response. Gastroenterology. 2019;156(8):2313–29.
Carrero A, Berenguer J, Hontanon V, Navarro J, Hernandez-Quero J, Galindo MJ, Quereda C, Santos I, Tellez MJ, Ortega E, et al. Effects of Eradication of HCV on cardiovascular risk and preclinical atherosclerosis in HIV/HCV-coinfected patients. J Acquir Immune Defic Syndr. 2020;83(3):292–300.
Wandrer F, Falk CS, John K, Skawran B, Manns MP, Schulze-Osthoff K, Bantel H. Interferon-mediated cytokine induction determines sustained virus control in chronic hepatitis C virus infection. J Infect Dis. 2016;213(5):746–54.
Irshad M, Gupta P, Irshad K. Immunopathogenesis of liver injury during hepatitis C virus infection. Viral Immunol. 2019;32(3):112–20.
Acknowledgements
We want to acknowledge the patients in this study for their participation and to the Spanish HIV HGM BioBank integrated into the Spanish AIDS Research Network (RIS) and collaborating Centers for the generous gifts of clinical samples used in this work. This study would not have been possible without the collaboration of all medical and nursing staff and data managers who have taken part in the project (see Additional file 7: Appendix S1, which shows all collaborators). We also want to thank the Bioinformatics Unit at the Institute of Health Carlos III for their valuable support for the bioinformatics analysis.
Funding
This study was supported by grants from Instituto de Salud Carlos III (ISCII; Grant Numbers PI20/00474 and PI17/00657 to JB, PI20/00507 and PI17/00903 to JGG, PI18CIII/00020 to AFR, and PI20CIII/00004 and PI17CIII/00003 to SR). The study was also funded by the RD16CIII/0002/0002 and RD16/0025/0017, and RD16/0025/0018 projects as part of the Plan Nacional R + D + I and co-funded by ISCIII- Subdirección General de Evaluación and the Fondo Europeo de Desarrollo Regional (FEDER). JB is an investigator from the Programa de Intensificación de la Actividad Investigadora en el Sistema Nacional de Salud (I3SNS), Ref. INT16/00100.
Author information
Authors and Affiliations
Consortia
Contributions
Conceptualization: SR, AFR. Data curation: JB, JGG, AC, VH, JN, JMG. Formal analysis: OBK, SR, IM, AFR. Funding acquisition: JB, JGG, and SR. Investigation and methodology: LMM, OBK, and AFR. Project Administration: JB. Supervision and visualization: SR. Writing – original draft preparation: IM, AFR, and SR. Writing – Review & Editing: MAJS. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The study was conducted following the Declaration of Helsinki, and patients gave their written consent. The Institutional Review Board and the Research Ethics Committee of the Instituto de Salud Carlos III (ISCIII) approved the study (CEI PI 23_2011).
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
: Table S1. Bioinformatics pipeline to analyze raw sequences from RNA-seq of HIV/HCV-infected and HIV-monoinfected patients.
Additional file 2
: Table S2. List of 4723 genes related to the immune system according to InnateDB, a public database of human immune genes (www.innatedb.com).
Additional file 3
: Table S3. Summary of significant differentially expressed genes (FC ≥1.5; FDR ≤0.05) in peripheral blood mononuclear cells: A) HIV/HCV-b versus HIV-mono. B) HIV/HCV-coinfected at week 24 (HIV/HCV-f) after SVR versus at baseline (HIV/HCV-b). C) HIV/HCV-f versus HIV-mono.
Additional file 4
: Table S4. Summary of significant KEGG pathways (FDR ≤0.05) in the HIV/HCV-f versus HIV/HCV-b comparison.
Additional file 5
: Table S5. Summary of significant KEGG pathways (FDR ≤0.05) in the HIV/HCV-f versus HIV-mono comparison.
Additional file 6
: Table S6. Summary of the most deregulated SDE genes (absolute FC ≥2 and FDR ≤0.05).
Additional file 7
: Appendix S1. The GESIDA 3603b Cohort Study Group.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Brochado, Ó., Martínez, I., Berenguer, J. et al. HCV eradication with IFN-based therapy does not completely restore gene expression in PBMCs from HIV/HCV-coinfected patients. J Biomed Sci 28, 23 (2021). https://doi.org/10.1186/s12929-021-00718-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12929-021-00718-6
Keywords
- HIV/HCV coinfection
- Gene expression
- Immune system
- HCV clearance
- Interferon therapy
- PBMCs