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

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.

Comprehensive contact tracing during an outbreak of alpha-variant SARS-CoV-2 in a rural community reveals less viral genomic diversity and higher household secondary attack rates than expected

Sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes throughout the COVID-19 pandemic has generated a wealth of data on viral evolution across populations, but only a few studies have so far explored SARS-CoV-2 evolution across large connected transmission networks. Here, we couple data from SARS-CoV-2 sequencing with contact tracing data from an outbreak with a single origin in a rural Norwegian community where samples from all exposed persons were collected prospectively. A total of 134 nasopharyngeal samples were positive by PCR. Among the 121 retrievable genomes, 81 were identical to the genome of the introductor, thus demonstrating that genomics beyond clustering genotypically similar viral genomes to confirm relatedness offers limited additional value to manual contact tracing. In the cases where mutations were discovered, five small genetic clusters were identified. We observed a household secondary attack rate of 77%, with 92% of household members infected among households with secondary transmission, suggesting that SARS-CoV-2 introduction into large families is likely to affect all household members.

IMPORTANCE

In outbreak investigations, obtaining a full overview of infected individuals within a population is seldom achieved. We here present an example where a single introduction of B1.1.7 SARS-CoV-2 within a rural community allowed for tracing of the virus from an introductor via dissemination through larger gatherings into households. The outbreak occurred before widespread vaccination, allowing for a "natural" outbreak development with community lockdown. We show through sequencing that the virus can infect up to five consecutive persons without gaining mutations, thereby showing that contact tracing seems more important than sequencing for local outbreak investigations in settings with few alternative introductory transmission pathways. We also show how larger households where a child introduced transmission appeared more likely to promote further spread of the virus compared to households with an adult as the primary introductor.

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

SARS-CoV-2 infection induces the generation of virus-specific CD4+ and CD8+ effector and memory T cells. However, the contribution of T cells in controlling SARS-CoV-2 during infection is not well understood. Following infection of C57BL/6 mice, 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 depleting antibodies, we found that T cells within the lungs play a minimal role in viral control, and viral clearance occurs in the absence of both CD4+ and CD8+ T cells through 28 days postinfection. In the nasal compartment, depletion of both CD4+ and CD8+ T cells, but not individually, results in persistent, culturable virus replicating in the nasal epithelial layer through 28 days postinfection. Viral sequencing analysis revealed adapted mutations across the SARS-CoV-2 genome, including a large 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.

Multiple DNA Viruses and HPV Integration in Inverted Papilloma and Associated Sinonasal Carcinoma

OBJECTIVES

Sinonasal inverted papilloma (IP) has a locally destructive growth pattern, can relapse, and can undergo malignant transformation (IP-associated sinonasal squamous cell carcinoma (IP-SNSCC)). Human papillomaviruses (HPV)-6 and -16 are frequently detected in IPs. To clarify the possible roles of other DNA viruses in IPs, we explored viruses not studied in this context before. With the setting of pre- and post-malignant transformation samples, we investigated HPV genomes in depth to assess the integration of HPV into the human genome and the presence of minor intratypic variants.

MATERIALS AND METHODS

We analyzed 35 IP samples representing 28 individuals, of which six had IP-SNSCC. For virus screening, we applied qPCR to detect 16 different DNA viruses in three virus families, comprising herpesviruses, parvoviruses, and polyomaviruses. In addition, targeted next generation sequencing (NGS) was used for detailed HPV analysis.

RESULTS

We detected herpes-, parvo-, and polyomaviruses in 13/28 (46%) patients, with codetections of multiple viruses in six (21%) patients. NGS revealed HPV16 DNA in 2/6 IP-SNSCC and in their respective earlier benign IP samples, as well as in a plasma sample from one of these patients. HPV6 was detected in two IP samples without subsequent malignant transformation. We identified sequence reads containing junctions of HPV6 and HPV16 and host genome suggestive of viral integration. HPV6 and HPV16 minor intratypic variants were present across pre- and post-malignant transformation, with mostly nonsynonymous mutations.

CONCLUSIONS

Multiple DNA viruses were present in IPs. HPV16 was detected only in IP-SNSCCs or in tumors that later underwent malignant transformation.

LEVEL OF EVIDENCE

3 Laryngoscope, 2024.

High-resolution epidemiological landscape from ~290,000 SARS-CoV-2 genomes from Denmark

Vast amounts of pathogen genomic, demographic and spatial data are transforming our understanding of SARS-CoV-2 emergence and spread. We examined the drivers of molecular evolution and spread of 291,791 SARS-CoV-2 genomes from Denmark in 2021. With a sequencing rate consistently exceeding 60%, and up to 80% of PCR-positive samples between March and November, the viral genome set is broadly whole-epidemic representative. We identify a consistent rise in viral diversity over time, with notable spikes upon the importation of novel variants (e.g., Delta and Omicron). By linking genomic data with rich individual-level demographic data from national registers, we find that individuals aged  < 15 and  > 75 years had a lower contribution to molecular change (i.e., branch lengths) compared to other age groups, but similar molecular evolutionary rates, suggesting a lower likelihood of introducing novel variants. Similarly, we find greater molecular change among vaccinated individuals, suggestive of immune evasion. We also observe evidence of transmission in rural areas to follow predictable diffusion processes. Conversely, urban areas are expectedly more complex due to their high mobility, emphasising the role of population structure in driving virus spread. Our analyses highlight the added value of integrating genomic data with detailed demographic and spatial information, particularly in the absence of structured infection surveys.

The application status of sequencing technology in global respiratory infectious disease diagnosis

