Molecular Epidemiology of Mayaro Virus among Febrile Patients, Roraima State, Brazil, 2018-2021

We detected Mayaro virus (MAYV) in 3.4% (28/822) of febrile patients tested during 2018-2021 from Roraima State, Brazil. We also isolated MAYV strains and confirmed that these cases were caused by genotype D. Improved surveillance is needed to better determine the burden of MAYV in the Amazon Region.

Pathophysiology of chikungunya virus infection associated with fatal outcomes

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that causes acute, subacute, and chronic human arthritogenic diseases and, in rare instances, can lead to neurological complications and death. Here, we combined epidemiological, virological, histopathological, cytokine, molecular dynamics, metabolomic, proteomic, and genomic analyses to investigate viral and host factors that contribute to chikungunya-associated (CHIK) death. Our results indicate that CHIK deaths are associated with multi-organ infection, central nervous system damage, and elevated serum levels of pro-inflammatory cytokines and chemokines compared with survivors. The histopathologic, metabolite, and proteomic signatures of CHIK deaths reveal hemodynamic disorders and dysregulated immune responses. The CHIKV East-Central-South-African lineage infecting our study population causes both fatal and survival cases. Additionally, CHIKV infection impairs the integrity of the blood-brain barrier, as evidenced by an increase in permeability and altered tight junction protein expression. Overall, our findings improve the understanding of CHIK pathophysiology and the causes of fatal infections.

Clinical and epidemiological characteristics of Madariaga and Venezuelan equine encephalitis virus infections

Madariaga virus (MADV) and Venezuelan equine encephalitis virus (VEEV) are emerging arboviruses affecting rural and remote areas of Latin America. However, there are limited clinical and epidemiological reports available, and outbreaks are occurring at an increasing frequency. We addressed this gap by analyzing all the available clinical and epidemiological data of MADV and VEEV infections recorded since 1961 in Panama. A total of 168 of human alphavirus encephalitis cases were detected in Panama from 1961 to 2023. Here we describe the clinical signs and symptoms and epidemiological characteristics of these cases, and also explored signs and symptoms as potential predictors of encephalitic alphavirus infection when compared to those of other arbovirus infections occurring in the region. Our results highlight the challenges clinical diagnosis of alphavirus disease in endemic regions with overlapping circulation of multiple arboviruses.

Chikungunya: a decade of burden in the Americas

In the Americas, one decade following its emergence in 2013, chikungunya virus (CHIKV) continues to spread and cause epidemics across the region. To date, 3.7 million suspected and laboratory-confirmed chikungunya cases have been reported in 50 countries or territories in the Americas. Here, we outline the current status and epidemiological aspects of chikungunya in the Americas and discuss prospects for future research and public health strategies to combat CHIKV in the region.

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.

Real-time RT-PCR for Venezuelan equine encephalitis complex, Madariaga, and Eastern equine encephalitis viruses: application in human and mosquito public health surveillance in Panama

Eastern equine encephalitis virus (EEEV), Madariaga virus (MADV), and Venezuelan equine encephalitis virus complex (VEEV) are New World alphaviruses transmitted by mosquitoes. They cause febrile and sometimes severe neurological diseases in human and equine hosts. Detecting them during the acute phase is hindered by non-specific symptoms and limited diagnostic tools. We designed and clinically assessed real-time reverse transcription polymerase chain reaction assays (rRT-PCRs) for VEEV complex, MADV, and EEEV using whole-genome sequences. Validation involved 15 retrospective serum samples from 2015 to 2017 outbreaks, 150 mosquito pools from 2015, and 118 prospective samples from 2021 to 2022 surveillance in Panama. The rRT-PCRs detected VEEV complex RNA in 10 samples (66.7%) from outbreaks, with one having both VEEV complex and MADV RNAs. VEEV complex RNA was found in five suspected dengue cases from disease surveillance. The rRT-PCR assays identified VEEV complex RNA in three Culex (Melanoconion) vomerifer pools, leading to VEEV isolates in two. Phylogenetic analysis revealed the VEEV ID subtype in positive samples. Notably, 11.9% of dengue-like disease patients showed VEEV infections. Together, our rRT-PCR validation in human and mosquito samples suggests that this method can be incorporated into mosquito and human encephalitic alphavirus surveillance programs in endemic regions.

Correction: Implementation of Genomic Surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings

[This corrects the article DOI: 10.1371/journal.pgph.0001455.].

Scale-free dynamics of COVID-19 in a Brazilian city

A common basis to address the dynamics of directly transmitted infectious diseases, such as COVID-19, are compartmental (or SIR) models. SIR models typically assume homogenous population mixing, a simplification that is convenient but unrealistic. Here we validate an existing model of a scale-free fractal infection process using high-resolution data on COVID-19 spread in São Caetano, Brazil. We find that transmission can be described by a network in which each infectious individual has a small number of susceptible contacts, of the order of 2-5. This model parameter correlated tightly with physical distancing measured by mobile phone data, such that in periods of greater distancing the model recovered a lower average number of contacts, and vice versa. We show that the SIR model is a special case of our scale-free fractal process model in which the parameter that reflects population structure is set at unity, indicating homogeneous mixing. Our more general framework better explained the dynamics of COVID-19 in São Caetano, used fewer parameters than a standard SIR model and accounted for geographically localized clusters of disease. Our model requires further validation in other locations and with other directly transmitted infectious agents.

Convergent trends and spatiotemporal patterns of Aedes-borne arboviruses in Mexico and Central America

BACKGROUND

Aedes-borne arboviruses cause both seasonal epidemics and emerging outbreaks with a significant impact on global health. These viruses share mosquito vector species, often infecting the same host population within overlapping geographic regions. Thus, comparative analyses of the virus evolutionary and epidemiological dynamics across spatial and temporal scales could reveal convergent trends.

METHODOLOGY/PRINCIPAL FINDINGS

Focusing on Mexico as a case study, we generated novel chikungunya and dengue (CHIKV, DENV-1 and DENV-2) virus genomes from an epidemiological surveillance-derived historical sample collection, and analysed them together with longitudinally-collected genome and epidemiological data from the Americas. Aedes-borne arboviruses endemically circulating within the country were found to be introduced multiple times from lineages predominantly sampled from the Caribbean and Central America. For CHIKV, at least thirteen introductions were inferred over a year, with six of these leading to persistent transmission chains. For both DENV-1 and DENV-2, at least seven introductions were inferred over a decade.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that CHIKV, DENV-1 and DENV-2 in Mexico share evolutionary and epidemiological trajectories. The southwest region of the country was determined to be the most likely location for viral introductions from abroad, with a subsequent spread into the Pacific coast towards the north of Mexico. Virus diffusion patterns observed across the country are likely driven by multiple factors, including mobility linked to human migration from Central towards North America. Considering Mexico's geographic positioning displaying a high human mobility across borders, our results prompt the need to better understand the role of anthropogenic factors in the transmission dynamics of Aedes-borne arboviruses, particularly linked to land-based human migration.

Madariaga and Venezuelan equine encephalitis virus seroprevalence in rodent enzootic hosts in Eastern and Western Panama

While rodents are primary reservoirs of Venezuelan equine encephalitis virus (VEEV), their role in Madariaga virus (MADV) transmission remains uncertain, particularly given their overlapping geographic distribution. This study explores the interplay of alphavirus prevalence, rodent diversity, and land use within Darien and Western Panama provinces. A total of three locations were selected for rodent sampling in Darien province: Los Pavitos, El Real de Santa Maria and Santa Librada. Two sites were selected in Western Panama province: El Cacao and Cirí Grande. We used plaque reduction neutralization tests to assess MADV and VEEV seroprevalences in 599 rodents of 16 species across five study sites. MADV seroprevalence was observed at higher rates in Los Pavitos (Darien province), 9.0%, 95% CI: 3.6-17.6, while VEEV seroprevalence was elevated in El Cacao (Western Panama province), 27.3%, 95% CI: 16.1-40.9, and El Real de Santa María (Darien province), 20.4%, 95% CI: 12.6-29.7. Species like Oryzomys coesi, 23.1%, 95% CI: 5.0-53.8, and Transandinomys bolivaris, 20.0%, 95% CI: 0.5-71.6 displayed higher MADV seroprevalences than other species, whereas Transandinomys bolivaris, 80.0%, 95% CI: 28.3-99.4, and Proechimys semispinosus, 27.3%, 95% CI: 17.0-39.6, exhibited higher VEEV seroprevalences. Our findings provide support to the notion that rodents are vertebrate reservoirs of MADV and reveal spatial variations in alphavirus seropositivity among rodent species, with different provinces exhibiting distinct rates for MADV and VEEV. Moreover, specific rodent species are linked to unique seroprevalence patterns for these viruses, suggesting that rodent diversity and environmental conditions might play a significant role in shaping alphavirus distribution.

Dispersal patterns and influence of air travel during the global expansion of SARS-CoV-2 variants of concern

The Alpha, Beta, and Gamma SARS-CoV-2 variants of concern (VOCs) co-circulated globally during 2020 and 2021, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, which, in turn, was displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC and identify countries that acted as global and regional hubs of dissemination. We demonstrate the declining role of presumed origin countries of VOCs in their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1%-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of its emergence, associated with accelerated passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants, with implications for genomic surveillance along the hierarchical airline network.

Spatiotemporal dynamics and recurrence of chikungunya virus in Brazil: an epidemiological study

BACKGROUND

Chikungunya virus (CHIKV) is an Aedes mosquito-borne virus that has caused large epidemics linked to acute, chronic, and severe clinical outcomes. Currently, Brazil has the highest number of chikungunya cases in the Americas. We aimed to investigate the spatiotemporal dynamics and recurrence pattern of chikungunya in Brazil since its introduction in 2013.

METHODS

In this epidemiological study, we used CHIKV genomic sequencing data, CHIKV vector information, and aggregate clinical data on chikungunya cases from Brazil. The genomic data comprised 241 Brazilian CHIKV genome sequences from GenBank (n=180) and the 2022 CHIKV outbreak in Ceará state (n=61). The vector data (Breteau index and House index) were obtained from the Brazilian Ministry of Health for all 184 municipalities in Ceará state and 116 municipalities in Tocantins state in 2022. Epidemiological data on laboratory-confirmed cases of chikungunya between 2013 and 2022 were obtained from the Brazilian Ministry of Health and Laboratory of Public Health of Ceará. We assessed the spatiotemporal dynamics of chikungunya in Brazil via time series, mapping, age-sex distribution, cumulative case-fatality, linear correlation, logistic regression, and phylogenetic analyses.

FINDINGS

Between March 3, 2013, and June 4, 2022, 253 545 laboratory-confirmed chikungunya cases were reported in 3316 (59·5%) of 5570 municipalities, mainly distributed in seven epidemic waves from 2016 to 2022. To date, Ceará in the northeast has been the most affected state, with 77 418 cases during the two largest epidemic waves in 2016 and 2017 and the third wave in 2022. From 2016 to 2022 in Ceará, the odds of being CHIKV-positive were higher in females than in males (odds ratio 0·87, 95% CI 0·85-0·89, p<0·0001), and the cumulative case-fatality ratio was 1·3 deaths per 1000 cases. Chikungunya recurrences in the states of Ceará, Tocantins (recurrence in 2022), and Pernambuco (recurrence in 2021) were limited to municipalities with few or no previously reported cases in the previous epidemic waves. The recurrence of chikungunya in Ceará in 2022 was associated with a new East-Central-South-African lineage. Population density metrics of the main CHIKV vector in Brazil, Aedes aegypti, were not correlated spatially with locations of chikungunya recurrence in Ceará and Tocantins.

INTERPRETATION

Spatial heterogeneity of CHIKV spread and population immunity might explain the recurrence pattern of chikungunya in Brazil. These results can be used to inform public health interventions to prevent future chikungunya epidemic waves in urban settings.

FUNDING

Global Virus Network, Burroughs Wellcome Fund, Wellcome Trust, US National Institutes of Health, São Paulo Research Foundation, Brazil Ministry of Education, UK Medical Research Council, Brazilian National Council for Scientific and Technological Development, and UK Royal Society.

TRANSLATION

For the Portuguese translation of the abstract see Supplementary Materials section.

Implementation of genomic surveillance of SARS-CoV-2 in the Caribbean: Lessons learned for sustainability in resource-limited settings

The COVID-19 pandemic highlighted the importance of global genomic surveillance to monitor the emergence and spread of SARS-CoV-2 variants and inform public health decision-making. Until December 2020 there was minimal capacity for viral genomic surveillance in most Caribbean countries. To overcome this constraint, the COVID-19: Infectious disease Molecular epidemiology for PAthogen Control & Tracking (COVID-19 IMPACT) project was implemented to establish rapid SARS-CoV-2 whole genome nanopore sequencing at The University of the West Indies (UWI) in Trinidad and Tobago (T&T) and provide needed SARS-CoV-2 sequencing services for T&T and other Caribbean Public Health Agency Member States (CMS). Using the Oxford Nanopore Technologies MinION sequencing platform and ARTIC network sequencing protocols and bioinformatics pipeline, a total of 3610 SARS-CoV-2 positive RNA samples, received from 17 CMS, were sequenced in-situ during the period December 5th 2020 to December 31st 2021. Ninety-one Pango lineages, including those of five variants of concern (VOC), were identified. Genetic analysis revealed at least 260 introductions to the CMS from other global regions. For each of the 17 CMS, the percentage of reported COVID-19 cases sequenced by the COVID-19 IMPACT laboratory ranged from 0·02% to 3·80% (median = 1·12%). Sequences submitted to GISAID by our study represented 73·3% of all SARS-CoV-2 sequences from the 17 CMS available on the database up to December 31st 2021. Increased staffing, process and infrastructural improvement over the course of the project helped reduce turnaround times for reporting to originating institutions and sequence uploads to GISAID. Insights from our genomic surveillance network in the Caribbean region directly influenced non-pharmaceutical countermeasures in the CMS countries. However, limited availability of associated surveillance and clinical data made it challenging to contextualise the observed SARS-CoV-2 diversity and evolution, highlighting the need for development of infrastructure for collecting and integrating genomic sequencing data and sample-associated metadata.

Global Expansion of SARS-CoV-2 Variants of Concern: Dispersal Patterns and Influence of Air Travel

In many regions of the world, the Alpha, Beta and Gamma SARS-CoV-2 Variants of Concern (VOCs) co-circulated during 2020-21 and fueled waves of infections. During 2021, these variants were almost completely displaced by the Delta variant, causing a third wave of infections worldwide. This phenomenon of global viral lineage displacement was observed again in late 2021, when the Omicron variant disseminated globally. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of SARS-CoV-2 VOCs worldwide. We find that the source-sink dynamics of SARS-CoV-2 varied substantially by VOC, and identify countries that acted as global hubs of variant dissemination, while other countries became regional contributors to the export of specific variants. We demonstrate a declining role of presumed origin countries of VOCs to their global dispersal: we estimate that India contributed <15% of all global exports of Delta to other countries and South Africa <1-2% of all global Omicron exports globally. We further estimate that >80 countries had received introductions of Omicron BA.1 100 days after its inferred date of emergence, compared to just over 25 countries for the Alpha variant. This increased speed of global dissemination was associated with a rebound in air travel volume prior to Omicron emergence in addition to the higher transmissibility of Omicron relative to Alpha. Our study highlights the importance of global and regional hubs in VOC dispersal, and the speed at which highly transmissible variants disseminate through these hubs, even before their detection and characterization through genomic surveillance.

Highlights

Global phylogenetic analysis reveals relationship between air travel and speed of dispersal of SARS-CoV-2 variants of concern (VOCs)Omicron VOC spread to 5x more countries within 100 days of its emergence compared to all other VOCsOnward transmission and dissemination of VOCs Delta and Omicron was primarily from secondary hubs rather than initial country of detection during a time of increased global air travelAnalysis highlights highly connected countries identified as major global and regional exporters of VOCs.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.

First Report of Wenzhou sobemo-like virus 4 in Aedes albopictus (Diptera: Culicidae) in Latin America

Insect-specific viruses (ISVs) are viruses that replicate exclusively in arthropod cells. Many ISVs have been studied in mosquitoes as many of them act as vectors for human etiological agents, such as arboviruses. Aedes (Stegomyia) albopictus is an important potential vector of several arboviruses in Brazil, such as dengue (DENV), Zika (ZIKV) and chikungunya (CHIKV). The development of next-generation sequencing metagenomics has enabled the discovery and characterization of new ISVs. Ae. albopictus eggs were collected using oviposition traps placed in two urban parks in the city of São Paulo, Brazil. The Aedes albopictus females were divided into pools and the genetic material was extracted and processed for sequencing by metagenomics. Complete genomes of ISV Wenzhou sobemo-like virus 4 (WSLV4) were obtained in three of the four pools tested. This is the first detection of ISV WSLV4 in Ae. albopictus females in Latin America. Further studies on ISVs in Ae. albopictus are needed to better understand the role of this species in the dynamics of arbovirus transmission in the Americas.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

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.

SARS-CoV-2 antibody dynamics in blood donors and COVID-19 epidemiology in eight Brazilian state capitals: A serial cross-sectional study

Background

The COVID-19 situation in Brazil is complex due to large differences in the shape and size of regional epidemics. Understanding these patterns is crucial to understand future outbreaks of SARS-CoV-2 or other respiratory pathogens in the country.

Methods

We tested 97,950 blood donation samples for IgG antibodies from March 2020 to March 2021 in 8 of Brazil's most populous cities. Residential postal codes were used to obtain representative samples. Weekly age- and sex-specific seroprevalence were estimated by correcting the crude seroprevalence by test sensitivity, specificity, and antibody waning.

Results

