key: cord-0079208-k2560aqw
authors: Han, Peiyu; Shen, Li; Nan, Nan; Zhou, Renwu; Li, Tianci; Dai, Xiaofeng
title: Cold Atmospheric Plasma Boosts Virus Multiplication via EGFR(Tyr1068) Phosphorylation-Mediated Control on Cell Mitophagy
date: 2022-05-09
journal: Int J Biol Sci
DOI: 10.7150/ijbs.71983
sha: 3e7d711e79ffd7f5ad2899a60bef0e2f4fc02db4
doc_id: 79208
cord_uid: k2560aqw

Objectives: Vaccination still remains as the most effective approach for preventing infectious diseases such as those caused by virus infection, with cell-based vaccine manufacturing being one flexible solution regarding the spectrum of infectious disorders it can prevent. Rapid cell-based virus propagation can enable high yield of vaccines against viral diseases that may offer critical values in the industry when handling emergent situations such as the ongoing viral disease pandemic. Methods: Through investigating the phenomenon and biological mechanism underlying redox-triggered cell survival towards enhanced viral particle production, this study explores novel strategies for improved yield of viral particles at a reduced cost to meet the increasing demand on cell-based vaccine manufacturing against viral diseases. Results: We found in this study that cold atmospheric plasma (CAP), composed of multiple reactive oxygen and nitrogen species including H(2)O(2), could effectively enhance virus replication via triggering cell mitophagy that was dynamically modulated by the p-EGFR(Tyr1068)/p-Drp1(Ser616) axis using IBRV and MDBK as the virus and cell models, respectively; and removing H(2)O(2) can further enhance virus yield via releasing cells from excessive G(0)/G(1) cell cycle arrest. The observed efficacy of CAP was extended to other viruses such as CDV and CPV. Conclusion: This study provides experimental evidences supporting the use of CAP as a modulator of cell survival including mitophagy and mitochondria dynamics, and makes CAP an interesting and promising tool for enhancing the yield of viral vaccines if translated into the industry.

Vaccination is an effective approach for infectious disease prevention and possibly the only solution for eradicating viral diseases that have low virulence but rapid transmission and with the potential of eventually evolving into symbiosis with human. Cell-based vaccine manufacturing that utilizes cultured mammalian cells as the virus production factory represents a flexible approach for viral vaccine production regarding the yield and spectrum of viral diseases it can control [1] . Rapid vaccine production can deliver critical benefits especially under urgent situations such as the outbreak of Coronavirus Disease 2019 (COVID-19) epidemic [2] . Enhancing virus titer is a practical strategy to meet the large and increasing demand on cell-based vaccine manufacturing against viral diseases.

Cold atmospheric plasma (CAP), belonging to the fourth state of matter, has been shown capable of preventing COVID-19 [3] , deactivating viruses, removing pathogens such as bacteria and fungi [4] , and selectively killing cancer cells via inducing, e.g.,

International Publisher cell cycle arrest [5] , senescence [6] , autophagy [7] , apoptosis [8] , ferroptosis [9, 10] , ICD [11] , and necrosis [12] in a dose-dependent fashion [13] . It is composed of long-lived and short-lived reactive oxygen and nitrogen species, the chemical properties of which have been intensively studied [14] [15] [16] [17] [18] [19] . We had previously shown that CAP could hijack host signaling towards enhanced virus production under controlled dosing and following specified protocol [7, 20] . It is thus possible that CAP can function as a tool to boost virus multiplication that ultimately leads to rapid vaccine production for viral disease control.

Motivated as such, we explored the efficacy of CAP in triggering increased virus titration and the underlying biological mechanism with the focus laid on cell survival including selective autophagy and cell cycle arrest. Using infectious bovine rhinotrachieitis virus (IBRV) and Madin-Darby Bovine Kidney (MDBK) cells as the virus and cell models, respectively, we found that CAP could boost virus propagation through enhancing cell mitophagy that was dynamically modulated by the activity and location of phosphorylated EGFR (Tyr1068). Interestingly, removing H2O2 from CAP-activated medium (PAM) further enhanced virus titer via releasing cells from H 2 O 2 -triggered halt at the G 0 /G 1 stage. Similar effect of CAP on boosting virus titer was observed for other viruses such as canine parvovirus (CPV) and canine distemper virus (CDV). Our study makes CAP an effective approach enabling rapid viral vaccine manufacturing and, also a potential dose-dependent tool for modulating cell mitophagy and mitochondria dynamics towards desirable cell outcomes for other research and clinical purposes.

The home-made CAP ejection device is composed of a power controller, a helium (He) gas cylinder, a rotor flow meter, and a CAP jet (Supplementary Figure 1) . The peak-to-peak voltage applied to the electrode, the sinusoidal wave frequency, the He flow rate, and the distance from the CAP jet to the medium surface were set to 0.96-1.24 kV, 10 kHz, 1 L/min, and 13 mm, respectively. Plasma activated medium (PAM) was prepared by exposing 1 mL medium to CAP for 2 or 4 min for each well in a 24-well plate.

IBRV (NCBI TaxId: 79889) was purchased from the Type Culture Collection of the Chinese Academy of Sciences, Shanghai, China and propagated using MDBK cells. The MDBK cell line was purchased from Cell Bank of the Chinese Academy of Sciences Cell Bank and cultured in RPMI medium containing 10% fetal bovine serum (FBS) and 1% penicillinstreptomycin.

