Distinct lineages of Ebola virus in Guinea during the 2014 West African epidemic

Molecular characterization of the 2022 Sudan virus disease outbreak in Uganda

IMPORTANCE

Ebola disease (EBOD) is a public health threat with a high case fatality rate. Most EBOD outbreaks have occurred in remote locations, but the 2013-2016 Western Africa outbreak demonstrated how devastating EBOD can be when it reaches an urban population. Here, the 2022 Sudan virus disease (SVD) outbreak in Mubende District, Uganda, is summarized, and the genetic relatedness of the new variant is evaluated. The Mubende variant exhibited 96% amino acid similarity with historic SUDV sequences from the 1970s and a high degree of conservation throughout the outbreak, which was important for ongoing diagnostics and highly promising for future therapy development. Genetic differences between viruses identified during the Mubende SVD outbreak were linked with epidemiological data to better interpret viral spread and contact tracing chains. This methodology should be used to better integrate discrete epidemiological and sequence data for future viral outbreaks.

Early warning of emerging infectious diseases based on multimodal data

The coronavirus disease 2019 (COVID-19) pandemic has dramatically increased the awareness of emerging infectious diseases. The advancement of multiomics analysis technology has resulted in the development of several databases containing virus information. Several scientists have integrated existing data on viruses to construct phylogenetic trees and predict virus mutation and transmission in different ways, providing prospective technical support for epidemic prevention and control. This review summarized the databases of known emerging infectious viruses and techniques focusing on virus variant forecasting and early warning. It focuses on the multi-dimensional information integration and database construction of emerging infectious viruses, virus mutation spectrum construction and variant forecast model, analysis of the affinity between mutation antigen and the receptor, propagation model of virus dynamic evolution, and monitoring and early warning for variants. As people have suffered from COVID-19 and repeated flu outbreaks, we focused on the research results of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and influenza viruses. This review comprehensively viewed the latest virus research and provided a reference for future virus prevention and control research.

High-throughput sequencing approaches applied to SARS-CoV-2

High-throughput sequencing is crucial for surveillance and control of viral outbreaks. During the ongoing coronavirus disease 2019 (COVID-19) pandemic, advances in the high-throughput sequencing technology resources have enhanced diagnosis, surveillance, and vaccine discovery. From the onset of the pandemic in December 2019, several genome-sequencing approaches have been developed and supported across the major sequencing platforms such as Illumina, Oxford Nanopore, PacBio, MGI DNBSEQTM and Ion Torrent. Here, we share insights from the sequencing approaches developed for sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between December 2019 and October 2022.

Distinct mutations and lineages of SARS-CoV-2 virus in the early phase of COVID-19 pandemic and subsequent 1-year global expansion

A novel coronavirus, SARS-CoV-2, has caused over 274 million cases and over 5.3 million deaths worldwide since it occurred in December 2019 in Wuhan, China. Here we conceptualized the temporospatial evolutionary and expansion dynamics of SARS-CoV-2 by taking a series of the cross-sectional view of viral genomes from early outbreak in January 2020 in Wuhan to the early phase of global ignition in early April, and finally to the subsequent global expansion by late December 2020. Based on the phylogenetic analysis of the early patients in Wuhan, Wuhan/WH04/2020 is supposed to be a more appropriate reference genome of SARS-CoV-2, instead of the first sequenced genome Wuhan-Hu-1. By scrutinizing the cases from the very early outbreak, we found a viral genotype from the Seafood Market in Wuhan featured with two concurrent mutations (i.e., M type) had become the overwhelmingly dominant genotype (95.3%) of the pandemic 1 year later. By analyzing 4013 SARS-CoV-2 genomes from different continents by early April, we were able to interrogate the viral genomic composition dynamics of the initial phase of global ignition over a time span of 14 weeks. Eleven major viral genotypes with unique geographic distributions were also identified. WE1 type, a descendant of M and predominantly witnessed in western Europe, consisted of half of all the cases (50.2%) at the time. The mutations of major genotypes at the same hierarchical level were mutually exclusive, which implies that various genotypes bearing the specific mutations were propagated during human-to-human transmission, not by accumulating hot-spot mutations during the replication of individual viral genomes. As the pandemic was unfolding, we also used the same approach to analyze 261 323 SARS-CoV-2 genomes from the world since the outbreak in Wuhan (i.e., including all the publicly available viral genomes) to recapitulate our findings over 1-year time span. By December 25, 2020, 95.3% of global cases were M type and 93.0% of M-type cases were WE1. In fact, at present all the five variants of concern (VOC) are the descendants of WE1 type. This study demonstrates that viral genotypes can be utilized as molecular barcodes in combination with epidemiologic data to monitor the spreading routes of the pandemic and evaluate the effectiveness of control measures. Moreover, the dynamics of viral mutational spectrum in the study may help the early identification of new strains in patients to reduce further spread of infection, guide the development of molecular diagnosis and vaccines against COVID-19, and help assess their accuracy and efficacy in real world at real time.

In Silico Analysis of Honeybee Venom Protein Interaction with Wild Type and Mutant (A82V + P375S) Ebola Virus Spike Protein

Venom from different organisms was used in ancient times to treat a wide range of diseases, and to combat a variety of enveloped and non-enveloped viruses. The aim of this in silico research was to investigate the impact of honeybee venom proteins and peptides against Ebola virus. In the current in silico study, different online and offline tools were used. RaptorX (protein 3D modeling) and PatchDock (protein–protein docking) were used as online tools, while Chimera and LigPlot + v2.1 were used for visualizing protein–protein interactions. We screened nine venom proteins and peptides against the normal Ebola virus spike protein and found that melittin, MCD and phospholipase A2 showed a strong interaction. We then screened these peptides and proteins against mutated strains of Ebola virus and found that the enzyme phospholipase A2 showed a strong interaction. According to the findings, phospholipase A2 found in honeybee venom may be an effective source of antiviral therapy against the deadly Ebola virus. Although the antiviral potency of phospholipase A2 has been recorded previously, this is the first in silico analysis of honeybee phospholipase A2 against the Ebola viral spike protein and its more lethal mutant strain.

Viral evolution sustains a dengue outbreak of enhanced severity

Compared to the previous 2013-2014 outbreak, dengue 2016-2017 outbreak in New Caledonia was characterized by an increased number of severe forms associated with hepatic presentations. In this study, we assessed the virological factors associated with this enhanced severity. Whole-genome sequences were retrieved from dengue virus (DENV)-1 strains collected in 2013-2014 and from severe and non-severe patients in 2016-2017. Fitness, hepatic tropism and cytopathogenicity of DENV 2016-2017 strains were compared to those of 2013-2014 strains using replication kinetics in the human hepatic cell line HuH7. Whole-genome sequencing identified four amino acid substitutions specific to 2016-2017 strains and absent from 2013-2014 strains. Three of these mutations occurred in predicted T cell epitopes, among which one was also a B cell epitope. Strains retrieved from severe forms did not exhibit specific genetic features. DENV strains from 2016-2017 exhibited a trend towards reduced replicative fitness and cytopathogenicity in vitro compared to strains from 2013-2014. Overall, the 2016-2017 dengue outbreak in New Caledonia was associated with a viral genetic evolution which had limited impact on DENV hepatic tropism and cytopathogenicity. These mutations, however, may have modified DENV strains antigenicity, altering the anti-DENV immune response in some patients, in turn favoring the development of severe forms.Trial registration: ClinicalTrials.gov identifier: NCT04615364.

