To gain a comprehensive view of the changes in host gene expression underlying Zika virus (ZIKV) pathogenesis, we performed whole-genome RNA sequencing (RNA-seq) of ZIKV-infected Drosophila adult flies. RNA-seq analysis revealed that ZIKV infection alters several and diverse biological processes, including stress, locomotion, lipid metabolism, imaginal disc morphogenesis and regulation of JAK/STAT signaling. To explore the interaction between ZIKV infection and JAK/STAT signaling regulation, we generated genetic constructs overexpressing ZIKV-specific non-structural proteins NS2A, NS2B, NS4A and NS4B. We found that ectopic expression of non-structural proteins in the developing Drosophila eye significantly restricts growth of the larval and adult eye and correlates with considerable repression of the in vivo JAK/STAT reporter, 10XStat92E-GFP. At the cellular level, eye growth defects are associated with reduced rate of proliferation without affecting the overall rate of apoptosis. In addition, ZIKV NS4A genetically interacts with the JAK/STAT signaling components; co-expression of NS4A along with the dominant-negative form of domeless or StatRNAi results in aggravated reduction in eye size, while co-expression of NS4A in HopTuml (also known as hopTum) mutant background partially rescues the hop-induced eye overgrowth phenotype. The function of ZIKV NS4A in regulating growth is maintained in the wing, where ZIKV NS4A overexpression in the pouch domain results in reduced growth linked with diminished expression of Notch targets, Wingless (Wg) and Cut, and the Notch reporter, NRE-GFP. Thus, our study provides evidence that ZIKV infection in Drosophila results in restricted growth of the developing eye and wing, wherein eye phenotype is induced through regulation of JAK/STAT signaling, whereas restricted wing growth is induced through regulation of Notch signaling. The interaction of ZIKV non-structural proteins with the conserved host signaling pathways further advance our understanding of ZIKV-induced pathogenesis.

This article has an associated First Person interview with the first author of the paper.

Zika virus (ZIKV) is an emerging pathogen of substantial public health concern, and belongs to the flavivirus family that also includes dengue, West Nile and yellow fever viruses. Although most infections are asymptomatic or associated with flu-like symptoms (Simpson, 1964), during the recent epidemic, ZIKV has been associated with multi-organ failure resulting in congenital abnormalities in fetuses of pregnant women and neurological complications (Guillain-Barre syndrome) characterized by progressive muscle weakness (Dos Santos et al., 2016; do Rosario et al., 2016). ZIKV can also infect the eye, resulting in conjunctivitis in up to 15% of the patients (Sun et al., 2016; Miner et al., 2016; Furtado et al., 2016). Although ZIKV is primarily transmitted by mosquitos, perinatal and congenital infections, infection through blood transfusion as well as sexual transmission have also been reported (Musso et al., 2014; Mead et al., 2018; Foy et al., 2011). In particular, it is the correlation of ZIKV-associated outbreak with microcephaly that makes ZIKA infection even more serious (Li et al., 2016a). Until now, no drug or vaccine is available to prevent or treat ZIKV infection (Fauci and Morens, 2016).

ZIKV is an obligate intracellular pathogen; therefore, understanding the interaction dynamics between ZIKV and the host and the resulting host pathology is valuable for developing anti-ZIKV therapeutic strategies (Shirasu-Hiza and Schneider, 2007). In particular, the symptoms or physiological consequences of infection prior to disease development reflect the physical state of the organism and contribute to the successful completion of the virus life cycle. These host-virus interactions and the resulting host response are equally important in defining the pathological outcome of an infection. Although efforts have been made to understand the intricacies of ZIKV replication, there have been very few pieces of evidence showing the ZIKV effects on host physiology (Wang et al., 2018; Liang et al., 2016).

Investigating the function of viral components is an elegant strategy to further our understanding of the molecular basis of viral diseases. Structurally, ZIKV is similar to the other flaviviruses and possesses a 25- to 30-nm nucleocapsid surrounded by a host membrane-derived lipid bilayer. ZIKV contains a positive-sense single-stranded RNA of ∼10 kb. Similar to other flaviviruses, host protease-mediated viral processing results in three structural (capsid, pre-membrane, envelope) and seven non-structural (NS) proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, NS5) (Chambers et al., 1990; Wang et al., 2017). In other flaviviruses, the NS proteins are required for viral replication and immune evasion, specifically through interference with RIG-I like receptor (RLR) signaling and type I interferon response (Ngono and Shresta, 2018; Chen et al., 2017). NS2 from dengue and Kunjin virus inhibits interferon-mediated response (Liu et al., 2004; Dalrymple et al., 2015). NS4A and NS4B from flaviviruses are able to inhibit Janus kinase/Signal transducer and activators of transcription (JAK/STAT) and RLR signaling through multiple mechanisms (Ding et al., 2013; Munoz-Jordan et al., 2003). However, to date, there is little information on the functional significance of ZIKV NS proteins and their role in ZIKV-induced pathogenesis. In an earlier study, overexpression of ZIKV NS proteins NS4A and NS4B in the fetal neuronal stem cells (fNSCs) reduced neurosphere formation and inhibited differentiation (Liang et al., 2016). ZIKV NS4A- and NS4B-mediated effect on neurogenesis was further linked to increased autophagy mediated by Akt-mTOR signaling (Liang et al., 2016). This result was specific to ZIKV, as NS4A and NS4B from the closely related dengue virus failed to show a similar effect (Liang et al., 2016).

The common fruit fly, Drosophila melanogaster, with a vast number of genetic and genomic tools available and highly conserved developmental signaling pathways, is widely recognized as an excellent model for studying host-pathogen interactions and human disease. The latter is demonstrated by the fact that 70% of Drosophila genes have human homologs and 75% of human disease-associated genes have homologs in the fly (Bernards and Hariharan, 2001; Edgar and Lehner, 1996; Pandey and Nichols, 2011). Evidence indicates that the fly is also a suitable model for dissecting pathologies related to human pathogenic viruses (Xu and Cherry, 2014; Hughes et al., 2012). In addition, Drosophila can be efficiently used to underscore the in vivo function of viral genes, which are further validated in mammalian models (Hughes et al., 2012). Using the Gal4/UAS system, viral transgenes can be expressed in a spatial and temporal manner followed by the analysis of the resulting phenotype. For instance, overexpression of SARS-CoV 3a and SARS-CoV membrane proteins induces apoptosis in the developing eye through the mitochondrial pathway via Cytochrome c and suppressing survival signaling pathways, respectively (Wong et al., 2005; Chan et al., 2007). Similarly, overexpression of human immunodeficiency virus (HIV) nef gene in the wing results in increased apoptosis without affecting the rate of proliferation (Lee et al., 2005). Overexpression of human cytomegalovirus immediate-early genes in Drosophila embryos results in abnormal embryonic development associated with disrupted adherens junctions (Steinberg et al., 2008). Recently, there have also been efforts to dissect host-ZIKV interactions using Drosophila as the model organism (Harsh et al., 2018; Link et al., 2019; Liu et al., 2018). Findings from these studies indicate that Drosophila can be a reliable model to analyze ZIKV tropism. ZIKV was shown to replicate in the fat body, midgut, crop and brain of the infected adult fly, and result in perturbation in lipid metabolism, intestinal homeostasis and autophagy (Liu et al., 2018; Harsh et al., 2018).

Here, we have challenged Drosophila adult flies with the MR766 strain of ZIKV and used RNA sequencing (RNA-seq) to gain a comprehensive understanding of the differentially regulated genes during infection. We report that ZIKV triggers a large number of biological processes, ranging from misregulation of developmental pathways to perturbed muscle development. Among the developmental signaling pathways, we have shown that ZIKV infection induces negative regulation of JAK/STAT signaling. Eye-specific expression of ZIKV transgenes NS2A, NS2B, NS4A and NS4B results in striking reduction in the size of the developing eye. We also show that the reduced eye size upon ZIKV NS4A protein expression, in particular, correlates with the reduced level of the JAK/STAT reporter, 10XStat92E-GFP. At the cellular level, the expression of ZIKV transgenes results in reduced rate of proliferation in the JAK/STAT-regulated anterior compartment of eye imaginal epithelia without affecting the level of apoptosis. We further demonstrate that ZIKV transgene expression has no effect on the differentiation of the photoreceptors of the eye, and that ZIKV NS4A interacts with the JAK/STAT signaling components. Co-expression of NS4A and the dominant-negative form of domeless or StatRNAi results in aggravated reduction in eye size, while co-expression of NS4A with activated Hop kinase partially rescues the eye enlargement. Finally, the ZIKV NS4A-mediated regulation of growth is also maintained in the wing, where NS4A overexpression restricts the size of the wing pouch. This effect is linked with reduced activity of Notch signaling, which is instrumental for several cell developmental processes, including proliferation and wing development. Altogether, our findings provide the first evidence linking ZIKV-induced pathogenesis and eye/wing development, and reveal a functional link between ZIKV and regulation of JAK/STAT signaling.

ZIKV infection induces distinct transcriptomic profiles in Drosophila

We generated complete transcriptomes from Drosophila wild-type female adult flies infected with the MR766 strain of ZIKV (Harsh et al., 2018; Liu et al., 2018) or injected with PBS (negative/sterile control) (Harsh et al., 2018). We examined gene transcript levels at two timepoints, 4 and 8 days post-injection (dpi) (Fig. 1A). The timepoints were chosen based on the ZIKV load, where 4 dpi corresponds to the early stage of infection, while 8 dpi corresponds to the stage when the infection reaches peak titers (Harsh et al., 2018). The numbers of sequence reads mapped to 96.37% of the D. melanogaster genome (Fig. 1B). Differences in gene expression levels across the two treatments, timepoints and technical replicates are illustrated in the form of a heat map (Fig. 1C). There was a substantial difference in the number of differentially expressed genes (DEGs) when the adjusted P-value was kept at 0.1 versus 0.05 (Fig. 1D). In particular, the number of upregulated genes in ZIKV-injected flies compared to PBS-injected flies increased from 158 to 336 at 4 dpi, and from 148 to 315 at 8 dpi (Fig. 1D). To determine the number of genes that were transcriptionally regulated upon infection with ZIKV at 4 and 8 dpi, we performed pairwise multiple comparison analyses with Limma (Liu et al., 2015; Ritchie et al., 2015). We found that the number of differentially regulated genes varied between 4 and 8 dpi (Fig. 1E). The number of upregulated genes at both 4 and 8 dpi was higher than the number of downregulated genes upon ZIKV infection (Fig. 1E). These results suggest that a large set of genes is differentially regulated in Drosophila adult flies during the early and late stages of infection with ZIKV.

Fig. 1.

Transcriptome analysis of ZIKV-infected flies through RNA-seq. (A) Overview of the experimental workflow. 5- to 6-day-old wild-type female flies (w1118 strain) were injected with the MR766 strain of ZIKV. PBS-injected flies served as negative controls. Total RNA was extracted at 4 and 8 days post-injection (dpi) for RNA-seq. (B) Transcriptome summary (number of reads and percentage mapped to D. melanogaster genome) from wild-type flies injected with PBS and ZIKV at 4 and 8 dpi. (C) Heat map showing the differentially expressed genes (DEGs) in PBS- and ZIKV-injected flies (in triplicate) at 4 and 8 dpi. (D) DEGs (downregulated/upregulated) in wild-type flies injected with ZIKV at 4 and 8 dpi when adjusted P-value is set at 0.05 and 0.1, respectively. (E) Venn diagrams showing the number of Drosophila genes that are differentially expressed (upregulated or downregulated) in wild-type flies injected with ZIKV at 4 and 8 dpi. Expression patterns are indicated (Up/Up, gene upregulation at both 4 and 8 dpi; Down/Down, gene downregulation at both timepoints; Up/Down, gene upregulation at 4 dpi and downregulation at 8 dpi; Down/Up, gene downregulation at 4 dpi and upregulation at 8 dpi). Significance was tested using a hypergeometric test: P-values for Up/Up, Down/Down, Up/Down and Down/Up trending genes were 2.065×10–33, 1.361×10–11, 3.050×10–9 and 5.485×10–11, respectively.

Fig. 1.

Transcriptome analysis of ZIKV-infected flies through RNA-seq. (A) Overview of the experimental workflow. 5- to 6-day-old wild-type female flies (w1118 strain) were injected with the MR766 strain of ZIKV. PBS-injected flies served as negative controls. Total RNA was extracted at 4 and 8 days post-injection (dpi) for RNA-seq. (B) Transcriptome summary (number of reads and percentage mapped to D. melanogaster genome) from wild-type flies injected with PBS and ZIKV at 4 and 8 dpi. (C) Heat map showing the differentially expressed genes (DEGs) in PBS- and ZIKV-injected flies (in triplicate) at 4 and 8 dpi. (D) DEGs (downregulated/upregulated) in wild-type flies injected with ZIKV at 4 and 8 dpi when adjusted P-value is set at 0.05 and 0.1, respectively. (E) Venn diagrams showing the number of Drosophila genes that are differentially expressed (upregulated or downregulated) in wild-type flies injected with ZIKV at 4 and 8 dpi. Expression patterns are indicated (Up/Up, gene upregulation at both 4 and 8 dpi; Down/Down, gene downregulation at both timepoints; Up/Down, gene upregulation at 4 dpi and downregulation at 8 dpi; Down/Up, gene downregulation at 4 dpi and upregulation at 8 dpi). Significance was tested using a hypergeometric test: P-values for Up/Up, Down/Down, Up/Down and Down/Up trending genes were 2.065×10–33, 1.361×10–11, 3.050×10–9 and 5.485×10–11, respectively.

