Genomic epidemiology reveals multiple introductions of Zika virus into the United States.

Zika virus (ZIKV) is causing an unprecedented epidemic linked to severe congenital abnormalities. In July 2016, mosquito-borne ZIKV transmission was reported in the continental United States; since then, hundreds of locally acquired infections have been reported in Florida. To gain insights into the timing, source, and likely route(s) of ZIKV introduction, we tracked the virus from its first detection in Florida by sequencing ZIKV genomes from infected patients and Aedes aegypti mosquitoes. We show that at least 4 introductions, but potentially as many as 40, contributed to the outbreak in Florida and that local transmission is likely to have started in the spring of 2016-several months before its initial detection. By analysing surveillance and genetic data, we show that ZIKV moved among transmission zones in Miami. Our analyses show that most introductions were linked to the Caribbean, a finding corroborated by the high incidence rates and traffic volumes from the region into the Miami area. Our study provides an understanding of how ZIKV initiates transmission in new regions.

Person-to-Person Transmission of Andes Virus in Hantavirus Pulmonary Syndrome, Argentina, 2014.

Andes virus is unique among hantaviruses because it can be transmitted from person to person. This mechanism was previously supported by epidemiologic data and genetic evidence based only on partial sequences. We used full-length virus sequencing to confirm person-to-person transmission of this virus in a cluster of 3 cases in Argentina in 2014.

Complete genomic sequences of Venezuelan equine encephalitis virus subtype IIID isolates from mosquitoes.

Venezuelan equine encephalitis virus (VEEV) is an important pathogen of medical and veterinary importance in the Americas. In this report, we present the complete genome sequences of five VEEV isolates obtained from pools of Culex (Melanoconion) gnomatos (4) or Culex (Melanoconion) pedroi (1) from Iquitos, Peru. Genetic and phylogenetic analyses showed that all five isolates grouped within the VEEV complex sister to VEEV IIIC and are members of subtype IIID. This is the first report of full-length genomic sequences of VEEV IIID.

Molecular Evidence of Sexual Transmission of Ebola Virus.

A suspected case of sexual transmission from a male survivor of Ebola virus disease (EVD) to his female partner (the patient in this report) occurred in Liberia in March 2015. Ebola virus (EBOV) genomes assembled from blood samples from the patient and a semen sample from the survivor were consistent with direct transmission. The genomes shared three substitutions that were absent from all other Western African EBOV sequences and that were distinct from the last documented transmission chain in Liberia before this case. Combined with epidemiologic data, the genomic analysis provides evidence of sexual transmission of EBOV and evidence of the persistence of infective EBOV in semen for 179 days or more after the onset of EVD. (Funded by the Defense Threat Reduction Agency and others.).

Monitoring of Ebola Virus Makona Evolution through Establishment of Advanced Genomic Capability in Liberia.

To support Liberia's response to the ongoing Ebola virus (EBOV) disease epidemic in Western Africa, we established in-country advanced genomic capabilities to monitor EBOV evolution. Twenty-five EBOV genomes were sequenced at the Liberian Institute for Biomedical Research, which provided an in-depth view of EBOV diversity in Liberia during September 2014-February 2015. These sequences were consistent with a single virus introduction to Liberia; however, shared ancestry with isolates from Mali indicated at least 1 additional instance of movement into or out of Liberia. The pace of change is generally consistent with previous estimates of mutation rate. We observed 23 nonsynonymous mutations and 1 nonsense mutation. Six of these changes are within known binding sites for sequence-based EBOV medical countermeasures; however, the diagnostic and therapeutic impact of EBOV evolution within Liberia appears to be low.

SUMOylation affects the interferon blocking activity of the influenza A nonstructural protein NS1 without affecting its stability or cellular localization.

Our pioneering studies on the interplay between the small ubiquitin-like modifier (SUMO) and influenza A virus identified the nonstructural protein NS1 as the first known SUMO target of influenza virus and one of the most abundantly SUMOylated influenza virus proteins. Here, we further characterize the role of SUMOylation for the A/Puerto Rico/8/1934 (PR8) NS1 protein, demonstrating that NS1 is SUMOylated not only by SUMO1 but also by SUMO2/3 and mapping the main SUMOylation sites in NS1 to residues K219 and K70. Furthermore, by using SUMOylatable and non-SUMOylatable forms of NS1 and an NS1-specific artificial SUMO ligase (ASL) that increases NS1 SUMOylation ~4-fold, we demonstrate that SUMOylation does not affect the stability or cellular localization of PR8 NS1. However, NS1's ability to be SUMOylated appears to affect virus multiplication, as indicated by the delayed growth of a virus expressing the non-SUMOylatable form of NS1 in the interferon (IFN)-competent MDCK cell line. Remarkably, while a non-SUMOylatable form of NS1 exhibited a substantially diminished ability to neutralize IFN production, increasing NS1 SUMOylation beyond its normal levels also exerted a negative effect on its IFN-blocking function. This observation indicates the existence of an optimal level of NS1 SUMOylation that allows NS1 to achieve maximal activity and suggests that the limited amount of SUMOylation normally observed for most SUMO targets may correspond to an optimal level that maximizes the contribution of SUMOylation to protein function. Finally, protein cross-linking data suggest that SUMOylation may affect NS1 function by regulating the abundance of NS1 dimers and trimers in the cell.

