Preprint
Article

This version is not peer-reviewed.

From Airport to School to Hospital – the Story of a Rare SARS-CoV-2 Virus Variant Import Told with the Help of Patient Swabs and Wastewater-Based Epidemiology

A peer-reviewed article of this preprint also exists.

Submitted:

12 February 2025

Posted:

13 February 2025

You are already at the latest version

Abstract
In December 2021, an outbreak of the SARS-CoV-2 B.1.640.2 variant, potentially originating from Cameroon, was investigated among schoolchildren in Germany. The index case, an adult who had recently returned from a three-week stay in the Republic of Congo, introduced the variant into a school setting via their children, resulting in subsequent transmission within the school and ultimately to a hospital ward. Whole-genome sequencing of viral samples identified both B.1.640.1 and B.1.640.2 lineages. This outbreak highlights the unpredictable nature of emerging SARS-CoV-2 variants and emphasizes the importance of early detection and containment to mitigate transmission to high-risk populations. Notably, wastewater surveillance detected the variant during the study peri-od, reinforcing the utility of wastewater-based epidemiology as a complementary tool for early warning and containment of novel variants. These findings underline the critical need for timely research and adherence to quarantine measures to enhance outbreak control efforts.
Keywords: 
;  ;  

1. Introduction

The COVID-19 pandemic has highlighted the importance of understanding viral variants and their potential introduction into new geographical areas and vulnerable institutions. Due to their impact on transmission rates, immune escape, etc., the emergence of new viral variants poses a significant challenge in containing and managing the COVID-19 pandemic. This study investigates a COVID-19 outbreak in a high school in Germany caused by the B.1.640.2 variant, a novel variant at a time when the Delta (B.1.617.2) variant represented 98% of infections in Germany. The study aims to shed light on the introduction and transmission dynamics of the B.1.640.2 variant in this specific setting and its impact on various institutions. The findings of this study contribute to our understanding of the unpredictability of SARS-CoV-2 variants and the importance of early detection and containment measures to prevent further spread and mitigate the impact on vulnerable populations. Monitoring emerging viral variants such as B.1.640.2 is an effective tool for responding to and managing the pandemics. Accurate surveillance and research on viral variants like B.1.640.2 are crucial for public health authorities to make informed decisions regarding testing, diagnostics, treatment, and vaccine development. It is also crucial to implement public health policies, strategies and interventions to prevent the spread of COVID-19, especially in vulnerable institutions such as schools and hospitals. In addition, it plays a pivotal role in understanding the transmission dynamics and potential impact of these variants, as well as informing targeted interventions and strategies to mitigate their spread and protect vulnerable populations. Furthermore, our study demonstrates that complementing clinical testing with wastewater-based epidemiology (WBE) will enhance our ability to observe emerging variants, building the basis for earlier and more targeted mitigation measures.

1.1. Background

SARS-CoV-2 variants have repeatedly emerged over the last years and have since become a major health concern. Thus, the WHO monitors all emerging SARS-CoV-2 virus variants in the world. While some assume that Viruses mutate and become more infectious and less deadly over time [1-3], current research demonstrates that the human immune status underlies this observation [4]. Novel variants might gain the ability to evade the host’s immune defenses and therefore need to be under surveillance, as variants with a high number of mutations are especially alarming due to their unexpected behavior, making them much more transmissible and capable of spreading more rapidly. Omicron followed this pattern. The original Omicron variant (BA.1/ B.1.1.529) was first reported by the WHO and uploaded to the GISAID disease variant database on November 23rd, 2024, and has more than 300,000 sequences in the database to date [5-7].
However, another variant of interest, Pangolin lineage B.1.640, predated Omicron. The world was alarmed because both variants had many mutations in common, making them much more rapidly transmissible [8,9]. The question was whether and to what extent these two variants would spread and shape the spread of infection. This second variant of interest was first identified in France in April 2021, and later in the Republic of Congo (Brazzaville) in September and was further involved in a cluster of cases in southeastern France from a vaccinated traveler returning from Cameroon around mid-October 2021 [9]. Until November 2021, a total of 12 cases were detected in the region, and the variant was named B.1.640 and uploaded to GISAID on Nov. 4th, 2021 [10].
In addition, French researchers detected an atypical combination of mutations within the B.1.640 variant at the beginning of December 2021 in samples from twelve people infected with Corona, including five children under the age of 15 [10]. The atypical variant was identified by screening with qPCR assays from the Delta (L452R-positive) and Omicron (L452R-negative and negative for S gene detection by the TaqPath COVID-19 assay) variants that co-circulated in Europe at that time [11].
This new lineage, known since December 12th, 2021, was named B.1.640.2. It is a phylogenetic sister group to the B.1.640 lineage renamed B.1.640.1. B.1.640.2 has up to 46 mutations that are already known from other virus variants. These include the mutation N501Y, which is known from the alpha, beta, gamma and (other) omicron variants. It enables the coronavirus to bind more strongly to its molecular target on human cells, making it more easily transmissible [12]. The mutations suggest that B.1.640.2 may be more contagious than the original virus discovered in Wuhan. However, at that time, it remained unclear whether this virus variant would prevail. By that time, B.1.640.2 had been detected in only a few sequencing runs [10].
At that time Omicron was about to take over as the predominant variant circulating worldwide but had not yet started to dominate as it would a month later, being responsible for over 75% of sequences reported to GISAID [13]. Even if only a few cases of B.1.640.2 were present, this case report outlines how the virus variant could spread to a high school in Western Germany causing a mid-scale outbreak cluster.

2. Materials and Methods

2.1. Setting

