Next Article in Journal
The Effects of Genetic Variation on H7N9 Avian Influenza Virus Pathogenicity
Previous Article in Journal
Metagenomic Approach with the NetoVIR Enrichment Protocol Reveals Virus Diversity within Ethiopian Honey Bees (Apis mellifera simensis)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 12
Issue 11
10.3390/v12111219
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Zika Virus Infection on Human Neural Stem Cell MicroRNA Signatures

by
Denna Tabari
1,
Catharina Scholl
1,*,
Michael Steffens
1,
Sandra Weickhardt
1,
Fabian Elgner
2,
Daniela Bender
2,
Marie-Luise Herrlein
2,
Catarina Sabino
2,
Vesselina Semkova
3,
Michael Peitz
3,4,
Andreas Till
3,
Oliver Brüstle
3,
Eberhard Hildt
2 and
Julia Stingl
5
1
Research Division, Federal Institute for Drugs and Medical Devices, 53175 Bonn, Germany
2
Department of Virology, Paul-Ehrlich-Institut, 63225 Langen, Germany
3
Institute of Reconstructive Neurobiology, LIFE & BRAIN Center, University of Bonn Medical Faculty & University Hospital Bonn, 53127 Bonn, Germany
4
Cell Programming Core Facility, Medical Faculty, University of Bonn, 53172 Bonn, Germany
5
Department of Clinical Pharmacology, University Hospital, RWTH Aachen University, 52074 Aachen, Germany
*
Author to whom correspondence should be addressed.
Viruses 2020, 12(11), 1219; https://doi.org/10.3390/v12111219
Submission received: 9 September 2020 / Revised: 9 October 2020 / Accepted: 23 October 2020 / Published: 27 October 2020
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Zika virus (ZIKV) is a mosquito-borne virus, which can cause brain abnormalities in newborns, including microcephaly. MicroRNAs (miRNAs) are small non-coding RNAs, which post- transcriptionally regulate gene expression. They are involved in various processes including neurological development and host responses to viral infection, but their potential role in ZIKV pathogenesis remains poorly understood. MiRNAs can be incorporated into extracellular vesicles (EVs) and mediate cell-to-cell communication. While it is well known that in viral infections EVs carrying miRNAs can play a crucial role in disease pathogenesis, ZIKV effects on EV-delivered miRNAs and their contribution to ZIKV pathogenesis have not been elucidated. In the present study, we profiled intracellular and EV-derived miRNAs by next generation sequencing and analyzed the host mRNA transcriptome of neural stem cells during infection with ZIKV Uganda and French Polynesia strains. We identified numerous miRNAs, including miR-4792, which were dysregulated at the intracellular level and had altered levels in EVs during ZIKV infection. Integrated analyses of differentially expressed genes and miRNAs showed that ZIKV infection had an impact on processes associated with neurodevelopment and oxidative stress. Our results provide insights into the roles of intracellular and EV-associated host miRNAs in ZIKV pathogenesis.

1. Introduction

Zika virus (ZIKV), a mosquito-borne virus and member of the Flaviviridae family, was initially isolated in 1947 from a rhesus monkey in the Zika forest in Uganda [1]. In 2007 the first human epidemic was recorded on Yap Island, followed by a major outbreak in French Polynesia in 2013 [2,3]. In the course of this outbreak, ZIKV infection was for the first time linked to Guillain-Barré Syndrome and severe neurodevelopmental defects, including microcephaly [4,5,6,7]. International awareness of the potential for severe disease outcomes grew in 2015 when ZIKV emerged in Brazil and subsequently spread to most of the Americas in 2016 [8,9,10]. During the global spread, the virus diversified and presently ZIKV strains are classified into African and Asian/American lineages. ZIKV has circulated for decades in Africa without any documented clinical relevance. It is unclear whether genetic change might have resulted in a greater capacity of neurovirulence or whether ZIKV has always been teratogenic and severe cases have not been documented before 2013 [11,12]. Due to the rapid spread of the virus and its severe effects on infants, on 1st February 2016, the World Health Organization (WHO) declared a Public Health Emergency of International Concern (PHEIC) [13]. Since then, significant progress has been made in elucidating the effects of ZIKV on the host transcriptome; however, the molecular mechanisms by which ZIKV leads to microcephaly are still not deeply understood. Previous studies have shown that ZIKV can infect human neural stem cells (hNSCs) and human neural progenitor cells (hNPCs), thus disrupting pathways involved in differentiation, apoptosis and cell cycle control [14,15,16,17].
MicroRNAs (miRNAs) are a class of small non-coding RNA molecules (containing about 22 nucleotides) that act as post-transcriptional regulators of mRNA. By binding to partially complementary sites of the target mRNA, miRNAs can suppress the translation [18]. MiRNAs play an important role in the regulation of genes involved in various processes including development, differentiation and apoptosis [19,20,21]. MiRNAs can further be incorporated into extracellular vesicles (EVs) and serve as signaling molecules to mediate intercellular communication [22]. Various Flaviviridae viruses modulate intracellular and EV-associated miRNA profiles, which can have an impact on the pathogenesis [23,24]. In the context of ZIKV infection, the relevance of EV-derived miRNAs has remained unexplored. Since only a few studies have analyzed the effects of ZIKV infection on global miRNA regulation at the intracellular level, their role in ZIKV pathogenesis is not fully understood [25,26]. Kozak et al. have found that ZIKV infection alters the host transcriptomic profiles in astrocytes and dysregulates miRNAs involved in the unfolded protein response pathway, which is a cellular response mechanism to endoplasmic reticulum (ER) stress [26].
Accumulating evidences have linked ER stress and oxidative stress to the pathogenesis of Flaviviridae viruses, including Dengue virus (DENV) and hepatitis C virus (HCV) [27]. An essential factor in the regulation of oxidative stress is the transcription factor Nrf2 (nuclear factor erythroid 2 (NF-E2)-related factor 2), which triggers the expression of numerous cytoprotective genes [28,29]. Under unstressed conditions, Nrf2 is associated with the actin binding protein Keap1 (Kelch-like ECH associating protein 1) and localized within the cytoplasm. Keap1 targets Nrf2 for rapid proteasomal degradation though ubiquitination, thereby maintaining it at relatively low levels in the cell [30]. Upon oxidative stress, Nrf2 dissociates from Keap1 and translocates into the nucleus, where it forms a heterodimer with small Maf (sMaf) proteins and subsequently binds to the antioxidant response element (ARE). Nrf2/ARE-regulated genes play a crucial role in cellular adaptation to reactive oxygen species (ROS) and some viruses manipulate the antioxidant system for their own benefit [28,31,32]. Little is currently known about the role of oxidative stress and the underlying molecular details during ZIKV infection in neuronal tissue on neurodevelopmental effects [33].
In this study, we investigated the role of miRNA-mediated gene expression during ZIKV infection of induced pluripotent stem cell (iPSC)-derived hNSCs. To gain more knowledge on the mechanisms of pathogenesis of neurodevelopmental damage, we profiled host miRNA and mRNA transcriptomes and performed integrative analyses.

2. Materials and Methods

2.1. Cell Culture

