key: cord-0976035-ajt3noh8
authors: Sarwar, Muhammad Bilal; Yasir, Muhammad; Alikhan, Nabil-Fareed; Afzal, Nadeem; de Oliveira Martins, Leonardo; Le Viet, Thanh; Trotter, Alexander J.; Prosolek, Sophie J.; Kay, Gemma L.; Foster-Nyarko, Ebenezer; Rudder, Steven; Baker, David J.; Muntaha, Sidra Tul; Roman, Muhammad; Webber, Mark A.; Shafiq, Almina; Shabbir, Bilquis; Akram, Javed; Page, Andrew J.; Jahan, Shah
title: SARS-CoV-2 variants of concern dominate in Lahore, Pakistan in April 2021
date: 2021-11-30
journal: Microb Genom
DOI: 10.1099/mgen.0.000693
sha: 5853030d36475c44182e421d5c3c73f68086c9fd
doc_id: 976035
cord_uid: ajt3noh8

The SARS-CoV-2 pandemic continues to expand globally, with case numbers rising in many areas of the world, including the Indian sub-continent. Pakistan has one of the world’s largest populations, of over 200 million people and is experiencing a severe third wave of infections caused by SARS-CoV-2 that began in March 2021. In Pakistan, during the third wave until now only 12 SARS-CoV-2 genomes have been collected and among these nine are from Islamabad. This highlights the need for more genome sequencing to allow surveillance of variants in circulation. In fact, more genomes are available among travellers with a travel history from Pakistan, than from within the country itself. We thus aimed to provide a snapshot assessment of circulating lineages in Lahore and surrounding areas with a combined population of 11.1 million. Within a week of April 2021, 102 samples were sequenced. The samples were randomly collected from two hospitals with a diagnostic PCR cutoff value of less than 25 cycles. Analysis of the lineages shows that the Alpha variant of concern (first identified in the UK) dominates, accounting for 97.9 % (97/99) of cases, with the Beta variant of concern (first identified in South Africa) accounting for 2.0 % (2/99) of cases. No other lineages were observed. In depth analysis of the Alpha lineages indicated multiple separate introductions and subsequent establishment within the region. Eight samples were identical to genomes observed in Europe (seven UK, one Switzerland), indicating recent transmission. Genomes of other samples show evidence that these have evolved, indicating sustained transmission over a period of time either within Pakistan or other countries with low-density genome sequencing. Vaccines remain effective against Alpha, however, the low level of Beta against which some vaccines are less effective demonstrates the requirement for continued prospective genomic surveillance.

The COVID-19 pandemic has spread rapidly throughout the world and continues to expand in many regions. It began with OPEN ACCESS an unknown case of pneumonia in the city of Wuhan, PR China [2] . The causative pathogen has since been named 'severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2)' . As of May 2021, there have been over 167 million reported cases and 3.4 million fatalities [3] . Genomic surveillance has assisted the pandemic response providing information for outbreak investigations and detecting possible epitope changes that would allow the virus to escape vaccines. Multiple classification systems have been developed to quickly communicate SARS-CoV-2 variants circulating in a community [4] [5] [6] . From these definitions, certain lineages have been designated as variants of concern (VOCs), which are defined as such due to indications of increased transmission patterns and/or possible resistance to vaccine and/or other treatments [7] . The World Health Organisation (WHO) has introduced a new nomenclature of these VOCs and variants of Interest (VOIs) based on Greek alphabets. B.1.1.7 (Alpha) and B.1.351 (Beta) are two VOCs that have circulated globally. The SARS-CoV-2 lineage B.1.1.7, designated variant of concern 202012/01 (VOC) by Public Health England, was first identified in the UK in late Summer to early Autumn 2020 [7] . B.1.351 is another VOC identified in South Africa and defined by eight mutations in the spike protein including (K417N, E484K and N501Y) [8] .

Pakistan is currently experiencing a severe third wave of infections caused by SARS-CoV-2, which began in March 2021. Vaccination rates nationally are under 2 % [9] , leaving large segments of the community at risk of serious illness from COVID-19. Currently very few SARS-CoV-2 genomes collected in Pakistan are available, with just 12 covering the third wave, nine of which are from one city, Islamabad. This highlights the need for more genome sequencing from Pakistan, particularly given the current situation and the very high population in order to allow surveillance of variants in circulation. Currently, more genomes are available in the GISAID database for travellers with a travel history from Pakistan, than from within the country itself [10] .

