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

High-resolution and real-time wastewater viral surveillance by Nanopore sequencing

Wastewater genomic sequencing stands as a pivotal complementary tool for viral surveillance in populations. While long-read Nanopore sequencing is a promising platform to provide real-time genomic data, concerns over the sequencing accuracy of the earlier Nanopore versions have somewhat restrained its widespread application in wastewater analysis. Here, we evaluate the latest improved version of Nanopore sequencing (R10.4.1), using SARS-CoV-2 as the model infectious virus, to demonstrate its effectiveness in wastewater viral monitoring. By comparing amplicon lengths of 400 bp and 1200 bp, we revealed that shorter PCR amplification is more suitable for wastewater samples due to viral genome fragmentation. Utilizing mock wastewater samples, we validated the reliability of Nanopore sequencing for variant identification by comparing it with Illumina sequencing results. The strength of Nanopore sequencing in generating real-time genomic data for providing early warning signals was also showcased, indicating that as little as 0.001 Gb of data can provide accurate results for variant prevalence. Our evaluation also identified optimal alteration frequency cutoffs (>50 %) for precise mutation profiling, achieving >99 % precision in detecting single nucleotide variants (SNVs) and insertions/deletions (indels). Monitoring two major wastewater treatment plants in Hong Kong from September 2022 to April 2023, covering over 4.5 million population, we observed a transition in dominant variants from BA.5 to XBB lineages, with XBB.1.5 being the most prevalent variants. Mutation detection also highlighted the potential of wastewater Nanopore sequencing in uncovering novel mutations and revealed links between signature mutations and specific variants. This study not only reveals the environmental implications of Nanopore sequencing in SARS-CoV-2 surveillance but also underscores its potential in broader applications including environmental health monitoring of other epidemic viruses, which could significantly enhance the field of wastewater-based epidemiology.

Variations in the persistence of 5'-end genomic and subgenomic SARS-CoV-2 RNAs in wastewater from aircraft, airports and wastewater treatment plants

Wastewater sequencing has become a powerful supplement to clinical testing in monitoring SARS-CoV-2 infections in the post-COVID-19 pandemic era. While its applications in measuring the viral burden and main circulating lineages in the community have proved their efficacy, the variations in sequencing quality and coverage across the different regions of the SARS-CoV-2 genome are not well understood. Furthermore, it is unclear how different sample origins, viral extraction and concentration methods and environmental factors impact the reads sequenced from wastewater. Using high-coverage, amplicon-based, paired-end read sequencing of viral RNA extracted from wastewater collected directly from aircraft, pooled from different aircraft and airport buildings or from regular wastewater plants, we assessed the genome coverage across the sample groups with a focus on the 5'-end region covering the leader sequence and investigated whether it was possible to detect subgenomic RNA from viral material recovered from wastewater. We identified distinct patterns in the persistence of the different genomic regions across the different types of wastewaters and the existence of chimeric reads mapping to non-amplified regions. Our findings suggest that preservation of the 5'-end of the genome and the ability to detect subgenomic RNA reads, though highly susceptible to environment and sample processing conditions, may be indicative of the quality and amount of the viral RNA present in wastewater.

Considerations and Opportunities for Probe Capture Enrichment Sequencing of Emerging Viruses from Wastewater

Until recently, wastewater-based monitoring for pathogens of public health concern primarily used PCR-based quantification methods and targeted sequencing for specific pathogens (e.g., SARS-CoV-2). In the past three years, researchers have expanded sequencing to monitor a broad range of pathogens, applying probe capture enrichment to wastewater. The goals of those studies included (1) monitoring and expanding fundamental knowledge of disease dynamics for known pathogens and (2) evaluating the potential for early detection of emerging diseases resulting from zoonotic spillover or novel viral variants. Several studies using off-the-shelf probe panels designed for clinical and environmental surveillance reported that enrichment increased virus relative abundance but did not recover complete genomes for most nonenteric viruses. Based on our experience and recent results reported by others using these panels for wastewater, clinical, and synthetic samples, we discuss challenges and technical factors that affect the rates of false positive and false negative results. We identify trade-offs and opportunities throughout the workflow, including in wastewater sample processing, probe panel design, and bioinformatic analysis. We suggest tailored methods of virus concentration and background removal, carefully designed probe panels, and multithresholded bioinformatics analysis.

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.

Optimisation and evaluation of viral genomic sequencing of SARS-CoV-2 rapid diagnostic tests: a laboratory and cohort-based study

BACKGROUND

Sequencing of SARS-CoV-2 from rapid diagnostic tests (RDTs) can bolster viral genomic surveillance efforts; however, approaches to maximise and standardise pathogen genome recovery from RDTs remain underdeveloped. We aimed to systematically optimise the elution of genetic material from RDT components and to evaluate the efficacy of RDT sequencing for outbreak investigation.

METHODS

In this laboratory and cohort-based study we seeded RDTs with inactivated SARS-CoV-2 to optimise the elution of genomic material from RDT lateral flow strips. We measured the effect of changes in buffer type, time in buffer, and rotation on PCR cycle threshold (Ct) value. We recruited individuals older than 18 years residing in the greater Boston area, MA, USA, from July 18 to Nov 5, 2022, via email advertising to students and staff at Harvard University, MA, USA, and via broad social media advertising. All individuals recruited were within 5 days of a positive diagnostic test for SARS-CoV-2; no other relevant exclusion criteria were applied. Each individual completed two RDTs and one PCR swab. On Dec 29, 2022, we also collected RDTs from a convenience sample of individuals who were positive for SARS-CoV-2 and associated with an outbreak at a senior housing facility in MA, USA. We extracted all returned PCR swabs and RDT components (ie, swab, strip, or buffer); samples with a Ct of less than 40 were subject to amplicon sequencing. We compared the efficacy of elution and sequencing across RDT brands and components and used RDT-derived sequences to infer transmission links within the outbreak at the senior housing facility. We conducted metagenomic sequencing of negative RDTs from symptomatic individuals living in the senior housing facility.

FINDINGS

Neither elution duration of greater than 10 min nor rotation during elution impacted viral titres. Elution in Buffer AVL (Ct=31·4) and Tris-EDTA Buffer (Ct=30·8) were equivalent (p=0·34); AVL outperformed elution in lysis buffer and 50% lysis buffer (Ct=40·0, p=0·0029 for both) as well as Universal Viral Transport Medium (Ct=36·7, p=0·079). Performance of RDT strips was poorer than that of matched PCR swabs (mean Ct difference 10·2 [SD 4·3], p<0·0001); however, RDT swabs performed similarly to PCR swabs (mean Ct difference 4·1 [5·2], p=0·055). No RDT brand significantly outperformed another. Across sample types, viral load predicted the viral genome assembly length. We assembled greater than 80% complete genomes from 12 of 17 RDT-derived swabs, three of 18 strips, and four of 11 residual buffers. We generated outbreak-associated SARS-CoV-2 genomes using both amplicon and metagenomic sequencing and identified multiple introductions of the virus that resulted in downstream transmission.

INTERPRETATION

RDT-derived swabs are a reasonable alternative to PCR swabs for viral genomic surveillance and outbreak investigation. RDT-derived lateral flow strips yield accurate, but significantly fewer, viral reads than matched PCR swabs. Metagenomic sequencing of negative RDTs can identify viruses that might underlie patient symptoms.

FUNDING

The National Science Foundation, the Hertz Foundation, the National Institute of General Medical Sciences, Harvard Medical School, the Howard Hughes Medical Institute, the US Centers for Disease Control and Prevention, the Broad Institute and the National Institute of Allergy and Infectious Diseases.

Dynamic diversity of SARS-CoV-2 genetic mutations in a lung transplantation patient with persistent COVID-19

Numerous SARS-CoV-2 variant strains with altered characteristics have emerged since the onset of the COVID-19 pandemic. Remdesivir (RDV), a ribonucleotide analogue inhibitor of viral RNA polymerase, has become a valuable therapeutic agent. However, immunosuppressed hosts may respond inadequately to RDV and develop chronic persistent infections. A patient with respiratory failure caused by interstitial pneumonia, who had undergone transplantation of the left lung, developed COVID-19 caused by Omicron BA.5 strain with persistent chronic viral shedding, showing viral fusogenicity. Genome-wide sequencing analyses revealed the occurrence of several viral mutations after RDV treatment, followed by dynamic changes in the viral populations. The C799F mutation in nsp12 was found to play a pivotal role in conferring RDV resistance, preventing RDV-triphosphate from entering the active site of RNA-dependent RNA polymerase. The occurrence of diverse mutations is a characteristic of SARS-CoV-2, which mutates frequently. Herein, we describe the clinical case of an immunosuppressed host in whom inadequate treatment resulted in highly diverse SARS-CoV-2 mutations that threatened the patient's health due to the development of drug-resistant variants.

Performance of amplicon and capture based next-generation sequencing approaches for the epidemiological surveillance of Omicron SARS-CoV-2 and other variants of concern

