Genomic Epidemiology as a Public Health Tool to Combat Mosquito-Borne Virus Outbreaks

National investment case development for pathogen genomics

Sustaining and expanding genomic surveillance capacity requires broader investment in genomics that target both novel and pandemic pathogens. Currently, there is no standardized methodology to evaluate the cost and benefit of a multi-pathogen surveillance system. We propose a framework for pathogen genomic surveillance that links public health and systems considerations to a stepwise approach.

Advancing West Nile virus monitoring through whole genome sequencing: Insights from a One Health genomic surveillance study in Romagna (Italy)

In the last 6 years, Italy accounted for 36 % of the total autochthonous European West Nile virus (WNV) cases reported to ECDC. Since 2001, the country put in place a multi-species national surveillance plan. The plan was enhanced in 2020 by adopting a fully integrated "One Health" approach, including human, wild bird, equine, and mosquito surveillance for the early detection of WNV. In this context, the systematic acquisition of whole viral genetic information from human patients and animals is fundamental to obtain an in-depth knowledge on the patterns of virus evolution and transmission and to gain insights on the role virus genetics in morbidity and mortality, The purpose of this pilot study was thus to design a One-Health surveillance framework based on the genomic surveillance of WNV circulating at the vector-human-animal interface, in the endemic territory of Romagna (North-Eastern Italy) during the 2023 transmission season. Whole genome sequencing (WGS) analyses confirmed the circulation of WNV lineage 2 showing high nucleotide and amino acid identity of 99.82 % and 99.92 % respectively among viral sequences from human patients, vectors and birds. All the sequences clustered with other Italian strains in the Central and Southern European clade with robust bootstrap support and BLASTn identity exceeding 99.7 %. The highest nucleotide identity was observed with sequences from Emilia-Romagna and Veneto regions (Italy), confirming a local virus circulation and overwintering of WNV lineage 2 with a confined virus spread and no (or limited) external introduction of viral strains. Our results, support the adoption of a One Health approach to WNV surveillance, based on WGS and integrating the clinical diagnosis, epidemiology, and genomic characterisation, to create a suitable operational process for the characterisation of autochthonous and imported Arboviruses circulating in Romagna to effectively integrate the already established surveillance plan.

Long-term co-circulation of multiple arboviruses in southeast Australia revealed by xeno-monitoring and viral whole-genome sequencing

Arbovirus surveillance of wild-caught mosquitoes is an affordable and sensitive means of monitoring virus transmission dynamics at various spatial-temporal scales, and emergence and re-emergence during epidemic and interepidemic periods. A variety of molecular diagnostics for arbovirus screening of mosquitoes (known as xeno-monitoring) are available, but most provide limited information about virus diversity. Polymerase chain reaction (PCR)-based screening coupled with RNA sequencing is an increasingly affordable and sensitive pipeline for integrating complete viral genome sequencing into surveillance programs. This enables large-scale, high-throughput arbovirus screening from diverse samples. We collected mosquitoes in CO2-baited light traps from five urban parks in Brisbane from March 2021 to May 2022. Mosquito pools of ≤200 specimens were screened for alphaviruses and flaviviruses using virus genus-specific primers and reverse transcription quantitative PCR (qRT-PCR). A subset of virus-positive samples was then processed using a mosquito-specific ribosomal RNA depletion method and then sequenced on the Illumina NextSeq. Overall, 54,670 mosquitoes representing 26 species were screened in 382 pools. Thirty detections of arboviruses were made in 28 pools. Twenty of these positive pools were further characterized using RNA sequencing generating 18 full-length genomes. These full-length sequences belonged to four medically relevant arboviruses: Barmah Forest, Ross River, Sindbis-like, and Stratford viruses. Phylogenetic and evolutionary analyses revealed the evolutionary progression of arbovirus lineages over the last 100 years, demonstrating that different epidemiological, immunological, and evolutionary processes may actively shape the evolution of Australian arboviruses. These results underscore the need for more genomic surveillance data to explore the complex evolutionary pressures acting on arboviruses. Overall, our findings highlight the effectiveness of our methodology, which can be applied broadly to enhance arbovirus surveillance in various ecological contexts and improve understanding of transmission dynamics.

Pathogen genomic surveillance status among lower resource settings in Asia

Asia remains vulnerable to new and emerging infectious diseases. Understanding how to improve next generation sequencing (NGS) use in pathogen surveillance is an urgent priority for regional health security. Here we developed a pathogen genomic surveillance assessment framework to assess capacity in low-resource settings in South and Southeast Asia. Data collected between June 2022 and March 2023 from 42 institutions in 13 countries showed pathogen genomics capacity exists, but use is limited and under-resourced. All countries had NGS capacity and seven countries had strategic plans integrating pathogen genomics into wider surveillance efforts. Several pathogens were prioritized for human surveillance, but NGS application to environmental and human-animal interface surveillance was limited. Barriers to NGS implementation include reliance on external funding, supply chain challenges, trained personnel shortages and limited quality assurance mechanisms. Coordinated efforts are required to support national planning, address capacity gaps, enhance quality assurance and facilitate data sharing for decision making.

Molecular epidemiology of dengue in Malaysia: 2015-2021

Dengue has been one of the major public health problems in Malaysia for decades. Over 600,000 dengue cases and 1,200 associated fatalities have been reported in Malaysia from 2015 to 2021, which was 100% increase from the cumulative total of dengue cases reported during the preceding 07-year period from 2008 to 2014. However, studies that describe the molecular epidemiology of dengue in Malaysia in recent years are limited. In the present study, we describe the genetic composition and dispersal patterns of Dengue virus (DENV) by using 4,004 complete envelope gene sequences of all four serotypes (DENV-1 = 1,567, DENV-2 = 1,417, DENV-3 = 762 and DENV-4 = 258) collected across Malaysia from 2015 to 2021. The findings revealed that DENV populations in Malaysia were highly diverse, and the overall heterogeneity was maintained through repetitive turnover of genotypes. Phylogeography analyses suggested that DENV dispersal occurred through an extensive network, mainly among countries in South and East Asia and Malaysian states, as well as among different states, especially within Peninsular Malaysia. The results further suggested Selangor and Johor as major hubs of DENV emergence and spread in Malaysia.

Leveraging insect-specific viruses to elucidate mosquito population structure and dynamics

Several aspects of mosquito ecology that are important for vectored disease transmission and control have been difficult to measure at epidemiologically important scales in the field. In particular, the ability to describe mosquito population structure and movement rates has been hindered by difficulty in quantifying fine-scale genetic variation among populations. The mosquito virome represents a possible avenue for quantifying population structure and movement rates across multiple spatial scales. Mosquito viromes contain a diversity of viruses, including several insect-specific viruses (ISVs) and "core" viruses that have high prevalence across populations. To date, virome studies have focused on viral discovery and have only recently begun examining viral ecology. While nonpathogenic ISVs may be of little public health relevance themselves, they provide a possible route for quantifying mosquito population structure and dynamics. For example, vertically transmitted viruses could behave as a rapidly evolving extension of the host's genome. It should be possible to apply established analytical methods to appropriate viral phylogenies and incidence data to generate novel approaches for estimating mosquito population structure and dispersal over epidemiologically relevant timescales. By studying the virome through the lens of spatial and genomic epidemiology, it may be possible to investigate otherwise cryptic aspects of mosquito ecology. A better understanding of mosquito population structure and dynamics are key for understanding mosquito-borne disease ecology and methods based on ISVs could provide a powerful tool for informing mosquito control programs.

Detection of Zika Virus in Aedes aegypti and Aedes albopictus Mosquitoes Collected in Urban Forest Fragments in the Brazilian Amazon