We have amplicon sequenced 102 samples, randomly chosen from a 1 week period in April 2021 from Lahore and surrounding areas, to get a snapshot assessment of the circulating lineage in the region. This has identified the Alpha variant of concern as the primary lineage circulating, found in 97.9 % of cases, with clear signals of repeated overseas introductions into the region.

RNA was extracted using Viral RNA Extraction kit, Favor-Prep (Cat. No. FAVNK 001-1) and and detection was carried out on IQ5five real time PCR using GenomeCoV19 Detection Kit by abm (cat: 628) in facility of resource lab for research in biomedical sciences and COVID-19 research lab of University of Health Sciences, Lahore, Pakistan.

Positive samples with a CT <25 were randomly selected for genome sequencing. Viral RNA was converted in cDNA (SuperScript III) and was amplified using the ARTIC protocol v3 (LoCost) [11] with sequencing libraries prepared using CoronaHiT [12] . We carried out genome sequencing using the Illumina NextSeq 500 platform.

The raw reads were demultiplexed using bcl2fastq (v2.20). The reads were used to generate a consensus sequence using the ARTIC bioinformatic pipeline [13] . Briefly, the reads had adapters trimmed with TrimGalore [14] and were aligned to the WuhanHu-1 reference genome (accession MN908947.3) using BWA-MEM (v0.7.17) [15] ; the ARTIC amplicons were trimmed and a consensus built using iVAR (v.1.2.3) [16] .

PANGO lineages assigned using Pangolin v2.4.2 and PangoLEARN model dated 12 May 2021 [6, 17] .

For the phylogenetic analysis, all sequences from GISAID (available on 14 May 2021) where Pakistan is the country of exposure were downloaded and added to the sequences from the current study. From these 268 samples, 131 were sequenced outside Pakistan, with 110 from the Airport Quarantine Station in Japan (travellers from Pakistan sampled at arrival) [18] . All remaining sequences from GISAID were then compared to this data set where we kept the closest ones -for each Pakistan sequence, the four closest from non-Pakistan were kept for subsequent analysis to provide context. The alignment and neighbour search were done with uvaia [19] , and problematic (homoplasic or difficult to sequence) sites were masked from the alignment [20] . For the B.1.1.7 (Alpha) and B.1.351 (Beta) dated phylogenetic inference, we further enriched each alignment with more distant neighbours, which maximized the phylogenetic diversity [21] , based on the neighbour-joining tree [22] of all sequences close to each Pakistan sequence, for each lineage. Sequences with more than 10 % N, or with incomplete date information were excluded from analysis. Clusters distant and unrelated to Pakistan sequences were reduced or removed by visual inspection of maximum-likelihood trees.

From the 102 UHS-PAK sequences, 90 samples were included in the phylogenetic analysis: 88 from Alpha and 2 samples

Genome sequencing of SARS-CoV-2 helps public health scientists to understand the dynamics of the pandemic in different parts of the world. Low-and middle-income countries have low levels of genomic surveillance, depriving them of effective data sources to manage their pandemic response and set public health policy. We provide a snapshot of samples taken over 1 week from Lahore, Pakistan and find that all of the samples sequenced were variants of concern, which have either increased transmissibility or reduced efficacy against vaccines. The 102 samples sequenced represented a fold increase in the number of genomes available for Pakistan (at the time of sequencing).