Cells were grown to 80%-90% confluence prior to incubation with different doses of PAM for 1 h. Cells were infected with viruses for 1 h. The medium was refreshed, and cells were cocultured with viruses for additional 46 h.

Tissue culture infective dose (TCID50) was measured as the virus dilution fold when 50% cells were infected. Cells were plated in a 96-well plate. Inoculating cells with a 10-gradient sequentially diluted solutions of viral fluid at a step-size of 10 folds, and monitoring cells for 5 to 7 days until the occurrence of cytopathic effect (CPE) in 50% cells. TCID50 was calculated following previously recommended protocol [21] .

Basic qRT-PCR protocol: 4 μL cDNA samples, 10 μL 2×UltraSYBR MixTure, 1 μL forward and 1 μL backward primers, 4 μL ddH 2 O were mixed and centrifuged before running the qRT-PCR program (pre-denature at 95 °C for 10 min; 95 °C for 10 s, 60 °C for 1 min, 72 °C for 20 s, for 40 cycles) in Roche LightCycler 480 RT-PCR. Each sample had 3 replicates.

For virus detection, supernatants were collected and boiled following the protocol of UltraSYBR MixTure Kit (Catalog Number: CW0957M, Cwbio Co. Ltd., China) prior to the aforementioned basic q-PCR protocol.

For EGFR mRNA quantification, intracellular mRNA was extracted followed by reverse transcription into cDNA prior to the aforementioned basic q-PCR protocol.

Primers for IBRV detection and EGFR quantification were listed in Supplementary Table 1 .

Cells were grown in 6-well plates, washed twice with PBS, lysed on ice for 20 min using RIPA lysis buffer containing protease inhibitor, and centrifuged at 12,000 g for 20 min before collecting the supernatants. Protein concentration was estimated using BCA analysis kit (Catalog Number: P0012 Beyotime Biotechnology, China). The protein per lane (20 μg) was resolved by SDS-PAGE and transferred to a PVDF membrane. After blocking with 5% skim milk powder in TBS and Tween-20 buffer, the membrane was incubated overnight at 4 °C with a suitable primary antibody followed by incubation with the secondary antibody for 1 h at the room temperature. After TBST cleaning, the membrane was detected by Bio-Rad Gel Doc XR. Details on antibodies used in this study were summarized in Supplementary Table 2 .

Cells were fixed with cell fixature and permeabilized with 1% Triton X-100. Cells were blocked in 1% bovine serum albumin (BSA) followed by incubation with the primary antibodies for 8 h. Cells were incubated with the secondary antibody and imaged using laser-scanning confocal microscope Zeiss Axio Imager Z2 (Zeiss, Germany).

For mitochondria morphology imaging, cells were incubated with 100 nM MitoTracker (Catalog Number: M7512 Thermo Scientific, USA) for 30 min at 37 °C prior to the performance of the aforementioned protocol.

The pRP[CRISPR]-EGFP/Puro-hCas9-U6> {EGFR[gRNA]} and Single-stranded Oligo Donor (ssODN) (CAGGGATCCTGCATGGGATGGTGCTTT GCTGATTACTTCACCTCTGATTTCTTTCCACTTTC AGAGTTCATTAATCAGTCCGTTCCCAAAAGGCC CGCTGGCTCTGTGCAGAATCCTGTCTATCACAA TCAGCCTCT) were purchased from Vector Builder. The ssODN-EGFR was ethanol precipitated and re-suspended (100 μM) in 10 mM Tris-HCl (pH=7.5). Oligos of EGFR(Tyr1068) (not including the signal peptide of 24 Aa at N-terminal) were diluted to 1 μM prior to transfection. When MDBK cells grew to 60% -70% confluence in 6-well plates, 1 μM ssODN, 1 μg pRP[CRISPR]-EGFP/Puro-hCas9-U6 and 6 μL lipofectamine 3000 were added followed by incubation for 8 h.

Cells were grown in 6-well plates, washed using PBS and digested with EDTA-free trypsin. Cells were centrifuged at 1000 rpm for 5 min to retain the pellets. Cell pellets were washed with 500 μL PBS, centrifuged at 1000 rpm/min for 5 min, and the supernatant was removed to retain the pellet. Cell pellets were re-suspended in 70% ethanol and placed in a refrigerator at 4 °C overnight for fixation. Fixed cells were centrifuged at 1000 rpm for 5 min to remove the supernatant, and cell pellets were suspended in 500 μL PBS. 5 μL propidium iodide (Catalog Number: ST511, Beyotime Biotechnology, China) was added and mixed with cells on ice for 30 min. Cell cycle detection was performed using a BD Accuri C6 flow cytometer, and data analysis was performed using Flowjo software.

