Information Flow in the White Matter During a Motor Task: A Structural Connectivity Driven Approach Guillermo Gallardo1 , Demian Wassermann 1 , Rachid Deriche 1 , Maxime Descoteaux2 and Samuel Deslauriers-Gauthier 2 Introduction: Cognitive tasks emerge from the interaction of functionally specialized cortical regions [START_REF] Verhagen | Hierarchical organization of parietofrontal circuits during goal-directed action[END_REF]. These interactions are supported by information flow through white matter fiber bundles connecting distant cortical regions. Estimating the information flow through white matter fiber bundles would therefore provide valuable information into the necessary cortical interactions to realize a task.

In this work, we build a Bayesian network representing cortical regions and their connections using a structural connectivity driven parcellation [START_REF] Gallardo | Groupwise Structural Parcellation of the Cortex: A Sound Approach Based on Logistic Models[END_REF] derived from diffusion MRI (dMRI). We then introduce Magnetoencephalography (MEG) measurements as evidence into this network to infer the information flow between cortical regions [START_REF] Deslauriers-Gauthier | Fibre directionality and information flow through the white matter: Preliminary results on the fusion of diffusion MRI and EEG[END_REF]. We show, for the first time, results on the interaction between the precentral, postcentral and occipital regions during a hand-movement task.

Methods: We applied our technique to subject 109123 (Male, aged 32-35) of the Human Connectome Project (HCP). First, we performed probabilistic tractography using 4000 seeds for each vertex of a pial mesh (8k vertices) representing the subject's cortex. We then obtained a connectivity driven cortical parcellation of 110 regions by using the parceling technique described in [START_REF] Gallardo | Groupwise Structural Parcellation of the Cortex: A Sound Approach Based on Logistic Models[END_REF]. The resulting parcellation, which follows the assumption of shared physical connections, was used to identify the number of streamlines connecting region pairs and their average length. As described in Deslauriers-Gauthier et al. ( 2016), the identified connections were then used to build a Bayesian network which associates a probability to all possible combinations of connections and cortical region states. Specifically, each connection and cortical region at every time point can be either active or inactive. By injecting MEG data as evidence into our Bayesian network and maximizing its entropy, we are able to obtain the posterior probability that a connection is active at any given time point. Here, the MEG data consisted of averaged evoked responses of a motor task where the subject was asked to move either his left or right hand following a visual cue.

Results: Our results are illustrated in the information flow diagrams of figures 1 and 2. Each row of the graph represents a cortical region and the circles correspond to individual time points. Information flow in connections is illustrated by lines connecting cortical regions at different time points. In both cases, green indicates a high probability that a region or connection is active whereas a white indicates a low probability.

Figure 1 illustrates information flow during the left hand movement task. The first activation at 100ms shows information transfer between the left occipital region and the left precentral gyrus and correspond to known visual evoked potential delays. The second activation which last from 200ms to 300ms presents information flow between the right occipital region and the right postcentral gyrus. This activation is consistent with motor activation contralateral to the hand movement. Furthermore, activity in the precentral gyrus within this time window is consistent with the initiation of a movement occurring at 300ms. Figure 2 illustrates information flow during the right hand movement task. The activation occurring at 100ms is similar to the results obtained with the left hand movement. However, during the 200 to 300ms window the information transfer is centralized in the left precentral gyrus and the occipital region. Again, this activation is consistent with motor activation contralateral to the hand movement.

Discussion and Conclusion:

We proposed a new method that combines structural information obtained from dMRI and functional information obtained from MEG. Our preliminary results obtained on HCP data show that we are able to recover information flow in the white matter that is consistent with expected cortical activation during a hand movement task.

Figure 1 :

 1 Figure 1: Information flow during the left hand movement task and corresponding MEG measurements. The cortical parcellation is connectivity driven, however cortical regions are identified using the overlapping Desikan-Killiany atlas for clarity.

Figure 2 :

 2 Figure 2: Information flow during the right hand movement task and corresponding MEG measurements. The cortical parcellation is connectivity driven, however cortical regions are identified using the overlapping Desikan-Killiany atlas for clarity.

Athena, Inria Sophia Antipolis, Mediterranée France

Université de Sherbrooke, Sherbrooke, Canada

This work has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (ERC Advanced Grant agreement No 694665 : CoBCoM -Computational Brain Connectivity Mapping). Data were provided by the Human Connectome Project, WU-Minn Consortium (Principal Investigators: David Van Essen and Kamil Ugurbil; 1U54MH091657) funded by the 16 NIH Institutes and Centers that support the NIH Blueprint for Neuroscience Research; and by the McDonnell Center for Systems Neuroscience at Washington University.


Audrey Alignier
email: audrey.alignier@toulouse.inra.fr

Helmut Meiss
email: helmut.meiss@dijon.inra.fr

Sandrine Petit

Xavier Reboud

Postdispersal weed seed predation ranged between 19 to 84% per week following a species preference rank

OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible.

Seed predation was significantly higher than removals of plastic globules and seed losses under the wire mesh.

Seed predation differed among the tested species (Wilcoxon χ2 = 381; df = 7; p < 0.0001). ("species preference rank", Figure)

  Double control for accidental seed losses:1) Seed cards placed under a narrow wire mesh 2) Presentation of dummy seeds (plastic globules) rates (19-84%) suggest that: [1] Pullaro, et al. 2006. Effects of killed cover crop mulch on weeds, weeds seeds, and herbivores. Agriculture, Ecosystems and Environment, 115, 97-104. [2] Westerman et al. 2003. Relative importance of vertebrates and invertebrates in epigeaic weed seed predation in organic cereals fields. Agriculture, Ecosystems and Environment, 95, 417-425.

H 1 )

 1 Weed seeds may constitute an important food resource for animals in the agro-ecosystem.2) Seed predation may be an important factor shaping weed population dynamics and community composition of arable weeds.So, to make optimal use of seed predation for reducing weed infestation in arable farming, research should focus on the mechanisms causing variability in predation.24th German Conference on Weed Biology and Weed Control. 4-6 March 2008, Hohenheim, Germany.

  Seed characteristics like weight, color, lipid or protein content are not linearly correlated to the observed predation rate.Spatial and temporal variation could not directly be explained by climatic variables or field border distance.

	Postdispersal Postdispersal weed weed seed seed predation predation ranged ranged between between
	19 to 84% per 19 to 84% per week week following following a a species species preference preference rank rank
	Seed Seed cards cards
	Wire Wire mesh mesh
	Dummy Dummy seeds seeds


Mirko Cortese

Sarah Goellner

Eliana Gisela Acosta

Christopher John Neufeldt

Olga Oleksiuk

Marko Lampe

Uta Haselmann

Charlotta Funaya

Nicole Schieber

Paolo Ronchi

Laurent Chatel-Chaix

Alessia Ruggieri

