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

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.

Human movement and environmental barriers shape the emergence of dengue

Understanding how emerging infectious diseases spread within and between countries is essential to contain future pandemics. Spread to new areas requires connectivity between one or more sources and a suitable local environment, but how these two factors interact at different stages of disease emergence remains largely unknown. Further, no analytical framework exists to examine their roles. Here we develop a dynamic modelling approach for infectious diseases that explicitly models both connectivity via human movement and environmental suitability interactions. We apply it to better understand recently observed (1995-2019) patterns as well as predict past unobserved (1983-2000) and future (2020-2039) spread of dengue in Mexico and Brazil. We find that these models can accurately reconstruct long-term spread pathways, determine historical origins, and identify specific routes of invasion. We find early dengue invasion is more heavily influenced by environmental factors, resulting in patchy non-contiguous spread, while short and long-distance connectivity becomes more important in later stages. Our results have immediate practical applications for forecasting and containing the spread of dengue and emergence of new serotypes. Given current and future trends in human mobility, climate, and zoonotic spillover, understanding the interplay between connectivity and environmental suitability will be increasingly necessary to contain emerging and re-emerging pathogens.

Hybrid-Capture Target Enrichment in Human Pathogens: Identification, Evolution, Biosurveillance, and Genomic Epidemiology

High-throughput sequencing (HTS) has revolutionised the field of pathogen genomics, enabling the direct recovery of pathogen genomes from clinical and environmental samples. However, pathogen nucleic acids are often overwhelmed by those of the host, requiring deep metagenomic sequencing to recover sufficient sequences for downstream analyses (e.g., identification and genome characterisation). To circumvent this, hybrid-capture target enrichment (HC) is able to enrich pathogen nucleic acids across multiple scales of divergences and taxa, depending on the panel used. In this review, we outline the applications of HC in human pathogens-bacteria, fungi, parasites and viruses-including identification, genomic epidemiology, antimicrobial resistance genotyping, and evolution. Importantly, we explored the applicability of HC to clinical metagenomics, which ultimately requires more work before it is a reliable and accurate tool for clinical diagnosis. Relatedly, the utility of HC was exemplified by COVID-19, which was used as a case study to illustrate the maturity of HC for recovering pathogen sequences. As we unravel the origins of COVID-19, zoonoses remain more relevant than ever. Therefore, the role of HC in biosurveillance studies is also highlighted in this review, which is critical in preparing us for the next pandemic. We also found that while HC is a popular tool to study viruses, it remains underutilised in parasites and fungi and, to a lesser extent, bacteria. Finally, weevaluated the future of HC with respect to bait design in the eukaryotic groups and the prospect of combining HC with long-read HTS.

MASCOT-Skyline integrates population and migration dynamics to enhance phylogeographic reconstructions

The spread of infectious diseases is shaped by spatial and temporal aspects, such as host population structure or changes in the transmission rate or number of infected individuals over time. These spatiotemporal dynamics are imprinted in the genome of pathogens and can be recovered from those genomes using phylodynamics methods. However, phylodynamic methods typically quantify either the temporal or spatial transmission dynamics, which leads to unclear biases, as one can potentially not be inferred without the other. Here, we address this challenge by introducing a structured coalescent skyline approach, MASCOT-Skyline that allows us to jointly infer spatial and temporal transmission dynamics of infectious diseases using Markov chain Monte Carlo inference. To do so, we model the effective population size dynamics in different locations using a non-parametric function, allowing us to approximate a range of population size dynamics. We show, using a range of different viral outbreak datasets, potential issues with phylogeographic methods. We then use these viral datasets to motivate simulations of outbreaks that illuminate the nature of biases present in the different phylogeographic methods. We show that spatial and temporal dynamics should be modeled jointly even if one seeks to recover just one of the two. Further, we showcase conditions under which we can expect phylogeographic analyses to be biased, particularly different subsampling approaches, as well as provide recommendations of when we can expect them to perform well. We implemented MASCOT-Skyline as part of the open-source software package MASCOT for the Bayesian phylodynamics platform BEAST2.

[High-throughput sequencing for infectious disease diagnoses: Example of shotgun metagenomics in central nervous system infections]

The advent of high-throughput sequencing in clinical microbiology is opening the way to new diagnostic and prognostic approaches in infectious diseases. Detection, identification and characterisation of pathogenic microorganisms are essential steps in diagnosis and implementation of appropriate antimicrobial therapy. However, standard methods of microbiological diagnosis are failing in some cases. In addition, the emergence of new infections, facilitated by international travel and global warming, requires the implementation of innovative diagnostic methods. Among the different strategies used in clinical microbiology and reviewed in this article, shotgun metagenomics is the only technique that allows today a panpathogenic and unbiased detection of all microorganisms potentially responsible for an infectious disease, including those still unknown. The aims of this article are to present the different possible strategies of high-throughput sequencing used in the microbiological diagnosis of infectious diseases and to highlight the diagnostic contribution of shotgun metagenomics in the field of central nervous system infections.

Zika Virus-A Reemerging Neurotropic Arbovirus Associated with Adverse Pregnancy Outcomes and Neuropathogenesis

Zika virus (ZIKV) is a reemerging flavivirus that is primarily spread through bites from infected mosquitos. It was first discovered in 1947 in sentinel monkeys in Uganda and has since been the cause of several outbreaks, primarily in tropical and subtropical areas. Unlike earlier outbreaks, the 2015-2016 epidemic in Brazil was characterized by the emergence of neurovirulent strains of ZIKV strains that could be sexually and perinatally transmitted, leading to the Congenital Zika Syndrome (CZS) in newborns, and Guillain-Barre Syndrome (GBS) along with encephalitis and meningitis in adults. The immune response elicited by ZIKV infection is highly effective and characterized by the induction of both ZIKV-specific neutralizing antibodies and robust effector CD8+ T cell responses. However, the structural similarities between ZIKV and Dengue virus (DENV) lead to the induction of cross-reactive immune responses that could potentially enhance subsequent DENV infection, which imposes a constraint on the development of a highly efficacious ZIKV vaccine. The isolation and characterization of antibodies capable of cross-neutralizing both ZIKV and DENV along with cross-reactive CD8+ T cell responses suggest that vaccine immunogens can be designed to overcome these constraints. Here we review the structural characteristics of ZIKV along with the evidence of neuropathogenesis associated with ZIKV infection and the complex nature of the immune response that is elicited by ZIKV infection.

Development of a quantitative NS1 antigen enzyme-linked immunosorbent assay (ELISA) for Zika virus detection using a novel virus-specific mAb

Viruses from the Flaviviridae family, such as Dengue virus (DENV), Yellow fever virus (YFV), and Zika virus (ZIKV) are notorious global public health problems. ZIKV emergence in Polynesia and the Americas from 2013 to 2016 raised concerns as new distinguishing features set it apart from previous outbreaks, including its association with neurological complications and heightened disease severity. Virus detection is impaired as cross-reactivity to other closely related orthoflaviviruses is common among commercially available diagnostic kits. While non-structural protein 1 (NS1) has been used as an early marker of DENV and West Nile virus (WNV) infection, little is known about NS1 expression during ZIKV infection. In the present work, we developed a NS1 capture ELISA using a novel ZIKV-specific monoclonal antibody to study NS1 expression dynamics in vitro in mosquito and human cell lines. While detectable in culture supernatants, higher concentrations of NS1 were predominantly cell-associated. To our knowledge, this is the first report of NS1 detection in human cells despite viral clearance over time. Tests with human samples need to be conducted to validate the applicability of NS1 detection for diagnosis, but overall, the tools developed in this work are promising for specific detection of acute ZIKV infection.

Tracking arboviruses, their transmission vectors and potential hosts by nanopore sequencing of mosquitoes

The risk to human health from mosquito-borne viruses such as dengue, chikungunya and yellow fever is increasing due to increased human expansion, deforestation and climate change. To anticipate and predict the spread and transmission of mosquito-borne viruses, a better understanding of the transmission cycle in mosquito populations is needed. We present a pathogen-agnostic combined sequencing protocol for identifying vectors, viral pathogens and their hosts or reservoirs using portable Oxford Nanopore sequencing. Using mosquitoes collected in São Paulo, Brazil, we extracted RNA for virus identification and DNA for blood meal and mosquito identification. Mosquitoes and blood meals were identified by comparing cytochrome c oxidase I (COI) sequences against a curated Barcode of Life Data System (BOLD). Viruses were identified using the SMART-9N protocol, which allows amplified DNA to be prepared with native barcoding for nanopore sequencing. Kraken 2 was employed to detect viral pathogens and Minimap2 and BOLD identified the contents of the blood meal. Due to the high similarity of some species, mosquito identification was conducted using blast after generation of consensus COI sequences using RACON polishing. This protocol can simultaneously uncover viral diversity, mosquito species and mosquito feeding habits. It also has the potential to increase understanding of mosquito genetic diversity and transmission dynamics of zoonotic mosquito-borne viruses.

Interactions between climate change, urban infrastructure and mobility are driving dengue emergence in Vietnam

Dengue is expanding globally, but how dengue emergence is shaped locally by interactions between climatic and socio-environmental factors is not well understood. Here, we investigate the drivers of dengue incidence and emergence in Vietnam, through analysing 23 years of district-level case data spanning a period of significant socioeconomic change (1998-2020). We show that urban infrastructure factors (sanitation, water supply, long-term urban growth) predict local spatial patterns of dengue incidence, while human mobility is a more influential driver in subtropical northern regions than the endemic south. Temperature is the dominant factor shaping dengue's distribution and dynamics, and using long-term reanalysis temperature data we show that warming since 1950 has expanded transmission risk throughout Vietnam, and most strongly in current dengue emergence hotspots (e.g., southern central regions, Ha Noi). In contrast, effects of hydrometeorology are complex, multi-scalar and dependent on local context: risk increases under either short-term precipitation excess or long-term drought, but improvements in water supply mitigate drought-associated risks except under extreme conditions. Our findings challenge the assumption that dengue is an urban disease, instead suggesting that incidence peaks in transitional landscapes with intermediate infrastructure provision, and provide evidence that interactions between recent climate change and mobility are contributing to dengue's expansion throughout Vietnam.

A hybridization target enrichment approach for pathogen genomics

IMPORTANCE

Emerging infectious diseases require continuous pathogen monitoring. Rapid clinical diagnosis by nucleic acid amplification is limited to a small number of targets and may miss target detection due to new mutations in clinical isolates. Whole-genome sequencing (WGS) identifies genome-wide variations that may be used to determine a pathogen's drug resistance patterns and phylogenetically characterize isolates to track disease origin and transmission. WGS is typically performed using DNA isolated from cultured clinical isolates. Culturing clinical specimens increases turn-around time and may not be possible for fastidious bacteria. To overcome some of these limitations, direct sequencing of clinical specimens has been attempted using expensive capture probes to enrich the entire genomes of target pathogens. We present a method to produce a cost-effective, time-efficient, and large-scale synthesis of probes for whole-genome enrichment. We envision that our method can be used for direct clinical sequencing of a wide range of microbial pathogens for genomic epidemiology.

Two invasions at once: update on the introduction of the invasive species Aedes aegypti and Aedes albopictus in Cyprus - a call for action in Europe

Aedes aegypti, the yellow fever mosquito and Aedes albopictus, the tiger mosquito, continue to expand their geographical distribution, reshaping the European epidemiological risks for mosquito-borne diseases. The reintroduction of Aedes aegypti near the airport and port in Larnaka and the detection of Aedes albopictus near the marina and old port of the Limassol area in Cyprus are reported herein. The measures taken to investigate these events included (i) communication to health authorities, (ii) expert on-site visits and verification of findings, (iii) enhanced active surveillance, and (iv) development of an Emergency Action Plan followed by a Contingency Plan. These emergency action plans were developed to delimitate the infested areas and to prevent the spreading of the mosquito populations into new areas. The general principles are presented along with their rationale to serve as guidelines for other geographical regions targeting suppression/eradication with a sterile insect technique component. In parallel, this manuscript serves as a call for action at the European level to impede the further spread of these species and support the activities being undertaken in Cyprus to combat the incursions of Aedes invasive species.

Highly multiplexed targeted sequencing strategy for infectious disease surveillance

BACKGROUND

Global efforts to characterize diseases of poverty are hampered by lack of affordable and comprehensive detection platforms, resulting in suboptimal allocation of health care resources and inefficient disease control. Next generation sequencing (NGS) can provide accurate data and high throughput. However, shotgun and metagenome-based NGS approaches are limited by low concentrations of microbial DNA in clinical samples, requirements for tailored sample and library preparations plus extensive bioinformatics analysis. Here, we adapted molecular inversion probes (MIPs) as a cost-effective target enrichment approach to characterize microbial infections from blood samples using short-read sequencing. We designed a probe panel targeting 2 bacterial genera, 21 bacterial and 6 fungi species and 7 antimicrobial resistance markers (AMRs).

RESULTS

Our approach proved to be highly specific to detect down to 1 in a 1000 pathogen DNA targets contained in host DNA. Additionally, we were able to accurately survey pathogens and AMRs in 20 out of 24 samples previously profiled with routine blood culture for sepsis.

CONCLUSIONS

Overall, our targeted assay identifies microbial pathogens and AMRs with high specificity at high throughput, without the need for extensive sample preparation or bioinformatics analysis, simplifying its application for characterization and surveillance of infectious diseases in medium- to low- resource settings.

Genomic epidemiology of dengue in Shantou, China, 2019

Objectives

Dengue has been endemic in Southeast Asian countries for decades. There are few reports tracing the dynamics of dengue in real time. In this study, we generated hundreds of pathogen genomes to understand the genomic epidemiology of an outbreak in a hyper-endemic area of dengue.

Methods

We leveraged whole-genome short-read sequencing (PE150) to generate genomes of the dengue virus and investigated the genomic epidemiology of a dengue virus transmission in a mesoscale outbreak in Shantou, China, in 2019.

Results

The outbreak was sustained from July to December 2019. The total accumulated number of laboratory-confirmed cases was 944. No gender bias or fatalities were recorded. Cambodia and Singapore were the main sources of imported dengue cases (74.07%, n = 20). A total of 284 dengue virus strains were isolated, including 259 DENV-1, 24 DENV-2, and 1 DENV-3 isolates. We generated the entire genome of 252 DENV isolates (229 DENV-1, 22 DENV-2, and 1 DENV-3), which represented 26.7% of the total cases. Combined epidemiological and phylogenetic analyses indicated multiple independent introductions. The internal transmission evaluations and transmission network reconstruction supported the inference of phylodynamic analysis, with high Bayes factor support in BSSVS analysis. Two expansion founders and transmission chains were detected in CCH and LG of Shantou.

Conclusions

We observed the instant effects of genomic epidemiology in monitoring the dynamics of DENV and highlighted its prospects for real-time tracing of outbreaks of other novel agents in the future.

Analyses of Early ZIKV Genomes Are Consistent with Viral Spread from Northeast Brazil to the Americas

The Americas, particularly Brazil, were greatly impacted by the widespread Zika virus (ZIKV) outbreak in 2015 and 2016. Efforts were made to implement genomic surveillance of ZIKV as part of the public health responses. The accuracy of spatiotemporal reconstructions of the epidemic spread relies on the unbiased sampling of the transmission process. In the early stages of the outbreak, we recruited patients exhibiting clinical symptoms of arbovirus-like infection from Salvador and Campo Formoso, Bahia, in Northeast Brazil. Between May 2015 and June 2016, we identified 21 cases of acute ZIKV infection and subsequently recovered 14 near full-length sequences using the amplicon tiling multiplex approach with nanopore sequencing. We performed a time-calibrated discrete phylogeographic analysis to trace the spread and migration history of the ZIKV. Our phylogenetic analysis supports a consistent relationship between ZIKV migration from Northeast to Southeast Brazil and its subsequent dissemination beyond Brazil. Additionally, our analysis provides insights into the migration of ZIKV from Brazil to Haiti and the role Brazil played in the spread of ZIKV to other countries, such as Singapore, the USA, and the Dominican Republic. The data generated by this study enhances our understanding of ZIKV dynamics and supports the existing knowledge, which can aid in future surveillance efforts against the virus.

A scoping review of waterborne and water-related disease in the Florida environment from 1999 to 2022

Florida's environments are suitable reservoirs for many disease-causing agents. Pathogens and toxins in Florida waterways have the potential to infect mosquito vectors, animals, and human hosts. Through a scoping review of the scientific literature published between 1999 and 2022, we examined the presence of water-related pathogens, toxins, and toxin-producers in the Florida environment and the potential risk factors for human exposure. Nineteen databases were searched using keywords relating to the waterborne, water-based toxins, and water-related vector-borne diseases which are reportable to the Florida Department of Health. Of the 10,439 results, 84 titles were included in the final qualitative analysis. The resulting titles included environmental samples of water, mosquitoes, algae, sand, soil/sediment, air, food, biofilm, and other media. Many of the waterborne, water-related vector-borne, and water-based toxins and toxin-producers of public health and veterinary importance from our search were found to be present in Florida environments. Interactions with Florida waterways can expose humans and animals to disease and toxins due to nearby human and/or animal activity, proximal animal or human waste, failing or inadequate water and/or sanitation, weather patterns, environmental events, and seasonality, contaminated food items, preference of agent for environmental media, high-risk populations, urban development and population movement, and unregulated and unsafe environmental activities. A One Health approach will be imperative to maintaining healthy waterways and shared environments throughout the state to protect the health of humans, animals, and our ecosystems.

Spatiotemporal distribution of vector mosquito species and areas at risk for arbovirus transmission in Maricopa County, Arizona

Mosquito-borne diseases are a major global public health concern and mosquito surveillance systems are essential for the implementation of effective mosquito control strategies. The objective of our study is to determine the spatiotemporal distribution of vector mosquito species in Maricopa County, AZ from 2011 to 2021, and to identify the hotspot areas for West Nile virus (WNV) and St. Louis Encephalitis virus (SLEV) transmission in 2021. The Maricopa County Mosquito Control surveillance system utilizes BG-Sentinel and EVS-CDC traps throughout the entire urban and suburban areas of the county. We estimated specific mosquito species relative abundance per unit area using the Kernel density estimator in ArcGIS 10.2. We calculated the distance between all traps in the surveillance system and created a 4 km buffer radius around each trap to calculate the extent to which each trap deviated from the mean number of Culex quinquefasciatus and Culex tarsalis collected in 2021. Our results show that vector mosquito species are widely distributed and abundant in the urban areas of Maricopa County. A total of 691,170Cx. quinquefasciatus, 542,733 Cx. tarsalis, and 292,305 Aedes aegypti were collected from 2011 to 2022. The relative abundance of Ae. aegypti was highly seasonal peaking in the third and fourth quarters of the year. Culex quinquefasciatus, on the other hand, was abundant throughout the year with several regions consistently yielding high numbers of mosquitoes. Culex tarsalis was abundant but it only reached high numbers in well-defined areas near irrigated landscapes. We also detected high levels of heterogeneity in the risk of WNV and SLEV transmission to humans disregarding traps geographical proximity. The well-defined species-specific spatiotemporal and geographical patterns found in this study can be used to inform vector control operations.

Outbreak.info genomic reports: scalable and dynamic surveillance of SARS-CoV-2 variants and mutations

In response to the emergence of SARS-CoV-2 variants of concern, the global scientific community, through unprecedented effort, has sequenced and shared over 11 million genomes through GISAID, as of May 2022. This extraordinarily high sampling rate provides a unique opportunity to track the evolution of the virus in near real-time. Here, we present outbreak.info , a platform that currently tracks over 40 million combinations of Pango lineages and individual mutations, across over 7,000 locations, to provide insights for researchers, public health officials and the general public. We describe the interpretable visualizations available in our web application, the pipelines that enable the scalable ingestion of heterogeneous sources of SARS-CoV-2 variant data and the server infrastructure that enables widespread data dissemination via a high-performance API that can be accessed using an R package. We show how outbreak.info can be used for genomic surveillance and as a hypothesis-generation tool to understand the ongoing pandemic at varying geographic and temporal scales.

High-throughput sequencing approaches applied to SARS-CoV-2

High-throughput sequencing is crucial for surveillance and control of viral outbreaks. During the ongoing coronavirus disease 2019 (COVID-19) pandemic, advances in the high-throughput sequencing technology resources have enhanced diagnosis, surveillance, and vaccine discovery. From the onset of the pandemic in December 2019, several genome-sequencing approaches have been developed and supported across the major sequencing platforms such as Illumina, Oxford Nanopore, PacBio, MGI DNBSEQTM and Ion Torrent. Here, we share insights from the sequencing approaches developed for sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) between December 2019 and October 2022.

Ethical Challenges Associated with Pathogen and Host Genetics in Infectious Disease

The Covid-19 pandemic has demonstrated the potential of genomic technologies for the detection and surveillance of infectious diseases. Pathogen genomics is likely to play a major role in the future of research and clinical implementation of genomic technologies. However, unlike human genetics, the specific ethical and social challenges associated with the implementation of infectious disease genomics has received comparatively little attention. In this paper, we contribute to this literature, focusing on the potential consequences for individuals and communities of the use of these technologies. We concentrate on areas of challenges related to privacy, stigma, discrimination and the return of results in the cases of the surveillance of known pathogens, metagenomics and host genomics.

Using Multiplex Amplicon PCR Technology to Efficiently and Timely Generate Rift Valley Fever Virus Sequence Data for Genomic Surveillance

Rift Valley fever (RVF) is a febrile vector-borne disease endemic in Africa and continues to spread in new territories. It is a climate-sensitive disease mostly triggered by abnormal rainfall patterns. The disease is associated with high mortality and morbidity in both humans and livestock. RVF is caused by the Rift Valley fever virus (RVFV) of the genus Phlebovirus in the family Phenuiviridae. It is a tripartite RNA virus with three genomic segments: small (S), medium (M) and large (L). Pathogen genomic sequencing is becoming a routine procedure and a powerful tool for understanding the evolutionary dynamics of infectious organisms, including viruses. Inspired by the utility of amplicon-based sequencing demonstrated in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and Ebola, Zika and West Nile viruses, we report an RVFV sample preparation based on amplicon multiplex polymerase chain reaction (amPCR) for template enrichment and reduction of background host contamination. The technology can be implemented rapidly to characterize and genotype RVFV during outbreaks in a near-real-time manner. To achieve this, we designed 74 multiplex primer sets covering the entire RVFV genome to specifically amplify the nucleic acid of RVFV in clinical samples from an animal tissue. Using this approach, we demonstrate achieving complete RVFV genome coverage even from samples containing a relatively low viral load. We report the first primer scheme approach of generating multiplex primer sets for a tripartite virus which can be replicated for other segmented viruses.

Leon K, Ontiveros CO, Agarwal R, Tsetsarkin KA, Harris E, Ebel GD, Wohl S, Grubaugh ND, Andersen KG. 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.

Rapid range shifts in African Anopheles mosquitoes over the last century

Facing a warming climate, many tropical species-including the arthropod vectors of several infectious diseases-will be displaced to higher latitudes and elevations. These shifts are frequently projected for the future, but rarely documented in the present day. Here, we use one of the most comprehensive datasets ever compiled by medical entomologists to track the observed range limits of African malaria mosquito vectors (Anopheles spp.) from 1898 to 2016. Using a simple regression approach, we estimate that these species' ranges gained an average of 6.5 m of elevation per year, and the southern limits of their ranges moved polewards 4.7 km per year. These shifts would be consistent with the local velocity of recent climate change, and might help explain the incursion of malaria transmission into new areas over the past few decades. Confirming that climate change underlies these shifts, and applying similar methods to other disease vectors, are important directions for future research.

Real-time digital data of international passengers will shine in the precaution of epidemics

International movement plays an important role in spatial spread of infectious diseases. Here, we share two successful COVID-19 interventions based on real-time digital information collected from international passengers, which have been performed in Greece and China respectively. Both of the interventions demonstrated good performance and showed the potential of real-time digital data in containing the spread. However, several key points should not be ignored when we promote similar strategies.

The Zika virus: an opportunity to revisit reproductive health needs and disparities

First isolated in 1947, the Zika virus was initially connected only to limited or sporadic human infections. In late 2015, the temporal clustering of a Zika outbreak and microcephaly in newborn babies from northeastern Brazil, and the identification of a causal link between the two, led to the characterization of the congenital Zika syndrome. In the wake of the epidemic, several countries from Latin America advised women to postpone pregnancies for periods ranging from six months to two years. These recommendations initiated critical conversations about the challenges of implementing them in societies with limited access to contraception, widespread socioeconomic inequalities, and high rates of unplanned and adolescent pregnancies. The messaging targeted exclusively women, despite a high prevalence of imbalances in the relationship power, and addressed all women as a group, failing to recognize that the decision to postpone pregnancies will impact different women in different ways, depending on their age at the time. Finally, in several countries affected by the Zika epidemic, due to restrictive reproductive policies, legally terminating a pregnancy is no longer an option even at the earliest time when brain malformations as part of the congenital Zika syndrome can be detected by ultrasonography. The virus continued to circulate after 2016 in several countries. Climate change models predict an expansion of the geographical area where local Zika transmission may occur, indicating that the interface between the virus, teratogenesis, and reproductive rights is a topic of considerable interest for medicine, social sciences, and public health for years to come.

Repeated Rapid Active Sampling Surveys Demonstrated a Rapidly Changing Zika Seroprevalence among Children in a Rural Dengue-endemic Region in Southwest Guatemala during the Zika Epidemic (2015-2016)

Although Central America is largely dengue virus (DENV)-endemic, the 2015-2016 Zika virus (ZIKV) pandemic brought new urgency to develop surveillance approaches capable of characterizing the rapidly changing disease burden in resource-limited settings. We conducted a pediatric DENV surveillance study in rural Guatemala, including serial cross-sectional surveys from April through September 2015 (Survey 1), in October-November 2015 (Survey 2), and January-February 2016 (Survey 3). Serum underwent DENV IgM MAC ELISA and polymerase chain reaction testing. Using banked specimens from Surveys 2 and 3, we expanded testing to include DENV 1-4 and ZIKV microneutralization (MN50), DENV NS1 IgG ELISA, and ZIKV anti-NS1 antibody Blockage of Binding (BoB) ELISA testing. Demographic risk factors for ZIKV BoB positivity were explored using multivariable generalized linear regression models. Of Survey 2 and 3 samples available (N = 382), DENV seroprevalence slightly increased (+1%-10% depending on the assay) during the surveillance period and increased with age. In contrast, ZIKV seroprevalence consistently increased over the 3-month period, including from 6% to 34% (P < 0.0001) and 10%-37% (P < 0.0001) using the MN50 ≥100 and BoB ELISA assays, respectively. Independent risk factors for ZIKV seropositivity included older age (prevalence ratio (PR)/year = 1.12, 95% confidence interval (CI) = 1.07-1.17) and primary caregiver literacy (PR = 2.80, CI = 1.30-6.06). Rapid active surveillance (RAS) surveys demonstrated a nearly 30% increase in ZIKV prevalence and a slight (≤ 10%) increase in DENV seroprevalence from October to November 2015 to January to February 2016 in rural southwest Guatemala, regardless of serologic assay used. RAS surveys may be a useful "off-the-shelf" tool to characterize arboviruses and other emerging pathogens rapidly in resource-limited settings.

Neutralizing Immunity Induced Against the Omicron BA.1 and BA.2 Variants in Vaccine Breakthrough Infections

BACKGROUND

As of early 2022, the Omicron variants are the predominant circulating lineages globally. Understanding neutralizing antibody responses against Omicron BA.1 and BA.2 after vaccine breakthrough infections will provide insights into BA.2 infectivity and susceptibility to subsequent reinfection.

METHODS

Live virus neutralization assays were used to study immunity against Delta and Omicron BA.1 and BA.2 variants in samples from 86 individuals, 24 unvaccinated (27.9%) and 62 vaccinated (72.1%), who were infected with Delta (n = 42, 48.8%) or BA.1 (n = 44, 51.2%). Among the 62 vaccinated individuals, 39 were unboosted (62.9%), whereas 23 were boosted (37.1%).

RESULTS

In unvaccinated infections, neutralizing antibodies (nAbs) against the three variants were weak or undetectable, except against Delta for Delta-infected individuals. Both Delta and BA.1 breakthrough infections resulted in strong nAb responses against ancestral wild-type and Delta lineages, but moderate nAb responses against BA.1 and BA.2, with similar titers between unboosted and boosted individuals. Antibody titers against BA.2 were generally higher than those against BA.1 in breakthrough infections.

CONCLUSIONS

These results underscore the decreased immunogenicity of BA.1 compared to BA.2, insufficient neutralizing immunity against BA.2 in unvaccinated individuals, and moderate to strong neutralizing immunity induced against BA.2 in Delta and BA.1 breakthrough infections.

Tracing the international arrivals of SARS-CoV-2 Omicron variants after Aotearoa New Zealand reopened its border

In the second quarter of 2022, there was a global surge of emergent SARS-CoV-2 lineages that had a distinct growth advantage over then-dominant Omicron BA.1 and BA.2 lineages. By generating 10,403 Omicron genomes, we show that Aotearoa New Zealand observed an influx of these immune-evasive variants (BA.2.12.1, BA.4, and BA.5) through the border. This is explained by the return to significant levels of international travel following the border's reopening in March 2022. We estimate one Omicron transmission event from the border to the community for every ~5,000 passenger arrivals at the current levels of travel and restriction. Although most of these introductions did not instigate any detected onward transmission, a small minority triggered large outbreaks. Genomic surveillance at the border provides a lens on the rate at which new variants might gain a foothold and trigger new waves of infection.

Analysis of Zika Virus Sequence Data Associated with a School Cohort in Haiti

Zika virus (ZIKV) infections occurred in epidemic form in the Americas in 2014-2016, with some of the earliest isolates in the region coming from Haiti. We isolated ZIKV from 20 children with acute undifferentiated febrile illness who were part of a cohort of children seen at a school clinic in the Gressier region of Haiti. The virus was also isolated from three pools of Aedes aegypti mosquitoes collected at the same location. On phylogenetic analysis, three distinct ZIKV clades were identified. Strains from all three clades were present in Haiti in 2014, making them among the earliest isolates identified in the Western Hemisphere. Strains from all three clades were also isolated in 2016, indicative of their persistence across the time period of the epidemic. Mosquito isolates were collected in 2016 and included representatives from two of the three clades; in one instance, ZIKV was isolated from a pool of male mosquitoes, suggestive of vertical transmission of the virus. The identification of multiple ZIKV clades in Haiti at the beginning of the epidemic suggests that Haiti served as a nidus for transmission within the Caribbean.

Using multiple sampling strategies to estimate SARS-CoV-2 epidemiological parameters from genomic sequencing data

The choice of viral sequences used in genetic and epidemiological analysis is important as it can induce biases that detract from the value of these rich datasets. This raises questions about how a set of sequences should be chosen for analysis. We provide insights on these largely understudied problems using SARS-CoV-2 genomic sequences from Hong Kong, China, and the Amazonas State, Brazil. We consider multiple sampling schemes which were used to estimate Rt and rt as well as related R0 and date of origin parameters. We find that both Rt and rt are sensitive to changes in sampling whilst R0 and the date of origin are relatively robust. Moreover, we find that analysis using unsampled datasets result in the most biased Rt and rt estimates for both our Hong Kong and Amazonas case studies. We highlight that sampling strategy choices may be an influential yet neglected component of sequencing analysis pipelines.

Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic

Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.

Culture-Free Phylogenetic Analysis of Legionella pneumophila Using Targeted CRISPR/Cas9 Next-Generation Sequencing

Currently available methods for the laboratory investigation of Legionella pneumophila outbreaks require organism culture. The ability to sequence L. pneumophila directly from clinical samples would significantly reduce delays. Here, we develop a method for targeted next-generation sequencing (NGS) of selected L. pneumophila genes utilizing a CRISPR/Cas9-based target enrichment system. We determine the method's utility by typing cultured L. pneumophila isolates and subsequently apply the method directly to patient samples. We sequenced 10 L. pneumophila isolates by 2 methods, (i) whole-genome sequencing (WGS) and (ii) targeted (CRISPR/Cas9-based) finding low-abundance sequences by hybridization (FLASH)-NGS, sequencing 57 selected genes. The targeted NGS of 57 genes was more efficient than WGS, and phylogenetic analysis of the 57 genes yielded the same classification of the L. pneumophila isolates as that based on analysis of whole-genome data. Furthermore, targeted NGS of L. pneumophila performed directly on patient respiratory samples correctly classified the patients according to their corresponding cultured isolates. This provides proof of concept that targeted NGS can be used to sequence L. pneumophila directly from patient samples. Studies on a larger number of patient samples will further validate this method. Nonetheless, CRISPR/Cas9 targeted NGS methods have the potential to be widely applicable to microbial-outbreak investigations in the future, particularly in the context of difficult and slow-growing organisms. IMPORTANCE The bacterium Legionella pneumophila is responsible for outbreaks of serious and life-threatening pneumonia called Legionnaires' disease. There is a need for new molecular methods that allow investigation of Legionella outbreaks directly from patient samples, without the need for prior microbiological culture, which causes delays. Our study aims to address this problem. We have utilized a CRISPR/Cas9-based targeted next-generation sequencing (NGS) method that can be applied directly on human specimens. Furthermore, we show that analysis of the sequences of a small number of targeted genes offers the same classification of L. pneumophila as that based on data derived from the whole genome. Given the rising interest globally in sequencing pathogens directly from human samples, CRISPR/Cas9 targeted NGS methods have the potential to be widely applicable to microbial-outbreak investigations in the future, particularly in the context of difficult and slow-growing organisms.

Outbreak.info genomic reports: scalable and dynamic surveillance of SARS-CoV-2 variants and mutations

The emergence of SARS-CoV-2 variants of concern has prompted the need for near real-time genomic surveillance to inform public health interventions. In response to this need, the global scientific community, through unprecedented effort, has sequenced and shared over 11 million genomes through GISAID, as of May 2022. This extraordinarily high sampling rate provides a unique opportunity to track the evolution of the virus in near real-time. Here, we present outbreak.info, a platform that currently tracks over 40 million combinations of PANGO lineages and individual mutations, across over 7,000 locations, to provide insights for researchers, public health officials, and the general public. We describe the interpretable and opinionated visualizations in the variant and location focussed reports available in our web application, the pipelines that enable the scalable ingestion of heterogeneous sources of SARS-CoV-2 variant data, and the server infrastructure that enables widespread data dissemination via a high performance API that can be accessed using an R package. We present a case study that illustrates how outbreak.info can be used for genomic surveillance and as a hypothesis generation tool to understand the ongoing pandemic at varying geographic and temporal scales. With an emphasis on scalability, interactivity, interpretability, and reusability, outbreak.info provides a template to enable genomic surveillance at a global and localized scale.

Mosquito surveillance in maritime entry ports in Miami-Dade County, Florida to increase preparedness and allow the early detection of invasive mosquito species

Invasive mosquito vector species have been inadvertently transported to new areas by humans for decades. Strong evidence supports that monitoring maritime, terrestrial, and aerial points of entry is an essential part of the effort to curb the invasion and establishment of invasive vector mosquito species. Miami-Dade County, Florida is an important operational hub for the cruise ship industry and leisure boats that routinely visit nearby areas in the Caribbean, and freight cargo ships transporting goods from Miami-Dade to Caribbean countries and vice versa. To deal with the increasing public health concern, we hypothesized that mosquito surveillance in small- and medium-sized maritime ports of entry in Miami-Dade is crucial to allow the early detection of invasive mosquito species. Therefore, we have selected 12 small- and medium-sized maritime ports of entry in Miami-Dade County with an increased flow of people and commodities that were not covered by the current mosquito surveillance system. Collection sites were comprised of two distinct environments, four marinas with international traffic of leisure boats, and eight maintenance and commercial freight cargo ship ports. Mosquitoes were collected weekly at each of the 12 collection sites for 24 hours for 6 weeks in the Spring and then for 6 additional weeks in the Summer using BG-Sentinel traps. A total of 32,590 mosquitoes were collected, with Culex quinquefasciatus and Aedes aegypti being the most abundant species totaling 19,987 and 11,247 specimens collected, respectively. Our results show that important mosquito vector species were present in great numbers in all of the 12 maritime ports of entry surveyed during this study. The relative abundance of Cx. quinquefasciatus and Ae. aegypti was substantially higher in the commercial freight cargo ship ports than in the marinas. These results indicate that even though both areas are conducive for the proliferation of vector mosquitoes, the port area in the Miami River is especially suitable for the proliferation of vector mosquitoes. Therefore, this potentially allows the establishment of invasive mosquito species inadvertently brought in by cargo freights.