Next-generation sequencing (NGS) has revolutionized clinical microbiology, particularly in diagnosing respiratory infectious diseases and conducting epidemiological investigations. This narrative review summarizes conventional methods for routine respiratory infection diagnosis, including culture, smear microscopy, immunological assays, image techniques as well as polymerase chain reaction(PCR). In contrast to conventional methods, there is a new detection technology, sequencing technology, and here we mainly focus on the next-generation sequencing NGS, especially metagenomic NGS(mNGS). NGS offers significant advantages over traditional methods. Firstly, mNGS eliminates assumptions about pathogens, leading to faster and more accurate results, thus reducing diagnostic time. Secondly, it allows unbiased identification of known and novel pathogens, offering broad-spectrum coverage. Thirdly, mNGS not only identifies pathogens but also characterizes microbiomes, analyzes human host responses, and detects resistance genes and virulence factors. It can complement targeted sequencing for bacterial and fungal classification. Unlike traditional methods affected by antibiotics, mNGS is less influenced due to the extended survival of pathogen DNA in plasma, broadening its applicability. However, barriers to full integration into clinical practice persist, primarily due to cost constraints and limitations in sensitivity and turnaround time. Despite these challenges, ongoing advancements aim to improve cost-effectiveness and efficiency, making NGS a cornerstone technology for global respiratory infection diagnosis.

Tracking the prevalence and emergence of SARS-CoV-2 variants of concern using a regional genomic surveillance program

SARS-CoV-2 molecular testing coupled with whole-genome sequencing is instrumental for real-time genomic surveillance. Genomic surveillance is critical for monitoring the spread of variants of concern (VOCs) as well as discovery of novel variants. Since the beginning of the pandemic, millions of SARS-CoV-2 genomes have been deposited into public sequence databases. This is the result of efforts of both national and regional diagnostic laboratories. In this study, we describe the results of SARS-CoV-2 genomic surveillance from February 2021 to June 2022 at a metropolitan hospital in the United States. We demonstrate that consistent daily sampling is sufficient to track the regional prevalence and emergence of VOCs and recapitulate national trends. Similar sampling efforts should be considered a viable option for local SARS-CoV-2 genomic surveillance at other regional laboratories.

IMPORTANCE

In our manuscript, we describe the results of SARS-CoV-2 genomic surveillance from February 2021 to June 2022 at a metropolitan hospital in the United States. We demonstrate that consistent daily sampling is sufficient to track the regional prevalence and emergence of variants of concern (VOCs). Similar sampling efforts should be considered a viable option for local SARS-CoV-2 genomic surveillance at other regional laboratories. While the SARS-CoV-2 pandemic has evolved into a more endemic form, we still believe that additional real-world information about sampling, procedures, and data interpretation is valuable for ongoing as well as future genomic surveillance efforts. Our study should be of substantial interest to clinical virologists.

Optimization data for an ARTIC-/Illumina-based whole-genome sequencing protocol and pipeline for SARS-CoV-2 analysis

In January 2021, Germany commenced surveillance of SARS-CoV-2 variants under the Corona Surveillance Act, which ceased in July 2023. The objective was to bolster pandemic control, as specific alterations in amino acids, particularly within the spike protein, were linked to heightened transmission and decreased vaccine effectiveness. Consequently, our team conducted whole genome sequencing using the commercially accessible ARTIC protocol on Illumina's NextSeq500 platform and MiSeq for SARS-CoV-2 positive samples obtained from patients at Heidelberg University Hospital, affiliated hospitals, and the public health office in the Rhine-Neckar/Heidelberg region. Throughout the pandemic, we refined the existing ARTIC V4 protocol as well as our bioinformatics pipeline, the details of which are outlined in this report. This report reflects the protocol for the MiSeq analysis, the protocol for the NextSeq500 can be found in our previous publication.

Wastewater surveillance of severe acute respiratory syndrome coronavirus-2 in open drains of two Indian megacities captures evolutionary lineage transitions: a zonation approach

Wastewater-based environmental surveillance (WBES) has been proven as proxy tool for monitoring nucleic acids of pathogens shed by infected population before clinical outcomes. The poor sewershed network of low to middle-income countries (LMICs) leads to most of the wastewater flow through open drains. We studied the effectiveness of WBES using open drain samples to monitor the emergence of the SARS-CoV-2 variants in 2 megacities of India having dense population through zonation approach. Samples from 28 locations spanned into 5 zones of Pune region, Maharashtra, India, were collected on a weekly basis during October 2021 to July 2022. Out of 1115 total processed samples, 303 (~ 27%) tested positive for SARS-CoV-2. The periodical rise and fall in the percentage positivity of the samples was found to be in sync with the abundance of SARS-CoV-2 RNA and the reported COVID-19 active cases for Pune city. Sequencing of the RNA obtained from wastewater samples confirmed the presence of SARS-CoV-2. Of 337 sequences, lineage identification for 242 samples revealed 265 distinct SARS-CoV-2 variants including 10 highly transmissible ones. Importantly, transition from Delta to Omicron variant could be detected in wastewater samples 2 weeks prior to any clinically reported Omicron cases in India. Thus, this study demonstrates the usefulness of open drain samples for real-time monitoring of a viral pathogen's evolutionary dynamics and could be implemented in LMICs.

SARS-CoV-2-related bat viruses evade human intrinsic immunity but lack efficient transmission capacity

Circulating bat coronaviruses represent a pandemic threat. However, our understanding of bat coronavirus pathogenesis and transmission potential is limited by the lack of phenotypically characterized strains. We created molecular clones for the two closest known relatives of SARS-CoV-2, BANAL-52 and BANAL-236. We demonstrated that BANAL-CoVs and SARS-CoV-2 have similar replication kinetics in human bronchial epithelial cells. However, BANAL-CoVs have impaired replication in human nasal epithelial cells and in the upper airway of mice. We also observed reduced pathogenesis in mice and diminished transmission in hamsters. Further, we observed that diverse bat coronaviruses evade interferon and downregulate major histocompatibility complex class I. Collectively, our study demonstrates that despite high genetic similarity across bat coronaviruses, prediction of pandemic potential of a virus necessitates functional characterization. Finally, the restriction of bat coronavirus replication in the upper airway highlights that transmission potential and innate immune restriction can be uncoupled in this high-risk family of emerging viruses.

An amplicon-based sequencing approach for Usutu virus characterization

Usutu virus (USUV), an arbovirus from the Flaviviridae family, genus Flavivirus, has recently gained increasing attention because of its potential for emergence. After his discovery in South Africa, USUV spread to other African countries, then emerged in Europe where it was responsible for epizootics. The virus has recently been found in Asia. USUV infection in humans is considered to be most often asymptomatic or to cause mild clinical signs. However, a few cases of neurological complications such as encephalitis or meningo-encephalitis have been reported in both immunocompromised and immunocompetent patients. USUV natural life cycle involves Culex mosquitoes as its main vector, and multiple bird species as natural viral reservoirs or amplifying hosts, humans and horses can be incidental hosts. Phylogenetic studies carried out showed eight lineages, showing an increasing genetic diversity for USUV. This work describes the development and validation of a novel whole-genome amplicon-based sequencing approach to Usutu virus. This study was carried out on different strains from Senegal and Italy. The new approach showed good coverage using samples derived from several vertebrate hosts and may be valuable for Usutu virus genomic surveillance to better understand the dynamics of evolution and transmission of the virus.

Genomic Epidemiology of Rift Valley Fever Virus Involved in the 2018 and 2022 Outbreaks in Livestock in Rwanda

Rift Valley fever (RVF), a mosquito-borne transboundary zoonosis, was first confirmed in Rwanda's livestock in 2012 and since then sporadic cases have been reported almost every year. In 2018, the country experienced its first large outbreak, which was followed by a second one in 2022. To determine the circulating virus lineages and their ancestral origin, two genome sequences from the 2018 outbreak, and thirty-six, forty-one, and thirty-eight sequences of small (S), medium (M), and large (L) genome segments, respectively, from the 2022 outbreak were generated. All of the samples from the 2022 outbreak were collected from slaughterhouses. Both maximum likelihood and Bayesian-based phylogenetic analyses were performed. The findings showed that RVF viruses belonging to a single lineage, C, were circulating during the two outbreaks, and shared a recent common ancestor with RVF viruses isolated in Uganda between 2016 and 2019, and were also linked to the 2006/2007 largest East Africa RVF outbreak reported in Kenya, Tanzania, and Somalia. Alongside the wild-type viruses, genetic evidence of the RVFV Clone 13 vaccine strain was found in slaughterhouse animals, demonstrating a possible occupational risk of exposure with unknown outcome for people working in meat-related industry. These results provide additional evidence of the ongoing wide spread of RVFV lineage C in Africa and emphasize the need for an effective national and international One Health-based collaborative approach in responding to RVF emergencies.

Wastewater surveillance together with metaviromic data revealed the unusual resurgence of infectious diseases after the first wave of the COVID-19 outbreak

How to address public health priorities after COVID-19 is becoming a critical task. To this end, we conducted wastewater surveillance for six leading pathogens, namely, SARS-CoV-2, norovirus, rotavirus, influenza A virus (IAV), enteroviruses and respiratory syncytial virus (RSV), in Nanchang city from January to April 2023. Metaviromic sequencing was conducted at the 1st, 4th, 7th, 9th, 12th and 14th weeks to reveal the dynamics of viral pathogens that were not covered by qPCR. Amplicon sequencing of the conserved region of norovirus GI and GII and the rotavirus and region encoding nonstructural protein of RSV was also conducted weekly. The results showed that after a rapid decrease in SARS-CoV-2 sewage concentrations occurred in January 2023, surges of norovirus, rotavirus, IAV and RSV started at the 6th, 7th, 8th and 11th weeks, respectively. The dynamics of the sewage concentrations of norovirus, rotavirus, IAV and RSV were consistent with the off-season resurgence of the above infectious diseases. Notably, peak sewage concentrations of norovirus GI, GII, rotavirus, IAV and RSV were found at the 6th, 3rd, 7th, 7th and 8th weeks, respectively. Astroviruses also resurge after the 7th week, as revealed by metaviromic data, suggesting that wastewater surveillance together with metaviromic data provides an essential early warning tool for revealing patterns of infectious disease resurgence.

Optimising nucleic acid recovery from rapid antigen tests for whole genome sequencing of respiratory viruses

BACKGROUND

Whole genome sequencing (WGS) of respiratory viruses from rapid antigen tests (RAT-WGS) is a novel approach to expanding genomic surveillance of respiratory infections. To date however, there are limited data on the genomic stability of these viruses on RATs. In this study, we investigated the effect of storage conditions and nucleic acid preservatives on the ability to enhance stability and improve recovery of respiratory virus genomes from RATs.

METHODS

A mixture of common respiratory viruses was used to inoculate RATs at different environmental temperatures (4°C, 20°C and 36°C), with two preservative reagents (RNALater and DNA/RNA shield) Nucleic acid was extracted from RATs at two different timepoints (72 h and seven days) and subject to real-time multiplex respiratory PCR to detect a range of respiratory viruses. WGS was performed using target-enrichment with the TWIST Comprehensive Viral Research Panel. Defined metrics from an automated in-house bioinformatic pipeline were used to assess and compare viral genome recovery under different conditions.

RESULTS

Nucleic acid degradation (indicated by relative change in PCR cycle threshold and WGS-based metrics) was most notable at 20 °C and 36 °C. Storage in either RNALater or DNA / RNA shield improved genome recovery for respiratory viruses across all temperature conditions, although this was most pronounced for RNALater. Subtyping of Influenza viruses demonstrated the applicability of RAT-WGS in downstream genomic epidemiological surveillance.

CONCLUSIONS

