Volume 29, Number 12—December 2023
Research Letter
SARS-CoV-2 Variants BQ.1 and XBB.1.5 in Wastewater of Aircraft Flying from China to Denmark, 2023
Abstract
We analyzed wastewater samples from 14 aircraft arriving in Denmark directly from China during January 9–February 12, 2023. Wastewater from 11 aircraft was SARS-CoV-2–positive by PCR; 6 predominantly contained BQ.1 and XBB.1 subvariants. Wastewater-based surveillance can contribute to public health monitoring of SARS-CoV-2 and other emerging infectious agents.
Relaxation of China’s zero-COVID policy in December 2022 led the European Centre for Disease Prevention and Control to recommend several nonpharmaceutical interventions to curb COVID-19 spread and monitor any emerging SARS-CoV-2 variants; those interventions included wastewater-based surveillance (1). We report results of subsequent wastewater surveillance of aircraft arriving at Copenhagen Airport in Copenhagen, Denmark, directly from Beijing or Shanghai, China.
During weeks 2–6 of 2023 (January 9–February 12), a total of 14 aircraft arrived at Copenhagen Airport from China. A service truck extracted waste from the aircraft by using vacuum pressure, after which a rinsing program was performed, and the disinfectant Idu-Flight (Brenntag Nordic A/S, https://www.brenntag.com) was added to the waste tank. Wastewater samples were collected as grab samples from the service truck and immediately transported to Statens Serum Institut in Copenhagen for analysis.
The pH value of the sample material ranged from 9–10 because of the addition of Idu-Flight. Idu-Flight contains the active ingredients glutaraldehyde and benzalkonium chloride; the disinfectant is expected to negatively affect the stability of virus particles and hinder amplification of RNA sequences. We adjusted the samples to pH 7.5–8.5 by using HCl and homogenized them by vigorous vortexing. We split the 14 samples into a total of 43 aliquots and then centrifuged those at either 4,000 × g or 10,000 × g for 10 min to pellet solid material. For the first aliquot from aircraft AC1, we analyzed 10 mL of sample material without any centrifugation; for all other samples, we analyzed 10 mL of supernatant after centrifugation. We purified viruses by using NanoTrap Microbiome A particles (Ceres Nanosciences Inc., https://www.ceresnano.com) and RNA by using Maxwell RSC Cartridges (Promega Corporation, https://www.promega.com). We performed quantitative reverse transcription PCR (qRT-PCR) in technical triplicate by using the GoTaq Enviro kit (Promega) and the US Centers for Disease Control and Prevention N2 primer/probe for SARS-CoV-2 detection (Table; Appendix Table).
Of the 43 qRT-PCR reactions, 31 (72%) were positive for SARS-CoV-2, representing 11 aircraft. We conducted whole-genome sequencing of samples from those 11 aircraft by using the Illumina MiSeq platform (https://www.illumina.com) according to the ARTIC protocol; we generated 2 × 150-bp paired-end reads by using the ARTIC 4.1 primer scheme (2). Wastewater raw reads are available from the European Nucleotide Archive (https://www.ebi.ac.uk/ena; accession no. PRJEB66221). We trimmed reads by using Trim Galore with default settings (3; https://zenodo.org/record/5127899). We removed human sequence reads by using the BWA-MEM alignment algorithm with default settings (H. Li, unpub. data, http://arxiv.org/abs/1303.3997) and the human genome reference build GRCh38. We then used BWA-MEM with default settings to map SARS-CoV-2 reads to the SARS-CoV-2 wild-type reference genome (GenBank accession no. MN908947.3). We performed primer trimming by using iVar with a minimum read length of 30 nt (4) and estimated SARS-CoV-2 lineage abundance in each sample by using Freyja; depth cutoff was 10×, and the lineage abundance filter was 5% (5). We used a 50% coverage minimum across the genome as the threshold for lineage calling. We obtained sequencing results for 13 (42%) of 31 SARS-CoV-2–positive samples (Appendix Table).
We analyzed sequence reads for each sample aliquot and also after combining raw reads for each aircraft (Table; Appendix Table). When reads were combined for each aircraft, we found that the SARS-CoV-2 BQ.1 variant was dominant in wastewater of 1 aircraft, and XBB.1 variants were dominant in wastewater of 5 aircraft; the XBB.1.5 subvariant was dominant in 4 of those 5 aircraft (Table). The discovery of predominant XBB subvariants (dominant in Europe and the United States during the study period) in the aircraft samples contrasts with variant data uploaded to the GISAID database (https://www.gisaid.org) from China within the same time frame, which were mainly subvariants BA.5.2.48 and BF.7.14 (6).
For wastewater-based surveillance, limited information is generally available regarding the persons who contributed to the samples, and, consequently, data related to travel history and place of residence are lacking. Because of the lack of supporting information for passengers, the SARS-CoV-2 variants observed in wastewater-based surveillance of aircraft arriving in Copenhagen might have come from passengers infected outside of China.
In conclusion, our findings indicate that the largely infection-naive population of China might not have comprised a strong selective force driving SARS-CoV-2 toward variants with immune evasive features, such as BQ.1 and XBB. Thus, if BQ.1 and XBB.1.5 subvariants were indeed circulating in China to the extent suggested by our analysis, their dominance in wastewater samples might have occurred because of a founder effect in selected communities instead of those variants arising in China. Since January and February 2023, XBB has become the dominant variant in sequence data from China (6). No new variants have been identified, but our study highlights the potential for wastewater-based surveillance to monitor virus spread among airline passengers in a cost-effective, anonymous, and noninvasive manner and to potentially identify circulating variants. This method can be rapidly modified to include other emerging infectious agents and can contribute substantially to future public health surveillance.
Ms. Qvesel is a bioinformatician at the Virus Surveillance and Research Section, Statens Serum Institut, Denmark. Her research interests focus on viral sequence analysis and molecular evolution.
Acknowledgments
We thank Man-Hung Eric Tang and Jannik Fonager for their insights regarding the analyses and the study in general; Mohammad El-Najjar, Hannibal Morten Schultz, Randi Thøgersen, Cecilie Muss, and Nadia Hamza for their thorough laboratory work and sample sequencing; and the staff at Copenhagen Airport for their cooperation.
National wastewater monitoring of SARS-CoV-2 in Denmark is supported by the Ministry of the Interior and Health of Denmark and by a grant from the European Commission Directorate-General Environment (agreement no. 060701/2021/864407/SUB/ENV.C2).
References
- European Centre for Disease Prevention and Control. EASA/ECDC provide guidelines for aviation as part of European response to COVID-19 developments in China. 2023 [cited 2023 Mar 13]. https://www.ecdc.europa.eu/en/news-events/easaecdc-provide-guidelines-aviation-part-european-response-covid-19-developments-china
- Quick J. nCoV-2019 sequencing protocol V.1. 2020 Jan 22 [cited 2023 Apr 17]. https://www.protocols.io/view/ncov-2019-sequencing-protocol-bbmuik6w
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2. DOIGoogle Scholar
- Grubaugh ND, Gangavarapu K, Quick J, Matteson NL, De Jesus JG, Main BJ, et al. An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar. Genome Biol. 2019;20:8. DOIPubMedGoogle Scholar
- Karthikeyan S, Levy JI, De Hoff P, Humphrey G, Birmingham A, Jepsen K, et al. Wastewater sequencing reveals early cryptic SARS-CoV-2 variant transmission. Nature. 2022;609:101–8. DOIPubMedGoogle Scholar
- Khare S, Gurry C, Freitas L, Schultz MB, Bach G, Diallo A, et al. GISAID’s role in pandemic response. China CDC Wkly. 2021;3:1049–51. DOIPubMedGoogle Scholar
Table
Cite This ArticleOriginal Publication Date: October 26, 2023
1These first authors contributed equally to this article.
2These senior authors contributed equally to this article.
Table of Contents – Volume 29, Number 12—December 2023
EID Search Options |
---|
Advanced Article Search – Search articles by author and/or keyword. |
Articles by Country Search – Search articles by the topic country. |
Article Type Search – Search articles by article type and issue. |
Please use the form below to submit correspondence to the authors or contact them at the following address:
Amanda Gammelby Qvesel, Virus Surveillance and Research Section, Department of Virus and Microbiological Special Diagnostics, Artillerivej 5, DK-2300 Copenhagen S, Denmark
Top