On December 6, 2021, routine screening of healthcare workers at a maximum care hospital in Western Germany led to the detection of the B.1.640.2 sequence. The variant was first identified in France in April 2021 and traced to Congo. According to the Robert Koch Institute, this was only the second identification of this variant in Germany.
The case involved an unvaccinated female nurse working in the neonatal ICU who was exempt from vaccination for medical reasons, necessitating strict adherence to hygiene protocols and regular PCR testing. She developed mild respiratory symptoms on December 5, one day before her positive test. Notably, her 12-year-old daughter had tested positive for SARS-CoV-2 three days prior, on December 2, as part of routine testing at her school in Western Germany, where saliva-based “lollipop” rapid tests were conducted multiple times per week. Several SARS-CoV-2 cases had already been identified in the daughter's class in November and December 2021.
From a virological standpoint, immediate analysis of this variant was essential. From epidemiological and clinical perspectives, tracing the transmission chain of this rare B.1.640.2 variant – previously detected in Europe only in France – was crucial to understanding its introduction and spread within the German school and hospital settings. The index case was retrospectively identified through extensive epidemiological investigations.
The primary setting of the first outbreak with variant B.1.640.2 in Germany was a high school in the Ruhr metropolitan region of North Rhine-Westphalia, Germany. Additionally, several sites, such as a maximum care hospital, have been affected. We monitored the outbreak from the first case we became aware of on December 6, 2021, until containment at the school in mid-January 2022. This period included recruitment, data collection and follow-up. Our main goal was to prevent the spread of this variant at the hospital and within the school setting.
We collected preliminary data, interviewed and informed the school’s principal of the urgent need for strict hygiene precautions and the exclusion of all symptomatic children and contact persons from school.

2.2. Data and Investigation Procedures

Recruitment was performed by passing on information to the school, teachers and parent groups. On this basis, parents and positive individuals contacted us independently and voluntarily.
The inclusion criteria were positivity and close contact with affected class members and/or healthcare workers. Close contacts were defined as persons having more than 15 minutes of face-to-face contact in and outdoors and family/household members with a confirmed COVID-19 case during the infectious period.
According to the principles of the Declaration of Helsinki, approved by the local institutional research ethics board (ethics vote 468/20), informed consent was given by all participants, in the case of minors, by the signature of at least one parent.
Participation criteria were recorded in personal interviews following a structured questionnaire. All participants who agreed to participate were visited at their homes. Data on sociodemographic factors, exposures, medical history, comorbidities, vaccination schedules, and epidemiological data were collected for all participants. Furthermore, study participants were asked to provide nasopharyngeal and deep throat swabs.
The participants were sampled up to 3 times within a fortnight at various time points during periods of illness. Initially, all children from the affected fifth grade who became ill and/or tested positive by rapid test were tested by the volunteer fire department on December 10, 2021. In this step, nasopharyngeal and deep throat swabs were collected using M 40 compliant Transwab systems (Medical Wire, Corsham, UK).
In the next testing steps i) one week later, December 17, ii) on December 23, and/or iii) on December 27 and/or iv) mid-January 2022, throat swabs for virus diagnostics were obtained using liquid medium swabs (UTM, Befl. swabs, Minitip m. 3 mL UTM (16x100mm) COPAN, Brescia, Italy). Swabs were stored in UTM Viral Stabilization Media at 4 °C at the study center for up to four hours. Swabs were resuspended at a maximum of 6 h after collection, and all laboratory analyses were performed within 48 h.

2.3. Laboratory Analyses

2.3.1. Detection of SARS-CoV-2 RNA

SARS-CoV-2 RNA was qualitatively detected by the real-time PCR SARS-CoV-2 AMP Kit (Abbott Diagnostics, Wiesbaden, Germany), a dual-target assay for the viral RdRp and N genes, which was run on the Alinity m system with automated sample preparation (Abbott Molecular, 2021). Alternatively, viral nucleic acids were extracted by MagNA Pure 96 DNA and Viral NA Small Volume Kit (Roche Diagnostics, Mannheim, Germany) and afterwards reverse-transcribed and amplified using RealStar SARS-CoV-2 RT-PCR Kit 1.0 (Altona Diagnostics, Hamburg, Germany) on the CFX 96 system (Bio-Rad, Feldkirchen, Germany).

2.3.2. Molecular Characterization of SARS-CoV-2 Isolates

Melting curve assays: Probe-based melting curve analyses were used to screen for common mutations in the viral genome and to provisionally characterize SARS-CoV-2 isolates. The following VirSNiP SARS-CoV-2 Spike assays, provided by TIB MOLBIOL (Berlin, Germany) and run on a Roche LC 480 System, were used: del HV69/70, S371L/S373P, E484K, and N501Y. In the resulting curves, a higher melting peak generally indicates the presence of the mutated variant, whereas wild-type sequences show a peak at lower melting temperatures.
Sanger sequencing and data analysis: Extracted viral RNA was amplified using the Super-Script™ III One-step-RT-PCR System with Platinum™ Taq DNA Polymerase (Invitrogen, Waltham, USA) with specific primers from two regions of the SARS-CoV-2 spike gene (Supplementary Table 1) and the following cycling conditions: initial RT step 30 min. at 55°C followed by an initial denaturation for 2 min. at 94°C, 35 cycles of denaturation 15 sec. at 94°C, annealing 30 sec. at 60°C and elongation 1 min. at 68°C, final elongation was 5 min. at 68 °C. In a second step, a nested-PCR with gene specific primers (Supplementary Table 1) was per-formed using the Platinum Taq DNA Polymerase HF (Invitrogen) with following cycling conditions: initial denaturation 2 min. at 94°C followed by 35 cycles of denaturation 15 sec. 94°C, annealing 30 sec. at 60°C and elongation 1 min. at 68°C, final elongation was 5 min. at 68°C. The PCR amplicons were quality checked on a 1 % agarose gel, purified using the QIAquick PCR purification Kit (Qiagen, Hilden, Germany) and sent for Sanger sequencing to Eurofins Genomics (Ebersberg, Germany).
Sequence editing, generation and translation of the consensus sequence into the corresponding amino acid sequence were done using Geneious Pro 5.1.7. A multiple sequence alignment was done by Clustal Omega using the SARS-CoV-2 isolate Wuhan-HU-1 (NC_045512.2) as a reference. Visualization of the Clustal Omega alignment was done by a Highlighter analysis using the Los Alamos National Laboratory pathogen database. Initial manual lineage prediction was performed based on mutations identified in a multiple sequence alignment using out-break.info 14. Whole genome sequencing (see below) later confirmed the initial lineage assignment.
A phylogenetic tree was inferred using the Maximum Likelihood method and the Tamura-Nei nucleotide substitution model based on the receptor binding domain of the viral spike gene (nucleotides 963 to 1737). The tree was drawn to scale, with branch lengths measured in the number of substitutions per site. Statistical robustness was tested using the bootstrap approach with 1000 replicates. Analysis was conducted in MEGA, version 11 15.
The partial SARS-CoV-2 spike gene sequences obtained from the infected individuals in this study were deposited in the GISAID database under the accession numbers EPI_ISL_10033340 - EPI_ISL_10033367.

