Genetic diversity and evolutionary dynamics of Ebola virus in Sierra Leone

Correlation of fluid-attenuated inversion recovery sequence vascular hyperintensity in magnetic resonance with collateral circulation and short-term clinical prognosis in acute ischemic stroke

Background

Accurately assessing the prognosis of patient with large-scale cerebral infarction caused by acute middle cerebral artery (MCA) occlusion in the early stages of onset can help clinicians to actively and effectively intervene, thus reducing mortality and disability rates. This study set out to investigate the predictive value of fluid-attenuated inversion recovery vascular hyperintensity (FVH) on collateral circulation and clinical prognosis.

Methods

The clinical data of 70 patients admitted to The First People's Hospital of Lianyungang from January 2018 to December 2021 with acute cerebral infarction due to occlusion of the proximal end of the M1 segment in the MCA were retrospectively collected. All patients had their first onset of disease and did not receive thrombolytic therapy at the time of onset. Subsequently, they underwent endovascular thrombectomy for treatment. The FVH and collateral vessel scores were derived according to patients' fluid-attenuated in version recovery (FLAIR) sequence and time-of-flight magnetic resonance angiography images. Based on the 90-day Modified Rankin Scale (mRS), patients were allocated to a good prognosis group (mRS ≤2) and a poor prognosis group (mRS =3-6). The correlation between the FVH and collateral vessel scores was assessed using the Spearman rank correlation test. Pearson correlation coefficient analysis was used to assess the correlation between FVH and the 90-day mRS together with the infarct size. Univariate analysis, multivariate binary logistic regression analysis, and receiver operating characteristic (ROC) curve analysis were adopted to identify those factors potentially. associated with the prognosis of patients with acute ischemic stroke (AIS).

Results

Out of 70 patients with acute unilateral MCA occlusion (MCAO) who met the inclusion criteria, 62 showed positive FVH sign. These 62 patients were divided into a good prognosis group (n=32) and a poor prognosis group (n=30) based on the mRS score 90 days after discharge. The Spearman rank correlation test indicated that FVH was positively correlated with collateral vessel grade (Spearman rho =0.865; P<0.001); meanwhile, Pearson correlation coefficient analysis indicated that FVH score had moderate negative correlation with 90-day mRS score (r=-0.605; P<0.001). The results of multivariate binary logistic regression analysis indicated that collateral vessel grade and FVH score may be associated with the prognosis of patients with AIS, and the area under the curve (AUC) of FVH score was larger than collateral vessel grade (AUC =0.738).

Conclusions

There was a positive correlation between FVH score and collateral vessel grade, and FVH score could indicate collateral circulation. FVH score was negatively correlated with 90-day mRS score and infarct volume and thus can predict clinical prognosis.

Natural history of Ebola virus disease in rhesus monkeys shows viral variant emergence dynamics and tissue-specific host responses

Ebola virus (EBOV) causes Ebola virus disease (EVD), marked by severe hemorrhagic fever; however, the mechanisms underlying the disease remain unclear. To assess the molecular basis of EVD across time, we performed RNA sequencing on 17 tissues from a natural history study of 21 rhesus monkeys, developing new methods to characterize host-pathogen dynamics. We identified alterations in host gene expression with previously unknown tissue-specific changes, including downregulation of genes related to tissue connectivity. EBOV was widely disseminated throughout the body; using a new, broadly applicable deconvolution method, we found that viral load correlated with increased monocyte presence. Patterns of viral variation between tissues differentiated primary infections from compartmentalized infections, and several variants impacted viral fitness in a EBOV/Kikwit minigenome system, suggesting that functionally significant variants can emerge during early infection. This comprehensive portrait of host-pathogen dynamics in EVD illuminates new features of pathogenesis and establishes resources to study other emerging pathogens.

Emerging Infectious Diseases Are Virulent Viruses-Are We Prepared? An Overview

The recent pandemic caused by SARS-CoV-2 affected the global population, resulting in a significant loss of lives and global economic deterioration. COVID-19 highlighted the importance of public awareness and science-based decision making, and exposed global vulnerabilities in preparedness and response systems. Emerging and re-emerging viral outbreaks are becoming more frequent due to increased international travel and global warming. These viral outbreaks impose serious public health threats and have transformed national strategies for pandemic preparedness with global economic consequences. At the molecular level, viral mutations and variations are constantly thwarting vaccine efficacy, as well as diagnostic, therapeutic, and prevention strategies. Here, we discuss viral infectious diseases that were epidemic and pandemic, currently available treatments, and surveillance measures, along with their limitations.

Oligonucleotide usage in coronavirus genomes mimics that in exon regions in host genomes

BACKGROUND

Viruses use various host factors for their growth, and efficient growth requires efficient use of these factors. Our previous study revealed that the occurrence frequency of oligonucleotides in the influenza virus genome is distinctly different among derived hosts, and the frequency tends to adapt to the host cells in which they grow. We aimed to study the adaptation mechanisms of a zoonotic virus to host cells.

METHODS

Herein, we compared the frequency of oligonucleotides in the genome of alpha- and betacoronavirus with those in the genomes of humans and bats, which are typical hosts of the viruses.

RESULTS

By comparing the oligonucleotide frequency in coronaviruses and their host genomes, we found a statistically tested positive correlation between the frequency of coronaviruses and that of the exon regions of the host from which the virus is derived. To examine the characteristics of early-stage changes in the viral genome, which are assumed to accompany the host change from non-humans to humans, we compared the oligonucleotide frequency between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at the beginning of the pandemic and the prevalent variants thereafter, and found changes towards the frequency of the host exon regions.

CONCLUSIONS

In alpha- and betacoronaviruses, the genome oligonucleotide frequency is thought to change in response to the cellular environment in which the virus is replicating, and actually the frequency has approached the frequency in exon regions in the host.

Biological characteristics, biosafety prevention and control strategies for the 2022 multi-country outbreak of monkeypox

Monkeypox is a zoonotic disease caused by the monkeypox virus (MPXV), which is a potential biological warfare agent of bioterrorism and poses the greatest threat to the world's public biosafety and health after variola virus (VARV). While the coronavirus disease 2019 (COVID-19) pandemic has not ended yet, monkeypox is spreading menacingly. The first case of monkeypox in a nonendemic country was confirmed on May 6^th, 2022, while the first imported case from Asia was found on June 21^st. There were more than 16 thousand reported cases as of July 23^rd, the day the World Health Organization (WHO) declared the global monkeypox outbreak a public health emergency of international concern (PHEIC) at the same level as smallpox and COVID-19; while there were more than 53 thousand cases as of September 1^st. Therefore, we will propose relevant biosafety prevention and control strategies after analyzing the etiology of the 2022 multi-country monkeypox outbreak from the biological feature, transmissibility, epidemic, and variability of MPXV.