To control the SARS-CoV-2 pandemic, healthcare systems have focused on ramping up their capacity for epidemiological surveillance through viral whole genome sequencing. In this paper, we tested the performance of two protocols of SARS-CoV-2 nucleic acid enrichment, an amplicon enrichment using different versions of the ARTIC primer panel and a hybrid-capture method using KAPA RNA Hypercap. We focused on the challenge of the Omicron variant sequencing, the advantages of automated library preparation and the influence of the bioinformatic analysis in the final consensus sequence. All 94 samples were sequenced using Illumina iSeq 100 and analysed with two bioinformatic pipelines: a custom-made pipeline and an Illumina-owned pipeline. We were unsuccessful in sequencing six samples using the capture enrichment due to low reads. On the other hand, amplicon dropout and mispriming caused the loss of mutation G21987A and the erroneous addition of mutation T15521A respectively using amplicon enrichment. Overall, we found high sequence agreement regardless of method of enrichment, bioinformatic pipeline or the use of automation for library preparation in eight different SARS-CoV-2 variants. Automation and the use of a simple app for bioinformatic analysis can simplify the genotyping process, making it available for more diagnostic facilities and increasing global vigilance.

A hepatitis B virus (HBV) sequence variation graph improves alignment and sample-specific consensus sequence construction

Nearly 300 million individuals live with chronic hepatitis B virus (HBV) infection (CHB), for which no curative therapy is available. As viral diversity is associated with pathogenesis and immunological control of infection, improved methods to characterize this diversity could aid drug development efforts. Conventionally, viral sequencing data are mapped/aligned to a reference genome, and only the aligned sequences are retained for analysis. Thus, reference selection is critical, yet selecting the most representative reference a priori remains difficult. We investigate an alternative pangenome approach which can combine multiple reference sequences into a graph which can be used during alignment. Using simulated short-read sequencing data generated from publicly available HBV genomes and real sequencing data from an individual living with CHB, we demonstrate alignment to a phylogenetically representative 'genome graph' can improve alignment, avoid issues of reference ambiguity, and facilitate the construction of sample-specific consensus sequences more genetically similar to the individual's infection. Graph-based methods can, therefore, improve efforts to characterize the genetics of viral pathogens, including HBV, and have broader implications in host-pathogen research.

INSaFLU-TELEVIR: an open web-based bioinformatics suite for viral metagenomic detection and routine genomic surveillance

BACKGROUND

Implementation of clinical metagenomics and pathogen genomic surveillance can be particularly challenging due to the lack of bioinformatics tools and/or expertise. In order to face this challenge, we have previously developed INSaFLU, a free web-based bioinformatics platform for virus next-generation sequencing data analysis. Here, we considerably expanded its genomic surveillance component and developed a new module (TELEVIR) for metagenomic virus identification.

RESULTS

The routine genomic surveillance component was strengthened with new workflows and functionalities, including (i) a reference-based genome assembly pipeline for Oxford Nanopore technologies (ONT) data; (ii) automated SARS-CoV-2 lineage classification; (iii) Nextclade analysis; (iv) Nextstrain phylogeographic and temporal analysis (SARS-CoV-2, human and avian influenza, monkeypox, respiratory syncytial virus (RSV A/B), as well as a "generic" build for other viruses); and (v) algn2pheno for screening mutations of interest. Both INSaFLU pipelines for reference-based consensus generation (Illumina and ONT) were benchmarked against commonly used command line bioinformatics workflows for SARS-CoV-2, and an INSaFLU snakemake version was released. In parallel, a new module (TELEVIR) for virus detection was developed, after extensive benchmarking of state-of-the-art metagenomics software and following up-to-date recommendations and practices in the field. TELEVIR allows running complex workflows, covering several combinations of steps (e.g., with/without viral enrichment or host depletion), classification software (e.g., Kaiju, Kraken2, Centrifuge, FastViromeExplorer), and databases (RefSeq viral genome, Virosaurus, etc.), while culminating in user- and diagnosis-oriented reports. Finally, to potentiate real-time virus detection during ONT runs, we developed findONTime, a tool aimed at reducing costs and the time between sample reception and diagnosis.

CONCLUSIONS

The accessibility, versatility, and functionality of INSaFLU-TELEVIR are expected to supply public and animal health laboratories and researchers with a user-oriented and pan-viral bioinformatics framework that promotes a strengthened and timely viral metagenomic detection and routine genomics surveillance. INSaFLU-TELEVIR is compatible with Illumina, Ion Torrent, and ONT data and is freely available at https://insaflu.insa.pt/ (online tool) and https://github.com/INSaFLU (code).

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.

Long-term monitoring of SARS-CoV-2 seroprevalence and variants in Ethiopia provides prediction for immunity and cross-immunity

Under-reporting of COVID-19 and the limited information about circulating SARS-CoV-2 variants remain major challenges for many African countries. We analyzed SARS-CoV-2 infection dynamics in Addis Ababa and Jimma, Ethiopia, focusing on reinfection, immunity, and vaccination effects. We conducted an antibody serology study spanning August 2020 to July 2022 with five rounds of data collection across a population of 4723, sequenced PCR-test positive samples, used available test positivity rates, and constructed two mathematical models integrating this data. A multivariant model explores variant dynamics identifying wildtype, alpha, delta, and omicron BA.4/5 as key variants in the study population, and cross-immunity between variants, revealing risk reductions between 24% and 69%. An antibody-level model predicts slow decay leading to sustained high antibody levels. Retrospectively, increased early vaccination might have substantially reduced infections during the delta and omicron waves in the considered group of individuals, though further vaccination now seems less impactful.

Monkeypox virus genomic accordion strategies

The 2023 monkeypox (mpox) epidemic was caused by a subclade IIb descendant of a monkeypox virus (MPXV) lineage traced back to Nigeria in 1971. Person-to-person transmission appears higher than for clade I or subclade IIa MPXV, possibly caused by genomic changes in subclade IIb MPXV. Key genomic changes could occur in the genome's low-complexity regions (LCRs), which are challenging to sequence and are often dismissed as uninformative. Here, using a combination of highly sensitive techniques, we determine a high-quality MPXV genome sequence of a representative of the current epidemic with LCRs resolved at unprecedented accuracy. This reveals significant variation in short tandem repeats within LCRs. We demonstrate that LCR entropy in the MPXV genome is significantly higher than that of single-nucleotide polymorphisms (SNPs) and that LCRs are not randomly distributed. In silico analyses indicate that expression, translation, stability, or function of MPXV orthologous poxvirus genes (OPGs), including OPG153, OPG204, and OPG208, could be affected in a manner consistent with the established "genomic accordion" evolutionary strategies of orthopoxviruses. We posit that genomic studies focusing on phenotypic MPXV differences should consider LCR variability.

Longitudinal dynamics of the nasopharyngal microbiome in response to SARS-CoV-2 Omicron variant and HIV infection in Kenyan women and their infants

The nasopharynx and its microbiota are implicated in respiratory health and disease. The interplay between viral infection and the nasopharyngeal microbiome is an area of increased interest and of clinical relevance. The impact of SARS-CoV-2, the etiological agent of the Coronavirus Disease 2019 (COVID-19) pandemic, on the nasopharyngeal microbiome, particularly among individuals living with HIV, is not fully characterized. Here we describe the nasopharyngeal microbiome before, during and after SARS-CoV-2 infection in a longitudinal cohort of Kenyan women (21 living with HIV and 14 HIV-uninfected) and their infants (18 HIV-exposed, uninfected and 18 HIV-unexposed, uninfected), followed between September 2021 through March 2022. We show using genomic epidemiology that mother and infant dyads were infected with the same strain of the SARS-CoV-2 Omicron variant that spread rapidly across Kenya. Additionally, we used metagenomic sequencing to characterize the nasopharyngeal microbiome of 20 women and infants infected with SARS-CoV-2, 6 infants negative for SARS-CoV-2 but experiencing respiratory symptoms, and 34 timepoint matched SARS-CoV-2 negative mothers and infants. Since individuals were sampled longitudinally before and after SARS-CoV-2 infection, we could characterize the short- and long-term impact of SARS-CoV-2 infection on the nasopharyngeal microbiome. We found that mothers and infants had significantly different microbiome composition and bacterial load (p-values <.0001). However, in both mothers and infants, the nasopharyngeal microbiome did not differ before and after SARS-CoV-2 infection, regardless of HIV-exposure status. Our results indicate that the nasopharyngeal microbiome is resilient to SARS-CoV-2 infection and was not significantly modified by HIV.

Lessons learned: overcoming common challenges in reconstructing the SARS-CoV-2 genome from short-read sequencing data via CoVpipe2

Background