Under simulated conditions, RAT-WGS demonstrated that (i) viral genomes were generally stable at 4°C at 72 h and 1 week, (ii) RNALater has a more significant preservation of nucleic acids compared to DNA/RNA Shield and (iii) genome recovery can be achieved using a sequencing depth of 500,000 reads per sample in RNALater, across all respiratory viruses and conditions.

Dispersal history of SARS-CoV-2 in Galicia, Spain

The dynamics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission are influenced by a variety of factors, including social restrictions and the emergence of distinct variants. In this study, we delve into the origins and dissemination of the Alpha, Delta, and Omicron-BA.1 variants of concern in Galicia, northwest Spain. For this, we leveraged genomic data collected by the EPICOVIGAL Consortium and from the GISAID database, along with mobility information from other Spanish regions and foreign countries. Our analysis indicates that initial introductions during the Alpha phase were predominantly from other Spanish regions and France. However, as the pandemic progressed, introductions from Portugal and the United States became increasingly significant. The number of detected introductions varied from 96 and 101 for Alpha and Delta to 39 for Omicron-BA.1. Most of these introductions left a low number of descendants (<10), suggesting a limited impact on the evolution of the pandemic in Galicia. Notably, Galicia's major coastal cities emerged as critical hubs for viral transmission, highlighting their role in sustaining and spreading the virus. This research emphasizes the critical role of regional connectivity in the spread of SARS-CoV-2 and offers essential insights for enhancing public health strategies and surveillance measures.

Entourage: all-in-one sequence analysis software for genome assembly, virus detection, virus discovery, and intrasample variation profiling

BACKGROUND

Pan-virus detection, and virome investigation in general, can be challenging, mainly due to the lack of universally conserved genetic elements in viruses. Metagenomic next-generation sequencing can offer a promising solution to this problem by providing an unbiased overview of the microbial community, enabling detection of any viruses without prior target selection. However, a major challenge in utilising metagenomic next-generation sequencing for virome investigation is that data analysis can be highly complex, involving numerous data processing steps.

RESULTS

Here, we present Entourage to address this challenge. Entourage enables short-read sequence assembly, viral sequence search with or without reference virus targets using contig-based approaches, and intrasample sequence variation quantification. Several workflows are implemented in Entourage to facilitate end-to-end virus sequence detection analysis through a single command line, from read cleaning, sequence assembly, to virus sequence searching. The results generated are comprehensive, allowing for thorough quality control, reliability assessment, and interpretation. We illustrate Entourage's utility as a streamlined workflow for virus detection by employing it to comprehensively search for target virus sequences and beyond in raw sequence read data generated from HeLa cell culture samples spiked with viruses. Furthermore, we showcase its flexibility and performance on a real-world dataset by analysing a preassembled Tara Oceans dataset. Overall, our results show that Entourage performs well even with low virus sequencing depth in single digits, and it can be used to discover novel viruses effectively. Additionally, by using sequence data generated from a patient with chronic SARS-CoV-2 infection, we demonstrate Entourage's capability to quantify virus intrasample genetic variations, and generate publication-quality figures illustrating the results.

CONCLUSIONS

Entourage is an all-in-one, versatile, and streamlined bioinformatics software for virome investigation, developed with a focus on ease of use. Entourage is available at https://codeberg.org/CENMIG/Entourage under the MIT license.

Direct whole-genome sequencing of HIV-1 for clinical drug-resistance analysis and public health surveillance

BACKGROUND

Human Immunodeficiency virus type 1 (HIV-1) remains a significant global health threat partly due to its ability to develop resistance to anti-retroviral therapies. HIV-1 genotype and drug resistance analysis of the polymerase (pol) sequence is a mainstay of its clinical and public health management. However, as new treatments and resistances evolve, analysis methods must change accordingly. In this study, we outline the development and implementation of a direct whole-genome sequencing approach (dWGS) using probe-capture target-enrichment for HIV-1 genotype and drug resistance analysis.

METHODS

We implemented dWGS and performed parallel pol Sanger sequencing for clinical samples, followed by comparative genotype and drug-resistance analysis. These HIV-1 WGS sequences were also utilised for a novel partitioned phylogenetic analysis.

RESULTS

Optimised nucleic acid extraction and DNAse I treatment significantly increased HIV-1 whole-genome coverage and depth, and improved recovery of high-quality genomes from low viral load clinical samples, enabling routine sequencing of viral loads as low as 1000 copies/mL. Overall, dWGS was robust, accurate and more sensitive for detecting low-frequency variants at drug-resistance sites compared to Sanger sequencing. Analysis of multiple sequence regions improved phylogenetic reconstruction for recombinant HIV-1 sequences compared to analysis of pol sequence alone.

CONCLUSIONS

These findings demonstrate dWGS enhances HIV-1 drug-resistance analysis by quantitative variant detection and improves reconstruction of HIV-1 phylogenies compared to traditional pol sequencing. This work supports that HIV-1 dWGS is a viable option to replace Sanger sequencing for clinical and public health applications.

Epidemiological Characterization and Genetic Variation of the SARS-CoV-2 Delta Variant in Palestine

The emergence of new SARS-CoV-2 variants in Palestine highlights the need for continuous genetic surveillance and accurate screening strategies. This case series study aimed to investigate the geographic distribution and genetic variation of the SARS-CoV-2 Delta Variant in Palestine in August 2021. Samples were collected at random in August 2021 (n = 571) from eight districts in the West Bank, Palestine. All samples were confirmed as positive for COVID-19 by RT-PCR. The samples passed the quality control test and were successfully sequenced using the ARTIC protocol. The Delta Variant was revealed to have four dominant lineages: B.1.617 (19%), AY.122 (18%), AY.106 (17%), and AY.121 (13%). The study revealed eight significant purely spatial clusters (p < 0.005) distributed in the northern and southern parts of Palestine. Phylogenetic analysis of SARS-CoV-2 genomes (n = 552) showed no geographically specific clades. The haplotype network revealed three haplogroups without any geographic distribution. Chronologically, the Delta Variant peak in Palestine was shortly preceded by the one in the neighboring Israeli community and shortly followed by the peak in Jordan. In addition, the study revealed an extremely intense transmission network of the Delta Variant circulating between the Palestinian districts as hubs (SHR ≈ 0.5), with Al-Khalil, the district with the highest prevalence of COVID-19, witnessing the highest frequency of transitions. Genetic diversity analysis indicated closely related haplogroups, as haplotype diversity (Hd) is high but has low nucleotide diversity (π). However, nucleotide diversity (π) in Palestine is still higher than the global figures. Neutrality tests were significantly (p < 0.05) low, including Tajima's D, Fu-Li's F, and Fu-Li's D, suggesting one or more of the following: population expansion, selective sweep, and natural negative selection. Wright's F-statistic (Fst) showed genetic differentiation (Fst > 0.25) with low to medium gene flow (Nm). Recombination events were minimal between clusters (Rm) and between adjacent sites (Rs). The study confirms the utility of the whole genome sequence as a surveillance system to track the emergence of new SARS-CoV-2 variants for any possible geographical association and the use of genetic variation analysis and haplotype networking to delineate any minimal change or slight deviation in the viral genome from a reference strain.

T cell-mediated Immune response and correlates of inflammation and their relationship with COVID-19 clinical severity: not an intuitive guess

BACKGROUND

Predictors of the outcome of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection remain to be fully determined. We evaluated selected viral characteristics and immunological responses that might predict and/or correlate to the clinical outcome of COVID-19.

METHODS

For individuals developing divergent clinical outcomes, the magnitude and breadth of T cell-mediated responses were measured within 36 h of symptom onset. Peripheral Blood Mononuclear Cells (PBMCs) were subjected to in vitro stimulation with SARS-CoV-2-based peptides. In addition, SARS-CoV-2 sequences were generated by metagenome, and HLA typing was performed using Luminex technology.

FINDINGS

CD4+ T cell activation was negatively correlated with SARS-CoV-2 basal viral load in patients with severe COVID-19 (p = 0·043). The overall cellular immune response, as inferred by the IFN-γ signal, was higher at baseline for patients who progressed to mild disease compared to patients who progressed to severe disease (p = 0·0044). Subjects with milder disease developed higher T cell responses for MHC class I and II-restricted peptides (p = 0·033).

INTERPRETATION

Mounting specific cellular immune responses in the first days after symptom onset, as inferred by IFN-γ magnitude in the ELISPOT assay, may efficiently favor a positive outcome. In contrast, progression to severe COVID-19 was accompanied by stronger cellular immune responses, higher CD4 + T cell activation, and a higher number of in silico predicted high-affinity class I HLA alleles.

Molecular epidemiology of recurrent zoonotic transmission of mpox virus in West Africa

Nigeria and Cameroon reported their first mpox cases in over three decades in 2017 and 2018 respectively. The outbreak in Nigeria is recognised as an ongoing human epidemic. However, owing to sparse surveillance and genomic data, it is not known whether the increase in cases in Cameroon is driven by zoonotic or sustained human transmission. Notably, the frequency of zoonotic transmission remains unknown in both Cameroon and Nigeria. To address these uncertainties, we investigated the zoonotic transmission dynamics of the mpox virus (MPXV) in Cameroon and Nigeria, with a particular focus on the border regions. We show that in these regions mpox cases are still driven by zoonotic transmission of a newly identified Clade IIb.1. We identify two distinct zoonotic lineages that circulate across the Nigeria-Cameroon border, with evidence of recent and historic cross border dissemination. Our findings support that the complex cross-border forest ecosystems likely hosts shared animal populations that drive cross-border viral spread, which is likely where extant Clade IIb originated. We identify that the closest zoonotic outgroup to the human epidemic circulated in southern Nigeria in October 2013. We also show that the zoonotic precursor lineage circulated in an animal population in southern Nigeria for more than 45 years. This supports findings that southern Nigeria was the origin of the human epidemic. Our study highlights the ongoing MPXV zoonotic transmission in Cameroon and Nigeria, underscoring the continuous risk of MPXV (re)emergence.

Genomic epidemiology uncovers the timing and origin of the emergence of mpox in humans

Five years before the 2022-2023 global mpox outbreak Nigeria reported its first cases in nearly 40 years, with the ongoing epidemic since driven by sustained human-to-human transmission. However, limited genomic data has left questions about the timing and origin of the mpox virus' (MPXV) emergence. Here we generated 112 MPXV genomes from Nigeria from 2021-2023. We identify the closest zoonotic outgroup to the human epidemic in southern Nigeria, and estimate that the lineage transmitting from human-to-human emerged around July 2014, circulating cryptically until detected in September 2017. The epidemic originated in Southern Nigeria, particularly Rivers State, which also acted as a persistent and dominant source of viral dissemination to other states. We show that APOBEC3 activity increased MPXV's evolutionary rate twenty-fold during human-to-human transmission. We also show how Delphy, a tool for near-real-time Bayesian phylogenetics, can aid rapid outbreak analytics. Our study sheds light on MPXV's establishment in West Africa before the 2022-2023 global outbreak and highlights the need for improved pathogen surveillance and response.

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.

Emergence of transmissible SARS-CoV-2 variants with decreased sensitivity to antivirals in immunocompromised patients with persistent infections