ZIKV infection leads to the enrichment of genes associated with diverse biological processes in Drosophila

To identify the molecular functions and biological activities regulated by ZIKV in Drosophila, we performed Database for Annotation, Visualization and Integrated Discovery (DAVID)-based analyses (Huang da et al., 2009a,b) (Fig. 2). At both 4 and 8 dpi, there was considerable enrichment in the number of genes controlling biological processes and molecular functions (Fig. 2). In order to retrieve a significant number of DEGs, we set a log2 fold change of 1.2 as an arbitrary cutoff threshold and adjusted P-value at 0.1 (Fig. 1D). This cutoff led to identification of 627 DEGs at 4 dpi, of which 336 were upregulated and 291 were downregulated (Tables S1 and S2). At 8 dpi, 706 genes were differentially expressed, of which 315 genes were upregulated and 391 genes were downregulated (Tables S3 and S4). Among the DEGs regulating immunity, Diedel and genes of the Turandot family of proteins were significantly enriched at both 4 and 8 dpi (Fig. 2A,E; Tables S1 and S3). Unlike for Drosophila C virus (DCV) infection (Merkling et al., 2015), we did not encounter enrichment of heat shock proteins at any stage of infection. However, gene ontology (GO) terms such as ‘response to heat, UV and hypoxia’ were significantly enriched at both 4 and 8 dpi (Fig. 2A,E). Another set of genes enriched at both 4 and 8 dpi included genes regulating muscle development (Fig. 2A,E). Interestingly, ZIKV infection also resulted in enrichment of genes regulating developmental pathways, and, in particular, JAK/STAT and Wnt signaling pathways (Fig. 2A,E). Genes that were downregulated at both 4 and 8 dpi were enriched for GO terms such as central nervous system development or synapse assembly (Fig. 2C,G). Genes regulating antibacterial humoral immune response and egg chorion assembly were specifically downregulated at 8 dpi (Fig. 2G). In the context of molecular function, the genes enriched at 4 and 8 dpi were mostly associated with signal transducer activity, actin binding, cadherin binding and structural constituent of chorion (Fig. 2B,D,F,H). Together, these data suggest that ZIKV infection triggers some of the conserved microbe-specific responses, including immune and stress response, and also activates some specific processes related to the regulation of developmental pathways, such as JAK/STAT signaling and vitelline assembly.

Fig. 2.

Infection of adult flies with ZIKV induces diverse biological processes and molecular functions. (A,E) Representative enrichment of upregulated (log>1.2 fold) biological processes using Database for Annotation, Visualization and Integrated Discovery (DAVID) classification database at 4 and 8 dpi. (C,G) Representative enrichment of downregulated (log<−1.2 fold) biological processes using DAVID at 4 and 8 dpi. (B,D,F,H) GO-based molecular functions regulated by the DEGs at 4 dpi (B,D) and 8 dpi (F,H).

Fig. 2.

Infection of adult flies with ZIKV induces diverse biological processes and molecular functions. (A,E) Representative enrichment of upregulated (log>1.2 fold) biological processes using Database for Annotation, Visualization and Integrated Discovery (DAVID) classification database at 4 and 8 dpi. (C,G) Representative enrichment of downregulated (log<−1.2 fold) biological processes using DAVID at 4 and 8 dpi. (B,D,F,H) GO-based molecular functions regulated by the DEGs at 4 dpi (B,D) and 8 dpi (F,H).

qRT-PCR-based transcriptome analysis reveals upregulation of negative regulators of JAK/STAT signaling upon ZIKV infection

To validate our RNA-seq results (Fig. 2), we performed quantitative real-time PCR (qRT-PCR) on Drosophila flies infected with ZIKV at 4 and 8 dpi (Fig. 3). Interestingly, we found that one of the prominent ZIKV infection-induced processes relates to the negative regulation of JAK/STAT signaling. First identified as a key regulator of interferon and cytokine signaling in mammals, JAK/STAT signaling regulates pleiotropic effects including growth and differentiation (Schindler et al., 1992; Shuai et al., 1994, 1992; Velazquez et al., 1992; Watling et al., 1993). JAK/STAT signaling is also known to regulate cellular proliferation, immune responses and maintenance, and proliferation of stem cells in the gonads (reviewed in Arbouzova and Zeidler, 2006). Given the range of biological roles played by JAK/STAT signaling, it is not surprising that this pathway is controlled through multiple regulatory mechanisms. The wide range of phenotypes associated with the expression of negative regulators of JAK/STAT signaling prompted us to validate the RNA-seq results.

Fig. 3.

Validation of the enriched negative regulators of JAK/STAT pathway by qRT-PCR. (A,B) Log2 fold change (RNA-seq) and mRNA levels (qRT-PCR) of Socs36E, E(bx) and Et in wild-type adult flies infected with ZIKV at 4 and 8 dpi. (C,D) Log2 fold change (RNA-seq) and mRNA levels (qRT-PCR) of Tep1, TotM, Diedel, CecA1 and Dipt in wild-type adult flies infected with ZIKV at 4 and 8 dpi. All data were normalized to the housekeeping gene RpL32 and are shown relative to wild-type flies injected with PBS (sterile control). Three independent experiments were carried out with ten flies per sample in triplicate (‘*P<0.05, **P<0.01, ***P<0.001, ****P<0.0001). Bars represent the mean±s.d. Statistical analysis was performed using unpaired two-tailed Student's t-test.

Fig. 3.

Validation of the enriched negative regulators of JAK/STAT pathway by qRT-PCR. (A,B) Log2 fold change (RNA-seq) and mRNA levels (qRT-PCR) of Socs36E, E(bx) and Et in wild-type adult flies infected with ZIKV at 4 and 8 dpi. (C,D) Log2 fold change (RNA-seq) and mRNA levels (qRT-PCR) of Tep1, TotM, Diedel, CecA1 and Dipt in wild-type adult flies infected with ZIKV at 4 and 8 dpi. All data were normalized to the housekeeping gene RpL32 and are shown relative to wild-type flies injected with PBS (sterile control). Three independent experiments were carried out with ten flies per sample in triplicate (‘*P<0.05, **P<0.01, ***P<0.001, ****P<0.0001). Bars represent the mean±s.d. Statistical analysis was performed using unpaired two-tailed Student's t-test.

We noticed that the established negative regulators of JAK/STAT signaling did not show significant upregulation (in terms of the adjusted P-value). Also, comparing the individual fold change in all three biological replicates revealed that, out of the three replicates, one replicate deviated from the other two replicates, resulting in the higher P-value and hence non-significant upregulation (Table S5). However, consistent with several pieces of evidence showing the interaction of flaviviruses with the regulators of JAK/STAT signaling (Laurent-Rolle et al., 2010; Kumthip et al., 2012; Best, 2017; Grant et al., 2016), we nevertheless decided to proceed with the validation of these negative regulators. qRT-PCR analysis indeed showed that the negative regulators of JAK/STAT signaling, E(bx) and Socs36E, were significantly upregulated in ZIKV-infected Drosophila flies at both 4 dpi (Fig. 3A) and 8 dpi (Fig. 3B). RNA-seq data were also validated through qRT-PCR of immunity-related genes, and, in particular, showed that Tep1, TotM and Diedel were strongly upregulated, whereas CecA1 and Dipt were downregulated (Fig. 3C,D). Thus, qRT-PCR analysis showed a good correlation with RNA-seq data and suggests that Drosophila can be a suitable model for studying host-ZIKV interaction, where ZIKV can trigger unique processes including negative regulation of the well-characterized JAK/STAT signaling pathway.

ZIKV transgene expression triggers retarded eye growth linked with reduced levels of JAK/STAT signaling

The demonstration of upregulation of Et, E(bx) and Socs36E through qRT-PCR analysis prompted us to probe the functional significance of the altered expression of these genes, particularly in the context of ZIKV transgene-induced eye development. In order to examine the effect of ZIKV on Drosophila eye development, we took advantage of genetic overexpression constructs. Gal4/UAS-based studies in Drosophila have led to identification of the functions of several viral genes, facilitating a comprehensive understanding of the outcome of viral infection (Hughes et al., 2012).

Along with the viral genes that encode structural proteins, the non-structural proteins are also critical because they interact with the host cells to promote viral pathogenesis. Non-structural proteins of flaviviruses play crucial roles in their replication (Bollati et al., 2010). ZIKV, like other flaviviruses, possesses seven non-structural proteins: NS1, NS2A, NS2B, NS3, NS4A, NS4B and NS5 (Shi and Gao, 2017; Hasan et al., 2018). Here, we would like to emphasize that although the perturbation of JAK/STAT signaling was revealed from the infection dynamics of the whole fly, we switched to the Gal4/UAS-based genetic system to validate this finding and further explore the functions of ZIKV NS proteins. We generated genetic constructs overexpressing four of these non-structural proteins, namely UAS-NS2A, UAS-NS2B, UAS-NS4A and UAS-NS4B. To score for the developmental defect upon ZIKV NS protein overexpression, we crossed these constructs with an eye-specific, eyeless-Gal4 (referred to as E1-Gal4), where Eyeless is one of the most critical transcription factors regulating eye development in Drosophila (Pallavi et al., 2012). Overexpression of genes encoding ZIKV NS2A, NS2B, NS4A and NS4B under E1-Gal4 did not result in any survival defects (Fig. S1A). Scoring for the time of pupariation further showed that these flies did not exhibit any developmental delay (Fig. S1B). However, compared to E1-Gal4, overexpression of genes encoding ZIKV non-structural proteins NS2A, NS2B, NS4A and NS4B driven under E1-Gal4 resulted in significantly reduced size of the eye imaginal disc, a phenotype also resulting from loss of function of hop (Perrimon and Mahowald, 1986; Mukherjee et al., 2005) (Fig. 4A,B).

Fig. 4.

ZIKV non-structural proteins induce restricted eye growth and are linked with downregulated JAK/STAT signaling. (A) Representative eye imaginal discs upon overexpression of ZIKV non-structural protein-coding genes driven under eye-specific eyeless-Gal4 (referred to as E1-Gal4) (E1>NS4A, E1>NS4B, E1>NS2B and E1>NS2A). Lower row shows the enlarged view of the eye imaginal discs. The cytoarchitecture was marked with Actin (gray). (B) Quantification of the size of eye imaginal discs upon E1-Gal4 driven overexpression of ZIKV non-structural protein coding genes, compared to E1-Gal4 alone. Bars show mean±s.d. (*P<0.05, **P<0.01; ns, not significant). (C) Representative images of the adult eye in the indicated genotypes. (D) Quantification of the size of adult eye upon overexpression of ZIKV non-structural protein-coding genes, compared to E1-Gal4. Bars show mean±s.d. (***P=0.0001, **P=0.0016, *P<0.05). (E) qRT-PCR analysis depicting the mRNA level of NS4A in the eye imaginal epithelia overexpressing ZIKV NS4A driven under E1-Gal4 (E1>NS4A), compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001). (F) qRT-PCR analysis depicting the mRNA level of chinmo, Mo25 and domeless in the eye imaginal disc where ZIKV NS4A was overexpressed. All data were normalized to the housekeeping gene RpL32 and are shown relative to E1-Gal4. Three independent experiments were carried out in triplicate and bars represent mean±s.d. (**P<0.05). (G) Expression of 10XStat92E-GFP (abbreviated to 10XStat-GFP, shown by red in merge and gray in separate channel) in eye imaginal epithelia carrying eyeless-specific overexpression of ZIKV NS4A compared to E1-Gal4 alone. The cytoarchitecture is marked with Actin (green). The 10XStat-GFP levels are marked by yellow arrows in both E1-Gal4- and E1>NS4A-carrying eye imaginal epithelia. (H) Quantification of 10XStat-GFP fluorescence intensity in the eye imaginal disc proper from E1-Gal4- and E1>NS4A-carrying larvae. Bars show mean±s.d. (**P=0.0015). Statistical analysis was performed using unpaired two-tailed Student's t-test. AU, arbitrary units. Scale bars: 100 μm.

Fig. 4.

ZIKV non-structural proteins induce restricted eye growth and are linked with downregulated JAK/STAT signaling. (A) Representative eye imaginal discs upon overexpression of ZIKV non-structural protein-coding genes driven under eye-specific eyeless-Gal4 (referred to as E1-Gal4) (E1>NS4A, E1>NS4B, E1>NS2B and E1>NS2A). Lower row shows the enlarged view of the eye imaginal discs. The cytoarchitecture was marked with Actin (gray). (B) Quantification of the size of eye imaginal discs upon E1-Gal4 driven overexpression of ZIKV non-structural protein coding genes, compared to E1-Gal4 alone. Bars show mean±s.d. (*P<0.05, **P<0.01; ns, not significant). (C) Representative images of the adult eye in the indicated genotypes. (D) Quantification of the size of adult eye upon overexpression of ZIKV non-structural protein-coding genes, compared to E1-Gal4. Bars show mean±s.d. (***P=0.0001, **P=0.0016, *P<0.05). (E) qRT-PCR analysis depicting the mRNA level of NS4A in the eye imaginal epithelia overexpressing ZIKV NS4A driven under E1-Gal4 (E1>NS4A), compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001). (F) qRT-PCR analysis depicting the mRNA level of chinmo, Mo25 and domeless in the eye imaginal disc where ZIKV NS4A was overexpressed. All data were normalized to the housekeeping gene RpL32 and are shown relative to E1-Gal4. Three independent experiments were carried out in triplicate and bars represent mean±s.d. (**P<0.05). (G) Expression of 10XStat92E-GFP (abbreviated to 10XStat-GFP, shown by red in merge and gray in separate channel) in eye imaginal epithelia carrying eyeless-specific overexpression of ZIKV NS4A compared to E1-Gal4 alone. The cytoarchitecture is marked with Actin (green). The 10XStat-GFP levels are marked by yellow arrows in both E1-Gal4- and E1>NS4A-carrying eye imaginal epithelia. (H) Quantification of 10XStat-GFP fluorescence intensity in the eye imaginal disc proper from E1-Gal4- and E1>NS4A-carrying larvae. Bars show mean±s.d. (**P=0.0015). Statistical analysis was performed using unpaired two-tailed Student's t-test. AU, arbitrary units. Scale bars: 100 μm.