The diagnostic value of metagenomic next-generation sequencing for identifying Pneumocystis jirovecii infection in non-HIV immunocompromised patients

Background

Pneumocystis jirovecii pneumonia (PJP) remains an important cause of morbidity and mortality in non-HIV immunocompromised patients especially in transplant recipients. But its diagnosis remains challenging due to the insuffificient performance of conventional methods for diagnosing Pneumocystis jirovecii(P. jirovecii) infection. Therefore, the auxiliary diagnostic function of metagenomics next-generation sequencing (mNGS) in clinical practice is worth of exploring.

Method

34 non-HIV immunocompromised patients who were diagnosed as PJP by clinical manifestations, imaging findings, immune status of the host, and Methenamine silver staining were tested by mNGS from October 2018 to December 2020 in Sichuan Provincial People’s Hospital. The clinical performances of mNGS for P. jirovecii infection diagnosis were also evaluated with genome reads abundance and comparing with other traditional diagnostic methods.

Results

We diagnosed a total of 34 non-HIV PJP patients by the clinical composite diagnosis. Our data shows that, compared with the clinical microbiological test, the detection rate of mNGS for P. jirovecii in non-HIV infected PJP patients is significantly higher than that of Methenamine silver staining and serum 1-3-β-D-glucan. mNGS can be used as an auxiliary diagnostic tool to help diagnosis. The number of reads mapped to the genome of P. jirovecii and the duration of patients from onset to sampling collection were statistically significant between the two groups (Reads>100 and Reads ≤ 100) (8days vs. 23days, p=0.020). In addition, univariate analysis showed that C-reactive protein (15.8mg/L vs.79.56mg/L, p=0.016), lactate dehydrogenase (696U/l vs. 494U/l, p=0.030) and procalcitonin (0.09ng/ml vs. 0.59ng/ml, p=0.028) was also statistically significant between the two groups.

Conclusions

An effective detection rate was achieved in PJP patients using mNGS testing of bronchoalveolar lavage fluid (BALF) or blood. The study also confirmed that the abundance of reads of P. jirovecii is related to the interval between the onset and sample collection. And the inflammation status during simultaneous mNGS detection might determine the abundance of pathogens. Hence, we conclude that the mNGS strategy could benefit disease diagnosis as well as treatment when complicated clinical infections appeared.

A microfluidic cell chip for virus isolation via rapid screening for permissive cells

Virus identification is a prerequisite not only for the early diagnosis of viral infectious diseases but also for the effective prevention of epidemics. Successful cultivation is the gold standard for identifying a virus, according to the Koch postulates. However, this requires screening for a permissive cell line, which is traditionally time-, reagent- and labor-intensive. Here, a simple and easy-to-operate microfluidic chip, formed by seeding a variety of cell lines and culturing them in parallel, is reported for use in virus cultivation and virus-permissive host-cell screening. The chip was tested by infection with two known viruses, enterovirus 71 (EV71) and influenza virus H1N1. Infection with EV71 and H1N1 caused significant cytopathic effects (CPE) in RD and MDCK cells, respectively, demonstrating that virus cultivation based on this microfluidic cell chip can be used as a substitute for the traditional plate-based culture method and reproduce the typical CPE caused by virus infection. Using this microfluidic cell chip method for virus cultivation could make it possible to identify an emerging virus in a high-throughput, automatic, and unprecedentedly fast way.

•

A simple microfluidic chip for tandem culture of different cell lines is achieved.

•

The cell chip has been used for permissive cell screening and culture of viruses.

•

The cell chip has advantages of being sample-, reagent-, and time-saving.

•

The cell chip system holds potential for high-throughput and automated screening.

Intra-host variation and evolutionary dynamics of adenoviruses correlate to neutrophils in infected patients

With deep sequencing of virus genomes within the hosts, intra-host single nucleotide variations (iSNVs) have been used for analyses of virus genome variation and evolution, which is indicated to correlate with viral pathogenesis and disease severity. Little is known about the features of iSNVs among DNA viruses. We performed the epidemiological and laboratory investigation of one outbreak of adenovirus. The whole genomes of viruses in both original oral swabs and cell-cultured virus isolates were deeply sequenced. We identified 737 iSNVs in the viral genomes sequenced from original samples and 46 viral iSNVs in cell cultured isolates, with 33 iSNVs shared by original samples and cultured isolates. Meanwhile, we found these 33 iSNVs were shared by different patients, among which, three hot-spot areas 6367-6401, 9213-9247 and 10584-10606 within the functional genes of the adenovirus genome were found. Notably, the substitution rates of iSNVs were closely correlated with the clinical and immune indicators of the patients. Especially a positive correlation to neutrophils was found, indicating a predictable biomarker of iSNV dynamics. Our findings demonstrated the neutrophil-correlated dynamic evolution features of the iSNVs within adenoviruses, which indicates a virus-host interaction during human infection of a DNA virus. This article is protected by copyright. All rights reserved.

Investigation of cross-regional spread and evolution of equine influenza H3N8 at US and global scales using Bayesian phylogeography based on balanced subsampling

Equine influenza virus (EIV) is a highly contagious pathogen of equids, and a well-known burden in global equine health. EIV H3N8 variants seasonally emerged and resulted in EIV outbreaks in the United States (US) and worldwide. The present study evaluated the pattern of cross-regional EIV H3N8 spread and evolutionary characteristics at US and global scales using Bayesian phylogeography with balanced subsampling based on regional horse population size. A total of 297 Haemagglutinin (HA) sequences of global EIV H3N8 were collected from 1963 to 2019 and subsampled to global subset (n = 67), raw US sequences (n = 100) and US subset (n = 44) datasets. Discrete trait phylogeography analysis was used to estimate the transmission history of EIV using four global and US genome datasets. The North American lineage was the major source of globally dominant EIV variants and spread to other global regions. The US EIV strains generally spread from the southern and midwestern regions to other regions. The EIV H3N8 accumulated approximately three nucleotide substitutions per year in the HA gene under heterogenous local positive selection. Our findings will guide better decision making of target intervention strategies of EIV H3N8 infection and provide the better scheme of genomic surveillance in the US and global equine health. This article is protected by copyright. All rights reserved.

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.

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.