DCFH-DA (Catalog Number: S0033S, Beyotime Biotechnology, China) was diluted with serum-free culture medium at a ratio of 1:1000 to make the final concentration 10 mM/L. Cell culture medium was removed followed by addition of DCFH-DA in an appropriate volume. Cells were incubated in a 37 °C cell incubator for 20 min. Cells were washed 3 times using serum-free cell culture medium to fully remove the rest of DCFH-DA. Cells were collected by trypsin without adding EDTA. Real-time fluorescence intensity detection was performed using Zeiss Axio Imager Z2 (Zeiss, Germany) before and after the signal stimulation using 488 nm excitation wavelength and 525 nm emission wavelength. After treating cells with CAP for 2 or 4 min followed by incubation for 1 h, the quencher of each CAP component was added, separately. Viruses were added to cells and incubated for 1 h. The medium was refreshed, and cells were cocultured with viruses for additional 46 h prior to subsequent analyses.

The optical emission spectrum, in the gas phase, of CAP ejected from our home-made He-based device (Supplementary Figure 1 ) was detected using a spectrometer (Andor Shamrock SR-500i-A-R, England), where the optical probe was placed at a distance of approximately 10 mm from the center of the ejection device.

Gene expression, represented as FPKM (Fragments Per Kilo Per Million reads) [22] , was computed using Htseq software (V 0.6.1) following Equation 1.

The mapped reads represent how much percent of all reads are mapped to a gene divided by the length of the gene. FPKM normalization normalizes both the sequencing depth and the gene length to make the expression of genes at different lengths and different sequencing depths comparable. In general, FPKM value 0.1 or 1 is considered as the threshold to determine whether a gene is expressed.

The correlation of gene expression levels between samples is an important index to test the reliability of the experiments. Pearson correlation was conducted among samples in a pair-wise fashion, with R 2 >0.8 being considered highly correlated.

Principle component analysis (PCA) was conducted using R packages 'gmodels', 'ggplot2', 'ggrepel'.

A hierarchical clustering analysis was performed based on the FPKM values of all genes in different samples, where 'euclidean' was used as the distance and 'complete linkage' was set as the linkage approach.

Differentially expressed genes were identified using the 'DESeq2' bioconductor package, where genes differentially expressed over 2 folds between the compared groups were considered as differentially expressed, and those with adjusted p value ≤ 0.05 and FDR ≤ 0.05 were selected.

Gene ontology (GO) analysis was conducted using GOseq [23] that is based on Wallenius non-central hyper-geometric distribution. Threshold used for enrichment filtering was set as p ≤ 0.05 for over represented genes.

KEGG (Kyoto Encyclopedia of Genes and Genomes) [24] analysis was conducted using the 'clusterProfiler' package [25] from R. The analysis was conducted using hypergeometric test where KEGG pathway was considered as the unit following Equation 2.

In Equation 2 , 'N' is the number of genes annotated in the pathways among all genes, 'n' is the number of differentially expressed genes among 'N', 'M' is the number of genes annotated to the specific pathway among all genes, and 'm' is the number genes annotated to the specific pathway among differentially expressed genes. Threshold for enrichment significance was set as p value ≤ 0.05.

Through activating the cell culture medium (Supplementary Figure 1) , CAP could substantially enhance the number of cells undergoing CPE after IBRV infection ( Figure 1A ) and effectively increase IBRV titer from 10 -2.57 /mL to 10 -3 /mL on 2 min CAP exposure and to 10 -4.5 /mL on 4 min exposure ( Figure  1B) . The qRT-PCR also revealed reduced PCR cycles required for IBRV detection on 2 min and 4 min CAP exposure ( Figure 1C ).

CAP is featured by dose-dependent cellular redox level modulation [8, 13, 26, 27] . Substantially enhanced level of ROS was observed on 4-min CAP exposure as compared with 2-min CAP exposure in IBRV-infected MDBK cells ( Figure 1D ).

Selective autophagy is one important process governing virus replication [20] and constitutes to an important part of cell stress responses including redox stress in mammalian cells [28] . P62 is the best characterized selective autophagy substrate in mammals that localizes at the autophagosome formation site [29] , and LC3B is responsible for selective degradation of p62 through selective autophagy [30] , both of which are typical markers of selective autophagy [31] . We observed considerably elevated LC3B and substantially reduced p62 expression on joint CAP exposure and IBRV infection ( Figure 1E) . The fluorescence imaging signal of LC3B in cells receiving both CAP exposure and IBRV infection was considerably higher than that in cells treated with CAP or IBRV alone ( Figure 1F) . These results collectively suggested the accompany of selective autophagy in CAP-triggered increase on virus titer.

Various evidences have supported the roles of EGFR in ROS-induced cellular response [32] [33] [34] [35] [36] [37] . EGFR level was considerably reduced when cells were exposed to IBRV infection and 4-min CAP treatment followed by 46 h culturing (p<0.001, Figure 2A ). EGFR phosphorylation sites (Tyr1068, Tyr1148, Tyr1086, Tyr845) are highly conserved between human and bovine protein sequences (Figure 2B) , enabling the use of human rather than bovine antibodies in the following experiments here.

For all the four phosphorylation sites examined, p-EGFR (Tyr1068) was synergistically enhanced by IBRV infection and CAP exposure ( Figure 2C) . With the increase of CAP dose (from 2 min to 4 min) that did not trigger cell death (Figure 2D) , p-EGFR(Tyr1068) started to accumulate in cell nucleus ( Figure 2E ). IJoint exposure of cells to CAP and IBRV infection substantially enhanced p-EGFR(Tyr1068) level that was co-localized with mitochondria ( Figure 2F) , associating observed selective autophagy to mitophagy.