Wastewater sequencing uncovers early, cryptic SARS-CoV-2 variant transmission

As SARS-CoV-2 continues to spread and evolve, detecting emerging variants early is critical for public health interventions. Inferring lineage prevalence by clinical testing is infeasible at scale, especially in areas with limited resources, participation, or testing/sequencing capacity, which can also introduce biases. SARS-CoV-2 RNA concentration in wastewater successfully tracks regional infection dynamics and provides less biased abundance estimates than clinical testing. Tracking virus genomic sequences in wastewater would improve community prevalence estimates and detect emerging variants. However, two factors limit wastewater-based genomic surveillance: low-quality sequence data and inability to estimate relative lineage abundance in mixed samples. Here, we resolve these critical issues to perform a high-resolution, 295-day wastewater and clinical sequencing effort, in the controlled environment of a large university campus and the broader context of the surrounding county. We develop and deploy improved virus concentration protocols and deconvolution software that fully resolve multiple virus strains from wastewater. We detect emerging variants of concern up to 14 days earlier in wastewater samples, and identify multiple instances of virus spread not captured by clinical genomic surveillance. Our study provides a scalable solution for wastewater genomic surveillance that allows early detection of SARS-CoV-2 variants and identification of cryptic transmission.

Investigation of cross-regional spread and evolution of equine influenza H3N8 at US and global scales using Bayesian phylogeography based on balanced subsampling

Equine influenza virus (EIV) is a highly contagious pathogen of equids, and a well-known burden in global equine health. EIV H3N8 variants seasonally emerged and resulted in EIV outbreaks in the United States (US) and worldwide. The present study evaluated the pattern of cross-regional EIV H3N8 spread and evolutionary characteristics at US and global scales using Bayesian phylogeography with balanced subsampling based on regional horse population size. A total of 297 Haemagglutinin (HA) sequences of global EIV H3N8 were collected from 1963 to 2019 and subsampled to global subset (n = 67), raw US sequences (n = 100) and US subset (n = 44) datasets. Discrete trait phylogeography analysis was used to estimate the transmission history of EIV using four global and US genome datasets. The North American lineage was the major source of globally dominant EIV variants and spread to other global regions. The US EIV strains generally spread from the southern and midwestern regions to other regions. The EIV H3N8 accumulated approximately three nucleotide substitutions per year in the HA gene under heterogenous local positive selection. Our findings will guide better decision making of target intervention strategies of EIV H3N8 infection and provide the better scheme of genomic surveillance in the US and global equine health. This article is protected by copyright. All rights reserved.

Intrinsic variation in the vertically transmitted core virome of the mosquito Aedes aegypti

Virome studies among metazoans have revealed the ubiquity of RNA viruses in animals, contributing to a fundamental re-thinking of the relationships between organisms and their microbiota. Mosquito viromes, often scrutinized due to their public health relevance, may also provide insight into broadly applicable concepts, such as a "core virome," a set of viruses consistently associated with a host species or population that may fundamentally impact its basic biology. A subset of mosquito-associated viruses (MAVs) could comprise such a core, and MAVs can be categorized as (i) arboviruses, which alternate between mosquito and vertebrate hosts, (ii) insect-specific viruses, which cannot replicate in vertebrate cells, and (iii) viruses with unknown specificity. MAVs have been widely characterized in the disease vector Aedes aegypti, and the occurrence of a core virome in this species has been proposed but remains unclear. Using a wild population previously surveyed for MAVs and a common laboratory strain, we investigated viromes in reproductive tissue via metagenomic RNA sequencing. Virome composition varied across samples, but four groups comprised >97% of virus sequences: a novel partiti-like virus (Partitiviridae), a toti-like virus (Totiviridae), unclassified Riboviria, and four orthomyxo-like viruses (Orthormyxoviridae). Whole or partial genomes for the partiti-like virus, toti-like virus, and one orthomyxo-like virus were assembled and analyzed phylogenetically. Multigenerational maintenance of these MAVs was confirmed by RT-PCR, indicating vertical transmission as a mechanism for persistence. This study provides fundamental information regarding MAV ecology and variability in A. aegypti and the potential for vertically maintained core viromes at the population level.

Urbanization favors the proliferation of Aedes aegypti and Culex quinquefasciatus in urban areas of Miami-Dade County, Florida

Urbanization processes are increasing globally. Anthropogenic alterations in the environment have profound effects on biodiversity. Decreased biodiversity due to biotic homogenization processes as a consequence of urbanization often result in increased levels of mosquito vector species and vector-borne pathogen transmission. Understanding how anthropogenic alterations in the environment will affect the abundance, richness, and composition of vector mosquito species is crucial for the implementation of effective and targeted mosquito control strategies. We hypothesized that anthropogenic alterations in the environment are responsible for increasing the abundance of mosquito species that are adapted to urban environments such as Aedesaegypti and Culexquinquefasciatus. Therefore, our objective was to survey mosquito relative abundance, richness, and community composition in Miami-Dade County, Florida, in areas with different levels of urbanization. We selected 24 areas, 16 remote areas comprised of natural and rural areas, and 8 urban areas comprised of residential and touristic areas in Miami-Dade County, Florida. Mosquitoes were collected weekly in each area for 24 h for 5 consecutive weeks from August to October 2020 using BG-Sentinel traps baited with dry ice. A total of 36,645 mosquitoes were collected, from which 34,048 were collected in the remote areas and 2,597 in the urban areas. Our results show a clear and well-defined pattern of abundance, richness, and community composition according to anthropogenic modifications in land use and land cover. The more urbanized a given area the fewer species were found and those were primary vectors of arboviruses, Ae.aegypti and Cx.quinquefasciatus.

Genomic diversity of SARS-CoV-2 in Oxford during United Kingdom's first national lockdown

Epidemiological efforts to model the spread of SARS-CoV-2, the virus that causes COVID-19, are crucial to understanding and containing current and future outbreaks and to inform public health responses. Mutations that occur in viral genomes can alter virulence during outbreaks by increasing infection rates and helping the virus evade the host immune system. To understand the changes in viral genomic diversity and molecular epidemiology in Oxford during the first wave of infections in the United Kingdom, we analyzed 563 clinical SARS-CoV-2 samples via whole-genome sequencing using Nanopore MinION sequencing. Large-scale surveillance efforts during viral epidemics are likely to be confounded by the number of independent introductions of the viral strains into a region. To avoid such issues and better understand the selection-based changes occurring in the SARS-CoV-2 genome, we utilized local isolates collected during the UK's first national lockdown whereby personal interactions, international and national travel were considerably restricted and controlled. We were able to track the short-term evolution of the virus, detect the emergence of several mutations of concern or interest, and capture the viral diversity of the region. Overall, these results demonstrate genomic pathogen surveillance efforts have considerable utility in controlling the local spread of the virus.

Genomic epidemiology reveals multiple introductions of SARS-CoV-2 followed by community and nosocomial spread, Germany, February to May 2020

Background

In the SARS-CoV-2 pandemic, viral genomes are available at unprecedented speed, but spatio-temporal bias in genome sequence sampling precludes phylogeographical inference without additional contextual data.

Aim

We applied genomic epidemiology to trace SARS-CoV-2 spread on an international, national and local level, to illustrate how transmission chains can be resolved to the level of a single event and single person using integrated sequence data and spatio-temporal metadata.

Methods

We investigated 289 COVID-19 cases at a university hospital in Munich, Germany, between 29 February and 27 May 2020. Using the ARTIC protocol, we obtained near full-length viral genomes from 174 SARS-CoV-2-positive respiratory samples. Phylogenetic analyses using the Auspice software were employed in combination with anamnestic reporting of travel history, interpersonal interactions and perceived high-risk exposures among patients and healthcare workers to characterise cluster outbreaks and establish likely scenarios and timelines of transmission.

Results

We identified multiple independent introductions in the Munich Metropolitan Region during the first weeks of the first pandemic wave, mainly by travellers returning from popular skiing areas in the Alps. In these early weeks, the rate of presumable hospital-acquired infections among patients and in particular healthcare workers was high (9.6% and 54%, respectively) and we illustrated how transmission chains can be dissected at high resolution combining virus sequences and spatio-temporal networks of human interactions.

Conclusions

Early spread of SARS-CoV-2 in Europe was catalysed by superspreading events and regional hotspots during the winter holiday season. Genomic epidemiology can be employed to trace viral spread and inform effective containment strategies.

Imidacloprid Crystal Polymorphs for Disease Vector Control and Pollinator Protection

Imidacloprid, the world's leading insecticide, has been approved recently for controlling infectious disease vectors; yet, in agricultural settings, it has been implicated in the frightening decline of pollinators. This argues for strategies that sharply reduce the environmental impact of imidacloprid. When used as a contact insecticide, the effectiveness of imidacloprid relies on physical contact between its crystal surfaces and insect tarsi. Herein, seven new imidacloprid crystal polymorphs are reported, adding to two known forms. Anticipating that insect uptake of imidacloprid molecules would depend on the respective free energies of crystal polymorph surfaces, measurements of insect knockdown times for the metastable crystal forms were as much as nine times faster acting than the commercial form against Aedes, Anopheles, and Culex mosquitoes as well as Drosophila (fruit flies). These results suggest that replacement of commercially available imidacloprid crystals (a.k.a. Form I) in space-spraying with any one of three new polymorphs, Forms IV, VI, IX, would suppress vector-borne disease transmission while reducing environmental exposure and harm to nontarget organisms.

Progress and challenges in virus genomic epidemiology

Genomic epidemiology, which links pathogen genomes with associated metadata to understand disease transmission, has become a key component of outbreak response. Decreasing costs of genome sequencing and increasing computational power provide opportunities to generate and analyse large viral genomic datasets that aim to uncover the spatial scales of transmission, the demographics contributing to transmission patterns, and to forecast epidemic trends. Emerging sources of genomic data and associated metadata provide new opportunities to further unravel transmission patterns. Key challenges include how to integrate genomic data with metadata from multiple sources, how to generate efficient computational algorithms to cope with large datasets, and how to establish sampling frameworks to enable robust conclusions.

Wave-wise comparative genomic study for revealing the complete scenario and dynamic nature of COVID-19 pandemic in Bangladesh

As the COVID-19 pandemic continues to ravage across the globe and take millions of lives and like many parts of the world, the second wave of the pandemic hit Bangladesh, this study aimed at understanding its causative agent, SARS-CoV-2 at the genomic and proteomic level and provide precious insights about the pathogenesis, evolution, strengths and weaknesses of the virus. As of Mid-June 2021, over 1500 SARS-CoV-2 genomesequences have been deposited in the GISAID database from Bangladesh which were extracted and categorized into two waves. By analyzing these genome sequences, it was discovered that the wave-2 samples had a significantly greater average rate of mutation/sample (30.79%) than the wave-1 samples (12.32%). Wave-2 samples also had a higher frequency of deletion, and transversion events. During the first wave, the GR clade was the most predominant but it was replaced by the GH clade in the latter wave. The B.1.1.25 variant showed the highest frequency in wave-1 while in case of wave-2, the B.1.351.3 variant, was the most common one. A notable presence of the delta variant, which is currently at the center of concern, was also observed. Comparison of the Spike protein found in the reference and the 3 most common lineages found in Bangladesh namely, B.1.1.7, B.1.351, B.1.617 in terms of their ability to form stable complexes with ACE2 receptor revealed that B.1.617 had the potential to be more transmissible than others. Importantly, no indigenous variants have been detected so far which implies that the successful prevention of import of foreign variants can diminish the outbreak in the country.

VGEA: an RNA viral assembly toolkit

Next generation sequencing (NGS)-based studies have vastly increased our understanding of viral diversity. Viral sequence data obtained from NGS experiments are a rich source of information, these data can be used to study their epidemiology, evolution, transmission patterns, and can also inform drug and vaccine design. Viral genomes, however, represent a great challenge to bioinformatics due to their high mutation rate and forming quasispecies in the same infected host, bringing about the need to implement advanced bioinformatics tools to assemble consensus genomes well-representative of the viral population circulating in individual patients. Many tools have been developed to preprocess sequencing reads, carry-out de novo or reference-assisted assembly of viral genomes and assess the quality of the genomes obtained. Most of these tools however exist as standalone workflows and usually require huge computational resources. Here we present (Viral Genomes Easily Analyzed), a Snakemake workflow for analyzing RNA viral genomes. VGEA enables users to map sequencing reads to the human genome to remove human contaminants, split bam files into forward and reverse reads, carry out de novo assembly of forward and reverse reads to generate contigs, pre-process reads for quality and contamination, map reads to a reference tailored to the sample using corrected contigs supplemented by the user's choice of reference sequences and evaluate/compare genome assemblies. We designed a project with the aim of creating a flexible, easy-to-use and all-in-one pipeline from existing/stand-alone bioinformatics tools for viral genome analysis that can be deployed on a personal computer. VGEA was built on the Snakemake workflow management system and utilizes existing tools for each step: fastp (Chen et al., 2018) for read trimming and read-level quality control, BWA (Li & Durbin, 2009) for mapping sequencing reads to the human reference genome, SAMtools (Li et al., 2009) for extracting unmapped reads and also for splitting bam files into fastq files, IVA (Hunt et al., 2015) for de novo assembly to generate contigs, shiver (Wymant et al., 2018) to pre-process reads for quality and contamination, then map to a reference tailored to the sample using corrected contigs supplemented with the user's choice of existing reference sequences, SeqKit (Shen et al., 2016) for cleaning shiver assembly for QUAST, QUAST (Gurevich et al., 2013) to evaluate/assess the quality of genome assemblies and MultiQC (Ewels et al., 2016) for aggregation of the results from fastp, BWA and QUAST. Our pipeline was successfully tested and validated with SARS-CoV-2 (n = 20), HIV-1 (n = 20) and Lassa Virus (n = 20) datasets all of which have been made publicly available. VGEA is freely available on GitHub at: https://github.com/pauloluniyi/VGEA under the GNU General Public License.

Discriminating Clonotypes of Influenza A Virus Genes by Nanopore Sequencing

Influenza viruses still pose a serious threat to humans, and we have not yet been able to effectively predict future pandemic strains and prepare vaccines in advance. One of the main reasons is the high genetic diversity of influenza viruses. We do not know the individual clonotypes of a virus population because some are the majority and others make up only a small fraction of the population. First-generation (FGS) and next-generation sequencing (NGS) technologies have inherent limitations that are unable to resolve a minority clonotype’s information in the virus population. Third-generation sequencing (TGS) technologies with ultra-long reads have the potential to solve this problem but have a high error rate. Here, we evaluated emerging direct RNA sequencing and cDNA sequencing with the MinION platform and established a novel approach that combines the high accuracy of Illumina sequencing technology and long reads of nanopore sequencing technology to resolve both variants and clonotypes of influenza virus. Furthermore, a new program was written to eliminate the effect of nanopore sequencing errors for the analysis of the results. By using this pipeline, we identified 47 clonotypes in our experiment. We conclude that this approach can quickly discriminate the clonotypes of virus genes, allowing researchers to understand virus adaptation and evolution at the population level.

Zika Virus Non-Structural Protein 1 Antigen-Capture Immunoassay

Infection with Zika virus (ZIKV), a member of the Flavivirus genus of the Flaviviridae family, typically results in mild self-limited illness, but severe neurological disease occurs in a limited subset of patients. In contrast, serious outcomes commonly occur in pregnancy that affect the developing fetus, including microcephaly and other major birth defects. The genetic similarity of ZIKV to other widespread flaviviruses, such as dengue virus (DENV), presents a challenge to the development of specific ZIKV diagnostic assays. Nonstructural protein 1 (NS1) is established for use in immunodiagnostic assays for flaviviruses. To address the cross-reactivity of ZIKV NS1 with proteins from other flaviviruses we used site-directed mutagenesis to modify putative epitopes. Goat polyclonal antibodies to variant ZIKV NS1 were affinity-purified to remove antibodies binding to the closely related NS1 protein of DENV. An antigen-capture ELISA configured with the affinity-purified polyclonal antibody showed a linear dynamic range between approximately 500 and 30 ng/mL, with a limit of detection of between 1.95 and 7.8 ng/mL. NS1 proteins from DENV, yellow fever virus, St. Louis encephalitis virus and West Nile virus showed significantly reduced reactivity in the ZIKV antigen-capture ELISA. Refinement of approaches similar to those employed here could lead to development of ZIKV-specific immunoassays suitable for use in areas where infections with related flaviviruses are common.

The impact of long-term non-pharmaceutical interventions on COVID-19 epidemic dynamics and control: the value and limitations of early models

Mathematical models of epidemics are important tools for predicting epidemic dynamics and evaluating interventions. Yet, because early models are built on limited information, it is unclear how long they will accurately capture epidemic dynamics. Using a stochastic SEIR model of COVID-19 fitted to reported deaths, we estimated transmission parameters at different time points during the first wave of the epidemic (March-June, 2020) in Santa Clara County, California. Although our estimated basic reproduction number ([Formula: see text]) remained stable from early April to late June (with an overall median of 3.76), our estimated effective reproduction number ([Formula: see text]) varied from 0.18 to 1.02 in April before stabilizing at 0.64 on 27 May. Between 22 April and 27 May, our model accurately predicted dynamics through June; however, the model did not predict rising summer cases after shelter-in-place orders were relaxed in June, which, in early July, was reflected in cases but not yet in deaths. While models are critical for informing intervention policy early in an epidemic, their performance will be limited as epidemic dynamics evolve. This paper is one of the first to evaluate the accuracy of an early epidemiological compartment model over time to understand the value and limitations of models during unfolding epidemics.

Assessing the impact of COVID-19 border restrictions on dengue transmission in Yunnan Province, China: an observational epidemiological and phylogenetic analysis

BACKGROUND

In response to the COVID-19 pandemic, China implemented strict restrictions on cross-border travel to prevent disease importation. Yunnan, a Chinese province that borders dengue-endemic countries in Southeast Asia, experienced unprecedented reduction in dengue, from 6840 recorded cases in 2019 to 260 in 2020.

METHODS

Using a combination of epidemiological and virus genomic data, collected from 2013 to 2020 in Yunnan and neighbouring countries, we conduct a series of analyses to characterise the role of virus importation in driving dengue dynamics in Yunnan and assess the association between recent international travel restrictions and the decline in dengue reported in Yunnan in 2020.

FINDINGS

We find strong evidence that dengue incidence between 2013-2019 in Yunnan was closely linked with international importation of cases. A 0-2 month lag in incidence not explained by seasonal differences, absence of local transmission in the winter, effective reproductive numbers < 1 (as estimated independently using genetic data) and diverse cosmopolitan dengue virus phylogenies all suggest dengue is non-endemic in Yunnan. Using a multivariate statistical model we show that the substantial decline in dengue incidence observed in Yunnan in 2020 but not in neighbouring countries is closely associated with the timing of international travel restrictions, even after accounting for other environmental drivers of dengue incidence.

INTERPRETATION

We conclude that Yunnan is a regional sink for DENV lineage movement and that border restrictions may have substantially reduced dengue burden in 2020, potentially averting thousands of cases. Targeted testing and surveillance of travelers returning from high-risk areas could help to inform public health strategies to minimise or even eliminate dengue outbreaks in non-endemic settings like southern China.

FUNDING

Funding for this study was provided by National Key Research and Development Program of China, Beijing Science and Technology Planning Project (Z201100005420010); Beijing Natural Science Foundation (JQ18025); Beijing Advanced Innovation Program for Land Surface Science; National Natural Science Foundation of China (82073616); Young Elite Scientist Sponsorship Program by CAST (YESS) (2018QNRC001); H.T., O.P.G. and M.U.G.K. acknowledge support from the Oxford Martin School. O.J.B was supported by a Wellcome Trust Sir Henry Wellcome Fellowship (206471/Z/17/Z). Chinese translation of the abstract (Appendix 2).

The Potential Role of Clinical Metagenomics in Infectious Diseases: Therapeutic Perspectives

Clinical metagenomics (CMg) is the process of sequencing nucleic acid of clinical samples to obtain clinically relevant information such as the identification of microorganisms and their susceptibility to antimicrobials. Over the last decades, sequencing and bioinformatic solutions supporting CMg have much evolved and an increasing number of case reports and series covering various infectious diseases have been published. Metagenomics is a new approach to infectious disease diagnosis that is currently being developed and is certainly one of the most promising for the coming years. However, most CMg studies are retrospective, and few address the potential impact CMg could have on patient management, including initiation, adaptation, or cessation of antimicrobials. In this narrative review, we have discussed the potential role of CMg in bacteriology, virology, mycology, and parasitology. Several reports and case-series confirm that CMg is an innovative tool with which one can (i) identify more microorganisms than with conventional methods in a single test, (ii) obtain results within hours, and (iii) tailor the antimicrobial regimen of patients. However, the cost-efficiency of CMg and its real impact on patient management are still to be determined.