2.3.3. Next Generation Sequencing and Data Analysis

After generating cDNA from the extracted viral RNA using the LunaScript RT SuperMix Kit (NEB), SARS-CoV-2 whole genome libraries were obtained using the EasySeq™ SARS-CoV-2 Whole Genome NGS Sequencing kit (Nimagen). Next-generation sequencing was performed on normalized pooled libraries on an Illumina MiSeq instrument employing the V2 chemistry (300 cycles).
Reads with an average Phred quality below 20 and a length below 30 base pairs were excluded, enabling downstream analysis. Data analysis was performed with UnCoVar 16. UnCoVar performs a series of QC steps, initially attempts de-novo assembly and then resorts to co-assembly for recalcitrant samples, it subsequently uses pangolin 17 and Kallisto 18 matching to GISAID 19 to obtain lineage calls. Additionally, Freebayes 20, Delly 21 and Varlociraptor 22 are utilized for variant calling.
Six full-length sequences were obtained for six samples, which were classified by pangolin (v3.1.20) as B.1.640.1. However, a deletion of 9 amino acids in the S protein at locations 135-144 was observed, which Colson et al. describe as indicative of lineage B.1.640.2. Kallisto also supported this assignment, as the majority of reads were quantified as B.1.640.2 in six samples.
The full-length SARS-CoV-2 sequences obtained from the infected individuals in this study were deposited in the GISAID database under the accession numbers EPI_ISL_9402306 – EPI_ISL_9402318.

2.3.4. Viral detection in Wastewater

During this study, wastewater samples were collected within the metropolitan area Ruhr and processed for SARS-CoV-2 detection, as described 23. The sampling area comprises both the high school and maximum care hospital described in this work. Viral variants were identified by a modified version of UnCoVar 16, based on the detection of variant-specific mutations without the attempt of genome assembly.

3. Results

3.1. Laboratory Results

According to the outcome of the melting curve analyses, the SARS-CoV-2 isolates harbored by the infected schoolchildren and their relatives did neither contain the spike gene mutations S371L/S373P (Figure 1 A) nor the amino acid exchange L452R (Figure 1 B). Therefore, they most likely did not represent the emerging Omicron or the dominant Delta variants that were the focus of attention in Germany in December 2021.
Furthermore, the isolates lacked the deletion HV69/70 and the mutation E484K (data not shown) but contained the N501Y mutation (Figure 1 C), so they could not be clearly assigned to the SARS-CoV-2 Alpha, Beta, or Gamma lineages. Consequently, the results of the melting curve assays remained inconclusive.
Direct sequencing of two 804 bp and 825 bp long nucleotide fragments from the SARS-CoV-2 spike gene and subsequent translation into the corresponding amino acid sequences revealed the presence of eight mutations when compared to the wild-type Wuhan-Hu-1 strain (Figure 2).
This characteristic pattern of amino acid exchanges had already been reported for variants from the former SARS-CoV-2 B.1.640 lineage, which was renamed to B.1.640.1 9. Analyses of full-lengths sequences obtained by next-generation sequencing confirmed and extended these findings and the formation of a so-called monophyletic group in a tree, which reconstructs the evolutionary history of the viral isolates (Figure 3) showing that the children and their relatives were infected in a single-source outbreak by a yet rare SARS-CoV-2 B.1.640.1 variant, which possibly had originated in the Republic of Congo 9.

3.2. Detailed Outbreak Course

The outbreak described here was first identified on December 6, 2021, following a positive SARS-CoV-2 B.1.640.2 test result in a nurse working in the neonatal ICU. Unvaccinated due to a chronic autoimmune disorder, she had developed mild symptoms one day prior to testing. The investigation suggested that she likely contracted the infection from her teenage daughter, who had shown symptoms three days earlier and tested positive at her school. Due to the nurse’s absence from the neonatal ward during this period, secondary hospital transmission was ruled out.
Following an interview with the initial case, we focused on tracing contacts within the school setting. Between December 6 and December 31, 2021, a total of 20 individuals, including classmates and household contacts, participated in the study, all providing informed consent. Data collection included face-to-face interviews conducted by a hygienist, covering symptom history, clinical signs, prior laboratory results, and potential COVID-19 exposure within familial and social networks. We also collected student schedules to map potential in-school trans-mission pathways.
The first follow-up testing, conducted on December 13, confirmed the isolation of all SARS-CoV-2-positive students and their contacts. Additional infection prevention measures were implemented within the school. Based on interviews with the affected students and their families, we hypothesized that the introduction of the B.1.640.2 variant stemmed from a classmate who was in close contact with the outbreak’s index case. This index case, an adult diagnosed with SARS-CoV-2 via qPCR on November 26, had recently returned from a three-week stay in the Republic of Congo, where he was likely infected with the IHU variant. His initial symptoms appeared on November 25, coinciding with his return journey.
Data from pharyngeal swabs and questionnaires were categorized into three groups for analysis: (i) Full Sets (sample, survey, follow-up interview), (ii) Half Sets (sample and survey), and (iii) Sample or Questionnaire Only. Detailed analyses were performed using only the Full-Set data. The outbreak course is summarized in a flow diagram illustrating transmission within different institutions.
Figure 4. Figure 4 provides a timeline and transmission pathway of the SARS-CoV-2 B.1.640.2 variant outbreak across multiple settings, including family units, a hospital, high school, primary school, kindergarten, and a retirement home. The outbreak originated on November 26, 2021, and was traced back to an individual returning from the Republic of Congo. Circles represent confirmed cases, with larger circles indicating clusters of infections. Transmission links between family members and within institutional settings are shown, demonstrating how the virus spreads through close contacts in schools, a healthcare setting, and a family network. Preventive measures and Christmas holidays achieved containment.
Figure 4. Figure 4 provides a timeline and transmission pathway of the SARS-CoV-2 B.1.640.2 variant outbreak across multiple settings, including family units, a hospital, high school, primary school, kindergarten, and a retirement home. The outbreak originated on November 26, 2021, and was traced back to an individual returning from the Republic of Congo. Circles represent confirmed cases, with larger circles indicating clusters of infections. Transmission links between family members and within institutional settings are shown, demonstrating how the virus spreads through close contacts in schools, a healthcare setting, and a family network. Preventive measures and Christmas holidays achieved containment.
Preprints 149220 g004

3.3. Symptoms and Vaccination Status

Clinical data were available for all confirmed COVID-19 cases linked to the school outbreak and affected healthcare workers. Eight tested positive via PCR at the school between December 1 and December 13, 2021. An additional 12 cases, initially negative, turned positive between December 13 and December 28 due to close contact with infected students. Altogether, 20 cases were confirmed, none of which required hospitalization or resulted in fatality.
Mild or asymptomatic courses were observed in 59% of cases. Among symptomatic individuals, the most common symptoms included fatigue (73%), rhinorrhea (68%), sore throat (38%), headache and cough (35%), fever, gastrointestinal symptoms (21%), and loss of taste or smell (12%). No complications, hospitalizations, or severe outcomes were recorded.
Regarding vaccination status, five students aged 12 had received one or two doses of the BNT162b2 vaccine or had previous SARS-CoV-2 infections, while others had no prior COVID-19 history.

3.4. Epidemiological Links

3.4.1. The Index Case and Entrance to Germany

The index case was likely infected with the B.1.640.2 variant during a professional stay in the Republic of Congo. Upon returning to Germany on November 26, 2021, the individual developed mild respiratory symptoms en route, which led to a SARS-CoV-2 diagnosis via real-time qPCR testing on November 26. Given the timeline and the case's travel history, alternative routes of B.1.640.2 introduction into Germany were carefully examined and excluded based on epidemiological and travel data. No other exposures were identified in Germany before symptom onset, reinforcing that the variant was likely contracted during his stay abroad.
This scenario emphasizes the potential for travel-related transmission and the need for vigilant surveillance of travelers returning from regions with emerging SARS-CoV-2 variants. The introduction of B.1.640.2 highlights how even rare variants can establish footholds in the new regions. It underscores the importance of early detection and rapid containment strategies to prevent further spread within high-contact settings like schools and healthcare facilities.

3.4.2. Exclusion of Cases and Transmission Events in a School

Initial transmission within the school cluster appeared to occur during break periods, with the earliest identified transmission event traced to a school break on November 29. According to student reports, nasal swab testing was conducted following this break, revealing subsequent cases. Potential transmission settings, such as staff meetings or parent-teacher events, were excluded from consideration, as no COVID-19 cases were linked to these gatherings.
Several cases were ruled out based on symptom onset and sequencing data. One infected student, an affected sibling and their household contacts without sequencing data were classified as part of a separate outbreak due to distinct symptom timelines. Another sibling, at-tending a kindergarten in a separate outbreak, was similarly excluded. Two students with early symptom onset were also ruled out as index cases. Sequence analysis confirmed that one of these students carried a different lineage (B.1.258), disqualifying them from association with the B.1.640.2 outbreak.
To further confirm the outbreak lineage, a phylogenetic tree was constructed using a custom pipeline and 9471 selected SARS-CoV-2 sequences from both German and international databases (see Methods section). This analysis helped delineate the specific lineage relationships and confirmed the outbreak’s association with the B.1.640.2 variant, supporting targeted containment efforts within the school setting.

3.4.3. Spill over from school to hospital

Infections in the hospital setting occurred before a comprehensive hygiene protocol was fully implemented, likely originating from an asymptomatic healthcare worker. While on weekend duty, this nurse performed blood draws across two different hospital wards where mask-wearing was not yet mandatory. The following Monday, the nurse developed symptoms, including fever and cough, and subsequently tested positive for COVID-19. Within 24 hours of diagnosis, the employee entered home isolation. However, 27 patients across the affected wards developed confirmed COVID-19 infections over the next ten days. After these initial cases, no further in-house transmissions were observed following the rigorous application of a standardized hygiene protocol, highlighting the effectiveness of rapid containment measures.
In the school setting, mask-wearing was only required indoors, and students spent most of their break times outdoors, where masks were not mandated. By December 13, 2021, all identified SARS-CoV-2-positive students, and their contacts had been isolated, with additional preventive measures promptly enforced. The Christmas holiday period contributed to out-break containment by reducing school-based transmission opportunities. Ultimately, the outbreak was contained, and all restrictions were lifted 35 days after the outbreak’s initial onset.
Figure 5. Transmission timeline and pathway of SARS-CoV-2 infections across different biotopes from November 26, 2021, to early January 2022. The diagram illustrates the spread of COVID-19 from an index family through various environments, including a retirement home, hospital, high school, and primary school. Circles represent infected individuals in each setting, with arrows indicating transmission pathways among family members and within institutional environments. Key nodes include the initial infection of a healthcare worker and subsequent spread to patients in the hospital, as well as secondary transmissions within a high school class and primary school.
Figure 5. Transmission timeline and pathway of SARS-CoV-2 infections across different biotopes from November 26, 2021, to early January 2022. The diagram illustrates the spread of COVID-19 from an index family through various environments, including a retirement home, hospital, high school, and primary school. Circles represent infected individuals in each setting, with arrows indicating transmission pathways among family members and within institutional environments. Key nodes include the initial infection of a healthcare worker and subsequent spread to patients in the hospital, as well as secondary transmissions within a high school class and primary school.
Preprints 149220 g005

3.4. Integration of clinical and wastewater monitoring