Pathogenicity and Virulence of Ebolaviruses with Species- and Variant-specificity

Ebola virus (EBOV), belonging to the species Zaire ebolavirus in the genus Ebolavirus, causes a severe febrile illness in humans with case fatality rates (CFRs) up to 90%. While there have been six virus species classified, which each have a single type virus in the genus Ebolavirus, CFRs of ebolavirus infections vary among viruses belonging to each distinct species. In this review, we aim to define the ebolavirus species-specific virulence on the basis of currently available laboratory and experimental findings. In addition, this review will also cover the variant-specific virulence of EBOV by referring to the unique biological and pathogenic characteristics of EBOV variant Makona, a new EBOV variant isolated from the 2013-2016 EBOV disease outbreak in West Africa. A better definition of species-specific and variant-specific virulence of ebolaviruses will facilitate our comprehensive knowledge on genus Ebolavirus biology, leading to the development of therapeutics against well-focused pathogenic mechanisms of each Ebola disease.

[The latest research findings on Ebola virus]

013-2016 Ebola virus disease (EVD) outbreak was the largest EVD outbreak ever documented that started earlier in Guinea and later widely spread throughout West Africa, ending up a total of > 28,000 human infections. In this review, we outline research findings on Ebola virus (EBOV) variant Makona, a new EBOV variant isolated from the 2013-2016 EVD outbreak, and introduce the unique biological and pathogenic characteristics of Makona variant. We also discuss about the relevance of persistent infection of EBOV in EVD survivors with resurgence of EVD outbreak in Guinea in 2021. Moreover, this review covers a recent case report of EVD relapse and deliberates new interpretations of EBOV biology and EVD outbreak.

A Systematic Evaluation of High-Throughput Sequencing Approaches to Identify Low-Frequency Single Nucleotide Variants in Viral Populations

High-throughput sequencing such as those provided by Illumina are an efficient way to understand sequence variation within viral populations. However, challenges exist in distinguishing process-introduced error from biological variance, which significantly impacts our ability to identify sub-consensus single-nucleotide variants (SNVs). Here we have taken a systematic approach to evaluate laboratory and bioinformatic pipelines to accurately identify low-frequency SNVs in viral populations. Artificial DNA and RNA "populations" were created by introducing known SNVs at predetermined frequencies into template nucleic acid before being sequenced on an Illumina MiSeq platform. These were used to assess the effects of abundance and starting input material type, technical replicates, read length and quality, short-read aligner, and percentage frequency thresholds on the ability to accurately call variants. Analyses revealed that the abundance and type of input nucleic acid had the greatest impact on the accuracy of SNV calling as measured by a micro-averaged Matthews correlation coefficient score, with DNA and high RNA inputs (107 copies) allowing for variants to be called at a 0.2% frequency. Reduced input RNA (105 copies) required more technical replicates to maintain accuracy, while low RNA inputs (103 copies) suffered from consensus-level errors. Base errors identified at specific motifs identified in all technical replicates were also identified which can be excluded to further increase SNV calling accuracy. These findings indicate that samples with low RNA inputs should be excluded for SNV calling and reinforce the importance of optimising the technical and bioinformatics steps in pipelines that are used to accurately identify sequence variants.

Discovery of sclareol and sclareolide as filovirus entry inhibitors

Filoviruses cause severe hemorrhagic fever in humans. Ebola virus (EBOV) is the most contagious filovirus. Although compassionate treatments have been used during the latest Ebola outbreak, novel anti-EBOV agents are still urgently needed. In this study, sclareol and sclareolide, two natural products in Salvia sclarea, were identified as EBOV entry inhibitors with EC₅₀s of 2.4 μmol/L and 8.0 μmol/L, respectively, through blocking the viral fusion process. Moreover, both compounds exhibited inhibitory effects on all tested filoviruses' entry, indicating their wide-spectrum activities against filoviruses. This study provides insights into the two natural products and their applications against filovirus infectious diseases.

Identification and molecular characterization of the first complete genome sequence of Human Parechovirus type 15

Using a metagenomics approach, we have determined the first full-length genome sequence of a human parechovirus type 15 (HPeV15) strain, isolated from a child with acute flaccid paralysis and co-infected with EV-A71. HPeV15 is a rarely reported type. To date, no full-length genome sequence of HPeV15 is available in the GenBank database, where only limited VP1 sequences of this virus are available. Pairwise comparisons of the complete VP1 nucleotide and deduced amino acid sequences revealed that the study strain belongs to type 15 as it displayed 79.6% nucleotide and 93.4% amino acid identity with the HPeV15 prototype strain. Comparative analysis of available genomic regions and phylogenetic analysis using the P2 and P3 coding regions revealed low nucleotide identity to HPeV reference genomes. Phylogenetic and similarity plot analyses showed that genomic recombination events might have occurred in the UTRs and nonstructural region during HPeV15 evolution. The study strain has high similarity features with different variants of HPeV3 suggesting intertypic recombination. Our data contributes to the scarce data available on HPeVs in Africa and provides valuable information for future studies that aim to understand the evolutionary history, molecular epidemiology or biological and pathogenic properties of HPeV15.

Mucin-Like Domain of Ebola Virus Glycoprotein Enhances Selective Oncolytic Actions against Brain Tumors

Given that the Ebola virus (EBOV) infects a wide array of organs and cells yet displays a relative lack of neurotropism, we asked whether a chimeric vesicular stomatitis virus (VSV) expressing the EBOV glycoprotein (GP) might selectively target brain tumors. The mucin-like domain (MLD) of the EBOV GP may enhance virus immune system evasion. Here, we compared chimeric VSVs in which EBOV GP replaces the VSV glycoprotein, thereby reducing the neurotoxicity associated with wild-type VSV. A chimeric VSV expressing the full-length EBOV GP (VSV-EBOV) containing the MLD was substantially more effective and safer than a parallel construct with an EBOV GP lacking the MLD (VSV-EBOVΔMLD). One-step growth, reverse transcription-quantitative PCR, and Western blotting assessments showed that VSV-EBOVΔMLD produced substantially more progeny faster than VSV-EBOV. Using immunodeficient SCID mice, we focused on targeting human brain tumors with these VSV-EBOVs. Similar to the findings of our previous study in which we used an attenuated VSV-EBOV with no MLD that expressed green fluorescent protein (GFP) (VSV-EBOVΔMLD-GFP), VSV-EBOVΔMLD without GFP targeted glioma but yielded only a modest extension of survival. In contrast, VSV-EBOV containing the MLD showed substantially better targeting and elimination of brain tumors after intravenous delivery and increased the survival of brain tumor-bearing mice. Despite the apparent destruction of most tumor cells by VSV-EBOVΔMLD, the virus remained active within the SCID mouse brain and showed widespread infection of normal brain cells. In contrast, VSV-EBOV eliminated the tumors and showed relatively little infection of normal brain cells. Parallel experiments with direct intracranial virus infection generated similar results. Neither VSV-EBOV nor VSV-EBOVΔMLD showed substantive infection of the brains of normal immunocompetent mice.IMPORTANCE The Ebola virus glycoprotein contains a mucin-like domain which may play a role in immune evasion. Chimeric vesicular stomatitis viruses with the EBOV glycoprotein substituted for the VSV glycoprotein show greater safety and efficacy in targeting brain tumors in immunodeficient mice when the MLD was expressed within the EBOV glycoprotein than when EBOV lacked the mucin-like domain.

Standardized phylogenetic and molecular evolutionary analysis applied to species across the microbial tree of life