We next examined the role of p-EGFR(Tyr1068) in CAP-triggered mitophagy after IBRV infection. Mitochondria accumulated to the nucleus surrounding area, and the signal was intensified with the increase of CAP treatment duration (Figure 3A) , implicating the occurrence of mitophagy. Mitofusin 1 (Mfn1) and dynamin-related protein 1 (Drp1) coordinate the processes of mitochondrial fusion and fission, which are canonical markers of mitophagy [38, 39] . On joint exposure to IBRV infection and CAP treatment, p-EGFR(Tyr1068) and p-Drp1(Ser616) expression were substantially increased, whereas Mfn1 level was reduced (Figure 3B) , suggesting the occurrence of cell mitophagy and involvement of p-EGFR(Tyr1068) in cells' response to both CAP exposure and IBRV infection. LC3B was co-localized with mitochondria and the signal was intensified on joint exposure to CAP and IBRV infection (Figure 3C) , confirming the presence of mitophagy. Co-immunoprecipitation suggested physical interactions of p-EGFR(Tyr1068) with p-Drp1(Ser616) but not with Mfn1 on joint exposure to IBRV infection and CAP treatment (Figure 3D) , indicating that the enhanced cell fission was mediated via the p-EGFR(Tyr1068)/ p-Drp1(Ser616) axis. 

To consolidate the role of p-EGFR(Tyr1068) in CAP-triggered increase on virus titer and mitophagy after IBRV infection, we generated an EGFR(Tyr1068) mutant plasmid where tyrosine (Tyr) at the 1068 site was mutated to phenylalanine (Phe) (Figure 4A) . The p-EGFR(Tyr1068) level was successfully reduced in cells transfected with the mutant plasmid (Figure 4E,  4F) . Cells transfected with the mutant plasmid produced considerably less CPE after IBRV infection with and without CAP exposure ( Figure 4B) . Accordingly, virus titer reduced substantially ( Figure  4C ) and the Ct values from the qRT-PCR assay were significantly enhanced in EGFR(Tyr1068) mutant cells ( Figure 4D) . Each primary marker used for selective autophagy and mitophagy characterization (i.e., LC3B, p-Drp1(Ser616), Mfn1) exhibited a relapsed pattern in EGFR(Tyr1068) mutant cells as compared with the wild type under the same treatment ( Figure  4E, 4G, 4H) .

CAP ejected from our home-made He-based device (Supplementary Figure 1) was primarily composed of OH, N 2 , N 2 + and O in the gas phase (Figure 5A) , which became H 2 O 2 , O 3 , •O 2 -, •OH, •NO, ein the liquid form [40] (Figure 5B) . By removing/quenching H 2 O 2 from PAM, cell amount ( Figure 5C ) and virus yield ( Figure 5D ) were substantially increased after PAM treatment. We next confirmed the roles of H 2 O 2 in decreasing virus titer using TCID50, i.e., it dropped from 10 -5 to 10 -5.33 in IBRV-infected cells and from 10 -6.625 to <10 -10 in CPV-infected cells on CAP exposure (Figure 5E, 5F) .

The p-EGFR(Tyr1068) mediated virus titer increase was primarily triggered via •NO and egiven the dramatic decrease on p-EGFR(Tyr1068) expression after removing each of these species (Figure 5G, 5H) . The effect of H 2 O 2 on virus titer was in-dependent of p-EGFR(Tyr1068) signaling as no visible alteration was observed regarding p-EGFR(Tyr1068) expression when H 2 O 2 was removed ( Figure 5G, 5H) , suggestive of the involvement of additional mechanisms in virus titer control. Autophagy as assessed via p62 and LC3B expression here was primarily affected by •O 2 -, •OH, •NO that, once being removed separately, showed opposite profiles on CAP exposure (Figure 5I, 5J) . Mitochondria fission through p-Drp1(Ser616) and fusion via Mfn1 were not primarily affected by any of these components as the fluctuations were not dramatic by adding any of these reactive species quenchers (Figure 5G, 5K, 5L) , suggestive of the combinatorial role of these species in triggering cell mitophagy.

We next conducted transcriptome sequencing to explore the mechanism of H 2 O 2 in decreasing virus titer and gain a broad view on CAP-triggered cell response. Gene expression did not considerably vary among sample groups (Supplementary Figure 2A) , and those of samples within the same group were highly correlated (Supplementary Figure 2B) . Samples of the same group were clustered together from both PCA and clustering analyses (Figures 6A,  6B) . CAP treatment caused rightward movement of samples along the first principal component (PC1) that explained 27.4% of the total variance; and removing H2O2 from CAP leftward shifted samples along PC1 and downward along PC2 that explained 15.8% of the total variance ( Figure 6A ). If PC1 reflected the effect of mitophagy, PC2 might represent another cell behavior (namely 'behavior X') that was affected by CAP and, in particular, H2O2. Additionally, removing H 2 O 2 decreased the expression of under-expressed (0-0.1) genes (p=0.0129 for IBRV+CAP-H 2 O 2 vs. IBRV; p=0.0019 for IBRV+CAP-H 2 O 2 vs. IBRV+CAP) and caused relapsed expression of 1-3 level genes (p=0.0046 for IBRV+CAP-H 2 O 2 vs. IBRV+CAP) with statistical significance that was oppositely regulated by CAP ( Figure 6C) . These are suggestive of the controversial roles of H 2 O 2 and CAP in modulating virus titer, as well as the association of these altered genes with the 'behavior X'. Specifically, 420, 168, 78 genes were differentially expressed in the pairwise comparisons 'IBRV vs IBRV+CAP', 'IBRV vs IBRV+CAP-H 2 O 2 ' and 'IBRV+CAP vs IBRV+CAP-H 2 O 2 '; and genes simultaneously modulated in all of these three group-wise comparisons showed a consistent regulatory direction in 'IBRV vs IBRV+CAP', 'IBRV vs IBRV+CAP-H 2 O 2 ' but an opposite direction in 'IBRV+CAP vs IBRV+CAP-H 2 O 2 ' (Figure 6D,  Supplementary Table 3, Supplementary Figure 3A) . These results suggested that CAP was the primary contributor to the altered transcriptome profile, and H 2 O 2 played an opposite role.