Ralf Bartenschlager
email: ralf.bartenschlager@med.uni-heidelberg.de

Martin Schorb

Priit Pruunsild

Yannick Schwab

Ultrastructural Characterization of Zika Virus Replication Factories

à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Ultrastructural Characterization of Zika Virus Replication Factories

INTRODUCTION

The pandemic outbreak of Zika virus (ZIKV), a member of the Flavivirus genus within the Flaviviridae family, and the association of ZIKV infections with severe disease raised intense concerns. Like other flaviviruses, ZIKV is primarily transmitted by mosquitoes, but more recently, unsuspected transmission modes and symptoms unique to ZIKV have been reported. Of note, infection of pregnant women with ZIKV can lead to congenital transmission and eventually to severe microcephaly in new-borns [START_REF] Pierson | Zika virus: immunity and vaccine development[END_REF]. Moreover, ZIKV is sexually transmissible through undefined mechanisms, and other neurological manifestations such as Guillain-Barre ´syndrome have been observed in infected individuals. Thus, intensive efforts have been undertaken to control this emerging disease, but neither antiviral therapies nor a prophylactic vaccine for ZIKV are currently available [START_REF] Pierson | Zika virus: immunity and vaccine development[END_REF]. This unmet medical need is exacerbated by the spread of the Aedes mosquito vector beyond tropical countries.

ZIKV was identified almost 70 years ago in Uganda, but severe symptoms were never reported. However, the recent outbreaks of the contemporary Asian lineage in the Yap Islands and French Polynesia in 2007 and 2013, respectively, have raised major interest in ZIKV [START_REF] Cao-Lormeau | Zika virus, French Polynesia, South Pacific[END_REF][START_REF] Duffy | Zika virus outbreak on Yap Island, Federated States of Micronesia[END_REF]. This Asian lineage is genetically closely related to the ''Brazilian strain'' currently spreading in the Americas and causing congenital microcephaly. Retrospective studies demonstrated that the French Polynesian strain also caused microcephaly [START_REF] Cauchemez | Association between Zika virus and microcephaly in French Polynesia, 2013-15: a retrospective study[END_REF], while it is hypothesized that the African lineage does not. More recently, the use of developing brain organoid models showed that in contrast to the African strain MR766, the Brazilian strain (Asian lineage) caused apoptosis of neural progenitor cells (NPCs) [START_REF] Cugola | The Brazilian Zika virus strain causes birth defects in experimental models[END_REF]. However, the viral determinants governing ZIKV pathogenesis are still unknown.

Upon entry into the host cell and release of the positive-strand RNA genome into the cytoplasm, the genome is used for the synthesis of a polyprotein that is cleaved co-and post-translationally into 10 proteins [START_REF] Acosta | Revisiting dengue virus-host cell interaction: new insights into molecular and cellular virology[END_REF]. These are the nonstructural proteins NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5 that are responsible for the replication of the viral RNA and the structural proteins Capsid, prM, and Envelope that assemble, together with the viral RNA genome, into infectious virions [START_REF] Acosta | Revisiting dengue virus-host cell interaction: new insights into molecular and cellular virology[END_REF]. The genome of ZIKV is replicated via a negative-strand RNA intermediate.

Flavivirus-infected cells undergo a massive remodeling of the endoplasmic reticulum (ER) to form membranous replication factories (RFs) [START_REF] Paul | Architecture and biogenesis of plusstrand RNA virus replication factories[END_REF]. Notably, flaviviruses, including ZIKV [START_REF] Hamel | Biology of Zika virus infection in human skin cells[END_REF], induce invaginations of the ER giving rise to clusters of vesicles designated vesicle packets (VPs), which are the presumed sites of viral RNA replication. Using electron tomography, several studies reported the 3D architecture of VPs formed during Dengue virus (DENV), West Nile virus (WNV), and tick-borne encephalitis virus (TBEV) infections [START_REF] Paul | Architecture and biogenesis of plusstrand RNA virus replication factories[END_REF]. In most of these studies, the vesicles are linked to the cytoplasm via an $10-nm-wide opening. This pore might allow the exit of the positive-strand RNA genome into the cytoplasm to be used for polyprotein synthesis (RNA translation) or packaging into virus particles budding into the ER lumen at sites juxtaposed to the VPs [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF]. In addition, in the case of WNV, vesicles were reported to be interconnected by a pore [START_REF] Gillespie | The endoplasmic reticulum provides the membrane platform for biogenesis of the flavivirus replication complex[END_REF]. Two additional ER-derived structures can be detected in flavivirusinfected cells: (1) convoluted membranes (CMs) accumulating in DENV-, TBEV-, and WNV-infected cells (Romero-Brey and Bartenschlager, 2014); and (2) DENV particles accumulating in regular arrays within enlarged ER cisternae. These ''virus bags'' might constitute a pool of assembled particles accumulating at this site prior to transport out of the cell. This transport is linked to proteolytic cleavage of the envelope glycoprotein by the cellular enzyme furin and it is required to render virions infectious [START_REF] Acosta | Revisiting dengue virus-host cell interaction: new insights into molecular and cellular virology[END_REF].

Although it has been reported that ZIKV induces VPs and CMs in infected cells [START_REF] Hamel | Biology of Zika virus infection in human skin cells[END_REF][START_REF] Hanners | Western Zika virus in human fetal neural progenitors persists long term with partial cytopathic and limited immunogenic effects[END_REF][START_REF] Offerdahl | Cytoarchitecture of Zika virus infection in human neuroblastoma and Aedes albopictus cell lines[END_REF], the high-resolution 3D architecture of ZIKV RFs and their biogenesis are unknown. Importantly, these analyses were conducted in cancer-derived hepatic and neuronal cell systems, which may be considered physiologically distant to ZIKV target cells such as human NPCs (hNPCs).

In addition to the ER, other cellular organelles such as mitochondria, lipid droplets, and autophagosomes are morphologically remodeled and functionally perturbed by flaviviruses [START_REF] Chatel-Chaix | Dengue virus perturbs mitochondrial morphodynamics to dampen innate immune responses[END_REF][START_REF] Heaton | Dengue virus and autophagy[END_REF][START_REF] You | Flavivirus infection impairs peroxisome biogenesis and early antiviral signaling[END_REF]. Furthermore, several studies identified a massive alteration of the cytoskeleton by some flaviviruses. For instance, intermediate filament and microtubule networks are reorganized in cells infected with DENV serotype 2 and their pharmacological disruption affects several steps of the viral replication cycle [START_REF] Chen | Vimentin is required for dengue virus serotype 2 infection but microtubules are not necessary for this process[END_REF][START_REF] Teo | Cellular vimentin regulates construction of dengue virus replication complexes through interaction with NS4A protein[END_REF]. Interestingly, the intermediate filament component vimentin associates with DENV NS4A and is important for the perinuclear localization of replication complexes [START_REF] Teo | Cellular vimentin regulates construction of dengue virus replication complexes through interaction with NS4A protein[END_REF]. Moreover, electron tomography-based 3D reconstructions using TBEV-infected human neurons revealed a connection between microtubules and virioncontaining vesicles [START_REF] Bı ´ly | Electron tomography analysis of tick-borne encephalitis virus infection in human neurons[END_REF]. Disrupting microtubules severely impeded virus production, highlighting the importance of this cytoskeleton network for the TBEV life cycle [START_REF] Bı ´ly | Electron tomography analysis of tick-borne encephalitis virus infection in human neurons[END_REF]. Whether ZIKV exploits similar or different strategies remains unknown.

