Postnatal Zika virus infection is associated with persistent abnormalities in brain structure, function, and behavior in infant macaques

See allHide authors and affiliations

Science Translational Medicine  04 Apr 2018:
Vol. 10, Issue 435, eaao6975
DOI: 10.1126/scitranslmed.aao6975

Postnatal perturbation by Zika virus

Much of the concern surrounding Zika virus infections focuses on fetuses infected in utero. Mavigner et al. reasoned that this neurotropic virus may have deleterious effects even after birth, so they set up a postnatal infection model to investigate. They found that infant rhesus macaques infected with Zika virus also had peripheral and central nervous system pathology. Longitudinal magnetic resonance imaging studies revealed that macaques that had been infected with Zika virus had structural and functional abnormalities and also altered emotional responses. These differences persisted months after the virus had been cleared. Although the work involved a small number of animals, their results suggest that infants and young children exposed to Zika virus should undergo more than just routine monitoring.


The Zika virus (ZIKV) epidemic is associated with fetal brain lesions and other serious birth defects classified as congenital ZIKV syndrome. Postnatal ZIKV infection in infants and children has been reported; however, data on brain anatomy, function, and behavioral outcomes following infection are absent. We show that postnatal ZIKV infection of infant rhesus macaques (RMs) results in persistent structural and functional alterations of the central nervous system compared to age-matched controls. We demonstrate ZIKV lymphoid tropism and neurotropism in infant RMs and histopathologic abnormalities in the peripheral and central nervous systems including inflammatory infiltrates, astrogliosis, and Wallerian degeneration. Structural and resting-state functional magnetic resonance imaging (MRI/rs-fMRI) show persistent enlargement of lateral ventricles, maturational changes in specific brain regions, and altered functional connectivity (FC) between brain areas involved in emotional behavior and arousal functions, including weakened amygdala-hippocampal connectivity in two of two ZIKV-infected infant RMs several months after clearance of ZIKV RNA from peripheral blood. ZIKV infection also results in distinct alterations in the species-typical emotional reactivity to acute stress, which were predicted by the weak amygdala-hippocampal FC. We demonstrate that postnatal ZIKV infection of infants in this model affects neurodevelopment, suggesting that long-term clinical monitoring of pediatric cases is warranted.


Zika virus (ZIKV) is a mosquito-borne flavivirus associated with severe birth defects when women are infected during pregnancy (13). ZIKV was first isolated in 1947 from rhesus macaques (RMs) (4), then from humans in 1952 (5), and its passage into the brain and ability to induce pathological changes have been reported since the late 1950s (68). Recent outbreaks resulting in congenital defects have caused ZIKV to be considered a global health emergency (9). Congenital infection with ZIKV occurs throughout gestation with resultant microcephaly and other brain malformations in fetuses and newborns (1013).

Whereas critical steps in the formation of the central nervous system take place in utero, postnatal development of the brain is a highly dynamic process. In the first year of life, drastic and rapid changes occur, including doubling of brain volume, increased synapse formation, pruning, and myelination (1417). Limbic regions that regulate behavior, such as the hippocampus, undergo substantial postnatal growth and maturation (18). The prolonged synaptic proliferation and neuronal maturation during the postnatal period and early childhood not only contribute to social learning and plasticity but also allow for environmental factors to affect the adult configuration (19). We hypothesized that infection with ZIKV during the first year of life may therefore have adverse neuroanatomic and functional consequences that lead to persistent behavioral abnormalities.

The potential for ZIKV to negatively affect brain development in the first year of life is suggested by other viral infections that can be acquired during infancy, such as cytomegalovirus (CMV) and HIV. These infections are transmitted not only in utero but also postnatally via breastfeeding and can lead to significant long-term neurologic complications, such as encephalopathy and cognitive impairment (2022). The route of postnatal transmission of ZIKV differs from CMV and HIV (arthropod versus breast milk). However, we posit that, when infection occurs during the critical period of brain growth in infancy, like CMV and HIV, postnatal ZIKV infection may cause sequelae similar to that observed in congenital infection. In addition, the demonstrated neurotropism of ZIKV in both fetuses and adults suggests that a similar pattern of infection and resultant neuronal damage may be seen in infants. Postnatal ZIKV of infants and children is reported (23, 24); however, detailed clinical data regarding the outcome of early postnatal ZIKV infection in humans are lacking. Without these data, the Centers for Disease Control and Prevention currently recommends only routine pediatric follow-up for infants and children with postnatal acquisition of ZIKV (25).

We sought to evaluate the long-term impact of postnatal ZIKV infection on infant brain development by performing a longitudinal study in infant RMs infected 5 weeks after birth. The model of ZIKV infection of adult macaques has been established to recapitulate key features of human infection (2629). Our evaluation of experimental postnatal ZIKV infection in infant RMs included structural and functional magnetic resonance imaging (MRI) brain scans at 3 and 6 months of age, as well as testing of socioemotional behavior at 6 and 12 months of age. We here report that postnatal ZIKV infection in infant RMs is associated with neuroinvasion, enlargement of lateral ventricles, and specific maturational alterations in the hippocampus and other brain regions. Furthermore, we observed altered functional connectivity (FC) between brain regions regulating emotional behavior and arousal functions (for example, between amygdala and hippocampus) that were predictive of distinct alterations in the species-typical emotional reactivity to acute stress. Together, this work demonstrates the pathologic features of postnatal ZIKV infection in the first year of life.


Postnatal ZIKV infection: Viral dynamics and tissue tropism

Eight infant RMs were examined, six of which were challenged subcutaneously at a median age of 37.5 days (29 to 38 days) with 105 plaque-forming units (pfu) of ZIKV PRVABC59 (30). Two of these ZIKV-infected infants and two age-matched, similarly reared, uninfected control infants (all females) were followed until 12 months of age. The remaining four ZIKV-infected infant RMs (two females and two males) were euthanized during the first 2 weeks after infection for tissue analysis (fig. S1). Clinical and hematological parameters remained stable in infant RMs exposed to ZIKV postnatally (fig. S2). ZIKV RNA in plasma peaked at days 2 and 3 after infection, with concentrations ranging from 1.7 × 106 to 3 × 107 copies of ZIKV RNA per milliliter of plasma (Fig. 1A). Consistent with previous reports in adult macaques (28, 29), ZIKV RNA was rapidly cleared from plasma, reaching undetectable concentrations by day 7 after infection (Fig. 1A). We did not detect ZIKV RNA in saliva or cerebrospinal fluid (CSF) from infant RMs, and ZIKV RNA was quantified in the urine of three infants at very low concentrations (15 to 37 copies/ml) only at days 2 and 3 after infection (table S1).

Fig. 1 ZIKV RNA in blood and tissues.

(A) Plasma ZIKV RNA concentration as measured by RT-PCR. (B) Assessment of ZIKV RNA concentration by RT-PCR in lymphoid (top) and nonlymphoid tissues (bottom) at day 3 or days 14 and 15 post-infection (p.i.). Low RNA indicates that <20 ng/μl was run in the PCR. LN, lymph node. (C) Representative detection of ZIKV RNA by RNAscope in situ hybridization in nervous system tissues. Arrowheads indicate some of the most intense RNAscope signals (red). No arrowheads were placed on the bottom left panel because the RNAscope signal was diffusely positive in the cerebellum at this magnification. Scale bars, 100 μm (left panels) and 50 μm (right panels). Four to five images per brain area were examined.

To assess ZIKV tissue tropism, we euthanized two ZIKV-infected infant RMs at predicted peak viremia (day 3 after infection) and another two after resolution of viremia (days 14 and 15 after infection). At days 3, 14, and 15 after infection, ZIKV RNA was detected by quantitative reverse transcription polymerase chain reaction (qRT-PCR) at concentrations ranging from 1 × 102 to 6 × 105 copies per microgram of RNA in lymph nodes located in the head, neck, axilla, abdomen, pelvis, and spleen (Fig. 1B). In addition, at days 14 and 15 after infection, ZIKV RNA was quantified by RT-PCR in several peripheral and central nervous system tissues including cauda equina, trigeminal ganglion, frontal cortex, parietal cortex, and occipital cortex (Fig. 1B). ZIKV RNA was also detected in situ by RNAscope hybridization in the cauda equina and in the granular cell layer of the cerebellum at days 14 and 15 after infection (Fig. 1C). Confocal microscopy analyses of these tissues demonstrated that productively infected ZIKV RNA–positive cells were CD68/CD163+ microglial cells in the cauda equina and cerebellum, as well as NeuN+ neurons in the granular cell layer of the cerebellum (fig. S3).

Humoral immune responses to postnatal ZIKV infection

Postnatal ZIKV infection of infant RMs resulted in the rapid generation of humoral immune responses, with ZIKV Envelope (E)–specific immunoglobulin M (IgM) detected in plasma of infant RMs by day 7 after infection [with concentrations declining by day 58 after infection (Fig. 2A)] and ZIKV E-specific IgG detected in infant RMs by day 10 after infection [with concentrations persisting for up to 142 days after infection (Fig. 2B)]. Using focus reduction neutralization test (FRNT) (31), anti-ZIKV neutralization activity was measured in plasma of ZIKV-infected infant RMs by day 10 after infection and persisted for up to 142 days after infection (Fig. 2C). These results suggest that infant RMs mount a similar humoral immune response to ZIKV as described in adult RMs (28).

Fig. 2 Humoral immune response to ZIKV infection.

(A) ZIKV E-specific IgM plasma concentrations measured by enzyme-linked immunosorbent assay (ELISA). (B) ZIKV E-specific IgG plasma concentrations measured by ELISA. (C) ZIKV neutralization activity of RM plasma measured by FRNT. The dashed line indicates the limit of detection for the respective assays (the lowest dilution tested). Undetectable values were set at half the lower limit of detection for visualization purposes.

Neurohistopathology of postnatal ZIKV infection

Central and peripheral nervous system tissues were further subjected to histopathological analysis and immunohistochemistry (IHC) to investigate the consequences of ZIKV neuroinvasion in postnatally infected infants as compared to age-matched controls. The entire brain and spinal cord were examined for pathologic lesions, and selected images were captured for further analysis. By day 14 after infection, consistent with ZIKV neuroinvasion, inflammatory infiltrates were observed in the brain and spinal cord of ZIKV-infected RM. A representative image of occipital cortex from a ZIKV-infected infant RM with multifocal meningeal mononuclear cell infiltrates is shown in Fig. 3A. The anterior thoracic spinal cord contained circumscribed inflammatory foci in white matter (WM), with evidence of Wallerian degeneration with multiple dilated axons and spheroid formation (Fig. 3B). Examination of the ependymal epithelial cell lining of the basal ganglia revealed hypercellularity of stromal glial cells that was not present in controls (Fig. 3C). Multifocal perivascular inflammatory cells composed predominately of plasma cells, other lymphocytes, and few macrophages were present in the cauda equina (Fig. 3D). Moreover, glial nodules were observed in WM of the occipital cortex and gray matter (GM) of the basal ganglia (Fig. 4A). IHC confirmed increased density of glial fibrillary acidic protein (GFAP)–positive reactive astrocytes as compared to controls, indicating astrogliosis in these brain regions (Fig. 4B). Caspase-3 staining revealed increased apoptosis in ZIKV-infected infant RMs compared to controls (Fig. 4C). Similar brain lesions including gliosis and axon injuries have been reported in the fetus of a pigtail macaque infected by ZIKV during gestation (26) and in autopsy findings from in utero ZIKV-infected human fetuses (2).

Fig. 3 Neurohistopathology of postnatal ZIKV infection in infant RMs.

Hematoxylin and eosin (H&E) staining of nervous tissue sections from ZIKV-infected infant RM (left) and age-matched controls (right). (A) Occipital cortex with multifocal meningeal infiltrates (arrowheads) in ZIKV-infected RM. (B) Anterior thoracic spinal cord with circumscribed inflammatory focus (arrowhead 1), dilated axons (arrowhead 2), and spheroid formation (arrowhead 3) in ZIKV-infected RM. (C) Basal ganglia showing hypercellularity with glial cells in the underlying stroma of the ependymal epithelial cell lining (arrowheads) in ZIKV-infected RM. (D) Cauda equina with multifocal perivascular inflammatory cells (arrowheads) in ZIKV-infected RM. Scale bars, 200 μm (10×) and 100 μm (20×). Four to five images per brain area were examined.

Fig. 4 Astrogliosis and apoptosis of brain regions in postnatally infected infant RMs.

Staining of central nervous system tissue sections from ZIKV-infected infant RM (left) and age-matched controls (right). (A) H&E staining with moderate astrogliosis in the WM in the periventricular region of the occipital cortex and glial nodule in the GM of the basal ganglia in ZIKV-infected RM. (B) GFAP staining of the occipital cortex and basal ganglia demonstrating reactive astrocytes in ZIKV-infected RM. (C) Caspase-3 staining of the occipital cortex and basal ganglia demonstrating increased apoptosis in ZIKV-infected RM. Scale bars, 200 μm (10×) and 100 μm (20×). Four to five images per brain area were examined.

Volumetric changes in specific brain areas: Structural MRI analyses

Having established the model of postnatal ZIKV infection of infant RMs and associated histopathology, our next objective was to investigate persistent brain structural consequences of ZIKV infection using longitudinal MRI. Structural MRI brain scans were acquired from two ZIKV-infected infant RMs and two age-matched uninfected controls at 3 and 6 months of age. Because of the small sample size, MRI findings were considered descriptive in nature; however, a consistent effect was considered potentially clinically relevant. T1-weighted structural MRI revealed no differences between ZIKV-infected RMs and controls for total brain volume (TBV) or intracranial volume, calculated as the sum of brain GM and WM and CSF (tables S2 and S3). However, the size of the lateral ventricles in the ZIKV-infected infant RMs was increased at both 3 and 6 months compared to controls (Fig. 5, A and B, and fig. S4), despite the similar increases in TBV in both groups over time. When corrected for TBV, smaller WM volume (Fig. 5B) paralleled by larger GM volume (table S2) was found in ZIKV-infected RMs at both ages compared to controls.

Fig. 5 Structural neuroimaging at 3 and 6 months of age after postnatal ZIKV infection of infant RMs.

(A) Sagittal T1-weighted structural MRI images through the corpus callosum (top) and hippocampus (bottom) at 6 months of age; yellow arrows illustrate the increased lateral ventricle volume in ZIKV-infected RMs. (B) Structural MRI measurements in cubic millimeters for TBV and lateral ventricle volume, as well as specific GM and WM areas corrected for TBV at 3 and 6 months of age in ZIKV-infected RMs and controls. The correction used was specific brain region/TBV × 100.

The protracted postnatal development of brain areas important for socioemotional behavior may make infants particularly vulnerable to ZIKV infection early in life. The hippocampus is one of the brain regions with extended postnatal development, doubling its volume from 1 week to 2 years of age in RMs (32, 33).Thus, the increase in total hippocampal volume between 3 and 6 months of age seen in the uninfected control infant RMs (Fig. 5B) is consistent with normative development in this species (18). However, ZIKV-infected infant RMs exhibited relatively blunted increases in hippocampal volume between 3 and 6 months of age (Fig. 5B). The amygdala volume, in contrast, did not differ in ZIKV-infected versus control infant RMs at either age (Fig. 5B). Postnatal ZIKV infection resulted in a greater reduction in GM volume between 3 and 6 months in the right occipital lobe (Fig. 5B) and right temporal-visual cortex (table S2) compared to age-matched uninfected controls. In summary, early postnatal ZIKV infection appears to be associated with specific brain maturational changes in cortical and subcortical regions evidenced by structural MRI.

Changes in brain FC: Resting-state functional MRI analyses

Resting-state functional MRI (rs-fMRI) can be used to measure the FC between specific brain regions of interest (ROIs) at rest, based on analysis of their temporal correlations of blood oxygen level–dependent (BOLD) signal fluctuations. Using this tool, we examined the impact of postnatal ZIKV infection on functional coupling/connectivity between ROIs during development, in addition to the brain structural impact reported in the previous section. The two ZIKV-infected infant RMs and two age-matched uninfected controls underwent rs-fMRI scans at 3 and 6 months of age. FC between ROIs was defined on the basis of published anatomical parcellations (34, 35) mapped onto the University of North Carolina (UNC)–Emory RM infant atlases registered to the F99 space (36). We first asked whether the blunted increase in hippocampal volume that we detected in ZIKV-infected RMs would have a functional effect on hippocampal connectivity. Uninfected control infant RMs exhibited strengthening in positive FC between amygdala and hippocampus (right and left hemisphere) from 3 to 6 months of age (Fig. 6, A and B). By contrast, ZIKV-infected infant RMs only exhibited this FC increase in the right hemisphere, whereas amygdala-hippocampus connectivity actually decreased in the left hemisphere from 3 to 6 months of age (Fig. 6B).

Fig. 6 Functional neuroimaging at 3 and 6 months of age after postnatal ZIKV infection of infant RMs.

(A) Depiction of amygdala-hippocampus resting-state FC measured by rs-fMRI at 6 months of age displayed in the UNC-Emory RM infant atlas (35). Images show the same two coronal slice series in each subject, with the seed (green) placed in the left amygdala and the FC r values in the hippocampus colored according to the scale provided. Results for the right amygdala seed are not shown for clarity but are provided in table S4. L, left hemisphere; R, right hemisphere; r values, raw FC correlations. (B) rs-fMRI longitudinal measurements of FC between specific ROIs at 3 and 6 months of age in ZIKV-infected RM and controls.

We next examined connections between the amygdala and hippocampus and the prefrontal cortex to better understand the impact of postnatal ZIKV infection on the limbic system network. Group differences in the directionality of FC were found between left amygdala and orbital frontal cortex (OFC) areas 11 and 13 (Fig. 6B), with ZIKV-infected infant RMs exhibiting positive connectivity at 6 months of age and uninfected controls demonstrating a negative connectivity at the same age. Weak connectivity was also observed between left hippocampus and OFC area 14 in ZIKV-infected RMs at both ages (with one of the infants showing uncoupling between these regions), whereas stronger negative connectivity developed in controls by 6 months (Fig. 6B). Last, ZIKV-infected infant RMs also showed changes in FC between visual cortical areas along the ventral motion pathway, which was of interest because of the potential for ZIKV to affect vision. FC between visual area 3 (V3) and middle temporal (MT) area in the right hemisphere increased from 3 to 6 months of age for ZIKV-infected RMs, but this trend was not detected in uninfected controls (Fig. 6B). All raw FC (r) values between ROIs are listed in table S4. These results suggest not only that postnatal ZIKV infection caused histopathological alterations in the brain and quantifiable structural brain changes by MRI scan but also that some of the structural effects resulted in functional changes, in particular alterations in FC between limbic regions important for socioemotional behavior.

Behavioral alterations after postnatal ZIKV infection: Human intruder paradigm

We next sought to assess the impact of postnatal ZIKV infection on infant behavior and emotional reactivity. The Human Intruder (HI) paradigm is a validated measure of the behavioral response to an acute stressor in nonhuman primates (Fig. 7A) (3740) similar to the task used for assessing dispositional anxiety and behavioral inhibition in children (41). The HI task robustly quantifies emotional reactivity based on the salience of threat presented by the presence and gaze direction of the intruder (table S5). Normally developing RM infants begin to modulate their emotional behavior responses based on the level of threat during the HI task at approximately 4 months of age (38, 40, 42); thus, we first conducted testing at 6 months of age. Normally developing infant RMs exhibit increased freezing and avoid vocalizations when faced with the mild threat of the intruder’s profile (Profile condition) as a defensive/anti-predator response. Subsequently, normally developing infant RMs exhibit increased hostility, anxiety, fear, and affiliative behaviors when faced with the more severe threat of the intruder’s direct gaze (Stare condition). This species-typical behavioral response to the increasing level of threat was observed in the control infant RMs (Fig. 7B). However, the ZIKV-infected infant RMs did not freeze during any condition (P = 0.039, ηp2 = 0.80; Fig. 7B), and during the Stare condition, they expressed less hostility (P = 0.029, ηp2 = 0.83) and fewer anxiety behaviors (P = 0.001, ηp2 = 0.97) than controls (Fig. 7B). In addition, ZIKV-infected RMs emitted more vocalizations and self-directed behaviors compared to controls during all conditions [P = 0.025 (ηp2 = 0.95) and P = 0.011 (ηp2 = 0.98), respectively; Fig. 7B].

Fig. 7 Assessment of the behavioral response to acute stress using the HI paradigm at 6 and 12 months of age after postnatal ZIKV infection of infant RMs.

(A) Description of the HI task measuring the ability to modulate emotional behavior based on the salience of the threat presented. Alone: the animal is placed alone in a novel room; Profile: an HI wearing a mask presents his profile to the animal; Stare: the intruder makes direct eye contact with the animal. (B) Comparison of behaviors using a standard ethogram (described in table S5) in response to the three stress levels in ZIKV-infected infant RMs and age-matched uninfected controls at 6 months of age. Bottom right graph shows pre- and post-stress plasma cortisol concentrations. n.s., not significant. (C) Correlations between the observed behavioral responses to stress with FC (r values) between specific brain ROIs in ZIKV-infected infant RMs and age-matched uninfected controls at 6 months of age. (D) Comparison of behaviors using a standard ethogram (described in table S5) in response to the three stress levels in ZIKV-infected infant RMs and age-matched uninfected controls at 12 months of age. Bottom right graph shows pre- and post-stress plasma cortisol concentration.

These data, together with the neuronal changes in temporal lobe structures that we observed by structural MRI and rs-fMRI, suggest that postnatal ZIKV infection of infant RMs may disrupt the normal maturation of brain areas important for socioemotional behavior, resulting in abnormal emotional processing. Correlations between rs-fMRI and emotional responses on the HI paradigm were conducted to examine these potential relationships. A strong negative correlation was found for left amygdala–hippocampal FC and self-directed behaviors (r = −0.98, P = 0.006; Fig. 7C), consistent with the previous finding that early hippocampal damage (and presumably decreased activity) results in increased self-directed behaviors on the HI paradigm (42). In addition, there was a strong negative correlation between left hippocampus–OFC area 14 FC and hostile behaviors (r = −0.93, P = 0.031; Fig. 7C), indicating that the weak connectivity between these regions in ZIKV-infected RM infants may predict the decreased hostile behaviors on the HI task because OFC 14 is involved in arousal and facial emotion recognition (43, 44).

Neonatal lesions to temporal lobe structures result in altered emotional behavior and changes in cortisol stress reactivity (39, 42). Because ZIKV-infected infant RMs exhibited structural and functional changes in temporal lobe structures, as well as altered emotional response, we examined stress reactivity using blood collected immediately before and after the acute stressor (HI task). Despite the structural and functional changes that we observed in the brain in ZIKV-infected infant RMs, they exhibited similar basal (pre-stress) cortisol concentrations and cortisol stress response (as measured by a percent change from pre-stressor to post-stressor) as controls (P = 0.81; Fig. 7B), suggesting that alterations after postnatal ZIKV infection occur in specific brain circuits regulating emotional responses but not in the hypothalamic-pituitary-adrenal axis at 6 months of age.

The HI paradigm was repeated in the same animals at 12 months of age to examine whether the behavioral alterations that we observed at 6 months were persistent. As shown in Fig. 7D, many of the atypical responses exhibited by the ZIKV-infected RMs continued, whereas the control RMs maintained a normal developmental response to this test of emotional reactivity. Specifically, at 12 months of age, ZIKV-infected RMs continued to demonstrate less freezing during the Profile condition (P = 0.018, ηp2 = 0.96; Fig. 7D) and to emit more vocalizations during all phases of the task compared to controls (P = 0.03, ηp2 = 0.93; Fig. 7D). In addition, ZIKV-infected RMs continued to express less hostility during the Stare condition (P = 0.042, ηp2 = 0.92; Fig. 7D). Whereas at 6 months ZIKV-infected RMs had been less anxious and exhibited more self-directed behaviors compared to controls, at 12 months ZIKV-infected RMs did not differ from controls in terms of anxiety or self-directed behaviors, largely because of wide variations in behavior between the two animals (Fig. 7, B and D). Both ZIKV-infected RMs demonstrated increased affiliative behaviors during the Stare condition compared to controls at 12 months (P = 0.035, ηp2 = 0.81; Fig. 7D), with a similar trend seen at 6 months of age (Fig. 7B). Examination of cortisol concentrations at 12 months of age revealed no difference between ZIKV-infected RMs and controls in the cortisol stress response, consistent with the findings at 6 months (Fig. 7, B and D). Mid-day basal (pre-stress) cortisol concentrations were lower in ZIKV-infected RMs compared to control RMs at 12 months of age (P = 0.022; Fig. 7D), a finding that has been previously noted in macaques with neonatal hippocampal lesions and humans with hippocampal damage (42, 45).

Sparing of visual recognition memory after postnatal ZIKV infection: Visual paired comparison task

Because we observed structural, functional, and behavioral alterations related to the hippocampus after ZIKV infection, we next interrogated another function of the hippocampus, namely, memory. The visual paired comparison (VPC) task has been widely used to study the development of visual recognition memory in infant humans and macaques (4649). Both ZIKV-infected infant RMs displayed the typical strong preference for novelty throughout the task (fig. S5), indicating intact vision and recognition memory. One control animal displayed a delay-independent, equally strong preference for either the familiar or the novel across trials. Thus, to avoid bias, we analyzed the data as different from chance, rather than preference for novelty. Overall, performance between the two groups was similar with no effect of group (P = 0.34), or of time (P = 0.106), nor was the interaction of group × time significant (P = 0.08).


This study represents a longitudinal analysis of brain structure, function, and behavior after postnatal ZIKV infection of young infant primates during an important period of brain development. Most research in newborns and infants has to date focused on the congenital ZIKV syndrome (13, 50), leaving the outcome of early postnatal infection on the developing infant brain uncertain. Reports of microcephaly developing after birth with congenital ZIKV infection (51, 52) suggest that the postnatal insult of ZIKV to the developing brain can be severe.

We found that ZIKV displays replication kinetics in infant macaques that are similar to adult macaques and humans (28, 29, 53), and initially seeds the lymph nodes and spleen. By 2 weeks after infection, ZIKV had disseminated to the peripheral and central nervous systems of postnatally infected infants, as evidenced by the presence of viral nucleic acid and histopathologic changes. Consistent with previous work, our results suggest that ZIKV can infect cerebellar neurons and microglial cells (8, 5462). Similar to the fetal brain injury observed after maternal ZIKV infection in a pigtail macaque model, as well as congenital ZIKV syndrome in humans (2, 26, 63, 64), our analysis of postnatal ZIKV infection in RM infants revealed inflammation, reactive astrocytes, gliosis in the GM and WM, axonal injury, and apoptosis in the brain and spinal cord. In particular, the consistent finding of astrogliosis between the infant infection model described here and fetal infection illustrates a common pathology of postnatal and in utero ZIKV infection (64). In our model, ZIKV RNA was not detected in the brain regions showing increased apoptosis and gliosis, suggesting that these abnormalities may be due to an indirect inflammatory mechanism rather than direct cytopathic effect of the virus.

In adult macaques, ZIKV has been detected by PCR in multiple tissues, including brain, at days 35 to 72 after infection (65, 66). In addition, ZIKV was found in the CSF of an adult RM at 42 days after inoculation (65). Although we did not specifically attempt to measure virus at these later time points, the neuroimaging and behavioral abnormalities described in this work suggest that ZIKV may lead to long-term neurobehavioral damage to primate infants. Whether these findings are due to immune events triggered during acute infection or to persistence of virus in the brain (or both) remains unknown.

Using structural MRI scans at 3 and 6 months of age, we observed a persistent reduction in WM volume and increase in ventricular size in ZIKV-infected infant RMs as compared to age-matched controls. To our knowledge, normative MRI data regarding the volume of lateral ventricles in a larger sample of control RM infants matching our study conditions are not available. However, 11 female RM at 3 months of age and 18 female RM at 6 months of age were imaged using similar protocols, as described in the UNC-Emory infant atlas (36), demonstrating an SD of 56 mm3 for lateral ventricle volumes across this dam-reared social group. The difference in the volume of the lateral ventricles in ZIKV-infected RM compared to control RM in the present work was greater, with a mean at 3 months of 341 mm3 versus 164 mm3 and at 6 months of 321 mm3 versus 184 mm3. Thus, the difference that we observed is not merely a variation from normal but rather likely indicates a true enlargement of the ventricle size after postnatal ZIKV infection. Similar but more extensive ventriculomegaly typically of the lateral ventricles has been reported in live-born infants and autopsied fetuses after in utero ZIKV infection (1, 2, 67, 68), supporting a potentially similar underlying mechanism for these effects.

We also found blunted hippocampal volume growth in ZIKV-infected infant RMs as compared to controls. Specific damage to the hippocampus has also been noted in models of ZIKV-infected adult rodents (56) and can be explained by the protracted postnatal development of this brain region (18), making it vulnerable to postnatal insults. The structural impact to the hippocampus was paralleled by alterations in its FC with other limbic regions such as the amygdala and OFC, which also demonstrated an abnormal interaction on rs-fMRI, in ZIKV-infected infant RM. The structural and functional brain changes that we observed corresponded with changes in emotional behavior. Although studies have demonstrated an impact of nursery-rearing on brain development and behavior (69, 70), rearing conditions and handling experiences were the same for control and ZIKV-infected RM infants in this study. ZIKV-infected infants were not impaired in their ability to produce emotional responses but, instead, exhibited an atypical emotional response based on the salience of the threat presented. For example, ZIKV-infected infants consistently vocalized more and froze less than controls during the mild threat of the intruder profile yet exhibited less hostility during the salient threat of the intruder stare. OFC area 14 has been implicated in arousal and facial emotion recognition (43, 44), and the weak connectivity between OFC 14 and the hippocampus seen in ZIKV-infected infant RMs may have resulted in less arousal and fewer hostile behaviors compared to uninfected controls. OFC areas 11 and 13 are involved in the ability to flexibly respond to rapidly changing social information, and cross-talk between the amygdala and OFC areas 11 and 13 guides optimal decision making (71). ZIKV-infected and control infant RM displayed opposite connectivity between these regions, possibly influencing the incorporation of new information into decisions, although this behavior was not specifically tested in this work. Neonatal lesions of the amygdala or hippocampus have resulted in similar behavioral changes exhibited by the ZIKV-infected infant RMs on the HI task in nursery-reared animals, as well as in adolescents and adults in a social cohort (39, 42, 7274). Damage to the hippocampus may also affect visual recognition memory (75). However, we found that recognition memory at 6 months of age, as measured by the VPC task, was similar in ZIKV-infected RM infants and uninfected controls, indicating selective sparing of this function. Although performance on the VPC task in adults is affected by hippocampal damage, a lack of memory deficit in infancy was also reported in animals with selective damage to the hippocampal formation (76) or exposure to anesthesia (46) that results in similar emotional changes (42, 77). RM infants with experimental selective hippocampal lesions evidenced a late-developing delay-dependent memory impairment as they reached adulthood (76), providing further rationale for longitudinal neurocognitive testing after postnatal ZIKV infection.

ZIKV has been linked with damage to the visual system, with ocular and brain visual system abnormalities seen in congenital ZIKV infection (51, 78, 79). In addition, adults with Guillain-Barre syndrome after ZIKV infection exhibit cytotoxic edema in the occipital lobe (80), the location of the visual cortex. We found that postnatal ZIKV infection in infant RM resulted in a greater reduction in GM volume between 3 and 6 months in the right occipital lobe and right temporal-visual cortex compared to age-matched uninfected controls, with abnormal inflammation, astrogliosis, and apoptosis apparent in the occipital cortex early in infection. Performance on the VPC task, however, confirmed intact visual processes necessary to encode and recognize a visual stimulus in ZIKV-infected RM infants.

There were limitations to this study. First, a small number of animals were studied; therefore, comparisons between groups are descriptive in nature. The strength of these results does not reside in statistical significance but rather in the effect size and potential clinical relevance of this work. The National Institutes of Health recently funded a study of neurodevelopmental outcomes of Guatemalan infants and children infected with ZIKV after birth (81), which may replicate our findings and/or provide new insights into postnatal ZIKV infection. Second, our results reflect infection with a single dose of a single ZIKV isolate. The PRVABC59 strain of ZIKV is an Asian genotype virus that has been widely used for RM studies and is most closely related to isolates from Brazil (30), where the largest outbreak of ZIKV infection has occurred to date. Additional work to substantiate our findings regarding postnatal ZIKV infection should use a larger cohort of infant RMs infected with multiple virus isolates that reflect the diversity of current infections. Third, it should be noted that rs-fMRI has several limitations: (i) FC was obtained under anesthesia, but we note that isoflurane doses were kept below those used in previous studies reporting patterns of coherent BOLD fluctuations similar to the awake state, including in sensory and visual systems (8285); (ii) only female subjects were studied, and there may be sex differences in the impact of ZIKV infection on FC brain development; and (iii) it is difficult to relate changes in strength or direction of FC between cortical areas with specific physiological, cellular, or molecular brain maturational processes, which will need to be addressed in future studies with larger sample sizes.

In summary, early postnatal ZIKV infection appears to be associated with histopathologic, behavioral, and brain structural and functional abnormalities during the first year of life. The impact of early postnatal ZIKV infection during later stages of life (for example, school age and adolescence) remains unknown. Previous research using macaques with hippocampal lesions made in the neonatal period demonstrated schizophrenic-like characteristics that only manifested during adolescence (86). Furthermore, studies of postnatally acquired CMV revealed that some of the cognitive sequelae do not become evident until school age (22). Our finding that postnatal ZIKV infection of infant RMs causes persistent changes consistent with delayed or dysfunctional brain maturation suggests that long-term clinical monitoring of all ZIKV-infected infants and children may be critical. Heightened surveillance of this population, with attention to neurodevelopment and behavior, is warranted. This RM model sheds insight into potential outcomes of human infants infected with ZIKV postnatally and provides a platform with which to test therapeutic approaches to alleviate the neurologic consequences of ZIKV infection.


Study design

The objective of this study was to evaluate the neurodevelopmental outcome of postnatal ZIKV infection of infant RMs. After infection, animals underwent virologic and immunologic measurements, as well as neurohistopathology, neuroimaging, and behavioral assessments. This was an exploratory pilot study of postnatal ZIKV infection in infant RMs; therefore, a formal power analysis was not performed. Infants were randomized to either the ZIKV or control groups based on the day of birth. Veterinary and animal resource staff were not blinded to study group due to the necessary environmental precautions that were taken with ZIKV. Laboratory staff were blinded to study group during the conduct of experiments and analyses.


This study was conducted in strict accordance with U.S. Department of Agriculture regulations and the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health, approved by the Emory University Institutional Animal Care and Use Committee (AWA# A3180-01), and conducted in an AAALAC (Association for Assessment and Accreditation of Laboratory Animal Care)–accredited facility. Eight infant Indian RM (Macaca mulatta) including two males and six females were used in this study. The infants were delivered naturally by their dams while housed in indoor/outdoor social groups. Infants were then removed from their dams at 5 to 10 days of age and transported to the Yerkes National Primate Research Center (YNPRC) nursery facility. All infants were singly housed in warming incubators for the first 3 to 4 weeks. Infants were hand-fed formula every 2 to 3 hours for the first month, then via self-feeders for the next 3 months as per standard YNPRC protocol. Soft blankets, plush toys, or fleece were provided and changed daily. Soft chow and fruits were introduced starting at 1 month of age, and by 4 months, formula was discontinued and all were fed a diet of Purina Primate Chow, Old World Monkey formulation, supplemented with daily fruits and vegetables. Water was provided ad libitum. At 4 weeks of age, infants transitioned into age-appropriate caging with visual and auditory contact with conspecifics. At 6 to 7 weeks of age, they began socialization consisting of protected contact housing supplemented with 2 to 3 hours daily of full contact with their age-matched peer. By 3 months of age, each pair was housed entirely together. ZIKV-infected infants were housed in an ABSL-2+ (Animal Biosafety Level 2) room until cleared of virus in blood and urine. Control infants were initially housed in a separate ABSL-2 room, and then all infant pairs were housed together after ZIKV clearance with visual contact to increase visual socialization and control for environmental factors.

Virus and infection

The ZIKV of Puerto Rican origin (PRVABC59, GenBank accession number: KU501215.1) used in this study had been passaged four times and titered on Vero cells. Experimental infections were performed via the subcutaneous route using 105 pfu of ZIKV PRVABC59.

ZIKV RNA detection by qRT-PCR

Total RNA was extracted from 140 μl of plasma and CSF samples using the QIAamp Viral RNA Mini Kit (Qiagen). RNA from urine and saliva was isolated using the QIAamp Circulating Nucleic Acid Kit (Qiagen) for the extraction of larger volumes (up to 3 ml). Tissues were suspended in TRI reagent and mechanically homogenized using beads. Total RNA was isolated from homogenized tissues using the Direct-zol RNA MiniPrep Plus Kit (Zymo Research) per the manufacturer’s instructions. ZIKV RNA standard was generated by annealing two oligonucleotides spanning the target ZIKV prM-E gene region and performing in vitro transcription using the MEGAscript T3 Transcription Kit (Ambion). Purified RNA was reverse-transcribed using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems) and random hexamers. For quantitation of viral RNA, a standard curve was generated using 10-fold dilutions of ZIKV RNA standard, and qRT-PCR was performed using TaqMan Gene Expression Master Mix (Applied Biosystems) and ZIKV primers (1 μM) and probe (250 nM). The ZIKV primer-probe set targeting the prM-E gene region included forward primer ZIKV/PR 907 (5′-TTGGTCATGATACTGCTGATTGC-3′), reverse primer ZIKV/PR 983c (5′-CCTTCCACAAAGTCCCTATTGC-3′), and probe ZIKV/PR 932 FAM (5′-CGGCATACAGCATCAGGTGCATAGGAG-3′). The probe was labeled with 6-carboxyfluorescein (FAM) at the 5′ end and two quenchers, an internal ZEN quencher and 3′ end Iowa Black FQ (Integrated DNA Technologies). The standard curve had an R2 value greater than 0.99. Viral RNA copies were interpolated from the standard curve using the sample CT value. For tissue samples, viral RNA copies were normalized to RNA input (per microgram of RNA). For plasma samples, viral RNA copies were represented as copies per milliliter of plasma.