We investigated the impact of antiviral treatment on the emergence of SARS-CoV-2 resistance during persistent infections in immunocompromised patients (n=15). All patients received remdesivir and some also received nirmatrelvir-ritonavir or monoclonal antibodies. Sequence analysis showed that nine patients carried viruses with mutations in the nsp12 (RNA dependent RNA polymerase), while four had viruses with nsp5 (3C protease) mutations. Infectious SARS-CoV-2 with a double mutation in nsp5 (T169I) and nsp12 (V792I) was recovered from respiratory secretions 77 days after initial COVID-19 diagnosis from a patient treated with remdesivir and nirmatrelvir-ritonavir. In vitro characterization confirmed its decreased sensitivity to remdesivir and nirmatrelvir, which was overcome by combined antiviral treatment. Studies in golden Syrian hamsters demonstrated efficient transmission to contact animals. This study documents the isolation of SARS-CoV-2 carrying resistance mutations to both nirmatrelvir and remdesivir from a patient and demonstrates its transmissibility in vivo.

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 were detected 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 for wildlife as well as commercial poultry operations have detected these viruses. In this study, 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 whole-genome sequencing. Sequencing 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.IMPORTANCEWhile surveillance programs for avian influenza viruses are often focused on migratory routes and their associated stop-over locations or commercial poultry operations, many bird species-including migratory birds-frequent or live in urban green spaces and wetlands. This brings them into contact with a highly dense population of humans and pets, providing an extensive urban animal-human interface in which the general public may have little awareness of circulating infectious diseases. This study focuses on virus surveillance of this interface, combined with culturally responsive science education and community outreach.

ViralFlow v1.0-a computational workflow for streamlining viral genomic surveillance

ViralFlow v1.0 is a computational workflow developed for viral genomic surveillance. Several key changes turned ViralFlow into a general-purpose reference-based genome assembler for all viruses with an available reference genome. New virus-agnostic modules were implemented to further study nucleotide and amino acid mutations. ViralFlow v1.0 runs on a broad range of computational infrastructures, from laptop computers to high-performance computing (HPC) environments, and generates standard and well-formatted outputs suited for both public health reporting and scientific problem-solving. ViralFlow v1.0 is available at: https://viralflow.github.io/index-en.html.

Mapping the Evolutionary Space of SARS-CoV-2 Variants to Anticipate Emergence of Subvariants Resistant to COVID-19 Therapeutics

New sublineages of SARS-CoV-2 variants-of-concern (VOCs) continuously emerge with mutations in the spike glycoprotein. In most cases, the sublineage-defining mutations vary between the VOCs. It is unclear whether these differences reflect lineage-specific likelihoods for mutations at each spike position or the stochastic nature of their appearance. Here we show that SARS-CoV-2 lineages have distinct evolutionary spaces (a probabilistic definition of the sequence states that can be occupied by expanding virus subpopulations). This space can be accurately inferred from the patterns of amino acid variability at the whole-protein level. Robust networks of co-variable sites identify the highest-likelihood mutations in new VOC sublineages and predict remarkably well the emergence of subvariants with resistance mutations to COVID-19 therapeutics. Our studies reveal the contribution of low frequency variant patterns at heterologous sites across the protein to accurate prediction of the changes at each position of interest.

Sequencing and characterization of human bocavirus genomes from patients diagnosed in Southern France between 2017 and 2022

The diversity and evolution of the genomes of human bocavirus (HBoV), which causes respiratory diseases, have been scarcely studied. Here, we aimed to obtain and characterize HBoV genomes from patients's nasopharyngeal samples collected between 2017 and 2022 period (5 years and 7 months). Next-generation sequencing (NGS) used Illumina technology after having implemented using GEMI an in-house multiplex PCR amplification strategy. Genomes were assembled and analyzed with CLC Genomics, Mafft, BioEdit, MeV, Nextclade, MEGA, and iTol. A total of 213 genomes were obtained. Phylogeny classified them all as of Bocavirus 1 (HBoV1) species. Five HBoV1 genotypic clusters determined by hierarchical clustering analysis of 27 variable genome positions were scattered over the study period although with differences in yearly prevalence. A total of 167 amino acid substitutions were detected. Besides, coinfection was observed for 52% of the samples, rhinoviruses then adenoviruses (HAdVs) being the most common viruses. Principal component analysis showed that HBoV1 genotypic cluster α tended to be correlated with HAdV co-infection. Subsequent HAdV typing for HBoV1-positive samples and negative controls demonstrated that HAdVC species predominated but HAdVB was that significantly HBoV1-associated. Overall, we described here the first HBoV1 genomes sequenced for France. HBoV1 and HAdVB association deserves further investigation.

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.

Dispersion of SARS-CoV-2 lineage BA.5.1.25 and its descendants in Peru during two COVID-19 waves in 2022

During the third year of the pandemic in Peru, the persistent transmission of SARS-CoV-2 led to the appearance of more transmissible and immune-evasive Omicron sublineages; in that context, the National Genomic Surveillance of SARS-CoV-2 performed by the Peruvian National Institute of Health detected spike mutations in the circulating Omicron BA.5.1.25 sublineage which was later designated as DJ.1 and increased during the fourth COVID-19 wave, this eventually branched into new sublineages. The introduction, emergence, and timing of the most recent common ancestor (tMRCA) of BA.5.1.25 and its descendants (DJ.1, DJ.1.1, DJ.1.2, and DJ.1.3) were investigated in this paper as well as the time lags between their emergence and identification by the Peruvian National Institute of Health. Our findings show that ongoing genomic surveillance of SARS-CoV-2 is critical for understanding its phylogenetic evolution and the emergence of novel variations.

Development and application of a whole genome amplicon sequencing method for infectious salmon anemia virus (ISAV)

Infectious salmon anemia (ISA) is an infectious disease primarily affecting farmed Atlantic salmon, Salmo salar, which is caused by the ISA virus (ISAV). ISAV belongs to the Orthomyxoviridae family. The disease is a serious condition resulting in reduced fish welfare and high mortality. In this study, we designed an amplicon-based sequencing protocol for whole genome sequencing of ISAV. The method consists of 80 ISAV-specific primers that cover 92% of the virus genome and was designed to be used on an Illumina MiSeq platform. The sequencing accuracy was investigated by comparing sequences with previously published Sanger sequences. The sequences obtained were nearly identical to those obtained by Sanger sequencing, thus demonstrating that sequences produced by this amplicon sequencing protocol had an acceptable accuracy. The amplicon-based sequencing method was used to obtain the whole genome sequence of 12 different ISAV isolates from a small local epidemic in the northern part of Norway. Analysis of the whole genome sequences revealed that segment reassortment took place between some of the isolates and could identify which segments that had been reassorted.