Long-term self-renewing neuroepithelial-like stem cells (lt-NES®) represent a reductionist in vitro model for the early stages of human neural stem cell differentiation [34]. These cells, which are derived from both human embryonal stem cells (hESCs) [34] and induced pluripotent stem cells (iPSCs) [35], retain a constant neuro- and gliogenic potential even after long-term proliferation. Lt-NES® cells used in the present study were generated from iPSC line iLB-C-31f [36], which is derived from a female donor in her mid-20s by retroviral reprogramming. Routine characterization of the resulting cell line was performed as depicted in the supplemental information (Figures S1 and S2). The cells were maintained on plastic 6-well cell culture dishes coated with 1x poly-L-ornithine and 10 mg/mL laminin (Merck, Darmstadt, Germany) in Dulbecco’s modified Eagle’s/F12 medium (DMEM-F12) (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with N2 Supplement (PAA Laboratories, Pasching, Austria) and containing 10 ng/mL basic fibroblast growth factor 2 (FGF2), 10 ng/mL epidermal growth factor (EGF) (both from R&D Systems, Minneapolis, MN, USA) and 1:1000 B-27 Supplement (Thermo Fisher Scientific, Waltham, MA, USA). Passaging was performed every 2–3 days by splitting at a 1:2 ratio. In short, cells were detached using incubation with 0.125% trypsin (Thermo Fisher Scientific, MA, USA) for 5–10 min at room temperature. Trypsin was deactivated with trypsin inhibitor (Thermo Fisher Scientific, MA, USA) and the cells were spun down at 300× g for 3 min. Supernatant was aspirated and the cell pellet was resuspended in maintenance medium and plated onto poly-L-ornithine/laminin pre-coated cell culture dishes.

2.2. Infection

Two ZIKV strains were used in the study: ZIKV Uganda 976 (kindly provided by the European Virus Archive), for which no associations with severe symptoms have been documented, and ZIKV PF13/251013-18 from French Polynesia (kindly provided by Professor Musso from the Institute Louis Marladé in Tahiti and by the European Virus Archive), which has been linked to severe pathological effects [2]. The lt-NES® cells were seeded into 6-well plates (1 million cells per well) and infected with a multiplicity of infection (MOI) of 0.1 or 1 using serum free virus stocks. MOI of 1 was chosen with the aim of seeing more significant effects. The medium of the infected and uninfected cells was changed daily (24, 48 and 72 h post infection) and the supernatants were collected 8 h after medium change for further analyses by plaque assay. Serum free virus stocks were generated by infection of Vero cells with the respective virus strain in Dulbecco’s Modified Eagle’s Medium (DMEM) high glucose (BioWest, Nuaillé, France) supplemented with 10% fetal bovine serum superior (Biochrom GmbH, Berlin, Germany), 2 mM L-glutamine (Biochrom GmbH, Berlin, Germany), 100 U/mL penicillin and 100 μg/mL streptomycin (Paul-Ehrlich-Institut facilities, Langen, Germany). On day 3 post infection, cells were washed once with phosphate-buffered saline (PBS) followed by a medium change to DMEM supplemented with 10% KnockOutTM serum replacement (Thermo Fisher Scientific, Waltham, MA, USA). On day 7 post infection the supernatant was harvested, centrifuged for 5 min at 1000× g and stored at −80 °C until determination of titer by plaque assay and infection of lt-NES® cells.

2.3. Plaque-Forming-Assay

Plaque assay was performed to quantify infectious viral particles. Vero cells were seeded in 6-well plates (3 × 105 per well) and infected with 100 µl of serial dilutions from cell culture supernatants. After 2 h, the medium was removed and the cells were layered with 0.4% agarose solution (SeaPlaque® Agarose; Lonza, Basel, Switzerland). Five days after infection, the agarose overlay was removed and cells were fixed with 4% formaldehyde in PBS for 20 min. Visualization of the plaques was achieved by staining with crystal violet (0.1% in 20% ethanol; Merck, Darmstadt, Germany) 15 min at room temperature. After washing with distilled water to remove excess crystal violet solution, the plates were dried, and plaques were counted.

2.4. MiRNAs: Sequencing and Data Analysis

For the analysis of intracellular miRNAs, total RNA was extracted from infected and uninfected lt-NES® cells at 24, 48 and 72 hpi using miRNeasy Mini Kit (217004, Qiagen, Hilden, Germany) according to the manufacturer’s instructions. For the analysis of EV-derived miRNAs, cell culture supernatants were collected at the indicated time points and EV-associated miRNAs were extracted using exoRNeasy Serum/Plasma Kit (77064, Qiagen, Hilden, Germany) as described by the manufacturer. RNA sequencing libraries were generated using the NEBNext Mulitplex Small RNA Library Prep Set for Illumina (E7300S, New England Biolabs, Frankfurt, Germany) following size selection by running the samples on a 10% polyacrylamide gel and isolating the 147 bp. Samples were purified using the QIAquick PCR Purification Kit (28104, Qiagen, Hilden, Germany). Libraries were validated on a high sensitivity DNA chip using the 2100 Bioanalyzer (5067-4626, Agilent, Santa Clara, CA, USA). Next generation sequencing was conducted on the Illumina® platform using MiSeq Reagent Kit v2 (50 cycles) (MS-102-2001, Illumina, CA, USA). To analyze the miRNA sequencing data, fastq files were generated and adapter trimming was performed with the program cutadapt. Reads were aligned to the whole human genome using bowtie2. The package DESeq2 was used to calculate the fold change, p-value (Wald test) and adjusted p-value (Benjamini Hochberg Adjustment) of differentially expressed miRNAs. To focus on miRNAs that potentially have a greater impact on host transcripts, miRNAs with copy numbers lower than 100 were removed from further data analysis [37]. Gene ontology (GO) analyses of biological processes and pathways of KEGG (Kyoto Encyclopedia of Genes and Genomes), which is a knowledge database for systematic analysis of gene functions, were performed using DIANA mirpath v3 [38,39]. Putative targets were determined using the homology search algorithm microT-CDS and TarBase, which is a database of experimentally validated miRNA-gene interactions (p-value threshold of 0.05) [40]. For the microT-CDS prediction, a microT threshold of 0.8 was set. The miRNA sequencing data sets have been deposited in Gene Expression Omnibus (GEO) database (GSE157532).

2.5. Microarray-Based Gene Expression Analysis

Gene expression profiling was performed using SurePrint G3 Human Gene Expression 8x60K Microarray Kit (Agilent One Color Microarray Technology, Agilent Technologies, Santa Clara, CA, USA). Of the total RNA, 100 ng was extracted using miRNeasy Mini Kit (217004, Qiagen, Hilden, Germany) according to the manufacturer’s instructions and was used for amplification.
Microarray results were extracted using Agilent Feature Extraction version 11.0.1.1 and analyzed using Genespring 14.9.1 GX software (Agilent Technologies, Santa Clara, CA, USA). To reduce artefacts, signal intensities below 50 were removed from the dataset. A two-way ANOVA test was used as statistical analysis and p-values were adjusted for multiple test correction using the Benjamini Hochberg algorithm. Genespring 14.9.1 GX software was further used for fold change analyses and to crosscheck the list of differentially expressed genes with putative miRNA targets obtained from bioinformatics databases (PiCTar, TarBase, TargetScan, miRDB) [41,42,43,44]. Pathways and network patterns associated with significantly dysregulated genes were analyzed using Reactome version 73 (https://reactome.org) [45]. The microarray data sets have been deposited in GEO (GSE157532).

2.6. Real-Time Quantitative PCR

Real-time quantitative PCR (qRT-PCR) was performed to validate selected genes and miRNAs. According to the manufacturer’s instructions, cDNA was generated from 500 ng of total RNA for the validation of miRNAs using miScript II RT Kit (218161, Qiagen, Hilden, Germany) and from 1000 ng of total RNA for the validation of mRNA using QuantiTect Rev. Transcription Kit (205311, Qiagen, Hilden, Germany). Quantitative PCR was performed with the SYBR Green PCR Kit (204145, Qiagen, Hilden, Germany) using a Roche LightCycler 480 (Roche, Basel, Switzerland). For miRNA expression analyses, miScript Primer Assays were used (Hs_miR-205_1; Hs-miR-182_2; Hs_miR-4792-1; Hs-miR-16_2, Qiagen, Hilden, Germany) and Ct values were normalized to miR-16 levels. Gene expression was analyzed using Quantitect Primer Assays (Hs_ASNS_1_SG; Hs_FOXC1_2_SG; Hs_SESN2_1_SG; Hs_GAPDGH_1_SG, Qiagen, Hilden, Germany) and Ct values were normalized to GAPDH. Fold changes were calculated with the ∆∆Ct method [46].

2.7. Oxyblot

Carbonylation of proteins by ROS was analyzed using the OxyBlotTM protein oxidation detection kit (Merck, Darmstadt, Germany) according to the manufacturer’s instructions. Briefly, ZIKV-infected cells and uninfected cells were harvested after 24 h, 48 h and 72 h, washed with PBS and lysed with RIPA buffer. For each sample the same amount of protein (15 µg) was used. Band intensity values were obtained using Image Studio Lite Software Version 5.2 (LI-COR Biosciences, Lincoln, NE, USA).

2.8. Immunofluorescence Microscopy

The cells were seeded on cover slides and fixed with 4% formaldehyde in PBS for 20 min at room temperature. Permeabilization and blocking were performed as described by Elgner et al. [47]. The cells were incubated with anti-ZIKV NS1 (1:1000; Biofront Technologies, Tallahassee, FL, USA) and anti-MafG (1.300; Abcam, Cambridge, UK) and afterwards with the secondary antibodies anti-mouse IgG-Alexa 488 (1:1000; Thermo Fisher Scientific, MA, USA) and anti-rabbit-Cy3 (1:400; Jackson Immunoresearch, West Grove, PA, USA). Nuclei were stained using 4′,6-diamidino-2-phenylindole (DAPI) (Carl Roth, Karlsruhe, Germany). Finally, the stained coverslips were mounted with Mowiol (Merck, Darmstadt, Germany) and analyzed using the confocal laser scanning microscope LSM 510 Meta and ZEN 2009 software (Carl Zeiss, Oberkochen, Germany).

2.9. Statistical Analysis

Statistically significant differences between group means were assessed by 2-tailed, unpaired Student’s t test and visualized using GraphPad Prism 8 software. Differentially expressed miRNAs from sequencing data were analyzed using DESeq2. A two-way ANOVA test was used as a statistic test for the microarray gene expression data. An adjusted p-value ≤ 0.05 was considered statistically significant.

3. Results

3.1. ZIKV Uganda and French Polynesia Alter Gene Expression in lt-NES® Cells

To investigate the impact of ZIKV infection on gene expression and host miRNA regulation in a well characterized model for human neural stem cells, we analyzed the mRNA transcriptome of infected lt-NES® cells at 24, 48 and 72 hpi. To this aim, lt-NES® cells were infected using a MOI of 0.1 with either the French Polynesia PF13/251013-18 (ZIKV Polynesia) or the Uganda 976 strain. The Uganda strain in general had a higher impact on gene expression in lt-NES® cells (Figure 1). By more than two-fold, ZIKV Uganda significantly upregulated 1280 genes and downregulated 1208 genes at 72 hpi, whereas only 258 genes were upregulated and 12 genes were downregulated by ZIKV Polynesia at the same time point. Plaque-forming assays revealed a significantly higher number of infectious particles produced by Uganda than by Polynesia strains at 24 and 48 hpi and no significant difference at 72 hpi (Figure 2).
To identify in which processes the differentially expressed genes were involved, we performed pathway analyses using the reactome pathway database. The top ranked terms for ZIKV Uganda and Polynesia primarily included antiviral response related pathways (Figure 3). Among the most affected processes by ZIKV Uganda, cell cycle related pathways including G1/S-specific transcription were also significantly enriched. Both strains induced upregulation of genes, which are mainly involved in processes related to antiviral interferon signaling, ER stress and the unfolded protein response pathway, which is consistent with previous reports [25,26]. Downregulated genes were largely involved in signal transduction, cell cycle and gene expression related pathways (Figure 4, Tables S1 and S2). Using the comparative toxicogenomics database [48], which provides information about gene–disease relationships, we screened all differentially expressed genes for association with microcephaly. Although for the Uganda strain no neuropathological symptoms are documented, the screening revealed five genes related to microcephaly, which were only dysregulated by ZIKV Uganda and one gene, namely ASNS (asparagine synthetase), which was upregulated by both strains (Table 1). RT-qPCR analyses showed a higher upregulation of ASNS by ZIKV Uganda than by the Polynesia strain (Figure S3).
Overall, our results indicate an increased expression of immune response related genes over time (Table 2 and Table 3) including genes that are known to act as restriction factors for other viruses, namely BST2 (tetherin, Figure S4) and IFIT2 (interferon-induced protein with tetratricopeptide repeats 2) [49,50].
Taken together, our results show that ZIKV Uganda and Polynesia modulate the gene expression pattern in lt-NES® cells by causing a dysregulation of genes mainly involved in interferon signaling and cell cycle, which may contribute to ZIKV pathogenesis.

3.2. Dysregulated miRNAs Play a Role in Neurodevelopment and Cell Cycle

Since miRNAs regulate posttranscriptional gene expression and can play an important role in pathological processes during viral infections, we profiled differential miRNA regulation upon ZIKV infection. Similar to our gene expression analysis, we analyzed miRNA profiles of lt-NES® cells infected with ZIKV Uganda and Polynesia at 24, 48 and 72 hpi at MOI of 0.1. Both virus strains induced dysregulation of miRNAs, whereby overall more miRNAs were upregulated than downregulated (Figure 5a). Interestingly, we found a large overlap between the ZIKV strains, both dysregulating the same 70 miRNAs (Figure 5b, Table S3), with 69 miRNAs upregulated and only one miRNA downregulated by both strains, namely miR-182-5p. To further characterize the differentially regulated miRNAs, potential target genes were predicted using the microT-CDS (v 5.0) database and gene ontology analysis was performed with a focus on biological processes. MiRNAs dysregulated by both strains were involved in the neurotrophin TRK receptor signaling pathway, synaptic transmission, and mitotic cell cycle, among others (Figure 6a). These processes were also enriched for the total of 30 miRNAs, which were dysregulated only by ZIKV Uganda (Table S4, Figure S9).
One specific miRNA—miR-7704—which was upregulated by ZIKV Polynesia and by ZIKV Uganda, showed up with potential roles in developmental processes including neural tube formation (Figure 6b). To further identify miRNA–mRNA interactions during ZIKV infection in lt-NES® cells, we crosschecked the list of potential miRNA targets with mRNAs that were differentially expressed after ZIKV infection of lt-NES® cells. This integrative analysis revealed various genes downregulated by ZIKV Uganda with roles in neurodevelopment, such as EPHB4 (Ephrin Receptor B4) and CNP (2′,3′-Cyclic Nucleotide 3′ Phosphodiesterase),which may be depressed by numerous miRNAs that were upregulated in lt-NES® cells. Thus, multiple upregulated miRNAs could potentially affect the same mRNA targets. Two miRNAs were significantly downregulated only by ZIKV Polynesia, namely miR-433-5p and miR-205-5p (Figure S5). Interestingly, miR-205-5p potentially targets the upregulated gene ASNS (Asparagine Synthetase) [51]. However, further experiments should be performed to assess the interaction between miR-205-5p and ASNS in lt-NES® cells.
Collectively, our findings indicate an alteration of miRNA regulation following ZIKV infection in lt-NES® cells, which was broadly similar between ZIKV Uganda and ZIKV Polynesia. We detected dysregulated miRNAs that might play a role in neurodevelopment and cell cycle related pathways and could be involved in ZIKV pathogenesis.

3.3. ZIKV Modifies the Host miRNA Cargo in Secreted EVs

Since exosome-derived miRNAs are crucial for cell–cell communication and can be a central player in disease pathogenesis during viral infection, our next aim was to study how ZIKV infection affects the level of host miRNAs incorporated into EVs. To this aim, we analyzed EV-derived miRNAs 24, 48 and 72 h post infection of lt-NES® cells at MOI of 1. Uninfected cells kept under the same conditions were used as a control. We found altered levels of EV-associated miRNAs by ZIKV Uganda, increasing over time. In contrast, for ZIKV Polynesia the effect on the signature of EV-derived miRNAs was higher at 48 hpi (Figure 7a). We detected five miRNAs within EVs with levels altered by both virus strains (Figure 7b). Using Diana mirPath v3 we performed KEGG pathway analyses of potential targets for EV-associated miRNAs with altered levels upon infection (listed in Tables S5 and S6), which revealed an enrichment of signaling pathways involved in neurodevelopmental processes including Wnt signaling, Pi3K-Akt signaling and signaling pathways regulating pluripotency of stem cell, among others (Figure 8). Both strains modulated the level of EV-associated miRNAs, which were involved in FoxO signaling and axon guidance. Interestingly, EV-derived miR-4792 was the only miRNA for which we found both an increased level following infection with ZIKV Polynesia (by 1.4-fold) and ZIKV Uganda (by 2.9-fold) at 48 hpi and additionally an upregulation by both strains at the intracellular level at any time point (Figure S6), which we validated by RT-qPCR. Using Diana mirPath v3 we further performed a gene ontology analysis of potential miR-4792 targets predicted by Tarbase, which uses experimentally validated data. Potential miR-4792 targets were involved in processes including in utero embryonic development and negative regulation of oxidoreductase activity, respectively (Table 4). Because of previously described interactions between miR-4792 and FOXC1 (Forkhead Box C1) in various cell lines [20,52] we analyzed the expression of FOXC1 mRNA in lt-NES® cells during ZIKV infection. RT-qPCR analyses revealed that FOXC1 was slightly downregulated by both ZIKV strains (Figure S7). However, further studies should be performed to investigate if FOXC1 is repressed due to the upregulation of miR-4792 in lt-NES® cells upon ZIKV infection.
Taken together, we observed that ZIKV infection induces multiple changes in secretion of host miRNAs in EVs. Overall, these EV-associated miRNAs, particularly miR-4792, were involved in oxidative and neurodevelopmental processes, suggesting that they could be a central player in ZIKV-induced neuropathogenesis given the important role of EVs in intercellular communication.

3.4. Dysregulated miRNAs and mRNAs Involved in Oxidative Stress

Since previous studies have linked ER stress and ROS increase to the pathogenesis of several Flaviviridae virus infections [33], we investigated how ZIKV infection affects miRNA-dependent regulation of genes related to oxidative stress and ER stress. To first assess if oxidative stress is induced by ZIKV infection, the carbonylation of proteins by ROS was quantified using oxyblots. Thereby the level of oxidative modified proteins was significantly higher after infection with either ZIKV Uganda or Polynesia compared with uninfected cells at 48 and 72 hpi, indicating an induction of oxidative stress during ZIKV infection (Figure 9 and Figure S10). We next studied the effect of ZIKV on Nrf2/ARE signaling, an antioxidant pathway known to be hijacked by other viruses [31]. To this aim, we analyzed the influence of ZIKV infection on the amount and subcellular distribution of sMaf proteins by immunofluorescence analysis. We could observe a colocalization of ZIKV Uganda and Polynesia envelope proteins with sMaf, as well as an enrichment of sMaf proteins in perinuclear regions after infection (Figure 10). In addition, RT-qPCR analysis showed a significant upregulation of the antioxidant gene Sestrin 2 (Figure S8a). Interestingly, miR-182-5p was downregulated upon ZIKV infection which potentially targets the gene Sestrin2, as described by Lin et al. [53] (Figure S8b).
We next performed a gene ontology analysis to generate a list of differentially expressed genes upon infection with ZIKV Polynesia, which were significantly annotated to the terms response to oxidative stress and response to endoplasmic reticulum stress. Our analysis revealed 19 upregulated genes, which are listed in Table 5 and Table 6, including DDIT3 (DNA-damage-inducible transcript 3), XBP1 (X-box binding protein 1) and EIF2AK2 (eukaryotic translation initiation factor 2-alpha kinase 2), among others. To identify if these genes could be targeted by miRNAs, which were downregulated during infection with ZIKV Polynesia, we performed an integrative analysis using Tarbase. Our analysis revealed that miR-205-5p, which was downregulated at the intracellular level, could potentially target the gene ATF3 (activating transcription factor 3), besides ASNS. Furthermore, we found multiple EV-associated miRNAs with decreased levels upon infection with ZIKV Polynesia, which potentially regulate the same mRNA targets implicated in oxidative stress or ER stress (Table 7).
Taken together, we observed that during ZIKV infection, oxidative stress is induced, and we identified a number of potential miRNA-mRNA interactions related to oxidative stress and ER stress, which were differentially regulated and may contribute to ZIKV pathogenesis.

4. Discussion

ZIKV infection induces severe neurodevelopmental damage with microcephaly observed in children of women infected during pregnancy. However, the pathogenesis of these viral induced effects on neurogenesis are not yet fully understood. In this study, we used a human neural stem cell (NSC) model system and integrative methods to assess whether infection of NSCs with ZIKV may affect miRNA signatures and gene expression profiles. We analyzed intracellular host miRNAs and measured changes in levels of miRNAs incorporated into EVs in response to ZIKV infection. We identified a variety of signaling pathways that were altered during infection and may contribute to the development of brain malformations.
Despite the important role of miRNAs in the pathogenesis of flavivirus infections [54,55], only a few studies have been conducted to analyze how ZIKV could alter global miRNA regulation, and their potential involvement in ZIKV pathogenesis and microcephaly has not yet been fully elucidated [25,26,56]. Through global miRNA transcriptome analyses, we found many overlaps in the intracellular miRNA profiles modulated by the ZIKV strains Uganda and French Polynesia. However, ZIKV Uganda had a greater impact on mRNA and miRNA dysregulation, which might be a result of its higher infection level [25].
Both virus strains caused a dysregulation of numerous pathways involved in neurodevelopment and cell cycle, suggesting that they may contribute to the microcephaly phenotype. Interestingly, we found several differentially expressed genes following infection with ZIKV Uganda, including COL4A1, MIR17HG, TUBA1A, SLC2A1 and STAMBP, which are associated with microcephaly according to the comparative toxicogenomics database [48]. Accumulating evidence reveals that African strains of ZIKV can cause equal or worse damage to neuronal cells compared to Asian lineages [11,57]. Our findings support the hypothesis that ZIKV Uganda could also influence neurodevelopment and that clinical signs have remained undetected due to insufficient capacity to detect rare cases in Africa before 2013 [12,58]. However, another hypothesis is that differences between Asian and African ZIKV lineages regarding cytopathic effects might influence pathogenicity and that Asian ZIKV may cause chronic infections within the central nervous system due to its lower ability to cause cell death [59,60,61].
We further identified the gene ASNS to be differentially expressed by both strains, which has previously been associated with developmental disorders including microcephaly [62,63]. ASNS encodes the protein asparagine synthetase, a metabolic enzyme that converts aspartate and glutamine to asparagine and glutamate in an ATP-dependent reaction [64]. A dysregulation of the balance of amino acids in the brain induced by ZIKV could be a potential contributing factor to the microcephaly phenotype [62,63]. The overexpression of ASNS following ZIKV infection could be related to the downregulation of miR-205-5p, which has previously been described to target ASNS [51], although further experiments should be performed to study this interaction in lt-NES® cells.
Recent studies have highlighted the role of EV-associated miRNAs in the pathogenesis of several viruses, however, changes in EV-derived miRNA levels during ZIKV infection have not yet been described [65,66]. Our findings revealed that ZIKV infection of lt-NES® cells altered the signature of EV-associated miRNAs. Differentially incorporated miRNAs in EVs produced in response to ZIKV infection were associated with processes related to neural and embryonic development, including Axon guidance and FoxO signaling.
We found miR-4792 to be highly upregulated by both ZIKV strains at the intracellular level and significantly enriched in EVs secreted from ZIKV infected cells. Previous studies by other groups have described that overexpression of miR-4792 can suppress the growth of nasopharyngeal tumors in mice and induce cell cycle arrest and activate apoptosis in A459 cells [20,52]. It is possible that overexpression of miR-4792 following ZIKV infection could induce cell death of lt-NES® cells and through transmission via EVs this could even affect uninfected neighboring cells and thus be a possible mechanism for impaired development of the brain. Furthermore, miR-4792 has been described to directly target the transcription factor FOXC1 in both, nasopharyngeal carcinoma cells and A549 cells [20,52]. According to previous studies, FOXC1 is required for normal cerebellar development and impairment of FOXC1 function contributes to the development of Dandy–Walker syndrome [67,68], which is a brain malformation associated with ZIKV infection [69]. Since FOXC1 was downregulated in lt-NES® cells following ZIKV infection, it could be a contributing factor to the development of congenital ZIKV syndrome.
FOXC1 further plays a role in oxidative stress response, whereby decreased FOXC1 transcript and protein levels were found in HTM cells exposed to oxidative stress [70]. Oxidative stress has been implicated in various neurodegenerative disorders [71] and several Flaviviridae viruses are known to trigger ROS increase and to interfere with antioxidant systems, which affects the virus life cycle and the cellular metabolism [72,73]. Hepatitis C virus (HCV) impairs the Nrf2 mediated antioxidant response, which in turn favors the release of viral particles. In HCV-positive cells, sMafs are bound to NS3 (nonstructural protein 3) and delocalized from the nucleus to the replicon surface, thereby trapping free Nrf2 and preventing its entry into the nucleus [31,32,72]. Impaired Nrf2/ARE-signaling contributes to elevated ROS levels, which in turn trigger the induction of autophagy, a pathway that is often manipulated by positive stranded RNA viruses to facilitate their replication and viral release. The crosstalk between oxidative stress and autophagy can support exosome-mediated release of viral particles, since autophagy and exocytosis both involve membrane trafficking and are functionally connected [72,74,75].
In this study, we found significantly higher levels of oxidative modified proteins after ZIKV infection, reflecting the induction of oxidative stress following ZIKV infection. Immunofluorescence analysis showed a colocalization of sMaf proteins with either ZIKV Uganda or ZIKV Polynesia envelope proteins and an accumulation of sMafs in perinuclear regions. It is possible that ZIKV, similar to HCV, interferes with sMafs to impair the induction of ARE-dependent genes, since major Nrf2 target genes including HO-1, NQO1, GCL, GST were not dysregulated upon ZIKV infection despite the induction of oxidative stress; however, further studies are needed to assess this [76,77,78]. Furthermore, in this study we investigated the colocalization of sMafs with the structural envelope proteins of ZIKV and not with non-structural proteins, whereas in the case of HCV the non-structural proteins NS3 interfere with sMafs. While the expression of these cytoprotective genes, which are crucial for the defense against oxidative stress, was not altered, the antioxidant gene Sestrin2 (SESN2) was significantly upregulated after infection. SESN2 expression can be induced in response to ER stress by the transcription factors ATF4 (Activating Transcription Factor 4) or ATF6 (Activating Transcription Factor 6) and via the Nrf2/ARE pathway [79,80,81]. Consistent with our finding, Carr et al. recently described that SESN2 was upregulated upon infection with Japanese encephalitis virus and ZIKV Uganda in neuroblastoma cells [82]. SESN2 is known to regulate ROS generation and to activate autophagy via inhibition of mTORC1, making it relevant as a therapeutic target and biomarker in various diseases [83]. Previous studies have shown that SESN2 levels are altered in patients with diseases including neurodegenerative disorders and plasma levels of SESN2 are positively correlated with the severity of those diseases [84,85,86]. It has been shown that suppression of miR-182-5p can lead to upregulation of SESN2 in glioma and adenocarcinoma cell lines [53]. Interestingly, miR-182-5p was downregulated in lt-NES® cells following ZIKV infection, which could be associated with the upregulation of SESN2. However, we acknowledge the limitation that further studies including gain and loss of function should be performed to further assess the interaction between miRNAs and mRNAs.
Collectively, our findings highlight multiple dysregulated host miRNAs and mRNAs upon ZIKV infection with roles in neurodevelopment and oxidative stress, which may contribute to a better understanding of the neurodegenerative phenotype induced by ZIKV infection. Furthermore, the lt-NES® cells resulted as a suitable cell model to study dysregulation of miRNAs and genes by ZIKV infection and could be used for comparative studies of drug effects on neurodevelopment as well.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/11/1219/s1, Figure S1: Characterization of iLB-C-31f-r1 lt-NES® cells, Figure S2: SNP analysis of iLB-C-31f-r1 lt-NES® cells, Figure S3: Dysregulation of ASNS by ZIKV Uganda and Polynesia, Figure S4: Dysregulation of BST2 by ZIKV Uganda and Polynesia, Figure S5: Dysregulation of miR-205 by ZIKV Uganda and Polynesia, Figure S6: Dysregulation of miR-4792 by ZIKV Uganda and Polynesia, Figure S7: Dysregulation of FOXC1 by ZIKV Uganda and Polynesia, Figure S8: Dysregulation of Sestrin2 and miR-182 by ZIKV Uganda and Polynesia, Figure S9: Gene ontology (GO) analysis of potential target genes of differentially expressed miRNAs, Table S1: Top downregulated genes by ZIKV Uganda at 72 hpi, Table S2: Downregulated genes by ZIKV Polynesia at 72 hpi, Table S3: Intracellular mature miRNAs dysregulated by ZIKV Polynesia and Uganda, Table S4: Intracellular mature miRNAs dysregulated by ZIKV Uganda, Table S5: EV-associated mature miRNAs dysregulated by ZIKV Polynesia, Table S6: EV-associated mature miRNAs dysregulated by ZIKV Uganda, Figure S10: Protein oxidation during ZIKV infection.

Author Contributions

Conceptualization, D.T., C.S. (Catharina Scholl), A.T., O.B., F.E., E.H. and J.S.; methodology, D.T., C.S. (Catharina Scholl)., D.B., F.E., M.-L.H., C.S. (Catarina Sabino), A.T., V.S., O.B., E.H. and J.S.; software, D.T. and M.S.; validation, D.T., S.W., M.P.; formal analysis, D.T., M.S., A.T.; investigation, D.T., C.S. (Catharina Scholl), A.T., O.B., F.E., E.H. and J.S.; resources, M.P., A.T., V.S., D.B., M.-L.H., C.S. (Catarina Sabino), F.E., O.B., J.S. and E.H.; data curation, D.T., M.S., A.T.; writing—Original draft preparation, D.T.; writing—Review and editing, C.S. (Catharina Scholl), M.-L.H., D.B., E.H., J.S.; visualization, D.T.; supervision, O.B., E.H. and J.S.; project administration, E.H. and J.S.; funding acquisition, E.H. and J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a grant from the Federal Ministry of Health (BMG). Furthermore, this study was supported by the European Virus Archive goes Global (EVAg) project that has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 653316. Based on this funding, EVAg provided the ZIKV Uganda 976 strain. In addition, this work was supported by grants from the German Federal Ministry of Education and Research (grant 01EK1603A-Neuro2D3) and the European Regional Development Fund (grants EFRE-0800978, EFRE-0800977, Stem Cell Factory III).

Acknowledgments

The authors thank Dider Musso, Director of the Research and Diagnosis Laboratory, Institute Louis Malardé in Papeete, Tahiti, for donating the ZIKV French Polynesia strain (PF 13/251013-18). The authors also thank Michaela Segschneider for expert technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dick, G.W.A.; Kitchen, S.F.; Haddow, A.J. Zika Virus (I). Isolations and serological specificity. Trans. R. Soc. Trop. Med. Hyg. 1952, 46, 509–520. [Google Scholar] [CrossRef]
  2. Cauchemez, S.; Besnard, M.; Bompard, P.; Dub, T.; Guillemette-Artur, P.; Eyrolle-Guignot, D.; Salje, H.; Van Kerkhove, M.D.; Abadie, V.; Garel, C.; et al. Association between Zika virus and microcephaly in French Polynesia, 2013–15: A retrospective study. Lancet 2016, 387, 2125–2132. [Google Scholar] [CrossRef] [Green Version]
  3. Lanciotti, R.S.; Kosoy, O.L.; Laven, J.J.; Velez, J.O.; Lambert, A.J.; Johnson, A.J.; Stanfield, S.M.; Duffy, M.R. Genetic and Serologic Properties of Zika Virus Associated with an Epidemic, Yap State, Micronesia 2007. CDC Emerg. Infect. Dis. J. 2008, 14, 8. [Google Scholar] [CrossRef] [PubMed]
  4. Calvet, G.; Aguiar, R.S.; Melo, A.; A Sampaio, S.; De Filippis, I.; Fabri, A.; Araujo, E.S.M.; De Sequeira, P.C.; De Mendonça, M.C.L.; De Oliveira, L.; et al. Detection and sequencing of Zika virus from amniotic fluid of fetuses with microcephaly in Brazil: A case study. Lancet Infect. Dis. 2016, 16, 653–660. [Google Scholar] [CrossRef] [Green Version]
  5. Driggers, R.W.; Ho, C.-Y.; Korhonen, E.M.; Kuivanen, S.; Jääskeläinen, A.J.; Smura, T.; Rosenberg, A.; Hill, D.A.; DeBiasi, R.L.; Vezina, G.; et al. Zika Virus Infection with Prolonged Maternal Viremia and Fetal Brain Abnormalities. N. Engl. J. Med. 2016, 374, 2142–2151. [Google Scholar] [CrossRef] [PubMed]
  6. Mlakar, J.; Korva, M.; Tul, N.; Popovic, M.; Poljšak-Prijatelj, M.; Mraz, J.; Kolenc, M.; Rus, K.R.; Vipotnik, T.V.; Vodušek, V.F.; et al. Zika Virus Associated with Microcephaly. N. Engl. J. Med. 2016, 374, 951–958. [Google Scholar] [CrossRef]
  7. Oehler, E.; Watrin, L.; Larre, P.; Leparc-Goffart, I.; Lastère, S.; Valour, F.; Baudouin, L.; Mallet, H.P.; Musso, D.; Ghawche, F. Zika virus infection complicated by Guillain-Barré syndrome—Case report, French Polynesia, December 2013. Eurosurveillance 2014, 19, 20720. [Google Scholar] [CrossRef] [Green Version]
  8. Schuler-Faccini, L.; Ribeiro, E.M.; Feitosa, I.M.; Horovitz, D.D.; Cavalcanti, D.P.; Pessoa, A.; Doriqui, M.J.R.; Neri, J.I.; Neto, J.M.D.P.; Wanderley, H.Y.; et al. Possible Association Between Zika Virus Infection and Microcephaly—Brazil, 2015. MMWR Morb. Mortal. Wkly. Rep. 2016, 65, 59–62. [Google Scholar] [CrossRef]
  9. Zanluca, C.; De Melo, V.C.A.; Mosimann, A.L.P.; Dos Santos, G.I.V.; Dos Santos, C.N.D.; Luz, K. First report of autochthonous transmission of Zika virus in Brazil. Memórias Inst. Oswaldo Cruz 2015, 110, 569–572. [Google Scholar] [CrossRef]
  10. Zhang, Q.; Sun, K.; Chinazzi, M.; Piontti, A.P.Y.; Dean, N.E.; Rojas, D.P.; Merler, S.; Mistry, D.; Poletti, P.; Rossi, L.; et al. Spread of Zika virus in the Americas. Proc. Natl. Acad. Sci. USA 2017, 114, E4334–E4343. [Google Scholar] [CrossRef] [Green Version]
  11. Jaeger, A.S.; Murrieta, R.A.; Goren, L.R.; Crooks, C.M.; Moriarty, R.V.; Weiler, A.M.; Rybarczyk, S.; Semler, M.R.; Huffman, C.; Mejia, A.; et al. Zika viruses of African and Asian lineages cause fetal harm in a mouse model of vertical transmission. PLoS Neglected Trop. Dis. 2019, 13, e0007343. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Meda, N.; Salinas, S.; Kagoné, T.; Simonin, Y.; Van De Perre, P. Zika virus epidemic: Africa should not be neglected. Lancet 2016, 388, 337–338. [Google Scholar] [CrossRef] [Green Version]
  13. Heymann, D.L.; Hodgson, A.; Sall, A.A.; O Freedman, D.; Staples, J.E.; Althabe, F.; Baruah, K.; Mahmud, G.; Kandun, N.; Vasconcelos, P.F.C.; et al. Zika virus and microcephaly: Why is this situation a PHEIC? Lancet 2016, 387, 719–721. [Google Scholar] [CrossRef] [Green Version]
  14. Dang, J.; Tiwari, S.K.; Lichinchi, G.; Qin, Y.; Patil, V.S.; Eroshkin, A.M.; Rana, T.M. Zika Virus Depletes Neural Progenitors in Human Cerebral Organoids through Activation of the Innate Immune Receptor TLR3. Cell Stem Cell 2016, 19, 258–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Li, C.; Xu, D.; Ye, Q.; Hong, S.; Jiang, Y.; Liu, X.; Zhang, N.; Shi, L.; Qin, C.-F.; Xu, Z. Zika Virus Disrupts Neural Progenitor Development and Leads to Microcephaly in Mice. Cell Stem Cell 2016, 19, 672. [Google Scholar] [CrossRef] [Green Version]
  16. Ming, G.-L.; Tang, H.; Song, H. Advances in Zika Virus Research: Stem Cell Models, Challenges, and Opportunities. Cell Stem Cell 2016, 19, 690–702. [Google Scholar] [CrossRef] [Green Version]
  17. Tang, H.; Hammack, C.; Ogden, S.C.; Wen, Z.; Qian, X.; Li, Y.; Yao, B.; Shin, J.; Zhang, F.; Lee, E.M.; et al. Zika Virus Infects Human Cortical Neural Progenitors and Attenuates Their Growth. Cell Stem Cell 2016, 18, 587–590. [Google Scholar] [CrossRef] [Green Version]
  18. Bartel, D.P. MicroRNAs: Target Recognition and Regulatory Functions. Cell 2009, 136, 215–233. [Google Scholar] [CrossRef] [Green Version]
  19. Fu, G.; Brkić, J.; Hayder, H.; Peng, C. MicroRNAs in Human Placental Development and Pregnancy Complications. Int. J. Mol. Sci. 2013, 14, 5519–5544. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, P.; Pu, J.; Zhang, J.; Chen, Z.; Wei, K.; Shi, L. Bioinformatic analysis of miR-4792 regulates Radix Tetrastigma hemsleyani flavone to inhibit proliferation, invasion, and induce apoptosis of A549 cells. Onco Targets Ther. 2019, 12, 1401–1412. [Google Scholar] [CrossRef] [Green Version]
  21. Yao, S. MicroRNA biogenesis and their functions in regulating stem cell potency and differentiation. Biol. Proced. Online 2016, 18, 8. [Google Scholar] [CrossRef] [Green Version]
  22. Gurunathan, S.; Kang, M.-H.; Jeyaraj, M.; Qasim, M.; Kim, J.-H. Review of the Isolation, Characterization, Biological Function, and Multifarious Therapeutic Approaches of Exosomes. Cells 2019, 8, 307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Devhare, P.B.; Sasaki, R.; Shrivastava, S.; Di Bisceglie, A.M.; Ray, R.; Ray, R.B. Exosome-Mediated Intercellular Communication between Hepatitis C Virus-Infected Hepatocytes and Hepatic Stellate Cells. J. Virol. 2017, 91, e02225-16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mukherjee, S.; Akbar, I.; Kumari, B.; Vrati, S.; Basu, A.; Banerjee, A. Japanese Encephalitis Virus-induced let-7a/b interacted with the NOTCH—TLR 7 pathway in microglia and facilitated neuronal death via caspase activation. J. Neurochem. 2019, 149, 518–534. [Google Scholar] [CrossRef] [PubMed]
  25. Dang, J.W.; Tiwari, S.K.; Qin, Y.; Rana, T.M. Genome-wide Integrative Analysis of Zika-Virus-Infected Neuronal Stem Cells Reveals Roles for MicroRNAs in Cell Cycle and Stemness. Cell Rep. 2019, 27, 3618–3628. [Google Scholar] [CrossRef] [Green Version]
  26. Kozak, R.A.; Majer, A.; Biondi, M.J.; Medina, S.J.; Goneau, L.W.; Sajesh, B.V.; Slota, J.A.; Zubach, V.; Severini, A.; Safronetz, D.; et al. MicroRNA and mRNA Dysregulation in Astrocytes Infected with Zika Virus. Viruses 2017, 9, 297. [Google Scholar] [CrossRef] [Green Version]
  27. Valadão, A.L.C.; Aguiar, R.S.; De Arruda, L.B. Interplay between Inflammation and Cellular Stress Triggered by Flaviviridae Viruses. Front. Microbiol. 2016, 7, 1233. [Google Scholar] [CrossRef] [Green Version]
  28. Kobayashi, M.; Li, L.; Iwamoto, N.; Nakajima-Takagi, Y.; Kaneko, H.; Nakayama, Y.; Eguchi, M.; Wada, Y.; Kumagai, Y.; Yamamoto, M. The Antioxidant Defense System Keap1-Nrf2 Comprises a Multiple Sensing Mechanism for Responding to a Wide Range of Chemical Compounds. Mol. Cell. Biol. 2008, 29, 493–502. [Google Scholar] [CrossRef] [Green Version]
  29. Ungvari, Z.; Csiszar, A. Resveratrol Confers Endothelial Protection in Insulin-Dependent Diabetes Mellitus. Cardiovasc. Drugs Ther. 2011, 25, 111–113. [Google Scholar] [CrossRef]
  30. Yamamoto, M.; Kensler, T.W.; Motohashi, H. The KEAP1-NRF2 System: A Thiol-Based Sensor-Effector Apparatus for Maintaining Redox Homeostasis. Physiol. Rev. 2018, 98, 1169–1203. [Google Scholar] [CrossRef] [Green Version]
  31. Bender, D.; Hildt, E. Effect of Hepatitis Viruses on the Nrf2/Keap1-Signaling Pathway and Its Impact on Viral Replication and Pathogenesis. Int. J. Mol. Sci. 2019, 20, 4659. [Google Scholar] [CrossRef] [Green Version]
  32. Carvajal-Yepes, M.; Himmelsbach, K.; Schaedler, S.; Ploen, D.; Krause, J.; Ludwig, L.; Weiss, T.S.; Klingel, K.; Hildt, E. Hepatitis C Virus Impairs the Induction of Cytoprotective Nrf2 Target Genes by Delocalization of Small Maf Proteins. J. Biol. Chem. 2011, 286. [Google Scholar] [CrossRef] [Green Version]
  33. Zhang, Z.; Rong, L.; Li, Y.-P. Flaviviridae Viruses and Oxidative Stress: Implications for Viral Pathogenesis. Oxid. Med. Cell. Longev. 2019, 2019, 1409582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Koch, P.; Opitz, T.; Steinbeck, J.A.; Ladewig, J.; Brustle, O. A rosette-type, self-renewing human ES cell-derived neural stem cell with potential for in vitro instruction and synaptic integration. Proc. Natl. Acad. Sci. USA 2009, 106, 3225–3230. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Falk, A.; Koch, P.; Kesavan, J.; Takashima, Y.; Ladewig, J.; Alexander, M.; Wiskow, O.; Tailor, J.; Trotter, M.; Pollard, S.; et al. Capture of Neuroepithelial-Like Stem Cells from Pluripotent Stem Cells Provides a Versatile System for In Vitro Production of Human Neurons. PLoS ONE 2012, 7, e29597. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Rehbach, K.; Kesavan, J.; Hauser, S.; Ritzenhofen, S.; Jungverdorben, J.; Schüle, R.; Schöls, L.; Peitz, M.; Brüstle, O. Multiparametric rapid screening of neuronal process pathology for drug target identification in HSP patient-specific neurons. Sci. Rep. 2019, 9, 9615. [Google Scholar] [CrossRef] [PubMed]
  37. Tenoever, B.R. RNA viruses and the host microRNA machinery. Nat. Rev. Genet. 2013, 11, 169–180. [Google Scholar] [CrossRef] [PubMed]
  38. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3.0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef]
  39. Kanehisa, M. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  40. Reczko, M.; Maragkakis, M.; Alexiou, P.; Grosse, I.; Hatzigeorgiou, A.G. Functional microRNA targets in protein coding sequences. Bioinformatics 2012, 28, 771–776. [Google Scholar] [CrossRef]
  41. Agarwal, V.; Bell, G.W.; Nam, J.-W.; Bartel, D.P. Predicting effective microRNA target sites in mammalian mRNAs. eLife 2015, 4, e05005. [Google Scholar] [CrossRef] [PubMed]
  42. Krek, A.; Grün, D.; Poy, M.N.; Wolf, R.; Rosenberg, L.; Epstein, E.J.; MacMenamin, P.; Da Piedade, I.; Gunsalus, K.C.; Stoffel, M.; et al. Combinatorial microRNA target predictions. Nat. Genet. 2005, 37, 495–500. [Google Scholar] [CrossRef]
  43. Sethupathy, P.; Corda, B.; Hatzigeorgiou, A.G. TarBase: A comprehensive database of experimentally supported animal microRNA targets. RNA 2005, 12, 192–197. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wang, X. miRDB: A microRNA target prediction and functional annotation database with a wiki interface. RNA 2008, 14, 1012–1017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Croft, D.; Mundo, A.F.; Haw, R.; Milacic, M.; Weiser, J.; Wu, G.; Caudy, M.; Garapati, P.; Gillespie, M.; Kamdar, M.R.; et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2013, 42, D472–D477. [Google Scholar] [CrossRef] [PubMed]
  46. Schmittgen, T.D.; Livak, K.J. Analyzing real-time PCR data by the comparative CT method. Nat. Protoc. 2008, 3, 1101–1108. [Google Scholar] [CrossRef]
  47. Elgner, F.; Sabino, C.; Basic, M.; Ploen, D.; Grünweller, A.; Hildt, E. Inhibition of Zika Virus Replication by Silvestrol. Viruses 2018, 10, 149. [Google Scholar] [CrossRef] [Green Version]
  48. Davis, A.P.; Grondin, C.J.; Johnson, R.J.; Sciaky, D.; King, B.L.; McMorran, R.; Wiegers, J.; Wiegers, T.C.; Mattingly, C.J. The Comparative Toxicogenomics Database: Update 2017. Nucleic Acids Res. 2016, 45, D972–D978. [Google Scholar] [CrossRef] [Green Version]
  49. Li, M.; Wang, P.; Zheng, Z.; Hu, K.; Zhang, M.; Guan, X.; Fu, M.; Zhang, D.; Wang, W.; Xiao, G.; et al. Japanese encephalitis virus counteracts BST2 restriction via its envelope protein E. Virology 2017, 510, 67–75. [Google Scholar] [CrossRef]
  50. Cho, H.; Shrestha, B.; Sen, G.C.; Diamond, M.S. A Role for Ifit2 in Restricting West Nile Virus Infection in the Brain. J. Virol. 2013, 87, 8363–8371. [Google Scholar] [CrossRef] [Green Version]
  51. Xie, H.; Zhao, Y.; Caramuta, S.; Larsson, C.; Lui, W.-O. miR-205 Expression Promotes Cell Proliferation and Migration of Human Cervical Cancer Cells. PLoS ONE 2012, 7, e46990. [Google Scholar] [CrossRef] [PubMed]
  52. Li, Y.; Chen, X. miR-4792 inhibits epithelial–mesenchymal transition and invasion in nasopharyngeal carcinoma by targeting FOXC1. Biochem. Biophys. Res. Commun. 2015, 468, 863–869. [Google Scholar] [CrossRef] [PubMed]
  53. Lin, L.-T.; Liu, S.-Y.; Leu, J.-D.; Chang, C.-Y.; Chiou, S.-H.; Lee, T.-C.; Lee, Y.-J. Arsenic trioxide-mediated suppression of miR-182-5p is associated with potent anti-oxidant effects through up-regulation of SESN2. Oncotarget 2018, 9, 16028–16042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Lowey, B.; Hertz, L.; Chiu, S.; Valdez, K.; Li, Q.; Liang, T.J.; Weinman, S.; Saito, T. Hepatitis C Virus Infection Induces Hepatic Expression of NF-κB-Inducing Kinase and Lipogenesis by Downregulating miR-122. mBio 2019, 10, 19. [Google Scholar] [CrossRef] [Green Version]
  55. Zhu, X.; He, Z.; Hu, Y.; Wen, W.; Lin, C.; Yu, J.; Pan, J.; Li, R.; Deng, H.; Liao, S.; et al. MicroRNA-30e* Suppresses Dengue Virus Replication by Promoting NF-κB–Dependent IFN Production. PLoS Negl. Trop. Dis. 2014, 8, e3088. [Google Scholar] [CrossRef]
  56. Azouz, F.; Arora, K.; Krause, K.; Nerurkar, V.R.; Kumar, M. Integrated MicroRNA and mRNA Profiling in Zika Virus-Infected Neurons. Viruses 2019, 11, 162. [Google Scholar] [CrossRef] [Green Version]
  57. Nutt, C.; Adams, P. Zika in Africa—The invisible epidemic? Lancet 2017, 389, 1595–1596. [Google Scholar] [CrossRef]
  58. Simonin, Y.; Loustalot, F.; Desmetz, C.; Foulongne, V.; Constant, O.; Fournier-Wirth, C.; Leon, F.; Molès, J.-P.; Goubaud, A.; Lemaitre, J.-M.; et al. Zika Virus Strains Potentially Display Different Infectious Profiles in Human Neural Cells. EBioMedicine 2016, 12, 161–169. [Google Scholar] [CrossRef] [Green Version]
  59. Sheridan, M.A.; Balaraman, V.; Schust, D.J.; Ezashi, T.; Roberts, R.M.; Franz, A.W.E. African and Asian strains of Zika virus differ in their ability to infect and lyse primitive human placental trophoblast. PLoS ONE 2018, 13, e0200086. [Google Scholar] [CrossRef] [Green Version]
  60. Simonin, Y.; van Riel, D.; Van de Perre, P.; Rockx, B.; Salinas, S. Differential virulence between Asian and African lineages of Zika virus. PLoS Negl. Trop. Dis. 2017, 11. [Google Scholar] [CrossRef]
  61. Anfasa, F.; Siegers, J.Y.; Van Der Kroeg, M.; Mumtaz, N.; Raj, V.S.; De Vrij, F.M.S.; Widagdo, W.; Gabriel, G.; Salinas, S.; Simonin, Y.; et al. Phenotypic Differences between Asian and African Lineage Zika Viruses in Human Neural Progenitor Cells. mSphere 2017, 2, e00292-17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Lomelino, C.L.; Andring, J.T.; McKenna, R.; Kilberg, M.S. Asparagine synthetase: Function, structure, and role in disease. J. Biol. Chem. 2017, 292, 19952–19958. [Google Scholar] [CrossRef] [Green Version]
  63. Schleinitz, D.; Seidel, A.; Stassart, R.; Klammt, J.; Hirrlinger, P.G.; Winkler, U.; Köhler, S.; Heiker, J.T.; Schönauer, R.; Bialek, J.; et al. Novel Mutations in the Asparagine Synthetase Gene (ASNS) Associated With Microcephaly. Front. Genet. 2018, 9, 245. [Google Scholar] [CrossRef] [PubMed]
  64. Richards, N.G.J.; Kilberg, M.S. Asparagine Synthetase Chemotherapy. Annu. Rev. Biochem. 2006, 75, 629–654. [Google Scholar] [CrossRef] [Green Version]
  65. Slonchak, A.; Clarke, B.; MacKenzie, J.; Amarilla, A.A.; Setoh, Y.X.; Khromykh, A.A. West Nile virus infection and interferon alpha treatment alter the spectrum and the levels of coding and noncoding host RNAs secreted in extracellular vesicles. BMC Genom. 2019, 20, 474. [Google Scholar] [CrossRef] [Green Version]
  66. Yoshikawa, F.S.Y.; Teixeira, F.M.E.; Sato, M.N.; Oliveira, L.M.S. Delivery of microRNAs by Extracellular Vesicles in Viral Infections: Could the News be Packaged? Cells 2019, 8, 611. [Google Scholar] [CrossRef] [Green Version]
  67. A Aldinger, K.; Lehmann, O.J.; Hudgins, L.; Chizhikov, V.V.; Bassuk, A.G.; Ades, L.C.; Krantz, I.D.; Dobyns, W.B.; Millen, K.J. FOXC1 is required for normal cerebellar development and is a major contributor to chromosome 6p25.3 Dandy-Walker malformation. Nat. Genet. 2009, 41, 1037–1042. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Haldipur, P.; Dang, D.; Aldinger, K.A.; Janson, O.K.; Guimiot, F.; Adle-Biasette, H.; Dobyns, W.B.; Siebert, J.R.; Russo, R.; Millen, K.J. Phenotypic outcomes in Mouse and Human Foxc1 dependent Dandy-Walker cerebellar malformation suggest shared mechanisms. eLife 2017, 6, 321. [Google Scholar] [CrossRef] [Green Version]
  69. Pool, K.-L.; Adachi, K.; Karnezis, S.; Salamon, N.; Romero, T.; Nielsen-Saines, K.; Pone, S.; Boechat, M.; Aibe, M.; Da Silva, T.G.; et al. Association Between Neonatal Neuroimaging and Clinical Outcomes in Zika-Exposed Infants From Rio de Janeiro, Brazil. JAMA Netw. Open 2019, 2, e198124. [Google Scholar] [CrossRef]
  70. Ito, Y.A.; Goping, I.S.; Berry, F.; Walter, M.A. Dysfunction of the stress-responsive FOXC1 transcription factor contributes to the earlier-onset glaucoma observed in Axenfeld-Rieger syndrome patients. Cell Death Dis. 2014, 5, e1069. [Google Scholar] [CrossRef] [Green Version]
  71. Singh, A.; Kukreti, R.; Saso, L.; Kukreti, S. Oxidative Stress: A Key Modulator in Neurodegenerative Diseases. Molecules 2019, 24, 1583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Medvedev, R.; Ploen, D.; Spengler, C.; Elgner, F.; Ren, H.; Bunten, S.; Hildt, E. HCV-induced oxidative stress by inhibition of Nrf2 triggers autophagy and favors release of viral particles. Free Radic. Biol. Med. 2017, 110, 300–315. [Google Scholar] [CrossRef] [PubMed]
  73. Olagnier, D.; Peri, S.; Steel, C.; Van Montfoort, N.; Chiang, C.; Beljanski, V.; Slifker, M.; He, Z.; Nichols, C.N.; Lin, R.; et al. Cellular Oxidative Stress Response Controls the Antiviral and Apoptotic Programs in Dengue Virus-Infected Dendritic Cells. PLoS Pathog. 2014, 10, e1004566. [Google Scholar] [CrossRef] [PubMed]
  74. Choi, Y.; Bowman, J.W.; Jung, J.U. Autophagy during viral infection—A double-edged sword. Nat. Rev. Genet. 2018, 16, 341–354. [Google Scholar] [CrossRef] [PubMed]
  75. Dreux, M.; Chisari, F.V. Viruses and the autophagy machinery. Cell Cycle 2010, 9, 1295–1307. [Google Scholar] [CrossRef] [PubMed]
  76. Itoha, K.; Chibabc, T.; Takahashia, S.; Ishiia, T.; Igarashia, K.; Katoha, Y.; Oyaked, T.; Hayashid, N.; Satohe, K.; Hatayamae, I.; et al. An Nrf2/Small Maf Heterodimer Mediates the Induction of Phase II Detoxifying Enzyme Genes through Antioxidant Response Elements. Biochem. Biophys. Res. Commun. 1997, 236, 313–322. [Google Scholar] [CrossRef]
  77. Rushworth, S.A.; MacEwan, D.; O’Connell, M.A. Lipopolysaccharide-induced expression of NAD(P)H:quinone oxidoreductase 1 and heme oxygenase-1 protects against excessive inflammatory responses in human monocytes. J. Immunol. 2008, 181, 6730–6737. [Google Scholar] [CrossRef] [Green Version]
  78. Vasiliou, V.; Ross, D.; Nebert, D.W. Update of the NAD(P)H:quinone oxidoreductase (NQO) gene family. Hum. Genom. 2006, 2, 329–335. [Google Scholar] [CrossRef] [Green Version]
  79. Shin, B.Y.; Jin, S.H.; Cho, I.J.; Ki, S.H. Nrf2-ARE pathway regulates induction of Sestrin-2 expression. Free Radic. Biol. Med. 2012, 53, 834–841. [Google Scholar] [CrossRef]
  80. Ding, B.; Parmigiani, A.; Divakaruni, A.S.; Archer, K.; Murphy, A.N.; Budanov, A.V. Sestrin2 is induced by glucose starvation via the unfolded protein response and protects cells from non-canonical necroptotic cell death. Sci. Rep. 2016, 6, 22538. [Google Scholar] [CrossRef] [Green Version]
  81. Jegal, K.H.; Park, S.M.; Cho, S.S.; Byun, S.H.; Ku, S.-K.; Kim, S.C.; Ki, S.H.; Cho, I.J. Activating transcription factor 6-dependent sestrin 2 induction ameliorates ER stress-mediated liver injury. Biochim. Biophys. Acta Bioenerg. 2017, 1864, 1295–1307. [Google Scholar] [CrossRef] [PubMed]
  82. Carr, M.; Gonzalez, G.; Martinelli, A.; Wastika, C.E.; Ito, K.; Orba, Y.; Sasaki, M.; Hall, W.W.; Sawa, H. Upregulated expression of the antioxidant sestrin 2 identified by transcriptomic analysis of Japanese encephalitis virus-infected SH-SY5Y neuroblastoma cells. Virus Genes 2019, 55, 630–642. [Google Scholar] [CrossRef] [PubMed]
  83. Pasha, M.; Eid, A.H.; Eid, A.A.; Gorin, Y.; Munusamy, S. Sestrin2 as a Novel Biomarker and Therapeutic Target for Various Diseases. Available online: https://www.hindawi.com/journals/omcl/2017/3296294/ (accessed on 3 June 2020).
  84. Soontornniyomkij, V.; Soontornniyomkij, B.; Moore, D.J.; Gouaux, B.; Masliah, E.; Tung, S.; Vinters, H.V.; Grant, I.; Achim, C.L. Antioxidant sestrin-2 redistribution to neuronal soma in human immunodeficiency virus-associated neurocognitive disorders. J. Neuroimmune Pharmacol. 2012, 7, 579–590. [Google Scholar] [CrossRef] [PubMed]
  85. Xiao, T.; Zhang, L.; Huang, Y.; Shi, Y.; Wang, J.; Ji, Q.; Ye, J.; Lin, Y.; Liu, H. Sestrin2 increases in aortas and plasma from aortic dissection patients and alleviates angiotensin II-induced smooth muscle cell apoptosis via the Nrf2 pathway. Life Sci. 2019, 218, 132–138. [Google Scholar] [CrossRef]
  86. Ye, J.; Wang, M.; Xu, Y.; Liu, J.; Jiang, H.; Wang, Z.; Lin, Y.; Wan, J. Sestrins increase in patients with coronary artery disease and associate with the severity of coronary stenosis. Clin. Chim. Acta 2017, 472, 51–57. [Google Scholar] [CrossRef]
Viruses 12 01219 g001 550
Figure 1. Differentially expressed genes (DEGs). Number of ≥ two-fold that were significant (p-value ≤ 0.05) and differentially expressed genes at 24, 48 and 72 h post infection with ZIKV Polynesia and ZIKV Uganda.
Figure 1. Differentially expressed genes (DEGs). Number of ≥ two-fold that were significant (p-value ≤ 0.05) and differentially expressed genes at 24, 48 and 72 h post infection with ZIKV Polynesia and ZIKV Uganda.
Viruses 12 01219 g001
Viruses 12 01219 g002 550
Figure 2. Plaque forming assays. Quantification was performed in Vero cells by plaque forming assays and the number of infectious viral particles is expressed in plaque forming units per mL (PFU/mL). Error bar represents the mean ± SEM of n = 3 biological replicates; * p < 0.05. Quantification of the number of extracellular infectious viral particles from lt-NES® cells infected with either the Polynesia or the Uganda strain at 24, 48 and 72 hpi at multiplicity of infection (MOI) of 0.1.
Figure 2. Plaque forming assays. Quantification was performed in Vero cells by plaque forming assays and the number of infectious viral particles is expressed in plaque forming units per mL (PFU/mL). Error bar represents the mean ± SEM of n = 3 biological replicates; * p < 0.05. Quantification of the number of extracellular infectious viral particles from lt-NES® cells infected with either the Polynesia or the Uganda strain at 24, 48 and 72 hpi at multiplicity of infection (MOI) of 0.1.
Viruses 12 01219 g002
Viruses 12 01219 g003 550
Figure 3. Top 20 enriched pathways identified by functional reactome analysis. Pathway analysis of ≥ two-fold up- and downregulated genes upon infection with either (a) ZIKV Polynesia or (b) ZIKV Uganda. Pathways ranked according to −log10 p-value.
Figure 3. Top 20 enriched pathways identified by functional reactome analysis. Pathway analysis of ≥ two-fold up- and downregulated genes upon infection with either (a) ZIKV Polynesia or (b) ZIKV Uganda. Pathways ranked according to −log10 p-value.
Viruses 12 01219 g003
Viruses 12 01219 g004a 550Viruses 12 01219 g004b 550
Figure 4. Reactome pathway diagram showing enriched signaling pathways as red lines. Genome-wide overview pathway analysis of ≥ two-fold downregulated genes for (a) ZIKV Polynesia and (b) ZIKV Uganda. The red lines denote over-representation of those pathways.
Figure 4. Reactome pathway diagram showing enriched signaling pathways as red lines. Genome-wide overview pathway analysis of ≥ two-fold downregulated genes for (a) ZIKV Polynesia and (b) ZIKV Uganda. The red lines denote over-representation of those pathways.
Viruses 12 01219 g004aViruses 12 01219 g004b
Viruses 12 01219 g005 550
Figure 5. Dysregulated miRNAs. (a) Number of more than two-fold, significantly (p-value ≤ 0.05) dysregulated miRNAs by ZIKV Polynesia and Uganda at 24, 48 and 72 hpi. (b) Significantly more than two-fold dysregulated miRNAs by ZIKV Polynesia and Uganda including all time points.
Figure 5. Dysregulated miRNAs. (a) Number of more than two-fold, significantly (p-value ≤ 0.05) dysregulated miRNAs by ZIKV Polynesia and Uganda at 24, 48 and 72 hpi. (b) Significantly more than two-fold dysregulated miRNAs by ZIKV Polynesia and Uganda including all time points.
Viruses 12 01219 g005
Viruses 12 01219 g006 550
Figure 6. Gene ontology (GO) analysis of potential target genes of differentially expressed miRNAs. (a) GO analysis of biological processes of putative targets of miRNAs which were dysregulated in lt-NES® cells during infection with ZIKV Polynesia and Uganda. (b) GO analysis of biological processes of putative targets of miR-7704, which was upregulated by ZIKV Polynesia and Uganda.
Figure 6. Gene ontology (GO) analysis of potential target genes of differentially expressed miRNAs. (a) GO analysis of biological processes of putative targets of miRNAs which were dysregulated in lt-NES® cells during infection with ZIKV Polynesia and Uganda. (b) GO analysis of biological processes of putative targets of miR-7704, which was upregulated by ZIKV Polynesia and Uganda.
Viruses 12 01219 g006
Viruses 12 01219 g007 550
Figure 7. Level of extracellular vesicles (EV)-derived miRNAs. (a) Significantly (p-value ≤ 0.05) ≥ two-fold dysregulated mature miRNAs, which were incorporated into EVs during infection of lt-NES® cells with either ZIKV Polynesia or Uganda. (b) EV-derived miRNAs with levels significantly (p-value ≤ 0.05) altered (≥two-fold) following infection with ZIKV Polynesia and Uganda at any time point.
Figure 7. Level of extracellular vesicles (EV)-derived miRNAs. (a) Significantly (p-value ≤ 0.05) ≥ two-fold dysregulated mature miRNAs, which were incorporated into EVs during infection of lt-NES® cells with either ZIKV Polynesia or Uganda. (b) EV-derived miRNAs with levels significantly (p-value ≤ 0.05) altered (≥two-fold) following infection with ZIKV Polynesia and Uganda at any time point.
Viruses 12 01219 g007
Viruses 12 01219 g008 550
Figure 8. KEGG Pathway analysis of potential targets of EV-derived miRNAs which were differentially incorporated following infection with (a) ZIKV Polynesia and (b) ZIKV Uganda.
Figure 8. KEGG Pathway analysis of potential targets of EV-derived miRNAs which were differentially incorporated following infection with (a) ZIKV Polynesia and (b) ZIKV Uganda.
Viruses 12 01219 g008
Viruses 12 01219 g009 550
Figure 9. ZIKV infection induces elevated levels of oxidized proteins. Oxyblot analysis of lysates derived from uninfected lt-NES® cells and lt-NES® cells, which were infected with ZIKV Uganda and ZIKV Polynesia at MOI of 1 at 24, 48 and 72 hpi. Protein oxidation was analyzed by incubation with 2,4,-dinitrophenylhydrazine (DNPH) for covalent modification of oxidative modified proteins. Error bar represents the mean ± SEM of n = 3 biological replicates; * p < 0.05; ** p < 0.01.
Figure 9. ZIKV infection induces elevated levels of oxidized proteins. Oxyblot analysis of lysates derived from uninfected lt-NES® cells and lt-NES® cells, which were infected with ZIKV Uganda and ZIKV Polynesia at MOI of 1 at 24, 48 and 72 hpi. Protein oxidation was analyzed by incubation with 2,4,-dinitrophenylhydrazine (DNPH) for covalent modification of oxidative modified proteins. Error bar represents the mean ± SEM of n = 3 biological replicates; * p < 0.05; ** p < 0.01.
Viruses 12 01219 g009
Viruses 12 01219 g010 550
Figure 10. Increased amount of sMaf proteins in ZIKV infected cells. Confocal immunofluorescence microscopy of ZIKV-positive cells (ZIKV Uganda at 48 hpi and ZIKV Polynesia at 24 hpi) or negative control cells (24 hpi). The immunofluorescence staining was performed using the polyclonal sMaf-specific serum and an envelope specific antibody. Nuclei were stained with DAPI.
Figure 10. Increased amount of sMaf proteins in ZIKV infected cells. Confocal immunofluorescence microscopy of ZIKV-positive cells (ZIKV Uganda at 48 hpi and ZIKV Polynesia at 24 hpi) or negative control cells (24 hpi). The immunofluorescence staining was performed using the polyclonal sMaf-specific serum and an envelope specific antibody. Nuclei were stained with DAPI.
Viruses 12 01219 g010
Table
Table 1. Comparative toxicogenomics database screening.
Table 1. Comparative toxicogenomics database screening.
GeneSymbolFold Change ZIKV Ugandap-Value
ZIKV Uganda		Fold Change ZIKV Polynesiap-Value
ZIKV Polynesia
COL4A1−2.10.0093−1.20.635
MIR17HG−2.30.0497−1.30.749
TUBA1A−2.22.56 × 10−4−1.50.256
SLC2A1−2.30.02391.00.806
ASNS9.60.03142.90.006
STAMBP2.51.07 × 10−41.40.975
Genes dysregulated by ZIKV Uganda > two-fold and by ZIKV Polynesia, which are associated with microcephaly according to the comparative toxicogenomics database [48].
Table
Table 2. Top upregulated genes by ZIKV Uganda.
Table 2. Top upregulated genes by ZIKV Uganda.
GeneSymbolFold Change
ZIKV Uganda 24 hpi		Fold Change
ZIKV Uganda 48 hpi		Fold Change
ZIKV Uganda 72 hpi
OASL1.649.9766.7
CH25H2.546.8611.3
IFNB12.337.9554.3
CXCL10−1.340.8521.3
IFIT21.030.7218.8
IFIH12.138.1176.4
PRDM168.123.5141.9
CCL5−1.310.3117.1
CXCL111.112.8104.1
GSC1.29.060.4
BST21.826.659.0
IFI442.420.855.5
PKD1L22.89.649.9
DHX581.38.346.3
ATP4A1.13.738.3
TMEM712.58.735.8
CEACAM11.46.835.8
NEURL31.210.534.0
IFI61.09.133.8
Table
Table 3. Top upregulated genes by ZIKV Polynesia.
Table 3. Top upregulated genes by ZIKV Polynesia.
GeneSymbolFold Change
ZIKV Polynesia 24 hpi		Fold Change
ZIKV Polynesia 48 hpi		Fold Change
ZIKV Polynesia 72 hpi
BST21.210.2459.5
OAS1−1.02.7261.2
IFIH11.511.8219.1
IFNB11.44.3211.2
OASL−1.16.3193.1
CXCL101.04.2147.1
IFITM11.44.8133.1
HERC61.13.2130.7
CH25H1.53.0126.7
IFI442.111.6124.7
IFI6−1.03.7118.3
SAMD9L−1.06.2105.6
IFIT21.35.196.8
MX11.14.390.2
ISG15−1.05.167.5
OAS2−1.12.454.4
DDX601.04.352.0
Table
Table 4. miR-4792: Diana mirPath v3 gne ontology analysis.
Table 4. miR-4792: Diana mirPath v3 gne ontology analysis.
GO Categoryp-Value#Genes
Gene expression2.49 × 10−512
Cellular nitrogen compound metabolic process0.0001337
Macromolecular complex assembly0.0016213
Protein complex assembly0.0101211
Cellular component assembly0.0198614
Biosynthetic process0.0267328
Biological process0.0290190
Negative regulation of oxidoreductase activity0.038212
In utero embryonic development0.038217
Biological processes, in which potential targets of miR-4792 are involved, which were predicted by TarBase.
Table
Table 5. Differentially expressed genes (DEGs) during infection with ZIKV Polynesia at 72 hpi, which were annotated to GO-term ”response to oxidative stress”.
Table 5. Differentially expressed genes (DEGs) during infection with ZIKV Polynesia at 72 hpi, which were annotated to GO-term ”response to oxidative stress”.
DEGFold Change ZIKV Polynesia
PNPT13.36
FOS5.18
XBP12.06
PPARGC1A2.06
KLF46.28
PML3.76
TNFAIP34.96
STAT117.09
SLC7A112.11
JUN2.10
Table
Table 6. Differentially expressed genes (DEGs) during infection with ZIKV Polynesia at 72 hpi, which were annotated to GO-term ”response to endoplasmic reticulum stress”.
Table 6. Differentially expressed genes (DEGs) during infection with ZIKV Polynesia at 72 hpi, which were annotated to GO-term ”response to endoplasmic reticulum stress”.
DEGFold Change ZIKV Polynesia
XBP12.06
EIF2AK24.72
ASNS2.13
TRIB32.84
DDIT35.17
ATF35.56
HERPUD12.48
PPP1R15A2.31
FBXO62.33
CHAC14.46
Table
Table 7. EV associated miRNAs with decreased levels and their potential mRNA targets.
Table 7. EV associated miRNAs with decreased levels and their potential mRNA targets.
miRNA IDGene Targets (ID)
hsa-miR-26b-5pTRIB3, FOS, HERPUD1, TNFAIP3, JUN, ASNS, CHAC1, ATF3
hsa-let-7g-5pTRIB3, FBXO6, CHAC1
hsa-miR-17-5pFOS, TNFAIP3, JUN, KLF4
hsa-miR-221-3pFOS, PPP1R15A
hsa-miR-16-5pHERPUD1, PML, EIF2AK2, JUN, ASNS, PPP1R15A, CHAC1, STAT1
hsa-miR-106b-5pTNFAIP3
hsa-miR-374a-5pTNFAIP3, EIF2AK2
hsa-miR-93-5pTNFAIP3, JUN
hsa-miR-20a-5pTNFAIP3
hsa-miR-186-5pEIF2AK2, JUN, ASNS, CHAC1
hsa-miR-32-5pJUN, KLF4
hsa-miR-218-5pCHAC1
hsa-miR-340-5pPPARGC1A
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tabari, D.; Scholl, C.; Steffens, M.; Weickhardt, S.; Elgner, F.; Bender, D.; Herrlein, M.-L.; Sabino, C.; Semkova, V.; Peitz, M.; et al. Impact of Zika Virus Infection on Human Neural Stem Cell MicroRNA Signatures. Viruses 2020, 12, 1219. https://doi.org/10.3390/v12111219

AMA Style

Tabari D, Scholl C, Steffens M, Weickhardt S, Elgner F, Bender D, Herrlein M-L, Sabino C, Semkova V, Peitz M, et al. Impact of Zika Virus Infection on Human Neural Stem Cell MicroRNA Signatures. Viruses. 2020; 12(11):1219. https://doi.org/10.3390/v12111219

Chicago/Turabian Style

Tabari, Denna, Catharina Scholl, Michael Steffens, Sandra Weickhardt, Fabian Elgner, Daniela Bender, Marie-Luise Herrlein, Catarina Sabino, Vesselina Semkova, Michael Peitz, and et al. 2020. "Impact of Zika Virus Infection on Human Neural Stem Cell MicroRNA Signatures" Viruses 12, no. 11: 1219. https://doi.org/10.3390/v12111219

APA Style

Tabari, D., Scholl, C., Steffens, M., Weickhardt, S., Elgner, F., Bender, D., Herrlein, M. -L., Sabino, C., Semkova, V., Peitz, M., Till, A., Brüstle, O., Hildt, E., & Stingl, J. (2020). Impact of Zika Virus Infection on Human Neural Stem Cell MicroRNA Signatures. Viruses, 12(11), 1219. https://doi.org/10.3390/v12111219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop