Spatiotemporal transmission of SARS-CoV-2 lineages during 2020-2021 in Pernambuco-Brazil.

SARS-CoV-2 virus emerged as a new threat to humans and spread around the world, leaving a large death toll. As of January 2023, Brazil is among the countries with the highest number of registered deaths. Nonpharmacological and pharmacological interventions have been heterogeneously implemented in the country, which, associated with large socioeconomic differences between the country regions, has led to distinct virus spread dynamics. Here, we investigate the spatiotemporal dispersion of SARS-CoV-2 lineages in the Pernambuco state (Northeast Brazil) throughout the distinct epidemiological scenarios that unfolded in the first 2 years of the pandemic. We generated a total of 1,389 new SARS-CoV-2 genomes from June 2020 to August 2021. This sampling captured the arrival, communitary transmission, and the circulation of the B1.1, B.1.1.28, and B.1.1.33 lineages; the emergence of the former variant of interest P.2; and the emergence and fast replacement of all previous variants by the more transmissible variant of concern P.1 (Gamma). Based on the incidence and lineage spread pattern, we observed an East-to-West to inner state pattern of transmission, which is in agreement with the transmission of more populous metropolitan areas to medium- and small-size country-side cities in the state. Such transmission patterns may be partially explained by the main routes of traffic across municipalities in the state. Our results highlight that the fine-grained intrastate analysis of lineages and incidence spread can provide actionable insights for planning future nonpharmacological intervention for air-borne transmissible human pathogens.IMPORTANCEDuring the COVID-19 pandemic, Brazil was one of the most affected countries, mainly due its continental-size, socioeconomic differences among regions, and heterogeneous implementation of intervention methods. In order to investigate SARS-CoV-2 dynamics in the state of Pernambuco, we conducted a spatiotemporal dispersion study, covering the period from June 2020 to August 2021, to comprehend the dynamics of viral transmission during the first 2 years of the pandemic. Throughout this study, we were able to track three significant epidemiological waves of transmission caused by B1.1, B.1.1.28, B.1.1.33, P.2, and P.1 lineages. These analyses provided valuable insights into the evolution of the epidemiological landscape, contributing to a deeper understanding of the dynamics of virus transmission during the early years of the pandemic in the state of Pernambuco.

Evidence of Zika Virus Reinfection by Genome Diversity and Antibody Response Analysis, Brazil.

We generated 238 Zika virus (ZIKV) genomes from 135 persons in Brazil who had samples collected over 1 year to evaluate virus persistence. Phylogenetic inference clustered the genomes together with previously reported ZIKV strains from northern Brazil, showing that ZIKV has been remained relatively stable over time. Temporal phylogenetic analysis revealed limited within-host diversity among most ZIKV-persistent infected associated samples. However, we detected unusual virus temporal diversity from >5 persons, uncovering the existence of divergent genomes within the same patient. All those patients showed an increase in neutralizing antibody levels, followed by a decline at the convalescent phase of ZIKV infection. Of interest, in 3 of those patients, titers of neutralizing antibodies increased again after 6 months of ZIKV infection, concomitantly with real-time reverse transcription PCR re-positivity, supporting ZIKV reinfection events. Altogether, our findings provide evidence for the existence of ZIKV reinfection events.

Chikungunya virus infection in the southernmost state of Brazil was characterised by self-limited transmission (2017-2019) and a larger 2021 outbreak.

BACKGROUND: Chikungunya is a mosquito-borne virus that has been causing large outbreaks in the Americas since 2014. In Brazil, Asian-Caribbean (AC) and East-Central-South-African (ECSA) genotypes have been detected and lead to large outbreaks in several Brazilian states. In Rio Grande do Sul (RS), the southernmost state of Brazil, the first cases were reported in 2016. OBJECTIVES AND METHODS: We employed genome sequencing and epidemiological investigation to characterise the Chikungunya fever (CHIKF) burden in RS between 2017-2021. FINDINGS: We detected an increasing CHIKF burden linked to travel associated introductions and communitary transmission of distinct lineages of the ECSA genotype during this period. MAIN CONCLUSIONS: Until 2020, CHIKV introductions were most travel associated and transmission was limited. Then, in 2021, the largest outbreak occurred in the state associated with the introduction of a new ECSA lineage. CHIKV outbreaks are likely to occur in the near future due to abundant competent vectors and a susceptible population, exposing more than 11 million inhabitants to an increasing infection risk.