IAVCP (Influenza A Virus Consensus and Phylogeny): Automatic Identification of the Genomic Sequence of the Influenza A Virus from High-Throughput Sequencing Data

Recently, high-throughput sequencing of influenza A viruses has become a routine test. It should be noted that the extremely high diversity of the influenza A virus complicates the task of determining the sequences of all eight genome segments. For a fast and accurate analysis, it is necessary to select the most suitable reference for each segment. At the same time, there is no standardized method in the field of decoding sequencing results that allows the user to update the sequence databases to which the reads obtained by virus sequencing are compared. The IAVCP (influenza A virus consensus and phylogeny) was developed with the goal of automatically analyzing high-throughput sequencing data of influenza A viruses. Its goals include the extraction of a consensus genome directly from paired raw reads. In addition, the pipeline enables the identification of potential reassortment events in the evolutionary history of the virus of interest by analyzing the topological structure of phylogenetic trees that are automatically reconstructed.

Crykey: Rapid identification of SARS-CoV-2 cryptic mutations in wastewater

Wastewater surveillance for SARS-CoV-2 provides early warnings of emerging variants of concerns and can be used to screen for novel cryptic linked-read mutations, which are co-occurring single nucleotide mutations that are rare, or entirely missing, in existing SARS-CoV-2 databases. While previous approaches have focused on specific regions of the SARS-CoV-2 genome, there is a need for computational tools capable of efficiently tracking cryptic mutations across the entire genome and investigating their potential origin. We present Crykey, a tool for rapidly identifying rare linked-read mutations across the genome of SARS-CoV-2. We evaluated the utility of Crykey on over 3,000 wastewater and over 22,000 clinical samples; our findings are three-fold: i) we identify hundreds of cryptic mutations that cover the entire SARS-CoV-2 genome, ii) we track the presence of these cryptic mutations across multiple wastewater treatment plants and over three years of sampling in Houston, and iii) we find a handful of cryptic mutations in wastewater mirror cryptic mutations in clinical samples and investigate their potential to represent real cryptic lineages. In summary, Crykey enables large-scale detection of cryptic mutations in wastewater that represent potential circulating cryptic lineages, serving as a new computational tool for wastewater surveillance of SARS-CoV-2.

Enhancing the epidemiological surveillance of SARS-CoV-2 using Sanger sequencing to identify circulating variants and recombinants

Since the emergence of SARS-CoV-2 in December 2019, more than 12,000 mutations in the virus have been identified. These could cause changes in viral characteristics and directly impact global public health. The emergence of variants is a great concern due to the chance of increased transmissibility and infectivity. Sequencing for surveillance and monitoring circulating strains is extremely necessary as the early identification of new variants allows public health agencies to make faster and more effective decisions to contain the spread of the virus. In the present study, we identified circulating variants in samples collected in Belo Horizonte, Brazil, and detected a recombinant lineage using the Sanger method. The identification of lineages was done through gene amplification of SARS-CoV-2 by Reverse Transcription-Polymerase Chain Reaction (RT-PCR). By using these specific fragments, we were able to differentiate one variant of interest and five circulating variants of concern. We were also able to detect recombinants. Randomly selected samples were sequenced by either Sanger or Next Generation Sequencing (NGS). Our findings validate the effectiveness of Sanger sequencing as a powerful tool for monitoring variants. It is easy to perform and allows the analysis of a larger number of samples in countries that cannot afford NGS.

Efficient SARS-CoV-2 variant detection and monitoring with Spike Screen next-generation sequencing

The emergence and rapid spread of SARS-CoV-2 prompted the global community to identify innovative approaches to diagnose infection and sequence the viral genome because at several points in the pandemic positive case numbers exceeded the laboratory capacity to characterize sufficient samples to adequately respond to the spread of emerging variants. From week 10, 2020, to week 13, 2023, Slovenian routine complete genome sequencing (CGS) surveillance network yielded 41 537 complete genomes and revealed a typical molecular epidemiology with early lineages gradually being replaced by Alpha, Delta, and finally Omicron. We developed a targeted next-generation sequencing based variant surveillance strategy dubbed Spike Screen through sample pooling and selective SARS-CoV-2 spike gene amplification in conjunction with CGS of individual cases to increase throughput and cost-effectiveness. Spike Screen identifies variant of concern (VOC) and variant of interest (VOI) signature mutations, analyses their frequencies in sample pools, and calculates the number of VOCs/VOIs at the population level. The strategy was successfully applied for detection of specific VOC/VOI mutations prior to their confirmation by CGS. Spike Screen complemented CGS efforts with an additional 22 897 samples sequenced in two time periods: between week 42, 2020, and week 24, 2021, and between week 37, 2021, and week 2, 2022. The results showed that Spike Screen can be applied to monitor VOC/VOI mutations among large volumes of samples in settings with limited sequencing capacity through reliable and rapid detection of novel variants at the population level and can serve as a basis for public health policy planning.

Outbreak of influenza and SARS-CoV-2 at the Armed Forces of the Philippines Health Service Education and Training Center, September 25-October 10, 2023

In the last week of September 2023, a surge of influenza-like illness was observed among students of the Armed Forces of the Philippines (AFP) Health Service Education and Training Center, where 48 (27 males and 21 females; age in years: mean 33, range 27-41) of 247 military students at the Center presented with respiratory symptoms. Between September 25 and October 10, 2023, all 48 symptomatic students were evaluated with real-time reverse transcription polymerase chain reaction and sequencing for both influenza and SARS-CoV-2. Thirteen (27%) students were found positive for influenza A/H3 only, 6 (13%) for SARS-CoV-2 only, and 4 (8%) were co-infected with influenza A/H3 and SARS-CoV-2. Seventeen influenza A/ H3N2 viruses belonged to the same clade, 3C.2a1b.2a.2a.3a, and 4 SARSCoV-2 sequences belonged to the JE1.1 lineage, indicating a common source outbreak for both. The influenza A/H3N2 circulating virus belonged to a different clade than the vaccine strain for 2023 (3C.2a1b.2a.2a). Only 4 students had received the influenza vaccine for 2023. In response, the AFP Surgeon General issued a memorandum to all military health institutions on October 19, 2023 that mandated influenza vaccination as a prerequisite for enrollment of students at all education and training centers, along with implementation of non-pharmaceutical interventions and early notification and testing of students exhibiting influenza-like-illness.

Genomic Context of SARS-CoV-2 Outbreaks in Farmed Mink in Spain during Pandemic: Unveiling Host Adaptation Mechanisms

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects various mammalian species, with farmed minks experiencing the highest number of outbreaks. In Spain, we analyzed 67 whole genome sequences and eight spike sequences from 18 outbreaks, identifying four distinct lineages: B.1, B.1.177, B.1.1.7, and AY.98.1. The potential risk of transmission to humans raises crucial questions about mutation accumulation and its impact on viral fitness. Sequencing revealed numerous not-lineage-defining mutations, suggesting a cumulative mutation process during the outbreaks. We observed that the outbreaks were predominantly associated with different groups of mutations rather than specific lineages. This clustering pattern by the outbreaks could be attributed to the rapid accumulation of mutations, particularly in the ORF1a polyprotein and in the spike protein. Notably, the mutations G37E in NSP9, a potential host marker, and S486L in NSP13 were detected. Spike protein mutations may enhance SARS-CoV-2 adaptability by influencing trimer stability and binding to mink receptors. These findings provide valuable insights into mink coronavirus genetics, highlighting both host markers and viral transmission dynamics within communities.

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.

Sequencing analysis of SARS-CoV-2 cases in Slovenian long-term care facilities to support outbreak control

Introduction

Residents of long-term care facilities (LTCFs) are at high risk of morbidity and mortality due to COVID-19, especially when new variants of concern (VOC) emerge. To provide intradisciplinary data in order to tailor public health interventions during future epidemics, available epidemiologic and genomic data from Slovenian LTCFs during the initial phases of the COVID-19 pandemic was analyzed.

Methods

The first part of the study included SARS-CoV-2 reverse-transcription Real-Time PCR (rtRT-PCR) positive LTCF residents, from 21 facilities with COVID-19 outbreaks occurring in October 2020. The second part of the study included SARS-CoV-2 rtRT-PCR positive LTCF residents and staff between January and April 2021, when VOC Alpha emerged in Slovenia. Next-generation sequencing (NGS) was used to acquire SARS-CoV-2 genomes, and lineage determination. In-depth phylogenetic and mutational profile analysis were performed and coupled with available field epidemiological data to assess the dynamics of SARS-CoV-2 introduction and transmission.

Results

370/498 SARS-CoV-2 positive residents as well as 558/699 SARS-CoV-2 positive residents and 301/358 staff were successfully sequenced in the first and second part of the study, respectively. In October 2020, COVID-19 outbreaks in the 21 LTCFs were caused by intra-facility transmission as well as multiple independent SARS-CoV-2 introductions. The Alpha variant was confirmed in the first LTCF resident approximately 1.5 months after the first Alpha case was identified in Slovenia. The data also showed a slower replacement of existing variants by Alpha in residents compared to staff and the general population.

Discussion

Multiple SARS CoV-2 introductions as well as intra-facility spreading impacted disease transmission in Slovenian LTCFs. Timely implementation of control measures aimed at limiting new introductions while controlling in-facility transmission are of paramount importance, especially as new VOCs emerge. Sequencing, in conjunction with epidemiological data, can facilitate the determination of the need for future improvements in control measures to protect LTCF residents from COVID-19 or other respiratory infections.

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.

Within-host influenza viral diversity in the pediatric population as a function of age, vaccine, and health status

Seasonal influenza virus predominantly evolves through antigenic drift, marked by the accumulation of mutations at antigenic sites. Because of antigenic drift, influenza vaccines are frequently updated, though their efficacy may still be limited due to strain mismatches. Despite the high levels of viral diversity observed across populations, most human studies reveal limited intrahost diversity, leaving the origin of population-level viral diversity unclear. Previous studies show host characteristics, such as immunity, might affect within-host viral evolution. Here we investigate influenza A viral diversity in children aged between 6 months and 18 years. Influenza virus evolution in children is less well characterized than in adults, yet may be associated with higher levels of viral diversity given the lower level of pre-existing immunity and longer durations of infection in children. We obtained influenza isolates from banked influenza A-positive nasopharyngeal swabs collected at the Children's Hospital of Philadelphia during the 2017-18 influenza season. Using next-generation sequencing, we evaluated the population of influenza viruses present in each sample. We characterized within-host viral diversity using the number and frequency of intrahost single-nucleotide variants (iSNVs) detected in each sample. We related viral diversity to clinical metadata, including subjects' age, vaccination status, and comorbid conditions, as well as sample metadata such as virus strain and cycle threshold. Consistent with previous studies, most samples contained low levels of diversity with no clear association between the subjects' age, vaccine status, or health status. Further, there was no enrichment of iSNVs near known antigenic sites. Taken together, these findings are consistent with previous observations that the majority of intrahost influenza virus infection is characterized by low viral diversity without evidence of diversifying selection.

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.