Genetic Diversity of SARS-CoV-2 among Travelers Arriving in Hong Kong

We sequenced 10% of imported severe acute respiratory syndrome coronavirus 2 infections detected in travelers to Hong Kong and revealed the genomic diversity of regions of origin, including lineages not previously reported from those countries. Our results suggest that international or regional travel hubs might be useful surveillance sites to monitor sequence diversity.


DISPATCHES
H ong Kong uses an elimination strategy to control coronavirus disease (COVID-19) that includes stringent travel restrictions to reduce the risk of introducing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into local communities (1). COVID-19 testing was mandated on departure and arrival for all inbound travelers. Compulsory 14-day home quarantine was put in place for all arrivals beginning March 19, 2020. Nonresidents were banned from entry after March 25. In subsequent months, persons arriving from high-risk locations were required to quarantine in hotels; by November, all arrivals had to quarantine in hotels. On December 25, the quarantine period was extended to 21 days. Predeparture COVID-19 testing was mandated for travelers inbound from high-risk locations. Furthermore, daily health declarations were required from all quarantined travelers and respiratory samples were collected on arrival, day 12, and day 19 (for 21-day quarantine) for reverse transcription PCR (RT-PCR) testing. As of April 25, 2021, authorities had recorded 11,731 RT-PCR-positive COVID-19 cases in Hong Kong. About 20% (2,350) of the laboratory-confi rmed COVID-19 cases were considered imported, detected in persons thought to have been infected outside of Hong Kong. Here, we report the analyses of 10% of these imported cases through genome sequencing.
To estimate the viral sequence diversity among these imported cases, we performed next-generation

Genetic Diversity of SARS-CoV-2 among Travelers Arriving in Hong Kong
We sequenced 10% of imported severe acute respiratory syndrome coronavirus 2 infections detected in travelers to Hong Kong and revealed the genomic diversity of regions of origin, including lineages not previously reported from those countries. Our results suggest that international or regional travel hubs might be useful surveillance sites to monitor sequence diversity.
sequencing on 10% (221) of clinical samples collected (2,3) (Appendix). We selected a greater proportion of samples (204) beginning in June 2020 when greater genetic diversity began to appear globally. The number of samples we sequenced by country of origin was proportional to all cases detected in travelers from that country (R = 0.91).
Using the Pangolin classification system  Table  2). We detected 2 variants of concern (VOC) and 3 variants of interest (VOI; Table 1) (5). VOC B.1.1.7 (Alpha variant), which began spreading rapidly in the United Kingdom in November 2020 (6,7), was the most common VOC (39) in our study. We first detected this lineage in a passenger arriving from the United Kingdom on December 13, 2020, and we subsequently detected it in another 38 travelers from other countries, predominantly from the Philippines and Pakistan (Table 1). This finding corresponds with data from global surveillance that indicate this lineage has been circulating over a wide geographic range beginning in December 2020. The second VOC, B.1.351 (Beta), which was first reported to circulate widely in South Africa beginning in November 2020 (8) Fifty percent of our cases were imported from 5 middle-income countries in Asia: India, Indonesia, Nepal, Pakistan, and the Philippines (https://databank.worldbank.org/data/download/site-content/ CLASS.xls; Appendix Table 1). We wanted to compare the genomic diversity of SARS-CoV-2 imported from these countries with those reported in the GI-SAID database (https://www.gisaid.org). However, the Philippines, Nepal, and Pakistan had limited SARS-CoV-2 sequence information in the GISAID database (Table 2) (9). Of the 3 VOC or VOI we identified in travelers from the Philippines (Table 2), B.1.351 was not among sequences the Philippines submitted to GISAID, but the March 6-20, 2021, arrival dates of the 5 case-patients with B.1.351 suggest unreported domestic circulation of that lineage. Similarly, Nepal had reported to GISAID only 15 of the 20 viral sequences from 8 lineages we had identified. Other countries also had not previously reported several lineages we identified to GISAID, including 3 from India and 1 each from Pakistan and Indonesia. We did not analyze samples from travelers from some countries, either because they had their own extensive domestic sequencing efforts or we had few samples from these countries (<5 per country).
We further compared GISAID data with our data from the Philippines, Nepal, and Pakistan. We retrieved the earliest collection date for each lineage we detected that these countries had also reported to GI-SAID; some of those dates were close to the first dates of arrival for case-patients with those lineages in our study. In fact, in over half of those lineages reported in both sources, we identified the lineage either before or <1 month after it was reported by the country (Appendix Table 3), highlighting the potential use of this method of surveillance to assess genomic diversity in regions with limited sequence information.
The emergence of VOC and VOI in different geographic locations highlights the need for global-level genomic surveillance of SARS-CoV-2 (10), but genomic sequencing information from some regions remains incomplete. Our findings suggest that SARS-CoV-2 in Travelers to Hong Kong travel hubs such as Hong Kong can be used as surveillance sites to identify infected travelers from regions with widespread circulation of lineages of interest. Such indirect surveillance might provide useful data to partially reveal virus diversity in countries with limited sequence information, leading to better preparedness for and response to newly emerging SARS-CoV-2 variants. However, findings from these indirect analyses are likely to be only partial and skewed by the level of passenger traffic to destination countries from various points of departure. Also, the extent of different virus lineages circulating in a country of departure may have affected our observations; lineages that circulate at a low level in a country of interest might be missed by our current strategy. Optimizing this approach, such as by directing sequencing efforts toward travelers departing from targeted countries or regions rather than at the points of arrival, might help overcome those limitations.