Here we explored how host cell functions are altered to establish the membranous ZIKV RF. We compared intracellular remodeling induced by both the historical MR766 African strain and the contemporary Asian H/PF/2013 strain. Using Huh7 human hepatic cells and hNPCs, a physiologically relevant target cell of ZIKV, we demonstrate that independent of the lineage, ZIKV induces a drastic reorganization of intermediate filament and microtubule networks and reorganizes the ER to form VPs, zippered ER, and CMs. Cytoskeleton-targeting drugs suppress ZIKV replication, arguing that microtubules and intermediate filaments are critical for robust virus amplification.

RESULTS

ZIKV Replication Depends on the Rearrangement of the Cytoskeletal Network

With the aim to study the reorganization of subcellular structures by ZIKV, we first determined the replication kinetics of two virus strains: the historical African strain MR766 and the contemporary Asian strain H/PF/2013. In addition to the brain and placenta, the liver is one of the organs targeted by ZIKV in vivo in both the dam and fetus [START_REF] Adams Waldorf | Fetal brain lesions after subcutaneous inoculation of Zika virus in a pregnant nonhuman primate[END_REF]. Therefore, we used Huh7 human hepatoma cells, because they are highly susceptible to infection with several flaviviruses. Using immunofluorescence-based visualization of double-stranded (ds) RNA, a well-established replication intermediate, virus replication was detected as early as 4 hr postinfection in a low percentage of cells, with the percentage increasing up to 100% of cells 48 hr post-infection (Figures 1A and1B). Interestingly, although the percentage of infected cells was similar for both strains, we retrieved up to 10-fold higher extracellular virus titers for the Asian H/PF/2013 strain (Figure 1C). This might be due to the stronger cytotoxic effects exerted by the African MR766 strain, as revealed by the higher number of dead cells at late time points after infection (data not shown).

To gain insight into the subcellular distribution of different organelles, we conducted immunofluorescence analyses using confocal microscopy of ZIKV-infected Huh7 cells (Figure 1D). Protein disulfide isomerase (PDI), GM130 protein, and TGN46 protein were used as markers of the ER, Golgi apparatus, and trans-Golgi network (TGN), respectively. Compared to mock cells, where the Golgi apparatus and the TGN showed the classic shape of stacked cisternae, these two organelles appeared fragmented in infected cells and localized toward the interior of the dsRNA-positive area. Additionally, the nuclear morphology was altered to form a kidney-shaped structure, with dense ER accumulating in the perinuclear region that was also enriched for dsRNA.

In addition to changes in the morphology of these organelles, we also observed alterations of the cytoskeleton network within infected cells. In mock-infected cells, both intermediate filaments and microtubules were tightly linked and were observed to completely enwrap the nucleus (Figures 2A and2B). However, in infected cells, the cytoskeleton network was displaced by a perinuclear inclusion enriched for viral dsRNA and viral proteins (as revealed by staining for ZIKV NS3), both of which accumulated on the concave surface of the kidney-shaped nucleus (Figures 2A and 2B). Thicker cytoskeletal filaments surrounded the virus-induced structure (Figures 2A and2B). Confocal microscopy and 3D reconstruction indicated the formation of a cage-like arrangement of cytoskeleton elements clustering in close proximity of viral RNA (Figure S1A). Of note, a-tubulin staining revealed a bright spot, corresponding to the microtubules organizing center (MTOC), frequently located in the center of the viral inclusion (Figure 2B).

To monitor the dynamics of this cytoskeleton rearrangement, we performed live-cell imaging analysis of ZIKV-infected cells. A Huh7-derived cell line stably expressing tubulin-GFP was generated and infected with ZIKV. In mock cells, no apparent changes of the cytoskeletal network were observed within the 2-day observation period (Figure S1B). In contrast, in cells infected with ZIKV MR766, microtubule displacement and cage formation was first observed at $20 hr post-infection and became more prominent as infection proceeded (Figure 2C; Movie S1). At later time points after infection, vacuolization and microtubule disassembly were noticed, right before cell death as a result of virus-induced cytopathic effects (Movie S1).

