Recurrent SARS-CoV-2 mutations at Spike D796 evade antibodies from pre-Omicron convalescent and vaccinated subjects

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages of the Omicron variant rapidly became dominant in early 2022 and frequently cause human infections despite vaccination or prior infection with other variants. In addition to antibody-evading mutations in the receptor-binding domain, Omicron features amino acid mutations elsewhere in the Spike protein; however, their effects generally remain ill defined. The Spike D796Y substitution is present in all Omicron sub-variants and occurs at the same site as a mutation (D796H) selected during viral evolution in a chronically infected patient. Here, we map antibody reactivity to a linear epitope in the Spike protein overlapping position 796. We show that antibodies binding this region arise in pre-Omicron SARS-CoV-2 convalescent and vaccinated subjects but that both D796Y and D796H abrogate their binding. These results suggest that D796Y contributes to the fitness of Omicron in hosts with pre-existing immunity to other variants of SARS-CoV-2 by evading antibodies targeting this site.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved substantially through the coronavirus disease 2019 (COVID-19) pandemic: understanding the drivers and consequences of this evolution is essential for projecting the course of the pandemic and developing new countermeasures. Here, we study the immunological effects of a particular mutation present in the Spike protein of all Omicron strains and find that it prevents the efficient binding of a class of antibodies raised by pre-Omicron vaccination and infection. These findings reveal a novel consequence of a poorly understood Omicron mutation and shed light on the drivers and effects of SARS-CoV-2 evolution.

Towards a post-pandemic future for global pathogen genome sequencing

Pathogen genome sequencing has become a routine part of our response to active outbreaks of infectious disease and should be an important part of our preparations for future epidemics. In this Essay, we discuss the innovations that have enabled routine pathogen genome sequencing, as well as how genome sequences can be used to understand and control the spread of infectious disease. We also explore the impact of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic on the field of pathogen genomics and outline the challenges we must address to further improve the utility of pathogen genome sequencing in the future.

Virome-wide detection of natural infection events and the associated antibody dynamics using longitudinal highly-multiplexed serology

Current methods for detecting infections either require a sample collected from an actively infected site, are limited in the number of agents they can query, and/or yield no information on the immune response. Here we present an approach that uses temporally coordinated changes in highly-multiplexed antibody measurements from longitudinal blood samples to monitor infection events at sub-species resolution across the human virome. In a longitudinally-sampled cohort of South African adolescents representing >100 person-years, we identify >650 events across 48 virus species and observe strong epidemic effects, including high-incidence waves of Aichivirus A and the D68 subtype of Enterovirus D earlier than their widespread circulation was appreciated. In separate cohorts of adults who were sampled at higher frequency using self-collected dried blood spots, we show that such events temporally correlate with symptoms and transient inflammatory biomarker elevations, and observe the responding antibodies to persist for periods ranging from ≤1 week to >5 years. Our approach generates a rich view of viral/host dynamics, supporting novel studies in immunology and epidemiology.

PepSeq: a fully in vitro platform for highly multiplexed serology using customizable DNA-barcoded peptide libraries

PepSeq is an in vitro platform for building and conducting highly multiplexed proteomic assays against customizable targets by using DNA-barcoded peptides. Starting with a pool of DNA oligonucleotides encoding peptides of interest, this protocol outlines a fully in vitro and massively parallel procedure for synthesizing the encoded peptides and covalently linking each to a corresponding cDNA tag. The resulting libraries of peptide/DNA conjugates can be used for highly multiplexed assays that leverage high-throughput sequencing to profile the binding or enzymatic specificities of proteins of interest. Here, we describe the implementation of PepSeq for fast and cost-effective epitope-level analysis of antibody reactivity across hundreds of thousands of peptides from <1 µl of serum or plasma input. This protocol includes the design of the DNA oligonucleotide library, synthesis of DNA-barcoded peptide constructs, binding of constructs to sample, preparation for sequencing and data analysis. Implemented in this way, PepSeq can be used for a number of applications, including fine-scale mapping of antibody epitopes and determining a subject's pathogen exposure history. The protocol is divided into two main sections: (i) design and synthesis of DNA-barcoded peptide libraries and (ii) use of libraries for highly multiplexed serology. Once oligonucleotide templates are in hand, library synthesis takes 1-2 weeks and can provide enough material for hundreds to thousands of assays. Serological assays can be conducted in 96-well plates and generate sequencing data within a further ~4 d. A suite of software tools, including the PepSIRF package, are made available to facilitate the design of PepSeq libraries and analysis of assay data.

Highly Multiplexed Serology for Nonhuman Mammals

Emerging infectious diseases represent a serious and ongoing threat to humans. Most emerging viruses are maintained in stable relationships with other species of animals, and their emergence within the human population results from cross-species transmission. Therefore, if we want to be prepared for the next emerging virus, we need to broadly characterize the diversity and ecology of viruses currently infecting other animals (i.e., the animal virosphere). High-throughput metagenomic sequencing has accelerated the pace of virus discovery. However, molecular assays can detect only active infections and only if virus is present within the sampled fluid or tissue at the time of collection. In contrast, serological assays measure long-lived antibody responses to infections, which can be detected within the blood, regardless of the infected tissues. Therefore, serological assays can provide a complementary approach for understanding the circulation of viruses, and while serological assays have historically been limited in scope, recent advancements allow thousands to hundreds of thousands of antigens to be assessed simultaneously using <1 μL of blood (i.e., highly multiplexed serology). The application of highly multiplexed serology for the characterization of the animal virosphere is dependent on the availability of reagents that can be used to capture or label antibodies of interest. Here, we evaluate the utility of commercial immunoglobulin-binding proteins (protein A and protein G) to enable highly multiplexed serology in 25 species of nonhuman mammals, and we describe a competitive fluorescence-linked immunosorbent assay (FLISA) that can be used as an initial screen for choosing the most appropriate capture protein for a given host species. IMPORTANCE Antibodies are generated in response to infections with viruses and other pathogens, and they help protect against future exposures. Mature antibodies are long lived, are highly specific, and can bind to their protein targets with high affinity. Thus, antibodies can also provide information about an individual's history of viral exposures, which has important applications for understanding the epidemiology and etiology of disease. In recent years, there have been large advances in the available methods for broadly characterizing antibody-binding profiles, but thus far, these have been utilized primarily with human samples only. Here, we demonstrate that commercial antibody-binding reagents can facilitate modern antibody assays for a wide variety of mammalian species, and we describe an inexpensive and fast approach for choosing the best reagent for each animal species. By studying antibody-binding profiles in captive and wild animals, we can better understand the distribution and prevalence of viruses that could spill over into humans.