Gene ontology (GO) and KEGG pathway analyses revealed that CAP-triggered cell behavior was associated with 'extracellular space' and 'EGFR tyrosine kinase inhibitor resistance', which were the top GO term and pathway, respectively, from both 'IBRV vs IBRV+CAP' and 'IBRV vs IBRV+CAP-H2O2' group-wise comparisons, but not from the 'IBRV+CAP vs IBRV+CAP-H 2 O 2 ' comparison ( Figure  6E, 6F, Supplementary Figure 3B, 3C) . These results are supportive to the mediating role of p-EGFR(Tyr1068) in CAP-triggered mitophagy (Figure 4) since the extracellular region of EGFR has been previously proposed as a potential target for anti-EGFR drug discovery [41] . The top GO term and KEGG pathway identified from the 'IBRV+CAP vs IBRV+CAP-H2O2' comparison but not from the other two group-wise comparisons was 'lysosome' ( Figure  6G, 6H, Supplementary Figure 3B, 3C) , which is known playing essential roles in autophagy [42] . Given the differential regulation of autophagy during cell cycle progression [43] and our knowledge on the roles of H2O2 in triggering cellular damage and consequently cell cycle arrest under certain concentrations [44] , we hypothesized H 2 O 2 -induced cellular 'behavior X' to be 'cell cycle alteration'.

We firstly explored CAP-imposed H 2 O 2 concentration in the medium (or precisely, quenched H 2 O 2 amount in the scavenger assay) via removing H 2 O 2 first using 10 mM sodium pyruvate followed by H 2 O 2 titration, the lowest Ct value of IBRV infected cells was achieved when 5 mM and 1 mM H 2 O 2 were supplemented, respectively, with and without adding sodium pyruvate (Figure 7A, 7B) ; thus, H 2 O 2 quenched by 10 mM sodium pyruvate in our (Figure 7C) . On the other hand, virus titer decreased when cellular H 2 O 2 exceeded 5 mM ( Figure  7B) , suggesting that the intracellular H 2 O 2 raised in response to virus infection was over 1 mM (5 mM -4 mM). Either virus-triggered H 2 O 2 (over 1 mM) or CAP-imposed H 2 O 2 (approximately 4 mM) led to moderate cell cycle arrest (around 40% G 0 /G 1 phase) that corresponded to slightly increased virus titer. Both sources of H 2 O 2 , once added together, dramatically decreased virus yield as a result of substantially increased number of cells arrested at the G 0 /G 1 stage (Figure 7D) . Cellular ROS increased with H 2 O 2 dose, and the intensity dramatically dropped at 10 mM (Figure 7E) , indicative of cell death due to exceeded ROS level. Given these results, together with the known functionality of ROS in triggering cell cycle arrest [45] , we assumed that the observed inhibition on virus titer was due to exceeded level of H2O2.

The levels of p-EGFR(Tyr1068) and p-Drp1(Ser616) were both slightly up-regulated when H 2 O 2 was low (< 1 mM) and dramatically dropped under high (1-10mM) level of H 2 O 2 ( Figure 7F,  Supplementary Figure 4) , which were supportive of the physical interactions of both phosphorylated proteins ( Figure 3D ) and the opposite role of high dose (1-10mM) of H 2 O 2 and CAP in triggering p-EGFR(Tyr1068)-mediated cell response. Autophagy markers p62 and LC3B showed a similar pattern with that of p-EGFR and p-Drp1 ( Figure 7F) . Cyclin B1 is primarily distributed in the M phase and serves as the master regulator of the G2/M transition [46] . CDK1 interacts with 9 different cyclins including cyclin B1 that becomes inactivated at the end of the M phase [47] . Both cyclin B1 and CDK1 expression declined with the increased dose of H 2 O 2 ( Figure 7G,  Supplementary Figure 4) , suggestive of a suppressive role of high H 2 O 2 dose on the M stage. Cyclin D1 is rapidly synthesized in the G1 phase and degraded as cells enter the S phase [48] , and CDK4/6 creates complexes with cyclin D1 to mediate cell progression through the G1 phase [49] . Cyclin D1 expression showed a positive correlation with H2O2 dose; CDK4 expression increased with H 2 O 2 concentration and declined when H 2 O 2 exceeded 1 mM (Figure 7G,  Supplementary Figure 4) , implicating that cells were accumulated and arrested at the G 0 /G 1 stage when H 2 O 2 exceeded 1 mM.

As integrin represents a common receptor mediating virus entry [50] and is linked to intracellular signaling via EGFR tyrosine activation, we focused on the four integrin-dependent EGFR tyrosine phosphorylation residues, i.e., Tyr1068, Tyr1086, Try845, Tyr1173 [51] in this study. We excluded Tyr1173 given its primary role in the MAPK cascade and cell proliferation [52] to narrow our focus down to cell selective autophagy/mitophagy. Among the three EGFR tyrosine phosphorylation sites, EGFR(Tyr1068) showed the highest phosphorylation level in response to joint IBRV infection and CAP treatment (Figure 2C) . The expression profile of p-EGFR(Tyr1068) under varied CAP doses matched well with that of selective autophagy biomarkers (LC3B, p62, Figure 1E ), in particular mitophagy markers (p-Drp1(Ser616), Mfn1, Figure 3B) , and virus titer (Figure 1A-1D) . In addition, p-EGFR(Tyr1068) co-localized with mitochondria ( Figure 2F ). These were, collectively, implicative of the regulatory role of p-EGFR(Tyr1068) in mitophagy. We, therefore, considered p-EGFR(Tyr1068) as the key EGFR phosphorylation site that drove the observed efficacy of CAP in boosting virus titer.

Selective autophagy degrades certain proteins, organelles to maintain cellular homeostasis and, on virus invasion, can be hijacked to recycle cell nutritional resources for virus propagation [20] . Elevated LC3B expression and decreased p62 level on CAP exposure after virus infection (Figure 1E, 1F) attributed the promotive role of CAP on virus multiplication to cell selective autophagy. Mitochondria accumulation at the nucleus surrounding region (Figure 3A) , enhanced Drp1(Ser616) phosphorylation and reduced Mfn1 level (Figure 3B, 3C) , providing additional evidence associating CAP-triggered selective autophagy to mitophagy. Mitochondrial fusion, fission, biogenesis and mitophagy constitute mitochondrial dynamics that determines mitochondrial morphology, quality and abundance [53, 54] . Activated mitochondrial fusion as represented by elevated Mfn1 [38] is protective of mitochondria on pressure. Mitochondrial fission machinery as featured by Drp1 [38] actively participates in apoptosis [55] and promotes or inhibits cell death depending on the initial lethal stimulus [56] . We observed physical interactions between phosphorylated p-EGFR(Tyr1086) and p-Drp1 (Ser616) on IBRV infection and/or CAP exposure (Figure 3D) , suggesting that mitochondria fission was enhanced due to elevated level of Drp1 that was possibly increased by activated EGFR(Tyr1068). Through introducing the EGFR(Tyr1068) point mutation to native EGFR genes using the CRISPR/Cas9 system (Figure 4) , we consolidated the role of p-EGFR(Tyr1068) in mediating CAPinduced virus titer increase. To this end, we could preliminarily conclude that CAP promoted IBRV titration via, at least partially, the trigger of cell mitophagy, and this process was driven by p-EGFR(Tyr1068)/p-Drp1(Ser616)-mediated mitochondria fission. We also observed elevated expression of the pro-apoptotic protein caspase 7 (Supplementary Figure 5) that was supportive to the pro-apoptotic role of mitochondria fission. Quenching each component of PAM using chemical scavengers, though being a qualitative approach, has been accepted as a standard measure to examine CAP-liquid interactions and the effect of CAP ingredients on cell behavior, based on the fact that each RONS scavenger has a higher reaction rate (often in an order of magnitude) with its targeted species than with other species [57, 58] . Here, we used tiron (20 mM), D-mannitol (200 mM), hemoglobin (20 µM), sodium pyruvate (10 mM), uric acid (100 µM) and monopotassium phosphate (1 mM) to trap ·O 2 -, ·OH, ·NO, H 2 O 2 , O 3 and e -, respectively, with the reaction rates being 2.9×10 9 mol/Ls, 1×10 9 mol/Ls, 1×10 8 mol/Ls, 2.4 mol/Ls, 1.4×10 6 mol/Ls [58, 59] . Except for H 2 O 2 whose removal led to a substantial increase on virus titer, eliminating each of the other ingredients unanimously decreased virus multiplication to some extent (Figure 5D ), suggestive of their indispensable and collective role in promoting virus multiplication.

Exceeded level of H2O2 (1-10 mM) triggered G 0 /G cell cycle arrest that was associated with recessed cell activities and reduced virus production (Figure 7) . Thus, cells underwent mitophagy in acute response to CAP-triggered redox stress, leading to the creation of a favorable environment for virus multiplication as a result of shrunken cell biomass usage. However, high H2O2 dose (1-10 mM) may result in cell cycle G 0 /G 1 arrest due to the trigger of moderate level of DNA damage signaling (not reaching the threshold of cell apoptosis) and thus counteract virus production. Therefore, a favorable environment for virus multiplication requires enhanced available cellular biomass materials and normal/accelerated cell cycle. Importantly, these results suggested the irreplaceable role of CAP in boosting virus titer that could not be surrogated by H2O2, an easily produced CAP component, and H 2 O 2 might function as an opposite regulator of CAP that made its effect on virus titer control reversible.