There is growing interest in reconstructing phylogenies from the copious amounts of genome sequencing projects that target related viral, bacterial or eukaryotic organisms. To facilitate the construction of standardized and robust phylogenies for disparate types of projects, we have developed a complete bioinformatic workflow, with a web-based component to perform phylogenetic and molecular evolutionary (PhaME) analysis from sequencing reads, draft assemblies or completed genomes of closely related organisms. Furthermore, the ability to incorporate raw data, including some metagenomic samples containing a target organism (e.g. from clinical samples with suspected infectious agents), shows promise for the rapid phylogenetic characterization of organisms within complex samples without the need for prior assembly.

Reporter Assays for Ebola Virus Nucleoprotein Oligomerization, Virion-Like Particle Budding, and Minigenome Activity Reveal the Importance of Nucleoprotein Amino Acid Position 111

For highly pathogenic viruses, reporter assays that can be rapidly performed are critically needed to identify potentially functional mutations for further study under maximal containment (e.g., biosafety level 4 [BSL-4]). The Ebola virus nucleoprotein (NP) plays multiple essential roles during the viral life cycle, yet few tools exist to study the protein under BSL-2 or equivalent containment. Therefore, we adapted reporter assays to measure NP oligomerization and virion-like particle (VLP) production in live cells and further measured transcription and replication using established minigenome assays. As a proof-of-concept, we examined the NP-R111C substitution, which emerged during the 2013‒2016 Western African Ebola virus disease epidemic and rose to high frequency. NP-R111C slightly increased NP oligomerization and VLP budding but slightly decreased transcription and replication. By contrast, a synthetic charge-reversal mutant, NP-R111E, greatly increased oligomerization but abrogated transcription and replication. These results are intriguing in light of recent structures of NP oligomers, which reveal that the neighboring residue, K110, forms a salt bridge with E349 on adjacent NP molecules. By developing and utilizing multiple reporter assays, we find that the NP-111 position mediates a complex interplay between NP's roles in protein structure, virion budding, and transcription and replication.

The ability of single genes vs full genomes to resolve time and space in outbreak analysis

Background

Inexpensive pathogen genome sequencing has had a transformative effect on the field of phylodynamics, where ever increasing volumes of data have promised real-time insight into outbreaks of infectious disease. As well as the sheer volume of pathogen isolates being sequenced, the sequencing of whole pathogen genomes, rather than select loci, has allowed phylogenetic analyses to be carried out at finer time scales, often approaching serial intervals for infections caused by rapidly evolving RNA viruses. Despite its utility, whole genome sequencing of pathogens has not been adopted universally and targeted sequencing of loci is common in some pathogen-specific fields.

Results

In this study we highlighted the utility of sequencing whole genomes of pathogens by re-analysing a well-characterised collection of Ebola virus sequences in the form of complete viral genomes (≈19 kb long) or the rapidly evolving glycoprotein (GP, ≈2 kb long) gene. We have quantified changes in phylogenetic, temporal, and spatial inference resolution as a result of this reduction in data and compared these to theoretical expectations.

Conclusions

We propose a simple intuitive metric for quantifying temporal resolution, i.e. the time scale over which sequence data might be informative of various processes as a quick back-of-the-envelope calculation of statistical power available to molecular clock analyses.

Recently Identified Mutations in the Ebola Virus-Makona Genome Do Not Alter Pathogenicity in Animal Models

Ebola virus (EBOV), isolate Makona, the causative agent of the West African EBOV epidemic, has been the subject of numerous investigations to determine the genetic diversity and its potential implication for virus biology, pathogenicity, and transmissibility. Despite various mutations that have emerged over time through multiple human-to-human transmission chains, their biological relevance remains questionable. Recently, mutations in the glycoprotein GP and polymerase L, which emerged and stabilized early during the outbreak, have been associated with improved viral fitness in cell culture. Here, we infected mice and rhesus macaques with EBOV-Makona isolates carrying or lacking those mutations. Surprisingly, all isolates behaved very similarly independent of the genotype, causing severe or lethal disease in mice and macaques, respectively. Likewise, we could not detect any evidence for differences in virus shedding. Thus, no specific biological phenotype could be associated with these EBOV-Makona mutations in two animal models.

Inferring time-dependent migration and coalescence patterns from genetic sequence and predictor data in structured populations

Population dynamics can be inferred from genetic sequence data by using phylodynamic methods. These methods typically quantify the dynamics in unstructured populations or assume migration rates and effective population sizes to be constant through time in structured populations. When considering rates to vary through time in structured populations, the number of parameters to infer increases rapidly and the available data might not be sufficient to inform these. Additionally, it is often of interest to know what predicts these parameters rather than knowing the parameters themselves. Here, we introduce a method to  infer the predictors for time-varying migration rates and effective population sizes by using a generalized linear model (GLM) approach under the marginal approximation of the structured coalescent. Using simulations, we show that our approach is able to reliably infer the model parameters and its predictors from phylogenetic trees. Furthermore, when simulating trees under the structured coalescent, we show that our new approach outperforms the discrete trait GLM model. We then apply our framework to a previously described Ebola virus dataset, where we infer the parameters and its predictors from genome sequences while accounting for phylogenetic uncertainty. We infer weekly cases to be the strongest predictor for effective population size and geographic distance the strongest predictor for migration. This approach is implemented as part of the BEAST2 package MASCOT, which allows us to jointly infer population dynamics, i.e. the parameters and predictors, within structured populations, the phylogenetic tree, and evolutionary parameters.

Active Ebola Virus Replication and Heterogeneous Evolutionary Rates in EVD Survivors

Following cessation of continuous Ebola virus (EBOV) transmission within Western Africa, sporadic EBOV disease (EVD) cases continued to re-emerge beyond the viral incubation period. Epidemiological and genomic evidence strongly suggests that this represented transmission from EVD survivors. To investigate whether persistent infections are characterized by ongoing viral replication, we sequenced EBOV from the semen of nine EVD survivors and a subset of corresponding acute specimens. EBOV evolutionary rates during persistence were either similar to or reduced relative to acute infection rates. Active EBOV replication/transcription continued during convalescence, but decreased over time, consistent with viral persistence rather than viral latency. Patterns of genetic divergence suggest a moderate relaxation of selective constraints within the sGP carboxy-terminal tail during persistent infections, but do not support widespread diversifying selection. Altogether, our data illustrate that EBOV persistence in semen, urine, and aqueous humor is not a quiescent or latent infection.

•

During persistence, EBOV exhibits heterogeneous evolutionary rates

•

Active EBOV transcription and replication occurs during persistence

•

RNA hyper-editing observed during viral persistence

•

No evidence for significant selective pressure during persistence

Longitudinal Analysis of the Human B Cell Response to Ebola Virus Infection

Ebola virus (EBOV) remains a public health threat. We performed a longitudinal study of B cell responses to EBOV in four survivors of the 2014 West African outbreak. Infection induced lasting EBOV-specific immunoglobulin G (IgG) antibodies, but their subclass composition changed over time, with IgG1 persisting, IgG3 rapidly declining, and IgG4 appearing late. Striking changes occurred in the immunoglobulin repertoire, with massive recruitment of naive B cells that subsequently underwent hypermutation. We characterized a large panel of EBOV glycoprotein-specific monoclonal antibodies (mAbs). Only a small subset of mAbs that bound glycoprotein by ELISA recognized cell-surface glycoprotein. However, this subset contained all neutralizing mAbs. Several mAbs protected against EBOV disease in animals, including one mAb that targeted an epitope under evolutionary selection during the 2014 outbreak. Convergent antibody evolution was seen across multiple donors, particularly among VH3-13 neutralizing antibodies specific for the GP1 core. Our study provides a benchmark for assessing EBOV vaccine-induced immunity.