To corroborate these findings in a noncarcinoma cell culture system that is physiologically relevant for ZIKV-associated neurological symptoms, we investigated ZIKV-induced subcellular alterations in hNPCs previously shown to be ZIKV target cells both in vitro and in vivo (Li et al., 2016a(Li et al., , 2016b;;[START_REF] Tang | Zika virus infects human cortical neural progenitors and attenuates their growth[END_REF]. Indeed, these cells were well permissive for both the MR766 and H/PF/2013 strains, with around 60% of cells becoming infected within 36 hr after infection (Figures S2A andS2B). Both strains replicated to comparable levels, as revealed by quantification of intracellular ZIKV RNA (Figure S2C), and virus titers were on the order of 10 6 PFU/mL (Figure S2D). Importantly, ZIKV also induced profound alterations of intermediate filaments and microtubules in this cell system, consistent with phenotypes detected in Huh7 cells (Figures S3A and S3B). Most striking was the effect on the intermediate filaments, as revealed by staining for nestin, an intermediate filament marker characteristic for neural stem and progenitor cells. As shown in Figure S3A, nestin filaments collapsed toward the perinuclear region, surrounding a dsRNA-positive structure presumably corresponding to the viral RFs (see below). Moreover, when visualized by stimulated emission depletion (STED) microscopy, which we used to discriminate between filamentous and non-filamentous structures, nestin filaments accumulated toward the periphery of the virus-induced structure, whereas the diffuse nestin signal present within the NS3-positive region suggested filament depolymerization within this particular area (Figure S3C).

Microtubule Stability Is Required for Efficient ZIKV Replication

Next we sought to define whether ZIKV-induced cytoskeleton remodeling was required for efficient replication. The observed microtubule displacement from the perinuclear region and the microtubule disassembly at the late time post-infection prompted us to evaluate the role of microtubules for ZIKV replication.

To avoid any potential effect on viral entry, we first identified the endpoint of entry by using bafilomycin A1, an inhibitor of endosome acidification that prevents the fusion step between the viral envelope and the endosomal membrane [START_REF] Yoshimori | Bafilomycin A1, a specific inhibitor of vacuolar-type H(+)-ATPase, inhibits acidification and protein degradation in lysosomes of cultured cells[END_REF]. Bafilomycin sensitivity was lost 3 hr post-infection, indicating that ZIKV entry and fusion were completed within this time span (Figure 2D). To examine the function of microtubules in ZIKV infection, we treated ZIKV-infected Huh7 cells with noncytotoxic concentrations of paclitaxel (Figure 2E), a microtubule-stabilizing drug. Immunofluorescence microscopy staining of a-tubulin demonstrated that the drug performed as expected, causing the formation of thick bundles of microtubules (Figure 2F). The compound was then provided 3 hr after infection with ZIKV strains, and virus titers were analyzed 24 hr post-infection (Figure 2G). Strikingly, paclitaxel treatment strongly reduced viral titers, suggesting that microtubule dynamics is required for efficient ZIKV infection. Of note, the same phenotype was observed in hNPCs, supporting the conclusion that cytoskeleton flexibility is also required for the ZIKV replication cycle in human neural progenitor cells (Figures S3D andS3E).

Ultrastructural Analysis of Membrane Rearrangements Induced in ZIKV-Infected Cells

Considering the profound ZIKV-induced rearrangements of the cytoskeletal network, we characterized the ultrastructural details of these alterations in infected Huh7 cells. To this end, we applied focus ion beam (FIB)-scanning electron microscopy to obtain a 3D overview of the cytoplasmic organization of the infected cell. To identify infected cells for this analysis, we used a reporter cell line expressing the fluorescently labeled ER marker Sec61b [START_REF] Chatel-Chaix | Dengue virus perturbs mitochondrial morphodynamics to dampen innate immune responses[END_REF]. Given the close similarity of DENV and ZIKV, we assumed that this protein accumulates in CMs, analogous to what we reported earlier for DENV-infected cells [START_REF] Chatel-Chaix | Dengue virus perturbs mitochondrial morphodynamics to dampen innate immune responses[END_REF]. Huh7 cells were infected with the MR766 strain and 24 hr later, cells containing Sec61b clusters were selected by correlative light and electron microscopy (EM) and processed for subsequent FIB-scanning electron microscopy analysis (Figure 3A). Following serial imaging and 3D reconstruction, the cytoplasm of the infected cells was found to display large areas covered by CMs with VPs located at the periphery of these structures and also accumulating in the perinuclear region (Figure 3A; Movie S2). Of note, the 3D reconstruction revealed cytoskeleton filaments surrounding most of the virus-induced VPs, while the large CM aggregates were located outside this region. This distribution nicely corroborated our immunofluorescence data described above and further implies that the perinuclear region containing both dsRNA and NS3 corresponds to the area where the majority of VPs accumulate.

Next, we focused our attention on the structure of ZIKV RFs and performed transmission electron microscopy (TEM) to achieve the resolution needed to analyze the ultrastructural details. Similar to our earlier observation with DENV, ZIKV-induced vesicles were found within the lumen of the rough ER. The average vesicle diameter was 80.82 ± 0.96 nm (n = 243) for the MR766 strain and 88.31 ± 1.25 nm (n = 219) for the H/PF/2013 strain (Figures 3B, 3D, andS4). Virus particles were found within the ER lumen or in isolated vesicles and had an outer diameter of around 40 nm (Figure 3D). Of note, CMs appeared to be crossed by microtubules (Figures 3C and3F). This organization was not observed in DENV-infected cells [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF]. In addition to these structures, we often observed tightly juxtaposed collapsed ER cisternae (Figures 3B and3E) with a restricted luminal space that was reminiscent of the zippered ER found in cells infected with infectious bronchitis virus (IBV) [START_REF] Maier | Infectious bronchitis virus generates spherules from zippered endoplasmic reticulum membranes[END_REF].

Time-course analyses revealed that virus-induced vesicles formed as early as 4 hr and 8 hr post-infection for the MR766 and H/PF/2013 strains, respectively (Figure 3G). Vesicles accumulated through the course of infection, reaching a peak at 36 hr post-infection (Figure 3G). A sharp increase in vesicle numbers was observed between 12 hr and 24 hr post-infection, correlating with the accumulation of viral RNA levels in the infected cells (Figure 3H). This result confirms previous findings regarding the close relationship between viral replication and VP formation [START_REF] Junjhon | Ultrastructural characterization and three-dimensional architecture of replication sites in dengue virus-infected mosquito cells[END_REF].

3D Architecture of ZIKV RFs

To determine the 3D organization of ZIKV RFs, we performed electron tomography on 250-nm-thick sections of Huh7 cells infected with the ZIKV H/PF/2013 strain (Figure 4). Consistent with the membranous alterations induced by the closely related DENV [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF] and other flaviviruses such as TBEV and WNV [START_REF] Gillespie | The endoplasmic reticulum provides the membrane platform for biogenesis of the flavivirus replication complex[END_REF][START_REF] Miorin | Three-dimensional architecture of tick-borne encephalitis virus replication sites and trafficking of the replicated RNA[END_REF], ZIKV induced an organelle-like compartment consisting of arrays of invaginated vesicles within dilated ER sheets; their decoration with ribosomes found in several instances (e.g., Figure 3B) argues for rough ER. Bundles of presumably intermediate filaments were observed in close proximity of the VPs (Figures 4A and4B; Movie S3). As previously reported for other flavirivuses (Romero-Brey and Bartenschlager, 2014), the lumen of the invaginated vesicles was connected to the cytosol via a narrow pore-like structure. The calculated pore width of 10.9 ± 1.6 nm (n = 8) was comparable to the pore size observed in DENV-infected cells (11.2 ± 1.6 nm) [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF]. Similar to DENV, ER sheets containing virus particles were often found adjacent to virus-induced invaginated vesicles on the opposite side of the pore (Figures 4D and4E), suggesting that ZIKV particles bud in close proximity to the genome replication site. Although the pore-like structure was readily visible, a connection between adjacent vesicles was not observed (Figures 4D and4E, yellow arrows). These results are similar to those previously described for DENV [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF] and TBEV [START_REF] Miorin | Three-dimensional architecture of tick-borne encephalitis virus replication sites and trafficking of the replicated RNA[END_REF] but differ from what has been reported for WNV [START_REF] Gillespie | The endoplasmic reticulum provides the membrane platform for biogenesis of the flavivirus replication complex[END_REF] or Langat virus [START_REF] Offerdahl | A three-dimensional comparison of tick-borne flavivirus infection in mammalian and tick cell lines[END_REF]. 3D reconstruction of regions encompassing the zippered ER revealed that the collapsed ER was contiguous with regions containing invaginated vesicles (Figures 4F-4H; Movie S4), suggesting a transition from a swollen ER area filled with virus-induced vesicles to a luminally restricted area.

Next, we sought to analyze the ultrastructural alterations induced by ZIKV infection in hNPCs. To this end, hNPCs were infected with either ZIKV strains; 24 hr later, cells were fixed and analyzed by TEM. Inspection of mock cells revealed stacked rough ER cisternae forming an organized network surrounding the nucleus, with mitochondria dispersed throughout the cytoplasm and often in contact with the ER (Figures 5A and5B). However, in infected cells, the rough ER appeared as a convoluted tangle of cisternae creating a large perinuclear inclusion residing on the concave side of the kidney-shaped nucleus (Figures 5C and5E). Mitochondria accumulated on the periphery of the convoluted rough (r)ER region (Figures 5C, 5E, and5F). Virus-induced vesicles were often observed to cluster within regions of more densely packed rER (Figures 5D, 5F, and5G). Such areas were more frequent and prominent in MR766-infected cells (Figures 5E-5H). Cytoskeletal filaments were observed on the edge of the ER region (Figures 5D and5G). The average diameters of ZIKV-induced vesicles was 62.80 ± 0.49 nm (n = 516) and 65.51 ± 0.56 nm (n = 457) for the African and the Asian strains, respectively. Compared to the vesicle diameters observed in Huh7 cells, ZIKV-induced vesicles in infected hNPCs appeared smaller (Figure S4), arguing that vesicle size might be cell-type dependent. Another difference among the two cell lines was the complete absence of convoluted smooth ER membranes and zippered rough ER in infected hNPCs.

With the aim of revealing the 3D architecture of the viral RFs, we employed electron tomography analysis of hNPCs infected for 24 hr with the ZIKV H/PF/2013 strain (Figure 6; Movie S5). Virus-induced vesicles appeared as tightly stacked invaginations within the rough ER with narrow pores connecting the luminal side of the vesicles with the cytosol (Figures 6C-6H). Comparable to Huh7 cells, no connections were observed among adjacent vesicles in hNPCs. Quantification of pore sizes revealed an average width of 11.4 ± 1.2 nm and 10.7 ± 1.9 nm in H/PF/2013and MR766-infected hNPCs, respectively. Interestingly, although vesicle diameters were different, pores sizes were comparable not only between hNPCs and Huh7 cells but also between the two virus strains, arguing that the vesicle pore is formed by viral proteins only or by a conserved cellular machinery. Consistent with previous observations, virions were found in close proximity to the invaginated vesicles, usually within tightly opposed ER sheets (Figures 6B and6F). Diffuse electron density was observed on the cytosolic side of the vesicle pore, together with frequent virion budding events occurring within the opposed ER sheet (Figure 6F). 3D reconstruction of virion accumulations revealed that ZIKV particles stacked in regular arrays inside enlarged ER cisternae, which are contiguous and connected to VPs within the same ER compartment. In addition, an intricate network of intermediate filaments was present in close proximity of the ER membranes, although single filaments rather than bundles were observed.

In conclusion, these results reveal commonalities and differences of the 3D architecture of presumed ZIKV RFs in a human hepatic cell and hNPCs. While the overall architecture is well comparable to the one we reported earlier for DENV [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF], the utilization and perturbation of the cytoskeletal network by ZIKV sets it apart from DENV.

DISCUSSION

In this study, we report the high-resolution 3D architecture of ZIKV RFs in Huh7 human hepatic cells (a standard cell system for flavi- viruses) and, more importantly, in hNPCs, one of the physiological targets of ZIKV. While [START_REF] Offerdahl | Cytoarchitecture of Zika virus infection in human neuroblastoma and Aedes albopictus cell lines[END_REF] recently reported low-resolution electron tomography of ZIKV-induced vesicles in infected neuroblastoma cells, our study unveils detailed features of ZIKV RFs in hNPCs. Similar to DENV, WNV, and TBEV [START_REF] Bı ´ly | Electron tomography analysis of tick-borne encephalitis virus infection in human neurons[END_REF][START_REF] Gillespie | The endoplasmic reticulum provides the membrane platform for biogenesis of the flavivirus replication complex[END_REF][START_REF] Junjhon | Ultrastructural characterization and three-dimensional architecture of replication sites in dengue virus-infected mosquito cells[END_REF][START_REF] Miorin | Three-dimensional architecture of tick-borne encephalitis virus replication sites and trafficking of the replicated RNA[END_REF][START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF], ZIKV-induced VPs are also located within the ER lumen, implying that the vesicles originate from this organelle. Hence, all flaviviruses appear to share similar hallmarks of the architecture of their respective RFs. Since DENV VPs were shown to be associated with dsRNA [START_REF] Junjhon | Ultrastructural characterization and three-dimensional architecture of replication sites in dengue virus-infected mosquito cells[END_REF][START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF], these structures are believed to be the site of RNA replication. This spatial segregation would allow concentrating metabolites required for RNA replication while protecting viral RNA from cellular nucleases and innate immunity-triggering cytosolic RNA sensors. Furthermore, we observed a pore-like opening toward the cytosol allowing the exit of neo-synthesized positive-strand viral RNA into the cytoplasm, where it could be used for protein synthesis or virion formation. Very interestingly, electron-dense structures resembling virions could be detected in ER regions that were juxtaposed to the pore of VPs (Figures 6F and6H). This suggests that particle assembly (including viral RNA encapsidation) and budding are spatially coordinated with RNA replication, reminiscent of what we reported for DENV [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF].

No obvious differences between the MR766 and H/PF/2013 strains were detected in terms of VP ultrastructural features, implying that the African and the Asian lineages use similar strategies to induce their RFs. The mean diameter of ZIKVinduced vesicles in Huh7 cells is 81 nm and 88 nm for MR766 and H/PF/2013, respectively, comparable with the reported size of DENV2 VPs (80-90 nm) [START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF]. Interestingly, the vesicle diameter was lower in NPCs, independent of the ZIKV lineage (63-66 nm), even though these cells were as permissive as Huh7 cells in terms of virus yield (Figures 1C and S2D). This result suggests that cell-specific determinants influence the morphogenesis of ZIKV RFs. Whether this is due to cell-specific differences in the stoichiometry of viral proteins, RNA, and host factors within ZIKV VPs remains to be determined.

In Huh7 cells, both MR766 and H/PF/2013 induced CMs comparable to DENV, WNV, and TBEV [START_REF] Gillespie | The endoplasmic reticulum provides the membrane platform for biogenesis of the flavivirus replication complex[END_REF][START_REF] Welsch | Composition and three-dimensional architecture of the dengue virus replication and assembly sites[END_REF]. Interestingly, arrays of microtubules passing through CMs could be detected, a phenotype that was never observed in DENV-infected Huh7 cells (M.C. and R.B., unpublished data). Using live-cell imaging, we recently showed that DENV CMs are highly dynamic structures, which also undergo merging and fission events [START_REF] Chatel-Chaix | Dengue virus perturbs mitochondrial morphodynamics to dampen innate immune responses[END_REF]. Hence, in the case of ZIKV infection, CM-associated microtubules might influence mobility and morphology of this dynamic membranous structure.

The zippered ER observed in ZIKV-infected cells was reminiscent of structures observed during infection with IBV, an avian gamma-coronavirus generating ER-derived double-membrane vesicles and spherules [START_REF] Maier | Infectious bronchitis virus generates spherules from zippered endoplasmic reticulum membranes[END_REF]. As for IBV, the precise role of zippered ER for ZIKV replication remains unknown, although they appear to be unique to ZIKV in the genus Flavivirus. Of note, CMs and zippered ER were not found in NPCs. This was probably not due to virus replication, since virus yields were comparable at the time point of TEM analysis (24 hr post-infection). This discrepancy might be due to cell-type-specific differences in the availability of host factors required for CM and zippered ER biogenesis. In any case, the physiological relevance of these structures for ZIKV-associated neurogenesis defects is questionable, although important roles for other aspects of viral pathogenesis cannot be excluded.

Taking advantage of various imaging techniques, we demonstrate that ZIKV induces major rearrangements of the cytoskeleton. Indeed, both intermediate filaments and microtubule networks are reorganized to form a cage-like structure in which dsRNA, NS3, and remodeled ER accumulate. Microtubule dynamics play a critical role in ZIKV replication, as deduced from the strong antiviral effect of the microtubule-stabilizing drug paclitaxel (Figures 2G andS3E). Interestingly, DENV, WNV, JEV, and TBEV also induce cytoskeleton remodeling [START_REF] Chen | Vimentin is required for dengue virus serotype 2 infection but microtubules are not necessary for this process[END_REF][START_REF] Chiou | Association of Japanese encephalitis virus NS3 protein with microtubules and tumour susceptibility gene 101 (TSG101) protein[END_REF][START_REF] Ng | Flavivirus infection: essential ultrastructural changes and association of Kunjin virus NS3 protein with microtubules[END_REF][START_REF] Vancova | Morphological changes in human neural cells following tick-borne encephalitis virus infection[END_REF][START_REF] Teo | Cellular vimentin regulates construction of dengue virus replication complexes through interaction with NS4A protein[END_REF]. In fact, treatment with a microtubule-destabilizing drug enhances DENV production [START_REF] Chen | Vimentin is required for dengue virus serotype 2 infection but microtubules are not necessary for this process[END_REF], suggesting that DENV and ZIKV have developed similar cytoskeleton remodeling strategies. Yeast two-hybridbased interaction screenings identified the intermediate filament vimentin as an interactant of flaviviral NS5 (Le [START_REF] Breton | Flavivirus NS3 and NS5 proteins interaction network: a highthroughput yeast two-hybrid screen[END_REF]. Moreover, vimentin was reported to interact with the N terminus of DENV NS4A and to colocalize with replication complexes [START_REF] Teo | Cellular vimentin regulates construction of dengue virus replication complexes through interaction with NS4A protein[END_REF]. In contrast, we found that intermediate filaments surround ZIKV RFs in the perinuclear (MTOCpositive) area, arguing for differences between ZIKV and DENV cytoskeleton remodeling.

The cytoskeleton cage observed here for ZIKV might contribute to concentrate the different ZIKV-induced membranous structures and to spatially organize replication and assembly complexes in a closed environment to which RIG-like receptors have limited access. Interestingly, such a model was recently proposed for the related hepatitis C virus (HCV). Indeed, the HCV RF excludes retinoic acid-inducible gene I (RIG-I) by the recruitment of specific nuclear pore proteins [START_REF] Neufeldt | The hepatitis C virus-induced membranous web and associated nuclear transport machinery limit access of pattern recognition receptors to viral replication sites[END_REF]. Considering the accepted role of the cytoskeleton in innate immunity [START_REF] Mostowy | The cytoskeleton in cell-autonomous immunity: structural determinants of host defence[END_REF], it is tempting to speculate that ZIKV-mediated remodeling of the cytoskeleton might interfere with antiviral responses to favor viral replication. This might be either by direct contribution of the cytoskeleton or indirectly by promoting RF formation. In support to the latter assumption, evidence has been presented that both TBEV and the flavivirus Kunjin virus shield viral dsRNA in ER membrane invaginations from sensing by RIG-like receptors [START_REF] Hoenen | West Nile virus-induced cytoplasmic membrane structures provide partial protection against the interferon-induced antiviral MxA protein[END_REF][START_REF] Overby | Tick-borne encephalitis virus delays interferon induction and hides its double-stranded RNA in intracellular membrane vesicles[END_REF]. Further work (e.g., by using pharmacological inhibitors) will be needed to unravel the possible role of the cytoskeleton in the biogenesis of the ZIKV RF and how ZIKV-induced cytoskeleton remodeling might interfere with the innate immune response.

Cytoskeleton dynamics is a crucial factor of neurogenesis. For instance, vimentin knockout mice have an altered morphology of Purkinje cells in the cerebellar cortex [START_REF] Colucci-Guyon | Cerebellar defect and impaired motor coordination in mice lacking vimentin[END_REF]. In addition, neurofilaments, a sub-class of intermediate filaments, determine the diameter of myelinated axons [START_REF] Elder | Agerelated atrophy of motor axons in mice deficient in the mid-sized neurofilament subunit[END_REF] and a mutation in b3-tubulin is associated with defects of commissural axons and cranial nerves [START_REF] Tischfield | Human TUBB3 mutations perturb microtubule dynamics, kinesin interactions, and axon guidance[END_REF]. Of note, intermediate filaments also possess cytoprotective functions such as apoptosis inhibition, organelle anchoring, and regulation of cell growth as well as motility [START_REF] Snider | Post-translational modifications of intermediate filament proteins: mechanisms and functions[END_REF][START_REF] Toivola | Intermediate filaments take the heat as stress proteins[END_REF]. It is therefore tempting to speculate that the remodeling of the cytoskeleton by ZIKV might contribute to the severe fetal neurogenesis defects following congenital viral transmission. Deciphering the molecular mechanisms regulating the biogenesis of ZIKV RFs might help to develop novel therapies suppressing virus replication and reducing the risk of congenital transmission and microcephaly.

EXPERIMENTAL PROCEDURES

Further details of experimental procedures used are given in the Supplemental Experimental Procedures.

Cell Lines Huh7, 293T, and VeroE6 cells were maintained in DMEM supplemented with 10% fetal bovine serum (FBS; Sigma-Aldrich), 100 U/mL penicillin, 100 mg/mL streptomycin, 2 mM L-glutamine, and 1% non-essential amino acids (all from Gibco, Life Technologies). C6/36 mosquito cells were grown in Leibovitz L-15 medium (Gibco, Life Technologies) supplemented with 10% FBS, 100 U/mL penicillin, 100 mg/mL streptomycin, 2 mM L-glutamine, and 10 mM HEPES (Gibco, Life Technologies). Huh7 yellow fluorescent protein (YFP)-Sec61b cells [START_REF] Chatel-Chaix | Dengue virus perturbs mitochondrial morphodynamics to dampen innate immune responses[END_REF] were additionally supplemented with 1 mg/mL puromycin. Generation and maintenance of hNPCs are specified in the Supplemental Experimental Procedures.

Immunofluorescence Microscopy

For immunofluorescence microscopy, cells grown on glass coverslips were washed twice in cold PBS before fixation in 100% ice-cold methanol for 10 min at À20 C. Methanol was removed and samples were air-dried at room temperature. Fixed cells were incubated for 1 hr in blocking buffer (5% BSA in PBS). Samples were incubated with primary antibodies diluted in blocking buffer for 2 hr. After washing thrice with PBS, samples were incubated for 1 hr with Alexa Fluor 488-, 568-, or 647-conjugated secondary antibodies (Molecular Probes, Life Technologies) in the dark and subsequently washed three times with PBS for 10 min each. Nuclear DNA was stained with 4,6-diamidino-2-phenylindole, dihydrochloride (DAPI; Sigma-Aldrich) for 10 min, washed with PBS, and rinsed with ultrapure water. Coverslips were mounted on microscopy slides using Fluoromount G (SouthernBiotech). Fluorescence images were acquired with a Nikon Ti Eclipse microscope and analyzed with the NIS-Element AR software package. For 3D reconstruction, 130-nm optical sections were acquired with a Nikon TE2000-E confocal microscope equipped with an Ultraview ERS spinning disc (PerkinElmer). Images were deconvoluted with Autoquant X3 software (Media Cybernetics) and 3D reconstruction was performed with Imaris 8 (Bitplane).

TEM and Tomography

Huh7 cells were fixed with 2.5% glutaraldehyde in cacodylate buffer (50 mM cacodylate, 50 mM KCl, 2.6 mM CaCl 2 , 2.6 mM MgCl 2 , and 2% sucrose) for at least 1 hr at room temperature. Samples were incubated with 2% osmium tetroxide in 50 mM cacodylate buffer for 40 min on ice, washed three times with EM-grade water and incubated for 30 min in 0.5% uranyl acetate in water. Samples were rinsed three times with water, dehydrated in a graded ethanol series (from 40% to 100%) at room temperature before being embedded in epoxy resin, and polymerized for at least 48 hr at 60 C. Ultrathin sections of 70 nm were obtained by sectioning with a Leica EM UC6 ultramicrotome (Leica Microsystems) and a diamond knife. Sections were collected on EM grids, counterstained using lead citrate and uranyl acetate, and examined with a Zeiss 10C transmission electron microscope (80 kV; Carl Zeiss).

For electron tomography, 250-nm-thick sections were prepared. Counterstained and gold-labeled grids were placed in a high-tilt holder and digital images were recoded as single-axis tilt series over a À53 to 53 tilt range (increment, 1 ) on an FEI Tecnai TF20 (Figure 4) or dual-axis tilt series over a À60 to 60 tilt range on an FEI Tecnai TF30 microscope (Figure 6). Tomograms were reconstructed using IMOD software [START_REF] Kremer | Computer visualization of three-dimensional image data using IMOD[END_REF] (http://bio3d.colorado.edu/imod).

Figure 1 .

 1 Figure 1. ZIKV Growth Kinetic in Huh7 Cells (A) Time course of the spread of ZIKV infection (MOI = 5) in Huh7 cells as detected by immunofluorescence using a dsRNA-specific antibody (green). Nuclear DNA was stained with DAPI (blue). Scale bars, 20 mm. (B) Quantification of the percentage of infected cells shown in (A). Means ± SD from two independent experiments (four fields of view each) are shown. (C) Titers of ZIKV released from infected Huh7 cells. Supernatants from infected cells (MOI = 5) were collected at the indicated time points and titers were determined by plaque assay. Means ± SD from two independent experiments are shown. (D) Huh7 cells were infected with ZIKV (MOI = 5) for 36 hr. Cells were fixed and cellular proteins specified in the left of each panel (red) as well as dsRNA (green) were detected by immunofluorescence microscopy. Colored images are merges of the corresponding individual panels on the left. hpi, hours post-infection; scale bars, 10 mm.

Figure 2 .

 2 Figure 2. Rearrangement of the Cytoskeleton in ZIKV-Infected Huh7 Cells (A and B) Huh7 cells were mock infected or infected for 36 hr with the MR766 or the H/PF/21013 ZIKV strain (MOI = 5) prior to methanol fixation. ZIKV infection was visualized by immunostaining for dsRNA and Cytokeratin 8 or NS3 and a-tubulin (red), respectively. Nuclear DNA was stained with DAPI (blue). Naı ¨ve cells (mock) served as the control. Representative view fields are shown. Yellow squares in the left panels indicate the magnified area shown in the corresponding right color panel. Scale bars, 10 mm. (C) Dynamics of tubulin in Huh7 GFPtubulin-expressing cells infected with ZIKV MR766 (MOI = 10). Cropped sections of the GFP-tubulin channel recorded 3, 10, 20, 30, and 40 hr postinfection (see Movie S1). A representative infected cell (a) is highlighted to facilitate tracking. Images were taken every 30 min. Scale bars, 20 mm. (D) Titers of ZIKV particles released from Huh7 cells pretreated 2 hr prior to infection or 3 hr post-infection with bafilomycin A1 (Baf; 2.5 nM). DMSO was used as the control. Supernatants were collected 24 hr post-infection and titers were determined using the plaque assay. Histogram bars represent the fold difference to DMSO-treated cells. Means ± SD from two replicates are shown. (E) Viability of Huh7 cells as determined by measurement of ATP level after a 24-hr treatment with 12.5 mM paclitaxel. Treatment with DMSO was used as the control. Means ± SD from three independent experiments are shown. (F) Huh7 cells were treated for 24 hr with 12.5 mM paclitaxel or with equal amounts of DMSO. a-tubulin was detected by immunofluorescence microscopy. Scale bars, 20 mm. (G) Titer of extracellular virus released from ZIKV-infected Huh7 cells treated with 12.5 mM paclitaxel for 20 hr. Cells were treated with the drug 3 hr after infection with ZIKV (MOI = 5). Histogram bars represent the fold difference to DMSO-treated cells. Means ± SD from three independent experiments are shown.

Figure 3 .

 3 Figure 3. Ultrastructural Analysis of ZIKV-Infected Cells (A) Huh7 YFP-Sec61b cells were infected with ZIKV MR766 (MOI = 5) for 24 hr and analyzed by FIB-scanning electron microscopy. Infected cells were identified by visual inspection for clusters of YFP-Sec61b corresponding to convoluted membranes (CMs; Chatel-Chaix et al., 2016) followed by correlative light and electron microscopy. Shown is a 3D-reconstruction of a representative infected cell: vesicle packets (Vp, yellow), cytoskeletal filaments (Fi, red), nucleus (NC, blue), convoluted membrane (green), and microtubuleorganizing center (MTOC, pink). Scale bar, 2 mm. (B-F) Transmission electron microscopy (TEM) images of 70-nm-thin sections of resin-embedded Huh7 cells infected for 24 hr with ZIKV MR766 (B and C) or H/PF/2013 (D-F) (MOI = 5). (B and D) Virus-induced vesicles (Ve) clustering within the rough ER lumen. (B and E) Tightly juxtaposed and collapsed ER cisternae (zippered ER [zER]). (C and F) Microtubules (MTs) are often associated with convoluted membranes. Red scale bars, 200 nm; yellow scale bars, 500 nm. (G) Time course of accumulation of ZIKV-induced vesicles. Huh7 cells were infected with MR766 or H/PF/2013 prior to fixation at given time points and processed for TEM; 70-nm-thick sections were acquired. For quantification of vesicle numbers per cell, 10 cell profiles were randomly taken for each time point and the average number of vesicles per profile was plotted. Means ± SEM are shown. (H) Time course of accumulation of intracellular ZIKV RNA in infected Huh7 cells (MOI = 5) measured by two-step qRT-PCR. Values for each sample were normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and are expressed relative to mock cells. Results from two independent experiments (three replicates each) are shown. Means ± SD are displayed.

Figure 4 .

 4 Figure 4. Electron Tomography of ZIKV-Induced Vesicle Packets in Huh7 Cells (A) Huh7 cells were fixed 24 hr post-infection, embedded in epoxy resin, and analyzed by electron tomography. A tomographic slice showing zippered ER, ZIKV virions (Vi), and virus-induced invaginated vesicles (Ve) within the ER is displayed. The yellow square represents the cropped section shown in (B) which is a 3D surface model of the ER lumen containing vesicles and virions. ER membranes are depicted in blue, virions in gold, and cytoskeletal filaments in red. (C) Side view of the 3D model showing the pore-like openings (red arrowheads) in the ER membrane. (D) Slice through the tomogram showing ER cisternae containing virions opposing ZIKV-induced vesicles. A close apposition of adjacent vesicles is visible (yellow arrowheads), but direct connections are not visible. Red arrowhead, vesicle pore. (E) Reconstruction of the area shown in (D). (F) Slice through a tomogram showing vesicle packets and details of the zippered ER (zER) (blue arrowheads). (G) Magnification of the yellow squared area shown in (F). (H) 3D reconstruction of the area in (F). Regions of zippered ER are indicated with blue arrowheads. Vesicles sharing the same ER lumen as a virion can be observed. Scale bars, 100 nm.

Figure 5 .

 5 Figure 5. Ultrastructural Analysis of ZIKV-Infected hNPCs (A-H) TEM images of 70-nm-thin sections of ZIKVinfected and resin-embedded hNPCs. Cells were mock infected (A and B) or infected (MOI = 5) with the H/PF/2013 (C and D) or MR766 (E-H) ZIKV strain for 24 hr. Yellow boxed areas are enlarged in the given panels. Overviews of ZIKV-infected cells (C and E) revealing incised nuclei and the ER inclusion. Mitochondria (mito) are located in the cell periphery adjacent to the outer rim of the ER (C and E). Virus-induced vesicles (Ve) and virions (Vi) are indicated in the magnified panels (D and H) and are marked with white and gold arrowheads, respectively. Red arrowheads, cytoskeletal filaments. Scale bars, 2 mm in (A), (C), and (E); 500 nm in (F) and (G); 200 nm in (B), (D), and (H).

Figure 6 .

 6 Figure 6. Electron Tomography of ZIKV-Induced Vesicle Packets in hNPCs hNPCs were fixed 24 hr post-infection, embedded in epoxy resin, and analyzed by electron tomography. (A and B) A slice through a tomogram is shown depicting ZIKV-induced vesicles (Ve) within the rough ER as well as virions (Vi). (C) 3D surface model of the boxed area in (A), showing virus-induced vesicles, virions, and intermediate filaments. ER membranes are depicted in light blue, vesicles in dark blue, virus particles in gold, and filaments in red. ER tubules connected the virions with the invaginated vesicles. (D-F) Slice through the tomogram showing the porelike openings (colored arrowheads) of ZIKV-induced vesicles toward the cytoplasm. A potential ZIKV budding event (Vi?) on the ER tubule opposing the vesicle pore can be observed in (F). (G and H) Reconstruction of the areas shown in (D) to (F). Colors of arrowheads refer to the vesicle pores marked in (D) and (E). Scale bars, 100 nm in (A), (B), and (D-H); 200 nm in (C).

Cell Reports 18, 2113-2123, February 28, 2017

Cell Reports18, 2113-2123, February 28, 2017 2115

Cell Reports18, 2113-2123, February 28, 2017 2117

Cell Reports18, 2113-2123, February 28, 2017 2119

ACKNOWLEDGMENTS

We gratefully acknowledge the expert support provided by the Electron Microscopy Core Facility at Heidelberg University, the Flow Cytometry and FACS Core Facility of the Center for Molecular Biology of Heidelberg University (ZMBH), and the Electron Microscopy Core Facility at the European Molecular Biology Laboratory (EMBL) Heidelberg. We thank the Advanced Light Microscopy Facility (ALMF) at the European Molecular Biology Laboratory (EMBL) and Leica Microsystems for support. We are grateful to Dr. Vibor Laketa and the Infectious Diseases Imaging Platform (IDIP) at the Department of Infectious Diseases, University Hospital Heidelberg, for the excellent microscopy support. We are grateful to Hilmar Bading, in whose laboratory the hNPCs were generated, to Dr. Gualtiero Alvisi for constructing the YFP-Sec61b-encoding plasmid, and to the European Virus Archive (EVAg) for provision of ZIKV strains. We are very grateful to Philipp Klein for support in cell sorting and live-cell imaging. We thank Thais Moraes for virus stock production. This work was supported by grants from the Deutsche Forschungsgemeinschaft (SFB 1129 to R.B and A.R). L.C.-C. holds a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada (RGPIN-2016-05584). P.P. holds a fellowship by the Excellence Cluster Cell Networks at Heidelberg University. C.J.N is funded by a European Molecular Biology Organization (EMBO) Long-Term Fellowship (ALTF 466-2016).

SUPPLEMENTAL INFORMATION

Supplemental Information includes Supplemental Experimental Procedures, four figures, and five movies and can be found with this article online at http://dx.doi.org/10.1016/j.celrep.2017.02.014.

AUTHOR CONTRIBUTIONS