Zika virus (ZIKV) is an RNA flavivirus (Flaviviridae family) endemic in tropical and subtropical regions that is transmitted to humans by Aedes (Stegomyia) species mosquitoes. The two main urban vectors of ZIKV are Aedes aegypti and Aedes albopictus, which can be found throughout Brazil. This study investigated ZIKV infection in mosquito species sampled from urban forest fragments in Manaus (Brazilian Amazon). A total of 905 non-engorged female Ae. aegypti (22 specimens) and Ae. albopictus (883 specimens) were collected using BG-Sentinel traps, entomological hand nets, and Prokopack aspirators during the rainy and dry seasons between 2018 and 2021. All pools were macerated and used to inoculate C6/36 culture cells. Overall, 3/20 (15%) Ae. aegypti and 5/241 (2%) Ae. albopictus pools screened using RT-qPCR were positive for ZIKV. No supernatants from Ae. aegypti were positive for ZIKV (0%), and 15 out of 241 (6.2%) Ae. albopictus pools were positive. In this study, we provide the first-ever evidence of Ae. albopictus naturally infected with ZIKV in the Amazon region.

Enhanced Arbovirus Surveillance with High-Throughput Metatranscriptomic Processing of Field-Collected Mosquitoes

Surveillance programs are essential for the prevention and control of mosquito-borne arboviruses that cause serious human and animal diseases. Viral metatranscriptomic sequencing can enhance surveillance by enabling untargeted, high-throughput arbovirus detection. We used metatranscriptomic sequencing to screen field-collected mosquitoes for arboviruses to better understand how metatranscriptomics can be utilised in routine surveillance. Following a significant flood event in 2016, more than 56,000 mosquitoes were collected over seven weeks from field traps set up in Victoria, Australia. The traps were split into samples of 1000 mosquitoes or less and sequenced on the Illumina HiSeq. Five arboviruses relevant to public health (Ross River virus, Sindbis virus, Trubanaman virus, Umatilla virus, and Wongorr virus) were detected a total of 33 times in the metatranscriptomic data, with 94% confirmed using reverse transcription quantitative PCR (RT-qPCR). Analysis of metatranscriptomic cytochrome oxidase I (COI) sequences enabled the detection of 12 mosquito and two biting midge species. Screening of the same traps by an established public health arbovirus surveillance program corroborated the metatranscriptomic arbovirus and mosquito species detections. Assembly of genome sequences from the metatranscriptomic data also led to the detection of 51 insect-specific viruses, both known and previously undescribed, and allowed phylogenetic comparison to past strains. We have demonstrated how metatranscriptomics can enhance surveillance by enabling untargeted arbovirus detection, providing genomic epidemiological data, and simultaneously identifying vector species from large, unsorted mosquito traps.

A unified global genotyping framework of dengue virus serotype-1 for a stratified coordinated surveillance strategy of dengue epidemics

BACKGROUND

Dengue is the fastest spreading arboviral disease, posing great challenges on global public health. A reproduceable and comparable global genotyping framework for contextualizing spatiotemporal epidemiological data of dengue virus (DENV) is essential for research studies and collaborative surveillance.

METHODS

Targeting DENV-1 spreading prominently in recent decades, by reconciling all qualified complete E gene sequences of 5003 DENV-1 strains with epidemiological information from 78 epidemic countries/areas ranging from 1944 to 2018, we established and characterized a unified global high-resolution genotyping framework using phylogenetics, population genetics, phylogeography, and phylodynamics.

RESULTS

The defined framework was discriminated with three hierarchical layers of genotype, subgenotype and clade with respective mean pairwise distances 2-6%, 0.8-2%, and ≤ 0.8%. The global epidemic patterns of DENV-1 showed strong geographic constraints representing stratified spatial-genetic epidemic pairs of Continent-Genotype, Region-Subgenotype and Nation-Clade, thereby identifying 12 epidemic regions which prospectively facilitates the region-based coordination. The increasing cross-transmission trends were also demonstrated. The traditional endemic countries such as Thailand, Vietnam and Indonesia displayed as persisting dominant source centers, while the emerging epidemic countries such as China, Australia, and the USA, where dengue outbreaks were frequently triggered by importation, showed a growing trend of DENV-1 diffusion. The probably hidden epidemics were found especially in Africa and India. Then, our framework can be utilized in an accurate stratified coordinated surveillance based on the defined viral population compositions. Thereby it is prospectively valuable for further hampering the ongoing transition process of epidemic to endemic, addressing the issue of inadequate monitoring, and warning us to be concerned about the cross-national, cross-regional, and cross-continental diffusions of dengue, which can potentially trigger large epidemics.

