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.</p

Zika virus infection among symptomatic patients from two healthcare centers in Sao Paulo State, Brazil: prevalence, clinical characteristics, viral detection in body fluids and serodynamics.

Zika virus (ZIKV) clinical presentation and frequency/duration of shedding need further clarification. Symptomatic ZIKV-infected individuals identified in two hospitals in Sao Paulo State, Brazil, were investigated regarding clinical characteristics, shedding in body fluids, and serodynamics. Ninety-four of 235 symptomatic patients (Site A: 58%; Site B: 16%) had Real-Time PCR-confirmed ZIKV infection; fever, headache and gastrointestinal symptoms were less frequent, and rash was more frequent compared to ZIKV-negative patients. Real-Time PCR in serum had worse performance compared to plasma, while urine had the highest sensitivity. Shedding in genital fluids and saliva was rare. IgM positivity was the highest 28 days (24%); IgG positivity increased >14 days (96%) remaining positive in 94% of patients >28 days. ZIKV prevalence varied importantly in two neighboring cities during the same transmission season. Urine Real-Time PCR can improve diagnostic sensitivity; serum testing is less useful. Accurate serological tests are needed to improve diagnosis and surveillance

Zika virus infection among symptomatic patients from two healthcare centers in Sao Paulo State, Brazil: prevalence, clinical characteristics, viral detection in body fluids and serodynamics

Zika virus (ZIKV) clinical presentation and frequency/duration of shedding need further clarification. Symptomatic ZIKV-infected individuals identified in two hospitals in Sao Paulo State, Brazil, were investigated regarding clinical characteristics, shedding in body fluids, and serodynamics. Ninety-four of 235 symptomatic patients (Site A: 58%; Site B: 16%) had Real-Time PCR-confirmed ZIKV infection; fever, headache and gastrointestinal symptoms were less frequent, and rash was more frequent compared to ZIKV-negative patients. Real-Time PCR in serum had worse performance compared to plasma, while urine had the highest sensitivity. Shedding in genital fluids and saliva was rare. IgM positivity was the highest &lt;14 days after the symptoms onset (86%), decreasing &gt;28 days (24%); IgG positivity increased &gt;14 days (96%) remaining positive in 94% of patients &gt;28 days. ZIKV prevalence varied importantly in two neighboring cities during the same transmission season. Urine Real-Time PCR can improve diagnostic sensitivity; serum testing is less useful. Accurate serological tests are needed to improve diagnosis and surveillance

Multiplex qPCR Discriminates Variants of Concern to Enhance Global Surveillance of SARS-CoV-2

With the emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants that may increase transmissibility and/or cause escape from immune responses, there is an urgent need for the targeted surveillance of circulating lineages. It was found that the B.1.1.7 (also 501Y.V1) variant, first detected in the United Kingdom, could be serendipitously detected by the Thermo Fisher TaqPath COVID-19 PCR assay because a key deletion in these viruses, spike Δ69-70, would cause a spike gene target failure (SGTF) result. However, a SGTF result is not definitive for B.1.1.7, and this assay cannot detect other variants of concern (VOC) that lack spike Δ69-70, such as B.1.351 (also 501Y.V2), detected in South Africa, and P.1 (also 501Y.V3), recently detected in Brazil. We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all 3 variants, which has not yet been widely detected in other SARS-CoV-2 lineages. Using ORF1a Δ3675-3677 as the primary target and spike Δ69-70 to differentiate, we designed and validated an open-source PCR assay to detect SARS-CoV-2 VOC. Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence and spread of B.1.1.7, B.1.351, and P.1

Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness

Multiplex PCR method for MinION and Illumina sequencing of Zika and other virus genomes directly from clinical samples

Genome sequencing has become a powerful tool for studying emerging infectious diseases; however, genome sequencing directly from clinical samples (i.e., without isolation and culture) remains challenging for viruses such as Zika, for which metagenomic sequencing methods may generate insufficient numbers of viral reads. Here we present a protocol for generating coding-sequence-complete genomes, comprising an online primer design tool, a novel multiplex PCR enrichment protocol, optimized library preparation methods for the portable MinION sequencer (Oxford Nanopore Technologies) and the Illumina range of instruments, and a bioinformatics pipeline for generating consensus sequences. The MinION protocol does not require an Internet connection for analysis, making it suitable for field applications with limited connectivity. Our method relies on multiplex PCR for targeted enrichment of viral genomes from samples containing as few as 50 genome copies per reaction. Viral consensus sequences can be achieved in 1-2 d by starting with clinical samples and following a simple laboratory workflow. This method has been successfully used by several groups studying Zika virus evolution and is facilitating an understanding of the spread of the virus in the Americas. The protocol can be used to sequence other viral genomes using the online Primal Scheme primer designer software. It is suitable for sequencing either RNA or DNA viruses in the field during outbreaks or as an inexpensive, convenient method for use in the lab

Genomic Surveillance of Yellow Fever Virus Epizootic in São Paulo, Brazil, 2016 – 2018

São Paulo, a densely inhabited state in southeast Brazil that contains the fourth most populated city in the world, recently experienced its largest yellow fever virus (YFV) outbreak in decades. YFV does not normally circulate extensively in São Paulo, so most people were unvaccinated when the outbreak began. Surveillance in non-human primates (NHPs) is important for determining the magnitude and geographic extent of an epizootic, thereby helping to evaluate the risk of YFV spillover to humans. Data from infected NHPs can give more accurate insights into YFV spread than when using data from human cases alone. To contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo, we generated and analysed virus genomic data and epizootic case data from NHPs in São Paulo. We report the occurrence of three spatiotemporally distinct phases of the outbreak in São Paulo prior to February 2018. We generated 51 new virus genomes from YFV positive cases identified in 23 different municipalities in São Paulo, mostly sampled from NHPs between October 2016 and January 2018. Although we observe substantial heterogeneity in lineage dispersal velocities between phylogenetic branches, continuous phylogeographic analyses of generated YFV genomes suggest that YFV lineages spread in São Paulo at a mean rate of approximately 1km per day during all phases of the outbreak. Viral lineages from the first epizootic phase in northern São Paulo subsequently dispersed towards the south of the state to cause the second and third epizootic phases there. This alters our understanding of how YFV was introduced into the densely populated south of São Paulo state. Our results shed light on the sylvatic transmission of YFV in highly fragmented forested regions in São Paulo state and highlight the importance of continued surveillance of zoonotic pathogens in sentinel species

Phylogenetics, Epidemiology and Temporal Patterns of Dengue Virus in Araraquara, São Paulo State

Dengue virus (DENV) is a prominent arbovirus with global spread, causing approximately 390 million infections each year. In Brazil, yearly epidemics follow a well-documented pattern of serotype replacement every three to four years on average. Araraquara, located in the state of São Paulo, has faced significant impacts from DENV epidemics since the emergence of DENV-1 in 2010. The municipality then transitioned from low to moderate endemicity in less than 10 years. Yet, there remains an insufficient understanding of virus circulation dynamics, particularly concerning DENV-1, in the region, as well as the genetic characteristics of the virus. To address this, we sequenced 37 complete or partial DENV-1 genomes sampled from 2015 to 2022 in Araraquara. Then, using also Brazilian and worldwide DENV-1 sequences we reconstructed the evolutionary history of DENV-1 in Araraquara and estimated the time to the most recent common ancestor (tMRCA) for serotype 1, for genotype V and its main lineages. Within the last ten years, there have been at least three introductions of genotype V in Araraquara, distributed in two main lineages (L Ia and L Ib, and L II). The tMRCA for the first sampled lineage (2015/2016 epidemics) was approximately 15 years ago (in 2008). Crucially, our analysis challenges existing assumptions regarding the emergence time of the DENV-1 genotypes, suggesting that genotype V might have diverged more recently than previously described. The presence of the two lineages of genotype V in the municipality might have contributed to the extended persistence of DENV-1 in the region

Genomic Surveillance of Yellow Fever Virus Epizootic in São Paulo, Brazil, 2016 – 2018

São Paulo, a densely inhabited state in southeast Brazil that contains the fourth most populated city in the world, recently experienced its largest yellow fever virus (YFV) outbreak in decades. YFV does not normally circulate extensively in São Paulo, so most people were unvaccinated when the outbreak began. Surveillance in non-human primates (NHPs) is important for determining the magnitude and geographic extent of an epizootic, thereby helping to evaluate the risk of YFV spillover to humans. Data from infected NHPs can give more accurate insights into YFV spread than when using data from human cases alone. To contextualise human cases, identify epizootic foci and uncover the rate and direction of YFV spread in São Paulo, we generated and analysed virus genomic data and epizootic case data from NHPs in São Paulo. We report the occurrence of three spatiotemporally distinct phases of the outbreak in São Paulo prior to February 2018. We generated 51 new virus genomes from YFV positive cases identified in 23 different municipalities in São Paulo, mostly sampled from NHPs between October 2016 and January 2018. Although we observe substantial heterogeneity in lineage dispersal velocities between phylogenetic branches, continuous phylogeographic analyses of generated YFV genomes suggest that YFV lineages spread in São Paulo at a mean rate of approximately 1km per day during all phases of the outbreak. Viral lineages from the first epizootic phase in northern São Paulo subsequently dispersed towards the south of the state to cause the second and third epizootic phases there. This alters our understanding of how YFV was introduced into the densely populated south of São Paulo state. Our results shed light on the sylvatic transmission of YFV in highly fragmented forested regions in São Paulo state and highlight the importance of continued surveillance of zoonotic pathogens in sentinel species.ISSN:1553-7374ISSN:1553-736