Rapid SARS-CoV-2 whole-genome sequencing and analysis for informed public health decision-making in the Netherlands

Prospective and systematic screening for invasive aspergillosis in the ICU during the COVID-19 pandemic, a proof of principle for future pandemics

The diagnostic performance of a prospective, systematic screening strategy for COVID-19 associated pulmonary aspergillosis (CAPA) during the COVID-19 pandemic was investigated. Patients with COVID-19 admitted to the ICU were screened for CAPA twice weekly by collection of tracheal aspirate (TA) for Aspergillus culture and PCR. Subsequently, bronchoalveolar lavage (BAL) sampling was performed in patients with positive screening results and clinical suspicion of infection. Patient data were collected from April 2020-February 2022. Patients were classified according to 2020 ECMM/ISHAM consensus criteria. In total, 126/370 (34%) patients were positive in screening and CAPA frequency was 52/370 (14%) (including 13 patients negative in screening). CAPA was confirmed in 32/43 (74%) screening positive patients who underwent BAL sampling. ICU mortality was 62% in patients with positive screening and confirmed CAPA, and 31% in CAPA cases who were screening negative. The sensitivity, specificity, positive and negative predictive value (PPV & NPV) of screening for CAPA were 0.71, 0.73, 0.27, and 0.95, respectively. The PPV was higher if screening was culture positive compared to PCR positive only, 0.42 and 0.12 respectively. CAPA was confirmed in 74% of screening positive patients, and culture of TA had a better diagnostic performance than PCR. Positive screening along with clinical manifestations appeared to be a good indication for BAL sampling since diagnosis of CAPA was confirmed in most of these patients. Prospective, systematic screening allowed to quickly gain insight into the epidemiology of fungal superinfections during the pandemic and could be applicable for future pandemics.

Clinical and Virological Outcome of Monoclonal Antibody Therapies Across Severe Acute Respiratory Syndrome Coronavirus 2 Variants in 245 Immunocompromised Patients: A Multicenter Prospective Cohort Study

BACKGROUND

Immunocompromised patients (ICPs) have an increased risk for a severe and prolonged COVID-19. SARS-CoV-2 monoclonal antibodies (mAbs) were extensively used in these patients, but data from randomized trials that focus on ICPs are lacking. We evaluated the clinical and virological outcome of COVID-19 in ICPs treated with mAbs across SARS-CoV-2 variants.

METHODS

In this multicenter prospective cohort study, we enrolled B-cell- and/or T-cell-deficient patients treated with casirivimab/imdevimab, sotrovimab, or tixagevimab/cilgavimab. SARS-CoV-2 RNA was quantified and sequenced weekly, and time to viral clearance, viral genome mutations, hospitalization, and death rates were registered.

RESULTS

Two hundred and forty five patients infected with the Delta (50%) or Omicron BA.1, 2, or 5 (50%) variant were enrolled. Sixty-seven percent were vaccinated; 78 treated as outpatients, of whom 2 required hospital admission, but both survived. Of the 159 patients hospitalized at time of treatment, 43 (27%) required mechanical ventilation or died. The median time to viral clearance was 14 days (interquartile range, 7-22); however, it took >30 days in 15%. Resistance-associated spike mutations emerged in 9 patients in whom the median time to viral clearance was 63 days (95% confidence interval, 57-69; P < .001). Spike mutations were observed in 1 of 42 (2.4%) patients after treatment with 2 active mAbs, in 5 of 34 (14.7%) treated with actual monotherapy (sotrovimab), and 3 of 20 (12%) treated with functional monotherapy (ie, tixagevimab/cilgavimab against tixagevimab-resistant variant).

CONCLUSIONS

Despite treatment with mAbs, morbidity and mortality of COVID-19 in ICPs remained substantial. Combination antiviral therapy should be further explored and may be preferred in severely ICPs.

Whole genome sequencing and phylogenetic analysis of SARS-CoV-2 strains isolated during the COVID-19 pandemic in Nigeria

Objectives

The emergence and spread of SARS-CoV-2 have stimulated ongoing research into the virus transmission dynamics, circulating variants, and potential mutations. This study was conducted to understand the genomic dynamics of the epidemic in Nigeria.

Design

Whole genome sequencing was conducted on SARS-CoV-2 samples collected during the first and second outbreaks using the Oxford Nanopore MinION sequencing platform. Phylogenetic analysis was conducted, and genomes were grouped into different pangolin lineages.

Results

The study revealed four circulating SARS-CoV-2 variants. The Alpha (B.1.1.7) variant was the most prevalent (32.7%), followed by Beta (B.1 B.1.1, L.3, and B.1.1.318) (30.8%), Eta (B.1.525) (28.9%), and Delta (B.1.617, AY.1, AY.109, and AY.36) (7.7%). Phylogenetic analysis revealed three clusters with four Nextstrain clades (20I, 20B, 21D, and 21J). The Alpha lineages (B.1.1.7) clustered with references from Italy. The Beta lineages (Clade 20B) (B.11, B.11318, and L3) and sub-lineage B.11 were distinct. Sub-lineage B.11318 is clustered with references from the USA, whereas sub-lineage L3 is clustered with references from Russia, the Philippines, Australia, and Japan. The 21D and 21J, belonging to two Pango lineages, Eta (B.1525) and Delta (B.1.617 and AY.109), showed high genetic similarity.

Conclusion

The phylogenetic relatedness of the lineages suggests multiple virus introduction, which could be a source of more virulent, locally adapted variants.

Virus sequencing performance during the SARS-CoV-2 pandemic: a retrospective analysis of data from multiple rounds of external quality assessment in Austria

Introduction: A notable feature of the 2019 coronavirus disease (COVID-19) pandemic was the widespread use of whole genome sequencing (WGS) to monitor severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections. Countries around the world relied on sequencing and other forms of variant detection to perform contact tracing and monitor changes in the virus genome, in the hopes that epidemic waves caused by variants would be detected and managed earlier. As sequencing was encouraged and rewarded by the government in Austria, but represented a new technicque for many laboratories, we designed an external quality assessment (EQA) scheme to monitor the accuracy of WGS and assist laboratories in validating their methods. Methods: We implemented SARS-CoV-2 WGS EQAs in Austria and report the results from 7 participants over 5 rounds from February 2021 until June 2023. The participants received sample material, sequenced genomes with routine methods, and provided the sequences as well as information about mutations and lineages. Participants were evaluated on the completeness and accuracy of the submitted sequence and the ability to analyze and interpret sequencing data. Results: The results indicate that performance was excellent with few exceptions, and these exceptions showed improvement over time. We extend our findings to infer that most publicly available sequences are accurate within ≤1 nucleotide, somewhat randomly distributed through the genome. Conclusion: WGS continues to be used for SARS-CoV-2 surveillance, and will likely be instrumental in future outbreak scenarios. We identified hurdles in building next-generation sequencing capacity in diagnostic laboratories. EQAs will help individual laboratories maintain high quality next-generation sequencing output, and strengthen variant monitoring and molecular epidemiology efforts.

Deciphering complex antibiotic resistance patterns in Helicobacter pylori through whole genome sequencing and machine learning

Introduction

Helicobacter pylori (H.pylori, Hp) affects billions of people worldwide. However, the emerging resistance of Hp to antibiotics challenges the effectiveness of current treatments. Investigating the genotype-phenotype connection for Hp using next-generation sequencing could enhance our understanding of this resistance.

Methods

In this study, we analyzed 52 Hp strains collected from various hospitals. The susceptibility of these strains to five antibiotics was assessed using the agar dilution assay. Whole-genome sequencing was then performed to screen the antimicrobial resistance (AMR) genotypes of these Hp strains. To model the relationship between drug resistance and genotype, we employed univariate statistical tests, unsupervised machine learning, and supervised machine learning techniques, including the development of support vector machine models.

Results

Our models for predicting Amoxicillin resistance demonstrated 66% sensitivity and 100% specificity, while those for Clarithromycin resistance showed 100% sensitivity and 100% specificity. These results outperformed the known resistance sites for Amoxicillin (A1834G) and Clarithromycin (A2147), which had sensitivities of 22.2% and 87%, and specificities of 100% and 96%, respectively.

Discussion

Our study demonstrates that predictive modeling using supervised learning algorithms with feature selection can yield diagnostic models with higher predictive power compared to models relying on single single-nucleotide polymorphism (SNP) sites. This approach significantly contributes to enhancing the precision and effectiveness of antibiotic treatment strategies for Hp infections. The application of whole-genome sequencing for Hp presents a promising pathway for advancing personalized medicine in this context.

Innovation and Patenting Activities During COVID-19 and Advancement of Biochemical and Molecular Diagnosis in the Post- COVID-19 Era

The COVID-19 pandemic is to escalate globally and acquire new mutations quickly, so accurate diagnostic technologies play a vital role in controlling and understanding the epidemiology of the disease. A plethora of technologies acquires diagnosis of individuals and informs clinical management of COVID. Some important biochemical parameters for COVID diagnosis are the elevation of liver enzymes, creatinine, and nonspecific inflammatory markers such as C-reactive protein (CRP) and Interleukin 6 (IL-6). The main progression predictors are lymphopenia, elevated D-dimer, and hyperferritinemia, although it is also necessary to consider LDH, CPK, and troponin in the marker panel of diagnosis. Owing to the greater sensitivity and accuracy, molecular technologies such as conventional polymerase chain reaction (PCR), reverse transcription (RT)-PCR, nested PCR, loop-mediated isothermal amplification (LAMP), and xMAP technology have been extensively used for COVID diagnosis for some time now. To make so many diagnostics accessible to general people, many techniques may be exploited, including point of care (POC), also called bedside testing, which is developing as a portable promising tool in pathogen identification. Some other lateral flow assay (LFA)-centered techniques like SHERLOCK, CRISPR-Cas12a (AIOD-CRISPR), and FNCAS9 editor limited uniform detection assay (FELUDA), etc. have shown auspicious results in the rapid detection of pathogens. More recently, low-cost sequencing and advancements in big data management have resulted in a slow but steady rise of next-generation sequencing (NGS)-based approaches for diagnosis that have potential relevance for clinical purposes and may pave the way toward a better future. Due to the COVID-19 pandemic, various institutions provided free, specialized websites and tools to promote research and access to critically needed advanced solutions by alleviating research and analysis of data within a substantial body of scientific and patent literature regarding biochemical and molecular diagnosis published since January 2020. This circumstance is unquestionably unique and difficult for anyone using patent information to find pertinent disclosures at a specific date in a trustworthy manner.

Advancing pathogen genomics in resource-limited settings

Genomic sequencing has emerged as a powerful tool to enhance early pathogen detection and characterization with implications for public health and clinical decision making. Although widely available in developed countries, the application of pathogen genomics among low-resource, high-disease burden settings remains at an early stage. In these contexts, tailored approaches for integrating pathogen genomics within infectious disease control programs will be essential to optimize cost efficiency and public health impact. We propose a framework for embedding pathogen genomics within national surveillance plans across a spectrum of surveillance and laboratory capacities. We adopt a public health approach to genomics and examine its application to high-priority diseases relevant in resource-limited settings. For each grouping, we assess the value proposition for genomics to inform public health and clinical decision-making, alongside its contribution toward research and development of novel diagnostics, therapeutics, and vaccines.

Advancing animal tuberculosis surveillance using culture-independent long-read whole-genome sequencing

Animal tuberculosis is a significant infectious disease affecting both livestock and wildlife populations worldwide. Effective disease surveillance and characterization of Mycobacterium bovis (M. bovis) strains are essential for understanding transmission dynamics and implementing control measures. Currently, sequencing of genomic information has relied on culture-based methods, which are time-consuming, resource-demanding, and concerning in terms of biosafety. This study explores the use of culture-independent long-read whole-genome sequencing (WGS) for a better understanding of M. bovis epidemiology in African buffaloes (Syncerus caffer). By comparing two sequencing approaches, we evaluated the efficacy of Illumina WGS performed on culture extracts and culture-independent Oxford Nanopore adaptive sampling (NAS). Our objective was to assess the potential of NAS to detect genomic variants without sample culture. In addition, culture-independent amplicon sequencing, targeting mycobacterial-specific housekeeping and full-length 16S rRNA genes, was applied to investigate the presence of microorganisms, including nontuberculous mycobacteria. The sequencing quality obtained from DNA extracted directly from tissues using NAS is comparable to the sequencing quality of reads generated from culture-derived DNA using both NAS and Illumina technologies. We present a new approach that provides complete and accurate genome sequence reconstruction, culture independently, and using an economically affordable technique.

Public Health Response Measures for COVID-19 in Kazakhstan

The article describes measures developed to counter the spread of coronavirus infection in the Republic of Kazakhstan. The first cases of coronavirus disease 2019 (COVID-19) in Kazakhstan were detected on March 13, 2020, among people who arrived from Germany. After declaring the state of emergency in the country, the Ministry of Healthcare of the Republic of Kazakhstan began to formulate and implement a comprehensive package of measures aimed at slowing down and stopping the transmission of infection, preventing outbreaks, ensuring optimal care for all patients, especially the seriously ill, minimizing the negative impact of the pandemic on health systems, social services, and economic activities. Developed set of restrictive measures was approved by the Country Office of Word Health Organization (WHO) in Kazakhstan, being later adapted by the European Union (EU) countries and applied in Kyrgyzstan. In addition, article identifies Kazakhstan's experience in creating epidemiological surveillance systems, studying virus mutations, and the clinical aspects of dealing with it to combat the infection. It also indicates the impact of the epidemic on health-care workers and the development of measures to protect them, strengthening infection prevention, and control in medical organizations.

Whole-genome sequencing to investigate transmission of SARS-CoV-2 in the acute healthcare setting: a systematic review

BACKGROUND

Whole-genome sequencing (WGS) has been used widely to elucidate transmission of SARS-CoV-2 in acute healthcare settings, and to guide infection, prevention, and control (IPC) responses.

AIM

To systematically appraise available literature, published between January 1st, 2020 and June 30th, 2022, describing the implementation of WGS in acute healthcare settings to characterize nosocomial SARS-CoV-2 transmission.

METHODS

Searches of the PubMed, Embase, Ovid MEDLINE, EBSCO MEDLINE, and Cochrane Library databases identified studies in English reporting the use of WGS to investigate SARS-CoV-2 transmission in acute healthcare environments. Publications involved data collected up to December 31st, 2021, and findings were reported in accordance with the Preferred Reporting Items for Systematic Reviews and Meta-Analyses statement.

FINDINGS

In all, 3088 non-duplicate records were retrieved; 97 met inclusion criteria, involving 62 outbreak analyses and 35 genomic surveillance studies. No publications from low-income countries were identified. In 87/97 (90%), WGS supported hypotheses for nosocomial transmission, while in 46 out of 97 (47%) suspected transmission events were excluded. An IPC intervention was attributed to the use of WGS in 18 out of 97 (18%); however, only three (3%) studies reported turnaround times ≤7 days facilitating near real-time IPC action, and none reported an impact on the incidence of nosocomial COVID-19 attributable to WGS.

CONCLUSION

WGS can elucidate transmission of SARS-CoV-2 in acute healthcare settings to enhance epidemiological investigations. However, evidence was not identified to support sequencing as an intervention to reduce the incidence of SARS-CoV-2 in hospital or to alter the trajectory of active outbreaks.

Genomic epidemiology of SARS-CoV-2 from Uttar Pradesh, India

The various strains and mutations of SARS-CoV-2 have been tracked using several forms of genomic classification systems. The present study reports high-throughput sequencing and analysis of 99 SARS-CoV-2 specimens from Western Uttar Pradesh using sequences obtained from the GISAID database, followed by phylogeny and clade classification. Phylogenetic analysis revealed that Omicron lineages BA-2-like (55.55%) followed by Delta lineage-B.1.617.2 (45.5%) were predominantly circulating in this area Signature substitution at positions S: N501Y, S: D614G, S: T478K, S: K417N, S: E484A, S: P681H, and S: S477N were commonly detected in the Omicron variant-BA-2-like, however S: D614G, S: L452R, S: P681R and S: D950N were confined to Delta variant-B.1.617.2. We have also identified three escape variants in the S gene at codon position 19 (T19I/R), 484 (E484A/Q), and 681 (P681R/H) during the fourth and fifth waves in India. Based on the phylogenetic diversification studies and similar changes in other lineages, our analysis revealed indications of convergent evolution as the virus adjusts to the shifting immunological profile of its human host. To the best of our knowledge, this study is an approach to comprehensively map the circulating SARS-CoV-2 strains from Western Uttar Pradesh using an integrated approach of whole genome sequencing and phylogenetic analysis. These findings will be extremely valuable in developing a structured approach toward pandemic preparedness and evidence-based intervention plans in the future.

Integrating full and partial genome sequences to decipher the global spread of canine rabies virus

Despite the rapid growth in viral genome sequencing, statistical methods face challenges in handling historical viral endemic diseases with large amounts of underutilized partial sequence data. We propose a phylogenetic pipeline that harnesses both full and partial viral genome sequences to investigate historical pathogen spread between countries. Its application to rabies virus (RABV) yields precise dating and confident estimates of its geographic dispersal. By using full genomes and partial sequences, we reduce both geographic and genetic biases that often hinder studies that focus on specific genes. Our pipeline reveals an emergence of the present canine-mediated RABV between years 1301 and 1403 and reveals regional introductions over a 700-year period. This geographic reconstruction enables us to locate episodes of human-mediated introductions of RABV and examine the role that European colonization played in its spread. Our approach enables phylogeographic analysis of large and genetically diverse data sets for many viral pathogens.

Evolution of Sequence and Structure of SARS-CoV-2 Spike Protein: A Dynamic Perspective

Novel coronavirus (SARS-CoV-2) enters its host cell through a surface spike protein. The viral spike protein has undergone several modifications/mutations at the genomic level, through which it modulated its structure-function and passed through several variants of concern. Recent advances in high-resolution structure determination and multiscale imaging techniques, cost-effective next-generation sequencing, and development of new computational methods (including information theory, statistical methods, machine learning, and many other artificial intelligence-based techniques) have hugely contributed to the characterization of sequence, structure, function of spike proteins, and its different variants to understand viral pathogenesis, evolutions, and transmission. Laying on the foundation of the sequence-structure-function paradigm, this review summarizes not only the important findings on structure/function but also the structural dynamics of different spike components, highlighting the effects of mutations on them. As dynamic fluctuations of three-dimensional spike structure often provide important clues for functional modulation, quantifying time-dependent fluctuations of mutational events over spike structure and its genetic/amino acidic sequence helps identify alarming functional transitions having implications for enhanced fusogenicity and pathogenicity of the virus. Although these dynamic events are more difficult to capture than quantifying a static, average property, this review encompasses those challenging aspects of characterizing the evolutionary dynamics of spike sequence and structure and their implications for functions.

SARS-CoV-2 in lions, gorillas and zookeepers in the Rotterdam Zoo, the Netherlands, a One Health investigation, November 2021

In November 2021, seven western lowland gorillas and four Asiatic lions were diagnosed with COVID-19 at Rotterdam Zoo. An outbreak investigation was undertaken to determine the source and extent of the outbreak and to identify possible transmission routes. Interviews were conducted with staff to identify human and animal contacts and cases, compliance with personal protective equipment (PPE) and potential transmission routes. Human and animal contacts and other animal species suspected to be susceptible to SARS-CoV-2 were tested for SARS-CoV-2 RNA. Positive samples were subjected to sequencing. All the gorillas and lions that could be tested (3/7 and 2/4, respectively) were RT-PCR positive between 12 November and 10 December 2021. No other animal species were SARS-CoV-2 RNA positive. Forty direct and indirect human contacts were identified. Two direct contacts tested RT-PCR positive 10 days after the first COVID-19 symptoms in animals. The zookeepers' viral genome sequences clustered with those of gorillas and lions. Personal protective equipment compliance was suboptimal at instances. Findings confirm transmission of SARS-CoV-2 among animals and between humans and animals but source and directionality could not be established. Zookeepers were the most likely source and should have periodic PPE training. Sick animals should promptly be tested and isolated/quarantined.

Tracking SARS-CoV-2 introductions in Mozambique using pandemic-scale phylogenies: a retrospective observational study

BACKGROUND

From the start of the SARS-CoV-2 outbreak, global sequencing efforts have generated an unprecedented amount of genomic data. Nonetheless, unequal sampling between high-income and low-income countries hinders the implementation of genomic surveillance systems at the global and local level. Filling the knowledge gaps of genomic information and understanding pandemic dynamics in low-income countries is essential for public health decision making and to prepare for future pandemics. In this context, we aimed to discover the timing and origin of SARS-CoV-2 variant introductions in Mozambique, taking advantage of pandemic-scale phylogenies.

METHODS

We did a retrospective, observational study in southern Mozambique. Patients from Manhiça presenting with respiratory symptoms were recruited, and those enrolled in clinical trials were excluded. Data were included from three sources: (1) a prospective hospital-based surveillance study (MozCOVID), recruiting patients living in Manhiça, attending the Manhiça district hospital, and fulfilling the criteria of suspected COVID-19 case according to WHO; (2) symptomatic and asymptomatic individuals with SARS-CoV-2 infection recruited by the National Surveillance system; and (3) sequences from SARS-CoV-2-infected Mozambican cases deposited on the Global Initiative on Sharing Avian Influenza Data database. Positive samples amenable for sequencing were analysed. We used Ultrafast Sample placement on Existing tRees to understand the dynamics of beta and delta waves, using available genomic data. This tool can reconstruct a phylogeny with millions of sequences by efficient sample placement in a tree. We reconstructed a phylogeny (~7·6 million sequences) adding new and publicly available beta and delta sequences.

FINDINGS

A total of 5793 patients were recruited between Nov 1, 2020, and Aug 31, 2021. During this time, 133 328 COVID-19 cases were reported in Mozambique. 280 good quality new SARS-CoV-2 sequences were obtained after the inclusion criteria were applied and an additional 652 beta (B.1.351) and delta (B.1.617.2) public sequences were included from Mozambique. We evaluated 373 beta and 559 delta sequences. We identified 187 beta introductions (including 295 sequences), divided in 42 transmission groups and 145 unique introductions, mostly from South Africa, between August, 2020 and July, 2021. For delta, we identified 220 introductions (including 494 sequences), with 49 transmission groups and 171 unique introductions, mostly from the UK, India, and South Africa, between April and November, 2021.

INTERPRETATION

The timing and origin of introductions suggests that movement restrictions effectively avoided introductions from non-African countries, but not from surrounding countries. Our results raise questions about the imbalance between the consequences of restrictions and health benefits. This new understanding of pandemic dynamics in Mozambique can be used to inform public health interventions to control the spread of new variants.

FUNDING

European and Developing Countries Clinical Trials, European Research Council, Bill & Melinda Gates Foundation, and Agència de Gestió d'Ajuts Universitaris i de Recerca.

Oxford nanopore technologies-a valuable tool to generate whole-genome sequencing data for in silico serotyping and the detection of genetic markers in Salmonella

Bacteria of the genus Salmonella pose a major risk to livestock, the food economy, and public health. Salmonella infections are one of the leading causes of food poisoning. The identification of serovars of Salmonella achieved by their diverse surface antigens is essential to gain information on their epidemiological context. Traditionally, slide agglutination has been used for serotyping. In recent years, whole-genome sequencing (WGS) followed by in silico serotyping has been established as an alternative method for serotyping and the detection of genetic markers for Salmonella. Until now, WGS data generated with Illumina sequencing are used to validate in silico serotyping methods. Oxford Nanopore Technologies (ONT) opens the possibility to sequence ultra-long reads and has frequently been used for bacterial sequencing. In this study, ONT sequencing data of 28 Salmonella strains of different serovars with epidemiological relevance in humans, food, and animals were taken to investigate the performance of the in silico serotyping tools SISTR and SeqSero2 compared to traditional slide agglutination tests. Moreover, the detection of genetic markers for resistance against antimicrobial agents, virulence, and plasmids was studied by comparing WGS data based on ONT with WGS data based on Illumina. Based on the ONT data from flow cell version R9.4.1, in silico serotyping achieved an accuracy of 96.4 and 92% for the tools SISTR and SeqSero2, respectively. Highly similar sets of genetic markers comparing both sequencing technologies were identified. Taking the ongoing improvement of basecalling and flow cells into account, ONT data can be used for Salmonella in silico serotyping and genetic marker detection.

SARS-CoV-2 transmission dynamics in bars, restaurants, and nightclubs

Background

In an attempt to control the spread of SARS-CoV-2, many governments decided to close public venues including bars, restaurants, and nightclubs during the pandemic, making it difficult to study how transmission occurs in these environments. In this study, we were able to gain insight into the transmission dynamics of SARS-CoV-2 in 16 venues in the city of Maastricht using a combination of epidemiological and whole-genome sequencing (WGS) data during a period of 2 weeks in 2021, when bars, restaurants, and nightclubs were temporarily reopened in the Netherlands. This led to a subsequent rise of SARS-CoV-2 cases in the community following the reopening.

Methods

WGS was performed on samples from 154/348 of selected cases and combined with epidemiological investigation (e.g., contact tracing and linking cases to specific venues) to identify SARS-CoV-2 transmission clusters. In addition, genomic surveillance data were used to investigate spillover of outbreak-associated genotypes into the community.

Results

Clustering was observed in 129/136 (95%) successfully genotyped samples. We established that most cases were linked to venues with dancing facilities and that specific genotypes of the Delta variant were more frequently spread within and from these venues compared to venues without dancing facilities. In addition, we show indications of spillover of certain genotypes from the bar and restaurant industry into the community, with the number of hospital admissions increasing in the weeks following peak cases in the community.

Conclusion

Lifting restrictions on bar and restaurant industry venues with a corona entree ticket in a largely unvaccinated population led to a surge in COVID-19 cases and promoted the spread of new (sub)variants. Nightclubs were identified as potential super-spreading locations.

Rise and fall of SARS-CoV-2 variants in Rotterdam: Comparison of wastewater and clinical surveillance

Monitoring of SARS-CoV-2 in wastewater (WW) is a promising tool for epidemiological surveillance, correlating not only viral RNA levels with the infection dynamics within the population, but also to viral diversity. However, the complex mixture of viral lineages in WW samples makes tracking of specific variants or lineages circulating in the population a challenging task. We sequenced sewage samples of 9 WW-catchment areas within the city of Rotterdam, used specific signature mutations from individual SARS-CoV-2 lineages to estimate their relative abundances in WW and compared them against those observed in clinical genomic surveillance of infected individuals between September 2020 and December 2021. We showed that especially for dominant lineages, the median of the frequencies of signature mutations coincides with the occurrence of those lineages in Rotterdam's clinical genomic surveillance. This, along with digital droplet RT-PCR targeting signature mutations of specific variants of concern (VOCs), showed that several VOCs emerged, became dominant and were replaced by the next VOC in Rotterdam at different time points during the study. In addition, single nucleotide variant (SNV) analysis provided evidence that spatio-temporal clusters can also be discerned from WW samples. We were able to detect specific SNVs in sewage, including one resulting in the Q183H amino acid change in the Spike gene, that was not captured by clinical genomic surveillance. Our results highlight the potential use of WW samples for genomic surveillance, increasing the set of epidemiological tools to monitor SARS-CoV-2 diversity.

Combining epidemiological data and whole genome sequencing to understand SARS-CoV-2 transmission dynamics in a large tertiary care hospital during the first COVID-19 wave in The Netherlands focusing on healthcare workers

BACKGROUND

Healthcare facilities have been challenged by the risk of SARS-CoV-2 transmission between healthcare workers (HCW) and patients. During the first wave of the COVID-19 pandemic, infections among HCW were observed, questioning infection prevention and control (IPC) measures implemented at that time.

AIM

This study aimed to identify nosocomial transmission routes of SARS-CoV-2 between HCW and patients in a tertiary care hospital.

METHODS

All SARS-CoV-2 PCR positive HCW and patients identified between 1 March and 19 May 2020, were included in the analysis. Epidemiological data were collected from patient files and HCW contact tracing interviews. Whole genome sequences of SARS-CoV-2 were generated using Nanopore sequencing (WGS). Epidemiological clusters were identified, whereafter WGS and epidemiological data were combined for re-evaluation of epidemiological clusters and identification of potential transmission clusters. HCW infections were further classified into categories based on the likelihood that the infection was acquired via nosocomial transmission. Secondary cases were defined as COVID-19 cases in our hospital, part of a transmission cluster, of which the index case was either a patient or HCW from our hospital.

FINDINGS

The study population consisted of 293 HCW and 245 patients. Epidemiological data revealed 36 potential epidemiological clusters, with an estimated 222 (75.7%) HCW as secondary cases. WGS results were available for 195 HCW (88.2%) and 20 patients (12.8%) who belonged to an epidemiological cluster. Re-evaluation of the epidemiological clusters, with the available WGS data identified 31 transmission clusters with 65 (29.4%) HCW as secondary cases. Transmission clusters were all part of 18 (50.0%) previously determined epidemiological clusters, demonstrating that several larger outbreaks actually consisted, of several smaller transmission clusters. A total of 21 (7.2%) HCW infections were classified as from confirmed nosocomial, of which 18 were acquired from another HCW and 3 from a patient.

CONCLUSION

The majority of SARS-CoV-2 infections among HCW could be attributed to community-acquired infection. Infections among HCW that could be classified as due to nosocomial transmission, were mainly caused by HCW-to-HCW transmission rather than patient-to-HCW transmission. It is important to recognize the uncertainties of cluster analyses based solely on epidemiological data.

Robust Performance of SARS-CoV-2 Whole-Genome Sequencing from Wastewater with a Nonselective Virus Concentration Method

The sequencing of human virus genomes from wastewater samples is an efficient method for tracking viral transmission and evolution at the community level. However, this requires the recovery of viral nucleic acids of high quality. We developed a reusable tangential-flow filtration system to concentrate and purify viruses from wastewater for genome sequencing. A pilot study was conducted with 94 wastewater samples from four local sewersheds, from which viral nucleic acids were extracted, and the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was sequenced using the ARTIC V4.0 primers. Our method yielded a high probability (0.9) of recovering complete or near-complete SARS-CoV-2 genomes (>90% coverage at 10× depth) from wastewater when the COVID-19 incidence rate exceeded 33 cases per 100 000 people. The relative abundances of sequenced SARS-CoV-2 variants followed the trends observed from patient-derived samples. We also identified SARS-CoV-2 lineages in wastewater that were underrepresented or not present in the clinical whole-genome sequencing data. The developed tangential-flow filtration system can be easily adopted for the sequencing of other viruses in wastewater, particularly those at low concentrations.

Optimised protocol for monitoring SARS-CoV-2 in wastewater using reverse complement PCR-based whole-genome sequencing

Monitoring the spread of viral pathogens in the population during epidemics is crucial for mounting an effective public health response. Understanding the viral lineages that constitute the infections in a population can uncover the origins and transmission patterns of outbreaks and detect the emergence of novel variants that may impact the course of an epidemic. Population-level surveillance of viruses through genomic sequencing of wastewater captures unbiased lineage data, including cryptic asymptomatic and undiagnosed infections, and has been shown to detect infection outbreaks and novel variant emergence before detection in clinical samples. Here, we present an optimised protocol for quantification and sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in influent wastewater, used for high-throughput genomic surveillance in England during the COVID-19 pandemic. This protocol utilises reverse compliment PCR for library preparation, enabling tiled amplification across the whole viral genome and sequencing adapter addition in a single step to enhance efficiency. Sequencing of synthetic SARS-CoV-2 RNA provided evidence validating the efficacy of this protocol, while data from high-throughput sequencing of wastewater samples demonstrated the sensitivity of this method. We also provided guidance on the quality control steps required during library preparation and data analysis. Overall, this represents an effective method for high-throughput sequencing of SARS-CoV-2 in wastewater which can be applied to other viruses and pathogens of humans and animals.

A comparison of five Illumina, Ion Torrent, and nanopore sequencing technology-based approaches for whole genome sequencing of SARS-CoV-2

Rapid identification of the rise and spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern remains critical for monitoring of the efficacy of diagnostics, therapeutics, vaccines, and control strategies. A wide range of SARS-CoV-2 next-generation sequencing (NGS) methods have been developed over the last years, but cross-sequence technology benchmarking studies have been scarce. In the current study, 26 clinical samples were sequenced using five protocols: AmpliSeq SARS-CoV-2 (Illumina), EasySeq RC-PCR SARS-CoV-2 (Illumina/NimaGen), Ion AmpliSeq SARS-CoV-2 (Thermo Fisher), custom primer sets (Oxford Nanopore Technologies (ONT)), and capture probe-based viral metagenomics (Roche/Illumina). Studied parameters included genome coverage, depth of coverage, amplicon distribution, and variant calling. The median SARS-CoV-2 genome coverage of samples with cycle threshold (Ct) values of 30 and lower ranged from 81.6 to 99.8% for, respectively, the ONT protocol and Illumina AmpliSeq protocol. Correlation of coverage with PCR Ct values varied per protocol. Amplicon distribution signatures differed across the methods, with peak differences of up to 4 log₁₀ at disbalanced positions in samples with high viral loads (Ct values ≤ 23). Phylogenetic analyses of consensus sequences showed clustering independent of the workflow used. The proportion of SARS-CoV-2 reads in relation to background sequences, as a (cost-)efficiency metric, was the highest for the EasySeq protocol. The hands-on time was the lowest when using EasySeq and ONT protocols, with the latter additionally having the shortest sequence runtime. In conclusion, the studied protocols differed on a variety of the studied metrics. This study provides data that assist laboratories when selecting protocols for their specific setting.

Looking back on the COVID-19 pandemic in an elite sports team using whole genome sequencing

Objectives

The aim of this study was to investigate the effectiveness of infection control measures to prevent transmission of SARS-CoV-2 within a professional sports team using whole genome sequencing (WGS).

Design

Prospective cohort study.

Methods

74 players and staff members of a Dutch professional male football team were followed from August 2020 until May 2021. A set of health and safety measures were introduced and all participants underwent regular SARS-CoV-2 RNA testing. All positive samples were subsequently sequenced (Nanopore sequencing) to assess whether infections were acquired within the training center or in the community.

Results

Throughout the study period, 13 participants tested positive for SARS-CoV-2. The phylogenetic analysis revealed 3 clusters (of 2 and 3 cases respectively), indicating that 3/13 cases (23%) acquired infection from another player or staff member. The first cluster was diagnosed upon enrolment, thus transmission had occurred prior to the implementation of health and safety protocols. Finally, 4 cases were diagnosed prior to symptom onset, emphasizing that frequent testing leads to early detection and isolation.

Conclusion