The inferred attack rate of SARS-CoV-2 in December 2020, before the Gamma variant of concern (VOC) was dominant, ranged from 19.3% (95% credible interval [CrI] 17.5-21.2%) in Curitiba to 75.0% (95% CrI 70.8-80.3%) in Manaus. Seroprevalence was consistently smaller in women and donors older than 55 years. The age-specific infection fatality rate (IFR) differed between cities and consistently increased with age. The infection hospitalisation rate increased significantly during the Gamma-dominated second wave in Manaus, suggesting increased morbidity of the Gamma VOC compared to previous variants circulating in Manaus. The higher disease penetrance associated with the health system's collapse increased the overall IFR by a minimum factor of 2.91 (95% CrI 2.43-3.53).

Conclusions

These results highlight the utility of blood donor serosurveillance to track epidemic maturity and demonstrate demographic and spatial heterogeneity in SARS-CoV-2 spread.

Funding

This work was supported by Itaú Unibanco 'Todos pela Saude' program; FAPESP (grants 18/14389-0, 2019/21585-0); Wellcome Trust and Royal Society Sir Henry Dale Fellowship 204311/Z/16/Z; the Gates Foundation (INV- 034540 and INV-034652); REDS-IV-P (grant HHSN268201100007I); the UK Medical Research Council (MR/S0195/1, MR/V038109/1); CAPES; CNPq (304714/2018-6); Fundação Faculdade de Medicina; Programa Inova Fiocruz-CE/Funcap - Edital 01/2020 Number: FIO-0167-00065.01.00/20 SPU N°06531047/2020; JBS - Fazer o bem faz bem.

Predicting SARS-CoV-2 Variant Spread in a Completely Seropositive Population Using Semi-Quantitative Antibody Measurements in Blood Donors

SARS-CoV-2 serologic surveys estimate the proportion of the population with antibodies against historical variants, which nears 100% in many settings. New approaches are required to fully exploit serosurvey data. Using a SARS-CoV-2 anti-Spike (S) protein chemiluminescent microparticle assay, we attained a semi-quantitative measurement of population IgG titers in serial cross-sectional monthly samples of blood donations across seven Brazilian state capitals (March 2021−November 2021). Using an ecological analysis, we assessed the contributions of prior attack rate and vaccination to antibody titer. We compared anti-S titer across the seven cities during the growth phase of the Delta variant and used this to predict the resulting age-standardized incidence of severe COVID-19 cases. We tested ~780 samples per month, per location. Seroprevalence rose to >95% across all seven capitals by November 2021. Driven by vaccination, mean antibody titer increased 16-fold over the study, with the greatest increases occurring in cities with the highest prior attack rates. Mean anti-S IgG was strongly correlated (adjusted R2 = 0.89) with the number of severe cases caused by Delta. Semi-quantitative anti-S antibody titers are informative about prior exposure and vaccination coverage and may also indicate the potential impact of future SARS-CoV-2 variants.

Understanding the Potential Impact of Different Drug Properties on Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Transmission and Disease Burden: A Modelling Analysis

BACKGROUND

The public health impact of the coronavirus disease 2019 (COVID-19) pandemic has motivated a rapid search for potential therapeutics, with some key successes. However, the potential impact of different treatments, and consequently research and procurement priorities, have not been clear.

METHODS

Using a mathematical model of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission, COVID-19 disease and clinical care, we explore the public-health impact of different potential therapeutics, under a range of scenarios varying healthcare capacity, epidemic trajectories; and drug efficacy in the absence of supportive care.

RESULTS

The impact of drugs like dexamethasone (delivered to the most critically-ill in hospital and whose therapeutic benefit is expected to depend on the availability of supportive care such as oxygen and mechanical ventilation) is likely to be limited in settings where healthcare capacity is lowest or where uncontrolled epidemics result in hospitals being overwhelmed. As such, it may avert 22% of deaths in high-income countries but only 8% in low-income countries (assuming R = 1.35). Therapeutics for different patient populations (those not in hospital, early in the course of infection) and types of benefit (reducing disease severity or infectiousness, preventing hospitalization) could have much greater benefits, particularly in resource-poor settings facing large epidemics.

CONCLUSIONS

Advances in the treatment of COVID-19 to date have been focused on hospitalized-patients and predicated on an assumption of adequate access to supportive care. Therapeutics delivered earlier in the course of infection that reduce the need for healthcare or reduce infectiousness could have significant impact, and research into their efficacy and means of delivery should be a priority.

Emergence and widespread circulation of a recombinant SARS-CoV-2 lineage in North America

Although recombination is a feature of coronavirus evolution, previously detected recombinant lineages of SARS-CoV-2 have shown limited circulation thus far. Here, we present a detailed phylogenetic analysis of four SARS-CoV-2 lineages to investigate the possibility of virus recombination among them. Our analyses reveal well-supported phylogenetic differences between the Orf1ab region encoding viral non-structural proteins and the rest of the genome, including Spike (S) protein and remaining reading frames. By accounting for several deletions in NSP6, Orf3a, and S, we conclude that the B.1.628 major cluster, now designated as lineage XB, originated from a recombination event between viruses of B.1.631 and B.1.634 lineages. This scenario is supported by the spatiotemporal distribution of these lineages across the USA and Mexico during 2021, suggesting that the recombination event originated in this geographical region. This event raises important questions regarding the role and potential effects of recombination on SARS-CoV-2 evolution.

Understanding Sabiá virus infections (Brazilian mammarenavirus)

BACKGROUND

Only two naturally occurring human Sabiá virus (SABV) infections have been reported, and those occurred over 20 years ago.

METHODS

We diagnosed two new cases of SABV infection using metagenomics in patients thought to have severe yellow fever and described new features of histopathological findings.

RESULTS

We characterized clinical manifestations, histopathology and analyzed possible nosocomial transmission. Patients presented with hepatitis, bleeding, neurological alterations and died. We traced twenty-nine hospital contacts and evaluated them clinically and by RT-PCR and neutralizing antibodies. Autopsies uncovered unique features on electron microscopy, such as hepatocyte "pinewood knot" lesions. Although previous reports with similar New-World arenavirus had nosocomial transmission, our data did not find any case in contact tracing.

CONCLUSIONS

Although an apparent by rare, Brazilian mammarenavirus infection is an etiology for acute hemorrhagic fever syndrome. The two fatal cases had peculiar histopathological findings not previously described. The virological diagnosis was possible only by contemporary techniques such as metagenomic assays. We found no subsequent infections when we used serological and molecular tests to evaluate close contacts.

Spatial and temporal fluctuations in COVID-19 fatality rates in Brazilian hospitals

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Gamma variant of concern has spread rapidly across Brazil since late 2020, causing substantial infection and death waves. Here we used individual-level patient records after hospitalization with suspected or confirmed coronavirus disease 2019 (COVID-19) between 20 January 2020 and 26 July 2021 to document temporary, sweeping shocks in hospital fatality rates that followed the spread of Gamma across 14 state capitals, during which typically more than half of hospitalized patients aged 70 years and older died. We show that such extensive shocks in COVID-19 in-hospital fatality rates also existed before the detection of Gamma. Using a Bayesian fatality rate model, we found that the geographic and temporal fluctuations in Brazil's COVID-19 in-hospital fatality rates were primarily associated with geographic inequities and shortages in healthcare capacity. We estimate that approximately half of the COVID-19 deaths in hospitals in the 14 cities could have been avoided without pre-pandemic geographic inequities and without pandemic healthcare pressure. Our results suggest that investments in healthcare resources, healthcare optimization and pandemic preparedness are critical to minimize population-wide mortality and morbidity caused by highly transmissible and deadly pathogens such as SARS-CoV-2, especially in low- and middle-income countries.

Tracking the turnover of SARS-CoV-2 VOCs Gamma to Delta in a Brazilian state (Minas Gerais) with a high-vaccination status

The emergence and global dissemination of Severe Acute Respiratory Syndrome virus 2 (SARS-CoV-2) variants of concern (VOCs) have been described as the main factor driving the Coronavirus Disease 2019 pandemic. In Brazil, the Gamma variant dominated the epidemiological scenario during the first period of 2021. Many Brazilian regions detected the Delta variant after its first description and documented its spread. To monitor the introduction and spread of VOC Delta, we performed Polymerase Chain Reaction (PCR) genotyping and genome sequencing in ten regional sentinel units from June to October 2021 in the State of Minas Gerais (MG). We documented the introduction and spread of Delta, comprising 70 per cent of the cases 8 weeks later. Comparing the viral loads of the Gamma and Delta dominance periods, we provide additional evidence that the latter is more transmissible. The spread and dominance of Delta did not culminate in the increase in cases and deaths, suggesting that the vaccination may have restrained the epidemic growth. Analysis of 224 novel Delta genomes revealed that Rio de Janeiro state was the primary source for disseminating this variant in the state of MG. We present the establishment of Delta, providing evidence of its enhanced transmissibility and showing that this variant shift did not aggravate the epidemiological scenario in a high immunity setting.

Shotgun metagenomic sequencing of the first case of monkeypox virus in Brazil, 2022

Monkeypox virus (MPXV), a zoonotic virus endemic to the African continent, has been reported in 33 non-endemic countries since May 2022. We report an almost complete genome of the first confirmed case of MPXV in Brazil. Shotgun metagenomic sequencing was completed in 18 hours, from DNA extraction to consensus sequence generation.

Clearance of Persistent SARS-CoV-2 RNA Detection in a NFκB-Deficient Patient in Association with the Ingestion of Human Breast Milk: A Case Report

Currently, there are no evidence-based treatment options for long COVID-19, and it is known that SARS-CoV-2 can persist in part of the infected patients, especially those with immunosuppression. Since there is a robust secretion of SARS-CoV-2-specific highly-neutralizing IgA antibodies in breast milk, and because this immunoglobulin plays an essential role against respiratory virus infection in mucosa cells, being, in addition, more potent in neutralizing SARS-CoV-2 than IgG, here we report the clinical course of an NFκB-deficient patient chronically infected with the SARS-CoV-2 Gamma variant, who, after a non-full effective treatment with plasma infusion, received breast milk from a vaccinated mother by oral route as treatment for COVID-19. After such treatment, the symptoms improved, and the patient was systematically tested negative for SARS-CoV-2. Thus, we hypothesize that IgA and IgG secreted antibodies present in breast milk could be useful to treat persistent SARS-CoV-2 infection in immunodeficient patients.

Epidemiology of COVID-19 after Emergence of SARS-CoV-2 Gamma Variant, Brazilian Amazon, 2020-2021

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Gamma variant has been hypothesized to cause more severe illness than previous variants, especially in children. Successive SARS-CoV-2 IgG serosurveys in the Brazilian Amazon showed that age-specific attack rates and proportions of symptomatic SARS-CoV-2 infections were similar before and after Gamma variant emergence.

Reinfection by the SARS-CoV-2 Gamma variant in blood donors in Manaus, Brazil

BACKGROUND

The city of Manaus, north Brazil, was stricken by a second epidemic wave of SARS-CoV-2 despite high seroprevalence estimates, coinciding with the emergence of the Gamma (P.1) variant. Reinfections were postulated as a partial explanation for the second surge. However, accurate calculation of reinfection rates is difficult when stringent criteria as two time-separated RT-PCR tests and/or genome sequencing are required. To estimate the proportion of reinfections caused by Gamma during the second wave in Manaus and the protection conferred by previous infection, we identified anti-SARS-CoV-2 antibody boosting in repeat blood donors as a mean to infer reinfection.

METHODS

We tested serial blood samples from unvaccinated repeat blood donors in Manaus for the presence of anti-SARS-CoV-2 IgG antibodies using two assays that display waning in early convalescence, enabling the detection of reinfection-induced boosting. Donors were required to have three or more donations, being at least one during each epidemic wave. We propose a strict serological definition of reinfection (reactivity boosting following waning like a V-shaped curve in both assays or three spaced boostings), probable (two separate boosting events) and possible (reinfection detected by only one assay) reinfections. The serial samples were used to divide donors into six groups defined based on the inferred sequence of infection and reinfection with non-Gamma and Gamma variants.

RESULTS

From 3655 repeat blood donors, 238 met all inclusion criteria, and 223 had enough residual sample volume to perform both serological assays. We found 13.6% (95% CI 7.0-24.5%) of all presumed Gamma infections that were observed in 2021 were reinfections. If we also include cases of probable or possible reinfections, these percentages increase respectively to 22.7% (95% CI 14.3-34.2%) and 39.3% (95% CI 29.5-50.0%). Previous infection conferred a protection against reinfection of 85.3% (95% CI 71.3-92.7%), decreasing to respectively 72.5% (95% CI 54.7-83.6%) and 39.5% (95% CI 14.1-57.8%) if probable and possible reinfections are included.

CONCLUSIONS

Reinfection by Gamma is common and may play a significant role in epidemics where Gamma is prevalent, highlighting the continued threat variants of concern pose even to settings previously hit by substantial epidemics.

Mapping environmental suitability of Haemagogus and Sabethes spp. mosquitoes to understand sylvatic transmission risk of yellow fever virus in Brazil

BACKGROUND

Yellow fever (YF) is an arboviral disease which is endemic to Brazil due to a sylvatic transmission cycle maintained by infected mosquito vectors, non-human primate (NHP) hosts, and humans. Despite the existence of an effective vaccine, recent sporadic YF epidemics have underscored concerns about sylvatic vector surveillance, as very little is known about their spatial distribution. Here, we model and map the environmental suitability of YF's main vectors in Brazil, Haemagogus spp. and Sabethes spp., and use human population and NHP data to identify locations prone to transmission and spillover risk.

METHODOLOGY/PRINCIPAL FINDINGS

We compiled a comprehensive set of occurrence records on Hg. janthinomys, Hg. leucocelaenus, and Sabethes spp. from 1991-2019 using primary and secondary data sources. Linking these data with selected environmental and land-cover variables, we adopted a stacked regression ensemble modelling approach (elastic-net regularized GLM, extreme gradient boosted regression trees, and random forest) to predict the environmental suitability of these species across Brazil at a 1 km x 1 km resolution. We show that while suitability for each species varies spatially, high suitability for all species was predicted in the Southeastern region where recent outbreaks have occurred. By integrating data on NHP host reservoirs and human populations, our risk maps further highlight municipalities within the region that are prone to transmission and spillover.

CONCLUSIONS/SIGNIFICANCE

Our maps of sylvatic vector suitability can help elucidate potential locations of sylvatic reservoirs and be used as a tool to help mitigate risk of future YF outbreaks and assist in vector surveillance. Furthermore, at-risk regions identified from our work could help disease control and elucidate gaps in vaccination coverage and NHP host surveillance.

Tracking the emergence of disparities in the subnational spread of COVID-19 in Brazil using an online application for real-time data visualisation: A longitudinal analysis

BACKGROUND

Brazil is one of the countries worst affected by the COVID-19 pandemic with over 20 million cases and 557,000 deaths reported by August 2021. Comparison of real-time local COVID-19 data between areas is essential for understanding transmission, measuring the effects of interventions, and predicting the course of the epidemic, but are often challenging due to different population sizes and structures.

METHODS

We describe the development of a new app for the real-time visualisation of COVID-19 data in Brazil at the municipality level. In the CLIC-Brazil app, daily updates of case and death data are downloaded, age standardised and used to estimate the effective reproduction number (Rt ). We show how such platforms can perform real-time regression analyses to identify factors associated with the rate of initial spread and early reproduction number. We also use survival methods to predict the likelihood of occurrence of a new peak of COVID-19 incidence.

FINDINGS

After an initial introduction in São Paulo and Rio de Janeiro states in early March 2020, the epidemic spread to northern states and then to highly populated coastal regions and the Central-West. Municipalities with higher metrics of social development experienced earlier arrival of COVID-19 (decrease of 11·1 days [95% CI:8.9,13.2] in the time to arrival for each 10% increase in the social development index). Differences in the initial epidemic intensity (mean Rt ) were largely driven by geographic location and the date of local onset.

INTERPRETATION

This study demonstrates that platforms that monitor, standardise and analyse the epidemiological data at a local level can give useful real-time insights into outbreak dynamics that can be used to better adapt responses to the current and future pandemics.

FUNDING

This project was supported by a Medical Research Council UK (MRC-UK) -São Paulo Research Foundation (FAPESP) CADDE partnership award (MR/S0195/1 and FAPESP 18/14389-0).

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases 1-3 . The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions 4,5 . Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta's invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta’s invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Global disparities in SARS-CoV-2 genomic surveillance

Genomic sequencing provides critical information to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments and vaccines, and guide public health responses. To investigate the spatiotemporal heterogeneity in the global SARS-CoV-2 genomic surveillance, we estimated the impact of sequencing intensity and turnaround times (TAT) on variant detection in 167 countries. Most countries submit genomes >21 days after sample collection, and 77% of low and middle income countries sequenced <0.5% of their cases. We found that sequencing at least 0.5% of the cases, with a TAT <21 days, could be a benchmark for SARS-CoV-2 genomic surveillance efforts. Socioeconomic inequalities substantially impact our ability to quickly detect SARS-CoV-2 variants, and undermine the global pandemic preparedness.

Interacting Epidemics in Amazonian Brazil: Prior Dengue Infection Associated With Increased Coronavirus Disease 2019 (COVID-19) Risk in a Population-Based Cohort Study

BACKGROUND

Immunity after dengue virus (DENV) infection has been suggested to cross-protect from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and mortality.

METHODS

We tested whether serologically proven prior DENV infection diagnosed in September-October 2019, before the coronavirus disease 2019 (COVID-19) pandemic, reduced the risk of SARS-CoV-2 infection and clinically apparent COVID-19 over the next 13 months in a population-based cohort in Amazonian Brazil. Mixed-effects multiple logistic regression analysis was used to identify predictors of infection and disease, adjusting for potential individual and household-level confounders. Virus genomes from 14 local SARS-CoV-2 isolates were obtained using whole-genome sequencing.

RESULTS

