Next Article in Journal
Pathogenicity and Transmissibility of North American H7 Low Pathogenic Avian Influenza Viruses in Chickens and Turkeys
Next Article in Special Issue
The Application and Interpretation of IgG Avidity and IgA ELISA Tests to Characterize Zika Virus Infections
Previous Article in Journal
Ebola Virus Isolation Using Huh-7 Cells has Methodological Advantages and Similar Sensitivity to Isolation Using Other Cell Types and Suckling BALB/c Laboratory Mice
Previous Article in Special Issue
The Roles of prM-E Proteins in Historical and Epidemic Zika Virus-mediated Infection and Neurocytotoxicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Integrated MicroRNA and mRNA Profiling in Zika Virus-Infected Neurons

1
Department of Tropical Medicine, Medical Microbiology and Pharmacology, Pacific Center for Emerging Infectious Diseases Research, John A. Burns School of Medicine, University of Hawai’i at Mānoa, Honolulu, HI 96813, USA
2
Department of Biology, College of Arts and Sciences, Georgia State University, Atlanta, GA 30303, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Viruses 2019, 11(2), 162; https://doi.org/10.3390/v11020162
Submission received: 10 January 2019 / Revised: 8 February 2019 / Accepted: 14 February 2019 / Published: 16 February 2019
(This article belongs to the Special Issue New Advances on Zika Virus Research)

Abstract

:
Zika virus (ZIKV) infections have caused a wide spectrum of neurological diseases, such as Guillain-Barré syndrome, myelitis, meningoencephalitis, and congenital microcephaly. No effective therapies currently exist for treating patients infected with ZIKV. MicroRNAs (miRNAs) are a group of small RNAs involved in the regulation of a wide variety of cellular and physiological processes. In this study, we analyzed digital miRNA and mRNA profiles in ZIKV-infected primary mouse neurons using the nCounter technology. A total of 599 miRNAs and 770 mRNAs were examined. We demonstrate that ZIKV infection causes global downregulation of miRNAs with only few upregulated miRNAs. ZIKV-modulated miRNAs including miR-155, miR-203, miR-29a, and miR-124-3p are known to play critical role in flavivirus infection, anti-viral immunity and brain injury. ZIKV infection also results in downregulation of miRNA processing enzymes. In contrast, ZIKV infection induces dramatic upregulation of anti-viral, inflammatory and apoptotic genes. Furthermore, our data demonstrate an inverse correlation between ZIKV-modulated miRNAs and target host mRNAs induced by ZIKV. Biofunctional analysis revealed that ZIKV-modulated miRNAs and mRNAs regulate the pathways related to neurological development and neuroinflammatory responses. Functional studies targeting specific miRNA are warranted to develop therapeutics for the management of ZIKV neurological disease.

1. Introduction

Zika virus (ZIKV) is an emerging mosquito-borne pathogen that is part of the Spondweni serocomplex of the genus Flavivirus, family Flaviviridae. ZIKV is closely related to other pathogens of public health importance including yellow fever virus (YFV), dengue virus (DENV), Japanese encephalitis virus (JEV), and West Nile virus (WNV). The ZIKV genome is comprised of a single-stranded, positive-sense 11-kb RNA that contains three structural and seven nonstructural genes [1,2]. ZIKV is highly neurotropic in human fetal infections and has been linked to the development of severe fetal abnormalities that include spontaneous abortion, stillbirth, hydranencephaly, microcephaly, and placental insufficiency that may cause intrauterine growth restriction [1]. An increased incidence of Guillain-Barré syndrome (GBS), neuropathy of the peripheral nervous system, has also been reported in ZIKV-infected patients [3]. No effective therapies currently exist for treating patients infected with ZIKV.
ZIKV has been shown to replicate and induce cell death in neuronal cells of fetal mice as well as in human neural progenitor cells and brain organoids, a mechanism thought to play an important role in the pathogenesis of ZIKV neurological disease [4,5,6,7]. It is known that immunocompetent adult mice are resistant to subcutaneous or intraperitoneal inoculation of ZIKV [8,9]. However, it has been demonstrated that intracerebral inoculation of ZIKV in adult immunocompetent mice results in neurological disease [10]. Additionally, several studies have reported that neonatal immunocompetent mice inoculated with ZIKV via subcutaneously or intracerebral route develop ZIKV disease, and ZIKV infection can be detected in the neurons [11,12,13]. It has also been demonstrated that ZIKV infection during the period of maximal brain growth causes microcephaly and corticospinal neuron apoptosis in wild-type mice [14]. However, to date the effect of ZIKV infection on microRNAs (miRNAs) expression in primary mouse neurons has not been examined.
miRNAs are a group of small RNAs involved in the regulation of a wide variety of cellular and physiological processes. miRNAs are considered novel diagnostic and interventional candidates due to their biochemical structure [15]. They function by directly binding to the 3’ untranslated regions (3’UTRs) of specific target mRNA, causing a block of translation or degradation of the target mRNA. miRNAs have been demonstrated to play a crucial regulatory role in neurodegenerative diseases such as Alzheimer and Parkinson [16]. miRNAs also play a critical role in the regulation of immune response; including differentiation, proliferation, cell fate determination, function of immune cells, and inflammatory mediator release as well as the intracellular signaling pathways [17,18].
miRNAs of infected cells can influence the ability of a virus to replicate or spread. It is known that endogenous miRNAs inhibit replication of a number of RNA viruses including HIV-1, Ebola virus and vesicular stomatitis virus [19,20,21,22,23]. For example, miR-28, miR-125b, miR-150, miR- 223, miR-198, and miR-382 inhibit HIV replication in CD4 T cells by directly targeting HIV RNA or by modulating cellular factors responsible for its replication [24]. Furthermore, miR-122 supports HCV replication by enhancing colony formation efficiency of HCV [25], whereas miR-196 and miR-296 substantially attenuate virus replication through type I interferon (IFN)-associated pathways in liver cells [26]. Over-expression of miRNA-30e, let-7c, and miRNA-126-5p inhibits DENV replication [27,28,29]. Cellular miR-532-5p inhibits WNV replication via suppression of host genes SESTD1 and TAB3 required for virus replication [30]. miRNA HS_154 contributes to WNV-mediated apoptosis in vitro in the human neuronal cell line, SK-N-MC [31]. Moreover, incorporation of a target sequence for cellular microRNAs expressed in the central nervous system (CNS) into the flavivirus genome alters the neurovirulence of the virus and prevents the development of lethal encephalitis in mice [32].
We have previously demonstrated the role of cellular miRNAs in the pathogenesis of WNV encephalitis [33,34]. In this study, we analyzed miRNA and mRNA profiles in ZIKV-infected neurons using the nCounter system. Unlike other platforms (such as microarray and next generation sequencing), the nCounter platform enables high throughput, sensitive, quantitative, and reproducible gene expression analysis without the need of enzymatic target amplification. This technology utilizes 100 nucleotide molecular bar codes which measure gene quantities without an amplification step [35]. To our knowledge, this is the first study to evaluate the modulation of miRNAs following ZIKV infection using relevant cells, primary neurons.

2. Materials and Methods

2.1. Neuronal Cultures

Mouse cortical neuron cultures were prepared from one-day old pups of either gender (approximately equivalent numbers) obtained from established colonies of wild-type C57BL/6J mice as described previously [36]. The neurons were plated onto poly-D-lysine-coated 6-well or 24-well plates in serum Neurobasal A medium. The cultures were maintained in serum-free Neurobasal A medium supplemented with B27 for seven days prior to infection [36,37]. This study was carried out in accordance with the recommendations of the National Institutes of Health and the Institutional Animal Care and Use Committee (IACUC). The protocol was approved by the University of Hawaii IACUC (protocol number 17-2721) and Georgia State University (protocol number A19005).

2.2. ZIKV Infection and Plaque Assay

In this study, we used a low cell-culture-passaged and sequence-verified ZIKV strain, PRVABC59 (BEI Resources, NR-50240). Virus strain was amplified once in Vero E6 cells and had titers of 5 × 106 plaque-forming units (PFU)/mL. Cells were infected with ZIKV or PBS (Mock) at multiplicity of infection (MOI)-1 and supernatants and cell lysates were harvested at 12, 24, 48, and 72 h after infection [38]. ZIKV titers were measured in cell supernatants using plaque assay [38,39]. Experiments were repeated to obtain four biological replicates of mock- and ZIKV-infected neurons at each time point for the NanoString analysis (n = 4 per group per time point).

2.3. Indirect-Immunofluorescence Microscopy

Neuronal cell monolayers were grown on coverslips in 24-well plates and infected with ZIKV or PBS at MOI-1. Cells were fixed in 4% paraformaldehyde (PFA) and immunostained using mouse anti-dsRNA (1:1000) antibody followed by secondary antibody conjugated with Alexa Fluor 555 (Millipore, Burlington, MA, USA) as described previously [40].