Wastewater viral genomic surveillance during the study led to the detection of the B.1.640 variant in the inlet of a municipal wastewater treatment plant in the Ruhr metropolitan area. The plant receives wastewater from the northern part of Essen and other nearby cities. Table 1 compares the clinical cases reported for this variant to observations in wastewater.
Surveillance data was obtained from the German Robert Koch Institute as part of the national genomic surveillance from symptomatic patients.
The table compares the national surveillance data and the Essen maximum care hospital, showing that the B.1.640 variant could be identified in wastewater one week before the first clinical case reported in the City of Essen. This suggests that the variant was already circulating in the population before its first clinically detected case, confirming the early warning capacity of WBE.

4. Discussion

This study presents a detailed examination of a COVID-19 outbreak caused by the rare B.1.640.2 variant in a German school. It highlights the complexities of variant containment in high-contact settings like schools and underscores the essential role of comprehensive, adaptive hygiene measures [14]. The outbreak’s progression emphasizes the need for rapid, robust preventive strategies to mitigate SARS-CoV-2 transmission, especially for emerging variants whose behaviors are still largely unknown.
The school’s pre-existing hygiene measures, including twice-weekly rapid antigen testing for staff and daily testing for symptomatic individuals, were aligned with standard educational guidelines and provided an effective framework for early detection [15]. However, the outbreak’s dynamics suggest that more intensive strategies, such as mandatory indoor mask-wearing within all grade cohorts, enhanced ventilation, and stricter distancing measures, might be necessary to manage highly transmissible variants. The combination of asymptomatic spreaders and minimal initial symptoms created an environment conducive to viral propagation, especially in poorly ventilated areas and during break times when distancing and mask-wearing were inconsistent [16-18].
Control of this outbreak heavily relied on immediate isolation protocols and extensive testing of close contacts, effectively curtailing further spread. These measures proved critical in limiting in-school transmissions, even after superspreading events in which mask-wearing and distancing were relaxed. This aligns with findings from other studies, which indicate that consistent mask use and improved ventilation significantly reduce transmission risk, particularly when coupled with comprehensive contact tracing and testing [15,19-21].
This outbreak also demonstrated that even non-variants of concern (VOC) can create substantial disruptions if hygiene protocols are not meticulously applied. Compared to previous variants like Delta, which rapidly impacted various sectors, the B.1.640.2 variant’s effect within the school setting illustrates that rigorous adherence to preventive protocols is necessary for all variant types [22]. In this case, inconsistent indoor ventilation, mask-wearing policies, and distancing contributed to the spread within the outbreak cohort, supporting the need for strict, consistent hygiene measures.
Furthermore, comparing the clinical testing with data from wastewater surveillance demonstrated the presence of this variant in the population before the first confirmed patient. Several WBE studies performed during the COVID-19 pandemic showed the early detection of viral variants in wastewater, as well as suggested higher anticipation of emerging variants when sampling in more minor, intra-urban scales [23-26]. The identification of the B.1.640 variant in a wastewater treatment plant with input from more than a million people places the question of a higher than-identified outbreak, with a large enough number of infected people to be detected by this large-scale sampling. This emphasizes the value of combining clinical surveillance with wastewater surveillance to identify potentially otherwise-overlooked outbreaks and to direct localized actions [27].
In addition, the study’s findings on the biphasic infection pattern, with a secondary viral load peak shortly after quarantine release, raise questions about variant-specific viral kinetics. This observation could have (policy) implications for quarantine duration and post-quarantine testing requirements to reduce reintroduction risks in schools and households. Furthermore, genomic analysis confirmed the variant’s origin, underscoring the utility of sequencing in distinguishing imported cases and informing containment strategies that can prevent wider community spread.
Strengthened protocols, including enforced mask use, optimized ventilation systems, and cohort isolation policies to mitigate close-contact exposure, would benefit future strategies for managing outbreaks in schools and similar settings [28]. These measures, particularly when coupled with ongoing genomic surveillance, as recommended by WHO and public health authorities, are vital in identifying and controlling new variant outbreaks and minimizing the risk of cross-border and community spread [29].
In conclusion, this outbreak illustrates that even rare SARS-CoV-2 variants can initiate significant transmission clusters in schools, underscoring the need for proactive, adaptive outbreak management. The findings highlight the critical role of enhanced hygiene protocols, regular testing, and rapid isolation measures in containing variant outbreaks, especially in educational environments with high contact rates and potential for asymptomatic transmission. Adopting such comprehensive measures as part of standard hygiene school protocols may enhance resilience against future SARS-CoV-2 variants, ultimately protecting both students and the broader community.

5. Conclusions

Based on this low number of cases, it is not possible to speculate on virological, epidemiological or clinical features of this B.1.640.2 variant. However, we need to learn from such outbreak scenarios. Overall, these observations again show the unpredictability of the emergence of new SARS-CoV-2 variants and their introduction from abroad, and they exemplify the difficulty of controlling such introduction and subsequent spread. These data also warrant implementing genomic surveillance of SARS-CoV-2 in geographical areas and at a country scale.
Continued monitoring is needed to further elucidate the role of young children in the ongoing COVID-19 pandemic, especially in environments such as DCCs where close contact between young children could facilitate the spread of SARS-CoV-2. Such surveillance is necessary to guide public health interventions such as screening, vaccination, and temporary closure of DCCs. Nonetheless, it is important to consider again the impossibility of predicting the emergence of a new variant introduction from abroad and then the difficulty of controlling its subsequent spread in a country. It is necessary to ensure rigorous genomic surveillance worldwide as an implementation of the genomic track system to prevent the eventual passage of this variant and its eventual spread in other regions. It is too early to speculate in this way, but we must be alert to this variant and other concerns that can be introduced in a country. Of course, vaccination and government advice must be considered by people to protect themselves and others.
COVID-19 prevention in schools involves studying in small groups and minimizing student mixing in activities and transportation. Teachers and parents should lead by wearing facemasks, maintaining hand hygiene, keeping physical distance, etc. School attendance should be avoided at any sign of illness. Learning from home may also reduce the need for class attendance. Outdoor classes should also be considered. COVID-19 prevention encompasses avoiding the ‘three Cs’: closed spaces with poor ventilation, crowded places and close-contact settings. The European Centre for Disease Prevention and Control’s report on air-conditioning and ventilation systems and COVID-19 recommends increasing air exchange rate and outdoor air use and decreasing air recirculation to reduce spread in indoor spaces [30,31]. Finally, appropriate planning for COVID-19 preparedness and prevention for the next school year is essential.

Author Contributions

Conceptualization, S.S., R.M.S., S.G., M.K., S.C.G.; Methodology, A.T., F.M., I. K., S.S., J.S., S.G., M.K., S.C.G., S.V., and R.M.S.; Software, A.T., F.M., I.K.; Validation, C.E., A.T., F.M., I.K., R.S.R.; Formal Analysis, C.E., U.D., A.T., F.M., I.K., R.S.R., S.C.G., O.E.A.; Investigation, C.E., U.D., A.T., F.M., I.K., R.M.S., S.S., R.S.R., S.R., S.V., S.C.G., O.E.A.; Resources, J.S., S.S., R.M.S., S.G., M.K., U.D.; Data Curation, A.T., F.M., I.K.; Writing – Original Draft Preparation, R.M.S., S.S., U.D., O.E.A.; Writing – Review & Editing, R.M.S., S.S., U.D., O.E.A., F.M.; Visualization, C.E., A.T., F.M., I.K., R.M.S., S.R., S.C.G., O.E.A.; Supervision, U.D., F.M., R.M.S.; Project Administration, S.S., R.M.S. The manuscript was written with input from all co-authors. All authors approved the last version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Ethical review and approval were waived for this study, as no formal ethics approval was re-quired for outbreak investigations under the applicable regulations. This study was conducted as part of a public health response to a SARS-CoV-2 outbreak, where ethical approval is not mandatory. No direct participant interventions were performed, and only anonymized epidemiological and sequencing data were analyzed. No personally identifiable information was collected, ensuring full compliance with ethical and data protection standards.

Informed Consent Statement

Informed consent was waived for this study, as it was conducted as part of an outbreak investigation in accordance with public health regulations. No direct interventions or clinical trials were performed, and all collected data were anonymized epidemiological and sequencing da-ta. No personally identifiable information was recorded or published, ensuring compliance with ethical and data protection standards.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. To ensure compliance with data protection policies, access may be granted to qualified re-searchers for non-commercial scientific purposes upon reasonable request.

Acknowledgments

The authors are grateful to B. Reimer, S. Dziubek, M. Hölper, B. Krellenberg, S. Sarr, S. Tschammer for excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest. OEA wishes to declare that she is currently working for the company Pfizer Hellas S.A. and that the content of the publication does not reflect Pfizer's position.

