Neuroanatomical abnormalities in a nonhuman primate model of congenital Zika virus infection
Abstract
We evaluated neuropathological consequences of fetal ZIKV exposure in rhesus monkeys, a translatable animal model for human neural development, by carrying out quantitative neuroanatomical analyses of the nearly full-term brains of fetuses infected with ZIKV and procedure-matched controls. For each animal, a complete cerebral hemisphere was evaluated using immunohistochemical (IHC) and neuroanatomical techniques to detect virus, identify affected cell types, and evaluate gross neuroanatomical abnormalities. IHC staining revealed the presence of ZIKV in the frontal lobe, which contained activated microglia and showed increased apoptosis of immature neurons. ZIKV-infected animals exhibited macrostructural changes within the visual pathway. Regional differences tracked with the developmental timing of the brain, suggesting inflammatory processes related to viral infiltration swept through the cortex, followed by a wave of cell death resulting in morphological changes. These findings may help explain why some infants born with normal sized heads during the ZIKV epidemic manifest developmental challenges as they age.
Editor's evaluation
This rigorous study provides compelling evidence that Zika virus infections in infants can markedly impact brain development through neuroinflammatory mechanisms. The work will have broad interest among developmental neurobiologists, as well as scientists whose work focuses on ZIKV pathogenesis.
https://doi.org/10.7554/eLife.64734.sa0Introduction
In late 2015, medical professionals in the northeast region of Brazil reported a surge in the number of children born with microcephaly, with an increase of 265% in the number of cases during 2015–2016 (Microcephaly Epidemic Research Group, 2016; Kleber de Oliveira et al., 2016). Infants born with microcephaly during the Zika outbreak also had a constellation of clinical presentations that ultimately became known as Congenital Zika Syndrome (CZS): skull deformation, abnormally small cerebral cortices, and retinal scarring, among others (Broussard et al., 2018; Moore et al., 2017; Musso et al., 2019). Clinicians have identified additional health issues as these children have continued to grow and develop (Wheeler et al., 2018; Cardoso et al., 2019), including visual impairment (Ventura et al., 2017), increased risk for autism (Vianna et al., 2018), atypical motor development, increased risk for cerebral palsy (Marques et al., 2019), and poor sleep quality (Pinato et al., 2018). Further screening has identified cases where children not diagnosed with CZS at birth exhibit CZS by their first birthday; up to 50% of children born to mothers infected with Zika virus during pregnancy infection exhibit anatomical or behavioral abnormalities (Cardoso et al., 2019; Wheeler, 2018). Taken together, these findings suggest that while microcephaly is an extreme manifestation of CZS, children born with normal sized heads may also have other significant neural deficits. Understanding the anatomy of those deficits is critical for predicting the developmental trajectories of children with CZS and ultimately for developing effective treatments and interventions, and the focus of the present report.
Since the onset of the ZIKV pandemic that began in 2015, a number of studies have revealed important information about how ZIKV targets specific cell types and proteins within the central nervous system (for reviews Khaiboullina et al., 2019; Sutarjono, 2019). The neurotropism of the virus may result from ZIKV using AXL, a tyrosine kinase receptor important in modulating the innate immune system (Rothlin et al., 2007; Hafizi and Dahlbäck, 2006), to enter the brain (Meertens et al., 2017; Rossi et al., 2020; Nowakowski et al., 2016). Single-cell RNA sequencing analyses have demonstrated that AXL is highly expressed in radial glia, endothelial cells, microglia, and astrocytes in the developing human cortex (Nowakowski et al., 2016). After breaching the placental–fetal barrier, ZIKV infects radial glial cells. Radial glia cells undergo symmetric and asymmetric cell division to generate the neurons that populate the cerebral cortex (Anthony et al., 2004; Falk and Götz, 2017), and the disruption of this process may ultimately lead to decreased neuronal numbers resulting in microcephaly. Although ZIKV might not be cytotoxic for microglia and astrocytes initially, once infected, these cells can propagate and spread the virus through the brain, maintaining a high viral load in the brain over time (Muffat et al., 2018).
An increasing number of studies point to the involvement of glial cells in ZIKV infection of the CNS. Microglia are well known for their role as cellular effectors of innate immunity; together with astrocytes, they drive the neuroinflammatory process through phagocytosis and cytokine release. They also play a critical role in neural development by promoting neuronal survival via release of neurotrophic factors to support neuronal circuit formation, phagocytosing immature neurons that fail to form proper neuronal circuits, and removing redundant or dysfunctional synapses in the developing brain (Thion et al., 2018; Reemst et al., 2016). During the late stages of brain development and after birth, glial cells develop and spread through the nervous system, constituting at least half of the cellular population within the brain (Zuchero and Barres, 2015). Two possible pathways exist by which ZIKV may interact with glial cells –either direct infection or induction of an inflammatory response that activates microglia and/or astrocytes. Congenitally acquired ZIKV has been shown to infect glial precursor cells, impairing their distribution in the brain, leading to reduced and delayed myelination (Chimelli and Avvad-Portari, 2018), ultimately contributing to neuroinflammation and microcephaly.
Nonhuman primate (NHP) models of ZIKV infection have been established and are able to bridge the translational divide between rodent models and humans (Nguyen et al., 2017; Coffey et al., 2018). Rodent models have significant limitations including the fact that adult mouse models of ZIKV must be immunocompromised or genetically modified (e.g., Winkler and Peterson, 2018). Fetal and neonatal wild-type mice do show some susceptibility to ZIKV infection, but they must be directly exposed to the virus, either through intracerebroventricular, intraperitoneal, or intra-amniotic (IA) inoculation, and as such are not a good model for vertical transmission of ZIKV from infected mothers to fetuses (e.g., Paul et al., 2018). In addition, the most severe symptoms of CZS include dysfunction in the neocortex, especially the prefrontal cortex (PFC), and the extent to which certain cortical regions are homologous in rodents and primates is not clear (Laubach et al., 2018). NHPs gestational development is also similar to humans, including its extended duration (relative to rodents), the prevalence of singleton pregnancies (rather than litters), and the structure of the placenta (Carter, 2007). The anatomical development and organization of the brain also proceed along a similar path in humans and NHPs (Workman et al., 2013). Similar to humans, pregnant NHPs do not need to be immunocompromised to be infected with ZIKV, experience viremia, and demonstrate transplacental transmission of the virus to the fetus (Koide et al., 2016; Dudley et al., 2016; Morrison and Diamond, 2017; Coffey et al., 2017). Although none had evidence of microcephaly, previous studies of ZIKV-exposed fetal and infant macaques have identified a number of pathological features, including calcifications, abnormal gliosis, and white matter hypoplasia (Coffey et al., 2018; Mavigner et al., 2018; Adams Waldorf et al., 2018; Adams Waldorf et al., 2016). Studies that employed immunohistochemical (IHC) analysis also identified loss of neural precursors, reduced neurogenesis, gliosis, and increased apoptosis in different brain regions (Coffey et al., 2018; Mavigner et al., 2018; Adams Waldorf et al., 2018). Neuroanatomical evaluations in these studies have adopted approaches that are standard for pathologic analyses, insofar as they typically evaluate small areas of brain tissue and use formalin fixation and paraffin embedding, procedures that can mask epitopes, reduce antigenicity, and make it challenging to understand the impact of infection on the whole brain (Cinar et al., 2006; Jones and Hartman, 1978).
The present study builds upon a previously described macaque model of CZS in which, to assure fetal infection at a defined time of gestation, fetuses were inoculated with ZIKV via the IA routes, mothers were concurrently inoculated with ZIKV intravenously (IV), and fetuses were harvested at the end of gestation (Coffey et al., 2018). We elected to use IA inoculation to ensure that the fetuses were infected so that we could study fetal outcomes, rather than inoculating the mother only and leaving open the possibility that fetuses did not become infected (i.e., at the time we initiated our study it was not clear whether vertical transmission rates were 100% although there is some subsequent evidence that is the case [Nguyen et al., 2017]). While RT-qPCR-based virology, immunology, and histology findings on these animals have been described previously (Coffey et al., 2018), the current report describes a more in-depth neuropathological assessment of the macro- and microstructural effects of ZIKV infection on the developing fetal brain using quantitative microscopy, comparing the effects of ZIKV infection in a developed cortical region (at the caudal extent of the brain, including the occipital lobe, Area 17) and an immature cortical region (at the rostral extent of the brain, including the frontal lobe, Area 46). Critically, we maintained the spatial integrity of the tissue and examined whole tissue sections across the complete hemisphere available for study; as such the analyses presented here differ from previous reports because we quantified specific anatomical regions of the brain (e.g., Areas 17 and 46) including features related to their macrostructure (e.g., cortical thickness) and microstructure (e.g., morphological analyses; identification of immature neurons in the process of apoptosis, etc.). These methods allowed us to identify where the virus was found, what cell types were affected, and whether gross neuroanatomical abnormalities were present.
Results
To identify and characterize neuroanatomical consequences of fetal ZIKV, we examined the brains of six near full-term monkeys, including three animals that had been inoculated with ZIKV at different gestational days (GD 50, 64, and 90; full term is approximately GD 165) and three procedure-matched controls (Coffey et al., 2018). The pregnant animal inoculated at GD 64 spontaneously gave birth to a small but viable baby on GD 151, at which point both mother and baby were euthanized. The remaining ZIKV-infected subjects and matched controls were monitored until GD 155 when fetectomies were performed followed by necropsy (Coffey et al., 2018). The left cerebral hemisphere was used for both analysis of viral load and histology (Coffey et al., 2018) while the right hemisphere was preserved for anatomical analyses. The right hemisphere was sectioned into four blocks to allow for more efficient immersion fixation and cryoprotection (see Figure 1A). Brains were sectioned on a freezing sliding microtome, and one series underwent Nissl staining following standard lab protocols (Kreutzberg, 1984; Lavenex et al., 2009). Nissl-stained sections from the most anterior block and most posterior block were evaluated to determine if there were gross anatomical abnormalities and if development had proceed as expected (from caudal to rostral; Workman et al., 2013). The frontal and occipital lobes were selected as initial targets for evaluation because neocortex develops along the caudal–rostral axis, and at the day of birth the occipital lobe exhibits a well-defined pattern of cortical lamination while the layers in the PFC are indistinct (Workman et al., 2013; Silbereis et al., 2016). Additionally, we visually inspected each Nissl-stained tissue section to see if there were additional macrostructural abnormalities.

ZIKV infection induce neuroinflammatory response associated with activated microglia in immature prefrontal cortex.
Unlike control animals (A), the tissue of the ZIKV-infected animals (B) evidenced clusters of round-shaped microglia (Iba1, red) and astrocytes ( Glial fibrillary acidic protein (GFAP), green) within selective cortical areas. Scale bar: 50 µM, zoom: 10 µM. In the Area 46 of the ZIKV-infected animals, the microglia marker Iba1 (red), colocalized with lysosomal protein CD68 (orange), a marker for phagocytic activity normally observed in activated microglia (C–G). Scale bar: 100 µM, zoom: 5 µM. Activated microglia surrounding apoptotic neurons (NeuN, orange) are shown within Area 46 in ZIKV animals (H, highlighted in I). 3D analysis of the micrographs shows apoptotic material (fragmented 4′,6-diamidino-2-phenylindole (DAPI) and NeuN) inside activated microglia (J, K, scale bar: 5 µM).
Neuroinflammatory abnormalities and ZIKV presence in brains from ZIKV-infected animals
IHC confocal microscopy revealed the presence of clusters of reactive glia within Area 46 of the PFC from ZIKV+ animals, but not in the controls (Figure 1A, B). Multiple labeling fluorescent microscopy allowed us to visualize that abnormal clusters of Iba1+ microglia which also express the lysosomal marker CD68 in the ZIKV-infected animals (Figure 1C–G). CD68 levels are known to be substantially upregulated in microglia during inflammatory processes and are thought to be involved in active phagocytosis (Jurga et al., 2020). We also observed that clusters of amoeboid-shaped microglia within Area 46, are directly engulfing neuronal debris in the ZIKV-infected monkeys (Figure 1H, I; 3D reconstruction: Figure 1J, K).
Next, we sought to detect the presence of ZIKV protein in the same regions we observed active inflammatory response. Interestingly, we observed the presence of ZIKV envelope protein within Area 46 of all the infected animals, with detection of viral protein inside clusters of microglia, but not astrocytes (GFAP, Figures 2A–G and 3D reconstruction: Figure 2H). Despite the fact that cells expressing ZIKV were detected in Area 46, we were not able to confirm the identity of the infected cells. As shown in Figure 2I, J, these ZIKV+ cells are thought to be in a direct process of phagocytosis, a common event occurring during viral encephalitis (Chhatbar and Prinz, 2021). While we could not detect the expression of other proteins in the ZIKV+ cells, they were surrounded and in direct contact with reactive microglia also expressing CD68, as shown in Figure 2K, L.