We further examined whether ZIKV-induced infection affects retinal differentiation and found that eye imaginal epithelia from E1-Gal4-driven overexpression of NS2A, NS2B, NS4A and NS4B did not display any defect in the expression of pan-neural marker, Elav, which marks the retinal neuron-specific fate (Fig. S2). In corroboration with the reduced size of eye imaginal epithelia, we noticed that the resulting adult flies harboring E1-Gal4-mediated ZIKV non-structural protein overexpression also displayed a significant reduction in eye size (Fig. 4C,D). We noticed that ZIKV-induced reduction in adult eye size was very similar to the ‘small eye’ phenotype triggered upon loss of function of unpaired (upd) (Tsai and Sun, 2004; Bach et al., 2003). To ascertain that the eye phenotype is indeed a consequence of ZIKV transgene expression, we examined the transcriptional status of NS4A in the eye imaginal epithelia carrying E1-Gal4-mediated NS4A overexpression. Assigning the arbitrary value of 1 to the level of NS4A in the control eye imaginal epithelia (E1-Gal4), there was 3000-fold enrichment of NS4A in the eye imaginal epithelia when NS4A was overexpressed (Fig. 4E).

Given the similarity between ZIKV- and JAK/STAT-induced eye phenotypes, we next aimed to confirm whether the ZIKV-induced eye phenotype is indeed linked to impaired JAK/STAT signaling. qRT-PCR analysis revealed that the eye imaginal epithelia carrying E1-Gal4-driven NS4A overexpression showed significantly reduced expression of the targets of JAK/STAT signaling – including chinmo, Mo25 and domeless – compared to E1-Gal4 (Flaherty et al., 2009) (Fig. 4F). For a more robust validation, we switched to the in vivo reporter, 10XStat92E-GFP, which accurately reflects the activation of the JAK/STAT pathway in several tissues (Bach et al., 2007). 10XStat92E-GFP is expressed throughout the development of Drosophila, including the embryo and larval stages, in a pattern overlapping with the expression pattern of Stat92E protein (Bach et al., 2007). Moreover, it is activated by ectopic JAK/STAT signaling and there is loss of reporter expression in Stat92E clones (Bach et al., 2007). The control eye imaginal disc showed 10XStat92E-GFP expression restricted to the antenna and the posterior compartment of the eye imaginal epithelia (Fig. 4G). Interestingly, overexpression of eye-specific ZIKV NS4A protein resulted in a substantial reduction of 10XStat92E-GFP expression (Fig. 4G). We also noticed that this reduction was autonomous, as NS4A overexpression had no effect on the expression pattern of 10XStat92E-GFP in the antenna of the eye imaginal disc. We further quantified the fluorescent intensity of the STAT reporter in the E1- and E1>NS4A-carrying eye imaginal epithelia and found a significant loss of STAT reporter activity where ZIKV NS4A was overexpressed (Fig. 4H). These findings suggest that ZIKV-regulated eye growth in Drosophila is linked with the substantial downregulation of the JAK/STAT signaling pathway.

ZIKV-induced eye phenotype is associated with a reduced rate of proliferation without affecting the rate of apoptosis

We next examined the eye imaginal epithelia overexpressing ZIKV non-structural proteins for the overall rate of proliferation. We used the pan mitosis-marker, anti-phosphohistone (PH3), for detecting the mitotically active cells. We noticed numerous PH3-positive cells in both the anterior (first mitotic wave) and posterior compartment (second mitotic wave) in the eye imaginal epithelia in the case of E1-Gal4 (Fig. 5A). Upon overexpression of eye-specific ZIKV non-structural proteins, the number of PH3-positive cells was reduced significantly (Fig. 5A). Quantification analysis showed that cells from both the first (Fig. 5B) and second (Fig. 5C) mitotic wave responded in a similar manner, and the total number of PH3-positive cells in the eye imaginal epithelia overexpressing ZIKV non-structural proteins was reduced to 50% of the PH3-positive cells in E1-Gal4 alone (Fig. 5B,C).

Fig. 5.

ZIKV non-structural protein overexpression results in reduced rate of proliferation in eye imaginal epithelia. (A,D) Representative eye imaginal discs overexpressing ZIKV non-structural protein-coding genes driven under eye-specific E1-Gal4 (E1>NS4A, E1>NS4B, E1>NS2B and E1>NS2A). (A) The rate of proliferation was marked with anti-phosphohistone (PH3) staining (green). The yellow lines mark the morphogenetic furrow. The compartments are marked: a, anterior; p, posterior. In all images, nuclei were stained with DAPI (blue). (B) Quantification of the number of PH3-marked cells in the first mitotic wave or anterior region of eye discs upon ZIKV non-structural protein overexpression compared to E1-Gal4 alone. Bars show mean±s.d. (*P<0.05, **P<0.01). (C) Quantification of the number of PH3-marked cells in the second mitotic wave or posterior region of eye discs upon ZIKV non-structural protein overexpression compared to E1-Gal4 alone. Bars show mean±s.d. (*P<0.05, **P<0.01). (D) Cell death was indicated by anti-Dcp-1 staining (green); Actin was used to mark the cytoarchitecture (red). In all images, nuclei were stained with DAPI (blue). Statistical analysis for the graphs were performed using unpaired two-tailed Student's t-test. Scale bars: 100 μm.

Fig. 5.

ZIKV non-structural protein overexpression results in reduced rate of proliferation in eye imaginal epithelia. (A,D) Representative eye imaginal discs overexpressing ZIKV non-structural protein-coding genes driven under eye-specific E1-Gal4 (E1>NS4A, E1>NS4B, E1>NS2B and E1>NS2A). (A) The rate of proliferation was marked with anti-phosphohistone (PH3) staining (green). The yellow lines mark the morphogenetic furrow. The compartments are marked: a, anterior; p, posterior. In all images, nuclei were stained with DAPI (blue). (B) Quantification of the number of PH3-marked cells in the first mitotic wave or anterior region of eye discs upon ZIKV non-structural protein overexpression compared to E1-Gal4 alone. Bars show mean±s.d. (*P<0.05, **P<0.01). (C) Quantification of the number of PH3-marked cells in the second mitotic wave or posterior region of eye discs upon ZIKV non-structural protein overexpression compared to E1-Gal4 alone. Bars show mean±s.d. (*P<0.05, **P<0.01). (D) Cell death was indicated by anti-Dcp-1 staining (green); Actin was used to mark the cytoarchitecture (red). In all images, nuclei were stained with DAPI (blue). Statistical analysis for the graphs were performed using unpaired two-tailed Student's t-test. Scale bars: 100 μm.

We next examined the rate of apoptosis as one of the cellular factors contributing to the reduced size of the eye imaginal epithelia overexpressing ZIKV non-structural proteins. Unlike proliferation, we did not notice any reduction in the level of apoptosis in the eye imaginal epithelia overexpressing ZIKV non-structural proteins compared to E1-Gal4 alone (Fig. 5D). Interestingly, however, ZIKV non-structural protein overexpression did trigger enhanced cell death in the optic stalk as shown by increased accumulation of anti-Dcp-1 (Fig. 5D). Together, these findings suggest that, at the cellular level, ZIKV infection triggers reduction in overall proliferation rate in the eye imaginal epithelia without significantly affecting apoptosis.

Knockdown of Stat along with NS4A overexpression leads to synergistic reduction in eye size, while NS4A overexpression rescues Hop-mediated eye overgrowth

After validating the status of JAK/STAT signaling in the developing eye tissue carrying overexpression of ZIKV NS4A protein, we next investigated the genetic interaction of ZIKV NS4A with the regulators of JAK/STAT signaling. We took advantage of estimating the adult as well as larval eye size as the assay for scoring genetic interaction. We genetically recombined the overexpression construct of ZIKV non-structural gene NS4A (UAS-NS4A) and dominant-negative form of Upd receptor, Domeless (UAS-domeDN) (Brown et al., 2001). The overexpression of the dominant-negative form of domeless driven under E1-Gal4 (E1>domeDN) resulted in a drastic reduction in the size of the eye imaginal epithelia (Tsai and Sun, 2004) (Fig. 6A). Introduction of the dominant-negative form of domeless (UAS-domeDN) in the background of E1-Gal4-driven NS4A (E1>NS4A) overexpression resulted in aggravated reduction in eye size (Fig. 6A,B).

Fig. 6.

ZIKV NS4A shows genetic interaction with different components of JAK/STAT signaling pathway. (A) Representative eye imaginal discs overexpressing ZIKV NS4A, a dominant negative form of domeless, and co-expression of the dominant-negative form of domeless and NS4A driven under the eye-specific driver, E1-Gal4 (E1>NS4A, E1>domeDN and E1>NS4A, domeDN, respectively), compared to E1-Gal4 alone. The lower row shows the enlarged view of the eye imaginal discs. Cytoarchitecture was marked with Actin (gray). (B) Quantification of the size of eye imaginal discs in the indicated genotypes. Bars show mean±s.d. (***P<0.0001, **P<0.05). (C) Representative images of the adult eye in E1>NS4A, E1>domeDN and E1>NS4A, domeDN, respectively. (D) Percentage of individuals displaying normal, moderate, small, very small or absence of adult retinal area in the indicated genotypes. (E) Quantification of the size of adult eye in E1>NS4A, E1>domeDN and E1>NS4A, domeDN, respectively, compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001, **P=0.0053). (F) Representative images of the adult eye overexpressing ZIKV NS4A, StatRNAi and co-expression of StatRNAi and NS4A driven under eye-specific E1-Gal4 (E1>NS4A, E1>StatRNAi and E1>NS4A, StatRNAi, respectively). (G) Penetrance of different eye phenotypes in the indicated genotypes. (H) Quantification of the size of adult eye in E1>NS4A, E1>StatRNAi and E1>NS4A, StatRNAi, respectively, compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001, ***P<0.001 and **P=0.0071). (I) Representative images of the adult eye upon E1-Gal4-driven overexpression of HopTuml and co-expression of NS4A and HopTuml (E1>HopTuml and E1>NS4A, HopTuml, respectively). (J) Quantification of the size of adult eye in E1>HopTuml and E1>NS4A, HopTuml, compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001 and ***P=0.0006). Statistical analysis was performed using unpaired two-tailed Student's t-test. Scale bars: 100 μm.

Fig. 6.

ZIKV NS4A shows genetic interaction with different components of JAK/STAT signaling pathway. (A) Representative eye imaginal discs overexpressing ZIKV NS4A, a dominant negative form of domeless, and co-expression of the dominant-negative form of domeless and NS4A driven under the eye-specific driver, E1-Gal4 (E1>NS4A, E1>domeDN and E1>NS4A, domeDN, respectively), compared to E1-Gal4 alone. The lower row shows the enlarged view of the eye imaginal discs. Cytoarchitecture was marked with Actin (gray). (B) Quantification of the size of eye imaginal discs in the indicated genotypes. Bars show mean±s.d. (***P<0.0001, **P<0.05). (C) Representative images of the adult eye in E1>NS4A, E1>domeDN and E1>NS4A, domeDN, respectively. (D) Percentage of individuals displaying normal, moderate, small, very small or absence of adult retinal area in the indicated genotypes. (E) Quantification of the size of adult eye in E1>NS4A, E1>domeDN and E1>NS4A, domeDN, respectively, compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001, **P=0.0053). (F) Representative images of the adult eye overexpressing ZIKV NS4A, StatRNAi and co-expression of StatRNAi and NS4A driven under eye-specific E1-Gal4 (E1>NS4A, E1>StatRNAi and E1>NS4A, StatRNAi, respectively). (G) Penetrance of different eye phenotypes in the indicated genotypes. (H) Quantification of the size of adult eye in E1>NS4A, E1>StatRNAi and E1>NS4A, StatRNAi, respectively, compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001, ***P<0.001 and **P=0.0071). (I) Representative images of the adult eye upon E1-Gal4-driven overexpression of HopTuml and co-expression of NS4A and HopTuml (E1>HopTuml and E1>NS4A, HopTuml, respectively). (J) Quantification of the size of adult eye in E1>HopTuml and E1>NS4A, HopTuml, compared to E1-Gal4 alone. Bars show mean±s.d. (****P<0.0001 and ***P=0.0006). Statistical analysis was performed using unpaired two-tailed Student's t-test. Scale bars: 100 μm.