CAP is a known cocktail of ROS and non-ROS [40] . Both CAP exposure and virus infection enhanced cellular redox level [26, 60] (Figure 1D) , suggesting that the promotive role of CAP in IBRV titer is due to its synergy with virus infection in enhancing cellular ROS level. We observed similar expression profiles between p-EGFR(Tyr1068) and FOXO1 (Supplementary Figure 5 ) on CAP exposure, where FOXO1 is a key measure of cells' anti-oxidative ability [61] . This implicated the mediating role of p-EGFR(Tyr1068) in enhancing cellular anti-oxidant ability in response to CAP treatment and/or virus infection.

We additionally found that blocking EGFR(Tyr1068) phosphorylation reduced caspase 3 expression that was elevated on CAP exposure or its joint CAP and IBRV treatments (Supplementary Figure 6) . This suggested the essential role of p-EGFR(Tyr1068) in regulating cell apoptosis, and that elevated cellular ROS, on CAP exposure, triggered mitophagy and concomitantly tilted cells towards the pro-apoptotic state. This warranted us to keep appropriate CAP dosing to achieve mild yet effective boost of virus titer as over-dosing might kill the host that ultimately led to reduced virus production. We also observed rescued p65 expression in p-EGFR(Tyr1068) mutant cells on CAP or CAP+IBRV treatments (Supplementary Figure 6) . This implicated that the cellular system was hijacked by viruses for maximum virus production and other cellular behaviors including cell migration (p65 is a key protein controlling cell migration [62] ) were reduced, whereas p-EGFR(Tyr1068) was the hub controlling these processes.

Though we claimed on the driving role of p-EGFR(Tyr1068) and cell mitophagy in CAPmediated boost of IBRV propagation, we could not exclude the possible involvement of p-EGFR(Try1086) and p-EGFR(Tyr845) in this process. However, these functionalities were more likely to be associated with cells' natural response to IBRV infection rather than CAP exposure given our results.

Though different viruses may attach and/or trigger cellular signaling using different cell surface receptors beyond EGFR that may limit the generality of the uncovered driving mechanism, the facts that virus infection imposes cells with redox pressure and CAP relies on ROS to deliver its efficacy do not vary. In consistent with this, we observed enhanced titer on CAP exposure when CPV and CDV were used (Supplementary Figure 7) . Therefore, the efficacy of CAP in boosting virus production may be generalizable to a broad spectrum of, if not all, viruses and hosts that deserves intensive experimental explorations and additional consolidations.

We report in this study that CAP can enhance IBRV propagation towards increased yield of vaccines against viral diseases through triggering cell selective autophagy and, in particular, mitophagy; this process is driven by the p-EGFR (Tyr1068)/p-Drp1(Ser616) signaling through physical interactions in response to CAP-induced cellular ROS level in a dose-dependent manner (Supplementary Figure 7A) . High level of H2O2 (1-10 mM) may decrease virus titer due to enhanced G 0 /G 1 cell cycle arrest. Thus, controlling H 2 O 2 below 1 mM through the use of H 2 O 2 quencher (such as sodium pyruvate) together with CAP may represent a good strategy towards maximal virus production (Supplementary Figure 7B) . Our results contribute in advancing our understandings on CAP-triggered cellular responses that make CAP a potential modulator on cell survival such as mitophagy and cell cycle arrest. Importantly, this study provides a novel strategy for boosting IBRV multiplication that can be potentially generalized to a broad spectrum of viruses towards enhanced cell-based viral vaccine production in the industry. The generality has been preliminarily examined using CDV and CPV, which needs to be validated using more virus and cell models.

Cell culture-derived influenza vaccines from Vero cells: a new horizon for vac-cine production

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and coronavirus disease-2019 (COVID-19): The epidemic and the challenges

Cold atmospheric plasma for preventing infection of viruses that use ACE2 for entry

Vulnerability of SARS-CoV-2 and PR8 H1N1 virus to cold atmospheric plasma activated media

Cold atmospheric plasma selectively induces G0/G1 cell cycle arrest and apoptosis in AR-independent prostate cancer cells

Cold atmospheric plasma causes a calcium influx in melanoma cells triggering CAP-induced senescence

Cold atmospheric plasma increases IBRV titer in MDBK cells by orchestrating the host cell network

Cold atmospheric plasma conveys selectivity on triple negative breast cancer cells both in vitro and in vivo

Non-thermal plasma as a simple ferroptosis inducer in cancer cells: A possible role of ferritin

) expression confers intrinsic resistance to physical plasma treatment in tumor cells

Non-Thermal Plasma as a Unique Delivery System of Short-Lived Reactive Oxygen and Nitrogen Species for Immunogenic Cell Death in Melanoma Cells

Cold Atmospheric Plasma Induces a Predominantly Necrotic Cell Death via the Microenvironment

Dosing: the key to precision plasma oncology

Gas-phase peroxynitrite generation using dielectric barrier discharge at atmospheric pressure: A prospective sterilizer

The chemical product mode transition of the air DBD driven by AC power: A plausible evaluation parameter and the chemical behaviors

