Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
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

Author affiliation: Author affiliation: Statens Serum Institut, Copenhagen, Denmark

Cite This Article

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.

Top

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).

Top

References

  1. 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
  2. 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
  3. Martin  M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:102. DOIGoogle Scholar
  4. 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
  5. 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:1018. DOIPubMedGoogle Scholar
  6. 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:104951. DOIPubMedGoogle Scholar

Top

Table

Top

Cite This Article

DOI: 10.3201/eid2912.230717

Original 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
presentation_01 Advanced Article Search – Search articles by author and/or keyword.
presentation_01 Articles by Country Search – Search articles by the topic country.
presentation_01 Article Type Search – Search articles by article type and issue.

Top

Comments

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

Send To

10000 character(s) remaining.

Top

Page created: October 11, 2023
Page updated: November 18, 2023
Page reviewed: November 18, 2023
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external