Andes Virus Genome Mutations That Are Likely Associated with Animal Model Attenuation and Human Person-to-Person Transmission.

We performed whole-genome sequencing with bait enrichment techniques to analyze Andes virus (ANDV), a cause of human hantavirus pulmonary syndrome. We used cryopreserved lung tissues from a naturally infected long-tailed colilargo, including early, intermediate, and late cell culture, passages of an ANDV isolate from that animal, and lung tissues from golden hamsters experimentally exposed to that ANDV isolate. The resulting complete genome sequences were subjected to detailed comparative genomic analysis against American orthohantaviruses. We identified four amino acid substitutions related to cell culture adaptation that resulted in attenuation of ANDV in the typically lethal golden hamster animal model of hantavirus pulmonary syndrome. Changes in the ANDV nucleocapsid protein, glycoprotein, and small nonstructural protein open reading frames correlated with mutations typical for ANDV strains associated with increased virulence in the small-animal model. Finally, we identified three amino acid substitutions, two in the small nonstructural protein and one in the glycoprotein, that were only present in the clade of viruses associated with efficient person-to-person transmission. Our results indicate that there are single-nucleotide polymorphisms that could be used to predict strain-specific ANDV virulence and/or transmissibility. IMPORTANCE Several orthohantaviruses cause the zoonotic disease hantavirus pulmonary syndrome (HPS) in the Americas. Among them, HPS caused by Andes virus (ANDV) is of great public health concern because it is associated with the highest case fatality rate (up to 50%). ANDV is also the only orthohantavirus associated with relatively robust evidence of person-to-person transmission. This work reveals nucleotide changes in the ANDV genome that are associated with virulence attenuation in an animal model and increased transmissibility in humans. These findings may pave the way to early severity predictions in future ANDV-caused HPS outbreaks.

Novel viruses in hard ticks collected in the Republic of Korea unveiled by metagenomic high-throughput sequencing analysis.

Ticks are vectors of a wide range of zoonotic viruses of medical and veterinary importance. Recently, metagenomics studies demonstrated that they are also the source of potentially pathogenic novel viruses. During the period from 2015 to 2017, questing ticks were collected by dragging the vegetation from geographically distant locations in the Republic of Korea (ROK) and a target-independent high-throughput sequencing method was utilized to study their virome. A total of seven viruses, including six putative novel viral entities, were identified. Genomic analysis showed that the novel viruses were most closely related to members in the orders Jingchuvirales and Bunyavirales. Phylogenetic reconstruction showed that the Bunyavirales-like viruses grouped in the same clade with other viruses within the Nairovirus and Phlebovirus genera, while the novel Jingchuvirales-like virus grouped together with other viruses within the family Chuviridae. Real-time RT-PCR was used to determine the geographic distribution and prevalence of these viruses in adult ticks. These novel viruses have a wide geographic distribution in the ROK with prevalences ranging from 2% to 18%. Our study expands the knowledge about the composition of the tick virome and highlights the wide diversity of viruses they harbor in the ROK. The discovery of novel viruses associated with ticks in the ROK highlights the need for an active tick-borne disease surveillance program to identify possible reservoirs of putative novel human pathogens.

Chikungunya viruses containing the A226V mutation detected retrospectively in Cameroon form a new geographical subclade.