The hydrophilicity improvement of polytetrafluoroethylene by Ar plasma jet: The relationship of hydrophilicity, ambient humidity and plasma parameters

Long-lived species in plasma-activated water generated by an AC multi-needle-to-water discharge: effects of gas flow on chemical reactions

Insights into generation of OH radicals in plasma jets with constant power: The effects of driving voltage and frequency

Study on the Physical and Chemical Characteristics of DBD: The Effect of N2/O2 Mixture Ratio on the Product Regulation

Orchestrated efforts on host network hijacking: Processes governing virus replication

Canine parvovirus induces G1/S cell cycle arrest that involves EGFR Tyr1086 phosphorylation

Mapping and quantifying mammalian transcriptomes by RNA-Seq

Gene ontology analysis for RNA-seq: accounting for selection bias

KEGG: kyoto encyclopedia of genes and genomes

clusterProfiler: an R package for comparing biological themes among gene clusters

The Emerging Role of Gas Plasma in Oncotherapy

Cold atmospheric plasma induces GSDME-dependent pyroptotic signaling pathway via ROS generation in tumor cells

Autophagy and the integrated stress response

The role of the selective adaptor p62 and ubiquitin-like proteins in autophagy

LC3B is indispensable for selective autophagy of p62 but not basal autophagy. Biochemical and biophysical research communications

Quantifying autophagy using novel LC3B and p62 TR-FRET assays

Inhibition of the ROS-EGFR Pathway Mediates the Protective Action of Nox1/4 Inhibitor GKT137831 against Hypertensive Cardiac Hypertrophy via Suppressing Cardiac Inflammation and Activation of Akt and ERK1/2

Epidermal Growth Factor (EGF) Augments the Invasive Potential of Human Glioblastoma Multiforme Cells via the Activation of Collaborative EGFR/ ROS-Dependent Signaling

Licochalcone D Induces ROS-Dependent Apoptosis in Gefitinib-Sensitive or Resistant Lung Cancer Cells by Targeting EGFR and MET

EGFR mutation decreases FDG uptake in nonsmall cell lung cancer via the NOX4/ROS/GLUT1 axis

ROS-Mediated Mitochondrial Pathway is Required for Manilkara Zapota (L.) P. Royen Leaf Methanol Extract Inducing Apoptosis in the Modulation of Caspase Activation and EGFR/NF-kappaB Activities of HeLa Human Cervical Cancer Cells

ROS mediated EGFR/MEK/ ERK/HIF-1alpha Loop Regulates Glucose metabolism in pancreatic cancer

MFN1 deacetylation activates adaptive mitochondrial fusion and protects metabolically challenged mitochondria

Zinc accumulation in mitochondria promotes ischemia-induced BBB disruption through Drp1-dependent mitochondria fission

Genome-Wide Comparison of the Target Genes of the Reactive Oxygen Species and Non-Reactive Oxygen Species Constituents of Cold Atmospheric Plasma in Cancer Cells

Extracellular region of epidermal growth factor receptor: a potential target for anti-EGFR drug discovery

Lysosome biology in autophagy

The Role of Hydrogen Peroxide and Peroxiredoxins throughout the Cell Cycle

Aflatoxin B1 Induces Neurotoxicity through Reactive Oxygen Species Generation, DNA Damage, Apoptosis, and S-Phase Cell Cycle Arrest

Quantitative Imaging of B1 Cyclin Expression Across the Cell Cycle Using Green Fluorescent Protein Tagging and Epifluorescence

An overview of Cdk1-controlled targets and processes

Cyclin D1 in Cancer: A Molecular Connection for Cell Cycle Control, Adhesion and Invasion in Tumor and Stroma

CDK4: A Key Player in the Cell Cycle, Development, and Cancer

Host Receptors: The Key to Establishing Cells with Broad Viral Tropism for Vaccine Production

The cooperative action of bacterial fibronectin-binding proteins and secreted proteins promote maximal Campylobacter jejuni invasion of host cells by stimulating membrane ruffling

Requirement of Tyr-992 and Tyr-1173 in phosphorylation of the epidermal growth factor receptor by ionizing radiation and modulation by SHP2

Dynamics of mitochondria in living cells: shape changes, dislocations, fusion, and fission of mitochondria

Altered fusion dynamics underlie unique morphological changes in mitochondria during hypoxia-reoxygenation stress

Mitochondrial fission in apoptosis

Mitochondrial fusion and fission in the control of apoptosis

Enhanced catalytic performance of graphene-TiO2 nanocomposites for synergetic degradation of fluoroquinolone antibiotic in pulsed discharge plasma system

Inactivation of virus in solution by cold atmospheric pressure plasma: identification of chemical inactivation pathways

The oxidation of tiron by superoxide anion. Kinetics of the reaction in aqueous solution in chloroplasts

Epstein-Barr virus infection induces bone resorption in apical periodontitis via increased production of reactive oxygen species

Aconitase 2 inhibits the proliferation of MCF-7 cells promoting mitochondrial oxidative metabolism and ROS/FoxO1-mediated autophagic response

TNFalpha-YAP/p65-HK2 axis mediates breast cancer cell migration

Supplementary figures and tables. https://www.ijbs.com/v18p3405s1.pdf 

The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.

The authors have declared that no competing interest exists.