•

Ebola virus infection causes massive recruitment of naive B cells

•

Virus-specific antibodies continue to class-switch and mutate for months after acute infection

•

Protective antibodies can be neutralizing or non-neutralizing and can appear early

•

Convergent, protective antibody rearrangements are seen in multiple donors

Intra-host Ebola viral adaption during human infection

The onsite next generation sequencing (NGS) of Ebola virus (EBOV) genomes during the 2013–2016 Ebola epidemic in Western Africa provides an opportunity to trace the origin, transmission, and evolution of this virus. Herein, we have diagnosed a cohort of EBOV patients in Sierra Leone in 2015, during the late phase of the outbreak. The surviving EBOV patients had a recovery process characterized by decreasing viremia, fever, and biochemical parameters. EBOV genomes sequenced through the longitudinal blood samples of these patients showed dynamic intra-host substitutions of the virus during acute infection, including the previously described short stretches of 13 serial T>C mutations. Remarkably, within individual patients, samples collected during the early phase of infection possessed Ts at these nucleotide sites, whereas they were replaced by Cs in samples collected in the later phase, suggesting that these short stretches of T>C mutations could emerge independently. In addition, up to a total of 35 nucleotide sites spanning the EBOV genome were mutated coincidently. Our study showed the dynamic intra-host adaptation of EBOV during patient recovery and gave more insight into the complex EBOV-host interactions.

Investigation of Outbreak-Specific Nonsynonymous Mutations on Ebolavirus GP in the Context of Known Immune Reactivity

The global response to the most recent EBOV outbreak has led to increased generation and availability of data, which can be globally analyzed to increase our understanding of immune responses to EBOV. We analyzed the published antibody epitope data to identify regions immunogenic for humans on the main GP antigenic target and determine sequence variance/nonsynonymous mutations between historical isolates and variants from the 2013-2016 outbreak. Approximately half of the GP sequence has been reported as targeted by antibody responses. Our results show an enrichment of nonsynonymous mutations (NSMs) within epitopic regions on GP (70%, p = 0.0133). Mapping NSMs to human epitope reactivity may be useful for future therapeutic and prophylaxis development as well as for our general understanding of immunity against EBOV.

Ebola virus, but not Marburg virus, replicates efficiently and without required adaptation in snake cells

Ebola virus (EBOV) disease is a viral hemorrhagic fever with a high case-fatality rate in humans. This disease is caused by four members of the filoviral genus Ebolavirus, including EBOV. The natural hosts reservoirs of ebolaviruses remain to be identified. Glycoprotein 2 of reptarenaviruses, known to infect only boa constrictors and pythons, is similar in sequence and structure to ebolaviral glycoprotein 2, suggesting that EBOV may be able to infect reptilian cells. Therefore, we serially passaged EBOV and a distantly related filovirus, Marburg virus (MARV), in boa constrictor JK cells and characterized viral infection/replication and mutational frequency by confocal imaging and sequencing. We observed that EBOV efficiently infected and replicated in JK cells, but MARV did not. In contrast to most cell lines, EBOV-infected JK cells did not result in an obvious cytopathic effect. Surprisingly, genomic characterization of serial-passaged EBOV in JK cells revealed that genomic adaptation was not required for infection. Deep sequencing coverage (>10,000×) demonstrated the existence of only a single nonsynonymous variant (EBOV glycoprotein precursor pre-GP T544I) of unknown significance within the viral population that exhibited a shift in frequency of at least 10 per cent over six serial passages. In summary, we present the first reptilian cell line that replicates a filovirus at high titers, and for the first time demonstrate a filovirus genus-specific restriction to MARV in a cell line. Our data suggest the possibility that there may be differences between the natural host spectra of ebolaviruses and marburgviruses.

Phylodynamic assessment of intervention strategies for the West African Ebola virus outbreak

Genetic analyses have provided important insights into Ebola virus spread during the recent West African outbreak, but their implications for specific intervention scenarios remain unclear. Here, we address this issue using a collection of phylodynamic approaches. We show that long-distance dispersal events were not crucial for epidemic expansion and that preventing viral lineage movement to any given administrative area would, in most cases, have had little impact. However, major urban areas were critical in attracting and disseminating the virus: preventing viral lineage movement to all three capitals simultaneously would have contained epidemic size to one-third. We also show that announcements of border closures were followed by a significant but transient effect on international virus dispersal. By quantifying the hypothetical impact of different intervention strategies, as well as the impact of barriers on dispersal frequency, our study illustrates how phylodynamic analyses can help to address specific epidemiological and outbreak control questions.

During the last Ebola virus outbreak in West Africa, a large amount of viral genomic data was obtained. Here, Dellicour et al. use phylodynamic approaches to assess effect of intervention strategies such as border closures.

Impact of the tree prior on estimating clock rates during epidemic outbreaks

Genetic sequencing data of pathogens allow one to quantify the evolutionary rate together with epidemiological dynamics using Bayesian phylodynamic methods. Such tools are particularly useful for obtaining a timely understanding of newly emerging epidemic outbreaks. During the West African Ebola virus disease epidemic, an unusually high evolutionary rate was initially estimated, promoting discussions regarding the potential danger of the strain quickly evolving into an even more dangerous virus. We show here that such high evolutionary rates are not necessarily real but can stem from methodological biases in the analyses. While most analyses of epidemic outbreak data are performed such that these biases may be present, we suggest a solution to overcome these biases in the future.

Bayesian phylogenetics aims at estimating phylogenetic trees together with evolutionary and population dynamic parameters based on genetic sequences. It has been noted that the clock rate, one of the evolutionary parameters, decreases with an increase in the sampling period of sequences. In particular, clock rates of epidemic outbreaks are often estimated to be higher compared with the long-term clock rate. Purifying selection has been suggested as a biological factor that contributes to this phenomenon, since it purges slightly deleterious mutations from a population over time. However, other factors such as methodological biases may also play a role and make a biological interpretation of results difficult. In this paper, we identify methodological biases originating from the choice of tree prior, that is, the model specifying epidemiological dynamics. With a simulation study we demonstrate that a misspecification of the tree prior can upwardly bias the inferred clock rate and that the interplay of the different models involved in the inference can be complex and nonintuitive. We also show that the choice of tree prior can influence the inference of clock rate on real-world Ebola virus (EBOV) datasets. While commonly used tree priors result in very high clock-rate estimates for sequences from the initial phase of the epidemic in Sierra Leone, tree priors allowing for population structure lead to estimates agreeing with the long-term rate for EBOV.

An all-atom molecular dynamics study of the anti-interferon signaling of Ebola virus: interaction mechanisms of EBOV VP24 binding to Karyopherin alpha5