BACKGROUND: Chikungunya virus (CHIKV) is a re-emerging arbovirus associated with sporadic outbreaks in Cameroon since 2006. Viral whole genomes were generated to analyze the origins of evolutionary lineages, the potential of emergence/re-emergence, and to infer transmission dynamics of recent Cameroon CHIKV outbreak strains. METHODS: Samples collected between 2016 and 2019 during CHIKV outbreaks in Cameroon were screened for CHIKV using reverse transcription PCR (RT-PCR), followed by whole genome sequencing of positive samples. RESULTS: Three coding-complete CHIKV genomes were obtained from samples, which belong to an emerging sub-lineage of the East/Central/South African genotype and formed a monophyletic taxon with previous Central African strains. This clade, which we have named the new Central African clade, appears to be evolving at 3.0 × 10-4 nucleotide substitutions per site per year (95% highest posterior density (HPD) interval of 1.94 × 10-4 to 4.1 × 10-4). Notably, mutations in the envelope proteins (E1-A226V, E2-L210Q, and E2-I211T), which are known to enhance CHIKV adaptability and infectious potential in Aedes albopictus, were present in all strains and mapped to established high-density Ae. albopictus populations. CONCLUSIONS: These new CHIKV strains constitute a conserved genomic pool of an emerging sub-lineage, reflecting a putative vector host adaptation to Ae. albopictus, which has practically displaced Aedes aegypti from select regions of Cameroon.

A Department of Defense Laboratory Consortium Approach to Next Generation Sequencing and Bioinformatics Training for Infectious Disease Surveillance in Kenya.

Epidemics of emerging and re-emerging infectious diseases are a danger to civilian and military populations worldwide. Health security and mitigation of infectious disease threats is a priority of the United States Government and the Department of Defense (DoD). Next generation sequencing (NGS) and Bioinformatics (BI) enhances traditional biosurveillance by providing additional data to understand transmission, identify resistance and virulence factors, make predictions, and update risk assessments. As more and more laboratories adopt NGS and BI technologies they encounter challenges in building local capacity. In addition to choosing the right sequencing platform and approach, considerations must also be made for the complexity of bioinformatics analyses, data storage, as well as personnel and computational requirements. To address these needs, a comprehensive training program was developed covering wet lab and bioinformatics approaches to NGS. The program is meant to be modular and adaptive to meet both common and individualized needs of medical research and public health laboratories across the DoD. The training program was first deployed internationally to the Basic Science Laboratory of the US Army Medical Research Directorate-Africa in Kisumu, Kenya, which is an overseas Lab of the Walter Reed Army Institute of Research (WRAIR). A week-long workshop with intensive focus on targeted sequencing and the bioinformatics of genome assembly (n = 24 participants) was held. Post-workshop self-assessment (completed by 21 participants) noted significant median gains in knowledge domains related to NGS targeted sequencing, bioinformatics for genome assembly, and sequence quality assessment. The participants also reported that the information on study design, sample preparation, sequencing quality control, data quality assessment, reporting, and basic and advanced bioinformatics analysis were the most useful information presented in the training. While longer-term evaluations are planned, the training resulted in significant short-term improvement of a laboratory's self-reported wet lab and bioinformatics capabilities. This framework can be used for future DoD laboratory development in the area of NGS and BI for infectious disease surveillance, ultimately enhancing this global DoD capability.

Prospective Cohort Study of Next-Generation Sequencing as a Diagnostic Modality for Unexplained Encephalitis in Children.

BACKGROUND: Encephalitis is an inflammatory condition of the brain associated with long-term neurologic sequelae and even death in children. Although viruses are often implicated, an etiology is not identified in the majority of cases. Metagenomics-based next-generation sequencing (mNGS) is a high-throughput sequencing technique that can enhance the detection of novel or low-frequency pathogens. METHODS: Hospitalized immunocompetent children aged 6 months to 18 years with encephalitis of unidentified etiology were eligible for enrollment. Demographic, historical, and clinical information was obtained, and residual blood and cerebrospinal fluid (CSF) samples were subjected to mNGS. Pathogens were identified by querying the sequence data against the NCBI GenBank database. RESULTS: Twenty children were enrolled prospectively between 2013 and 2017. mNGS of CSF identified 7 nonhuman nucleic acid sequences of significant frequency in 6 patients, including that of Mycoplasma bovis, parvovirus B19, Neisseria meningitidis, and Balamuthia mandrillaris. mNGS also detected Cladophialophora species, tobacco mosaic virus, and human bocavirus, which were presumed to be contaminants or nonpathogenic organisms. One patient was found to have positive serology results for California encephalitis virus, but mNGS did not detect it. Patients for whom mNGS identified a diagnosis had a significantly higher CSF white blood cell count, a higher CSF protein concentration, and a lower CSF glucose level than patients for whom mNGS did not identify a diagnosis. CONCLUSION: We describe here the results of a prospective cohort analysis to evaluate mNGS as a diagnostic tool for children with unexplained encephalitis. Although mNGS detected multiple nonpathogenic organisms, it also identified multiple pathogens successfully and was most useful in patients with a CSF abnormality.

