# Hepatitis C virus alters the morphology and function of peroxisomes

Esther Martin de Fourchambault, Nathalie Callens, Jean-Michel Saliou, Marie Fourcot, Oceane Delos, Nicolas Barois, Quentin Thorel, Santseharay Ramirez, Jens Bukh, Laurence Cocquerel, et al.

► To cite this version:

Esther Martin de Fourchambault, Nathalie Callens, Jean-Michel Saliou, Marie Fourcot, Oceane Delos, et al.. Hepatitis C virus alters the morphology and function of peroxisomes. Frontiers in Microbiology, 2023, 14, pp.1254728. 10.3389/fmicb.2023.1254728. hal-04245832

HAL Id: hal-04245832

https://hal.science/hal-04245832v1

Submitted on 16 Nov 2023

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

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

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

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

Check for updates

OPEN ACCESS

EDITED BY

Leiliang Zhang,
Shandong First Medical University and
Shandong Academy of Medical Sciences, China

REVIEWED BY

Michael Schindler,
University of Tübingen, Germany
Kai Xu,
Nanjing Normal University, China
Masahiko Ito,
Hamamatsu University School of Medicine,
Japan

*CORRESPONDENCE

Yves Rouillé

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

yves.rouille@ibl.cnrs.fr

RECEIVED 07 July 2023

ACCEPTED 07 September 2023

PUBLISHED 21 September 2023

CITATION

Martin de Fourchambault E, Callens N, Saliou J-M, Fourcot M, Delos O, Barois N, Thorel Q, Ramirez S, Bukh J, Cocquerel L, Bertrand-Michel J, Marot G, Sebti Y, Dubuisson J and Rouillé Y (2023) Hepatitis C virus alters the morphology and function of peroxisomes. *Front. Microbiol.* 14:1254728. doi: 10.3389/fmicb.2023.1254728

COPYRIGHT

© 2023 Martin de Fourchambault, Callens, Saliou, Fourcot, Delos, Barois, Thorel, Ramirez, Bukh, Cocquerel, Bertrand-Michel, Marot, Sebti, Dubuisson and Rouillé. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

# Hepatitis C virus alters the morphology and function of peroxisomes

Esther Martin de Fourchambault<fnref n="1" />, Nathalie Callens<fnref n="1" />,
Jean-Michel Saliou², Marie Fourcot², Oceane Delos³,⁴,
Nicolas Barois<fnref n="1" />,², Quentin Thorel⁵, Santseharay Ramirez⁶,
Jens Bukh⁶, Laurence Cocquerel<fnref n="1" />, Justine Bertrand-Michel³,⁴,
Guillemette Marot⁷, Yasmine Sebti⁵, Jean Dubuisson<fnref n="1" /> and
Yves Rouillé<fnref n="1" />*

Despite the introduction of effective treatments for hepatitis C in clinics, issues remain regarding the liver disease induced by chronic hepatitis C virus (HCV) infection. HCV is known to disturb the metabolism of infected cells, especially lipid metabolism and redox balance, but the mechanisms leading to HCV-induced pathogenesis are still poorly understood. In an APEX2-based proximity biotinylation screen, we identified ACBD5, a peroxisome membrane protein, as located in the vicinity of HCV replication complexes. Confocal microscopy confirmed the relocation of peroxisomes near HCV replication complexes and indicated that their morphology and number are altered in approximately 30% of infected Huh-7 cells. Peroxisomes are small versatile organelles involved among other functions in lipid metabolism and ROS regulation. To determine their importance in the HCV life cycle, we generated Huh-7 cells devoid of peroxisomes by inactivating the PEX5 and PEX3 genes using CRISPR/Cas9 and found that the absence of peroxisomes had no impact on replication kinetics or infectious titers of HCV strains JFH1 and DBN3a. The impact of HCV on peroxisomal functions was assessed using sub-genomic replicons. An increase of ROS was measured in peroxisomes of replicon-containing cells, correlated with a significant decrease of catalase activity with the DBN3a strain. In contrast, HCV replication had little to no impact on cytoplasmic and mitochondrial ROS, suggesting that the redox balance of peroxisomes is specifically impaired in cells replicating HCV. Our study provides evidence that peroxisome function and morphology are altered in HCV-infected cells.

KEYWORDS

hepatitis C virus, HCV genotype, peroxisome, proximity biotinylation, APEX2, CRISPR-Cas9, ROS

<footnote n="1">¹Université de Lille, CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, U 1019 – UMR9017 – CIIL – Center for Infection and Immunity of Lille, Lille, France, ²Université de Lille, CNRS, Inserm, CHU Lille, Institut Pasteur de Lille, UAR CNRS 2014 - US Inserm 41 - PLBS, Lille, France, ³MetaToul-MetaboHUB, National Infrastructure of Metabolomics and Fluxomics, Toulouse, France, ⁴I2MC, Université de Toulouse, Inserm, Université Toulouse III – Paul Sabatier (UPS), Toulouse, France, ⁵Université de Lille, Inserm, CHU Lille, Institut Pasteur de Lille, U1011 - EGID, Lille, France, ⁶Faculty of Health and Medical Sciences, Copenhagen Hepatitis C Program (CO-HEP), Department of Infectious Diseases, Copenhagen University Hospital Hvidovre and Department of Immunology and Microbiology, University of Copenhagen, Copenhagen, Denmark, ⁷Université de Lille, Inria, CHU Lille, ULR 2694 - METRICS: Évaluation des technologies de santé et des pratiques médicales, Lille, France</footnote>
# 1. Introduction

Hepatitis C is one of the leading causes of viral chronic hepatitis worldwide. According to WHO estimates, some 58 million people are chronically infected with the hepatitis C virus (HCV), causing about 0.3 million deaths every year (WHO, 2022). About 60-80% of people infected with HCV will develop a chronic infection and are at risk of progressive liver disease, including fibrosis, steatosis, cirrhosis and hepatocellular carcinoma. The infection remains asymptomatic until decades after the initial contact with the virus, when symptoms of liver damage appear. About 10 years ago, efficient antiviral treatments were introduced in the clinics, which cure about 95% of treated patients from their infection. Nevertheless, the risk of HCV-associated disease is reduced but not eradicated among patients with fibrosis who achieve sustained virological response following DAA treatment (van der Meer et al., 2017). Moreover, access to diagnosis remains limited, especially in low-income countries, and a vast majority of HCV-infected persons are unaware of their infected status, therefore HCV keeps spreading. Indeed, it is estimated that 1.5 million new infections occur each year worldwide (WHO, 2022). These data indicate that hepatitis C remains a global health problem, with a severe impact on the quality of life of chronically infected people.

HCV is a small enveloped RNA virus that infects human hepatocytes. Its genome encodes a single open reading frame. It is translated into a precursor polyprotein, which is cleaved into 10 mature structural and non-structural (NS) proteins by host and viral proteases (Lindenbach et al., 2007). All HCV proteins are initially inserted in the membrane of the endoplasmic reticulum (ER). A subset of these proteins (NS3/4A, NS4B, NS5A and NS5B), which are involved in the replication of the genomic RNA, induce membrane rearrangements, including single and double membrane vesicles collectively known as 'membranous web', in which the HCV genome is replicated (Ferraris et al., 2010; Romero-Brey et al., 2012). The membranous web likely originates from the ER membrane and includes proteins and lipids normally found in the ER or in other sub-cellular compartments in non-infected cells (Tu et al., 1999; Hamamoto et al., 2005; Stone et al., 2007; Tai et al., 2009; Reiss et al., 2011; Khan et al., 2014; Wang et al., 2014). In addition to these morphological changes, HCV infection alters different cellular pathways, notably metabolic pathways such as glucose metabolism, lipid metabolism and oxidative stress (Hui et al., 2003; Negro, 2010; Virzi et al., 2018). HCV infection also induces neutral lipid accumulation in lipid droplets (Filipe and McLauchlan, 2015; Galli et al., 2021). Recently peroxisome defects were also documented in cell culture and patient samples (Lupberger et al., 2019). As a consequence of these metabolic alterations, HCV infection is often associated with high prevalence of hepatic metabolic disorders such as insulin resistance and/or steatosis (Negro, 2010). HCV-induced alterations in cellular metabolism likely originate in interactions engaged between viral proteins and cellular partners. The HCV genome encodes only a limited number of proteins and relies on host factors for many steps of its life cycle. During HCV replication, proviral host factors are recruited to replication complexes by viral proteins (Reiss et al., 2011; Wang et al., 2014). The recruitment of these proteins may in turn alter their normal cellular functions. Thus, a better understanding of these interactions is crucial to deciphering the alterations in cellular metabolism caused by HCV infection, which are at the origin of pathophysiological disorders.

In order to get more insight into cellular proteins that are recruited in the vicinity of HCV replication complexes, we performed a proximity biotinylation screen based on the use of APEX2 (Rhee et al., 2013). Interestingly, we found that peroxisomes are recruited in the vicinity of HCV replication complexes and that their morphology and function are altered in HCV infected cells. Peroxisomes are small membrane-bound organelles with unique functions in the cell metabolism. Notably, they are involved in lipid metabolism and in the regulation of reactive oxygen species (ROS), two pathways that are critical for HCV infection.

# 2. Materials and methods

## 2.1. Chemicals

Dulbecco's modified Eagle's medium (DMEM), Dulbecco's phosphate-buffered saline (D-PBS), OptiMEM and geneticin were purchased from ThermoFischer Scientific. 4,6-Diamidino-2-phenylindole (DAPI) was from Molecular Probes. Mowiol 3-88 was from Calbiochem. TransIT-LT1 was from Mirus Bio. All other chemicals were from Sigma. Biotin-phenol (BP) was synthesized by mixing equimolar amounts of tyramine (from Sigma-Aldrich) and biotin-OSu (D-Biotin N-hydroxysuccinimide ester, from Iris Biotech) in anhydrous dimethylformamide. The mixture was incubated for 2 h at room temperature. BP synthesis was verified by LC-MS, purified by HPLC and lyophilized.

## 2.2. Antibodies

Mouse monoclonal antibodies (mAb) 9E10 against HCV NS5A (Lindenbach et al., 2005) and 486D39 against NS3 were kindly provided by C. M. Rice (The Rockefeller University) and by J.-F. Delagneau (Bio-Rad), respectively. Rabbit anti-PEX14 and anti-PEX5 antibodies were from ProteinTech. Rabbit anti-ACBD5 and anti-PMP-70 antibodies and mouse mAb TUB2.1 anti-β-tubulin were from Sigma. Anti-catalase rabbit mAb was from Cell Signaling Technologies. Mouse mAb14F3.2 anti-PEX3 was from Millipore. Mouse mAb anti-TOM20 was from BD Biosciences. Mouse mAb anti-myc was produced *in vitro* using a MiniPerm apparatus (Heraeus) following the manufacturer's protocol. All HRP-conjugated and fluorescent A488-, Cy3-, or Cy5-conjugated secondary antibodies and HRP- or Alexa 488-conjugated streptavidin were from Jackson Immunoresearch.

## 2.3. Plasmids

