Molecular analysis of the 2012 Bundibugyo virus disease outbreak.

Bundibugyo virus (BDBV) is one of four ebolaviruses known to cause disease in humans. Bundibugyo virus disease (BVD) outbreaks occurred in 2007-2008 in Bundibugyo District, Uganda, and in 2012 in Isiro, Province Orientale, Democratic Republic of the Congo. The 2012 BVD outbreak resulted in 38 laboratory-confirmed cases of human infection, 13 of whom died. However, only 4 BDBV specimens from the 2012 outbreak have been sequenced. Here, we provide BDBV sequences from seven additional patients. Analysis of the molecular epidemiology and evolutionary dynamics of the 2012 outbreak with these additional isolates challenges the current hypothesis that the outbreak was the result of a single spillover event. In addition, one patient record indicates that BDBV's initial emergence in Isiro occurred 50 days earlier than previously accepted. Collectively, this work demonstrates how retrospective sequencing can be used to elucidate outbreak origins and provide epidemiological contexts to a medically relevant pathogen.

Asymptomatic Infection of Marburg Virus Reservoir Bats Is Explained by a Strategy of Immunoprotective Disease Tolerance.

Marburg virus (MARV) is among the most virulent pathogens of primates, including humans. Contributors to severe MARV disease include immune response suppression and inflammatory gene dysregulation ("cytokine storm"), leading to systemic damage and often death. Conversely, MARV causes little to no clinical disease in its reservoir host, the Egyptian rousette bat (ERB). Previous genomic and in vitro data suggest that a tolerant ERB immune response may underlie MARV avirulence, but no significant examination of this response in vivo yet exists. Here, using colony-bred ERBs inoculated with a bat isolate of MARV, we use species-specific antibodies and an immune gene probe array (NanoString) to temporally characterize the transcriptional host response at sites of MARV replication relevant to primate pathogenesis and immunity, including CD14+ monocytes/macrophages, critical immune response mediators, primary MARV targets, and skin at the inoculation site, where highest viral loads and initial engagement of antiviral defenses are expected. Our analysis shows that ERBs upregulate canonical antiviral genes typical of mammalian systems, such as ISG15, IFIT1, and OAS3, yet demonstrate a remarkable lack of significant induction of proinflammatory genes classically implicated in primate filoviral pathogenesis, including CCL8, FAS, and IL6. Together, these findings offer the first in vivo functional evidence for disease tolerance as an immunological mechanism by which the bat reservoir asymptomatically hosts MARV. More broadly, these data highlight factors determining disparate outcomes between reservoir and spillover hosts and defensive strategies likely utilized by bat hosts of other emerging pathogens, knowledge that may guide development of effective antiviral therapies.

"Super-Spreaders" and Person-to-Person Transmission of Andes Virus in Argentina.

BACKGROUND: From November 2018 through February 2019, person-to-person transmission of Andes virus (ANDV) hantavirus pulmonary syndrome occurred in Chubut Province, Argentina, and resulted in 34 confirmed infections and 11 deaths. Understanding the genomic, epidemiologic, and clinical characteristics of person-to-person transmission of ANDV is crucial to designing effective interventions. METHODS: Clinical and epidemiologic information was obtained by means of patient report and from public health centers. Serologic testing, contact-tracing, and next-generation sequencing were used to identify ANDV infection as the cause of this outbreak of hantavirus pulmonary syndrome and to reconstruct person-to-person transmission events. RESULTS: After a single introduction of ANDV from a rodent reservoir into the human population, transmission was driven by 3 symptomatic persons who attended crowded social events. After 18 cases were confirmed, public health officials enforced isolation of persons with confirmed cases and self-quarantine of possible contacts; these measures most likely curtailed further spread. The median reproductive number (the number of secondary cases caused by an infected person during the infectious period) was 2.12 before the control measures were enforced and decreased to 0.96 after the measures were implemented. Full genome sequencing of the ANDV strain involved in this outbreak was performed with specimens from 27 patients and showed that the strain that was present (Epuyén/18-19) was similar to the causative strain (Epilink/96) in the first known person-to-person transmission of hantavirus pulmonary syndrome caused by ANDV, which occurred in El Bolsón, Argentina, in 1996. Clinical investigations involving patients with ANDV hantavirus pulmonary syndrome in this outbreak revealed that patients with a high viral load and liver injury were more likely than other patients to spread infection. Disease severity, genomic diversity, age, and time spent in the hospital had no clear association with secondary transmission. CONCLUSIONS: Among patients with ANDV hantavirus pulmonary syndrome, high viral titers in combination with attendance at massive social gatherings or extensive contact among persons were associated with a higher likelihood of transmission. (Funded by the Ministerio de Salud y Desarrollo Social de la Nación Argentina and others.).

Transcriptomics Reveal Antiviral Gene Induction in the Egyptian Rousette Bat Is Antagonized In Vitro by Marburg Virus Infection.

The Egyptian rousette bat (ERB) is the only known Marburg virus (MARV) reservoir host. ERBs develop a productive MARV infection with low viremia and shedding but no overt disease, suggesting this virus is efficiently controlled by ERB antiviral responses. This dynamic would contrast with humans, where MARV-mediated interferon (IFN) antagonism early in infection is thought to contribute to the severe, often fatal disease. The newly-annotated ERB genome and transcriptome have now enabled us to use a custom-designed NanoString nCounter ERB CodeSet in conjunction with RNA-seq to investigate responses in a MARV-infected ERB cell line. Both transcriptomic platforms correlated well and showed that MARV inhibited the antiviral program in ERB cells, while an IFN antagonism-impaired MARV was less efficient at suppressing the response gene induction, phenotypes previously reported for primate cells. Interestingly, and despite the expansion of IFN loci in the ERB genome, neither MARV showed specific induction of almost any IFN gene. However, we detected an upregulation of putative, unannotated ERB antiviral paralogs, as well as an elevated basal expression in uninfected ERB cells of key antiviral genes.

The Egyptian Rousette Genome Reveals Unexpected Features of Bat Antiviral Immunity.

Bats harbor many viruses asymptomatically, including several notorious for causing extreme virulence in humans. To identify differences between antiviral mechanisms in humans and bats, we sequenced, assembled, and analyzed the genome of Rousettus aegyptiacus, a natural reservoir of Marburg virus and the only known reservoir for any filovirus. We found an expanded and diversified KLRC/KLRD family of natural killer cell receptors, MHC class I genes, and type I interferons, which dramatically differ from their functional counterparts in other mammals. Such concerted evolution of key components of bat immunity is strongly suggestive of novel modes of antiviral defense. An evaluation of the theoretical function of these genes suggests that an inhibitory immune state may exist in bats. Based on our findings, we hypothesize that tolerance of viral infection, rather than enhanced potency of antiviral defenses, may be a key mechanism by which bats asymptomatically host viruses that are pathogenic in humans.

Error baseline rates of five sample preparation methods used to characterize RNA virus populations.

Individual RNA viruses typically occur as populations of genomes that differ slightly from each other due to mutations introduced by the error-prone viral polymerase. Understanding the variability of RNA virus genome populations is critical for understanding virus evolution because individual mutant genomes may gain evolutionary selective advantages and give rise to dominant subpopulations, possibly even leading to the emergence of viruses resistant to medical countermeasures. Reverse transcription of virus genome populations followed by next-generation sequencing is the only available method to characterize variation for RNA viruses. However, both steps may lead to the introduction of artificial mutations, thereby skewing the data. To better understand how such errors are introduced during sample preparation, we determined and compared error baseline rates of five different sample preparation methods by analyzing in vitro transcribed Ebola virus RNA from an artificial plasmid-based system. These methods included: shotgun sequencing from plasmid DNA or in vitro transcribed RNA as a basic "no amplification" method, amplicon sequencing from the plasmid DNA or in vitro transcribed RNA as a "targeted" amplification method, sequence-independent single-primer amplification (SISPA) as a "random" amplification method, rolling circle reverse transcription sequencing (CirSeq) as an advanced "no amplification" method, and Illumina TruSeq RNA Access as a "targeted" enrichment method. The measured error frequencies indicate that RNA Access offers the best tradeoff between sensitivity and sample preparation error (1.4-5) of all compared methods.

Therapeutic efficacy of the small molecule GS-5734 against Ebola virus in rhesus monkeys.

The most recent Ebola virus outbreak in West Africa, which was unprecedented in the number of cases and fatalities, geographic distribution, and number of nations affected, highlights the need for safe, effective, and readily available antiviral agents for treatment and prevention of acute Ebola virus (EBOV) disease (EVD) or sequelae. No antiviral therapeutics have yet received regulatory approval or demonstrated clinical efficacy. Here we report the discovery of a novel small molecule GS-5734, a monophosphoramidate prodrug of an adenosine analogue, with antiviral activity against EBOV. GS-5734 exhibits antiviral activity against multiple variants of EBOV and other filoviruses in cell-based assays. The pharmacologically active nucleoside triphosphate (NTP) is efficiently formed in multiple human cell types incubated with GS-5734 in vitro, and the NTP acts as an alternative substrate and RNA-chain terminator in primer-extension assays using a surrogate respiratory syncytial virus RNA polymerase. Intravenous administration of GS-5734 to nonhuman primates resulted in persistent NTP levels in peripheral blood mononuclear cells (half-life, 14 h) and distribution to sanctuary sites for viral replication including testes, eyes, and brain. In a rhesus monkey model of EVD, once-daily intravenous administration of 10 mg kg(-1) GS-5734 for 12 days resulted in profound suppression of EBOV replication and protected 100% of EBOV-infected animals against lethal disease, ameliorating clinical disease signs and pathophysiological markers, even when treatments were initiated three days after virus exposure when systemic viral RNA was detected in two out of six treated animals. These results show the first substantive post-exposure protection by a small-molecule antiviral compound against EBOV in nonhuman primates. The broad-spectrum antiviral activity of GS-5734 in vitro against other pathogenic RNA viruses, including filoviruses, arenaviruses, and coronaviruses, suggests the potential for wider medical use. GS-5734 is amenable to large-scale manufacturing, and clinical studies investigating the drug safety and pharmacokinetics are ongoing.