The adult fly carrying the overexpression of dominant-negative form of domeless under E1-Gal4 displayed variable phenotype, ranging from complete absence of the eye to small and very small eye size (Fig. 6D). Compared to NS4A overexpression alone, where most of the flies showed moderately small eye size (Fig. 6C,D), introduction of the dominant-negative form of domeless in the background of E1-Gal4-mediated NS4A overexpression resulted in a range of eye phenotypes (Fig. 6C,D). The majority of the flies showed either small or very small eye size, while the flies bearing moderately sized eye reduced to 12% of the total number of flies tested (Fig. 6D). Quantification analysis further validated that the co-expression of UAS-domeDN and ZIKV NS4A overexpression aggravates the reduction in eye size compared to expression of E1>UAS-NS4A or E1>UAS-domeDN alone (Fig. 6E). We further validated this synergistic reduction by co-expressing StatRNAi and ZIKV NS4A overexpression using eye-specific E1-Gal4. Similar to domeDN, E1>StatRNAi resulted in reduced size of the adult eye and also demonstrated penetrance of a range of adult eye phenotypes ranging from very small to total absence of eye (Fig. 6F,G). Consistent with the role of Stat in eye development, clonal loss of Stat92E results in a small or ablated eye (Ekas et al., 2006). Quantification of the size of adult eye showed that eyeless-specific overexpression of ZIKV NS4A and StatRNAi resulted in significantly smaller eyes than E1>NS4A or E1>StatRNAi alone (Fig. 6H).

To further validate the genetic interaction between ZIKV NS4A and JAK/STAT signaling, we then overexpressed JAK/STAT signaling components and asked if co-expression of ZIKV NS4A could rescue the STAT-mediated phenotype. Ectopic misexpression of the JAK/STAT ligand, Upd, using E1-Gal4 (E1>upd) results in enlarged eye phenotype (Bach et al., 2003) (Fig. S3). However, after co-expression of ZIKV NS4A and upd together in the developing eye, we still observed an enlarged eye (Fig. S3). We then switched to a less-severe eye phenotype to examine the effect of ZIKV NS4A overexpression. Overexpression of activated Hop kinase, HopTuml (also known as hopTum), also results in enlargement of the eye (Bach et al., 2003) (Fig. 6I,J). When we co-expressed ZIKV NS4A and HopTuml under E1-Gal4, we noticed that there was a significant, although modest, reduction in the size of the developing eye (Fig. 6I,J). Thus, the aggravation of the StatRNAi- and domeDN-induced eye phenotype and the partial rescue of the HopTuml-mediated overgrowth, confirms the genetic interaction between JAK/STAT signaling components and ZIKV NS4A.

ZIKV transgene expression fails to regulate the growth of posterior compartment of larval wing epithelia but engrailed-specific NS4A overexpression triggers thickening of veins in adult wing

Apart from regulating eye development, JAK/STAT signaling also regulates wing development and patterning. The mutation in upd resulting in ‘outstretched’ phenotype is characterized by the outstretched wing posture in the adult fly. The outstretched wings phenotype was also associated with the upregulation of a negative regulator of JAK/STAT signaling, Socs36E (Callus and Mathey-Prevot, 2002). Driven under a wing-specific Gal4, directed expression of Socs36E further resulted in loss of the wing anterior cross vein, formation of wing vein deltas and humeral outgrowths (Callus and Mathey-Prevot, 2002). In a recent study, JAK/STAT was also shown to regulate the posterior compartment of wing imaginal epithelia and promote cycling and survival of these cells (Recasens-Alvarez et al., 2017). The JAK/STAT-mediated reduced size of the posterior compartment was further attributed to the decreased expression of cell cycle marker, CycA, and increased apoptosis (Recasens-Alvarez et al., 2017).

We next investigated if ZIKV transgene expression triggers wing defects in Drosophila. Based on JAK/STAT signaling-induced wing phenotypes, ZIKV non-structural proteins were overexpressed using a wing-specific engrailed-Gal4 [Engrailed (En) marks the posterior compartment of the wing]. We noticed that there was no difference in the size of the posterior compartment (tagged with RFP), compared to the anterior compartment of the wing imaginal disc (shown by the absence of RFP), when engrailed-Gal4-specific ZIKV non-structural proteins were expressed (Fig. 7A). We also did not encounter any noticeable change in the level of CycA in the two compartments, compared to the anterior compartment of wing imaginal disc, when engrailed-Gal4-specific ZIKV non-structural proteins were expressed (Fig. 7A). We further examined apoptotic activity in the posterior compartment and found that ZIKV non-structural protein overexpression failed to trigger any change in the level of apoptosis (Fig. 7B). These findings in the larval stage were further consistent in the adult stage, when the eclosed flies did not show any defect in the size of the adult wing (Fig. 7C). Interestingly, overexpression of ZIKV NS4A, however, resulted in a very reproducible and highly penetrant (90% of the flies, n=50) phenotype marked by thickening of veins and formation of wing deltas (Fig. 7C), a phenotype characteristic of Socs36E overexpression (Callus and Mathey-Prevot, 2002). These findings suggest that while ZIKV non-structural proteins unambiguously regulate eye growth, the wing phenotype is mainly restricted to the overexpression of ZIKV NS4A.

Fig. 7.

ZIKV NS4A overexpression results in thickening of veins in adult wing. (A) Representative wing imaginal discs upon overexpression of ZIKV non-structural protein-coding genes driven under wing-specific, engrailed-Gal4 (en>NS4A, en>NS2A, en>NS2B and en>NS4B). The engrailed marked compartment was tagged with RFP (enGal4>UAS-RFP). Wing disc-specific JAK/STAT target, CycA, is marked in green. (B) Anti-Dcp-1 was used to show cell death (green) in the indicated genotypes. (C) Representative images of the adult wing in en>NS4A, en>NS2A, en>NS2B and en>NS4B. The thickening and branching of veins are depicted with red arrows. Scale bars: 100 μm.

Fig. 7.

ZIKV NS4A overexpression results in thickening of veins in adult wing. (A) Representative wing imaginal discs upon overexpression of ZIKV non-structural protein-coding genes driven under wing-specific, engrailed-Gal4 (en>NS4A, en>NS2A, en>NS2B and en>NS4B). The engrailed marked compartment was tagged with RFP (enGal4>UAS-RFP). Wing disc-specific JAK/STAT target, CycA, is marked in green. (B) Anti-Dcp-1 was used to show cell death (green) in the indicated genotypes. (C) Representative images of the adult wing in en>NS4A, en>NS2A, en>NS2B and en>NS4B. The thickening and branching of veins are depicted with red arrows. Scale bars: 100 μm.

ZIKV NS4A regulates wing growth linked with downregulated Notch signaling

We noticed that ZIKV NS4A-induced wing vein thickening is also a characteristic phenotype of mutation in Notch signaling. This prompted us to examine other possible wing defects upon overexpression of ZIKV non-structural proteins. We used nubbin-Gal4, where Nubbin is specific to the pouch domain of the wing imaginal disc, to specify the proximo-distal axis in Drosophila wing (Ng et al., 1995; Neumann and Cohen, 1998; Cifuentes and Garcia-Bellido, 1997). Nubbin-driven overexpression of NS4A resulted in a classic notched wing phenotype marked by loss of wing margin (Fig. 8A). Furthermore, this wing phenotype was highly reproducible and penetrant (85% of the flies, n=55). Similar to the engrailed-Gal4 findings, we found that overexpression of NS2A, NS2B and NS4B had no noticeable defect in the adult wing (Table S6). The classical notched wing phenotype is attributed to downregulated Notch signaling. We next examined the underlying mechanism for NS4A-induced notched phenotype. We found that nub-Gal4-driven NS4A expression resulted in drastic reduction in the size of the pouch region of the wing imaginal epithelia (Fig. 8B,C). In addition, the wing imaginal epithelia displayed increased apoptosis as marked by increased expression of Dcp-1 (Fig. 8B).

Fig. 8.

ZIKV NS4A overexpression results in notching of wing linked with downregulated Notch signaling. (A) Representative images of the adult wing upon overexpression of ZIKV NS4A driven under wing-specific, nubbin-Gal4 (nub>NS4A), compared to nub-Gal4 alone. (B-E) Representative wing imaginal discs displaying overexpression of ZIKV NS4A driven under nub-Gal4 (nub>NS4A), compared to nub-Gal4 alone. For clarity, the nubbin region is marked with a yellow outline. (B) Cell death was indicated with anti-Dcp-1 staining (blue in merge panel and shown in gray in separate channel). (C) Quantification of the size of the wing pouch upon overexpression of ZIKV non-structural protein NS4A under nub-Gal4 (nub>NS4A), compared to nub-Gal4 alone. The size of the wing pouch in each genotype was normalized with respect to the size of the whole wing disc. Bars show mean±s.d. (***P=0.0001). (D,E) Wingless and Cut, the targets of Notch signaling, are marked in red. Reduced expression of Wg and Cut upon NS4A overexpression (nub>NS4A) is indicated by yellow arrows. (F) Overexpression of NS4A in the wing pouch (nub>NS4A) significantly downregulated Notch signaling reporter NRE-GFP expression at the dorsal/ventral boundary, compared to nubbin-Gal4 alone. NRE-GFP is shown in red in the merge images and in gray in the separate channel. (G) The amount of Notch intracellular domain (NICD) was significantly reduced when ZIKV NS4A was overexpressed using nubbin-Gal4 (nub>NS4A), compared to nubbin-Gal4 alone. Notch protein level was marked with anti-NICD staining (red in merge channel and gray in the separate channel). In all images, cytoarchitecture was marked with Actin (green). Scale bars: 100 μm.

Fig. 8.

ZIKV NS4A overexpression results in notching of wing linked with downregulated Notch signaling. (A) Representative images of the adult wing upon overexpression of ZIKV NS4A driven under wing-specific, nubbin-Gal4 (nub>NS4A), compared to nub-Gal4 alone. (B-E) Representative wing imaginal discs displaying overexpression of ZIKV NS4A driven under nub-Gal4 (nub>NS4A), compared to nub-Gal4 alone. For clarity, the nubbin region is marked with a yellow outline. (B) Cell death was indicated with anti-Dcp-1 staining (blue in merge panel and shown in gray in separate channel). (C) Quantification of the size of the wing pouch upon overexpression of ZIKV non-structural protein NS4A under nub-Gal4 (nub>NS4A), compared to nub-Gal4 alone. The size of the wing pouch in each genotype was normalized with respect to the size of the whole wing disc. Bars show mean±s.d. (***P=0.0001). (D,E) Wingless and Cut, the targets of Notch signaling, are marked in red. Reduced expression of Wg and Cut upon NS4A overexpression (nub>NS4A) is indicated by yellow arrows. (F) Overexpression of NS4A in the wing pouch (nub>NS4A) significantly downregulated Notch signaling reporter NRE-GFP expression at the dorsal/ventral boundary, compared to nubbin-Gal4 alone. NRE-GFP is shown in red in the merge images and in gray in the separate channel. (G) The amount of Notch intracellular domain (NICD) was significantly reduced when ZIKV NS4A was overexpressed using nubbin-Gal4 (nub>NS4A), compared to nubbin-Gal4 alone. Notch protein level was marked with anti-NICD staining (red in merge channel and gray in the separate channel). In all images, cytoarchitecture was marked with Actin (green). Scale bars: 100 μm.

We then investigated the status of Notch signaling and found that ZIKV NS4A overexpression (nub>NS4A) substantially reduces Wg and Cut expression in the pouch domain of the wing imaginal epithelia, where wg and cut are the morphogens and are considered the targets of Notch signaling (Neumann and Cohen, 1996; Giraldez and Cohen, 2003) (Fig. 8D,E). We further validated the status of Notch signaling using the reporter NRE-GFP, which consists of a Notch responsive element (NRE) fused with a GFP reporter (Saj et al., 2010). The expression of GFP is restricted to the dorsal/ventral (DV) boundary and its expression thus reliably marks the cells with active Notch signaling. Corroborating with reduced expression of Wg and Cut, we found that NRE-GFP was significantly reduced when ZIKV NS4A was overexpressed in the nubbin region (marks the wing pouch encompassing the DV boundary) of the wing imaginal epithelia (Fig. 8F). To determine if Notch is reduced by overexpression of NS4A, we examined the expression level of Notch protein in wing disc carrying nubbin-specific overexpression of ZIKV NS4A. Antibody staining against the Notch receptor, Notch intracellular domain (NICD), showed that, similar to the expression pattern of Cut, Wg and NRE-GFP, there was reduced expression of Notch protein when ZIKV NS4A was overexpressed compared to nubbin-Gal4 alone (Fig. 8G).

Host-virus interaction and its effect on host physiology has so far remained mostly unexplored. Recently, studies based on transcriptomic and proteomic analyses have been instrumental in advancing our understanding of how viruses can affect host physiology (Merkling et al., 2015; Sun et al., 2017; Dong et al., 2017). Our results further emphasize the applicability of such transcriptomic analyses in understanding host-virus interactions and their molecular regulators. Here, we have shown that infection of adult Drosophila with the MR766 strain of ZIKV results in differential expression of genes implicated in crucial biological processes. Furthermore, analyzing the RNA-seq by putting a stringent filter cutoff resulted in just a handful of biological processes, whereas the number of biological processes increased considerably only when a modest filter cutoff was employed.

Our results show that ZIKV infection induces some of the generic stress-mediated responses, including enrichment of genes related to response to oxidative stress, upregulation of the Turandot family of genes and induction of Diedel (Harsh et al., 2018). However, ZIKV infection does not trigger the activity of heat shock proteins, unlike DCV infection (Merkling et al., 2015).