Accurate genome sequences form the basis for genomic surveillance programs, the added value of which was impressively demonstrated during the COVID-19 pandemic by tracing transmission chains, discovering new viral lineages and mutations, and assessing them for infectiousness and resistance to available treatments. Amplicon strategies employing Illumina sequencing have become widely established for variant detection and reference-based reconstruction of SARS-CoV-2 genomes, and are routine bioinformatics tasks. Yet, specific challenges arise when analyzing amplicon data, for example, when crucial and even lineage-determining mutations occur near primer sites.

Methods

We present CoVpipe2, a bioinformatics workflow developed at the Public Health Institute of Germany to reconstruct SARS-CoV-2 genomes based on short-read sequencing data accurately. The decisive factor here is the reliable, accurate, and rapid reconstruction of genomes, considering the specifics of the used sequencing protocol. Besides fundamental tasks like quality control, mapping, variant calling, and consensus generation, we also implemented additional features to ease the detection of mixed samples and recombinants.

Results

We highlight common pitfalls in primer clipping, detecting heterozygote variants, and dealing with low-coverage regions and deletions. We introduce CoVpipe2 to address the above challenges and have compared and successfully validated the pipeline against selected publicly available benchmark datasets. CoVpipe2 features high usability, reproducibility, and a modular design that specifically addresses the characteristics of short-read amplicon protocols but can also be used for whole-genome short-read sequencing data.

Conclusions

CoVpipe2 has seen multiple improvement cycles and is continuously maintained alongside frequently updated primer schemes and new developments in the scientific community. Our pipeline is easy to set up and use and can serve as a blueprint for other pathogens in the future due to its flexibility and modularity, providing a long-term perspective for continuous support. CoVpipe2 is written in Nextflow and is freely accessible from \href{https://github.com/rki-mf1/CoVpipe2}{github.com/rki-mf1/CoVpipe2} under the GPL3 license.

Systematic SARS-CoV-2 S-gene sequencing in wastewater samples enables early lineage detection and uncovers rare mutations in Portugal

As the COVID-19 pandemic reached its peak, many countries implemented genomic surveillance systems to track the evolution and transmission of SARS-CoV-2. Transition from the pandemic to the endemic phase prioritized alternative testing strategies to maintain effective epidemic surveillance at the population level, with less intensive sequencing efforts. One such promising approach was Wastewater-Based Surveillance (WBS), which offers non-invasive, cost-effective means for analysing virus trends at the sewershed level. From 2020 onwards, wastewater has been recognized as an instrumental source of information for public health, with national and international authorities exploring options to implement national wastewater surveillance systems and increasingly relying on WBS as early warning of potential pathogen outbreaks. In Portugal, several pioneer projects joined the academia, water utilities and Public Administration around WBS. To validate WBS as an effective genomic surveillance strategy, it is crucial to collect long term performance data. In this work, we present one year of systematic SARS-CoV-2 wastewater surveillance in Portugal, representing 35 % of the mainland population. We employed two complementary methods for lineage determination - allelic discrimination by RT-PCR and S-gene sequencing. This combination allowed us to monitor variant evolution in near-real-time and identify low-frequency mutations. Over the course of this year-long study, spanning from May 2022 to April 2023, we successfully tracked the dominant Omicron sub-lineages, their progression and evolution, which aligned with concurrent clinical surveillance data. Our results underscore the effectiveness of WBS as a tracking system for virus variants, with the ability to unveil mutations undetected via massive sequencing of clinical samples from Portugal, demonstrating the ability of WBS to uncover new mutations and detect rare genetic variants. Our findings emphasize that knowledge of the genetic diversity of SARS-CoV-2 at the population level can be extended far beyond via the combination of routine clinical genomic surveillance with wastewater sequencing and genotyping.

Scalable search of massively pooled nucleic acid samples enabled by a molecular database query language

The surge in nucleic acid analytics requires scalable storage and retrieval systems akin to electronic databases used to organize digital data. Such a system could transform disease diagnosis, ecological preservation, and molecular surveillance of biothreats. Current storage systems use individual containers for nucleic acid samples, requiring single-sample retrieval that falls short compared with digital databases that allow complex and combinatorial data retrieval on aggregated data. Here, we leverage protective microcapsules with combinatorial DNA labeling that enables arbitrary retrieval on pooled biosamples analogous to Structured Query Languages. Ninety-six encapsulated pooled mock SARS-CoV-2 genomic samples barcoded with patient metadata are used to demonstrate queries with simultaneous matches to sample collection date ranges, locations, and patient health statuses, illustrating how such flexible queries can be used to yield immunological or epidemiological insights. The approach applies to any biosample database labeled with orthogonal barcodes, enabling complex post-hoc analysis, for example, to study global biothreat epidemiology.

Metagenomic Analysis for Diagnosis of Hemorrhagic Fever in Minas Gerais, Brazil

Viral hemorrhagic fever poses a significant public health challenge due to its severe clinical presentation and high mortality rate. The diagnostic process is hindered by similarity of symptoms across different diseases and the broad spectrum of pathogens that can cause hemorrhagic fever. In this study, we applied viral metagenomic analysis to 43 serum samples collected by the Public Health Laboratory (Fundação Ezequiel Dias, FUNED) in Minas Gerais State, Brazil, from patients diagnosed with hemorrhagic fever who had tested negative for the standard local hemorrhagic disease testing panel. This panel includes tests for Dengue virus (DENV) IgM, Zika virus IgM, Chikungunya virus IgM, yellow fever IgM, Hantavirus IgM, Rickettsia rickettsii IgM/IgG, and Leptospira interrogans IgM, in addition to respective molecular tests for these infectious agents. The samples were grouped into 18 pools according to geographic origin and analyzed through next-generation sequencing on the NextSeq 2000 platform. Bioinformatic analysis revealed a prevalent occurrence of commensal viruses across all pools, but, notably, a significant number of reads corresponding to the DENV serotype 2 were identified in one specific pool. Further verification via real-time PCR confirmed the presence of DENV-2 RNA in an index case involving an oncology patient with hemorrhagic fever who had initially tested negative for anti-DENV IgM antibodies, thereby excluding this sample from initial molecular testing. The complete DENV-2 genome isolated from this patient was taxonomically classified within the cosmopolitan genotype that was recently introduced into Brazil. These findings highlight the critical role of considering the patient's clinical condition when deciding upon the most appropriate testing procedures. Additionally, this study showcases the potential of viral metagenomics in pinpointing the viral agents behind hemorrhagic diseases. Future research is needed to assess the practicality of incorporating metagenomics into standard viral diagnostic protocols.

Detection of clade 2.3.4.4b highly pathogenic H5N1 influenza virus in New York City

Highly pathogenic avian influenza viruses of the H5N1 clade 2.3.4.4b arrived in North America in the winter of 2021/2022. These viruses have spread across the Americas causing morbidity and mortality in both wild and domestic birds as well as some mammalian species, including cattle. Many surveillance programs in wildlife as well as commercial poultry operations have detected these viruses. Here we conducted surveillance of avian species in the urban environment in New York City. We detected highly pathogenic H5N1 viruses in six samples from four different bird species and performed full genome sequencing. Sequence analysis showed the presence of multiple different genotypes. Our work highlights that the interface between animals and humans that may give rise to zoonotic infections or even pandemics is not limited to rural environments and commercial poultry operations but extends into the heart of our urban centers.

Sequence signatures within the genome of SARS-CoV-2 can be used to predict host source

We conducted an in silico analysis to better understand the potential factors impacting host adaptation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in white-tailed deer, humans, and mink due to the strong evidence of sustained transmission within these hosts. Classification models trained on single nucleotide and amino acid differences between samples effectively identified white-tailed deer-, human-, and mink-derived SARS-CoV-2. For example, the balanced accuracy score of Extremely Randomized Trees classifiers was 0.984 ± 0.006. Eighty-eight commonly identified predictive mutations are found at sites under strong positive and negative selective pressure. A large fraction of sites under selection (86.9%) or identified by machine learning (87.1%) are found in genes other than the spike. Some locations encoded by these gene regions are predicted to be B- and T-cell epitopes or are implicated in modulating the immune response suggesting that host adaptation may involve the evasion of the host immune system, modulation of the class-I major-histocompatibility complex, and the diminished recognition of immune epitopes by CD8+ T cells. Our selection and machine learning analysis also identified that silent mutations, such as C7303T and C9430T, play an important role in discriminating deer-derived samples across multiple clades. Finally, our investigation into the origin of the B.1.641 lineage from white-tailed deer in Canada discovered an additional human sequence from Michigan related to the B.1.641 lineage sampled near the emergence of this lineage. These findings demonstrate that machine-learning approaches can be used in combination with evolutionary genomics to identify factors possibly involved in the cross-species transmission of viruses and the emergence of novel viral lineages.IMPORTANCESevere acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a highly transmissible virus capable of infecting and establishing itself in human and wildlife populations, such as white-tailed deer. This fact highlights the importance of developing novel ways to identify genetic factors that contribute to its spread and adaptation to new host species. This is especially important since these populations can serve as reservoirs that potentially facilitate the re-introduction of new variants into human populations. In this study, we apply machine learning and phylogenetic methods to uncover biomarkers of SARS-CoV-2 adaptation in mink and white-tailed deer. We find evidence demonstrating that both non-synonymous and silent mutations can be used to differentiate animal-derived sequences from human-derived ones and each other. This evidence also suggests that host adaptation involves the evasion of the immune system and the suppression of antigen presentation. Finally, the methods developed here are general and can be used to investigate host adaptation in viruses other than SARS-CoV-2.

Detection of SARS-COV-2 variants and their proportions in wastewater samples using next-generation sequencing in Finland

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) variants may have different characteristics, e.g., in transmission, mortality, and the effectiveness of vaccines, indicating the importance of variant detection at the population level. Wastewater-based surveillance of SARS-CoV-2 RNA fragments has been shown to be an effective way to monitor the COVID-19 pandemic at the population level. Wastewater is a complex sample matrix affected by environmental factors and PCR inhibitors, causing insufficient coverage in sequencing, for example. Subsequently, results where part of the genome does not have sufficient coverage are not uncommon. To identify variants and their proportions in wastewater over time, we utilized next-generation sequencing with the ARTIC Network's primer set and bioinformatics pipeline to evaluate the presence of variants in partial genome data. Based on the wastewater data from November 2021 to February 2022, the Delta variant was dominant until mid-December in Helsinki, Finland's capital, and thereafter in late December 2022 Omicron became the most common variant. At the same time, the Omicron variant of SARS-CoV-2 outcompeted the previous Delta variant in Finland in new COVID-19 cases. The SARS-CoV-2 variant findings from wastewater are in agreement with the variant information obtained from the patient samples when visually comparing trends in the sewerage network area. This indicates that the sequencing of wastewater is an effective way to monitor temporal and spatial trends of SARS-CoV-2 variants at the population level.