Ebola virus (EBOV) is highly lethal due to virally encoded immune antagonists, and the combination of EBOV VP24 with karyopherin alpha (KPNA) will trigger anti-interferon (IFN) signaling. The crystal structure of VP24-KPNA5 has been proposed in recent studies, but the precise binding mechanisms are still unclear. In order to explore the VP24-KPNA5 protein binding micro-mechanisms, Molecular Dynamic (MD) simulations and Molecular Mechanics Generalized Born Surface Area (MM-GB/SA) energy calculation are performed. The obtained results show that EBOV VP24 binding to KPNA5 will rigidify their binding-face, and both proteins will be compacted during binding. According to the analyses of binding free energies of WT and the eight mutant systems, MUT3 makes the most effective contributions to the interaction; additionally MUT4, R398A and the double mutant have the second most effective influence. Hydrogen bond analysis demonstrates that inhibitors which can interfere with the formation of hydrogen bonds D480-T138, E483-R137 and D205-R396 will prevent the anti-IFN effect. Meanwhile, by combining the decomposition of binding free energies (DC) with computational alanine scanning (CAS) results, it is shown that VP24 residues R137 and T138 will be potential targets for EBOV VP24 inhibitors, and KPNA5 residues R396, R398, R480, Y477 and F484 will be potential targets to prevent KPNA5 binding to VP24, which will ultimately block anti-IFN signaling. Our investigations provide theoretical data to understand the binding modes of VP24-KPNA5. The precise binding mechanisms of the complex may shed light on the development of potential novel inhibitors against EBOV infection.

Experimental Adaptive Evolution of Simian Immunodeficiency Virus SIVcpz to Pandemic Human Immunodeficiency Virus Type 1 by Using a Humanized Mouse Model

Human immunodeficiency virus type 1 (HIV-1), the causative agent of AIDS, originated from simian immunodeficiency virus from chimpanzees (SIVcpz), the precursor of the human virus, approximately 100 years ago. This indicates that HIV-1 has emerged through the cross-species transmission of SIVcpz from chimpanzees to humans. However, it remains unclear how SIVcpz has evolved into pandemic HIV-1 in humans. To address this question, we inoculated three SIVcpz strains (MB897, EK505, and MT145), four pandemic HIV-1 strains (NL4-3, NLCSFV3, JRCSF, and AD8), and two nonpandemic HIV-1 strains (YBF30 and DJO0131). Humanized mice infected with SIVcpz strain MB897, a virus phylogenetically similar to pandemic HIV-1, exhibited a peak viral load comparable to that of mice infected with pandemic HIV-1, while peak viral loads of mice infected with SIVcpz strain EK505 or MT145 as well as nonpandemic HIV-1 strains were significantly lower. These results suggest that SIVcpz strain MB897 is preadapted to humans, unlike the other SIVcpz strains. Moreover, viral RNA sequencing of MB897-infected humanized mice identified a nonsynonymous mutation in env, a G413R substitution in gp120. The infectivity of the gp120 G413R mutant of MB897 was significantly higher than that of parental MB897. Furthermore, we demonstrated that the gp120 G413R mutant of MB897 augments the capacity for viral replication in both in vitro cell cultures and humanized mice. Taken together, this is the first experimental investigation to use an animal model to demonstrate a gain-of-function evolution of SIVcpz into pandemic HIV-1.IMPORTANCE From the mid-20th century, humans have been exposed to the menace of infectious viral diseases, such as severe acute respiratory syndrome coronavirus, Ebola virus, and Zika virus. These outbreaks of emerging/reemerging viruses can be triggered by cross-species viral transmission from wild animals to humans, or zoonoses. HIV-1, the causative agent of AIDS, emerged by the cross-species transmission of SIVcpz, the HIV-1 precursor in chimpanzees, around 100 years ago. However, the process by which SIVcpz evolved to become HIV-1 in humans remains unclear. Here, by using a hematopoietic stem cell-transplanted humanized-mouse model, we experimentally recapitulate the evolutionary process of SIVcpz to become HIV-1. We provide evidence suggesting that a strain of SIVcpz, MB897, preadapted to infect humans over other SIVcpz strains. We further demonstrate a gain-of-function evolution of SIVcpz in infected humanized mice. Our study reveals that pandemic HIV-1 has emerged through at least two steps: preadaptation and subsequent gain-of-function mutations.

Private collection: high correlation of sample collection and patient admission date in clinical microbiological testing complicates sharing of phylodynamic metadata

Infectious pathogens are known for their rapid evolutionary rates with new mutations arising over days to weeks. The ability to rapidly recover whole genome sequences and analyze the spread and evolution of pathogens using genetic information and pathogen collection dates has lead to interest in real-time tracking of infectious transmission and outbreaks. However, the level of temporal resolution afforded by these analyses may conflict with definitions of what constitutes protected health information (PHI) and privacy requirements for de-identification for publication and public sharing of research data and metadata. In the United States, dates and locations associated with patient care that provide greater resolution than year or the first three digits of the zip code are generally considered patient identifiers. Admission and discharge dates are specifically named as identifiers in Department of Health and Human Services guidance. To understand the degree to which one can impute admission dates from specimen collection dates, we examined sample collection dates and patient admission dates associated with more than 270,000 unique microbiological results from the University of Washington Laboratory Medicine Department between 2010 and 2017. Across all positive microbiological tests, the sample collection date exactly matched the patient admission date in 68.8% of tests. Collection dates and admission dates were identical from emergency department and outpatient testing 86.7% and 96.5% of the time, respectively, with >99% of tests collected within 1 day from the patient admission date. Samples from female patients were significantly more likely to be collected closer to admission date that those from male patients. We show that PHI-associated dates such as admission date can confidently be imputed from deposited collection date. We suggest that publicly depositing microbiological collection dates at greater resolution than the year may not meet routine Safe Harbor-based requirements for patient de-identification. We recommend the use of Expert Determination to determine PHI for a given study and/or direct patient consent if clinical laboratories or phylodynamic practitioners desire to make these data available.

A Single Amino Acid Change in the Marburg Virus Glycoprotein Arises during Serial Cell Culture Passages and Attenuates the Virus in a Macaque Model of Disease

Marburg virus (MARV) causes disease with high case fatality rates, and there are no approved vaccines or therapies. Licensing of MARV countermeasures will likely require approval via the FDA's Animal Efficacy Rule, which requires well-characterized animal models that recapitulate human disease. This includes selection of the virus used for exposure and ensuring that it retains the properties of the original isolate. The consequences of amplification of MARV for challenge studies are unknown. Here, we serially passaged and characterized MARV through 13 passes from the original isolate. Surprisingly, the viral genome was very stable, except for a single nucleotide change that resulted in an amino acid substitution in the hydrophobic region of the signal peptide of the glycoprotein (GP). The particle/PFU ratio also decreased following passages, suggesting a role for the amino acid in viral infectivity. To determine if amplification introduces a phenotype in an animal model, cynomolgus macaques were exposed to either 100 or 0.01 PFU of low- and high-passage-number MARV. All animals succumbed when exposed to 100 PFU of either passage 3 or 13 viruses, although animals exposed to the high-passage-number virus survived longer. However, none of the passage 13 MARV-exposed animals succumbed to 0.01-PFU exposure compared to 75% of passage 3-exposed animals. This is consistent with other filovirus studies that show some particles that are unable to yield a plaque in cell culture can cause lethal disease in vivo. These results have important consequences for the design of experiments that investigate MARV pathogenesis and that test the efficacy of MARV countermeasures. IMPORTANCE Marburg virus (MARV) causes disease with a high case fatality rate, and there are no approved vaccines or therapies. Serial amplification of viruses in cell culture often results in accumulation of mutations, but the effect of such cell culture passage on MARV is unclear. Serial passages of MARV resulted in a single mutation in the region encoding the glycoprotein (GP). This is a region where mutations can have important consequences on outbreaks and human disease [S. Mahanty and M. Bray, Lancet Infect Dis 4:487-498, 2004, https://doi.org/10.1016/S1473-3099(04)01103-X]. We thus investigated whether this mutation impacted disease by using a cynomolgus macaque model of MARV infection. Monkeys exposed to virus containing the mutation had better clinical outcomes than monkeys exposed to virus without the mutation. We also observed that a remarkably low number of MARV particles was sufficient to cause death. Our results could have a significant impact on how future studies are designed to model MARV disease and test vaccines and therapeutics.