Uncovering neglected subtypes and zoonotic transmission of Hepatitis E virus (HEV) in Brazil.

Hepatitis E virus (HEV) circulation in humans and swine has been extensively studied in South America over the last two decades. Nevertheless, only 2.1% of reported HEV strains are available as complete genome sequences. Therefore, many clinical, epidemiological, and evolutionary aspects of circulating HEV in the continent still need to be clarified. Here, we conducted a retrospective evolutionary analysis of one human case and six swine HEV strains previously reported in northeastern, southern, and southeastern Brazil. We obtained two complete and four nearly complete genomic sequences. Evolutionary analysis comparing the whole genomic and capsid gene sequences revealed high genetic variability. This included the circulation of at least one unrecognized unique South American subtype. Our results corroborate that sequencing the whole capsid gene could be used as an alternative for HEV subtype assignment in the absence of complete genomic sequences. Moreover, our results substantiate the evidence for zoonotic transmission by comparing a larger genomic fragment recovered from the sample of the autochthonous human hepatitis E case. Further studies should continuously investigate HEV genetic diversity and zoonotic transmission of HEV in South America.

Chikungunya virus infection in the southernmost state of Brazil was characterised by self-limited transmission (2017-2019) and a larger 2021 outbreak

BACKGROUND Chikungunya is a mosquito-borne virus that has been causing large outbreaks in the Americas since 2014. In Brazil, Asian-Caribbean (AC) and East-Central-South-African (ECSA) genotypes have been detected and lead to large outbreaks in several Brazilian states. In Rio Grande do Sul (RS), the southernmost state of Brazil, the first cases were reported in 2016. OBJECTIVES AND METHODS We employed genome sequencing and epidemiological investigation to characterise the Chikungunya fever (CHIKF) burden in RS between 2017-2021. FINDINGS We detected an increasing CHIKF burden linked to travel associated introductions and communitary transmission of distinct lineages of the ECSA genotype during this period. MAIN CONCLUSIONS Until 2020, CHIKV introductions were most travel associated and transmission was limited. Then, in 2021, the largest outbreak occurred in the state associated with the introduction of a new ECSA lineage. CHIKV outbreaks are likely to occur in the near future due to abundant competent vectors and a susceptible population, exposing more than 11 million inhabitants to an increasing infection risk.

Multiple introductions and country-wide spread of DENV-2 genotype II (Cosmopolitan) in Brazil.

Dengue virus serotype 2, genotype Cosmopolitan (DENV-2-GII), is one of the most widespread DENV strains globally. In the USA, DENV-2 epidemics have been dominated by DENV-2 genotype Asian-American (DENV-2-GIII), and the first cases of DENV-2-GII were only described in 2019, in Peru, and in 2021 in Brazil. To gain new information about the circulation of DENV-2-GII in Brazil, we sequenced 237 DENV-2 confirmed cases sampled between March 2021 and March 2023 and revealed that DENV-2-GII is already present in all geographic regions of Brazil. The phylogeographic analysis inferred that DENV-2-GII was introduced at least four times in Brazil, between May 2020 and August 2022, generating multiple clades that spread throughout the country with different success. Despite multiple introductions of DENV-2-GII, analysis of the country-wide laboratory surveillance data showed that the Brazilian dengue epidemic in 2022 was dominated by DENV-1 in most states. We hypothesize that massive circulation of DENV-2-GIII in previous years in Brazil might have created a population immune barrier against symptomatic homotypic reinfections by DENV-2-GII, leading to sustained cryptic circulation in asymptomatic cases and localized outbreaks of this new genotype. In summary, our study stresses the importance of arboviral genomic surveillance to close monitoring and better understanding the potential impact of DENV-2-GII in the coming years.

Identification of a Virulent Newcastle Disease Virus Strain Isolated from Pigeons (Columbia livia) in Northeastern Brazil Using Next-Generation Genome Sequencing.

Newcastle disease virus (NDV), also known as avian paramyxoviruses 1 (APMV-1) is among the most important viruses infecting avian species. Given its widespread circulation, there is a high risk for the reintroduction of virulent strains into the domestic poultry industry, making the surveillance of wild and domestic birds a crucial process to appropriately respond to novel outbreaks. In the present study, we investigated an outbreak characterized by the identification of sick pigeons in a large municipality in Northeastern Brazil in 2018. The affected pigeons presented neurological signs, including motor incoordination, torticollis, and lethargy. Moribund birds were collected, and through a detailed histopathological analysis we identified severe lymphoplasmacytic meningoencephalitis with perivascular cuffs and gliosis in the central nervous system, and lymphoplasmacytic inflammation in the liver, kidney, and intestine. A total of five pigeons tested positive for NDV, as assessed by rRT-PCR targeted to the M gene. Laboratory virus isolation on Vero E6 cells confirmed infection, after the recovery of infectious NVD from brain and kidney tissues. We next characterized the isolated NDV/pigeon/PE-Brazil/MP003/2018 by next-generation sequencing (NGS). Phylogenetic analysis grouped the virus with other NDV class II isolates from subgenotype VI.2.1.2, including two previous NDV isolates from Brazil in 2014 and 2019. The diversity of aminoacid residues at the fusion F protein cleavage site was analyzed identifying the motif RRQKR↓F, typical of virulent strains. Our results all highlight the importance of virus surveillance in wild and domestic birds, especially given the risk of zoonotic NDV.

Unusual SARS-CoV-2 intrahost diversity reveals lineage superinfection.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has infected almost 200 million people worldwide by July 2021 and the pandemic has been characterized by infection waves of viral lineages showing distinct fitness profiles. The simultaneous infection of a single individual by two distinct SARS-CoV-2 lineages may impact COVID-19 disease progression and provides a window of opportunity for viral recombination and the emergence of new lineages with differential phenotype. Several hundred SARS-CoV-2 lineages are currently well phylogenetically defined, but two main factors have precluded major coinfection/codetection and recombination analysis thus far: (i) the low diversity of SARS-CoV-2 lineages during the first year of the pandemic, which limited the identification of lineage defining mutations necessary to distinguish coinfecting/recombining viral lineages; and the (ii) limited availability of raw sequencing data where abundance and distribution of intrasample/intrahost variability can be accessed. Here, we assembled a large sequencing dataset from Brazilian samples covering a period of 18 May 2020 to 30 April 2021 and probed it for unexpected patterns of high intrasample/intrahost variability. This approach enabled us to detect nine cases of SARS-CoV-2 coinfection with well characterized lineage-defining mutations, representing 0.61 % of all samples investigated. In addition, we matched these SARS-CoV-2 coinfections with spatio-temporal epidemiological data confirming its plausibility with the cocirculating lineages at the timeframe investigated. Our data suggests that coinfection with distinct SARS-CoV-2 lineages is a rare phenomenon, although it is certainly a lower bound estimate considering the difficulty to detect coinfections with very similar SARS-CoV-2 lineages and the low number of samples sequenced from the total number of infections.

Lying in wait: the resurgence of dengue virus after the Zika epidemic in Brazil.

After the Zika virus (ZIKV) epidemic in the Americas in 2016, both Zika and dengue incidence declined to record lows in many countries in 2017-2018, but in 2019 dengue resurged in Brazil, causing ~2.1 million cases. In this study we use epidemiological, climatological and genomic data to investigate dengue dynamics in recent years in Brazil. First, we estimate dengue virus force of infection (FOI) and model mosquito-borne transmission suitability since the early 2000s. Our estimates reveal that DENV transmission was low in 2017-2018, despite conditions being suitable for viral spread. Our study also shows a marked decline in dengue susceptibility between 2002 and 2019, which could explain the synchronous decline of dengue in the country, partially as a result of protective immunity from prior ZIKV and/or DENV infections. Furthermore, we performed phylogeographic analyses using 69 newly sequenced genomes of dengue virus serotype 1 and 2 from Brazil, and found that the outbreaks in 2018-2019 were caused by local DENV lineages that persisted for 5-10 years, circulating cryptically before and after the Zika epidemic. We hypothesize that DENV lineages may circulate at low transmission levels for many years, until local conditions are suitable for higher transmission, when they cause major outbreaks.

A Sanger-based approach for scaling up screening of SARS-CoV-2 variants of interest and concern.

The global spread of new SARS-CoV-2 variants of concern underscore an urgent need of simple deployed molecular tools that can differentiate these lineages. Several tools and protocols have been shared since the beginning of the COVID-19 pandemic, but they need to be timely adapted to cope with SARS-CoV-2 evolution. Although whole-genome sequencing (WGS) of the virus genetic material has been widely used, it still presents practical difficulties such as high cost, shortage of available reagents in the global market, need of a specialized laboratorial infrastructure and well-trained staff. These limitations result in SARS-CoV-2 surveillance blackouts across several countries. Here we propose a rapid and accessible protocol based on Sanger sequencing of a single PCR fragment that is able to identify and discriminate all SARS-CoV-2 variants of concern (VOCs) identified so far, according to each characteristic mutational profile at the Spike-RBD region (K417N/T, E484K, N501Y, A570D). Twelve COVID-19 samples from Brazilian patients were evaluated for both WGS and Sanger sequencing: three P.2, two P.1, six B.1.1 and one B.1.1.117 lineage. All results from the Sanger sequencing method perfectly matched the mutational profile of VOCs and non-VOCs RBD's characterized by WGS. In summary, this approach allows a much broader network of laboratories to perform molecular surveillance of SARS-CoV-2 VOCs and report results within a shorter time frame, which is of utmost importance in the context of rapid public health decisions in a fast evolving worldwide pandemic.

Metatranscriptomic analysis identifies different viral-like sequences in two neotropical Mansoniini mosquito species.

Mosquitoes interact with a wide range of viruses including both arboviruses and insect-specific viruses. This study aimed to characterize the RNA viruses that are interacting with Mansonia wilsoni and Coquillettidia hermanoi mosquito species. The total RNA extracted from mosquito pools were sequenced on a Ion torrent platform. Viral contigs were identified against viral databases and their evolutionary relationship were reconstructed. We identified a total of 107 viral sequences, 11 of which were assigned as endogenous viral elements, and at least six known viral families were identified. Phylogenetic reconstructions were performed for 4 viral families. All Mansoniini viruses investigated through phylogenetic analysis are closely related to insect-specific viruses found in other mosquito species although with considerable divergence at the amino acid level, suggesting that we have detected new viral lineages. This study enhanced our understanding about the virome of two sylvatic Mansoniini mosquitoes.

Multiple Introductions Followed by Ongoing Community Spread of SARS-CoV-2 at One of the Largest Metropolitan Areas of Northeast Brazil.

Multiple epicenters of the SARS-CoV-2 pandemic have emerged since the first pneumonia cases in Wuhan, China, such as Italy, USA, and Brazil. Brazil is the third-most affected country worldwide, but genomic sequences of SARS-CoV-2 strains are mostly restricted to states from the Southeast region. Pernambuco state, located in the Northeast region, is the sixth most affected Brazilian state, but very few genomic sequences from the strains circulating in this region are available. We sequenced 101 strains of SARS-CoV-2 from patients presenting Covid-19 symptoms that reside in Pernambuco. Phylogenetic reconstructions revealed that all genomes belong to the B lineage and most of the samples (88%) were classified as lineage B.1.1. We detected multiple viral introductions from abroad (likely from Europe) as well as six local B.1.1 clades composed by Pernambuco only strains. Local clades comprise sequences from the capital city (Recife) and other country-side cities, corroborating the community spread between different municipalities of the state. These findings demonstrate that different from Southeastern Brazilian states where the epidemics were majorly driven by one dominant lineage (B.1.1.28 or B.1.1.33), the early epidemic phase at the Pernambuco state was driven by multiple B.1.1 lineages seeded through both national and international traveling.