Informing the Historical Record of Experimental Nonhuman Primate Infections with Ebola Virus: Genomic Characterization of USAMRIID Ebola Virus/H.sapiens-tc/COD/1995/Kikwit-9510621 Challenge Stock "R4368" and Its Replacement "R4415".

The creation of licensed medical countermeasures against Select Agents such as Ebola virus (EBOV) is critically dependent on the use of standardized reagents, assays, and animal models. We performed full genome reconstruction, population genomics, contaminant analysis, and characterization of the glycoprotein gene editing site of historical United States Army Medical Research Institute of Infectious Diseases (USAMRIID) nonhuman-primate challenge stock Ebola virus Kikwit "R4368" and its 2014 replacement "R4415." We also provide characterization of the master stock used to create "R4415." The obtained data are essential to understanding the quality of the seed stock reagents used in pivotal animal studies that have been used to inform medical countermeasure development. Furthermore, these data might add to the understanding of the influence of EBOV variant populations on pathogenesis and disease outcome and inform attempts to avoid the evolution of EBOV escape mutants in response to current therapeutics. Finally, as the primary challenge stocks have changed over time, these data will provide a baseline for understanding and correlating past and future animal study results.

De novo transcriptome reconstruction and annotation of the Egyptian rousette bat.

BACKGROUND: The Egyptian Rousette bat (Rousettus aegyptiacus), a common fruit bat species found throughout Africa and the Middle East, was recently identified as a natural reservoir host of Marburg virus. With Ebola virus, Marburg virus is a member of the family Filoviridae that causes severe hemorrhagic fever disease in humans and nonhuman primates, but results in little to no pathological consequences in bats. Understanding host-pathogen interactions within reservoir host species and how it differs from hosts that experience severe disease is an important aspect of evaluating viral pathogenesis and developing novel therapeutics and methods of prevention. RESULTS: Progress in studying bat reservoir host responses to virus infection is hampered by the lack of host-specific reagents required for immunological studies. In order to establish a basis for the design of reagents, we sequenced, assembled, and annotated the R. aegyptiacus transcriptome. We performed de novo transcriptome assembly using deep RNA sequencing data from 11 distinct tissues from one male and one female bat. We observed high similarity between this transcriptome and those available from other bat species. Gene expression analysis demonstrated clustering of expression profiles by tissue, where we also identified enrichment of tissue-specific gene ontology terms. In addition, we identified and experimentally validated the expression of novel coding transcripts that may be specific to this species. CONCLUSION: We comprehensively characterized the R. aegyptiacus transcriptome de novo. This transcriptome will be an important resource for understanding bat immunology, physiology, disease pathogenesis, and virus transmission.

Evolution and Spread of Ebola Virus in Liberia, 2014-2015.

The 2013-present Western African Ebola virus disease (EVD) outbreak is the largest ever recorded with >28,000 reported cases. Ebola virus (EBOV) genome sequencing has played an important role throughout this outbreak; however, relatively few sequences have been determined from patients in Liberia, the second worst-affected country. Here, we report 140 EBOV genome sequences from the second wave of the Liberian outbreak and analyze them in combination with 782 previously published sequences from throughout the Western African outbreak. While multiple early introductions of EBOV to Liberia are evident, the majority of Liberian EVD cases are consistent with a single introduction, followed by spread and diversification within the country. Movement of the virus within Liberia was widespread, and reintroductions from Liberia served as an important source for the continuation of the already ongoing EVD outbreak in Guinea. Overall, little evidence was found for incremental adaptation of EBOV to the human host.

Emergence of Ebola Virus Escape Variants in Infected Nonhuman Primates Treated with the MB-003 Antibody Cocktail.

MB-003, a plant-derived monoclonal antibody cocktail used effectively in treatment of Ebola virus infection in non-human primates, was unable to protect two of six animals when initiated 1 or 2 days post-infection. We characterized a mechanism of viral escape in one of the animals, after observation of two clusters of genomic mutations that resulted in five nonsynonymous mutations in the monoclonal antibody target sites. These mutations were linked to a reduction in antibody binding and later confirmed to be present in a viral isolate that was not neutralized in vitro. Retrospective evaluation of a second independent study allowed the identification of a similar case. Four SNPs in previously identified positions were found in this second fatality, suggesting that genetic drift could be a potential cause for treatment failure. These findings highlight the importance selecting different target domains for each component of the cocktail to minimize the potential for viral escape.