DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus.

BACKGROUND: The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes. RESULTS: We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 10-100 RNA copies/µL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments. CONCLUSIONS: DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean.

Dengue is the most prevalent mosquito-borne viral disease in humans, and cases are continuing to rise globally. In particular, islands in the Caribbean have experienced more frequent outbreaks, and all four dengue virus (DENV) serotypes have been reported in the region, leading to hyperendemicity and increased rates of severe disease. However, there is significant variability regarding virus surveillance and reporting between islands, making it difficult to obtain an accurate understanding of the epidemiological patterns in the Caribbean. To investigate this, we used travel surveillance and genomic epidemiology to reconstruct outbreak dynamics, DENV serotype turnover, and patterns of spread within the region from 2009-2022. We uncovered two recent DENV-3 introductions from Asia, one of which resulted in a large outbreak in Cuba, which was previously under-reported. We also show that while outbreaks can be synchronized between islands, they are often caused by different serotypes. Our study highlights the importance of surveillance of infected travelers to provide a snapshot of local introductions and transmission in areas with limited local surveillance and suggests that the recent DENV-3 introductions may pose a major public health threat in the region.

Contribution of climate change to the spatial expansion of West Nile virus in Europe.

West Nile virus (WNV) is an emerging mosquito-borne pathogen in Europe where it represents a new public health threat. While climate change has been cited as a potential driver of its spatial expansion on the continent, a formal evaluation of this causal relationship is lacking. Here, we investigate the extent to which WNV spatial expansion in Europe can be attributed to climate change while accounting for other direct human influences such as land-use and human population changes. To this end, we trained ecological niche models to predict the risk of local WNV circulation leading to human cases to then unravel the isolated effect of climate change by comparing factual simulations to a counterfactual based on the same environmental changes but a counterfactual climate where long-term trends have been removed. Our findings demonstrate a notable increase in the area ecologically suitable for WNV circulation during the period 1901-2019, whereas this area remains largely unchanged in a no-climate-change counterfactual. We show that the drastic increase in the human population at risk of exposure is partly due to historical changes in population density, but that climate change has also been a critical driver behind the heightened risk of WNV circulation in Europe.

Genomic epidemiology of West Nile virus in Europe.