Intra-amniotic inoculation of ZIKV results in active brain infections of fetuses up to 105-day postinoculation.
ZIKV protein (antiflavivirus, green) (A) is identified within Area 46 of ZIKV inoculated animals. ZIKV+ cells are surrounded by activated microglia (Iba1, red) (B) and astrocytes (GFAP, purple) (C), with DAPI (D) staining all nuclei. Scale bar: 10 µM. 3D analysis and reconstruction of confocal micrographs obtained in Area 46, highlight the presence of viral protein within the microglia (E–H). Scale bar: 10 µM, 3D: 5 µM. ZIKV+ cells observed in the same brain region are surrounded and actively contacted by reactive microglia (I, J). Scale bar: 25 µM, zoom: 3 µM. Multiple labeling combining both microglia markers and viral protein shows the presence of viral protein inside phagocytic microglia in the prefrontal cortex (PFC) of the ZIKV-infected fetuses (K, L). Scale bar: 25 µM, zoom: 3 µM.
ZIKV-induced changes to the occipital and frontal lobe
Macrostructural changes
The most prominent manifestation of CZS is microcephaly (Moore et al., 2017), but none of the ZIKV-infected animals studied here exhibited a disproportionately (compared to their bodies) smaller biparietal diameter during gestation or head sizes at necropsy (Coffey et al., 2018). However, in less severe cases of CZS there have been documented changes to gross anatomical features of the brain including gyral simplification (defined as a reduction in the number of gyri associated with shallow sulci) and lissencephaly. Given those findings in humans, an analysis of the surface-area-to-volume ratio (or gyrencephality) of the occipital and frontal lobes from our infected and control animals allowed for an evaluation of these more subtle (compared to microcephaly) changes related to infection (Figure 3). These analyses revealed important differences in macrostructural features in Zika-infected versus control brains.

Fetal ZIKV infection results in variation in cortical development.
(A) Each fetal brain was divided into four blocks to ensure complete rapid fixation (top). Representative Nissl-stained sections from each block (middle), and the distributions of gray matter (red), white matter (blue), and subcortical tissue (gray) are shown for each section (bottom). For (B)–(D), data from control animals are depicted in blue circles and data from ZIKV-infected animals are depicted in red squares. Gestational day (GD) of infection is indicated with different intensities of blue and red. The mean value for each group is depicted with a horizontal line (blue for control, red for ZIKV). (B) The gyrencephality index of the occipital lobe (Block 1; left) and frontal lobe (Block 4; right) is shown for each individual subject. Areas 17 and 46 were analyzed from Blocks 1 and 4, respectively. (C) Proportion of gray matter (left) and white matter (right) in Block 1. (D) Proportion of gray matter (left) and white matter (right) in Block 4.
Gyrification of the occipital lobe was significantly reduced in ZIKV-infected animals as compared with control animals (t4 = 3.33, p = 0.029; d = 2.72) (Figure 3B). Within the occipital lobes there were no significant differences between ZIKV and control animals in the total lobe volumes (t4 = 0.46, p = 0.67; d = 0.37) or white and gray matter volumes (t4 = 1.33, p = 0.25; d = 1.09 and t4 = 0.264, p = 0.805; d = 0.22, respectively). However, differences between the proportions of white and gray matter varied by condition. The occipital lobes from ZIKV animals had a significantly higher proportion of white matter and lower proportion of gray matter than the control animals (t4 = 2.93, p = 0.043; d = 2.39; because the occipital lobe included only gray and white matter, the proportions are inverses of each other and the statistics identical) (Figure 3C). To evaluate cortical thickness, we selected one clearly defined anatomical region completely contained within the occipital lobe – Brodmann’s Area 17. Cortical thickness in that area did not differ significantly between ZIKV and control subjects (t4 = 1.42, p = 0.230; d = 1.16; Figure 3A), although all the data points for the ZIKV group were below the mean of the control group, suggesting the possibility of cortical thinning within the primary visual cortex.
There was no difference between ZIKV and control animals in the gyrification of the frontal lobes (t4 = 0.56, p = 0.60; d = 0.46) (Figure 3B). The total volume of the frontal lobes did not differ between ZIKV and control animals (t4 = 1.20, p = 0.30; d = 0.98), and neither did the white and gray matter volumes (t4 = 0.35, p = 0.75; d = 0.28 and t4 = 1.53, p = 0.20; d = 1.25, respectively), nor relative proportions of white and gray matter (t4 = 1.80, p = 0.15; d = 1.47 and t4 = 0.59, p = 0.59; d = 0.48, respectively; Figure 3D). As in the occipital lobe analyses, we selected one anatomical region completely contained within the frontal lobe – Brodmann’s Area 46 – and carried out a cortical thickness evaluation of that entire area. There were no significant differences between cortical thickness in Area 17 (Figure 4A) and in Area 46 (Figure 4B) of ZIKV and control animals (t4 = 0.57, p = 0.60; d = 0.46; Figure 4B). Taken together, these findings suggest that ZIKV infection did not cause gross morphological changes to frontal cortex during development.

Cortical thickness does not differ across immature versus mature cortical regions of the brain in ZIKV-infected and noninfected animals.
The cortical thickness of a mature cortical region (Area 17, primary visual cortex, V1) (A) did not differ between the control (blue circles) and ZIKV-infected (red squares) animals. The anatomical localization of Area 17 is shown in a photomicrograph of a Nissl-stained section stained (dashed lines = the extent of Area 17). (B) Cortical thickness of Area 46 did not differ between the control (blue circles) and ZIKV-infected animals, although there was more individual variation than was seen in Area 17. The anatomical localization of Area 46 (dorsolateral prefrontal cortex, DLPFC) is shown in a photomicrograph of a Nissl-stained section (dashed line = the extent of Area 46). Conventions as in previous figures.
M1 cortical thickness, and volumes of the amygdala, hippocampus, and optic tract
We carried out a few additional gross level analyses of one additional cortical target and three subcortical targets. We evaluated the cortical thickness of one additional cortical area, Area 4, primary motor cortex, was included because it is roughly in the middle of the brain (on the rostral–caudal axis) and thus represents a cortex at a different point in development compared to Areas 17 and 46. Three thickness measurements were taken on each of three different slides and then the value were averaged. There was not a difference in cortical thickness between the Zika and control animals in Area 4 (t4 = 0.80, p = 0.47; d = 0.65; MeanZIKV = 2427.76 μm, SDZIKV = 95.64 μm; MeanControl = 2361.72 μm, SDControl = 106.62 μm).
Given the senior author’s lab’s focus on the amygdala and hippocampus, we also analyzed volumes of those structures. There were no differences between the ZIKV and control groups in the volumes of the amygdala, hippocampus. Amygdala: t4 = 0.095, p = 0.93; d = 0.078; MeanZika = 79.85 mm3, SDZika = 3.88 mm3; MeanControl = 79.28 mm3, SDControl = 9.45 mm3. Hippocampus: t4 = 0.58, p = 0.59; d = 0.47; MeanZika = 135.82 mm3, SDZika = 10.65 mm3; MeanControl = 143.97 mm3, SDControl = 21.84 mm3.
Finally, given the impact of ZIKV on the lateral geniculate nucleus (LGN) and observations about its impact on the retina, we also carried out a volumetric analysis of the optic tract. There were no group differences between the ZIKV and control animals in the volume of their optic tracts: t4 = 1.061, p = 0.35; d = 0.87; MeanZika = 31.24 mm3, SDZika = 6.54 mm3; MeanControl = 25.09 mm3, SDControl = 7.62 mm3.
Microstructural changes to glia
Following macrostructural evaluations, we carried out a series of IHC analyses to quantify cell-level features that might be impacted by ZIKV infection with a specific focus on glia (microglia and astrocytes) in Brodmann’s Areas 17 and 46. The three-dimensional space glia occupy, especially microglia, change constantly, because these cells are highly dynamic and retract and expand their processes in response to subtle changes in the surrounding environment. Several studies have shown that glial cells have a strong form-to-function mapping and investigating morphological alterations in these cells allows inferences about their activation state and their inflammatory status across different brain regions (Karperien et al., 2013; Morrison and Filosa, 2013; Young and Morrison, 2018). For example, using the combination of high-resolution confocal microscopy with three-dimensional reconstruction of individual glia volumes, we recently described how morphological changes in these cell types are connected with spine loss and neuron death in two different monkey models of Alzheimer’s disease (Seelke et al., 2020; Beckman et al., 2021). Here, we use a similar approach of analysis by randomly selecting 28 cells from each anatomical area, for each animal. Iba1 and GFAP were used as general microglia and astrocytes markers, respectively, and each cell was exported and analyzed individually for total and cell body volumes, and for number of terminal point (branching ramification).
Analysis of the microglia in Area 17 revealed no significant differences in the cell body size (t4 = 0.36, p = 0.74; d = 0.29; Figure 5C). The whole cell volumes were marginally smaller in ZIKV-infected animals than in control animals and while the p value did not reach conventional levels of significance, the effect size was very large (t4 = 2.71, p = 0.054; d = 2.11; Figure 5D) There were no group differences in the number of terminal points (t4 = 0.83, p = 0.45; d = 0.68; Figure 5E). In contrast, while there were no significant differences between the cell body size of microglia in Area 46 of ZIKV-infected compared to control animals (t4 = 2.07, p = 0.11; d = 1.69; Figure 5M), ZIKV-infected animals, compared to controls, had significantly smaller whole cell volumes (t4 = 9.41, p < 0.001; d = 7.68; Figure 5N) and significantly fewer terminal points (t4 = 4.18, p = 0.01; d = 3.41; Figure 5O). The smaller whole cell volumes and terminal points suggest that the microglia were in their activated state in frontal cortex.

Glial cell morphological changes induced by ZIKV infection of the brain.
For each cell type (microglia and astrocytes), 28 cells were randomly selected from Area 17 and from Area 46 for each animal. 3D confocal images of each individual cell were exported to Imaris software and analyzed for cell body and whole cell volume, and total number of terminal points. Representative microglia from Area 17 in control and ZIKV-infected animals are depicted in (A) and (B). The average for each condition analyzed is represented for each animal in both groups. There were no group differences in cell body volume (C), whole cell volume (D), or number of terminal points (E), a measurement for glia complexity. Similarly, there were no group differences in astrocyte anatomy (F and G for representative anatomy, H–J for data). In contrast, there were group differences in Area 46 (representative images of Area 46 microglia, K and L; astrocytes, P and Q). The cell bodies of ZIKV-infected animals’ microglia were larger in size (M) than control animals. Compared to control animals’ microglia, ZIKV-infected animals’ microglia whole cell volumes were smaller (N) and had fewer of terminal points (O). Compared to controls animals’ astrocytes, ZIKV-infected animals astrocytes tended to have smaller cell body volumes (R), no differences in whole cell volumes (S), and a greater number of terminal points (T).
There were no significant group differences in astrocyte morphology in either Area 17 or 46. In Area 17, astrocyte cell body size did not differ between ZIKV-infected and control animals (t4 = 0.34, p = 0.75; d = 0.28; Figure 5H), whole cell volume (t4 = 2.13, p = 0.09; d = 1.81; Figure 5I), or number of terminal points (t4 = 0.44, p = 0.68; d = 0.34; Figure 5J). In Area 46, astrocyte cell body size did not differ between ZIKV-infected and control animals (t4 = 0.10, p = 0.93; d = 0.81; Figure 5R), whole cell volume (t4 = 1.03, p = 0.36; d = 0.84; Figure 5S), or number of terminal points (t4 = 1.22, p = 0.29; d = 0.99; Figure 5T).
Microstructural changes related to neural development
A neuroinflammatory response, mainly driven by CD68/Iba+ microglia, in combination with the presence of ZIKV envelope protein, indicated that the frontal lobe for all ZIKV-infected subjects was a site of persistent ZIKV infection for a minimum of 60-day postinoculation. However, the consequences of the infection persist beyond the window of virus replication. Previous studies have demonstrated that ZIKV infection is associated with increased apoptosis of neural progenitor cells (Li et al., 2016; Souza et al., 2016). IHC fluorescent microscopy analyses for DAPI, cleaved caspase-3 (CC3), and SATB2, which are markers of nuclei, apoptosis onset, and immature neurons, respectively, were carried out in order to determine if the increased inflammatory response or active viral presence could induce neuronal loss across the regions analyzed 3D analysis of the average total number of cells (DAPI+), average number of immature neurons (SATB2+), and total number of apoptotic cells (CC3+), did not significantly differ between ZIKV and control animals in both regions analyzed (Figure 6F – Area 17, Figure 6R – Area 46). The number of cells in which CC3 and SATB2 was colocalized did not differ across the ZIKV-infected and control groups and the numbers were fairly low (t4 = 0.65, p = 0.55; d = 0.53; Figure 5A–K), indicating low frequency of immature neuron death. In contrast, there was significantly more frequent death of immature neurons in Area 46 of the ZIKV-infected compared to control animals, as indicated by a greater number of cells that expressed both CC3 and SATB2 (t2.20 = 5.37, p = 0.027; d = 4.39; Figure 6L–Z). Many CC3+ cells in the Area 46 of ZIKV animals are constituted by apoptotic immature neurons, higher resolution 3D analysis showed that activated microglia (HLA-Dr+) also present high expression of cell death marker CC3+ (Figure 6—figure supplement 1). The higher frequency of immature neuron death, the active immune response (activated microglia), and the presence of ZIKV protein in Area 46 suggest that frontal cortex remained a site of active ZIKV-induced neuronal remodeling at the time the brains were analyzed. Fragmented DNA was also observed in this area (Figure 6—figure supplement 2).

Increased death of immature neurons in Area 46 of ZIKV-infected animals.
Representative tissue samples from Area 17 (A–L) and Area 46 (M–Z) in control and ZIKV-infected animals. Immature neuronal marker SATB2 (A, G, M, S) colocalization with CC3 (B, H, N, T), an apoptotic marker, and DAPI (C, I, O, U) were analyzed in Areas 17 and 46 of control and ZIKV-infected animals. Merged images (D, J, P, V). 3D surface rendering reconstruction of SATB2 and CC3 are shown in E, K, Q, and X. Average total cell counting (F, R) and quantification of SATB2/CC3 colocalization (L, Z) indicated no group differences in the number of immature neurons undergoing apoptosis in Area 17 (L), but a significantly greater immature neurons undergoing apoptosis in ZIKV-infected compared to control animals in Area 46 (R). Scale bar = 25 µm.
Additional neuronal abnormalities
Visual inspection of the Nissl-stained sections of each brain revealed other ZIKV infection-related abnormalities within the visual pathway in addition to those observed in Area 17 (discussed above). These abnormalities included structural changes to the LGN. The LGN is a six-layered structure that is found in the ventral thalamus that receives retinal input from both eyes via the optic tract (Jones, 1985; e.g., Figure 7A). In the ZIKV-infected subjects, we observed multiple instances of laminar discontinuities within the LGN (Figure 7B–D). The LGNs of the control subjects were all normal.