These data show that a combination of regular testing and basic control measures can prevent outbreaks of COVID-19 in a professional sports team. WGS is an important tool to distinguish between infections introduced from the community and infections transmitted between athletes.

Recent Advances in the Lipid Nanoparticle-Mediated Delivery of mRNA Vaccines

Lipid nanoparticles (LNPs) have recently emerged as one of the most advanced technologies for the highly efficient in vivo delivery of exogenous mRNA, particularly for COVID-19 vaccine delivery. LNPs comprise four different lipids: ionizable lipids, helper or neutral lipids, cholesterol, and lipids attached to polyethylene glycol (PEG). In this review, we present recent the advances and insights for the design of LNPs, as well as their composition and properties, with a subsequent discussion on the development of COVID-19 vaccines. In particular, as ionizable lipids are the most critical drivers for complexing the mRNA and in vivo delivery, the role of ionizable lipids in mRNA vaccines is discussed in detail. Furthermore, the use of LNPs as effective delivery vehicles for vaccination, genome editing, and protein replacement therapy is explained. Finally, expert opinion on LNPs for mRNA vaccines is discussed, which may address future challenges in developing mRNA vaccines using highly efficient LNPs based on a novel set of ionizable lipids. Developing highly efficient mRNA delivery systems for vaccines with improved safety against some severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants remains difficult.

A high-resolution melt curve toolkit to identify lineage-defining SARS-CoV-2 mutations

The emergence of severe acute respiratory syndrome 2 (SARS-CoV-2) variants of concern (VOCs), with mutations linked to increased transmissibility, vaccine escape and virulence, has necessitated the widespread genomic surveillance of SARS-CoV-2. This has placed a strain on global sequencing capacity, especially in areas lacking the resources for large scale sequencing activities. Here we have developed three separate multiplex high-resolution melting assays to enable the identification of Alpha, Beta, Delta and Omicron VOCs. The assays were evaluated against whole genome sequencing on upper-respiratory swab samples collected during the Alpha, Delta and Omicron [BA.1] waves of the UK pandemic. The sensitivities of the eight individual primer sets were all 100%, and specificity ranged from 94.6 to 100%. The multiplex HRM assays have potential as a tool for high throughput surveillance of SARS-CoV-2 VOCs, particularly in areas with limited genomics facilities.

Accurate Detection of SARS-CoV-2 by Next-Generation Sequencing in Low Viral Load Specimens

As new SARS-CoV-2 variants emerge, there is an urgent need to increase the efficiency and availability of viral genome sequencing, notably to detect the lineage in samples with a low viral load. SARS-CoV-2 genome next-generation sequencing (NGS) was performed retrospectively in a single center on 175 positive samples from individuals. An automated workflow used the Ion AmpliSeq SARS-CoV-2 Insight Research Assay on the Genexus Sequencer. All samples were collected in the metropolitan area of the city of Nice (France) over a period of 32 weeks (from 19 July 2021 to 11 February 2022). In total, 76% of cases were identified with a low viral load (Ct ≥ 32, and ≤200 copies/µL). The NGS analysis was successful in 91% of cases, among which 57% of cases harbored the Delta variant, and 34% the Omicron BA.1.1 variant. Only 9% of cases had unreadable sequences. There was no significant difference in the viral load in patients infected with the Omicron variant compared to the Delta variant (Ct values, p = 0.0507; copy number, p = 0.252). We show that the NGS analysis of the SARS-CoV-2 genome provides reliable detection of the Delta and Omicron SARS-CoV-2 variants in low viral load samples.

Magnetic hydrogel particles improve nanopore sequencing of SARS-CoV-2 and other respiratory viruses

Presented here is a magnetic hydrogel particle enabled workflow for capturing and concentrating SARS-CoV-2 from diagnostic remnant swab samples that significantly improves sequencing results using the Oxford Nanopore Technologies MinION sequencing platform. Our approach utilizes a novel affinity-based magnetic hydrogel particle, circumventing low input sample volumes and allowing for both rapid manual and automated high throughput workflows that are compatible with Nanopore sequencing. This approach enhances standard RNA extraction protocols, providing up to 40 × improvements in viral mapped reads, and improves sequencing coverage by 20-80% from lower titer diagnostic remnant samples. Furthermore, we demonstrate that this approach works for contrived influenza virus and respiratory syncytial virus samples, suggesting that it can be used to identify and improve sequencing results of multiple viruses in VTM samples. These methods can be performed manually or on a KingFisher automation platform.

Evaluation of Nanopore sequencing for Mycobacterium tuberculosis drug susceptibility testing and outbreak investigation: a genomic analysis

BACKGROUND

Mycobacterium tuberculosis whole-genome sequencing (WGS) has been widely used for genotypic drug susceptibility testing (DST) and outbreak investigation. For both applications, Illumina technology is used by most public health laboratories; however, Nanopore technology developed by Oxford Nanopore Technologies has not been thoroughly evaluated. The aim of this study was to determine whether Nanopore sequencing data can provide equivalent information to Illumina for transmission clustering and genotypic DST for M tuberculosis.

METHODS

In this genomic analysis, we analysed 151 M tuberculosis isolates from Madagascar, South Africa, and England, which were collected between 2011 and 2018, using phenotypic DST and matched Illumina and Nanopore data. Illumina sequencing was done with the MiSeq, HiSeq 2500, or NextSeq500 platforms and Nanopore sequencing was done on the MinION or GridION platforms. Using highly reliable PacBio sequencing assemblies and pairwise distance correlation between Nanopore and Illumina data, we optimise Nanopore variant filters for detecting single-nucleotide polymorphisms (SNPs; using BCFtools software). We then used those SNPs to compare transmission clusters identified by Nanopore with the currently used UK Health Security Agency Illumina pipeline (COMPASS). We compared Illumina and Nanopore WGS-based DST predictions using the Mykrobe software and mutation catalogue.

FINDINGS

The Nanopore BCFtools pipeline identified SNPs with a median precision of 99·3% (IQR 99·1-99·6) and recall of 90·2% (88·1-94·2) compared with a precision of 99·6% (99·4-99·7) and recall of 91·9% (87·6-98·6) using the Illumina COMPASS pipeline. Using a threshold of 12 SNPs for putative transmission clusters, Illumina identified 98 isolates as unrelated and 53 as belonging to 19 distinct clusters (size range 2-7). Nanopore reproduced 15 out of 19 clusters perfectly; two clusters were merged into one cluster, one cluster had a single sample missing, and one cluster had an additional sample adjoined. Illumina-based clusters were also closely replicated using a five SNP threshold and clustering accuracy was maintained using mixed Illumina and Nanopore datasets. Genotyping resistance variants with Nanopore was highly concordant with Illumina, having zero discordant SNPs across more than 3000 SNPs and four insertions or deletions (indels), across 60 000 indels.

INTERPRETATION

Illumina and Nanopore technologies can be used independently or together by public health laboratories performing M tuberculosis genotypic DST and outbreak investigations. As a result, clinical and public health institutions making decisions on which sequencing technology to adopt for tuberculosis can base the choice on cost (which varies by country), batching, and turnaround time.

FUNDING

Academy for Medical Sciences, Oxford Wellcome Institutional Strategic Support Fund, and the Swiss South Africa Joint Research Award (Swiss National Science Foundation and South African National Research Foundation).

[Spreading Dynamics of the Alpha and Beta Virus Variants of SARS-CoV-2 at the Beginning of Their Appearance in the City of Cologne in January and February 2021]

BACKGROUND AND AIM OF THE STUDY

Constantly changing virus variants of the novel coronavirus SARS-CoV-2 pose major challenges to the healthcare system. The aim of the present study was to analyse major outbreaks of the alpha and beta variants in Cologne in order to enable effective and rapid response to new virus variants in future pandemics as well as to derive targeted measures to combat the pandemic.

METHODS

In the observation period from January 22 to February 23, 2021, all individuals testing positive for SARS-CoV-2 and their contact persons who were reported to the Cologne Public Health Department were interviewed by employees of the Public Health Department over the telephone. On the one hand, these data formed the basis for the epidemiological and descriptive comparison of the alpha and beta variants to the previously dominant wild type. On the other hand, they were also the basis for the graphical processing of clusters formed by the two virus variants in the form of so-called timelines. For the present work, all clusters with ≥10 individuals were taken into account for the period under consideration.

RESULTS

Of the 3780 individuals that tested positive for SARS-CoV-2 in Cologne during the observation period, 818 cases were due to the virus variants alpha and beta. The alpha versus the beta variant spread quickly in Cologne despite strict non-pharmaceutical interventions. As part of the cluster analysis, five major outbreak were identified in Cologne during the observation period. The alpha variant clusters included two daycare centers and one monastery, while the beta variant clusters included a communal accommodation for refugees and an old people's and nursing home. With the help of cluster analysis, the core role of the spread of the virus variants examined was shown, especially in the context of the home setting. In addition, a high proportion of cases of unknown infection site/contact was found for the wild type and alpha variant.