Infertility is considered another infection-induced pathology where in order to sustain the fitness of an individual, the reproducing ability is reduced (Shirasu-Hiza and Schneider, 2007). Infection with Flock house virus causes oocyte destruction, marked by degenerating egg chambers, disorganized posterior follicle cells and reduction in fecundity (Thomson et al., 2012). Our RNA-seq analysis shows that the genes associated with vitelline membrane and eggshell assembly – such as Cp36, nudel, Femcoat, dec-1 and Cp7Fc – are downregulated upon ZIKV infection, while others have reported that mutation in Cp36 and dec-1 results in female sterility (Hawley and Waring, 1988; Mauzy-Melitz and Waring, 2003).

In mammals, interferon-mediated JAK/STAT signaling is considered to be antiviral in nature. Various strategies of IFN-induced JAK/STAT antagonism have been shown for flaviviruses. Dengue virus infection results in loss of STAT2 expression, while West Nile virus infection results in failed JAK activation in an attempt to evade JAK/STAT-mediated immunity (Guo et al., 2005; Jones et al., 2005). In vitro studies indicate that, similar to dengue virus, ZIKV NS5 binds to STAT2 and targets it for degradation (Grant et al., 2016). Accumulating evidence indicates that the contribution of JAK/STAT signaling is virus specific in in vivo models like Drosophila and mosquitos. Flies mutant for JAK are more sensitive to cricket paralysis virus (CrPV) and DCV; however, these mutant flies show a rather weak phenotype to Drosophila X virus (DXV), Sindbis virus (SINV) and vesicular stomatitis virus (VSV) (Kemp et al., 2013; Dostert et al., 2005). Whereas activation of JAK/STAT restricts dengue virus infection, it does not impart resistance in Aedes aegypti to ZIKV or Chikungunya virus infection (Jupatanakul et al., 2017; Souza-Neto et al., 2009).

In the case of ZIKV infection in Drosophila, a recent finding indicates that one of the negative regulators of JAK/STAT signaling, Diedel, is significantly enriched in ZIKV-infected flies (Harsh et al., 2018). The upregulation of Diedel throughout the stages of infection could be a potential strategy of ZIKV to evade JAK/STAT-mediated immune response. Our findings showing the upregulation of negative regulators Eye transformer (Et) and Socs36E, in addition to reduced level of JAK/STAT reporter, 10XStat92E-GFP, highlight the complex nature of this interaction and indicate that ZIKV infection in Drosophila also inhibits the development-induced JAK/STAT signaling. Et is structurally related to Drosophila JAK/STAT receptor, Domeless, and Et knockdown triggers hyperactivation of septic injury-induced JAK/STAT targets and enhances unpaired-induced eye overgrowth (Kallio et al., 2010). The Suppressor of cytokine signaling (Socs) genes are one of the best-characterized JAK/STAT pathway negative regulators and constitute one of the first feedback loops identified in JAK/STAT signaling (Kile and Alexander, 2001). Overexpression of Socs36E, in particular, phenocopies the outstretched phenotype of upd mutants and the venation defects associated with Stat allele (Rawlings et al., 2004).

The intracellular nature of infection and the complexities linked with host-virus interaction has made it difficult to decipher the pathways responsible for symptoms or disease. Drosophila and its strong genetic tools enable a spatial- and temporal-based genetic manipulation and therefore have emerged as an excellent model system to dissect host-virus interactions (Hughes et al., 2012). Expression of a viral transgene in a tissue-specific manner has led to the deciphering of the functions of several viral components in terms of interaction with host signaling pathways or factors (Hughes et al., 2012). Flavivirus non-structural proteins are considered to be the most critical for their replication and assembly. The host and viral protease-mediated cleavage of flaviviral polyprotein results in seven non-structural proteins, which share structural and functional similarity among the different members of the family (Apte-Sengupta et al., 2014; Bollati et al., 2010). fNSCs expressing ZIKV NS4A and NS4B display impaired neurogenesis and elevated autophagy, which in turn has been linked with inhibited Akt-mTOR signaling (Liang et al., 2016). In line with these findings, our results show that eye-specific expression of ZIKV non-structural proteins, specifically NS4A, results in markedly restricted size of the developing eye. The consistency of the eye phenotype throughout the larval and the adult stage strongly indicates that ZIKV non-structural proteins interact with the regulatory proteins involved in Drosophila eye development. Indeed, our findings validating the reduced level of JAK/STAT signaling and genetic interaction between ZIKV NS4A and JAK/STAT signaling components further imply that ZIKV pathogenesis is linked with eye development, which in turn is associated with regulation of JAK/STAT signaling. These observations also become more relevant in the light of evidence that, unlike other flaviviruses, ZIKV can also be found in the eyes, reproductive organs and bodily secretions, including saliva and urine (Gourinat et al., 2015; Musso et al., 2015; Sun et al., 2016).

As intracellular pathogens, viruses devise strategies to ensure their own proliferation and persistence at the cost of host growth. Some viruses facilitate cell proliferation while others are known to inhibit cell proliferation and in turn result in infection-induced pathologies. Human T cell leukemia virus and human papilloma viruses, for example, encode proteins that promote cell proliferation and result in oncogenic transformation (Op De Beeck and Caillet-Fauquet, 1997). HIV viral protein R (Vpr) is one of the best-characterized proteins known to inhibit T-cell clonal expansion and support the replication of the virus (Poon et al., 1998; Zhao and Elder, 2005; Zhao et al., 2011). Emerging evidence indicates that ZIKV can infect neural precursor cells and subsequently inhibit proliferation and induce apoptosis (Wang et al., 2018; Li et al., 2016a,b). The effects of ZIKV on cell proliferation have also been suggested to be the underlying mechanism for ZIKV-induced pathologies, including microcephaly. Our data depicting the impaired proliferation in the developing eye overexpressing ZIKV transgene further stress the interaction between ZIKV and host cell proliferation. The reduced rate of proliferation could also be the cellular basis for the ZIKV-induced restricted growth of the eye. Future investigations will focus on how reduced cell proliferation benefit ZIKV replication and result in adverse effects to the host.

In the context of neuropathology and ZIKV infection, a recent study demonstrated that overexpression of ZIKV NS4A results in reduced size of larval brain in Drosophila (Shah et al., 2018). It was further shown that ZIKV-induced microcephaly is mediated through human Ankryin repeat and LEM domain-containing 2 (ANKLE2). Interestingly, Drosophila eye epithelia are attached to the optic lobes of the brain. During larval development, while the morphogenetic furrow is developing in the eye, the receptor neurons are specified (Freeman, 1997; Silies et al., 2007). These neurons project their axons to the posterior end of the eye disc, resulting in the optic stalk, which in turn connects the eye disc with optic lobes of the brain. Mutation in the gene disco, responsible for proper formation of the optic stalk, results in defective innervation of neurons from the eye disc to the brain (Steller et al., 1987; Holmes et al., 1998). Given the close proximity of the eye and brain lobes, and our observation that there is an increased accumulation of cell death in the optic stalk (connects the brain and developing eye) upon overexpression of ZIKV NS4A, further suggests that ZIKV transgene expression in Drosophila can also result in neurological abnormalities. Future investigations could provide insight into the neurological pathologies of ZIKV transgene expression.

Our results further show that, in addition to the eye, ZIKV transgene expression can also affect Drosophila wing and regulate its development and patterning. Our findings indicate that ZIKV infection-induced growth regulation is maintained in Drosophila wing as well. However, unlike the eye, the outcome of ZIKV infection is compartment specific in wing. While ZIKV NS4A overexpression does not affect the growth of the posterior compartment of the wing, ZIKV NS4A expression in the pouch domain (encompassing the dorsal/ventral compartment) results in the restricted growth of this domain. The subdivision of the wing into anterior/posterior (a/p) compartments occurs in the embryonic stage, while the DV boundary is established during proliferation of the imaginal disc (Garcia-Bellido et al., 1973). ZIKV effect on the pouch domain without affecting the posterior compartment thus emphasizes the regulation of cell proliferation by ZIKV-induced infection. In Drosophila wing, Notch signaling is crucial in regulating cell differentiation and patterning. Notch signaling, in particular, regulates expression of several genes involved in the DV boundary and elicits wing cell proliferation (Irvine and Vogt, 1997). Our finding that ZIKV transgene-mediated effect on wing growth is regulated via reduced Notch signaling is suggestive of the role of ZIKV-induced infection in regulating cell proliferation in the wing imaginal epithelia as well.

Taken altogether, our results illustrate the efficacy of Drosophila as a model for understanding the pathology resulting from host-ZIKV interactions. Our findings reveal that ZIKV plays a crucial role in regulating the growth of Drosophila wing and eye and results in their reduced size. Whereas ZIKV-induced pathogenesis in the eye is linked with reduced JAK/STAT signaling, in the wing it correlates with reduced level of Notch signaling. The key findings presented in this study will assist in unraveling ZIKV-induced pathogenesis and will advance our understanding of host-ZIKV interactions.

Fly stocks

The following fly lines were used: w1118 (wild-type control), E1-Gal4 (Pallavi et al., 2012), engrailed Gal4UAS-RFP (Bloomington no. 30557), UAS-NS4A, UAS-NS4B, UAS-NS2A, UAS-NS2B, UAS-StatRNAi (Vienna Drosophila Resource Center no. 43866), w1118; P(10XStat92E-GFP) (Bloomington no. 26197), w1118; P(NRE-EGFP.S)5A (Bloomington no. 30727), nubbin-Gal4 (Bloomington no. 25754). UAS-upd and UAS-HopTuml were generously provided by Heinrich Jasper (Buck Institute for Research and Aging, Novato, CA, USA), and UAS-domeDN was provided by Dan Hultmark (Umeå University, Umeå, Sweden). Flies were reared on instant Drosophila diet (Formula 4-24 Drosophila medium) supplemented with yeast (Carolina Biological Supply), and maintained at 25°C and in a 12:12-h light:dark photoperiodic cycle. Female adult flies aged 4-6 days old were used in infection assays with ZIKV.

DNA cloning and generation of transgenic fly strains

Four ZIKV open reading frames (ORFs) (NS2A, NS2B, NS4A and NS4B) in mammalian expression vectors were gifts from Dr Hongjun Liu's laboratory in Johns Hopkins University (Yoon et al., 2017). To generate Flag-tagged UAS-ZIKV-ORF constructs, each ZIKV ORF was PCR amplified with specific primers using mammalian expression vector as template. Each amplicon was inserted into the pUASTattB vector with EcoRI and KpnI restriction enzymes and introduced into the germ cells of flies by standard P element-mediated germ line transformation. Primers used for PCR were as follows: ZIKV-NS2A: forward, 5′-CCGGAATTCATGGATTACAAGGATGACGACGATAAGGGGTCAACCGATCATATGGAC-3′; reverse, 5′-CGGGGTACCCTACCGCTTCCCACTCCTTGTGAG-3′. ZIKV-NS2B: forward, 5′-CCGGAATTCATGGATTACAAGGATGACGACGATAAGAGCTGGCCCCCTAGTGAAGTTC-3′; reverse, 5′-CGGGGTACCCTACCTTTTCCCAGTCTTCACATAC-3′. ZIKV-NS4A: forward, prime;-CCGGAATTCATGGATTACAAGGATGACGACGATAAGGGAGCGGCTTTGGGAGTAATG-3′; reverse, 5′-CGGGGTACCCTATCTTTGCTTCTCTGGCTCGGG-3′. ZIKV-NS4B: forward, 5′-CCGGAATTCATGGATTACAAGGATGACGACGATAAGAACGAACTTGGATGGCTGGAAAG-3′; reverse, 5′-CGGGGTACCCTAACGTCTCTTAACCAGGCCAGC-3′. Bold italic sequences are restriction enzyme cutting sites and underlined sequences are Flag tag.

ZIKV infection

The ZIKV infection protocol and propagation have been described in detail previously. Briefly, the African strain of ZIKV, MR766, was propagated in Vero cells followed by determination of ZIKV titers using plaque assay on Vero cells as described in Delvecchio et al. (2016). For fly infection, adult female flies were intrathoracically injected with 100 nl live ZIKV solution [11,000 plaque-forming units (PFU)/fly] using a nanoinjector (Nanoject III, Drummond Scientific). Injection of the same volume of PBS served as a negative control. Injected flies were kept at 25°C and transferred to fresh vials every third day. They were collected at the different timepoints and directly processed for RNA analysis.

RNA isolation

Total RNA was extracted from ten adult female flies injected with ZIKV and PBS at 4 and 8 dpi in triplicate using Trizol according to the manufacturer's protocol. Total RNA was re-suspended in 30 μl sterile nuclease-free water. RNA concentration was measured using Nanodrop. RNA integrity and quality were assessed on formaldehyde agarose gel.

Library preparation and RNA-seq

TruSeq Total RNA with Ribo-Zero Gold Kit (Illumina, San Diego, CA, USA) was used to generate strand-specific RNA-seq libraries following the manufacturer's recommendations. The libraries were pooled and sequenced in a single lane on an Illumina HiSeq 2500 sequencing instrument with single-end 50 bp reads. RNA-seq was performed at The Genome Technology Access Center (GTAC) based in Washington University in St Louis, MO, USA.