Target-independent high-throughput sequencing methods provide evidence that already known human viral pathogens play a main role in respiratory infections with unexplained etiology.

Despite the advanced PCR-based assays available, a fraction of the pediatric respiratory infections remain unexplained every epidemic season, and there is a perception that novel viruses might be present in these specimens. We systematically collected samples from a prospective cohort of pediatric patients with respiratory infections, that returned negative results by validated molecular RT-PCR assays, and studied them with a target-independent, high-throughput sequencing-based approach. We also included a matched cohort of children with no symptoms of respiratory infection, as a contrast study population. More than fifty percent of the specimens from the group of patients with unexplained respiratory infections were resolved. However, the higher rate of detection was not due to the presence of novel viruses, but to the identification of well-known viral respiratory pathogens. Our results show that already known viral pathogens are responsible for the majority of cases that remain unexplained after the epidemic season. High-throughput sequencing approaches that use pathogen-specific probes are easier to standardize because they ensure reproducible library enrichment and sequencing. In consequence, these techniques might be desirable from a regulatory standpoint for diagnostic laboratories seeking to benefit from the many advantages of these sequencing technologies.

Lassa virus circulating in Liberia: a retrospective genomic characterisation.

BACKGROUND: An alarming rise in reported Lassa fever cases continues in west Africa. Liberia has the largest reported per capita incidence of Lassa fever cases in the region, but genomic information on the circulating strains is scarce. The aim of this study was to substantially increase the available pool of data to help foster the generation of targeted diagnostics and therapeutics. METHODS: Clinical serum samples collected from 17 positive Lassa fever cases originating from Liberia (16 cases) and Guinea (one case) within the past decade were processed at the Liberian Institute for Biomedical Research using a targeted-enrichment sequencing approach, producing 17 near-complete genomes. An additional 17 Lassa virus sequences (two from Guinea, seven from Liberia, four from Nigeria, and four from Sierra Leone) were generated from viral stocks at the US Centers for Disease Control and Prevention (Atlanta, GA) from samples originating from the Mano River Union (Guinea, Liberia, and Sierra Leone) region and Nigeria. Sequences were compared with existing Lassa virus genomes and published Lassa virus assays. FINDINGS: The 23 new Liberian Lassa virus genomes grouped within two clades (IV.A and IV.B) and were genetically divergent from those circulating elsewhere in west Africa. A time-calibrated phylogeographic analysis incorporating the new genomes suggests Liberia was the entry point of Lassa virus into the Mano River Union region and estimates the introduction to have occurred between 300-350 years ago. A high level of diversity exists between the Liberian Lassa virus genomes. Nucleotide percent difference between Liberian Lassa virus genomes ranged up to 27% in the L segment and 18% in the S segment. The commonly used Lassa Josiah-MGB assay was up to 25% divergent across the target sites when aligned to the Liberian Lassa virus genomes. INTERPRETATION: The large amount of novel genomic diversity of Lassa virus observed in the Liberian cases emphasises the need to match deployed diagnostic capabilities with locally circulating strains and underscores the importance of evaluating cross-lineage protection in the development of vaccines and therapeutics. FUNDING: Defense Biological Product Assurance Office of the US Department of Defense and the Armed Forces Health Surveillance Branch and its Global Emerging Infections Surveillance and Response Section.

Medical countermeasures during the 2018 Ebola virus disease outbreak in the North Kivu and Ituri Provinces of the Democratic Republic of the Congo: a rapid genomic assessment.

BACKGROUND: The real-time generation of information about pathogen genomes has become a vital goal for transmission analysis and characterisation in rapid outbreak responses. In response to the recently established genomic capacity in the Democratic Republic of the Congo, we explored the real-time generation of genomic information at the start of the 2018 Ebola virus disease (EVD) outbreak in North Kivu Province. METHODS: We used targeted-enrichment sequencing to produce two coding-complete Ebola virus genomes 5 days after declaration of the EVD outbreak in North Kivu. Subsequent sequencing efforts yielded an additional 46 genomes. Genomic information was used to assess early transmission, medical countermeasures, and evolution of Ebola virus. FINDINGS: The genomic information demonstrated that the EVD outbreak in the North Kivu and Ituri Provinces was distinct from the 2018 EVD outbreak in Équateur Province of the Democratic Republic of the Congo. Primer and probe mismatches to Ebola virus were identified in silico for all deployed diagnostic PCR assays, with the exception of the Cepheid GeneXpert GP assay. INTERPRETATION: The first two coding-complete genomes provided actionable information in real-time for the deployment of the rVSVΔG-ZEBOV-GP Ebola virus envelope glycoprotein vaccine, available therapeutics, and sequence-based diagnostic assays. Based on the mutations identified in the Ebola virus surface glycoprotein (GP12) observed in all 48 genomes, deployed monoclonal antibody therapeutics (mAb114 and ZMapp) should be efficacious against the circulating Ebola virus variant. Rapid Ebola virus genomic characterisation should be included in routine EVD outbreak response procedures to ascertain efficacy of medical countermeasures. FUNDING: Defense Biological Product Assurance Office.

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.