CONCLUSIONS

The framework and its utilization in quantitatively assessing DENV-1 epidemics has laid a foundation and re-unveiled the urgency for establishing a stratified coordinated surveillance platform for blocking global spreading of dengue. This framework is also expected to bridge classical DENV-1 genotyping with genomic epidemiology and risk modeling. We will promote it to the public and update it periodically.

The human "contaminome": bacterial, viral, and computational contamination in whole genome sequences from 1000 families

The unmapped readspace of whole genome sequencing data tends to be large but is often ignored. We posit that it contains valuable signals of both human infection and contamination. Using unmapped and poorly aligned reads from whole genome sequences (WGS) of over 1000 families and nearly 5000 individuals, we present insights into common viral, bacterial, and computational contamination that plague whole genome sequencing studies. We present several notable results: (1) In addition to known contaminants such as Epstein-Barr virus and phiX, sequences from whole blood and lymphocyte cell lines contain many other contaminants, likely originating from storage, prep, and sequencing pipelines. (2) Sequencing plate and biological sample source of a sample strongly influence contamination profile. And, (3) Y-chromosome fragments not on the human reference genome commonly mismap to bacterial reference genomes. Both experiment-derived and computational contamination is prominent in next-generation sequencing data. Such contamination can compromise results from WGS as well as metagenomics studies, and standard protocols for identifying and removing contamination should be developed to ensure the fidelity of sequencing-based studies.

Effect of low-passage number on dengue consensus genomes and intra-host variant frequencies

Intra-host single nucleotide variants (iSNVs) have been increasingly used in genomic epidemiology to increase phylogenetic resolution and reconstruct fine-scale outbreak dynamics. These analyses are preferably done on sequence data from direct clinical samples, but in many cases due to low viral loads, there might not be enough genetic material for deep sequencing and iSNV determination. Isolation of the virus from clinical samples with low-passage number increases viral load, but few studies have investigated how dengue virus (DENV) culture isolation from a clinical sample impacts the consensus sequence and the intra-host virus population frequencies. In this study, we investigate consensus and iSNV frequency differences between DENV sequenced directly from clinical samples and their corresponding low-passage isolates. Twenty five DENV1 and DENV2 positive sera and their corresponding viral isolates (T. splendens inoculation and C6/36 passage) were obtained from a prospective cohort study in the Philippines. These were sequenced on MiSeq with minimum nucleotide depth of coverage of 500×, and iSNVs were detected using LoFreq. For both DENV1 and DENV2, we found a maximum of one consensus nucleotide difference between clinical sample and isolate. Interestingly, we found that iSNVs with frequencies ≥5 % were often preserved between the samples, and that the number of iSNV positions, and sample diversity, at this frequency cutoff did not differ significantly between the sample pairs (clinical sample and isolate) in either DENV1 or DENV2 data. Our results show that low-passage DENV isolate consensus genomes are largely representative of their direct sample parental viruses, and that low-passage isolates often mirror high frequency within-host variants from direct samples.

Revealing fine-scale spatiotemporal differences in SARS-CoV-2 introduction and spread

Evidence-based public health approaches that minimize the introduction and spread of new SARS-CoV-2 transmission clusters are urgently needed in the United States and other countries struggling with expanding epidemics. Here we analyze 247 full-genome SARS-CoV-2 sequences from two nearby communities in Wisconsin, USA, and find surprisingly distinct patterns of viral spread. Dane County had the 12th known introduction of SARS-CoV-2 in the United States, but this did not lead to descendant community spread. Instead, the Dane County outbreak was seeded by multiple later introductions, followed by limited community spread. In contrast, relatively few introductions in Milwaukee County led to extensive community spread. We present evidence for reduced viral spread in both counties following the statewide "Safer at Home" order, which went into effect 25 March 2020. Our results suggest patterns of SARS-CoV-2 transmission may vary substantially even in nearby communities. Understanding these local patterns will enable better targeting of public health interventions.

SARS-CoV-2 sequencing: The technological initiative to strengthen early warning systems for public health emergencies in Latin America and the Caribbean

The COVID-19 pandemic caused by SARS-CoV-2 is a public health problem on a scale unprecedented in the last 100 years, as has been the response focused on the rapid genomic characterization of SARS-CoV-2 in virtually all regions of the planet. This pandemic emerged during the era of genomic epidemiology, a science fueled by continued advances in next-generation sequencing. Since its recent appearance, genomic epidemiology included the precise identification of new lineages or species of pathogens and the reconstruction of their genetic variability in real time, evidenced in past outbreaks of influenza H1N1, MERS, and SARS. However, the global and uncontrolled scale of this pandemic created a scenario where genomic epidemiology was put into practice en masse, from the rapid identification of SARS-CoV-2 to the registration of new lineages and their active surveillance throughout the world. Prior to the COVID-19 pandemic, the availability of genomic data on circulating pathogens in several Latin America and the Caribbean countries was scarce or nil. With the arrival of SARS-CoV-2, this scenario changed significantly, although the amount of available information remains scarce and, in countries such as Colombia, Brazil, Argentina, and Chile, the genomic information of SARS-CoV-2 was obtained mainly by research groups in genomic epidemiology rather than the product of a public health surveillance policy or program. This indicates the need to establish public health policies aimed at implementing genomic epidemiology as a tool to strengthen surveillance and early warning systems against threats to public health in the region.

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

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

Highlights in Medical Entomology, 2019: Familiar Foes and New Frontiers

The 2019 Entomological Society of America annual meeting was held in St. Louis, Missouri, just blocks away from the iconic Gateway Arch. Representing a 'gateway to the West', this inspired the theme of the Highlights in Medical Entomology to reflect on the accomplishments of the past year as we move into a 'new frontier' of vector biology research. Papers were selected broadly across arthropods that influence public health, focusing on topics ranging from West Nile virus transmission, ticks and tick-borne disease, to advances in genetics and 'big data' studies. This included current perspectives on West Nile virus ecology and epidemiology, which has now been endemic in the United States for 20 yr. Additional topics such as the advantages of citizen science and the importance of scientific communication were also discussed. Together, these papers demonstrate the achievements of the vector community while emphasizing the challenges that we collectively face to reduce the burden of vector-borne disease.

Distinct patterns of SARS-CoV-2 transmission in two nearby communities in Wisconsin, USA

Evidence-based public health approaches that minimize the introduction and spread of new SARS-CoV-2 transmission clusters are urgently needed in the United States and other countries struggling with expanding epidemics. Here we analyze 247 full-genome SARS-CoV-2 sequences from two nearby communities in Wisconsin, USA, and find surprisingly distinct patterns of viral spread. Dane County had the 12th known introduction of SARS-CoV-2 in the United States, but this did not lead to descendant community spread. Instead, the Dane County outbreak was seeded by multiple later introductions, followed by limited community spread. In contrast, relatively few introductions in Milwaukee County led to extensive community spread. We present evidence for reduced viral spread in both counties, and limited viral transmission between counties, following the statewide Safer-at-Home public health order, which went into effect 25 March 2020. Our results suggest that early containment efforts suppressed the spread of SARS-CoV-2 within Wisconsin.

Phylogenetic interpretation during outbreaks requires caution

How viruses are related, and how they have evolved and spread over time, can be investigated using phylogenetics. Here, we set out how genomic analyses should be used during an epidemic and propose that phylogenetic insights from the early stages of an outbreak should heed all of the available epidemiological information.