Air travel-related outbreak of multiple SARS-CoV-2 variants =========================================================== * Vijaykrishna Dhanasekaran * Kimberly M Edwards * Ruopeng Xie * Haogao Gu * Dillon C Adam * Lydia DJ Chang * Sammi SY Cheuk * Shreya Gurung * Pavithra Krishnan * Daisy YM Ng * Gigi YZ Liu * Carrie KC Wan * Samuel SM Cheng * Dominic NC Tsang * Benjamin Cowling * Malik Peiris * Leo LM Poon ## ABSTRACT **Background** A large cluster of 59 cases were linked to a single flight with 146 passengers from New Delhi to Hong Kong in April 2021. This outbreak coincided with early reports of exponential pandemic growth in New Delhi, which reached a peak of >400,000 newly confirmed cases on 7 May 2021. **Methods** Epidemiological information including date of symptom onset, date of positive-sample detection, and travel and contact history for individual cases from this flight were collected. Whole genome sequencing was performed, and sequences were classified based on the dynamic Pango nomenclature system. Maximum-likelihood phylogenetic analysis compared sequences from this flight alongside other cases imported from India to Hong Kong on 26 flights between June 2020 and April 2021, as well as sequences from India or associated with India-related travel from February to April 2021, and 1,217 reference sequences. **Results** Sequence analysis identified six lineages of SARS-CoV-2 belonging to two variants of concern (Alpha and Delta) and one variant of public health interest (Kappa) involved in this outbreak. Phylogenetic analysis confirmed at least three independent sub-lineages of Alpha with limited onward transmission, a superspreading event comprising 37 cases of Kappa, and transmission of Delta to only one passenger. Additional analysis of another 26 flights from India to Hong Kong confirmed widespread circulation of all three variants in India since early March 2021. **Conclusions** The broad spectrum of disease severity and long incubation period of SARS-CoV-2 pose a challenge for surveillance and control. As illustrated by this particular outbreak, opportunistic infections of SARS-CoV-2 can occur irrespective of variant lineage, and requiring a nucleic acid test within 72 hours of departure may be insufficient to prevent importation or in-flight transmission. KEYWORDS * Genomic epidemiology * air travel-related outbreak * whole genome sequencing * SARS-CoV-2 ## Introduction Since March 2020, increasingly stringent travel measures were implemented in Hong Kong to reduce the risk of SARS-CoV-2 importation1, 2, including the testing pre-departure, on-arrival, and post-arrival, and quarantine of all passenger arrivals outside the home3. Detection of cases on arrival or during quarantine results in contact tracing and additional testing of contacts. In recognition of the emergence of SARS-CoV-2 variants of concern (VOC) with public health significance4, on-arrival quarantine was extended from 14 to 21 days on 25 December 20205. Travellers departing from high-risk countries are required to produce a negative nucleic acid test from an accredited laboratory within 72 hours of departure, and passenger flights from countries with widespread VOC circulation are classified as extremely high risk and temporarily banned. Multiple cases of SARS-CoV-2 were detected upon arrival of a flight from New Delhi to Hong Kong (DEL-HKG). The arrival of this flight in early April 2021 coincided with early reports of exponential growth of SARS-CoV-2 in India6 and immediately prior to widespread reports of healthcare infrastructure failure. Because India was at that time classified by the Hong Kong government as high-risk, passengers were required to show documentation of a 21-day quarantine hotel reservation and proof of a negative nucleic acid test taken within 72 hours of departure. After six passengers were denied boarding, the six-hour flight was at 85% capacity. Upon arrival, passengers were again tested by RT-PCR, at which point five positive cases were detected. Cases among passengers continued to be detected throughout the duration of the 21-day quarantine period. For passengers sharing accommodations with positive cases, quarantine periods were extended to 30 days post-arrival. Representing the largest known air travel-related outbreak of SARS-CoV-27, a total of 59 confirmed cases were detected from the 146 passengers onboard, though only 20% (12/59) of cases were reportedly symptomatic. Through whole genome sequencing and phylogenetic analysis, we aimed to infer the dynamics of SARS-CoV-2 transmission associated with this air travel-related related outbreak. ## Methods Epidemiological information including date of symptom onset, date of positive-sample detection, and travel and contact history for individual cases are provided in Table 1. Of the 59 RT-PCR confirmed SARS-CoV-2 cases, samples from 48 individuals contained adequate viral load for sequencing using previously described methods8, 9. Of those, 46 sequences were of sufficient coverage to be classified based on the dynamic Pango nomenclature system ([https://pangolin.cog-uk.io/](https://pangolin.cog-uk.io/), accessed 10-May-2021)10, however those with <70% genome coverage were excluded from subsequent phylogenetic analysis, resulting in 43 complete or nearly complete SARS-CoV-2 genomes. We performed phylogenetic analysis to compare sequences from this outbreak (Table 1) alongside cases imported on 26 additional flights from India to Hong Kong between July 2020 and May 2021 (n=51) (Table 2), as well as GISAID sequences from India or associated with India-related travel from February to April 2021 (n=2,660), and global reference sequences (n=1,217) including Wuhan-Hu-1 (MN908947.3) (accessed 20-May-2021) (Table S1). After removing outliers based on a root-to-tip regression analysis in TempEst v.1.5.311, a maximum likelihood (ML) tree was generated for the final dataset (n=3,971) in IQ-TREE v.212, 13 and dated using the least square dating (LSD2) method14. View this table: [Table 1.](http://medrxiv.org/content/early/2021/07/23/2021.07.22.21260854/T1) Table 1. SARS-CoV-2 cases detected from a New Delhi to Hong Kong flight in April 2021. View this table: [Table 2.](http://medrxiv.org/content/early/2021/07/23/2021.07.22.21260854/T2) Table 2. Sequenced SARS-CoV-2 genomes from cases imported to Hong Kong on 26 flights from July 2020 to May 2021. ## Results and Discussion Whole genome sequencing verified the presence of three distinct variants onboard the DEL-HKG flight, Alpha, Kappa, and Delta4, covering three Pango lineages10 B.1.1.7, B.1.617.1, and B.1.617.2, respectively. Kappa (B.1.617.1) was detected in the majority of passenger sequences (n=39/46, 84.8%), followed by Alpha (B.1.1.7) in five passengers including a crew member who tested positive on arrival (n=5/46, 10.9%) and Delta (B.1.617.2) in two passengers (n=2/46, 4.3%). In a maximum likelihood phylogeny, 37 of the Kappa sequences clustered together with zero or near-zero branch lengths suggesting a single transmission cluster (inset in Fig. 1A), however one sequence (case 3) was genetically distinct enough to infer pre-flight infection without further onboard transmission. High genome coverage (>98%) of three of the five Alpha variant sequences are consistent with three separate introductions by (i) an infected crew member, (ii) an infected infant, and (iii) an infected but asymptomatic adult in economy seating. Any of these three could represent the index case for onward transmission of the Alpha variant to two individuals from which only partial sequence could be obtained (identified as travel group O in Table 1), testing positive on days 14 and 15 post-arrival. The two Delta variant cases have nearly identical genomes (with >95% coverage), and the two individuals sat in close proximity but tested positive 12 days apart, suggestive of a single onboard transmission event (Table 2 and Fig. 1B). However, there were additional nearby cases which could not be sequenced due to insufficient viral load, and those cases may represent missing links in the inferred transmission chains. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/07/23/2021.07.22.21260854/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/07/23/2021.07.22.21260854/F1) Figure 1. SARS-CoV-2 transmission associated with a flight from New Delhi to Hong Kong in April 2021. Colors represent WHO VOC/VOI and Pango lineage designations10. (A) Maximum likelihood phylogeny of 43 genomes sequenced from 59 cases. Wuhan-Hu-1 (MN908947.3, black tip) was used to root the tree. Diamond tip shapes indicate index cases and circles represent secondary-transmission cases; tips are labelled with case numbers and travel groups. (*) denotes sequences with lower coverage (∼73%). Nodes are labelled with observed amino acid mutations. The Kappa (B.1.617.1) transmission cluster is magnified to show nucleotide substitutions among samples. (B) Case detection colored by SARS-CoV-2 lineage in the 30 days post-arrival. Please contact the corresponding authors to request access to data on seat assignments of passengers testing positive. Using the temporal window from arrival-to-detection as a proxy for SARS-CoV-2 incubation, we estimate that seven individual cases were likely infected prior to travel (Fig. 1A and Table 1). Forty one cases were detected between days 5 and 14 post arrival, suggesting they were likely infected in transit. Of the cases infected in transit that were sequenced, 94.1% (n=32/34) were related to a single Kappa variant transmission cluster, demonstrating potential for superspreading of SARS-CoV-2 on congested flights and reinforcing the importance of stringent travel-related testing. Another 11 cases were detected after the 14th day post arrival, of which only five could be sequenced. At least two of these Kappa variant cases were likely infected by co-quarantined family members (cases 51 and 53, detected on days 16 and 19, respectively, Fig. 1A and Table 1). Case 59, detected 30 days post-arrival, may be an example of an unusually long incubation period as this person travelled alone and was the only known case from this flight in their quarantine hotel (Table 1). However, this was one of 13 cases not sequenced, and quarantine hotel information was missing for 15 cases in total (n=15/59, 25.4%). Given that instances of SARS-CoV-2 transmission among epidemiologically unlinked cases during quarantine have been documented15, we cannot conclusively determine whether transmission occurred in transit or quarantine. The age distribution of positive cases ranged from four months to 63 years (median age = 33 years), and notably, eight of the positive cases were detected in children under the age of two, who are generally exempt from masking requirements. Interestingly, members of at least two families were infected with independent transmission lineages, and sequencing revealed two instances of children infected with different variants than their family members (Travel groups A and B in Fig. 1A and Table 1). In a family of three individuals (Travel group A), one parent tested positive for Alpha, and the infant tested positive for Kappa (the other parent tested positive, but no sequence was recovered from this patient). As before, using arrival-to-detection as a proxy for incubation, it is likely that Travel group A was infected several days prior to departure. In a second family group (n=6, Travel group B), one child was infected with Alpha, while the other five family members, two children and three adults, were infected with Kappa. The child who tested positive for the Alpha variant at 5 days post arrival was likely infected prior to departure, while the other five family members tested positive on days 7-8 post arrival, suggesting they were most likely infected in-transit. As a major travel hub, Hong Kong’s practice of testing inbound passengers has proved a valuable source for sentinel SARS-CoV-2 surveillance and containment. Between June 2020 and April 2021, 51 cases of SARS-CoV-2 were detected on 26 other passenger flights from India to Hong Kong (Table 2)16. Among those flights, the second largest air-travel related SARS-CoV-2 cluster (n=8 cases) was detected just five days after the DEL-HKG outbreak on another flight arriving from New Delhi (Table 2 and Supplementary Fig. S1). All imported cases detected since mid-March 2021 belonged to Alpha (n=8), Delta (n=13), or Kappa (n=2) lineages and originated from Kolkata, New Delhi, and Mumbai, reflecting the predominant circulation of these variants across major Indian cities during the 2021 resurgence. Similar to the DEL-HKG outbreak, a majority (n=33/51, 64.7%) of imported cases sequenced from the other 26 flights were asymptomatic. ## Conclusions The broad spectrum of disease severity and potentially long incubation period of SARS-CoV-2 pose a challenge for surveillance and control. Based on the mean incubation period of 5.7 days for SARS-CoV-217 and that infectiousness peaks prior to symptom onset18, requiring a nucleic acid test within 72 hours of departure is clearly insufficient to eliminate the risk of SARS-CoV-2 importation or in-flight transmission. Thus, these results highlight the value of quarantine and testing on arrival to prevent virus introduction into the community. A recent review of in-flight transmission of SARS-CoV-2 shows secondary attack rates are highly variable7, though, admittedly so are adherence to and enforcement of mask mandates and distancing policies19. Furthermore, as illustrated by this particular outbreak, opportunistic infections of SARS-CoV-2 can occur irrespective of variant lineage. ## Supporting information Figure S1 [[supplements/260854_file02.pdf]](pending:yes) Figure S2 [[supplements/260854_file03.pdf]](pending:yes) Table S1 [[supplements/260854_file04.pdf]](pending:yes) ## Data Availability Hong Kong SARS-CoV-2 genome sequences and associated metadata is deposited at gisaid.org. All anonymized data and sequence accession numbers are available within the manuscript. ## Data Availability Hong Kong SARS-CoV-2 genome sequences and associated metadata is deposited at gisaid.org. All anonymized data and sequence accession numbers are available within the manuscript. ## Data Availability Hong Kong SARS-CoV-2 genome sequences and associated metadata is deposited at gisaid.org. All anonymized data and sequence accession numbers are available within the manuscript. ## Author Contributions V.D. and L.L.M.P. designed the study and were responsible for project supervision. K.M.E, R.X., H.G., D.C.A. and D.N.C.T. collected and analyzed the data. L.D.J.C., S.S.Y.C., S.G., P.K., D.Y.M.N., G.Y.Z.L., C.K.C.W., and S.S.M.C performed genome sequencing. V.D., K.M.E and D.C.A. wrote the original draft, and all authors reviewed and edited the manuscript. ## Funding This work was supported by the Health and Medical Research Fund, Food and Health Bureau of the Hong Kong SAR Government [grant number COVID190205]; and the National Institute of Allergy and Infectious Diseases, National Institutes of Health, Department of Health and Human Services, under contract numbers U01AI151810 and 75N93021C00016. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and writing of the manuscript. ## Conflict of interests The authors declare no conflict of interests. ## Supplementary Data **Figure S1**. Time-scaled phylogeny of global SARS-CoV-2 genomes related to a flight from New Delhi to Hong Kong (DEL-HKG) in April 2021. Tips are coloured by DEL-HKG flight (n=43, red); global Pango reference sequences (n=1217, purple); cases imported from India to Hong Kong on 26 flights between July 2020 and May 2021 (n=51, green); and GISAID sequences from cases in India or associated with India-related travel (n=2660, blue). **Figure S2**. Maximum-likelihood phylogeny of the tree in Figure S1 with branch lengths measured in nucleotide substitutions. **Table S1**. Acknowledgements table. ## Acknowledgements We gratefully acknowledge the staff from the originating laboratories responsible for obtaining the specimens and from the submitting laboratories where the genome data were generated and shared via GISAID (Supplementary Table S1). We acknowledge the technical support provided by colleagues from the Centre for PanorOmic Sciences of the University of Hong Kong. We also acknowledge the Centre for Health Protection of the Department of Health for providing epidemiological data for the study. * Received July 22, 2021. * Revision received July 22, 2021. * Accepted July 23, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. BJ Cowling, S. Ali, TWY Ng, et al. Impact assessment of non-pharmaceutical interventions against coronavirus disease 2019 and influenza in Hong kong: An observational study. Lancet Public Health 2020; 5(5):e279–e288. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 2. DC Adam, P Wu, JY Wong, et al. Clustering and superspreading potential of SARS-CoV-2 infections in Hong Kong. Nat Med 2020; 26(11):1714–1719. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 3. GM Leung, BJ Cowling, JT Wu. From a sprint to a marathon in Hong Kong. N Engl J Med 2020; 382(18):e45. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2009790&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32294373&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 4. PR Krause, TR Fleming, IM Longini, et al. SARS-CoV-2 variants and vaccines. N Engl J Med 2021. 5. The Government of the Hong Kong SAR Press Releases. Government tightens restrictions and compulsory quarantine requirements for persons arriving at Hong kong who have stayed in places outside China. [www.info.gov.hk](http://www.info.gov.hk): xGovHK, 2020.12.25. 6. S Kumar. Second wave of COVID-19: Emergency situation in India. J Travel Med 2021. 7. DO Freedman, A Wilder-Smith. In-flight transmission of SARS-CoV-2: A review of the attack rates and available data on the efficacy of face masks. J Travel Med 2020; 27(8). 8. THC Sit, CJ Brackman, SM Ip, et al. Infection of dogs with SARS-CoV-2. Nature 2020; 586(7831):776–778. 9. H Gu, R Xie, DC Adam, et al. SARS-CoV-2 under an elimination strategy in Hong kong. medRxiv 2021:2021.06.19.21259169. 10. A Rambaut, EC Holmes, A O’Toole, et al. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat Microbiol 2020; 5(11):1403–1407. 11. A Rambaut, TT Lam, L Max Carvalho, OG Pybus. Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen). Virus Evol 2016; 2(1):vew007. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ve/vew007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27774300&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 12. LT Nguyen, HA Schmidt, A von Haeseler, BQ Minh. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 2015; 32(1):268–74. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/molbev/msu300&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25371430&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 13. BQ Minh, HA Schmidt, O Chernomor, et al. IQ-TREE 2:New models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol 2020; 37(5):1530–1534. 14. TH To, M Jung, S Lycett, O Gascuel. Fast dating using least-squares criteria and algorithms. Syst Biol 2016; 65(1):82–97. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/sysbio/syv068&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26424727&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 15. N Eichler, C Thornley, T Swadi, et al. Transmission of severe acute respiratory syndrome coronavirus 2 during border quarantine and air travel, New Zealand (Aotearoa). Emerg Infect Dis 2021; 27(5):1274–1278. 16. H Gu., Dkw Chu., Ldj Chang., et al. Genetic diversity of SARS-CoV-2 in travelers arriving in Hong Kong. Emerging Infectious Diseases 2021; 27(10). 17. H Xin, JY Wong, C Murphy, et al. The incubation period distribution of coronavirus disease 2019 (COVID-19): A systematic review and meta-analysis. Clin Infect Dis 2021. 18. X He, EHY Lau, P Wu, et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med 2020; 26(5):672–675. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-3012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F07%2F23%2F2021.07.22.21260854.atom) 19. Government of India Office of the Director General of Civil Aviation. Strict compliance of COVID-19 protocols during air travel. [https://bit.ly/3wGt6I2](https://bit.ly/3wGt6I2): Government of India, 2021.03.13.