CONCLUSION

Cluster analyses are an extremely useful tool in the determination of infection sites/contacts and transmission paths as well as in determining existing protective measures and hygiene concepts. Since clusters are to be regarded as the most unfavorable spread scenario, cluster analyses provide important suggestions for modifying further action for both this, as well as for future pandemics.

How public health authorities can use pathogen genomics in health protection practice: a consensus-building Delphi study conducted in the United Kingdom

Pathogen sequencing guided understanding of SARS-CoV-2 evolution during the COVID-19 pandemic. Many health systems developed pathogen genomics services to monitor SARS-CoV-2. There are no agreed guidelines about how pathogen genomic information should be used in public health practice. We undertook a modified Delphi study in three rounds to develop expert consensus statements about how genomic information should be used. Our aim was to inform health protection policy, planning and practice. Participants were from organisations that produced or used pathogen genomics information in the United Kingdom. The first round posed questions derived from a rapid literature review. Responses informed statements for the subsequent rounds. Consensus was accepted when 70 % or more of the responses were strongly agree/agree, or 70 % were disagree/strongly disagree on the five-point Likert scale. Consensus was achieved in 26 (96 %) of 27 statements. We grouped the statements into six categories: monitoring the emergence of new variants; understanding the epidemiological context of genomic data; using genomic data in outbreak risk assessment and risk management; prioritising the use of limited sequencing capacity; sequencing service performance; and sequencing service capability. The expert consensus statements will help guide public health authorities and policymakers to integrate pathogen genomics in health protection practice.

Rapid System to Detect Variants of SARS-CoV-2 in Nasopharyngeal Swabs

Currently, the reference method for identifying the presence of variants of SARS-CoV-2 is whole genome sequencing. Although it is less expensive than in the past, it is still time-consuming, and interpreting the results is difficult, requiring staff with specific skills who are not always available in diagnostic laboratories. The test presented in this study aimed to detect, using traditional real-time PCR, the presence of the main variants described for the spike protein of the SARS-CoV-2 genome. The primers and probes were designed to detect the main deletions that characterize the different variants. The amplification targets were deletions in the S gene: 25-27, 69-70, 241-243, and 157-158. In the ORF1a gene, the deletion 3675-3677 was chosen. Some of these mutations can be considered specific variants, while others can be identified by the simultaneous presence of one or more deletions. We avoided using point mutations in order to improve the speed of the test. Our test can help clinical and medical microbiologists quickly recognize the presence of variants in biological samples (particularly nasopharyngeal swabs). The test can also be used to identify variants of the virus that could potentially be more diffusive as well as not responsive to the vaccine.

Comparative Analysis of Library Preparation Approaches for SARS-CoV-2 Genome Sequencing on the Illumina MiSeq Platform

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been responsible for over two years of the COVID-19 pandemic and a global health emergency. Genomic surveillance plays a key role in overcoming the ongoing COVID-19 pandemic despite its relative successive waves and the continuous emergence of new variants. Many technological approaches are currently applied for the whole genome sequencing (WGS) of SARS-CoV-2. They differ in key stages of the process, and they feature some differences in genomic coverage, sequencing depth, and in the accuracy of variant-calling options. In this study, three different protocols for SARS-CoV-2 WGS library construction are compared: an amplicon-based protocol with a commercial primer panel; an amplicon-based protocol with a custom panel; and a hybridization capture protocol. Specific differences in sequencing depth and genomic coverage as well as differences in SNP number were found. The custom panel showed suitable results and a predictable output applicable for the epidemiological surveillance of SARS-CoV-2 variants.

Towards precision medicine: Omics approach for COVID-19

The coronavirus disease 2019 (COVID-19) pandemic had a devastating impact on human society. Beginning with genome surveillance of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the development of omics technologies brought a clearer understanding of the complex SARS-CoV-2 and COVID-19. Here, we reviewed how omics, including genomics, proteomics, single-cell multi-omics, and clinical phenomics, play roles in answering biological and clinical questions about COVID-19. Large-scale sequencing and advanced analysis methods facilitate COVID-19 discovery from virus evolution and severity risk prediction to potential treatment identification. Omics would indicate precise and globalized prevention and medicine for the COVID-19 pandemic under the utilization of big data capability and phenotypes refinement. Furthermore, decoding the evolution rule of SARS-CoV-2 by deep learning models is promising to forecast new variants and achieve more precise data to predict future pandemics and prevent them on time.

Reduced Seasonal Coronavirus Antibody Responses in Children Following COVID-19 Mitigation Measures, The Netherlands

SARS-CoV-2 prevention and control measures did not only impact SARS-CoV-2 circulation, but also the timing and prevalence of other seasonal respiratory viruses. Especially in children, information on exposure and infections to seasonal coronaviruses as well as SARS-CoV-2 in the first year of the pandemic is largely lacking. Therefore, we set up a one-year serological survey in a large tertiary hospital in the Netherlands. We show that seasonal coronavirus seroprevalence significantly decreased in 2021 in children less than one year, most likely due to COVID-19 control measures. The SARS-CoV-2 seroprevalence in children and adolescents increased from 0.4% to 11.3%, the highest in adolescents. This implies higher exposure rates in adolescents as compared to the general population (>18 years old). It is clear that there have been significant changes in the circulation and subsequent immunity against most respiratory pathogens as a result of the mitigation measures. The implications on shorter as well as longer term are still largely unknown, but the impact of the SARS-CoV-2 pandemic and subsequent control measures will continue to affect the dynamics of other pathogens.

An Automated Strategy to Handle Antigenic Variability in Immunisation Protocols, Part I: Nanopore Sequencing of Infectious Agent Variants

Infectious agents often challenge therapeutics, from antibiotics resistance to antigenic variability affecting inoculation measures. Over the last decades, genome sequencing arose as an important ally to address such challenges. In bacterial infection, whole-genome-sequencing (WGS) supports tracking pathogenic alterations affecting the human microbiome. In viral infection, the analysis of the relevant sequence of nucleotides helps with determining historical variants of a virus and elucidates details about infection clusters and their distribution. Additionally, genome sequencing is now an important step in inoculation protocols, isolating target genes to design more robust immunisation assays. Ultimately, genetic engineering has empowered repurposing at scale, allowing long-lasting repeating clinical trials to be automated within a much shorter time-frame, by adjusting existing protocols. This is particularly important during sanitary emergencies as the ones caused by the 2014 West African Ebola outbreak, the Zika virus rapid spread in both South and North America in 2015, followed by Asia in 2016, and the pandemic caused by the SARS-CoV-2, which has infected more than 187 million people and caused more than 4 million deaths, worldwide, as per July 2021 statistics. In this scenery, this chapter presents a novel fully automated strategy to handle antigenic variability in immunisation protocols. The methodology comprises of two major steps (1) nanopore sequencing of infectious agent variants - the focus is on the SARS-CoV-2 and its variants; followed by (2) mRNA vector design for immunotherapy. This chapter presents the nanopore sequencing step and Chapter 17 introduces a protocol for mRNA vector design.

SARS-CoV-2 outbreaks in secondary school settings in the Netherlands during fall 2020; silent circulation

BACKGROUND

In fall 2020 when schools in the Netherlands operated under a limited set of COVID-19 measures, we conducted outbreaks studies in four secondary schools to gain insight in the level of school transmission and the role of SARS-CoV-2 transmission via air and surfaces.

METHODS

Outbreak studies were performed between 11 November and 15 December 2020 when the wild-type variant of SARS-CoV-2 was dominant. Clusters of SARS-CoV-2 infections within schools were identified through a prospective school surveillance study. All school contacts of cluster cases, irrespective of symptoms, were invited for PCR testing twice within 48 h and 4-7 days later. Combined NTS and saliva samples were collected at each time point along with data on recent exposure and symptoms. Surface and active air samples were collected in the school environment. All samples were PCR-tested and sequenced when possible.

RESULTS

Out of 263 sampled school contacts, 24 tested SARS-CoV-2 positive (secondary attack rate 9.1%), of which 62% remained asymptomatic and 42% had a weakly positive test result. Phylogenetic analysis on 12 subjects from 2 schools indicated a cluster of 8 and 2 secondary cases, respectively, but also other distinct strains within outbreaks. Of 51 collected air and 53 surface samples, none were SARS-CoV-2 positive.

CONCLUSION

Our study confirmed within school SARS-CoV-2 transmission and substantial silent circulation, but also multiple introductions in some cases. Absence of air or surface contamination suggests environmental contamination is not widespread during school outbreaks.

Potential antiviral peptides targeting the SARS-CoV-2 spike protein

BACKGROUND

The coronavirus disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection became an international pandemic and created a public health crisis. The binding of the viral Spike glycoprotein to the human cell receptor angiotensin-converting enzyme 2 (ACE2) initiates viral infection. The development of efficient treatments to combat coronavirus disease is considered essential.

METHODS

An in silico approach was employed to design amino acid peptide inhibitor against the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein. The designed inhibitor (SARS-CoV-2 PEP 49) consists of amino acids with the α1 helix and the β4 - β5 sheets of ACE2. The PEP-FOLD3 web tool was used to create the 3D structures of the peptide amino acids. Analyzing the interaction between ACE2 and the RBD of the Spike protein for three protein data bank entries (6M0J, 7C8D, and 7A95) indicated that the interacting amino acids were contained inside two regions of ACE2: the α1 helical protease domain (PD) and the β4 - β5 sheets.

RESULTS

Molecular docking analysis of the designed inhibitor demonstrated that SARS-CoV-2 PEP 49 attaches directly to the ACE2 binding site of the Spike protein with a binding affinity greater than the ACE2, implying that the SARS-CoV-2 PEP 49 model may be useful as a potential RBD binding blocker.

Genomic analysis of SARS-CoV-2 variants of concern circulating in Hawai'i to facilitate public-health policies

Using genomics, bioinformatics and statistics, herein we demonstrate the effect of statewide and nationwide quarantine on the introduction of SARS-CoV-2 variants of concern (VOC) in Hawai'i. To define the origins of introduced VOC, we analyzed 260 VOC sequences from Hawai'i, and 301,646 VOC sequences worldwide, deposited in the GenBank and global initiative on sharing all influenza data (GISAID), and constructed phylogenetic trees. The trees define the most recent common ancestor as the origin. Further, the multiple sequence alignment used to generate the phylogenetic trees identified the consensus single nucleotide polymorphisms in the VOC genomes. These consensus sequences allow for VOC comparison and identification of mutations of interest in relation to viral immune evasion and host immune activation. Of note is the P71L substitution within the E protein, the protein sensed by TLR2 to produce cytokines, found in the B.1.351 VOC may diminish the efficacy of some vaccines. Based on the phylogenetic trees, the B.1.1.7, B.1.351, B.1.427, and B.1.429 VOC have been introduced in Hawai'i multiple times since December 2020 from several definable geographic regions. From the first worldwide report of VOC in GenBank and GISAID, to the first arrival of VOC in Hawai'i, averages 320 days with quarantine, and 132 days without quarantine. As such, the effect of quarantine is shown to significantly affect the time to arrival of VOC in Hawai'i. Further, the collective 2020 quarantine of 43-states in the United States demonstrates a profound impact in delaying the arrival of VOC in states that did not practice quarantine, such as Utah. Our data demonstrates that at least 76% of all definable SARS-CoV-2 VOC have entered Hawai'i from California, with the B.1.351 variant in Hawai'i originating exclusively from the United Kingdom. These data provide a foundation for policy-makers and public-health officials to apply precision public health genomics to real-world policies such as mandatory screening and quarantine.

A comparison of high-throughput SARS-CoV-2 sequencing methods from nasopharyngeal samples

The COVID-19 pandemic caused by the new Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) continues to threaten public health and burden healthcare systems worldwide. Whole SARS-CoV-2 genome sequencing has become essential for epidemiological monitoring and identification of new variants, which could represent a risk of increased transmissibility, virulence, or resistance to vaccines or treatment. Different next-generation sequencing approaches are used in SARS-CoV-2 sequencing, although with different ability to provide whole genome coverage without gaps and to reliably detect new variants. In this study, we compared the performance of three target enrichment methods (two multiplex amplification methods and one hybridization capture) using nasopharyngeal swabs from infected individuals. We applied these target enrichment methods to the same set of nasopharyngeal samples (N = 93) in high-throughput mode. SARS-CoV-2 genome was obtained using short-read next-generation sequencing. We observed that each method has some advantages, such as high mapping rate (CleanPlex and COVIDSeq) or absence of systematic variant calling error (SureSelect) as well as their limitations such as suboptimal uniformity of coverage (CleanPlex), high cost (SureSelect) or supply shortages (COVIDSeq). Nevertheless, each of the three target enrichment kits tested in this study yielded acceptable results of whole SARS-CoV-2 genome sequencing and either of them can therefore be used in prospective programs of genomic surveillance of SARS-CoV-2. Genomic surveillance will be crucial to overcoming the ongoing pandemic of COVID-19, despite its successive waves and continually emerging variants.

Optimized conditions for Listeria, Salmonella and Escherichia whole genome sequencing using the Illumina iSeq100 platform with point-and-click bioinformatic analysis

Whole-genome sequencing (WGS) data have become an integral component of public health investigations and clinical diagnostics. Still, many veterinary diagnostic laboratories cannot afford to implement next generation sequencing (NGS) due to its high cost and the lack of bioinformatic knowledge of the personnel to analyze NGS data. Trying to overcome these problems, and make NGS accessible to every diagnostic laboratory, thirteen veterinary diagnostic laboratories across the United States (US) initiated the assessment of Illumina iSeq100 sequencing platform for whole genome sequencing of important zoonotic foodborne pathogens Escherichia coli, Listeria monocytogenes, and Salmonella enterica. The work presented in this manuscript is a continuation of this multi-laboratory effort. Here, seven AAVLD accredited diagnostic laboratories explored a further reduction in sequencing costs and the usage of user-friendly platforms for genomic data analysis. Our investigation showed that the same genomic library quality could be achieved by using a quarter of the recommended reagent volume and, therefore a fraction of the actual price, and confirmed that Illumina iSeq100 is the most affordable sequencing technology for laboratories with low WGS demand. Furthermore, we prepared step-by-step protocols for genomic data analysis in three popular user-friendly software (BaseSpace, Geneious, and GalaxyTrakr), and we compared the outcomes in terms of genome assembly quality, and species and antimicrobial resistance gene (AMR) identification. No significant differences were found in assembly quality, and the three analysis methods could identify the target bacteria species. However, antimicrobial resistance genes were only identified using BaseSpace and GalaxyTrakr; and GalaxyTrakr was the best tool for this task.

Global disparities in SARS-CoV-2 genomic surveillance

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

Population-based screening in a municipality after a primary school outbreak of the SARS-CoV-2 Alpha variant, the Netherlands, December 2020-February 2021

An outbreak of SARS-CoV-2 Alpha variant (Pango lineage B.1.1.7) was detected at a primary school (School X) in Lansingerland, the Netherlands, in December 2020. The outbreak was studied retrospectively, and population-based screening was used to assess the extent of virus circulation and decelerate transmission. Cases were SARS-CoV-2 laboratory confirmed and were residents of Lansingerland (November 16th 2020 until February 22th 2021), or had an epidemiological link with School X or neighbouring schools. The SARS-CoV-2 variant was determined using variant PCR or whole genome sequencing. A questionnaire primarily assessed clinical symptoms. A total of 77 Alpha variant cases were found with an epidemiological link to School X, 16 Alpha variant cases linked to the neighbouring schools, and 146 Alpha variant cases among residents of Lansingerland without a link to the schools. The mean number of self-reported symptoms was not significantly different among Alpha variant infected individuals compared to non-Alpha infected individuals. The secondary attack rate (SAR) among Alpha variant exposed individuals in households was 52% higher compared to non-Alpha variant exposed individuals (p = 0.010), with the mean household age, and mean number of children and adults per household as confounders. Sequence analysis of 60 Alpha variant sequences obtained from cases confirmed virus transmission between School X and neighbouring schools, and showed that multiple introductions of the Alpha variant had already taken place in Lansingerland at the time of the study. The alpha variant caused a large outbreak at both locations of School X, and subsequently spread to neighbouring schools, and households. Population-based screening (together with other public health measures) nearly stopped transmission of the outbreak strain, but did not prevent variant replacement in the Lansingerland municipality.

Dynamic SARS-CoV-2 emergence algorithm for rationally-designed logical next-generation vaccines

SARS-CoV-2 worldwide spread and evolution has resulted in variants containing mutations resulting in immune evasive epitopes that decrease vaccine efficacy. We acquired SARS-CoV-2 positive clinical samples and compared the worldwide emerged spike mutations from Variants of Concern/Interest, and developed an algorithm for monitoring the evolution of SARS-CoV-2 in the context of vaccines and monoclonal antibodies. The algorithm partitions logarithmic-transformed prevalence data monthly and Pearson's correlation determines exponential emergence of amino acid substitutions (AAS) and lineages. The SARS-CoV-2 genome evaluation indicated 49 mutations, with 44 resulting in AAS. Nine of the ten most worldwide prevalent (>70%) spike protein changes have Pearson's coefficient r > 0.9. The tenth, D614G, has a prevalence >99% and r-value of 0.67. The resulting algorithm is based on the patterns these ten substitutions elucidated. The strong positive correlation of the emerged spike protein changes and algorithmic predictive value can be harnessed in designing vaccines with relevant immunogenic epitopes. Monitoring, next-generation vaccine design, and mAb clinical efficacy must keep up with SARS-CoV-2 evolution, as the virus is predicted to remain endemic.