Robust and sensitive amplicon-based whole-genome sequencing assay of respiratory syncytial virus subtype A and B

Prevention of respiratory syncytial virus (RSV) infection is now a global health priority, with a long-acting monoclonal antibody and two RSV vaccines recently licenced for clinical use. Most licenced and candidate interventions target the RSV fusion (RSV-F) protein. New interventions may be associated with the spread of mutations, reducing susceptibility to antibody neutralization in RSV-F. There is a need for ongoing longitudinal global surveillance of circulating RSV strains. To achieve this large-scale genomic surveillance, a reliable, high-throughput RSV sequencing assay is required. Here we report an improved high-throughput RSV whole-genome sequencing (WGS) assay performed directly on clinical samples without additional enrichment, using a 4-primer-pool, short-amplicon PCR-tiling approach that is suitable for short-read sequencing platforms. Using upper respiratory tract (URT) RSV-positive clinical samples obtained from a sentinel network of primary care providers and from hospital patients (29.7% and 70.2%, respectively; n = 1,037), collected over the period 2019 to 2023, this assay had a threshold of approximately 4 × 103 to 8 × 103 copies/mL (RSV-B and RSV-A sub-types, respectively) as the lowest amount of virus needed in the sample to achieve >96% of whole-genome coverage at a high-quality level. Using a Ct value of 31 as an empirical cut-off, the overall assay success rate of obtaining >90% genome coverage at a read depth minimum of 20 was 96.83% for clinical specimens successfully sequenced from a total of 1,071. The RSV WGS approach described in this study has increased sensitivity compared to previous approaches and can be applied to clinical specimens without the requirement for enrichment. The updated approach produces sequences of high quality consistently and cost-effectively, suitable for implementation to underpin national programs for the surveillance of RSV genomic variation.

IMPORTANCE

In this paper, we report an improved high-throughput respiratory syncytial virus (RSV) whole-genome sequencing (WGS) assay performed directly on clinical samples, using a 4-primer-pool, short-amplicon PCR-tiling approach that is suitable for short-read sequencing platforms. The RSV WGS approach described in this study has increased sensitivity compared to previous approaches and can be applied to clinical specimens without the requirement for enrichment. The updated approach produces sequences of high quality consistently and cost-effectively, suitable for implementation to underpin national and global programs for the surveillance of RSV genomic variation. The quality of sequence produced is essential for preparedness for new interventions in monitoring antigenic escape, where a single point mutation might lead to a reduction in antibody binding effectiveness and neutralizing activity, or indeed in the monitoring of retaining susceptibility to neutralization by existing and new interventions.

Evolution of SARS-CoV-2 in the Rhine-Neckar/Heidelberg Region 01/2021 - 07/2023

In January 2021, the monitoring of circulating variants of SARS-CoV-2 was initiated in Germany under the Corona Surveillance Act, which was discontinued after July 2023. This initiative aimed to enhance pandemic containment, as specific amino acid changes, particularly in the spike protein, were associated with increased transmission and reduced vaccine efficacy. Our group conducted whole genome sequencing using the ARTIC protocol (currently V4) on Illumina's NextSeq 500 platform (and, starting in May 2023, on the MiSeq DX platform) for SARS-CoV-2 positive specimen from patients at Heidelberg University Hospital, associated hospitals, and the public health office in the Rhine-Neckar/Heidelberg region. In total, we sequenced 26,795 SARS-CoV-2-positive samples between January 2021 and July 2023. Valid sequences, meeting the requirements for upload to the German electronic sequencing data hub (DESH) operated by the Robert Koch Institute (RKI), were determined for 24,852 samples, and the lineage/clade could be identified for 25,912 samples. The year 2021 witnessed significant dynamics in the circulating variants in the Rhine-Neckar/Heidelberg region, including A.27.RN, followed by the emergence of B.1.1.7 (Alpha), subsequently displaced by B.1.617.2 (Delta), and the initial occurrences of B.1.1.529 (Omicron). By January 2022, B.1.1.529 had superseded B.1.617.2, dominating with over 90%. The years 2022 and 2023 were then characterized by the dominance of B.1.1.529 and its sublineages, particularly BA.5 and BA.2, and more recently, the emergence of recombinant variants like XBB.1.5. Since the global dominance of B.1.617.2, the identified variant distribution in our local study, apart from a time delay in the spread of new variants, can be considered largely representative of the global distribution. om a time delay in the spread of new variants, can be considered largely representative of the global distribution.

An Amplicon-Based Application for the Whole-Genome Sequencing of GI-19 Lineage Infectious Bronchitis Virus Directly from Clinical Samples

Infectious bronchitis virus (IBV) causes a highly contagious respiratory disease in chickens, leading to significant economic losses in the poultry industry worldwide. IBV exhibits a high mutation rate, resulting in the continuous emergence of new variants and strains. A complete genome analysis of IBV is crucial for understanding its characteristics. However, it is challenging to obtain whole-genome sequences from IBV-infected clinical samples due to the low abundance of IBV relative to the host genome. Here, we present a novel approach employing next-generation sequencing (NGS) to directly sequence the complete genome of IBV. Through in silico analysis, six primer pairs were designed to match various genotypes, including the GI-19 lineage of IBV. The primer sets successfully amplified six overlapping fragments by long-range PCR and the size of the amplicons ranged from 3.7 to 6.4 kb, resulting in full coverage of the IBV genome. Furthermore, utilizing Illumina sequencing, we obtained the complete genome sequences of two strains belonging to the GI-19 lineage (QX genotype) from clinical samples, with 100% coverage rates, over 1000 × mean depth coverage, and a high percentage of mapped reads to the reference genomes (96.63% and 97.66%). The reported method significantly improves the whole-genome sequencing of IBVs from clinical samples; thus, it can improve understanding of the epidemiology and evolution of IBVs.

Whole-Genome Sequencing and Mutation Analyses of SARS-CoV-2 Isolates from Indonesia

The SARS-CoV-2 infection that caused the COVID-19 pandemic has become a significant public health concern. New variants with distinct mutations have emerged, potentially impacting its infectivity, immune evasion capacity, and vaccine response. A whole-genome sequencing study of 292 SARS-CoV-2 isolates collected from selected regions of Indonesia between January and October 2021 was performed to identify the distribution of SARS-CoV-2 variants and common mutations in Indonesia. During January-April 2021, Indonesian lineages B.1.466.2 and B.1.470 dominated, but from May 2021, Delta's AY.23 lineage outcompeted them. An analysis of 7515 published sequences from January 2021 to June 2022 revealed a decline in Delta in November 2021, followed by the emergence of Omicron variants in December 2021. We identified C241T (5'UTR), P314L (NSP12b), F106F (NSP3), and D614G (Spike) mutations in all sequences. The other common substitutions included P681R (76.4%) and T478K (60%) in Spike, D377Y in Nucleocapsid (61%), and I82T in Membrane (60%) proteins. Breakthrough infection and prolonged viral shedding cases were associated with Delta variants carrying the Spike T19R, G142D, L452R, T478K, D614G, P681R, D950N, and V1264L mutations. The dynamic of SARS-CoV-2 variants in Indonesia highlights the importance of continuous genomic surveillance in monitoring and identifying potential strains leading to disease outbreaks.

Rapid response to a COVID-19 outbreak at a nightclub in Kagoshima prefecture, Japan, in the early phase of the COVID-19 pandemic, June and July 2020: A descriptive epidemiological study

INTRODUCTION

During COVID-19 pandemic in Japan, nightclubs were identified as high-risk locations for COVID-19 outbreaks, but an outbreak investigation in this setting is challenging because of the anonymous and opportunistic nature of interactions.

METHODS

The joint rapid response team collected epidemiological data, conducted descriptive epidemiology to determine the characteristics of cases associated with the nightclub, and implemented countermeasures. Polymerase chain reaction (PCR) tests were performed by the Local Institute of Public Health, Kagoshima University, and several commercial laboratories.

RESULTS

Between June 15 and July 20, 2020, 121 individuals tested positive for SARS-CoV-2 (59 confirmed and 62 asymptomatic) of whom 8 were nightclub staff who had no travel history of outside Kagoshima, 66 were guests, and 47 were subsequent contacts. The median age was 32 years (interquartile range: 24-43 years). One individual showed severe symptoms but there were no fatal. The epidemic curve showed one peak on June 30 and July 1 with a limited number of cases subsequently. Of the 121 cases, 116 and 5 were in individuals living in and outside Kagoshima Prefecture, respectively. Haplotype network analysis showed 5 genome-wide single-nucleotide variants between the isolates before and during this outbreak.

CONCLUSIONS

There is a possibility that unidentified guests from outside Kagoshima Prefecture could infect staff who could subsequently spread the virus to guests and other staff, who were mainly a younger population. The rapid outbreak response enabled onward transmission in the community to be minimized. This outbreak investigation could provide insights for effective responses to challenging situations in future pandemic.

Immune and behavioral correlates of protection against symptomatic post-vaccination SARS-CoV-2 infection

Introduction

We sought to determine pre-infection correlates of protection against SARS-CoV-2 post-vaccine inzfections (PVI) acquired during the first Omicron wave in the United States.

Methods

Serum and saliva samples from 176 vaccinated adults were collected from October to December of 2021, immediately before the Omicron wave, and assessed for SARS-CoV-2 Spike-specific IgG and IgA binding antibodies (bAb). Sera were also assessed for bAb using commercial assays, and for neutralization activity against several SARS-CoV-2 variants. PVI duration and severity, as well as risk and precautionary behaviors, were assessed by questionnaires.

Results

Serum anti-Spike IgG levels assessed by research assay, neutralization titers against Omicron subvariants, and low home risk scores correlated with protection against PVIs after multivariable regression analysis. Commercial assays did not perform as well as research assay, likely due to their lower dynamic range.

Discussion

In the 32 participants that developed PVI, anti-Spike IgG bAbs correlated with lower disease severity and shorter duration of illness.

Natural variation in neuraminidase activity influences the evolutionary potential of the seasonal H1N1 lineage hemagglutinin

The antigenic evolution of the influenza A virus hemagglutinin (HA) gene poses a major challenge for the development of vaccines capable of eliciting long-term protection. Prior efforts to understand the mechanisms that govern viral antigenic evolution mainly focus on HA in isolation, ignoring the fact that HA must act in concert with the viral neuraminidase (NA) during replication and spread. Numerous studies have demonstrated that the degree to which the receptor binding avidity of HA and receptor cleaving activity of NA are balanced with each other influences overall viral fitness. We recently showed that changes in NA activity can significantly alter the mutational fitness landscape of HA in the context of a lab-adapted virus strain. Here, we test whether natural variation in relative NA activity can influence the evolutionary potential of HA in the context of the seasonal H1N1 lineage (pdmH1N1) that has circulated in humans since the 2009 pandemic. We observed substantial variation in the relative activities of NA proteins encoded by a panel of H1N1 vaccine strains isolated between 2009 and 2019. We comprehensively assessed the effect of NA background on the HA mutational fitness landscape in the circulating pdmH1N1 lineage using deep mutational scanning and observed pronounced epistasis between NA and residues in or near the receptor binding site of HA. To determine whether NA variation could influence the antigenic evolution of HA, we performed neutralizing antibody selection experiments using a panel of monoclonal antibodies targeting different HA epitopes. We found that the specific antibody escape profiles of HA were highly contingent upon NA background. Overall, our results indicate that natural variation in NA activity plays a significant role in governing the evolutionary potential of HA in the currently circulating pdmH1N1 lineage.

SARS-CoV-2 evolution during prolonged infection in immunocompromised patients

Prolonged infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in immunocompromised patients provides an opportunity for viral evolution, potentially leading to the generation of new pathogenic variants. To investigate the pathways of viral evolution, we carried out a study on five patients experiencing prolonged SARS-CoV-2 infection (quantitative polymerase chain reaction-positive for 79-203 days) who were immunocompromised due to treatment for lymphoma or solid organ transplantation. For each timepoint analyzed, we generated at least two independent viral genome sequences to assess the heterogeneity and control for sequencing error. Four of the five patients likely had prolonged infection; the fifth apparently experienced a reinfection. The rates of accumulation of substitutions in the viral genome per day were higher in hospitalized patients with prolonged infection than those estimated for the community background. The spike coding region accumulated a significantly greater number of unique mutations than other viral coding regions, and the mutation density was higher. Two patients were treated with monoclonal antibodies (bebtelovimab and sotrovimab); by the next sampled timepoint, each virus population showed substitutions associated with monoclonal antibody resistance as the dominant forms (spike K444N and spike E340D). All patients received remdesivir, but remdesivir-resistant substitutions were not detected. These data thus help elucidate the trends of emergence, evolution, and selection of mutational variants within long-term infected immunocompromised individuals.

IMPORTANCE

SARS-CoV-2 is responsible for a global pandemic, driven in part by the emergence of new viral variants. Where do these new variants come from? One model is that long-term viral persistence in infected individuals allows for viral evolution in response to host pressures, resulting in viruses more likely to replicate efficiently in humans. In this study, we characterize replication in several hospitalized and long-term infected individuals, documenting efficient pathways of viral evolution.

Remote surveillance and detection of SARS-CoV-2 transmission among household members in King County, Washington

BACKGROUND

Early during the COVID-19 pandemic, it was important to better understand transmission dynamics of SARS-CoV-2, the virus that causes COVID-19. Household contacts of infected individuals are particularly at risk for infection, but delays in contact tracing, delays in testing contacts, and isolation and quarantine posed challenges to accurately capturing secondary household cases.

METHODS

In this study, 346 households in the Seattle region were provided with respiratory specimen collection kits and remotely monitored using web-based surveys for respiratory illness symptoms weekly between October 1, 2020, and June 20, 2021. Symptomatic participants collected respiratory specimens at symptom onset and mailed specimens to the central laboratory in Seattle. Specimens were tested for SARS-CoV-2 using RT-PCR with whole genome sequencing attempted when positive. SARS-CoV-2-infected individuals were notified, and their household contacts submitted specimens every 2 days for 14 days.

RESULTS

In total, 1371 participants collected 2029 specimens that were tested; 16 individuals (1.2%) within 6 households tested positive for SARS-CoV-2 during the study period. Full genome sequences were generated from 11 individuals within 4 households. Very little genetic variation was found among SARS-CoV-2 viruses sequenced from different individuals in the same household, supporting transmission within the household.

CONCLUSIONS

This study indicates web-based surveillance of respiratory symptoms, combined with rapid and longitudinal specimen collection and remote contact tracing, provides a viable strategy to monitor households and detect household transmission of SARS-CoV-2.

TRIAL REGISTRATION IDENTIFIER

NCT04141930, Date of registration 28/10/2019.

Genome sequence of bubaline herpesvirus-1 associated with pustular vulvovaginitis in Indian water buffalo

We report here the genome sequence of a bubaline herpesvirus 1 isolated from Indian water buffalo. The bubaline herpesvirus 1 strain S102_1 was isolated in 2021 from a Murrah buffalo heifer with clinical presentation of pustular vulvovaginitis.

Recommendations for Uniform Variant Calling of SARS-CoV-2 Genome Sequence across Bioinformatic Workflows

Genomic sequencing of clinical samples to identify emerging variants of SARS-CoV-2 has been a key public health tool for curbing the spread of the virus. As a result, an unprecedented number of SARS-CoV-2 genomes were sequenced during the COVID-19 pandemic, which allowed for rapid identification of genetic variants, enabling the timely design and testing of therapies and deployment of new vaccine formulations to combat the new variants. However, despite the technological advances of deep sequencing, the analysis of the raw sequence data generated globally is neither standardized nor consistent, leading to vastly disparate sequences that may impact identification of variants. Here, we show that for both Illumina and Oxford Nanopore sequencing platforms, downstream bioinformatic protocols used by industry, government, and academic groups resulted in different virus sequences from same sample. These bioinformatic workflows produced consensus genomes with differences in single nucleotide polymorphisms, inclusion and exclusion of insertions, and/or deletions, despite using the same raw sequence as input datasets. Here, we compared and characterized such discrepancies and propose a specific suite of parameters and protocols that should be adopted across the field. Consistent results from bioinformatic workflows are fundamental to SARS-CoV-2 and future pathogen surveillance efforts, including pandemic preparation, to allow for a data-driven and timely public health response.

