ORIGINAL RESEARCH article

Front. Microbiol., 26 March 2019

Sec. Virology

Volume 10 - 2019 | https://doi.org/10.3389/fmicb.2019.00596

Zika Virus Infection Disrupts Astrocytic Proteins Involved in Synapse Control and Axon Guidance

  • 1. Department of Medical Microbiology and Infectious Diseases, University of Manitoba, Winnipeg, MB, Canada

  • 2. Manitoba Centre for Proteomics and Systems Biology, Winnipeg, MB, Canada

  • 3. Children’s Hospital Research Institute of Manitoba, Winnipeg, MB, Canada

Abstract

The first human Zika virus (ZIKV) outbreak was reported in Micronesia in 2007, followed by one in Brazil in 2015. Recent studies have reported cases in Europe, Oceania and Latin America. In 2016, ZIKV transmission was also reported in the US and the World Health Organization declared it a Public Health Emergency of International Concern. Because various neurological conditions are associated with ZIKV, such as microcephaly, Guillain-Barré syndrome, and other disorders of both the central and peripheral nervous systems, including encephalopathy, (meningo)encephalitis and myelitis, and because of the lack of reliable patient diagnosis, numerous ongoing studies seek to understand molecular mechanisms underlying ZIKV pathogenesis. Astrocytes are one of the most abundant cells in the CNS. They control axonal guidance, synaptic signaling, neurotransmitter trafficking and maintenance of neurons, and are targeted by ZIKV. In this study, we used a newly developed multiplexed aptamer-based technique (SOMAScan) to examine > 1300 human astrocyte cell proteins. We identified almost 300 astrocyte proteins significantly dysregulated by ZIKV infection that span diverse functions and signaling pathways, including protein translation, synaptic control, cell migration and differentiation.

Introduction

More than 3 billion people are at risk for infection by arthropod-borne viruses (arboviruses). Zika virus (ZIKV) is a newly re-emerging mosquito-borne virus that is known to cause rare but serious birth defects (De Carvalho et al., 2016; Teixeira et al., 2016). ZIKV is a flavivirus that contains a positive single-stranded RNA genome that encodes three structural (capsid, pre-membrane and envelope) and seven non-structural proteins (Peng et al., 2018). ZIKV was first isolated in 1947 from sentinel rhesus monkeys in the Zika forest near Entebbe, Uganda and has since spread to the Pacific and the Americas (Ioos et al., 2014; Barzon et al., 2016). ZIKV is transmitted by the mosquito vectors Aedes aegypti and A. albopictus. A. aegypti feeds on the blood of many hosts and in the process is able to spread ZIKV to numerous individuals, making it the main mosquito species involved in ZIKV transmission (Glover and Coombs, 2017). Sexual transmission also has been reported (Moreira et al., 2017).

The serious brain malformations in newborn babies infected with ZIKV in utero include microcephaly, cortical thinning and blindness, and reflect the ability of ZIKV to readily infect the developing human brain. ZIKV can successfully target and replicate in cells in different brain regions, including brain stem/progenitor cells of the proliferating ventricular zone (Gabriel et al., 2017), oligodendroglia and microglia, whereas neurons are less susceptible (Retallack et al., 2016). Recently, it was shown that ZIKV efficiently infects numerous cell types ranging from astrocytes, neural progenitor cells and microglia-like cells to Sertoli cells and human umbilical vein endothelial cells (Kozak et al., 2017; Siemann et al., 2017; Muffat et al., 2018; Peng et al., 2018). Astrocytes are a major brain cell population that are considered an essential first site of ZIKV infection (van den Pol et al., 2017). Astrocytes act as initiators and drivers of ZIKV infection within the developing brain. Thus, given the ability of this virus to infect not only brain cells and brain cell precursors, and its ability to penetrate the placental barrier, ZIKV infection results in various neurological conditions associated with the central nervous system (CNS) [(meningo)encephalitis, myelitis and acute disseminated encephalomyetis (ADEM)] and with the peripheral nervous system (PNS) (Guillain-Barré syndrome, chronic inflammatory demyelinating polyneuropathy and acute transient polyneuritis) (Zhang et al., 2017; Leonhard et al., 2018; Mehta et al., 2018). Many of these clinical abnormalities are associated with problems occurring during the development of the brain in the fetus or in adult CNS and PNS. Of these two, CNS abnormalities have been associated with problems involving astrocytes (Wishart et al., 2006). Since astrocytes are highly involved in the survival of neurons, in axonal and synaptic guidance, and in coordination with other microglia in postnatal synaptogenesis, astrocytic death or functional disruption can result in numerous neurological conditions (Jacobsen and Miller, 2003; Gonzalez-Perez and Quinones-Hinojosa, 2012; Reemst et al., 2016). Astrocytes also migrate and assist in the migration of new neurons via the organization of the rostral migratory stream during early developmental stages (Bozoyan et al., 2012).

Since viruses are obligate intracellular parasites, they depend upon host cells and processes for replication. Numerous studies have demonstrated profound virus-induced alterations in the transcriptomic and proteomic profiles of virus-infected cells and organisms (for examples: Kobasa et al., 2007; Vester et al., 2009). Similar recent studies have identified ZIKV-induced changes in expression of large numbers of cellular genes and proteins. For example, Garcez et al. (2016) examined proteomic and transcriptomic responses to ZIKV infection in neurospheres and Kumar et al. (2018) measured transcriptomic alterations induced in Sertoli cells. Other global proteomic non-biased mass spectrometry-based approaches have been performed on ZIKV-infected C6/36 mosquito cells (Xin et al., 2017) and on neural progenitor cells (Jiang et al., 2018). Scaturro et al. (2018) recently performed an integrated affinity and phosphoproteomic analysis of human neural progenitor cells. Despite mass spectrometry being a powerful tool to identify and measure large numbers of proteins (reviewed in Yates et al., 2009; Coombs, 2011), drawbacks include that low abundant proteins may not be detected and there often is poor overlap between experimental replicates. We recently used a complementary newly developed multiplexed aptamer-based technique to measure > 1300 proteins in monkey Vero cells that had been infected with ZIKV for different amounts of time to determine temporal patterns of ZIKV-induced host protein dysregulation (Glover et al., 2019). Here, we used these same slow off-rate modified aptamers (SOMAmers®; SomaLogics, Inc., Denver, CO, United States; SOMAScan, 2017) to assay alterations in the levels of ZIKV-induced protein expression in human glioblastoma U251 cells as a model for these more physiologically relevant human astrocyte cells. We identified 170 unique U-251 proteins that were significantly over- or under-expressed at least 1.5-fold compared to non-infected that belong to several broad cellular networks, including apoptosis and pathways important for astrocytic cell synaptic control, cell growth and migration.

Materials and Methods

Cells and Viruses

Cells

Human U-251 glioblastoma astrocytoma [U-251 MG (formerly known as U-373 MG; European Collection of Authenticated Cell Cultures – ECACC 09063001)] were cultured in Dulbecco’s modified Eagle’s/Nutrient Mixture F-12 (DMEM/F-12) medium. Monkey Vero cells (ATCC® # CCL-81TM) were cultured in Dulbecco’s modified Eagle’s medium (DMEM). Both media were supplemented with 10% Fetal Bovine Serum (FBS), non-essential amino acids, 2 mM L-glutamine, and sodium pyruvate, respectively, and cells were maintained at 37°C in 5% CO2. Cells were trypsinized and passaged every 2 – 3 days.

Virus

The Asian strain of ZIKV was a gift from Dr. David Safronetz, Chief of Special Pathogens, the National Microbiology Laboratory, Public Health Agency of Canada. Virus was propagated in Vero cells in DMEM supplemented with 2.5% FBS, non-essential amino acids, 2 mM L-glutamine, sodium pyruvate, 2× gentamicin and 2× Amphotericin B at 37°C in 5% CO2.

ZIKV Infections

Stock Virus Growth

Vero cells were grown to 60–70% confluency in fully supplemented DMEM media and infected with ZIKV at MOI ∼ 0.001 PFU/cell. Media were replaced after virus adsorption with DMEM supplemented as above but containing 2.5% FBS. Supernatants were collected at different time-points to harvest virus, centrifuged at 600 ×g for 8 min at 4°C to remove cell debris, supplemented to 20% FBS and frozen at -80°C until used later for titration and U251 infections.

ZIKV Titration

Vero cells were grown to 60–70% confluency in 6-well plates. Serial 10-fold dilutions of viral stocks were made in gel saline [137 mM NaCl, 0.2 mM CaCl2, 0.8 mM MgCl2, 19 mM HBO3, 0.1 mM Na2B4O7, 0.3% (wt/vol) gelatin], and adsorbed in duplicate on the Vero cell monolayers. After viral attachment (2 h in a 37°C CO2 incubator with periodic rocking), 2.5 ml per well of a 50:50 (vol/vol) mixture of [1.2% Type I Agarose (Difco Laboratories, Detroit, Mich.)] and [2 × completed Medium 199 (Medium 199 with 6% FBS, 4 mM glutamine, 10 μg of gentamicin sulfate per ml, and 3 μg of amphotericin B per ml)] was added to each well to overlay cells. Cells were incubated at 37°C for 5 days and then overlaid with 1% Bacto Agar in 1× PBS that contained 0.04% neutral red. Plaques were counted 18–21 h after neutral red staining.

Experimental Infections

U-251 cells were grown to 60-70% confluency in fully supplemented DMEM/F-12 media and infected at MOI = 3, which, by Poisson calculation, suggests > 95% of cells are infected. Virus was adsorbed in gel saline as described above for 2 h in a 37°C CO2 incubator with periodic rocking. Mock-infected cells were treated similarly in gel saline that contained no virus. After adsorption, cells were overlaid with DMEM/F12 supplemented as described above but containing 2.5% FBS and incubated at 37°C for various periods of time as indicated below.

Photomicroscopy

Bright-Field Microscopy

Infected- and mock-treated U-251 cells were microscopically observed for cytopathic effect (CPE) at 6, 12, 24, 48, and 72 h post-infection (hpi) using a Nikon TE-2000 inverted microscope. All images were taken with a Canon A-700 digital camera and exported into PowerPoint for minimal adjustment of brightness and contrast.

Immunofluorescence Microscopy

U-251 cells were seeded on spotted slides at 20,000 cells/spot, infected at an MOI = 3 as described above and processed at 12, 24, and 48 hpi. Cells were fixed using 4% paraformaldehyde for 15 min, followed by permeabilization with 0.1% Triton X-100 in PBS for 5 min. Cells were then blocked in 3% BSA blocking solution for 90 min, followed by overnight incubation at 4°C with mouse anti-ZIKV NS1 primary antibody (BioFront Technologies; cat # BF-1225-06) in 1% BSA and PBS. Cells were washed with PBT buffer (PBS with 0.2% Tween-20) and secondary anti-mouse antibody tagged with Alexa FluorTM 488 (Invitrogen #A-11001) was added for 1 h. Slides were washed 3× with PBT buffer followed by mounting using DAPI-Prolong® Gold AntiFade to stain nuclei. Slide edges were sealed using nail polish and the slides kept at 4°C. until visualized. Images were taken with a Zeiss Axio Observer ZI inverted microscope using 20× objective. AxioVision 4.8.2 software was used to capture and edit the images before exporting them.

Cell Viability

Cell viability was measured using the WST-1 assay according to manufacturer’s instructions (Roche), but adding 8 μL of WST-1 reagent instead of 10 μL into each well of a 96-well plate at desired times. The cells were incubated for 2 h at 37°C and then absorbance was measured at 440 and 610 nm. Absorbance at 610 nm was subtracted from that at 440 nm and the values of the ZIKV-infected samples were normalized with their respective mock samples at each time point. A minimum of triplicate experiments were analyzed.

Protein Quantification

Mock- and ZIKV-infected U-251 samples were harvested at each of various time points (12, 24, and 48 hpi) and washed 3× with >50-volumes of PBS to remove media and culture FBS. Washed cells were lysed with MPER® (Pierce; Rockford, IL, United States) supplemented with 1× HALT® Protease inhibitor (Pierce; Rockford, IL, United States). After cell lysis, lysates were centrifuged at 14,000 ×g for 15 min at 11°C to remove insoluble cellular components. Cell lysate protein quantities were determined by BCATM Protein Assay (Pierce; Rockford, IL, United States) and normalized to bovine serum albumin standards.

SOMAScan® Analyses

Protein concentrations of BCA-determined cell lysates were adjusted to 200 ng/μL, and 70 μL of each sample submitted for SOMAScan® analysis in-house on a SomaLogics®-licensed platform in the Manitoba Centre for Proteomics and Systems Biology as described (SOMAScan, 2017 #14453; Glover et al., 2019). Briefly, the SOMAScan assay is a novel proteomic tool that uses single-stranded DNA-based Slow Off-rate Modified Aptamer reagents (SOMAmers). These chemically modified nucleotides were selected based upon their capacity to bind to specific human proteins. The SOMAmers capture proteins in their native state, and, after a series of washing steps, are released and their quantities measured on DNA microarray chips. When the SOMAmers are used to probe a range of sample concentrations, they are capable of measuring femtomolar to micromolar quantities of proteins. We used the SOMAScan version 1.3, capable of simultaneously measuring 1,307 distinct proteins in each of up to 92 samples (SOMAScan, 2017). Three biologic replicates of infected samples collected at 12, 24, and 48 hpi, and of each time-matched mock-infected control (=18 total samples) were simultaneously analyzed in a single SOMAScan 96-well plate. Results were reported in Relative fluorescent units (RFU) for each sample, which are directly proportional to the amounts of target protein quantities in the initial samples, as confirmed by a standard curve generated for each protein-SOMAmer pair (SOMAScan, 2017). RFU differences between Mock and ZIKV-infected were analyzed as described below in Statistical and bioinformatics analyses.

Immunoblotting

BCA-quantified protein samples were adjusted to load 10 μg of protein per gel lane. Samples were heated to 95°C for 5 min and resolved by mini-10% sodium dodecyl sulfate polyacrylamide electrophoresis (SDS-PAGE) until the loading dye had just run off the bottom edge of the gel. The proteins were transferred to PVDF [Immobilon-P polyvinylidene difluoride membrane (Millipore)] for 2 h in ice-cold buffer, followed by overnight blocking of the membrane in 5% skim milk in 1× TBST. Primary antibodies were added to each blot at 1:1000 dilution in 1% milk/TBST overnight. Primary antibodies used were: mouse monoclonal anti-ZIKV NS1 (BioFront Technologies # BF-1225-06); anti-ZIKV NS3 (GeneTex Inc. # GTX133309), anti-ZIKV envelope protein (GeneTex Inc. # GTX133314), anti-STAT3 (Cell Signaling # 124H6), anti-β-actin (Cell Signaling # 8H10D10); rabbit polyclonal anti-fibronectin (Abcam # ab2413), anti-SPARC (Cell Signaling # 5420S), and anti-cystatin C (Abcam # ab109508). After overnight binding, membranes were washed 3× with TBST and appropriate goat HRP-conjugated anti-rabbit (Cell Signaling # 7074) or anti-mouse (Cell Signaling # 7076) secondary antibodies were added for 1 h. The blots were washed 3 additional times with 1× TBST, developed with ECL western blotting peroxidase substrate for chemiluminescence and imaged with an enhanced chemiluminescence (ECL) detection machine (Amersham-Pharmacia Biotech); ImageJ was used to analyze each blot and each band in each blot was normalized to its respective actin control and to its time-matched mock-infected band intensity.

Statistical and Bioinformatics Analyses

Relative fluorescent units values for each of the 1322 analytes (1307 unique proteins + 15 internal SOMA control analytes) in each of 3 biologic replicates, each consisting of a ZIKV-infected sample and a time-matched non-infected mock sample at 12, 24, and 48 hpi were imported into Excel and converted to Log2 values. Fold-changes were determined for each of the nine infected samples compared to their time-matched mock samples. The fold-changes were analyzed for significance by both Students T-test with 2 tails, and by Z-score analysis, as described (Coombs et al., 2010; Glover et al., 2019). Briefly, all fold-changes not deemed to be significant by T-test were examined by Z-score, expressing each value as its number of standard deviations away from the population mean. Each protein’s Z-score was considered significant if the average Z-score for that protein was > 1.96σ or <-1.96σ; and if the Z-score was < -1.96σ in each of 2 or more replicates and was < -0.98σ in no more than a single replicate, or if the Z-score was > 1.96σ in each of ≥ 2 replicates and was > 0.98σ in the remaining replicate. For increased stringency, we also applied a fold-change cut-off of 1.50-fold dysregulation (= ≥ 1.50-fold if up-regulated, or ≤0.667-fold if down-regulated) to those proteins considered significantly dysregulated. Fold changes and P-values were imported into DAVID, Network Analyst and Ingenuity Pathway Analysis (IPA®) for additional bioinformatics and pathway analyses. Western blot data were examined for significance by one-way ANOVA, using a significance cut-off of 0.05.

Results

ZIKV Successfully Establishes Infection in U-251 Cells

Unlike the mock-treated astrocytic U-251 cells, ZIKV-infected astrocytes show a visible cytopathic effect (CPE) over time as they start to circularize and lift from the surface (Figure 1A), and increasing CPE with time post-infection was measured by cell viability assays (Figure 1B). ZIKV also promotes increasing production of its viral non-structural proteins, NS1 and NS3, along with structural protein, envelope, as time progresses (Figure 1C). Densitometric analyses of each of these viral proteins shows the first signal by 12 hpi (Figure 1D). NS3 protein gives the strongest signal at 24 hpi (Figure 1D). By 48 hpi, all tested viral proteins show 4 – 10-fold increase in expression (Figure 1D). We also used immunofluorescence microscopy to probe expression of viral non-structural protein NS1, which is only expressed during successful infection, to confirm that >90% of cells were infected by 48 hpi (Figure 1E), which allows subsequent meaningful proteomic analyses.

FIGURE 1

ZIKV Alters the Proteomic State of U-251 Cells With Increasing Amounts of Dysregulation as the Time After Infection Increases

We analyzed the cellular proteome of ZIKV-infected U-251 astrocytes at 12, 24, and 48 hpi and compared them to time-matched mock-treated cells using the novel SOMAScan platform. A total of 296 cellular proteins were significantly over- or under-expressed compared to mock, non-infected, with 17 over-expressed and 17 under-expressed at 12 hpi, 23 over- and 40 under-expressed at 24 hpi, and 45 over- and 209 under-expressed at 48 hpi (Table 1). For further analysis, we also consider each protein’s fold-change along with its significance. We set the cut-offs as P-value < 0.05 and fold-change of ≥ 1.50 [over-expressed to ≥ 1.50 (≥0.585 Log2) or under-expressed to ≤0.667 (≤-0.585 Log2)]. This cut-off value was chosen to keep the numbers of significant differentially expressed molecules large enough to perform meaningful bioinformatics analyses, while applying sufficient stringency to maintain physiologic relevance. Using these criteria, no proteins were significantly over- or under-expressed at 12 hpi, 2 proteins were significantly up-regulated and 17 were significantly down-regulated at 24 hpi, and 18 proteins were significantly up-regulated and 151 were significantly down-regulated at 48 hpi (Table 1 and Figure 2A). The 170 unique proteins significantly over- or under-expressed ≥ 1.50-fold at any of the tested time points (12, 24, or 48 hpi) are listed and bolded in Table 2.

Table 1

Number that are significantTotal Unique12 hpi24 hpi48 hpi
and fold-change > 1.0002961723  45
and fold-change < 0.99991740209
and fold-change > 1.100272  221  43
and fold-change < 0.9091  840209
and fold-change > 1.250244  010  30
and fold-change < 0.8000  337204
and fold-change > 1.5001700218
and fold-change < 0.6667017151
and fold-change > 1.750113  0  1    9
and fold-change < 0.5714  0  4102
and fold-change > 2.000  80  0  1    7
and fold-change < 0.5000  0  2  73

Numbers of significantly dysregulated ZIKV-infected U-251 proteins.

Significance was determined by T-test and Z-score as detailed in Section “Materials and Methods.” Values that are bolded correspond to significant differential expression and fold-change ≥ 1.50 that are listed in Table 2 and that were used for subsequent bioinformatics analyses.

FIGURE 2

Table 2

12 hpi
24 hpi
48 hpi
Entrez Gene SymbolProteinFold ChangeP-valueFold ChangeP-valueFold ChangeP-value
Significantly up-regulated proteins
CXCL11C-X-C motif chemokine 110.8650.0052.7630.0058.6144.4E-05
NACANascent polypeptide-associated complex subunit alpha1.0400.2811.5810.0060.9740.761
CCL5C-C motif chemokine 51.0410.1781.3160.0324.7212.1E-04
CKBCreatine kinase B-type0.9890.8600.9880.9313.8470.001
SSRP1FACT complex subunit SSRP10.9950.9681.4220.0402.9513.9E-04
ISG15Ubiquitin-like protein ISG150.9980.9821.2610.0392.7490.006
HMGN1Non-histone chromosomal protein HMG-141.0160.8151.1690.1432.3120.002
NUDCD3NudC domain-containing protein 30.9860.7381.1350.1252.0550.001
PLCG11-phosphatidylinositol 4,5-bisphosphate phosphodiesterase gamma-11.0150.7751.1830.0431.8180.004
METAP2Methionine aminopeptidase 20.9920.9471.1180.0761.8120.003
ADSLAdenylosuccinate lyase1.0940.1011.1920.1431.6940.002
PARK7Protein deglycase DJ-11.0280.3261.2250.0571.6650.009
HSPA1AHeat shock 70 kDa protein 1A1.1600.0421.2240.0191.6420.011
GDI2Rab GDP dissociation inhibitor beta1.0650.5121.1140.1111.5930.002
PEBP1Phosphatidylethanolamine-binding protein 10.9930.9421.1690.1891.5640.003
UBE2NUbiquitin-conjugating enzyme E2 N1.0470.5851.1010.1971.5570.001
PPA1Inorganic pyrophosphatase1.0350.6841.1390.0751.5340.002
IFNL1Interferon lambda-10.9080.0601.2390.0171.5120.016
EIF5AEukaryotic translation initiation factor 5A-11.0400.6531.1430.0391.5050.003
Significantly down-regulated proteins
STC1Stanniocalcin-10.9020.4370.3630.0020.0692.3E-04
CST3Cystatin-C0.9130.4450.4820.0010.2320.001
IGFBP2Insulin-like growth factor-binding protein 20.9760.6890.5370.0010.2134.6E-04
SFRP1Secreted frizzled-related protein 10.9040.2110.5620.0070.1373.5E-04
APPAmyloid beta A4 protein0.9760.7860.5850.0030.1541.5E-04
MFGE8Lactadherin0.9460.4100.5870.0030.2290.001
MDKMidkine1.0190.6290.5880.0040.1990.002
IGFBP7Insulin-like growth factor-binding protein 70.9190.2560.6110.0100.2370.001
IGFBP5Insulin-like growth factor-binding protein 50.9820.7670.6300.0060.1452.8E-05
CTSALysosomal protective protein0.9460.3530.6370.0050.2810.001
LRIG3Leucine-rich repeats and immunoglobulin-like domains protein 30.9780.6880.6380.0090.2032.2E-04
FGFR1Fibroblast growth factor receptor 10.9370.2910.6480.0110.2003.0E-04
PCSK9Proprotein convertase subtilisin/kexin type 91.0120.8400.6530.0170.2212.2E-04
CFIComplement factor I0.9110.3420.6550.0120.1700.001
TIMP2Metalloproteinase inhibitor 20.9990.9860.6570.0040.3273.3E-04
FSTL1Follistatin-related protein 11.0120.7690.6580.0100.2610.001
GRNGranulins0.9610.5780.6620.0040.1981.7E-04
HS6ST1Heparan-sulfate 6-O-sulfotransferase 10.9640.2410.7100.0070.1962.5E-04
C4A C4BComplement C40.9160.2651.0440.7110.2460.004
PLAUUrokinase-type plasminogen activator0.9720.5300.7790.0410.2530.001
TNCTenascin0.9730.4800.7320.0100.2574.4E-04
GPNMBTransmembrane glycoprotein NMB1.0140.7970.9140.4190.2580.002
GFRA1GDNF family receptor alpha-10.9710.6380.7440.0160.2670.001
LAMA1 LAMB1 LAMC1Laminin0.9490.3190.7620.0190.2770.001
THBS2Thrombospondin-20.9300.1330.7100.0160.2900.001
ACVR1BActivin receptor type-1B1.0010.9901.1420.1940.2930.002
XRCC6X-ray repair cross-complementing protein 60.9360.2220.9060.3570.2933.9E-04
GNSN-acetylglucosamine-6-sulfatase1.0040.9480.8780.1790.3090.001
CCL13C-C motif chemokine 131.0030.9581.0840.3310.3260.001
CCL20C-C motif chemokine 200.9690.6861.1000.4210.3270.003
L1CAMNeural cell adhesion molecule L10.9000.1870.8800.0700.3360.001
IGF1RInsulin-like growth factor 1 receptor1.0180.6240.9720.6690.3380.005
ICAM5Intercellular adhesion molecule 50.9710.4940.9390.4280.3380.002
PLXNB2Plexin-B20.9740.6040.8260.0310.3420.001
GSNGelsolin0.9130.3120.8990.2940.3440.003
C4A C4BComplement C4b0.9580.5701.1150.3360.3520.001
IL17RAInterleukin-17 receptor A0.9510.5410.9300.4450.3610.001
LYNTyrosine-protein kinase Lyn, isoform B0.9560.6350.9040.4500.3640.002
PF4Platelet factor 40.9810.7021.1180.1920.3740.002
PROK1Prokineticin-10.9860.8061.1390.1360.3850.003
HIST1H1CHistone H1.20.9230.2890.8210.1790.3900.005
RTN4RReticulon-4 receptor0.9930.7720.9060.3190.3920.002
UNC5CNetrin receptor UNC5C1.0040.8750.8370.0860.3920.002
LIFRLeukemia inhibitory factor receptor0.9550.5930.6860.0120.3960.003
LYNTyrosine-protein kinase Lyn0.9250.3870.8970.4290.3980.003
FGF12Fibroblast growth factor 120.9810.7091.0610.5120.4040.002
HSD17B103-hydroxyacyl-CoA dehydrogenase type-20.9840.6880.9910.9170.4230.003
CMA1Chymase0.9890.8041.1260.1410.4240.003
CRISP3Cysteine-rich secretory protein 30.9940.9341.1220.2080.4250.002
ALBSerum albumin0.9980.9651.1040.2950.4340.001
TNFRSF1ATumor necrosis factor receptor superfamily member 1A1.0050.9040.7830.0460.4410.003
SRCProto-oncogene tyrosine-protein kinase Src0.9770.7060.9070.3530.4420.005
FERTyrosine-protein kinase Fer0.9820.7901.2260.5270.4520.006
TFPITissue factor pathway inhibitor0.9480.2160.7150.0130.4580.002
DYNLL1Dynein light chain 1, cytoplasmic0.9900.8231.0950.3080.4590.003
MMP272 kDa type IV collagenase0.9830.7290.8770.1030.4600.003
RNASEH1Ribonuclease H10.9760.7091.0060.9570.4630.003
FN1Fibronectin1.0130.7500.8510.0750.4640.003
SERPINA3Alpha-1-antichymotrypsin0.9450.1761.0530.5960.4670.003
NRXN3Neurexin-3-beta0.9830.4790.9180.2810.4680.006
PON1Serum paraoxonase/arylesterase 10.9710.5981.0980.3470.4720.004
FSTL3Follistatin-related protein 30.9100.3700.7270.0460.4730.001
NID1Nidogen-10.9770.6410.8430.0660.4750.001
B2MBeta-2-microglobulin0.9690.5800.7210.0140.4850.007
UNC5DNetrin receptor UNC5D0.9350.3620.8690.2390.4860.002
AGERAdvanced glycosylation end product-specific receptor, soluble1.0120.8541.0920.3510.4890.0046
MICBMHC class I polypeptide-related sequence B1.0090.8750.9780.7740.4890.005
ECM1Extracellular matrix protein 10.9810.6550.7630.0190.4890.003
LY9T-lymphocyte surface antigen Ly-91.0090.9251.1390.0800.4930.003
ADAM12Disintegrin and metalloproteinase domain-containing protein 120.9880.7640.9520.6160.4930.004
PTNPleiotrophin1.0310.5480.7690.0260.4960.004
HIST1H3AHistone H3.10.9940.9290.7610.0100.4981.8E-04
TFSerotransferrin0.9030.2061.1000.3860.4990.001
NOTCH1Neurogenic locus notch homolog protein 10.9870.8540.8770.2150.5060.005
CDONCell adhesion molecule-related/down-regulated by oncogenes0.9540.2631.0100.9060.5070.008
STAT3Signal transducer and activator of transcription 31.0180.7700.9030.1780.5080.003
CTSVCathepsin L20.9900.7770.7140.0490.5110.009
C5C5a anaphylatoxin1.0080.9031.1220.1720.5110.007
KIF23Kinesin-like protein KIF230.9300.1820.8840.0660.5200.006
BOCBrother of CDO0.9310.3141.0770.6190.5260.004
IL15RAInterleukin-15 receptor subunit alpha0.9920.8621.1370.1690.5280.006
GPC2Glypican-21.0050.8951.1480.3030.5320.015
SERPINE2Glia-derived nexin0.9980.9671.0580.5000.5320.006
ANGPT2Angiopoietin-20.9710.5011.1370.1360.5330.007
F11Coagulation Factor XI1.0060.8991.1100.2350.5340.006
CTSSCathepsin S1.0070.8941.0930.2590.5350.006
COL18A1Endostatin0.9920.7590.7810.0240.5360.003
PRKACAcAMP-dependent protein kinase catalytic subunit alpha1.0600.0711.0460.6080.5360.002
SPARCSPARC1.0920.1160.8670.1050.5370.004
PLATTissue-type plasminogen activator0.9920.7800.8850.1450.5400.005
PRKCAProtein kinase C alpha type0.9460.2961.0180.9250.5410.004
THBS1Thrombospondin-11.0140.7920.9050.4020.5520.010
LGALS8Galectin-80.9790.6400.8130.0310.5550.008
AKT1 AKT2 AKT3RAC-alpha/beta/gamma serine/threonine-protein kinase0.9600.5260.9730.7800.5570.003
KLK13Kallikrein-131.0110.8331.1110.2000.5600.009
POSTNPeriostin0.9880.8521.1610.2300.5610.007
MRC2C-type mannose receptor 20.9750.4551.1200.3890.5660.014
ICOSInducible T-cell costimulator1.0100.8491.0990.3050.5690.008
PDK1[Pyruvate dehydrogenase (acetyl-transferring)] kinase isozyme 1, mitochondrial1.0250.5871.0750.2440.5690.008
DKK1Dickkopf-related protein 10.9160.0200.7480.0170.5690.010
METAP1Methionine aminopeptidase 10.9940.9561.0300.8050.5700.003
DDCAromatic-L-amino-acid decarboxylase0.9500.2811.0790.3490.5700.003
MAP2K3Dual specificity mitogen-activated protein kinase kinase 31.0110.6091.0870.2690.5730.012
CAMK2DCalcium/calmodulin-dependent protein kinase type II subunit delta0.9970.9631.0860.5770.5750.009
METHepatocyte growth factor receptor1.0600.2901.0740.4120.5840.012
ALCAMCD166 antigen1.0300.5400.9730.5740.5860.003
EGFREpidermal growth factor receptor1.0530.2361.0260.5790.5900.002
LGALS3BPGalectin-3-binding protein0.9830.8160.8260.1200.5900.013
TGFBITransforming growth factor-beta-induced protein ig-h31.0160.7200.8360.0440.5900.003
PSMA2Proteasome subunit alpha type-20.9880.7821.0350.7320.5910.017
CTSHCathepsin H1.0000.9921.0710.4250.5950.013
EFNA2Ephrin-A21.0360.1950.9840.8410.6000.010
CAMK2BCalcium/calmodulin-dependent protein kinase type II subunit beta0.9890.8731.0750.6440.6000.017
EPHA3Ephrin type-A receptor 31.0050.9141.0650.3560.6000.016
FYNTyrosine-protein kinase Fyn0.9700.4921.0300.6430.6020.010
TBPTATA-box-binding protein0.9860.8190.9770.7970.6030.013
EPHB2Ephrin type-B receptor 20.9950.9141.0280.6830.6030.016
DKK4Dickkopf-related protein 40.8970.0530.7790.0320.6040.016
CSKTyrosine-protein kinase CSK1.0110.9121.0510.8220.6090.008
PLAURUrokinase plasminogen activator surface receptor0.9880.7350.9250.2770.6150.008
ADRBK1beta-adrenergic receptor kinase 10.9890.8650.9950.9410.6160.021
CTSBCathepsin B0.9780.6221.1230.1200.6160.022
HSPB1Heat shock protein beta-11.0100.8191.0070.9510.6170.011
CAMK2ACalcium/calmodulin-dependent protein kinase type II subunit alpha0.9720.5471.1200.5280.6170.020
MCL1Induced myeloid leukemia cell differentiation protein Mcl-10.9840.6790.9410.4010.6180.014
ROBO2Roundabout homolog 21.0440.1640.9730.6590.6260.011
MATN2Matrilin-20.9920.8251.0370.6490.6270.010
CSRP3Cysteine and glycine-rich protein 30.9900.8880.9020.2710.6280.011
PDE11ADual 3′,5′-cyclic-AMP and -GMP phosphodiesterase 11A1.0220.7160.9970.9200.6310.001
TK1Thymidine kinase, cytosolic0.9660.5720.9940.9620.6330.028
MAPK14Mitogen-activated protein kinase 140.9850.7991.0160.8740.6350.004
NRP1Neuropilin-10.9520.3610.8390.0940.6360.018
PIK3CA PIK3R1PIK3CA/PIK3R10.9820.8280.9720.8020.6370.024
PDPK13-phosphoinositide-dependent protein kinase 10.9670.6381.0310.8070.6390.029
APOA1Apolipoprotein A-I1.0000.9911.1060.2280.6430.020
SNX4Sorting nexin-41.0230.4771.1470.5030.6450.047
SLAMF7SLAM family member 71.0300.4981.0810.2640.6460.016
PPP3CA PPP3R1Calcineurin0.9820.7601.0050.9510.6510.009
PRTN3Myeloblastin0.9880.7851.0930.2960.6510.019
CCL15C-C motif chemokine 151.0000.9961.1210.2050.6510.014
GAS1Growth arrest-specific protein 11.0120.6751.0070.9170.6520.032
IL6STInterleukin-6 receptor subunit beta1.0210.4680.9070.2410.6550.016
LRPAP1alpha-2-macroglobulin receptor-associated protein1.0470.0840.9910.7940.6570.001
LCN2Neutrophil gelatinase-associated lipocalin0.9430.3221.1420.1380.6570.016
TNFRSF21Tumor necrosis factor receptor superfamily member 211.0230.4350.9680.7200.6570.027
CXCL6C-X-C motif chemokine 61.0290.4480.9750.7880.6570.006
MPLThrombopoietin Receptor0.9700.4811.1110.2760.6590.017
COL23A1Collagen alpha-1(XXIII) chain0.9710.5631.1500.2160.6590.042
NLGN4XNeuroligin-4, X-linked0.9780.5871.0220.8220.6640.027
SERPINE1Plasminogen activator inhibitor 11.0280.7811.1410.3410.6640.029
CCDC80Coiled-coil domain-containing protein 801.0630.1450.9690.6890.6660.020

U-251 proteins significantly dysregulated by ZIKV infection.

Values represent averages of all three biologic replicates and represent infected alterations compared to time-matched mock controls. Significance was determined by T-test and by Z-Score as detailed in Section “Materials and Methods” and values < 0.05 are bolded. Significant fold-changes > 1.50-fold at 12, 24, and 48 hpi are indicated in larger bold font in Fold-Change columns. Proteins are sorted by time point (early to late) with up-regulated at top and down-regulated below.

To validate the SOMAScan results, we selected four cellular proteins predicted to have dysregulated expression by ZIKV infection and assessed their quantities at various times post infection. Fibronectin is involved in astrocytic cell growth, migration and actin polymerization, and STAT3 and SPARC are involved in astrocytic synaptic control. Commercial antibodies are available for each of these, and also for CST3 and were used to probe ZIKV-infected astrocytes (Figure 2B). Western blot confirmed that SPARC was significantly under-expressed at 24 hpi and that all proteins were under-expressed by 48 hpi (Figure 2B,C). There were apparent differences in the degree of SPARC expression as measured by both immunoblot and SOMAScan, but the Western blot and SOMAScan results trended in the same directions for all tested proteins. Minor differences in the extent of over/under-expression were observed between the two techniques (Figure 2C), but, overall, the immunoblot data validated the SOMAScan-measured protein dysregulation.

A global overview of these proteins with altered levels of expression, using IPA, indicates that cell death functions are highly activated at 48 hpi, while cell survival, cell outgrowth, cell migration, cell spreading and cell-to-cell contact functions are inhibited at 48 hpi (Figure 2D).

ZIKV Differentially Disrupts Cellular Protein Levels and Pathways in Human Astrocytes and Monkey Kidney Cells

Lists of the ZIKV-induced differentially expressed U-251 proteins were uploaded to Ingenuity Pathway Analysis (IPA) for network analysis. Multiple networks that contain ≥ 12 focus molecules and with scores ≥ 20 were identified. The IPA-generated scores for each network are based on the fit between the network and that of the set of focus genes identified in that network; this score is a measure of likelihood of finding those identified genes in that particular network by random chance (Long et al., 2004). A number of proteins involved in post-translational modification, protein degradation and protein synthesis were significantly under-expressed, with many more significantly under-expressed at later time points (Figure 3A). Molecules significantly over- or under-expressed in this network included transmembrane receptors [(AGER; receptor for advanced glycosylation end products), PLAUR], transporters (ECM1), various enzymes like peptidases (CMA1, F11, MMP2, PLAT, PLAU, PCSK9, and PRTN3) and phosphatases (PON1) and other molecules like RNASEH1, SERPINES 1 and 2 and TFPI (Figure 3A). Most molecules involved in cellular assembly, organization, maintenance and development are under-expressed except for GDI2, SSRP1 and Rho gdi (Figure 3B). Molecules significantly under-expressed in this network included kinases (PDPK1, EPHA3, EFNA2, PRKCA, and SRC), transmembrane receptors (GFRA1 and ROBO2) and other molecules (CSRP3, GDI2, GNS, GNS, ICAM5, and SSRP1) (Figure 3B). Cellular development, cellular growth and proliferation and cell cycle networks also had numerous over-/under-expressed proteins, with the majority being under-expressed (Figure 3C). Only IFNL1 and EIF5A were significantly over-expressed by 48 hpi. Molecules significantly under-expressed included kinases (CSK and EGFR), transmembrane receptors (ICOS, IL6ST and MICB), translation regulator (EIF5A), cytokine (IFNL1) and other molecules (LY9, MATN2, SLAMF7, SNX4, and TF) (Figure 3C).

FIGURE 3

We have previously reported SOMAScan® analyses of ZIKV-infected monkey Vero cells (Glover et al., 2019). We therefore compared our differentially expressed protein lists using IPA. Three of the top U-251 networks (with ≥ 12 focus molecules and with scores ≥ 20) were: (1) Post-translational modification, protein degradation, protein synthesis; (2) Cellular assembly and organization, cellular function and maintenance, cellular development; and (3) Cellular development, cellular growth and proliferation, cell cycle. Numerous molecules within these networks were significantly under-expressed by ZIKV infection by 48 hpi. By contrast, while a few molecules (i.e., SSRP1, EIF5A, GSN, and SNX4) were similarly over- or under-expressed in both U-251 and Vero cells by 48 hpi, most of the differentially expressed U-251 proteins were not significantly affected in Vero cells (Supplementary Figure 1). Similarly, protein expression patterns impacted by ZIKV were different in the U-251 cells at 48 hpi compared to proteins present in the top two Vero cell networks (Supplementary Figure 2). For example, far more cell cycle, cell death and survival and cancer network proteins were over-expressed in Vero cell than in U-251 cells. In addition, while many proteins associated with the cancer, organismal injury and abnormalities and cellular movement network were similarly expressed in both cell types (i.e., IGFBP, CCL5 and TNFRSF21), several proteins (i.e., IGFBP3, ITGA1, and C4A/C4B) were dissimilarly impacted in the two cell types.

ZIKV Impairs Numerous Astrocytic Canonical Pathways Including Axonal Guidance Signaling, FGF Signaling, STAT3 Signaling, AMPK and ERK/MAPK Signaling

Six canonical pathways; Axonal guidance signaling, FGF signaling, Glioma Signaling, STAT3 signaling, ERK/MAPK and AMPK signaling, which are highly important in astrocytic functions, were predicted to be dysregulated by ZIKV infection by 48 hpi (Figure 4A), with Z-scores exceeding ± 1.96. Most canonical pathways were predicted by IPA to be inhibited (blue bars), based upon the program considering the expression values of all identified proteins within each pathway. However, the most significantly dysregulated canonical pathway (axonal guidance signaling), could not be assigned by the program as overall activated or inhibited despite the large numbers of differentially expressed proteins (Table 3 and Supplementary Figure 3). Many of the molecules in these pathways control two major astrocyte functions; synaptic control and astrocytic cell growth and migration (actin polymerization) (Figure 4B). Protein Kinase C alpha/beta-II type (PRKCA/PRKCB-II), CKB, Nidogen-1, NCAM1 and fibronectin are involved in astrocytic cell growth, migration and actin polymerization (Figure 4B). HSD17B10, STAT3, Thrombospondin 1/2, CAMK2A/B/D, SPARC and EPHA3 are involved in astrocytic synaptic control (Figure 4B). SOMAScan indicated that CKB was over-expressed whereas PRKCA/PRKCB-II, Nidogen-1, NCAM1 and fibronectin were under-expressed by ZIKV by 48 hpi (Figure 4B). However, HSD17B10, STAT3, Thrombospondin 1/2, CAMK2A/B/D, SPARC and EPHA3 were under-expressed by ZIKV at 48 hpi (Figure 4B). Numerous differentially expressed astrocytic proteins also are predicted to inhibit differentiation of nervous system, differentiation of astrocytes, growth of neurites, maintain radial glial cells, neuron development and neuron migration by 48 hpi (Figure 4C).

FIGURE 4

Table 3

Entrez Gene SymbolProteinFold ChangeP-valueLocationType
ADAM12ADAM metallopeptidase domain 120.490.004Plasma MembranePeptidase
EFNA2Ephrin A20.600.010Plasma MembraneKinase
EPHA3EPH receptor A30.600.016Plasma MembraneKinase
EPHB2EPH receptor B20.600.016Plasma MembraneKinase
FGFR1Fibroblast growth factor receptor 10.200.000Plasma MembraneKinase
FYNFYN proto-oncogene, Src family tyrosine kinase0.600.010Plasma MembraneKinase
L1CAML1 cell adhesion molecule0.340.001Plasma MembraneOther
METMET proto-oncogene, receptor tyrosine kinase0.580.012Plasma MembraneKinase
MMP2Matrix metallopeptidase 20.460.003Extracellular SpacePeptidase
NRP1Neuropilin 10.640.018Plasma MembraneTransmembrane receptor
PLCG1Phospholipase C gamma 11.810.004CytoplasmEnzyme
PLXNB2Plexin B20.340.001Plasma MembraneTransmembrane receptor
PRKACAProtein kinase cAMP-activated catalytic subunit alpha0.540.002CytoplasmKinase
PRKCAProtein kinase C alpha0.540.004CytoplasmKinase
PROK1Prokineticin 10.390.003Extracellular SpaceGrowth factor
ROBO2Roundabout guidance receptor 20.630.011Plasma MembraneTransmembrane receptor
RTN4RReticulon 4 receptor0.390.002Plasma MembraneTransmembrane receptor
UNC5CUNC-5 netrin receptor C0.390.002Plasma MembraneTransmembrane receptor
UNC5DUNC-5 netrin receptor D0.490.002Plasma MembraneOther

U-251 axonal guidance proteins significantly dysregulated by ZIKV by 48 hpi.

Values represent averages of all three biologic replicates and represent infected alterations compared to time-matched mock controls. Significance was determined by T-test and by Z-Score as detailed in Section “Materials and Methods.” Only those proteins whose significant fold-change was ≥ 1.50-fold in either direction (i.e., ≥1.50 or ≤0.666) are shown.

ZIKV Infection Leads to the Dysregulation of Numerous Upstream Cellular Regulator Molecules in U251 Cells by 48 hpi

Based on the over- and under-expression of a number of downstream molecules, 56 upstream molecules were identified and their activity was predicted to be significantly inhibited, and 46 upstream molecules were identified and predicted to be significantly activated by ZIKV by 48 hpi (Figure 5). These upstream molecules are predicted by the IPA program and therefore not necessarily probed by SOMAScan. Based on the over-/under-expression of each of their downstream molecules, that are probed by SOMAScan, IPA generates a prediction of the activation status of these upstream molecules. Therefore, even though the expression of some of these upstream molecules was not directly examined by the SOMAScan, IPA was able to predict their activation status (with high significance; i.e., p-values < 10-7 for inhibited and <10-6 for activated). These molecules belong to numerous categories including transmembrane receptors, transcription regulators, ligand-dependent nuclear receptors, growth factors, cytokines and G protein coupled receptors. A subset of these molecules, and the implications of their dysregulation on astrocytic regulation of nervous system organization and differentiation (Landstrom et al., 2000; Mazars et al., 2001) will be further discussed below in the Section “Discussion.”

FIGURE 5

Discussion

We employed a novel proteomic approach to assess ZIKV-induced alterations in levels of U-251 astrocyte protein expression during a time course after infection. We showed that ZIKV can successfully infect and replicate in U-251 cells and that under our infection conditions (MOI = 3), >90% of the cells were infected by 48 hpi (Figure 1E), allowing meaningful proteomic analyses. The vast majority of proteins were significantly under-expressed, especially by late time points, and many of these differentially expressed proteins belong to specific astrocyte cell functions, upstream complexes/group of molecules/transcriptional regulators, and several individual key molecules (Figure 2). A few of these cellular functions and molecules, such as inhibition of differentiation, cell growth dysregulation and modulation of stress responses, have been identified in other ZIKV-infected cell types, including Vero cells, C6/36 cells and human neural stem cells (Devhare et al., 2017; Hou et al., 2017; Glover et al., 2019).

One of the dysregulated cellular systems was the Post-translational modification, protein degradation and synthesis network (Figure 3A). Molecules significantly over- or under-expressed in this network included transmembrane receptors [AGER (Receptor for advanced glycation end products), PLAUR], transporters (ECM1), various enzymes like peptidases (CMA1, F11, MMP2, PLAT, PLAU, PCSK9, and PRTN3) and phosphatases (PON1) and other molecules like RNASEH1, SERPINES 1 and 2 and TFPI (Figure 3A). They localize in the nucleus, plasma membrane and extracellular space. Their differential expression is known to impair a number of astrocytic functions like apoptosis and astrocytic reactivity. For instance, down-regulation in the expression of AGER is associated with modulation of cell apoptosis, localization of NF-κB, phosphorylation of ERK1/2, and astrocyte reactivity leading toward neurodegeneration (Wang et al., 2002; Ohtaki et al., 2007; Choi et al., 2014; Kim et al., 2015; Gu et al., 2016). Differential expression of other molecules like MMP-2 (Matrix Matalloproteinase-2) and SERPINES 1/2 (Serpin family E member 1 and 2) results in the proteolytic cleavage of pericellular proteins and breakdown of extra-cellular matrix, stalls cell migration and neurite remodeling in the brain, which in turn can impact brain development (Ogier et al., 2006; Hagemann et al., 2012; Kamat et al., 2014; Wu et al., 2018; Yang et al., 2018).

Another dysregulated cellular system was the Cellular assembly and organization, cellular function and maintenance and cellular development network. Molecules significantly affected in this network included kinases (PDPK1, EPHA3, EFNA2, PRKCA, and SRC), transmembrane receptors (GFRA1 and ROBO2) and other molecules (CSRP3, GDI2, GNS, GNS, ICAM5, and SSRP1) (Figure 3B). These are localized in the nucleus, cytoplasm, plasma membrane and extra cellular space. They control various astrocytic cellular functions. For example, PDPK1 (3-phosphoinositide dependent protein kinase 1) dysregulated expression is known to damage glucose metabolism, cell migration, cell survival and neural plasticity (Scheid et al., 2005; Park et al., 2013). Under-expression of PRKCA (Protein kinase C alpha) is linked with modulation of apoptotic signaling and cell proliferation while Src (SRC proto-oncogene, non-receptor tyrosine kinase), which has a role in actin dynamics in U-251 cells, is associated with the disruption of cell migration (Angers-Loustau et al., 2004; Marquina-Sanchez et al., 2017). In addition, the under-expression of GFRA1 (GDNF family receptor alpha 1) by ZIKV can disrupt its binding with GDNF (glial cell-line derived neurotrophic factor), dysregulating the process of cell differentiation and proliferation indirectly (Takashima et al., 2015; Ayanlaja et al., 2018).

A third major dysregulated cellular system was the Cell development, cell growth, proliferation and cell cycle network. The molecules significantly affected in this network included kinases (CSK and EGFR), transmembrane receptors (ICOS, IL6ST, and MICB), translation regulator (EIF5A), cytokine (IFNL1) and other molecules (LY9, MATN2, SLAMF7, SNX4, and TF) (Figure 3C). These are localized in the cytoplasm, plasma membrane and extracellular space. Their abnormal expression leads to disruption of numerous astrocytic cellular functions. For example, under-expression of EGFR (Epidermal growth factor receptor) is known to result in dysregulation of the formation of reactive astrocytes and negatively impact EGFR-mediated cell migration (Liu and Neufeld, 2007). Since EGFR controls astrocytic survival and differential growth in cortical and midbrain astrocytes, its under-expression by ZIKV may result in differential apoptotic responses in various regions of the brain (Wagner et al., 2006; Liu and Neufeld, 2007) Disruption in EGFR is also linked with the disorganization of axons, fewer cell-to-cell contacts between neurons and astrocytes and apoptosis induction, resulting in the state of neurodegeneration (Wagner et al., 2006; Liu and Neufeld, 2007). In addition, EGFR serves as a coreceptor for facilitating entry of influenza A virus (IAV) and other flaviviruses, such as Hepatitis C virus (HCV) (Lupberger et al., 2013; Zheng et al., 2014), and EGFR also regulates interferons via STAT1 (Lupberger et al., 2013; Zheng et al., 2014). ZIKV also caused under-expression of STAT3 in our U-251 cells. Collectively, these imply ZIKV may modulate the EGFR signaling pathway and host proinflammatory response to promote its successful replication. The eukaryotic translation factor eIF5A is another up-regulated molecule, which is involved in translation elongation and termination beyond the ribosome pause state due to polyproline motifs (Pelechano and Alepuz, 2017; Schuller et al., 2017). This seems plausible since ZIKV needs to hijack the cellular protein translation machinery for production of its progeny proteins. Moreover, other eukaryotic translation factors that were significantly up-regulated by ZIKV infection, although only ≈30%, include eIF4H and eIF5A-1. Another aspect of eIF5A that should be considered is the hypusination process that it undergoes as part of its unique post-translational modification; this process is critical for its function and is important for mRNA stability (Olsen and Connor, 2017). Furthermore, since other flaviviruses like dengue virus are sensitive to utilizing polyamines, needed for the hypusination process (Olsen and Connor, 2017; Smirnova et al., 2018), ZIKV-induced over-expression of eIF5A in U-251 cells seems plausible. Comparative analyses with the Vero cells at 48 hpi also showed a few of these same differentially expressed proteins; namely eIF5A, C4, SSRP1, EPHA3 and PDPK1 (Supplementary Figure 1). In addition to overlaying ZIKV-infected Vero cell data on our top ZIKV-infected U-251 networks, we also overlaid U-251 data onto the top identified dysregulated Vero cell networks (Supplementary Figure 2). Once again, a number of molecules, including HSPA1A/B, PARK7, HMGN1, IGF1R, IGFBP2/5, L1CAM and CCL5 were similarly over- or under-expressed in both systems (Supplementary Figure 2). All of these impairments of astrocytic cellular differentiation, migration and effect on neurite growth via dysregulated astrocytic functions were also significantly predicted by the IPA program (Figure 4). In addition, a few of the highly crucial canonical pathways needed in order for astrocytes to function normally were also predicted to be inhibited by IPA (Figure 4).

In addition to astrocyte cell migration, reactivity and differentiation, IPA also predicted ZIKV impacting their interaction with neural cells (Reemst et al., 2016). The formation of contact between astrocytes and neural synapses is necessary for the establishment of the “tripartite synapse,” brain homoeostasis, synaptic transmission regulation and supplying substrates to the neurons (Reemst et al., 2016; Becerra-Calixto and Cardona-Gomez, 2017). Astrocytic cell movement and outgrowth are important during the developmental stage when astrocytes need to migrate to various sites in the CNS which later become their final destination (Chaboub and Deneen, 2013). Both of these functions were predicted to be downregulated by ZIKV (Figure 2D,4C). Several important individual molecules highly involved in these processes and that are highly under-expressed include protein kinase C (PRKCA/PRKCB-II), which aids adhesion of astrocytes to the local neurons; thrombospondins 1 and 2 (TSP1, TSP2), which drive synapse assembly, aid recovery of synaptic plasticity and increase excitatory synapse density in the CNS; and SPARC, which is a negative regulator of maturation of presynaptic termini and promotes synapse elimination (Figure 5B) (Hama et al., 2004; Allen, 2013; Terni et al., 2017). PRKACA and PRKCA also are involved in axonal guidance (Figure 4A, Table 3, and Supplementary Figure 3). Numerous other proteins like MMP2, EPHA3, EFNA2, and FGFR1 which are involved in cellular assembly, organization, astrocytic reactivity, and adhesion of astrocytes to local neurons also regulate the axonal guidance pathways (Table 3). Despite the fact that IPA was unable to generate a prediction for this canonical pathway as either activated or inhibited, the list of 19 highly under-expressed proteins points toward dysregulation in the overall pathway (Figure 4A). Further detailed prediction patterns were also generated based on the expression profiles of proteins involved in the axonal guidance pathway and they point toward both inhibition and activation of axon dynamics, repulsion, outgrowth, and growth cone spreading (Supplementary Figure 3). Overall, this ZIKV-induced dysregulation at both the expression level and the predicted activity level highlights another way ZIKV may disrupt astrocytic function, leading toward overall brain conditions like microcephaly. FN1 (fibronectin) and laminin are another set of ZIKV-induced differentially expressed proteins in astrocytes that are important in the elongation of axons and pathfinding via contact-mediated attraction of the growth cone (Figure 5B) (Reemst et al., 2016). Moreover, CKB (creatine kinase B-type), which was highly over-expressed by ZIKV infection, is crucial in maintaining ion homeostasis in the brain, and glutamate uptake in neurotransmitter trafficking between astrocytes and oligodendrocytes and between axons and the contact points with the nodes of Ranvier (Steen et al., 2010). This ZIKV-induced differential expression in CKB levels could also disrupt the astrocyte function of maintaining glutamate toxicity for local neurons (Steen et al., 2010). Since CKB also co-localizes with actin and helps in astrocytic cell spreading, its under-expression by ZIKV could explain impairment in cell spreading and migration (Kuiper et al., 2009). Glutamate uptake by astrocytes is also controlled by CAMKII and our SOMAScan data showed its under-expression by ZIKV infection at 48 hpi, in part explaining the dysregulation in neurotransmitter trafficking via modulation of CAMKII (Figure 5B) (Chawla et al., 2014). These CAMKII fluctuations are also known to induce neuronal death by disrupting calcium homoeostasis (Ashpole et al., 2012, 2013). There were numerous other differentially expressed targets shown in Figure 4 along with the proposed astrocytic function that could be impaired by ZIKV-induced dysregulation (Yang et al., 2009, 2014; Tyzack et al., 2014). Thus, these multiple impaired functions help explain effects which may lead toward neuronal damage, neuronal and astrocytic degeneration and overall cell death resulting in conditions like microcephaly or other disorders (Wen et al., 2017).

IPA analyses also predicted activation/inhibition of various upstream regulators that could stall a variety of cellular canonical pathways (Figure 5). They belong to numerous broad classes like transmembrane receptors, transcription regulators, growth factors, cytokines, kinases and others (Figure 5C). These upstream molecules are generated by the IPA program and are not necessarily examined by SOMAScan. Based on the over/under-expressed of each of their downstream molecules, that are probed by SOMAScan, IPA generates a prediction on the activation status of these upstream molecules. Therefore, even though the expression of some of these upstream molecules was not directly examined by the SOMAScan, IPA was able to predict their activation status (with high significance; i.e., p-values < 10-7 for inhibited and <10-6 for activated). Signal transducer and activator of transcription-3 (STAT3) is a transcriptional regulator that was predicted to be inhibited (Figure 5A) in addition to being under-expressed by ZIKV (Figure 4C). Since it is involved in various key aspects of astrogliosis (promoting formation of reactive astrocytes) and astrocyte migration, it helps explain ZIKV-induced dysregulation leading to the abnormal production of reactive astrocytes, most of which are related to glial scar formation (Okada et al., 2006; Ceyzeriat et al., 2016; Renault-Mihara et al., 2017). In addition, STAT3 also has a role in type 1 interferon signaling pathway; its under-expression can result in modulation of the antiviral activity/response due to IFN type 1 (Ho and Ivashkiv, 2006; Hong and Song, 2014; Suarez et al., 2018). Interferon type I (especially IFNα) activates the STAT3 gene which in turn attenuates type I IFN inflammatory responses (Ho and Ivashkiv, 2006). This STAT3-induced negative regulation of type I IFN response has been shown in mice models where a STAT3 KO resulted in an increased IFN type I production upon viral infection (Wang et al., 2011). This enhanced immune response stimulates expression of a number of ISGs, including ISG15 (Wang et al., 2011), which we found was over-expressed 2.75-fold by 48 hpi (Table 2). Deficiency in IFN-α/β receptor (IFNAR) is also known to make mice highly susceptible to ZIKV infection (Yockey et al., 2018). Since interferon type I activates STAT3, the decrease in STAT3 expression in our data set can in part be explained with the under-expression of IFNA2 (Table 2). To better understand the functions of STAT3, analyses of STAT3 activity and its modulation in response to ZIKV-induced IFNA2 under-expression seems like a next plausible approach. Given recent research showing induction of alpha/beta interferon (IFN-α/β) in response to coronavirus infection in the CNS, which is followed by up-regulation of various interferon stimulating genes (ISGs), the impact of ZIKV infection on the expression of ISGs can in part be explained by modulation in the IFN type 1 responses (Hwang and Bergmann, 2018). The observed increase in ISG15 due to ZIKV infection at 48 hpi in U251 cells is consistent with the findings on ISG15 expression modulation due to other flaviviruses like Dengue virus (DENV) and West Nile Virus (WNV) (Dai et al., 2011). Since ISG15 has such an important immunomodulatory role, its over-expression and antiviral roles have also been linked with other viruses like Hepatitis E virus (HEV), influenza A/B, herpes simplex virus type 1, and HIV (Dai et al., 2011; Sooryanarain et al., 2017). Furthermore, other interferons, such as interferon lambda 1 (IFNL1) were also over-expressed due to ZIKV infection by 48 hpi (Table 2). IFN-λ is already known to inhibit ZIKV infection in the female reproductive tract (Caine et al., 2019). Additional studies focusing on each of these types of interferon I, II, and III, and their interplay with each other and with STAT3 in response to ZIKV infection, could help explain viral-induced host responses.

The 26S proteasome complex was another upstream regulator predicted to be activated by ZIKV infection (Figure 5B). The 26S proteasome complex is necessary for degradation of target proteins tagged by multiple ubiquitin markers, utilizing the ubiquitin conjugation cascade and the proteolytic properties of the proteasome core (Myung et al., 2001; Baert et al., 2007). Thus, it is likely that the proteasome, which would normally use this protein degradation process to remove misfolded proteins from the cellular system, is turned on by ZIKV to aid its replication. Further study to fully understand the relationship between the activated 26S signaling pathway and the production of viral proteins is warranted. The PI3K complex also was predicted to be inhibited by ZIKV. The predicted inhibition of PI3K complex results in dysregulation of various aspects of cell function like cell survival, cell proliferation, cellular trafficking, protein translation and RNA processing, since the PI3K complex lies upstream of most of these cellular functions’ respective molecular pathways (Wang L. L et al., 2017; Ahmad and Siddiqui, 2018). One of its downstream molecules is AKT, RAC-alpha/beta/gamma serine/threonine-protein kinase, and its downregulation by ZIKV is known to modulate glucose metabolism, cell differentiation, regulation of gene transduction and apoptosis (Wang L. L et al., 2017; Ahmad and Siddiqui, 2018). The ERK1/2 (Extracellular signal-regulated kinase 1/2) and MAP2K1 (Dual specificity mitogen-activated protein kinase 1) molecules were another set of molecules predicted to be inhibited by ZIKV (Figure 4, 5). These proteins also are modulated by other flaviviruses (Cheng et al., 2018). For example, HCV modulates the Ras/Raf/MEK pathways to attenuate the JAK-STAT pathway and induce viral replication, YFV and DENV-2 both induce MAPK and ERK1/2 pathways respectively to promote their replication (Cheng et al., 2018). In addition, MEK/ERK signaling is involved in modulating astrocytic differentiation and an inhibition of these molecules is associated with the impairment of cell differentiation and modulation of astrocyte reactivity (Mandell and VandenBerg, 1999; Li et al., 2012). This predicted inhibition also is linked with disruption in neuronal support and synaptic functioning, two crucial functions of normally differentiated astrocytic cells (Mandell and VandenBerg, 1999; Li et al., 2012). Some of these ZIKV-mediated alterations on the infected U-251 cellular proteome also were observed in our recent study of ZIKV-infected Vero cells (Glover et al., 2019). For example, ZIKV down-regulated STAT3 and up-regulated eIF5A expression in both cell systems (Glover et al., 2019). In Vero cells, ZIKV induced activation of the Tec Kinase pathway, linking it to Ca+2 influx, cell proliferation, differentiation, migration and apoptosis, showing that even though the underlying effect ZIKV has on the molecular profile of Vero cells varies from that of U-251 cells, the effect at the functional level has some similarities (Glover et al., 2019). Recognizing these similarities and differences in the protein expression levels in response to ZIKV infection in both U251 and Vero cells are very important and relevant. Since the Vero cell line still is the most commonly used cell line to test a number of antiviral drugs, it is crucial to understand how well the proteomic profile of ZIKV infection in U251 cells, which is a more physiologically relevant cell type, compares to that in the Vero cell type (Perera et al., 2008; Farias et al., 2013; Zamora et al., 2016; Wang S. B et al., 2017; Haviernik et al., 2018). This way, eventual translation of the large number of antiviral drugs and vaccines that have already been tested in Vero cells against a number of other flaviviruses like Dengue Virus Type 2 and West Nile Virus, into the ZIKV system can be made more efficiently (Perera et al., 2008; Farias et al., 2013; Zamora et al., 2016).

Conclusion

In conclusion, ZIKV successfully infects and replicates in astrocytic U-251 cells. A detailed proteomic analysis of the ZIKV-induced dysregulated cellular proteins revealed numerous dysregulated networks and cellular functions, especially those pertaining to astrocytic functions in the brain, and upstream regulator molecules that provide more explanation of how ZIKV can cause neurodegeneration and impairment of neural circuit and network development resulting in conditions like microcephaly. Extension of similar studies to other cell types infectable by ZIKV are warranted to develop a fuller comprehension of ZIKV-induced proteomic alterations that could be targeted for therapeutic intervention.

Statements

Data availability statement

All datasets generated for this study are included in the manuscript and/or the Supplementary Files.

Author contributions

AS performed most of the experiments, interpreted the results, and wrote the manuscript first draft. KG assisted in experimental set-up, in results interpretation, and edited the manuscript. KC provided overall guidance, interpreted the results, and edited the manuscript.

Funding

This research was supported by an operating grant from the Children’s Hospital Research Institute of Manitoba, which was not involved in data generation or interpretation. The authors would like to thank members of their laboratory group for helpful suggestions, and Gao Ang and Victor Spicer for technical support. KG was the recipient of a Graduate Enhancement Training Studentship from the University of Manitoba.

Conflict of interest

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

Supplementary material

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

References

  • 1

    AhmadF.SiddiquiA. (2018). Zika virus deregulate AKT signaling pathway: could be reason of microcephaly.Virol. Immunol. J.2:000183.

  • 2

    AllenN. J. (2013). Role of glia in developmental synapse formation.Curr. Opin. Neurobiol.2310271033. 10.1016/j.conb.2013.06.004

  • 3

    Angers-LoustauA.HeringR.WerbowetskiT. E.KaplanD. R.Del MaestroR. F. (2004). Src regulates actin dynamics and invasion of malignant glial cells in three dimensions.Mol. Cancer Res.2595605.

  • 4

    AshpoleN. M.ChawlaA. R.MartinM. P.BrustovetskyT.BrustovetskyN.HudmonA. (2013). Loss of calcium/calmodulin-dependent protein kinase II activity in cortical astrocytes decreases glutamate uptake and induces neurotoxic release of ATP.J. Biol. Chem.2881459914611. 10.1074/jbc.M113.466235

  • 5

    AshpoleN. M.SongW. H.BrustovetskyT.EnglemanE. A.BrustovetskyN.CumminsT. R.et al (2012). Calcium/calmodulin-dependent protein kinase II (CaMKII) inhibition induces neurotoxicity via dysregulation of glutamate/calcium signaling and hyperexcitability.J. Biol. Chem.28784958506. 10.1074/jbc.M111.323915

  • 6

    AyanlajaA. A.ZhangB.JiG.GaoY.WangJ.KanworeK.et al (2018). The reversible effects of glial cell line-derived neurotrophic factor (GDNF) in the human brain.Semin. Cancer Biol.53212222. 10.1016/j.semcancer.2018.07.005

  • 7

    BaertJ. L.BeaudoinC.MonteD.DegernyC.MauenS.de LaunoitY. (2007). The 26S proteasome system degrades the ERM transcription factor and regulates its transcription-enhancing activity.Oncogene26415424. 10.1038/sj.onc.1209801

  • 8

    BarzonL.TrevisanM.SinigagliaA.LavezzoE.PaluG. (2016). Zika virus: from pathogenesis to disease control.FEMS Microbiol. Lett.363:202. 10.1093/femsle/fnw202

  • 9

    Becerra-CalixtoA.Cardona-GomezG. P. (2017). The role of astrocytes in neuroprotection after brain stroke: potential in cell therapy.Front. Mol. Neurosci.10:88. 10.3389/fnmol.2017.00088

  • 10

    BozoyanL.KhlghatyanJ.SaghatelyanA. (2012). Astrocytes control the development of the migration-promoting vasculature scaffold in the postnatal brain via VEGF signaling.J. Neurosci.3216871704. 10.1523/JNEUROSCI.5531-11.2012

  • 11

    CaineE. A.ScheafferS. M.AroraN.ZaitsevK.ArtyomovM. N.CoyneC. B.et al (2019). Interferon lambda protects the female reproductive tract against Zika virus infection.Nat. Commun.10:280. 10.1038/s41467-018-07993-2

  • 12

    CeyzeriatK.AbjeanL.Carrillo-de SauvageM. A.Ben HaimL.EscartinC. (2016). The complex STATes of astrocyte reactivity: How are they controlled by the JAK-STAT3 pathway?Neuroscience330205218. 10.1016/j.neuroscience.2016.05.043

  • 13

    ChaboubL. S.DeneenB. (2013). Astrocyte form and function in the developing central nervous system.Semin. Pediatr. Neurol.20230235. 10.1016/j.spen.2013.10.003

  • 14

    ChawlaA.JohnsonD.HudmonA. (2014). Mechanisms underlying CaMKII regulation of astrocytic glutamate transporters.Biophys. J.106:104a. 10.1016/j.bpj.2013.11.646

  • 15

    ChengF.da SilvaS. R.HuangI. C.JungJ. U.GaoS. J. (2018). Suppression of Zika virus infection and replication in endothelial cells and astrocytes by PKA inhibitor PKI 14-22.J. Virol.92:e02019-17. 10.1128/JVI.02019-17

  • 16

    ChoiB. R.ChoW. H.KimJ.LeeH. J.ChungC.JeonW. K.et al (2014). Increased expression of the receptor for advanced glycation end products in neurons and astrocytes in a triple transgenic mouse model of Alzheimer’s disease.Exp. Mol. Med.46:e75. 10.1038/emm.2013.147

  • 17

    CoombsK. M. (2011). Quantitative proteomics of complex mixtures.Expert Rev. Proteomics8659677. 10.1586/epr.11.55

  • 18

    CoombsK. M.BerardA.XuW.KrokhinO.MengX.CortensJ. P.et al (2010). Quantitative proteomic analyses of influenza virus-infected cultured human lung cells.J. Virol.841088810906. 10.1128/JVI.00431-10

  • 19

    DaiJ.PanW.WangP. (2011). ISG15 facilitates cellular antiviral response to dengue and west nile virus infection in vitro.Virol J.8:468. 10.1186/1743-422X-8-468

  • 20

    De CarvalhoN. S.De CarvalhoB. F.FugacaC. A.DorisB.BiscaiaE. S. (2016). Zika virus infection during pregnancy and microcephaly occurrence: a review of literature and Brazilian data.Braz. J. Infect Dis.20282289. 10.1016/j.bjid.2016.02.006

  • 21

    DevhareP.MeyerK.SteeleR.RayR. B.RayR. (2017). Zika virus infection dysregulates human neural stem cell growth and inhibits differentiation into neuroprogenitor cells.Cell Death Dis.8:e3106. 10.1038/cddis.2017.517

  • 22

    FariasK. J. S.MachadoP. R. L.da FonsecaB. A. L. (2013). Chloroquine inhibits Dengue virus type 2 replication in Vero cells but not in C6/36 cells.Sci. World J.2013:282734. 10.1155/2013/282734

  • 23

    GabrielE.RamaniA.KarowU.GottardoM.NatarajanK.GooiL. M. (2017). Recent Zika virus isolates induce premature differentiation of neural progenitors in human brain organoids.Cell Stem Cell20397406. 10.1016/j.stem.2016.12.005

  • 24

    GarcezP. P.LoiolaE. C.Da CostaR. M.HigaL. M.TrindadeP.DelvecchioR.et al (2016). Zika virus impairs growth in human neurospheres and brain organoids.Science352816818. 10.1126/science.aaf6116

  • 25

    GloverK. K. M.CoombsK. M. (2017). Zika virus infection: a review of available techniques towards early detection.J. Mol. Biochem.62632.

  • 26

    GloverK. K. M.GaoA.CoombsK. M. (2019). Vero cell proteomic changes induced by Zika virus infection.Proteomics19:e1800309. 10.1002/pmic.201800309

  • 27

    Gonzalez-PerezO.Quinones-HinojosaA. (2012). Astrocytes as neural stem cells in the adult brain.J Stem Cells7181188.

  • 28

    GuJ. J.XuR.LiY. X.ZhangJ. H.WangS. S. (2016). MicroRNA-218 modulates activities of glioma cells by targeting HMGB1.Am. J. Transl. Res.837803790.

  • 29

    HagemannC.AnackerJ.ErnestusR.VinceG. H. (2012). A complete compilation of matrix metalloproteinase expression in human malignant gliomas.World J. Clin. Oncol.36779. 10.5306/wjco.v3.i5.67

  • 30

    HamaH.HaraC.YamaguchiK.MiyawakiA. (2004). PKC signaling mediates global enhancement of excitatory synaptogenesis in neurons triggered by local contact with astrocytes.Neuron41405415. 10.1016/S0896-6273(04)00007-8

  • 31

    HaviernikJ.StefanikM.FojtikovaM.KaliS.TordoN.RudolfI.et al (2018). Arbidol (Umifenovir): a broad-spectrum antiviral drug that inhibits medically important arthropod-borne Flaviviruses.Viruses Basel10:184. 10.3390/v10040184

  • 32

    HoH. H.IvashkivL. B. (2006). Role of STAT3 in type I interferon responses - Negative regulation of STAT1-dependent inflammatory gene activation.J. Biol. Chem.2811411114118. 10.1074/jbc.M511797200

  • 33

    HongS.SongM. R. (2014). STAT3 but not STAT1 is required for astrocyte differentiation.PLoS One9:e86851. 10.1371/journal.pone.0086851

  • 34

    HouS. M.KumarA.XuZ. K.AiroA. M.StryapuninaI.WongC. P.et al (2017). Zika virus hijacks stress granule proteins and modulates the host stress response.J. Virol.91:e00474-17. 10.1128/JVI.00474-17

  • 35

    HwangM.BergmannC. C. (2018). Alpha/beta interferon (IFN-alpha/beta) signaling in astrocytes mediates protection against viral encephalomyelitis and regulates IFN-gamma-dependent responses.J. Virol.92:e01901-17. 10.1128/JVI.01901-17

  • 36

    IoosS.MalletH. P.GoffartI. L.GauthierV.CardosoT.HeridaM. (2014). Current Zika virus epidemiology and recent epidemics.Med. Maladies Infect.44302307. 10.1016/j.medmal.2014.04.008

  • 37

    JacobsenC. T.MillerR. H. (2003). Control of astrocyte migration in the developing cerebral cortex.Dev. Neurosci.25207216. 10.1159/000072269

  • 38

    JiangX.DongX.LiS. H.ZhouY. P.RaynerS.XiaH. M.et al (2018). Proteomic analysis of Zika virus infected primary human fetal neural progenitors suggests a role for Doublecortin in the pathological consequences of infection in the cortex.Front. Microbiol.9:1067. 10.3389/fmicb.2018.01067

  • 39

    KamatP. K.SwarnkarS.RaiS.KumarV.TyagiN. (2014). Astrocyte mediated MMP-9 activation in the synapse dysfunction: an implication in Alzheimer disease.Ther. Targets Neurol. Dis.1:e243.

  • 40

    KimJ.WaldvogelH. J.FaullR. L. M.CurtisM. A.NicholsonL. F. B. (2015). The RAGE receptor and its ligands are highly expressed in astrocytes in a grade-dependant manner in the striatum and subependymal layer in Huntington’s disease.J. Neurochem.134927942. 10.1111/jnc.13178

  • 41

    KobasaD.JonesS. M.ShinyaK.KashJ. C.CoppsJ.EbiharaH.et al (2007). Aberrant innate immune response in lethal infection of macaques with the 1918 influenza virus.Nature445319323. 10.1038/nature05495

  • 42

    KozakR. A.MajerA.BiondiM. J.MedinaS. J.GoneauL. W.SajeshB. V.et al (2017). MicroRNA and mRNA dysregulation in astrocytes infected with Zika virus.Viruses9:E297. 10.3390/v9100297

  • 43

    KuiperJ. W. P.van HorssenR.OerlemansF.PetersW.van DommelenM. M. T.LindertM. M. T.et al (2009). Local ATP generation by brain-type creatine kinase (CK-B) facilitates cell motility.PLoS One4:e5030. 10.1371/journal.pone.0005030

  • 44

    KumarA.JovelJ.Lopez-OrozcoJ.LimontaD.AiroA. M.HouS. M.et al (2018). Human Sertoli cells support high levels of Zika virus replication and persistence.Sci. Rep.8:5477. 10.1038/s41598-018-23899-x

  • 45

    LandstromM.HeldinN. E.BuS. Z.HermanssonA.ItohS.ten DijkeP.et al (2000). Smad7 mediates apoptosis induced by transforming growth factor beta in prostatic carcinoma cells.Curr. Biol.10535538. 10.1016/S0960-9822(00)00470-X

  • 46

    LeonhardS. E.LantS.JacobsB. C.Wilder-SmithA.FerreiraM. L. B.SolomanT.et al (2018). Zika virus infection in the returning traveller: what every neurologist should know.Pract. Neurol.18271277. 10.1136/practneurol-2017-001789

  • 47

    LiX. Y.NewbernJ. M.WuY. H.Morgan-SmithM.ZhongJ.CharronJ.et al (2012). MEK is a key regulator of gliogenesis in the developing brain.Neuron7510351050. 10.1016/j.neuron.2012.08.031

  • 48

    LiuB.NeufeldA. H. (2007). Activation of epidermal growth factor receptors in astrocytes: from development to neural injury.J. Neurosci. Res.8535233529. 10.1002/jnr.21364

  • 49

    LongF.LiuH.HahnC.SumazinP.ZhangM. Q.ZilbersteinA. (2004). Genome-wide prediction and analysis of function-specific transcription factor binding sites.In Silico Biol.4395410.

  • 50

    LupbergerJ.DuongF. H. T.FofanaI.ZonaL.XiaoF.ThumannC.et al (2013). Epidermal growth factor receptor signaling impairs the antiviral activity of interferon-alpha.Hepatology5812251235. 10.1002/hep.26404

  • 51

    MandellJ. W.VandenBergS. R. (1999). ERK/MAP kinase is chronically activated in human reactive astrocytes.Neuroreport1035673572. 10.1097/00001756-199911260-00019

  • 52

    Marquina-SanchezB.Gonzalez-JorgeJ.Hansberg-PastorV.Wegman-OstroskyT.Baranda-AvilaN.Mejia-PerezS.et al (2017). The interplay between intracellular progesterone receptor and PKC plays a key role in migration and invasion of human glioblastoma cells.J. Steroid Biochem. Mol. Biol.172198206. 10.1016/j.jsbmb.2016.10.001

  • 53

    MazarsA.LallemandF.PrunierC.MaraisJ.FerrandN.PessahM.et al (2001). Evidence for a role of the JNK cascade in Smad7-mediated apoptosis.J. Biol. Chem.2763679736803. 10.1074/jbc.M101672200

  • 54

    MehtaR.SoaresC. N.Medialdea-CarreraR.EllulM.da SilvaM. T. T.Rosala-HallasA.et al (2018). The spectrum of neurological disease associated with Zika and chikungunya viruses in adults in Rio de Janeiro, Brazil: A case series.PLoS Neglect. Trop. Dis.12:e0006212. 10.1371/journal.pntd.0006212

  • 55

    MoreiraJ.PeixotoT. M.SiqueiraA. M.LamasC. C. (2017). Sexually acquired Zika virus: a systematic review.Clin. Microbiol. Infect.23296305. 10.1016/j.cmi.2016.12.027

  • 56

    MuffatJ.LiY.OmerA.DurbinA.BoschI.BakiasiG.et al (2018). Human induced pluripotent stem cell-derived glial cells and neural progenitors display divergent responses to Zika and dengue infections.Proc. Natl. Acad. Sci. U.S.A.11571177122. 10.1073/pnas.1719266115

  • 57

    MyungJ.KimK. B.CrewsC. M. (2001). The ubiquitin-proteasome pathway and proteasome inhibitors.Med. Res. Rev.21245273. 10.1002/med.1009

  • 58

    OgierC.BernardA.CholletA. M.Le DiguardherT.HanessianS.ChartonG.et al (2006). Matrix metalloproteinase-2 (MMP-2) regulates astrocyte motility in connection with the actin cytoskeleton and integrins.Glia54272284. 10.1002/glia.20349

  • 59

    OhtakiN.KamitaniW.WatanabeY.HayashiY.YanaiH.IkutaK.et al (2007). Downregulation of an astrocyte-derived inflammatory protein, S100B, reduces vascular inflammatory responses in brains persistently infected with Borna disease virus.J. Virol.8159405948. 10.1128/JVI.02137-06

  • 60

    OkadaS.NakamuraM.KatohH.MiyaoT.ShimazakiT.IshiiK.et al (2006). Conditional ablation of Stat3 or Socs3 discloses a dual role for reactive astrocytes after spinal cord injury.Nat. Med.12829834. 10.1038/nm1425

  • 61

    OlsenM. E.ConnorJ. H. (2017). Hypusination of eIF5A as a target for antiviral therapy.DNA Cell Biol.36198201. 10.1089/dna.2016.3611

  • 62

    ParkJ.KwonK.KimS. H.YiM. H.ZhangE.KongG.et al (2013). Astrocytic phosphorylation of PDK1 on Tyr9 following an excitotoxic lesion in the mouse hippocampus.Brain Res.15333743. 10.1016/j.brainres.2013.08.028

  • 63

    PelechanoV.AlepuzP. (2017). eIF5A facilitates translation termination globally and promotes the elongation of many non polyproline-specific tripeptide sequences.Nucleic Acids Res.4573267338. 10.1093/nar/gkx479

  • 64

    PengH. P.LiuB.YvesT. D.HeY. H.WangS. J.TangH. L.et al (2018). Zika virus induces autophagy in human umbilical vein endothelial cells.Viruses10:E259. 10.3390/v10050259

  • 65

    PereraR.KhaliqM.KuhnR. J. (2008). Closing the door on flaviviruses: entry as a target for antiviral drug design.Antiviral Res.801122. 10.1016/j.antiviral.2008.05.004

  • 66

    ReemstK.NoctorS. C.LucassenP. J.HolE. M. (2016). The indispensable roles of microglia and astrocytes during brain development.Front. Hum. Neurosci.10:566. 10.3389/fnhum.2016.00566

  • 67

    Renault-MiharaF.MukainoM.ShinozakiM.KumamaruH.KawaseS.BaudouxM.et al (2017). Regulation of RhoA by STAT3 coordinates glial scar formation.J. Biol. Cell2162533255010.1083/jcb.201610102

  • 68

    RetallackH.Di LulloE.AriasC.KnoppK. A.LaurieM. T.Sandoval-EspinosaC.et al (2016). Zika virus cell tropism in the developing human brain and inhibition by azithromycin.Proc. Natl. Acad. Sci. U.S.A.1131440814413. 10.1073/pnas.1618029113

  • 69

    ScaturroP.StukalovA.HaasD. A.CorteseM.DraganovK.PlaszczycaA.et al (2018). An orthogonal proteomic survey uncovers novel Zika virus host factors.Nature561253257. 10.1038/s41586-018-0484-5

  • 70

    ScheidM. P.ParsonsM.WoodgettJ. R. (2005). Phosphoinositide-dependent phosphorylation of PDK1 regulates nuclear translocation.Mol. Cell. Biol.2523472363. 10.1128/MCB.25.6.2347-2363.2005

  • 71

    SchullerA. P.WuC. C. C.DeverT. E.BuskirkA. R.GreenR. (2017). eIF5A functions globally in translation elongation and termination.Mol. Cell66194205. 10.1016/j.molcel.2017.03.003

  • 72

    SiemannD. N.StrangeD. P.MaharajP. N.ShiP. Y.VermaS. (2017). Zika virus infects human Sertoli cells and modulates the integrity of the in vitro blood-testis barrier model.J. Virol.91:e00623-17. 10.1128/JVI.00623-17

  • 73

    SmirnovaO. A.BartoschB.ZakirovaN. E.KochetkovS. N.IvanovA. V. (2018). Polyamine metabolism and oxidative protein folding in the ER as ROS-producing systems neglected in virology.Int. J. Mol. Sci.19:1219. 10.3390/ijms19041219

  • 74

    SOMAScan (2017). Proteomic Assay Technical White Paper, SSM-002 Rev. Available at: http://somalogic.com/wp-content/uploads/2017/06/SSM-002-Technical-White-Paper_010916_LSM1.pdf.

  • 75

    SooryanarainH.RogersA. J.CaoD. J.HaacM. E. R.KarpeY. A.MengX. J. (2017). ISG15 modulates type I interferon signaling and the antiviral response during hepatitis E virus replication.J. Virol.91:e00621-17. 10.1128/JVI.00621-17

  • 76

    SteenC.WilczakN.HoogduinJ. M.KochM.De KeyserJ. (2010). Reduced creatine kinase B activity in multiple sclerosis normal appearing white matter.PLoS One5:e10811. 10.1371/journal.pone.0010811

  • 77

    SuarezA. A. R.Van RenneN.BaumertT. F.LupbergerJ. (2018). Viral manipulation of STAT3: evade, exploit, and injure.PLoS Pathog.14:e1006839. 10.1371/journal.ppat.1006839

  • 78

    TakashimaS.Kanatsu-ShinoharaM.TanakaT.MorimotoH.InoueK.OgonukiN.et al (2015). Functional differences between GDNF-dependent and FGF2-dependent mouse spermatogonial stem cell self-renewal.Stem Cell Rep.4489502. 10.1016/j.stemcr.2015.01.010

  • 79

    TeixeiraM. G.CostaM. D. N.de OliveiraW. K.NunesM. L.RodriguesL. C. (2016). The epidemic of Zika virus-related microcephaly in Brazil: detection, control, etiology, and future scenarios.Am. J. Public Health106601605. 10.2105/AJPH.2016.303113

  • 80

    TerniB.Lopez-MurciaF. J.LlobetA. (2017). Role of neuron-glia interactions in developmental synapse elimination.Brain Res. Bull.1297481. 10.1016/j.brainresbull.2016.08.017

  • 81

    TyzackG. E.SitnikovS.BarsonD.Adams-CarrK. L.LauN. K.KwokJ. C.et al (2014). Astrocyte response to motor neuron injury promotes structural synaptic plasticity via STAT3-regulated TSP-1 expression.Nat. Commun.5:4294. 10.1038/ncomms5294

  • 82

    van den PolA. N.MaoG.YangY.OrmaghiS.DavisJ. N. (2017). Zika virus targeting in the developing brain.J. Neurosci.3721612175. 10.1523/JNEUROSCI.3124-16.2017

  • 83

    VesterD.RappE.GadeD.GenzelY.ReichlU. (2009). Quantitative analysis of cellular proteome alterations in human influenza A virus-infected mammalian cell lines.Proteomics933163327. 10.1002/pmic.200800893

  • 84

    WagnerB.NatarajanA.GrunaugS.KroismayrR.WagnerE. F.SibiliaM. (2006). Neuronal survival depends on EGFR signaling in cortical but not midbrain astrocytes.EMBO J.25752762. 10.1038/sj.emboj.7600988

  • 85

    WangL. L.YangL.FikrigE.WangP. H. (2017). An essential role of PI3K in the control of West Nile virus infection.Sci. Rep.7:3724. 10.1038/s41598-017-03912-5

  • 86

    WangS. B.LiuY.GuoJ.WangP. L.ZhangL. K.XiaoG. F.et al (2017). Screening of FDA-approved drugs for inhibitors of Japanese Encephalitis virus infection.J. Virol.91:e01055-17. 10.1128/JVI.01055-17

  • 87

    WangW. B.LevyD. E.LeeC. K. (2011). STAT3 negatively regulates type I IFN-mediated antiviral response.J. Immunol.18725782585. 10.4049/jimmunol.1004128

  • 88

    WangZ.LiD. D.LiangY. Y.WangD. S.CaiN. S. (2002). Activation of astrocytes by advanced glycation end products: cytokines induction and nitric oxide release.Acta Pharmacol. Sin.23974980.

  • 89

    WenZ. X.SongH. J.MingG. L. (2017). How does Zika virus cause microcephaly?Genes Dev.31849861. 10.1101/gad.298216.117

  • 90

    WishartT. M.ParsonS. H.GillingwaterT. H. (2006). Synaptic vulnerability in neurodegenerative disease.J. Neuropathol. Exp. Neurol.65733739. 10.1097/01.jnen.0000228202.35163.c4

  • 91

    WuD. M.WangS.WenX.HanX. R.WangY. J.FanS. H.et al (2018). MircoRNA-1275 promotes proliferation, invasion and migration of glioma cells via SERPINE1.J. Cell. Mol. Med.2249634974. 10.1111/jcmm.13760

  • 92

    XinQ. L.DengC. L.ChenX.WangJ.WangS. B.WangW.et al (2017). Quantitative proteomic analysis of mosquito C6/36 cells reveals host proteins involved in Zika virus infection.J. Virol.91:e00554-17. 10.1128/JVI.00554-17

  • 93

    YangS. Y.HeX. Y.IsaacsC.DobkinC.MillerD.PhilippM. (2014). Roles of 17 beta-hydroxysteroid dehydrogenase type 10 in neurodegenerative disorders.J. Steroid Biochem. Mol. Biol.143460472. 10.1016/j.jsbmb.2014.07.001

  • 94

    YangS. Y.HeX. Y.OlpinS. E.SuttonV. R.McMenaminJ.PhilippM.et al (2009). Mental retardation linked to mutations in the HSD17B10 gene interfering with neurosteroid and isoleucine metabolism.Proc. Natl. Acad. Sci. U.S.A.1061482014824. 10.1073/pnas.0902377106

  • 95

    YangY.XinX. K.FuX.XuD. M. (2018). Expression pattern of human SERPINE2 in a variety of human tumors.Oncol. Lett.1545234530. 10.3892/ol.2018.7819

  • 96

    YatesJ. R.RuseC. I.NakorchevskyA. (2009). Proteomics by mass spectrometry: approaches, advances, and applications.Ann. Rev. Biomed. Eng.114979. 10.1146/annurev-bioeng-061008-124934

  • 97

    YockeyL. J.JuradoK. A.AroraN.MilletA.RakibT.MilanoK. M.et al (2018). Type I interferons instigate fetal demise after Zika virus infection.Sci. Immunol.3:eaao1680

  • 98

    ZamoraA. P.EdmondsJ. H.ReynoldsM. J.KhromykhA. A.RalphS. J. (2016). The in vitro and in vivo antiviral properties of combined monoterpene alcohols against West Nile virus infection.Virology4951832. 10.1016/j.virol.2016.04.021

  • 99

    ZhangZ. W.LiZ. L.YuanS. (2017). The role of secretory autophagy in Zika virus transfer through the placental barrier.Front. Cell Infect. Microbiol.6:206. 10.3389/fcimb.2016.00206

  • 100

    ZhengK.KitazatoK.WangY. F. (2014). Viruses exploit the function of epidermal growth factor receptor.Rev. Med. Virol.24274286. 10.1002/rmv.1796

Summary

Keywords

RNA virus infection, proteomics, aptamers, SOMAScan, astrocytic dysregulation, virus-host interaction

Citation

Sher AA, Glover KKM and Coombs KM (2019) Zika Virus Infection Disrupts Astrocytic Proteins Involved in Synapse Control and Axon Guidance. Front. Microbiol. 10:596. doi: 10.3389/fmicb.2019.00596

Received

08 January 2019

Accepted

08 March 2019

Published

26 March 2019

Volume

10 - 2019

Edited by

Juliano Bordignon, Fundação Oswaldo Cruz (Fiocruz), Brazil

Reviewed by

Florencia Meyer, Mississippi State University, United States; Yhojan Rodríguez, Universidad del Rosario, Colombia

Updates

Copyright

*Correspondence: Kevin M. Coombs, ;

This article was submitted to Virology, a section of the journal Frontiers in Microbiology

Disclaimer

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

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics