Molecularly barcoded Zika virus libraries to probe in vivo evolutionary dynamics

Shedding dynamics of a DNA virus population during acute and long-term persistent infection

Although much is known of the molecular mechanisms of virus infection within cells, substantially less is understood about within-host infection. Such knowledge is key to understanding how viruses take up residence and transmit infectious virus, in some cases throughout the life of the host. Here, using murine polyomavirus (muPyV) as a tractable model, we monitor parallel infections of thousands of differentially barcoded viruses within a single host. In individual mice, we show that numerous viruses (>2600) establish infection and are maintained for long periods post-infection. Strikingly, a low level of many different barcodes is shed in urine at all times post-infection, with a minimum of at least 80 different barcodes present in every sample throughout months of infection. During the early acute phase, bulk shed virus genomes derive from numerous different barcodes. This is followed by long term persistent infection detectable in diverse organs. Consistent with limited productive exchange of virus genomes between organs, each displays a unique pattern of relative barcode abundance. During the persistent phase, constant low-level shedding of typically hundreds of barcodes is maintained but is overlapped with rare, punctuated shedding of high amounts of one or a few individual barcodes. In contrast to the early acute phase, these few infrequent highly shed barcodes comprise the majority of bulk shed genomes observed during late times of persistent infection, contributing to a stark decrease in bulk barcode diversity that is shed over time. These temporally shifting patterns, which are conserved across hosts, suggest that polyomaviruses balance continuous transmission potential with reservoir-driven high-level reactivation. This offers a mechanistic basis for polyomavirus ubiquity and long-term persistence, which are typical of many DNA viruses.

Author Summary / Importance

Polyomavirus infections, mostly benign but potentially fatal for immunocompromised individuals, undergo acute and long-term persistent infections. Typically, polyomavirus-associated diseases arise due to virus infection occurring in the context of a persistently infected individual. However, little is understood regarding the mechanisms of how polyomaviruses establish, maintain, and reactivate from persistent infection. We developed a non-invasive virus shedding assay combining barcoded murine polyomavirus, massively parallel sequencing technology, and novel computational approaches to track long-term infections in mice. We expect these methods to be of use not only to the study of DNA viruses but also for understanding persitent infection of diverse microbes. The study revealed organ-specific virus reservoirs and two distinct shedding patterns: constant low-level shedding of numerous barcodes and episodic high-level shedding of few barcodes. Over time, the diversity of shed barcodes decreased substantially. These findings suggest a persistent low-level infection in multiple reservoirs, with occasional bursts of replication in a small subset of infected cells. This combination of broad reservoirs and varied shedding mechanisms may contribute to polyomavirus success in transmission and maintaining long-term infections.

Barcoded SARS-CoV-2 viruses define the impact of duration and route of exposure on the transmission bottleneck in a hamster model

The transmission bottleneck, defined as the number of viruses shed from one host to infect another, is an important determinant of the rate of virus evolution and the level of immunity required to protect against virus transmission. Despite its importance, the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission bottleneck remains poorly characterized. We adapted a SARS-CoV-2 reverse genetics system to generate a pool of >200 isogenic SARS-CoV-2 viruses harboring specific 6-nucleotide barcodes, infected donor hamsters with this pool, and exposed contact hamsters to paired infected donors, varying the duration and route of exposure. Following exposure, the nasal turbinates, trachea, and lungs were collected and the number of barcodes in each tissue was enumerated. We found that longer and more direct exposures increased the transmission bottleneck and that the upper airway is the primary source of transmitted virus in this model. Together, these findings highlight the utility of barcoded viruses as tools to rigorously study virus transmission.

Barcoded SARS-CoV-2 viruses define the impact of time and route of transmission on the transmission bottleneck in a Syrian hamster model

The transmission bottleneck, defined as the number of viruses that transmit from one host to infect another, is an important determinant of the rate of virus evolution and the level of immunity required to protect against virus transmission. Despite its importance, SARS-CoV-2's transmission bottleneck remains poorly characterized, in part due to a lack of quantitative measurement tools. To address this, we adapted a SARS-CoV-2 reverse genetics system to generate a pool of >200 isogenic SARS-CoV-2 viruses harboring specific 6-nucleotide barcodes inserted in ORF10, a non-translated ORF. We directly inoculated donor Syrian hamsters intranasally with this barcoded virus pool and exposed a paired naïve contact hamster to each donor. Following exposure, the nasal turbinates, trachea, and lungs were collected, viral titers were measured, and the number of barcodes in each tissue were enumerated to quantify the transmission bottleneck. The duration and route (airborne, direct contact, and fomite) of exposure were varied to assess their impact on the transmission bottleneck. In airborne-exposed hamsters, the transmission bottleneck increased with longer exposure durations. We found that direct contact exposure produced the largest transmission bottleneck (average 27 BCs), followed by airborne exposure (average 16 BCs) then fomite exposure (average 8 BCs). Interestingly, we detected unique BCs in both the upper and lower respiratory tract of contact animals from all routes of exposure, suggesting that SARS-CoV-2 can directly infect hamster lungs. Altogether, these findings highlight the utility of barcoded viruses as tools to rigorously study virus transmission. In the future, barcoded SARS-CoV-2 will strengthen studies of immune factors that influence virus transmission.

Asian Zika virus can acquire generic African-lineage mutations during in utero infection

The Zika virus 2015 epidemic showed an unusual phenotype for human flaviviruses, specifically fetal infection. We previously showed that in utero inoculation with the Asian Zika virus isolated from the human sample causes persistent infection in porcine fetuses. Here, we characterized the evolution of the Asian Zika virus in the fetal brain and placenta. Interestingly, the Asian Zika virus acquired generic African lineage K101R (A408G) and R1609 K (G4932A) mutations during in utero infection. Both African mutations were nonsynonymous and had a high frequency of nearly 100% in the fetal brain. Then, we synthetically generated the wild-type Asian variant and fetal brain-specific variant with generic African-lineage K101R and R1609 K mutations. In mosquito C6/36 cells, but not in human and pig cells, the fetal brain-specific variant showed higher virus loads compared to the Asian wild-type prototype. While in utero infection with both variants caused comparable virus loads in the placenta and amniotic fluids, fetuses injected with the fetal brain-specific variant had the trend to higher virus loads in lymph nodes. Also, introduced K101R and R1609 K mutations were stable and had high nearly 100% frequency at 28 days after in utero inoculation in both directly injected and trans-infected fetuses. These findings evoke concerns because Zika persists in pig herds and mosquitoes on farms in Mexico. It will be essential to identify how persistent in utero infection affects virus evolution and whether in utero-emerged Zika variants have the potential for shedding into the environment, more efficient transmission, and more aggressive infection phenotypes.

Gain without pain: adaptation and increased virulence of Zika virus in vertebrate host without fitness cost in mosquito vector

IMPORTANCE

Previously, we modeled direct transmission chains of Zika virus (ZIKV) by serially passaging ZIKV in mice and mosquitoes and found that direct mouse transmission chains selected for viruses with increased virulence in mice and the acquisition of non-synonymous amino acid substitutions. Here, we show that these same mouse-passaged viruses also maintain fitness and transmission capacity in mosquitoes. We used infectious clone-derived viruses to demonstrate that the substitution in nonstructural protein 4A contributes to increased virulence in mice.

Zika virus diversity in mice is maintained during early vertical transmission from placenta to fetus, but reduced in fetal bodies and brains at late stages of infection

Since emerging in French Polynesia and Brazil in the 2010s, Zika virus (ZIKV) has been associated with fetal congenital disease. Previous studies have compared ancestral and epidemic ZIKV strains to identify strain differences that may contribute to vertical transmission and fetal disease. However, within-host diversity in ZIKV populations during vertical transmission has not been well studied. Here, we used the established anti-interferon treated Rag1-/- mouse model of ZIKV vertical transmission to compare genomic variation within ZIKV populations in matched placentas, fetal bodies, and fetal brains via RNASeq. At early stages of vertical transmission, the ZIKV populations in the matched placentas and fetal bodies were similar. Most ZIKV single nucleotide variants were present in both tissues, indicating little to no restriction in transmission of ZIKV variants from placenta to fetus. In contrast, at later stages of fetal infection there was a sharp reduction in ZIKV diversity in fetal bodies and fetal brains. All fetal brain ZIKV populations were comprised of one of two haplotypes, containing either a single variant or three variants together, as largely homogenous populations. In most cases, the dominant haplotype present in the fetal brain was also the dominant haplotype present in the matched fetal body. However, in two of ten fetal brains the dominant ZIKV haplotype was undetectable or present at low frequencies in the matched placenta and fetal body ZIKV populations, suggesting evidence of a strict selective bottleneck and possible selection for certain variants during neuroinvasion of ZIKV into fetal brains.

Wolbachia-mediated resistance to Zika virus infection in Aedes aegypti is dominated by diverse transcriptional regulation and weak evolutionary pressures

A promising candidate for arbovirus control and prevention relies on replacing arbovirus-susceptible Aedes aegypti populations with mosquitoes that have been colonized by the intracellular bacterium Wolbachia and thus have a reduced capacity to transmit arboviruses. This reduced capacity to transmit arboviruses is mediated through a phenomenon referred to as pathogen blocking. Pathogen blocking has primarily been proposed as a tool to control dengue virus (DENV) transmission, however it works against a range of viruses, including Zika virus (ZIKV). Despite years of research, the molecular mechanisms underlying pathogen blocking still need to be better understood. Here, we used RNA-seq to characterize mosquito gene transcription dynamics in Ae. aegypti infected with the wMel strain of Wolbachia that are being released by the World Mosquito Program in Medellín, Colombia. Comparative analyses using ZIKV-infected, uninfected tissues, and mosquitoes without Wolbachia revealed that the influence of wMel on mosquito gene transcription is multifactorial. Importantly, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to pathogen blocking. Therefore, to understand the influence of Wolbachia on within-host ZIKV evolution, we characterized the genetic diversity of molecularly barcoded ZIKV virus populations replicating in Wolbachia-infected mosquitoes and found that within-host ZIKV evolution was subject to weak purifying selection and, unexpectedly, loose anatomical bottlenecks in the presence and absence of Wolbachia. Together, these findings suggest that there is no clear transcriptional profile associated with Wolbachia-mediated ZIKV restriction, and that there is no evidence for ZIKV escape from this restriction in our system.

Wolbachia -mediated resistance to Zika virus infection in Aedes aegypti is dominated by diverse transcriptional regulation and weak evolutionary pressures

A promising candidate for arbovirus control and prevention relies on replacing arbovirus-susceptible Aedes aegypti populations with mosquitoes that have been colonized by the intracellular bacterium Wolbachia and thus have a reduced capacity to transmit arboviruses. This reduced capacity to transmit arboviruses is mediated through a phenomenon referred to as pathogen blocking. Pathogen blocking has primarily been proposed as a tool to control dengue virus (DENV) transmission, however it works against a range of viruses, including Zika virus (ZIKV). Despite years of research, the molecular mechanisms underlying pathogen blocking still need to be better understood. Here, we used RNA-seq to characterize mosquito gene transcription dynamics in Ae. aegypti infected with the w Mel strain of Wolbachia that are being released by the World Mosquito Program in Medellín, Colombia. Comparative analyses using ZIKV-infected, uninfected tissues, and mosquitoes without Wolbachia revealed that the influence of w Mel on mosquito gene transcription is multifactorial. Importantly, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to pathogen blocking. Therefore, to understand the influence of Wolbachia on within-host ZIKV evolution, we characterized the genetic diversity of molecularly barcoded ZIKV virus populations replicating in Wolbachia -infected mosquitoes and found that within-host ZIKV evolution was subject to weak purifying selection and, unexpectedly, loose anatomical bottlenecks in the presence and absence of Wolbachia . Together, these findings suggest that there is no clear transcriptional profile associated with Wolbachia -mediated ZIKV restriction, and that there is no evidence for ZIKV escape from this restriction in our system.

Author Summary

When Wolbachia bacteria infect Aedes aegypti mosquitoes, they dramatically reduce the mosquitoes' susceptibility to infection with a range of arthropod-borne viruses, including Zika virus (ZIKV). Although this pathogen-blocking effect has been widely recognized, its mechanisms remain unclear. Furthermore, because Wolbachia limits, but does not completely prevent, replication of ZIKV and other viruses in coinfected mosquitoes, there is a possibility that these viruses could evolve resistance to Wolbachia -mediated blocking. Here, we use host transcriptomics and viral genome sequencing to examine the mechanisms of ZIKV pathogen blocking by Wolbachia and viral evolutionary dynamics in Ae. aegypti mosquitoes. We find complex transcriptome patterns that do not suggest a single clear mechanism for pathogen blocking. We also find no evidence that Wolbachia exerts detectable selective pressures on ZIKV in coinfected mosquitoes. Together our data suggest that it may be difficult for ZIKV to evolve Wolbachia resistance, perhaps due to the complexity of the pathogen blockade mechanism.

Gain without pain: Adaptation and increased virulence of Zika virus in vertebrate host without fitness cost in mosquito vector

Zika virus (ZIKV) is now in a post-pandemic period, for which the potential for re-emergence and future spread is unknown. Adding to this uncertainty is the unique capacity of ZIKV to directly transmit between humans via sexual transmission. Recently, we demonstrated that direct transmission of ZIKV between vertebrate hosts leads to rapid adaptation resulting in enhanced virulence in mice and the emergence of three amino acid substitutions (NS2A-A117V, NS2A-A117T, and NS4A-E19G) shared among all vertebrate-passaged lineages. Here, we further characterized these host-adapted viruses and found that vertebrate-passaged viruses also have enhanced transmission potential in mosquitoes. To understand the contribution of genetic changes to the enhanced virulence and transmission phenotype, we engineered these amino acid substitutions, singly and in combination, into a ZIKV infectious clone. We found that NS4A-E19G contributed to the enhanced virulence and mortality phenotype in mice. Further analyses revealed that NS4A-E19G results in increased neurotropism and distinct innate immune signaling patterns in the brain. None of the substitutions contributed to changes in transmission potential in mosquitoes. Together, these findings suggest that direct transmission chains could enable the emergence of more virulent ZIKV strains without compromising mosquito transmission capacity, although the underlying genetics of these adaptations are complex.

The Interaction between Oxidative Stress Biomarkers and Gut Microbiota in the Antioxidant Effects of Extracts from Sonchus brachyotus DC. in Oxazolone-Induced Intestinal Oxidative Stress in Adult Zebrafish

Oxidative stress is a phenomenon caused by an imbalance between the production and accumulation of reactive oxygen species in cells and tissues that eventually leads to the production of various diseases. Here, we investigated the antioxidant effects of the extract from Sonchus brachyotus DC. (SBE) based on the 0.2% oxazolone-induced intestinal oxidative stress model of zebrafish. Compared to the model group, the treatment group alleviated oxazolone-induced intestinal tissue damage and reduced the contents of malondialdehyde, reactive oxygen species, IL-1β, and TNF-α and then increased the contents of superoxide dismutase, glutathione peroxidase, and IL-10. The 16s rDNA gene sequencing findings demonstrated that SBE could increase the relative abundance of Fusobacteriota, Actinobacteriota, and Firmicutes and decrease the relative abundance of Proteobacteria. Based on the correlation analysis between the oxidative stress biomarkers and intestinal flora, we found that the trends of oxidative stress biomarkers were significantly correlated with intestinal microorganisms, especially at the genus level. The correlations of MDA, IL-1β, and TNF-α were significantly negative with Shewanella, while SOD, GSH-Px, and IL-10 were significantly positive with Cetobacterium, Gemmobacter, and Flavobacterium. Consequently, we concluded that the antioxidant effect of SBE was realized through the interaction between oxidative stress biomarkers and gut microbiota.

Interpreting and de-noising genetically engineered barcodes in a DNA virus

The concept of a nucleic acid barcode applied to pathogen genomes is easy to grasp and the many possible uses are straightforward. But implementation may not be easy, especially when growing through multiple generations or assaying the pathogen long-term. The potential problems include: the barcode might alter fitness, the barcode may accumulate mutations, and construction of the marked pathogens may result in unintended barcodes that are not as designed. Here, we generate approximately 5,000 randomized barcodes in the genome of the prototypic small DNA virus murine polyomavirus. We describe the challenges faced with interpreting the barcode sequences obtained from the library. Our Illumina NextSeq sequencing recalled much greater variation in barcode sequencing reads than the expected 5,000 barcodes-necessarily stemming from the Illumina library processing and sequencing error. Using data from defined control virus genomes cloned into plasmid backbones we develop a vetted post-sequencing method to cluster the erroneous reads around the true virus genome barcodes. These findings may foreshadow problems with randomized barcodes in other microbial systems and provide a useful approach for future work utilizing nucleic acid barcoded pathogens.

Human immune globulin treatment controls Zika viremia in pregnant rhesus macaques

There are currently no approved drugs to treat Zika virus (ZIKV) infection during pregnancy. Hyperimmune globulin products such as VARIZIG and WinRho are FDA-approved to treat conditions during pregnancy such as Varicella Zoster virus infection and Rh-incompatibility. We administered ZIKV-specific human immune globulin as a treatment in pregnant rhesus macaques one day after subcutaneous ZIKV infection. All animals controlled ZIKV viremia following the treatment and generated robust levels of anti-Zika virus antibodies in their blood. No adverse fetal or infant outcomes were identified in the treated animals, yet the placebo control treated animals also did not have signs related to congenital Zika syndrome (CZS). Human immune globulin may be a viable prophylaxis and treatment option for ZIKV infection during pregnancy, however, more studies are required to fully assess the impact of this treatment to prevent CZS.

Zika virus persistence in the male macaque reproductive tract

Zika virus (ZIKV) is unique among mosquito-borne flaviviruses in that it is also vertically and sexually transmitted by humans. The male reproductive tract is thought to be a ZIKV reservoir; however, the reported magnitude and duration of viral persistence in male genital tissues vary widely in humans and non-human primate models. ZIKV tissue and cellular tropism and potential effects on male fertility also remain unclear. The objective of this study was to resolve these questions by analyzing archived genital tissues from 51 ZIKV-inoculated male macaques and correlating data on plasma viral kinetics, tissue tropism, and ZIKV-induced pathological changes in the reproductive tract. We hypothesized that ZIKV would persist in the male macaque genital tract for longer than there was detectable viremia, where it would localize to germ and epithelial cells and associate with lesions. We detected ZIKV RNA and infectious virus in testis, epididymis, seminal vesicle, and prostate gland. In contrast to prepubertal males, sexually mature macaques were significantly more likely to harbor persistent ZIKV RNA or infectious virus somewhere in the genital tract, with detection as late as 60 days post-inoculation. ZIKV RNA localized primarily to testicular stem cells/sperm precursors and epithelial cells, including Sertoli cells, epididymal duct epithelium, and glandular epithelia of the seminal vesicle and prostate gland. ZIKV infection was associated with microscopic evidence of inflammation in the epididymis and prostate gland of sexually mature males, pathologies that were absent in uninfected controls, which could have significant effects on male fertility. The findings from this study increase our understanding of persistent ZIKV infection which can inform risk of sexual transmission during assisted reproductive therapies as well as potential impacts on male fertility.

Zika virus (ZIKV) spread since 2015 led to establishment of urban epidemic cycles involving humans and Aedes mosquitoes. ZIKV is also sexually and vertically transmitted and causes congenital Zika syndrome. Together, these features show that ZIKV poses significant global public health risks. By virtue of similar reproductive anatomy and physiology to humans, macaques serve as a useful model for ZIKV infection. However, macaque studies to date have been limited by small sample size, typically 1 to 5 animals. Although mounting evidence identifies the male reproductive tract as a significant ZIKV reservoir, data regarding the duration of ZIKV persistence, potential for sexual transmission, and male genitourinary sequelae remain sparse. Here, we analyzed archived genital tissues from more than 50 ZIKV-inoculated male macaques. Our results show that ZIKV can persist in the male macaque reproductive tract after the resolution of viremia, with virus localization to sperm precursors and epithelial cells, and microscopic evidence of inflammation in the epididymis and prostate gland. Our findings help explain cases of sexual transmission of ZIKV in humans, which also carries a risk for transmission via assisted fertility procedures, even after resolution of detectable viremia.

Reinventing positive-strand RNA virus reverse genetics

Reverse genetics is the prospective analysis of how genotype determines phenotype. In a typical experiment, a researcher alters a viral genome, then observes the phenotypic outcome. Among RNA viruses, this approach was first applied to positive-strand RNA viruses in the mid-1970s and over nearly 50 years has become a powerful and widely used approach for dissecting the mechanisms of viral replication and pathogenesis. During this time the global health importance of two virus groups, flaviviruses (genus Flavivirus, family Flaviviridae) and betacoronaviruses (genus Betacoronavirus, subfamily Orthocoronavirinae, family Coronaviridae), have dramatically increased, yet these viruses have genomes that are technically challenging to manipulate. As a result, several new techniques have been developed to overcome these challenges. Here I briefly review key historical aspects of positive-strand RNA virus reverse genetics, describe some recent reverse genetic innovations, particularly as applied to flaviviruses and coronaviruses, and discuss their benefits and limitations within the larger context of rigorous genetic analysis.

The mucosal barrier and anti-viral immune responses can eliminate portions of the viral population during transmission and early viral growth

Little is known about how specific individual viral lineages replicating systemically during acute Human Immunodeficiency Virus or Simian Immunodeficiency Virus (HIV/SIV) infection persist into chronic infection. In this study, we use molecularly barcoded SIV (SIVmac239M) to track distinct viral lineages for 12 weeks after intravenous (IV) or intrarectal (IR) challenge in macaques. Two Mafa-A1*063+ cynomolgus macaques (Macaca fascicularis, CM) were challenged IV, and two Mamu-A1*001+ rhesus macaques (Macaca mulatta, RM) were challenged IR with 200,000 Infectious Units (IU) of SIVmac239M. We sequenced the molecular barcode of SIVmac239M from all animals over the 12 weeks of the study to characterize the diversity and persistence of virus lineages. During the first three weeks post-infection, we found ~70-560 times more unique viral lineages circulating in the animals challenged IV compared to those challenged IR, which is consistent with the hypothesis that the challenge route is the primary driver restricting the transmission of individual viral lineages. We also characterized the sequences of T cell epitopes targeted during acute SIV infection, and found that the emergence of escape variants in acutely targeted epitopes can occur on multiple virus templates simultaneously, but that elimination of some of these templates is likely a consequence of additional host factors. These data imply that virus lineages present during acute infection can still be eliminated from the systemic virus population even after initial selection.

Use of hemagglutinin and neuraminidase amplicon-based high-throughput sequencing with variant analysis to detect co-infection and resolve identical consensus sequences of seasonal influenza in a university setting

Background

Local transmission of seasonal influenza viruses (IVs) can be difficult to resolve. Here, we study if coupling high-throughput sequencing (HTS) of hemagglutinin (HA) and neuraminidase (NA) genes with variant analysis can resolve strains from local transmission that have identical consensus genome. We analyzed 24 samples collected over four days in January 2020 at a large university in the US. We amplified complete hemagglutinin (HA) and neuraminidase (NA) genomic segments followed by Illumina sequencing. We identified consensus complete HA and NA segments using BLASTn and performed variant analysis on strains whose HA and NA segments were 100% similar.

Results

Twelve of the 24 samples were PCR positive, and we detected complete HA and/or NA segments by de novo assembly in 83.33% (10/12) of them. Similarity and phylogenetic analysis showed that 70% (7/10) of the strains were distinct while the remaining 30% had identical consensus sequences. These three samples also had IAV and IBV co-infection. However, subsequent variant analysis showed that they had distinct variant profiles. While the IAV HA of one sample had no variant, another had a T663C mutation and another had both C1379T and C1589A.

Conclusion

In this study, we showed that HTS coupled with variant analysis of only HA and NA genes can help resolve variants that are closely related. We also provide evidence that during a short time period in the 2019–2020 season, co-infection of IAV and IBV occurred on the university campus and both 2020/2021 and 2021/2022 WHO recommended H1N1 vaccine strains were co-circulating.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06526-5.

African-Lineage Zika Virus Replication Dynamics and Maternal-Fetal Interface Infection in Pregnant Rhesus Macaques

Following the Zika virus (ZIKV) outbreak in the Americas, ZIKV was causally associated with microcephaly and a range of neurological and developmental symptoms, termed congenital Zika syndrome (CZS). The viruses responsible for this outbreak belonged to the Asian lineage of ZIKV. However, in vitro and in vivo studies assessing the pathogenesis of African-lineage ZIKV demonstrated that African-lineage isolates often replicated to high titers and caused more-severe pathology than Asian-lineage isolates. To date, the pathogenesis of African-lineage ZIKV in a translational model, particularly during pregnancy, has not been rigorously characterized. Here, we infected four pregnant rhesus macaques with a low-passage-number strain of African-lineage ZIKV and compared its pathogenesis to those for a cohort of four pregnant rhesus macaques infected with an Asian-lineage isolate and a cohort of mock-inoculated controls. The viral replication kinetics for the two experimental groups were not significantly different, and both groups developed robust neutralizing antibody titers above levels considered to be protective. There was no evidence of significant fetal head growth restriction or gross fetal harm at delivery (1 to 1.5 weeks prior to full term) in either group. However, a significantly higher burden of ZIKV viral RNA (vRNA) was found in the maternal-fetal interface tissues of the macaques exposed to an African-lineage isolate. Our findings suggest that ZIKV of any genetic lineage poses a threat to pregnant individuals and their infants.

IMPORTANCE ZIKV was first identified in 1947 in Africa, but most of our knowledge of ZIKV is based on studies of the distinct Asian genetic lineage, which caused the outbreak in the Americas in 2015 to 2016. In its most recent update, the WHO stated that improved understanding of African-lineage ZIKV pathogenesis during pregnancy must be a priority. The recent detection of African-lineage isolates in Brazil underscores the need to understand the impact of these viruses. Here, we provide the first comprehensive assessment of African-lineage ZIKV infection during pregnancy in a translational nonhuman primate model. We show that African-lineage isolates replicate with kinetics similar to those of Asian-lineage isolates and can infect the placenta. However, there was no evidence of more-severe outcomes with African-lineage isolates. Our results highlight both the threat that African-lineage ZIKV poses to pregnant individuals and their infants and the need for epidemiological and translational in vivo studies with African-lineage ZIKV.

Previous exposure to dengue virus is associated with increased Zika virus burden at the maternal-fetal interface in rhesus macaques

Concerns have arisen that pre-existing immunity to dengue virus (DENV) could enhance Zika virus (ZIKV) disease, due to the homology between ZIKV and DENV and the observation of antibody-dependent enhancement (ADE) among DENV serotypes. To date, no study has examined the impact of pre-existing DENV immunity on ZIKV pathogenesis during pregnancy in a translational non-human primate model. Here we show that macaques with a prior DENV-2 exposure had a higher burden of ZIKV vRNA in maternal-fetal interface tissues as compared to DENV-naive macaques. However, pre-existing DENV immunity had no detectable impact on ZIKV replication kinetics in maternal plasma, and all pregnancies progressed to term without adverse outcomes or gross fetal abnormalities detectable at delivery. Understanding the risks of ADE to pregnant women worldwide is critical as vaccines against DENV and ZIKV are developed and licensed and as DENV and ZIKV continue to circulate.

Rapid evolution of enhanced Zika virus virulence during direct vertebrate transmission chains

Zika virus (ZIKV) has the unusual capacity to circumvent natural alternating mosquito-human transmission and be directly transmitted human-to-human via sexual and vertical routes. The impact of direct transmission on ZIKV evolution and adaptation to vertebrate hosts is unknown. Here we show that molecularly barcoded ZIKV rapidly adapted to a mammalian host during direct transmission chains in mice, coincident with the emergence of an amino acid substitution previously shown to enhance virulence. In contrast, little to no adaptation of ZIKV to mice was observed following chains of direct transmission in mosquitoes or alternating host transmission. Detailed genetic analyses revealed that ZIKV evolution in mice was generally more convergent and subjected to more relaxed purifying selection than in mosquitoes or alternate passages. These findings suggest that prevention of direct human transmission chains may be paramount to resist gains in ZIKV virulence.Importance We used experimental evolution to model chains of direct and indirect Zika virus (ZIKV) transmission by serially passaging a synthetic swarm of molecularly barcoded ZIKV within and between mosquitoes and mice. We observed that direct mouse transmission chains facilitated a rapid increase in ZIKV replication and enhanced virulence in mice. These findings demonstrate that ZIKV is capable of rapid adaptation to a vertebrate host and indicate that direct human-to-human transmission could pose a greater threat to public health than currently realized.

Population bottlenecks and founder effects: implications for mosquito-borne arboviral emergence

Transmission of arthropod-borne viruses (arboviruses) involves infection and replication in both arthropod vectors and vertebrate hosts. Nearly all arboviruses are RNA viruses with high mutation frequencies, which leaves them vulnerable to genetic drift and fitness losses owing to population bottlenecks during vector infection, dissemination from the midgut to the salivary glands and transmission to the vertebrate host. However, despite these bottlenecks, they seem to avoid fitness declines that can result from Muller's ratchet. In addition, founder effects that occur during the geographic introductions of human-amplified arboviruses, including chikungunya virus and Zika virus, can affect epidemic and endemic circulation, as well as virulence. In this Review, we discuss the role of genetic drift following population bottlenecks and founder effects in arboviral evolution and spread, and the emergence of human disease.

Genome Number and Size Polymorphism in Zika Virus Infectious Units

Zika virus (ZIKV; Flaviviridae, Flavivirus) is an arthropod-borne infection that can result in severe outcomes, particularly in fetuses infected in utero It has been assumed that infection by ZIKV, as well as other viruses, is largely initiated by individual virus particles binding to and entering a cell. However, recent studies have demonstrated that multiple virus particles are frequently delivered to a cell simultaneously and that this collective particle delivery enhances infection. ZIKV is maintained in nature between Aedes aegypti mosquitos and vertebrate hosts, including humans. Human infection is initiated through the injection of a relatively small initial inoculum comprised of a genetically complex virus population. Since most mutations decrease virus fitness, collective particle transmission could benefit ZIKV and other arthropod-borne diseases by facilitating the maintenance of genetic complexity and adaptability during infection or through other mechanisms. Therefore, we utilized a barcoded ZIKV to quantify the number of virus genomes that initiate a plaque. We found that individual plaques contain a mean of 10 infecting viral genomes (range, 1 to 212). Few plaques contained more than two dominant genomes. To determine whether multigenome infectious units consist of collectively transmitting virions, infectious units of ZIKV were then separated mechanically by centrifugation, and heavier fractions were found to contain more genomes per plaque-forming unit, with larger diameters. Finally, larger/heavier infectious units reformed after removal. These data suggest that ZIKV populations consist of a variety of infectious unit sizes, likely mostly made up of aggregates, and only rarely begin with a single virus genome.IMPORTANCE The arthropod-borne Zika virus (ZIKV) infects humans and can cause severe neurological sequelae, particularly in fetuses infected in utero How this virus has been able to spread across vast geological ranges and evolve in new host populations is not yet understood. This research demonstrates a novel mechanism of ZIKV transmission through multigenome aggregates, providing insight into ZIKV evolution, immunologic evasion, and better future therapeutic design. This study shows that ZIKV plaques result from collections of genomes rather than individual genomes, increasing the potential for interactions between ZIKV genotypes.

Mapping the evolutionary landscape of Zika virus infection in immunocompromised mice

The fundamental basis of how arboviruses evolve in nature and what regulates the adaptive process remain unclear. To address this problem, we established a Zika virus (ZIKV) vector-borne transmission system in immunocompromised mice to study the evolutionary characteristics of ZIKV infection. Using this system, we defined factors that influence the evolutionary landscape of ZIKV infection and show that transmission route and specific organ microenvironments impact viral diversity and defective viral genome production. In addition, we identified in mice the emergence of ZIKV mutants previously seen in natural infections, including variants present in currently circulating Asian and American strains, as well as mutations unique to the mouse infections. With these studies, we have established an insect-to-mouse transmission model to study ZIKV evolution in vivo. We also defined how organ microenvironments and infection route impact the ZIKV evolutionary landscape, providing a deeper understanding of the factors that regulate arbovirus evolution and emergence.

Rolling circle amplification: A high fidelity and efficient alternative to plasmid preparation for the rescue of infectious clones

Alphaviruses (genus Alphavirus; family Togaviridae) are a medically relevant family of viruses that include chikungunya virus and Mayaro virus. Infectious cDNA clones of these viruses are necessary molecular tools to understand viral biology. Traditionally, rescuing virus from an infectious cDNA clone requires propagating plasmids in bacteria, which can result in mutations in the viral genome due to bacterial toxicity or recombination and requires specialized equipment and knowledge to propagate the bacteria. Here, we present an alternative- rolling circle amplification (RCA), an in vitro technology. We demonstrate that the viral yield of transfected RCA product is comparable to midiprepped plasmid, albeit with a slight delay in kinetics. RCA, however, is cheaper and less time-consuming. Further, sequential RCA did not introduce mutations into the viral genome, subverting the need for glycerol stocks and retransformation. These results indicate that RCA is a viable alternative to traditional plasmid-based approaches to viral rescue.

Evolving viral and serological stages of Zika virus RNA-positive blood donors and estimation of incidence of infection during the 2016 Puerto Rican Zika epidemic: an observational cohort study

BACKGROUND

Puerto Rico began screening blood donations for Zika virus RNA with nucleic acid amplification tests (NAATs) on April 3, 2016, because of an emerging Zika virus outbreak. We followed up positive donors to assess the dynamics of viral and serological markers during the early stages of Zika virus infection and update the estimate of infection incidence in the Puerto Rican population during the outbreak.

METHODS

Blood donations from volunteer donors in Puerto Rico were screened for the presence of Zika virus RNA using the cobas Zika NAAT. Positive donations were further tested to confirm infection, estimate viral load, and identify Zika virus-specific IgM antibodies. Individuals with positive blood donations were invited to attend follow-up visits. Donations with confirmed infection (defined as detection of Zika virus RNA or IgM on additional testing of index or follow-up samples) were assessed for stage of infection according to Zika virus RNA detectability in simulated minipools, viral load, and Zika virus IgM status. A three-step process was used to estimate the mean duration of NAAT reactivity of Zika virus in human plasma from individuals identified pre-seroconversion with at least one follow up visit and to update the 2016 incidence estimate of Zika virus infection.

FINDINGS

Between April 3 and Dec 31, 2016, 53 112 blood donations were screened for Zika virus, of which 351 tested positive, 339 had confirmed infections, and 319 could be staged. Compared with IgM-positive index donations (n=110), IgM-negative index donations (n=209) had higher mean viral loads (1·1 × 10⁶vs 8·3 × 10⁴ international units per mL) and were more likely to be detected in simulated minipools (93% [n=194] vs 26% [n=29]). The proportions of donations with confirmed infections that had viral RNA detected only in individual-donation NAATs (ie, not in simulated minipools) and were IgM positive increased as the epidemic evolved. The estimated mean duration of NAAT detectability in the 140 donors included in the follow-up study was 11·70 days (95% CI 10·06-14·36). Applying this detection period to the observed proportion of donations that were confirmed NAAT positive yielded a Zika virus seasonal incidence estimate of 21·1% (95% CI 18·1-24·1); 768 101 infections in a population of 3 638 773 in 2016.

INTERPRETATION

Characterisation of early Zika virus infection has implications for blood safety because infectivity of blood donations and utility of screening methods likely correlate with viral load and serological stage of infection. Our findings also have implications for diagnostic testing, public health surveillance, and epidemiology, and we estimate that around 21% of the Puerto Rican population was infected during the 2016 outbreak.

FUNDING

Biomedical Advanced Research and Development Authority, National Heart, Lung, and Blood Institute.

Within-Host Viral Diversity: A Window into Viral Evolution

The evolutionary dynamics of a virus can differ within hosts and across populations. Studies of within-host evolution provide an important link between experimental studies of virus evolution and large-scale phylodynamic analyses. They can determine the extent to which global processes are recapitulated on local scales and how accurately experimental infections model natural ones. They may also inform epidemiologic models of disease spread and reveal how host-level dynamics contribute to a virus's evolution at a larger scale. Over the last decade, advances in viral sequencing have enabled detailed studies of viral genetic diversity within hosts. I review how within-host diversity is sampled, measured, and expressed, and how comparative studies of viral diversity can be leveraged to elucidate a virus's evolutionary dynamics. These concepts are illustrated with detailed reviews of recent research on the within-host evolution of influenza virus, dengue virus, and cytomegalovirus.

Rapid Dissemination and Monopolization of Viral Populations in Mice Revealed Using a Panel of Barcoded Viruses

The gastrointestinal tract presents a formidable barrier for pathogens to initiate infection. Despite this barrier, enteroviruses, including coxsackievirus B3 (CVB3), successfully penetrate the intestine to initiate infection and spread systemically prior to shedding in stool. However, the effect of the gastrointestinal barrier on CVB3 population dynamics is relatively unexplored, and the selective pressures acting on CVB3 in the intestine are not well characterized. To examine viral population dynamics in orally infected mice, we produced over 100 CVB3 clones harboring nine unique nucleotide "barcodes." Using this collection of barcoded viruses, we found diverse viral populations throughout each mouse within the first day postinfection, but by 48 h the viral populations were dominated by fewer than three barcoded viruses in intestinal and extraintestinal tissues. Using light-sensitive viruses to track replication status, we found that diverse viruses had replicated prior to loss of diversity. Sequencing whole viral genomes from samples later in infection did not reveal detectable viral adaptations. Surprisingly, orally inoculated CVB3 was detectable in pancreas and liver as soon as 20 min postinoculation, indicating rapid systemic dissemination. These results suggest rapid dissemination of diverse viral populations, followed by a major restriction in population diversity and monopolization in all examined tissues. These results underscore a complex dynamic between dissemination and clearance for an enteric virus.IMPORTANCE Enteric viruses initiate infection in the gastrointestinal tract but can disseminate to systemic sites. However, the dynamics of viral dissemination are unclear. In this study, we created a library of 135 barcoded coxsackieviruses to examine viral population diversity across time and space following oral inoculation of mice. Overall, we found that the broad population of viruses disseminates early, followed by monopolization of mouse tissues with three or fewer pool members at later time points. Interestingly, we detected virus in systemic tissues such as pancreas and liver just 20 min after oral inoculation. These results suggest rapid dissemination of diverse viral populations, followed by a major restriction in population diversity and monopolization in all examined tissues.

Mutations present in a low-passage Zika virus isolate result in attenuated pathogenesis in mice

Zika virus (ZIKV) infection can result in neurological disorders including Congenital Zika Syndrome in infants exposed to the virus in utero. Pregnant women can be infected by mosquito bite as well as by sexual transmission from infected men. Herein, the variants of ZIKV within the male reproductive tract and ejaculates were assessed in inoculated mice. We identified two non-synonymous variants at positions E-V330L and NS1-W98G. These variants were also present in the passage three PRVABC59 isolate and infectious clone relative to the patient serum PRVABC59 sequence. In subsequent studies, ZIKV E-330L was less pathogenic in mice than ZIKV E-330V as evident by increased average survival times. In Vero cells, ZIKV E-330L/NS1-98G outcompeted ZIKV E-330V/NS1-98W within 3 passages. These results suggest that the E-330L/NS1-98G variants are attenuating in mice and were enriched during cell culture passaging. Cell culture propagation of ZIKV could significantly affect animal model development and vaccine efficacy studies.

Dengue and Zika Virus Cross-Reactive Human Monoclonal Antibodies Protect against Spondweni Virus Infection and Pathogenesis in Mice

Spondweni virus (SPOV) is the flavivirus that is most closely related to Zika virus (ZIKV). Although SPOV causes sporadic human infections in Africa, recently it was found in Culex mosquitoes in Haiti. To investigate the pathogenic spectrum of SPOV, we developed infection models in mice. Although two SPOV strains failed to cause disease in immunocompetent mice, each accumulated in the brain, spleen, eye, testis, and kidney when type I interferon signaling was blocked and unexpectedly caused infection, immune cell infiltration, and swelling in the ankle. In pregnant mice, SPOV replicated in the placenta and fetus but did not cause placental insufficiency or microcephaly. We identified human antibodies from ZIKV or DENV immune subjects that neutralized SPOV infection and protected against lethal challenge. Our experiments describe similarities and differences in clinical syndromes between SPOV and ZIKV and suggest that their serological relatedness has implications for antibody therapeutics and flavivirus vaccine development.

Antibody responses to Zika virus proteins in pregnant and non-pregnant macaques

The specificity of the antibody response against Zika virus (ZIKV) is not well-characterized. This is due, in part, to the antigenic similarity between ZIKV and closely related dengue virus (DENV) serotypes. Since these and other similar viruses co-circulate, are spread by the same mosquito species, and can cause similar acute clinical syndromes, it is difficult to disentangle ZIKV-specific antibody responses from responses to closely-related arboviruses in humans. Here we use high-density peptide microarrays to profile anti-ZIKV antibody reactivity in pregnant and non-pregnant macaque monkeys with known exposure histories and compare these results to reactivity following DENV infection. We also compare cross-reactive binding of ZIKV-immune sera to the full proteomes of 28 arboviruses. We independently confirm a purported ZIKV-specific IgG antibody response targeting ZIKV nonstructural protein 2B (NS2B) that was recently reported in ZIKV-infected people and we show that antibody reactivity in pregnant animals can be detected as late as 127 days post-infection (dpi). However, we also show that these responses wane over time, sometimes rapidly, and in one case the response was elicited following DENV infection in a previously ZIKV-exposed animal. These results suggest epidemiologic studies assessing seroprevalence of ZIKV immunity using linear epitope-based strategies will remain challenging to interpret due to susceptibility to false positive results. However, the method used here demonstrates the potential for rapid profiling of proteome-wide antibody responses to a myriad of neglected diseases simultaneously and may be especially useful for distinguishing antibody reactivity among closely related pathogens.

Animal Models of Zika Virus Infection during Pregnancy

Zika virus (ZIKV) emerged suddenly in the Americas in 2015 and was associated with a widespread outbreak of microcephaly and other severe congenital abnormalities in infants born to mothers infected during pregnancy. Vertical transmission of ZIKV in humans was confirmed when viral RNA was detected in fetal and placental tissues, and this outcome has been recapitulated experimentally in animals. Unlike other flaviviruses, ZIKV is both arthropod- and sexually-transmitted, and has a broad tissue tropism in humans, including multiple tissues of the reproductive tract. The threats posed by ZIKV have prompted the development of multiple in vivo models to better understand the pathogenesis of ZIKV, particularly during pregnancy. Here, we review the progress on animal models of ZIKV infection during pregnancy. These studies have generated a foundation of insights into the biology of ZIKV, and provide a means for evaluating vaccines and therapeutics.

Mathematical modeling of within-host Zika virus dynamics

Recent Zika virus outbreaks have been associated with severe outcomes, especially during pregnancy. A great deal of effort has been put toward understanding this virus, particularly the immune mechanisms responsible for rapid viral control in the majority of infections. Identifying and understanding the key mechanisms of immune control will provide the foundation for the development of effective vaccines and antiviral therapy. Here, we outline a mathematical modeling approach for analyzing the within-host dynamics of Zika virus, and we describe how these models can be used to understand key aspects of the viral life cycle and to predict antiviral efficacy.

Fetal demise and failed antibody therapy during Zika virus infection of pregnant macaques

Zika virus (ZIKV) infection of pregnant women is associated with pathologic complications of fetal development. Here, we infect pregnant rhesus macaques (Macaca mulatta) with a minimally passaged ZIKV isolate from Rio de Janeiro, where a high rate of fetal development complications was observed. The infection of pregnant macaques with this virus results in maternal viremia, virus crossing into the amniotic fluid (AF), and in utero fetal deaths. We also treated three additional ZIKV-infected pregnant macaques with a cocktail of ZIKV-neutralizing human monoclonal antibodies (nmAbs) at peak viremia. While the nmAbs can be effective in clearing the virus from the maternal sera of treated monkeys, it is not sufficient to clear ZIKV from AF. Our report suggests that ZIKV from Brazil causes fetal demise in non-human primates (NHPs) without additional mutations or confounding co-factors. Treatment with a neutralizing anti-ZIKV nmAb cocktail is insufficient to fully stop vertical transmission.

Zika virus (ZIKV) infection in pregnant women has been associated with fetal developmental defects. Here, the authors show that a Brazilian ZIKV isolate causes fetal demise in non-human primates and that antibody treatment at time of peak viremia is insufficient to clear ZIKV replication from amniotic fluid.