Culicidae evolutionary history focusing on the Culicinae subfamily based on mitochondrial phylogenomics.

Mosquitoes are insects of medical importance due their role as vectors of different pathogens to humans. There is a lack of information about the evolutionary history and phylogenetic positioning of the majority of mosquito species. Here we characterized the mitogenomes of mosquito species through low-coverage whole genome sequencing and data mining. A total of 37 draft mitogenomes of different species were assembled from which 16 are newly-sequenced species. We datamined additional 49 mosquito mitogenomes, and together with our 37 mitogenomes, we reconstructed the evolutionary history of 86 species including representatives from 15 genera and 7 tribes. Our results showed that most of the species clustered in clades with other members of their own genus with exception of Aedes genus which was paraphyletic. We confirmed the monophyletic status of the Mansoniini tribe including both Coquillettidia and Mansonia genus. The Aedeomyiini and Uranotaeniini were consistently recovered as basal to other tribes in the subfamily Culicinae, although the exact relationships among these tribes differed between analyses. These results demonstrate that low-coverage sequencing is effective to recover mitogenomes, establish phylogenetic knowledge and hence generate basic fundamental information that will help in the understanding of the role of these species as pathogen vectors.

Sequencing of ZIKV genomes directly from Ae. aegypti and Cx. quinquefasciatus mosquitoes collected during the 2015-16 epidemics in Recife.

Zika virus (ZIKV) is a negative sense RNA virus from the Flaviviridae family, which was relatively unknown until the first human epidemic in Micronesia, in 2007. Since then, it spread to French Polynesia and the Americas. Recife, the capital of Pernambuco state and epicenter of the Zika epidemic in Brazil, experienced a large number of microcephaly cases and other congenital abnormalities associated to the ZIKV infection from, 2015 to 16. Evidences suggest that both Aedes aegypti and Culex quinquefasciatus mosquitoes from Recife are capable of replicating and transmitting the virus. Here, we conducted high throughput sequencing of ZIKV genomes directly from Ae. aegypti and Cx. quinquefasciatus mosquitoes collected during the ZIKV epidemics in Recife, in order to investigate the variability and evolution of the virus. We obtained 11 draft ZIKV genomes derived from 5 pools from each Ae. aegypti and Cx. quinquefasciatus species. Genome coverage breadth ranged from 16 to 100% and average depth from 45 to 46,584×. Two of these genomes were obtained from pools of Cx. quinquefasciatus females with no sign of blood in the abdomen. Amino acid substitutions found here were not species-specific. In addition, molecular clock dating estimated that ZIKV draft genomes obtained here were co-circulating in other regions of the country during the epidemics. Overall results highlight that viral mutations and even minor variants can be detected in genomes directly sequenced from mosquito samples and insights about natural viral genomic variability and viral evolution can be useful when designing tools for mosquito control programs.

Spread of two Zika virus lineages in Midwest Brazil.

Zika virus (ZIKV) has been intensively studied in South America and across the globe since 2015-2016 epidemics. However, in Brazil - the largest and the most affected country in terms of human infection by this virus, most of the viral molecular information is restricted to metropolitan centers distributed along the Brazilian coast and almost no information is known about the virus spread in most difficult access areas such as the Midwest region of the country. Here, we report two ZIKV complete genomes from samples obtained during arboviral surveillance at the Sinop city, southern border of the Amazonian forest, Midwest Brazil in 2015. Our results show that the virus was introduced in this region through two independent introductions: one occurred at the end of 2014, around the period that the virus was already distributed in other regions of the country and abroad, and a second at the end of 2015. Moreover, these genomes were clustered with other viral strains sampled at distant Brazilian states in line with other findings about the rapid spread of the virus throughout the country.

Genome sequencing reveals coinfection by multiple chikungunya virus genotypes in a recent outbreak in Brazil.