ZIKV RNA in situ hybridization

ZIKV in situ hybridization was done on cerebellum and cauda equina using riboprobes that target multiple genes of ZIKV strains including PRVABC59 (Advanced Cell Diagnostics). RNAscope was performed as per the manufacturer’s recommended protocol. Briefly, prepared slides were baked for 1 hour at 60°C before use. After deparaffinization and hydration, tissue slides were air-dried and treated with a peroxidase blocker before heating in a target retrieval solution for 20 min at 95° to 100°C. Protease was then applied for 30 min at 40°C. Target probes were hybridized for 2 hours at 40°C, followed by a series of signal amplification and washing steps. Red chromogen development was performed following the RNAscope 2.5 HD detection protocol and reagents. The chromogenic detection was performed using Fast Red as substrate to generate red signal, followed by counterstaining with hematoxylin.


ZIKV-infected RMs sacrificed on day 3 and days 14 and 15 after infection were subjected to histopathologic analyses. For histopathologic examination, various tissue samples were fixed in 10% neutral buffered formalin, paraffin-embedded, sectioned at 5 μm, and stained with H&E. To define specific cell types, we performed IHC analysis by using antibodies against GFAP (M0761, Agilent Technologies) and caspase-3 (9662, Cell Signaling). IHC staining on sections of brain from both ZIKV-infected and uninfected animals was performed using a biotin-free polymer system (Biocare Medical). The paraffin-embedded sections were subjected to deparaffinization in xylene and rehydration in graded series of ethanol and rinsed with double-distilled water. The brain sections were incubated with rabbit anti-human GFAP and caspase-3 antibodies for 1 hour after heat-induced epitope retrieval. Antibody labeling was visualized by development of the chromogen (Warp Red or DAB Chromogen Kits, Biocare Medical). Digital images of brain were captured at ×100 and ×200 magnification with an Olympus BX43 microscope equipped with a digital camera (DP26, Olympus) and evaluated using cellSens digital imaging software 1.15 (Olympus).

MRI studies

Neuroimaging data were collected longitudinally at 3 and 6 months of age using a Siemens 3T Tim Trio system (Siemens Medical Solutions) with a Tx/Rx 8-channel volume coil. Data were acquired in a single session, which included T1- and T2-weighted structural scans, and two 15-min rs-fMRI (T2*-weighted) scans to measure temporal changes in regional BOLD signal. Animals were scanned in the supine position in the same orientation, achieved by placement and immobilization of the head in a custom-made head holder via ear bars and a mouthpiece. After initial telazol induction and intubation, scans were collected under isoflurane anesthesia (0.8 to 1% to effect), kept at the lowest dose possible to minimize effects of anesthesia on BOLD signal (8285). End-tidal CO2, inhaled CO2, O2 saturation, heart rate, respiratory rate, blood pressure, and body temperature were monitored continuously and maintained during each MRI session.

Structural MRI acquisition. High-resolution T1-weighted scans were acquired for volumetrics and registration of the functional scans using a three-dimensional magnetization-prepared rapid gradient echo (3D-MPRAGE) parallel image sequence [repetition time/echo time (TR/TE) = 2600/3.46 ms; field of view (FOV), 116 mm × 116 mm; voxel size, 0.5 mm3 isotropic; eight averages] with GeneRalized Autocalibrating Partially Parallel Acquisitions (GRAPPA) acceleration factor of R = 2. T2-weighted MRI scans were collected using a 3D fast spin-echo sequence (TR/TE = 3200/373 ms; FOV, 128 mm × 128 mm; voxel size, 0.5 mm3 isotropic; three averages) with GRAPPA (R = 2) to aid with registration and delineation of anatomical borders.

Structural MRI data processing and analysis. Data sets were processed using AutoSeg_3.3.2 segmentation package (87) to get the volumes of brain WM and GM, CSF, and cortical (temporal visual area, temporal auditory area, and prefrontal, frontal, parietal, and occipital lobes) and subcortical (hippocampus and amygdala) brain areas. Image processing steps included the following: (i) averaging T1 and T2 images to improve signal-to-noise ratio, (ii) intensity inhomogeneity correction, (iii) rigid-body registration of the subject MRI to the 3- or 6-month UNC-Emory infant RM atlases (36), (iv) tissue segmentation and skull stripping, (v) registration of the atlas to the subject’s brain to generate cortical parcellations (affine followed by deformable SyN/ANTS registration), and (vi) manual editing of the amygdala and hippocampus done on all scans using previously published neuroanatomical boundaries (32). For more detailed description of the structural MRI analysis, see (87, 88). Cortical and subcortical volume measurements (in cubic millimeters) were corrected by TBV and then multiplied by 100 (for example, hippocampal volume/TBV × 100).

Rs-fMRI acquisition. BOLD-weighted functional images were collected using a single-shot echo-planar imaging (EPI) sequence (400 volumes; TR/TE = 2060/25 ms, 2 × 15 min; voxel size, 1.5 mm3 isotropic) after the T1-MRI scan. An additional short, reverse-phase encoding scan was also acquired for unwarping susceptibility-induced distortions in the EPI images using previously validated methods (89).

Functional MRI data preprocessing. Imaging data were preprocessed with an in-house pipeline in Nipype (90) with modifications (9194), including adaptations for the RM brain (95). After the EPI functional time series were concatenated and rigid-body registered to the subject’s T1-weighted structural image, this was transformed to conform to age-specific T1-weighted RM infant atlases (36) using nonlinear registration methods in FSL (FNIRT). These atlases were registered to the 112RM-SL atlas in F99 space (96, 97). We registered the scans to the age-specific appropriate infant atlas (3 or 6 months). Frames with displacement values greater than 0.2 mm were removed (98, 99), and data were visually inspected upon preprocessing completion to exclude series with unsatisfactory co-registration or significant BOLD signal dropout.

FC analysis. FC between ROIs including the amygdala, hippocampus, OFC, and visual areas (V1, V3, and MT area) was analyzed. ROIs were defined on the basis of published anatomical parcellations (34, 35) mapped onto the UNC-Emory RM infant atlases registered to the F99 space (36). The left and right amygdala label maps were manually drawn by experts using cytoarchitectonic maps in the existing UNC-Wisconsin adolescent atlas (100), and other ROIs were manually edited in the infant atlases following established anatomical landmarks (101, 102). The time course of the BOLD signal was averaged across all voxels within each ROI using a single ROI as the seed and then correlated with the time course of the other ROIs. FC values are calculated as correlation coefficients (r values) between ROI BOLD time courses, which were extracted from correlation matrices using MATLAB (MathWorks Inc.).

HI paradigm

Infants were previously trained to quickly transfer from their home cage to a transport box. On the testing day, the subject was separated from the cagemate and transported from the nursery to a novel testing room where a basal blood sample (1 ml, awake femoral venipuncture) was collected within 10 min of initial disturbance, reflecting baseline concentrations of plasma biomarkers. The RM was then placed in a modified housing wire-mesh cage, with one side in clear acrylic to allow video recording of animal’s behavior. The paradigm consisted of three conditions: Alone condition, animal remained alone in the cage for 9 min to acclimate to the environment and obtain a baseline measure of behavior. Profile condition: An HI (researcher wearing a rubber mask) entered the room, sat on a stool 2 m from the test cage, while presenting his profile (without eye contact) to the animal for 9 min. Then, the intruder left the room, giving the animal a 3-min break. Stare condition: The HI reentered the room but this time stared directly at the animal for 9 min. A second blood draw was collected (1 ml, awake femoral draw) immediately after the task, reflecting concentrations of plasma biomarkers in response to acute stress. Blood samples were collected in prechilled plastic 2-ml vacutainer tubes containing EDTA and immediately placed on ice. Samples were centrifuged at 3000 rpm for 15 min in a refrigerated centrifuge (at 4°C). Plasma was pipetted into sterile tubes and stored at −80°C until assayed.

Statistical analysis

HI paradigm analyses were conducted using SPSS 24 for Windows (IBM Corporation), and significance was set at P < 0.05. Group differences in emotional behaviors were examined with repeated-measures analysis of variance (ANOVA), with Group (control, ZIKV) as the between-subjects factor and Condition (Alone, Profile, Stare) as the within-subject repeated measures. Effect sizes were calculated using partial eta squared (ηp2). Emotional reactivity to the intruder was assessed via videotape recording for later coding using the Observer XT 10 software (Noldus Inc.) and a detailed ethogram (table S4). Two blinded experimenters coded all of the videotapes but had a high degree of interrater reliability Cohen’s κ = 0.86 and an average intrarater reliability of Cohen’s κ = 0.97. Pearson correlation coefficient was used to assess for association between behavior on HI task and FC. Cortisol assay was analyzed by independent t tests, with Group (control, ZIKV) as the between-subject factor and percent change in cortisol concentration as the dependent variable. All primary data for this study are presented in table S6.


Materials and Methods

Fig. S1. Schematic of the study design.

Fig. S2. Clinical parameters in ZIKV-infected and uninfected control infant RMs.

Fig. S3. Identification of ZIKV target cells in infant RMs.

Fig. S4. Series of coronal T1-weighted structural MRI images through the level of the lateral ventricle at 6 months of age in control and ZIKV-infected RMs.

Fig. S5. VPC task.

Table S1. ZIKV RNA detection by RT-PCR in body fluids.

Table S2. Brain regional volumes corrected for TBV.

Table S3. Brain regional volumes corrected for intracranial volume.

Table S4. FC (r values) between brain ROIs.

Table S5. Behavioral ethogram.

Table S6. Primary data.

Reference (103)


Acknowledgments: We thank S. Ehnert and the YNPRC Division of Research Resources for expert assistance with animal procedures, M. Thompson and the Emory Environmental Health and Safety Office for supervising environmental infection control in the YNPRC nursery, and A. McElroy for training and helpful discussions. We acknowledge K. De Paris from UNC and K. Van Rompay from the California National Primate Research Center (CNPRC) for providing control infant brains for histopathology (grant numbers R01 DE022287 and R01 DE019064 to K. De Paris and P51 OD011107 to CNPRC). Funding: This work was supported by the Pilot Grant Program of the YNPRC (P51 OD011132), as well as funding from the Emory Center for Childhood Infections and Vaccines, and Children’s Healthcare of Atlanta (all to A.C.). Research reported in this publication was supported by the National Library of Medicine of the NIH under award number T15LM007088 (to E.F.). Author contributions: A.C. designed the study with input from M.M., M.C.A., M.S.S., J.W., D.H.O., T.H.V., R.P.J., D.I.W., and G.S. M.S.S. provided the virus. J.R. and M.C.A. performed the HI test and interpreted the data. Z.K.-B., M.M.S., and X.Z. interpreted the neuroimaging. D.F., E.F., E.E., and M.S. provided tools for imaging analysis. S.G. and M.W.B. performed histopathology. S.K.B. and J.W. performed and analyzed ELISA and neutralization assays. J.T.O., C.E.M., M.S.S., V.K.B., D.M.M., and D.I.W. performed and analyzed data from virologic assays with help from C.M. J.H., C.M., and M.M. processed samples. S.M.J. and J.K.C. provided animal care and collected samples. A.C., M.M., and V.A. wrote the manuscript with help from the other authors. Competing interests: D.H.O. is a paid a member of the Institutional Biosafety Committee for Takeda Vaccines, which is involved in ZIKV research. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH.
View Abstract

Navigate This Article