COVID-19 vaccination elicits an evolving, cross-reactive antibody response to epitopes conserved with endemic coronavirus spike proteins

The COVID-19 pandemic has triggered the first widespread vaccination campaign against a coronavirus. Many vaccinated subjects are previously naive to SARS-CoV-2; however, almost all have previously encountered other coronaviruses (CoVs), and the role of this immunity in shaping the vaccine response remains uncharacterized. Here, we use longitudinal samples and highly multiplexed serology to identify mRNA-1273 vaccine-induced antibody responses against a range of CoV Spike epitopes, in both phylogenetically conserved and non-conserved regions. Whereas reactivity to SARS-CoV-2 epitopes shows a delayed but progressive increase following vaccination, we observe distinct kinetics for the endemic CoV homologs at conserved sites in Spike S2: these become detectable sooner and decay at later time points. Using homolog-specific antibody depletion and alanine-substitution experiments, we show that these distinct trajectories reflect an evolving cross-reactive response that can distinguish rare, polymorphic residues within these epitopes. Our results reveal mechanisms for the formation of antibodies with broad reactivity against CoVs.

COVID-19 vaccination recruits and matures cross-reactive antibodies to conserved epitopes in endemic coronavirus Spike proteins.. [PMID: 35118479 PMCID: PMC8811912 DOI: 10.1101/2022.01.24.22269542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The COVID-19 pandemic has triggered the first widespread vaccination campaign against a coronavirus. Most vaccinated subjects are naïve to SARS-CoV-2, however almost all have previously encountered other coronaviruses (CoVs) and the role of this immunity in shaping the vaccine response remains uncharacterized. Here we use longitudinal samples and highly-multiplexed serology to identify mRNA-1273 vaccine-induced antibody responses against a range of CoV Spike epitopes and in both phylogenetically conserved and non-conserved regions. Whereas reactivity to SARS-CoV-2 epitopes showed a delayed but progressive increase following vaccination, we observed distinct kinetics for the endemic CoV homologs at two conserved sites in Spike S2: these became detectable sooner, and decayed at later timepoints. Using homolog-specific depletion and alanine-substitution experiments, we show that these distinctly-evolving specificities result from cross-reactive antibodies as they mature against rare, polymorphic residues within these epitopes. Our results reveal mechanisms for the formation of antibodies with broad reactivity against CoVs.

Comparative Genomics Analyses Support the Reclassification of Bisgaard Taxon 40 as Mergibacter gen. nov., With Mergibacter septicus sp. nov. as Type Species: Novel Insights Into the Phylogeny and Virulence Factors of a Pasteurellaceae Family Member Associated With Mortality Events in Seabirds

The Pasteurellaceae family has been associated with fatal diseases in numerous avian species. Several new taxa within this family, including Bisgaard taxon 40, have been recently described in wild birds, but their genomic characteristics and pathogenicity are not well understood. We isolated Bisgaard taxon 40 from four species of seabirds, including one sampled during a mass, multi-species mortality event in Florida, United States. Here, we present a comprehensive phenotypic and genetic characterization of Bisgaard taxon 40 and comparative genomic analysis with reference strains from the Pasteurellaceae family, aiming at determining its phylogenetic position, antimicrobial susceptibility profile, and identifying putative virulence factors. In silico multilocus sequence-based and whole-genome-based phylogenetic analysis clustered all Bisgaard taxon 40 strains together on a distinct branch separated from the other members of the Pasteurellaceae family, indicating that Bisgaard taxon 40 could represent a new genus. These findings were further supported by protein similarity analyses using the concatenation of 31 conserved proteins and other taxonomic approaches such as the percentage of conserved protein test. Additionally, several putative virulence factors were identified, including those associated with adhesion (capsule, ompA, ompH) and colonization (exbD, fur, galU, galE, lpxA, lpxC, and kdsA) of the host and a cytolethal distending toxin (cdt), which may have played a role in disease development leading to the mortality event. Considerably low minimum inhibitory concentrations (MICs) were found for all the drugs tested, in concordance with the absence of antimicrobial resistance genes in these genomes. The novel findings of this study highlight genomic and phenotypic characteristics of this bacterium, providing insights into genome evolution and pathogenicity. We propose a reclassification of these organisms within the Pasteurellaceae family, designated as Mergibacter gen. nov., with Mergibacter septicus sp. nov. as the type species. The type strain is Mergibacter septicus A25201^T (=DSM 112696).

Multiplexed Simian Immunodeficiency Virus-Specific Paired RNA-Guided Cas9 Nickases Inactivate Proviral DNA

Human and simian immunodeficiency virus (HIV and SIV) infections establish lifelong reservoirs of cells harboring an integrated proviral genome. Genome editing CRISPR-associated Cas9 nucleases, combined with SIV-specific guiding RNA (gRNA) molecules, inactivate integrated provirus DNA in vitro and in animal models. We generated RNA-guided Cas9 nucleases (RGNu) and nickases (RGNi) targeting conserved SIV regions with no homology in the human or rhesus macaque genome. Assays in cells cotransfected with SIV provirus and plasmids coding for RGNus identified SIV long terminal repeat (LTR), trans-activation response (TAR) element, and ribosome slip site (RSS) regions as the most effective at virus suppression; RGNi targeting these regions inhibited virus production significantly. Multiplex plasmids that coexpressed these three RGNu (Nu3), or six (three pairs) RGNi (Ni6), were more efficient at virus suppression than any combination of individual RGNu and RGNi plasmids. Both Nu3 and Ni6 plasmids were tested in lymphoid cells chronically infected with SIVmac239, and whole-genome sequencing was used to determine on- and off-target mutations. Treatment with these all-in-one plasmids resulted in similar levels of mutations of viral sequences from the cellular genome; Nu3 induced indels at the 3 SIV-specific sites, whereas for Ni6 indels were present at the LTR and TAR sites. Levels of off-target effects detected by two different algorithms were indistinguishable from background mutations. In summary, we demonstrate that Cas9 nickase in association with gRNA pairs can specifically eliminate parts of the integrated provirus DNA; also, we show that careful design of an all-in-one plasmid coding for 3 gRNAs and Cas9 nuclease inhibits SIV production with undetectable off-target mutations, making these tools a desirable prospect for moving into animal studies. IMPORTANCE Our approach to HIV cure, utilizing the translatable SIV/rhesus macaque model system, aims at provirus inactivation and its removal with the least possible off-target side effects. We developed single molecules that delivered either three truncated SIV-specific gRNAs along with Cas9 nuclease or three pairs of SIV-specific gRNAs (six individual gRNAs) along with Cas9 nickase to enhance efficacy of on-target mutagenesis. Whole-genome sequencing demonstrated effective SIV sequence mutation and inactivation and the absence of demonstrable off-target mutations. These results open the possibility to employ Cas9 variants that introduce single-strand DNA breaks to eliminate integrated proviral DNA.

Genomic signatures for predicting the zoonotic potential of novel viruses

Powered by metagenomics, viral discovery is outpacing our capacity for the downstream characterization needed to fully assess zoonotic potential. A study published in PLOS Biology uses machine learning to prioritize novel viruses based only on genomic signatures.

Epitope-resolved profiling of the SARS-CoV-2 antibody response identifies cross-reactivity with endemic human coronaviruses

The SARS-CoV-2 proteome shares regions of conservation with endemic human Coronaviruses (CoVs), but it remains unknown to what extent these may be cross-recognized by the antibody response. Here, we study cross-reactivity using a highly-multiplexed peptide assay (‘PepSeq’) to generate an epitope-resolved view of IgG reactivity across all human CoVs in both COVID-19 convalescent and negative donors. PepSeq resolves epitopes across the SARS-CoV-2 Spike and Nucleocapsid proteins that are commonly targeted in convalescent donors, including several sites also recognized in some uninfected controls. By comparing patterns of homologous reactivity between CoVs and using targeted antibody-depletion experiments, we demonstrate that SARS-CoV-2 elicits antibodies that cross-recognize pandemic and endemic CoV antigens at two Spike S2 subunit epitopes. We further show that these cross-reactive antibodies preferentially bind endemic homologs. Our findings highlight sites at which the SARS-CoV-2 response appears to be shaped by previous CoV exposures and which have potential to raise broadly-neutralizing responses.

Inference of Nipah virus evolution, 1999-2015

Despite near-annual human outbreaks of Nipah virus (NiV) disease in Bangladesh, typically due to individual spillover events from the local bat population, only twenty whole-genome NiV sequences exist from humans and ten from bats. NiV whole-genome sequences from annual outbreaks have been challenging to generate, primarily due to the low viral load in human throat swab and serum specimens. Here, we used targeted enrichment with custom NiV-specific probes and generated thirty-five additional unique full-length genomic sequences directly from human specimens and viral isolates. We inferred the temporal and geographic evolutionary history of NiV in Bangladesh and expanded a tool to visualize NiV spatio-temporal spread from a Bayesian continuous diffusion analysis. We observed that strains from Bangladesh segregated into two distinct clades that have intermingled geographically in Bangladesh over time and space. As these clades expanded geographically and temporally, we did not observe evidence for significant branch and site-specific selection, except for a single site in the Henipavirus L polymerase. However, the Bangladesh 1 and 2 clades are differentiated by mutations initially occurring in the polymerase, with additional mutations accumulating in the N, G, F, P, and L genes on external branches. Modeling the historic geographical and temporal spread demonstrates that while widespread, NiV does not exhibit significant genetic variation in Bangladesh. Thus, future public health measures should address whether NiV within in the bat population also exhibits comparable genetic variation, if zoonotic transmission results in a genetic bottleneck and if surveillance techniques are detecting only a subset of NiV.

Epitope-resolved profiling of the SARS-CoV-2 antibody response identifies cross-reactivity with endemic human coronaviruses

The SARS-CoV-2 proteome shares regions of conservation with endemic human coronaviruses (CoVs), but it remains unknown to what extent these may be cross-recognized by the antibody response. Here, we study cross-reactivity using a highly multiplexed peptide assay (PepSeq) to generate an epitope-resolved view of IgG reactivity across all human CoVs in both COVID-19 convalescent and negative donors. PepSeq resolves epitopes across the SARS-CoV-2 Spike and Nucleocapsid proteins that are commonly targeted in convalescent donors, including several sites also recognized in some uninfected controls. By comparing patterns of homologous reactivity between CoVs and using targeted antibody-depletion experiments, we demonstrate that SARS-CoV-2 elicits antibodies that cross-recognize pandemic and endemic CoV antigens at two Spike S2 subunit epitopes. We further show that these cross-reactive antibodies preferentially bind endemic homologs. Our findings highlight sites at which the SARS-CoV-2 response appears to be shaped by previous CoV exposures and which have the potential to raise broadly neutralizing responses.

PepSeq enables fully in vitro, highly multiplexed peptide-based antibody assays

Epitope mapping shows preexisting antibody reactivity to SARS-CoV-2 antigens

Antibodies cross-recognize endemic and pandemic antigens in the Spike S2 subunit

Cross-reactive antibodies raised by SARS-CoV-2 preferentially bind endemic homologs

Evolution of Antibiotic Resistance in Surrogates of Francisella tularensis (LVS and Francisella novicida): Effects on Biofilm Formation and Fitness

Francisella tularensis, the causative agent of tularemia, is capable of causing disease in a multitude of mammals and remains a formidable human pathogen due to a high morbidity, low infectious dose, lack of a FDA approved vaccine, and ease of aerosolization. For these reasons, there is concern over the use of F. tularensis as a biological weapon, and, therefore, it has been classified as a Tier 1 select agent. Fluoroquinolones and aminoglycosides often serve as the first line of defense for treatment of tularemia. However, high levels of resistance to these antibiotics has been observed in gram-negative bacteria in recent years, and naturally derived resistant Francisella strains have been described in the literature. The acquisition of antibiotic resistance, either natural or engineered, presents a challenge for the development of medical countermeasures. In this study, we generated a surrogate panel of antibiotic resistant F. novicida and Live Vaccine Strain (LVS) by selection in the presence of antibiotics and characterized their growth, biofilm capacity, and fitness. These experiments were carried out in an effort to (1) assess the fitness of resistant strains; and (2) identify new targets to investigate for the development of vaccines or therapeutics. All strains exhibited a high level of resistance to either ciprofloxacin or streptomycin, a fluoroquinolone and aminoglycoside, respectively. Whole genome sequencing of this panel revealed both on-pathway and off-pathway mutations, with more mutations arising in LVS. For F. novicida, we observed decreased biofilm formation for all ciprofloxacin resistant strains compared to wild-type, while streptomycin resistant isolates were unaffected in biofilm capacity. The fitness of representative antibiotic resistant strains was assessed in vitro in murine macrophage-like cell lines, and also in vivo in a murine model of pneumonic infection. These experiments revealed that mutations obtained by these methods led to nearly all ciprofloxacin resistant Francisella strains tested being completely attenuated while mild attenuation was observed in streptomycin resistant strains. This study is one of the few to examine the link between acquired antibiotic resistance and fitness in Francisella spp., as well as enable the discovery of new targets for medical countermeasure development.

An Early Pandemic Analysis of SARS-CoV-2 Population Structure and Dynamics in Arizona

In December of 2019, a novel coronavirus, SARS-CoV-2, emerged in the city of Wuhan, China, causing severe morbidity and mortality. Since then, the virus has swept across the globe, causing millions of confirmed infections and hundreds of thousands of deaths. To better understand the nature of the pandemic and the introduction and spread of the virus in Arizona, we sequenced viral genomes from clinical samples tested at the TGen North Clinical Laboratory, the Arizona Department of Health Services, and those collected as part of community surveillance projects at Arizona State University and the University of Arizona. Phylogenetic analysis of 84 genomes from across Arizona revealed a minimum of 11 distinct introductions inferred to have occurred during February and March. We show that >80% of our sequences descend from strains that were initially circulating widely in Europe but have since dominated the outbreak in the United States. In addition, we show that the first reported case of community transmission in Arizona descended from the Washington state outbreak that was discovered in late February. Notably, none of the observed transmission clusters are epidemiologically linked to the original travel-related case in the state, suggesting successful early isolation and quarantine. Finally, we use molecular clock analyses to demonstrate a lack of identifiable, widespread cryptic transmission in Arizona prior to the middle of February 2020.IMPORTANCE As the COVID-19 pandemic swept across the United States, there was great differential impact on local and regional communities. One of the earliest and hardest hit regions was in New York, while at the same time Arizona (for example) had low incidence. That situation has changed dramatically, with Arizona now having the highest rate of disease increase in the country. Understanding the roots of the pandemic during the initial months is essential as the pandemic continues and reaches new heights. Genomic analysis and phylogenetic modeling of SARS-COV-2 in Arizona can help to reconstruct population composition and predict the earliest undetected introductions. This foundational work represents the basis for future analysis and understanding as the pandemic continues.

Epitope-resolved profiling of the SARS-CoV-2 antibody response identifies cross-reactivity with an endemic human CoV

A high-resolution understanding of the antibody response to SARS-CoV-2 is important for the design of effective diagnostics, vaccines and therapeutics. However, SARS-CoV-2 antibody epitopes remain largely uncharacterized, and it is unknown whether and how the response may cross-react with related viruses. Here, we use a multiplexed peptide assay ('PepSeq') to generate an epitope-resolved view of reactivity across all human coronaviruses. PepSeq accurately detects SARS-CoV-2 exposure and resolves epitopes across the Spike and Nucleocapsid proteins. Two of these represent recurrent reactivities to conserved, functionally-important sites in the Spike S2 subunit, regions that we show are also targeted for the endemic coronaviruses in pre-pandemic controls. At one of these sites, we demonstrate that the SARS-CoV-2 response strongly and recurrently cross-reacts with the endemic virus hCoV-OC43. Our analyses reveal new diagnostic and therapeutic targets, including a site at which SARS-CoV-2 may recruit common pre-existing antibodies and with the potential for broadly-neutralizing responses.

Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity

The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2's unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.

Reproducibly sampling SARS-CoV-2 genomes across time, geography, and viral diversity

The COVID-19 pandemic has led to a rapid accumulation of SARS-CoV-2 genomes, enabling genomic epidemiology on local and global scales. Collections of genomes from resources such as GISAID must be subsampled to enable computationally feasible phylogenetic and other analyses. We present genome-sampler, a software package that supports sampling collections of viral genomes across multiple axes including time of genome isolation, location of genome isolation, and viral diversity. The software is modular in design so that these or future sampling approaches can be applied independently and combined (or replaced with a random sampling approach) to facilitate custom workflows and benchmarking. genome-sampler is written as a QIIME 2 plugin, ensuring that its application is fully reproducible through QIIME 2’s unique retrospective data provenance tracking system. genome-sampler can be installed in a conda environment on macOS or Linux systems. A complete default pipeline is available through a Snakemake workflow, so subsampling can be achieved using a single command. genome-sampler is open source, free for all to use, and available at https://caporasolab.us/genome-sampler. We hope that this will facilitate SARS-CoV-2 research and support evaluation of viral genome sampling approaches for genomic epidemiology.

Persistence of Brucella abortus lineages revealed by genomic characterization and phylodynamic analysis

Brucellosis, caused by Brucella abortus, is a major disease of cattle and humans worldwide distributed. Eradication and control of the disease has been difficult in Central and South America, Central Asia, the Mediterranean and the Middle East. Epidemiological strategies combined with phylogenetic methods provide the high-resolution power needed to study relationships between surveillance data and pathogen population dynamics, using genetic diversity and spatiotemporal distributions. This information is crucial for prevention and control of disease spreading at a local and worldwide level. In Costa Rica (CR), the disease was first reported at the beginning of the 20th century and has not been controlled despite many efforts. We characterized 188 B. abortus isolates from CR recovered from cattle, humans and water buffalo, from 2003 to 2018, and whole genome sequencing (WGS) was performed in 95 of them. They were also assessed based on geographic origin, date of introduction, and phylogenetic associations in a worldwide and national context. Our results show circulation of five B. abortus lineages (I to V) in CR, phylogenetically related to isolates from the United States, United Kingdom, and South America. Lineage I was dominant and probably introduced at the end of the 19th century. Lineage II, represented by a single isolate from a water buffalo, clustered with a Colombian sample, and was likely introduced after 1845. Lineages III and IV were likely introduced during the early 2000s. Fourteen isolates from humans were found within the same lineage (lineage I) regardless of their geographic origin within the country. The main CR lineages, introduced more than 100 years ago, are widely spread throughout the country, in contrast to new introductions that seemed to be more geographically restricted. Following the brucellosis prevalence and the farming practices of several middle- and low-income countries, similar scenarios could be found in other regions worldwide.

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.

Travel Surveillance and Genomics Uncover a Hidden Zika Outbreak during the Waning Epidemic

The Zika epidemic in the Americas has challenged surveillance and control. As the epidemic appears to be waning, it is unclear whether transmission is still ongoing, which is exacerbated by discrepancies in reporting. To uncover locations with lingering outbreaks, we investigated travel-associated Zika cases to identify transmission not captured by reporting. We uncovered an unreported outbreak in Cuba during 2017, a year after peak transmission in neighboring islands. By sequencing Zika virus, we show that the establishment of the virus was delayed by a year and that the ensuing outbreak was sparked by long-lived lineages of Zika virus from other Caribbean islands. Our data suggest that, although mosquito control in Cuba may initially have been effective at mitigating Zika virus transmission, such measures need to be maintained to be effective. Our study highlights how Zika virus may still be "silently" spreading and provides a framework for understanding outbreak dynamics. VIDEO ABSTRACT.

A novel host-adapted strain of Salmonella Typhimurium causes renal disease in olive ridley turtles (Lepidochelys olivacea) in the Pacific

Salmonella spp. are frequently shed by wildlife including turtles, but S. enterica subsp. enterica serovar Typhimurium or lesions associated with Salmonella are rare in turtles. Between 1996 and 2016, we necropsied 127 apparently healthy pelagic olive ridley turtles (Lepidochelys olivacea) that died from drowning bycatch in fisheries and 44 live or freshly dead stranded turtles from the west coast of North and Central America and Hawaii. Seven percent (9/127) of pelagic and 47% (21/44) of stranded turtles had renal granulomas associated with S. Typhimurium. Stranded animals were 12 times more likely than pelagic animals to have Salmonella-induced nephritis suggesting that Salmonella may have been a contributing cause of stranding. S. Typhimurium was the only Salmonella serovar detected in L. olivacea, and phylogenetic analysis from whole genome sequencing showed that the isolates from L. olivacea formed a single clade distinct from other S. Typhimurium. Molecular clock analysis revealed that this novel clade may have originated as recently as a few decades ago. The phylogenetic lineage leading to this group is enriched for non-synonymous changes within the genomic area of Salmonella pathogenicity island 1 suggesting that these genes are important for host adaptation.

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

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 (GP₁₂) 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.

Precision epidemiology for infectious disease control

Advances in genomics and computing are transforming the capacity for the characterization of biological systems, and researchers are now poised for a precision-focused transformation in the way they prepare for, and respond to, infectious diseases. This includes the use of genome-based approaches to inform molecular diagnosis and individual-level treatment regimens. In addition, advances in the speed and granularity of pathogen genome generation have improved the capability to track and understand pathogen transmission, leading to potential improvements in the design and implementation of population-level public health interventions. In this Perspective, we outline several trends that are driving the development of precision epidemiology of infectious disease and their implications for scientists' ability to respond to outbreaks.

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.

Persistence of Ebola virus after the end of widespread transmission in Liberia: an outbreak report

BACKGROUND

Outbreak response efforts for the 2014-15 Ebola virus disease epidemic in west Africa brought widespread transmission to an end. However, subsequent clusters of infection have occurred in the region. An Ebola virus disease cluster in Liberia in November, 2015, that was identified after a 15-year-old boy tested positive for Ebola virus infection in Monrovia, raised the possibility of transmission from a persistently infected individual.

METHODS

Case investigations were done to ascertain previous contact with cases of Ebola virus disease or infection with Ebola virus. Molecular investigations on blood samples explored a potential linkage between Ebola virus isolated from cases in this November, 2015, cluster and epidemiologically linked cases from the 2014-15 west African outbreak, according to the national case database.

FINDINGS

The cluster investigated was the family of the index case (mother, father, three siblings). Ebola virus genomes assembled from two cases in the November, 2015, cluster, and an epidemiologically linked Ebola virus disease case in July, 2014, were phylogenetically related within the LB5 sublineage that circulated in Liberia starting around August, 2014. Partial genomes from two additional individuals, one from each cluster, were also consistent with placement in the LB5 sublineage. Sequencing data indicate infection with a lineage of the virus from a former transmission chain in the country. Based on serology and epidemiological and genomic data, the most plausible scenario is that a female case in the November, 2015, cluster survived Ebola virus disease in 2014, had viral persistence or recurrent disease, and transmitted the virus to three family members a year later.

INTERPRETATION

Investigation of the source of infection for the November, 2015, cluster provides evidence of Ebola virus persistence and highlights the risk for outbreaks after interruption of active transmission. These findings underscore the need for focused prevention efforts among survivors and sustained capacity to rapidly detect and respond to new Ebola virus disease cases to prevent recurrence of a widespread outbreak.

FUNDING

US Centers for Disease Control and Prevention, Defense Threat Reduction Agency, and WHO.

Countering Zika Virus: The USAMRIID Response

The United States Army Medical Research Institute of Infectious Diseases (USAMRIID) possesses an array of expertise in diverse capabilities for the characterization of emerging infectious diseases from the pathogen itself to human or animal infection models. The recent Zika virus (ZIKV) outbreak was a challenge and an opportunity to put these capabilities to work as a cohesive unit to quickly respond to a rapidly developing threat. Next-generation sequencing was used to characterize virus stocks and to understand the introduction and spread of ZIKV in the United States. High Content Imaging was used to establish a High Content Screening process to evaluate antiviral therapies. Functional genomics was used to identify critical host factors for ZIKV infection. An animal model using the temporal blockade of IFN-I in immunocompetent laboratory mice was investigated in conjunction with Positron Emission Tomography to study ZIKV. Correlative light and electron microscopy was used to examine ZIKV interaction with host cells in culture and infected animals. A quantitative mass spectrometry approach was used to examine the protein and metabolite type or concentration changes that occur during ZIKV infection in blood, cells, and tissues. Multiplex fluorescence in situ hybridization was used to confirm ZIKV replication in mouse and NHP tissues. The integrated rapid response approach developed at USAMRIID presented in this review was successfully applied and provides a new template pathway to follow if a new biological threat emerges. This streamlined approach will increase the likelihood that novel medical countermeasures could be rapidly developed, evaluated, and translated into the clinic.

An attenuated Machupo virus with a disrupted L-segment intergenic region protects guinea pigs against lethal Guanarito virus infection

Machupo virus (MACV) is a New World (NW) arenavirus and causative agent of Bolivian hemorrhagic fever (HF). Here, we identified a variant of MACV strain Carvallo termed Car91 that was attenuated in guinea pigs. Infection of guinea pigs with an earlier passage of Carvallo, termed Car68, resulted in a lethal disease with a 63% mortality rate. Sequencing analysis revealed that compared to Car68, Car91 had a 35 nucleotide (nt) deletion and a point mutation within the L-segment intergenic region (IGR), and three silent changes in the polymerase gene that did not impact amino acid coding. No changes were found on the S-segment. Because it was apathogenic, we determined if Car91 could protect guinea pigs against Guanarito virus (GTOV), a distantly related NW arenavirus. While naïve animals succumbed to GTOV infection, 88% of the Car91-exposed guinea pigs were protected. These findings indicate that attenuated MACV vaccines can provide heterologous protection against NW arenaviruses. The disruption in the L-segment IGR, including a single point mutant and 35 nt partial deletion, were the only major variance detected between virulent and avirulent isolates, implicating its role in attenuation. Overall, our data support the development of live-attenuated arenaviruses as broadly protective pan-arenavirus vaccines.

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.

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.

Signals of selection in conditionally expressed genes in the diversification of three horned beetle species

Species radiations may be facilitated by phenotypic differences already present within populations, such as those arising through sex-specific development or developmental processes biased towards particular reproductive or trophic morphs. We sought to test this hypothesis by utilizing a comparative transcriptomic approach to contrast among- and within-species differentiation using three horned beetle species in the genus Onthophagus. These three species exhibit differences along three phenotypic axes reflective of much of the interspecific diversity present within the genus: horn location, polarity of sexual dimorphism and degree of nutritional sensitivity. Our approach combined de novo transcript assembly, assessment of amino acid substitutions (dN/dS) across orthologous gene pairs and integration of gene function and conditional gene expression data. We identified 17 genes across the three species pairs related to axis patterning, development and metabolism with dN/dS > 1 and detected elevated dN/dS in genes related to metabolism and biosynthesis in the most closely related species pair, which is characterized by a loss of nutritional polyphenism and a reversal of sexual dimorphism. Further, we found that genes that are conditionally expressed (i.e. as a function of sex, nutrition or body region) within one of our focal species also showed significantly stronger signals of positive or relaxed purifying selection between species divergent along the same morphological axis (i.e. polarity of sexual dimorphism, degree of nutritional sensitivity or location of horns). Our findings thus reveal a positive relationship between intraspecific differentiation due to condition-specific development and genetic divergences among species.

Virus genomes reveal factors that spread and sustained the Ebola epidemic

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

A Multicomponent Animal Virus Isolated from Mosquitoes

RNA viruses exhibit a variety of genome organization strategies, including multicomponent genomes in which each segment is packaged separately. Although multicomponent genomes are common among viruses infecting plants and fungi, their prevalence among those infecting animals remains unclear. We characterize a multicomponent RNA virus isolated from mosquitoes, designated Guaico Culex virus (GCXV). GCXV belongs to a diverse clade of segmented viruses (Jingmenvirus) related to the prototypically unsegmented Flaviviridae. The GCXV genome comprises five segments, each of which appears to be separately packaged. The smallest segment is not required for replication, and its presence is variable in natural infections. We also describe a variant of Jingmen tick virus, another Jingmenvirus, sequenced from a Ugandan red colobus monkey, thus expanding the host range of this segmented and likely multicomponent virus group. Collectively, this study provides evidence for the existence of multicomponent animal viruses and their potential relevance for animal and human health.

Standard finishing categories for high-throughput sequencing of viral genomes

Viral genome sequencing has become the cornerstone of almost all aspects of virology. In particular, high-throughput, next-generation viral genome sequencing has become an integral part of molecular epidemiological investigations into outbreaks of viral disease, such as the recent outbreaks of Middle Eastern respiratory syndrome, Ebola virus disease and Zika virus infection. Multiple institutes have acquired the expertise and necessary infrastructure to perform such investigations, as evidenced by the accumulation of thousands of novel viral sequences over progressively shorter time periods. The authors recently proposed a nomenclature comprised of five high-throughput sequencing standard categories to describe the quality of determined viral genome sequences. These five categories (standard draft, high quality, coding complete, complete and finished) cover all levels of viral genome finishing and can be applied to sequences determined by any technology platform or assembly technique.

Ebola Virus Epidemiology and Evolution in Nigeria

Containment limited the 2014 Nigerian Ebola virus (EBOV) disease outbreak to 20 reported cases and 8 fatalities. We present here clinical data and contact information for at least 19 case patients, and full-length EBOV genome sequences for 12 of the 20. The detailed contact data permits nearly complete reconstruction of the transmission tree for the outbreak. The EBOV genomic data are consistent with that tree. It confirms that there was a single source for the Nigerian infections, shows that the Nigerian EBOV lineage nests within a lineage previously seen in Liberia but is genetically distinct from it, and supports the conclusion that transmission from Nigeria to elsewhere did not occur.

Genomic Characterization of the Genus Nairovirus (Family Bunyaviridae)

Nairovirus, one of five bunyaviral genera, includes seven species. Genomic sequence information is limited for members of the Dera Ghazi Khan, Hughes, Qalyub, Sakhalin, and Thiafora nairovirus species. We used next-generation sequencing and historical virus-culture samples to determine 14 complete and nine coding-complete nairoviral genome sequences to further characterize these species. Previously unsequenced viruses include Abu Mina, Clo Mor, Great Saltee, Hughes, Raza, Sakhalin, Soldado, and Tillamook viruses. In addition, we present genomic sequence information on additional isolates of previously sequenced Avalon, Dugbe, Sapphire II, and Zirqa viruses. Finally, we identify Tunis virus, previously thought to be a phlebovirus, as an isolate of Abu Hammad virus. Phylogenetic analyses indicate the need for reassignment of Sapphire II virus to Dera Ghazi Khan nairovirus and reassignment of Hazara, Tofla, and Nairobi sheep disease viruses to novel species. We also propose new species for the Kasokero group (Kasokero, Leopards Hill, Yogue viruses), the Ketarah group (Gossas, Issyk-kul, Keterah/soft tick viruses) and the Burana group (Wēnzhōu tick virus, Huángpí tick virus 1, Tǎchéng tick virus 1). Our analyses emphasize the sister relationship of nairoviruses and arenaviruses, and indicate that several nairo-like viruses (Shāyáng spider virus 1, Xīnzhōu spider virus, Sānxiá water strider virus 1, South Bay virus, Wǔhàn millipede virus 2) require establishment of novel genera in a larger nairovirus-arenavirus supergroup.

Reduced evolutionary rate in reemerged Ebola virus transmission chains

On 29 June 2015, Liberia's respite from Ebola virus disease (EVD) was interrupted for the second time by a renewed outbreak ("flare-up") of seven confirmed cases. We demonstrate that, similar to the March 2015 flare-up associated with sexual transmission, this new flare-up was a reemergence of a Liberian transmission chain originating from a persistently infected source rather than a reintroduction from a reservoir or a neighboring country with active transmission. Although distinct, Ebola virus (EBOV) genomes from both flare-ups exhibit significantly low genetic divergence, indicating a reduced rate of EBOV evolution during persistent infection. Using this rate of change as a signature, we identified two additional EVD clusters that possibly arose from persistently infected sources. These findings highlight the risk of EVD flare-ups even after an outbreak is declared over.

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.

Pathosphere.org: pathogen detection and characterization through a web-based, open source informatics platform

BACKGROUND

The detection of pathogens in complex sample backgrounds has been revolutionized by wide access to next-generation sequencing (NGS) platforms. However, analytical methods to support NGS platforms are not as uniformly available. Pathosphere (found at Pathosphere.org) is a cloud - based open - sourced community tool that allows for communication, collaboration and sharing of NGS analytical tools and data amongst scientists working in academia, industry and government. The architecture allows for users to upload data and run available bioinformatics pipelines without the need for onsite processing hardware or technical support.

RESULTS

The pathogen detection capabilities hosted on Pathosphere were tested by analyzing pathogen-containing samples sequenced by NGS with both spiked human samples as well as human and zoonotic host backgrounds. Pathosphere analytical pipelines developed by Edgewood Chemical Biological Center (ECBC) identified spiked pathogens within a common sample analyzed by 454, Ion Torrent, and Illumina sequencing platforms. ECBC pipelines also correctly identified pathogens in human samples containing arenavirus in addition to animal samples containing flavivirus and coronavirus. These analytical methods were limited in the detection of sequences with limited homology to previous annotations within NCBI databases, such as parvovirus. Utilizing the pipeline-hosting adaptability of Pathosphere, the analytical suite was supplemented by analytical pipelines designed by the United States Army Medical Research Insititute of Infectious Diseases and Walter Reed Army Institute of Research (USAMRIID-WRAIR). These pipelines were implemented and detected parvovirus sequence in the sample that the ECBC iterative analysis previously failed to identify.

CONCLUSIONS

By accurately detecting pathogens in a variety of samples, this work demonstrates the utility of Pathosphere and provides a platform for utilizing, modifying and creating pipelines for a variety of NGS technologies developed to detect pathogens in complex sample backgrounds. These results serve as an exhibition for the existing pipelines and web-based interface of Pathosphere as well as the plug-in adaptability that allows for integration of newer NGS analytical software as it becomes available.

Ebola Virus Infections in Nonhuman Primates Are Temporally Influenced by Glycoprotein Poly-U Editing Site Populations in the Exposure Material

Recent experimentation with the variants of the Ebola virus that differ in the glycoprotein's poly-uridine site, which dictates the form of glycoprotein produced through a transcriptional stutter, has resulted in questions regarding the pathogenicity and lethality of the stocks used to develop products currently undergoing human clinical trials to combat the disease. In order to address these concerns and prevent the delay of these critical research programs, we designed an experiment that permitted us to intramuscularly challenge statistically significant numbers of naïve and vaccinated cynomolgus macaques with either a 7U or 8U variant of the Ebola virus, Kikwit isolate. In naïve animals, no difference in survivorship was observed; however, there was a significant delay in the disease course between the two groups. Significant differences were also observed in time-of-fever, serum chemistry, and hematology. In vaccinated animals, there was no statistical difference in survivorship between either challenge groups, with two succumbing in the 7U group compared to 1 in the 8U challenge group. In summary, survivorship was not affected, but the Ebola virus disease course in nonhuman primates is temporally influenced by glycoprotein poly-U editing site populations.

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

Characterization of the murine macrophage response to infection with virulent and avirulent Burkholderia species

BACKGROUND

Burkholderia pseudomallei (Bp) and Burkholderia mallei (Bm) are Gram-negative facultative intracellular pathogens, which are the causative agents of melioidosis and glanders, respectively. Depending on the route of exposure, aerosol or transcutaneous, infection by Bp or Bm can result in an extensive range of disease - from acute to chronic, relapsing illness to fatal septicemia. Both diseases are associated with difficult diagnosis and high fatality rates. About ninety five percent of patients succumb to untreated septicemic infections and the fatality rate is 50 % even when standard antibiotic treatments are administered.

RESULTS

The goal of this study is to profile murine macrophage-mediated phenotypic and molecular responses that are characteristic to a collection of Bp, Bm, Burkholderia thailandensis (Bt) and Burkholderia oklahomensis (Bo) strains obtained from humans, animals, environment and geographically diverse locations. Burkholderia spp. (N = 21) were able to invade and replicate in macrophages, albeit to varying degrees. All Bp (N = 9) and four Bm strains were able to induce actin polymerization on the bacterial surface following infection. Several Bp and Bm strains showed reduced ability to induce multinucleated giant cell (MNGC) formation, while Bo and Bp 776 were unable to induce this phenotype. Measurement of host cytokine responses revealed a statistically significant Bm mediated IL-6 and IL-10 production compared to Bp strains. Hierarchical clustering of transcriptional data from 84 mouse cytokines, chemokines and their corresponding receptors identified 29 host genes as indicators of differential responses between the Burkholderia spp. Further validation confirmed Bm mediated Il-1b, Il-10, Tnfrsf1b and Il-36a mRNA expressions were significantly higher when compared to Bp and Bt.

CONCLUSIONS

These results characterize the phenotypic and immunological differences in the host innate response to pathogenic and avirulent Burkholderia strains and provide insight into the phenotypic alterations and molecular targets underlying host-Burkholderia interactions.

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

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

Development of real-time PCR assays for the detection of Moraxella macacae associated with bloody nose syndrome in rhesus (Macaca mulatta) and cynomolgus (Macaca fascicularis) macaques

BACKGROUND

Moraxella macacae is a recently described bacterial pathogen that causes epistaxis or so-called bloody nose syndrome in captive macaques. The aim of this study was to develop specific molecular diagnostic assays for M. macacae and to determine their performance characteristics.

METHODS

We developed six real-time PCR assays on the Roche LightCycler. The accuracy, precision, selectivity, and limit of detection (LOD) were determined for each assay, in addition to further validation by testing nasal swabs from macaques presenting with epistaxis at the Tulane National Primate Research Center.

RESULTS

All assays exhibited 100% specificity and were highly sensitive with an LOD of 10 fg for chromosomal assays and 1 fg for the plasmid assay. Testing of nasal swabs from 10 symptomatic macaques confirmed the presence of M. macacae in these animals.

CONCLUSIONS

We developed several accurate, sensitive, and species-specific real-time PCR assays for the detection of M. macacae in captive macaques.

Ebola Virus Epidemiology, Transmission, and Evolution during Seven Months in Sierra Leone

The 2013–2015 Ebola virus disease (EVD) epidemic is caused by the Makona variant of Ebola virus (EBOV). Early in the epidemic, genome sequencing provided insights into virus evolution and transmission and offered important information for outbreak response. Here, we analyze sequences from 232 patients sampled over 7 months in Sierra Leone, along with 86 previously released genomes from earlier in the epidemic. We confirm sustained human-to-human transmission within Sierra Leone and find no evidence for import or export of EBOV across national borders after its initial introduction. Using high-depth replicate sequencing, we observe both host-to-host transmission and recurrent emergence of intrahost genetic variants. We trace the increasing impact of purifying selection in suppressing the accumulation of nonsynonymous mutations over time. Finally, we note changes in the mucin-like domain of EBOV glycoprotein that merit further investigation. These findings clarify the movement of EBOV within the region and describe viral evolution during prolonged human-to-human transmission.

•

In Sierra Leone, transmission has primarily been within-country, not between-country

•

Infectious doses are large enough for intrahost variants to transmit between hosts

•

A prolonged epidemic removes deleterious mutations from the viral population

•

There is preliminary evidence for human RNA editing effects on the Ebola genome

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.