West Nile virus is one of the most widespread mosquito-borne zoonotic viruses, with unique transmission dynamics in various parts of the world. Genomic surveillance has provided important insights in the global patterns of West Nile virus emergence and spread. In Europe, multiple West Nile virus lineages have been isolated, with lineage 1a and 2 being the main lineages responsible for human infections. In contrast to North America, where a single introduction of lineage 1a resulted in the virus establishing itself in a new continent, at least 13 introductions of lineages 1a and 2 have occurred into Europe, which is likely a vast underestimation of the true number of introductions. Historically, lineage 1a was the main lineage circulating in Europe, but since the emergence of lineage 2 in the early 2000s, the latter has become the predominant lineage. This shift in West Nile virus lineage prevalence has been broadly linked to the expansion of the virus into northerly temperate regions, where autochthonous cases in animals and humans have been reported in Germany and The Netherlands. Here, we discuss how genomic analysis has increased our understanding of the epidemiology of West Nile virus in Europe, and we present a global Nextstrain build consisting of publicly available West Nile virus genomes (https://nextstrain.org/community/grubaughlab/WNV-Global). Our results elucidate recent insights in West Nile virus lineage dynamics in Europe, and discuss how expanded programs can fill current genomic surveillance gaps.

Introduction and Spread of Dengue Virus 3, Florida, USA, May 2022-April 2023.

During May 2022-April 2023, dengue virus serotype 3 was identified among 601 travel-associated and 61 locally acquired dengue cases in Florida, USA. All 203 sequenced genomes belonged to the same genotype III lineage and revealed potential transmission chains in which most locally acquired cases occurred shortly after introduction, with little sustained transmission.

Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean.

Dengue is the most prevalent mosquito-borne viral disease in humans, and cases are continuing to rise globally. In particular, islands in the Caribbean have experienced more frequent outbreaks, and all four dengue virus (DENV) serotypes have been reported in the region, leading to hyperendemicity and increased rates of severe disease. However, there is significant variability regarding virus surveillance and reporting between islands, making it difficult to obtain an accurate understanding of the epidemiological patterns in the Caribbean. To investigate this, we used travel surveillance and genomic epidemiology to reconstruct outbreak dynamics, DENV serotype turnover, and patterns of spread within the region from 2009-2022. We uncovered two recent DENV-3 introductions from Asia, one of which resulted in a large outbreak in Cuba, which was previously under-reported. We also show that while outbreaks can be synchronized between islands, they are often caused by different serotypes. Our study highlights the importance of surveillance of infected travelers to provide a snapshot of local introductions and transmission in areas with limited local surveillance and suggests that the recent DENV-3 introductions may pose a major public health threat in the region.

Survey of white-footed mice in Connecticut, USA reveals low SARS-CoV-2 seroprevalence and infection with divergent betacoronaviruses.

Diverse mammalian species display susceptibility to and infection with SARS-CoV-2. Potential SARS-CoV-2 spillback into rodents is understudied despite their host role for numerous zoonoses and human proximity. We assessed exposure and infection among white-footed mice (Peromyscus leucopus) in Connecticut, USA. We observed 1% (6/540) wild-type neutralizing antibody seroprevalence among 2020-2022 residential mice with no cross-neutralization of variants. We detected no SARS-CoV-2 infections via RT-qPCR, but identified non-SARS-CoV-2 betacoronavirus infections via pan-coronavirus PCR among 1% (5/468) of residential mice. Sequencing revealed two divergent betacoronaviruses, preliminarily named Peromyscus coronavirus-1 and -2. Both belong to the Betacoronavirus 1 species and are ~90% identical to the closest known relative, Porcine hemagglutinating encephalomyelitis virus. Low SARS-CoV-2 seroprevalence suggests white-footed mice may not be sufficiently susceptible or exposed to SARS-CoV-2 to present a long-term human health risk. However, the discovery of divergent, non-SARS-CoV-2 betacoronaviruses expands the diversity of known rodent coronaviruses and further investigation is required to understand their transmission extent.

DengueSeq: A pan-serotype whole genome amplicon sequencing protocol for dengue virus.

Background: The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes. Results: We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 101-102 RNA copies/µL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments. Conclusions: DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

The known unknowns of Powassan virus ecology.

Powassan virus (POWV; Family: Flaviviridae, Genus: Flavivirus) is the sole North American member of the tick-borne encephalitis sero-complex. While associated with high rates of morbidity and mortality, POWV has historically been of little public health concern due to low incidence rates. However, over the last 20 yr, incidence rates have increased highlighting the growing epidemiological threat. Currently, there are no vaccines or therapeutics with tick habitat reduction, acaricide application, and public awareness programs being our primary means of intervention. The effectiveness of these control strategies is dependent on having a sound understanding of the virus's ecology. In this Forum, we review what is currently known about POWV ecology, identify gaps in our knowledge, and discuss prevailing and alternative hypotheses about transmission dynamics, reservoir hosts, and spatial focality.

Dynamics of eastern equine encephalitis virus during the 2019 outbreak in the Northeast United States.

Eastern equine encephalitis virus (EEEV) causes a rare but severe disease in horses and humans and is maintained in an enzootic transmission cycle between songbirds and Culiseta melanura mosquitoes. In 2019, the largest EEEV outbreak in the United States for more than 50 years occurred, centered in the Northeast. To explore the dynamics of the outbreak, we sequenced 80 isolates of EEEV and combined them with existing genomic data. We found that, similar to previous years, cases were driven by multiple independent but short-lived virus introductions into the Northeast from Florida. Once in the Northeast, we found that Massachusetts was important for regional spread. We found no evidence of any changes in viral, human, or bird factors which would explain the increase in cases in 2019, although the ecology of EEEV is complex and further data is required to explore these in more detail. By using detailed mosquito surveillance data collected by Massachusetts and Connecticut, however, we found that the abundance of Cs. melanura was exceptionally high in 2019, as was the EEEV infection rate. We employed these mosquito data to build a negative binomial regression model and applied it to estimate early season risks of human or horse cases. We found that the month of first detection of EEEV in mosquito surveillance data and vector index (abundance multiplied by infection rate) were predictive of cases later in the season. We therefore highlight the importance of mosquito surveillance programs as an integral part of public health and disease control.

Phylogeographic reconstruction of the emergence and spread of Powassan virus in the northeastern United States.

Powassan virus is an emerging tick-borne virus of concern for public health, but very little is known about its transmission patterns and ecology. Here, we expanded the genomic dataset by sequencing 279 Powassan viruses isolated from Ixodes scapularis ticks from the northeastern United States. Our phylogeographic reconstructions revealed that Powassan virus lineage II was likely introduced or emerged from a relict population in the Northeast between 1940 and 1975. Sequences strongly clustered by sampling location, suggesting a highly focal geographical distribution. Our analyses further indicated that Powassan virus lineage II emerged in the northeastern United States mostly following a south-to-north pattern, with a weighted lineage dispersal velocity of ~3 km/y. Since the emergence in the Northeast, we found an overall increase in the effective population size of Powassan virus lineage II, but with growth stagnating during recent years. The cascading effect of population expansion of white-tailed deer and I. scapularis populations likely facilitated the emergence of Powassan virus in the northeastern United States.

Nonsystematic Reporting Biases of the SARS-CoV-2 Variant Mu Could Impact Our Understanding of the Epidemiological Dynamics of Emerging Variants.

Developing a timely and effective response to emerging SARS-CoV-2 variants of concern (VOCs) is of paramount public health importance. Global health surveillance does not rely on genomic data alone to identify concerning variants when they emerge. Instead, methods that utilize genomic data to estimate the epidemiological dynamics of emerging lineages have the potential to serve as an early warning system. However, these methods assume that genomic data are uniformly reported across circulating lineages. In this study, we analyze differences in reporting delays among SARS-CoV-2 VOCs as a plausible explanation for the timing of the global response to the former VOC Mu. Mu likely emerged in South America in mid-2020, where its circulation was largely confined. In this study, we demonstrate that Mu was designated as a VOC ∼1 year after it emerged and find that the reporting of genomic data for Mu differed significantly than that of other VOCs within countries, states, and individual laboratories. Our findings suggest that nonsystematic biases in the reporting of genomic data may have delayed the global response to Mu. Until they are resolved, the surveillance gaps that affected the global response to Mu could impede the rapid and accurate assessment of future emerging variants.

Toward a global virus genomic surveillance network.

The COVID-19 pandemic galvanized the field of virus genomic surveillance, demonstrating its utility for public health. Now, we must harness the momentum that led to increased infrastructure, training, and political will to build a sustainable global genomic surveillance network for other epidemic and endemic viruses. We suggest a generalizable modular sequencing framework wherein users can easily switch between virus targets to maximize cost-effectiveness and maintain readiness for new threats. We also highlight challenges associated with genomic surveillance and when global inequalities persist. We propose solutions to mitigate some of these issues, including training and multilateral partnerships. Exploring alternatives to clinical sequencing can also reduce the cost of surveillance programs. Finally, we discuss how establishing genomic surveillance would aid control programs and potentially provide a warning system for outbreaks, using a global respiratory virus (RSV), an arbovirus (dengue virus), and a regional zoonotic virus (Lassa virus) as examples.

Dynamics of Eastern equine encephalitis virus during the 2019 outbreak in the Northeast United States.

Eastern equine encephalitis virus (EEEV) causes a rare but severe disease in horses and humans, and is maintained in an enzootic transmission cycle between songbirds and Culiseta melanura mosquitoes. In 2019, the largest EEEV outbreak in the United States for more than 50 years occurred, centered in the Northeast. To explore the dynamics of the outbreak, we sequenced 80 isolates of EEEV and combined them with existing genomic data. We found that, like previous years, cases were driven by frequent short-lived virus introductions into the Northeast from Florida. Once in the Northeast, we found that Massachusetts was important for regional spread. We found no evidence of any changes in viral, human, or bird factors which would explain the increase in cases in 2019. By using detailed mosquito surveillance data collected by Massachusetts and Connecticut, however, we found that the abundance of Cs. melanura was exceptionally high in 2019, as was the EEEV infection rate. We employed these mosquito data to build a negative binomial regression model and applied it to estimate early season risks of human or horse cases. We found that the month of first detection of EEEV in mosquito surveillance data and vector index (abundance multiplied by infection rate) were predictive of cases later in the season. We therefore highlight the importance of mosquito surveillance programs as an integral part of public health and disease control.

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.

Age-dependent impairment in antibody responses elicited by a homologous CoronaVac booster dose.

The emergence of the SARS-CoV-2 Omicron sublineages resulted in increased transmission rates and reduced protection from vaccines. To counteract these effects, multiple booster strategies were used in different countries, although data comparing their efficiency in improving protective immunity remain sparse, especially among vulnerable populations, including older adults. The inactivated CoronaVac vaccine was among the most widely distributed vaccine worldwide and was essential in the early control of SARS-CoV-2-related hospitalizations and deaths. However, it is not well understood whether homologous versus heterologous booster doses in those fully vaccinated with CoronaVac induce distinct humoral responses or whether these responses vary across age groups. We analyzed plasma antibody responses from CoronaVac-vaccinated younger or older individuals who received a homologous CoronaVac or heterologous BNT162b2 or ChAdOx1 booster vaccine. All three evaluated boosters resulted in increased virus-specific IgG titers 28 days after the booster dose. However, we found that both IgG titers against SARS-CoV-2 Spike or RBD and neutralization titers against Omicron sublineages were substantially reduced in participants who received homologous CoronaVac compared with the heterologous BNT162b2 or ChAdOx1 booster. This effect was specifically prominent in recipients >50 years of age. In this group, the CoronaVac booster induced low virus-specific IgG titers and failed to elevate neutralization titers against any Omicron sublineage. Our results point to the notable inefficiency of CoronaVac immunization and boosting in mounting protective antiviral humoral immunity, particularly among older adults, during the Omicron wave. These observations also point to benefits of heterologous regimens in high-risk populations fully vaccinated with CoronaVac.

Powassan Virus Lineage I in Field-Collected Dermacentor variabilis Ticks, New York, USA.

Powassan virus is a tickborne flavivirus that can cause lethal or debilitating neurologic illness. It is canonically transmitted by Ixodes spp. ticks but might spill over to sympatric Dermacentor species. We detected Powassan virus lineage I from a pool of field-collected D. variabilis ticks in New York, USA.

Identification and characterization of novel lineage 1 Powassan virus strains in New York State.

Powassan virus (POWV, family Flaviviridae) is a reemerging tick-borne virus endemic in North America and Russia. In 1997, a POWV-like agent was isolated from Ixodes scapularis in New England and determined to be genetically distinct from the original POWV isolate. This revealed the existence of two lineages: lineage 1, prototype Powassan virus (POWV-1) and lineage 2, deer tick virus (DTV). POWV-1 is thought to be primarily maintained in a cycle between I. cookei and woodchucks and I. marxi and squirrels, while DTV is primarily maintained in a cycle between I. scapularis and small mammal hosts. Recent tick, mammalian, and human isolates from New York State (NYS) have been identified as DTV, but for the first time in 45 years, we detected four POWV-1 isolates, including the first reported isolation of POWV-1 from I. scapularis. We aimed to investigate genotypic and phenotypic characteristics of recent NYS isolates through sequence analysis and evaluation of replication kinetics in vitro and in vivo. Our sequencing revealed genetic divergence between NYS POWV-1 isolates, with two distinct foci. We found that POWV-1 isolates displayed variable replication kinetics in nymphal ticks but not in cell culture. POWV-1 isolated from I. scapularis displayed increased fitness in experimentally infected I. scapularis as compared to historic and recent POWV-1 isolates from I. cookei. These data suggest the emergence of divergent POWV-1 strains in alternate tick hosts and maintenance of genetically and phenotypically discrete POWV-1 foci.

Nasal host response-based screening for undiagnosed respiratory viruses: a pathogen surveillance and detection study.

BACKGROUND: Symptomatic patients who test negative for common viruses are an important possible source of unrecognised or emerging pathogens, but metagenomic sequencing of all samples is inefficient because of the low likelihood of finding a pathogen in any given sample. We aimed to determine whether nasopharyngeal CXCL10 screening could be used as a strategy to enrich for samples containing undiagnosed viruses. METHODS: In this pathogen surveillance and detection study, we measured CXCL10 concentrations from nasopharyngeal swabs from patients in the Yale New Haven health-care system, which had been tested at the Yale New Haven Hospital Clinical Virology Laboratory (New Haven, CT, USA). Patients who tested negative for a panel of respiratory viruses using multiplex PCR during Jan 23-29, 2017, or March 3-14, 2020, were included. We performed host and pathogen RNA sequencing (RNA-Seq) and analysis for viral reads on samples with CXCL10 higher than 1 ng/mL or CXCL10 testing and quantitative RT-PCR (RT-qPCR) for SARS-CoV-2. We used RNA-Seq and cytokine profiling to compare the host response to infection in samples that were virus positive (rhinovirus, seasonal coronavirus CoV-NL63, or SARS-CoV-2) and virus negative (controls). FINDINGS: During Jan 23-29, 2017, 359 samples were tested for ten viruses on the multiplex PCR respiratory virus panel (RVP). 251 (70%) were RVP negative. 60 (24%) of 251 samples had CXCL10 higher than 150 pg/mL and were identified for further analysis. 28 (47%) of 60 CXCL10-high samples were positive for seasonal coronaviruses. 223 (89%) of 251 samples were PCR negative for 15 viruses and, of these, CXCL10-based screening identified 32 (13%) samples for further analysis. Of these 32 samples, eight (25%) with CXCL10 concentrations higher than 1 ng/mL and sufficient RNA were selected for RNA-Seq. Microbial RNA analysis showed the presence of influenza C virus in one sample and revealed RNA reads from bacterial pathobionts in four (50%) of eight samples. Between March 3 and March 14, 2020, 375 (59%) of 641 samples tested negative for 15 viruses on the RVP. 32 (9%) of 375 samples had CXCL10 concentrations ranging from 100 pg/mL to 1000 pg/mL and four of those were positive for SARS-CoV-2. CXCL10 elevation was statistically significant, and a distinguishing feature was found in 28 (8%) of 375 SARS-CoV-2-negative samples versus all four SARS-CoV-2-positive samples (p=4·4 × 10-5). Transcriptomic signatures showed an interferon response in virus-positive samples and an additional neutrophil-high hyperinflammatory signature in samples with high amounts of bacterial pathobionts. The CXCL10 cutoff for detecting a virus was 166·5 pg/mL for optimal sensitivity and 1091·0 pg/mL for specificity using a clinic-ready automated microfluidics-based immunoassay. INTERPRETATION: These results confirm CXCL10 as a robust nasopharyngeal biomarker of viral respiratory infection and support host response-based screening followed by metagenomic sequencing of CXCL10-high samples as a practical approach to incorporate clinical samples into pathogen discovery and surveillance efforts. FUNDING: National Institutes of Health, the Hartwell Foundation, the Gruber Foundation, Fast Grants for COVID-19 research from the Mercatus Center, and the Huffman Family Donor Advised Fund.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques.

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.