Genomic epidemiology of West Nile virus in Europe

West Nile virus is one of the most widespread mosquito-borne zoonotic viruses, with unique transmission dynamics in various parts of the world. Genomic surveillance has provided important insights in the global patterns of West Nile virus emergence and spread. In Europe, multiple West Nile virus lineages have been isolated, with lineage 1a and 2 being the main lineages responsible for human infections. In contrast to North America, where a single introduction of lineage 1a resulted in the virus establishing itself in a new continent, at least 13 introductions of lineages 1a and 2 have occurred into Europe, which is likely a vast underestimation of the true number of introductions. Historically, lineage 1a was the main lineage circulating in Europe, but since the emergence of lineage 2 in the early 2000s, the latter has become the predominant lineage. This shift in West Nile virus lineage prevalence has been broadly linked to the expansion of the virus into northerly temperate regions, where autochthonous cases in animals and humans have been reported in Germany and The Netherlands. Here, we discuss how genomic analysis has increased our understanding of the epidemiology of West Nile virus in Europe, and we present a global Nextstrain build consisting of publicly available West Nile virus genomes (https://nextstrain.org/community/grubaughlab/WNV-Global). Our results elucidate recent insights in West Nile virus lineage dynamics in Europe, and discuss how expanded programs can fill current genomic surveillance gaps.

DengueSeq: a pan-serotype whole genome amplicon sequencing protocol for dengue virus

BACKGROUND

The increasing burden of dengue virus on public health due to more explosive and frequent outbreaks highlights the need for improved surveillance and control. Genomic surveillance of dengue virus not only provides important insights into the emergence and spread of genetically diverse serotypes and genotypes, but it is also critical to monitor the effectiveness of newly implemented control strategies. Here, we present DengueSeq, an amplicon sequencing protocol, which enables whole-genome sequencing of all four dengue virus serotypes.

RESULTS

We developed primer schemes for the four dengue virus serotypes, which can be combined into a pan-serotype approach. We validated both approaches using genetically diverse virus stocks and clinical specimens that contained a range of virus copies. High genome coverage (>95%) was achieved for all genotypes, except DENV2 (genotype VI) and DENV 4 (genotype IV) sylvatics, with similar performance of the serotype-specific and pan-serotype approaches. The limit of detection to reach 70% coverage was 10-100 RNA copies/μL for all four serotypes, which is similar to other commonly used primer schemes. DengueSeq facilitates the sequencing of samples without known serotypes, allows the detection of multiple serotypes in the same sample, and can be used with a variety of library prep kits and sequencing instruments.

CONCLUSIONS

DengueSeq was systematically evaluated with virus stocks and clinical specimens spanning the genetic diversity within each of the four dengue virus serotypes. The primer schemes can be plugged into existing amplicon sequencing workflows to facilitate the global need for expanded dengue virus genomic surveillance.

Travel surveillance uncovers dengue virus dynamics and introductions in the Caribbean

Dengue is the most prevalent mosquito-borne viral disease in humans, and cases are continuing to rise globally. In particular, islands in the Caribbean have experienced more frequent outbreaks, and all four dengue virus (DENV) serotypes have been reported in the region, leading to hyperendemicity and increased rates of severe disease. However, there is significant variability regarding virus surveillance and reporting between islands, making it difficult to obtain an accurate understanding of the epidemiological patterns in the Caribbean. To investigate this, we used travel surveillance and genomic epidemiology to reconstruct outbreak dynamics, DENV serotype turnover, and patterns of spread within the region from 2009-2022. We uncovered two recent DENV-3 introductions from Asia, one of which resulted in a large outbreak in Cuba, which was previously under-reported. We also show that while outbreaks can be synchronized between islands, they are often caused by different serotypes. Our study highlights the importance of surveillance of infected travelers to provide a snapshot of local introductions and transmission in areas with limited local surveillance and suggests that the recent DENV-3 introductions may pose a major public health threat in the region.

Contribution of climate change to the spatial expansion of West Nile virus in Europe

West Nile virus (WNV) is an emerging mosquito-borne pathogen in Europe where it represents a new public health threat. While climate change has been cited as a potential driver of its spatial expansion on the continent, a formal evaluation of this causal relationship is lacking. Here, we investigate the extent to which WNV spatial expansion in Europe can be attributed to climate change while accounting for other direct human influences such as land-use and human population changes. To this end, we trained ecological niche models to predict the risk of local WNV circulation leading to human cases to then unravel the isolated effect of climate change by comparing factual simulations to a counterfactual based on the same environmental changes but a counterfactual climate where long-term trends have been removed. Our findings demonstrate a notable increase in the area ecologically suitable for WNV circulation during the period 1901-2019, whereas this area remains largely unchanged in a no-climate-change counterfactual. We show that the drastic increase in the human population at risk of exposure is partly due to historical changes in population density, but that climate change has also been a critical driver behind the heightened risk of WNV circulation in Europe.

Introduction and Spread of Dengue Virus 3, Florida, USA, May 2022-April 2023

During May 2022-April 2023, dengue virus serotype 3 was identified among 601 travel-associated and 61 locally acquired dengue cases in Florida, USA. All 203 sequenced genomes belonged to the same genotype III lineage and revealed potential transmission chains in which most locally acquired cases occurred shortly after introduction, with little sustained transmission.

Identification and characterization of novel lineage 1 Powassan virus strains in New York State

Powassan virus (POWV, family Flaviviridae) is a reemerging tick-borne virus endemic in North America and Russia. In 1997, a POWV-like agent was isolated from Ixodes scapularis in New England and determined to be genetically distinct from the original POWV isolate. This revealed the existence of two lineages: lineage 1, prototype Powassan virus (POWV-1) and lineage 2, deer tick virus (DTV). POWV-1 is thought to be primarily maintained in a cycle between I. cookei and woodchucks and I. marxi and squirrels, while DTV is primarily maintained in a cycle between I. scapularis and small mammal hosts. Recent tick, mammalian, and human isolates from New York State (NYS) have been identified as DTV, but for the first time in 45 years, we detected four POWV-1 isolates, including the first reported isolation of POWV-1 from I. scapularis. We aimed to investigate genotypic and phenotypic characteristics of recent NYS isolates through sequence analysis and evaluation of replication kinetics in vitro and in vivo. Our sequencing revealed genetic divergence between NYS POWV-1 isolates, with two distinct foci. We found that POWV-1 isolates displayed variable replication kinetics in nymphal ticks but not in cell culture. POWV-1 isolated from I. scapularis displayed increased fitness in experimentally infected I. scapularis as compared to historic and recent POWV-1 isolates from I. cookei. These data suggest the emergence of divergent POWV-1 strains in alternate tick hosts and maintenance of genetically and phenotypically discrete POWV-1 foci.

The known unknowns of Powassan virus ecology

Powassan virus (POWV; Family: Flaviviridae, Genus: Flavivirus) is the sole North American member of the tick-borne encephalitis sero-complex. While associated with high rates of morbidity and mortality, POWV has historically been of little public health concern due to low incidence rates. However, over the last 20 yr, incidence rates have increased highlighting the growing epidemiological threat. Currently, there are no vaccines or therapeutics with tick habitat reduction, acaricide application, and public awareness programs being our primary means of intervention. The effectiveness of these control strategies is dependent on having a sound understanding of the virus's ecology. In this Forum, we review what is currently known about POWV ecology, identify gaps in our knowledge, and discuss prevailing and alternative hypotheses about transmission dynamics, reservoir hosts, and spatial focality.

Survey of white-footed mice in Connecticut, USA reveals low SARS-CoV-2 seroprevalence and infection with divergent betacoronaviruses

Diverse mammalian species display susceptibility to and infection with SARS-CoV-2. Potential SARS-CoV-2 spillback into rodents is understudied despite their host role for numerous zoonoses and human proximity. We assessed exposure and infection among white-footed mice (Peromyscus leucopus) in Connecticut, USA. We observed 1% (6/540) wild-type neutralizing antibody seroprevalence among 2020-2022 residential mice with no cross-neutralization of variants. We detected no SARS-CoV-2 infections via RT-qPCR, but identified non-SARS-CoV-2 betacoronavirus infections via pan-coronavirus PCR among 1% (5/468) of residential mice. Sequencing revealed two divergent betacoronaviruses, preliminarily named Peromyscus coronavirus-1 and -2. Both belong to the Betacoronavirus 1 species and are ~90% identical to the closest known relative, Porcine hemagglutinating encephalomyelitis virus. Low SARS-CoV-2 seroprevalence suggests white-footed mice may not be sufficiently susceptible or exposed to SARS-CoV-2 to present a long-term human health risk. However, the discovery of divergent, non-SARS-CoV-2 betacoronaviruses expands the diversity of known rodent coronaviruses and further investigation is required to understand their transmission extent.

Toward a global virus genomic surveillance network

The COVID-19 pandemic galvanized the field of virus genomic surveillance, demonstrating its utility for public health. Now, we must harness the momentum that led to increased infrastructure, training, and political will to build a sustainable global genomic surveillance network for other epidemic and endemic viruses. We suggest a generalizable modular sequencing framework wherein users can easily switch between virus targets to maximize cost-effectiveness and maintain readiness for new threats. We also highlight challenges associated with genomic surveillance and when global inequalities persist. We propose solutions to mitigate some of these issues, including training and multilateral partnerships. Exploring alternatives to clinical sequencing can also reduce the cost of surveillance programs. Finally, we discuss how establishing genomic surveillance would aid control programs and potentially provide a warning system for outbreaks, using a global respiratory virus (RSV), an arbovirus (dengue virus), and a regional zoonotic virus (Lassa virus) as examples.

Development of an amplicon-based sequencing approach in response to the global emergence of mpox

The 2022 multicountry mpox outbreak concurrent with the ongoing Coronavirus Disease 2019 (COVID-19) pandemic further highlighted the need for genomic surveillance and rapid pathogen whole-genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical specimens that tested presumptively positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (Ct) (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR Ct below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon-based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole-genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Dynamics of eastern equine encephalitis virus during the 2019 outbreak in the Northeast United States

Eastern equine encephalitis virus (EEEV) causes a rare but severe disease in horses and humans and is maintained in an enzootic transmission cycle between songbirds and Culiseta melanura mosquitoes. In 2019, the largest EEEV outbreak in the United States for more than 50 years occurred, centered in the Northeast. To explore the dynamics of the outbreak, we sequenced 80 isolates of EEEV and combined them with existing genomic data. We found that, similar to previous years, cases were driven by multiple independent but short-lived virus introductions into the Northeast from Florida. Once in the Northeast, we found that Massachusetts was important for regional spread. We found no evidence of any changes in viral, human, or bird factors which would explain the increase in cases in 2019, although the ecology of EEEV is complex and further data is required to explore these in more detail. By using detailed mosquito surveillance data collected by Massachusetts and Connecticut, however, we found that the abundance of Cs. melanura was exceptionally high in 2019, as was the EEEV infection rate. We employed these mosquito data to build a negative binomial regression model and applied it to estimate early season risks of human or horse cases. We found that the month of first detection of EEEV in mosquito surveillance data and vector index (abundance multiplied by infection rate) were predictive of cases later in the season. We therefore highlight the importance of mosquito surveillance programs as an integral part of public health and disease control.

Nonsystematic Reporting Biases of the SARS-CoV-2 Variant Mu Could Impact Our Understanding of the Epidemiological Dynamics of Emerging Variants

Developing a timely and effective response to emerging SARS-CoV-2 variants of concern (VOCs) is of paramount public health importance. Global health surveillance does not rely on genomic data alone to identify concerning variants when they emerge. Instead, methods that utilize genomic data to estimate the epidemiological dynamics of emerging lineages have the potential to serve as an early warning system. However, these methods assume that genomic data are uniformly reported across circulating lineages. In this study, we analyze differences in reporting delays among SARS-CoV-2 VOCs as a plausible explanation for the timing of the global response to the former VOC Mu. Mu likely emerged in South America in mid-2020, where its circulation was largely confined. In this study, we demonstrate that Mu was designated as a VOC ∼1 year after it emerged and find that the reporting of genomic data for Mu differed significantly than that of other VOCs within countries, states, and individual laboratories. Our findings suggest that nonsystematic biases in the reporting of genomic data may have delayed the global response to Mu. Until they are resolved, the surveillance gaps that affected the global response to Mu could impede the rapid and accurate assessment of future emerging variants.

Dynamics of Eastern equine encephalitis virus during the 2019 outbreak in the Northeast United States

Eastern equine encephalitis virus (EEEV) causes a rare but severe disease in horses and humans, and is maintained in an enzootic transmission cycle between songbirds and Culiseta melanura mosquitoes. In 2019, the largest EEEV outbreak in the United States for more than 50 years occurred, centered in the Northeast. To explore the dynamics of the outbreak, we sequenced 80 isolates of EEEV and combined them with existing genomic data. We found that, like previous years, cases were driven by frequent short-lived virus introductions into the Northeast from Florida. Once in the Northeast, we found that Massachusetts was important for regional spread. We found no evidence of any changes in viral, human, or bird factors which would explain the increase in cases in 2019. By using detailed mosquito surveillance data collected by Massachusetts and Connecticut, however, we found that the abundance of Cs. melanura was exceptionally high in 2019, as was the EEEV infection rate. We employed these mosquito data to build a negative binomial regression model and applied it to estimate early season risks of human or horse cases. We found that the month of first detection of EEEV in mosquito surveillance data and vector index (abundance multiplied by infection rate) were predictive of cases later in the season. We therefore highlight the importance of mosquito surveillance programs as an integral part of public health and disease control.

Age-dependent impairment in antibody responses elicited by a homologous CoronaVac booster dose

The emergence of the SARS-CoV-2 Omicron sublineages resulted in increased transmission rates and reduced protection from vaccines. To counteract these effects, multiple booster strategies were used in different countries, although data comparing their efficiency in improving protective immunity remain sparse, especially among vulnerable populations, including older adults. The inactivated CoronaVac vaccine was among the most widely distributed vaccine worldwide and was essential in the early control of SARS-CoV-2-related hospitalizations and deaths. However, it is not well understood whether homologous versus heterologous booster doses in those fully vaccinated with CoronaVac induce distinct humoral responses or whether these responses vary across age groups. We analyzed plasma antibody responses from CoronaVac-vaccinated younger or older individuals who received a homologous CoronaVac or heterologous BNT162b2 or ChAdOx1 booster vaccine. All three evaluated boosters resulted in increased virus-specific IgG titers 28 days after the booster dose. However, we found that both IgG titers against SARS-CoV-2 Spike or RBD and neutralization titers against Omicron sublineages were substantially reduced in participants who received homologous CoronaVac compared with the heterologous BNT162b2 or ChAdOx1 booster. This effect was specifically prominent in recipients >50 years of age. In this group, the CoronaVac booster induced low virus-specific IgG titers and failed to elevate neutralization titers against any Omicron sublineage. Our results point to the notable inefficiency of CoronaVac immunization and boosting in mounting protective antiviral humoral immunity, particularly among older adults, during the Omicron wave. These observations also point to benefits of heterologous regimens in high-risk populations fully vaccinated with CoronaVac.

Genomic and phenotypic analyses suggest moderate fitness differences among Zika virus lineages

RNA viruses have short generation times and high mutation rates, allowing them to undergo rapid molecular evolution during epidemics. However, the extent of RNA virus phenotypic evolution within epidemics and the resulting effects on fitness and virulence remain mostly unknown. Here, we screened the 2015-2016 Zika epidemic in the Americas for lineage-specific fitness differences. We engineered a library of recombinant viruses representing twelve major Zika virus lineages and used them to measure replicative fitness within disease-relevant human primary cells and live mosquitoes. We found that two of these lineages conferred significant in vitro replicative fitness changes among human primary cells, but we did not find fitness changes in Aedes aegypti mosquitoes. Additionally, we found evidence for elevated levels of positive selection among five amino acid sites that define major Zika virus lineages. While our work suggests that Zika virus may have acquired several phenotypic changes during a short time scale, these changes were relatively moderate and do not appear to have enhanced transmission during the epidemic.

Powassan Virus Lineage I in Field-Collected Dermacentor variabilis Ticks, New York, USA

Powassan virus is a tickborne flavivirus that can cause lethal or debilitating neurologic illness. It is canonically transmitted by Ixodes spp. ticks but might spill over to sympatric Dermacentor species. We detected Powassan virus lineage I from a pool of field-collected D. variabilis ticks in New York, USA.

Nasal host response-based screening for undiagnosed respiratory viruses: a pathogen surveillance and detection study

BACKGROUND

Symptomatic patients who test negative for common viruses are an important possible source of unrecognised or emerging pathogens, but metagenomic sequencing of all samples is inefficient because of the low likelihood of finding a pathogen in any given sample. We aimed to determine whether nasopharyngeal CXCL10 screening could be used as a strategy to enrich for samples containing undiagnosed viruses.

METHODS

In this pathogen surveillance and detection study, we measured CXCL10 concentrations from nasopharyngeal swabs from patients in the Yale New Haven health-care system, which had been tested at the Yale New Haven Hospital Clinical Virology Laboratory (New Haven, CT, USA). Patients who tested negative for a panel of respiratory viruses using multiplex PCR during Jan 23-29, 2017, or March 3-14, 2020, were included. We performed host and pathogen RNA sequencing (RNA-Seq) and analysis for viral reads on samples with CXCL10 higher than 1 ng/mL or CXCL10 testing and quantitative RT-PCR (RT-qPCR) for SARS-CoV-2. We used RNA-Seq and cytokine profiling to compare the host response to infection in samples that were virus positive (rhinovirus, seasonal coronavirus CoV-NL63, or SARS-CoV-2) and virus negative (controls).

FINDINGS

During Jan 23-29, 2017, 359 samples were tested for ten viruses on the multiplex PCR respiratory virus panel (RVP). 251 (70%) were RVP negative. 60 (24%) of 251 samples had CXCL10 higher than 150 pg/mL and were identified for further analysis. 28 (47%) of 60 CXCL10-high samples were positive for seasonal coronaviruses. 223 (89%) of 251 samples were PCR negative for 15 viruses and, of these, CXCL10-based screening identified 32 (13%) samples for further analysis. Of these 32 samples, eight (25%) with CXCL10 concentrations higher than 1 ng/mL and sufficient RNA were selected for RNA-Seq. Microbial RNA analysis showed the presence of influenza C virus in one sample and revealed RNA reads from bacterial pathobionts in four (50%) of eight samples. Between March 3 and March 14, 2020, 375 (59%) of 641 samples tested negative for 15 viruses on the RVP. 32 (9%) of 375 samples had CXCL10 concentrations ranging from 100 pg/mL to 1000 pg/mL and four of those were positive for SARS-CoV-2. CXCL10 elevation was statistically significant, and a distinguishing feature was found in 28 (8%) of 375 SARS-CoV-2-negative samples versus all four SARS-CoV-2-positive samples (p=4·4 × 10-5). Transcriptomic signatures showed an interferon response in virus-positive samples and an additional neutrophil-high hyperinflammatory signature in samples with high amounts of bacterial pathobionts. The CXCL10 cutoff for detecting a virus was 166·5 pg/mL for optimal sensitivity and 1091·0 pg/mL for specificity using a clinic-ready automated microfluidics-based immunoassay.

INTERPRETATION

These results confirm CXCL10 as a robust nasopharyngeal biomarker of viral respiratory infection and support host response-based screening followed by metagenomic sequencing of CXCL10-high samples as a practical approach to incorporate clinical samples into pathogen discovery and surveillance efforts.

FUNDING

National Institutes of Health, the Hartwell Foundation, the Gruber Foundation, Fast Grants for COVID-19 research from the Mercatus Center, and the Huffman Family Donor Advised Fund.

Lineage abundance estimation for SARS-CoV-2 in wastewater using transcriptome quantification techniques

Effectively monitoring the spread of SARS-CoV-2 mutants is essential to efforts to counter the ongoing pandemic. Predicting lineage abundance from wastewater, however, is technically challenging. We show that by sequencing SARS-CoV-2 RNA in wastewater and applying algorithms initially used for transcriptome quantification, we can estimate lineage abundance in wastewater samples. We find high variability in signal among individual samples, but the overall trends match those observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in mutant prevalence in situations where clinical sequencing is unavailable.

Interactions between seasonal temperature variation and temporal synchrony drive increased arbovirus co-infection incidence

Though instances of arthropod-borne (arbo)virus co-infection have been documented clinically, the overall incidence of arbovirus co-infection and its drivers are not well understood. Now that dengue, Zika and chikungunya viruses are all in circulation across tropical and subtropical regions of the Americas, it is important to understand the environmental and biological conditions that make co-infections more likely to occur. To understand this, we developed a mathematical model of co-circulation of two arboviruses, with transmission parameters approximating dengue, Zika and/or chikungunya viruses, and co-infection possible in both humans and mosquitoes. We examined the influence of seasonal timing of arbovirus co-circulation on the extent of co-infection. By undertaking a sensitivity analysis of this model, we examined how biological factors interact with seasonality to determine arbovirus co-infection transmission and prevalence. We found that temporal synchrony of the co-infecting viruses and average temperature were the most influential drivers of co-infection incidence. Our model highlights the synergistic effect of co-transmission from mosquitoes, which leads to more than double the number of co-infections than would be expected in a scenario without co-transmission. Our results suggest that appreciable numbers of co-infections are unlikely to occur except in tropical climates when the viruses co-occur in time and space.

Rapid emergence of SARS-CoV-2 Omicron variant is associated with an infection advantage over Delta in vaccinated persons

Background

The SARS-CoV-2 Omicron variant became a global concern due to its rapid spread and displacement of the dominant Delta variant. We hypothesized that part of Omicron's rapid rise was based on its increased ability to cause infections in persons that are vaccinated compared to Delta.

Methods

We analyzed nasal swab PCR tests for samples collected between December 12 and 16, 2021, in Connecticut when the proportion of Delta and Omicron variants was relatively equal. We used the spike gene target failure (SGTF) to classify probable Delta and Omicron infections. We fitted an exponential curve to the estimated infections to determine the doubling times for each variant. We compared the test positivity rates for each variant by vaccination status, number of doses, and vaccine manufacturer. Generalized linear models were used to assess factors associated with odds of infection with each variant among persons testing positive for SARS-CoV-2.

Findings

For infections with high virus copies (Ct < 30) among vaccinated persons, we found higher odds that they were infected with Omicron compared to Delta, and that the odds increased with increased number of vaccine doses. Compared to unvaccinated persons, we found significant reduction in Delta positivity rates after two (43.4%-49.1%) and three vaccine doses (81.1%), while we only found a significant reduction in Omicron positivity rates after three doses (62.3%).

Conclusion

The rapid rise in Omicron infections was likely driven by Omicron's escape from vaccine-induced immunity.

Funding

This work was supported by the Centers for Disease Control and Prevention (CDC).

Rapid emergence of SARS-CoV-2 Omicron variant is associated with an infection advantage over Delta in vaccinated persons

BACKGROUND

The SARS-CoV-2 Omicron variant became a global concern due to its rapid spread and displacement of the dominant Delta variant. We hypothesized that part of Omicron's rapid rise was based on its increased ability to cause infections in persons that are vaccinated compared to Delta.

METHODS

We analyzed nasal swab PCR tests for samples collected between December 12 and 16, 2021, in Connecticut when the proportion of Delta and Omicron variants was relatively equal. We used the spike gene target failure (SGTF) to classify probable Delta and Omicron infections. We fitted an exponential curve to the estimated infections to determine the doubling times for each variant. We compared the test positivity rates for each variant by vaccination status, number of doses, and vaccine manufacturer. Generalized linear models were used to assess factors associated with odds of infection with each variant among persons testing positive for SARS-CoV-2.

FINDINGS

For infections with high virus copies (Ct < 30) among vaccinated persons, we found higher odds that they were infected with Omicron compared to Delta, and that the odds increased with increased number of vaccine doses. Compared to unvaccinated persons, we found significant reduction in Delta positivity rates after two (43.4%-49.1%) and three vaccine doses (81.1%), while we only found a significant reduction in Omicron positivity rates after three doses (62.3%).

CONCLUSION

The rapid rise in Omicron infections was likely driven by Omicron's escape from vaccine-induced immunity.

FUNDING

This work was supported by the Centers for Disease Control and Prevention (CDC).

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA

The SARS-CoV-2 Delta variant rose to dominance in mid-2021, likely propelled by an estimated 40%-80% increased transmissibility over Alpha. To investigate if this ostensible difference in transmissibility is uniform across populations, we partner with public health programs from all six states in New England in the United States. We compare logistic growth rates during each variant's respective emergence period, finding that Delta emerged 1.37-2.63 times faster than Alpha (range across states). We compute variant-specific effective reproductive numbers, estimating that Delta is 63%-167% more transmissible than Alpha (range across states). Finally, we estimate that Delta infections generate on average 6.2 (95% CI 3.1-10.9) times more viral RNA copies per milliliter than Alpha infections during their respective emergence. Overall, our evidence suggests that Delta's enhanced transmissibility can be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on underlying population attributes and sequencing data availability.

Insights into the limited global spread of the immune evasive SARS-CoV-2 variant Mu

SARS-CoV-2 'Variants of Concern' (VOCs) continue to reshape the trajectory of the COVID-19 pandemic. However, why some VOCs, like Omicron, become globally dominant while the spread of others is limited is not fully understood. To address this question, we investigated the VOC Mu, which was first identified in Colombia in late 2020. Our study demonstrates that, although Mu is less sensitive to neutralization compared to variants that preceded it, it did not spread significantly outside of South and Central America. Additionally, we find evidence that the response to Mu was impeded by reporting delays and gaps in the global genomic surveillance system. Our findings suggest that immune evasion alone was not sufficient to outcompete highly transmissible variants that were circulating concurrently with Mu. Insights into the complex relationship between genomic and epidemiological characteristics of previous variants should inform our response to variants that are likely to emerge in the future.

Neutralizing antibodies against the SARS-CoV-2 Delta and Omicron variants following heterologous CoronaVac plus BNT162b2 booster vaccination

The recent emergence of the SARS-CoV-2 Omicron variant is raising concerns because of its increased transmissibility and its numerous spike mutations, which have the potential to evade neutralizing antibodies elicited by COVID-19 vaccines. Here we evaluated the effects of a heterologous BNT162b2 mRNA vaccine booster on the humoral immunity of participants who had received a two-dose regimen of CoronaVac, an inactivated vaccine used globally. We found that a heterologous CoronaVac prime vaccination of two doses followed by a BNT162b2 booster induces elevated virus-specific antibody levels and potent neutralization activity against the ancestral virus and the Delta variant, resembling the titers obtained after two doses of mRNA vaccines. Although neutralization of Omicron was undetectable in participants who had received a two-dose regimen of CoronaVac, the BNT162b2 booster resulted in a 1.4-fold increase in neutralization activity against Omicron compared with the two-dose mRNA vaccine. Despite this increase, neutralizing antibody titers were reduced by 7.1-fold and 3.6-fold for Omicron compared with the ancestral strain and the Delta variant, respectively. These findings have immediate implications for multiple countries that previously used a CoronaVac regimen and reinforce the idea that the Omicron variant is associated with immune escape from vaccines or infection-induced immunity, highlighting the global need for vaccine boosters to combat the impact of emerging variants.

Single-cell multi-omics reveals dyssynchrony of the innate and adaptive immune system in progressive COVID-19

Dysregulated immune responses against the SARS-CoV-2 virus are instrumental in severe COVID-19. However, the immune signatures associated with immunopathology are poorly understood. Here we use multi-omics single-cell analysis to probe the dynamic immune responses in hospitalized patients with stable or progressive course of COVID-19, explore V(D)J repertoires, and assess the cellular effects of tocilizumab. Coordinated profiling of gene expression and cell lineage protein markers shows that S100Ahi/HLA-DRlo classical monocytes and activated LAG-3hi T cells are hallmarks of progressive disease and highlights the abnormal MHC-II/LAG-3 interaction on myeloid and T cells, respectively. We also find skewed T cell receptor repertories in expanded effector CD8+ clones, unmutated IGHG+ B cell clones, and mutated B cell clones with stable somatic hypermutation frequency over time. In conclusion, our in-depth immune profiling reveals dyssynchrony of the innate and adaptive immune interaction in progressive COVID-19.

Sequencing SARS-CoV-2 Genomes from Saliva

Genomic sequencing is crucial to understanding the epidemiology and evolution of SARS-CoV-2. Often, genomic studies rely on remnant diagnostic material, typically nasopharyngeal swabs, as input into whole genome SARS-CoV-2 next-generation sequencing pipelines. Saliva has proven to be a safe and stable specimen for the detection of SARS-CoV-2 RNA via traditional diagnostic assays, however saliva is not commonly used for SARS-CoV-2 sequencing. Using the ARTIC Network amplicon-generation approach with sequencing on the Oxford Nanopore MinION, we demonstrate that sequencing SARS-CoV-2 from saliva produces genomes comparable to those from nasopharyngeal swabs, and that RNA extraction is necessary to generate complete genomes from saliva. In this study, we show that saliva is a useful specimen type for genomic studies of SARS-CoV-2.

Impact of circulating SARS-CoV-2 variants on mRNA vaccine-induced immunity

The emergence of SARS-CoV-2 variants with mutations in major neutralizing antibody-binding sites can affect humoral immunity induced by infection or vaccination1-6. Here we analysed the development of anti-SARS-CoV-2 antibody and T cell responses in individuals who were previously infected (recovered) or uninfected (naive) and received mRNA vaccines to SARS-CoV-2. While individuals who were previously infected sustained higher antibody titres than individuals who were uninfected post-vaccination, the latter reached comparable levels of neutralization responses to the ancestral strain after the second vaccine dose. T cell activation markers measured upon spike or nucleocapsid peptide in vitro stimulation showed a progressive increase after vaccination. Comprehensive analysis of plasma neutralization using 16 authentic isolates of distinct locally circulating SARS-CoV-2 variants revealed a range of reduction in the neutralization capacity associated with specific mutations in the spike gene: lineages with E484K and N501Y/T (for example, B.1.351 and P.1) had the greatest reduction, followed by lineages with L452R (for example, B.1.617.2). While both groups retained neutralization capacity against all variants, plasma from individuals who were previously infected and vaccinated displayed overall better neutralization capacity than plasma from individuals who were uninfected and also received two vaccine doses, pointing to vaccine boosters as a relevant future strategy to alleviate the effect of emerging variants on antibody neutralizing activity.

Longitudinal immune profiling of a SARS-CoV-2 reinfection in a solid organ transplant recipient

BACKGROUND

The underlying immunologic deficiencies enabling SARS-CoV-2 reinfection are currently unknown. We describe deep longitudinal immune profiling of a transplant recipient hospitalized twice for COVID-19.

METHODS

A 66-year-old male renal transplant recipient was hospitalized with COVID-19 March 2020 then readmitted to the hospital with COVID-19 233 days after initial diagnosis. Virologic and immunologic investigation were performed on samples from the primary and secondary infections.

RESULTS

Whole viral genome sequencing and phylogenic analysis revealed that viruses causing both infections were caused by distinct genetic lineages without evidence of immune escape mutations. Longitudinal comparison of cellular and humoral responses during primary SARS-CoV-2 infection revealed that this patient responded to the primary infection with low neutralization titer anti-SARS-CoV-2 antibodies that were likely present at the time of reinfection.

DISCUSSION

The development of neutralizing antibodies and humoral memory responses in this patient failed to confer protection against reinfection, suggesting that they were below a neutralizing titer threshold or that additional factors may be required for efficient prevention of SARS-CoV-2 reinfection. Development of poorly neutralizing antibodies may have been due to profound and relatively specific reduction in naïve CD4 T-cell pools. Seropositivity alone may not be a perfect correlate of protection in immunocompromised patients.

COVID-19 Outcomes and Genomic Characterization of SARS-CoV-2 Isolated From Veterans in New England States: Retrospective Analysis

BACKGROUND

Clinical and virologic characteristics of COVID-19 infections in veterans in New England have not been described. The average US veteran is a male older than the general US population. SARS-CoV-2 infection is known to cause poorer outcomes among men and older adults, making the veteran population an especially vulnerable group for COVID-19.

OBJECTIVE

This study aims to evaluate clinical and virologic factors impacting COVID-19 outcomes.

METHODS

This retrospective chart review included 476 veterans in six New England states with confirmed SARS-CoV-2 infection between April and September 2020. Whole genome sequencing was performed on SARS-CoV-2 RNA isolated from these veterans, and the correlation of genomic data to clinical outcomes was evaluated. Clinical and demographic variables were collected by manual chart review and were correlated to the end points of peak disease severity (based on oxygenation requirements), hospitalization, and mortality using multivariate regression analyses.

RESULTS

Of 476 veterans, 274 had complete and accessible charts. Of the 274 veterans, 92.7% (n=254) were men and 83.2% (n=228) were White, and the mean age was 63 years. In the multivariate regression, significant predictors of hospitalization (C statistic 0.75) were age (odds ratio [OR] 1.05, 95% CI 1.03-1.08) and non-White race (OR 2.39, 95% CI 1.13-5.01). Peak severity (C statistic 0.70) also varied by age (OR 1.07, 95% CI 1.03-1.11) and O2 requirement on admission (OR 45.7, 95% CI 18.79-111). Mortality (C statistic 0.87) was predicted by age (OR 1.06, 95% CI 1.01-1.11), dementia (OR 3.44, 95% CI 1.07-11.1), and O2 requirement on admission (OR 6.74, 95% CI 1.74-26.1). Most (291/299, 97.3%) of our samples were dominated by the spike protein D614G substitution and were from SARS-CoV-2 B.1 lineage or one of 37 different B.1 sublineages, with none representing more than 8.7% (26/299) of the cases.

CONCLUSIONS

In a cohort of veterans from the six New England states with a mean age of 63 years and a high comorbidity burden, age was the largest predictor of hospitalization, peak disease severity, and mortality. Non-White veterans were more likely to be hospitalized, and patients who required oxygen on admission were more likely to have severe disease and higher rates of mortality. Multiple SARS-CoV-2 lineages were distributed in patients in New England early in the COVID-19 era, mostly related to viruses from New York State with D614G mutation.

Variant abundance estimation for SARS-CoV-2 in wastewater using RNA-Seq quantification

Effectively monitoring the spread of SARS-CoV-2 variants is essential to efforts to counter the ongoing pandemic. Wastewater monitoring of SARS-CoV-2 RNA has proven an effective and efficient technique to approximate COVID-19 case rates in the population. Predicting variant abundances from wastewater, however, is technically challenging. Here we show that by sequencing SARS-CoV-2 RNA in wastewater and applying computational techniques initially used for RNA-Seq quantification, we can estimate the abundance of variants in wastewater samples. We show by sequencing samples from wastewater and clinical isolates in Connecticut U.S.A. between January and April 2021 that the temporal dynamics of variant strains broadly correspond. We further show that this technique can be used with other wastewater sequencing techniques by expanding to samples taken across the United States in a similar timeframe. We find high variability in signal among individual samples, and limited ability to detect the presence of variants with clinical frequencies <10%; nevertheless, the overall trends match what we observed from sequencing clinical samples. Thus, while clinical sequencing remains a more sensitive technique for population surveillance, wastewater sequencing can be used to monitor trends in variant prevalence in situations where clinical sequencing is unavailable or impractical.

Viral dynamics of acute SARS-CoV-2 infection and applications to diagnostic and public health strategies

SARS-CoV-2 infections are characterized by viral proliferation and clearance phases and can be followed by low-level persistent viral RNA shedding. The dynamics of viral RNA concentration, particularly in the early stages of infection, can inform clinical measures and interventions such as test-based screening. We used prospective longitudinal quantitative reverse transcription PCR testing to measure the viral RNA trajectories for 68 individuals during the resumption of the 2019-2020 National Basketball Association season. For 46 individuals with acute infections, we inferred the peak viral concentration and the duration of the viral proliferation and clearance phases. According to our mathematical model, we found that viral RNA concentrations peaked an average of 3.3 days (95% credible interval [CI] 2.5, 4.2) after first possible detectability at a cycle threshold value of 22.3 (95% CI 20.5, 23.9). The viral clearance phase lasted longer for symptomatic individuals (10.9 days [95% CI 7.9, 14.4]) than for asymptomatic individuals (7.8 days [95% CI 6.1, 9.7]). A second test within 2 days after an initial positive PCR test substantially improves certainty about a patient's infection stage. The effective sensitivity of a test intended to identify infectious individuals declines substantially with test turnaround time. These findings indicate that SARS-CoV-2 viral concentrations peak rapidly regardless of symptoms. Sequential tests can help reveal a patient's progress through infection stages. Frequent, rapid-turnaround testing is needed to effectively screen individuals before they become infectious.

Diverse functional autoantibodies in patients with COVID-19

COVID-19 manifests with a wide spectrum of clinical phenotypes that are characterized by exaggerated and misdirected host immune responses1-6. Although pathological innate immune activation is well-documented in severe disease1, the effect of autoantibodies on disease progression is less well-defined. Here we use a high-throughput autoantibody discovery technique known as rapid extracellular antigen profiling7 to screen a cohort of 194 individuals infected with SARS-CoV-2, comprising 172 patients with COVID-19 and 22 healthcare workers with mild disease or asymptomatic infection, for autoantibodies against 2,770 extracellular and secreted proteins (members of the exoproteome). We found that patients with COVID-19 exhibit marked increases in autoantibody reactivities as compared to uninfected individuals, and show a high prevalence of autoantibodies against immunomodulatory proteins (including cytokines, chemokines, complement components and cell-surface proteins). We established that these autoantibodies perturb immune function and impair virological control by inhibiting immunoreceptor signalling and by altering peripheral immune cell composition, and found that mouse surrogates of these autoantibodies increase disease severity in a mouse model of SARS-CoV-2 infection. Our analysis of autoantibodies against tissue-associated antigens revealed associations with specific clinical characteristics. Our findings suggest a pathological role for exoproteome-directed autoantibodies in COVID-19, with diverse effects on immune functionality and associations with clinical outcomes.

A stem-loop RNA RIG-I agonist confers prophylactic and therapeutic protection against acute and chronic SARS-CoV-2 infection in mice

As SARS-CoV-2 continues to cause morbidity and mortality around the world, there is an urgent need for the development of effective medical countermeasures. Here, we assessed the antiviral capacity of a minimal RIG-I agonist, stem-loop RNA 14 (SLR14), in viral control, disease prevention, post-infection therapy, and cross-variant protection in mouse models of SARS-CoV-2 infection. A single dose of SLR14 prevented viral replication in the lower respiratory tract and development of severe disease in a type I interferon (IFN-I) dependent manner. SLR14 demonstrated remarkable protective capacity against lethal SARS-CoV-2 infection when used prophylactically and retained considerable efficacy as a therapeutic agent. In immunodeficient mice carrying chronic SARS-CoV-2 infection, SLR14 elicited near-sterilizing innate immunity by inducing IFN-I responses in the absence of the adaptive immune system. In the context of infection with variants of concern (VOC), SLR14 conferred broad protection and uncovered an IFN-I resistance gradient across emerging VOC. These findings demonstrate the therapeutic potential of SLR14 as a host-directed, broad-spectrum antiviral for early post-exposure treatment and for treatment of chronically infected immunosuppressed patients.

Maternal respiratory SARS-CoV-2 infection in pregnancy is associated with a robust inflammatory response at the maternal-fetal interface

BACKGROUND

Pregnant women are at increased risk for severe outcomes from coronavirus disease 2019 (COVID-19), but the pathophysiology underlying this increased morbidity and its potential effect on the developing fetus is not well understood.

METHODS

We assessed placental histology, ACE2 expression, and viral and immune dynamics at the term placenta in pregnant women with and without respiratory severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection.

FINDINGS

The majority (13 of 15) of placentas analyzed had no detectable viral RNA. ACE2 was detected by immunohistochemistry in syncytiotrophoblast cells of the normal placenta during early pregnancy but was rarely seen in healthy placentas at full term, suggesting that low ACE2 expression may protect the term placenta from viral infection. Using immortalized cell lines and primary isolated placental cells, we found that cytotrophoblasts, the trophoblast stem cells and precursors to syncytiotrophoblasts, rather than syncytiotrophoblasts or Hofbauer cells, are most vulnerable to SARS-CoV-2 infection in vitro. To better understand potential immune mechanisms shielding placental cells from infection in vivo, we performed bulk and single-cell transcriptomics analyses and found that the maternal-fetal interface of SARS-CoV-2-infected women exhibited robust immune responses, including increased activation of natural killer (NK) and T cells, increased expression of interferon-related genes, as well as markers associated with pregnancy complications such as preeclampsia.

CONCLUSIONS

SARS-CoV-2 infection in late pregnancy is associated with immune activation at the maternal-fetal interface even in the absence of detectable local viral invasion.

FUNDING

NIH (T32GM007205, F30HD093350, K23MH118999, R01AI157488, U01DA040588) and Fast Grant funding support from Emergent Ventures at the Mercatus Center.

Diverse Functional Autoantibodies that Underlie Immune Perturbations in COVID-19

Infection with SARS-CoV-2 results in a broad spectrum of immunological and clinical outcomes for patients with COVID-19. While dysregulated innate immunity has been extensively explored in severe disease, autoantibody responses and their impact on disease trajectories remain uncharted. Here, we enrolled a cohort of 194 SARS-CoV-2 infected patients and uninfected individuals, and employed a high-throughput autoantibody discovery platform, Rapid Extracellular Antigen Profiling (REAP), to assess the prevalence of autoantibodies against 2,770 extracellular and secreted proteins in COVID-19 disease. We found that COVID-19 patients harbor dramatically expanded autoantibody repertoires compared to uninfected controls. We further established that immunomodulatory proteins including cytokines, chemokines, complements, and cell surface receptors, are frequent targets of autoantibodies. Longitudinal analysis revealed the existence of both pre-existing autoantibodies, as well as a broad subset of autoantibodies induced upon SARS-CoV-2 infection. We uncovered distinct mechanisms by which immune-targeting autoantibodies perturb immunological functions, including interference of immunoreceptor signaling, depletion of circulating leukocytes, and modulation of antiviral antibody responses. Concordantly, murine surrogates of these autoantibodies hinder immune activation and exacerbate disease severity in a mouse model of SARS-CoV-2 infection. In summary, through the lens of an unbiased proteome-scale survey for autoantibody targets, these findings implicate humoral immunopathology as an integral aspect of COVID-19 pathogenesis with diverse impacts on immune functionality and clinical outcomes.

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.

Early introductions and transmission of SARS-CoV-2 variant B.1.1.7 in the United States

The emergence and spread of SARS-CoV-2 lineage B.1.1.7, first detected in the United Kingdom, has become a global public health concern because of its increased transmissibility. Over 2,500 COVID-19 cases associated with this variant have been detected in the United States (US) since December 2020, but the extent of establishment is relatively unknown. Using travel, genomic, and diagnostic data, we highlight that the primary ports of entry for B.1.1.7 in the US were in New York, California, and Florida. Furthermore, we found evidence for many independent B.1.1.7 establishments starting in early December 2020, followed by interstate spread by the end of the month. Finally, we project that B.1.1.7 will be the dominant lineage in many states by mid- to late March. Thus, genomic surveillance for B.1.1.7 and other variants urgently needs to be enhanced to better inform the public health response.

Stability of SARS-CoV-2 RNA in Nonsupplemented Saliva

The expense of saliva collection devices designed to stabilize severe acute respiratory syndrome coronavirus 2 RNA is prohibitive to mass testing. However, virus RNA in nonsupplemented saliva is stable for extended periods and at elevated temperatures. Simple plastic tubes for saliva collection will make large-scale testing and continued surveillance easier.

Case Study: Longitudinal immune profiling of a SARS-CoV-2 reinfection in a solid organ transplant recipient

Prior to the emergence of antigenically distinct SARS-CoV-2 variants, reinfections were reported infrequently - presumably due to the generation of durable and protective immune responses. However, case reports also suggested that rare, repeated infections may occur as soon as 48 days following initial disease onset. The underlying immunologic deficiencies enabling SARS-CoV-2 reinfections are currently unknown. Here we describe a renal transplant recipient who developed recurrent, symptomatic SARS-CoV-2 infection - confirmed by whole virus genome sequencing - 7 months after primary infection. To elucidate the immunological mechanisms responsible for SARS-CoV-2 reinfection, we performed longitudinal profiling of cellular and humoral responses during both primary and recurrent SARS-CoV-2 infection. We found that the patient responded to the primary infection with transient, poor-quality adaptive immune responses. The patient's immune system was further compromised by intervening treatment for acute rejection of the renal allograft prior to reinfection. Importantly, we also identified the development of neutralizing antibodies and the formation of humoral memory responses prior to SARS-CoV-2 reinfection. However, these neutralizing antibodies failed to confer protection against reinfection, suggesting that additional factors are required for efficient prevention of SARS-CoV-2 reinfection. Further, we found no evidence supporting viral evasion of primary adaptive immune responses, suggesting that susceptibility to reinfection may be determined by host factors rather than pathogen adaptation in this patient. In summary, our study suggests that a low neutralizing antibody presence alone is not sufficient to confer resistance against reinfection. Thus, patients with solid organ transplantation, or patients who are otherwise immunosuppressed, who recover from infection with SARS-CoV-2 may not develop sufficient protective immunity and are at risk of reinfection.

MOG-associated encephalitis following SARS-COV-2 infection

A variety of neurologic manifestations of COVID-19 infections have been reported. Here, we present a case of steroid-responsive MOG-antibody associated encephalitis, characterized by cognitive decline, headaches, fever, unilateral FLAIR-hyperintensities, and leptomeningeal enhancement, that occurred in the setting of recent COVID-19 infection.

Diverse Functional Autoantibodies in Patients with COVID-19

COVID-19 manifests with a wide spectrum of clinical phenotypes that are characterized by exaggerated and misdirected host immune responses1-8. While pathological innate immune activation is well documented in severe disease1, the impact of autoantibodies on disease progression is less defined. Here, we used a high-throughput autoantibody discovery technique called Rapid Extracellular Antigen Profiling (REAP) to screen a cohort of 194 SARS-CoV-2 infected COVID-19 patients and healthcare workers for autoantibodies against 2,770 extracellular and secreted proteins (the "exoproteome"). We found that COVID-19 patients exhibit dramatic increases in autoantibody reactivities compared to uninfected controls, with a high prevalence of autoantibodies against immunomodulatory proteins including cytokines, chemokines, complement components, and cell surface proteins. We established that these autoantibodies perturb immune function and impair virological control by inhibiting immunoreceptor signaling and by altering peripheral immune cell composition, and found that murine surrogates of these autoantibodies exacerbate disease severity in a mouse model of SARS-CoV-2 infection. Analysis of autoantibodies against tissue-associated antigens revealed associations with specific clinical characteristics and disease severity. In summary, these findings implicate a pathological role for exoproteome-directed autoantibodies in COVID-19 with diverse impacts on immune functionality and associations with clinical outcomes.

Reproducing the Rift Valley fever virus mosquito-lamb-mosquito transmission cycle

Rift Valley fever virus (RVFV) is a mosquito-borne bunyavirus that is pathogenic to ruminants and humans. The virus is endemic to Africa and the Arabian Peninsula where outbreaks are characterized by abortion storms and mortality of newborns, particularly in sheep herds. Vector competence experiments in laboratory settings have suggested that over 50 mosquito species are capable of transmitting RVFV. Transmission of mosquito-borne viruses in the field is however influenced by numerous factors, including population densities, blood feeding behavior, extrinsic incubation period, longevity of vectors, and viremia levels in vertebrate hosts. Animal models to study these important aspects of RVFV transmission are currently lacking. In the present work, RVFV was transmitted to European (Texel-swifter cross-breed) lambs by laboratory-reared Aedes aegypti mosquitoes that were infected either by membrane feeding on a virus-spiked blood meal or by feeding on lambs that developed viremia after intravenous inoculation of RVFV. Feeding of mosquitoes on viremic lambs resulted in strikingly higher infection rates as compared to membrane feeding. Subsequent transmission of RVFV from lamb to lamb by infected mosquitoes was highly efficient in both models. The animal models described here can be used to study mosquito-mediated transmission of RVFV among the major natural target species and to evaluate the efficacy of vaccines against mosquito-mediated RVFV infection.

Saliva viral load is a dynamic unifying correlate of COVID-19 severity and mortality

While several clinical and immunological parameters correlate with disease severity and mortality in SARS-CoV-2 infection, work remains in identifying unifying correlates of coronavirus disease 2019 (COVID-19) that can be used to guide clinical practice. Here, we examine saliva and nasopharyngeal (NP) viral load over time and correlate them with patient demographics, and cellular and immune profiling. We found that saliva viral load was significantly higher in those with COVID-19 risk factors; that it correlated with increasing levels of disease severity and showed a superior ability over nasopharyngeal viral load as a predictor of mortality over time (AUC=0.90). A comprehensive analysis of immune factors and cell subsets revealed strong predictors of high and low saliva viral load, which were associated with increased disease severity or better overall outcomes, respectively. Saliva viral load was positively associated with many known COVID-19 inflammatory markers such as IL-6, IL-18, IL-10, and CXCL10, as well as type 1 immune response cytokines. Higher saliva viral loads strongly correlated with the progressive depletion of platelets, lymphocytes, and effector T cell subsets including circulating follicular CD4 T cells (cTfh). Anti-spike (S) and anti-receptor binding domain (RBD) IgG levels were negatively correlated with saliva viral load showing a strong temporal association that could help distinguish severity and mortality in COVID-19. Finally, patients with fatal COVID-19 exhibited higher viral loads, which correlated with the depletion of cTfh cells, and lower production of anti-RBD and anti-S IgG levels. Together these results demonstrated that viral load - as measured by saliva but not nasopharyngeal - is a dynamic unifying correlate of disease presentation, severity, and mortality over time.

SalivaDirect: A simplified and flexible platform to enhance SARS-CoV-2 testing capacity

Background

Scaling SARS-CoV-2 testing to meet demands of safe reopenings continues to be plagued by assay costs and supply chain shortages. In response, we developed SalivaDirect, which received Emergency Use Authorization (EUA) from the U.S. Food and Drug Administration (FDA).

Methods

We simplified our saliva-based diagnostic test by (1) not requiring collection tubes with preservatives, (2) replacing nucleic acid extraction with a simple enzymatic and heating step, and (3) testing specimens with a dualplex qRT-PCR assay. Moreover, we validated SalivaDirect with reagents and instruments from multiple vendors to minimize supply chain issues.

Findings

From our hospital cohort, we show a high positive agreement (94%) between saliva tested with SalivaDirect and nasopharyngeal swabs tested with a commercial qRT-PCR kit. In partnership with the National Basketball Association (NBA) and National Basketball Players Association (NBPA), we tested 3,779 saliva specimens from healthy individuals and detected low rates of invalid (0.3%) and false-positive (<0.05%) results.

Conclusions

We demonstrate that saliva is a valid alternative to swabs for SARS-CoV-2 screening and that SalivaDirect can make large-scale testing more accessible and affordable. Uniquely, we can designate other laboratories to use our sensitive, flexible, and simplified platform under our EUA (https://publichealth.yale.edu/salivadirect/).

Funding

This study was funded by the NBA and NBPA (N.D.G.), the Huffman Family Donor Advised Fund (N.D.G.), a Fast Grant from Emergent Ventures at the Mercatus Center at George Mason University (N.D.G.), the Yale Institute for Global Health (N.D.G.), and the Beatrice Kleinberg Neuwirth Fund (A.I.K.). C.B.F.V. is supported by NWO Rubicon 019.181EN.004.

Frequent testing is critical to limit SARS-CoV-2 transmission. In response to this need, we developed SalivaDirect, a sensitive, simplified, and flexible testing framework, which received Emergency Use Authorization (EUA) from the U.S. Food and Drug Administration (FDA). We tested saliva collected from a hospital cohort and showed a high positive agreement (94%) as compared to paired nasopharyngeal swabs tested with a commercial diagnostic kit. Then, we partnered with the National Basketball Association (NBA) to test a large cohort of mostly healthy individuals, and we detected low rates of invalid (0.3%) and false-positive (0.03%–0.05%) results. Our study shows that SalivaDirect can help to increase testing capacity by providing access to an affordable framework that is less prone to supply chain shortages.

Kinetics of antibody responses dictate COVID-19 outcome

Recent studies have provided insights into innate and adaptive immune dynamics in coronavirus disease 2019 (COVID-19). Yet, the exact feature of antibody responses that governs COVID-19 disease outcomes remain unclear. Here, we analysed humoral immune responses in 209 asymptomatic, mild, moderate and severe COVID-19 patients over time to probe the nature of antibody responses in disease severity and mortality. We observed a correlation between anti-Spike (S) IgG levels, length of hospitalization and clinical parameters associated with worse clinical progression. While high anti-S IgG levels correlated with worse disease severity, such correlation was time-dependent. Deceased patients did not have higher overall humoral response than live discharged patients. However, they mounted a robust, yet delayed response, measured by anti-S, anti-RBD IgG, and neutralizing antibody (NAb) levels, compared to survivors. Delayed seroconversion kinetics correlated with impaired viral control in deceased patients. Finally, while sera from 89% of patients displayed some neutralization capacity during their disease course, NAb generation prior to 14 days of disease onset emerged as a key factor for recovery. These data indicate that COVID-19 mortality does not correlate with the cross-sectional antiviral antibody levels per se, but rather with the delayed kinetics of NAb production.

Kinetics of antibody responses dictate COVID-19 outcome

Recent studies have provided insights into innate and adaptive immune dynamics in coronavirus disease 2019 (COVID-19). Yet, the exact feature of antibody responses that governs COVID-19 disease outcomes remain unclear. Here, we analysed humoral immune responses in 209 asymptomatic, mild, moderate and severe COVID-19 patients over time to probe the nature of antibody responses in disease severity and mortality. We observed a correlation between anti-Spike (S) IgG levels, length of hospitalization and clinical parameters associated with worse clinical progression. While high anti-S IgG levels correlated with worse disease severity, such correlation was time-dependent. Deceased patients did not have higher overall humoral response than live discharged patients. However, they mounted a robust, yet delayed response, measured by anti-S, anti-RBD IgG, and neutralizing antibody (NAb) levels, compared to survivors. Delayed seroconversion kinetics correlated with impaired viral control in deceased patients. Finally, while sera from 89% of patients displayed some neutralization capacity during their disease course, NAb generation prior to 14 days of disease onset emerged as a key factor for recovery. These data indicate that COVID-19 mortality does not correlate with the cross-sectional antiviral antibody levels per se , but rather with the delayed kinetics of NAb production.

Impact of Gut Bacteria on the Infection and Transmission of Pathogenic Arboviruses by Biting Midges and Mosquitoes

Tripartite interactions among insect vectors, midgut bacteria, and viruses may determine the ability of insects to transmit pathogenic arboviruses. Here, we investigated the impact of gut bacteria on the susceptibility of Culicoides nubeculosus and Culicoides sonorensis biting midges for Schmallenberg virus, and of Aedes aegypti mosquitoes for Zika and chikungunya viruses. Gut bacteria were manipulated by treating the adult insects with antibiotics. The gut bacterial communities were investigated using Illumina MiSeq sequencing of 16S rRNA, and susceptibility to arbovirus infection was tested by feeding insects with an infectious blood meal. Antibiotic treatment led to changes in gut bacteria for all insects. Interestingly, the gut bacterial composition of untreated Ae. aegypti and C. nubeculosus showed Asaia as the dominant genus, which was drastically reduced after antibiotic treatment. Furthermore, antibiotic treatment resulted in relatively more Delftia bacteria in both biting midge species, but not in mosquitoes. Antibiotic treatment and subsequent changes in gut bacterial communities were associated with a significant, 1.8-fold increased infection rate of C. nubeculosus with Schmallenberg virus, but not for C. sonorensis. We did not find any changes in infection rates for Ae. aegypti mosquitoes with Zika or chikungunya virus. We conclude that resident gut bacteria may dampen arbovirus transmission in biting midges, but not so in mosquitoes. Use of antimicrobial compounds at livestock farms might therefore have an unexpected contradictory effect on the health of animals, by increasing the transmission of viral pathogens by biting midges.