Persistence and Intra-Host Genetic Evolution of Zika Virus Infection in Symptomatic Adults: A Special View in the Male Reproductive System.

We followed the presence of Zika virus (ZIKV) in four healthy adults (two men and two women), for periods ranging from 78 to 298 days post symptom onset. The patients were evaluated regarding the presence of the virus in different body fluids (blood, saliva, urine and semen), development of immune responses (including antibodies, cytokines and chemokines), and virus genetic variation within samples collected from semen and urine during the infection course. The analysis was focused primarily on the two male patients who shed the virus for up to 158 days after the initial symptoms. ZIKV particles were detected in the spermatozoa cytoplasm and flagella, in immature sperm cells and could also be isolated from semen in cell culture, confirming that the virus is able to preserve integrity and infectivity during replication in the male reproductive system (MRS). Despite the damage caused by ZIKV infection within the MRS, our data showed that ZIKV infection did not result in infertility at least in one of the male patients. This patient was able to conceive a child after the infection. We also detected alterations in the male genital cytokine milieu, which could play an important role in the replication and transmission of the virus which could considerably increase the risk of ZIKV sexual spread. In addition, full genome ZIKV sequences were obtained from several samples (mainly semen), which allowed us to monitor the evolution of the virus within a patient during the infection course. We observed genetic changes over time in consensus sequences and lower frequency intra-host single nucleotide variants (iSNV), that suggested independent compartmentalization of ZIKV populations in the reproductive and urinary systems. Altogether, the present observations confirm the risks associated with the long-term replication and shedding of ZIKV in the MRS and help to elucidate patterns of intra-host genetic evolution during long term replication of the virus.

Beatrice Hill Virus Represents a Novel Species in the Genus Tibrovirus (Mononegavirales: Rhabdoviridae).

The rhabdoviral genus Tibrovirus currently has three official members assigned to two species: Bivens Arm virus and Tibrogargan virus (species Tibrogargan tibrovirus) and Coastal Plains virus (species Coastal Plains tibrovirus). Here, we report the complete genome sequence of a new putative member of this genus, Beatrice Hill virus. Although relatively closely related to the three classified viruses, Beatrice Hill virus represents a novel Tibrovirus species.

First report of naturally infected Aedes aegypti with chikungunya virus genotype ECSA in the Americas.

BACKGROUND: The worldwide expansion of new emergent arboviruses such as Chikungunya and Zika reinforces the importance in understanding the role of mosquito species in spreading these pathogens in affected regions. This knowledge is essential for developing effective programs based on species specificity to avoid the establishment of endemic transmission cycles sustained by the identified local vectors. Although the first autochthonous transmission of Chikungunya virus was described in 2014 in the north of Brazil, the main outbreaks were reported in 2015 and 2016 in the northeast of Brazil. METHODOLOGY/PRINCIPAL FINDINGS: During 5 days of February 2016, we collected mosquitoes in homes of 6 neighborhoods of Aracaju city, the capital of Sergipe state. Four mosquito species were identified but Culex quinquefasciatus and Aedes aegypti were the most abundant. Field-caught mosquitoes were tested for Chikungunya (CHIKV), Zika (ZIKV) and Dengue viruses (DENV) by qRT-PCR and one CHIKV-infected Ae. aegypti female was detected. The complete sequence of CHIKV genome was obtained from this sample and phylogenetic analysis revealed that this isolate belongs to the East-Central-South-African (ECSA) genotype. CONCLUSIONS: Our study describes the first identification of a naturally CHIKV-infected Ae. aegypti in Brazil and the first report of a CHIKV from ECSA genotype identified in this species in the Americas. These findings support the notion of Ae. aegypti being a vector involved in CHIKV outbreaks in northeast of Brazil.

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

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

Complete Genome Sequences of Five Zika Virus Isolates.

Zika virus is an emerging human pathogen of great concern due to putative links to microcephaly and Guillain-Barre syndrome. Here, we report the complete genomes, including the 5' and 3' untranslated regions, of five Zika virus isolates, one from the Asian lineage and four from the African lineage.