The Ebola outbreak, 2013-2016: old lessons for new epidemics

Ebola virus causes a severe haemorrhagic fever in humans with high case fatality and significant epidemic potential. The 2013-2016 outbreak in West Africa was unprecedented in scale, being larger than all previous outbreaks combined, with 28 646 reported cases and 11 323 reported deaths. It was also unique in its geographical distribution and multicountry spread. It is vital that the lessons learned from the world's largest Ebola outbreak are not lost. This article aims to provide a detailed description of the evolution of the outbreak. We contextualize this outbreak in relation to previous Ebola outbreaks and outline the theories regarding its origins and emergence. The outbreak is described by country, in chronological order, including epidemiological parameters and implementation of outbreak containment strategies. We then summarize the factors that led to rapid and extensive propagation, as well as highlight the key successes, failures and lessons learned from this outbreak and the response.This article is part of the themed issue 'The 2013-2016 West African Ebola epidemic: data, decision-making and disease control'.

Evolutionary Analysis of HIV-1 Pol Proteins Reveals Representative Residues for Viral Subtype Differentiation

RNA viruses have been used as model systems to understand the patterns and processes of molecular evolution because they have high mutation rates and are genetically diverse. Human immunodeficiency virus 1 (HIV-1), the etiological agent of acquired immune deficiency syndrome, is highly genetically diverse, and is classified into several groups and subtypes. However, it has been difficult to use its diverse sequences to establish the overall phylogenetic relationships of different strains or the trends in sequence conservation with the construction of phylogenetic trees. Our aims were to systematically characterize HIV-1 subtype evolution and to identify the regions responsible for HIV-1 subtype differentiation at the amino acid level in the Pol protein, which is often used to classify the HIV-1 subtypes. In this study, we systematically characterized the mutation sites in 2,052 Pol proteins from HIV-1 group M (144 subtype A; 1,528 subtype B; 380 subtype C), using sequence similarity networks. We also used spectral clustering to group the sequences based on the network graph structures. A stepwise analysis of the cluster hierarchies allowed us to estimate a possible evolutionary pathway for the Pol proteins. The subtype A sequences also clustered according to when and where the viruses were isolated, whereas both the subtype B and C sequences remained as single clusters. Because the Pol protein has several functional domains, we identified the regions that are discriminative by comparing the structures of the domain-based networks. Our results suggest that sequence changes in the RNase H domain and the reverse transcriptase (RT) connection domain are responsible for the subtype classification. By analyzing the different amino acid compositions at each site in both domain sequences, we found that a few specific amino acid residues (i.e., M357 in the RT connection domain and Q480, Y483, and L491 in the RNase H domain) represent the differences among the subtypes. These residues were located on the surface of the RT structure and in the vicinity of the amino acid sites responsible for RT enzymatic activity or function.

Mapping the clinical outcomes and genetic evolution of Ebola virus in Sierra Leone

Sierra Leone was the most severely affected country in Western Africa during the 2013-2016 outbreak of Ebola virus disease (EVD). Previous genome surveillance studies have revealed the origin, diversity, and evolutionary dynamics of the Ebola virus (EBOV); however, the information regarding EBOV sequences is insufficient, especially the clinical outcomes, given that the correlation between the clinical outcomes and the genetic evolution of EBOV is still not clear. Here, we collected and curated a comprehensive data set that includes 514 EBOV genome sequences from patients with confirmed EVD (including 60 sequences not previously studied), >87.5% of which have residence information and definitive clinical outcomes. Phylogenetic reconstruction revealed 11 lineages of EBOV in Sierra Leone. The median-joining haplotype network showed that haplotypes that are associated with lethal outcomes tend to contribute more to the spread of the EBOV in Sierra Leone than those with live outcomes. Analyses of the spatial-temporal distribution unraveled the lineage-distinctive distribution patterns. Different viral lineages have different case fatality rates (CFRs) during the same stage of the outbreak, implying that several lineages featuring SNPs may correlate with increased/decreased CFRs. This study provides invaluable data sets of EBOV infection and highlights the potential SNPs for further in-depth investigation.

Flavivirus and Filovirus EvoPrinters: New alignment tools for the comparative analysis of viral evolution

Background

Flavivirus and Filovirus infections are serious epidemic threats to human populations. Multi-genome comparative analysis of these evolving pathogens affords a view of their essential, conserved sequence elements as well as progressive evolutionary changes. While phylogenetic analysis has yielded important insights, the growing number of available genomic sequences makes comparisons between hundreds of viral strains challenging. We report here a new approach for the comparative analysis of these hemorrhagic fever viruses that can superimpose an unlimited number of one-on-one alignments to identify important features within genomes of interest.

Methodology/Principal finding

We have adapted EvoPrinter alignment algorithms for the rapid comparative analysis of Flavivirus or Filovirus sequences including Zika and Ebola strains. The user can input a full genome or partial viral sequence and then view either individual comparisons or generate color-coded readouts that superimpose hundreds of one-on-one alignments to identify unique or shared identity SNPs that reveal ancestral relationships between strains. The user can also opt to select a database genome in order to access a library of pre-aligned genomes of either 1,094 Flaviviruses or 460 Filoviruses for rapid comparative analysis with all database entries or a select subset. Using EvoPrinter search and alignment programs, we show the following: 1) superimposing alignment data from many related strains identifies lineage identity SNPs, which enable the assessment of sublineage complexity within viral outbreaks; 2) whole-genome SNP profile screens uncover novel Dengue2 and Zika recombinant strains and their parental lineages; 3) differential SNP profiling identifies host cell A-to-I hyper-editing within Ebola and Marburg viruses, and 4) hundreds of superimposed one-on-one Ebola genome alignments highlight ultra-conserved regulatory sequences, invariant amino acid codons and evolutionarily variable protein-encoding domains within a single genome.

Conclusions/Significance

EvoPrinter allows for the assessment of lineage complexity within Flavivirus or Filovirus outbreaks, identification of recombinant strains, highlights sequences that have undergone host cell A-to-I editing, and identifies unique input and database SNPs within highly conserved sequences. EvoPrinter’s ability to superimpose alignment data from hundreds of strains onto a single genome has allowed us to identify unique Zika virus sublineages that are currently spreading in South, Central and North America, the Caribbean, and in China. This new set of integrated alignment programs should serve as a useful addition to existing tools for the comparative analysis of these viruses.

Flaviviruses, including Zika and Dengue viruses, and Filoviruses, including Ebola and Marburg viruses, are significant global public health threats. Genetic surveillance of viral isolates provides important insights into the origin of outbreaks, reveals lineage heterogeneity and diversification, and facilitates identification of novel recombinant strains and host cell modified viral genomes. We report the development of EvoPrinter, a web-accessed alignment tool for the rapid comparative analysis of viral genomes. EvoPrinter superimposes alignment data from multiple pairwise comparisons onto a single reference sequence of interest, to reveal both similarities and differences detected in hundreds of selected viral isolates. Evoprinter databases provide easy access to hundreds of non-redundant Flavivirus and Filovirus genomes. allowing the user to distinguish between sublineage identity SNPs and unique strain-specific SNPs, thus facilitating analysis of the history of viral diversification during an epidemic. EvoPrinter also proves useful in identifying recombinant strains and their parental lineages and detecting host-cell genomic editing. EvoPrinter should serve as a useful addition to existing tools for the comparative analysis of these viruses.