ADAR Editing in Viruses: An Evolutionary Force to Reckon with

Adenosine Deaminases that Act on RNA (ADARs) are RNA editing enzymes that play a dynamic and nuanced role in regulating transcriptome and proteome diversity. This editing can be highly selective, affecting a specific site within a transcript, or nonselective, resulting in hyperediting. ADAR editing is important for regulating neural functions and autoimmunity, and has a key role in the innate immune response to viral infections, where editing can have a range of pro- or antiviral effects and can contribute to viral evolution. Here we examine the role of ADAR editing across a broad range of viral groups. We propose that the effect of ADAR editing on viral replication, whether pro- or antiviral, is better viewed as an axis rather than a binary, and that the specific position of a given virus on this axis is highly dependent on virus- and host-specific factors, and can change over the course of infection. However, more research needs to be devoted to understanding these dynamic factors and how they affect virus-ADAR interactions and viral evolution. Another area that warrants significant attention is the effect of virus-ADAR interactions on host-ADAR interactions, particularly in light of the crucial role of ADAR in regulating neural functions. Answering these questions will be essential to developing our understanding of the relationship between ADAR editing and viral infection. In turn, this will further our understanding of the effects of viruses such as SARS-CoV-2, as well as many others, and thereby influence our approach to treating these deadly diseases.

On-Demand Patient-Specific Phenotype-to-Genotype Ebola Virus Characterization

Biosafety, biosecurity, logistical, political, and technical considerations can delay or prevent the wide dissemination of source material containing viable virus from the geographic origin of an outbreak to laboratories involved in developing medical countermeasures (MCMs). However, once virus genome sequence information is available from clinical samples, reverse-genetics systems can be used to generate virus stocks de novo to initiate MCM development. In this study, we developed a reverse-genetics system for natural isolates of Ebola virus (EBOV) variants Makona, Tumba, and Ituri, which have been challenging to obtain. These systems were generated starting solely with in silico genome sequence information and have been used successfully to produce recombinant stocks of each of the viruses for use in MCM testing. The antiviral activity of MCMs targeting viral entry varied depending on the recombinant virus isolate used. Collectively, selecting and synthetically engineering emerging EBOV variants and demonstrating their efficacy against available MCMs will be crucial for answering pressing public health and biosecurity concerns during Ebola disease (EBOD) outbreaks.

Assessment of two-pool multiplex long-amplicon nanopore sequencing of SARS-CoV-2

Genomic surveillance of severe acute respiratory syndrome coronavirus 2 (SARS‐CoV‐2) plays an important role in COVID‐19 pandemic control and elimination efforts, especially by elucidating its global transmission network and illustrating its viral evolution. The deployment of multiplex PCR assays that target SARS‐CoV‐2 followed by either massively parallel or nanopore sequencing is a widely‐used strategy to obtain genome sequences from primary samples. However, multiplex PCR‐based sequencing carries an inherent bias of sequencing depth among different amplicons, which may cause uneven coverage. Here we developed a two‐pool, long‐amplicon 36‐plex PCR primer panel with ~1000‐bp amplicon lengths for full‐genome sequencing of SARS‐CoV‐2. We validated the panel by assessing nasopharyngeal swab samples with a <30 quantitative reverse transcription PCR cycle threshold value and found that ≥90% of viral genomes could be covered with high sequencing depths (≥20% mean depth). In comparison, the widely‐used ARTIC panel yielded 79%–88% high‐depth genome regions. We estimated that ~5 Mbp nanopore sequencing data may ensure a >95% viral genome coverage with a ≥10‐fold depth and may generate reliable genomes at consensus sequence levels. Nanopore sequencing yielded false‐positive variations with frequencies of supporting reads <0.8, and the sequencing errors mostly occurred on the 5′ or 3′ ends of reads. Thus, nanopore sequencing could not elucidate intra‐host viral diversity.

Molecular evolutionary characteristics of SARS-CoV-2 emerging in the United States

SARS‐CoV‐2 is a newly discovered beta coronavirus at the end of 2019, which is highly pathogenic and poses a serious threat to human health. In this paper, 1875 SARS‐CoV‐2 whole genome sequences and the sequence coding spike protein (S gene) sampled from the United States were used for bioinformatics analysis to study the molecular evolutionary characteristics of its genome and spike protein. The MCMC method was used to calculate the evolution rate of the whole genome sequence and the nucleotide mutation rate of the S gene. The results showed that the nucleotide mutation rate of the whole genome was 6.677 × 10⁻⁴ substitution per site per year, and the nucleotide mutation rate of the S gene was 8.066 × 10⁻⁴ substitution per site per year, which was at a medium level compared with other RNA viruses. Our findings confirmed the scientific hypothesis that the rate of evolution of the virus gradually decreases over time. We also found 13 statistically significant positive selection sites in the SARS‐CoV‐2 genome. In addition, the results showed that there were 101 nonsynonymous mutation sites in the amino acid sequence of S protein, including seven putative harmful mutation sites. This paper has preliminarily clarified the evolutionary characteristics of SARS‐CoV‐2 in the United States, providing a scientific basis for future surveillance and prevention of virus variants.

3' untranslated regions of Marburg and Ebola virus mRNAs possess negative regulators of translation that are modulated by ADAR1 editing

The filovirus family includes deadly pathogens such as Ebola virus (EBOV) and Marburg virus (MARV). A substantial portion of filovirus genomes encode 5′ and 3′ untranslated regions (UTRs) of viral mRNAs. Select viral genomic RNA sequences corresponding to 3′ UTRs are prone to editing by adenosine deaminase acting on RNA 1 (ADAR1). A reporter mRNA approach, in which different 5′ or 3′ UTRs were inserted into luciferase-encoding mRNAs, demonstrates that MARV 3′ UTRs yield different levels of reporter gene expression, suggesting modulation of translation. The modulation occurs in cells unable to produce microRNAs (miRNAs) and can be recapitulated in a MARV minigenome assay. Deletion mutants identified negative regulatory regions at the ends of the MARV nucleoprotein (NP) and large protein (L) 3′ UTRs. Apparent ADAR1 editing mutants were previously identified within the MARV NP 3′ UTR. Introduction of these changes into the MARV nucleoprotein (NP) 3′ UTR or deletion of the region targeted for editing enhances translation, as indicated by reporter assays and polysome analysis. In addition, the parental NP 3′ UTR, but not the edited or deletion mutant NP 3′ UTRs, induces a type I interferon (IFN) response upon transfection into cells. Because some EBOV isolates from the West Africa outbreak exhibited ADAR1 editing of the viral protein of 40 kDa (VP40) 3′ UTR, VP40 3′ UTRs with parental and edited sequences were similarly assayed. The EBOV VP40 3′ UTR edits also enhanced translation, but neither the wild-type nor the edited 3′ UTRs induced IFN. These findings implicate filoviral mRNA 3′ UTRs as negative regulators of translation that can be inactivated by innate immune responses that induce ADAR1.