2.4. NanoString nCounter® Gene Expression

Total RNA was isolated using miRNeasy Mini Kit (Qiagen, Hilden, Germany) as described previously [33]. Genomic DNA contamination was eliminated by digesting the RNA with RNase-free DNase (Ambion, Cambridge, MA, USA). RNA was quantitated using Nanodrop (Thermo Scientific, Waltham, MA, USA), and the 28S/18S RNA ratios of all RNA samples were between 1.8 and 2.0. RNA quality was analyzed using the Bioanalyzer [33,41]. For miRNA analysis, we used the nCounter® Mouse miRNA Expression Panel (NanoString, Seattle, WA, USA, Cat: CSO-MMIR15-12). Raw data was normalized using the geometric mean values of the top 100 expressed miRNA in each sample using the nSolver Analysis Software (NanoString), according to the manufacturer’s guidelines. For mRNA analysis, we utilized the nCounter® Mouse PanCancer Immune Profiling Panel to count 770 immune-related genes (NanoString, Cat: XT-CSO-MIP1-12). Raw data was normalized with a set of housekeeping genes and analyzed using the nSolver Analysis Software (NanoString), according to the manufacturer’s guidelines.

2.5. qRT-PCR

Total RNA was isolated using miRNeasy Mini Kit, and cDNA prepared using a miScript II RT Kit (Qiagen) [33,34]. qRT-PCR was performed using specific miRNA primer (Qiagen), and the miScript SYBR green PCR kit (containing Universal reverse primer) [33]. For mRNA analysis, cDNA was prepared using iScript™ cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA), and qRT-PCR was conducted as described previously [41,42]. The primer sequences used for qRT-PCR are listed in Table 1.

2.6. Measurement of Cytokines and Chemokines

The levels of cytokines and chemokines were measured in the cell supernatants by multiplex immunoassay using MILLIPLEX MAP mouse Cytokine/Chemokine magnetic panel as per manufacturer’s instructions (Millipore) [43,44].

2.7. Ingenuity Pathways Analysis (IPA)

Target prediction and pathway analysis were conducted using IPA (Ingenuity Systems Inc., Redwood City, CA, USA) as described previously [33,34,41]. Fisher’s exact test, using IPA, was used to calculate the cut-off point of significance. p < 0.05 is considered significant. We also conducted correlation pairing with mRNA expression data and significantly modulated miRNAs using IPA. For multiplex immunoassay analysis, unpaired Student’s t-test using Graph Pad was used to calculate p values.

3. Results

3.1. ZIKV Can Infect Primary Mouse Cortical Neurons

ZIKV infection of neuronal cells plays an important role in the pathogenesis of ZIKV neurological disease. Therefore, we used primary neurons for our experiments. We first determined ZIKV infection and replication kinetics in primary mouse cortical neurons by plaque assay. Mouse cortical neuron cultures were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1 and supernatants were collected at 12, 24, 48, and 72 h after infection. High ZIKV replication was observed as early as 12 h after infection. Viral titers peaked at 48 h (log 7–8 PFU/mL) followed by a slight decline at 72 h after infection (Figure 1A). Immunofluorescence staining of ZIKV-infected neurons demonstrated robust dsRNA staining in the cytoplasm. Based on a total of 5,000 cells counted in 10 independent fields, dsRNA was detected in approximately 60% of cells at 48 h after infection (Figure 1B,C).

3.2. ZIKV Infection Modulates Cellular miRNA Expression

Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1 and cell lysates were harvested at 24 and 48 h after infection. We used nCounter® Technology (NanoString) to evaluate the global miRNA expression profiles in the mock-and ZIKV-infected neurons at 24 and 48 h after infection (n = 4 per group per time point). Among 599 miRNAs present on the array, 67 and 45 miRNAs were significantly modulated at 24 and 48 h, respectively. miRNAs that were altered at least 2-fold were considered significant. While 62 miRNAs were significantly downregulated (between 2- to 4-fold), only five miRNAs were upregulated (between 2 to 18-fold) in ZIKV-infected neurons when compared to mock-infected neurons at 24 h (Table 2 and Table 3). Similarly, 40 miRNAs were significantly downregulated (between 2- to 6-fold), and only five were upregulated (between 2 to 10-fold) at 48 h (Table 2 and Table 4). As depicted in the Venn diagram (Figure 2A), three upregulated and 26 downregulated miRNAs were common in both 24 and 48 h-infected neurons. Among upregulated miRNAs; miR-155, miR-29a, and miR-29b were induced at both 24 and 48 h. miR-3471 and miR-2145 were upregulated only at 24 h, and miR-203 and miR-1902 were upregulated only at 48 h. The miRNAs with the highest induction during the course of infection were miR-155 (18.2-fold), miR-3471 (4.9-fold), and miR-203 (4.6-fold) (Table 2). Among downregulated miRNAs, miR-124-3p was the highest downregulated miRNA at both 24 (4-fold) and 48 (6.6-fold) h. Other miRNAs with high downregulation were miR-883a-3p (3.9-fold) and miR-2137 (2.9-fold). We also conducted qRT-PCR to confirm the expression changes of a selected number of differentially expressed miRNAs. Similar to the NanoString data, miR-155, miR-29a, and miR-203 were significantly upregulated, and miR-124-3p was downregulated in ZIKV-infected neurons (Figure 2B).

3.3. ZIKV Infection Results in Downregulation of miRNA Processing Enzymes

Flaviviruses have been shown to induce downregulation in the expression of cellular miRNAs by targeting miRNA processing enzymes [45,46,47,48]. Our data also demonstrated a trend toward decrease in miRNA expression in ZIKA-infected cells. Therefore, we next evaluated the expression of miRNA processing enzymes in ZIKV-infected neurons. mRNA expression levels of Dicer-1, Drosha, AGO1, and AGO2 were downregulated in ZIKV-infected neurons as compared to mock-infected neurons at both 24 and 48 h (Figure 2C). mRNA expression of DGCR8 increased slightly at 24 h followed by a decrease in the expression at 48 h.

3.4. Functional Analysis of ZIKV-Modulated miRNAs and Their Predicted Targets

Biofunctional analysis of ZIKV-modulated miRNAs and their targets revealed organismal injury and abnormalities, immunological disease, inflammatory response, neurological disease, and nervous system development and function as the top pathways in signaling pathways category (Table 5). Since ZIKV infection is associated with a wide spectrum of neurological and immunological diseases, these miRNAs may play an important role in the development of brain abnormalities following ZIKV infection.

3.5. ZIKV Infection Induces Dramatic Upregulation of Anti-Viral, Inflammatory, and Apoptotic Genes in Neurons

We next analyzed the mRNA expression of key anti-viral, inflammatory, and apoptotic genes in the ZIKV-infected neurons to examine whether differentially expressed miRNAs could regulate their target mRNAs. Mouse cortical neuron cultures were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1 and cell lysates were harvested at 24 h and 48 h after infection. We utilized the nCounter® Mouse PanCancer Immune Profiling Panel to count 770 immune-related genes. mRNA expression profiles for ZIKV-infected neurons were compared with mock-infected neurons. Only mRNAs that were altered at least 2-fold were considered significant. Figure 3A demonstrates total numbers of up- and downregulated differentially expressed mRNAs in ZIKV-infected neurons at 24 and 48 h. The number of differentially expressed mRNAs was higher at 48 h after infection, which correlates with significantly higher viral load observed at 48 h as compared to 24 h (Figure 1A). At 24 h, 65 mRNAs were upregulated with fold change values ranging from 2.0 to 26 (Table 6). mRNAs were not downregulated at 24 h. At 48 h, 116 mRNAs were upregulated with fold change values ranging from 2.0 to 244 (Table 6 and Table 7) and 12 mRNAs were downregulated (Table 8). 62 upregulated mRNAs were common in both 24 and 48 h- infected neurons (Table 6). MMP9, SOCS3 and IFI27 were upregulated only at 24 h. 54 mRNAs were upregulated only at 48 h (Table 7).
Genes associated with virus sensing and type 1 IFN signaling were the most upregulated genes after ZIKV infection. ZIKV infection also induced a strong upregulation of multiple cytokines and chemokines in the neurons. Most of the chemokines and cytokines were significantly upregulated after ZIKV infection, including CXCL10, CCL5, CCL7, CCL2, CXCL1, CCL12, CXCL9, CXCL2, CCL4, CXCL11, CXCL13, CCL11, CXCL5, IL6, IL7, PTGS2, and LIF. Several genes involved in cell death (CASP1, TNFSF10, RIPK2, and BID) were also activated upon ZIKV infection (Table 6 and Table 7). We also validated the expression of selected upregulated mRNAs using qRT-PCR. Similar to the NanoString data, IFIT1, IFIT3, IL6, and Caspase1 were significantly upregulated in the ZIKV-infected neurons as compared to mock-infected neurons at both 24 and 48 h (Figure 3B). To examine that this increase in mRNA expression also lead to increased protein levels, we measured protein levels of key chemokines and cytokines in the cell culture supernatants using multiplex immunoassay. Similar to the mRNA expression, protein levels of the key chemokines such as CCL2, CCL4, CCL5, CCL11, CXCL1, CXCL2, CXCL9, and CXCL10, and cytokines such as IL6 and LIF were significantly increased in ZIKV-infected neurons as compared to mock-infected neurons (Figure 4).

