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.

Power and sample size calculations for testing the ratio of reproductive values in phylogenetic samples.

The quality of the inferences we make from pathogen sequence data is determined by the number and composition of pathogen sequences that make up the sample used to drive that inference. However, there remains limited guidance on how to best structure and power studies when the end goal is phylogenetic inference. One question that we can attempt to answer with molecular data is whether some people are more likely to transmit a pathogen than others. Here we present an estimator to quantify differential transmission, as measured by the ratio of reproductive numbers between people with different characteristics, using transmission pairs linked by molecular data, along with a sample size calculation for this estimator. We also provide extensions to our method to correct for imperfect identification of transmission linked pairs, overdispersion in the transmission process, and group imbalance. We validate this method via simulation and provide tools to implement it in an R package, phylosamp.

Combining genomics and epidemiology to track mumps virus transmission in the United States.

Unusually large outbreaks of mumps across the United States in 2016 and 2017 raised questions about the extent of mumps circulation and the relationship between these and prior outbreaks. We paired epidemiological data from public health investigations with analysis of mumps virus whole genome sequences from 201 infected individuals, focusing on Massachusetts university communities. Our analysis suggests continuous, undetected circulation of mumps locally and nationally, including multiple independent introductions into Massachusetts and into individual communities. Despite the presence of these multiple mumps virus lineages, the genomic data show that one lineage has dominated in the US since at least 2006. Widespread transmission was surprising given high vaccination rates, but we found no genetic evidence that variants arising during this outbreak contributed to vaccine escape. Viral genomic data allowed us to reconstruct mumps transmission links not evident from epidemiological data or standard single-gene surveillance efforts and also revealed connections between apparently unrelated mumps outbreaks.

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.

Zika virus evolution and spread in the Americas.

Although the recent Zika virus (ZIKV) epidemic in the Americas and its link to birth defects have attracted a great deal of attention, much remains unknown about ZIKV disease epidemiology and ZIKV evolution, in part owing to a lack of genomic data. Here we address this gap in knowledge by using multiple sequencing approaches to generate 110 ZIKV genomes from clinical and mosquito samples from 10 countries and territories, greatly expanding the observed viral genetic diversity from this outbreak. We analysed the timing and patterns of introductions into distinct geographic regions; our phylogenetic evidence suggests rapid expansion of the outbreak in Brazil and multiple introductions of outbreak strains into Puerto Rico, Honduras, Colombia, other Caribbean islands, and the continental United States. We find that ZIKV circulated undetected in multiple regions for many months before the first locally transmitted cases were confirmed, highlighting the importance of surveillance of viral infections. We identify mutations with possible functional implications for ZIKV biology and pathogenesis, as well as those that might be relevant to the effectiveness of diagnostic tests.

Sample size calculation for phylogenetic case linkage.

Sample size calculations are an essential component of the design and evaluation of scientific studies. However, there is a lack of clear guidance for determining the sample size needed for phylogenetic studies, which are becoming an essential part of studying pathogen transmission. We introduce a statistical framework for determining the number of true infector-infectee transmission pairs identified by a phylogenetic study, given the size and population coverage of that study. We then show how characteristics of the criteria used to determine linkage and aspects of the study design can influence our ability to correctly identify transmission links, in sometimes counterintuitive ways. We test the overall approach using outbreak simulations and provide guidance for calculating the sensitivity and specificity of the linkage criteria, the key inputs to our approach. The framework is freely available as the R package phylosamp, and is broadly applicable to designing and evaluating a wide array of pathogen phylogenetic studies.

Repeated Coronavirus Disease 2019 Molecular Testing: Correlation of Severe Acute Respiratory Syndrome Coronavirus 2 Culture With Molecular Assays and Cycle Thresholds.

BACKGROUND: Repeated coronavirus disease 2019 (COVID-19) molecular testing can lead to positive test results after negative results and to multiple positive results over time. The association between positive test results and infectious virus is important to quantify. METHODS: A 2-month cohort of retrospective data and consecutively collected specimens from patients with COVID-19 or patients under investigation were used to understand the correlation between prolonged viral RNA positive test results, cycle threshold (Ct) values and growth of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in cell culture. Whole-genome sequencing was used to confirm virus genotype in patients with prolonged viral RNA detection. Droplet digital polymerase chain reaction was used to assess the rate of false-negative COVID-19 diagnostic test results. RESULTS: In 2 months, 29 686 specimens were tested and 2194 patients underwent repeated testing. Virus recovery in cell culture was noted in specimens with a mean Ct value of 18.8 (3.4) for SARS-CoV-2 target genes. Prolonged viral RNA shedding was associated with positive virus growth in culture in specimens collected up to 21 days after the first positive result but mostly in individuals symptomatic at the time of sample collection. Whole-genome sequencing provided evidence the same virus was carried over time. Positive test results following negative results had Ct values >29.5 and were not associated with virus culture. Droplet digital polymerase chain reaction results were positive in 5.6% of negative specimens collected from patients with confirmed or clinically suspected COVID-19. CONCLUSIONS: Low Ct values in SARS-CoV-2 diagnostic tests were associated with virus growth in cell culture. Symptomatic patients with prolonged viral RNA shedding can also be infectious.

Genomic Analysis of Lassa Virus during an Increase in Cases in Nigeria in 2018.

During 2018, an unusual increase in Lassa fever cases occurred in Nigeria, raising concern among national and international public health agencies. We analyzed 220 Lassa virus genomes from infected patients, including 129 from the 2017-2018 transmission season, to understand the viral populations underpinning the increase. A total of 14 initial genomes from 2018 samples were generated at Redeemer's University in Nigeria, and the findings were shared with the Nigerian Center for Disease Control in real time. We found that the increase in cases was not attributable to a particular Lassa virus strain or sustained by human-to-human transmission. Instead, the data were consistent with ongoing cross-species transmission from local rodent populations. Phylogenetic analysis also revealed extensive viral diversity that was structured according to geography, with major rivers appearing to act as barriers to migration of the rodent reservoir.

Evidence of Ebola Virus Replication and High Concentration in Semen of a Patient During Recovery.

In one patient over time, we found that concentration of Ebola virus RNA in semen during recovery is remarkably higher than blood at peak illness. Virus in semen is replication-competent with no change in viral genome over time. Presence of sense RNA suggests replication in cells present in semen.

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.

New Vibrio cholerae sequences from Eastern and Southern Africa alter our understanding of regional cholera transmission.

Despite ongoing containment and vaccination efforts, cholera remains prevalent in many countries in sub-Saharan Africa. Part of the difficulty in containing cholera comes from our lack of understanding of how it circulates throughout the region. To better characterize regional transmission, we generated and analyzed 118 Vibrio cholerae genomes collected between 2007-2019 from five different countries in Southern and Eastern Africa. We showed that V. cholerae sequencing can be successful from a variety of sample types and filled in spatial and temporal gaps in our understanding of circulating lineages, including providing some of the first sequences from the 2018-2019 outbreaks in Uganda, Kenya, Tanzania, Zambia, and Malawi. Our results present a complex picture of cholera transmission in the region, with multiple lineages found to be co-circulating within several countries. We also find evidence that previously identified sporadic cases may be from larger, undersampled outbreaks, highlighting the need for careful examination of sampling biases and underscoring the need for continued and expanded cholera surveillance across the African continent.

Sample size calculations for pathogen variant surveillance in the presence of biological and systematic biases.

Tracking the emergence and spread of pathogen variants is an important component of monitoring infectious disease outbreaks. To that end, accurately estimating the number and prevalence of pathogen variants in a population requires carefully designed surveillance programs. However, current approaches to calculating the number of pathogen samples needed for effective surveillance often do not account for the various processes that can bias which infections are detected and which samples are ultimately characterized as a specific variant. In this article, we introduce a framework that accounts for the logistical and epidemiological processes that may bias variant characterization, and we demonstrate how to use this framework (implemented in a publicly available tool) to calculate the number of sequences needed for surveillance. Our framework is designed to be easy to use while also flexible enough to be adapted to various pathogens and surveillance scenarios.

Single-cell profiling of lncRNA expression during Ebola virus infection in rhesus macaques.

Long non-coding RNAs (lncRNAs) are involved in numerous biological processes and are pivotal mediators of the immune response, yet little is known about their properties at the single-cell level. Here, we generate a multi-tissue bulk RNAseq dataset from Ebola virus (EBOV) infected and not-infected rhesus macaques and identified 3979 novel lncRNAs. To profile lncRNA expression dynamics in immune circulating single-cells during EBOV infection, we design a metric, Upsilon, to estimate cell-type specificity. Our analysis reveals that lncRNAs are expressed in fewer cells than protein-coding genes, but they are not expressed at lower levels nor are they more cell-type specific when expressed in the same number of cells. In addition, we observe that lncRNAs exhibit similar changes in expression patterns to those of protein-coding genes during EBOV infection, and are often co-expressed with known immune regulators. A few lncRNAs change expression specifically upon EBOV entry in the cell. This study sheds light on the differential features of lncRNAs and protein-coding genes and paves the way for future single-cell lncRNA studies.

Genomic and phenotypic analyses suggest moderate fitness differences among Zika virus lineages.