from Beta; the final alignments have 723 and 107 sequences, respectively. The maximum-likelihood trees were inferred with IQTREE2 v2.1.2 [23] with 1000 ultrafast bootstrap replicates [24] . The divergence times were estimated by marginalization under a relaxed clock using TreeTime [25] and assuming a coalescent prior on three lengths with a constant effective population size. Polytomies were resolved by treetime when possible by using the time information. Ancestral state reconstruction of the country of exposure was done under parsimony over the timed tree with castor [26] and trees were plotted with ggtree [27] , both for R. By weighting reconstructed states (countries of exposure) at the internal nodes by the number of most parsimonious scenarios, castor gives the probability of each node belonging to Pakistan or abroad. The number of introductions into Pakistan was estimated by counting branches where the reconstructed probability of being exposed in Pakistan increased to one. Equivalently, the exportations out of Pakistan were given by branches with a parental node most probably from Pakistan (prob >0.5) and a child node certainly not from Pakistan (prob=0). We assume that the date of the event is given by the parental node, to discount for 'importation' and 'detection' lags [28] , and thus may be biassed since the true transition happened somewhere along the branch (i.e. more recent than our estimate). For this analysis we excluded the 21 Pakistani samples sequenced outside Pakistan, which might have been circulating abroad (i.e. except those from the Japanese Airport Quarantine Station). Notice that if the parsimony assumption fails, e.g. if we expect frequent international transmission events hidden along nodes in the tree, then this should be modelled under a probabilistic (likelihood or Bayesian) approach. Since we included the closest sequences from abroad, the reconstruction should not change by adding more available data. However, due to unequal sampling and sequencing, the number of transitions is an underestimate and their dates are subject to selection bias (e.g. previous introductions were not sequenced due to regional differences or severity of infection).

The distances between all leaf pairs along the maximumlikelihood tree were calculated with the package ape for R. These so-called patristic (or cophenetic) distances describe the path length along all branches separating two samples, and by taking the phylogeny into account are a better descriptor of the evolutionary differences between the sequences. By multiplying the patristic distances by genome length we obtain the expected number of SNP differences between samples, since branch lengths in a maximum-likelihood tree are given in 'expected substitutions per site' .

In this study, from 6 April 2021 to 11 April 2021, 102 SARS-CoV-2 samples in viral transport medium (VTM) were randomly selected from SARS-CoV-2 diagnostic positives by the Mayo Hospital and Sheikh Zayed Hospital in Lahore, Pakistan. All samples had a Ct <25. These public hospitals serve a predominantly middle-income population. The cases were between 21 and 91 years of age with a broad distribution of ages, with 74 % (n=74) male, 26 % (n=26) female and two unknown. Four samples were from cases who died (all aged over 60), while all others recovered. Cases had no history of international travel, with 77 cases explicitly reporting that they and their household did not recently travel overseas (Table S1 ). Viral RNA was amplified using the ARTIC protocol [11] with sequencing libraries prepared using CoronaHiT [12] . Resulting sequenced reads were used to generate consensus sequences with the ARTIC bioinformatics protocol (see Methods). SARS-CoV-2 was first identified in Pakistan on 26 February 2020 [29] . The first Alpha genome was submitted to GISAID on 25 December 2020. There were previously 268 SARS-CoV-2 genomes available on GISAID where Pakistan was listed as the country of exposure with 90 assigned as Alpha or Beta (Tables 1  and S1 ). Most of these samples were associated with international travel from the country and were collected outside of Pakistan (Table S2) . A list of submitting authors of these data can be found in Table S3 . We combined all Alpha and Beta public data with the genomes presented in this study to provide a snapshot of Alpha and Beta dissemination (Table S4) .

By reconstructing the exposure history over the phylogenetic tree, we estimate the number of Alpha importations as 112 using the dated tree, as well as 11 exports out of the country (Fig. 1) . From the histogram at the bottom of Fig. 1 , we observe that the first introduction of Alpha into Pakistan may have happened as early as mid October 2020. This is after the date of emergence of Alpha in the UK, which has been calculated as September 2020 elsewhere [30] . There was an increasing number of independent introductions more recently, such that most Alpha samples can be traced back to a recent ancestor abroad, while a few may have been circulating for longer (or belong to unsampled groups).

Genomes from Pakistan were intermingled with genomes from elsewhere, as shown in Fig. 2 where we calculated the patristic distance of genomes in this study to their closest neighbour and found that they were generally closer to samples from abroad, i.e. samples without a travel history to Pakistan (see also Fig. S5 , with the maximum-ikelihood tree). The number of substitutions to the nearest neighbour also varied (0-8 substitutions) (Fig. 2) . All samples in this study were sourced from cases with no travel history. The virus is detected and sequenced from the source, usually the UK, where 2 % of the population describe their ethnicity as Pakistani [22] and where there are many daily direct flights, and then in some cases sampled in Pakistan immediately or after several months, where the number of substitutions were lower or higher, respectively.

Beta was likely introduced into Pakistan recently. Beta samples were found in four clades in the phylogeny of Beta, suggesting at least four separate introductions into the country (Fig. 3) . The time to the most recent common ancestor of Beta globally was calculated here as during September 2020, whereas the clades containing Beta genomes from Pakistan date no earlier than February 2021 (Fig. 3) . The date of emergence of Beta is before these dates and has been calculated as early August 2020 elsewhere [8] .

It has been demonstrated [31] that vaccines, in particular BNT162b2, are effective against Alpha. Therefore, if this remains the dominant lineage in the region, public health vaccination policy can be implemented accordingly. However, whilst the prevalence of Beta, which has been linked to lower vaccine efficacy, is very low, this reinforces the need for prospective surveillance of SARS-CoV-2 using genome sequencing to inform public health interventions in a continual manner.

The most regular source of genomic sequencing is from travellers from Pakistan being sequenced by their destination countries, and being annotated as such in the public databases (GISAID). Japan is the largest contributor of genomes from travellers originating from Pakistan (111 out of 269). This indirect surveillance is useful but unreliable as travel restrictions and pre-flight testing can bias the results.

This report shows the critical importance of whole-genome sequencing of SARS-CoV-2 to determine the prevalence and changing epidemiology of different variants of the virus. This data is crucial to inform public health decision makers as well as to allow global epidemiology to be understood. Building capacity for sequencing and analysis of genomes in countries with high infection rates will be crucial for the global response to COVID-19.

SARS-cov-2 variants of concern dominate in lahore

Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China

An interactive web-based dashboard to track COVID-19 in real time

Geographical and temporal distribution of SARS-CoV-2 clades in the WHO European Region

Year-letter Genetic Clade Naming for SARS-CoV-2 on Nextstrain. org

A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology

Transmission of SARS-COV-2 lineage b.1.1.7 in England: Insights from linking epidemiological and genetic data. medRxiv

Detection of a SARS-CoV-2 variant of concern in South Africa

Pakistan Cases Details

Global initiative on sharing all influenza data -from vision to reality

NCOV-2019 sequencing protocol v2 v1

CoronaHiT: high-throughput sequencing of SARS-CoV-2 genomes

Artic network

UK: The Babraham Institute

An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar

Assignment of epidemiological lineages in an emerging pandemic using the pangolin tool

COVID-19 genome surveillance at international airport quarantine stations in Japan

UK: Quadram Institute Bioscience

Stability of SARS-CoV-2 phylogenies

Phylogenetic diversity within seconds

Rapid computation of distance estimators from nucleotide and amino acid alignments

IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era

UFBoot2: Improving the Ultrafast Bootstrap Approximation

TreeTime: Maximum-likelihood phylodynamic analysis

Efficient comparative phylogenetics on large trees

ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data

Establishment and lineage dynamics of the SARS-CoV-2 epidemic in the UK

Is Pakistan's response to coronavirus (SARS-CoV-2) adequate to prevent an outbreak

Emergence of SARS-CoV-2 B.1.1.7 Lineage -United States

National Study Group for COVID-19 Vaccination. Effectiveness of the BNT162b2 Covid-19 Vaccine against the B.1.1.7 and B.1.351 Variants

The Quadram Institute authors gratefully acknowledge the support of the Biotechnology and Biological Sciences Research Council (BBSRC); their research was funded by the BBSRC Institute Strategic Programme Microbes in the Food Chain BB/R012504/1 and its constituent project BBS/E/F/000PR10352, also Quadram Institute Bioscience BBSRC funded Core Capability Grant (project number BB/CCG1860/1). The University of Health Sciences authors acknowledge the support provided by the Higher Education Commission (HEC) of Pakistan under the project RRG-211.

The authors declare that there are no conflicts of interest.

This project was conducted under approval number UHS/REG-20/ ERC/1758 from the University of Health Sciences Lahore Ethical Review Committee.

Five reasons to publish your next article with a Microbiology Society journal 1 . The Microbiology Society is a not-for-profit organization. 2. We offer fast and rigorous peer review -average time to first decision is 4-6 weeks. 3. Our journals have a global readership with subscriptions held in research institutions around the world. 4. 80% of our authors rate our submission process as 'excellent' or 'very good'. 5. Your article will be published on an interactive journal platform with advanced metrics.Find out more and submit your article at microbiologyresearch.org.