References

  1. Samieefar, N.; Rashedi, R.; Akhlaghdoust, M.; Mashhadi, M.; Darzi, P.; Rezaei, N. Delta Variant: The New Challenge of COVID-19 Pandemic, an Overview of Epidemiological, Clinical, and Immune Characteristics. Acta Biomed 2022, 93, e2022179. [CrossRef]
  2. Butt, A.A.; Dargham, S.R.; Chemaitelly, H.; Al Khal, A.; Tang, P.; Hasan, M.R.; Coyle, P.V.; Thomas, A.G.; Borham, A.M.; Concepcion, E.G.; et al. Severity of Illness in Persons Infected With the SARS-CoV-2 Delta Variant vs Beta Variant in Qatar. JAMA Intern Med 2022, 182, 197-205. [CrossRef]
  3. Kuehn, B.M. Delta Variant Wasn't Linked With More Severe Disease. JAMA 2021, 326, 2251. [CrossRef]
  4. Winkel, A.; Kozanli, E.; Haverkort, M.E.; Euser, S.M.; Sluiter-Post, J.G.C.; Mariman, R.; Vogelzang, A.; de Bakker, J.; Lap, C.R.; van Houten, M.A.; et al. Lower levels of household transmission of SARS-CoV-2 VOC Omicron compared to Wild-type: an interplay between transmissibility and immune status. J Infect Dis 2024. [CrossRef]
  5. Nyberg, T.; Ferguson, N.M.; Nash, S.G.; Webster, H.H.; Flaxman, S.; Andrews, N.; Hinsley, W.; Bernal, J.L.; Kall, M.; Bhatt, S.; et al. Comparative analysis of the risks of hospitalisation and death associated with SARS-CoV-2 omicron (B.1.1.529) and delta (B.1.617.2) variants in England: a cohort study. Lancet 2022, 399, 1303-1312. [CrossRef]
  6. Veneti, L.; Bøås, H.; Bråthen Kristoffersen, A.; Stålcrantz, J.; Bragstad, K.; Hungnes, O.; Storm, M.L.; Aasand, N.; Rø, G.; Starrfelt, J.; et al. Reduced risk of hospitalisation among reported COVID-19 cases infected with the SARS-CoV-2 Omicron BA.1 variant compared with the Delta variant, Norway, December 2021 to January 2022. Euro Surveill 2022, 27. [CrossRef]
  7. (WHO), W.H.O. Classification of Omicron (B.1.1.529): SARS-CoV-2 Variant of Concern. 2021.
  8. Greene, S.K.; Levin-Rector, A.; Kyaw, N.T.T.; Luoma, E.; Amin, H.; McGibbon, E.; Mathes, R.W.; Ahuja, S.D. Comparative hospitalization risk for SARS-CoV-2 Omicron and Delta variant infections, by variant predominance periods and patient-level sequencing results, New York City, August 2021-January 2022. Influenza Other Respir Viruses 2023, 17, e13062. [CrossRef]
  9. Colson, P.; Delerce, J.; Burel, E.; Beye, M.; Fournier, P.E.; Levasseur, A.; Lagier, J.C.; Raoult, D. Occurrence of a substitution or deletion of SARS-CoV-2 spike amino acid 677 in various lineages in Marseille, France. Virus Genes 2022, 58, 53-58. [CrossRef]
  10. Picard, G.; Fournier, L.; Maisa, A.; Grolhier, C.; Chent, S.; Huchet-Kervalla, C.; Sudour, J.; Pretet, M.; Josset, L.; Behillil, S.; et al. Emergence, spread and characterisation of the SARS-CoV-2 variant B.1.640 circulating in France, October 2021 to February 2022. Euro Surveill 2023, 28. [CrossRef]
  11. Català, M.; Coma, E.; Alonso, S.; Andrés, C.; Blanco, I.; Antón, A.; Bordoy, A.E.; Cardona, P.J.; Fina, F.; Martró, E.; et al. Transmissibility, hospitalization, and intensive care admissions due to omicron compared to delta variants of SARS-CoV-2 in Catalonia: A cohort study and ecological analysis. Front Public Health 2022, 10, 961030. [CrossRef]
  12. Liu, Y.; Liu, J.; Plante, K.S.; Plante, J.A.; Xie, X.; Zhang, X.; Ku, Z.; An, Z.; Scharton, D.; Schindewolf, C.; et al. The N501Y spike substitution enhances SARS-CoV-2 infection and transmission. Nature 2022, 602, 294-299. [CrossRef]
  13. Liddor Naim, M.; Fu, Y.; Shagan, M.; Bar-Or, I.; Marks, R.; Sun, Q.; Granek, R.; Kushmaro, A. The Rise and Fall of Omicron BA.1 Variant as Seen in Wastewater Supports Epidemiological Model Predictions. Viruses 2023, 15. [CrossRef]
  14. Zimmerman, K.O.; Akinboyo, I.C.; Brookhart, M.A.; Boutzoukas, A.E.; McGann, K.A.; Smith, M.J.; Maradiaga Panayotti, G.; Armstrong, S.C.; Bristow, H.; Parker, D.; et al. Incidence and Secondary Transmission of SARS-CoV-2 Infections in Schools. Pediatrics 2021, 147. [CrossRef]
  15. Young, B.C.; Eyre, D.W.; Kendrick, S.; White, C.; Smith, S.; Beveridge, G.; Nonnenmacher, T.; Ichofu, F.; Hillier, J.; Oakley, S.; et al. Daily testing for contacts of individuals with SARS-CoV-2 infection and attendance and SARS-CoV-2 transmission in English secondary schools and colleges: an open-label, cluster-randomised trial. Lancet 2021, 398, 1217-1229. [CrossRef]
  16. Yuan, H.; Reynolds, C.; Ng, S.; Yang, W. Factors affecting the transmission of SARS-CoV-2 in school settings. Influenza Other Respir Viruses 2022, 16, 643-652. [CrossRef]
  17. Allen, J.G.; Ibrahim, A.M. Indoor Air Changes and Potential Implications for SARS-CoV-2 Transmission. JAMA 2021, 325, 2112-2113. [CrossRef]
  18. Pavilonis, B.; Ierardi, A.M.; Levine, L.; Mirer, F.; Kelvin, E.A. Estimating aerosol transmission risk of SARS-CoV-2 in New York City public schools during reopening. Environ Res 2021, 195, 110805. [CrossRef]
  19. Chen, X.; Ran, L.; Liu, Q.; Hu, Q.; Du, X.; Tan, X. Hand Hygiene, Mask-Wearing Behaviors and Its Associated Factors during the COVID-19 Epidemic: A Cross-Sectional Study among Primary School Students in Wuhan, China. International journal of environmental research and public health 2020, 17. [CrossRef]
  20. McLeod, R.S.; Hopfe, C.J.; Bodenschatz, E.; Moriske, H.J.; Poschl, U.; Salthammer, T.; Curtius, J.; Helleis, F.; Niessner, J.; Herr, C.; et al. A multi-layered strategy for COVID-19 infection prophylaxis in schools: A review of the evidence for masks, distancing, and ventilation. Indoor Air 2022, 32, e13142. [CrossRef]
  21. Zivelonghi, A.; Lai, M. Mitigating aerosol infection risk in school buildings: the role of natural ventilation, volume, occupancy and CO2 monitoring. Building and Environment 2021, 204, 108139. [CrossRef]
  22. Foster, A.; Kinzel, M. SARS-CoV-2 transmission in classroom settings: Effects of mitigation, age, and Delta variant. Phys Fluids (1994) 2021, 33, 113311. [CrossRef]
  23. Wilhelm, A.; Widera, M.; Grikscheit, K.; Toptan, T.; Schenk, B.; Pallas, C.; Metzler, M.; Kohmer, N.; Hoehl, S.; Marschalek, R.; et al. Limited neutralisation of the SARS-CoV-2 Omicron subvariants BA.1 and BA.2 by convalescent and vaccine serum and monoclonal antibodies. EBioMedicine 2022, 82, 104158. [CrossRef]
  24. Westhaus, S.; Weber, F.A.; Schiwy, S.; Linnemann, V.; Brinkmann, M.; Widera, M.; Greve, C.; Janke, A.; Hollert, H.; Wintgens, T.; et al. Detection of SARS-CoV-2 in raw and treated wastewater in Germany - Suitability for COVID-19 surveillance and potential transmission risks. Sci Total Environ 2021, 751, 141750. [CrossRef]
  25. Schmiege, D.; Haselhoff, T.; Thomas, A.; Kraiselburd, I.; Meyer, F.; Moebus, S. Small-scale wastewater-based epidemiology (WBE) for infectious diseases and antibiotic resistance: A scoping review. International Journal of Hygiene and Environmental Health 2024, 259. [CrossRef]
  26. Schmiege, D.; Kraiselburd, I.; Haselhoff, T.; Thomas, A.; Doerr, A.; Gosch, J.; Schoth, J.; Teichgraber, B.; Moebus, S.; Meyer, F. Analyzing community wastewater in sub-sewersheds for the small-scale detection of SARS-CoV-2 variants in a German metropolitan area. Sci Total Environ 2023, 898, 165458. [CrossRef]
  27. Smith, T.; Holm, R.H.; Yeager, R.; Moore, J.B.t.; Rouchka, E.C.; Sokoloski, K.J.; Elliott, E.M.; Talley, D.; Arora, V.; Moyer, S.; et al. Combining Community Wastewater Genomic Surveillance with State Clinical Surveillance: A Framework for SARS-CoV-2 Public Health Practice. Food Environ Virol 2022, 14, 410-416. [CrossRef]
  28. Villers, J.; Henriques, A.; Calarco, S.; Rognlien, M.; Mounet, N.; Devine, J.; Azzopardi, G.; Elson, P.; Andreini, M.; Tarocco, N.; et al. SARS-CoV-2 aerosol transmission in schools: the effectiveness of different interventions. Swiss Med Wkly 2022, 152, w30178. [CrossRef]
  29. Diamond, M.B.; Keshaviah, A.; Bento, A.I.; Conroy-Ben, O.; Driver, E.M.; Ensor, K.B.; Halden, R.U.; Hopkins, L.P.; Kuhn, K.G.; Moe, C.L.; et al. Wastewater surveillance of pathogens can inform public health responses. Nat Med 2022, 28, 1992-1995. [CrossRef]
  30. Elsaid, A.M.; Mohamed, H.A.; Abdelaziz, G.B.; Ahmed, M.S. A critical review of heating, ventilation, and air conditioning (HVAC) systems within the context of a global SARS-CoV-2 epidemic. Process Saf Environ Prot 2021, 155, 230-261. [CrossRef]
  31. 2020., E.C.f.D.P.a.C.S.E. Heating, ventilation and air-conditioning systems in the context of COVID-19.; ECDC: 10. November 2020 2020.
Figure 1. Results of probe-based melting curve analyses to detect SARS-CoV-2 spike gene mutations S371L/S373P (A), L452R (B), and N501Y (C). Viral RNA extracted from five selected schoolchildren and their relatives (GISAID database accession numbers given) was run in parallel with materials obtained from epidemiologically unlinked patients.
Figure 1. Results of probe-based melting curve analyses to detect SARS-CoV-2 spike gene mutations S371L/S373P (A), L452R (B), and N501Y (C). Viral RNA extracted from five selected schoolchildren and their relatives (GISAID database accession numbers given) was run in parallel with materials obtained from epidemiologically unlinked patients.
Preprints 149220 g001
Figure 2. (A) Schematic overview of the SARS-CoV-2 spike protein (blue) and location of selected protein domains such as N-terminal domain (NTD, orange), receptor binding domain (RBD, red), subdomain 1 and 2 (SD1, SD2, grey) and the S1/S2 cleavage site (yellow) (adapted from https://covdb.stanford.edu/. Last accessed: February 17th, 2022). (B) and (C) Protein alignment of the sequences resulting from the SARS-CoV-2 B.1.640.1 outbreak in a secondary school. Horizontal grey lines correspond to the consensus sequence of each patient (GISAID database accession numbers given) and the SARS-CoV-2 Wuhan-Hu-1 (NC_045512.2) reference sequence, respectively. Amino acid exchanges are depicted in pink. (B) Alignment of about 268 amino acids covering the RBD of the SARS-CoV-2 spike protein (red dashed line). (C) Alignment of about 275 amino acids covering the SD2 and part of the C-terminus of the protein (grey dashed line).
Figure 2. (A) Schematic overview of the SARS-CoV-2 spike protein (blue) and location of selected protein domains such as N-terminal domain (NTD, orange), receptor binding domain (RBD, red), subdomain 1 and 2 (SD1, SD2, grey) and the S1/S2 cleavage site (yellow) (adapted from https://covdb.stanford.edu/. Last accessed: February 17th, 2022). (B) and (C) Protein alignment of the sequences resulting from the SARS-CoV-2 B.1.640.1 outbreak in a secondary school. Horizontal grey lines correspond to the consensus sequence of each patient (GISAID database accession numbers given) and the SARS-CoV-2 Wuhan-Hu-1 (NC_045512.2) reference sequence, respectively. Amino acid exchanges are depicted in pink. (B) Alignment of about 268 amino acids covering the RBD of the SARS-CoV-2 spike protein (red dashed line). (C) Alignment of about 275 amino acids covering the SD2 and part of the C-terminus of the protein (grey dashed line).
Preprints 149220 g002
Figure 3. Phylogenetic tree based on SARS-CoV-2 spike gene RBD sequences (nucleotides 963 to 1737) from 14 infected schoolchildren and their relatives (GISAID database accession numbers given). Five epidemiologically unlinked isolates were included as references (EPI_ISL_3115625 [Alpha], EPI_ISL_1282471 [Beta], EPI_ISL_1849483 [Gamma], EPI_ISL_8038881 [Delta], and NC_045512.2 [SARS-CoV-2 Wuhan-HU-1]). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. Codon positions included were 1st, with 258 amino acid positions in the final dataset.
Figure 3. Phylogenetic tree based on SARS-CoV-2 spike gene RBD sequences (nucleotides 963 to 1737) from 14 infected schoolchildren and their relatives (GISAID database accession numbers given). Five epidemiologically unlinked isolates were included as references (EPI_ISL_3115625 [Alpha], EPI_ISL_1282471 [Beta], EPI_ISL_1849483 [Gamma], EPI_ISL_8038881 [Delta], and NC_045512.2 [SARS-CoV-2 Wuhan-HU-1]). The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. Codon positions included were 1st, with 258 amino acid positions in the final dataset.
Preprints 149220 g003
Table 1. A comparison of clinical cases associated with this variant detected in wastewater samples from the University Hospital Essen and national data from Germany (WWTP: wastewater treatment plant). In 2021, no cases of this variant had been reported in Germany.
Table 1. A comparison of clinical cases associated with this variant detected in wastewater samples from the University Hospital Essen and national data from Germany (WWTP: wastewater treatment plant). In 2021, no cases of this variant had been reported in Germany.
Wastewater Patients
Week (2021) WWTP Hospital Number of B.1.640 cases In Germany
43 - B.1.617.2 no cases
44 B.1.617.2 B.1.617.2 no cases
45 B.1.617.2 B.1.617.2 no cases
46 - B.1.617.2 no cases
47 B.1.617.2
B.1.640
B.1.617.2 1 case
48 B.1.617.2
B.1.640
B.1.617.2
B.1.640
5 cases
49 B.1.617.2
BA.1
B.1.617.2
B.1.640.*
14 cases
50 not analyzed B.1.617.2
B.1.640.*
BA.1
7 cases
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated