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.

Viral kinetics of sequential SARS-CoV-2 infections.

The impact of a prior SARS-CoV-2 infection on the progression of subsequent infections has been unclear. Using a convenience sample of 94,812 longitudinal RT-qPCR measurements from anterior nares and oropharyngeal swabs, we identified 71 individuals with two well-sampled SARS-CoV-2 infections between March 11th, 2020, and July 28th, 2022. We compared the SARS-CoV-2 viral kinetics of first vs. second infections in this group, adjusting for viral variant, vaccination status, and age. Relative to first infections, second infections usually featured a faster clearance time. Furthermore, a person's relative (rank-order) viral clearance time, compared to others infected with the same variant, was roughly conserved across first and second infections, so that individuals who had a relatively fast clearance time in their first infection also tended to have a relatively fast clearance time in their second infection (Spearman correlation coefficient: 0.30, 95% credible interval (0.12, 0.46)). These findings provide evidence that, like vaccination, immunity from a prior SARS-CoV-2 infection shortens the duration of subsequent acute SARS-CoV-2 infections principally by reducing viral clearance time. Additionally, there appears to be an inherent element of the immune response, or some other host factor, that shapes a person's relative ability to clear SARS-CoV-2 infection that persists across sequential infections.

Phylogeographic reconstruction of the emergence and spread of Powassan virus in the northeastern United States.

Powassan virus is an emerging tick-borne virus of concern for public health, but very little is known about its transmission patterns and ecology. Here, we expanded the genomic dataset by sequencing 279 Powassan viruses isolated from Ixodes scapularis ticks from the northeastern United States. Our phylogeographic reconstructions revealed that Powassan virus lineage II was likely introduced or emerged from a relict population in the Northeast between 1940 and 1975. Sequences strongly clustered by sampling location, suggesting a highly focal geographical distribution. Our analyses further indicated that Powassan virus lineage II emerged in the northeastern United States mostly following a south-to-north pattern, with a weighted lineage dispersal velocity of ~3 km/y. Since the emergence in the Northeast, we found an overall increase in the effective population size of Powassan virus lineage II, but with growth stagnating during recent years. The cascading effect of population expansion of white-tailed deer and I. scapularis populations likely facilitated the emergence of Powassan virus in the northeastern United States.

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.

How the Commonwealth of the Northern Mariana Islands stalled COVID-19 for 22 months and managed its first significant community transmission.

Objective: The Commonwealth of the Northern Mariana Islands (CNMI) is a remote Pacific island territory with a population of 47 329 that successfully prevented the significant introduction of coronavirus disease (COVID-19) until late 2021. This study documents how the response to the introduction of COVID-19 in CNMI in 2021 was conducted with limited resources without overwhelming local clinical capacity or compromising health service delivery for the population. Methods: Data from COVID-19 case investigations, contact tracing, the Commonwealth's immunization registry and whole genome sequencing were collated and analysed as part of this study. Results: Between 26 March 2020 and 31 December 2021, 3281 cases and 14 deaths due to COVID-19 were reported in CNMI (case fatality rate, 0.4%). While notification rates were highest among younger age groups, hospitalization and mortality rates were disproportionately greater among those aged > 50 years and among the unvaccinated. The first widespread community transmission in CNMI was detected in October 2021, with genomic epidemiology and contact tracing data indicating a single introduction event involving the AY.25 lineage and subsequent rapid community spread. Vaccination coverage was high before widespread transmission occurred in October 2021 and increased further over the study period. Discussion: Robust preparedness and strong leadership generated resilience within the public health sector such that COVID-19 did not overwhelm CNMI's health system as it did in other jurisdictions and countries around the world. At no point was hospital capacity exceeded, and all patients received adequate care without the need for health-care rationing.

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.

Quantifying the impact of immune history and variant on SARS-CoV-2 viral kinetics and infection rebound: A retrospective cohort study.

Background: The combined impact of immunity and SARS-CoV-2 variants on viral kinetics during infections has been unclear. Methods: We characterized 1,280 infections from the National Basketball Association occupational health cohort identified between June 2020 and January 2022 using serial RT-qPCR testing. Logistic regression and semi-mechanistic viral RNA kinetics models were used to quantify the effect of age, variant, symptom status, infection history, vaccination status and antibody titer to the founder SARS-CoV-2 strain on the duration of potential infectiousness and overall viral kinetics. The frequency of viral rebounds was quantified under multiple cycle threshold (Ct) value-based definitions. Results: Among individuals detected partway through their infection, 51.0% (95% credible interval [CrI]: 48.3-53.6%) remained potentially infectious (Ct <30) 5 days post detection, with small differences across variants and vaccination status. Only seven viral rebounds (0.7%; N=999) were observed, with rebound defined as 3+days with Ct <30 following an initial clearance of 3+days with Ct ≥30. High antibody titers against the founder SARS-CoV-2 strain predicted lower peak viral loads and shorter durations of infection. Among Omicron BA.1 infections, boosted individuals had lower pre-booster antibody titers and longer clearance times than non-boosted individuals. Conclusions: SARS-CoV-2 viral kinetics are partly determined by immunity and variant but dominated by individual-level variation. Since booster vaccination protects against infection, longer clearance times for BA.1-infected, boosted individuals may reflect a less effective immune response, more common in older individuals, that increases infection risk and reduces viral RNA clearance rate. The shifting landscape of viral kinetics underscores the need for continued monitoring to optimize isolation policies and to contextualize the health impacts of therapeutics and vaccines. Funding: Supported in part by CDC contract #200-2016-91779, a sponsored research agreement to Yale University from the National Basketball Association contract #21-003529, and the National Basketball Players Association.

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.

Combining genomic and epidemiological data to compare the transmissibility of SARS-CoV-2 variants Alpha and Iota.

SARS-CoV-2 variants shaped the second year of the COVID-19 pandemic and the discourse around effective control measures. Evaluating the threat posed by a new variant is essential for adapting response efforts when community transmission is detected. In this study, we compare the dynamics of two variants, Alpha and Iota, by integrating genomic surveillance data to estimate the effective reproduction number (Rt) of the variants. We use Connecticut, United States, in which Alpha and Iota co-circulated in 2021. We find that the Rt of these variants were up to 50% larger than that of other variants. We then use phylogeography to show that while both variants were introduced into Connecticut at comparable frequencies, clades that resulted from introductions of Alpha were larger than those resulting from Iota introductions. By monitoring the dynamics of individual variants throughout our study period, we demonstrate the importance of routine surveillance in the response to COVID-19.

Sequencing SARS-CoV-2 genomes from saliva.

Genomic sequencing is crucial to understanding the epidemiology and evolution of Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Often, genomic studies rely on remnant diagnostic material, typically nasopharyngeal (NP) 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 NP 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.

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.

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.

Severe Acute Respiratory Syndrome Coronavirus 2 Reinfection: A Case Series From a 12-Month Longitudinal Occupational Cohort.

Findings are described in 7 patients with severe acute respiratory syndrome coronavirus 2 reinfection from the National Basketball Association 2020-2021 occupational testing cohort, including clinical details, antibody test results, genomic sequencing, and longitudinal reverse-transcription polymerase chain reaction results. Reinfections were infrequent and varied in clinical presentation, viral dynamics, and immune response.

Longitudinal Immune Profiling of a Severe Acute Respiratory Syndrome Coronavirus 2 Reinfection in a Solid Organ Transplant Recipient.

BACKGROUND: The underlying immunologic deficiencies enabling severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) reinfection are currently unknown. We describe deep longitudinal immune profiling of a transplant recipient hospitalized twice for coronavirus disease 2019 (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 investigations were performed on samples from the primary and secondary infections. RESULTS: Whole viral genome sequencing and phylogenetic 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. CONCLUSIONS: 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 naive CD4 T-cell pools. Seropositivity alone may not be a perfect correlate of protection in immunocompromised patients.

Impact of extrinsic incubation temperature on natural selection during Zika virus infection of Aedes aegypti and Aedes albopictus.

Arthropod-borne viruses (arboviruses) require replication across a wide range of temperatures to perpetuate. While vertebrate hosts tend to maintain temperatures of approximately 37°C-40°C, arthropods are subject to ambient temperatures which can have a daily fluctuation of > 10°C. Temperatures impact vector competence, extrinsic incubation period, and mosquito survival unimodally, with optimal conditions occurring at some intermediate temperature. In addition, the mean and range of daily temperature fluctuations influence arbovirus perpetuation and vector competence. The impact of temperature on arbovirus genetic diversity during systemic mosquito infection, however, is poorly understood. Therefore, we determined how constant extrinsic incubation temperatures of 25°C, 28°C, 32°C, and 35°C control Zika virus (ZIKV) vector competence and population dynamics within Aedes aegypti and Aedes albopictus mosquitoes. We also examined fluctuating temperatures which better mimic field conditions in the tropics. We found that vector competence varied in a unimodal manner for constant temperatures peaking between 28°C and 32°C for both Aedes species. Transmission peaked at 10 days post-infection for Aedes aegypti and 14 days for Aedes albopictus. Conversely, fluctuating temperature decreased vector competence. Using RNA-seq to characterize ZIKV population structure, we identified that temperature alters the selective environment in unexpected ways. During mosquito infection, constant temperatures more often elicited positive selection whereas fluctuating temperatures led to strong purifying selection in both Aedes species. These findings demonstrate that temperature has multiple impacts on ZIKV biology, including major effects on the selective environment within mosquitoes.

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

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Delta variant quickly rose to dominance in mid-2021, displacing other variants, including Alpha. Studies using data from the United Kingdom and India estimated that Delta was 40-80% more transmissible than Alpha, allowing Delta to become the globally dominant variant. However, it was unclear if the ostensible difference in relative transmissibility was due mostly to innate properties of Delta's infectiousness or differences in the study populations. To investigate, we formed a partnership with SARS-CoV-2 genomic surveillance programs from all six New England US states. By comparing logistic growth rates, we found that Delta emerged 37-163% faster than Alpha in early 2021 (37% Massachusetts, 75% New Hampshire, 95% Maine, 98% Rhode Island, 151% Connecticut, and 163% Vermont). We next computed variant-specific effective reproductive numbers and estimated that Delta was 58-120% more transmissible than Alpha across New England (58% New Hampshire, 68% Massachusetts, 76% Connecticut, 85% Rhode Island, 98% Maine, and 120% Vermont). Finally, using RT-PCR data, we estimated that Delta infections generate on average ∼6 times more viral RNA copies per mL than Alpha infections. Overall, our evidence indicates that Delta's enhanced transmissibility could be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on the underlying immunity and behavior of distinct populations.

Tracing the Origin, Spread, and Molecular Evolution of Zika Virus in Puerto Rico, 2016-2017.

We reconstructed the 2016-2017 Zika virus epidemic in Puerto Rico by using complete genomes to uncover the epidemic's origin, spread, and evolutionary dynamics. Our study revealed that the epidemic was propelled by multiple introductions that spread across the island, intricate evolutionary patterns, and ≈10 months of cryptic transmission.

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.