3.6. Functional Analysis of ZIKV-Modulated mRNAs

To investigate the biological interactions of differentially expressed mRNAs and identify important functional networks, significantly modulated mRNAs were imported into the IPA tool. The highest activated networks (high z-score) were identified using IPA. Figure 5 depicts the top 10 activated canonical pathways after ZIKV infection. The topmost activated canonical pathway after ZIKV infection was ‘Neuroinflammation Signaling’. In addition, key players in inflammation and innate immunity, such as ‘Activation of IRFs by Cytosolic Pattern Recognition Receptors (PRR)’, ‘Role of PRR in Recognition of Viruses and Bacteria’, and ‘IFN Signaling’ were also highly activated (positive z-score). To further understand the role of these differentially activated canonical pathways, we generated the network maps of ‘PRR in Recognition of Viruses and Bacteria’ and ‘IFN Signaling’.

3.6.1. Pattern Recognition Receptors

Our data demonstrate that ZIKV infection induces the expression of all three major PRR; retinoic-acid-inducible gene-I (RIG-I)-like receptors (RLR), toll-like receptors (TLR), and the nucleotide oligomerization domain (Nod)-like receptors (NLR) (Figure 6A). The RLR are a family of cytosolic RNA helicase proteins comprised of three members: RIG-I, myeloma differentiation antigen 5 (MDA5), and LGP2 [49]. Both RIG-I and MDA5 were upregulated after ZIKV infection. The importance of the RLR signaling pathway in protection against flaviviruses has been validated by several studies in vivo [49,50,51]. The TLR family is composed of more than 10 members, with each acting as a sensor of conserved microbial component, that drive the induction of immune response [52]. Our data show that ZIKV infection induces the expression of TLR2, TLR3, and adaptor molecule MYD88. NLR are soluble or cytosolic receptors in the mammalian cell cytoplasm [53]. Activation of NOD1 and NLRC5 was observed following ZIKV infection. In addition to RLR, TLR, and NLR; PKR and OAS are classes of IFN-inducible PRR that can recognize dsRNA and restrict a number of viruses. Our data demonstrate the activation of OAS after infection with ZIKV. Furthermore, we observed increased levels of type 1 interferons and several pro-inflammatory mediators after ZIKV infection in neurons, which correlate with PRR activation in these cells.

3.6.2. IFN Signaling

ZIKV infection induced strong upregulation of genes associated with IFN signaling such as IFNα, IFNβ, STAT1, and STAT2 (Figure 6B). Interferon-stimulated genes (ISG) such as MX2, IFIT1, IFIT3, IFITM1, IFI35, and RSAD2 (Viperin) were significantly upregulated after ZIKV infection. The IFN response is central to the innate defense mechanisms of the host against flavivirus infection. The paracrine and autocrine secretion of IFN creates an anti-viral state by inducing several genes including ISG [49,54].

3.7. Network Analysis of Expression of miRNAs and mRNAs From ZIKV-Infected Neurons

We next sought to determine whether mRNA expression changes might be influenced by the differential expression of cellular miRNAs during ZIKV infection. To analyze the direct and indirect miRNA–mRNA interactions, we conducted IPA expression pairing analysis with mRNA expression data and significantly modulated miRNAs. Several miRNAs were found to directly or indirectly target multiple mRNAs analyzed in our study and demonstrated an inverse correlation with mRNA expression induced by ZIKV. These targets are predicted by TargetScan. miR-124-3p was the highest downregulated miRNA. Our data demonstrated increased expression of all the predicted targets of miR-124-3p including IL7, CCL2, LITAF, IRF1, and SBNO2 in ZIKV-infected neurons (Figure 7A). Other known targets of miR-124 includes SOCS5, TLR6, STAT3, TNF, and NF-kB [55]. Similarly, miR-654-3p was downregulated and its predicted targets—CD69, FLT3LG, IFIT1, TLR2, ZBP1, LITAF—were upregulated in ZIKV-infected neurons (Figure 7B). These targets involve genes belonging to anti-viral and inflammatory response signaling pathways. Furthermore, downregulation of miR-331-5p and miR-509-5p was inversely correlated to their predicted targets IL7, CD274, HLA-DRB5, XAF1, IL13RA1, TAP1, and CD80 (Figure 7C,D). Our data further indicate that miR-335-3p may play an important role in regulating inflammatory and apoptotic genes in the neurons following ZIKV infection. miR-335-3p targets Caspase1, CCL5, CXCL3, PTPRC, and GBP5. These genes play significant roles in mediating inflammation and cell death (Figure 7E). We also show an increase in the protein levels of the target genes such as CCL2, CCL5, and CXCL10 (Figure 4). It is interesting to see the correlation between ZIKV-modulated miRNAs and target genes at both mRNA and protein level, which is consistent with the miRNA–mRNA–protein triad and demonstrate the functional importance of our results.

4. Discussion

In this study, we used nCounter technology to identify miRNAs and mRNAs modulated by ZIKV infection in neurons. Our data demonstrate that ZIKV infection causes global downregulation of miRNAs with only few upregulated miRNAs. ZIKV infection also results in downregulation of miRNA processing enzymes. Upregulated miRNAs including miR-155, miR-203, and miR-29a have been previously shown to play critical role in flavivirus infection and anti-viral immunity [22,56,57,58,59,60]. ZIKV infection in the neurons significantly induced the expression of anti-viral, inflammatory and apoptotic genes. Biofunctional analysis revealed that ZIKV-modulated miRNAs and their target genes regulate the pathways related to neurological development and neuroinflammatory responses. Furthermore, our data show an inverse correlation between ZIKV-modulated miRNAs and target host immune mRNAs induced by ZIKV infection.
Flaviviruses have been shown to induce downregulation in the expression of cellular miRNAs by targeting miRNA processing enzymes [33,45,46,47,48,61]. It has been reported that depletion of Dicer and Drosha by siRNA-mediated silencing results in increase in flavivirus replication [46,47,61]. We also observed a trend towards decrease in miRNAs expression in ZIKA-infected cells. Similar to our data, ZIKV infection in human astrocytes induces global downregulation of miRNAs [62]. It has been demonstrated that ZIKV infection in Aedes aegypti mosquitoes downregulates expression of host miRNAs [63]. In this study, we also show that ZIKV infection results in downregulation of miRNA processing enzymes. Downregulation of genes involved in miRNA processing enzymes including Dicer-1 was also reported in ZIKV-infected astrocytes [62]. Similarly, ZIKV infection in HepG2 cells led to a decrease in DGCR8, Ago1, and Ago3 expression [64]. It is known that non-coding, subgenomic RNA (sfRNA) of WNV and DENV can interfere with miRNA biogenesis pathway. WNV sfRNA is processed by Dicer and suppresses RNAi [33,45,46,47,48,61]. Future experiments are warranted to elucidate whether this sfRNA function holds for all flaviviruses including ZIKV. Since miRNAs negatively regulate mRNA expression, decrease in cellular miRNA expression results in increase in apoptotic and inflammatory genes associated with flavivirus infection [27,30,31,33,41]. Our data also demonstrate that ZIKV infection is associated with dramatic upregulation of several anti-viral, apoptotic and inflammatory genes.
We observed significant upregulation of miR-155, miR-203, miR-29a, and miR-29b in ZIKV-infected neurons. These miRNAs have been shown to play an important role in viral infection [19,20,21,65,66]. miR-155 is multifunctional and widely reported to modulate different stages of innate immune response during inflammation and infection [21,22,56]. miR-155 not only modulates TLR-mediated innate immune response, but also targets complement regulatory proteins and facilitate complement activation [56,57]. This phenomenon is critical to eliminate the virus from infected cells. Several published studies have demonstrated the essential role of miR-155 in viral infections caused by Epstein–Barr, Borna disease, and reticuloendotheliosis viruses [19,20,56,65]. For example, overexpression of miR-155 significantly suppressed human HIV infection in activated macrophages [21]. Similarly, miR-155 suppresses JEV replication in microglial cells and regulates JEV-induced inflammatory response in mice brain [22,67]. Studies have shown that miR-203 is also involved in regulating the anti-viral immune response [59,60]. It has been demonstrated that miR-29b regulates JEV-induced neuroinflammation [58]. Since these miRNAs are known to have a role in anti-viral immunity, upregulation of these miRNAs might play a role in controlling ZIKV infection and associated neurotoxicity. A recent study identified upregulation of miR-30e-3p, miR-30e-5p, and miR-17-5p in ZIKV infection of human astrocytes [62]. However, we did not observe upregulation of these miRNAs in ZIKV-infected mouse neurons. This could be due to the difference in the neural cell type and species studied.
In our study, miRNAs with the highest downregulation were miR-124-3p, miR-883a-5p and miR-2137. miR-124 is the most abundant miRNA in the brain and affects a broad spectrum of biological functions in the CNS [55,68,69,70]. miR-124 has been reported to participate in chronic stress, neurodegeneration, alcohol/cocaine neuroadaptation, synapse morphology, neurotransmission, long-term potentiation, neurodevelopment, myeloid cell function, and hematopoiesis. In mammalian neurons, miR-124 suppresses the levels of many non-neural genes, which contributes to the acquisition and maintenance of neuronal identity [69]. Furthermore, when miR-124 is aberrantly expressed, it contributes to pathological conditions involving the CNS [68,69,70]. It has also been shown to be useful as a diagnostic and prognostic indicator of CNS disorders, such as brain tumors and stroke [68]. miR-883a has been demonstrated to be involved in inflammatory pathways, and target genes belonging to the TLR signaling pathway as well as the VEGF and chemokine signaling pathways [71]. miR-2137 is involved in the pathogenesis of traumatic brain injury [72]. IPA analysis also revealed that individual miRNA including miR-124a-3p, miR-654-3p, miR-331-5p, miR-335-3p, and miR-509-5p can modulate multiple upstream regulatory genes. There was an inverse correlation between these miRNAs and target host immune mRNAs induced by ZIKV infection. Since these miRNAs are known to have a role in brain injury and inflammatory responses, downregulation of these miRNAs following ZIKV infection possibly resulted in upregulation of neuroinflammatory and apoptotic genes. However, further experiments are required to validate their biological functions in ZIKV infection.
This study has few limitations. First, we only analyzed the mRNA expression of genes primarily involved in immune response, which introduce a bias in the analysis towards anti-viral response. Second, additional experimental studies are warranted to validate the interplay between miRNA, mRNA, and ZIKV replication, which includes knockdown of individual miRNA and its effect on the expression levels of its target mRNAs and virus titers.
In conclusion, this is the first study to evaluate the modulation of miRNAs following ZIKV infection in neurons. ZIKV-modulated miRNAs in neurons are also known to play a role in the pathogenesis of other flaviviruses and anti-viral immune response. In this study, the utility of the nCounter system enabling rapid miRNA and mRNA expression analysis was also demonstrated. Collectively, these data suggest that miRNAs regulate downstream gene expression, important in ZIKV disease pathogenesis, and can be targeted in the future to develop therapeutics for the management of ZIKV neurological disease.