RNA viruses have short generation times and high mutation rates, allowing them to undergo rapid molecular evolution during epidemics. However, the extent of RNA virus phenotypic evolution within epidemics and the resulting effects on fitness and virulence remain mostly unknown. Here, we screened the 2015-2016 Zika epidemic in the Americas for lineage-specific fitness differences. We engineered a library of recombinant viruses representing twelve major Zika virus lineages and used them to measure replicative fitness within disease-relevant human primary cells and live mosquitoes. We found that two of these lineages conferred significant in vitro replicative fitness changes among human primary cells, but we did not find fitness changes in Aedes aegypti mosquitoes. Additionally, we found evidence for elevated levels of positive selection among five amino acid sites that define major Zika virus lineages. While our work suggests that Zika virus may have acquired several phenotypic changes during a short time scale, these changes were relatively moderate and do not appear to have enhanced transmission during the epidemic.

Prediction of Antibiotic Susceptibility Profiles of Vibrio cholerae Isolates From Whole Genome Illumina and Nanopore Sequencing Data: CholerAegon.

During the last decades, antimicrobial resistance (AMR) has become a global public health concern. Nowadays multi-drug resistance is commonly observed in strains of Vibrio cholerae, the etiological agent of cholera. In order to limit the spread of pathogenic drug-resistant bacteria and to maintain treatment options the analysis of clinical samples and their AMR profiles are essential. Particularly, in low-resource settings a timely analysis of AMR profiles is often impaired due to lengthy culturing procedures for antibiotic susceptibility testing or lack of laboratory capacity. In this study, we explore the applicability of whole genome sequencing for the prediction of AMR profiles of V. cholerae. We developed the pipeline CholerAegon for the in silico prediction of AMR profiles of 82 V. cholerae genomes assembled from long and short sequencing reads. By correlating the predicted profiles with results from phenotypic antibiotic susceptibility testing we show that the prediction can replace in vitro susceptibility testing for five of seven antibiotics. Because of the relatively low costs, possibility for real-time data analyses, and portability, the Oxford Nanopore Technologies MinION sequencing platform-especially in light of an upcoming less error-prone technology for the platform-appears to be well suited for pathogen genomic analyses such as the one described here. Together with CholerAegon, it can leverage pathogen genomics to improve disease surveillance and to control further spread of antimicrobial resistance.

Regional connectivity drove bidirectional transmission of SARS-CoV-2 in the Middle East during travel restrictions.

Regional connectivity and land travel have been identified as important drivers of SARS-CoV-2 transmission. However, the generalizability of this finding is understudied outside of well-sampled, highly connected regions. In this study, we investigated the relative contributions of regional and intercontinental connectivity to the source-sink dynamics of SARS-CoV-2 for Jordan and the Middle East. By integrating genomic, epidemiological and travel data we show that the source of introductions into Jordan was dynamic across 2020, shifting from intercontinental seeding in the early pandemic to more regional seeding for the travel restrictions period. We show that land travel, particularly freight transport, drove introduction risk during the travel restrictions period. High regional connectivity and land travel also drove Jordan's export risk. Our findings emphasize regional connectedness and land travel as drivers of transmission in the Middle East.

Regional sequencing collaboration reveals persistence of the T12 Vibrio cholerae O1 lineage in West Africa.

Background: Despite recent insights into cholera transmission patterns in Africa, regional and local dynamics in West Africa-where cholera outbreaks occur every few years-are still poorly understood. Coordinated genomic surveillance of Vibrio cholerae in the areas most affected may reveal transmission patterns important for cholera control. Methods: During a regional sequencing workshop in Nigeria, we sequenced 46 recent V. cholerae isolates from Cameroon, Niger, and Nigeria (37 from 2018 to 2019) to better understand the relationship between the V. cholerae bacterium circulating in these three countries. Results: From these isolates, we generated 44 whole Vibrio cholerae O1 sequences and analyzed them in the context of 1280 published V. cholerae O1 genomes. All sequences belonged to the T12 V. cholerae seventh pandemic lineage. Conclusions: Phylogenetic analysis of newly generated and previously published V. cholerae genomes suggested that the T12 lineage has been continuously transmitted within West Africa since it was first observed in the region in 2009, despite lack of reported cholera in the intervening years. The results from this regional sequencing effort provide a model for future regionally coordinated surveillance efforts. Funding: Funding for this project was provided by Bill and Melinda Gates Foundation OPP1195157.

Genomic diversity of SARS-CoV-2 during early introduction into the Baltimore-Washington metropolitan area.