Longitudinal Household Assessment of Respiratory Illness in Children and Parents During the COVID-19 Pandemic

IMPORTANCE

In the early COVID-19 pandemic, SARS-CoV-2 testing was only accessible and recommended for symptomatic persons or adults. This restriction hampered assessment of the true incidence of SARS-CoV-2 infection in children as well as detailed characterization of the SARS-CoV-2 disease spectrum and how this spectrum compared with that of other common respiratory illnesses.

OBJECTIVE

To estimate the community incidence of SARS-CoV-2 infection in children and parents and to assess the symptoms and symptom severity of respiratory illness episodes involving SARS-CoV-2-positive test results relative to those with SARS-CoV-2-negative test results.

DESIGN, SETTING, AND PARTICIPANTS

This cohort study randomly selected Dutch households with at least 1 child younger than 18 years. A total of 1209 children and adults from 307 households were prospectively followed up between August 25, 2020, and July 29, 2021, covering the second and third waves of the COVID-19 pandemic. Participation included SARS-CoV-2 screening at 4- to 6-week intervals during the first 23 weeks of participation (core study period; August 25, 2020, to July 29, 2021). Participants in all households finishing the core study before July 1, 2021, were invited to participate in the extended follow-up and to actively report respiratory symptoms using an interactive app until July 1, 2021. At new onset of respiratory symptoms or a SARS-CoV-2 positive test result, a household outbreak study was initiated, which included daily symptom recording, repeated polymerase chain reaction testing (nose-throat swabs and saliva and fecal samples), and SARS-CoV-2 antibody measurement (paired dried blood spots) in all household members. Outbreaks, households, and episodes of respiratory illness were described as positive or negative depending on SARS-CoV-2 test results. Data on participant race and ethnicity were not reported because they were not uniformly collected in the original cohorts and were therefore not representative or informative.

EXPOSURES

SARS-CoV-2-positive and SARS-CoV-2-negative respiratory illness episodes.

MAIN OUTCOMES AND MEASURES

Age-stratified incidence rates, symptoms, and symptom severity for SARS-CoV-2-positive and SARS-CoV-2-negative respiratory illness episodes.

RESULTS

Among 307 households including 1209 participants (638 female [52.8%]; 403 [33.3%] aged <12 years, 179 [14.8%] aged 12-17 years, and 627 [51.9%] aged ≥18 years), 183 household outbreaks of respiratory illness were observed during the core study and extended follow-up period, of which 63 (34.4%) were SARS-CoV-2 positive (59 outbreaks [32.2%] during the core study and 4 outbreaks [2.2%] during follow-up). SARS-CoV-2 incidence was similar across all ages (0.24/person-year [PY]; 95% CI, 0.21-0.28/PY). Overall, 33 of 134 confirmed SARS-CoV-2 episodes (24.6%) were asymptomatic. The incidence of SARS-CoV-2-negative respiratory illness episodes was highest in children younger than 12 years (0.94/PY; 95% CI, 0.89-0.97/PY). When comparing SARS-CoV-2-positive vs SARS-CoV-2-negative respiratory illness episodes in children younger than 12 years, no differences were observed in number of symptoms (median [IQR], 2 [2-4] for both groups), symptom severity (median [IQR] maximum symptom severity score, 6 [4-9] vs 7 [6-13]), or symptom duration (median [IQR], 6 [5-12] days vs 8 [4-13] days). However, among adults, SARS-CoV-2-positive episodes had a significantly higher number (median [IQR], 6 [4-8] vs 3 [2-4]), severity (median [IQR] maximum symptom severity score, 15 [9-19] vs 7 [6-11]), and duration (median [IQR] 13 [8-29] days vs 5 [3-11] days; P < .001 for all comparisons) of symptoms vs SARS-CoV-2-negative episodes.

CONCLUSIONS AND RELEVANCE

In this cohort study, during the first pandemic year when mostly partial or full in-person learning occurred, the SARS-CoV-2 incidence rate in children was substantially higher than estimated from routine testing or seroprevalence data and was similar to that of adult household members. Unlike in unvaccinated adults, SARS-CoV-2 symptoms and symptom severity in children were similar to other common respiratory illnesses. These findings may prove useful when developing pediatric COVID-19 vaccine recommendations.

COVID-19 Whole-Genome Resequencing with Redundant Tiling PCR and Subtract-Based Amplicon Normalization Successfully Characterized SARS-CoV-2 Variants in Clinical Specimens

With an increasing number of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) sequences gathered worldwide, we recognize that deletion mutants and nucleotide substitutions that may affect whole-genome sequencing are accumulating. Here, we propose an additional strategy for tiling PCR for whole-genome resequencing, which can make the pipeline robust for mutations at the primer annealing site by a redundant amplicon scheme. We further demonstrated that subtracting overrepresented amplicons from the multiplex PCR products reduced the bias of the next-generation sequencing (NGS) library, resulting in decreasing required sequencing reads per sample. We applied this sequencing strategy to clinical specimens collected in Bangladesh. More than 80% out of the 304 samples were successfully sequenced. Less than 5% were ambiguous nucleotides, and several known variants were detected. With the additional strategies presented here, we believe that whole-genome resequencing of SARS-CoV-2 from clinical samples can be optimized.

Regional importation and asymmetric within-country spread of SARS-CoV-2 variants of concern in the Netherlands

Background:

Variants of concern (VOCs) of SARS-CoV-2 have caused resurging waves of infections worldwide. In the Netherlands, the Alpha, Beta, Gamma, and Delta VOCs circulated widely between September 2020 and August 2021. We sought to elucidate how various control measures, including targeted flight restrictions, had impacted the introduction and spread of these VOCs in the Netherlands.

Methods:

We performed phylogenetic analyses on 39,844 SARS-CoV-2 genomes collected under the Dutch national surveillance program.

Results:

We found that all four VOCs were introduced before targeted flight restrictions were imposed on countries where the VOCs first emerged. Importantly, foreign introductions, predominantly from other European countries, continued during these restrictions. After their respective introductions into the Netherlands, the Alpha and Delta VOCs largely circulated within more populous regions of the country with international connections before asymmetric bidirectional transmissions occurred with the rest of the country and the VOC became the dominant circulating lineage.

Conclusions:

Our findings show that flight restrictions had limited effectiveness in deterring VOC introductions due to the strength of regional land travel importation risks. As countries consider scaling down SARS-CoV-2 surveillance efforts in the post-crisis phase of the pandemic, our results highlight that robust surveillance in regions of early spread is important for providing timely information for variant detection and outbreak control.

Funding:

None.

Highly Divergent SARS-CoV-2 Alpha Variant in Chronically Infected Immunocompromised Person

We detected a highly divergent SARS-CoV-2 Alpha variant in an immunocompromised person several months after the latest detection of the Alpha variant in the Netherlands. The patient was infected for 42 weeks despite several treatment regimens and disappearance of most clinical symptoms. We identified several potential immune escape mutations in the spike protein.

Increased transmissibility of SARS-CoV-2 alpha variant (B.1.1.7) in children: three large primary school outbreaks revealed by whole genome sequencing in the Netherlands

Background

Variant of concern (VOC) SARS-CoV-2 alpha variant (B.1.1.7) was the dominant strain in the Netherlands between March 2021–June 2021. We describe three primary school outbreaks due to the alpha variant using whole genome sequencing with evidence of large-scale transmission among children, teachers and their household contacts.

Method

All outbreaks described were investigated by the South Limburg Public Health Service, the Netherlands. A case was defined as an individual with a real-time polymerase chain reaction test or antigen test positive for SARS-CoV-2. Whole genome sequencing was performed on random samples from at least one child and one teacher of each affected class.

Results

Peak attack rates in classes were 53%, 33% and 39%, respectively. Specific genotypes were identified for each school across a majority of affected classes. Attack rates were high among staff members, likely to promote staff-to-children transmission. Cases in some classes were limited to children, indicating child-to-child transmission. At 39%, the secondary attack rate (SAR) in household contacts of infected children was remarkably high, similar to SAR in household contacts of staff members (42%). SAR of household contacts of asymptomatic children was only 9%.

Conclusion

Our findings suggest increased transmissibility of the alpha variant in children compared to preceding non-VOC variants, consistent with a substantial rise in the incidence of cases observed in primary schools and children aged 5–12 since the alpha variant became dominant in March 2021. Lack of mandatory masking, insufficient ventilation and lack of physical distancing also probably contributed to the school outbreaks. The rise of the delta variant (B.1.617.2) since July 2021 which is estimated to be 55% more transmissible than the alpha variant, provides additional urgency to adequate infection prevention in school settings.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07623-9.