Anatomical abnormalities within the lateral geniculate nucleus (LGN) of ZIKV-infected subjects.
(A) An example of normal caudal LGN from the control animal who was sham inoculated at GD 90. Images (B–D) illustrate LGN abnormalities from each of the three subjects infected with ZIKV. Black ovals indicate areas of anatomical abnormalities within each section. (B) Two instances of anatomical abnormalities within the animal infected with ZIKV at GD 50. The top oval highlights an area of undifferentiated lamination, where layers 1 and 2, and 3 appear to blend together. The bottom oval highlights an area where the normal laminar structure is interrupted by a region of low cell density. (C) Two instances of anatomical abnormalities within the animal infected with ZIKV at GD 64. The top oval highlights an area of undifferentiated lamination, where layers 2 and 3 appear to blend together. The bottom oval highlights an area where the normal laminar structure is interrupted by a region of low cell density. (D) An instance of laminar interruption of layers 1, 2, and 3 by a region of low cell density in the third ZIKV-infected animal (infected at GD 90). (E) The proportion of the histological section occupied by the LGN did not significantly differ between control animals and the ZIKV-infected animals, t4 = 1.00, p = 0.37, d = 0.81. There were no group differences in the proportion the (F) magnocellular (t2.08 = 1.02, p = 0.36, d = 0.84), (G) parvocellular (t4 = 0.14, p = 0.90, d = 0.11), and (H) koniocellular layers (t4 = 0.57, p = 0.60, d = 0.47). Scale = 250 µm.
The first type of LGN pathology consisted of a blurring of the expected boundaries between cell layers. In healthy tissue, the layers of the LGN are separated by a thin layer of koniocellular tissue. However, each ZIKV-infected animal had discrete instances where the koniocellular layer was absent, resulting in indistinct boundaries between layers (Figure 7B–D). The location of the restricted absence of the koniocellular layer varied across cases. These abnormalities occurred between both magnocelluar and parvocellular layers and could be found in both the dorsomedial and ventrolateral aspects of the structure. To investigate the differences observed in the LGN of ZIKV-infected animals, fluorescent microscopy was performed in adjacent sections of the Nissl-stained LGN (Figure 8A–L). In comparison to the control animals (Figure 8B–F), in the ZIKV-infected animals, there was an increase in the frequency and intensity of the markers for microglia (Iba1) and astrocytes (GFAP) across different layers of the LGN (Figure 8H–L). Although these data suggest the possibility of active inflammatory response occurring, no active virus was observed in this region. The LGN sends direct projections to layer 4 of the primary visual cortex (V1, Area 17) in the occipital lobe. We identified two categories of anatomical abnormalities in the occipital lobe of ZIKV-infected animals (but not controls): laminar disruptions and cortical heterotopias. Laminar disruptions were identified within Area 17 and consisted of a break in the continuity of layer 4a (Figure 9A–R). These disruptions appeared to be specific to layer 4a and not associated with the normal transition from Area 17 to 18, because layers 4b and 4c were still present immediately proximate to the abnormal layer 4a. This disruption is likely functionally significant, as layer 4a receives direct projections from the LGN (Thomson and Bannister, 2003; Miller, 2003), and the absence of cells that make up layer 4a likely represent lost visual information. Cortical heterotopias were identified as ectopic islands of neurons extending into cortical layer 1 (Figure 9—figure supplement 1).

Neuroinflammatory abnormalities in the lateral geniculate nucleus (LGN) of ZIKV-infected animals.
The LGNs of control animals were all normal (A–F), but structural abnormalities were detected in the LGNs of all three ZIKV-infected fetuses in the Nissl-stained sections (an example depicted in G–L). A tissue section adjacent to each Nissl stained was analyzed using confocal microscopy. Neuronal layers and glia were analyzed using NeuN (B, H), Iba1 (C, D; red), and GFAP (E, K; green). Zoomed NeuN images in the same area are in D and J. Merged images are found in F and L. Scale = 250 µm in G–K, scale = 25 µm in B–F and L–P, scale = 10 µm in Q–U.