The early COVID-19 pandemic was characterized by rapid global spread. In Maryland and Washington, DC, United States, more than 2500 cases were reported within 3 weeks of the first COVID-19 detection in March 2020. We aimed to use genomic sequencing to understand the initial spread of SARS-CoV-2 - the virus that causes COVID-19 - in the region. We analyzed 620 samples collected from the Johns Hopkins Health System during March 11-31, 2020, comprising 28.6% of the total cases in Maryland and Washington, DC. From these samples, we generated 114 complete viral genomes. Analysis of these genomes alongside a subsampling of over 1000 previously published sequences showed that the diversity in this region rivaled global SARS-CoV-2 genetic diversity at that time and that the sequences belong to all of the major globally circulating lineages, suggesting multiple introductions into the region. We also analyzed these regional SARS-CoV-2 genomes alongside detailed clinical metadata and found that clinically severe cases had viral genomes belonging to all major viral lineages. We conclude that efforts to control local spread of the virus were likely confounded by the number of introductions into the region early in the epidemic and the interconnectedness of the region as a whole.

Successive epidemic waves of cholera in South Sudan between 2014 and 2017: a descriptive epidemiological study.

BACKGROUND: Between 2014 and 2017, successive cholera epidemics occurred in South Sudan within the context of civil war, population displacement, flooding, and drought. We aim to describe the spatiotemporal and molecular features of the three distinct epidemic waves and explore the role of vaccination campaigns, precipitation, and population movement in shaping cholera spread in this complex setting. METHODS: In this descriptive epidemiological study, we analysed cholera linelist data to describe the spatiotemporal progression of the epidemics. We placed whole-genome sequence data from pandemic Vibrio cholerae collected throughout these epidemics into the global phylogenetic context. Using whole-genome sequence data in combination with other molecular attributes, we characterise the relatedness of strains circulating in each wave and the region. We investigated the association of rainfall and the instantaneous basic reproduction number using distributed lag non-linear models, compared county-level attack rates between those with early and late reactive vaccination campaigns, and explored the consistency of the spatial patterns of displacement and suspected cholera case reports. FINDINGS: The 2014 (6389 cases) and 2015 (1818 cases) cholera epidemics in South Sudan remained spatially limited whereas the 2016-17 epidemic (20 438 cases) spread among settlements along the Nile river. Initial cases of each epidemic were reported in or around Juba soon after the start of the rainy season, but we found no evidence that rainfall modulated transmission during each epidemic. All isolates analysed had similar genotypic and phenotypic characteristics, closely related to sequences from Uganda and Democratic Republic of the Congo. Large-scale population movements between counties of South Sudan with cholera outbreaks were consistent with the spatial distribution of cases. 21 of 26 vaccination campaigns occurred during or after the county-level epidemic peak. Counties vaccinated on or after the peak incidence week had 2·2 times (95% CI 2·1-2·3) higher attack rates than those where vaccination occurred before the peak. INTERPRETATION: Pandemic V cholerae of the same clonal origin was isolated throughout the study period despite interepidemic periods of no reported cases. Although the complex emergency in South Sudan probably shaped some of the observed spatial and temporal patterns of cases, the full scope of transmission determinants remains unclear. Timely and well targeted use of vaccines can reduce the burden of cholera; however, rapid vaccine deployment in complex emergencies remains challenging. FUNDING: The Bill & Melinda Gates Foundation.

Genomic Diversity of SARS-CoV-2 During Early Introduction into the United States National Capital Region.

BACKGROUND: The early COVID-19 pandemic has been characterized by rapid global spread. In the United States National Capital Region, over 2,000 cases were reported within three weeks of its first detection in March 2020. We aimed to use genomic sequencing to understand the initial spread of SARS-CoV-2, the virus that causes COVID-19, in the region. By correlating genetic information to disease phenotype, we also aimed to gain insight into any correlation between viral genotype and case severity or transmissibility. METHODS: We performed whole genome sequencing of clinical SARS-CoV-2 samples collected in March 2020 by the Johns Hopkins Health System. We analyzed these regional SARS-CoV-2 genomes alongside detailed clinical metadata and the global phylogeny to understand early establishment of the virus within the region. RESULTS: We analyzed 620 samples from the Johns Hopkins Health System collected between March 11-31, 2020, comprising 37.3% of the total cases in Maryland during this period. We selected 143 of these samples for sequencing, generating 114 complete viral genomes. These genomes belong to all five major Nextstrain-defined clades, suggesting multiple introductions into the region and underscoring the diversity of the regional epidemic. We also found that clinically severe cases had genomes belonging to all of these clades. CONCLUSIONS: We established a pipeline for SARS-CoV-2 sequencing within the Johns Hopkins Health system, which enabled us to capture the significant viral diversity present in the region as early as March 2020. Efforts to control local spread of the virus were likely confounded by the number of introductions into the region early in the epidemic and interconnectedness of the region as a whole.