Author Contributions

Conceptualization, M.K. and V.R.N.; Methodology, M.K., K.A., F.A., and K.K.; Software, M.K., F.A., and K.A.; Validation, F.A. and K.A.; Formal analysis, M.K., K.A., F.A., and V.R.N.; Writing—original draft preparation, M.K. and K.A.; Writing—review and editing, M.K., K.A., F.A., K.K., and V.R.N.; Funding acquisition, M.K.

Funding

This work was supported by a grant (P30GM114737) from the Centers of Biomedical Research Excellence, National Institute of General Medical Sciences, grant (R21NS099838) from National Institute of Neurological Disorders and Stroke, grant (R21OD024896) from the Office of the Director, National Institutes of Health, and grant (18CON-90812) from Hawaii Community Foundation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Musso, D.; Gubler, D.J. Zika virus. Clin. Microbiol. Rev. 2016, 29, 487–524. [Google Scholar] [CrossRef] [PubMed]
  2. Klase, Z.A.; Khakhina, S.; Schneider Ade, B.; Callahan, M.V.; Glasspool-Malone, J.; Malone, R. Zika fetal neuropathogenesis: Etiology of a viral syndrome. PLoS Negl. Trop. Dis. 2016, 10, e0004877. [Google Scholar] [CrossRef] [PubMed]
  3. Cao-Lormeau, V.M.; Blake, A.; Mons, S.; Lastere, S.; Roche, C.; Vanhomwegen, J.; Dub, T.; Baudouin, L.; Teissier, A.; Larre, P.; et al. Guillain-Barré syndrome outbreak associated with zika virus infection in french polynesia: A case-control study. Lancet 2016, 387, 1531–1539. [Google Scholar] [CrossRef]
  4. 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] [PubMed]
  5. Anfasa, F.; Siegers, J.Y.; van der Kroeg, M.; Mumtaz, N.; Stalin Raj, V.; 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. [Google Scholar] [CrossRef] [PubMed]
  6. Gaburro, J.; Bhatti, A.; Sundaramoorthy, V.; Dearnley, M.; Green, D.; Nahavandi, S.; Paradkar, P.N.; Duchemin, J.B. Zika virus-induced hyper excitation precedes death of mouse primary neuron. Virol. J. 2018, 15, 79. [Google Scholar] [CrossRef]
  7. Rosenfeld, A.B.; Doobin, D.J.; Warren, A.L.; Racaniello, V.R.; Vallee, R.B. Replication of early and recent zika virus isolates throughout mouse brain development. Proc. Natl. Acad. Sci. USA 2017, 114, 12273–12278. [Google Scholar] [CrossRef]
  8. Lazear, H.M.; Govero, J.; Smith, A.M.; Platt, D.J.; Fernandez, E.; Miner, J.J.; Diamond, M.S. A mouse model of zika virus pathogenesis. Cell host microbe 2016, 19, 720–730. [Google Scholar] [CrossRef]
  9. Rossi, S.L.; Tesh, R.B.; Azar, S.R.; Muruato, A.E.; Hanley, K.A.; Auguste, A.J.; Langsjoen, R.M.; Paessler, S.; Vasilakis, N.; Weaver, S.C. Characterization of a novel murine model to study zika virus. Am. J. Trop Med. Hyg. 2016, 94, 1362–1369. [Google Scholar] [CrossRef]
  10. Duggal, N.K.; Ritter, J.M.; McDonald, E.M.; Romo, H.; Guirakhoo, F.; Davis, B.S.; Chang, G.J.; Brault, A.C. Differential neurovirulence of african and asian genotype zika virus isolates in outbred immunocompetent mice. Am. J. Trop Med. Hyg. 2017, 97, 1410–1417. [Google Scholar] [CrossRef]
  11. Manangeeswaran, M.; Ireland, D.D.; Verthelyi, D. Zika (prvabc59) infection is associated with t cell infiltration and neurodegeneration in cns of immunocompetent neonatal c57bl/6 mice. PLoS pathog. 2016, 12, e1006004. [Google Scholar] [CrossRef] [PubMed]
  12. Van den Pol, A.N.; Mao, G.; Yang, Y.; Ornaghi, S.; Davis, J.N. Zika virus targeting in the developing brain. J. Neurosci. 2017, 37, 2161–2175. [Google Scholar] [CrossRef] [PubMed]
  13. Fernandes, N.C.; Nogueira, J.S.; Ressio, R.A.; Cirqueira, C.S.; Kimura, L.M.; Fernandes, K.R.; Cunha, M.S.; Souza, R.P.; Guerra, J.M. Experimental zika virus infection induces spinal cord injury and encephalitis in newborn swiss mice. Exp. Toxicol. Pathol. 2017, 69, 63–71. [Google Scholar] [CrossRef] [PubMed]
  14. Huang, W.C.; Abraham, R.; Shim, B.S.; Choe, H.; Page, D.T. Zika virus infection during the period of maximal brain growth causes microcephaly and corticospinal neuron apoptosis in wild type mice. Sci. Rep. 2016, 6, 34793. [Google Scholar] [CrossRef] [PubMed]
  15. Garofalo, M.; Condorelli, G.; Croce, C.M. Micrornas in diseases and drug response. Curr. Opin. Pharmacol. 2008, 8, 661–667. [Google Scholar] [CrossRef] [PubMed]
  16. Junn, E.; Mouradian, M.M. Micrornas in neurodegenerative diseases and their therapeutic potential. Pharmacol. Ther. 2012, 133, 142–150. [Google Scholar] [CrossRef] [PubMed]
  17. Pauley, K.M.; Chan, E.K. Micrornas and their emerging roles in immunology. Ann. NY Acad. Sci. 2008, 1143, 226–239. [Google Scholar] [CrossRef]
  18. O’Connell, R.M.; Rao, D.S.; Baltimore, D. Microrna regulation of inflammatory responses. Annu. Rev. Immunol. 2012, 30, 295–312. [Google Scholar] [CrossRef]
  19. Zhai, A.; Qian, J.; Kao, W.; Li, A.; Li, Y.; He, J.; Zhang, Q.; Song, W.; Fu, Y.; Wu, J.; et al. Borna disease virus encoded phosphoprotein inhibits host innate immunity by regulating mir-155. Antiviral Res. 2013, 98, 66–75. [Google Scholar] [CrossRef]
  20. Bolisetty, M.T.; Dy, G.; Tam, W.; Beemon, K.L. Reticuloendotheliosis virus strain t induces Mir-155, which targets JARID2 and promotes cell survival. J. Virol. 2009, 83, 12009–12017. [Google Scholar] [CrossRef]
  21. Swaminathan, G.; Rossi, F.; Sierra, L.J.; Gupta, A.; Navas-Martin, S.; Martin-Garcia, J. A role for microrna-155 modulation in the anti-HIV-1 effects of toll-like receptor 3 stimulation in macrophages. PLoS Pathog. 2012, 8, e1002937. [Google Scholar] [CrossRef] [PubMed]
  22. Pareek, S.; Roy, S.; Kumari, B.; Jain, P.; Banerjee, A.; Vrati, S. Mir-155 induction in microglial cells suppresses japanese encephalitis virus replication and negatively modulates innate immune responses. J. Neuroinflammation 2014, 11, 97. [Google Scholar] [CrossRef] [PubMed]
  23. Fabozzi, G.; Nabel, C.S.; Dolan, M.A.; Sullivan, N.J. Ebolavirus proteins suppress the effects of small interfering rna by direct interaction with the mammalian rna interference pathway. J. Virol. 2011, 85, 2512–2523. [Google Scholar] [CrossRef]
  24. Huang, J.; Wang, F.; Argyris, E.; Chen, K.; Liang, Z.; Tian, H.; Huang, W.; Squires, K.; Verlinghieri, G.; Zhang, H. Cellular micrornas contribute to hiv-1 latency in resting primary CD4+ T lymphocytes. Nat. Med. 2007, 13, 1241–1247. [Google Scholar] [CrossRef] [PubMed]
  25. Jopling, C.L.; Yi, M.; Lancaster, A.M.; Lemon, S.M.; Sarnow, P. Modulation of hepatitis c virus rna abundance by a liver-specific microrna. Science 2005, 309, 1577–1581. [Google Scholar] [CrossRef] [PubMed]
  26. Pedersen, I.M.; Cheng, G.; Wieland, S.; Volinia, S.; Croce, C.M.; Chisari, F.V.; David, M. Interferon modulation of cellular micrornas as an antiviral mechanism. Nature 2007, 449, 919–922. [Google Scholar] [CrossRef] [PubMed]
  27. Bavia, L.; Mosimann, A.L.; Aoki, M.N.; Duarte Dos Santos, C.N. A glance at subgenomic flavivirus rnas and micrornas in flavivirus infections. Virol. J. 2016, 13, 84. [Google Scholar] [CrossRef] [PubMed]
  28. 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]
  29. Escalera-Cueto, M.; Medina-Martinez, I.; del Angel, R.M.; Berumen-Campos, J.; Gutierrez-Escolano, A.L.; Yocupicio-Monroy, M. Let-7c overexpression inhibits dengue virus replication in human hepatoma Huh-7 cells. Virus Res. 2015, 196, 105–112. [Google Scholar] [CrossRef]
  30. Slonchak, A.; Shannon, R.P.; Pali, G.; Khromykh, A.A. Human microrna MiR-532-5p exhibits antiviral activity against west nile virus via suppression of host genes SESTD1 and TAB3 required for virus replication. J. virol. 2015, 90, 2388–2402. [Google Scholar] [CrossRef]
  31. Smith, J.L.; Grey, F.E.; Uhrlaub, J.L.; Nikolich-Zugich, J.; Hirsch, A.J. Induction of the cellular microrna, hs_154, by west nile virus contributes to virus-mediated apoptosis through repression of antiapoptotic factors. J. Virol. 2012, 86, 5278–5287. [Google Scholar] [CrossRef] [PubMed]
  32. Heiss, B.L.; Maximova, O.A.; Pletnev, A.G. Insertion of microrna targets into the flavivirus genome alters its highly neurovirulent phenotype. J. virol. 2011, 85, 1464–1472. [Google Scholar] [CrossRef] [PubMed]
  33. Kumar, M.; Nerurkar, V.R. Integrated analysis of micrornas and their disease related targets in the brain of mice infected with west nile virus. Virology 2014, 452–453, 143–151. [Google Scholar] [CrossRef] [PubMed]
  34. Shin, O.S.; Kumar, M.; Yanagihara, R.; Song, J.W. Hantaviruses induce cell type- and viral species-specific host microrna expression signatures. Virology 2013, 446, 217–224. [Google Scholar] [CrossRef] [PubMed]
  35. Geiss, G.K.; Bumgarner, R.E.; Birditt, B.; Dahl, T.; Dowidar, N.; Dunaway, D.L.; Fell, H.P.; Ferree, S.; George, R.D.; Grogan, T.; et al. Direct multiplexed measurement of gene expression with color-coded probe pairs. Nat. Biotechnol. 2008, 26, 317–325. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Wakayama, I.; Song, K.J.; Nerurkar, V.R.; Yoshida, S.; Garruto, R.M. Slow dendritic transport of dissociated mouse hippocampal neurons exposed to aluminum. Brain Res. 1997, 748, 237–240. [Google Scholar] [CrossRef]
  37. Forest, K.H.; Alfulaij, N.; Arora, K.; Taketa, R.; Sherrin, T.; Todorovic, C.; Lawrence, J.L.M.; Yoshikawa, G.T.; Ng, H.L.; Hruby, V.J.; et al. Protection against beta-amyloid neurotoxicity by a non-toxic endogenous n-terminal beta-amyloid fragment and its active hexapeptide core sequence. J. Neurochem. 2018, 144, 201–217. [Google Scholar] [CrossRef] [PubMed]
  38. Kim, J.A.; Seong, R.K.; Kumar, M.; Shin, O.S. Favipiravir and ribavirin inhibit replication of asian and african strains of zika virus in different cell models. Viruses 2018, 10, 72. [Google Scholar] [CrossRef] [PubMed]
  39. Kumar, M.; Krause, K.K.; Azouz, F.; Nakano, E.; Nerurkar, V.R. A guinea pig model of zika virus infection. Virol. J. 2017, 14, 75. [Google Scholar] [CrossRef] [PubMed]
  40. Verma, S.; Lo, Y.; Chapagain, M.; Lum, S.; Kumar, M.; Gurjav, U.; Luo, H.; Nakatsuka, A.; Nerurkar, V.R. West nile virus infection modulates human brain microvascular endothelial cells tight junction proteins and cell adhesion molecules: Transmigration across the in vitro blood-brain barrier. Virology 2009, 385, 425–433. [Google Scholar] [CrossRef]
  41. Kumar, M.; Belcaid, M.; Nerurkar, V.R. Identification of host genes leading to west nile virus encephalitis in mice brain using RNA-Seq analysis. Sci. Rep. 2016, 6, 26350. [Google Scholar] [CrossRef] [PubMed]
  42. Kumar, M.; Verma, S.; Nerurkar, V.R. Pro-inflammatory cytokines derived from west nile virus (WNV)-infected SK-N-SH cells mediate neuroinflammatory markers and neuronal death. J. Neuroinflammation 2010, 7, 73. [Google Scholar] [CrossRef] [PubMed]
  43. Kumar, M.; Roe, K.; Nerurkar, P.V.; Namekar, M.; Orillo, B.; Verma, S.; Nerurkar, V.R. Impaired virus clearance, compromised immune response and increased mortality in type 2 diabetic mice infected with west nile virus. PLoS ONE 2012, 7, e44682. [Google Scholar] [CrossRef] [PubMed]
  44. Kumar, M.; Roe, K.; Orillo, B.; Muruve, D.A.; Nerurkar, V.R.; Gale, M., Jr.; Verma, S. Inflammasome adaptor protein apoptosis-associated speck-like protein containing card (ASC) is critical for the immune response and survival in west nile virus encephalitis. J. Virol. 2013, 87, 3655–3667. [Google Scholar] [CrossRef] [PubMed]
  45. Pijlman, G.P. Flavivirus rnai suppression: Decoding non-coding rna. Curr. Opin. Virol. 2014, 7, 55–60. [Google Scholar] [CrossRef]
  46. Schnettler, E.; Sterken, M.G.; Leung, J.Y.; Metz, S.W.; Geertsema, C.; Goldbach, R.W.; Vlak, J.M.; Kohl, A.; Khromykh, A.A.; Pijlman, G.P. Noncoding flavivirus rna displays rna interference suppressor activity in insect and mammalian cells. J. Virol. 2012, 86, 13486–13500. [Google Scholar] [CrossRef]
  47. Kakumani, P.K.; Ponia, S.S.; S., R.K.; Sood, V.; Chinnappan, M.; Banerjea, A.C.; Medigeshi, G.R.; Malhotra, P.; Mukherjee, S.K.; Bhatnagar, R.K. Role of RNA interference (RNAi) in dengue virus replication and identification of ns4b as an rnai suppressor. J. Virol. 2013, 87, 8870–8883. [Google Scholar] [CrossRef]
  48. Moon, S.L.; Dodd, B.J.; Brackney, D.E.; Wilusz, C.J.; Ebel, G.D.; Wilusz, J. Flavivirus sfrna suppresses antiviral rna interference in cultured cells and mosquitoes and directly interacts with the rnai machinery. Virology 2015, 485, 322–329. [Google Scholar] [CrossRef]
  49. Fredericksen, B.L. The neuroimmune response to west nile virus. J. Neurovirol. 2014, 20, 113–121. [Google Scholar] [CrossRef]
  50. Fredericksen, B.L.; Keller, B.C.; Fornek, J.; Katze, M.G.; Gale, M., Jr. Establishment and maintenance of the innate antiviral response to west nile virus involves both rig-i and mda5 signaling through IPS-1. J. Virol. 2008, 82, 609–616. [Google Scholar] [CrossRef]
  51. Errett, J.S.; Suthar, M.S.; McMillan, A.; Diamond, M.S.; Gale, M., Jr. The essential, nonredundant roles of rig-i and mda5 in detecting and controlling west nile virus infection. J. Virol. 2013, 87, 11416–11425. [Google Scholar] [CrossRef] [PubMed]
  52. Akira, S.; Uematsu, S.; Takeuchi, O. Pathogen recognition and innate immunity. Cell 2006, 124, 783–801. [Google Scholar] [CrossRef] [PubMed]
  53. Kaushik, D.K.; Gupta, M.; Kumawat, K.L.; Basu, A. Nlrp3 inflammasome: Key mediator of neuroinflammation in murine japanese encephalitis. PLoS ONE 2012, 7, e32270. [Google Scholar] [CrossRef] [PubMed]
  54. Diamond, M.S.; Gale, M., Jr. Cell-intrinsic innate immune control of west nile virus infection. Trends Immunol. 2012, 33, 522–530. [Google Scholar] [CrossRef] [PubMed]
  55. Qin, Z.; Wang, P.Y.; Su, D.F.; Liu, X. Mirna-124 in immune system and immune disorders. Front. Immunol. 2016, 7, 406. [Google Scholar] [CrossRef] [PubMed]
  56. Jiang, M.; Broering, R.; Trippler, M.; Wu, J.; Zhang, E.; Zhang, X.; Gerken, G.; Lu, M.; Schlaak, J.F. Microrna-155 controls toll-like receptor 3- and hepatitis C virus-induced immune responses in the liver. J. Viral. Hepat. 2014, 21, 99–110. [Google Scholar] [CrossRef] [PubMed]
  57. Dickey, L.L.; Hanley, T.M.; Huffaker, T.B.; Ramstead, A.G.; O’Connell, R.M.; Lane, T.E. MicroRNA 155 and viral-induced neuroinflammation. J. Neuroimmunol. 2017, 308, 17–24. [Google Scholar] [CrossRef]
  58. Thounaojam, M.C.; Kaushik, D.K.; Kundu, K.; Basu, A. MicroRNA-29b modulates japanese encephalitis virus-induced microglia activation by targeting tumor necrosis factor alpha-induced protein 3. J. Neurochem. 2014, 129, 143–154. [Google Scholar] [CrossRef] [PubMed]
  59. Buggele, W.A.; Horvath, C.M. Microrna profiling of sendai virus-infected a549 cells identifies MiR-203 as an interferon-inducible regulator of IFIT1/ISG56. J. Virol. 2013, 87, 9260–9270. [Google Scholar] [CrossRef] [PubMed]
  60. Zhang, S.; Li, J.; Li, J.; Yang, Y.; Kang, X.; Li, Y.; Wu, X.; Zhu, Q.; Zhou, Y.; Hu, Y. Upregulation of microrna-203 in influenza A virus infection inhibits viral replication by targeting dr1. Sci. Rep. 2018, 8, 6797. [Google Scholar] [CrossRef]
  61. Song, M.S.; Rossi, J.J. Molecular mechanisms of dicer: Endonuclease and enzymatic activity. Biochem. J. 2017, 474, 1603–1618. [Google Scholar] [CrossRef] [PubMed]
  62. 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] [PubMed]
  63. Saldana, M.A.; Etebari, K.; Hart, C.E.; Widen, S.G.; Wood, T.G.; Thangamani, S.; Asgari, S.; Hughes, G.L. Zika virus alters the microrna expression profile and elicits an rnai response in aedes aegypti mosquitoes. PLoS Negl. Trop Dis. 2017, 11, e0005760. [Google Scholar] [CrossRef] [PubMed]
  64. Ferreira, R.N.; Holanda, G.M.; Pinto Silva, E.V.; Casseb, S.M.M.; Melo, K.F.L.; Carvalho, C.A.M.; Lima, J.A.; Vasconcelos, P.F.C.; Cruz, A.C.R. Zika virus alters the expression profile of microrna-related genes in liver, lung, and kidney cell lineages. Viral. Immunol. 2018, 31, 583–588. [Google Scholar] [CrossRef] [PubMed]
  65. Lu, F.; Weidmer, A.; Liu, C.G.; Volinia, S.; Croce, C.M.; Lieberman, P.M. Epstein–Barr virus-induced mir-155 attenuates NF-KB signaling and stabilizes latent virus persistence. J.Virol. 2008, 82, 10436–10443. [Google Scholar] [CrossRef] [PubMed]
  66. Tili, E.; Croce, C.M.; Michaille, J.J. Mir-155: On the crosstalk between inflammation and cancer. Int. Rev. Immunol. 2009, 28, 264–284. [Google Scholar] [CrossRef] [PubMed]
  67. Thounaojam, M.C.; Kundu, K.; Kaushik, D.K.; Swaroop, S.; Mahadevan, A.; Shankar, S.K.; Basu, A. MicroRNA 155 regulates Japanese encephalitis virus-induced inflammatory response by targeting src homology 2-containing inositol phosphatase 1. J. Virol. 2014, 88, 4798–4810. [Google Scholar] [CrossRef] [PubMed]
  68. Sun, Y.; Luo, Z.M.; Guo, X.M.; Su, D.F.; Liu, X. An updated role of microrna-124 in central nervous system disorders: A review. Front. Cell Neurosci. 2015, 9, 193. [Google Scholar] [CrossRef] [PubMed]
  69. Papagiannakopoulos, T.; Kosik, K.S. MicroRNA-124: Micromanager of neurogenesis. Cell Stem Cell 2009, 4, 375–376. [Google Scholar] [CrossRef]
  70. Roy, B.; Dunbar, M.; Shelton, R.C.; Dwivedi, Y. Identification of microRNA-124-3p as a putative epigenetic signature of major depressive disorder. Neuropsychopharmacology 2017, 42, 864–875. [Google Scholar] [CrossRef]
  71. Ge, Q.; Gerard, J.; Noel, L.; Scroyen, I.; Brichard, S.M. MicroRNAs regulated by adiponectin as novel targets for controlling adipose tissue inflammation. Endocrinology 2012, 153, 5285–5296. [Google Scholar] [CrossRef] [PubMed]
  72. Meissner, L.; Gallozzi, M.; Balbi, M.; Schwarzmaier, S.; Tiedt, S.; Terpolilli, N.A.; Plesnila, N. Temporal profile of microRNA expression in contused cortex after traumatic brain injury in mice. J. Neurotrauma. 2016, 33, 713–720. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Zika virus (ZIKV) infection of the primary mouse neurons. Mouse cortical neuron cultures were prepared from one-day old pups. Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at multiplicity of infection (MOI)-1. (A) ZIKV titers in culture supernatant were determined by plaque assay. Viral titers are expressed as plaque forming units (PFU)/mL of supernatant. Data represents the mean ± SEM. Neurons grown and fixed on coverslips at 48 h after infection were stained with anti-dsRNA antibody (red) and counterstained with DAPI (blue). (B) Mock-infected cells. 20× magnification. (C) ZIKV-infected cells demonstrate robust virus staining in the cytoplasm. 20× magnification.
