Derivation of brain cells from equine induced pluripotent stem cells (eiPSCs) and their permissivity to WNV.
Equine iPSCs generated previously [15] were differentiated into the neural lineage as described (Fig. 1A). Neural induction led to a stable population of proliferating eNPCs that homogeneously expressed the neural nuclear marker SOX-2 (Fig. 1B). Further differentiation of eNPCs, achieved by bFGF and EGF withdrawal, allowed the progressive generation of equine neurons (eNe) over time, as revealed by immunostaining for two neuronal markers, βIII-Tubulin and huC/D (Fig. 1D, 1E). Fourteen days after growth factors withdrawal, however, the differentiation process was incomplete, as progenitor cells (SOX-2 positive) were still numerous in the vicinity of young neurons (SOX-2/βIII-Tubulin-positive) (Fig. 1E), and glial cells were not detected upon immunostaining for markers of astrocytes (GFAP) and oligodendrocytes (OLIG-2) (not shown). Also, extensive cell death, as revealed by the presence of numerous floating cells, occurred at that time, therefore precluding prolongation of the differentiation process.
We thus next assessed the permissivity of equine brain cells at day 0 (eNPCs) and day 14 (equine differentiated neural cells-eDNCs) of differentiation. Equine NPCs and eDNCs were infected with WNVNY99 at MOI 10− 1 for 48h and 96h, respectively. Immunostaining with anti-WNV-E3 antibody revealed clusters of massively infected eNPCs (Fig. 2A, left) whereas infection of eNe, albeit detected, was rare (Fig. 2A, right). We therefore decided to use eNPCs for antiviral screening and aimed at further characterizing their infection while defining optimal conditions for screening. Kinetic studies performed between 24h and 72h with WNVNY99 at MOI 10− 1 demonstrated virus replication and spread in eNPCs, as shown by immunofluorescence labeling with WNV-E3 antibody (Fig. 2B) and quantification of viral RNA and viral titers in supernatants by RT-qPCR (Fig. 2C) and end-point dilution (Fig. 2D), respectively. Dose response studies performed with WNVNY99 at MOI ranging from 10− 4 to 1, associated to automated quantification of the percentage of infected cells and total cell number, showed dose-dependent changes in viral infection (Fig. 2E) and survival (Fig. 2F), with deep alteration of eNPCs growth and survival at the highest MOI. Of note, at MOI 10− 1, increase in viral infection from 24hpi to 72hpi was similarly detected by RT-qPCR, titration and image analysis (Fig. 2C-E), establishing that image-based analysis is a suitable method to quantify viral infection in eNPCs. Based on these results, we established the screening conditions (MOI 10− 2 for 48h) such as there is a substantial percentage of infected cells (approximately 45%) with no impact on cell growth or survival.
Phenotypic screen using eNPCs identifies compounds with antiviral and proviral properties
WNVNY99-infected eNPCs were used to screen a library of 41 chemical compounds selected for their antiviral activity against human and equine viruses (Supplementary Table 1), as schematically represented (Fig. 3A). Toxicity and antiviral effect were determined by automated quantification of total cells labeled with DAPI and infected cells immuno-labelled with WNV-E3 antibody, respectively. A hit was arbitrarily defined as a compound inducing a reduction of at least 25% of infected cells and less than 20% cell loss. Of all compounds tested at 10 µM (Fig. 3B), 39% (16/41) were toxic, suggesting that eNPCs were particularly sensitive to drugs. Among the non-toxic molecules (25/41), 21 (amantadine, capecitabine, DMXAA,eflornithin, favipiravir, herpes virus and reverse transcriptase inhibitors, isatin, maribavir, nelarabine and sofosbuvir) had no antiviral activity against WNV. Three coumpounds, 2’C-methylcytidine (2’-CMC), arbidol and ribavirine, reduced the percentage of infected cells to 35.3+/-7.5%, 69.7+/-21.3% and 56+/-20.1%, respectively, and were thus considered to be hits, and surprisingly, one of the compounds, atorvastatin, induced an increase in the percentage of infected cells (279.2+/-60.6%), revealing a pro-viral effect. For compounds exerting toxicity at 10µM, an additional screen was performed at 1µM (Fig. 3C). Although cytotoxicity was generally reduced, cellular loss remained above 20% for all but two compounds, fludarabine and 25-hydroxycholesterol, which nonetheless showed no antiviral activity. Of note, mycophenolic acid and brequinar displayed strong antiviral activity, albeit inducing 50% or greater cell loss (Fig. 3C). These results are summarised (Table 1). In order to verify the effect of the 4 compounds identified as modulators of WNV infection, we next quantified viral RNA and infectious viral particles in supernatants of eNPCs treated or not with 10 µM of 2’-CMC, arbidol, ribavirin or atorvastatin (Fig. 4). 2’-CMC and arbidol induced a decrease in both viral RNA (Fig. 4A, B) and viral titers (Fig. 4E, F), confirming their antiviral impact. This was not the case of ribavirin, for which a decrease in viral RNA was not detected (Fig. 4C), despite a decrease in virus titer (Fig. 4G). The proviral effect of atorvastatin was confirmed when quantifying both viral RNA and infectious viral particles (Fig. 4D, H). Thus, our newly developed screen based on image analyses allowed efficient identification of molecules that inhibit or promote WNV replication in equine brain cells, as well as simultaneous assessment of their toxicity.
Table 1
Antiviral and toxicity activities of compounds on WNV-infected eNPCs.
Compounds | Type | Viral replication | Cytotoxicity (10µM) |
2’-C-methylcytidine | DAA-NA | AV | |
25-hydroxycholesterol | HDA | | Tx |
Abacavir | DAA-RTI | NoA | |
Adefovir dipivoxil | DAA-HVI | | Tx |
Amantadine | HDA | NoA | |
Arbidol | HDA | AV | |
Atorvastatin | HDA | PV | |
Brequinar | HDA | | Tx |
Capecitabine | N/A | NoA | |
Cidofovir | DAA-HVI | | Tx |
Cladribine | N/A | | Tx |
Clofarabine | N/A | | Tx |
Cytarabine | N/A | | Tx |
Decitabine | DAA-NA | | Tx |
Didanosine | DAA-RTI | NoA | |
DMXAA | HDA | NoA | |
Eflornithin (dfmo) | HDA | NoA | |
Emtricitabine | DAA-RTI | NoA | |
F83233 | HDA | | Tx |
F83233RS | HDA | | Tx |
Famciclovir | DAA-HVI | NoA | |
Favipiravir | DAA-NA | NoA | |
Fludarabine | DAA-NA | | Tx |
Fluorouracile | DAA-NA | | Tx |
Gemcitabine | DAA-NA | | Tx |
Isatin | HDA | NoA | |
Lamivudine | DAA-RTI | NoA | |
Maribavir | DAA | NoA | |
Mercaptopurine | DAA-NA | | Tx |
Mycophenolic acid | HDA | | Tx |
Nelarabine | N/A | NoA | |
Penciclovir | DAA-HVI | NoA | |
Proguanil | HDA | NoA | |
Ribavirin | DAA-NA, HDA | AV | |
Sofosbuvir | DAA-NA | NoA | |
Stavudine | DAA-RTI | NoA | |
Telbivudine | DAA-RTI | NoA | |
Tenofovir disoproxil | DAA-RTI | NoA | |
Thioguanine | DAA-NA | | Tx |
Valaciclovir | DAA-HVI | NoA | |
Valganciclovir | DAA-HVI | NoA | |
DAA, direct-antiviral activity. HDA, host-directed antiviral. NA, nucleoside analog. RTI, reverse transcriptase inhibitor. HVI, herpes virus inhibitor. Tx, toxic. AV, antiviral activity. NoA, no activity. PV, pro-viral activity. |
Selectivity index for 2’-CMC, arbidol and ribavirin
Using the experimental design described (Fig. 3A), dose responses were evaluated, and IC50, CC50 and SI (CC50/IC50) determined for 2’-CMC, arbidol and ribavirin. As shown (Fig. 5), each drug was effective in the 10 micromolar range, with IC50 being 11 ± 1.7 µM, 15 ± 0.3 µM and 11.1 ± 1.8 µM for 2’-CMC (Fig. 5A), arbidol (Fig. 5B), and ribavirin (Fig. 5C), respectively. 2’-CMC presented the highest SI (5.3) and arbidol the lowest (1.2), revealing for the latter a toxicity in the same range of concentrations as antiviral activity.
Atorvastatin, simvastatin, lovastatin and fluvastatin have a pro-viral effect on WNV-infected eNPCs.
The proviral effect of atorvastatin in WNV-infected eNPCs raised the question of whether other statins may act similarly. We thus infected eNPCs with WNVNY99 at MOI 5x10− 3 for 48hpi (approximately 25% of cells were infected) and treated them as previously described with 3 additional statins: fluvastatin, simvastatin and lovastatin, all at 10µM. All statins induced an increase in the percentage of infected eNPCs compared with non-treated cells (Fig. 6A, B). Confirmation of their proviral role was obtained for all by quantification of viral RNA (Fig. 6C) and infectious viral particles (Fig. 6D) in the supernatants. In these latter experiments, fluvastatin was used at 1 µM (Fig. 6C, D), in order to avoid any bias due to toxicity when used at 10 µM (Fig. 6B).
Statins have no effect or an anti-viral effect on WNV-infected VERO, A549 and human NPCs.
Given that statins have been described to have broad spectrum antiviral activity [24], the observation of a pro-viral effect in WNV-infected eNPCs was striking. To clarify their role, we further assessed statin activity in 2 cell lines (VERO and A549) and in NPCs of human origin (hNPCs) using an experimental design similar to that described in Fig. 3A. Of all statins tested (at 10µM), an effect in VERO (Fig. 7A-C) and A549 (Fig. 7D-F) cells was observed only for fluvastatin, which induced a decrease of 61% and 37% in the percentage of infected cells, respectively (Fig. 7A,E). This was confirmed by quantification of virus titer in supernatant (Fig. 7C,F). In hNPCs, all four statins exerted an antiviral effect, with atorvastatin having the strongest impact, as determined by fluorescent microscopy and quantification of infected cells (Fig. 7G, H). For all statins except simvastatin, these findings were confirmed by quantification of virus titers in supernatants (Fig. 7I). Thus, our results revealed differential effects of statins depending on cellular type and species, with a proviral effect that is specific to NPCs of equine origin.