The plasmid pEGFP-ACBD5 (Costello et al., 2017a) was kindly provided by M. Schrader (University of Exeter). The plasmid pcDNA-myc-CD2AP (Monzo et al., 2005) was kindly provided by M. Cormont (Université Côte d'Azur). Plasmids expressing roGFP2 constructs pMF1707, pMF1706 and pMF1762 (Ivashchenko et al., 2011) were kindly provided by M. Fransen (Katholieke Universiteit Leuven) through Addgene. The plasmid pX330-U6-Chimeric_BB-CBh-hSpCas9 (Cong et al., 2013) was a gift from F. Zhang (Addgene plasmid # 42230). To construct the plasmid pmCherry-SKL, synthetic
oligonucleotides 5'-GATCTAAGCTTTG-3' and 5'-AATTCAAAGCTTA-3' were phosphorylated, annealed and inserted between BglII and EcoRI sites of the plasmid pmCherry-C1 (Clonetech).

## 2.4. Cell culture

Huh-7 (Nakabayashi et al., 1982) and CD81-deficient Huh-7w7 (Rocha-Perugini et al., 2009) cells were cultured at 37°C with 5% CO₂ in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% heat-inactivated fetal bovine serum (FBS) and glutamax.

## 2.5. Viruses

The viruses used in this study were the original genotype 2a JFH1 strain (Wakita et al., 2005), or a derivative of this strain bearing titer-enhancing mutations and the reconstitution of the A4 epitope in E1 (Goueslain et al., 2010), and the cell culture-adapted genotype 3a DBN3a strain (Ramirez et al., 2016). To generate genomic HCV RNA, plasmids were linearized by XbaI, treated with mung bean nuclease, and used as templates for in vitro transcription with the MEGAScript kit from Ambion. In vitro-transcribed RNA was delivered to Huh-7 cells by electroporation as previously described (Kato et al., 2003). Viral stocks were obtained as previously reported (Delgrange et al., 2007) and infectious titers were measured by the TCID50 method. For kinetics analyses, infected cells were trypsinized the day before, counted and 1.5×10⁵ cells were dispensed in a P-24 well in 1 mL of medium. The infectious titer was measured after 24 h of culture. For experiments with PEX3-KO and PEX5-KO Huh-7 cells, a Mann-Whitney test was used to determine statistical differences in infectious titers against the Huh-7 control group.

## 2.6. Replicons

The plasmid pSGR-JFH1 encoding a subgenomic replicon (SGR) of the JFH1 strain was obtained from T. Wakita (Kato et al., 2003). The plasmid pSGR-JFH1-GFP encoding an SGR containing the coding sequence of EGFP in domain III of NS5A has been described (Sahuc et al., 2019). The plasmid pSGR-JFH1-APEX2-A40 encoding an SGR containing the coding sequence of APEX2 in domain III and a deletion in domain II of NS5A has been described (Riva et al., 2021). The plasmid pSGR-JFH1-Rluc expressing a Renilla luciferase (Rluc) has been previously described (Riva et al., 2021). To construct an Rluc-expressing DBN3a replicon, the 3'-UTR of DBN3a was amplified by PCR using primer pair 5'-GCAGGAATTCTAATACGA CTCACTATAGCCTGCC-3' and 5'-TCATACGCGTTGG GCGGCGGATGGTGTTTCTTTTTGG-3' and inserted into EcoRI and MluI sites of pSGR-JFH1-Rluc, yielding the plasmid pSGR-JFH1-Rluc/3'UTR3a. Then, a DNA fragment containing the EMCV IRES and the coding sequence of the N-terminal part of DBN3a NS3 with no KpnI site was assembled by tripartite fusion PCR using primer pairs 5'-CCAGAAGGTACCCCATGTATGGG-3' / 5'-TGTGATCGGGG CCATGGTATCATCGTGTTTTTC-3', 5'-GATGATACCATGGGCC CCGATCACAGCATACACCC-3' / 5'-ATGAAGGTAgCCGAC TTGATAGCT-3', and 5'-AGCTATCAAGTCGGcTACCTTCAT-3' /

5'-CTGCTCTAGAGCACTCACAGAGAACAACCGAG-3' and inserted into KpnI and BsiWI sites of pSGR-JFH1-Rluc/3'UTR3a, yielding the plasmid pSGR-JFH1/3a-Rluc/3'UTR3a. Finally, the BsiWI-XbaI fragment of DBN3a was inserted into the same sites of pSGR-JFH1/3a-Rluc/3'UTR3a, yielding the plasmid pSGR-DBN3a-Rluc. The GND mutation was introduced by fusion PCR using primer pairs 5'-AGGTCGACTCTAGACATGATCTGC-3' / 5'-TTTCTTGT CTGCGGAaATGATCTG-3' and 5'-CACCACGACCAGATCATiT CCGCA-3' / 5'-ATGTGGACCTCAAAGAAGACCCCC-3'. The PCR product was digested with NotI and XbaI and inserted at the corresponding sites in pSGR-DBN3a-Rluc. To construct a DBN3a-based SGR expressing a neomycin resistance marker, the coding sequence of neomycin phosphotransferase was obtained by PCR using the primers 5'-GGGAGAGGGTTTAAACTCAGAAG-3' and 5'-CCCAACGCGTATGATTGAACAAGATGGATTG-3', and the PCR product was digested by PmeI and MluI and inserted at the corresponding sites in pSGR-DBN3a-Rluc, yielding the plasmid pSGR-DBN3a. All constructs were verified by sequencing. In vitro transcribed SGR RNA was introduced in Huh-7 cells by electroporation (Kato et al., 2006). Electroporated cells were kept in culture for up to 4 days post electroporation for Rluc-expressing SGRs, or were selected using 0.5 mg/mL geneticin to generate Huh-7 cell lines permanently expressing JFH1 or DBN3a SGR.

## 2.7. APEX2-mediated biotinylation

Biotinylation of control and SGR-APEX-containing cells were carried out as previously described (Hung et al., 2016). Briefly, cells were incubated for 30 min at 37°C/5% CO₂ in DMEM containing 0.5 mM BP. H₂O₂ was then added to a final concentration of 1 mM and the cells were incubated for 1 min at room temperature with gentle agitation. The reaction was quenched, and the cells were washed 3 times with D-PBS containing 5 mM Trolox, 10 mM sodium ascorbate and 10 mM NaN₃ before proceeding to immunoblotting or immunofluorescence.

For purification of biotinylated proteins, Huh-7 cells with or without SGR-APEX were biotinylated as described above. Each cell line was treated with or without H₂O₂ after BP incubation. Cell pellets were then lysed in RIPA buffer (50 mM TrisCl pH 7.5, 150 mM NaCl, 0.1% SDS, 0.5% Na deoxycholate, 1% Triton X-100) containing 5 mM Trolox, 10 mM sodium ascorbate and 10 mM NaN₃ and the lysates were centrifuged at 14,000 x g. Streptavidin-coated magnetic beads (Pierce) were washed twice with RIPA buffer, and 2.5 mg of each sample was incubated with 100 µL of magnetic beads slurry with rotation for 2 h at room temperature. The beads were subsequently washed 3 times with RIPA lysis buffer, twice with 1 M KCl, once with 0.1 M Na₂CO₃, once with 2 M urea in 10 mM Tris-HCl (pH 8.0), and twice with RIPA lysis buffer. Biotinylated proteins were then eluted from the beads at 95°C for 10 min in 100 µL of 2 X protein loading buffer containing 40 mM DTT and 2 mM biotin.

## 2.8. Generation of knockout cells

Gene invalidation using CRISPR-Cas9 was performed as described (Ran et al., 2013). Briefly, guide RNA coding sequences were ordered as synthetic oligonucleotides, phosphorylated with T4-polynucleotide kinase, annealed and inserted between BbsI sites
of plasmid pX330 (Ran et al., 2013). The coding sequences of the guide RNAs were CTCCAGCTGCAAGAAACTCC, AGATGCTGTGG ATGTAACTC and CTTGTAGGGGTACCAGTTTG to generate PEX5-KO cells lines and CACCTCCAAGGACCGTGCCC, TAAAGATGCTGAGGTCTGTA and CATATTTCTAGCTTGTTTGA to generate PEX3-KO cells lines. Transfection of Huh-7 cells and selection of KO cell lines were performed as previously described (Ferlin et al., 2018).

## 2.9. Immunofluorescence

Cells were processed for immunofluorescent detection of viral and cellular proteins as previously described (Rouillé et al., 2006). Coverslips were mounted on glass slides using a Mowiol-based medium. For colocalization experiments, confocal microscopy was performed with an LSM880 confocal microscope (Zeiss) using a 63×/1.4 numerical aperture oil immersion objective. Signals were sequentially collected by using single fluorescence excitation and acquisition settings to avoid crossover. For some experiments, high-resolution images (pixel 0.035×0.035 µm²) were acquired using an Airyscan detector used in SR mode and processed for high-resolution reconstruction using Zen software version 2 (Carl Zeiss Microscopy). Images were assembled by using Adobe Photoshop software. The Pearson's correlation coefficients (PCC) were calculated using the JACoP plugin of ImageJ software. For each calculation, at least 30 cells were analyzed.

## 2.10. Peroxisome morphology assessment

Huh-7 cells were infected with HCV strains JFH1 or cell culture adapted DBN3a (MOI=0.1). Infected cells were kept in culture and when they were passaged, approximately 5×10⁴ cells were seeded on glass coverslips. The day after, cells were fixed with 3% PFA, transferred out of the BSL3 laboratory and immunostained as described above using antibodies to NS5A and Pex14. Image stacks of cells of interest were acquired using an inverted Eclipse Ti inverted confocal microscope (Nikon) equipped with a CSU-W1 spinning-disk (Yokogawa, Roper Scientific). A live-SR module (Gataca Systems) was added to the system to improve the obtained resolutions. Observations were done with a 60× oil immersion objective (Nikon Plan Apo 60× NA 1.4). Signals were sequentially collected by using single fluorescence excitation and acquisition settings to avoid crossover. Image stacks were then processed using the Imaris software to model a 3D-image of peroxisomes. Objects smaller than 0.1 µm were filtered out and the volume of each object was measured, enabling the volume of each peroxisome and the average number of peroxisomes per cell to be calculated. For each measurement, 30 cells from 3 independent infections were analyzed. Statistical differences were determined with a Kruskall-Wallis test.

## 2.11. ROS measurements

Plasmids expressing roGFP2 sensors localized in the cytosol, mitochondria or peroxisomes were transfected into control, SGR2a- and SGR3a-containing Huh-7 cells grown in glass-bottom 4 well µ-slide (Ibidi). At 18 h post-transfection, cells were rinsed three times

with sterile PBS and placed in phenol red-free DMEM. As a positive control, H₂O₂ was added to naïve Huh-7 cells at a final concentration of 0.4 mM. Two hours later, cells were imaged at 37°C and 5% CO₂ with an LSM880 confocal microscope (Zeiss) using a 63×/1.4 numerical aperture oil immersion objective. Once the acquisition conditions had been set for each roGFP2 sensor, all cells expressing this sensor were acquired using the same parameters. Images were acquired by sequentially exciting with 405 nm and a 488 nm lasers, and then processed using ImageJ software. A region of interest (ROI) was drawn around the whole cell and for both wavelengths the mean intensity of fluorescence was measured. Ratios of fluorescence intensity obtained with each excitation wavelength (405/488) were then calculated. A Kruskall-Wallis test was used to determine statistical differences.

## 2.12. Immunoblotting

Cells were lysed in 50 mM Tris-Cl buffer (pH 7.5) containing 100 mM NaCl, 1 mM EDTA, 1% Triton X-100, 0.1% sodium dodecyl sulfate (SDS), and protease inhibitors for 20 min on ice. Cells were collected, and the nuclei were pelleted. Protein concentration in post-nuclear supernatants was determined by the bicinchoninic acid (BCA) method as recommended by the manufacturer (Sigma), using bovine serum albumin as a standard. Proteins were separated by SDS-polyacrylamide gel electrophoresis and transferred to nitrocellulose membranes (Hybond-ECL; Amersham) by using a Trans-Blot apparatus (Bio-Rad). Blocking was carried out in PBS containing 5% dried skim milk and 0.1% NP-40. Proteins of interest were revealed with specific primary antibodies, followed by species-specific secondary antibodies conjugated to HRP, and enhanced chemiluminescence detection (Super Signal West Pico), as recommended by the manufacturer (ThermoFischer Scientific). The signals were recorded using a LAS 3000 apparatus (Fujifilm). For quantification, unsaturated signals were measured using the gel quantification function of ImageJ (Fiji).

## 2.13. Catalase activity

Catalase activity was measured using a colorimetric catalase assay kit, as recommended by the supplier (Abcam), on cell lysates and normalized by protein content, measured by the BCA method (Sigma).

## 2.14. Lipidomics

For extraction, cells were solubilized in 1 mL of Water/EGTA 5 mM: Methanol (1:2, v/v). Lipids corresponding to total fatty acids in cells were extracted according to Bligh and Dyer (1959) in dichloromethane/water/methanol 2% acetic acid (2.5:2.5:2, v/v/v), in the presence of the internal standards glyceryl trinonadecanoate (4 µg). Extracts were centrifuged during 6 min at 2500 rpm. The lipid extract was hydrolyzed in KOH (0.5 M in methanol, 1 mL) at 55°C for 30 min. A second extraction was realized with 1.5 mL of Methanol, 2 mL of water, and 2.5 mL of dichloromethane. The lipids extract was dried under nitrogen and derivatized with PentafluoroBenzylBromide
in ACN (1%) and Diisopropylethylamine in ACN (1%) (1/1, v/v). Final extracts (1 µL) were analyzed on a gas chromatography Thermo Trace 1,310 with mass spectrometer Thermo TSQ 8000 EVO. Chemical ionization is used as ionization, helium as carrier gas, and a HP-5MS column (30 m x 0.25 mm, 0.25 µm) is used to separation. The oven is set at 180°C during 1 min. Then, the temperature is increased at 315°C at 6°C/min during 2 min. SIM method is used to detect the component (m/z 255 to C16:0, m/z 297 to internal standard and m/z 395 to C26:0) and the quantification was obtained through external calibration curves.

## 2.15. Mass spectrometry

APEX2-biotinylated proteins were loaded on SDS-PAGE to perform one gel slice trypsin digestion for each sample. Peptides were extracted with 0.1% formic acid in acetonitrile, evaporated to reduce volume at 8 µL and injected on an UltiMate 3,000 RSLCnano System (Thermo Fisher Scientific). Peptides were automatically fractionated onto a commercial C18 reversed-phase column (75 µm x 250 mm, 2-µm particle, PepMap100 RSLC column, Thermo Fisher Scientific (Waltham, USA), temperature 35°C). Trapping was performed for 4 min at 5 µL/min with solvent A (98% H₂O, 2% acetonitrile and 0.1% formic acid). Elution was performed using two solvents, A (0.1% formic acid in water) and B (0.1% formic acid in acetonitrile), at a flow rate of 300 nL/min. Gradient separation was 3 min at 3% B, 110 min from 3 to 20% B, 10 min from 20 to 80% B, and maintained for 15 min. The column was equilibrated for 6 min with 3% buffer B prior to the next sample analysis. The eluted peptides from the C18 column were analyzed by Q-Exactive instruments (Thermo Fisher Scientific, Waltham, USA). The electrospray voltage was 1.9 kV, and the capillary temperature was 275°C. Full mass spectrometry (MS) scans were acquired in the Orbitrap mass analyzer over the m/z 400–1,200 range with a resolution of 70,000 (m/z 200). The target value was 3.00E+06. The 15 most intense peaks with charge states between 2 and 5 were fragmented in the high-energy collisional dissociation cell (HCD) with a normalized collision energy of 27%, and the tandem mass spectrum was acquired using an Orbitrap mass analyzer with a resolution of 17,500 at m/z 200. The target value was 1.00E+05. The ion selection threshold was 5.0E+04 counts, and the maximum allowed ion accumulation times were 250 ms for full MS scans and 100 ms for tandem mass spectra. The dynamic exclusion was set to 30 s. The full proteome analysis was performed in same conditions with a column of 500 mm length, heated at 55°C and a 3 to 20% B gradient separation of 170 min. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (Perez-Riverol et al., 2022) partner repository with the dataset identifier PXD044994 and 10.6019/PXD044994.

## 2.16. Proteomic data analysis

Raw data collected during nano LC-MS/MS analyses were processed and converted into an \*.mgf peak list format with Proteome Discoverer 1.4 (Thermo Fisher Scientific). MS/MS data were analyzed using the search engine Mascot (version 2.4.0, Matrix Science, London, United Kingdom) installed on a local server. Searches were performed with a tolerance on mass measurement of 0.02 Da for

precursor and 10 ppm for fragment ions against a composite target-decoy database built with a Human Swissprot database (taxonomy 9,606) fused with the sequences of recombinant trypsin, NS3, NS4A, NS4B, NS5A, NS5B, and a list of classical contaminants. Cysteine carbamidomethylation, methionine oxidation, protein N-terminal acetylation, and cysteine propionamidation were searched as variable modifications. Up to one missed trypsin cleavage was allowed.

For the APEX2-biotinylated protein samples, the identification results were imported into Scaffold software 4.5.3 and peptides were filtered out according to the cutoff set for proteins hits with 1 or more peptides taller than 9 residues, ion score > 20, identity score > -7, corresponding to a 1% false positive rate.

For the full proteome analysis, the identification results were imported into Proline software<fnref n="1" /> (Bouyssié et al., 2020). Peptide spectrum matches taller than 9 residues and ion scores >10 were retained. The false discovery rate was then optimized to be below 1% at the protein level using the Mascot Modified Mudpit score. Extracted Ion Current based quantification was performed with Proline 2.0. At peptide level, a single charge level per peptide was selected. Then, for peptides with 0 spectra, those with abundance above the upper bound of the 99% confidence interval of the distribution were filtered out. The abundance of a protein was obtained by summing the abundances of its component peptides.

Statistical analysis was performed by using R (version 4.1.2). Data were log2 transformed and distributions were centered on the upper quartile. Only proteins with at least three measurements (non missing values) in at least one group were kept. Following recommendations of Tyanova et al. (2016), missing values were then imputed using a normal distribution X ~ (μ-1.8σ, 0.3σ). Differential expression analysis between control Huh-7 and SGR-JFH1 and between control Huh-7 and SGR-DBN3a was performed using limma R package (Ritchie et al., 2015), which uses an empirical Bayesian approach to estimate variances in moderated t tests. Raw *p* values were adjusted for multiple testing using the Benjamini–Hochberg procedure (Benjamini and Hochberg, 1995). Proteins were considered differentially expressed if their adjusted *p*-values were below 0.05.

In order to evaluate the ability of peroxisomal proteins to discriminate groups when considered altogether, a multivariate analysis was performed with the mixOmics R package, using the partial least squares – discriminant approach (PLS-DA) (Rohart et al., 2017). Only the peroxisomal proteins were kept in this multivariate analysis.

# 3. Results

## 3.1. Peroxisomes are recruited close to HCV replication complexes

To identify cellular proteins associated to HCV replication complexes, we used an SGR of JFH1 strain containing the small peroxidase APEX2 inserted into domain III of NS5A (Riva et al., 2021). The catalytic activity of APEX2 was used to biotinylate proteins that are in close proximity (less than 20 nm) to the NS5A-APEX2 fusion protein

<footnote n="1">1 http://proline.profiproteomics.fr</footnote>
using BP as a substrate (Rhee et al., 2013). Cells were incubated with APEX2 substrates BP and H₂O₂ for 1 min and the presence of biotinylated material was verified by immunofluorescence and immunoblotting. SGR-APEX-containing cells displayed double-labeled dot-like structures when they were probed with Alexa488-conjugated streptavidin to reveal the localization of biotinylated material and with an anti-NS5A antibody as a marker of HCV replication complexes (Figure 1A). These double-labeled structures were found in SGR-APEX-containing, but not in non-tagged SGR-containing Huh-7 cells (Figure 1B). When the biotinylated proteins were analyzed by immunoblot, several specific bands were detected in SGR-containing cells incubated with APEX2 substrates (Figure 1C, left blot), in addition to a few bands present in control samples, which correspond to cellular proteins that are endogenously biotinylated. When the same samples were analyzed with an anti-NS5A antibody, the NS5A-APEX2 fusion protein was detected at a molecular size very similar to that of a major band of the biotinylated protein (Figure 1C, right blot), indicating that, as expected, NS5A-APEX2 was one of the most efficiently biotinylated proteins.

Biotinylated proteins were purified with streptavidin-coupled agarose beads. Streptavidin blot analysis of the eluted proteins revealed a large range of biotinylated proteins, whereas the controls again showed only endogenously biotinylated proteins (Figure 1D, left blot). The same membrane was then re-probed with an anti-NS3 antibody to verify that a known component of replication complexes was actually biotinylated. As expected, HCV NS3 was detected among the biotinylated proteins purified from the lysate of SGR-containing cells incubated in the presence of APEX2 substrates, but not under control conditions with no substrate and/or no peroxidase (Figure 1D, right blot). These results validate the use of this method to biotinylate and purify constituent proteins of HCV replication complexes.

Next, we analyzed the proteins purified after APEX2-mediated biotinylation in SGR-containing cells. Affinity purified proteins derived from biotinylated and non-biotinylated control and SGR-APEX-containing Huh-7 cells were subjected to mass spectrometry analysis. Liquid chromatography tandem-mass spectrometry (LC-MS/MS) was performed from in-gel digested samples and spectral counting was used to compare protein enrichment between biotinylated and non-biotinylated SGR-APEX samples. Forty-eight cellular proteins were identified with at least 2 unique peptides in the biotinylated SGR-APEX sample, no peptide in the non-biotinylated SGR-APEX sample and in both control Huh-7 samples (Table 1). Of these 48 proteins, CD2AP (Igloi et al., 2015; Zhang et al., 2018) and PACSIN2 (Nguyen et al., 2020) have previously been shown to interact with NS5A and Scribble with NS4B (Hu et al., 2016) in HCV-infected cells. Another host factor potentially recruited to HCV replication complexes is Ago2, which interacts with HCV RNA 5'-UTR in association with miR122 (Shimakami et al., 2012). In addition, 11 other proteins have been reported to interact with NS5A and 5 proteins with NS3, NS4B or NS5B in proteomics or yeast two-hybrid screens (de Chassey et al., 2008; Dolan et al., 2013; Germain et al., 2014; Ramage et al., 2015; Goonawardane et al., 2017), but their recruitment to HCV replication complexes has not been established. It should be noted that proteins identified with this approach include proteins that are very close to each other, without necessarily being interactors.

We verified by immunofluorescence microscopy the recruitment to HCV-replication complexes of four out of the top five proteins in Table 1, which correspond to proteins identified with a spectral ratio at least equal to 10. Tagged constructs of these proteins were expressed

in control and SGR-containing Huh-7 cells, and their intracellular localization was compared to that of HCV NS5A in SGR-containing cells. CD2-associated protein (CD2AP), a protein of the plasma membrane, was massively recruited to NS5A-positive dot-like structures (Supplementary Figure S1), as expected (Igloi et al., 2015; Zhang et al., 2018). The intracellular distribution of Acyl-CoA-binding domain-containing protein 5 (ACBD5, a peroxisomal membrane protein) was also altered in SGR-containing cells, as compared to control Huh-7 cells (Supplementary Figure S2). In contrast, no visible change of intracellular localization was observed for Protein kinase C and casein kinase substrate in neurons protein 2 (PACSIN2, also known as syndapin II), or for filamin B (data not shown). The reason for this apparent absence of recruitment is unclear, but it could represent proteins for which only a minor fraction of the cellular pool may associate with HCV replication complexes. Alternatively, they may represent proteins interacting with a pool of NS5A located in the ER membrane, outside the replication complexes.

Interestingly, in addition to ACBD5, another peroxisome membrane protein, PEX14, was present in the list of proteins identified in this screen (Table 1). This suggested that peroxisomes could be present in the vicinity of HCV replication complexes. To confirm that ACBD5 was actually biotinylated in cells with an APEX2-containing SGR, biotinylated proteins were analyzed by immunoblotting using an anti-ACBD5 antibody. ACBD5 was identified in the pool of proteins biotinylated by APEX2 in the presence of both of its substrates, but not in the absence of BP or if APEX2 was replaced by GFP in the replicon (Figure 2A). This result confirmed that ACBD5 is actually biotinylated in cells containing an APEX2-SGR, and suggested that it could be recruited in the proximity of HCV-replication complexes.

To test whether ACBD5 relocalized from peroxisomes to HCV replication complexes or whole peroxisomes were recruited to the proximity of HCV-replication complexes, we expressed a GFP-tagged version of ACBD5 in naïve or SGR-containing Huh-7 cells and labeled them with antibodies to NS5A, a marker of HCV replication complexes, and PEX14, a marker of peroxisomes. In both naïve and SGR-containing cells, GFP-ACBD5 was co-distributed with PEX14, indicating that it was correctly localized in peroxisomes in naïve (Figure 2B) and SGR-containing cells (Figures 2C,D). The localization of GFP-ACBD5 to peroxisomes was supported by Pearson correlation coefficients of 0.80±0.07 measured in SGR-containing cells and 0.81±0.07 in control Huh-7 cells. In SGR-containing cells, GFP-ACBD5 was also associated with NS5A-positive structures (Figures 2C,D), confirming the recruitment of peroxisomes to the proximity of HCV-replication complexes. It is worth noting that the proximity of peroxisomes and HCV-replication complexes was not induced by GFP-ACBD5 expression, because a colocalization of NS5A and PEX14 was also observed in GFP-transfected cells (Figure 2E).

## 3.2. Alteration of peroxisomes morphology in replicon-containing Huh-7 cells

Interestingly, not only the localization of peroxisomes but also their morphology appeared to be altered in a fraction of JFH1 SGR-containing cells expressing GFP-ACBD5 (Figure 2D). To assess if the change of peroxisome morphology was induced or not by GFP-ACBD5 expression in transfected cells, we analyzed it in untransfected SGR-containing cells using PEX14 as an endogenous
<figure><img src="image_4.png" /><figcaption>FIGURE 1</figcaption></figure>

(A,B) Immunofluorescence microscopy of biotinylated cells expressing the SGR-APEX (A) or non-tagged SGR (B) double labeled with A488-streptavidin (green) and with anti-NS5A antibody (red). Nuclei were stained with DAPI (blue). Bars, 20 µm. (C) Immunoblot analysis of biotinylated proteins. HCV SGR-APEX-containing and control Huh-7 cells pre-incubated in the presence of BP were treated in the presence (+) or in the absence (-) of H₂O₂ for 1 min, quenched and lysed. Cell lysates were probed with streptavidin-HRP (left side of the blot) or with an anti-NS5A antibody (right side of the blot). (D) Immunoblot analysis of proteins purified with strepatvidin beads. HCV SGR-APEX-containing and control Huh-7 cells were processed as described in (C). Biotinylated proteins were purified using streptavidin beads and analyzed by immunoblot with streptavidin-HRP (left side of the blot) and with an antibody to HCV NS3 (right side of the blot). Note that mammalian cells have four endogenously biotinylated proteins that are detected in control conditions.

marker. Again, a partial relocalization of peroxisomes in the vicinity of HCV replication complexes was observed (Figure 3B). When compared with peroxisomes of control cells (Figure 3A), larger peroxisomes were present in some but not all SGR-containing cells. Altered peroxisomal structures included both larger structures and long thin tubular structures. Cells displaying peroxisomes with an altered morphology were also observed in cells containing a sub-genomic replicon of strain DBN3a (data not shown). These altered peroxisome structures were not always associated with NS5A, suggesting that the alteration of their morphology does not directly depend on their association with HCV replication complexes. In other cells, near-normal sized peroxisomes colocalize with HCV replication complexes, again indicating that the morphological alteration of peroxisomes does not uniquely depend on their association with HCV replication complexes.

## 3.3. Association of peroxisomes and HCV-replication complexes

The association of peroxisomes and HCV replication complexes could also be observed using ACBD5 as an endogenous peroxisome marker (Supplementary Figure S3A) and in cells expressing mCherry-SKL, a red fluorescent peroxisomal marker (Supplementary Figure S3B). In any case, peroxisomes and HCV replication complexes often appeared juxtaposed rather than completely colocalized, when visualized at high resolution (Figure 3C). Pearson correlation coefficient values obtained when the co-distribution of ACBD5 and NS5A or NS3 was analyzed (Figure 3D) are consistent with a partial colocalization or a close proximity of peroxisome and HCV
<table><caption>TABLE 1 Proteins identified in the proximity biotinylation screen.</caption><thead><tr><th rowspan="2">Protein name</th><th rowspan="2">Accession number</th><th rowspan="2">Size (Da)</th><th colspan="4">Unique peptides</th><th rowspan="2">coverage</th></tr><tr><th colspan="2">SGR</th><th colspan="2">Hu</th></tr><tr><th>Sample</th><th></th><th></th><th>+</th><th>-</th><th>+</th><th>-</th><th></th></tr></thead><tbody><tr><td>Biotinylation</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>CD2-associated protein</td><td>Q9Y5K6</td><td>71,454</td><td>16</td><td>0</td><td>0</td><td>0</td><td>32%</td></tr><tr><td>Acyl-CoA-binding domain-containing protein 5</td><td>Q5T8D3</td><td>60,092</td><td>6</td><td>0</td><td>0</td><td>0</td><td>12%</td></tr><tr><td>Sorbin and SH3 domain-containing protein 1</td><td>Q9BX66</td><td>142,517</td><td>5</td><td>0</td><td>0</td><td>0</td><td>6%</td></tr><tr><td>Filamin-B</td><td>O75369</td><td>278,159</td><td>6</td><td>0</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Protein kinase C and casein kinase substrate in neurons protein 2</td><td>Q9UNF0</td><td>55,739</td><td>5</td><td>0</td><td>0</td><td>0</td><td>12%</td></tr><tr><td>Serine/threonine-protein kinase Nek9</td><td>Q8TD19</td><td>107,170</td><td>4</td><td>0</td><td>0</td><td>0</td><td>6%</td></tr><tr><td>Arsenite methyltransferase</td><td>Q9HBK9</td><td>41,748</td><td>3</td><td>0</td><td>0</td><td>0</td><td>10%</td></tr><tr><td>Hepatocyte growth factor-regulated tyrosine kinase substrate</td><td>O14964</td><td>86,192</td><td>2</td><td>0</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Protein transport protein Sec16A</td><td>O15027</td><td>233,515</td><td>3</td><td>0</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Vinexin</td><td>O60504</td><td>75,342</td><td>3</td><td>0</td><td>0</td><td>0</td><td>8%</td></tr><tr><td>Peroxisomal membrane protein PEX14</td><td>O75381</td><td>41,237</td><td>2</td><td>0</td><td>0</td><td>0</td><td>6%</td></tr><tr><td>Integrin beta-1</td><td>P05556</td><td>88,415</td><td>3</td><td>0</td><td>0</td><td>0</td><td>5%</td></tr><tr><td>Alpha actinin-1</td><td>P12814</td><td>103,061</td><td>2</td><td>0</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Nuclear factor NF-kappa-B p100 subunit</td><td>Q00653</td><td>96,751</td><td>3</td><td>0</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>Deubiquitinating protein VCIP135</td><td>Q96JH7</td><td>134,322</td><td>3</td><td>0</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>Tumor susceptibility gene 101 protein</td><td>Q99816</td><td>43,946</td><td>3</td><td>0</td><td>0</td><td>0</td><td>9%</td></tr><tr><td>NSFL1 cofactor p47</td><td>Q9UNZ2</td><td>40,573</td><td>3</td><td>0</td><td>0</td><td>0</td><td>14%</td></tr><tr><td>Protein phosphatase 1 regulatory subunit 12A</td><td>O14974</td><td>115,283</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>U2 snRNP-associated SURP motif-containing protein</td><td>O15042</td><td>118,296</td><td>2</td><td>0</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Heterogeneous nuclear ribonucleoprotein R</td><td>O43390</td><td>70,944</td><td>2</td><td>0</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>Transforming acidic coiled-coil-containing protein 1</td><td>O75410</td><td>87,796</td><td>2</td><td>0</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>NADH dehydrogenase [ubiquinone] 1 beta subcomplex subunit 4</td><td>O95168</td><td>15,210</td><td>2</td><td>0</td><td>0</td><td>0</td><td>22%</td></tr><tr><td>Integrin alpha-2</td><td>P17301</td><td>129,298</td><td>2</td><td>0</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Heterogeneous nuclear ribonucleoproteins A2</td><td>P22626</td><td>37,430</td><td>2</td><td>0</td><td>0</td><td>0</td><td>8%</td></tr><tr><td>Nuclear pore complex protein Nup214</td><td>P35658</td><td>213,619</td><td>2</td><td>0</td><td>0</td><td>0</td><td>1%</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 (Continued)</caption><thead><tr><th>Protein name</th><th>Accession number</th><th>Size (Da)</th><th>Unique peptides</th><th>coverage</th></tr></thead><tbody><tr><td>Peripherin</td><td>P41219</td><td>53,652</td><td>2</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>E3 ubiquitin-protein ligase NEDD4</td><td>P50402</td><td>149,118</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Emerin</td><td>P50402</td><td>28,995</td><td>2</td><td>0</td><td>0</td><td>9%</td></tr><tr><td>14-3-3 protein gamma</td><td>P61981</td><td>28,303</td><td>2</td><td>0</td><td>0</td><td>10%</td></tr><tr><td>Focal adhesion kinase 1</td><td>Q05397</td><td>119,236</td><td>2</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Golgin subfamily A member 2</td><td>Q08379</td><td>113,086</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Protein scribble homolog</td><td>Q14160</td><td>174,887</td><td>2</td><td>0</td><td>0</td><td>1%</td></tr><tr><td>Proteasome activator complex subunit 4</td><td>Q14997</td><td>211,341</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Leucine-rich repeat flightless-interacting protein 1</td><td>Q32MZ4</td><td>89,253</td><td>2</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>Striatin-interacting protein 1</td><td>Q5VSL9</td><td>95,579</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Ankyrin repeat and LEM domain-containing protein 2</td><td>Q86XL3</td><td>104,117</td><td>2</td><td>0</td><td>0</td><td>4%</td></tr><tr><td>Cell division cycle and apoptosis regulator protein 1</td><td>Q8IX12</td><td>132,824</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Proline-, glutamic acid- and leucine-rich protein 1</td><td>Q8IZL8</td><td>119,701</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>EH domain-binding protein 1</td><td>Q8NDI1</td><td>140,021</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Pre-B-cell leukemia transcription factor-interacting protein 1</td><td>Q96AQ6</td><td>80,643</td><td>2</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Dedicator of cytokinesis protein 7</td><td>Q96N67</td><td>242,566</td><td>2</td><td>0</td><td>0</td><td>1%</td></tr><tr><td>182 kDa tankyrase-1-binding protein</td><td>Q9C0C2</td><td>181,795</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Tyrosine-protein phophatase non-receptor type 23</td><td>Q9H3S7</td><td>178,975</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Exocyst complex component 1</td><td>Q9NV70</td><td>101,984</td><td>2</td><td>0</td><td>0</td><td>2%</td></tr><tr><td>Prefoldin subunit 2</td><td>Q9UHV9</td><td>16,648</td><td>2</td><td>0</td><td>0</td><td>17%</td></tr><tr><td>Protein argonaute-2</td><td>Q9UKV8</td><td>97,208</td><td>2</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Melanoma-associated antigen D2</td><td>Q9UNF1</td><td>64,956</td><td>2</td><td>0</td><td>0</td><td>3%</td></tr><tr><td>Nuclear migration protein nudC</td><td>Q9Y266</td><td>38,244</td><td>2</td><td>0</td><td>0</td><td>8%</td></tr></tbody></table>

replication complex markers. Fluorescence intensity plots also indicated that the two markers are close from one another but not colocalized (Figure 3E).

ACBD5 interacts with VAP-A and VAP-B via a FFAT motif and has been reported to be involved in the interaction between peroxisomes and the membrane of the endoplasmic reticulum (Hua et al., 2017; Costello et al., 2017a). Since VAP-A and VAP-B also

interact with NS5A in HCV-infected cells and are recruited to HCV-replication complexes (Tu et al., 1999; Hamamoto et al., 2005), we investigated whether ACBD5 could be involved in the association of peroxisomes and HCV replication complexes through its interaction with VAP-A/B. The gene encoding ACBD5 was inactivated by CRISPR-Cas9 technology (Supplementary Figure S4) and a HCV SGR was introduced in ACBD5-KO cells. No difference
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>FIGURE 2</figcaption></figure>

(A) APEX2-mediated biotinylation of ACBD5. Huh-7 cells containing APEX2- or GFP-tagged HCV replicons (SGR) were pre-incubated in the presence (+) or the absence (−) of BP, treated with H₂O₂ for 1 min, quenched and lysed in RIPA buffer. Biotinylated proteins were purified using streptavidin beads and analyzed by immunoblot, together with 8% of the inputs (lysates), with antibodies to ACBD5, HCV NS3 and tubulin. (B–E) Intracellular localization of GFP-ACBD5 transfected in naïve Huh-7 cells (B) or in HCV SGR-containing Huh-7 cells (C,D). As a control, GFP was expressed in SGR-containing cells (E). Transfected cells were fixed and processed for immunofluorescence staining of HCV NS5A (left panels, red in merge panels) and PEX14 (third panels from the left, blue in merge panels), and imaged with a confocal microscope. GFP-ACBD5 and GFP (second panels from the left, green in merge panels) were detected using GFP fluorescence and nuclei were stained with DAPI (grey in merge panels). Yellow arrowheads indicate examples of regions where markers are co-distributed. Bar, 10 µm.

of replication efficiency was observed between ACBD5-KO and control Huh-7 cells (data not shown), indicating that ACBD5 is not essential for HCV replication. No major differences in peroxisome phenotypes regarding both association to HCV-replication complexes and morphological alterations was observed in ACBD5-KO cells as

compared to control Huh-7 cells (Supplementary Figure S4). This result suggests that ACBD5 is not involved in peroxisome association to HCV replication complexes. Alternatively, ACBD5 could be part of a complex involved in this association, but its inactivation would not be sufficient to inhibit it.
<figure><img src="image_7.png" /><figcaption>FIGURE 3</figcaption></figure>

(A–C) Intracellular localization of endogenous PEX14 in naïve (A) or HCV SGR-containing Huh-7 cells (B,C). Cells were fixed and processed for immunofluorescence staining of HCV NS5A (left panels, green in merge panels) and PEX14 (middle panels, red in merge panels), and imaged with a high-resolution (Airyscan) confocal microscope. Nuclei were stained with DAPI (blue in merge panels). Bars, 10 µm. (D) Analysis of ACBD5 colocalization with HCV NS3 and NS5A. HCV SGR-containing Huh-7 cells were double labeled with antibodies to ACBD5 and NS3 or NS5A, and the Pearson correlation coefficient (PCC) was calculated on at least 30 individual cells imaged with a confocal microscope. (E) Fluorescence intensity plot along the double arrowed line indicated in B, merge panel.

## 3.4. Alteration of peroxisomes morphology in HCV-infected cells

Next, we investigated peroxisome morphology in HCV-infected cells. Huh-7 cells were infected with HCV strains JFH1 and DBN3a and fixed at different dpi to analyze peroxisome morphology by immunofluorescence microscopy. For this series of experiments, we used the culture adapted DBN3a strain and the original, non-adapted JFH1 strain, because the cell culture-adapted JFH1 strain did not allow

to keep infected cells for up to 16 dpi. Like in SGR-containing cells, some infected Huh-7 cells displayed peroxisome/replication complex association and altered peroxisome morphology for both strains (Figures 4A–C). The number of cells displaying altered peroxisomes increased over time and reached 20–30% at 16 dpi, a number similar to what was observed in cells permanently expressing an SGR (Figure 4D).

We characterized these morphological alterations using confocal microscopy and 3D-reconstitution in naïve and infected cells over time. The average volume of a peroxisome increased steadily from
0.15±0.01 µm³ in naïve Huh-7 cells to 0.35±0.02 µm³ at 16 dpi in DBN3a-infected cells and to 0.33±0.02 µm³ in JFH1-infected cells (Figure 4E). At the same time, the average number of peroxisomes per cell dropped from 547±94 in naïve Huh-7 cells to 242±59 in DBN3a-infected cells and 204±95 in JFH1-infected cells (Figure 4F). As a consequence of this increase in mean volume coupled to a decrease in peroxisome number per cell, the total volume occupied by peroxisomes in a cell displaying altered peroxisome morphology did not significantly change over time from their total volume in non-infected cells.

## 3.5. Analysis of peroxisomal proteins expression by mass spectrometry

To further assess if these morphological alterations are linked to biochemical alterations of peroxisomes, their protein content was analyzed by mass spectrometry. Total cellular proteins of control Huh-7 cells, of JFH1 replicon-containing cells and of DBN3a replicon-containing cells (50 µg per sample, 4 samples per cell line) were subjected to analysis. Overall, 4,072 proteins were identified, and the expression levels of proteins identified in at least 3 samples from at least one cell line (3,768 proteins), of which 63 proteins with a GO annotation: peroxisome, were used to compare protein expression between samples. Overall, 14 proteins were found to be differentially

expressed in JFH1 replicon-containing cells (Figure 5A) and 97 proteins in DBN3a replicon-containing cells (Figure 5B), of which the only protein with a “peroxisome” GO annotation was vimentin in the DBN3a replicon. None of the peroxisomal membrane and matrix proteins was dysregulated in a significant manner. These data indicate that the variation of expression of peroxisomal proteins in SGR-containing cells is too moderate to be appreciated by differential analysis with only four replicates per group. In particular, some proteins, such as PECR for example, show a fold change comparable to differentially expressed proteins (Figure 5B) but the variability between replicates is too high to observe an adjusted value of p below 0.05. Moreover, these differential analyses were performed protein by protein. We observed that multivariate analyses separated Huh-7 from DBN3a and JFH1 when looking at only peroxisome proteins (Figure 5C). This suggests that even if individual peroxisome proteins present a moderate change of expression levels, HCV replication has a significant effect on their expression when being taken together.

## 3.6. Role of peroxisomes during the HCV life cycle

To assess the importance of peroxisomes in the HCV life cycle, we generated a series of Huh-7-derived cells with non-functional peroxisomes or no peroxisomes at all. Cells with non-functional

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

FIGURE 4

Alteration of peroxisome morphology in HCV-infected and replicon-containing Huh-7 cells. (A–C) Immunofluorescence analysis of NS5A (green) and PEX14 (red) at 7 days after infection with JFH1 (A) or DBN3a (B) strains, or in naïve Huh-7 cells (C). (D) Number of cells with altered (larger and/or tubule-like) peroxisomes in HCV infected cells at different time points post infection and in Huh-7 cells permanently expressing a JFH1 or DBN3a SGR. (E) Mean volume of peroxisomes. (F) Mean number of peroxisomes per cell. Bars represent the mean +/- SD of 30 cells from 3 different infections.
*, p < 0.05; **, p < 0.01; ***, p < 0.001; ****, p < 0.0001 (Kruskall-Wallis test).
<figure><img src="image_9.png" /><figcaption>FIGURE 5</figcaption></figure>

Differential expression analysis of proteins in JFH1 (A) and DBN3a (B) SGR-containing cells. Proteins of total cell lysates were analyzed by LC-MS/MS (4 replicates per sample). Peroxisome proteins are displayed in black and other identified proteins either in pink (those with a significant difference of expression) or in grey (non differentially expressed). (C) Multivariate analysis of peroxisomal proteins expression using PLS-DA (Partial Least-Squares Discriminant Analysis). Groups to be discriminated were control Huh-7 cells (Hu) and DBN3a and JFH1 SGR-containing cells.

peroxisomes were generated by inactivating the PEX5 gene. The Pex5 protein is involved in the import of peroxisomal matrix proteins containing a type 1 peroxisomal targeting signal (PTS1), consisting of a carboxy-terminal -serine-lysine-leucine (-SKL) tripeptide (Ma et al., 2011) and the inactivation of the PEX5 gene leads to non-functional peroxisomes (Baes et al., 1997). Three Huh-7 derived PEX5-KO cell lines were generated with 3 different guide RNAs using the CRISPR-Cas9 technology. These cell lines had less than 5% Pex5 protein expression levels (Figure 6A), as estimated by immunoblotting. Because cells were not cloned after the transfection of the CRISPR-Cas9 plasmid, we can expect that a percentage of them express a functional Pex5 protein, due to the random mode of repair exerted by the cell's NHEJ pathway after the genomic DNA has been cut by the Cas9 nuclease. Catalase is imported in peroxisomes in a Pex5-dependent manner (Freitas et al., 2011). To evaluate the number of cells with non-functional peroxisomes, we labeled cells by immunofluorescence with an anti-catalase antibody to visualize its intracellular localization. In control cells, catalase pattern of labeling was made of dot-like structures, corresponding to peroxisomes, while this pattern was lost in PEX5-KO cells and catalase appeared merely cytosolic (Figure 6B). Using this assay, the number of cells with no functional peroxisomes was quantified as approximately 90% for the 3 PEX5-KO cell lines (Figure 6C).

Since peroxisomes are partially recruited near replication complexes, we first studied the impact of the absence of functional peroxisomes on HCV replication using luciferase-expressing SGRs of HCV strains JFH1 and DBN3a. We monitored their replication during 4 days and did not observe any significant difference in HCV replication (Figures 6D,E). We also titered infectious particles secreted by PEX5-KO cell lines. Cells were infected and infectious particles were measured in cell culture media at 3 dpi for JFH1. For DBN3a, which is less toxic than the cell culture-adapted JFH1 strain, infectious titers could be measured up to 16 dpi (Figure 4). Huh-7w7 cells, which do not express CD81 (Rocha-Perugini et al., 2009), were used as controls. Again, the absence of functional peroxisomes had no impact on HCV infectious titers (Figures 6F,G). These results indicate that functional peroxisomes are not required for HCV to establish an infection.

To further assess the role of peroxisomes during the HCV life cycle, we generated PEX3-KO cell lines. The Pex3 protein is a component of the peroxisome membrane protein insertion machinery (Ma et al., 2011), and its inactivation results in cells completely devoid of peroxisomes (Hettema et al., 2000). Huh-7-derived PEX3-KO cell lines were generated with 3 different guide RNAs. However, only about 50% of cells displayed a catalase pattern consistent with a complete absence of peroxisome. Therefore, we isolated by limited dilution cellular clones completely devoid of peroxisomes using the
<figure><img src="image_10.png" /></figure>

catalase pattern as a screening method. The 4 selected clones contained 100% PEX3 KO cells as estimated by catalase staining. No Pex3 protein was detected by immunoblotting in these clones, and the peroxisome membrane protein PMP70 was also absent, as expected for a protein inserted in peroxisome membrane by a PEX3-dependent mechanism (Hettema et al., 2000). In contrast, the catalase expression levels were normal, in keeping with its cytosolic localization in the absence of peroxisomes (Figure 7A). As clones 1.7 and 1.6C originated from the same PEX3-KO cell line, we only used one of these clones to analyze HCV replication. In the luciferase-based replication assay, no differences in replication were observed between PEX3 KO clones and control Huh-7 cells both for JFH1 (Figure 7B) and DBN3a (Figure 7C). Surprisingly, when infectious titers were measured in the same PEX3-KO clones infected with JFH1, two clones yielded higher titers and the third one a lower titer than control Huh-7 cells at 3 dpi (Figure 7D). These differences are not correlated to Pex3 expression levels and most likely result from clonal differences unrelated to PEX3 gene inactivation. When the same clones were infected with DBN3a, no significant differences from control Huh-7 cells were observed up to 16 dpi, except for the clone 3.2E, which displayed a lower titer than the other two clones and Huh-7 cells at 16 dpi only (Figure 7E). It is worth noting that clone 3.2E is one of the clones that yielded higher titers with JFH1 at 3 dpi (Figure 7D). Again, this difference was attributed to a clonal effect rather than to PEX3 gene inactivation. Altogether, these results confirmed that HCV is able to replicate to near normal level in the absence of peroxisomes.

## 3.7. Impact of HCV infection on peroxisome functions

As peroxisomes do not apparently play any important role in the HCV life cycle, we sought to assess whether HCV infection would have an impact on peroxisome functions. Peroxisomes play a series of unique functions in the cell metabolism, of which we chose to investigate their role in lipid metabolism and in ROS production and degradation, which are two functions potentially related to HCV replication and pathogenesis.

To probe potential peroxisome dysfunction, we assessed whether the lipid metabolism of peroxisomes is altered in cells replicating HCV. The *β*-oxidation of very long chain fatty acids (VLCFA) occurs exclusively in peroxisomes, whereas mitochondria use shorter fatty acids. Therefore, an excess of very long chain fatty acids is indicative of peroxisome dysfunction (Moser et al., 1999). Moreover, the *β*-oxidation of VLCFA is considered the main source of H₂O₂ in peroxisomes. To assess peroxisome *β*-oxidation, we measured the total cellular content of C16 and C26 fatty acids in control and in SGR-containing Huh-7 cells, and in PEX3-KO cells as a control of cells with no active peroxisomes. As shown in Supplementary Figure S5, the ratio of C26 fatty acid over C16 fatty acid did not significantly change in cells containing sub-genomic replicons of JFH1 or DBN3a strains, as compared with control Huh-7 cells. In contrast, the C26:C16 ratio was significantly higher in PEX3-KO cells, as expected. These results suggest that the
A

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

B

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

C

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

D

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

E

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

FIGURE 7

HCV infection in PEX3-KO Huh-7 cells. (A) Immunoblot analysis of Pex3 protein, Pmp-70, catalase and tubulin expression in control and selected clones of PEX3-KO cells. (B,C) Huh-7 cells were electroporated with JFH1 SGR (B) or a DBN3a SGR (C) expressing Renilla luciferase or with a non-replicative SGR (GND). Samples were harvested for luciferase assay at 4, 24, 48, 72, and 96 h post electroporation. The luciferase activity at 4 h post electroporation is expressed as 1. Error bars indicate standard errors of the means for 3 independent experiments performed in quadruplicate samples. (D) Huh-7, CD81-deficient Huh-7w7 and PEX3-KO cells were infected with JFH1 and cultured for 3 days. Supernatants were collected, filtered and infectious titers were determined. (E) Huh-7, CD81-deficient Huh-7w7 and PEX3-KO cells were infected with DBN3a and passed at 3, 6 and 15 dpi. At each passage, 1.5x10⁵ cells were seeded in a well of a P-24 plate. Twenty-four hours later, supernatants were collected and infectious titers were determined. *, *p* < 0.05; **, *p* < 0.01 (Mann–Whitney test).

peroxisomal β-oxidation is not significantly affected by HCV infection.

To specifically survey peroxisomal ROS content, we used a peroxisome-localized roGFP2-based probe (Ivashchenko et al., 2011). For comparison, we used cytosolic and mitochondrial roGFP2 constructs. The 3 probes were expressed in Huh-7 cells and their correct intracellular localizations were verified by immunofluorescence (Figures 8A–C). The roGFP-2 sensor is based on the presence in GFP of 2 surface-exposed cysteine residues at positions 147 and 204. These 2 Cys residues can form a disulfide bond in an oxidative environment, and the presence of this disulfide bond alters the excitation spectrum of the probe (Dooley et al., 2004). Relative to reduced roGFP-2, oxidized roGFP-2 fluorescence excitation is increased at 405 nm and reduced at 488 nm, with no change of emission at 514 nm. Therefore, the ratio of roGFP2 fluorescence signals measured at 405 nm versus 488 nm (ratio 405/488) increases in relation to the oxidative content in the sub-cellular environment containing the probe. When this ratiometric redox sensor was expressed in control and SGR-containing Huh-7 cells, an increased fluorescence ratio was detected with the peroxisomal probe (Figure 8D). A small difference was also observed with the mitochondrial probe for the JFH1-SGR, but not for the DBN3a-SGR (Figure 8E). No difference of fluorescence ratio was

observed for the cytosolic probe both for DBN3a-based replicon (Figure 8F). These results indicate that the redox potential of peroxisomes is increased in cells replicating HCV.

To assess whether this increase of ROS in peroxisomes could result from a defect in H₂O₂ degradation, we measured the cellular activity of catalase in cell lysates of Huh-7 and SGR-containing cells. As shown in Figure 8G, catalase activity was reduced in cells containing the DBN3a-based replicon as compared with Huh-7 cells, and also appear to be reduced in cells containing the JFH1-based replicon, although the difference was not significant for this strain. This lower activity was not correlated with lower levels of catalase expression in SGR-containing cells, as determined by immunoblotting (Figure 8H). These results suggest that the higher peroxisomal ROS levels in HCV-infected cells could at least in part result from a defect in catalase activity.

## 4. Discussion

Viral infections often result in alterations of the cellular metabolism and of the morphology of intracellular compartments. For positive-strand RNA viruses, which replicate in association with
A

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

B

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

C

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

D

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

E

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

F

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

G

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

H

<figure><img src="image_23.png" /><figcaption>FIGURE 8</figcaption></figure>

Alteration of ROS metabolism in HCV SGR-containing cells. (A) Immunofluorescence of peroxisomal roGFP2 sensor (green) transfected in Huh-7 cells co-stained with PEX14 (red). (B) Immunofluorescence of mitochondrial roGFP2 sensor (green) transfected in Huh-7 cells co-stained with TOM20 (red). (C) Immunofluorescence of cytosolic roGFP2 sensor (green) transfected in Huh-7 cells. (D–F) *In situ* measurements of the oxidation state of roGFP2 sensors in control and HCV SGR-containing Huh-7 cells transiently transfected with a peroxisomal (D), a mitochondrial (E) or a cytosolic (F) roGFP2-based sensor. Fluorescence images of live cells were acquired and the ratio of fluorescence emitted from 405- and 488-nm excitation were used to calculate the 405/488-nm excitation ratio. At least 90 cells from 3 independent transfections were used for each sample. Ratio values are presented as violin plots with the median and quartiles indicated as solid and dotted lines, respectively. (G) Catalase activity in cell lysates of control and SGR-containing Huh-7 cells. (H) Immunoblot analysis of catalase and tubulin expression in control and HCV SGR-containing Huh-7 cells. \*, *p* < 0.05; \*,\*, *p* < 0.01; \*,\*,\*, *p* < 0.001; \*,\*,\*,\*, *p* < 0.0001 (Kruskall-Wallis test).

rearranged cytoplasmic membranes, the stress imposed by the
establishment of replication complexes can result in an intense
reorganization of various intracellular compartments. ER stress and
alterations of mitochondria and lipid droplets have previously been
documented in HCV-infected cells (Tardif et al., 2004; McLauchlan,
2009; Arciello et al., 2013). In this study, we found that peroxisome
morphology and functions are also altered in HCV-infected cells. In
line with our results, two recently published multi-omics studies

reported peroxisome dysfunctions in cells infected with HCV *in vitro*
and in patients samples (Lupberger et al., 2019; Ali et al., 2023).

Using a proximity-based biotinylation approach, we unexpectedly
uncovered a close association of peroxisomes and HCV replication
complexes coupled to an alteration of peroxisome morphology. Two
peroxisome membrane proteins, ACBD5 and PEX14, were detected
in the set of proteins biotinylated by HCV NS5A-APEX2 fusion in the
context of cells containing a sub-genomic HCV replicon. The
biotinylation of these proteins did not represent a relocalization of individual peroxisomal proteins from peroxisomes to replication complexes in infected cells, but rather the recruitment of whole peroxisomes in the vicinity of HCV replication complexes. ACBD5, a protein involved in ER/peroxisome membrane contact sites and in the transfer of lipids between peroxisomes and the ER (Hua et al., 2017; Costello et al., 2017a), ranked second in the list of proteins biotinylated by NS5A-APEX2. This is in keeping with its association to ER proteins VAPA and VAPB, which interact with NS5A in HCV-infected cells. We hypothesized that, for this reason, ACBD5 could be in close proximity to NS5A, through its VAPA/B interaction and that it could be instrumental in the association of peroxisomes and replication complexes. However, no difference of association between these two compartments could be observed in cells lacking ACBD5. This indicates that ACBD5 is either not involved in peroxisome-replication complex association, or that it is part of a larger set of proteins collectively responsible for this association. This second possibility could explain the lack of impact of ACBD5 gene inactivation on peroxisome-HCV replication complex association. Other membrane contact sites, such as ER-plasma membrane or ER-mitochondria membrane contact sites are known to involve numerous proteins (Wu et al., 2018). Not much is known about the protein composition of ER-peroxisome membrane contact sites in human cells, except that ACBD5 (Hua et al., 2017; Costello et al., 2017a) and a peroxisomal ACBD4 isoform (Costello et al., 2017b) are involved.

In addition to peroxisomes recruitment to the vicinity of HCV replication complexes, an alteration of peroxisome morphology was also observed both in HCV-infected and SGR-containing cells. These morphological alterations included enlarged peroxisomes and also extended tubules. Tubules were mostly visible in SGR-containing cells, whereas enlarged peroxisomes were present in both conditions. The appearance of extended tubules in SGR-containing cells but not very frequently in infected cells could result from a longer time of HCV replication in SGR-containing cells, as compared to infected cells, and could therefore represent a more advanced phenotype than the one visible in infected cells, which were kept in culture for a maximum of 16 days in our experiments. Alternatively, the difference of phenotype could also result from the difference of intracellular localization of HCV replication complexes in infected cells, in which they are recruited to lipid droplet-associated virion assembly sites via NS5A interaction with core (Miyanari et al., 2007). In contrast, replication complexes are not associated to lipid droplets in SGR-containing cells. The presence of peroxisomes with an altered morphology did not completely correlate with their association with replication complexes, suggesting that their recruitment is not likely to be the only trigger of their change of morphology. The metabolic state of cells can influence peroxisome size and number (Farré et al., 2019). HCV infection increases glucose metabolism and the STAT3 signaling pathway and it has been suggested that these metabolic changes reduce peroxisome function (Lupberger et al., 2019). We could hypothesize that these metabolic changes could also impact peroxisome morphology.

The turnover of peroxisomes is quite short (Islinger et al., 2018; Farré et al., 2019). Peroxisomes are selectively degraded by a special mode of autophagy named pexophagy, and new peroxisomes are mostly generated by growth and fission. An alternative mode of *de novo* production of peroxisomes also exists by fusion of two pre-peroxisomal vesicles originating from ER and mitochondria

(Sugiura et al., 2017). Our results indicate a reduction in the number of peroxisomes in infected cells, as previously reported (Lupberger et al., 2019). However, we also observed that the mean volume of peroxisomes is increased in HCV-infected cells and the total volume occupied by peroxisomes per cell does not appear to significantly change over time during infection. This suggests that the change of peroxisome morphology is not likely to result from increased or decreased peroxisome turnover in HCV-infected cells. A defect in peroxisome fission could result in a phenotype similar to that observed in infected cells. The presence of elongated tubule-like peroxisomal structures in SGR-containing cells also points to a possible defect in peroxisome fission. In yeast, the deletion of genes involved in peroxisome fission results in an increased size and a reduced number of peroxisomes (Yofe et al., 2017). Alternatively, larger structures could result from peroxisome clustering and aggregation. It would be interesting to image peroxisome alterations by electron microscopy. However, peroxisomes are small single membrane-bound organelles ranging in size from 0.1 to 1 µm, that are not easy to identify in electron microscopy in human cells, in contrast to rodent peroxisomes, easily identifiable by the presence of a prominent crystalloid core, which is absent in human cells.

Peroxisomes are not essential for cell viability, which made it possible to evaluate their implication in the establishment of HCV infection. Surprisingly, no impact on the HCV life cycle was observed in cells devoid of peroxisomes for the two strains we used up to 16 dpi, despite the alterations of their morphology and of peroxisomal ROS metabolism during the same period in normal peroxisome-containing cells. Other viruses have been reported to interact with peroxisomes. Rotavirus VP4 protein contains a PST1 motif and is localized in peroxisomes (Mohan et al., 2002), but the importance of peroxisomes in rotavirus infection remains undetermined. Viruses of the *Herpesviridae* family have been shown to increase peroxisome numbers in infected cells (Sychev et al., 2017; Jean Beltran et al., 2018). Interestingly, for human cytomegalovirus, it has been shown that the biogenesis of peroxisomes by fission is enhanced and coupled to alterations in peroxisome morphology, protein composition and lipid metabolism, which result in an increased production of plasmalogens that are required for virion assembly (Jean Beltran et al., 2018). This increase in plasmalogens production is linked to a reorganization of ER/peroxisome contact sites and an increased association of peroxisomes to the ER membrane (Cook et al., 2022). Other viruses, such as Dengue virus and West Nile virus, two viruses of the *Flaviviridae* family, cause a decrease in peroxisome number by targeting their capsid protein to the peroxisomal membrane (You et al., 2015), a mechanism different from that of HCV, for which the core protein, which is absent from SGRs, does not appear to be involved in peroxisome alteration. This decrease in peroxisomes in flavivirus-infected cells results in an inhibition of type III interferon expression. Indeed, it has been shown that the innate immunity adaptor MAVS is located both on mitochondrial and peroxisomal membranes and promotes different antiviral responses at each location (Dixit et al., 2010). Moreover, mitochondrial MAVS preferentially induces a type I interferon response, whereas peroxisomal MAVS exclusively induces a type III interferon response (Odendall et al., 2014). HCV uses another mechanism to counteract MAVS signaling. In HCV-infected cells, MAVS is cleaved by NS3 both on mitochondrial and peroxisomal membranes (Bender et al., 2015), and this cleavage inhibits all interferon responses (Foy et al., 2005; Meylan et al., 2005).
In HCV-infected cells, we did not detect any alteration of the lipid
metabolism part, which deals with the β-oxidation of VLCFA, a
metabolic pathway exclusively performed by peroxisomes. This finding
is at odds with lipidomics studies of patient samples, in which an excess
of VLCFAs was reported (Lupberger et al., 2019; Ali et al., 2023). The
reason for this discrepancy is still unclear, but may have to do with the
different system used (infectious virus and SGR). It would be interesting
to investigate if any of the HCV proteins absent in SGR could impact
peroxisomal β-oxidation, either by itself or in combination with
proteins of the HCV replication complex. In contrast, we found an
alteration of peroxisomal redox potential during HCV infection. Our
results indicate that ROS levels are specifically increased in
peroxisomes. This result suggests that peroxisomes are the major site
of ROS increased production in HCV SGR-containing cells. HCV core
has been reported to induce the production of ROS (Okuda et al., 2002)
through mitochondrial electron transport inhibition (Korenaga et al.,
2005; T. Wang et al., 2010). This has led to the idea that mitochondria
are the major site of ROS production in HCV-infected cells. The use of
SGRs, which do not express core, allowed us to uncover the additional
role of non-structural proteins in ROS production by peroxisomes. The
β-oxidation of VLCFA by peroxisomes is a source of H₂O₂ production
in peroxisomes. However, our results do not indicate that VLCFA
accumulate or are more actively degraded in HCV SGR-containing
cells. Rather, a lower rate of ROS elimination by catalase (and possibly
other detoxifying enzymes) appears to be a possible cause of ROS
accumulation in peroxisomes, and this loss of catalase activity is not
related to a decreased catalase expression. It would be interesting to
decipher the mechanism of catalase inhibition and its impact on ROS
accumulation and pathogenesis induced. Interestingly, peroxisomes are
altered and catalase activity is reduced in hepatocellular tumors (Litwin
et al., 1999). Given the role of oxidative stress in cancer development
and progression (Cheung and Vousden, 2022), peroxisomes may have
an as yet underestimated contribution to these processes.

In conclusion, our study indicates that peroxisome function and
morphology are altered in HCV-infected cells. These perturbations are
likely to contribute to liver dysfunction and pathogenesis in
chronically infected patients.

Data availability statement

The original contributions presented in the study are included in
the article/Supplementary material, further inquiries can be directed
to the corresponding author.

Author contributions

EM: Conceptualization, Investigation, Methodology, Writing – review & editing. NC: Investigation, Methodology, Writing – review & editing. J-MS: Investigation, Writing – review & editing. MF: Writing – review & editing, Formal analysis. OD: Writing – review & editing, Investigation. NB: Investigation, Writing – review & editing. QT: Investigation, Writing – review & editing. SR: Writing – review & editing, Resources. JB: Resources, Writing – review & editing. LC: Writing

- review & editing. JB-M: Writing – review & editing, Investigation. GM: Writing – review & editing, Data curation. YS: Conceptualization, Funding acquisition, Investigation, Methodology, Writing – review & editing. JD: Conceptualization, Writing – review & editing. YR: Conceptualization, Funding acquisition, Investigation, Methodology, Writing – original draft.

Funding

This work was supported in part by grants from the CPER-CTRL (Centre Transdisciplinaire de Recherche sur la Longévité) program (to YR and YS) and from the French Agence Nationale de Recherche sur le Sida et les Hépatites Virales (ANRS) (to YR). EM was supported by a pre-doctoral fellowship from ANRS.

Acknowledgments

We thank M. Schrader, M. Cormont, M. Fransen, and F. Zhang
for providing plasmids, Takaji Wakita for JFH1 virus and replicon,
Charles M. Rice and Jean-François Delagneau for antibodies, Hervé
Drobecq for help in BP synthesis, Sophie Salomé-Desnoulez and
Elizabeth Werkmeister for help with image acquisition and analysis,
and Sophana Ung for help in assembling figures. Immunofluorescence
and electron microscopy analyses were performed with the help of
the imaging core facility of the BioImaging Center Lille (BiCeL).
Lipidomic analyses were performed in MetaToul facility (Toulouse
metabolomics and fluxomics facilities, www.metatoul.fr), which is
part of the French National Infrastructure for Metabolomics and
Fluxomics MetaboHUB-ANR-11-INBS-0010. We are grateful to the
Institut Français de Bioinformatique (IFB) (ANR-11-INBS-0013) and
the Bilille platform (UAR 2014 - US 41 - PLBS) for providing help in
statistical analysis of proteomics data.

Conflict of interest

The authors declare that the research was conducted in the
absence of any commercial or financial relationships that could
be construed as a potential conflict of interest.

Publisher's note

All claims expressed in this article are solely those of the authors
and do not necessarily represent those of their affiliated organizations,
or those of the publisher, the editors and the reviewers. Any product
that may be evaluated in this article, or claim that may be made by its
manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online
at: https://www.frontiersin.org/articles/10.3389/fmicb.2023.1254728/
full#supplementary-material
## References

Ali, R. O., Quinn, G. M., Umarova, R., Haddad, J. A., Zhang, G. Y., Townsend, E. C., et al. (2023). Longitudinal multi-omics analyses of the gut-liver axis reveals metabolic dysregulation in hepatitis C infection and cirrhosis. *Nat. Microbiol.* 8, 12–27. doi: 10.1038/s41564-022-01273-y

Arciello, M., Gori, M., and Balsano, C. (2013). Mitochondrial dysfunctions and altered metals homeostasis: new weapons to counteract HCV-related oxidative stress. *Oxidative Med. Cell Longev.* 2013:971024. doi: 10.1155/2013/971024

Baes, M., Gressens, P., Baumgart, E., Carmeliet, P., Casteels, M., Fransen, M., et al. (1997). A mouse model for Zellweger syndrome. *Nat. Genet.* 17, 49–57. doi: 10.1038/ng0997-49

Bender, S., Reuter, A., Eberle, F., Einhorn, E., Binder, M., and Bartenschlager, R. (2015). Activation of type I and III interferon response by mitochondrial and Peroxisomal MAVS and inhibition by hepatitis C virus. *PLoS Pathog.* 11:e1005264. doi: 10.1371/journal.ppat.1005264

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. *J. R. Stat. Soc. Ser. B* 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x

Bligh, E. G., and Dyer, W. J. (1959). A rapid method of total lipid extraction and purification. *Can. J. Biochem. Physiol.* 37, 911–917. doi: 10.1139/o59-099

Bouyssié, D., Hesse, A.-M., Mouton-Barbosa, E., Rompais, M., Macron, C., Carapito, C., et al. (2020). Proline: an efficient and user-friendly software suite for large-scale proteomics. *Bioinformatics* 36, 3148–3155. doi: 10.1093/bioinformatics/btaa118

Cheung, E. C., and Vousden, K. H. (2022). The role of ROS in tumour development and progression. *Nat. Rev. Cancer* 22, 280–297. doi: 10.1038/s41568-021-00435-0

Cong, L., Ran, F. A., Cox, D., Lin, S., Barretto, R., Habib, N., et al. (2013). Multiplex genome engineering using CRISPR/Cas systems. *Science* 339, 819–823. doi: 10.1126/science.1231143

Cook, K. C., Tsopurashvili, E., Needham, J. M., Thompson, S. R., and Cristea, I. M. (2022). Restructured membrane contacts rewrite organelles for human cytomegalovirus infection. *Nat. Commun.* 13:4720. doi: 10.1038/s41467-022-32488-6

Costello, J. L., Castro, I. G., Hacker, C., Schrader, T. A., Metz, J., Zeuschner, D., et al. (2017a). ACBD5 and VAPB mediate membrane associations between peroxisomes and the ER. *J. Cell Biol.* 216, 331–342. doi: 10.1083/jcb.201607055

Costello, J. L., Castro, I. G., Schrader, T. A., Islinger, M., and Schrader, M. (2017b). Peroxisomal ACBD4 interacts with VAPB and promotes ER-peroxisome associations. *Cell Cycle* 16, 1039–1045. doi: 10.1080/15384101.2017.1314422

de Chassey, B., Navratil, V., Tafforeau, L., Hiet, M. S., Aublin-Gex, A., Agaugué, S., et al. (2008). Hepatitis C virus infection protein network. *Mol. Syst. Biol.* 4:230. doi: 10.1038/msb.2008.66

Delgrange, D., Pillez, A., Castelain, S., Cocquerel, L., Rouillé, Y., Dubuisson, J., et al. (2007). Robust production of infectious viral particles in Huh-7 cells by introducing mutations in hepatitis C virus structural proteins. *J. Gen. Virol.* 88, 2495–2503. doi: 10.1099/vir.0.82872-0

Dixit, E., Boulant, S., Zhang, Y., Lee, A. S. Y., Odendall, C., Shum, B., et al. (2010). Peroxisomes are signaling platforms for antiviral innate immunity. *Cells* 141, 668–681. doi: 10.1016/j.cell.2010.04.018

Dolan, P. T., Zhang, C., Khadka, S., Arumugaswami, V., Vangeloff, A. D., Heaton, N. S., et al. (2013). Identification and comparative analysis of hepatitis C virus-host cell protein interactions. *Mol. BioSyst.* 9, 3199–3209. doi: 10.1039/c3mb70343f

Dooley, C. T., Dore, T. M., Hanson, G. T., Jackson, W. C., Remington, S. J., and Tsien, R. Y. (2004). Imaging dynamic redox changes in mammalian cells with green fluorescent protein indicators. *J. Biol. Chem.* 279, 22284–22293. doi: 10.1074/jbc.M312847200

Farré, J. C., Mahalingam, S. S., Proietto, M., and Subramani, S. (2019). Peroxisome biogenesis, membrane contact sites, and quality control. *EMBO Rep.* 20:e46864. doi: 10.15252/embr.201846864

Ferlin, J., Farhat, R., Belouzard, S., Cocquerel, L., Bertin, A., Hober, D., et al. (2018). Investigation of the role of GBF1 in the replication of positive-sense single-stranded RNA viruses. *J. Gen. Virol.* 99, 1086–1096. doi: 10.1099/jgv.0.001099

Ferraris, P., Blanchard, E., and Roingead, P. (2010). Ultrastructural and biochemical analyses of hepatitis C virus-associated host cell membranes. *J. Gen. Virol.* 91, 2230–2237. doi: 10.1099/vir.0.022186-0

Filipe, A., and McLauchlan, J. (2015). Hepatitis C virus and lipid droplets: finding a niche. *Trends Mol. Med.* 21, 34–42. doi: 10.1016/j.molmed.2014.11.003

Foy, E., Li, K., Sumpter, R., Loo, Y.-M., Johnson, C. L., Wang, C., et al. (2005). Control of antiviral defenses through hepatitis C virus disruption of retinoic acid-inducible gene-I signaling. *Proc. Natl. Acad. Sci. U. S. A.* 102, 2986–2991. doi: 10.1073/pnas.0408707102

Freitas, M. O., Francisco, T., Rodrigues, T. A., Alencastre, I. S., Pinto, M. P., Grou, C. P., et al. (2011). PEX5 protein binds monomeric catalase blocking its tetramerization and releases it upon binding the N-terminal domain of PEX14. *J. Biol. Chem.* 286, 40509–40519. doi: 10.1074/jbc.M111.287201

Galli, A., Ramirez, S., and Bukh, J. (2021). Lipid droplets accumulation during hepatitis C virus infection in cell-culture varies among genotype 1-3 strains and does not correlate with virus replication. *Viruses* 13:389. doi: 10.3390/v13030389

Germain, M.-A., Chatel-Chaix, L., Gagné, B., Bonneil, É., Thibault, P., Pradezynski, F., et al. (2014). Elucidating novel hepatitis C virus-host interactions using combined mass spectrometry and functional genomics approaches. *Mol. Cell. Proteomics* 13, 184–203. doi: 10.1074/mcp.M113.030155

Goonawardane, N., Gebhardt, A., Bartlett, C., Pichlmair, A., and Harris, M. (2017). Phosphorylation of serine 225 in hepatitis C virus NS5A regulates protein-protein interactions. *J. Virol.* 91, e00805–17. doi: 10.1128/JVI.00805-17

Goueslain, L., Alsaleh, K., Horellou, P., Roingead, P., Descamps, V., Duverlie, G., et al. (2010). Identification of GBF1 as a cellular factor required for hepatitis C virus RNA replication. *J. Virol.* 84, 773–787. doi: 10.1128/JVI.01190-09

Hamamoto, I., Nishimura, Y., Okamoto, T., Aizaki, H., Liu, M., Mori, Y., et al. (2005). Human VAP-B is involved in hepatitis C virus replication through interaction with NS5A and NS5B. *J. Virol.* 79, 13473–13482. doi: 10.1128/JVI.79.21.13473-13482.2005

Hettema, E. H., Giralzalsky, W., van Den Berg, M., Erdmann, R., and Distel, B. (2000). *Saccharomyces cerevisiae* pex3p and pex19p are required for proper localization and stability of peroxisomal membrane proteins. *EMBO J.* 19, 223–233. doi: 10.1093/emboj/19.2.223

Hu, B., Li, S., Zhang, Z., Xie, S., Hu, Y., Huang, X., et al. (2016). HCV NS4B targets scribble for proteasome-mediated degradation to facilitate cell transformation. *Tumour Biol.* 37, 12387–12396. doi: 10.1007/s13277-016-5100-4

Hua, R., Cheng, D., Coyaud, É., Freeman, S., Di Pietro, E., Wang, Y., et al. (2017). VAPs and ACBD5 tether peroxisomes to the ER for peroxisome maintenance and lipid homeostasis. *J. Cell Biol.* 216, 367–377. doi: 10.1083/jcb.201608128

Hui, J. M., Sud, A., Farrell, G. C., Bandara, P., Byth, K., Kench, J. G., et al. (2003). Insulin resistance is associated with chronic hepatitis C virus infection and fibrosis progression [corrected]. *Gastroenterology* 125, 1695–1704. doi: 10.1053/j.gastro.2003.08.032

Hung, V., Udeshi, N. D., Lam, S. S., Loh, K. H., Cox, K. J., Pedram, K., et al. (2016). Spatially resolved proteomic mapping in living cells with the engineered peroxidase APEX2. *Nat. Protoc.* 11, 456–475. doi: 10.1038/nprot.2016.018

Igloi, Z., Kazlauskas, A., Saksela, K., Macdonald, A., Mankouri, J., and Harris, M. (2015). Hepatitis C virus NS5A protein blocks epidermal growth factor receptor degradation via a proline motif-dependent interaction. *J. Gen. Virol.* 96, 2133–2144. doi: 10.1099/vir.0.000145

Islinger, M., Voelkl, A., Fahimi, H. D., and Schrader, M. (2018). The peroxisome: an update on mysteries 2.0. *Histochem. Cell Biol.* 150, 443–471. doi: 10.1007/s00418-018-1722-5

Ivashchenko, O., Van Veldhoven, P. P., Brees, C., Ho, Y.-S., Terlecky, S. R., and Fransen, M. (2011). Intraperoxisomal redox balance in mammalian cells: oxidative stress and interorganellar cross-talk. *Mol. Biol. Cell* 22, 1440–1451. doi: 10.1091/mbc.E10-11-0919

Jean Beltran, P. M., Cook, K. C., Hashimoto, Y., Galitzine, C., Murray, L. A., Vitek, O., et al. (2018). Infection-induced peroxisome biogenesis is a metabolic strategy for herpesvirus replication. *Cell Host Microbe* 24, 526–541.e7. doi: 10.1016/j.chom.2018.09.002

Kato, T., Date, T., Miyamoto, M., Furusaka, A., Tokushige, K., Mizokami, M., et al. (2003). Efficient replication of the genotype 2a hepatitis C virus subgenomic replicon. *Gastroenterology* 125, 1808–1817. doi: 10.1053/j.gastro.2003.09.023

Kato, T., Date, T., Murayama, A., Morikawa, K., Akazawa, D., and Wakita, T. (2006). Cell culture and infection system for hepatitis C virus. *Nat. Protoc.* 1, 2334–2339. doi: 10.1038/nprot.2006.395

Khan, I., Katikaneni, D. S., Han, Q., Sanchez-Felipe, L., Hanada, K., Ambrose, R. L., et al. (2014). Modulation of hepatitis C virus genome replication by glycosphingolipids and four-phosphate adaptor protein 2. *J. Virol.* 88, 12276–12295. doi: 10.1128/JVI.00970-14

Korenaga, M., Wang, T., Li, Y., Showalter, L. A., Chan, T., Sun, J., et al. (2005). Hepatitis C virus core protein inhibits mitochondrial electron transport and increases reactive oxygen species (ROS) production. *J. Biol. Chem.* 280, 37481–37488. doi: 10.1074/jbc.M506412200

Lindenbach, B. D., Evans, M. J., Syder, A. J., Wölk, B., Tellinghuisen, T. L., Liu, C. C., et al. (2005). Complete replication of hepatitis C virus in cell culture. *Science* 309, 623–626. doi: 10.1126/science.1114016

Lindenbach, B. D., Thiel, H.-J., and Rice, C. M. (2007). Flaviviridae: the viruses and their replication. *Fields Virol.* 1, 1101–1152.

Litwin, J. A., Beier, K., Vökl, A., Hofmann, W. J., and Fahimi, H. D. (1999). Immunocytochemical investigation of catalase and peroxisomal lipid beta-oxidation enzymes in human hepatocellular tumors and liver cirrhosis. *Virchows Arch.* 435, 486–495. doi: 10.1007/s004280050432

Lupberger, J., Croonenborghs, T., Roca Suarez, A. A., Van Renne, N., Jühling, F., Oudot, M. A., et al. (2019). Combined analysis of metabolomes, proteomes, and transcriptomes of hepatitis C virus-infected cells and liver to identify pathways associated with disease development. *Gastroenterology* 157, 537–551.e9. doi: 10.1053/j.gastro.2019.04.003
Ma, C., Agrawal, G., and Subramani, S. (2011). Peroxisome assembly: matrix and membrane protein biogenesis. *J. Cell Biol.* 193, 7–16. doi: 10.1083/jcb.201010022

McLauchlan, J. (2009). Lipid droplets and hepatitis C virus infection. *Biochim. Biophys. Acta* 1791, 552–559. doi: 10.1016/j.bbalip.2008.12.012

Meylan, E., Curran, J., Hofmann, K., Moradpour, D., Binder, M., Bartenschlager, R., et al. (2005). Cardif is an adaptor protein in the RIG-I antiviral pathway and is targeted by hepatitis C virus. *Nature* 437, 1167–1172. doi: 10.1038/nature04193

Miyanari, Y., Atsuzawa, K., Usuda, N., Watashi, K., Hishiki, T., Zayas, M., et al. (2007). The lipid droplet is an important organelle for hepatitis C virus production. *Nat. Cell Biol.* 9, 1089–1097. doi: 10.1038/ncb1631

Mohan, K. V. K., Som, I., and Atreya, C. D. (2002). Identification of a type 1 Peroxisomal targeting signal in a viral protein and demonstration of its targeting to the organelle. *J. Virol.* 76, 2543–2547. doi: 10.1128/jvi.76.5.2543-2547.2002

Monzo, P., Gauthier, N. C., Keslair, F., Loubat, A., Field, C. M., Le Marchand-Brustel, Y., et al. (2005). Clues to CD2-associated protein involvement in cytokinesis. *Mol. Biol. Cell* 16, 2891–2902. doi: 10.1091/mbc.e04-09-0773

Moser, A. B., Kreiter, N., Bezman, L., Lu, S., Raymond, G. V., Naidu, S., et al. (1999). Plasma very long chain fatty acids in 3,000 peroxisome disease patients and 29,000 controls. *Ann. Neurol.* 45, 100–110. doi: 10.1002/1531-8249(199901)45:1<100::aid-art16>3.0.co;2-u

Nakabayashi, H., Taketa, K., Miyano, K., Yamane, T., and Sato, J. (1982). Growth of human hepatoma cells lines with differentiated functions in chemically defined medium. *Cancer Res.* 42, 3858–3863.

Negro, F. (2010). Abnormalities of lipid metabolism in hepatitis C virus infection. *Gut* 59, 1279–1287. doi: 10.1136/gut.2009.192732

Nguyen, L. P., Tran, S. C., Suetsugu, S., Lim, Y.-S., and Hwang, S. B. (2020). PACSIN2 interacts with nonstructural protein 5A and regulates hepatitis C virus assembly. *J. Virol.* 94, e01531-19. doi: 10.1128/JVI.01531-19

Odendall, C., Dixit, E., Stavru, F., Bierne, H., Franz, K. M., Durbin, A. F., et al. (2014). Diverse intracellular pathogens activate type III interferon expression from peroxisomes. *Nat. Immunol.* 15, 717–726. doi: 10.1038/ni.2915

Okuda, M., Li, K., Beard, M. R., Showalter, L. A., Scholle, F., Lemon, S. M., et al. (2002). Mitochondrial injury, oxidative stress, and antioxidant gene expression are induced by hepatitis C virus core protein. *Gastroenterology* 122, 366–375. doi: 10.1053/gast.2002.30983

Perez-Riverol, Y., Bai, J., Bandla, C., García-Seisdedos, D., Hewapathirana, S., Kamatchinathan, S., et al. (2022). The PRIDE database resources in 2022: a hub for mass spectrometry-based proteomics evidences. *Nucleic Acids Res.* 50, D543–D552. doi: 10.1093/nar/gkab1038

Ramage, H. R., Kumar, G. R., Verschueren, E., Johnson, J. R., Dollen, V. z. J., Johnson, T., et al. (2015). A combined proteomics/genomics approach links hepatitis C virus infection with nonsense-mediated mRNA decay. *Mol. Cell* 57, 329–340. doi: 10.1016/j.molcel.2014.12.028

Ramirez, S., Mikkelsen, L. S., Gottwein, J. M., and Bukh, J. (2016). Robust HCV genotype 3a infectious cell culture system permits identification of escape variants with resistance to Sofosbuvir. *Gastroenterology* 151, 973–985.e2. doi: 10.1053/j.gastro.2016.07.013

Ran, F. A., Hsu, P. D., Wright, J., Agarwala, V., Scott, D. A., and Zhang, F. (2013). Genome engineering using the CRISPR-Cas9 system. *Nat. Protoc.* 8, 2281–2308. doi: 10.1013/nprot.2013.143

Reiss, S., Rebhan, I., Backes, P., Romero-Brey, I., Erfle, H., Matula, P., et al. (2011). Recruitment and activation of a lipid kinase by hepatitis C virus NS5A is essential for integrity of the membranous replication compartment. *Cell Host Microbe* 9, 32–45. doi: 10.1016/j.chom.2010.12.002

Rhee, H.-W., Zou, P., Udeshi, N. D., Martell, J. D., Mootha, V. K., Carr, S. A., et al. (2013). Proteomic mapping of mitochondria in living cells via spatially restricted enzymatic tagging. *Science* 339, 1328–1331. doi: 10.1126/science.1230593

Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.* 43:e47. doi: 10.1093/nar/gkv007

Riva, L., Spriet, C., Barois, N., Popescu, C.-I., Dubuisson, J., and Rouillé, Y. (2021). Comparative analysis of hepatitis C virus NS5A dynamics and localization in assembly-deficient mutants. *Pathogens* 10:172. doi: 10.3390/pathogens10020172

Rocha-Perugini, V., Lavie, M., Delgrange, D., Canton, J., Pillez, A., Potel, J., et al. (2009). The association of CD81 with tetraspanin-enriched microdomains is not essential for hepatitis C virus entry. *BMC Microbiol.* 9:111. doi: 10.1186/1471-2180-9-111

Rohart, F., Gautier, B., Singh, A., and Lê Cao, K.-A. (2017). mixOmics: an R package for 'omics feature selection and multiple data integration. *PLoS Comput. Biol.* 13:e1005752. doi: 10.1371/journal.pcbi.1005752

Romero-Brey, I., Merz, A., Chiramel, A., Lee, J. Y., Chlanda, P., Haselman, U., et al. (2012). Three-dimensional architecture and biogenesis of membrane structures associated with hepatitis C virus replication. *PLoS Pathog.* 8:e1003056. doi: 10.1371/journal.ppat.1003056

Rouillé, Y., Helle, F., Delgrange, D., Roingeard, P., Voisset, C., Blanchard, E., et al. (2006). Subcellular localization of hepatitis C virus structural proteins in a cell culture system that efficiently replicates the virus. *J. Virol.* 80, 2832–2841. doi: 10.1128/JVI.80.6.2832-2841.2006

Sahuc, M.-E., Sahli, R., Rivière, C., Pène, V., Lavie, M., Vandeputte, A., et al. (2019). Dehydrojuncusol, a natural Phenanthrene compound extracted from *Juncus maritimus* is a new inhibitor of hepatitis C virus RNA replication. *J. Virol.* 93, e02009–18. doi: 10.1128/JVI.02009-18

Shimakami, T., Yamane, D., Jangra, R. K., Kempf, B. J., Spaniel, C., Barton, D. J., et al. (2012). Stabilization of hepatitis C virus RNA by an Ago2-miR-122 complex. *Proc. Natl. Acad. Sci. U. S. A.* 109, 941–946. doi: 10.1073/pnas.1112263109

Stone, M., Jia, S., Heo, W. D., Meyer, T., and Konan, K. V. (2007). Participation of rab5, an early endosome protein, in hepatitis C virus RNA replication machinery. *J. Virol.* 81, 4551–4563. doi: 10.1128/JVI.01366-06

Sugiura, A., Mattie, S., Prudent, J., and McBride, H. M. (2017). Newly born peroxisomes are a hybrid of mitochondrial and ER-derived pre-peroxisomes. *Nature* 542, 251–254. doi: 10.1038/nature21375

Sychev, Z. E., Hu, A., DiMaio, T. A., Gitter, A., Camp, N. D., Noble, W. S., et al. (2017). Integrated systems biology analysis of KSHV latent infection reveals viral induction and reliance on peroxisome mediated lipid metabolism. *PLoS Pathog.* 13:e1006256. doi: 10.1371/journal.ppat.1006256

Tai, A. W., Benita, Y., Peng, L. F., Kim, S.-S., Sakamoto, N., Xavier, R. J., et al. (2009). A functional genomic screen identifies cellular cofactors of hepatitis C virus replication. *Cell Host Microbe* 5, 298–307. doi: 10.1016/j.chom.2009.02.001

Tardif, K. D., Mori, K., Kaufman, R. J., and Siddiqui, A. (2004). Hepatitis C virus suppresses the IRE1-XBP1 pathway of the unfolded protein response. *J. Biol. Chem.* 279, 17158–17164. doi: 10.1074/jbc.M312144200

Tu, H., Gao, L., Shi, S. T., Taylor, D. R., Yang, T., Mircheff, A. K., et al. (1999). Hepatitis C virus RNA polymerase and NS5A complex with a SNARE-like protein. *Virology* 263, 30–41. doi: 10.1006/viro.1999.9893

Tyanova, S., Temu, T., Sinitcyn, P., Carlson, A., Hein, M. Y., Geiger, T., et al. (2016). The Perseus computational platform for comprehensive analysis of (prote)omics data. *Nat. Methods* 13, 731–740. doi: 10.1038/nmeth.3901

van der Meer, A. J., Feld, J. J., Hofer, H., Almasio, P. L., Calvaruso, V., Fernández-Rodríguez, C. M., et al. (2017). Risk of cirrhosis-related complications in patients with advanced fibrosis following hepatitis C virus eradication. *J. Hepatol.* 66, 485–493. doi: 10.1016/j.jhep.2016.10.017

Virzi, A., Roca Suarez, A. A., Baumert, T. F., and Lupberger, J. (2018). Oncogenic signaling induced by HCV infection. *Viruses* 10:538. doi: 10.3390/v10100538

Wakita, T., Pietschmann, T., Kato, T., Date, T., Miyamoto, M., Zhao, Z., et al. (2005). Production of infectious hepatitis C virus in tissue culture from a cloned viral genome. *Nat. Med.* 11, 791–796. doi: 10.1038/nm1268

Wang, T., Campbell, R. V., Yi, M. K., Lemon, S. M., and Weinman, S. A. (2010). Role of hepatitis C virus core protein in viral-induced mitochondrial dysfunction. *J. Viral Hepat.* 17, 784–793. doi: 10.1111/j.1365-2893.2009.01238.x

Wang, H., Perry, J. W., Lauring, A. S., Neddermann, P., De Francesco, R., and Tai, A. W. (2014). Oxysterol-binding protein is a phosphatidylinositol 4-kinase effector required for HCV replication membrane integrity and cholesterol trafficking. *Gastroenterology* 146, e1–e11. doi: 10.1053/j.gastro.2014.02.002

WHO. (2022). *Hepatitis C*. Available at: https://www.who.int/news-room/fact-sheets/detail/hepatitis-c (Accessed June 28, 2023).

Wu, H., Carvalho, P., and Voeltz, G. K. (2018). Here, there, and everywhere: the importance of ER membrane contact sites. *Science* 361:e aan5835. doi: 10.1126/science.aan5835

Yofe, I., Soliman, K., Chuartzman, S. G., Morgan, B., Weill, U., Yifrach, E., et al. (2017). Pex35 is a regulator of peroxisome abundance. *J. Cell Sci.* 130, 791–804. doi: 10.1242/jcs.187914

You, J., Hou, S., Malik-Soni, N., Xu, Z., Kumar, A., Rachubinski, R. A., et al. (2015). Flavivirus infection impairs peroxisome biogenesis and early antiviral signaling. *J. Virol.* 89, 12349–12361. doi: 10.1128/JVI.01365-15

Zhang, H., Zhang, C., Tang, H., Gao, S., Sun, F., Yang, Y., et al. (2018). CD2-associated protein contributes to hepatitis C, virus propagation and steatosis by disrupting insulin signaling. *Hepatology* 68, 1710–1725. doi: 10.1002/hep.30073