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.

Lessons from a Multilaboratorial Task Force for Diagnosis of a Fatal Toxoplasmosis Outbreak in Captive Primates in Brazil.

Toxoplasmosis is an important zoonotic disease caused by the parasite Toxoplasma gondii and is especially fatal for neotropical primates. In Brazil, the Ministry of Health is responsible for national epizootic surveillance, but some diseases are still neglected. Here, we present an integrated investigation of an outbreak that occurred during the first year of the COVID-19 pandemic among eleven neotropical primates housed at a primatology center in Brazil. After presenting non-specific clinical signs, all animals died within four days. A wide range of pathogens were evaluated, and we successfully identified T. gondii as the causative agent within four days after necropsies. The liver was the most affected organ, presenting hemorrhage and hepatocellular necrosis. Tachyzoites and bradyzoite cysts were observed in histological examinations and immunohistochemistry in different organs; in addition, parasitic DNA was detected through PCR in blood samples from all specimens evaluated. A high prevalence of Escherichia coli was also observed, indicating sepsis. This case highlights some of the obstacles faced by the current Brazilian surveillance system. A diagnosis was obtained through the integrated action of researchers since investigation for toxoplasmosis is currently absent in national guidelines. An interdisciplinary investigation could be a possible model for future epizootic investigations in animals.

Genetic differences of dengue virus 2 in patients with distinct clinical outcome.

The genetic diversity of the dengue virus is characterized by four circulating serotypes, several genotypes, and an increasing number of existing lineages that may have differences in the potential to cause epidemics and disease severity. Accurate identification of the genetic variability of the virus is essential to identify lineages responsible for an epidemic and understanding the processes of virus spread and virulence. Here, we characterize, using portable nanopore genomic sequencing, different lineages of dengue virus 2 (DENV-2) detected in 22 serum samples from patients with and without dengue warning signs attended at Hospital de Base of São José do Rio Preto (SJRP) in 2019, during a DENV-2 outbreak. Demographic, epidemiological, and clinical data were also analyzed. The phylogenetic reconstruction and the clinical data showed that two lineages belonging to the American/Asian genotype of DENV-2-BR3 and BR4 (BR4L1 and BR4L2)-were co-circulating in SJRP. Although preliminary, these results indicate no specific association between clinical form and phylogenetic clustering at the virus consensus sequence level. Studies with larger sample sizes and which explore single nucleotide variants are needed. Therefore, we showed that portable nanopore genome sequencing could generate quick and reliable sequences for genomic surveillance to monitor viral diversity and its association with disease severity as an epidemic unfolds.

Dynamics of Early Establishment of SARS-CoV-2 VOC Omicron Lineages in Minas Gerais, Brazil.

Brazil is one of the nations most affected by Coronavirus disease 2019 (COVID-19). The introduction and establishment of new virus variants can be related to an increase in cases and fatalities. The emergence of Omicron, the most modified SARS-CoV-2 variant, caused alarm for the public health of Brazil. In this study, we examined the effects of the Omicron introduction in Minas Gerais (MG), the second-most populous state of Brazil. A total of 430 Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) samples from November 2021 to June 2022 from Belo Horizonte (BH) city were sequenced. These newly sequenced genomes comprise 72% of all previously available SARS-CoV-2 genomes for the city. Evolutionary analysis of novel viral genomes reveals that a great diversity of Omicron sublineages have circulated in BH, a pattern in-keeping with observations across Brazil more generally. Bayesian phylogeographic reconstructions indicate that this diversity is a product of a large number of international and national importations. As observed previously, São Paulo state is shown as a significant hub for viral spread throughout the country, contributing to around 70% of all viral Omicron introductions detected in MG.

Molecular characterization and sequecing analysis of SARS-CoV-2 genome in Minas Gerais, Brazil.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), first identified in Wuhan, China, is the causative agent of the coronavirus disease 2019 (COVID-19). Since its first notification in São Paulo state (SP) on 26th February 2020, more than 22,300,000 cases and 619,000 deaths were reported in Brazil. In early pandemic, SARS-CoV-2 spread locally, however, over time, this virus was disseminated to other regions of the country. Herein, we performed genomic sequencing and phylogenetic analysis of SARS-CoV-2 using 20 clinical samples of COVID-19 confirmed cases from 9 cities of Minas Gerais state (MG), in order to evaluate the molecular properties of circulating viral strains in this locality from March to May 2020. Our analyses demonstrated the circulation of B.1 lineage isolates in the investigated locations and nucleotide substitutions were observed into the genomic regions related to important viral structures. Additionally, sequences generated in this study clustered with isolates from SP, suggesting a dissemination route between these two states. Alternatively, monophyletic groups of sequences from MG and other states or country were observed, indicating independent events of virus introduction. These results reinforce the need of genomic surveillance for understand the ongoing spread of emerging viral pathogens.

Retrospective Genomic Surveillance of Chikungunya Transmission in Minas Gerais State, Southeast Brazil.

Brazil accounted for a total number of 1,276,194 reported cases of chikungunya fever between 2014 and 2022. Additionally, since 2015, the country has experienced an increasing death toll, in which the Northeast and Southeast regions appear to report the worst scenarios. Although the CHIKV transmission dynamics have been studied in many parts of the country since its introduction in 2014, little is still known about chikungunya virus (CHIKV) transmission and genetic diversity in the state of Minas Gerais, located in southeast Brazil. Moreover, no studies have been published characterizing CHIKV genomic surveillance in this state. Thus, to retrospectively explore the CHIKV epidemic in Minas Gerais, we generated 40 genomes from clinical samples using Nanopore sequencing. Phylogenetic analysis indicated that multiple introductions of CHIKV occurred, likely from the northeastern Brazilian states, with the most recent common ancestral strain dating to early March 2016, which is in agreement with local epidemiological reports. Additionally, epidemiological data reveals a decline in the number of reported cases from 2017 to 2021, indicating that population immunity or changes in vector activity may have contributed to the decreasing waves of CHIKV infection. Together, our results shed light on the dispersion dynamics of CHIKV and show that infections decreased from March 2017 to January 2021 despite multiple introductions into Minas Gerais State. In conclusion, our study highlights the importance of combining genomic and epidemiological data in order to assist public health laboratories in monitoring and understanding the patterns and diversity of mosquito-borne viral epidemics. IMPORTANCE Arbovirus infections in Brazil, including chikungunya, dengue, yellow fever, and Zika, result in considerable morbidity and mortality and are pressing public health concerns. However, our understanding of these outbreaks is hampered by the limited availability of genomic data. In this study, we combine epidemiological analysis and portable genome sequencing to retrospectively describe the CHIKV epidemic in Minas Gerais between 2017 and 2021. Our results indicate that the East/Central/South African (ECSA) CHIKV lineage was introduced into Minas Gerais by three distinct events, likely from the North and Northeast regions of Brazil. Our study provides an understanding of how CHIKV initiates transmission in the region and illustrates that genomics in the field can augment traditional approaches to infectious disease surveillance and control.

Epidemiologic History and Genetic Diversity Origins of Chikungunya and Dengue Viruses, Paraguay.

Paraguay has been severely affected by emergent Zika and chikungunya viruses, and dengue virus is endemic. To learn more about the origins of genetic diversity and epidemiologic history of these viruses in Paraguay, we deployed portable sequencing technologies to strengthen genomic surveillance and determine the evolutionary and epidemic history of arthropod-borne viruses (arboviruses). Samples stored at the Paraguay National Central Laboratory were sequenced and subjected to phylogenetic analysis. Among 33 virus genomes generated, we identified 2 genotypes of chikungunya and 2 serotypes of dengue virus that circulated in Paraguay during 2014-2018; the main source of these virus lineages was estimated to be Brazil. The evolutionary history inferred by our analyses precisely matched the available travel history of the patients. The genomic surveillance approach used was valuable for describing the epidemiologic history of arboviruses and can be used to determine the origins and evolution of future arbovirus outbreaks.

PCR assay to enhance global surveillance for SARS-CoV-2 variants of concern.

With the emergence of SARS-CoV-2 variants that may increase transmissibility and/or cause escape from immune responses 1-3 , 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 UK 4,5 could be serendipitously detected by the ThermoFisher 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 that lack spike Δ69-70, such as B.1.351 (also 501Y.V2) detected in South Africa 6 and P.1 (also 501Y.V3) recently detected in Brazil 7 . We identified a deletion in the ORF1a gene (ORF1a Δ3675-3677) in all three 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 variants of concern 8 . Our assay can be rapidly deployed in laboratories around the world to enhance surveillance for the local emergence spread of B.1.1.7, B.1.351, and P.1.

Local Transmission of SARS-CoV-2 Lineage B.1.1.7, Brazil, December 2020.

In December 2020, research surveillance detected the B.1.1.7 lineage of severe acute respiratory syndrome coronavirus 2 in São Paulo, Brazil. Rapid genomic sequencing and phylogenetic analysis revealed 2 distinct introductions of the lineage. One patient reported no international travel. There may be more infections with this lineage in Brazil than reported.

First case of SARS-COV-2 sequencing in cerebrospinal fluid of a patient with suspected demyelinating disease.

The association between coronaviruses and central nervous system (CNS) demyelinating lesions has been previously shown. However, no case has been described of an association between the novel coronavirus (SARS-COV-2) and CNS demyelinating disease so far. SARS-COV-2 was previously detected in cerebrospinal fluid (CSF) sample of a patient with encephalitis. However, the virus identity was not confirmed by deep sequencing of SARS-COV-2 detected in the CSF. Here, we report a case of a patient with mild respiratory symptoms and neurological manifestations compatible with clinically isolated syndrome. The viral genome of SARS-COV-2 was detected and sequenced in CSF with 99.74-100% similarity between the patient virus and worldwide sequences. This report suggests a possible association of SARS-COV-2 infection with neurological symptoms of demyelinating disease, even in the absence of relevant upper respiratory tract infection signs.

Genomic detection of a virus lineage replacement event of dengue virus serotype 2 in Brazil, 2019.

BACKGROUND Despite efforts to mitigate the impact of dengue virus (DENV) epidemics, the virus remains a public health problem in tropical and subtropical regions around the world. Most DENV cases in the Americas between January and July 2019 were reported in Brazil. São Paulo State in the southeast of Brazil has reported nearly half of all DENV infections in the country. OBJECTIVES To understand the origin and dynamics of the 2019 DENV outbreak. METHODS Here using portable nanopore sequencing we generated20 new DENV genome sequences from viremic patients with suspected dengue infection residing in two of the most-affected municipalities of São Paulo State, Araraquara and São José do Rio Preto. We conducted a comprehensive phylogenetic analysis with 1,630 global DENV strains to better understand the evolutionary history of the DENV lineages that currently circulate in the region. FINDINGS The new outbreak strains were classified as DENV2 genotype III (American/Asian genotype). Our analysis shows that the 2019 outbreak is the result of a novel DENV lineage that was recently introduced to Brazil from the Caribbean region. Dating phylogeographic analysis suggests that DENV2-III BR-4 was introduced to Brazil in or around early 2014, possibly from the Caribbean region. MAIN CONCLUSIONS Our study describes the early detection of a newly introduced and rapidly-expanding DENV2 virus lineage in Brazil.

Author Correction: Nanopore native RNA sequencing of a human poly(A) transcriptome.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Genomic detection of a virus lineage replacement event of dengue virus serotype 2 in Brazil, 2019

BACKGROUND Despite efforts to mitigate the impact of dengue virus (DENV) epidemics, the virus remains a public health problem in tropical and subtropical regions around the world. Most DENV cases in the Americas between January and July 2019 were reported in Brazil. São Paulo State in the southeast of Brazil has reported nearly half of all DENV infections in the country. OBJECTIVES To understand the origin and dynamics of the 2019 DENV outbreak. METHODS Here using portable nanopore sequencing we generated20 new DENV genome sequences from viremic patients with suspected dengue infection residing in two of the most-affected municipalities of São Paulo State, Araraquara and São José do Rio Preto. We conducted a comprehensive phylogenetic analysis with 1,630 global DENV strains to better understand the evolutionary history of the DENV lineages that currently circulate in the region. FINDINGS The new outbreak strains were classified as DENV2 genotype III (American/Asian genotype). Our analysis shows that the 2019 outbreak is the result of a novel DENV lineage that was recently introduced to Brazil from the Caribbean region. Dating phylogeographic analysis suggests that DENV2-III BR-4 was introduced to Brazil in or around early 2014, possibly from the Caribbean region. MAIN CONCLUSIONS Our study describes the early detection of a newly introduced and rapidly-expanding DENV2 virus lineage in Brazil.

Nanopore native RNA sequencing of a human poly(A) transcriptome.

High-throughput complementary DNA sequencing technologies have advanced our understanding of transcriptome complexity and regulation. However, these methods lose information contained in biological RNA because the copied reads are often short and modifications are not retained. We address these limitations using a native poly(A) RNA sequencing strategy developed by Oxford Nanopore Technologies. Our study generated 9.9 million aligned sequence reads for the human cell line GM12878, using thirty MinION flow cells at six institutions. These native RNA reads had a median length of 771 bases, and a maximum aligned length of over 21,000 bases. Mitochondrial poly(A) reads provided an internal measure of read-length quality. We combined these long nanopore reads with higher accuracy short-reads and annotated GM12878 promoter regions to identify 33,984 plausible RNA isoforms. We describe strategies for assessing 3' poly(A) tail length, base modifications and transcript haplotypes.

Yellow Fever Virus Reemergence and Spread in Southeast Brazil, 2016-2019.

The recent reemergence of yellow fever virus (YFV) in Brazil has raised serious concerns due to the rapid dissemination of the virus in the southeastern region. To better understand YFV genetic diversity and dynamics during the recent outbreak in southeastern Brazil, we generated 18 complete and nearly complete genomes from the peak of the epidemic curve from nonhuman primates (NHPs) and human infected cases across the Espírito Santo and Rio de Janeiro states. Genomic sequencing of 18 YFV genomes revealed the estimated timing, source, and likely routes of yellow fever virus transmission and dispersion during one of the largest outbreaks ever registered in Brazil. We showed that during the recent epidemic, YFV was reintroduced from Minas Gerais to the Espírito Santo and Rio de Janeiro states multiple times between 2016 and 2019. The analysis of data from portable sequencing could identify the corridor of spread of YFV. These findings reinforce the idea that continued genomic surveillance strategies can provide information on virus genetic diversity and transmission dynamics that might assist in understanding arbovirus epidemics.IMPORTANCE Arbovirus infections in Brazil, including yellow fever, dengue, zika, and chikungunya, result in considerable morbidity and mortality and are pressing public health concerns. However, our understanding of these outbreaks is hampered by the limited availability of genomic data. In this study, we investigated the genetic diversity and spatial distribution of YFV during the current outbreak by analyzing genomic data from areas in southeastern Brazil not covered by other previous studies. To gain insights into the routes of YFV introduction and dispersion, we tracked the virus by sequencing YFV genomes sampled from nonhuman primates and infected patients from the southeastern region. Our study provides an understanding of how YFV initiates transmission in new Brazilian regions and illustrates that genomics in the field can augment traditional approaches to infectious disease surveillance and control.

Acute Vector-Borne Viral Infection: Zika and MinION Surveillance.

The MinION sequencer was launched by the Oxford Nanopore Technologies start-up as a disruptive technology for genome sequencing based on single-molecule synthesis. Its characteristics as a portable device, low cost, and simple library preparation have made it a good candidate for field researchers. MinION has been used to sequence a number of microorganisms, such as bacteria, viruses, and fungi. Based on the experience that characterized the Ebola virus genetic diversity in Guinea during the 2014-2015 outbreak, the ZiBRA (Zika in Brazil Real-time Analysis) project aimed to sequence a large number of Zika virus genomes during a mobile laboratory trip in northeast Brazil to provide important epidemiological information about the spread of this disease in this country. In response to the positive and rapid results obtained by the ZiBRA project, the Brazilian Ministry of Health and many leading institutions, such as the Pan American Health Organization and WHO, have shown interest in expanding the strategy used in this project to other countries dealing with arbovirus infection. *This article is part of a curated collection.

Lipidomic Analysis Reveals Serum Alteration of Plasmalogens in Patients Infected With ZIKA Virus.

Zika virus (ZIKV) is an arthropod-borne virus (arbovirus) in the Flavivirus genus of the Flaviviridae family. Since the large outbreaks in French Polynesia in 2013-2014 and in Brazil in 2015, ZIKV has been considered a new public health threat. Similar to other related flavivirus, ZIKV is associated with mild and self-limiting symptoms such as rash, pruritus, prostration, headache, arthralgia, myalgia, conjunctivitis, lower back pain and, when present, a short-term low grade fever. In addition, ZIKV has been implicated in neurological complications such as neonatal microcephaly and Guillain-Barré syndrome in adults. Herein, serum lipidomic analysis was used to identify possible alterations in lipid metabolism triggered by ZIKV infection. Patients who presented virus-like symptoms such as fever, arthralgia, headache, exanthema, myalgia and pruritus were selected as the control group. Our study reveals increased levels of several phosphatidylethanolamine (PE) lipid species in the serum of ZIKV patients, the majority of them plasmenyl-phosphatidylethanolamine (pPE) (or plasmalogens) linked to polyunsaturated fatty acids. Constituting up to 20% of total phospholipids in humans, plasmalogens linked to polyunsaturated fatty acids are particularly enriched in neural membranes of the brain. The biosynthesis of plasmalogens requires functional peroxisomes, which are important sites for viral replication, including ZIKV. Thus, increased levels of plasmalogens in serum of ZIKV infected subjects suggest a link between ZIKV life cycle and peroxisomes. Our data provide important insights into specific host cellular lipids that are likely associated with ZIKV replication and may serve as platform for antiviral strategy against ZIKV.

Genomic, epidemiological and digital surveillance of Chikungunya virus in the Brazilian Amazon.

BACKGROUND: Since its first detection in the Caribbean in late 2013, chikungunya virus (CHIKV) has affected 51 countries in the Americas. The CHIKV epidemic in the Americas was caused by the CHIKV-Asian genotype. In August 2014, local transmission of the CHIKV-Asian genotype was detected in the Brazilian Amazon region. However, a distinct lineage, the CHIKV-East-Central-South-America (ECSA)-genotype, was detected nearly simultaneously in Feira de Santana, Bahia state, northeast Brazil. The genomic diversity and the dynamics of CHIKV in the Brazilian Amazon region remains poorly understood despite its importance to better understand the epidemiological spread and public health impact of CHIKV in the country. METHODOLOGY/PRINCIPAL FINDINGS: We report a large CHIKV outbreak (5,928 notified cases between August 2014 and August 2018) in Boa vista municipality, capital city of Roraima's state, located in the Brazilian Amazon region. We generated 20 novel CHIKV-ECSA genomes from the Brazilian Amazon region using MinION portable genome sequencing. Phylogenetic analyses revealed that despite an early introduction of the Asian genotype in 2015 in Roraima, the large CHIKV outbreak in 2017 in Boa Vista was caused by an ECSA-lineage most likely introduced from northeastern Brazil. Epidemiological analyses suggest a basic reproductive number of R0 of 1.66, which translates in an estimated 39 (95% CI: 36 to 45) % of Roraima's population infected with CHIKV-ECSA. Finally, we find a strong association between Google search activity and the local laboratory-confirmed CHIKV cases in Roraima. CONCLUSIONS/SIGNIFICANCE: This study highlights the potential of combining traditional surveillance with portable genome sequencing technologies and digital epidemiology to inform public health surveillance in the Amazon region. Our data reveal a large CHIKV-ECSA outbreak in Boa Vista, limited potential for future CHIKV outbreaks, and indicate a replacement of the Asian genotype by the ECSA genotype in the Amazon region.

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar.

How viruses evolve within hosts can dictate infection outcomes; however, reconstructing this process is challenging. We evaluate our multiplexed amplicon approach, PrimalSeq, to demonstrate how virus concentration, sequencing coverage, primer mismatches, and replicates influence the accuracy of measuring intrahost virus diversity. We develop an experimental protocol and computational tool, iVar, for using PrimalSeq to measure virus diversity using Illumina and compare the results to Oxford Nanopore sequencing. We demonstrate the utility of PrimalSeq by measuring Zika and West Nile virus diversity from varied sample types and show that the accumulation of genetic diversity is influenced by experimental and biological systems.

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.