VIPERA: Viral Intra-Patient Evolution Reporting and Analysis

Viral mutations within patients nurture the adaptive potential of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) during chronic infections, which are a potential source of variants of concern. However, there is no integrated framework for the evolutionary analysis of intra-patient SARS-CoV-2 serial samples. Herein, we describe Viral Intra-Patient Evolution Reporting and Analysis (VIPERA), a new software that integrates the evaluation of the intra-patient ancestry of SARS-CoV-2 sequences with the analysis of evolutionary trajectories of serial sequences from the same viral infection. We have validated it using positive and negative control datasets and have successfully applied it to a new case, which revealed population dynamics and evidence of adaptive evolution. VIPERA is available under a free software license at https://github.com/PathoGenOmics-Lab/VIPERA.

High pathogenicity avian influenza A (H5N1) clade 2.3.4.4b virus infection in a captive Tibetan black bear (Ursus thibetanus): investigations based on paraffin-embedded tissues, France, 2022

High pathogenicity avian influenza viruses (HPAIVs) H5Nx of clade 2.3.4.4b have been circulating increasingly in both wild and domestic birds in recent years. In turn, this has led to an increase in the number of spillover events affecting mammals. In November 2022, an HPAIV H5N1 caused an outbreak in a zoological park in the south of France, resulting in the death of a Tibetan black bear (Ursus thibetanus) and several captive and wild bird species. We detected the virus in various tissues of the bear and a wild black-headed gull (Chroicocephalus ridibundus) found dead in its enclosure using histopathology, two different in situ detection techniques, and next-generation sequencing, all performed on formalin-fixed paraffin-embedded tissues. Phylogenetic analysis performed on the hemagglutinin gene segment showed that bear and gull strains shared 99.998% genetic identity, making the bird strain the closest related strain. We detected the PB2 E627K mutation in minute quantities in the gull, whereas it predominated in the bear, which suggests that this mammalian adaptation marker was selected during the bear infection. Our results provide the first molecular and histopathological characterization of an H5N1 virus infection in this bear species.

IMPORTANCE

Avian influenza viruses are able to cross the species barrier between birds and mammals because of their high genetic diversity and mutation rate. Using formalin-fixed paraffin-embedded tissues, we were able to investigate a Tibetan black bear's infection by a high pathogenicity H5N1 avian influenza virus at the molecular, phylogenetic, and histological levels. Our results highlight the importance of virological surveillance programs in mammals and the importance of raising awareness among veterinarians and zookeepers of the clinical presentations associated with H5Nx virus infection in mammals.

Evolutionary trajectory of diverse SARS-CoV-2 variants at the beginning of COVID-19 outbreak

Despite extensive scientific efforts directed toward the evolutionary trajectory of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans at the beginning of the COVID-19 epidemic, it remains unclear how the virus jumped into and evolved in humans so far. Herein, we recruited almost all adult coronavirus disease 2019 (COVID-19) cases appeared locally or imported from abroad during the first 8 months of the outbreak in Shanghai. From these patients, SARS-CoV-2 genomes occupying the important phylogenetic positions in the virus phylogeny were recovered. Phylogenetic and mutational landscape analyses of viral genomes recovered here and those collected in and outside of China revealed that all known SARS-CoV-2 variants exhibited the evolutionary continuity despite the co-circulation of multiple lineages during the early period of the epidemic. Various mutations have driven the rapid SARS-CoV-2 diversification, and some of them favor its better adaptation and circulation in humans, which may have determined the waxing and waning of various lineages.

PySNV for complex intra-host variation detection

MOTIVATION

Intra-host variants refer to genetic variations or mutations that occur within an individual host organism. These variants are typically studied in the context of viruses, bacteria, or other pathogens to understand the evolution of pathogens. Moreover, intra-host variants are also explored in the field of tumor biology and mitochondrial biology to characterize somatic mutations and inherited heteroplasmic mutations. Intra-host variants can involve long insertions, deletions, and combinations of different mutation types, which poses challenges in their identification. The performance of current methods in detecting of complex intra-host variants is unknown.

RESULTS

First, we simulated a dataset comprising 10 samples with 1869 intra-host variants involving various mutation patterns and benchmarked current variant detection software. The results indicated that though current software can detect most variants with F1-scores between 0.76 and 0.97, their performance in detecting long indels and low frequency variants was limited. Thus, we developed a new software, PySNV, for the detection of complex intra-host variations. On the simulated dataset, PySNV successfully detected 1863 variant cases (F1-score: 0.99) and exhibited the highest Pearson correlation coefficient (PCC: 0.99) to the ground truth in predicting variant frequencies. The results demonstrated that PySNV delivered promising performance even for long indels and low frequency variants, while maintaining computational speed comparable to other methods. Finally, we tested its performance on SARS-CoV-2 replicate sequencing data and found that it reported 21% more variants compared to LoFreq, the best-performing benchmarked software, while showing higher consistency (62% over 54%) within replicates. The discrepancies mostly exist in low-depth regions and low frequency variants.

AVAILABILITY AND IMPLEMENTATION

https://github.com/bnuLyndon/PySNV/.

Spike-Seq: An amplicon-based high-throughput sequencing approach for the sensitive detection and characterization of SARS-CoV-2 genetic variations in environmental samples

Ever since the outbreak of COVID-19 disease in Wuhan, China, different variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been identified. Wastewater-based epidemiology (WBE), an approach that has been successfully applied in numerous case studies worldwide, offers a cost-effective and rapid way for monitoring trends of SARS-Cov-2 in the community level without selection bias. Despite being a gold-standard procedure, WBE is a challenging approach due to the sample instability and the moderate efficiency of SARS-CoV-2 concentration in wastewater. In the present study, we introduce Spike-Seq, a custom amplicon-based approach for the S gene sequencing of SARS-CoV-2 in wastewater samples, which enables not only the accurate identification of the existing Spike-related genetic markers, but also the estimation of their frequency in the investigated samples. The implementation of Spike-Seq involves the combination of nested PCR-based assays that efficiently amplify the entire nucleotide sequence of the S gene and next-generation sequencing, which enables the variant detection and the estimation of their frequency. In the framework of the current work, Spike-Seq was performed to investigate the mutational profile of SARS-CoV-2 in samples from the Wastewater Treatment Plant (WWTP) of Athens, Greece, which originated from multiple timepoints, ranging from March 2021 until July 2022. Our findings demonstrate that Spike-Seq efficiently detected major genetic markers of B.1.1.7 (Alpha), B.1.617.2 (Delta) as well as B.1.1.529 (Omicron) variants in wastewater samples and provided their frequency levels, showing similar variant distributions with the published clinical data from the National Public Health organization. The presented approach can prove to be a useful tool for the detection of SARS-CoV-2 in challenging wastewater samples and the identification of the existing genetic variants of S gene.

Exploration of low-frequency allelic variants of SARS-CoV-2 genomes reveals coinfections in Mexico occurred during periods of VOCs turnover

A total of 14 973 alleles in 29 661 sequenced samples collected between March 2021 and January 2023 by the Mexican Consortium for Genomic Surveillance (CoViGen-Mex) and collaborators were used to construct a thorough map of mutations of the Mexican SARS-CoV-2 genomic landscape containing Intra-Patient Minor Allelic Variants (IPMAVs), which are low-frequency alleles not ordinarily present in a genomic consensus sequence. This additional information proved critical in identifying putative coinfecting variants included alongside the most common variants, B.1.1.222, B.1.1.519, and variants of concern (VOCs) Alpha, Gamma, Delta, and Omicron. A total of 379 coinfection events were recorded in the dataset (a rate of 1.28 %), resulting in the first such catalogue in Mexico. The most common putative coinfections occurred during the spread of Delta or after the introduction of Omicron BA.2 and its descendants. Coinfections occurred constantly during periods of variant turnover when more than one variant shared the same niche and high infection rate was observed, which was dependent on the local variants and time. Coinfections might occur at a higher frequency than customarily reported, but they are often ignored as only the consensus sequence is reported for lineage identification.

Molecular and phylogenetic characterization of the monkeypox outbreak in the South of Spain

Background and Aim

Until the May 2022 Monkeypox (MPXV) outbreak, which spread rapidly to many non-endemic countries, the virus was considered a viral zoonosis limited to some African countries. The Andalusian circuit of genomic surveillance was rapidly applied to characterize the MPXV outbreak in the South of Spain.

Methods

Whole genome sequencing was used to obtain the genomic profiles of samples collected across the south of Spain, representative of all the provinces of Andalusia. Phylogenetic analysis was used to study the relationship of the isolates and the available sequences of the 2022 outbreak.

Results

Whole genome sequencing of a total of 160 MPXV viruses from the different provinces that reported cases were obtained. Interestingly, we report the sequences of MPXV viruses obtained from two patients who died. While one of the isolates bore no noteworthy mutations that explain a potential heightened virulence, in another patient the second consecutive genome sequence, performed after the administration of tecovirimat, uncovered a mutation within the A0A7H0DN30 gene, known to be a prime target for tecovirimat in its Vaccinia counterpart. In general, a low number of mutations were observed in the sequences reported, which were very similar to the reference of the 2022 outbreak (OX044336), as expected from a DNA virus. The samples likely correspond to several introductions of the circulating MPXV viruses from the last outbreak. The virus sequenced from one of the two patients that died presented a mutation in a gene that bears potential connections to drug resistance. This mutation was absent in the initial sequencing before treatment.

Utilizing a university testing program to estimate relative effectiveness of monovalent COVID-19 mRNA booster vaccine versus two-dose primary series against symptomatic SARS-CoV-2 infection

Vaccine effectiveness (VE) studies utilizing the test-negative design are typically conducted in clinical settings, rather than community populations, leading to bias in VE estimates against mild disease and limited information on VE in healthy young adults. In a community-based university population, we utilized data from a large SARS-CoV-2 testing program to estimate relative VE of COVID-19 mRNA vaccine primary series and monovalent booster dose versus primary series only against symptomatic SARS-CoV-2 infection from September 2021 to July 2022. We used the test-negative design and logistic regression implemented via generalized estimating equations adjusted for age, calendar time, prior SARS-CoV-2 infection, and testing frequency (proxy for test-seeking behavior) to estimate relative VE. Analyses included 2,218 test-positive cases (59 % received monovalent booster dose) and 9,615 test-negative controls (62 %) from 9,066 individuals, with median age of 21 years, mostly students (71 %), White (56 %) or Asian (28 %), and with few comorbidities (3 %). More cases (23 %) than controls (6 %) had COVID-19-like illness. Estimated adjusted relative VE of primary series and monovalent booster dose versus primary series only against symptomatic SARS-CoV-2 infection was 40 % (95 % CI: 33-47 %) during the overall analysis period and 46 % (39-52 %) during the period of Omicron circulation. Relative VE was greater for those without versus those with prior SARS-CoV-2 infection (41 %, 34-48 % versus 33 %, 9 %-52 %, P < 0.001). Relative VE was also greater in the six months after receiving a booster dose (41 %, 33-47 %) compared to more than six months (27 %, 8-42 %), but this difference was not statistically significant (P = 0.06). In this relatively young and healthy adult population, an mRNA monovalent booster dose provided increased protection against symptomatic SARS-CoV-2 infection, overall and with the Omicron variant. University testing programs may be utilized for estimating VE in healthy young adults, a population that is not well-represented by routine VE studies.

Randomized controlled trial of molnupiravir SARS-CoV-2 viral and antibody response in at-risk adult outpatients

Viral clearance, antibody response and the mutagenic effect of molnupiravir has not been elucidated in at-risk populations. Non-hospitalised participants within 5 days of SARS-CoV-2 symptoms randomised to receive molnupiravir (n = 253) or Usual Care (n = 324) were recruited to study viral and antibody dynamics and the effect of molnupiravir on viral whole genome sequence from 1437 viral genomes. Molnupiravir accelerates viral load decline, but virus is detectable by Day 5 in most cases. At Day 14 (9 days post-treatment), molnupiravir is associated with significantly higher viral persistence and significantly lower anti-SARS-CoV-2 spike antibody titres compared to Usual Care. Serial sequencing reveals increased mutagenesis with molnupiravir treatment. Persistence of detectable viral RNA at Day 14 in the molnupiravir group is associated with higher transition mutations following treatment cessation. Viral viability at Day 14 is similar in both groups with post-molnupiravir treated samples cultured up to 9 days post cessation of treatment. The current 5-day molnupiravir course is too short. Longer courses should be tested to reduce the risk of potentially transmissible molnupiravir-mutated variants being generated. Trial registration: ISRCTN30448031.

Genomic surveillance of SARS-CoV-2 using long-range PCR primers

Introduction

Whole Genome Sequencing (WGS) of the SARS-CoV-2 virus is crucial in the surveillance of the COVID-19 pandemic. Several primer schemes have been developed to sequence nearly all of the ~30,000 nucleotide SARS-CoV-2 genome, using a multiplex PCR approach to amplify cDNA copies of the viral genomic RNA. Midnight primers and ARTIC V4.1 primers are the most popular primer schemes that can amplify segments of SARS-CoV-2 (400 bp and 1200 bp, respectively) tiled across the viral RNA genome. Mutations within primer binding sites and primer-primer interactions can result in amplicon dropouts and coverage bias, yielding low-quality genomes with 'Ns' inserted in the missing amplicon regions, causing inaccurate lineage assignments, and making it challenging to monitor lineage-specific mutations in Variants of Concern (VoCs).

Methods

In this study we used a set of seven long-range PCR primer pairs to sequence clinical isolates of SARS-CoV-2 on Oxford Nanopore sequencer. These long-range primers generate seven amplicons approximately 4500 bp that covered whole genome of SARS-CoV-2. One of these regions includes the full-length S-gene by using a set of flanking primers. We also evaluated the performance of these long-range primers with Midnight primers by sequencing 94 clinical isolates in a Nanopore flow cell.

Results and discussion

Using a small set of long-range primers to sequence SARS-CoV-2 genomes reduces the possibility of amplicon dropout and coverage bias. The key finding of this study is that long range primers can be used in single-molecule sequencing of RNA viruses in surveillance of emerging variants. We also show that by designing primers flanking the S-gene, we can obtain reliable identification of SARS-CoV-2 variants.

Improved detection of low-frequency within-host variants from deep sequencing: A case study with human papillomavirus

High-coverage sequencing allows the study of variants occurring at low frequencies within samples, but is susceptible to false-positives caused by sequencing error. Ion Torrent has a very low single nucleotide variant (SNV) error rate and has been employed for the majority of human papillomavirus (HPV) whole genome sequences. However, benchmarking of intrahost SNVs (iSNVs) has been challenging, partly due to limitations imposed by the HPV life cycle. We address this problem by deep sequencing three replicates for each of 31 samples of HPV type 18 (HPV18). Errors, defined as iSNVs observed in only one of three replicates, are dominated by C→T (G→A) changes, independently of trinucleotide context. True iSNVs, defined as those observed in all three replicates, instead show a more diverse SNV type distribution, with particularly elevated C→T rates in CCG context (CCG→CTG; CGG→CAG) and C→A rates in ACG context (ACG→AAG; CGT→CTT). Characterization of true iSNVs allowed us to develop two methods for detecting true variants: (1) VCFgenie, a dynamic binomial filtering tool which uses each variant's allele count and coverage instead of fixed frequency cut-offs; and (2) a machine learning binary classifier which trains eXtreme Gradient Boosting models on variant features such as quality and trinucleotide context. Each approach outperforms fixed-cut-off filtering of iSNVs, and performance is enhanced when both are used together. Our results provide improved methods for identifying true iSNVs in within-host applications across sequencing platforms, specifically using HPV18 as a case study.

CD4+ and CD8+ T cells are required to prevent SARS-CoV-2 persistence in the nasal compartment

SARS-CoV-2 is the causative agent of COVID-19 and continues to pose a significant public health threat throughout the world. Following SARS-CoV-2 infection, virus-specific CD4+ and CD8+ T cells are rapidly generated to form effector and memory cells and persist in the blood for several months. However, the contribution of T cells in controlling SARS-CoV-2 infection within the respiratory tract are not well understood. Using C57BL/6 mice infected with a naturally occurring SARS-CoV-2 variant (B.1.351), we evaluated the role of T cells in the upper and lower respiratory tract. Following infection, SARS-CoV-2-specific CD4+ and CD8+ T cells are recruited to the respiratory tract and a vast proportion secrete the cytotoxic molecule Granzyme B. Using antibodies to deplete T cells prior to infection, we found that CD4+ and CD8+ T cells play distinct roles in the upper and lower respiratory tract. In the lungs, T cells play a minimal role in viral control with viral clearance occurring in the absence of both CD4+ and CD8+ T cells through 28 days post-infection. In the nasal compartment, depletion of both CD4+ and CD8+ T cells, but not individually, results in persistent and culturable virus replicating in the nasal compartment through 28 days post-infection. Using in situ hybridization, we found that SARS-CoV-2 infection persisted in the nasal epithelial layer of tandem CD4+ and CD8+ T cell-depleted mice. Sequence analysis of virus isolates from persistently infected mice revealed mutations spanning across the genome, including a deletion in ORF6. Overall, our findings highlight the importance of T cells in controlling virus replication within the respiratory tract during SARS-CoV-2 infection.