Chikungunya virus (CHIKV) is an RNA virus from the Togaviridae family transmitted by mosquitoes in both sylvatic and urban cycles. In humans, CHIKV infection leads to a febrile illness, denominated Chikungunya fever (CHIKF), commonly associated with more intense and debilitating outcomes. CHIKV arrived in Brazil in 2014 through two independent introductions: the Asian/Caribbean genotype entered through the North region and the African ECSA genotype was imported through the Northeast region. Following their initial introduction, both genotypes established their urban cycle among large naive human populations causing several outbreaks in the Americas. Here, we sequenced CHIKV genomes from a recent outbreak in the Northeast region of Brazil, employing an in-house developed Next-Generation Sequencing (NGS) protocol capable of directly detecting multiple known CHIKV genotypes from clinical positive samples. Our results demonstrate that both Asian/Caribbean and ECSA genotypes expanded their ranges, reaching cocirculation in the Northeast region of Brazil. In addition, our NGS data supports the findings of simultaneous infection by these two genotypes, suggesting that coinfection might be more common than previously thought in highly endemic areas. Future efforts to understand CHIKV epidemiology should thus take into consideration the possibility of coinfection by different genotypes in the human population.

Zika virus detection, isolation and genome sequencing through Culicidae sampling during the epidemic in Vitória, Espírito Santo, Brazil.

BACKGROUND: Zika virus (ZIKV) has been isolated from many mosquito species in nature, but it is believed that the main vectors in urban environments are species of the genus Aedes. Here, we detected and isolated ZIKV in samples from Aedes aegypti, Aedes taeniorhynchus and Culex quinquefasciatus, collected during the Zika epidemic in Vitória, southeast Brazil. Using quantitative real-time polymerase chain reaction, ZIKV detection was performed in mosquito samples collected from February to April 2016. RESULTS: Overall, six pools of mosquitoes were positive for ZIKV: four of Cx. quinquefasciatus, one of Ae. aegypti and one of Ae. taeniorhynchus. Their genomes were sequenced. CONCLUSIONS: These results support and strengthen the hypothesis that other mosquito species can also be involved in ZIKV transmission.

Phylogenetic positioning of the Antarctic alga Prasiola crispa (Trebouxiophyceae) using organellar genomes and their structural analysis.

Antarctica is one of the most difficult habitats for sustaining life on earth; organisms that live there have developed different strategies for survival. Among these organisms is the green alga Prasiola crispa, belonging to the class Trebouxiophyceae. The literature on P. crispa taxonomy is scarce, and many gaps in the evolutionary relationship with its closest relatives remain. The goal of this study was to analyze the evolutionary relationships between P. crispa and other green algae using plastid and mitochondrial genomes. In addition, we analyzed the synteny conservation of these genomes of P. crispa with those of closely related species. Based on the plastid genome, P. crispa grouped with Prasiolopsis sp. SAG 84.81, another Trebouxiophyceaen species from the Prasiola clade. Based on the mitochondrial genome analysis, P. crispa grouped with other Trebouxiophyceaen species but had a basal position. The structure of the P. crispa chloroplast genome had low synteny with Prasiolopsis sp. SAG 84.81, despite some conserved gene blocks. The same was observed in the mitochondrial genome compared with Coccomyxa subellipsoidea C-169. We were able to establish the phylogenetic position of P. crispa with other species of Trebouxiophyceae using its genomes. In addition, we described the plasticity of these genomes using a structural analysis. The plastid and mitochondrial genomes of P. crispa will be useful for further genetic studies, phylogenetic analysis and resource protection of P. crispa as well as for further phylogenetic analysis of Trebouxiophyceaen green algae.

Draft Plastid and Mitochondrial Genome Sequences from Antarctic Alga Prasiola crispa.

The organelle genomes of the Antarctic alga Prasiola crispa (Lightfoot) Kützing have been sequenced. The plastid and mitochondrial genomes have a total length of 196,502 bp and 89,819 bp, respectively. These genomes have 19 putative photosynthesis-related genes and 17 oxidative metabolism-related genes, respectively.