Anatomical abnormalities within the primary and secondary visual cortex of a ZIKV-infected animals.
Compared to normal transitions between primary (V1) and secondary (V2) visual cortex (A), cellular abnormalities in the transition were detected in ZIKV-infected monkeys (B, C) on the Nissl-stained sections. A deeper investigation using quadruple fluorescent staining (merged images in D–F) was performed in adjacent sections for SATB2 (G–I), Iba1 (J–L), and GFAP (M– O). A normal transition between V1 and V2 is characterized by the uniform contraction of layer 4 and the absence of laminar discontinuity (A, C, E, G, I, K). We identified a region of laminar discontinuity within the boundaries of V1 (B) that appears to be the caused by a disruption to layer 4a that resulted in the displacement of layers 4b and 4c toward the pial surface (seen in M–R). An additional abnormality was found below layer six in area V2 (seen in A–K and magnified in B–L). Disrupted neuronal layers present are surrounded by robust presence of activated microglia and astrocytes, as highlighted in the zoom area N–R. Scale = 250 µm in A–L, scale = 25 µm in M–R.
Discussion
Our results demonstrate that ZIKV infection during fetal development results in central nervous system abnormalities that track with temporal patterns of brain development. Importantly, our previous report on these animals (Coffey et al., 2018) documented evidence of maternal and fetal virus replication consistent with evidence presented here that ZIKV directly infected the brains of the fetuses. In the current study, virus replication and local neuroinflammatory response were observed in the neocortex of all three ZIKV-infected fetal brains. The types and magnitudes of these abnormalities varied across the cortex, with altered macrostructural features in the occipital cortex and microstructural abnormalities found in the frontal regions of the cortex. These regional differences likely reflect the graded manner in which the neocortex develops, with ZIKV infection tracking temporally with neural development. Importantly, the neuroanatomical abnormalities that we observed in all three ZIKV-infected fetal brains, but not in control brains, occurred even in the absence of the extreme ZIKV infection phenotype documented in human infants, microcephaly. These findings may help explain why some infants born with normal sized heads during the ZIKV epidemic have developed significant sensory, motor, and sociocognitive challenges as they grow up (Wheeler et al., 2018; Cardoso et al., 2019; Ventura et al., 2017; Vianna et al., 2018; Marques et al., 2019; Pinato et al., 2018).
The neuroanatomical abnormalities that we identified here appeared to follow the pattern of otherwise normal neuronal development. In mammals the neocortex develops in a gradient. Neural progenitor cells undergo migration in a caudal-to-rostral pattern, resulting in occipital regions having mature laminar organization before mature laminar organization is present in frontal regions (Workman et al., 2013; Charvet and Finlay, 2014). Rhesus monkeys, like humans and NHPs, have frontal cortices that do not reach maturity until well after birth (Bourgeois et al., 1994; Petanjek et al., 2011; Mrzljak et al., 1991). We observed macrostructural differences (e.g., relative gray/white matter volumes) in the occipital lobe which develops first and microstructural and immune-related differences in frontal lobe which develops later. Our analyses revealed evidence that microglia were in their activated state in Area 46 of ZIKV-infected but not control animals, as indicated by smaller whole cell volumes and fewer terminal points. Microglia play a critical role in neural development via a number of mechanisms (including promoting neuronal survival and removing redundant or dysfunctional synapses) and morphological differences across the groups may be evidence of pathological neuronal development. Supporting the idea that ZIKV infection interferes with neural development was our finding that ZIKV-infected fetuses had a greater number of immature cortical neurons undergoing apoptosis compared to control animals. That we found these differences in Area 46, which was still undergoing development at the time of tissue harvest, and not Area 17, which was fully laminated at the time of tissue harvest, underscores the ability of ZIKV to interfere with normative neural development. One interpretation of these data is that at the point at which the tissue was harvested, damage had already been done to the part of the brain that was differentiated but the processes of damage were underway in the tissue that was actively developing. These results suggest that infants infected with ZIKV during fetal development and born with normal sized brains may have neurodevelopmental abnormalities that only become apparent later in development. This finding is of particular importance because most infants born during the 2015–2016 ZIKV epidemic in the ZIKV hotspots were born with normal sized heads and microcephaly, while an extreme phenotype, is thought to be relatively low frequency (Cardoso et al., 2019; Aragao et al., 2017; Faiçal et al., 2019).
One seemingly consistent anatomical abnormality found in both human patients and animal models of ZIKV is an increased frequency of abnormalities in the visual system (Aleman et al., 2017; Fernandez et al., 2017; Miner et al., 2016b; Simonin et al., 2019; Singh et al., 2018; Ventura et al., 2019; Ventura et al., 2016; Verçosa et al., 2017), especially the retina. Patients present with severe eye disease, including congenital cataracts, optic nerve abnormalities, chorioretinal atrophy, intraretinal hemorrhages, mottled retinal pigment epithelium, and blindness (Ventura et al., 2017; Aleman et al., 2017; Fernandez et al., 2017; Ventura et al., 2019; Ventura et al., 2016; Verçosa et al., 2017; Miner et al., 2016a; Ventura et al., 2018). Further, there is some evidence that ZIKV can be transmitted from the retina to the rest of the visual system via the optic nerve, although this has yet to be confirmed in humans or NHPs (van den Pol et al., 2017). Anatomical abnormalities that we observed in the LGN and Area 17/V1 are consistent with infection and developmental disruption of the eye and optic nerve, insofar as projections from the retina enervate the LGN in a retinotopic manner, and projections from the LGN in turn enervate the primary visual cortex (Brodmann’s Area 17) (Shatz, 2018; Rakic, 1977). The anatomical and functional retinotopic organization of both LGN and V1 is dependent upon activity during developmental critical periods (LeVay et al., 1980; Ghosh and Shatz, 1994). Insults during early development that disrupt this activity can result in damage to brain regions downstream of the initial insult (Ito et al., 2008; Eysel and Wolfhard, 1984; Yücel et al., 2003). The lesions and laminar dysregulations observed in the LGN and V1, respectively, have significant functional consequences.
While routine H&E histopathology of the eyes did not reveal any significant lesions in these animals (Coffey et al., 2018), subsequent work from our group has identified retinal pathology in one infant that was infected with ZIKV IA as a fetus and monitored for 2 years after birth (Yiu et al., 2020). The retinal tissue was used for pathological analysis and was unavailable for use in this study. However, the retinotopic organization and activity-dependent development of LGN layers suggests that the dysregulation we see within the LGN is likely related to lesions of the retina caused by the Zika virus. Future studies would benefit from examination of the entire visual pathway, from retina to LGN to V1, and visual association regions targeted by V1.
While our sample size was small (N = 3 in each group), only one hemisphere of the brain was available for our analyses, and the timing of ZIKV infection varied across animals (with one ZIKV and one procedure-matched control at each of three infection dates, GD 50, 64, or 90), all three infected subjects showed the same pattern of pathology: active inflammation in the still-developing regions of the frontal cortex paired with anatomical abnormalities in the already-developed regions of the caudal cortex. Additionally, each infected fetus had evidence of active virus replication as late as 105-day postinoculation, which could explain the ongoing neurologic damage. Although the sample size here is not atypically small for neuroanatomical analyses of NHP tissues or initial infectious diseases studies of this sort, it is small and thus important to note that we carried out both experimental and statistical procedures to address potential sample size concerns. First, each infected animal was paired with a procedure-matched control that underwent identical procedures without viral inoculation. Second, all conducted statistical analyses are reported (i.e., there was no p-hacking). Third, raw data but not statistics are displayed in figures encouraging the reader to draw conclusions from the patterns of data within and across groups rather than relying solely on p values. While replication in a larger sample is a laudable future goal, the present study also sets the stage for targeted hypothesis testing about the timing of infection and gestational development which seems to matter less in the case of ZIKV than potentially anticipated.
Previous studies from our group (Coffey et al., 2018) and others (Adams Waldorf et al., 2018; Martinot et al., 2018; Mohr et al., 2018; Steinbach et al., 2020) have identified discrete histopathological abnormalities in the brain tissue of macaque fetuses infected with ZIKV, including calcifications with gliosis and a loss of neural precursor cells. In addition to previous neuropathological evaluations, previous evaluation of the other hemisphere of each of these brains identified ZIKV RNA in Area 46 in all infected animals, and in occipital pole (Area 17) of animals infected at GD 64 and 90 (Coffey et al., 2018). The histological analyses presented here differ in that we measured specific cytoarchitectonic regions of the brain as well as morphological analyses of glial cells in the frontal and occipital cortex. Together, these analyses paint a picture of inflammatory processes related to viral infiltration sweeping through the cortex and a wave of cell death. Given our analyses, it is not possible to tell whether inflammatory processes are caused by or are causing cell death, but they certainly appear to co-occur given their proximity. Importantly, in addition to underscoring the importance of adopting quantitative neuroanatomical analyses to understand viral infections of the developing brain, these findings underscore the importance of using a monkey model because we would not expect to see such patterns in nonprimate models, given that rodents do not have a homolog of Area 46 and a significant period of their brain development occurs after birth.
While we focus heavily on the importance of homologs between humans and rhesus monkeys both in the course of ZIKV infection and the development of PFC (including Area 46), the rhesus monkey model opens the possibility of linking neural development to tractable, translatable behavioral development, and also developing translatable interventions to improve neural and behavioral function. Rhesus monkeys are a robust model for human infant development across a variety of behavioral domains – including cognition, affect, and social behavior (Phillips et al., 2014; Clancy et al., 2013; Fox, 1984; Cassidy and Shaver, 1999). Moving forward, this will allow us to understand how the neural abnormalities associated with fetal ZIKV infection, such as those documented here, impact infant behavioral development. Because rhesus monkeys develop three to four times faster than humans (Workman et al., 2013), a separate, currently ongoing cohort of infants who were infected with ZIKV as fetuses will soon out-age human infants who were infected in utero, thus providing a prospective look at the challenges those infants may face. In this context, our long-term developmental studies which have documented significant neural and behavioral plasticity following early brain damage (e.g., Grayson et al., 2017; Bliss-Moreau et al., 2017) provide hope that fairly simple behavioral interventions may be possible to improve the outcomes for impacted children. This is particularly important given that the ZIKV epidemic disproportionally impacted people of color living in socioeconomically disadvantaged communities (Moreno-Madriñán et al., 2017; Kapiriri and Ross, 2020; Gómez et al., 2018) and that ZIKV is expected to return seasonally in the future (Tesla et al., 2018).
Materials and methods
All procedures were approved by the University of California, Davis Institutional Animal Care and Use Committee (Protocol # 19211) which is accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC). Animal care was performed in compliance with the 2011 Guide for the Care and Use of Laboratory Animals provided by the Institute for Laboratory Animal Research.
Study design
Subjects
All adult female rhesus macaques (Macaca mulatta) enrolled in the study were from the conventional breeding colony and were born at CNPRC. Animals were determined to be free of West Nile virus (WNV) via simian WNV ELISA (Xpress Bio), SIV-free, type D retrovirus-free, and simian lymphocyte tropic virus type one antibody negative. Animals were selected from the timed breeding program based on the date of their mating and ultimately the GD of their fetus determined by comparing data obtained on ultrasound to CNPRC colony data for fetal development. All animals had at least one and up to six previous successful pregnancies with live births. ZIKV-infected animals were enrolled first, followed by control animals. The animals were housed indoors in standard stainless-steel cages (Lab Products, Inc) that met or exceed the sizing required by NIH standards, in rooms that maintained 12 hr light/dark cycles (0600 lights on to 1800 lights off), with temperature controlled between 65–75°F and 30–70% humidity. Monkeys had free access to water and received commercial chow (high protein diet, Ralston Purina Co.) twice a day and fresh produce multiple times a week. All but one female was paired with a compatible social partner.
Experimental infection with Zika virus
Request a detailed protocolFour pregnant females and their fetuses were inoculated with a 2015 Brazilian virus isolate (strain Zika virus/H. sapiens-tc/BRA/2015/Brazil_SPH2015; genbank accession number KU321639.1) by injecting 5.0 log10 PFU in 1 ml RPMI-1640 medium both IV and IA (Coffey et al., 2018). Each animal was inoculated only once at either 41, 50, 64, or 90 GDs. The fetus of the female inoculated at GD 41 was found dead 1 week later and its brain was autolyzed (Coffey et al., 2018). Three mother–fetus pairs were sham inoculated with 1 ml of RPMI-1640 medium both IV and IA at GD 50, 64, and 90 to serve as procedure-matched controls.
No formal sample size estimation was carried out because this was an initial pilot study to determine whether rhesus monkey fetuses could be infected with ZIKV (Coffey et al., 2018). When the fetus inoculated at GD 41 died, we did not include a control fetus with a sham inoculation at GD 41, resulting in the final sample of N = 3 in each experimental group.
Histological procedures
Tissue harvesting and fixation
Request a detailed protocolAs experimentally planned, the fetuses were harvested at GD 155 via fetectomy and immediately euthanized for detailed tissue collection for five of six pairs. The fetus that was infected at GD 64 was born alive at approximately GD 151 and euthanized the same day for tissue collection. Only evaluation of fetal tissues is reported here.
The left hemisphere of each brain was used for viral load determination or fixed in 10% buffered formalin for histology as described earlier (Coffey et al., 2018). The right hemisphere was fixed in 4% paraformaldehyde in 0.1 M sodium phosphate buffer for an hour, cut into four blocks to allow sufficient rapid penetration of the fixative into the tissues, and then fixed for 48 hr at 4°C with rotation. Blocks were numbered 1–4, starting at the caudal extent of the brain (Block 1, including the occipital lobe) and ending at the rostral extent of the brain (Block 4, including the frontal pole). The tissue was cryoprotected in 10% glycerin with 2% DMSO in 0.1 M sodium phosphate buffer overnight then 20% glycerin with 2% DMSO in 0.1 M sodium phosphate buffer for 72 hr (Rosene et al., 1986). Blocks were then frozen in isopentane following standard laboratory procedures (Bliss-Moreau et al., 2017).
Tissue sectioning
The tissue was sectioned on a sliding freezing microtome (Thermo Scientific Microm HM430, Waltham, MA, USA) into eight series (7 at 30 µm and 1 at 60 µm). 30 µm tissue sections were placed in a cryoprotectant solution and stored at −20°C. The 60 µm tissue sections were postfixed for 2 weeks in 10% formalin and stored at 4°C.
Nissl staining
Request a detailed protocolThe 60 µm sections were mounted on gelatin subbed slides and Nissl stained using 0.25% thionin (according to our standard protocols [Lavenex et al., 2009; Bliss-Moreau et al., 2017]). Slides were coverslipped using DPX mounting medium (Millipore Sigma, St. Louis, MO, USA). These sections were scanned (TissueScope LE; Huron Digital Pathology; St. Jacobs, ON, Canada) and digital images were used for analyses. Each image was coded to keep evaluators blind to experimental condition. Anatomical boundaries were determined by comparing Nissl-stained tissue to reference atlases (Saleem and Logothetis, 2012; Rohlfing et al., 2012), and all anatomical analyses were performed using StereoInvestigator (MBF Bioscience, Williston, VT, USA).
Immunohistochemistry
Request a detailed protocol30-µm-thick free-floating sections were incubated in an antigen retrieval solution (Wako, S1700) at 60°C for 30 min. Sections were incubated in blocking solution: 5% donkey serum, 5% goat serum, 5% bovine serum albumin in phosphate-buffered saline (PBS) 0.3% Triton for 2 hr at room temperature (RT) under agitation. Sections were then incubated overnight with antibodies to Iba1 (Wako, code 019-19741, 1:500), NeuN (Synaptic Systems, code 266-004, 1:1000), GFAP (Abcam, code ab4674, 1:1000), antiflavivirus group (Millipore, code MAB10216, 1:400), nonphosphorylated neurofilament (SMI32, Biolegend, code SMI-32P, 1:500), CC3 (Cell Signaling, code 9661, 1:500), CD68 (Abcam, code ab53444, 1:300), HLA-Dr (ThermoFisher, code MA5-11966), and SATB2 (Cell Signaling, code 39229, 1:500). Tissue was washed thoroughly with PBS and incubated with Alexa Fluor secondary antibodies (Invitrogen, 1:500) for 2 hr, at room temperature. Finally, sections were incubated with DAPI for 10 min after washing the secondary antibodies. Sections were mounted on microscope slides and coverslipped with Prolong Diamond Antifade (Invitrogen).
Neuroanatomical evaluations
We first reviewed each Nissl slide from each case to qualitatively identify potential differences across brains. Given hypotheses about ZIKV’s particular impact on developing neurons, we elected to quantitatively compare Areas 17 and 46, delineated from Blocks 1 and 4, respectively, because of the caudal–rostral developmental gradient of cortical development but also report evaluations of cortical and subcortical structures and those that appeared to be compromised by ZIKV infection (e.g., LGN). Macrolevel features of the tissue were quantified on the Nissl sections (including gyrencephalization, cortical thickness, etc., described below). Finally, a series of IHC analyses on targeted regions of tissue were completed. All people analyzing tissue were blind to condition.
Analysis of Nissl-stained tissues
Request a detailed protocolTo compute an index of gyrencephalization, we measured the area and perimeter of each section of tissue in StereoInvestigator. These areas were chosen because the cortex develops in a caudal to rostral pattern (Charvet and Finlay, 2014; Rakic, 2002; Colby et al., 2011); brain development at the time of tissue harvest should have been most complete in the occipital lobe (caudal) and least complete in the frontal lobe (rostral). Thus, if fetal ZIKV infection disrupted development in a timing specific way, we would expect to see the greatest difference in neuroanatomical patterning when these lobes were compared. The distance between each tissue section (section thickness + intersection interval) was then used to determine the surface area and volume of each section as follows:
Surface area and volume values were summed for each subject to determine the total surface area and total volume of the frontal and occipital lobes. The gyrencephality ratio for both lobes was computed by dividing the surface area (in mm2) by the volume (in mm3).
Gray matter and white matter volumes were computed on the same Nissl-stained sections in the occipital and frontal lobes. Gray matter was defined as cortical areas that exhibited clear lamination; white matter was defined as fiber tracks adjacent to the cortex. The outline of the section, areas of gray matter, and areas of white matter were traced in StereoInvestigator. Subcortical tissue, including the caudate and putamen and the rostral extent of the hippocampus, was present in a few sections, but these areas were not included in this analysis. The total area and volume of gray and white matter for each section were computed (as above). The proportion of each lobe that was composed of white matter and gray matter was computed by dividing the total volume of gray and white matter in each lobe by the total volume of the lobe. The occipital lobe included no subcortical structures and so the total area was occupied by either gray or white matter. The frontal lobe included subcortical structures and so the percentage of gray matter plus the percentage of white matter did not sum to 100.
Cortical thickness was computed for three areas: Brodmann’s Area 17 (primary visual cortex, V1), Brodmann’s Area 4 (primary motor cortex, M1), and Brodmann’s Area 46 (dorsolateral PFC). Each region was chosen because of its connection with a sensory input or motor output that is known to be negatively impacted by ZIKV infection and because Area 17 is in Block 1 at the caudal extent of the brain and Area 46 is in Block 4 at the rostral extent of the brain. Humans with CZS have been documented to have vision deficiencies (Ventura et al., 2017) as well as muscle spasms and hypertonia (Pessoa et al., 2018), implicating V1 and M1, respectively. Similarly, some evidence exists that children with CZS have symptoms consistent with intellectual developmental disorders and executive function disorders that implicate Area 46 (Cardoso et al., 2019; van der Linden et al., 2016).
Boundaries of the cortical regions of interest were identified through patterns of cortical lamination and other architectonic landmarks (Saleem and Logothetis, 2012; Rohlfing et al., 2012). Cortical thickness measurements were taken perpendicular to the pial surface. In Areas 17 and 4, three thickness measurements were taken from each of three adjacent sections. In Area 46, four thickness measurements were taken from each of two adjacent sections. Measurements were first averaged within a section, then between sections, for each subject.
Four sections from each subject that included the full LGN were evaluated. This procedure was required because the brains were blocked through the LGN in five of our six subjects. For each subject, the four largest complete sections were selected from the available complete sections. Total LGN area was measured on the Nissl-stained sections in StereoInvestigator and then the area of magnocellular (layers 1 and 2) and parvocellular (layers 3–6) was measured. Koniocellular area for each section was computed as the total area minus the magnocellular and parvocellular areas.
Analysis of tissues prepared via IHC procedures
Request a detailed protocolAll analyses of fluorescently stained tissue were conducted using a Zeiss LSM 800 confocal microscope (Carl Zeiss AG, Oberkochen, Germany).
Microglia and astrocyte quantification
Request a detailed protocolFor microglia and astrocyte morphological analyses, 28 microglia and 28 astrocytes were selected randomly from Areas 17 and 46 of each animal. For each cell, a z-stack at ×63 magnification was made, and the image was exported to Imaris software (Oxford Instruments, Bitplane Inc, Concord, MA, USA) to create 3D volume surface rendering. Cell body volume, whole cell volume, and terminal points were quantified for each cell and averaged within each subject.
Total DAPI, SATB2, and CC3 counting and SATB2 colocalization with CC3
Request a detailed protocolThree z-stack images at ×20 magnification were obtained within both Areas 46 and 17 for each subject. Images that contained both immature neurons (SATB2) and apoptotic cells (CC3) were exported to Imaris software to quantify colocalization.
DNA fragmentation
Request a detailed protocolApoptotic cells were analyzed using DeadEnd Fluorometric TUNEL System, from Promega (cat no. G3250). Four slides from each animal were incubated with equilibration buffer for 10 min. Slides were then incubated with TUNEL reaction mix for 60 min at 37°C in a humidifier chamber, avoiding exposure to light. The reaction consists in a catalytical incorporation of fluorescein-12-dUTP at 3-OH DNA end using Terminal Deoxynucleotidyl Transferase, forming a polymeric tail using the TUNEL principle. The reaction was stopped using SSC (NaCl and sodium citrate) and slides were washed three times with PBS, incubated with DAPI solution for 10 min and then mounted using Vectashield (Vector Lab). For each animal, one slide was incubated with RQ1 RNase-Free DNase (Prometa cat no. M6101) used as a positive control. Three images from each area analyzed were acquired at ×20 magnification for each slide for each animal.
Data analysis
Request a detailed protocolData were analyzed in SPSS 26.0 (IBM Corp). Data were checked for normality via the Kolmogorov–Smirnov and Shapiro–Wilk tests. Data were then analyzed using independent sample t-tests with experimental group as the between subject variable. Cohen’s d was computed as a measure of effect sizes. No individual data points were excluded from analysis (e.g., because they were outliers) and all conducted analyses are reported in the main or supplementary text. We elected to report each individual animals’ data and not to report statistics on the figures because the sample size is small.
Data availability
Data is available on OSF.
-
Open Science FrameworkID 9gdqw. Neuroanatomical abnormalities in a nonhuman primate model of congenital Zika virus infection.
References
-
Radial Glia Serve as Neuronal Progenitors in All Regions of the Central Nervous System has shown that they can be divided into multiple antigen-ically distinct subpopulations that differ in their cell cycle kinetics, with each subpopulation changing in aNeuron 41:881–890.https://doi.org/10.1016/S0896-6273(04)00140-0
-
A novel tau-based rhesus monkey model of Alzheimer’s pathogenesisAlzheimer’s & Dementia 17:933–945.https://doi.org/10.1002/alz.12318
-
The effects of neonatal amygdala or hippocampus lesions on adult social behaviorBehavioural Brain Research 322:123–137.https://doi.org/10.1016/j.bbr.2016.11.052
-
Congenital Zika infection: neurology can occur without microcephalyArchives of Disease in Childhood 104:199–200.https://doi.org/10.1136/archdischild-2018-314782
-
Animal models of human placentation--a reviewPlacenta 28 Suppl A:S41–S47.https://doi.org/10.1016/j.placenta.2006.11.002
-
BookSuomi, in Handbook of Attachment: Theory, Research, and Clinical ApplicationsGuilford Press.
-
Evo-devo and the primate isocortex: the central organizing role of intrinsic gradients of neurogenesisBrain, Behavior and Evolution 84:81–92.https://doi.org/10.1159/000365181
-
The roles of microglia in viral encephalitis: from sensome to therapeutic targetingCellular & Molecular Immunology 18:250–258.https://doi.org/10.1038/s41423-020-00620-5
-
Congenital Zika virus infection: a neuropathological reviewChild’s Nervous System 34:95–99.https://doi.org/10.1007/s00381-017-3651-3
-
A rhesus macaque model of Asian-lineage Zika virus infectionNature Communications 7:12204.https://doi.org/10.1038/ncomms12204
-
The effects of partial retinal lesions on activity and size of cells in the dorsal lateral geniculate nucleusThe Journal of Comparative Neurology 229:301–309.https://doi.org/10.1002/cne.902290214
-
Glial control of neurogenesisCurrent Opinion in Neurobiology 47:188–195.https://doi.org/10.1016/j.conb.2017.10.025
-
Ocular Histopathologic Features of Congenital Zika SyndromeJAMA Ophthalmology 135:1163–1169.https://doi.org/10.1001/jamaophthalmol.2017.3595
-
BookAffective Development: A Psychobiological PerspectiveLawrence Erlbaum Associates Publishers.
-
Segregation of geniculocortical afferents during the critical period: a role for subplate neuronsThe Journal of Neuroscience 14:3862–3880.https://doi.org/10.1523/JNEUROSCI.14-06-03862.1994
-
Signalling and functional diversity within the Axl subfamily of receptor tyrosine kinasesCytokine & Growth Factor Reviews 17:295–304.https://doi.org/10.1016/j.cytogfr.2006.04.004
-
Recent advances in neuroanatomical methodologyAnnual Review of Neuroscience 1:215–296.https://doi.org/10.1146/annurev.ne.01.030178.001243
-
Overview of General and Discriminating Markers of Differential Microglia PhenotypesFrontiers in Cellular Neuroscience 14:198.https://doi.org/10.3389/fncel.2020.00198
-
Quantitating the subtleties of microglial morphology with fractal analysisFrontiers in Cellular Neuroscience 7:3.https://doi.org/10.3389/fncel.2013.00003
-
Zika Virus Transmission Through Blood Tissue BarriersFrontiers in Microbiology 10:1465.https://doi.org/10.3389/fmicb.2019.01465
-
Development of a Zika Virus Infection Model in Cynomolgus MacaquesFrontiers in Microbiology 7:2028.https://doi.org/10.3389/fmicb.2016.02028
-
100 years of nissl stainingTrends in Neurosciences 7:236–237.https://doi.org/10.1016/S0166-2236(84)80213-1
-
What, If Anything, Is Rodent Prefrontal Cortex?ENeuro 5:ENEURO.0315-18.2018.https://doi.org/10.1523/ENEURO.0315-18.2018
-
Postmortem changes in the neuroanatomical characteristics of the primate brain: hippocampal formationThe Journal of Comparative Neurology 512:27–51.https://doi.org/10.1002/cne.21906
-
The development of ocular dominance columns in normal and visually deprived monkeysThe Journal of Comparative Neurology 191:1–51.https://doi.org/10.1002/cne.901910102
-
Microcephaly Epidemic Research Group, Microcephaly in Infants, Pernambuco StateEmerg Infect Dis 22:1090–1093.https://doi.org/10.3201/eid2206.160062
-
Factors of Concern Regarding Zika and Other Aedes aegypti-Transmitted Viruses in the United StatesJournal of Medical Entomology 54:251–257.https://doi.org/10.1093/jme/tjw212
-
Animal Models of Zika Virus Infection, Pathogenesis, and ImmunityJournal of Virology 91:e00009-17.https://doi.org/10.1128/JVI.00009-17
-
Neuronal development in human prefrontal cortex in prenatal and postnatal stagesProgress in Brain Research 85:185–222.https://doi.org/10.1016/s0079-6123(08)62681-3
-
Zika Virus Infection - After the PandemicThe New England Journal of Medicine 381:1444–1457.https://doi.org/10.1056/NEJMra1808246
-
Why primate models matterAmerican Journal of Primatology 76:801–827.https://doi.org/10.1002/ajp.22281
-
Prenatal development of the visual system in rhesus monkeyPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 278:245–260.https://doi.org/10.1098/rstb.1977.0040
-
Pre- and post-developmental neurogenesis in primatesClinical Neuroscience Research 2:29–39.https://doi.org/10.1016/S1566-2772(02)00005-1
-
The Indispensable Roles of Microglia and Astrocytes during Brain DevelopmentFrontiers in Human Neuroscience 10:566.https://doi.org/10.3389/fnhum.2016.00566
-
The INIA19 Template and NeuroMaps Atlas for Primate Brain Image Parcellation and Spatial NormalizationFrontiers in Neuroinformatics 6:00027.https://doi.org/10.3389/fninf.2012.00027
-
A cryoprotection method that facilitates cutting frozen sections of whole monkey brains for histological and histochemical processing without freezing artifactThe Journal of Histochemistry and Cytochemistry 34:1301–1315.https://doi.org/10.1177/34.10.3745909
-
BookA Combined MRI and Histology Atlas of the Rhesus Monkey Brain in Stereotaxic CoordinatesElsevier.
-
Impulse activity and the patterning of connections during CNS developmentPattern Form. Phys. Biol. Sci 5:299–332.https://doi.org/10.1201/9780429493362
-
Inner retinal vasculopathy in Zika virus diseaseAmerican Journal of Ophthalmology Case Reports 10:6–7.https://doi.org/10.1016/j.ajoc.2018.01.023
-
Can We Better Understand How Zika Leads to Microcephaly? A Systematic Review of the Effects of the Zika Virus on Human Brain OrganoidsThe Journal of Infectious Diseases 219:734–745.https://doi.org/10.1093/infdis/jiy572
-
Temperature drives Zika virus transmission: evidence from empirical and mathematical modelsProceedings. Biological Sciences 285:20180795.https://doi.org/10.1098/rspb.2018.0795
-
Interlaminar connections in the neocortexCerebral Cortex 13:5–14.https://doi.org/10.1093/cercor/13.1.5
-
Zika Virus Targeting in the Developing BrainThe Journal of Neuroscience 37:2161–2175.https://doi.org/10.1523/JNEUROSCI.3124-16.2017
-
Description of 13 Infants Born During October 2015-January 2016 With Congenital Zika Virus Infection Without Microcephaly at Birth - BrazilMMWR. Morbidity and Mortality Weekly Report 65:1343–1348.https://doi.org/10.15585/mmwr.mm6547e2
-
Visual impairment in children with congenital Zika syndromeJournal of AAPOS 21:295–299.https://doi.org/10.1016/j.jaapos.2017.04.003
-
Ocular Manifestations and Visual Outcome in Children With Congenital Zika SyndromeTopics in Magnetic Resonance Imaging 28:23–27.https://doi.org/10.1097/RMR.0000000000000192
-
Zika Virus as a Possible Risk Factor for Autism Spectrum Disorder: Neuroimmunological AspectsNeuroimmunomodulation 25:320–327.https://doi.org/10.1159/000495660
-
Modeling transformations of neurodevelopmental sequences across mammalian speciesThe Journal of Neuroscience 33:7368–7383.https://doi.org/10.1523/JNEUROSCI.5746-12.2013
-
Quantifying Microglia Morphology from Photomicrographs of Immunohistochemistry Prepared Tissue Using ImageJJournal of Visualized Experiments 1:57648.https://doi.org/10.3791/57648
-
Effects of retinal ganglion cell loss on magno-, parvo-, koniocellular pathways in the lateral geniculate nucleus and visual cortex in glaucomaProgress in Retinal and Eye Research 22:465–481.https://doi.org/10.1016/s1350-9462(03)00026-0
-
Glia in mammalian development and diseaseDevelopment 142:3805–3809.https://doi.org/10.1242/dev.129304
Decision letter
-
Rebecca ShanskyReviewing Editor; Northeastern University, United States
-
Michael A TaffeSenior Editor; University of California, San Diego, United States
-
Rebecca ShanskyReviewer; Northeastern University, United States
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Decision letter after peer review:
Thank you for submitting your article "Neuroanatomical abnormalities in a nonhuman primate model of congenital Zika virus infection" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, including Rebecca Shansky as Reviewing Editor and Reviewer #1, and the evaluation has been overseen by Michael Taffe as the Senior Editor.
The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.
The reviewers all agree that the work is novel, timely, and important. Reviewers 2 and 3 had a number of issues that need to be resolved before they consider the manuscript suitable for publication. Some involve additional immunohistochemical analyses or re-analysis of existing images (Reviewer 3). Please review these comments carefully in deciding whether to submit a revised version of this manuscript.
Summary:
This paper will be of interest to those that study the neuropathological effects of viral replication on brain development in pediatric populations and for those aiming to understand ZIKV pathogenesis. A major strength is the detailed histopathological analyses on the developing brain of non-human primates that include rigorous age-matched and procedure matched-controls.
Reviewer #1:
Children born during the Zika outbreak are known to have marked changes in neurological function and in some cases exhibit microcephaly, but the anatomical bases for these deficits is unknown. The goal of this study was to identify neuropathologies in near-term monkeys whose mothers were infected with Zika (ZIKV) during gestation. Using a combination of histological and immunofluorescent techniques, the authors demonstrate that compared to control monkeys, the ZIKV-infected fetuses exhibit a number of anatomically discrete abnormalities that are accompanied by evidence of localized neuro-immune responses and cell death. These abnormalities are shown through a comprehensive series of fluorescent micrographs that make the differences between the infected and control monkeys clear to the reader. The quantification methodologies are explained in detail, and the authors take care to note statistical considerations for the small sample size. Importantly, the observations map onto areas whose functions are affected in reports of children infected with ZIKV in utero, and therefore this work could have a major impact on public health.
No recommendations, this is a very nicely-prepared paper that needs no improvements.
Reviewer #2:
This is an excellent pathological description of neuropathological changes associated with congenital ZIKV infection in the non-human primate model. Major findings include regional variation in apoptosis, neuronal and glial morphology, and structural changes in nuclei of the visual system only seen in ZIKV challenged infants. The inclusion of three procedure control animals is exceptional and although the sample numbers are small the inclusion of these controls have uncovered clear pathological features of ZIKV infection that have not been reported convincingly prior to this manuscript. A major limitation to earlier ZIKV NHP pregnancy studies was the absence of procedure controlled uninfected animals.
This manuscript provides a detailed assessment of three different anatomic regions of the brain. The use of quantitative image analysis is a strength although the analysis could have been strengthened by evaluating increased sections from each region per animal. The small number of animals does hinder interpretation of the biological (and statistical) significance -especially since the different inoculation timepoints are grouped together which contributes to quantitative spread of the data in the image analysis.
The authors' data is supportive of a model of direct effects of viral replication on neurons in regions of the brain undergoing maturation (namely glial activation and possibly dysfunction and neuronal cell death). The rhesus macaque model is ideal to demonstrate this since regions of the macaque brain (frontal cortex) are still undergoing maturation at birth allowing for an internal control for each animal (developing versus mature brain).
There are a few stretches in interpretation of the data related to the effects of ZIKV induced immune activation and inflammation as the cause of neuronal death. Glial cells are immune cells in the brain, but there needs to be more evidence supporting the model presented that ZIKV induced inflammation and glial cell activation leads to neuronal apoptosis.
The fixation method used and sectioning method requires the use of a specialized histopathology core but the analytical methods should be able to be performed with standard image analysis tools. This paper will be useful for those doing NHP brain studies by providing a method for consistent trimming and evaluation of the brain.
The manuscript could benefit from some stylistic changes that would make the data more digestible for the non-neuropathologist.
Results:
Figure 1 and (Ln 166)-please include procedure controls for the IHC. ZIKV RNAscope for viral RNA should be included.
(Ln 213) Figure 4: 20 cells per animal, per region? This is a very small number and the findings of microstructural glial changes are very interesting-is it possible to increase the n to flesh out this data? The exact number of cells needs to be listed in the figure legend (with clarification). I think this data would benefit from showing all the cells in the graph in addition to the averages.
(Ln 234) – use of "active immune" response is inappropriate here unless there is a quantification of Iba-1 signal across groups. I do not think morphological changes can be reported as reflective of an active immune response unless this has been previously described and reported as a direct correlate to infiltrates in the brain. If so, please reference and describe better.
(Ln 237) – please indicate clearly which inoculation group had active virus in brain (all?).
Figure 7 – please include some normal for comparison here-otherwise can't interpret result without looking at supplemental figure. Would suggest including a column with just the merged from control that is in the Suppl and place next to ZIKV merged as a simple solution.
Figure 8 – see Figure 7 comments, include a high mag of normal next to high mag of Area V2 abnormal.
Figure S2 – this is exciting data and the images beautiful-can it be moved to main figures? The Iba-1 signal here is the most compelling of an "immune response" (see previous comment). Also, need to delineate the blood vessel location somehow otherwise difficult to interpret (for pathologists and non-pathologists). Maybe include H and E here? Also modify title to explain Area 46 is immature cortex. Also, please show controls at same magnification as ZIKV.
Reviewer #3:
Seelke, Beckman et al., describe and quantify neuroanatomical and histological abnormalities in the perinatal rhesus macaque brain following fetal infection with Zika virus (ZIKV) at three distinct gestational ages. The authors further show regional differences in pathological features in the cerebral cortex that are consistent with the caudal-to-rostral developmental timing of this structure and the known tropism of ZIKV for neural stem and progenitor cells. Specifically, macroanatomical defects are found in the occipital lobe/visual cortex, where the peak infection would have had more dramatic effect on cortical progenitor cell proliferation and survival, whereas phenotypes described in the frontal cortex are consistent with the presence of residual active infection, as nicely demonstrated by immunostaining for the virus envelope itself. The observation of combined cytoarchitectural abnormalities in both the visual thalamus (lateral geniculate nucleus, LGN) and primary visual cortex (area 17) is intriguing, suggestive of how disruptions in the retina or other more-sensitive structures/cell types might propagate up a circuit. In addition, the observations of ZIKV+ and apoptotic cells in the frontal cortex as many as 15 weeks after acute infection constitute remarkable and valuable data. Lastly, the authors provide evidence of the response of glial cell types in the brain to ZIKV, by quantitative morphological analyses of microglia and astrocytes.
Overall the authors' observations in these valuable tissues, although limited in statistical power by the small n, provide new data on the effects of ZIKV infection on brain development. However, the study is descriptive and aspects of the study, particularly with respect to affected cell types and glial cell responses, would benefit from additional cell type-specific stainings, improved imaging and presentation, and additional analysis in some places. The biggest weakness is the small sample size. Although this is a challenge with NHP studies in general, it nevertheless a weakness as reproducibility and animal to animal variability is a significant concern.
1) Additional immunofluorescent analysis of the frontal cortex in conjunction with the ZIKV staining shown in Figure 1 and caspase staining shown in Figure 4, to identify the cell types that show persistent infection many weeks after inoculation and clarify which cell types are most likely to be undergoing apoptosis at this late stage of cortical development, would further support the central argument that areal differences in pathology across the cortex are related to the caudorostral temporal developmental gradient. If the thesis is simply that ZIKV persists longer in frontal cortex because neurogenesis does as well, this could be definitively supported by such co-immunostaining.
2) Morphological changes observed in microglia of the frontal cortex are consistent with these immune cells' response to the active infection and cell death shown by the presence of ZIKV+ and caspase+ cells in this brain area. However, such morphology is quite nonspecific and the present data show correlation only, rather than suggesting any mediating, causal role for glial cell types in either the persistence or the pathological effects of Zika infection in this brain region. Without more detailed analysis of microglia state changes in the presence of Zika, it is impossible to infer whether these cells are performing their expected beneficial functions in responding to infection and tissue damage, or are additionally contributing to pathogenesis through aberrant activation or inflammatory responses.
3) The central assertion that areal differences in cortical phenotypes following Zika infection earlier in gestation are related to the temporal caudal-to-rostral gradient of cortical development is intriguing, but the manuscript would benefit from more detailed discussion and a clear model of this hypothesis, how it would arise from changes in the neural progenitor cell type composition and behaviors over cortical development, and how such interactions might be expected to result in differential structural and functional deficits in individuals infected at earlier vs later stages of fetal development.
In general this small (n=3 per group, if we ignore the different gestational ages of infection) pilot study provides additional evidence that fetal ZIKV infection disrupts brain development in ways that are influenced by the normal gradient of regional and cortical areal neurogenesis -- e.g. caudal-to-rostral in the cortex. Combined cytoarchitectural disruptions of the LGN, V1, and occipital gyrification are intriguing, suggesting that disruptions in the retina or other more-sensitive structures/cell types might propagate up a circuit, although no mechanistic interpretation is provided and further discussion of how such propagation might occur would be nice. In addition, the observation of ZIKV+ and apoptotic cells in the frontal cortex so long after the acute infection is a remarkable and valuable piece of data. The study would be greatly strengthened by additional immunohistochemical analyses of the frontal cortex to determine the cellular identity of ZIKV+ and apoptotic cells, and/or the spatial distribution of ZIKV+ cells relative to apoptotic neurons or other cell types. The authors' assertion that their methods reveal "where the virus was found, [and] what cell types were affected" is not fully supported by the data provided, given that ZIKV+ cells are not co-stained for any markers and their location/distribution is not described beyond the fact that they are present in both grey and white matter.
The glial angle of the work is lacking in strong phenotypes or clear mechanistic insight. Microglia morphological changes in the frontal cortex are a clear but non-specific observation, and not terribly surprising in proximity to active virus and/or apoptosis (although how proximal, again, is not clear without additional images or co-stainings). Although a relationship of the microglial phenotype to the increased cell death is implied, no mechanistic hypothesis is explicitly stated, and it is just as likely that microglia would become activated in response to cell death rather than cause it. Macroscopic phenotypes as observed in occipital cortex and LGN are even further removed from any potential glial mechanisms and in fact are quite similar to phenotypes described to arise from various non-viral/genetic/non-inflammatory disruptions of NPC survival, proliferation, and differentiation.
In general, the authors suggest or imply mechanisms which could stand to be more explicitly stated as conjectures or hypotheses supported by the current data and/or subject to future study. In addition to the specific role of microglia, the impact of the temporal caudal-rostral developmental gradient on the final neuropathology following infection at distinct gestational ages is not clearly explained or discussed. One gets the sense that reduced gyrification in occipital but not frontal cortex, or absence of changes to cortical thickness, could be related to the earlier development of occipital cortex, but exactly how this would occur -- e.g. through differential effects on subtypes of neural stem/progenitor cells and their relationships to tangential vs radial expansion of cortical plate. An explicit model of how regional phenotypic differences in cortex could arise given infection at different fetal stages would be a nice addition to the manuscript. (It's not clear why gyrification should not also be reduced in frontal lobe if infection was present at the earlier gestational ages when ventricular progenitors are highly proliferative -- possibly because the frontal lobe has not completed folding at birth in macaques?)
Even with further and clearer discussion of interpretations and model, I do not believe the paper merits publication in eLife without additional experiments, primarily co-immunostaining of ZIKV and CC3 with each other and with NPC and neuronal cell type markers.
Specific comments and suggestions to strengthen study (using existing tissue/samples).
Figure 1. Main finding: ZIKV is found in frontal but not occipital ctx nor in subcortical areas at near-term, after inoculation weeks earlier.
– What are the infected cell types? Sox2/NeuN/ZIKV co-stain would go far in supporting the argument that different phenotypes observed between frontal and occipital lobes relates to the late development of frontal lobe and persistence of NPC in this area, allowing ZIKV to continue to proliferate.
– A lower magnification view, or some quantifications, illustrating the frequency/density/laminar distribution of ZIKV+ cells in the frontal lobe would be very helpful in inferring the potential downstream impacts of persistent infection in this brain region.
– The term "activated" in the Figure 1 legend and associated Results section to describe astrocytes and microglia should be defined here, or simply removed given later discussion at Figure 4. In fact, text and figure legend here refer to "activated astrocytes" and "reactive astrocytes," but the analysis shown in Figure 4 revealed no significant differences in frontal lobe astrocytes between ZIKV+ and control, a discrepancy that should be addressed.
Figure 2. Main finding: ZIKV reduced gyrification in occipital but not frontal lobe, and caused no significant changes in total, GM, or WM volumes.
–To my knowledge a more commonly reported (though not necessarily better) local gyrification measurement is the ratio of outer cortical surface area to the "convex hull" area.
- the curve that encloses the exposed cortex, excluding sulci. See for example https://pubmed.ncbi.nlm.nih.gov/20176115/
Since this is a more commonly reported GI and quite straightforward/quick to implement with the existing images it would be nice to report this measurement in addition to what's already included.
– Given the abstract and introduction's emphasis on the rostrocaudal developmental timing of the cortex, the authors should further discuss why they think such timing might result in reduced gyrification in occipital lobe, but not frontal.
Figure 3. Main finding: no changes in cortical thickness. Again, there could perhaps be some discussion of the motivation for taking this measurement, or placing this finding in the context of the larger argument about developmental timing.
Figure 4. Main finding: frontal lobe microglia in ZIKV+ brains show smaller whole cell volume and fewer terminal branch points, suggestive of an activated state.
– Following on the comment to Figure 1 above, the main text and figure 4 legend need to be aligned - e.g. the figure legend states "astrocytes tended to have smaller cell body volumes (R)" but the main text reports this comparison as having a p-value of 0.37 and states astrocytes "did not differ between ZIKV-infected and control animals."
– If there is no statistically significant difference in astrocyte morphology, density, or distribution in ZIKV+ cortex, the terms "activated" or "reactive" astrocytes should be removed.
Figure 5. Main finding: increased number of caspase-positive neurons in cortical plate of area 46 of ZIKV+ brains.
– The following sentence/section should be rewritten to clarify the hypothesis being tested: "IHC fluorescent microscopy analyses for cleaved caspase-3 and SATB2, which are markers of apoptosis onset and immature cortical neurons, respectively, were carried out in order to determine whether NPCs underwent higher rates of apoptosis in immature neurons induced by the presence of ZIKV in the brain in both Area 17 and Area 46." Satb2 is a neuronal marker whereas previously apoptosis was observed in NPCs (progenitors, e.g. Sox2+ cells). ZIKV was present in area 46 but not 17 at the time of analysis.
– Are any of the caspase-positive SATB2+ neurons also ZIKV+? Are there any SATB2-negative caspase+ cells? Are there more total CC3+ cells in ZIKV+ vs control? Please be more specific in the units of measurement shown on the y-axis -- number of cells per what area/volume? It would also be helpful to report these data as total density of caspase+ cells regardless of SATB2 immunoreactivity, as well as the proportions of SATB2+ cells that are caspase+, and vice versa. Is there a preferential increase in apoptosis of neurons, or a general increase in apoptosis affecting multiple cell types?
– Sox2+ glial progenitor cells derived from earlier VZ/SVZ NPC are abundant across the cortical plate at late gestation. Given the previous reports of NPC infection by ZIKV, did the authors also observe caspase+ Sox2+ cells? (Again, what is the identity of the ZIKV+ cells shown in Figure 1?)
– If the apoptotic neurons are not themselves ZIKV+, what is the proposed mechanism/cause of their apoptosis?
Figure 6. Main finding: blurred laminar boundaries/absence of koniocellular layers in ZIKV+ LGN.
– Please include an un-pseudocolored version of the control section to allow the reader to better see what the normal layering pattern should look like across the whole LGN.
Figure 7. Main finding: reduced neuronal density, increased "glia activation" in dysmorphic areas of LGN.
– The phenotype asserted in the main text, "the density of the neuronal population (stained with NeuN) in these areas was also substantially reduced," is not immediately obvious from the images. It would help to indicate more precisely on the Nissl image the area shown at higher mag in IFL panels. This should be roughly achievable by comparing the Nissl and DAPI on the adjacent sections. Better yet would be to show a box on panels G-K to indicate the areas shown in B-F and L-P.
– Most convincing would be some quantification of NeuN+ cell density or cumulative signal in the dysmorphic ROI's, compared either to neighboring unaffected regions of LGN or to matched ROI's from the control LGNs.
– The assertion of "glia activation" is far from obvious given only the images in Figure 7. Quantitative morphological analysis of microglia similar to what was performed for cortex (Figure 4) would be a big help. Just improving the clarity of the images in panels C, D, M, and N would also be helpful; it is difficult to make out for example the increase in GFAP content in the NeuN-poor region of panel N.
https://doi.org/10.7554/eLife.64734.sa1Author response
Reviewer #1:
Children born during the Zika outbreak are known to have marked changes in neurological function and in some cases exhibit microcephaly, but the anatomical bases for these deficits is unknown. The goal of this study was to identify neuropathologies in near-term monkeys whose mothers were infected with Zika (ZIKV) during gestation. Using a combination of histological and immunofluorescent techniques, the authors demonstrate that compared to control monkeys, the ZIKV-infected fetuses exhibit a number of anatomically discrete abnormalities that are accompanied by evidence of localized neuro-immune responses and cell death. These abnormalities are shown through a comprehensive series of fluorescent micrographs that make the differences between the infected and control monkeys clear to the reader. The quantification methodologies are explained in detail, and the authors take care to note statistical considerations for the small sample size. Importantly, the observations map onto areas whose functions are affected in reports of children infected with ZIKV in utero, and therefore this work could have a major impact on public health.
No recommendations, this is a very nicely-prepared paper that needs no improvements.
Thank you so much for the positive feedback about our manuscript.
Reviewer #2:
[…]
The manuscript could benefit from some stylistic changes that would make the data more digestible for the non-neuropathologist.
Results:
Figure 1 and (Ln 166)-please include procedure controls for the IHC.
We thank the reviewer for the commentaries and suggestions to improve our manuscript. To better demonstrate that ZIKV induced glia abnormalities in Area 46 in comparison with control monkeys, we completely modified Figure 1. The first panel of Figure 1 shows now representative micrographs of the same region in Area 46 from a control and a ZIKV-infected animal. A second and a third panel were also added in figure 1, highlighting that Iba1+ microglia involved in the neuroinflammatory abnormalities observed in the ZIKV-infected monkeys, also expresses CD68, a lysosomal protein highly expressed in activated microglia. A new section titled: ‘Neuroinflammatory abnormalities and ZIKV presence in brains from ZIKV-infected animals”, was also added in the “Results” section, starting on Line 161.
It reads:
“Immunohistochemical (IHC) confocal microscopy revealed the presence of clusters of reactive glia within Area 46 of the PFC from ZIKV+ animals, but not in the controls (Figure 1A-B). Multiple labelling fluorescent microscopy allowed us to visualize that abnormal clusters of Iba1+ microglia which also express the lysosomal marker CD68 in the ZIKV-infected animals (Figure 1C-G). CD68 levels are known to be substantially upregulated in microglia during inflammatory processes and are thought to be involved in active phagocytosis (46). We also observed that clusters of amoeboid-shaped microglia within Area 46, are directly engulfing neuronal debris in the ZIKV-infected monkeys (Figure 1H-I, 3D reconstruction: J, K).
Next, we sought to detect the presence of ZIKV protein in the same regions we observed active inflammatory response. Interestingly, we observed the presence of ZIKV envelope protein within Area 46 of all the infected animals, with detection of viral protein inside clusters of microglia, but not astrocytes (GFAP, Figure 2A-G, 3D reconstruction: 2H). Despite the fact that cells expressing ZIKV were detected in Area 46, we were not able to confirm the identity of the infected cells. As shown in Figure 2I-J, these ZIKV+ cells are thought to be in a direct process of phagocytosis, a common event occurring during viral encephalitis (47). While we could not detect the expression of other proteins in the ZIKV+ cells, they were surrounded and in direct contact with reactive microglia also expressing CD68, as shown in Figure 2K-L.”
ZIKV RNAscope for viral RNA should be included.
Unfortunately, due to the protocol used for tissue fixation and preparation, we were unable to detect viable RNA for performing ZIKV RNAscope investigation. While there are groups who have used RNAscope in fixed tissue, many groups have been unsuccessful in processing fixed tissue in this way. Given this, we are electing to use evidence from immunohistochemistry demonstrating presence of Zika antibodies.
(Ln 213) Figure 4: 20 cells per animal, per region? This is a very small number and the findings of microstructural glial changes are very interesting-is it possible to increase the n to flesh out this data? The exact number of cells needs to be listed in the figure legend (with clarification). I think this data would benefit from showing all the cells in the graph in addition to the averages.
All the images for glia morphological analysis were acquired and processed using 3D high-resolution volumetric analysis. The number of 20 cells per animal per region was described recently by our group to be enough to detect subtle changes in glia morphology induced by inflammatory events (Beckman et al., PNAS 2019, http://doi.org/10.1073/pnas.1902301116, Beckman et al., Alzheimer’s and Dementia, 2021, https://doi.org/10.1002/alz.12318). Due to the fact that the mentioned studies were performed using middle-aged rhesus monkey brain, instead of the perinatal brains analyzed in the current study, we increased the number of glial cells analyzed to 28 cells per animal, per region. The exact number of cells analyzed and details about image processing were added in the legend of the figure. While we appreciate that the reviewer would prefer to see all of the cells this would create a huge figure (28 x 2 (astrocytes, microglia) X 2 brain areas (area 17 and 46) X 6 animals) and we think that the main story communicated by the example cells and statistical figure would not be clear from that figure. The graphs from all other figures only present a summary of the data collected from each animal, instead of all cells or micrographs analyzed. That way, we believe keeping the graph with this same visual identity in figure 4 will facilitate for the reader, comparisons between groups and figures. In addition, we included, as suggested by reviewer 3, a new image, more representative, for reactive microglia within area 46 of the ZIKV infected animals.
(Ln 234) – use of "active immune" response is inappropriate here unless there is a quantification of Iba-1 signal across groups. I do not think morphological changes can be reported as reflective of an active immune response unless this has been previously described and reported as a direct correlate to infiltrates in the brain. If so, please reference and describe better.
We agree with the reviewer that the term “active immune response” might induce the establishment of a connection between our observation and direct findings of infiltrates in the brain, which we cannot conclude based on the current data. We modified the beginning of sentence for “A neuroinflammatory response, mainly driven by CD68/Iba+ microglia” (Line 252). As mentioned in the previous reply, we and others have shown that alterations in microglia morphology can indicate a direct response of these cells to inflammatory events in the brain. We believe that by combining now the general marker Iba-1 with the marker for activation CD68, we were able to show more clearly the direct involvement of activated microglia in phagocytizing ZIKV protein in the frontal lobe.
That section that begins on Line 252 now reads:
“Microstructural changes related to neural development. A neuroinflammatory response, mainly driven by CD68/Iba+ microglia, in combination with the presence of ZIKV envelope protein, indicated that the frontal lobe for all ZIKV-infected subjects was a site of persistent ZIKV infection for a minimum of 60 days post-inoculation. […] The higher frequency of immature neuron death, the active immune response (activated microglia), and the presence of ZIKV protein in Area 46 suggest that frontal cortex remained a site of active ZIKV-induced neuronal remodeling at the time the brains were analyzed.”
(Ln 237) – please indicate clearly which inoculation group had active virus in brain (all?).
All three Zika infected animals analyzed in the study had evidence of active virus in their brains. We have now indicated that on Line 172 and 313 of the manuscript.
Figure 7 – please include some normal for comparison here-otherwise can't interpret result without looking at supplemental figure. Would suggest including a column with just the merged from control that is in the Suppl and place next to ZIKV merged as a simple solution.
We agree with the reviewer that showing side-by-side normal vs. zika infected images is important for understanding and interpreting the results. In the mentioned figure (previously Figure 7, now Figure 8), we have side-by-side representative images of Nissl and confocal images of LGN from both CTR and ZIKV groups.
Figure 8 – see Figure 7 comments, include a high mag of normal next to high mag of Area V2 abnormal.
High magnification representative images of normal V1/V2 transition vs. layer discontinuity and V2 abnormality, are now included side-by-side in the new figure (current Figure 9).
Figure S2 – this is exciting data and the images beautiful-can it be moved to main figures? The Iba-1 signal here is the most compelling of an "immune response" (see previous comment). Also, need to delineate the blood vessel location somehow otherwise difficult to interpret (for pathologists and non-pathologists). Maybe include H and E here? Also modify title to explain Area 46 is immature cortex. Also, please show controls at same magnification as ZIKV.
We thank the reviewer for the suggestion of moving Figure S2 (ZIKV neuroinflammatory response in area 46) to the main manuscript. As mentioned in the reply to comment 1A, we completely modified main figures 1 and 2 in the manuscript to highlight the details of how ZIKV infection is inducing a neuroinflammatory response driven by activated microglia (Iba+/CD68+), that is not present in the control animals. The title explaining Area 46 is immature cortex was also added in the legend of Figure 1.
Reviewer #3:
[…]
In general this small (n=3 per group, if we ignore the different gestational ages of infection) pilot study provides additional evidence that fetal ZIKV infection disrupts brain development in ways that are influenced by the normal gradient of regional and cortical areal neurogenesis -- e.g. caudal-to-rostral in the cortex. Combined cytoarchitectural disruptions of the LGN, V1, and occipital gyrification are intriguing, suggesting that disruptions in the retina or other more-sensitive structures/cell types might propagate up a circuit, although no mechanistic interpretation is provided and further discussion of how such propagation might occur would be nice. In addition, the observation of ZIKV+ and apoptotic cells in the frontal cortex so long after the acute infection is a remarkable and valuable piece of data. The study would be greatly strengthened by additional immunohistochemical analyses of the frontal cortex to determine the cellular identity of ZIKV+ and apoptotic cells, and/or the spatial distribution of ZIKV+ cells relative to apoptotic neurons or other cell types. The authors' assertion that their methods reveal "where the virus was found, [and] what cell types were affected" is not fully supported by the data provided, given that ZIKV+ cells are not co-stained for any markers and their location/distribution is not described beyond the fact that they are present in both grey and white matter.
We appreciate the comment from the reviewer and understand how useful it would be to know what cell types are particularly subject to apoptosis during fetal Zika virus infection. We did carry out some of these analyses, attempting to co-localize markers of apoptosis with other markers – for example our CC3 and SATB2 co-localization analyses give us information about the death of immature neurons. Unfortunately, in many cases, once apoptosis begins and certainly once it is underway, identification of the cells using IHC (via surface markers) does not paint a comprehensive picture of the system because the cells literally fall apart and so we were limited in terms of what approaches were successful. We agree that this is a potentially fruitful avenue for future research where tissue is processed at a different time point of infection (i.e., our time point is ~100 days post infection).
The glial angle of the work is lacking in strong phenotypes or clear mechanistic insight. Microglia morphological changes in the frontal cortex are a clear but non-specific observation, and not terribly surprising in proximity to active virus and/or apoptosis (although how proximal, again, is not clear without additional images or co-stainings). Although a relationship of the microglial phenotype to the increased cell death is implied, no mechanistic hypothesis is explicitly stated, and it is just as likely that microglia would become activated in response to cell death rather than cause it. Macroscopic phenotypes as observed in occipital cortex and LGN are even further removed from any potential glial mechanisms and in fact are quite similar to phenotypes described to arise from various non-viral/genetic/non-inflammatory disruptions of NPC survival, proliferation, and differentiation.
It is the case that our data do not allow us to determine whether inflammatory processes are causing or are caused by cell death, and we have now indicated this in the Discussion (Lines 398-400). It reads:
“The histological analyses presented here differ in that we measured specific cytoarchitectonic regions of the brain as well as morphological analyses of glial cells in the frontal and occipital cortex. Together, these analyses paint a picture of inflammatory processes related to viral infiltration sweeping through the cortex and a wave of cell death. Given our analyses, it is not possible to tell whether inflammatory processes are caused by or are causing cell death, but they certainly appear to co-occur given their proximity.”
We have also added a discussion of the relationship between morphology and activation in the Results section, beginning at Line 216. It reads:
“Microstructural Changes to Glia. Following macrostructural evaluations, we carried out a series of immunohistochemical analyses to quantify cell-level features that might be impacted by ZIKV infection with a specific focus on glia (microglia and astrocytes) in Brodmann’s Area 17 and Area 46. […] Iba1 and GFAP were used as general microglia and astrocytes markers, respectively, and each cell was exported and analyzed individually for total and cell body volumes, and for number of terminal point (branching ramification).”
In general, the authors suggest or imply mechanisms which could stand to be more explicitly stated as conjectures or hypotheses supported by the current data and/or subject to future study. In addition to the specific role of microglia, the impact of the temporal caudal-rostral developmental gradient on the final neuropathology following infection at distinct gestational ages is not clearly explained or discussed. One gets the sense that reduced gyrification in occipital but not frontal cortex, or absence of changes to cortical thickness, could be related to the earlier development of occipital cortex, but exactly how this would occur -- e.g. through differential effects on subtypes of neural stem/progenitor cells and their relationships to tangential vs radial expansion of cortical plate. An explicit model of how regional phenotypic differences in cortex could arise given infection at different fetal stages would be a nice addition to the manuscript. (It's not clear why gyrification should not also be reduced in frontal lobe if infection was present at the earlier gestational ages when ventricular progenitors are highly proliferative -- possibly because the frontal lobe has not completed folding at birth in macaques?)
Even with further and clearer discussion of interpretations and model, I do not believe the paper merits publication in eLife without additional experiments, primarily co-immunostaining of ZIKV and CC3 with each other and with NPC and neuronal cell type markers.
We appreciate the comments from the reviewer, and in order to address these questions we ran a new set of IHC experiments with multi-labeling fluorescence combinations. A completely modified Figure 1 and a whole new Figure 2 were added. The new micrographs included in Figure 1 highlight: 1 – The general glia alterations observed in the prefrontal cortex of ZIKV infected animals but not in control animals (Figure 1A-B); 2 – Inflammatory clusters observed in the ZIKV tissue were formed by reactive astrocytes and microglia, with the latest also co-expressing CD68, a lysosomal protein expressed in phagocytic cells (Figure 1C-G); 3 – 3D analysis of the reactive microglia from the ZIKV animals show these cells are actively engulfing apoptotic material formed by fragmented Dapi and NeuN (Figure 1H-K).
To address the identity of the ZIKV+ cells, we ran several multi-label antibodies combinations targeting glia, neurons, and apoptotic cells. Figure 2 panels now show that reactive microglia and astrocytes are surrounding ZIKV+ cells, and that ZIKV protein is observed inside Iba+/CD68+ microglia. Importantly, panel I-J shows in detail the 3D reconstruction of how microglia is interacting with ZIKV+ cells. After several rounds of IHC with different multi-labeling combinations, we cannot clearly point to an identity for the main ZIKV+ cells within these regions. We believe these cells are in an advanced cell death/phagocytosis stage, when protein synthesis is impaired inside the cell, and as consequence, labeling them becomes a real challenge, especially with limitations in the tissue fixation after necropsy. In addition, we do not observe a clear pattern of location/distribution of these infected cells within Area 46, besides the fact that they are observed in both, white and gray matter. Clear imaging of these cells can only be done using high-resolution microscopy, as often these cells are not expressing other markers, including DAPI. Imaging the entire Area 46 of the monkeys used in this study using the 63 or 100x objective to quantify the distribution and location of these cells would generate hundreds of images and an extensive amount of time to analyze it, that is beyond the scope of this project. But we expect that the new data added in the figures provides more evidence of the persistent of the virus within specific areas of the brain, and how this can reflect in neuroinflammatory abnormalities and neuronal loss, even after longer periods of viral incubation.
Specific comments and suggestions to strengthen study (using existing tissue/samples).
Figure 1. Main finding: ZIKV is found in frontal but not occipital ctx nor in subcortical areas at near-term, after inoculation weeks earlier.
– What are the infected cell types? Sox2/NeuN/ZIKV co-stain would go far in supporting the argument that different phenotypes observed between frontal and occipital lobes relates to the late development of frontal lobe and persistence of NPC in this area, allowing ZIKV to continue to proliferate.
– A lower magnification view, or some quantifications, illustrating the frequency/density/laminar distribution of ZIKV+ cells in the frontal lobe would be very helpful in inferring the potential downstream impacts of persistent infection in this brain region.
– The term "activated" in the Figure 1 legend and associated Results section to describe astrocytes and microglia should be defined here, or simply removed given later discussion at Figure 4. In fact, text and figure legend here refer to "activated astrocytes" and "reactive astrocytes," but the analysis shown in Figure 4 revealed no significant differences in frontal lobe astrocytes between ZIKV+ and control, a discrepancy that should be addressed.
Related to the cell type of the ZIKV infected cells, as we mentioned in the previous reply, we ran several IHC experiments with different antibodies combination, including Sox2/NeuN/ZIKV, as suggested by the reviewer. Fragments of both neuronal markers were found inside phagocytic microglia in direct contact with ZIKV+ cells, but not inside the infected cells the ZIKV+ cell (Figures 1 and 2).
We agree with the reviewer that understanding the frequency/density/laminar distribution of ZIKV+ cells in the frontal lobe would provide deeper understanding of the viral impact in the normal cortical development. Our intention is to do a detailed stereological quantification of the Zika viral distribution in the frontal lobe of the rhesus monkeys, but that is beyond the scope of this paper, this paper is intended to be an initial description of the pathology observed in a recent developed monkey model. We believe that for this type of analysis to be done correctly, substantial amount of time and work will be required, especially due to the size and complexity of this cerebral region in the monkey brain.
In relation to the use of the term “activated” in Figure 1 legend, we agree that there is no indication of astrocyte activation, and the term should not be used in relation to the astrocyte’s morphology. All mentions to the term “activated” in relation to astrocytes were removed from the text.
Figure 2. Main finding: ZIKV reduced gyrification in occipital but not frontal lobe, and caused no significant changes in total, GM, or WM volumes.
–To my knowledge a more commonly reported (though not necessarily better) local gyrification measurement is the ratio of outer cortical surface area to the "convex hull" area.
- the curve that encloses the exposed cortex, excluding sulci. See for example https://pubmed.ncbi.nlm.nih.gov/20176115/
Since this is a more commonly reported GI and quite straightforward/quick to implement with the existing images it would be nice to report this measurement in addition to what's already included.
– Given the abstract and introduction's emphasis on the rostrocaudal developmental timing of the cortex, the authors should further discuss why they think such timing might result in reduced gyrification in occipital lobe, but not frontal.
We thank the reviewer for the suggestion. We analyzed the local gyrification index in the frontal and occipital lobes using the method described by Rogers at al 2010 (the paper mentioned by the reviewer). We analyzed between 14 and 19 sections for each animal in each group but found no difference in the cerebral gyrification index between the groups for both regions of the cortex. As such, we’ve maintained our original presentation of the gyrification data.
Figure 3. Main finding: no changes in cortical thickness. Again, there could perhaps be some discussion of the motivation for taking this measurement, or placing this finding in the context of the larger argument about developmental timing.
Changes of cortical thickness during development might reflect spines/dendrites loss and ultimately neuronal death. In microencephaly induced by ZIKV infection, the loss of cortical white and gray matter can be so substantial that cortical thickness measurement can be used to indicative the degree of the pathology. The investigation of cortical thickness in the monkeys using Nissl was performed in combination with the fluorescent microscopy experiments investigating increased cell death of immature neurons in the frontal lobe of ZIKV animals. We believe is important to show that despite the fact that there is no microencephaly cases between the infected monkeys, nor a change in the cortical thickness, there are significant micro-structural changes such as increased death of immature neurons and robust neuroinflammatory response, that ultimately can contribute to long-term, chronic impairment related to early-development ZIKV infection.
Figure 4. Main finding: frontal lobe microglia in ZIKV+ brains show smaller whole cell volume and fewer terminal branch points, suggestive of an activated state.
– Following on the comment to Figure 1 above, the main text and figure 4 legend need to be aligned - e.g. the figure legend states "astrocytes tended to have smaller cell body volumes (R)" but the main text reports this comparison as having a p-value of 0.37 and states astrocytes "did not differ between ZIKV-infected and control animals."
– If there is no statistically significant difference in astrocyte morphology, density, or distribution in ZIKV+ cortex, the terms "activated" or "reactive" astrocytes should be removed.
We agree with the reviewer, and as mentioned before, all mentions to “activation” in relation to astrocytes, were removed from the main text and figure legends.
Figure 5. Main finding: increased number of caspase-positive neurons in cortical plate of area 46 of ZIKV+ brains.
– The following sentence/section should be rewritten to clarify the hypothesis being tested: "IHC fluorescent microscopy analyses for cleaved caspase-3 and SATB2, which are markers of apoptosis onset and immature cortical neurons, respectively, were carried out in order to determine whether NPCs underwent higher rates of apoptosis in immature neurons induced by the presence of ZIKV in the brain in both Area 17 and Area 46." Satb2 is a neuronal marker whereas previously apoptosis was observed in NPCs (progenitors, e.g. Sox2+ cells). ZIKV was present in area 46 but not 17 at the time of analysis.
– Are any of the caspase-positive SATB2+ neurons also ZIKV+? Are there any SATB2-negative caspase+ cells? Are there more total CC3+ cells in ZIKV+ vs control? Please be more specific in the units of measurement shown on the y-axis -- number of cells per what area/volume? It would also be helpful to report these data as total density of caspase+ cells regardless of SATB2 immunoreactivity, as well as the proportions of SATB2+ cells that are caspase+, and vice versa. Is there a preferential increase in apoptosis of neurons, or a general increase in apoptosis affecting multiple cell types?
– Sox2+ glial progenitor cells derived from earlier VZ/SVZ NPC are abundant across the cortical plate at late gestation. Given the previous reports of NPC infection by ZIKV, did the authors also observe caspase+ Sox2+ cells? (Again, what is the identity of the ZIKV+ cells shown in Figure 1?)
– If the apoptotic neurons are not themselves ZIKV+, what is the proposed mechanism/cause of their apoptosis?
We agree with the reviewer, and the sentence mentioned was rewritten to address that IHC fluorescent microscopy for cleaved caspase 3 was conducted in order to investigate if the increase in reactive microglia population and active presence of ZIKV can induce an increase in the number of immature neurons (SATB2) undergoing apoptosis (Ln 271-276). We carried out additional IHC combinations in order address the reviewer questions, and an updated version of the figure mentioned (currently Figure 6) and new Supplementary Figure 3 were added to the manuscript. The updated Figure 6 includes now the total cell counting for each of the markers analyzed, DAPI (total number of cells), SATB2 (total number of immature neurons) and CC3 (total number of apoptotic cells). A separate graph analyzing SATB2/CC3 colocalization was also added in the figure, for each region analyzed. Despite the fact that there is no significant difference in the total number of CC3 cells across the two analyzed regions between control and ZIKV animals, an increase in SATB2/CC3 was observed in the ZIKV animals in Area 46, but not 17. A deeper investigation of these regions with higher-resolution image acquisition showed that activated microglia (HLA-Dr+ cells) also express CC3 in the region, especially the cells in direct contact with neurons, possibly during a phagocytosis process. Images for this analysis are show in the new supplementary figure 3. Finally, once again we carry out different experiments to address if ZIKV+ is present in these SATB2/CC3 or HLA-Dr/CC3 cells, but once again, the virus could not be detected in these cells. Whether this is caused by loss of viral epitopes during the brain fixation process, or due to the stage of the viral infection/engulfment by microglia per se, could not be address at the time for this manuscript.
Figure 6. Main finding: blurred laminar boundaries/absence of koniocellular layers in ZIKV+ LGN.
– Please include an un-pseudocolored version of the control section to allow the reader to better see what the normal layering pattern should look like across the whole LGN.
Figure 7. Main finding: reduced neuronal density, increased "glia activation" in dysmorphic areas of LGN.
We thank the reviewer for the suggestion. To better allow the reader to compare between normal and altered layers of the LGN, we removed the coloring patterns from the micrographs. We concluded that the coloring pattern applied before was not necessarily and removed from the figure.
– The phenotype asserted in the main text, "the density of the neuronal population (stained with NeuN) in these areas was also substantially reduced," is not immediately obvious from the images. It would help to indicate more precisely on the Nissl image the area shown at higher mag in IFL panels. This should be roughly achievable by comparing the Nissl and DAPI on the adjacent sections. Better yet would be to show a box on panels G-K to indicate the areas shown in B-F and L-P.
– Most convincing would be some quantification of NeuN+ cell density or cumulative signal in the dysmorphic ROI's, compared either to neighboring unaffected regions of LGN or to matched ROI's from the control LGNs.
– The assertion of "glia activation" is far from obvious given only the images in Figure 7. Quantitative morphological analysis of microglia similar to what was performed for cortex (Figure 4) would be a big help. Just improving the clarity of the images in panels C, D, M, and N would also be helpful; it is difficult to make out for example the increase in GFAP content in the NeuN-poor region of panel N.
We agree with the reviewer, and the mentioned Figure (currently Figure 8) was redone to better show how the Nissl images correlate with the confocal images between controls and ZIKV animals. A “zoom” area was also added for both groups to highlight the differences observed between the groups. We agree with the reviewer that quantification of NeuN and Iba1 count would be informative, but require a substantial amount of time, and are beyond the scope of this manuscript now. Finally, we also agree that the assertion “glia activation” is not obvious, so besides improving the clarity of the images in the figure, we also removed the term “glia activation” from the figure legend and related text in the Result sections.
https://doi.org/10.7554/eLife.64734.sa2Article and author information
Author details
Funding
National Institute of Neurological Disorders and Stroke (R21NS104692)
- Eliza Bliss-Moreau
National Institute of Allergy and Infectious Diseases (R21AI129479)
- Lark LA Coffey
California National Primate Research Center (Pilot Grant P51OD011107)
- Lark LA Coffey
Eunice Kennedy Shriver National Institute of Child Health and Human Development (R01HD096436)
- Eliza Bliss-Moreau
National Institute of Allergy and Infectious Diseases (R21 AI129479)
- Koen KA Van Rompay
University of California, Davis (Murray B. Gardner Jr. Fellowship)
- Eliza Bliss-Moreau
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
Thank you to Jennifer Watanabe, Jodie Usachenko, Christina Cruzen, Kari Christe, Rachel Reader, Wilhelm von Morgenland, Anil Singapuri, Gilda Moadab, and members of the Bliss-Moreau Laboratory who contributed to this work. R21NS104692, R01HD096436, and the Murray B Gardner Junior Faculty Research Fellowship from the UC Davis Center of Comparative Medicine to EBM; R21AI129479 and a pilot grant from the California National Primate Research Center via P51OD011107 to LLC.
Ethics
All procedures were approved by the University of California, Davis Institutional Animal Care and Use Committee (Protocol # 19211) which is accredited by the Association for Assessment and Accreditation of Laboratory Animal Care International (AAALAC). Animal care was performed in compliance with the 2011 Guide for the Care and Use of Laboratory Animals provided by the Institute for Laboratory Animal Research.
Senior Editor
- Michael A Taffe, University of California, San Diego, United States
Reviewing Editor
- Rebecca Shansky, Northeastern University, United States
Reviewer
- Rebecca Shansky, Northeastern University, United States
Version history
- Received: November 9, 2020
- Accepted: December 14, 2021
- Version of Record published: March 9, 2022 (version 1)
- Version of Record updated: April 26, 2022 (version 2)
Copyright
© 2022, Beckman et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 785
- Page views
-
- 115
- Downloads
-
- 5
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Microbiology and Infectious Disease
- Plant Biology
Purinergic signaling activated by extracellular nucleotides and their derivative nucleosides trigger sophisticated signaling networks. The outcome of these pathways determine the capacity of the organism to survive under challenging conditions. Both extracellular ATP (eATP) and Adenosine (eAdo) act as primary messengers in mammals, essential for immunosuppressive responses. Despite the clear role of eATP as a plant damage-associated molecular pattern, the function of its nucleoside, eAdo, and of the eAdo/eATP balance in plant stress response remain to be fully elucidated. This is particularly relevant in the context of plant-microbe interaction, where the intruder manipulates the extracellular matrix. Here, we identify Ado as a main molecule secreted by the vascular fungus Fusarium oxysporum. We show that eAdo modulates the plant's susceptibility to fungal colonization by altering the eATP-mediated apoplastic pH homeostasis, an essential physiological player during the infection of this pathogen. Our work indicates that plant pathogens actively imbalance the apoplastic eAdo/eATP levels as a virulence mechanism.
-
- Microbiology and Infectious Disease
SARS-CoV-2 spike protein plays a key role in mediating viral entry and inducing host immune responses. It can adopt either an open or closed conformation based on the position of its receptor-binding domain (RBD). It is yet unclear what cause these conformational changes or how they influence the spike's functions. Here we show that Lys417 in the RBD plays dual roles in the spike's structure: it stabilizes the closed conformation of the trimeric spike by mediating inter-spike-subunit interactions; it also directly interacts with ACE2 receptor. Hence, a K417V mutation has opposing effects on the spike's function: it opens up the spike for better ACE2 binding while weakening the RBD's direct binding to ACE2. The net outcomes of this mutation are to allow the spike to bind ACE2 with higher probability, mediate viral entry more efficiently, but become more exposed to neutralizing antibodies. Given that residue 417 has been a viral mutational hotspot, SARS-CoV-2 may have been evolving to strike a balance between infection potency and immune evasion, contributing to its pandemic spread.