1. INTRODUCTION
The genus
Ebolavirus incudes Ebola virus (EBOV), Bundibugyo virus (BDBV), Sudan virus (SUDV), Taï Forest virus (TAFV), Reston virus (RESTV) and the most recently discovered Bombali virus (BOMV) which differ widely in the pathogenicity in humans and non-human primates [
1]. Among these, EBOV, SUDV and BDBV have caused periodic outbreaks of human disease, and EBOV caused an epidemic in West Africa in 2013-2016 [
2]. Severity of human disease caused by ebolaviruses varies greatly with the mean case fatality rates of 44% for EBOV, 34% for BDBV, 30% for SUDV, no documented fatalities for TAFV and no severe disease for RESTV; pathogenicity of BOMV for human is unknown [
2,
3]. Although human infections can differ greatly in routes of infection, dose of the virus, and healthcare conditions, differences in virulence of individual ebolaviruses can be reliably compared in experimental infections of non-human primates, which demonstrate 100% fatality rates for EBOV and only 50-75% for BDBV (reviewed in ref. [
4]). Multiple mechanisms of immune evasion explain the high pathogenicity of ebolaviruses [
5]. Two ebolavirus proteins mediate interferon antagonism: eVP35 inhibits retinoic acid-inducible gene I (RIG-I)-like receptor signaling and type I interferon (IFN-I) induction [6-10], and eVP24 prevents IFN induced Jak-STAT signaling [11-15]. Interestingly, in the closely related marburgviruses, VP40 rather than VP24 inhibits IFN-I signaling (reviewed in ref. [
5]). While wild type (WT) ebolaviruses do not cause a disease in rodents, rodent-adapted filoviruses do [
4]. Mutational profiling of eVP35 and eVP24 have identified amino acids that are critical for virulence phenotypes in mice [
16,
17], guinea pigs [
18,
19] and also non-human primates [
20].
eVP24 interaction with karyopherin alphas (KPNA) prevents translocation of tyrosine phosphorylated STAT1 to the nucleus, resulting in blockade of IFN-induced gene expression [
15]. eVP24 expression also damages the nuclear envelope [
21]. Comparison of the binding affinities of ebolaviruses with different pathogenicity demonstrated that BDBV VP24 (bVP24) binds to KPNA weakly compared to both EBOV and RESTV VP24 (eVP24 or rVP24). Prior studies that defined the binding interface between eVP24 and KPNA also defined three distinct clusters of amino acids, with each displaying varying impact on the interaction. These findings, coupled with sequence comparison of these three ebolaviruses identified the following four putative residues within cluster 1 (defined in ref. [
15]) in eVP24: P83, N135, R140 and V141. Importantly, these residues are located outside the proposed VP24 STAT-1 binding site [
22] and therefore are unlikely to significantly affect binding to STAT1. These amino acids, when mutated to alanine or substituted with corresponding BDBV residues resulted in a reduced inhibition of IFN-I signaling [
23]. In the meantime, multiple alanine substitutions were necessary to obtain near-complete loss of binding [
15] suggesting that they do not have a significant effect on VP24 expression, folding and stability. Collectively these findings further highlight the significance of this interface to modulate the host responses. However, the significance of KPNA binding residues has not yet been defined in the context of infectious viruses. Therefore, the rationale of this study was to investigate the role of individual amino acids in the VP24 KPNA interface, including those which differ between EBOV and BDBV, in virus replication and inhibition of IFN-I response in the context of infectious virus. Using a panel of recombinant EBOVs with mutations in eVP24, we further define the impact of KPNA binding residues on EBOV replication. Our studies provide further mechanistic insights into how individual amino acids in eVP24 KPNA interface affect eVP24 functions, including those dependent and independent on IFN-I. Our results further support the critical role of amino acids involved in binding to KPNA as a factor affecting ebolavirus replication.
2. MATERIALS AND METHODS
2.1. Cells
The Vero-E6 (African green monkey kidney epithelium) and 769-P (human kidney adenocarcinoma) cell lines were obtained from ATCC. The STAT1-deficient U3A cell line which is unresponsive to both IFN-α and IFN-γ [
24] was obtained from George R. Stark (Lerner Research Institute). Huh7 cells were a generous gift from the Gordan lab at the University of California in San Francisco. Vero-E6 cells were maintained in Modified Eagle Medium (MEM), supplemented with 10% fetal bovine serum (FBS), 1% sodium pyruvate, 1% MEM non-essential amino acids (Sigma) and 0.1% gentamicin sulfate, and plated on tissue culture treated 24-well plates. 769-P cells were maintained in RPMI1640 medium, whereas U3A cells were maintained in DMEM and supplemented with 10% FBS and antibiotics. Huh7 cells were maintained in Dulbecco’s modified Eagle’s medium supplemented with 10% fetal bovine serum. All cell culture reagents were from ThermoFisher Scientific unless indicated otherwise.
2.2. Generation of mutated viruses
To generate the recombinant viruses with mutations in eVP24, we used the EBOV reverse genetics system which included the full length clone, representing the plasmid encoding the genomic RNA of wt EBOV under the control of T7 polymerase, previously modified by addition of the transcriptional cassette encoding enhanced green fluorescent protein (eGFP) between the NP and VP35 genes [
25], generously provided by Dr. J. Towner and Dr. S. Nichol (CDC). The system also included the five support plasmids encoding the EBOV proteins involved in replication and transcription (NP, VP35, L, VP30), and also T7 polymerase [
26] kindly provided by Dr. Y. Kawaoka (University of Wisconsin) and Dr. H. Feldmann (NIH). The full-length clone was used to introduce the mutations shown in Figure 1. DNA fragments containing specific amino acid substitutions in eVP24 gene were synthesized by Integrated DNA Technologies (IDT; Coralville, IA). The synthesized DNA fragments were digested with NcoI and SacI restriction endonucleases and cloned into a pUC19 subclone containing the SalI-SacI fragment of the pEBOV plasmid that included the VP24 gene. Subsequently, the pUC19 subclone containing the mutations were digested with SalI and SacI restriction endonucleases and cloned into the full-length clone. The recovered recombinant viruses were amplified by two passages in Vero-E6 cells, and the presence of the mutations was confirmed by sequencing the entire viral genomes; no reversions to wt genotype or occurrence of any other mutations were observed in any part of the genome of any of the viruses.
The growth kinetics of the viruses were compared in triplicate Vero-E6 and 769-P monolayers; cells were infected at an MOI of 0.01 PFU/cell, incubated for 2 h, washed 3 times with 1x phosphate buffered saline (PBS), covered in MEM or RPMI1640 medium with 5% FBS for Vero-E6 and 769-P cells, respectively, and incubated for 6 days. To assess the growth kinetics in the presence of IFN-I, cells were treated with 100 U/ml of IFN-β (PBL Assay Science) at four different timepoints, -6, 0, 24 and 48 hours post infection (hpi). eGFP fluorescence at different timepoints in the presence and absence of IFNs was photographed. Daily aliquots were collected, flash frozen and used for quantitation of the viruses by plaque titration. Plaques positive for eGFP were counted under a UV microscope. The viral growth was quantified using the area under the curve (AUC) using the multistep growth curve and calculated by the trapezoidal rule [
27], using exact viral titers at 2-6 days post-infection, as determined by plaque assay in Vero-E6 and 769-P cells.
2.3. EBOV tetracistronic minigenome reporter assay
Huh7 cells were seeded at a density of 3 x 10
4 cells/well in a 96-well format. TransIt-LT1 transfection reagent (Mirus Bio) was used at a 3:1 ratio of reagent to DNA to transfect Huh7 cells with L at 166.7 ng, VP35 at 20.8 ng, NP at 20.8 ng, VP30 at 12.5 ng, T7 at 41.7 ng, and the EBOV tetracistronic minigenome at 41.7 ng per well. The EBOV tetracistronic minigenome encodes for Renilla luciferase, VP40, GP, and VP24 under a T7 promoter for initial transcription of the minigenome [
28]. As a transfection control, 0.8 ng per well of firefly luciferase was included. In addition to the plasmids described above, VP24 WT or mutant plasmids were also transfected at 5, 10 or 20.8 ng per well, with empty pCAGGS vector as a control. At 72 h post transfection, luciferase activity was measured using a dual-luciferase assay (Promega). Luciferase assay results were read using a BioTek Cytation 10 Confocal Imaging Reader.
2.4. ISG54 luciferase assay
Vero-E6 cells were transfected with an ISG54 firefly luciferase reporter plasmid and a constitutively active Renilla luciferase reporter plasmid (pRL-tk) (Promega, Madison, WI). At 24 h post-transfection, the cells were infected with WT or the mutant viruses at an MOI of 2 PFU/cell. The infected cells were treated or not with 1,000 U/ml of human IFN-β at either 24 or 48 hpi. At 24 h post-IFN-β treatment, cells were lysed, and a dual luciferase reporter assay (Promega) was performed. Firefly luciferase luminescence values were normalized to Renilla luciferase luminescence values.
2.5. STAT1 nuclear translocation assays
The assay was performed as described previously [
11]. Briefly, Vero-E6 cells were seeded onto 8-well chamber slides and were infected with either WT or the mutant viruses expressing eGFP at an MOI of 2 PFU/cell. The infected cells were treated with 1,000 U/ml human IFN-β for 30 min at 37°C. Cells were washed twice with PBS containing 0.9 mM CaCl
2 and 0.5 mM MgCl
2 (PBS-CM), fixed with 4% paraformaldehyde for 30 min, permeabilized for 10 min with PBS containing 0.3% Triton X-100 and blocked for 45 min at room temperature with 4% normal goat serum in PBS containing 0.5% BSA, 0.15% glycine (PBG) and 0.3% Triton X-100. The fixed cells were incubated with rabbit anti-STAT1 (1:400) (5 mg/ml; Santa Cruz) and Alexa Fluor 594 conjugated goat anti-rabbit IgG diluted 1:1,000 in PBS containing 1% BSA and 0.3% Triton X-100 for 1 h at room temperature. Slides were washed 3 times with PBS-CM, fixed with 10% formalin for 48 h and removed from the BSL-4 containment. Slides were washed with PBS 3 times, incubated with 6-diamin-2-phenylindole-dihydrochloride (DAPI) (ThermoFisher Scientific) at the final concentration of 1 µg/ml for 3 min and washed 3 times in PBS. Slides were analyzed by laser scanning confocal microscopy using an Olympus FV1000 confocal microscope. For DAPI, a 405 nm laser and 425/50 filter were used, for eGFP, a 488 nm laser and a 515/30 filter were used, and for AlexaFluor 594, a 543 nm laser and a 610/50 filter were used. All images were acquired using a 60x oil objective.
2.6. Viral replication kinetics in U3A cells
U3A cells were plated on 12-well plates at 1 × 105 cells per well and were infected with WT or mutant viruses at an MOI of 2 PFU/cell, adsorbed for one h, washed three times with PBS, covered in 1 ml of DMEM containing 5% FBS, and incubated for 2 days. eGFP fluorescence was photographed 48 h after infection. Supernatants were harvested, and the viruses were quantified by plaque titration in Vero-E6 cells as described above.
2.7. RNA isolation
Both supernatants and cell lysates of infected cells were inactivated by adding Trizol (ThermoFisher Scientific) with subsequent incubation for 10 min at room temperature, and removed from the BSL-4 facility. RNA isolation was performed using a Direct-Zol RNA miniprep kit (Zymo Research) according to the manufacturer’s recommendations, with on-column DNase treatment for cell lysates. The concentrations of the extracted RNA were measured on NanoDrop 2000 (ThermoFisher Scientific), and 200 ng was used in qPCR assays.
2.8. Quantification of gRNA in supernatants
Viral titers in supernatants of infected Vero-E6 cells were measured by plaque assay as described above. EBOV gRNA in supernatants was measured by quantitative RT-PCR as previously described [
29] with some modifications. First strand cDNA synthesis was performed with a sense-strand, nucleoprotein gene-specific primer and SuperScript IV reverse transcriptase (ThermoFisher Scientific) as per manufacturer’s protocols. qRT-PCR was performed with QuantStudio 6 Real-Time PCR system by combining 4 ng of cDNA, TaqManFast Advanced Master Mix and custom primer-probe (FAM) based qPCR assays (IDT). Absolute quantification was achieved using a standard curve generated by serial dilution of a virus stock template with known titer (PFU/ml).
2.9. Quantification of gRNA and mRNA in cell lysates
Strand-specific qRT-PCR assays to distinguish EBOV gRNA and VP35 mRNA were performed as previously described [
30] using reverse transcription primers that contain unique 18-20 nucleotide tag besides the strand specific sequence. Absolute quantification was performed using amplicons purchased from IDT.
2.10. Quantification of IFNλ1 expression in 769-P cells
cDNA synthesis was performed using oligo(dT)20 primers and SuperScript IV reverse transcriptase (ThermoFisher Scientific) as per manufacturer’s instructions. The following qPCR primers were used to amplify IFNλ1: forward, GGAGCTAGCGAGCTTCAAGA; reverse, ACTCCAGTTTTTCAGCTTGAGTG; FAM probe, GCCAGGGACGCCTTGGAAGAG. Predesigned PrimeTime qPCR assays were purchased for ACTB control (IDT Assay ID: Hs.PT.39a.22214847). Multiplex qPCR assay was performed by combining primer-probes (IDT) for both the target gene (IFNλ1-FAM probe) and housekeeping gene (ACTB-Cy5 probe) in the same tube, and the relative expression level of the IFNλ1 gene was analyzed using the 2−ΔΔCT method [
31]. gRNA quantitation in the 769-P cell lysates was also performed as described above and was used for further normalization of IFNλ1 gene expression for each virus.
2.11. Quantification of IFNβ and ISG15 expression in 769-P cells
0.2 µg of the RNA samples were reverse transcribed using the SuperScript™ IV First-Strand Synthesis System (Invitrogen) based on the manufacturer’s protocol with appropriate primer pairs (2 µM). The oligonucleotide primer sequences used in the experiment are as follows - IFNβ F: GTCAGAGTGGAAATCCTAAG, IFNβ R: ACAGCATCTGCTGGTTGAAG, ISG15 F: TCCTGGTGAGGAATAACAAGGG, ISG15 R: GTCAGCCAGAACAGGTCGTC, β-actin F: ACTGGAACGGTGAAGGTGAC and β-actin R: GTGGACTTGGGAGAGGACTG. Real-time qPCR was conducted using PerfeCTa® SYBR® Green FastMix® (QuantaBio). The PCR reaction was performed on a Bio-Rad CFX Opus 96 qPCR instrument under the following conditions: 95°C for 10 min, 40 cycles of 95°C for 15 sec, 60°C for 1 min. IFNβ and ISG15 were normalized using β-actin as the internal control and the fold difference in gene expression was calculated using the threshold cycle (ΔΔCT) method.
2.12. Analysis of VP24 by Western blotting
Vero-E6 cells were infected with EBOV WT or mutant viruses at an MOI of 2 PFU/cell for 3 days at 37°C, lysed in 4x SDS Laemmli Buffer, diluted to 1x in RIPA buffer, boiled for 15 min and loaded on 4-12% gradient gels (all materials from ThermoFisher Scientific). Separated proteins were blotted onto nitrocellulose membranes using the I-Blot2 system (ThermoFisher Scientific), blocked for 1 h with Odyssey PBS blocking buffer and incubated with primary antibodies overnight at 4°C. The following primary antibodies from ThermoFisher Scientific were used: rabbit polyclonal EBOV VP24 and mouse monoclonal actin antibody. Blots were washed 3 times for 5 min each with PBS buffer with 0.05% Tween-20 (PBST) on a shaker at room temperature. Anti-rabbit and anti-mouse fluorescent-conjugated secondary antibodies were purchased from LI-COR. Secondary antibodies were used at 1:15,000 dilution in PBST. Following 1 h-long incubation, blots were washed 3 times for 5 min each with PBST. Fluorescent bands were detected using an Odyssey Fc imaging system (LI-COR).
2.13. Statistical Methods
Statistical analyses were performed using GraphPad Prism version 6 (GraphPad Software). The following statistical tests were used: one-way ANOVA with multiple comparisons – Dunnett’s multiple comparison test (for the data shown in Figures 1, 5–8), two-way ANOVA with Tukey’s multiple comparison test (for the data shown in Figure 4), and two-way ANOVA with Dunnett’s multiple comparison test (for the data shown in
Table 1).
3. RESULTS
3.1. eVP24 KPNA interface amino acids are required for effective replication of the virus
eVP24 binds to the NPI-1 subfamily of KPNAs: KPNA1, KPNA5, and KPNA6 [
13,
14]. Interestingly, a previous study demonstrated that bVP24 binds to KPNA1, KPNA5, and KPNA6 with lower affinities compared to eVP24, resulting in a reduced inhibition of IFN-I signaling [12-15,22,23]. The amino acids on the eVP24-KPNA5 binding surface associated with this reduced binding affinity were identified: P83, N135, R140 and V141 [
23]. To examine the role of these eVP24 amino acids in a context of EBOV infection, we designed a panel of five mutant viruses, each carrying either one or more amino acid mutations in VP24 and each expressing eGFP (
Figure 1A,B). Specifically, we introduced point mutations N135A or R140A or both. We also designed two more recombinant viruses with substitutions in all four identified amino acids to that in bVP24 termed as eVP24 4x (eVP24 P83S/N135Q/R140H/V141A) or to alanine/glycine termed as eVP24 4xA (eVP24 P83A/N135A/R140A/V141G). All five mutated viruses were successfully recovered and demonstrated no significant difference in eVP24 expression (
Figure 1C). Over-expression of eVP24 can inhibit EBOV transcription [
32]. To ensure that the mutations do not disrupt the overall folding and biological properties of the protein, we over-expressed the eVP24 constructs in a context of an EBOV tetracistronic minigenome assay that was developed based on a published system [
28,
33]. All of the VP24 mutants caused a dose-dependent reduction of the
Renilla luciferase signal similar to WT VP24 (
Figure 1D). Retention of this function suggests that the mutations do not disrupt the overall structure of VP24.
As eVP24 antagonizes IFN-I response, we next compared the multistep growth kinetics of these mutants both in Vero-E6 cells, which are deficient in IFN-I [
34] and 769-P cells, which are IFN-I competent, either in the presence or absence of exogenously added IFN-β (
Figure 2, 3,
Table 1). Overall, the mutated viruses demonstrated decreased growth kinetics as compared to EBOV WT at 2-6 days post infection (dpi) in both cell lines with or without IFN-I, with the R140A mutant being the most attenuated (
Table 1). Specifically, in Vero-E6 cells, a significant decrease in peak viral titer relative to WT virus was observed for all mutant viruses in the presence of IFN-I. Unexpectedly, the R140A and N135/R140A mutants appeared to be also attenuated in the absence of exogenously added IFN-I (
Figure 2 and
Table 1). In 769-P cells, a significant decrease in peak viral titer was observed for the R140A, N135/R140A and 4xA mutants in presence of added IFN-I, and the R140A mutant showed a significant decrease in the absence of added IFN-I (
Figure 3 and
Table 1). These results demonstrate that all the identified amino acids play roles in viral replication, with the most prominent role of R140 whose substitution resulted in a reduced viral growth in both cell lines even in the absence of added IFN-I.
3.3. eVP24 R140 is required for effective EBOV replication independent of IFN-I antagonism
To determine the impact of STAT signaling on EBOV replication, we determined the growth kinetics of these mutant viruses in U3A cells, a human cell line deficient in STAT1 [
24]. Many viruses replicate well in STAT1 deficient cell lines due to deficient IFN responsiveness [
35,
36]. Although EBOV was able to replicate in U3A cells, the titers were lower than in Vero-E6 cells. We did not observe decreased viral titers for four of the mutant viruses at 2 dpi relative to WT (
Figure 5A,B). The R140A mutation, however, led to a significant decrease in the viral titer compared to WT virus (
Figure 5A,B). The reduced growth of the R140A mutant in STAT1 deficient cells suggests that the R140A mutation at least in part acts via an IFN-I-independent pathway, while the growth phenotypes of the remaining mutant viruses appear to be IFN-I-dependent. These data suggest that R140, despite being located in the KPNA interface, is required for effective viral replication at least in part not due to IFN-I antagonism.
3.4. Mutations in eVP24 KPNA interface reduce the amounts of both viral genome and mRNA
eVP24 has been known to regulate both replication and transcription of EBOV genome [
32,
37], but at present the role of KPNA interface in replication and transcription have yet to be described. We used the panel of eVP24 mutated viruses to compare the amounts of their genome and mRNA to WT virus. We infected Vero-E6 cells with EBOV WT or the mutant viruses at an MOI of 2 PFU/cell, harvested cells at 8 and 24 hpi and performed strand-specific qRT-PCR for both viral genomic RNA (gRNA) and mRNA (
Figure 6). Two single mutants, N135A and R140A, demonstrated decrease in gRNA levels at 8 hpi (
Figure 6A), while all five mutants demonstrated significantly reduced gRNA relative to WT at 24 hpi (
Figure 6C).
Further, the R140A mutant and the double mutant N135/R140A demonstrated a significantly reduced amounts of mRNA at 8 hpi (
Figure 6B), and three mutants, N135A, R140A and the double mutant N135/R140A demonstrated a significantly reduced amounts of mRNA at 24 hpi (
Figure 6D). The reduction in both gRNA/mRNA for the mutants N135A, R140A support their role in genome replication and/or transcription. It has been shown that eVP24 plays a role in regulation of viral replication and transcription [
32,
37,
38], and therefore the residues N135 and R140 are likely to be involved in these processes. Thus, the reduced growth of the R140A virus in STAT1 KO cells (
Figure 5B) can be attributed to the involvement of R140 in a yet poorly defined process that involves EBOV genome replication and/or transcription.
Next, we compared the amounts of viral particles released in the media by measuring viral titers by plaque assay, as well as quantifying gRNA by qRT-PCR. At 8 hpi, some of the mutants demonstrated reduced titers which appeared not to be significant due to the high sample-to-sample variation (
Figure 7A,B). However, at 24 hpi, a significant reduction in viral titers (57%, 63%, 44% and 73% for N135A, R140A, N135/R140A and 4xA mutants, respectively) was observed in both assays for all mutants except the 4x mutant relative to WT (
Figure 7C,D). For the 4x mutant, which has all four mutated residues analogous to BDBV, a modest but statistically significant reduction of 48% was observed only by the plaque assay. These results demonstrate the involvement of all four amino acids in EBOV replication.
3.5. Unlike wt eVP24, eVP24 with bVP24 interface residues does not inhibit IFNλ1 or IFNβ expression
EBOV infection is characterized by a very effective inhibition of IFN-I (IFNα, IFNβ) production and signaling, and both eVP35 and eVP24 viral proteins play a critical role in this inhibition [
39]. Recently eVP24 has also been shown to inhibit IFNλ1 gene expression, and this inhibition was lost when KPNA binding residues were mutated [
40]. To determine if any of the mutations in this study affect induction of IFNλ response, we measured IFNλ1 in 769-P cells infected with EBOV WT or the mutant viruses by qRT-PCR. Strikingly, we observed a 5.5- and 6.7-fold increase in IFNλ1 level at 2 and 3 dpi, respectively, for the 4x mutant which has all 4 residues identical to that of BDBV, and no increase for any of the other four mutants (
Figure 8A,B). Our data suggest that these 4 residues on eVP24 are critical for its IFNλ antagonistic activity.
To determine whether we see a similar expression pattern for type I IFN in these mutant viruses, we quantified gene expression of both IFN-β, as well as ISG15 which is induced by type I IFN. For both targets we did not observe any differences for any of the mutants relative to WT at 2 dpi (
Figure 8C,E). However, at 3 dpi we observed a significant decrease in both IFN-β and ISG15 levels for 2 mutants – N135/R140A and 4xA (
Figure 8D,F). Additionally, the two single mutants – N135A and R40A showed a decrease in ISG15 expression only. Conversely, 4x mutant showed a significant increase in both I IFN-β and ISG15 expression. Thus, the 4x mutant causes an increase in both type-I and type III IFN response. Taken together, these data suggest a novel role for the KPNA interacting cluster 1 residues and a possible mechanism of the lower pathogenicity of BDBV as compared to EBOV.
4. DISCUSSION
We performed this study to determine if the KPNA binding amino acids in bVP24 can alter EBOV infection. Incorporation of the selected mutations in the viral genome and successful recovery of infectious mutated viruses suggest that the mutations did not significantly affect VP24-specific functions. We observed a significant attenuation of all five mutated viruses in the presence of exogenous IFN-I in Vero-E6 cells, but only three of the four alanine mutants were attenuated in IFN-competent 769-P cells. Similarly, all five mutant viruses inhibited IFN-induced ISG54 gene expression less efficiently than WT EBOV due to increased STAT1 nuclear translocation. Additionally, only the EBOV mutant containing all four amino acids identical to that in BDBV caused a significant increase in IFNλ1, IFNβ and ISG15 gene expression. Furthermore, the mutant R140A was significantly attenuated in STAT1-KO cells and demonstrated a significant reduction in the amounts of viral gRNA and mRNA. Thus, the VP24 residues binding KPNAs are likely to be responsible for the reduced replication of BDBV compared to EBOV and may influence ebolavirus pathogenicity both by IFN-I-dependent and IFN-I-independent mechanisms.
We observed a reduced growth kinetics for all eVP24 mutant viruses in Vero-E6 cells mainly in the presence of IFN-I. Most importantly, one of the mutants, R140A, was severely attenuated in both Vero-E6 and 769-P cells irrespective of the presence or absence of exogenous IFN-I, suggesting a critical role of this amino acid in viral pathogenesis in an IFN-I-independent manner. The paradoxical lack of attenuation in the presence of additional mutations suggests possible compensatory effects.
It was demonstrated that the eVP24 residues that bind to karyopherin alpha are necessary for IFN-I antagonistic activity in a plasmid-based transfection system [
23]. Moreover, one study provided evidence that eVP24 can also directly bind STAT1 and identified residues on VP24 responsible for this interaction [
41]. In the current study, we demonstrated that all five mutated viruses, generated based on these structural studies and, prior comparison of eVP24 and bVP24, did not block STAT1 translocation to the nucleus resulting in enhanced expression of IFN stimulated genes. To delineate STAT1-independent effects, we also examined the growth of these viruses in U3A cells, a STAT1-deficient cell line. The single mutants N135A and R140A showed opposite effects, while the double mutant demonstrated no change compared to the WT virus. Mutations in the IFN-I antagonist C protein of Sendai virus led to reduced growth kinetics in IFN-I-responsive 2fTGH cells but similar growth in STAT1-deficient U3A cells suggesting that the IFN-inhibiting capacity of the C protein was critical for effective replication [
35]. In this study, only the R140A mutant, which was attenuated in STAT1-deficient U3A cells, was also attenuated in IFN-I-responsive 769-P cells, suggesting a STAT1- and IFN-I-independent mechanism of attenuation.
It is interesting that the ebolavirus VP24 residue 140, which is located in the cluster 1 of the KPNA binding residues, is not conserved between EBOV, BDBV and RESTV [
41]. The location of the residue adjacent to the KPNA1 binding domain of VP24 that resides in the amino acid 142-146 loop [
12] led to speculation that it is responsible for the different pathogenicity of the three ebolaviruses [
42,
43]. Additionally, the VP24 residues 137 and 140 are also critical for KPNA5 binding and had the greatest impact on IFN-I-sensitive response element activity [
15]. Multiple residues in this region have been identified as specificity determining positions (SDPs) which are positions that are conserved within protein subfamilies but differ between them. It was suggested that SDPs distinguish between the different functional specificities of proteins from the different
Ebolavirus species based on differences between the pathogenic and non-pathogenic ebolaviruses by molecular dynamics and structural analysis [
42,
44,
45]. However, this is the first study besides the one with the K142A mutant [
46] wherein multiple eVP24 residues have been mutated and characterized in context of recombinant viruses.
Unlike eVP35 that plays a major role in both IFN antagonism and viral transcription, the role of the eVP24 in viral transcription is less clear. Initial studies with plasmid-based monocistronic minigenome systems showed that eVP24 inhibits both genome replication and transcription [
32,
37]. But later studies using a tetracistronic minigenome, as well as infectious EBOV, demonstrated that the inhibiting effect of eVP24 on genome replication and transcription is modest [
33]. Recently it has been demonstrated that residues in the YxxL motif (amino acids 172-175) are critical for this inhibitory effect of eVP24 on transcription and replication [
47]. The current study addressed the effects of eVP24, with single and combined mutations N135A and R140A in a region outside the YxxL region, caused a significant reduction in the amounts of gRNA and mRNA relative to WT. As the R140A mutant demonstrated a reduced growth in STAT1-deficient U3A cells, the attenuation is likely to be attributed to the role of R140 in viral transcription and/or replication. Although none of the mutations significantly affected the VP24 expression and had no impact on the capacity of eVP24 to inhibit transcription of a tetracistronic minigenome, we cannot completely rule out that these mutations could have subtle effects on protein folding.
Recent studies have demonstrated that binding of eVP24 to KPNA is also critical for inhibiting of IFNλ1 gene expression [
40]. They confirmed that the mutant L201A/E203A/P204A/D205A/S207A, which did not bind KPNA [
15], was unable to inhibit IFNλ1 gene expression. In our study, we evaluated other KPNA binding residues and found that only the 4x mutant (P83S/N135Q/R140H/V141A) which had BDBV residues at all four positions induced significantly increased expression of IFN-λ1 as compared to WT virus. Further, only this 4x mutant caused a significant induction of both IFN-β as well as ISG15. These results support prior studies wherein mutation of KPNA binding sites in eVP24 reduces nuclear localization thereby reducing its ability to efficiently inhibit IFN gene expression [
15,
48,
49]. Thus, the reduced inhibition of both IFN-I and IFN-III activity may contribute to the reduced pathogenicity of BDBV relative to EBOV.
The limitations of the study include that, at least theoretically, the mutations in eVP24-karyopherin alpha interface may affect viral entry and uncoating. However, given the lack of any significant effect on VP24 conformation (
Figure 1B, C) and transcription inhibition function (
Figure 1D) we believe this scenario is unlikely.
Taken together, our studies for the first time evaluate the effect of a key KPNA binding interface of VP24 with infectious recombinant viruses. More importantly, our studies identify IFN-I-dependent and IFN-I-independent effects of eVP24, including virus replication. Future studies are required to fully define the mechanisms by which the KPNA binding interface impacts virus replication, to investigate if these mutated viruses are attenuated in animal models and to determine if they can induce a protective immune response. Collectively, these studies further highlight the significant role played by the KPNA binding interface in a key virulence factor and a critical therapeutic target, eVP24.
Author Contributions
Conceptualization, C.F.B., G.K.A. and A.B.; methodology: P.R., B.T., R.I.S., P.A.I., N.K., O.A.V., N.T., C.F.B., G.K.A. and A.B.; validation: P.R., B.T., R.I.S., P.A.I., N.K., O.A.V. and N.T.; formal analysis: P.R., B.T., R.I.S., P.A.I., N.K., O.A.V. and N.T.; investigation: P.R., B.T., R.I.S., P.A.I., N.K., O.A.V., N.T., H.A., C.W., C.F.B., G.K.A. and A.B.; resources: A.B.; data curation: P.R., B.T., R.I.S.; writing – original draft preparation: P.R.; writing – review and editing: P.R., B.T., P.A.I., C.F.B., G.K.A. and A.B.; visualization: P.R., B.T. and R.I.S.; supervision: A.B.; project administration: A.B.; funding acquisition: C.F.B., G.K.A. and A.B. All authors have read and agreed to the published version of the manuscript.
Funding
The authors were supported in part by U.S. National Institutes of Health grants U19 AI109664 to C.F.B., G.K.A. and A.B, P01AI120943 and R01AI143292 to G.K.A and C.F.B, R01AI123926, and R01AI161374 to G.K.A. and R01AI148663 to C.F.B. O.A.V. was supported in part by Public Health Service institutional research training award AI07647.
Conflicts of Interests
The authors declare no conflict of interests.
References
- Goldstein, T.; Anthony, S.J.; Gbakima, A.; Bird, B.H.; Bangura, J.; Tremeau-Bravard, A.; Belaganahalli, M.N.; Wells, H.L.; Dhanota, J.K.; Liang, E.; et al. The discovery of Bombali virus adds further support for bats as hosts of ebolaviruses. Nature microbiology 2018, 3, 1084–1089. [Google Scholar] [CrossRef] [PubMed]
- Jacob, S.T.; Crozier, I.; Fischer, W.A., 2nd; Hewlett, A.; Kraft, C.S.; Vega, M.A.; Soka, M.J.; Wahl, V.; Griffiths, A.; Bollinger, L.; et al. Ebola virus disease. Nat Rev Dis Primers 2020, 6, 13. [Google Scholar] [CrossRef]
- Martell, H.J.; Masterson, S.G.; McGreig, J.E.; Michaelis, M.; Wass, M.N. Is the Bombali virus pathogenic in humans? Bioinformatics (Oxford, England) 2019, 35, 3553–3558. [Google Scholar] [CrossRef] [PubMed]
- Yamaoka, S.; Ebihara, H. Pathogenicity and Virulence of Ebolaviruses with Species- and Variant-specificity. Virulence 2021, 12, 885–901. [Google Scholar] [CrossRef] [PubMed]
- Messaoudi, I.; Amarasinghe, G.K.; Basler, C.F. Filovirus pathogenesis and immune evasion: insights from Ebola virus and Marburg virus. Nat Rev Microbiol 2015, 13, 663–676. [Google Scholar] [CrossRef]
- Basler, C.F.; Mikulasova, A.; Martinez-Sobrido, L.; Paragas, J.; Muhlberger, E.; Bray, M.; Klenk, H.D.; Palese, P.; Garcia-Sastre, A. The Ebola virus VP35 protein inhibits activation of interferon regulatory factor 3. J Virol 2003, 77, 7945–7956. [Google Scholar] [CrossRef]
- Cardenas, W.B.; Loo, Y.M.; Gale, M., Jr.; Hartman, A.L.; Kimberlin, C.R.; Martinez-Sobrido, L.; Saphire, E.O.; Basler, C.F. Ebola Virus VP35 Protein Binds Double-Stranded RNA and Inhibits Alpha/Beta Interferon Production Induced by RIG-I Signaling. J Virol 2006, 80, 5168–5178. [Google Scholar] [CrossRef]
- Prins, K.C.; Cardenas, W.B.; Basler, C.F. Ebola virus protein VP35 impairs the function of interferon regulatory factor-activating kinases IKKepsilon and TBK-1. Journal of virology 2009, 83, 3069–3077. [Google Scholar] [CrossRef]
- Leung, D.W.; Ginder, N.D.; Fulton, D.B.; Nix, J.; Basler, C.F.; Honzatko, R.B.; Amarasinghe, G.K. Structure of the Ebola VP35 interferon inhibitory domain. Proceedings of the National Academy of Sciences of the United States of America 2009, 106, 411–416. [Google Scholar] [CrossRef]
- Leung, D.W.; Prins, K.C.; Borek, D.M.; Farahbakhsh, M.; Tufariello, J.M.; Ramanan, P.; Nix, J.C.; Helgeson, L.A.; Otwinowski, Z.; Honzatko, R.B.; et al. Structural basis for dsRNA recognition and interferon antagonism by Ebola VP35. Nat Struct Mol Biol 2010, 17, 165–172. [Google Scholar] [CrossRef]
- Reid, S.P.; Leung, L.W.; Hartman, A.L.; Martinez, O.; Shaw, M.L.; Carbonnelle, C.; Volchkov, V.E.; Nichol, S.T.; Basler, C.F. Ebola virus VP24 binds karyopherin alpha1 and blocks STAT1 nuclear accumulation. Journal of virology 2006, 80, 5156–5167. [Google Scholar] [CrossRef] [PubMed]
- Mateo, M.; Reid, S.P.; Leung, L.W.; Basler, C.F.; Volchkov, V.E. Ebolavirus VP24 Binding to Karyopherins Is Required for Inhibition of Interferon Signaling. J Virol 2010, 84, 1169–1175. [Google Scholar] [CrossRef] [PubMed]
- Reid, S.P.; Leung, L.W.; Hartman, A.L.; Martinez, O.; Shaw, M.L.; Carbonnelle, C.; Volchkov, V.E.; Nichol, S.T.; Basler, C.F. Ebola Virus VP24 Binds Karyopherin {alpha}1 and Blocks STAT1 Nuclear Accumulation. J Virol 2006, 80, 5156–5167. [Google Scholar] [CrossRef]
- Reid, S.P.; Valmas, C.; Martinez, O.; Sanchez, F.M.; Basler, C.F. Ebola virus VP24 proteins inhibit the interaction of NPI-1 subfamily karyopherin alpha proteins with activated STAT1. Journal of virology 2007, 81, 13469–13477. [Google Scholar] [CrossRef]
- Xu, W.; Edwards, M.R.; Borek, D.M.; Feagins, A.R.; Mittal, A.; Alinger, J.B.; Berry, K.N.; Yen, B.; Hamilton, J.; Brett, T.J.; et al. Ebola virus VP24 targets a unique NLS binding site on karyopherin alpha 5 to selectively compete with nuclear import of phosphorylated STAT1. Cell Host Microbe 2014, 16, 187–200. [Google Scholar] [CrossRef] [PubMed]
- Ebihara, H.; Takada, A.; Kobasa, D.; Jones, S.; Neumann, G.; Theriault, S.; Bray, M.; Feldmann, H.; Kawaoka, Y. Molecular determinants of Ebola virus virulence in mice. PLoS Pathog 2006, 2, e73. [Google Scholar] [CrossRef]
- Hartman, A.L.; Bird, B.H.; Towner, J.S.; Antoniadou, Z.A.; Zaki, S.R.; Nichol, S.T. Inhibition of IRF-3 activation by VP35 is critical for the high level of virulence of ebola virus. J Virol 2008, 82, 2699–2704. [Google Scholar] [CrossRef]
- Prins, K.C.; Delpeut, S.; Leung, D.W.; Reynard, O.; Volchkova, V.A.; Reid, S.P.; Ramanan, P.; Cardenas, W.B.; Amarasinghe, G.K.; Volchkov, V.E.; et al. Mutations abrogating VP35 interaction with double-stranded RNA render Ebola virus avirulent in guinea pigs. Journal of virology 2010, 84, 3004–3015. [Google Scholar] [CrossRef]
- Mateo, M.; Carbonnelle, C.; Reynard, O.; Kolesnikova, L.; Nemirov, K.; Page, A.; Volchkova, V.A.; Volchkov, V.E. VP24 is a molecular determinant of Ebola virus virulence in guinea pigs. J Infect Dis 2011, 204 Suppl 3, S1011–1020. [Google Scholar] [CrossRef]
- Woolsey, C.; Menicucci, A.R.; Cross, R.W.; Luthra, P.; Agans, K.N.; Borisevich, V.; Geisbert, J.B.; Mire, C.E.; Fenton, K.A.; Jankeel, A.; et al. A VP35 Mutant Ebola Virus Lacks Virulence but Can Elicit Protective Immunity to Wild-Type Virus Challenge. Cell Rep 2019, 28, 3032–3046. [Google Scholar] [CrossRef]
- Vidal, S.; Sánchez-Aparicio, M.; Seoane, R.; El Motiam, A.; Nelson, E.V.; Bouzaher, Y.H.; Baz-Martínez, M.; García-Dorival, I.; Gonzalo, S.; Vázquez, E.; et al. Expression of the Ebola Virus VP24 Protein Compromises the Integrity of the Nuclear Envelope and Induces a Laminopathy-Like Cellular Phenotype. mBio 2021, 12, e0097221. [Google Scholar] [CrossRef] [PubMed]
- Zhang, A.P.; Bornholdt, Z.A.; Liu, T.; Abelson, D.M.; Lee, D.E.; Li, S.; Woods, V.L., Jr.; Saphire, E.O. The Ebola Virus Interferon Antagonist VP24 Directly Binds STAT1 and Has a Novel, Pyramidal Fold. PLoS pathogens 2012, 8, e1002550. [Google Scholar] [CrossRef] [PubMed]
- Schwarz, T.M.; Edwards, M.R.; Diederichs, A.; Alinger, J.B.; Leung, D.W.; Amarasinghe, G.K.; Basler, C.F. VP24-Karyopherin Alpha Binding Affinities Differ between Ebolavirus Species, Influencing Interferon Inhibition and VP24 Stability. Journal of virology 2017, 91. [Google Scholar] [CrossRef] [PubMed]
- McKendry, R.; John, J.; Flavell, D.; Muller, M.; Kerr, I.M.; Stark, G.R. High-frequency mutagenesis of human cells and characterization of a mutant unresponsive to both alpha and gamma interferons. Proc Natl Acad Sci U S A 1991, 88, 11455–11459. [Google Scholar] [CrossRef] [PubMed]
- Towner, J.S.; Paragas, J.; Dover, J.E.; Gupta, M.; Goldsmith, C.S.; Huggins, J.W.; Nichol, S.T. Generation of eGFP expressing recombinant Zaire ebolavirus for analysis of early pathogenesis events and high-throughput antiviral drug screening. Virology 2005, 332, 20–27. [Google Scholar] [CrossRef] [PubMed]
- Neumann, G.; Feldmann, H.; Watanabe, S.; Lukashevich, I.; Kawaoka, Y. Reverse genetics demonstrates that proteolytic processing of the Ebola virus glycoprotein is not essential for replication in cell culture. J Virol 2002, 76, 406–410. [Google Scholar] [CrossRef] [PubMed]
- Wang, G.P.; Bushman, F.D. A statistical method for comparing viral growth curves. Journal of virological methods 2006, 135, 118–123. [Google Scholar] [CrossRef]
- El-Sayed, N.S.; Jureka, A.S.; Edwards, M.R.; Lohan, S.; Williams, C.G.; Keiser, P.T.; Davey, R.A.; Totonchy, J.; Tiwari, R.K.; Basler, C.F.; et al. Synthesis and antiviral activity of fatty acyl conjugates of remdesivir against severe acute respiratory syndrome coronavirus 2 and Ebola virus. Eur J Med Chem 2021, 226, 113862. [Google Scholar] [CrossRef]
- Meyer, M.; Garron, T.; Lubaki, N.M.; Mire, C.E.; Fenton, K.A.; Klages, C.; Olinger, G.G.; Geisbert, T.W.; Collins, P.L.; Bukreyev, A. Aerosolized Ebola vaccine protects primates and elicits lung-resident T cell responses. J Clin Invest 2015, 125, 3241–3255. [Google Scholar] [CrossRef]
- Prasad, A.N.; Ronk, A.J.; Widen, S.G.; Wood, T.G.; Basler, C.F.; Bukreyev, A. Ebola Virus Produces Discrete Small Noncoding RNAs Independently of the Host MicroRNA Pathway Which Lack RNA Interference Activity in Bat and Human Cells. Journal of virology 2020, 94. [Google Scholar] [CrossRef]
- Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods (San Diego, Calif.) 2001, 25, 402–408. [Google Scholar] [CrossRef]
- Watanabe, S.; Noda, T.; Halfmann, P.; Jasenosky, L.; Kawaoka, Y. Ebola virus (EBOV) VP24 inhibits transcription and replication of the EBOV genome. J Infect Dis 2007, 196 Suppl 2, S284–290. [Google Scholar] [CrossRef]
- Watt, A.; Moukambi, F.; Banadyga, L.; Groseth, A.; Callison, J.; Herwig, A.; Ebihara, H.; Feldmann, H.; Hoenen, T. A novel life cycle modeling system for Ebola virus shows a genome length-dependent role of VP24 in virus infectivity. J Virol 2014, 88, 10511–10524. [Google Scholar] [CrossRef] [PubMed]
- Desmyter, J.; Melnick, J.L.; Rawls, W.E. Defectiveness of interferon production and of rubella virus interference in a line of African green monkey kidney cells (Vero). J Virol 1968, 2, 955–961. [Google Scholar] [CrossRef] [PubMed]
- Kato, A.; Kiyotani, K.; Kubota, T.; Yoshida, T.; Tashiro, M.; Nagai, Y. Importance of the anti-interferon capacity of Sendai virus C protein for pathogenicity in mice. Journal of virology 2007, 81, 3264–3271. [Google Scholar] [CrossRef] [PubMed]
- Cuevas, R.A.; Ghosh, A.; Wallerath, C.; Hornung, V.; Coyne, C.B.; Sarkar, S.N. MOV10 Provides Antiviral Activity against RNA Viruses by Enhancing RIG-I-MAVS-Independent IFN Induction. J Immunol 2016, 196, 3877–3886. [Google Scholar] [CrossRef] [PubMed]
- Hoenen, T.; Jung, S.; Herwig, A.; Groseth, A.; Becker, S. Both matrix proteins of Ebola virus contribute to the regulation of viral genome replication and transcription. Virology 2010, 403, 56–66. [Google Scholar] [CrossRef]
- Mateo, M.; Carbonnelle, C.; Martinez, M.J.; Reynard, O.; Page, A.; Volchkova, V.A.; Volchkov, V.E. Knockdown of Ebola virus VP24 impairs viral nucleocapsid assembly and prevents virus replication. J Infect Dis 2011, 204 Suppl 3, S892–896. [Google Scholar] [CrossRef]
- Olejnik, J.; Hume, A.J.; Leung, D.W.; Amarasinghe, G.K.; Basler, C.F.; Mühlberger, E. Filovirus Strategies to Escape Antiviral Responses. Curr Top Microbiol Immunol 2017, 411, 293–322. [Google Scholar] [CrossRef]
- He, F.; Melén, K.; Maljanen, S.; Lundberg, R.; Jiang, M.; Österlund, P.; Kakkola, L.; Julkunen, I. Ebolavirus protein VP24 interferes with innate immune responses by inhibiting interferon-λ1 gene expression. Virology 2017, 509, 23–34. [Google Scholar] [CrossRef]
- Zhang, A.P.; Bornholdt, Z.A.; Liu, T.; Abelson, D.M.; Lee, D.E.; Li, S.; Woods, V.L., Jr.; Saphire, E.O. The ebola virus interferon antagonist VP24 directly binds STAT1 and has a novel, pyramidal fold. PLoS Pathog 2012, 8, e1002550. [Google Scholar] [CrossRef] [PubMed]
- Pappalardo, M.; Collu, F.; Macpherson, J.; Michaelis, M.; Fraternali, F.; Wass, M.N. Investigating Ebola virus pathogenicity using molecular dynamics. BMC Genomics 2017, 18, 566. [Google Scholar] [CrossRef] [PubMed]
- Chakraborty, S.; Rao, B.; Asgeirsson, B.; Dandekar, A. Correlating the ability of VP24 protein from Ebola and Marburg viruses to bind human karyopherin to their immune suppression mechanism and pathogenicity using computational methods [version 2; peer review: 2 approved with reservations]. F1000Research 2014, 3. [Google Scholar] [CrossRef]
- Pappalardo, M.; Juliá, M.; Howard, M.J.; Rossman, J.S.; Michaelis, M.; Wass, M.N. Conserved differences in protein sequence determine the human pathogenicity of Ebolaviruses. Scientific reports 2016, 6, 23743. [Google Scholar] [CrossRef] [PubMed]
- Ding, J.N.; Zhang, Y.J.; Zhong, H.; Ao, C.C.; Li, J.; Han, J.G. An all-atom molecular dynamics study of the anti-interferon signaling of Ebola virus: interaction mechanisms of EBOV VP24 binding to Karyopherin alpha5. Molecular bioSystems 2017, 13, 1031–1045. [Google Scholar] [CrossRef]
- Lubaki, N.M.; Ilinykh, P.; Pietzsch, C.; Tigabu, B.; Freiberg, A.N.; Koup, R.A.; Bukreyev, A. The lack of maturation of Ebola virus-infected dendritic cells results from the cooperative effect of at least two viral domains. Journal of virology 2013, 87, 7471–7485. [Google Scholar] [CrossRef]
- Takamatsu, Y.; Kolesnikova, L.; Schauflinger, M.; Noda, T.; Becker, S. The Integrity of the YxxL Motif of Ebola Virus VP24 Is Important for the Transport of Nucleocapsid-Like Structures and for the Regulation of Viral RNA Synthesis. Journal of virology 2020, 94, e02170–02119. [Google Scholar] [CrossRef]
- He, F.; Melen, K.; Maljanen, S.; Lundberg, R.; Jiang, M.; Osterlund, P.; Kakkola, L.; Julkunen, I. Ebolavirus protein VP24 interferes with innate immune responses by inhibiting interferon-lambda1 gene expression. Virology 2017, 509, 23–34. [Google Scholar] [CrossRef]
- He, F.B.; Khan, H.; Huttunen, M.; Kolehmainen, P.; Melen, K.; Maljanen, S.; Qu, M.; Jiang, M.; Kakkola, L.; Julkunen, I. Filovirus VP24 Proteins Differentially Regulate RIG-I and MDA5-Dependent Type I and III Interferon Promoter Activation. Front Immunol 2021, 12, 694105. [Google Scholar] [CrossRef]
Figure 1.
The EBOV VP24 mutants. A. The virus names and the mutations. B. Structure of VP24 displaying the mutated amino acids in yellow (PyMOL). C. Western blot analysis of VP24 in lysates of Vero-E6 cells infected with EBOV-eGFP WT or its mutated derivatives. D. EBOV tetracistronic minigenome reporter assay. Huh7 cells were transfected with the minigenome and support plasmids. Also included in the transfection was empty vector (No VP24) or plasmids expressing WT and mutated versions of VP24 at indicated amounts. On day 3, luciferase activity was quantified. Mean values based on triplicate samples ±SE. Statistical difference was assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the groups to no VP24 control: *P < 0.01, **P < 0.001, ** P < 0.0001.
Figure 1.
The EBOV VP24 mutants. A. The virus names and the mutations. B. Structure of VP24 displaying the mutated amino acids in yellow (PyMOL). C. Western blot analysis of VP24 in lysates of Vero-E6 cells infected with EBOV-eGFP WT or its mutated derivatives. D. EBOV tetracistronic minigenome reporter assay. Huh7 cells were transfected with the minigenome and support plasmids. Also included in the transfection was empty vector (No VP24) or plasmids expressing WT and mutated versions of VP24 at indicated amounts. On day 3, luciferase activity was quantified. Mean values based on triplicate samples ±SE. Statistical difference was assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the groups to no VP24 control: *P < 0.01, **P < 0.001, ** P < 0.0001.
Figure 2.
Characterization of the VP24 mutants in Vero-E6 cells. A. UV microscopy of cells infected with WT or mutant viruses and treated or not with IFN-β at 2, 4 and 6 dpi, 20× magnification. B. Viral growth kinetics: titers in supernatants of Vero-E6 cells infected at MOl of 0.01 PFU/cell and treated or not with IFN-ß determined by plaque assay. Mean values based on triplicate samples ±SE.
Figure 2.
Characterization of the VP24 mutants in Vero-E6 cells. A. UV microscopy of cells infected with WT or mutant viruses and treated or not with IFN-β at 2, 4 and 6 dpi, 20× magnification. B. Viral growth kinetics: titers in supernatants of Vero-E6 cells infected at MOl of 0.01 PFU/cell and treated or not with IFN-ß determined by plaque assay. Mean values based on triplicate samples ±SE.
Figure 3.
Characterization of the VP24 mutants in 769-P cells. A. UV microscopy of cells infected with WT or mutant viruses and treated or not with IFN-β at 2, 4 and 6 dpi, 20x magnification. B. Viral growth kinetics: titers in supernatants of Vero-E6 cells infected at MOl of 0.01 PFU/cell and treated or not with IFN-ß determined by plaque assay. Mean values based on triplicate samples ±SE.
Figure 3.
Characterization of the VP24 mutants in 769-P cells. A. UV microscopy of cells infected with WT or mutant viruses and treated or not with IFN-β at 2, 4 and 6 dpi, 20x magnification. B. Viral growth kinetics: titers in supernatants of Vero-E6 cells infected at MOl of 0.01 PFU/cell and treated or not with IFN-ß determined by plaque assay. Mean values based on triplicate samples ±SE.
Figure 4.
eVP24 mutants have reduced IFN-I inhibitory activity. A, B. STAT1 nuclear translocation after 30 min-long IFN-ß treatment. A. In uninfected Vero-E6 cells, STAT1 (red) translocates to the nucleus (DAPI staining, blue). B. In infected Vero-E6 cells (green), STAT1 nuclear staining is reduced for EBOVeGFP-WT (the top row) but increased for the mutant viruses. C, D. ISG-54 reporter assay. Inhibition of IFN-induced ISG54 promoter activity in Vero-E6 cells at 48 h (C) and 72 h (D) after infection with EBOVeGFP at MOl of 2 PFU/cell. Mean values based on triplicate samples ±SE. Statistical difference was assessed by two-way ANOVA-Tukey's multiple comparison test comparing the mutants to WT EBov: * P < 0.05, ** P < 0.01; **** P < 0.0001.
Figure 4.
eVP24 mutants have reduced IFN-I inhibitory activity. A, B. STAT1 nuclear translocation after 30 min-long IFN-ß treatment. A. In uninfected Vero-E6 cells, STAT1 (red) translocates to the nucleus (DAPI staining, blue). B. In infected Vero-E6 cells (green), STAT1 nuclear staining is reduced for EBOVeGFP-WT (the top row) but increased for the mutant viruses. C, D. ISG-54 reporter assay. Inhibition of IFN-induced ISG54 promoter activity in Vero-E6 cells at 48 h (C) and 72 h (D) after infection with EBOVeGFP at MOl of 2 PFU/cell. Mean values based on triplicate samples ±SE. Statistical difference was assessed by two-way ANOVA-Tukey's multiple comparison test comparing the mutants to WT EBov: * P < 0.05, ** P < 0.01; **** P < 0.0001.
Figure 5.
The R140A mutation in eVP24 inhibits viral growth in U3A (STAT1-KO) cells. Cells wereinfected with the indicated viruses at MOl of 2 PFU/cell and analyzed at 48 hpi. A. UV microscopy, 20x magnification. B. Viral titers in the supernatants. Mean values based on triplicate samples ±SE. Statistical difference was assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the mutants ' to WT EBOV: *P < 0.05, ** P < 0.01, **** P < 0.0001.
Figure 5.
The R140A mutation in eVP24 inhibits viral growth in U3A (STAT1-KO) cells. Cells wereinfected with the indicated viruses at MOl of 2 PFU/cell and analyzed at 48 hpi. A. UV microscopy, 20x magnification. B. Viral titers in the supernatants. Mean values based on triplicate samples ±SE. Statistical difference was assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the mutants ' to WT EBOV: *P < 0.05, ** P < 0.01, **** P < 0.0001.
Figure 6.
Effect of VP24 mutations on the amounts of EBOV genomic RNA and mRNA. Vero-E6 cells were infected with the panel of viruses at MOl of 2 PFU/cell, and the cells were harvested at 8 hpi (A, B) or 24 hpi (C, D) and analyzed by qRT-PCR for viral genomic RNA (gRNA) (A, C) or VP35 mRNA (B, D). The values show the mean numbers of copies in a 2 ng reaction based on triplicate samples ±SE. Statistical differences were assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the mutants to WT EBOV: * P < 0.05. ** P < 0.01, ***P < 0.001; **** P < 0.0001.
Figure 6.
Effect of VP24 mutations on the amounts of EBOV genomic RNA and mRNA. Vero-E6 cells were infected with the panel of viruses at MOl of 2 PFU/cell, and the cells were harvested at 8 hpi (A, B) or 24 hpi (C, D) and analyzed by qRT-PCR for viral genomic RNA (gRNA) (A, C) or VP35 mRNA (B, D). The values show the mean numbers of copies in a 2 ng reaction based on triplicate samples ±SE. Statistical differences were assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the mutants to WT EBOV: * P < 0.05. ** P < 0.01, ***P < 0.001; **** P < 0.0001.
Figure 7.
Effect of VP24 mutations on EBOV replication. Vero-E6 cells were infected with the panel of viruses at MOI of 2 PFU/cell, the supernatants were collected at 8 hpi (A, B) or 24 hpi (C, D), and the viral titers were determined by plaque assay (A, C) and qRT-PCR (B, D). Mean values ±SE based on triplicate samples. Statistical differences were assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the mutants to WT EBOV: *P < 0.05; ** P < 0.01.
Figure 7.
Effect of VP24 mutations on EBOV replication. Vero-E6 cells were infected with the panel of viruses at MOI of 2 PFU/cell, the supernatants were collected at 8 hpi (A, B) or 24 hpi (C, D), and the viral titers were determined by plaque assay (A, C) and qRT-PCR (B, D). Mean values ±SE based on triplicate samples. Statistical differences were assessed by one-way ANOVA-Dunnett's multiple comparison test comparing the mutants to WT EBOV: *P < 0.05; ** P < 0.01.
Figure 8.
Effect of mutants on IFNA1, IFNβ, ISG15 expression. Gene expression - IFNA1 (A,B), IFNβ (C,D) and ISG15 (E,F)-was quantified by qRT-PCR in 769-P cell lysates infected with EBOV-eGFP at MOl of 2 PFU/cell after 2 dpi (A,C,E) and 3 dpi (B,D,F). Gene expression was normalized to gRNA and expressed relative to WT (100%, shown with the dotted line). Statistical differences were assessed by one-way ANOVA - Dunnett's multiple comparison test comparing the mutants to WT EBOV: **** p < 0.0001; ***p < 0.001;**p < 0.01 and*p < 0.05.
Figure 8.
Effect of mutants on IFNA1, IFNβ, ISG15 expression. Gene expression - IFNA1 (A,B), IFNβ (C,D) and ISG15 (E,F)-was quantified by qRT-PCR in 769-P cell lysates infected with EBOV-eGFP at MOl of 2 PFU/cell after 2 dpi (A,C,E) and 3 dpi (B,D,F). Gene expression was normalized to gRNA and expressed relative to WT (100%, shown with the dotted line). Statistical differences were assessed by one-way ANOVA - Dunnett's multiple comparison test comparing the mutants to WT EBOV: **** p < 0.0001; ***p < 0.001;**p < 0.01 and*p < 0.05.
Table 1.
Growth of the mutated viruses in Vero-E6 and 769-P cells in the presence and absence of exogenously added IFN-β.
Table 1.
Growth of the mutated viruses in Vero-E6 and 769-P cells in the presence and absence of exogenously added IFN-β.
|
Vero-E6 cells |
769-P cells |
Virus |
No IFN-β added
|
IFN-β added
|
No IFN-β added
|
IFN-β added
|
N135A |
19.98 (6.5ns) |
14.48 (5.1****) |
18.78 (6ns) |
13.12 (4.2ns) |
R140A |
17.42 (6.2*) |
11.32 (4.1****) |
14.18 (4.1****) |
8.9 (3.0****) |
135/140 |
18.83 (6.1*) |
14.85 (5.0****) |
16.72 (5.6ns) |
11.97 (3.4***) |
4XA |
19.8 (6.4ns) |
14.82 (5.1****) |
17.07 (5.7ns) |
10.9 (3.2****) |
4X |
19.85 (6.3ns) |
15.18 (5.1****) |
18.17 (5.6ns) |
12.6 (3.9ns) |
WT |
21.68 (6.5) |
16.52 (5.7) |
19.12 (5.4) |
14.28 (4.5) |
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).