Mathematical models for devising the optimal Ebola virus disease eradication

BACKGROUND

The 2014-2015 epidemic of Ebola virus disease (EVD) in West Africa defines an unprecedented health threat for human.

METHODS

We construct a mathematical model to devise the optimal Ebola virus disease eradication plan. We used mathematical model to investigate the numerical spread of Ebola and eradication pathways, further fit our model against the real total cases data and calculated infection rate as 1.754.

RESULTS

With incorporating hospital isolation and application of medication in our model and analyzing their effect on resisting the spread, we demonstrate the second peak of 10,029 total cases in 23 days, and expect to eradicate EVD in 285 days. Using the regional spread of EVD with our transmission model analysis, we analyzed the numbers of new infections through four important transmission paths including household, community, hospital and unsafe funeral.

CONCLUSIONS

Based on the result of the model, we find out the key paths in different situations and propose our suggestion to control regional transmission. We fully considers Ebola characteristics, economic and time optimization, dynamic factors and local condition constraints, and to make our plan realistic, sensible and feasible.

Virus genomes reveal factors that spread and sustained the Ebola epidemic

The 2013-2016 West African epidemic caused by the Ebola virus was of unprecedented magnitude, duration and impact. Here we reconstruct the dispersal, proliferation and decline of Ebola virus throughout the region by analysing 1,610 Ebola virus genomes, which represent over 5% of the known cases. We test the association of geography, climate and demography with viral movement among administrative regions, inferring a classic 'gravity' model, with intense dispersal between larger and closer populations. Despite attenuation of international dispersal after border closures, cross-border transmission had already sown the seeds for an international epidemic, rendering these measures ineffective at curbing the epidemic. We address why the epidemic did not spread into neighbouring countries, showing that these countries were susceptible to substantial outbreaks but at lower risk of introductions. Finally, we reveal that this large epidemic was a heterogeneous and spatially dissociated collection of transmission clusters of varying size, duration and connectivity. These insights will help to inform interventions in future epidemics.

A Polymorphism within the Internal Fusion Loop of the Ebola Virus Glycoprotein Modulates Host Cell Entry

The large scale of the Ebola virus disease (EVD) outbreak in West Africa in 2013-2016 raised the question whether the host cell interactions of the responsible Ebola virus (EBOV) strain differed from those of other ebolaviruses. We previously reported that the glycoprotein (GP) of the virus circulating in West Africa in 2014 (EBOV2014) exhibited reduced ability to mediate entry into two nonhuman primate (NHP)-derived cell lines relative to the GP of EBOV1976. Here, we investigated the molecular determinants underlying the differential entry efficiency. We found that EBOV2014-GP-driven entry into diverse NHP-derived cell lines, as well as human monocyte-derived macrophages and dendritic cells, was reduced compared to EBOV1976-GP, although entry into most human- and all bat-derived cell lines tested was comparable. Moreover, EBOV2014 replication in NHP but not human cells was diminished relative to EBOV1976, suggesting that reduced cell entry translated into reduced viral spread. Mutagenic analysis of EBOV2014-GP and EBOV1976-GP revealed that an amino acid polymorphism in the receptor-binding domain, A82V, modulated entry efficiency in a cell line-independent manner and did not account for the reduced EBOV2014-GP-driven entry into NHP cells. In contrast, polymorphism T544I, located in the internal fusion loop in the GP2 subunit, was found to be responsible for the entry phenotype. These results suggest that position 544 is an important determinant of EBOV infectivity for both NHP and certain human target cells.IMPORTANCE The Ebola virus disease outbreak in West Africa in 2013 entailed more than 10,000 deaths. The scale of the outbreak and its dramatic impact on human health raised the question whether the responsible virus was particularly adept at infecting human cells. Our study shows that an amino acid exchange, A82V, that was acquired during the epidemic and that was not observed in previously circulating viruses, increases viral entry into diverse target cells. In contrast, the epidemic virus showed a reduced ability to enter cells of nonhuman primates compared to the virus circulating in 1976, and a single amino acid exchange in the internal fusion loop of the viral glycoprotein was found to account for this phenotype.

Human Adaptation of Ebola Virus during the West African Outbreak

The 2013–2016 outbreak of Ebola virus (EBOV) in West Africa was the largest recorded. It began following the cross-species transmission of EBOV from an animal reservoir, most likely bats, into humans, with phylogenetic analysis revealing the co-circulation of several viral lineages. We hypothesized that this prolonged human circulation led to genomic changes that increased viral transmissibility in humans. We generated a synthetic glycoprotein (GP) construct based on the earliest reported isolate and introduced amino acid substitutions that defined viral lineages. Mutant GPs were used to generate a panel of pseudoviruses, which were used to infect different human and bat cell lines. These data revealed that specific amino acid substitutions in the EBOV GP have increased tropism for human cells, while reducing tropism for bat cells. Such increased infectivity may have enhanced the ability of EBOV to transmit among humans and contributed to the wide geographic distribution of some viral lineages.

•

EBOV adapted to humans during the West African outbreak

•

Amino acid substitutions in the EBOV glycoprotein increase human cell tropism

•

The same glycoprotein amino acid substitutions decrease tropism for bat cells

Between Securitisation and Neglect: Managing Ebola at the Borders of Global Health

In 2014 the World Health Organization (WHO) was widely criticised for failing to anticipate that an outbreak of Ebola in a remote forested region of south-eastern Guinea would trigger a public health emergency of international concern (pheic). In explaining the WHO's failure, critics have pointed to structural restraints on the United Nations organisation and a leadership 'vacuum' in Geneva, among other factors. This paper takes a different approach. Drawing on internal WHO documents and interviews with key actors in the epidemic response, I argue that the WHO's failure is better understood as a consequence of Ebola's shifting medical identity and of triage systems for managing emerging infectious disease (EID) risks. Focusing on the discursive and non-discursive practices that produced Ebola as a 'problem' for global health security, I argue that by 2014 Ebola was no longer regarded as a paradigmatic EID and potential biothreat so much as a neglected tropical disease. The result was to relegate Ebola to the fringes of biosecurity concerns just at the moment when the virus was crossing international borders in West Africa and triggering large urban outbreaks for the first time. Ebola's fluctuating medical identity also helps explain the prominence of fear and rumours during the epidemic and social resistance to Ebola control measures. Contrasting the WHO's delay over declaring a pheic in 2014, with its rapid declaration of pheics in relation to H1N1 swine flu in 2009 and polio in 2014, I conclude that such 'missed alarms' may be an inescapable consequence of pandemic preparedness systems that seek to rationalise responses to the emergence of new diseases.

Biochemical Basis for Increased Activity of Ebola Glycoprotein in the 2013-16 Epidemic

Ebola virus (EBOV) infection is characterized by sporadic outbreaks caused by zoonotic transmission. Fixed changes in amino acid sequence, such as A82V in the EBOV glycoprotein (GP) that occurred early in the 2013-16 epidemic, are suspected to confer a selective advantage to the virus. We used biochemical assays of GP function to show that A82V, as well as a polymorphism in residue 544 identified in other outbreaks, enhances infection by decreasing the threshold for activation of membrane fusion activity triggered by the host factors cathepsin B and Niemann-Pick C1. Importantly, the increase in infectivity comes with the cost of decreased virus stability. Thus, emergence of a virus GP with altered properties that can affect transmission and virulence may have contributed to the severity and scope of the 2013-16 EBOV epidemic.

Functional Characterization of Adaptive Mutations during the West African Ebola Virus Outbreak

The Ebola virus (EBOV) outbreak in West Africa started in December 2013, claimed more than 11,000 lives, threatened to destabilize a whole region, and showed how easily health crises can turn into humanitarian disasters. EBOV genomic sequences of the West African outbreak revealed nonsynonymous mutations, which induced considerable public attention, but their role in virus spread and disease remains obscure. In this study, we investigated the functional significance of three nonsynonymous mutations that emerged early during the West African EBOV outbreak. Almost 90% of more than 1,000 EBOV genomes sequenced during the outbreak carried the signature of three mutations: a D759G substitution in the active center of the L polymerase, an A82V substitution in the receptor binding domain of surface glycoprotein GP, and an R111C substitution in the self-assembly domain of RNA-encapsidating nucleoprotein NP. Using a newly developed virus-like particle system and reverse genetics, we found that the mutations have an impact on the functions of the respective viral proteins and on the growth of recombinant EBOVs. The mutation in L increased viral transcription and replication, whereas the mutation in NP decreased viral transcription and replication. The mutation in the receptor binding domain of the glycoprotein GP improved the efficiency of GP-mediated viral entry into target cells. Recombinant EBOVs with combinations of the three mutations showed a growth advantage over the prototype isolate Makona C7 lacking the mutations. This study showed that virus variants with improved fitness emerged early during the West African EBOV outbreak.

IMPORTANCE

The dimension of the Ebola virus outbreak in West Africa was unprecedented. Amino acid substitutions in the viral L polymerase, surface glycoprotein GP, and nucleocapsid protein NP emerged, were fixed early in the outbreak, and were found in almost 90% of the sequences. Here we showed that these mutations affected the functional activity of viral proteins and improved viral growth in cell culture. Our results demonstrate emergence of adaptive changes in the Ebola virus genome during virus circulation in humans and prompt further studies on the potential role of these changes in virus transmissibility and pathogenicity.

West Africa 2013 Ebola: From Virus Outbreak to Humanitarian Crisis

The 2013 outbreak of Ebola virus disease (EVD) in West Africa constituted a major humanitarian crisis. The outbreak numbered over 28,500 cases, more than 10 times the number cumulatively registered from all previous EVD outbreaks combined, with at least 11,000 deaths, and resulted in billions of dollars of lost economic growth to an already impoverished region. The unprecedented scale of West Africa 2013 took the world by surprise and laid bare deficiencies in our response capacity to complex humanitarian disasters of highly infectious and lethal pathogens. However, the magnitude of West Africa 2013 also provided a unique opportunity and obligation to better understand not only the biology and epidemiology of EVD, but also the many scientific, economic, social, political, ethical, and logistical challenges in confronting emerging infectious diseases in the modern era.

Understanding the complex evolution of rapidly mutating viruses with deep sequencing: Beyond the analysis of viral diversity

With the advent of affordable deep sequencing technologies, detection of low frequency variants within genetically diverse viral populations can now be achieved with unprecedented depth and efficiency. The high-resolution data provided by next generation sequencing technologies is currently recognised as the gold standard in estimation of viral diversity. In the analysis of rapidly mutating viruses, longitudinal deep sequencing datasets from viral genomes during individual infection episodes, as well as at the epidemiological level during outbreaks, now allow for more sophisticated analyses such as statistical estimates of the impact of complex mutation patterns on the evolution of the viral populations both within and between hosts. These analyses are revealing more accurate descriptions of the evolutionary dynamics that underpin the rapid adaptation of these viruses to the host response, and to drug therapies. This review assesses recent developments in methods and provide informative research examples using deep sequencing data generated from rapidly mutating viruses infecting humans, particularly hepatitis C virus (HCV), human immunodeficiency virus (HIV), Ebola virus and influenza virus, to understand the evolution of viral genomes and to explore the relationship between viral mutations and the host adaptive immune response. Finally, we discuss limitations in current technologies, and future directions that take advantage of publically available large deep sequencing datasets.

Ebola Virus Glycoprotein with Increased Infectivity Dominated the 2013-2016 Epidemic

The magnitude of the 2013-2016 Ebola virus disease (EVD) epidemic enabled an unprecedented number of viral mutations to occur over successive human-to-human transmission events, increasing the probability that adaptation to the human host occurred during the outbreak. We investigated one nonsynonymous mutation, Ebola virus (EBOV) glycoprotein (GP) mutant A82V, for its effect on viral infectivity. This mutation, located at the NPC1-binding site on EBOV GP, occurred early in the 2013-2016 outbreak and rose to high frequency. We found that GP-A82V had heightened ability to infect primate cells, including human dendritic cells. The increased infectivity was restricted to cells that have primate-specific NPC1 sequences at the EBOV interface, suggesting that this mutation was indeed an adaptation to the human host. GP-A82V was associated with increased mortality, consistent with the hypothesis that the heightened intrinsic infectivity of GP-A82V contributed to disease severity during the EVD epidemic.

The 2014 Ebola virus outbreak in West Africa highlights no evidence of rapid evolution or adaptation to humans

Following its immergence in December 2013, the recent Zaire Ebola virus (EBOV) outbreak in West Africa has spread and persisted for more than two years, making it the largest EBOV epidemic in both scale and geographical region to date. In this study, a total of 726 glycoprotein (GP) gene sequences of the EBOV full-length genome obtained from West Africa from the 2014 outbreak, combined with 30 from earlier outbreaks between 1976 and 2008 were used to investigate the genetic divergence, evolutionary history, population dynamics, and selection pressure of EBOV among distinct epidemic waves. Results from our dataset showed that no non-synonymous substitutions occurred on the GP gene coding sequences of EBOV that were likely to have affected protein structure or function in any way. Furthermore, the significantly different dN/dS ratios observed between the 2014 West African outbreak and earlier outbreaks were more likely due to the confounding presence of segregating polymorphisms. Our results highlight no robust evidence that the 2014 EBOV outbreak is fast-evolving and adapting to humans. Therefore, the unprecedented nature of the 2014 EBOV outbreak might be more likely related to non-virological elements, such as environmental and sociological factors.

The evolution of Ebola virus: Insights from the 2013-2016 epidemic

The 2013-2016 epidemic of Ebola virus disease in West Africa was of unprecedented magnitude and changed our perspective on this lethal but sporadically emerging virus. This outbreak also marked the beginning of large-scale real-time molecular epidemiology. Here, we show how evolutionary analyses of Ebola virus genome sequences provided key insights into virus origins, evolution and spread during the epidemic. We provide basic scientists, epidemiologists, medical practitioners and other outbreak responders with an enhanced understanding of the utility and limitations of pathogen genomic sequencing. This will be crucially important in our attempts to track and control future infectious disease outbreaks.

Understanding Ebola: the 2014 epidemic

Near the end of 2013, an outbreak of Zaire ebolavirus (EBOV) began in Guinea, subsequently spreading to neighboring Liberia and Sierra Leone. As this epidemic grew, important public health questions emerged about how and why this outbreak was so different from previous episodes. This review provides a synthetic synopsis of the 2014-15 outbreak, with the aim of understanding its unprecedented spread. We present a summary of the history of previous epidemics, describe the structure and genetics of the ebolavirus, and review our current understanding of viral vectors and the latest treatment practices. We conclude with an analysis of the public health challenges epidemic responders faced and some of the lessons that could be applied to future outbreaks of Ebola or other viruses.