The most exposed regions of SARS-CoV-2 structural proteins are subject to strong positive selection and gene overlap may locally modify this behavior

The SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) pandemic that emerged in 2019 has been an unprecedented event in international science, as it has been possible to sequence millions of genomes, tracking their evolution very closely. This has enabled various types of secondary analyses of these genomes, including the measurement of their sequence selection pressure. In this work, we have been able to measure the selective pressure of all the described SARS-CoV-2 genes, even analyzed by sequence regions, and we show how this type of analysis allows us to separate the genes between those subject to positive selection (usually those that code for surface proteins or those exposed to the host immune system) and those subject to negative selection because they require greater conservation of their structure and function. We have also seen that when another gene with an overlapping reading frame appears within a gene sequence, the overlapping sequence between the two genes evolves under a stronger purifying selection than the average of the non-overlapping regions of the main gene. We propose this type of analysis as a useful tool for locating and analyzing all the genes of a viral genome when an adequate number of sequences are available.IMPORTANCEWe have analyzed the selection pressure of all severe acute respiratory syndrome coronavirus 2 genes by means of the nonsynonymous (Ka) to synonymous (Ks) substitution rate. We found that protein-coding genes are exposed to strong positive selection, especially in the regions of interaction with other molecules (host receptor and genome of the virus itself). However, overlapping coding regions are more protected and show negative selection. This suggests that this measure could be used to study viral gene function as well as overlapping genes.

Within-host evolutionary dynamics and tissue compartmentalization during acute SARS-CoV-2 infection

The global evolution of SARS-CoV-2 depends in part upon the evolutionary dynamics within individual hosts with varying immune histories. To characterize the within-host evolution of acute SARS-CoV-2 infection, we sequenced saliva and nasal samples collected daily from vaccinated and unvaccinated individuals early during infection. We show that longitudinal sampling facilitates high-confidence genetic variant detection and reveals evolutionary dynamics missed by less-frequent sampling strategies. Within-host dynamics in both unvaccinated and vaccinated individuals appeared largely stochastic; however, in rare cases, minor genetic variants emerged to frequencies sufficient for forward transmission. Finally, we detected significant genetic compartmentalization of viral variants between saliva and nasal swab sample sites in many individuals. Altogether, these data provide a high-resolution profile of within-host SARS-CoV-2 evolutionary dynamics.IMPORTANCEWe detail the within-host evolutionary dynamics of SARS-CoV-2 during acute infection in 31 individuals using daily longitudinal sampling. We characterized patterns of mutational accumulation for unvaccinated and vaccinated individuals, and observed that temporal variant dynamics in both groups were largely stochastic. Comparison of paired nasal and saliva samples also revealed significant genetic compartmentalization between tissue environments in multiple individuals. Our results demonstrate how selection, genetic drift, and spatial compartmentalization all play important roles in shaping the within-host evolution of SARS-CoV-2 populations during acute infection.

Dynamic clade transitions and the influence of vaccine rollout on the spatiotemporal circulation of SARS-CoV-2 variants in São Paulo, Brazil

Since 2021, the emergence of variants of concern (VOC) has led Brazil to experience record numbers of in COVID-19 cases and deaths. The expanded spread of the SARS-CoV-2 combined with a low vaccination rate has contributed to the emergence of new mutations that may enhance viral fitness, leading to the persistence of the disease. Due to limitations in the real-time genomic monitoring of new variants in some Brazilian states, we aimed to investigate whether genomic surveillance, coupled with epidemiological data and SARS-CoV-2 variants spatiotemporal spread in a smaller region, can reflect the pandemic progression at a national level. Our findings revealed three SARS-CoV-2 variant replacements from 2021 to early 2022, corresponding to the introduction and increase in the frequency of Gamma, Delta, and Omicron variants, as indicated by peaks of the Effective Reproductive Number (Reff). These distinct clade replacements triggered two waves of COVID-19 cases, influenced by the increasing vaccine uptake over time. Our results indicated that the effectiveness of vaccination in preventing new cases during the Delta and Omicron circulations was six and eleven times higher, respectively, than during the period when Gamma was predominant, and it was highly efficient in reducing the number of deaths. Furthermore, we demonstrated that genomic monitoring at a local level can reflect the national trends in the spread and evolution of SARS-CoV-2.

An amplicon-based nanopore sequencing workflow for rapid tracking of avian influenza outbreaks, France, 2020-2022

During the recent avian influenza epizootics that occurred in France in 2020/21 and 2021/22, the virus was so contagiousness that it was impossible to control its spread between farms. The preventive slaughter of millions of birds consequently was the only solution available. In an effort to better understand the spread of avian influenza viruses (AIVs) in a rapid and innovative manner, we established an amplicon-based MinION sequencing workflow for the rapid genetic typing of circulating AIV strains. An amplicon-based MinION sequencing workflow based on a set of PCR primers targeting primarily the hemagglutinin gene but also the entire influenza virus genome was developed. Thirty field samples from H5 HPAIV outbreaks in France, including environmental samples, were sequenced using the MinION MK1C. A real-time alignment of the sequences with MinKNOW software allowed the sequencing run to be stopped as soon as enough data were generated. The consensus sequences were then generated and a phylogenetic analysis was conducted to establish links between the outbreaks. The whole sequence of the hemagglutinin gene was obtained for the 30 clinical samples of H5Nx HPAIV belonging to clade 2.3.4.4b. The consensus sequences comparison and the phylogenetic analysis demonstrated links between some outbreaks. While several studies have shown the advantages of MinION for avian influenza virus sequencing, this workflow has been applied exclusively to clinical field samples, without any amplification step on cell cultures or embryonated eggs. As this type of testing pipeline requires only a short amount of time to link outbreaks or demonstrate a new introduction, it could be applied to the real-time management of viral epizootics.

Integration of individualized and population-level molecular epidemiology data to model COVID-19 outcomes

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants with enhanced transmissibility and immune escape have emerged periodically throughout the coronavirus disease 2019 (COVID-19) pandemic, but the impact of these variants on disease severity has remained unclear. In this single-center, retrospective cohort study, we examined the association between SARS-CoV-2 clade and patient outcome over a two-year period in Chicago, Illinois. Between March 2020 and March 2022, 14,252 residual diagnostic specimens were collected from SARS-CoV-2-positive inpatients and outpatients alongside linked clinical and demographic metadata, of which 2,114 were processed for viral whole-genome sequencing. When controlling for patient demographics and vaccination status, several viral clades were associated with risk for hospitalization, but this association was negated by the inclusion of population-level confounders, including case count, sampling bias, and shifting standards of care. These data highlight the importance of integrating non-virological factors into disease severity and outcome models for the accurate assessment of patient risk.

Distinct chikungunya virus polymerase palm subdomains contribute to virus replication and virion assembly

Alphaviruses encode an error-prone RNA-dependent RNA polymerase (RdRp), nsP4, required for genome synthesis, yet how the RdRp functions in the complete alphavirus life cycle is not well-defined. Previous work using chikungunya virus (CHIKV) has established the importance of the nsP4 residue cysteine 483 in maintaining viral genetic fidelity. Given the location of residue C483 in the nsP4 palm domain, we hypothesized that other residues within this domain and surrounding subdomains would also contribute to polymerase function. To test this hypothesis, we designed a panel of nsP4 variants via homology modeling based on the Coxsackievirus B3 3 polymerase. We rescued each variant in both mammalian and mosquito cells and discovered that the palm domain and ring finger subdomain contribute to polymerase host-specific replication and genetic stability. Surprisingly, in mosquito cells, these variants in the ring finger and palm domain were replication competent and produced viral structural proteins, but they were unable to produce infectious progeny, indicating a yet uncharacterized role for the polymerase in viral assembly. Finally, we have identified additional residues in the nsP4 palm domain that influence the genetic diversity of the viral progeny, potentially via an alteration in NTP binding and/or discrimination by the polymerase. Taken together, these studies highlight that distinct nsP4 subdomains regulate multiple processes of the alphavirus life cycle, placing nsP4 in a central role during the switch from RNA synthesis to packaging and assembly.

High-performance enrichment-based genome sequencing to support the investigation of hepatitis A virus outbreaks

IMPORTANCE

This proof-of-concept study introduces a hybrid capture oligo panel for whole-genome sequencing of all six human pathogenic hepatitis A virus (HAV) subgenotypes, exhibiting a higher sensitivity than some conventional genotyping assays. The ability of hybrid capture to enrich multiple targets allows for a single, streamlined workflow, thus facilitating the potential harmonization of molecular surveillance of HAV with other enteric viruses. Even challenging sample matrices can be accommodated, making them suitable for broad implementation in clinical and public health laboratories. This innovative approach has significant implications for enhancing multijurisdictional outbreak investigations as well as our understanding of the global diversity and transmission dynamics of HAV.