Figure 1. Zika virus (ZIKV) infection of the primary mouse neurons. Mouse cortical neuron cultures were prepared from one-day old pups. Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at multiplicity of infection (MOI)-1. (A) ZIKV titers in culture supernatant were determined by plaque assay. Viral titers are expressed as plaque forming units (PFU)/mL of supernatant. Data represents the mean ± SEM. Neurons grown and fixed on coverslips at 48 h after infection were stained with anti-dsRNA antibody (red) and counterstained with DAPI (blue). (B) Mock-infected cells. 20× magnification. (C) ZIKV-infected cells demonstrate robust virus staining in the cytoplasm. 20× magnification.
Viruses 11 00162 g001
Figure 2. ZIKV infection of the primary mouse neurons causes changes in cellular miRNA expression. Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1. (A) Venn diagram showing the number of differentially expressed miRNAs at 24 and 48 h after infection. Sets of upregulated miRNAs are represented by upward red arrows and sets of downregulated miRNAs are represented by downward green arrows. Pairs of arrows in the intersection refer to the number of miRNAs upregulated (double red arrows) or down regulated (double green arrows) at both 24 and 48 h after infection. (B) qRT-PCR was conducted on RNA extracted from mock and ZIKV-infected neurons to determine fold-change in miR-155, miR-203, miR-29a, and miR-124-3p expression. Changes in the levels of each miRNA were first normalized to the snoRNA and then the fold-change in ZIKV-infected cells was calculated in comparison to corresponding mock-infected cells. Data represents the mean ± SEM. (C) qRT-PCR was conducted on RNA extracted from mock and ZIKV-infected neurons to determine fold-change in Dicer-1, Drosha, DGCR8, AGO1, and AGO2 expression. Changes in the levels of each mRNA were first normalized to the β-actin and then the fold-change in ZIKV-infected cells was calculated in comparison to corresponding mock-infected cells. Data represents the mean ± SEM.
Figure 2. ZIKV infection of the primary mouse neurons causes changes in cellular miRNA expression. Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1. (A) Venn diagram showing the number of differentially expressed miRNAs at 24 and 48 h after infection. Sets of upregulated miRNAs are represented by upward red arrows and sets of downregulated miRNAs are represented by downward green arrows. Pairs of arrows in the intersection refer to the number of miRNAs upregulated (double red arrows) or down regulated (double green arrows) at both 24 and 48 h after infection. (B) qRT-PCR was conducted on RNA extracted from mock and ZIKV-infected neurons to determine fold-change in miR-155, miR-203, miR-29a, and miR-124-3p expression. Changes in the levels of each miRNA were first normalized to the snoRNA and then the fold-change in ZIKV-infected cells was calculated in comparison to corresponding mock-infected cells. Data represents the mean ± SEM. (C) qRT-PCR was conducted on RNA extracted from mock and ZIKV-infected neurons to determine fold-change in Dicer-1, Drosha, DGCR8, AGO1, and AGO2 expression. Changes in the levels of each mRNA were first normalized to the β-actin and then the fold-change in ZIKV-infected cells was calculated in comparison to corresponding mock-infected cells. Data represents the mean ± SEM.
Viruses 11 00162 g002
Figure 3. ZIKV infection of the primary mouse neurons causes changes in cellular mRNA expression. Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1. (A) Venn diagram showing the number of differentially expressed mRNAs at 24 and 48 h after infection. Sets of upregulated mRNAs are represented by upward red arrows and sets of downregulated mRNAs are represented by downward green arrows. Pairs of arrows in the intersection refer to the number of mRNAs upregulated (double red arrows) or down regulated (double green arrows) at both 24 and 48 h after infection. (B) qRT-PCR was conducted on RNA extracted from mock and ZIKV-infected neurons to determine fold-change in IFIT1, IFIT3, IL6, and Caspase1 expression. Changes in the levels of each mRNA were first normalized to the β-actin and then the fold-change in ZIKV-infected cells was calculated in comparison to corresponding mock-infected cells. Data represents the mean ± SEM.
Figure 3. ZIKV infection of the primary mouse neurons causes changes in cellular mRNA expression. Neurons were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1. (A) Venn diagram showing the number of differentially expressed mRNAs at 24 and 48 h after infection. Sets of upregulated mRNAs are represented by upward red arrows and sets of downregulated mRNAs are represented by downward green arrows. Pairs of arrows in the intersection refer to the number of mRNAs upregulated (double red arrows) or down regulated (double green arrows) at both 24 and 48 h after infection. (B) qRT-PCR was conducted on RNA extracted from mock and ZIKV-infected neurons to determine fold-change in IFIT1, IFIT3, IL6, and Caspase1 expression. Changes in the levels of each mRNA were first normalized to the β-actin and then the fold-change in ZIKV-infected cells was calculated in comparison to corresponding mock-infected cells. Data represents the mean ± SEM.
Viruses 11 00162 g003
Figure 4. Enhanced production of cytokines and chemokines in ZIKV-infected neurons. Mouse cortical neuron cultures were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1 and supernatants were collected at 24 and 48 h after infection. Levels of chemokines and cytokines as noted in the figure were measured in cell supernatants using multiplex immunoassay and are expressed as the mean concentration (pg/mL) ± SEM. *p < 0.05, **p < 0.001.
Figure 4. Enhanced production of cytokines and chemokines in ZIKV-infected neurons. Mouse cortical neuron cultures were infected with ZIKV (PRVABC59 strain) or PBS (Mock) at MOI-1 and supernatants were collected at 24 and 48 h after infection. Levels of chemokines and cytokines as noted in the figure were measured in cell supernatants using multiplex immunoassay and are expressed as the mean concentration (pg/mL) ± SEM. *p < 0.05, **p < 0.001.
Viruses 11 00162 g004
Figure 5. Core functional pathway analysis of ZIKV-modulated mRNAs using IPA. Top canonical signaling pathways regulated by significantly modulated mRNAs. Threshold bar indicates cut-off point of significance p < 0.05, using Fisher’s exact test. Range of activation z-score is also depicted in the figure. The color of the bars indicates predicted pathway activation based on z-score (orange = activation; blue = inhibition; gray = no prediction can be made; white = z-score close to 0). Orange line represents the ratio = number of genes in dataset/total number of genes that compose that pathway.
Figure 5. Core functional pathway analysis of ZIKV-modulated mRNAs using IPA. Top canonical signaling pathways regulated by significantly modulated mRNAs. Threshold bar indicates cut-off point of significance p < 0.05, using Fisher’s exact test. Range of activation z-score is also depicted in the figure. The color of the bars indicates predicted pathway activation based on z-score (orange = activation; blue = inhibition; gray = no prediction can be made; white = z-score close to 0). Orange line represents the ratio = number of genes in dataset/total number of genes that compose that pathway.
Viruses 11 00162 g005
Figure 6. Pathway analysis for PRR and IFN signaling. Genes associated with (A) PRR and (B) IFN signaling activated by ZIKV infection are shown. Differentially expressed mRNAs are highlighted in color. Color intensity indicates the degree of upregulation (red) relative to the mock-infected neurons. Solid lines represent direct interactions and dashed lines indirect interactions. Shading intensity indicates the degree that each mRNA was upregulated.
Figure 6. Pathway analysis for PRR and IFN signaling. Genes associated with (A) PRR and (B) IFN signaling activated by ZIKV infection are shown. Differentially expressed mRNAs are highlighted in color. Color intensity indicates the degree of upregulation (red) relative to the mock-infected neurons. Solid lines represent direct interactions and dashed lines indirect interactions. Shading intensity indicates the degree that each mRNA was upregulated.
Viruses 11 00162 g006
Figure 7. Networks of the interactions of the miRNA target genes. IPA tool was used to generate the miRNA-mRNA interaction network of (A) miR-124-3p, (B) miR-654-3p, (C) miR-331-5p, (D) miR-509-5p, and (E) miR-335-3p and mRNAs significantly modulated in neurons after ZIKV infection. Red (increased expression) and green (decreased expression).
Figure 7. Networks of the interactions of the miRNA target genes. IPA tool was used to generate the miRNA-mRNA interaction network of (A) miR-124-3p, (B) miR-654-3p, (C) miR-331-5p, (D) miR-509-5p, and (E) miR-335-3p and mRNAs significantly modulated in neurons after ZIKV infection. Red (increased expression) and green (decreased expression).
Viruses 11 00162 g007
Table 1. Primer sequences used for qRT-PCR.
Table 1. Primer sequences used for qRT-PCR.
Gene (Accession No.)Primer Sequence (5’‒3’)
IFIT1 (NM_008331)
ForwardGTTGTTGTTGTTGTTCGT
ReverseCAGCAGGAATCAGTTGTG
IL6 (NM_031168)
ForwardATCCAGTTGCCTTCTTGGGACTGA
ReverseTAAGCCTCCGACTTGTGAAGTGGT
IFIT3 (NM_010501)
ForwardGTCCTCTCTACTCTTTGG
ReverseCATCCTCTGTCTTCTCTC
Caspase1 (NM_009807)
ForwardGGAAGCAATTTATCAACTCAGTG
ReverseGCCTTGTCCATAGCAGTAATG
Dicer-1 (NM_148948)
ForwardTGTCATCTTGCGATTCTA
ReverseTCTCTTCCAATTCCTCTG
DROSHA (NM_001130149)
ForwardCTTCAACAGTTACCAGAAC
ReverseCCTTTGGGAGTGAGTATG
AGO1 (NM_001317174)
ForwardCCTGTGTATGATGGAAAGA
ReverseCACTTGATGGAGACCTTAA
AGO2 (NM_153178)
ForwardGGAGAACAATCAAACTACAG
ReverseCAGATTCTTCCTTCCATCA
DGCR8 (NM_033324)
ForwardCAGATAAGAAGGATGAGGAA
ReverseGCTCCAAATTGTCAGTAAA
miR-155 (MIMAT0000165)
ForwardUUAAUGCUAAUUGUGAUAGGGGUA
miR-203 (MIMAT0000236)
ForwardGUGAAAUGUUUAGGACCACUAG
miR-29a (MIMAT0000535)
ForwardUAGCACCAUCUGAAAUCGGUUA
miR-124-3p (MIMAT0000134)
ForwardUAAGGCACGCGGUGAAUGCC
Table 2. Upregulated miRNAs in Zika virus (ZIKV)-infected neurons at 24 h and 48 h.
Table 2. Upregulated miRNAs in Zika virus (ZIKV)-infected neurons at 24 h and 48 h.
miRNAFold-change (24 h)miRNAFold-change (48 h)
miR-15518.29miR-15510.61
miR-34714.96miR-2034.61
miR-21452.33miR-19022.95
miR-29a2.09miR-29b2.1
miR-29b2.05miR-29a2.08
Table 3. Downregulated miRNAs at 24 h.
Table 3. Downregulated miRNAs at 24 h.
miRNAFold-changemiRNAFold-change
miR-124-3p−4.02miR-144−2.15
miR-M1-1−2.97miR-201−2.15
miR-1892−2.85miR-764-5p−2.15
miR-883a-3p−2.7miR-1895−2.13
miR-879−2.69miR-871−2.12
miR-669g−2.65miR-m108-1−2.12
miR-654-3p−2.62miR-1928−2.11
miR-339-3p−2.52miR-1941-5p−2.11
miR-1960−2.51miR-759−2.11
miR-207−2.51miR-1187−2.1
miR-2861−2.51miR-666-3p−2.1
miR-412−2.48miR-1966−2.09
miR-741−2.47miR-1942−2.08
miR-770-5p−2.41miR-1946a−2.08
miR-673-5p−2.38miR-1967−2.08
miR-1941-3p−2.37miR-m01-1−2.07
miR-493−2.33miR-383−2.06
miR-465a-3p−2.31miR-186−2.05
miR-1956−2.26miR-323-5p−2.05
miR-1194−2.25miR-M1-3−2.05
miR-709−2.25miR-433−2.04
miR-m107-1-5p−2.25miR-767−2.04
miR-877−2.23miR-1188−2.03
miR-331-5p−2.22miR-694−2.03
miR-483−2.22miR-665−2.02
miR-675-5p−2.2miR-2139−2.02
miR-1957−2.19miR-883a-5p−2.02
miR-M23-1-3p−2.19miR-1943−2
miR-1898−2.18miR-346−2
miR-1940−2.17miR-764-3p−2
miR-1894-5p−2.16miR-710−2
Table 4. Downregulated miRNAs at 48 h.
Table 4. Downregulated miRNAs at 48 h.
miRNAFold-changemiRNAFold-change
miR-124-3p−6.68miR-m21-1−2.15
miR-883a-3p−3.93miR-335-3p−2.14
miR-2137−2.92miR-1957−2.12
miR-2133−2.74miR-764-5p−2.11
miR-714−2.56miR-1194−2.08
miR-669g−2.54miR-683−2.08
miR-467g−2.5miR-509-5p−2.07
miR-1188−2.41miR-463−2.06
miR-879−2.33miR-741−2.06
miR-760−2.25miR-761−2.05
miR-298−2.21miR-710−2.03
miR-m01-3−2.21miR-346−2.01
miR-764-3p−-2.2miR-882−2.01
miR-370−2.18miR-1898−2.01
miR-1892−2.16miR-759−2
miR-m01-2−2.16miR-433−2
miR-133b−2.15miR-709−2
miR-1894-5p−2.15miR-1956−2
miR-666-3p−2.15miR-483−2
miR-877−2.15miR-1941-5p−2
Table 5. Top biological functions regulated by significantly modulated miRNAs.
Table 5. Top biological functions regulated by significantly modulated miRNAs.
Biological Process/Pathwayp-ValueNumber of miRNAs
Cancer4.02 × 10−1018
Organismal Injury and Abnormalities4.02 × 10−1023
Reproductive System Disease4.02 × 10−1018
Immunological Disease6.98 × 10−1013
Inflammatory Disease6.98 × 10−1013
Inflammatory Response6.98 × 10−1011
Neurological Disease6.98 × 10−1010
Connective Tissue Disorders1.46 × 10−810
Respiratory Disease1.46 × 10−87
Nervous System Development and Function7.20 × 10−76
Table 6. Upregulated mRNAs in ZIKV-infected neurons at 24 and 48 h.
Table 6. Upregulated mRNAs in ZIKV-infected neurons at 24 and 48 h.
mRNAFold-Change (24 h)Fold-Change (48 h)mRNAFold-Change (24 h)Fold-Change (48 h)
Rsad226.44244.24Ccl23.9114.32
Cxcl1022.31113.73Ifna42.8514.02
Ifit321.63103.85Nlrc53.4213.39
Ifi4418.36103.12Psmb103.1312.78
Irf712.5396.42Stat25.5512.37
Ifit124.5275.06H2-K1212.34
Isg1520.5869.08Cfb2.0412.04
Lcn227.0768.88Ifi354.0511.3
Zbp113.5755.1Ddx583.8410.56
Gbp512.5144.43Psmb92.410.03
Ifit22.5341.91Tap12.779.85
Ccl524.1137.65C33.299.22
Oas210.6437.13Casp12.767.63
Oasl15.1534.55Socs12.136.63
Usp189.732.04H2-Ab12.215.82
Isg205.5730.07Cxcl13.495.5
H2-T233.6128.33Ifitm13.945.39
Ifih16.1125.37Ccl123.685.38
Ddx607.8723.88Pml25.38
Stat16.9221.99Vcam12.235.24
Bst24.9221.62Myd882.384.87
Mx27.5520.85Cxcl92.834.7
Psmb84.3120.39Il62.073.63
Ccl72.7119.98Ptgs22.023.32
Cd2744.0219.34Cxcl22.683.27
Cmpk24.0118.95Sbno23.052.86
Xaf18.0218.06Ccl42.522.84
Herc63.6817.33Lif2.752.47
Ifnb12.9716.75Il13ra12.012.33
Tlr33.8416.4Runx122.15
Irgm26.2915.13Litaf22.01
Table 7. Upregulated mRNAs in ZIKV-infected neurons at 48 h only.
Table 7. Upregulated mRNAs in ZIKV-infected neurons at 48 h only.
mRNAFold-ChangemRNAFold-Change
Cxcl1117.69Slamf72.81
C29.88Tnfrsf142.76
Tnfsf109.53Fcgr42.6
Cxcl139.29Mill22.59
Cd748.49Ptprc2.59
C1ra7.96H2-Dma2.57
H2-D17.64Nfkbia2.52
H2-Aa7.35Tnfaip32.51
Tap26.53Cxcl52.47
Ccl115.49A2m2.46
Fcgr15.32Ifna12.38
H2-M35.25Il3ra2.36
Ccrl24.9Ctss2.34
Cd474.2Ripk22.34
Tapbp3.97H2-DMb12.28
C4b3.96Lck2.27
Ifna23.96Cd802.26
Irf13.88Cxcl162.19
Tlr23.87Cybb2.19
C1s13.56Icosl2.18
Lbp3.44Cfi2.14
Il73.3Irf22.1
Irf53.03Nod12.02
Cd692.99Atm2
Serping12.91Axl2
Flt3l2.87H2-Eb12
Bid2.85Relb2
Table 8. Downregulated mRNAs in ZIKV-infected neurons.
Table 8. Downregulated mRNAs in ZIKV-infected neurons.
mRNAFold-Change (48 h)mRNAFold-Change (48 h)
Cd36−2Cd207−2.19
Elane−2Timd4−2.22
Il17f−2Xcl1−2.31
Ticam2−2Card9−2.65
Pax5−2.02Sh2d1b1−2.76
Il1rapl2−2.15Mpped1−2.87

Share and Cite

MDPI and ACS Style

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. https://doi.org/10.3390/v11020162

AMA Style

Azouz F, Arora K, Krause K, Nerurkar VR, Kumar M. Integrated MicroRNA and mRNA Profiling in Zika Virus-Infected Neurons. Viruses. 2019; 11(2):162. https://doi.org/10.3390/v11020162

Chicago/Turabian Style

Azouz, Francine, Komal Arora, Keeton Krause, Vivek R. Nerurkar, and Mukesh Kumar. 2019. "Integrated MicroRNA and mRNA Profiling in Zika Virus-Infected Neurons" Viruses 11, no. 2: 162. https://doi.org/10.3390/v11020162

APA Style

Azouz, F., Arora, K., Krause, K., Nerurkar, V. R., & Kumar, M. (2019). Integrated MicroRNA and mRNA Profiling in Zika Virus-Infected Neurons. Viruses, 11(2), 162. https://doi.org/10.3390/v11020162

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