RNA-seq reads were aligned to the Ensembl release 76 top-level assembly with STAR version 2.0.4b (Dobin et al., 2013). Gene counts were derived from the number of uniquely aligned unambiguous reads by Subread: feature Count version 1.4.5 (Liao et al., 2014). Transcript counts were produced by Sailfish version 0.6.3 (Patro et al., 2014). Sequencing performance was assessed for total number of aligned reads, total number of uniquely aligned reads, genes and transcripts detected, ribosomal fraction known junction saturation and read distribution over known gene models with RSeQC version 2.3 (Wang et al., 2012). All RNA-seq data have been deposited in NCBI Gene Expression Omnibus (GEO) with the accession number GSE130108.

Bioinformatics analysis

All gene counts were then imported into the R/Bioconductor package EdgeR (Robinson et al., 2010) and TMM normalization size factors were calculated to adjust for samples for differences in library size across samples after ribosomal genes. Genes not expressed greater than 1 count per million in at least two samples were removed from further analysis. The TMM size factors and the matrix of counts were then imported into R/Bioconductor package Limma and weighted likelihoods based on the observed mean-variance relationship of every gene/transcript and sample were calculated for all samples with the voomWithQualityWeights function (Liu et al., 2015). Generalized linear models were then created to test for gene/transcript level differential expression. DEGs and transcripts were then filtered for FDR adjusted P-values less than or equal to 0.1.

GO analysis