IMPORTANCE UTRs comprise a large percentage of filovirus genomes and are apparent targets of editing by ADAR1, an enzyme with pro- and antiviral activities. However, the functional significance of the UTRs and ADAR1 editing has been uncertain. This study demonstrates that MARV and EBOV 3′ UTRs can modulate translation, in some cases negatively. ADAR1 editing or deletion of select regions within the translation suppressing 3′ UTRs relieves the negative effects of the UTRs. These data indicate that filovirus 3′ UTRs contain translation regulatory elements that are modulated by activation of ADAR1, suggesting a complex interplay between filovirus gene expression and innate immunity.

Simultaneous Detection of Ebola Virus and Pathogens Associated With Hemorrhagic Fever by an Oligonucleotide Microarray

Ebola virus infection causes severe hemorrhagic fever, and its mortality rates varied from 25 to 90% in the previous outbreaks. The highly infectious and lethal nature of this virus highlights the need for reliable and sensitive diagnostic methods to distinguish it from other diseases present with similar clinical symptoms. Based on multiplex polymerase chain reaction (PCR) and oligonucleotide microarray technology, a cost-effective, multipathogen and high-throughput method was developed for simultaneous detection of Ebola virus and other pathogens associated with hemorrhagic fever, including Marburg virus, Lassa fever virus, Junin virus, Machupo virus, Rift Valley fever virus, Crimean-Congo hemorrhagic fever virus, malaria parasite, hantavirus, severe fever with thrombocytopenia syndrome virus, dengue virus, yellow fever virus, Chikungunya virus, influenza A virus, and influenza B virus. This assay had an excellent specificity for target pathogens, without overlap signal between the probes. The limit of detection was approximately 10³ pathogen copies/μl. A total of 60 positive nucleic acid samples for different pathogens were detected, a concordance of 100% was observed between microarray assay and real-time PCR analysis. Consequently, the described oligonucleotide microarray may be specific and sensitive assay for diagnosis and surveillance of infections caused by Ebola virus and other species of hemorrhagic fever pathogens.

[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.

Species-Specific Evolution of Ebola Virus during Replication in Human and Bat Cells

Ebola virus (EBOV) causes a severe, often fatal disease in humans and nonhuman primates. Within the past decade, EBOV has caused two large and difficult-to-control outbreaks, one of which recently ended in the Democratic Republic of the Congo. Bats are the likely reservoir of EBOV, but little is known of their relationship with the virus. We perform serial passages of EBOV in human and bat cells and use circular sequencing to compare the short-term evolution of the virus. Virus populations passaged in bat cells have sequence markers indicative of host RNA editing enzyme activity, including evidence for ADAR editing of the EBOV glycoprotein. Multiple regions in the EBOV genome appear to have undergone adaptive evolution when passaged in bat and human cells. Individual mutated viruses are rescued and characterized. Our results provide insight into the host species-specific evolution of EBOV and highlight the adaptive flexibility of the virus.

Forty-two years of responding to Ebola virus outbreaks in Sub-Saharan Africa: a review

Introduction

Ebola virus disease (EVD) is one of the deadliest haemorrhagic fevers affecting humans and non-human primates. Thirty-four outbreaks have been reported in Africa since it was first recognised in 1976. This review analysed 42 years of EVD outbreaks and identified various challenges and opportunities for its control and prevention in Sub-Saharan Africa.

Methods

A literature search of relevant articles on EVD was done in PubMed, Web of Science and Google Scholar electronic databases. Articles published from 1976 to 2019 were reviewed to document reports of EVD outbreaks in Sub-Saharan Africa. Data extraction focused on the year of outbreak, geographical spread, virus strain involved, number of cases and deaths, case fatality, and outbreak management. Analyses of trends in case fatality were performed by calculating ORs between times.

Results

In the past four decades, a total of 34 EVD outbreaks affecting 34 356 cases and causing 14 823 deaths were reported in 11 countries in Sub-Saharan Africa. The overall case fatality rate (95% CI) was 66% (62 to 71) and did not change substantially over time (OR in 2019 vs 1976=1.6 (95% CI 1.5 to 1.8), p<0.001). The results of this review indicate that challenges to control EVD outbreaks are related to epidemiological, sociocultural and health system factors.

Conclusions

Sub-Saharan Africa continues to face considerable challenges in EVD control, whereby there are no significant changes in case fatality rates observed during the past four decades. Socioeconomic and cultural processes need to be critically considered to shape the community behaviours that lead to exposure to EVD outbreaks. Areas that need to be addressed to prevent future EVD outbreaks include a broad-based, one-health approach, effective communication, social mobilisation, and strengthening of the health systems.

Ebola virus disease

Ebola virus disease (EVD) is a severe and frequently lethal disease caused by Ebola virus (EBOV). EVD outbreaks typically start from a single case of probable zoonotic transmission, followed by human-to-human transmission via direct contact or contact with infected bodily fluids or contaminated fomites. EVD has a high case-fatality rate; it is characterized by fever, gastrointestinal signs and multiple organ dysfunction syndrome. Diagnosis requires a combination of case definition and laboratory tests, typically real-time reverse transcription PCR to detect viral RNA or rapid diagnostic tests based on immunoassays to detect EBOV antigens. Recent advances in medical countermeasure research resulted in the recent approval of an EBOV-targeted vaccine by European and US regulatory agencies. The results of a randomized clinical trial of investigational therapeutics for EVD demonstrated survival benefits from two monoclonal antibody products targeting the EBOV membrane glycoprotein. New observations emerging from the unprecedented 2013-2016 Western African EVD outbreak (the largest in history) and the ongoing EVD outbreak in the Democratic Republic of the Congo have substantially improved the understanding of EVD and viral persistence in survivors of EVD, resulting in new strategies toward prevention of infection and optimization of clinical management, acute illness outcomes and attendance to the clinical care needs of patients.

Pre-emptive genomic surveillance of emerging ebolaviruses

Genomic surveillance during ebolavirus outbreaks to elucidate transmission chains and develop diagnostic tests is delayed by the laborious development of variant-specific laboratory assays. We developed a new protocol combining 31 parallel PCR assays with Illumina/MinION-based sequencing, allowing generic ebolavirus genomic surveillance, validated using cell culture-derived Ebola, Reston, Sudan and Taï Forest virus at concentrations compatible with patient viral loads. Our approach enables pre-emptive genomic surveillance of ongoing and future ebolavirus outbreaks irrespective of variant divergence.

Phylogeographic reconstruction using air transportation data and its application to the 2009 H1N1 influenza A pandemic

Influenza A viruses cause seasonal epidemics and occasional pandemics in the human population. While the worldwide circulation of seasonal influenza is at least partly understood, the exact migration patterns between countries, states or cities are not well studied. Here, we use the Sankoff algorithm for parsimonious phylogeographic reconstruction together with effective distances based on a worldwide air transportation network. By first simulating geographic spread and then phylogenetic trees and genetic sequences, we confirmed that reconstructions with effective distances inferred phylogeographic spread more accurately than reconstructions with geographic distances and Bayesian reconstructions with BEAST that do not use any distance information, and led to comparable results to the Bayesian reconstruction using distance information via a generalized linear model. Our method extends Bayesian methods that estimate rates from the data by using fine-grained locations like airports and inferring intermediate locations not observed among sampled isolates. When applied to sequence data of the pandemic H1N1 influenza A virus in 2009, our approach correctly inferred the origin and proposed airports mainly involved in the spread of the virus. In case of a novel outbreak, this approach allows to rapidly analyze sequence data and infer origin and spread routes to improve disease surveillance and control.

Influenza A viruses infect up to 5 million people in recurring epidemics every year. Further, viruses of zoonotic origin constantly pose a pandemic risk. Understanding the geographical spread of these viruses, including the origin and the main spread routes between cities, states or countries, could help to monitor or contain novel outbreaks. Based on genetic sequences and sampling locations, the geographic spread can be reconstructed along a phylogenetic tree. Our approach uses a parsimonious reconstruction with air transportation data and was verified using a simulation of the 2009 H1N1 influenza A pandemic. Applied to real sequence data of the outbreak, our analysis gave detailed insights into spread patterns of influenza A viruses, highlighting the origin as well as airports mainly involved in the spread.

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.

An evolutionary insight into emerging Ebolavirus strains isolated in Africa

On July 19, 2019, the World Health Organization declared the current Ebolavirus (EBOV) outbreak in Congo Democratic Republic (COD) a public health emergency of international concern. To address the potential threat of EBOV evolution outpacing antibody treatment and vaccine efforts, a detailed evolutionary analysis of EBOV strains circulating in different African countries was performed. Genome composition of EBOV strains was studied using multivariate statistical analysis. To investigate the patterns of evolution of EBOV strains, a Bayesian Markov Chain Monte Carlo approach was used. Two different genetic lineages, with a distinct genome composition gave rise to the recent EBOV outbreaks in central and western Africa. Strains isolated in COD in 2018 fall into two different genetic clusters, according to their geographical location of isolation. Different amino acid substitutions among strains from these two clusters have been found, particularly in NP, GP, and L proteins. Significant differences in codon and amino acid usage among clusters were found. Strains isolated in COD in 2018 belong to two distinct genetic clusters, with distinct codon and amino acid usage. Geographical diversity plays an important role in shaping the molecular evolution of EBOV populations.

A new twenty-first century science for effective epidemic response

With rapidly changing ecology, urbanization, climate change, increased travel and fragile public health systems, epidemics will become more frequent, more complex and harder to prevent and contain. Here we argue that our concept of epidemics must evolve from crisis response during discrete outbreaks to an integrated cycle of preparation, response and recovery. This is an opportunity to combine knowledge and skills from all over the world-especially at-risk and affected communities. Many disciplines need to be integrated, including not only epidemiology but also social sciences, research and development, diplomacy, logistics and crisis management. This requires a new approach to training tomorrow's leaders in epidemic prevention and response.

Chasing Intracellular Zika Virus Using Proteomics

Flaviviruses are the most medically relevant group of arboviruses causing a wide range of diseases in humans and are associated with high mortality and morbidity, as such posing a major health concern. Viruses belonging to this family can be endemic (e.g., dengue virus), but can also cause fulminant outbreaks (e.g., West Nile virus, Japanese encephalitis virus and Zika virus). Intense research efforts in the past decades uncovered shared fundamental strategies used by flaviviruses to successfully replicate in their respective hosts. However, the distinct features contributing to the specific host and tissue tropism as well as the pathological outcomes unique to each individual flavivirus are still largely elusive. The profound footprint of individual viruses on their respective hosts can be investigated using novel technologies in the field of proteomics that have rapidly developed over the last decade. An unprecedented sensitivity and throughput of mass spectrometers, combined with the development of new sample preparation and bioinformatics analysis methods, have made the systematic investigation of virus-host interactions possible. Furthermore, the ability to assess dynamic alterations in protein abundances, protein turnover rates and post-translational modifications occurring in infected cells now offer the unique possibility to unravel complex viral perturbations induced in the infected host. In this review, we discuss the most recent contributions of mass spectrometry-based proteomic approaches in flavivirus biology with a special focus on Zika virus, and their basic and translational potential and implications in understanding and characterizing host responses to arboviral infections.

Emergence and adaptive evolution of Nipah virus

Since its first emergence in 1998 in Malaysia, Nipah virus (NiV) has become a great threat to domestic animals and humans. Sporadic outbreaks associated with human-to-human transmission caused hundreds of human fatalities. Here, we collected all available NiV sequences and combined phylogenetics, molecular selection, structural biology and receptor analysis to study the emergence and adaptive evolution of NiV. NiV can be divided into two main lineages including the Bangladesh and Malaysia lineages. We formly confirmed a significant association with geography which is probably the result of long-term evolution of NiV in local bat population. The two NiV lineages differ in many amino acids; one change in the fusion protein might be involved in its activation via binding to the G protein. We also identified adaptive and positively selected sites in many viral proteins. In the receptor-binding G protein, we found that sites 384, 386 and especially 498 of G protein might modulate receptor-binding affinity and thus contribute to the host jump from bats to humans via the adaption to bind the human ephrin-B2 receptor. We also found that site 1645 in the connector domain of L was positive selected and involved in adaptive evolution; this site might add methyl groups to the cap structure present at the 5'-end of the RNA and thus modulate its activity. This study provides insight to assist the design of early detection methods for NiV to assess its epidemic potential in humans.

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

Recent advances in the development and evaluation of molecular diagnostics for Ebola virus disease

INTRODUCTION

The 2014-16 outbreak of ebola virus disease (EVD) in West Africa resulted in 11,308 deaths. During the outbreak only 60% of patients were laboratory confirmed and global health authorities have identified the need for accurate and readily deployable molecular diagnostics as an important component of the ideal response to future outbreaks, to quickly identify and isolate patients. Areas covered: Currently PCR-based techniques and rapid diagnostic tests (RDTs) that detect antigens specific to EVD infections dominate the diagnostic landscape, but recent advances in biosensor technologies have led to novel approaches for the development of EVD diagnostics. This review summarises the literature and available performance data of currently available molecular diagnostics for ebolavirus, identifies knowledge gaps and maps out future priorities for research in this field. Expert opinion: While there are now a plethora of diagnostic tests for EVD at various stages of development, there is an acute need for studies to compare their clinical performance, but the sporadic nature of EVD outbreaks makes this extremely challenging, demanding pragmatic new modalities of research funding and ethical/institutional approval, to enable responsive research in outbreak settings. Retrospective head-to-head diagnostic comparisons could also be implemented using biobanked specimens, providing this can be done safely.

Expanding the watch list for potential Ebola virus antibody escape mutations

The 2014 outbreak of Ebola virus disease (EVD) in Western Africa is the largest recorded filovirus disease outbreak and led to the death of over 11,000 people. The recent EVD outbreaks (since May 2018) in the Democratic Republic of the Congo has already claimed the lives of over 250 people. Tackling Ebola virus (EBOV) outbreaks remains a challenge. Over the years, significant efforts have been put into developing vaccines or antibody therapies which rely on an envelope glycoprotein (GP) of Zaire ebolavirus (strain Mayinga-76). Therefore, one key approach for combating EVD epidemics is to predict mutations that may diminish the effectiveness of the treatment. In a previous study we generated a watch list of potential antibody escape mutations of EBOV GP against the monoclonal antibody KZ52. Molecular modeling methods were applied to the three-dimensional experimental structure of EBOV GP bound to KZ52 to predict the effect of every possible single mutation in EBOV GP. The final watch list contained 34 mutations that were predicted to destabilize binding of KZ52 to EBOV GP but did not affect EBOV GP folding and its ability to form trimers. In this study, we expand our watch list by including three more monoclonal antibodies with distinct epitopes on GP, namely Antibody 100 (Ab100), Antibody 114 (Ab114) and 13F6-1-2. Our updated watch list contains 127 mutations, three of which have been seen in humans or are experimentally associated with reduced efficacy of antibody treatment. We believe mutations on this watch list require attention since they provide information about circumstances in which interventions could lose the effectiveness.

Bayesian inference of transmission chains using timing of symptoms, pathogen genomes and contact data

There exists significant interest in developing statistical and computational tools for inferring 'who infected whom' in an infectious disease outbreak from densely sampled case data, with most recent studies focusing on the analysis of whole genome sequence data. However, genomic data can be poorly informative of transmission events if mutations accumulate too slowly to resolve individual transmission pairs or if there exist multiple pathogens lineages within-host, and there has been little focus on incorporating other types of outbreak data. We present here a methodology that uses contact data for the inference of transmission trees in a statistically rigorous manner, alongside genomic data and temporal data. Contact data is frequently collected in outbreaks of pathogens spread by close contact, including Ebola virus (EBOV), severe acute respiratory syndrome coronavirus (SARS-CoV) and Mycobacterium tuberculosis (TB), and routinely used to reconstruct transmission chains. As an improvement over previous, ad-hoc approaches, we developed a probabilistic model that relates a set of contact data to an underlying transmission tree and integrated this in the outbreaker2 inference framework. By analyzing simulated outbreaks under various contact tracing scenarios, we demonstrate that contact data significantly improves our ability to reconstruct transmission trees, even under realistic limitations on the coverage of the contact tracing effort and the amount of non-infectious mixing between cases. Indeed, contact data is equally or more informative than fully sampled whole genome sequence data in certain scenarios. We then use our method to analyze the early stages of the 2003 SARS outbreak in Singapore and describe the range of transmission scenarios consistent with contact data and genetic sequence in a probabilistic manner for the first time. This simple yet flexible model can easily be incorporated into existing tools for outbreak reconstruction and should permit a better integration of genomic and epidemiological data for inferring transmission chains.

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.

From "A"IV to "Z"IKV: Attacks from Emerging and Re-emerging Pathogens

100 years after the infamous “Spanish flu” pandemic, the 2017–2018 flu season has been severe, with numerous infections worldwide. In between, there have been continuous, relentless attacks from (re-)emerging viruses. To fully understand viral pathogenesis and develop effective medical countermeasures, we must strengthen current surveillance and basic research efforts.

Phylodynamic Analysis of Ebola Virus Disease Transmission in Sierra Leone

We generated genome sequences from 218 cases of Ebola virus disease (EVD) in Sierra Leone (SLE) during 2014–2015 to complement available datasets, particularly by including cases from a period of low sequence coverage during peak transmission of Ebola virus (EBOV) in the highly-affected Western Area division of SLE. The combined dataset was utilized to produce phylogenetic and phylodynamic inferences, to study sink–source dynamics and virus dispersal from highly-populated transmission hotspots. We identified four districts in SLE where EBOV was introduced and transmission occurred without onward exportation to other districts. We also identified six districts that substantially contributed to the dispersal of the virus and prolonged the EVD outbreak: five of these served as major hubs, with lots of movement in and out, and one acted primarily as a source, exporting the virus to other areas of the country. Positive correlations between case numbers, inter-district transition events, and district population sizes reaffirm that population size was a driver of EBOV transmission dynamics in SLE. The data presented here confirm the role of urban hubs in virus dispersal and of a delayed laboratory response in the expansion and perpetuation of the EVD outbreak in SLE.

Advances in Molecular Diagnostic Approaches for Biothreat Agents

The advancement in Molecular techniques has been implicated in the development of sophisticated, high-end diagnostic platform and point-of-care (POC) devices for the detection of biothreat agents. Different molecular and immunological approaches such as Immunochromatographic and lateral flow assays, Enzyme-linked Immunosorbent assays (ELISA), Biosensors, Isothermal amplification assays, Nucleic acid amplification tests (NAATs), Next Generation Sequencers (NGS), Microarrays and Microfluidics have been used for a long time as detection strategies of the biothreat agents. In addition, several point of care (POC) devices have been approved by FDA and commercialized in markets. The high-end molecular platforms like NGS and Microarray are time-consuming, costly, and produce huge amount of data. Therefore, the future prospects of molecular based technique should focus on developing quick, user-friendly, cost-effective and portable devices against biological attacks and surveillance programs.

Exact limits of inference in coalescent models

Recovery of population size history from molecular sequence data is an important problem in population genetics. Inference commonly relies on a coalescent model linking the population size history to genealogies. The high computational cost of estimating parameters from these models usually compels researchers to select a subset of the available data or to rely on insufficient summary statistics for statistical inference. We consider the problem of recovering the true population size history from two possible alternatives on the basis of coalescent time data previously considered by Kim et al. (2015). We improve upon previous results by giving exact expressions for the probability of correctly distinguishing between the two hypotheses as a function of the separation between the alternative size histories, the number of individuals, loci, and the sampling times. In more complicated settings we estimate the exact probability of correct recovery by Monte Carlo simulation. Our results give considerably more pessimistic inferential limits than those previously reported. We also extended our analyses to pairwise SMC and SMC’ models of recombination. This work is relevant for optimal design when the inference goal is to test scientific hypotheses about population size trajectories in coalescent models with and without recombination.

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.

outbreaker2: a modular platform for outbreak reconstruction

BACKGROUND

Reconstructing individual transmission events in an infectious disease outbreak can provide valuable information and help inform infection control policy. Recent years have seen considerable progress in the development of methodologies for reconstructing transmission chains using both epidemiological and genetic data. However, only a few of these methods have been implemented in software packages, and with little consideration for customisability and interoperability. Users are therefore limited to a small number of alternatives, incompatible tools with fixed functionality, or forced to develop their own algorithms at considerable personal effort.

RESULTS

Here we present outbreaker2, a flexible framework for outbreak reconstruction. This R package re-implements and extends the original model introduced with outbreaker, but most importantly also provides a modular platform allowing users to specify custom models within an optimised inferential framework. As a proof of concept, we implement the within-host evolutionary model introduced with TransPhylo, which is very distinct from the original genetic model in outbreaker, and demonstrate how even complex model results can be successfully included with minimal effort.

CONCLUSIONS

outbreaker2 provides a valuable starting point for future outbreak reconstruction tools, and represents a unifying platform that promotes customisability and interoperability. Implemented in the R software, outbreaker2 joins a growing body of tools for outbreak analysis.

An unusually high substitution rate in transplant-associated BK polyomavirus in vivo is further concentrated in HLA-C-bound viral peptides

Infection with human BK polyomavirus, a small double-stranded DNA virus, potentially results in severe complications in immunocompromised patients. Here, we describe the in vivo variability and evolution of the BK polyomavirus by deep sequencing. Our data reveal the highest genomic evolutionary rate described in double-stranded DNA viruses, i.e., 10⁻³–10⁻⁵ substitutions per nucleotide site per year. High mutation rates in viruses allow their escape from immune surveillance and adaptation to new hosts. By combining mutational landscapes across viral genomes with in silico prediction of viral peptides, we demonstrate the presence of significantly more coding substitutions within predicted cognate HLA-C-bound viral peptides than outside. This finding suggests a role for HLA-C in antiviral immunity, perhaps through the action of killer cell immunoglobulin-like receptors. The present study provides a comprehensive view of viral evolution and immune escape in a DNA virus.

Little is known about the mechanisms of evolution and viral immune escape in double-stranded DNA (dsDNA) viruses. Here, we study the evolution of BK polyomavirus and observe the highest genomic evolutionary rate described so far for a dsDNA virus, in the range of RNA viruses, which usually evolve rapidly. Furthermore, the prediction of viral peptides to determine immune escape suggests a specific role of HLA-C in antiviral immunity. These findings are helpful for future advances in antiviral therapies and provide a step forward in our understanding of in vivo viral evolution in humans.

Growth-Adaptive Mutations in the Ebola Virus Makona Glycoprotein Alter Different Steps in the Virus Entry Pathway

The Zaire ebolavirus (EBOV) glycoprotein (GP) is cleaved into two subunits (GP1 and GP2) that are both required for virus attachment and entry into cells. Sequence changes in the GP have been proposed to increase pathogenesis and to alter virus growth properties. Mutations in GP acquired during EBOV tissue culture passage have also been reported to change virus growth properties. Here, we report the isolation of six amino acid mutations in EBOV GP that spontaneously appeared during recovery and passage of an EBOV-Makona GP-pseudotyped vesicular stomatitis virus (VSV), two of which also occur during passage of EBOV clinical isolates in tissue culture. Each of the six mutations resulted in increased virus growth in monkey and human cell lines. All mutations are located in the GP2 fusion subunit and increase entry kinetics of EBOV virus-like particles (VLPs). The gain-of-entry function mapped to two mechanistic phenotypes. Mutations in heptad repeat 1 (HR1) decreased the requirement for cathepsin B activity for viral infection. Mutations directly within the fusion loop increased entry kinetics without altering the cathepsin B dependence. Several mutations in the fusion loop were substitutions of residues present in other ebolavirus glycoproteins, illustrating the evolutionary paths for maintaining an optimally functioning fusion loop under selection pressure.IMPORTANCEZaire ebolavirus (EBOV) is the causative agent of the highly lethal Ebola virus disease and poses a significant threat to the global health community. Approved antivirals against EBOV are lacking; however, promising therapies targeting the EBOV glycoprotein are being developed. Efficacy testing of these candidate therapeutics relies on EBOV laboratory stocks, which when grown in tissue culture may acquire mutations in the glycoprotein. These mutations can produce inaccurate results in therapeutic testing. Until recently, distinguishing between tissue culture mutations and naturally occurring polymorphisms in EBOV GP was difficult in the absence of consensus clinical GP sequences. Here, we utilize recombinant VSV (rVSV) pseudotyped with the consensus clinical EBOV Makona GP to identify several mutations that have emerged or have potential to emerge in EBOV GP during tissue culture passage. Identifying these mutations informs the EBOV research community as to which mutations may arise during preparation of laboratory virus stocks.

Identifying spatio-temporal dynamics of Ebola in Sierra Leone using virus genomes

Containing the recent West African outbreak of Ebola virus (EBOV) required the deployment of substantial global resources. Despite recent progress in analysing and modelling EBOV epidemiological data, a complete characterization of the spatio-temporal spread of Ebola cases remains a challenge. In this work, we offer a novel perspective on the EBOV epidemic in Sierra Leone that uses individual virus genome sequences to inform population-level, spatial models. Calibrated to phylogenetic linkages of virus genomes, these spatial models provide unique insight into the disease mobility of EBOV in Sierra Leone without the need for human mobility data. Consistent with other investigations, our results show that the spread of EBOV during the beginning and middle portions of the epidemic strongly depended on the size of and distance between populations. Our phylodynamic analysis also revealed a change in model preference towards a spatial model with power-law characteristics in the latter portion of the epidemic, correlated with the timing of major intervention campaigns. More generally, we believe this framework, pairing molecular diagnostics with a dynamic model selection procedure, has the potential to be a powerful forecasting tool along with offering operationally relevant guidance for surveillance and sampling strategies during an epidemic.

Mechanisms and Concepts in RNA Virus Population Dynamics and Evolution

RNA viruses are unique in their evolutionary capacity, exhibiting high mutation rates and frequent recombination. They rapidly adapt to environmental changes, such as shifts in immune pressure or pharmacological challenge. The evolution of RNA viruses has been brought into new focus with the recent developments of genetic and experimental tools to explore and manipulate the evolutionary dynamics of viral populations. These studies have uncovered new mechanisms that enable viruses to overcome evolutionary challenges in the environment and have emphasized the intimate relationship of viral populations with evolution. Here, we review some of the emerging viral and host mechanisms that underlie the evolution of RNA viruses. We also discuss new studies that demonstrate that the relationship between evolutionary dynamics and virus biology spans many spatial and temporal scales, affecting transmission dynamics within and between hosts as well as pathogenesis.

Portable sequencing, genomic data, and scale in global emerging infectious disease surveillance

Emerging infectious diseases (EIDs) occur when pathogens unpredictably spread into new contexts. EID surveillance systems seek to rapidly identify EID outbreaks to contain spread and improve public health outcomes. Sequencing data has historically not been integrated into real-time responses, but portable DNA sequencing technology has prompted optimism among epidemiologists. Specifically, attention has focused on the goal of a "sequencing singularity": the integration of portable sequencers in a worldwide event-based surveillance network with other digital data (Gardy & Loman, Nature Reviews Genetics, 19, 2018, p. 9). The sequencing singularity vision is a powerful socio-technical imaginary, shaping the discourse around the future of portable sequencing. Ethical and practical issues are bound by the vision in two ways: they are framed only as obstacles, and they are formulated only at the scales made visible by its implicit geography. This geography privileges two extremes of scale - the genomic and the global - and leaves intermediate scales comparatively unmapped. We explore how widespread portable sequencing could challenge this geography. Portable sequencers put the ability to produce genomic data in the hands of the individual. The explicit assertion of rights over data may therefore become a matter disputed more at an interpersonal scale than an international one. Portable sequencers also promise ubiquitous, indiscriminate sequencing of the total metagenomic content of samples, raising the question of what (or who) is under surveillance and inviting consideration of the human microbiome and more-than-human geographies. We call into question a conception of a globally integrated stream of sequencing data as composed mostly of "noise," within which signals of pathogen "emergence" are "hidden," considering it instead from the perspective of recent work into more-than-human geographies. Our work highlights a practical need for researchers to consider both the alternative possibilities they foreclose as well as the exciting opportunities they move towards when they deploy their visions of the future.

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.

A Rapid Screen for Host-Encoded miRNAs with Inhibitory Effects against Ebola Virus Using a Transcription- and Replication-Competent Virus-Like Particle System

MicroRNAs (miRNAs) may become efficient antiviral agents against the Ebola virus (EBOV) targeting viral genomic RNAs or transcripts. We previously conducted a genome-wide search for differentially expressed miRNAs during viral replication and transcription. In this study, we established a rapid screen for miRNAs with inhibitory effects against EBOV using a tetracistronic transcription- and replication-competent virus-like particle (trVLP) system. This system uses a minigenome comprising an EBOV leader region, luciferase reporter, VP40, GP, VP24, EBOV trailer region, and three noncoding regions from the EBOV genome and can be used to model the life cycle of EBOV under biosafety level (BSL) 2 conditions. Informatic analysis was performed to select up-regulated miRNAs targeting the coding regions of the minigenome with the highest binding energy to perform inhibitory effect screening. Among these miRNAs, miR-150-3p had the most significant inhibitory effect. Reverse transcription polymerase chain reaction (RT-PCR), Western blot, and double fluorescence reporter experiments demonstrated that miR-150-3p inhibited the reproduction of trVLPs via the regulation of GP and VP40 expression by directly targeting the coding regions of GP and VP40. This novel, rapid, and convenient screening method will efficiently facilitate the exploration of miRNAs against EBOV under BSL-2 conditions.

Novel cyclo-peptides inhibit Ebola pseudotyped virus entry by targeting primed GP protein

Ebola virus (EBOV) causes fatal hemorrhagic fever with high death rates in human. Currently, there are no available clinically-approved prophylactic or therapeutic treatments. The recently solved crystal structure of cleavage-primed EBOV glycoprotein (GPcl) in complex with the C domain of endosomal protein Niemann-Pick C1 (NPC1) provides a new target for the development of EBOV entry inhibitors. In this work, a computational approach using docking and molecular dynamic simulations is carried out for the rational design of peptide inhibitors. A novel cyclo-peptide (Pep-3.3) was identified to target at the late stage of EBOV entry and exhibit specific inhibitory activity against EBOV-GP pseudotyped viruses, with 50% inhibitory concentration (IC50) of 5.1 μM. In vitro binding assay and molecular simulations revealed that Pep-3.3 binds to GPcl with a KD value of 69.7 μM, through interacting with predicted residues in the hydrophobic binding pocket of GPcl. Mutation of predicted residues T83 caused resistance to Pep-3.3 inhibition in viral infectivity, providing preliminary support for the model of the peptide binding to GPcl. This study demonstrates the feasibility of inhibiting EBOV entry by targeting GPcl with peptides.