Anti-DENV immunoglobulin G (IgG) was found in 37.0% of 1285 cohort participants (95% confidence interval [CI]: 34.3% to 39.7%) in 2019, with 10.4 (95% CI: 6.7-15.5) seroconversion events per 100 person-years during the follow-up. In 2020, 35.2% of the participants (95% CI: 32.6% to 37.8%) had anti-SARS-CoV-2 IgG and 57.1% of the 448 SARS-CoV-2 seropositives (95% CI: 52.4% to 61.8%) reported clinical manifestations at the time of infection. Participants aged >60 years were twice more likely to have symptomatic COVID-19 than children under 5 years. Locally circulating SARS-CoV-2 isolates were assigned to the B.1.1.33 lineage. Contrary to the cross-protection hypothesis, prior DENV infection was associated with twice the risk of clinically apparent COVID-19 upon SARS-CoV-2 infection, with P values between .025 and .039 after adjustment for identified confounders.

CONCLUSIONS

Higher risk of clinically apparent COVID-19 among individuals with prior dengue has important public health implications for communities sequentially exposed to DENV and SARS-CoV-2 epidemics.

Report 46: Factors driving extensive spatial and temporal fluctuations in COVID-19 fatality rates in Brazilian hospitals

The SARS-CoV-2 Gamma variant spread rapidly across Brazil, causing substantial infection and death waves. We use individual-level patient records following hospitalisation with suspected or confirmed COVID-19 to document the extensive shocks in hospital fatality rates that followed Gamma's spread across 14 state capitals, and in which more than half of hospitalised patients died over sustained time periods. We show that extensive fluctuations in COVID-19 in-hospital fatality rates also existed prior to Gamma's detection, and were largely transient after Gamma's detection, subsiding with hospital demand. Using a Bayesian fatality rate model, we find that the geographic and temporal fluctuations in Brazil's COVID-19 in-hospital fatality rates are primarily associated with geographic inequities and shortages in healthcare capacity. We project that approximately half of Brazil's COVID-19 deaths in hospitals could have been avoided without pre-pandemic geographic inequities and without pandemic healthcare pressure. Our results suggest that investments in healthcare resources, healthcare optimization, and pandemic preparedness are critical to minimize population wide mortality and morbidity caused by highly transmissible and deadly pathogens such as SARS-CoV-2, especially in low- and middle-income countries.

NOTE

The following manuscript has appeared as 'Report 46 - Factors driving extensive spatial and temporal fluctuations in COVID-19 fatality rates in Brazilian hospitals' at https://spiral.imperial.ac.uk:8443/handle/10044/1/91875 .

ONE SENTENCE SUMMARY

COVID-19 in-hospital fatality rates fluctuate dramatically in Brazil, and these fluctuations are primarily associated with geographic inequities and shortages in healthcare capacity.

Neutralisation of SARS-CoV-2 lineage P.1 by antibodies elicited through natural SARS-CoV-2 infection or vaccination with an inactivated SARS-CoV-2 vaccine: an immunological study

BACKGROUND

Mutations accrued by SARS-CoV-2 lineage P.1-first detected in Brazil in early January, 2021-include amino acid changes in the receptor-binding domain of the viral spike protein that also are reported in other variants of concern, including B.1.1.7 and B.1.351. We aimed to investigate whether isolates of wild-type P.1 lineage SARS-CoV-2 can escape from neutralising antibodies generated by a polyclonal immune response.

METHODS

We did an immunological study to assess the neutralising effects of antibodies on lineage P.1 and lineage B isolates of SARS-CoV-2, using plasma samples from patients previously infected with or vaccinated against SARS-CoV-2. Two specimens (P.1/28 and P.1/30) containing SARS-CoV-2 lineage P.1 (as confirmed by viral genome sequencing) were obtained from nasopharyngeal and bronchoalveolar lavage samples collected from patients in Manaus, Brazil, and compared against an isolate of SARS-CoV-2 lineage B (SARS.CoV2/SP02.2020) recovered from a patient in Brazil in February, 2020. Isolates were incubated with plasma samples from 21 blood donors who had previously had COVID-19 and from a total of 53 recipients of the chemically inactivated SARS-CoV-2 vaccine CoronaVac: 18 individuals after receipt of a single dose and an additional 20 individuals (38 in total) after receipt of two doses (collected 17-38 days after the most recent dose); and 15 individuals who received two doses during the phase 3 trial of the vaccine (collected 134-230 days after the second dose). Antibody neutralisation of P.1/28, P.1/30, and B isolates by plasma samples were compared in terms of median virus neutralisation titre (VNT50, defined as the reciprocal value of the sample dilution that showed 50% protection against cytopathic effects).

FINDINGS

In terms of VNT50, plasma from individuals previously infected with SARS-CoV-2 had an 8·6 times lower neutralising capacity against the P.1 isolates (median VNT50 30 [IQR <20-45] for P.1/28 and 30 [<20-40] for P.1/30) than against the lineage B isolate (260 [160-400]), with a binominal model showing significant reductions in lineage P.1 isolates compared with the lineage B isolate (p≤0·0001). Efficient neutralisation of P.1 isolates was not seen with plasma samples collected from individuals vaccinated with a first dose of CoronaVac 20-23 days earlier (VNT50s below the limit of detection [<20] for most plasma samples), a second dose 17-38 days earlier (median VNT50 24 [IQR <20-25] for P.1/28 and 28 [<20-25] for P.1/30), or a second dose 134-260 days earlier (all VNT50s below limit of detection). Median VNT50s against the lineage B isolate were 20 (IQR 20-30) after a first dose of CoronaVac 20-23 days earlier, 75 (<20-263) after a second dose 17-38 days earlier, and 20 (<20-30) after a second dose 134-260 days earlier. In plasma collected 17-38 days after a second dose of CoronaVac, neutralising capacity against both P.1 isolates was significantly decreased (p=0·0051 for P.1/28 and p=0·0336 for P.1/30) compared with that against the lineage B isolate. All data were corroborated by results obtained through plaque reduction neutralisation tests.

INTERPRETATION

SARS-CoV-2 lineage P.1 might escape neutralisation by antibodies generated in response to polyclonal stimulation against previously circulating variants of SARS-CoV-2. Continuous genomic surveillance of SARS-CoV-2 combined with antibody neutralisation assays could help to guide national immunisation programmes.

FUNDING

São Paulo Research Foundation, Brazilian Ministry of Science, Technology and Innovation and Funding Authority for Studies, Medical Research Council, National Council for Scientific and Technological Development, National Institutes of Health.

TRANSLATION

For the Portuguese translation of the abstract see Supplementary Materials section.

Spatiotemporal invasion dynamics of SARS-CoV-2 lineage B.1.1.7 emergence

Understanding the causes and consequences of the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern is crucial to pandemic control yet difficult to achieve because they arise in the context of variable human behavior and immunity. We investigated the spatial invasion dynamics of lineage B.1.1.7 by jointly analyzing UK human mobility, virus genomes, and community-based polymerase chain reaction data. We identified a multistage spatial invasion process in which early B.1.1.7 growth rates were associated with mobility and asymmetric lineage export from a dominant source location, enhancing the effects of B.1.1.7's increased intrinsic transmissibility. We further explored how B.1.1.7 spread was shaped by nonpharmaceutical interventions and spatial variation in previous attack rates. Our findings show that careful accounting of the behavioral and epidemiological context within which variants of concern emerge is necessary to interpret correctly their observed relative growth rates.

Relax, Keep Walking - A Practical Guide to Continuous Phylogeographic Inference with BEAST

Spatially explicit phylogeographic analyses can be performed with an inference framework that employs relaxed random walks to reconstruct phylogenetic dispersal histories in continuous space. This core model was first implemented 10 years ago and has opened up new opportunities in the field of phylodynamics, allowing researchers to map and analyze the spatial dissemination of rapidly evolving pathogens. We here provide a detailed and step-by-step guide on how to set up, run, and interpret continuous phylogeographic analyses using the programs BEAUti, BEAST, Tracer, and TreeAnnotator.

Respiratory Viral Shedding in Healthcare Workers Reinfected with SARS-CoV-2, Brazil, 2020

We documented 4 cases of severe acute respiratory syndrome coronavirus 2 reinfection by non-variant of concern strains among healthcare workers in Campinas, Brazil. We isolated infectious particles from nasopharyngeal secretions during both infection episodes. Improved and continued protection measures are necessary to mitigate the risk for reinfection among healthcare workers.

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

A Phylodynamic Workflow to Rapidly Gain Insights into the Dispersal History and Dynamics of SARS-CoV-2 Lineages

Since the start of the COVID-19 pandemic, an unprecedented number of genomic sequences of SARS-CoV-2 have been generated and shared with the scientific community. The unparalleled volume of available genetic data presents a unique opportunity to gain real-time insights into the virus transmission during the pandemic, but also a daunting computational hurdle if analyzed with gold-standard phylogeographic approaches. To tackle this practical limitation, we here describe and apply a rapid analytical pipeline to analyze the spatiotemporal dispersal history and dynamics of SARS-CoV-2 lineages. As a proof of concept, we focus on the Belgian epidemic, which has had one of the highest spatial densities of available SARS-CoV-2 genomes. Our pipeline has the potential to be quickly applied to other countries or regions, with key benefits in complementing epidemiological analyses in assessing the impact of intervention measures or their progressive easement.

Genomic epidemiology of SARS-CoV-2 transmission lineages in Ecuador

Characterisation of SARS-CoV-2 genetic diversity through space and time can reveal trends in virus importation and domestic circulation, and permit the exploration of questions regarding the early transmission dynamics. Here we present a detailed description of SARS-CoV-2 genomic epidemiology in Ecuador, one of the hardest hit countries during the early stages of the COVID-19 pandemic. We generate and analyse 160 whole genome sequences sampled from all provinces of Ecuador in 2020. Molecular clock and phylgeographic analysis of these sequences in the context of global SARS-CoV-2 diversity enable us to identify and characterise individual transmission lineages within Ecuador, explore their spatiotemporal distributions, and consider their introduction and domestic circulation. Our results reveal a pattern of multiple international importations across the country, with apparent differences between key provinces. Transmission lineages were mostly introduced before the implementation of non-pharmaceutical interventions (NPIs), with differential degrees of persistence and national dissemination.

Higher risk of death from COVID-19 in low-income and non-White populations of São Paulo, Brazil

INTRODUCTION

Little evidence exists on the differential health effects of COVID-19 on disadvantaged population groups. Here we characterise the differential risk of hospitalisation and death in São Paulo state, Brazil, and show how vulnerability to COVID-19 is shaped by socioeconomic inequalities.

METHODS

We conducted a cross-sectional study using hospitalised severe acute respiratory infections notified from March to August 2020 in the Sistema de Monitoramento Inteligente de São Paulo database. We examined the risk of hospitalisation and death by race and socioeconomic status using multiple data sets for individual-level and spatiotemporal analyses. We explained these inequalities according to differences in daily mobility from mobile phone data, teleworking behaviour and comorbidities.

RESULTS

Throughout the study period, patients living in the 40% poorest areas were more likely to die when compared with patients living in the 5% wealthiest areas (OR: 1.60, 95% CI 1.48 to 1.74) and were more likely to be hospitalised between April and July 2020 (OR: 1.08, 95% CI 1.04 to 1.12). Black and Pardo individuals were more likely to be hospitalised when compared with White individuals (OR: 1.41, 95% CI 1.37 to 1.46; OR: 1.26, 95% CI 1.23 to 1.28, respectively), and were more likely to die (OR: 1.13, 95% CI 1.07 to 1.19; 1.07, 95% CI 1.04 to 1.10, respectively) between April and July 2020. Once hospitalised, patients treated in public hospitals were more likely to die than patients in private hospitals (OR: 1.40%, 95% CI 1.34% to 1.46%). Black individuals and those with low education attainment were more likely to have one or more comorbidities, respectively (OR: 1.29, 95% CI 1.19 to 1.39; 1.36, 95% CI 1.27 to 1.45).

CONCLUSIONS

Low-income and Black and Pardo communities are more likely to die with COVID-19. This is associated with differential access to quality healthcare, ability to self-isolate and the higher prevalence of comorbidities.

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.