GO analysis was performed using DAVID 6.8 bioinformatics resources (http://david.abcc.ncifcrf.gov/). In all cases, analyses were performed using the list of DEGs. The P-value cutoff to determine enriched pathways was 0.1.

qRT-PCR validation

To validate DEGs, we selected eight candidate genes based on significant fold differences and analyzed their mRNA levels using qRT-PCR. Ten adult female flies injected with PBS or ZIKV were frozen at 4 and 8 days dpi. Total RNA was extracted using Trizol according to the manufacturer's protocol. Total RNA (350-500 ng) was used to synthesize cDNA using a High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). qRT-PCR experiments were performed in technical triplicates and with gene-specific primers (Table 1) using iQ SYBR Green Supermix (Bio-Rad) and a CFX96 Real-Time PCR detection system (Bio-Rad). Quantification was performed from three biological replicates for both test and control treatments. Fold changes were calculated with the delta-delta Ct method using RpL32 as a housekeeping gene.

Table 1.

List of the gene-specific primers used in the study

List of the gene-specific primers used in the study
List of the gene-specific primers used in the study

Fly survival

For each fly strain, three groups of 20 female flies carrying eyeless-specific overexpression of non-structural proteins were used. Flies were maintained at a constant temperature of 25°C with a 12-h light/dark cycle, and mortality was recorded daily. Log-rank (Mantel–Cox) test was used to analyze the survival curves.

Pupariation delay

Two-hour timed egg collections were carried out and L1 larvae were collected 24 h after egg deposition and reared at 30 animals per vial. The number of larvae that had pupariated at a given time after egg deposition was scored every 8 h.

Immunostaining and antibodies

Anti-PH3 (1:500; Abcam), anti-CycA [1:50; Developmental Studies Hybridoma Bank (DSHB)], anti-Dcp-1 (1:100; Cell Signaling Technology), anti-Cut (1:10; DSHB), anti-Wg (1:1000; DSHB), anti-Elav (1:50; DSHB) and anti-NICD (1:10; DSHB) were used as the source for primary antibodies. Secondary antibodies included Alexa Fluor 488, 555 and 633 (Invitrogen). For nuclear staining, 4′,6-diamidino-2-phenylindole (DAPI; Invitrogen) was used. Standard procedures were followed for immunostaining. The crosses were properly synchronized for the stages and third-instar larvae were collected for dissections. Briefly, fly tissues were dissected and fixed in PBS containing 4% formaldehyde for 30 min. Following double rinsing in PBS containing 0.1% Triton X-100, the samples were incubated overnight at 4°C with the primary antibody. The samples were then blocked with 1% BSA for 2 h followed by a 2-h incubation with secondary antibody at room temperature. Finally, the samples were mounted with Vectashield medium (Vector Laboratories). Images were acquired with a Zeiss LSM 510 confocal microscope and processed using Adobe Photoshop CS6.

Fluorescence quantification

Relative amounts of fluorescence were measured with ImageJ software by using Shanbhag thresholding on images and calculating the integrated density, resulting area and mean fluorescence of the background. The following equation was used: corrected total fluorescence=integrated density−(area×mean fluorescence of background).

Larval wing imaginal disc and adult wing size

The size of the third-instar wing imaginal discs and adult wings were measured using ImageJ. In the case of the wing pouch, the size in different genotypes was normalized to the total size of the wing disc.

Adult eye imaging

The photographs and measurements were acquired using a Keyence VHX 5000 Digital Microscope.

Statistical analysis

qRT-PCR results represent the means±s.d. of relative values from three biological replicates. An unpaired two-tailed Student's t-test was used for statistical analysis of data using Prism (GraphPad Software).

We thank the Genome Technology Access Center in the Department of Genetics at Washington University School of Medicine for help with genomic analysis. The Center is partially supported by NCI Cancer Center Support Grant #P30 CA91842 to the Siteman Cancer Center and by ICTS/CTSA Grant #UL1TR000448 from the National Center for Research Resources (NCRR), a component of the National Institutes of Health (NIH), and NIH Roadmap for Medical Research. This publication is solely the responsibility of the authors and does not necessarily represent the official view of the NCRR or NIH. We thank Dr Alex Jeremic's laboratory [The George Washington University (GWU), Department of Biological Sciences] for assistance with confocal microscopy. We thank Dr Arnaud Martin (GWU, Department of Biological Sciences) for assistance with adult fly eye imaging. We also thank members of the Department of Biological Sciences at GWU for critical reading of the manuscript.

Author contributions

Conceptualization: S.H., I.E.; Methodology: S.H., Y.F., E.K., Z.H.; Software: S.H., E.K.; Validation: S.H.; Formal analysis: S.H.; Investigation: S.H.; Resources: Y.F., Z.H.; Data curation: S.H., E.K.; Writing - original draft: S.H., I.E.; Writing - review & editing: S.H., I.E.; Supervision: I.E.; Project administration: I.E.; Funding acquisition: I.E.

Funding

This work was supported by the National Institute of Allergy and Infectious Diseases (1R01AI110675–01A1 to I.E.).

Data availability

RNA-seq data are available at GEO with the accession number GSE130108.

Agaisse
,
H.
,
Petersen
,
U.-M.
,
Boutros
,
M.
,
Mathey-Prevot
,
B.
and
Perrimon
,
N.
(
2003
).
Signaling role of hemocytes in Drosophila JAK/STAT-dependent response to septic injury
.
Dev. Cell
5
,
441
-
450
.
Apte-Sengupta
,
S.
,
Sirohi
,
D.
and
Kuhn
,
R. J.
(
2014
).
Coupling of replication and assembly in flaviviruses
.
Curr. Opin. Virol.
9
,
134
-
142
.
Arbouzova
,
N. I.
and
Zeidler
,
M. P.
(
2006
).
JAK/STAT signalling in Drosophila: insights into conserved regulatory and cellular functions
.
Development
133
,
2605
-
2616
.
Bach
,
E. A.
,
Vincent
,
S.
,
Zeidler
,
M. P.
and
Perrimon
,
N.
(
2003
).
A sensitized genetic screen to identify novel regulators and components of the Drosophila janus kinase/signal transducer and activator of transcription pathway
.
Genetics
165
,
1149
-
1166
.
Bach
,
E. A.
,
Ekas
,
L. A.
,
Ayala-Camargo
,
A.
,
Flaherty
,
M. S.
,
Lee
,
H.
,
Perrimon
,
N.
and
Baeg
,
G.-H.
(
2007
).
GFP reporters detect the activation of the Drosophila JAK/STAT pathway in vivo
.
Gene Expr. Patterns
7
,
323
-
331
.
Bernards
,
A.
and
Hariharan
,
I. K.
(
2001
).
Of flies and men--studying human disease in Drosophila
.
Curr. Opin. Genet. Dev.
11
,
274
-
278
.
Best
,
S. M.
(
2017
).
The many faces of the flavivirus NS5 protein in antagonism of type I interferon signaling
.
J. Virol.
91
,
e01970-16
.
Bollati
,
M.
,
Alvarez
,
K.
,
Assenberg
,
R.
,
Baronti
,
C.
,
Canard
,
B.
,
Cook
,
S.
,
Coutard
,
B.
,
Decroly
,
E.
,
De Lamballerie
,
X.
,
Gould
,
E. A.
, et al. 
(
2010
).
Structure and functionality in flavivirus NS-proteins: perspectives for drug design
.
Antiviral Res.
87
,
125
-
148
.
Brown
,
S.
,
Hu
,
N.
and
Hombría
,
J. C.-G.
(
2001
).
Identification of the first invertebrate interleukin JAK/STAT receptor, the Drosophila gene domeless
.
Curr. Biol.
11
,
1700
-
1705
.
Callus
,
B. A.
and
Mathey-Prevot
,
B.
(
2002
).
SOCS36E, a novel Drosophila SOCS protein, suppresses JAK/STAT and EGF-R signalling in the imaginal wing disc
.
Oncogene
21
,
4812
-
4821
.
Chambers
,
T. J.
,
Hahn
,
C. S.
,
Galler
,
R.
and
Rice
,
C. M.
(
1990
).
Flavivirus genome organization, expression, and replication
.
Annu. Rev. Microbiol.
44
,
649
-
688
.
Chan
,
C.-M.
,
Ma
,
C.-W.
,
Chan
,
W.-Y.
and
Chan
,
H. Y. E.
(
2007
).
The SARS-coronavirus membrane protein induces apoptosis through modulating the Akt survival pathway
.
Arch. Biochem. Biophys.
459
,
197
-
207
.
Chen
,
S.
,
Wu
,
Z.
,
Wang
,
M.
and
Cheng
,
A.
(
2017
).
Innate immune evasion mediated by flaviviridae non-structural proteins
.
Viruses
9
,
E291
.
Cifuentes
,
F. J.
and
Garcia-Bellido
,
A.
(
1997
).
Proximo-distal specification in the wing disc of Drosophila by the nubbin gene
.
Proc. Natl. Acad. Sci. USA
94
,
11405
-
11410
.
Dalrymple
,
N. A.
,
Cimica
,
V.
and
Mackow
,
E. R.
(
2015
).
Dengue virus NS proteins inhibit RIG-I/MAVS signaling by blocking TBK1/IRF3 phosphorylation: dengue virus serotype 1 NS4A is a unique interferon-regulating virulence determinant
.
MBio
6
,
e00553-15
.
Delvecchio
,
R.
,
Higa
,
L. M.
,
Pezzuto
,
P.
,
Valadao
,
A. L.
,
Garcez
,
P. P.
,
Monteiro
,
F. L.
,
Loiola
,
E. C.
,
Dias
,
A. A.
,
Silva
,
F. J.
,
Aliota
,
M. T.
, et al. 
(
2016
).
Chloroquine, an endocytosis blocking agent, inhibits Zika virus infection in different cell models
.
Viruses
8
,
E322
.
Ding
,
Q.
,
Cao
,
X.
,
Lu
,
J.
,
Huang
,
B.
,
Liu
,
Y.-J.
,
Kato
,
N.
,
Shu
,
H.-B.
and
Zhong
,
J.
(
2013
).
Hepatitis C virus NS4B blocks the interaction of STING and TBK1 to evade host innate immunity
.
J. Hepatol.
59
,
52
-
58
.
do Rosario
,
M. S.
,
De Jesus
,
P. A.
,
Vasilakis
,
N.
,
Farias
,
D. S.
,
Novaes
,
M. A.
,
Rodrigues
,
S. G.
,
Martins
,
L. C.
,
Vasconcelos
,
P. F.
,
Ko
,
A. I.
,
Alcantara
,
L. C.
and et al.  (
2016
).
Guillain-barre syndrome after Zika virus infection in Brazil
.
Am. J. Trop. Med. Hyg.
95
,
1157
-
1160
.
Dobin
,
A.
,
Davis
,
C. A.
,
Schlesinger
,
F.
,
Drenkow
,
J.
,
Zaleski
,
C.
,
Jha
,
S.
,
Batut
,
P.
,
Chaisson
,
M.
and
Gingeras
,
T. R.
(
2013
).
STAR: ultrafast universal RNA-seq aligner
.
Bioinformatics
29
,
15
-
21
.
Dong
,
S.
,
Behura
,
S. K.
and
Franz
,
A. W. E.
(
2017
).
The midgut transcriptome of Aedes aegypti fed with saline or protein meals containing chikungunya virus reveals genes potentially involved in viral midgut escape
.
BMC Genomics
18
,
382
.
dos Santos
,
T.
,
Rodriguez
,
A.
,
Almiron
,
M.
,
Sanhueza
,
A.
,
Ramon
,
P.
,
De Oliveira
,
W. K.
,
Coelho
,
G. E.
,
Badaró
,
R.
,
Cortez
,
J.
,
Ospina
,
M.
, et al. 
(
2016
).
Zika virus and the guillain-barre syndrome - case series from seven countries
.
N. Engl. J. Med.
375
,
1598
-
1601
.
Dostert
,
C.
,
Jouanguy
,
E.
,
Irving
,
P.
,
Troxler
,
L.
,
Galiana-Arnoux
,
D.
,
Hetru
,
C.
,
Hoffmann
,
J. A.
and
Imler
,
J.-L.
(
2005
).
The Jak-STAT signaling pathway is required but not sufficient for the antiviral response of drosophila
.
Nat. Immunol.
6
,
946
-
953
.
Edgar
,
B. A.
and
Lehner
,
C. F.
(
1996
).
Developmental control of cell cycle regulators: a fly's perspective
.
Science
274
,
1646
-
1652
.
Ekas
,
L. A.
,
Baeg
,
G.-H.
,
Flaherty
,
M. S.
,
Ayala-Camargo
,
A.
and
Bach
,
E. A.
(
2006
).
JAK/STAT signaling promotes regional specification by negatively regulating wingless expression in Drosophila
.
Development
133
,
4721
-
4729
.
Fauci
,
A. S.
and
Morens
,
D. M.
(
2016
).
Zika virus in the americas--yet another arbovirus threat
.
N. Engl. J. Med.
374
,
601
-
604
.
Flaherty
,
M. S.
,
Zavadil
,
J.
,
Ekas
,
L. A.
and
Bach
,
E. A.
(
2009
).
Genome-wide expression profiling in the Drosophila eye reveals unexpected repression of notch signaling by the JAK/STAT pathway
.
Dev. Dyn.
238
,
2235
-
2253
.
Foy
,
B. D.
,
Kobylinski
,
K. C.
,
Foy
,
J. L. C.
,
Blitvich
,
B. J.
,
Travassos Da Rosa
,
A.
,
Haddow
,
A. D.
,
Lanciotti
,
R. S.
and
Tesh
,
R. B.
(
2011
).
Probable non-vector-borne transmission of Zika virus, Colorado, USA
.
Emerg. Infect. Dis.
17
,
880
-
882
.
Freeman
,
M.
(
1997
).
Cell determination strategies in the Drosophila eye
.
Development
124
,
261
-
270
.
Furtado
,
J. M.
,
Espósito
,
D. L.
,
Klein
,
T. M.
,
Teixeira-Pinto
,
T.
and
da Fonseca
,
B. A.
(
2016
).
Uveitis associated with Zika virus infection
.
N. Engl. J. Med.
375
,
394
-
396
.
Garcia-Bellido
,
A.
,
Ripoll
,
P.
and
Morata
,
G.
(
1973
).
Developmental compartmentalisation of the wing disk of Drosophila
.
Nat. New Biol.
245
,
251
-
253
.
Giraldez
,
A. J.
and
Cohen
,
S. M.
(
2003
).
Wingless and Notch signaling provide cell survival cues and control cell proliferation during wing development
.
Development
130
,
6533
-
6543
.
Gourinat
,
A.-C.
,
O'connor
,
O.
,
Calvez
,
E.
,
Goarant
,
C.
and
Dupont-Rouzeyrol
,
M.
(
2015
).
Detection of Zika virus in urine
.
Emerg. Infect. Dis.
21
,
84
-
86
.
Grant
,
A.
,
Ponia
,
S. S.
,
Tripathi
,
S.
,
Balasubramaniam
,
V.
,
Miorin
,
L.
,
Sourisseau
,
M.
,
Schwarz
,
M. C.
,
Sánchez-Seco
,
M. P.
,
Evans
,
M. J.
,
Best
,
S. M.
, et al. 
(
2016
).
Zika virus targets human STAT2 to inhibit type I interferon signaling
.
Cell Host Microbe.
19
,
882
-
890
.
Guo
,
J.-T.
,
Hayashi
,
J.
and
Seeger
,
C.
(
2005
).
West Nile virus inhibits the signal transduction pathway of alpha interferon
.
J. Virol.
79
,
1343
-
1350
.
Harsh
,
S.
,
Ozakman
,
Y.
,
Kitchen
,
S. M.
,
Paquin-Proulx
,
D.
,
Nixon
,
D. F.
and
Eleftherianos
,
I.
(
2018
).
Dicer-2 regulates resistance and maintains homeostasis against Zika virus infection in Drosophila
.
J. Immunol.
201
,
3058
-
3072
.
Hasan
,
S. S.
,
Sevvana
,
M.
,
Kuhn
,
R. J.
and
Rossmann
,
M. G.
(
2018
).
Structural biology of Zika virus and other flaviviruses
.
Nat. Struct. Mol. Biol.
25
,
13
-
20
.
Hawley
,
R. J.
and
Waring
,
G. L.
(
1988
).
Cloning and analysis of the dec-1 female-sterile locus, a gene required for proper assembly of the Drosophila eggshell
.
Genes Dev.
2
,
341
-
349
.
Holmes
,
A. L.
,
Raper
,
R. N.
and
Heilig
,
J. S.
(
1998
).
Genetic analysis of Drosophila larval optic nerve development
.
Genetics
148
,
1189
-
1201
.
Huang da
,
W.
,
Sherman
,
B. T.
and
Lempicki
,
R. A.
(
2009a
).
Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists
.
Nucleic Acids Res.
37
,
1
-
13
.
Huang da
,
W.
,
Sherman
,
B. T.
and
Lempicki
,
R. A.
(
2009b
).
Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources
.
Nat. Protoc.
4
,
44
-
57
.
Hughes
,
T. T.
,
Allen
,
A. L.
,
Bardin
,
J. E.
,
Christian
,
M. N.
,
Daimon
,
K.
,
Dozier
,
K. D.
,
Hansen
,
C. L.
,
Holcomb
,
L. M.
and
Ahlander
,
J.
(
2012
).
Drosophila as a genetic model for studying pathogenic human viruses
.
Virology
423
,
1
-
5
.
Irvine
,
K. D.
and
Vogt
,
T. F.
(
1997
).
Dorsal-ventral signaling in limb development
.
Curr. Opin. Cell Biol.
9
,
867
-
876
.
Jones
,
M.
,
Davidson
,
A.
,
Hibbert
,
L.
,
Gruenwald
,
P.
,
Schlaak
,
J.
,
Ball
,
S.
,
Foster
,
G. R.
and
Jacobs
,
M.
(
2005
).
Dengue virus inhibits alpha interferon signaling by reducing STAT2 expression
.
J. Virol.
79
,
5414
-
5420
.
Jupatanakul
,
N.
,
Sim
,
S.
,
Angleró-Rodriguez
,
Y. I.
,
Souza-Neto
,
J.
,
Das
,
S.
,
Poti
,
K. E.
,
Rossi
,
S. L.
,
Bergren
,
N.
,
Vasilakis
,
N.
and
Dimopoulos
,
G.
(
2017
).
Engineered aedes aegypti JAK/STAT pathway-mediated immunity to dengue virus
.
PLoS Negl. Trop. Dis.
11
,
e0005187
.
Kallio
,
J.
,
Myllymäki
,
H.
,
Grönholm
,
J.
,
Armstrong
,
M.
,
Vanha-Aho
,
L. M.
,
Mäkinen
,
L.
,
Silvennoinen
,
O.
,
Valanne
,
S.
and
Rämet
,
M.
(
2010
).
Eye transformer is a negative regulator of Drosophila JAK/STAT signaling
.
FASEB J.
24
,
4467
-
4479
.
Kemp
,
C.
,
Mueller
,
S.
,
Goto
,
A.
,
Barbier
,
V.
,
Paro
,
S.
,
Bonnay
,
F.
,
Dostert
,
C.
,
Troxler
,
L.
,
Hetru
,
C.
,
Meignin
,
C.
, et al. 
(
2013
).
Broad RNA interference-mediated antiviral immunity and virus-specific inducible responses in Drosophila
.
J. Immunol.
190
,
650
-
658
.
Kile
,
B. T.
and
Alexander
,
W. S.
(
2001
).
The suppressors of cytokine signalling (SOCS)
.
Cell Mol. Life Sci.
58
,
1627
-
1635
.
Kumthip
,
K.
,
Chusri
,
P.
,
Jilg
,
N.
,
Zhao
,
L.
,
Fusco
,
D. N.
,
Zhao
,
H.
,
Goto
,
K.
,
Cheng
,
D.
,
Schaefer
,
E. A.
,
Zhang
,
L.
, et al. 
(
2012
).
Hepatitis C virus NS5A disrupts STAT1 phosphorylation and suppresses type I interferon signaling
.
J. Virol.
86
,
8581
-
8591
.
Laurent-Rolle
,
M.
,
Boer
,
E. F.
,
Lubick
,
K. J.
,
Wolfinbarger
,
J. B.
,
Carmody
,
A. B.
,
Rockx
,
B.
,
Liu
,
W.
,
Ashour
,
J.
,
Shupert
,
W. L.
,
Holbrook
,
M. R.
, et al. 
(
2010
).
The NS5 protein of the virulent West Nile virus NY99 strain is a potent antagonist of type I interferon-mediated JAK-STAT signaling
.
J. Virol.
84
,
3503
-
3515
.
Lee
,
S. B.
,
Park
,
J.
,
Jung
,
J. U.
and
Chung
,
J.
(
2005
).
Nef induces apoptosis by activating JNK signaling pathway and inhibits NF-kappaB-dependent immune responses in Drosophila
.
J. Cell Sci.
118
,
1851
-
1859
.
Li
,
C.
,
Xu
,
D.
,
Ye
,
Q.
,
Hong
,
S.
,
Jiang
,
Y.
,
Liu
,
X.
,
Zhang
,
N.
,
Shi
,
L.
,
Qin
,
C.-F.
and
Xu
,
Z.
(
2016a
).
Zika virus disrupts neural progenitor development and leads to microcephaly in mice
.
Cell Stem Cell
19
,
120
-
126
.
Li
,
H.
,
Saucedo-Cuevas
,
L.
,
Regla-Nava
,
J. A.
,
Chai
,
G.
,
Sheets
,
N.
,
Tang
,
W.
,
Terskikh
,
A. V.
,
Shresta
,
S.
and
Gleeson
,
J. G.
(
2016b
).
Zika virus infects neural progenitors in the adult mouse brain and alters proliferation
.
Cell Stem Cell
19
,
593
-
598
.
Liang
,
Q.
,
Luo
,
Z.
,
Zeng
,
J.
,
Chen
,
W.
,
Foo
,
S.-S.
,
Lee
,
S.-A.
,
Ge
,
J.
,
Wang
,
S.
,
Goldman
,
S. A.
,
Zlokovic
,
B. V.
, et al. 
(
2016
).
Zika virus NS4A and NS4B proteins deregulate Akt-mTOR signaling in human fetal neural stem cells to inhibit neurogenesis and induce autophagy
.
Cell Stem Cell
19
,
663
-
671
.
Liao
,
Y.
,
Smyth
,
G. K.
and
Shi
,
W.
(
2014
).
featureCounts: an efficient general purpose program for assigning sequence reads to genomic features
.
Bioinformatics
30
,
923
-
930
.
Link
,
N.
,
Chung
,
H.
,
Jolly
,
A.
,
Withers
,
M.
,
Tepe
,
B.
,
Arenkiel
,
B. R.
,
Shah
,
P. S.
,
Krogan
,
N. J.
,
Aydin
,
H.
,
Geckinli
,
B. B.
, et al. 
(
2019
).
Mutations in ANKLE2, a ZIKA virus target, disrupt an asymmetric cell division pathway in Drosophila neuroblasts to cause microcephaly
.
Dev. Cell
51
,
713
-
729.e6
.
Liu
,
W. J.
,
Chen
,
H. B.
,
Wang
,
X. J.
,
Huang
,
H.
and
Khromykh
,
A. A.
(
2004
).
Analysis of adaptive mutations in Kunjin virus replicon RNA reveals a novel role for the flavivirus nonstructural protein NS2A in inhibition of beta interferon promoter-driven transcription
.
J. Virol.
78
,
12225
-
12235
.
Liu
,
R.
,
Holik
,
A. Z.
,
Su
,
S.
,
Jansz
,
N.
,
Chen
,
K.
,
Leong
,
H. S.
,
Blewitt
,
M. E.
,
Asselin-Labat
,
M.-L.
,
Smyth
,
G. K.
and
Ritchie
,
M. E.
(
2015
).
Why weight? Modelling sample and observational level variability improves power in RNA-seq analyses
.
Nucleic Acids Res.
43
,
e97
.
Liu
,
Y.
,
Gordesky-Gold
,
B.
,
Leney-Greene
,
M.
,
Weinbren
,
N. L.
,
Tudor
,
M.
and
Cherry
,
S.
(
2018
).
Inflammation-induced, STING-dependent autophagy restricts Zika virus infection in the Drosophila brain
.
Cell Host Microbe.
24
,
57
-
68.e3
.
Mauzy-Melitz
,
D.
and
Waring
,
G. L.
(
2003
).
fc177, a minor dec-1 proprotein, is necessary to prevent ectopic aggregation of the endochorion during eggshell assembly in Drosophila
.
Dev. Biol
255
,
193
-
205
.
Mead
,
P. S.
,
Hills
,
S. L.
and
Brooks
,
J. T.
(
2018
).
Zika virus as a sexually transmitted pathogen
.
Curr. Opin. Infect. Dis.
31
,
39
-
44
.
Merkling
,
S. H.
,
Overheul
,
G. J.
,
van Mierlo
,
J. T.
,
Arends
,
D.
,
Gilissen
,
C.
and
van Rij
,
R. P.
(
2015
).
The heat shock response restricts virus infection in Drosophila
.
Sci. Rep.
5
,
12758
.
Miner
,
J. J.
,
Sene
,
A.
,
Richner
,
J. M.
,
Smith
,
A. M.
,
Santeford
,
A.
,
Ban
,
N.
,
Weger-Lucarelli
,
J.
,
Manzella
,
F.
,
Rückert
,
C.
,
Govero
,
J.
, et al. 
(
2016
).
Zika virus infection in mice causes panuveitis with shedding of virus in tears
.
Cell Rep.
16
,
3208
-
3218
.
Mukherjee
,
T.
,
Hombría
,
J. C.-G.
and
Zeidler
,
M. P.
(
2005
).
Opposing roles for Drosophila JAK/STAT signalling during cellular proliferation
.
Oncogene
24
,
2503
-
2511
.
Munoz-Jordan
,
J. L.
,
Sanchez-Burgos
,
G. G.
,
Laurent-Rolle
,
M.
and
Garcia-Sastre
,
A.
(
2003
).
Inhibition of interferon signaling by dengue virus
.
Proc. Natl. Acad. Sci. USA
100
,
14333
-
14338
.
Musso
,
D.
,
Nhan
,
T.
,
Robin
,
E.
,
Roche
,
C.
,
Bierlaire
,
D.
,
Zisou
,
K.
,
Shan Yan
,
A.
,
Cao-Lormeau
,
V. M.
and
Broult
,
J.
(
2014
).
Potential for Zika virus transmission through blood transfusion demonstrated during an outbreak in French Polynesia, November 2013 to February 2014
.
Euro Surveill.
19
,
20761
.
Musso
,
D.
,
Roche
,
C.
,
Nhan
,
T.-X.
,
Robin
,
E.
,
Teissier
,
A.
and
Cao-Lormeau
,
V.-M.
(
2015
).
Detection of Zika virus in saliva
.
J. Clin. Virol.
68
,
53
-
55
.
Neumann
,
C. J.
and
Cohen
,
S. M.
(
1996
).
A hierarchy of cross-regulation involving Notch, wingless, vestigial and cut organizes the dorsal/ventral axis of the Drosophila wing
.
Development
122
,
3477
-
3485
.
Neumann
,
C. J.
and
Cohen
,
S. M.
(
1998
).
Boundary formation in Drosophila wing: Notch activity attenuated by the POU protein Nubbin
.
Science
281
,
409
-
413
.
Ng
,
M.
,
Diaz-Benjumea
,
F. J.
and
Cohen
,
S. M.
(
1995
).
Nubbin encodes a POU-domain protein required for proximal-distal patterning in the Drosophila wing
.
Development
121
,
589
-
599
.
Ngono
,
A. E.
and
Shresta
,
S.
(
2018
).
Immune response to dengue and Zika
.
Annu. Rev. Immunol.
36
,
279
-
308
.
Op De Beeck
,
A.
and
Caillet-Fauquet
,
P
. (
1997
).
Viruses and the cell cycle
.
Prog. Cell Cycle Res.
3
,
1
-
19
.
Pallavi
,
S. K.
,
Ho
,
D. M.
,
Hicks
,
C.
,
Miele
,
L.
and
Artavanis-Tsakonas
,
S.
(
2012
).
Notch and Mef2 synergize to promote proliferation and metastasis through JNK signal activation in Drosophila
.
EMBO J.
31
,
2895
-
2907
.
Pandey
,
U. B.
and
Nichols
,
C. D.
(
2011
).
Human disease models in Drosophila melanogaster and the role of the fly in therapeutic drug discovery
.
Pharmacol. Rev.
63
,
411
-
436
.
Patro
,
R.
,
Mount
,
S. M.
and
Kingsford
,
C.
(
2014
).
Sailfish enables alignment-free isoform quantification from RNA-seq reads using lightweight algorithms
.
Nat. Biotechnol.
32
,
462
-
464
.
Perrimon
,
N.
and
Mahowald
,
A. P.
(
1986
).
l(1)hopscotch, A larval-pupal zygotic lethal with a specific maternal effect on segmentation in Drosophila
.
Dev. Biol.
118
,
28
-
41
.
Poon
,
B.
,
Grovit-Ferbas
,
K.
,
Stewart
,
S. A.
and
Chen
,
I. S.
(
1998
).
Cell cycle arrest by Vpr in HIV-1 virions and insensitivity to antiretroviral agents
.
Science
281
,
266
-
269
.
Rawlings
,
J. S.
,
Rennebeck
,
G.
,
Harrison
,
S. M.
,
Xi
,
R.
and
Harrison
,
D. A.
(
2004
).
Two Drosophila suppressors of cytokine signaling (SOCS) differentially regulate JAK and EGFR pathway activities
.
BMC Cell Biol.
5
,
38
.
Recasens-Alvarez
,
C.
,
Ferreira
,
A.
and
Milán
,
M.
(
2017
).
JAK/STAT controls organ size and fate specification by regulating morphogen production and signalling
.
Nat. Commun.
8
,
13815
.
Ritchie
,
M. E.
,
Phipson
,
B.
,
Wu
,
D.
,
Hu
,
Y.
,
Law
,
C. W.
,
Shi
,
W.
and
Smyth
,
G. K.
(
2015
).
limma powers differential expression analyses for RNA-sequencing and microarray studies
.
Nucleic Acids Res.
43
,
e47
.
Robinson
,
M. D.
,
Mccarthy
,
D. J.
and
Smyth
,
G. K.
(
2010
).
edgeR: a Bioconductor package for differential expression analysis of digital gene expression data
.
Bioinformatics
26
,
139
-
140
.
Saj
,
A.
,
Arziman
,
Z.
,
Stempfle
,
D.
,
Van Belle
,
W.
,
Sauder
,
U.
,
Horn
,
T.
,
Dürrenberger
,
M.
,
Paro
,
R.
,
Boutros
,
M.
and
Merdes
,
G.
(
2010
).
A combined ex vivo and in vivo RNAi screen for notch regulators in Drosophila reveals an extensive notch interaction network
.
Dev. Cell
18
,
862
-
876
.
Schindler
,
C.
,
Shuai
,
K.
,
Prezioso
,
V. R.
and
Darnell
,
J. E.
Jr.
(
1992
).
Interferon-dependent tyrosine phosphorylation of a latent cytoplasmic transcription factor
.
Science
257
,
809
-
813
.
Shah
,
P. S.
,
Link
,
N.
,
Jang
,
G. M.
,
Sharp
,
P. P.
,
Zhu
,
T.
,
Swaney
,
D. L.
,
Johnson
,
J. R.
,
Von Dollen
,
J.
,
Ramage
,
H. R.
,
Satkamp
,
L.
, et al. 
(
2018
).
Comparative flavivirus-host protein interaction mapping reveals mechanisms of dengue and Zika virus pathogenesis
.
Cell
175
,
1931
-
1945.e18
.
Shi
,
Y.
and
Gao
,
G. F.
(
2017
).
Structural biology of the Zika virus
.
Trends Biochem. Sci.
42
,
443
-
456
.
Shirasu-Hiza
,
M. M.
and
Schneider
,
D. S.
(
2007
).
Confronting physiology: how do infected flies die?
Cell. Microbiol.
9
,
2775
-
2783
.
Shuai
,
K.
,
Schindler
,
C.
,
Prezioso
,
V. R.
and
Darnell
,
J. E.
Jr.
(
1992
).
Activation of transcription by IFN-gamma: tyrosine phosphorylation of a 91-kD DNA binding protein
.
Science
258
,
1808
-
1812
.
Shuai
,
K.
,
Horvath
,
C. M.
,
Huang
,
L. H. T.
,
Qureshi
,
S. A.
,
Cowburn
,
D.
and
Darnell
,
J. E.
Jr.
(
1994
).
Interferon activation of the transcription factor Stat91 involves dimerization through SH2-phosphotyrosyl peptide interactions
.
Cell
76
,
821
-
828
.
Silies
,
M.
,
Yuva
,
Y.
,
Engelen
,
D.
,
Aho
,
A.
,
Stork
,
T.
and
Klambt
,
C.
(
2007
).
Glial cell migration in the eye disc
.
J. Neurosci.
27
,
13130
-
13139
.
Simpson
,
D. I. H.
(
1964
).
Zika virus infection in man
.
Trans. R. Soc. Trop. Med. Hyg.
58
,
335
-
338
.
Souza-Neto
,
J. A.
,
Sim
,
S.
and
Dimopoulos
,
G.
(
2009
).
An evolutionary conserved function of the JAK-STAT pathway in anti-dengue defense
.
Proc. Natl. Acad. Sci. USA
106
,
17841
-
17846
.
Steinberg
,
R.
,
Shemer-Avni
,
Y.
,
Adler
,
N.
and
Neuman-Silberberg
,
S.
(
2008
).
Human cytomegalovirus immediate-early-gene expression disrupts embryogenesis in transgenic Drosophila
.
Transgenic Res.
17
,
105
-
119
.
Steller
,
H.
,
Fischbach
,
K.-F.
and
Rubin
,
G. M.
(
1987
).
Disconnected: a locus required for neuronal pathway formation in the visual system of Drosophila
.
Cell
50
,
1139
-
1153
.
Sun
,
J.
,
Wu
,
D.
,
Zhong
,
H.
,
Guan
,
D.
,
Zhang
,
H.
,
Tan
,
Q.
and
Ke
,
C
. (
2016
).
Presence of Zika virus in conjunctival fluid
.
JAMA Ophthalmol.
134
,
1330
-
1332
.
Sun
,
X.
,
Hua
,
S.
,
Chen
,
H.-R.
,
Ouyang
,
Z.
,
Einkauf
,
K.
,
Tse
,
S.
,
Ard
,
K.
,
Ciaranello
,
A.
,
Yawetz
,
S.
,
Sax
,
P.
, et al. 
(
2017
).
Transcriptional changes during naturally acquired Zika virus infection render dendritic cells highly conducive to viral replication
.
Cell Rep.
21
,
3471
-
3482
.
Thomson
,
T. C.
,
Schneemann
,
A.
and
Johnson
,
J.
(
2012
).
Oocyte destruction is activated during viral infection
.
Genesis
50
,
453
-
465
.
Tsai
,
Y.-C.
and
Sun
,
Y. H.
(
2004
).
Long-range effect of upd, a ligand for Jak/STAT pathway, on cell cycle in Drosophila eye development
.
Genesis
39
,
141
-
153
.
Velazquez
,
L.
,
Fellous
,
M.
,
Stark
,
G. R.
and
Pellegrini
,
S.
(
1992
).
A protein tyrosine kinase in the interferon αβ signaling pathway
.
Cell
70
,
313
-
322
.
Wang
,
L.
,
Wang
,
S.
and
Li
,
W.
(
2012
).
RSeQC: quality control of RNA-seq experiments
.
Bioinformatics
28
,
2184
-
2185
.
Wang
,
A.
,
Thurmond
,
S.
,
Islas
,
L.
,
Hui
,
K.
and
Hai
,
R.
(
2017
).
Zika virus genome biology and molecular pathogenesis
.
Emerg. Microbes. Infect.
6
,
e13
.
Wang
,
J.
,
Liu
,
J.
,
Zhou
,
R.
,
Ding
,
X.
,
Zhang
,
Q.
,
Zhang
,
C.
and
Li
,
L.
(
2018
).
Zika virus infected primary microglia impairs NPCs proliferation and differentiation
.
Biochem. Biophys. Res. Commun.
497
,
619
-
625
.
Watling
,
D.
,
Guschin
,
D.
,
Müller
,
M.
,
Silvennoinen
,
O.
,
Witthuhn
,
B. A.
,
Quelle
,
F. W.
,
Rogers
,
N. C.
,
Schindler
,
C.
,
Stark
,
G. R.
,
Ihle
,
J. N.
, et al. 
(
1993
).
Complementation by the protein tyrosine kinase JAK2 of a mutant cell line defective in the interferon-gamma signal transduction pathway
.
Nature
366
,
166
-
170
.
Wong
,
S. L. A.
,
Chen
,
Y.
,
Chan
,
C. M.
,
Chan
,
C. S. M.
,
Chan
,
P. K. S.
,
Chui
,
Y. L.
,
Fung
,
K. P.
,
Waye
,
M. M. Y.
,
Tsui
,
S. K. W.
and
Chan
,
H. Y. E.
(
2005
).
In vivo functional characterization of the SARS-Coronavirus 3a protein in Drosophila
.
Biochem. Biophys. Res. Commun.
337
,
720
-
729
.
Xu
,
J.
and
Cherry
,
S.
(
2014
).
Viruses and antiviral immunity in Drosophila
.
Dev. Comp. Immunol.
42
,
67
-
84
.
Yoon
,
K. J.
,
Song
,
G.
,
Qian
,
X.
,
Pan
,
J.
,
Xu
,
D.
,
Rho
,
H. S.
,
Kim
,
N. S.
,
Habela
,
C.
,
Zheng
,
L.
,
Jacob
,
F.
, et al. 
(
2017
).
Zika-virus-encoded NS2A disrupts mammalian cortical neurogenesis by degrading adherens junction proteins
.
Cell Stem Cell
21
,
349
-
358.e6
.
Zhao
,
R. Y.
and
Elder
,
R. T.
(
2005
).
Viral infections and cell cycle G2/M regulation
.
Cell Res.
15
,
143
-
149
.
Zhao
,
R. Y.
,
Li
,
G.
and
Bukrinsky
,
M. I.
(
2011
).
Vpr-host interactions during HIV-1 viral life cycle
.
J. Neuroimmune Pharmacol.
6
,
216
-
229
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information