Detection and genome sequencing of SARS-CoV-2 belonging to the B.1.1.7 lineage in the Philippines ================================================================================================= * Francis A. Tablizo * Cynthia P. Saloma * Marc Jerrone R. Castro * Kenneth M. Kim * Maria Sofia L. Yangzon * Carlo M. Lapid * Benedict A. Maralit * Marc Edsel C. Ayes * Jan Michael C. Yap * Jo-Hannah S. Llames * Shiela Mae M. Araiza * Kris P. Punayan * Irish Coleen A. Asin * Candice Francheska B. Tambaoan * Asia Louisa U. Chong * Karol Sophia Agape R. Padilla * Rianna Patricia S. Cruz * El King D. Morado * Joshua Gregor A. Dizon * Eva Maria Cutiongco-de la Paz * Alethea R. de Guzman * Razel Nikka M. Hao * Arianne A. Zamora * Devon Ray Pacial * Juan Antonio R. Magalang * Marissa Alejandria * Celia Carlos * Anna Ong-Lim * Edsel Maurice Salvaña * John Q. Wong * Jaime C. Montoya * Maria Rosario Singh-Vergeire ## ABSTRACT We report the sequencing and detection of 39 Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) samples containing lineage-defining mutations specific to viruses belonging to the B.1.1.7 lineage (UK variant) in the Philippines. ## ANNOUNCEMENT The Coronavirus Disease 2019 (COVID-19) is an infectious disease that has gained pandemic status from the World Health Organization, with millions of cases and deaths recorded worldwide. This global health crisis is caused by viruses referred to as Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), a member of the genus *Betacoronavirus* (Coronaviridae) together with the causative agents of the first SARS outbreak in 2003 and the Middle East Respiratory Syndrome (MERS) in 2012. In this study, we present the genome sequences of 39 cases of COVID-19 in the Philippines caused by viruses belonging to SARS-CoV-2 lineage B.1.1.7, more commonly known as the United Kingdom (UK) variant. This particular SARS-CoV-2 variant was initially identified in the United Kingdom and has been reported to cause a surge of COVID-19 infections in the said country (Kirby, 2021). Initial studies also suggest that the B.1.1.7 viruses appear to have replicative advantage (Grabowski et al., 2021) and are more transmissible (Leung et al., 2021). A total of 39 Philippine SARS-CoV-2 samples were classified under the B.1.1.7 lineage. Table 1 shows the primary consensus assembly metrics for these samples. The average depth of coverage across all the sequences is 1,128x, with 33 of the samples carrying all 17 hallmark mutations of the B.1.1.7 lineage as listed in the PANGO Lineages report for the B.1.1.7 variant of concern ([https://cov-lineages.org/global\_report_B.1.1.7.html](https://cov-lineages.org/global_report_B.1.1.7.html)). View this table: [Table 1.](http://medrxiv.org/content/early/2021/03/08/2021.03.04.21252557.1/T1) Table 1. Primary consensus sequence assembly metrics. The detection of B.1.1.7 from returning overseas Filipino workers and in the community highlights the need for genomic surveillance at the country’s ports of entry and the general population to monitor the importation and local transmission of emerging variants of concern that may impact the public health response to the SARS-CoV-2 pandemic in the Philippines. ## METHODS In order to detect the entry of B.1.1.7 infection in the Philippines, nasopharyngeal swabs were collected between December 10, 2020 and January 31, 2021 from COVID-19 cases detected in returning overseas Filipinos, as well as from local case clusters mainly from the Cordillera Administrative Region of the country, among others. The collected swab samples were subjected to RNA extraction using the QIAamp Viral RNA Mini Kit, the product of which were then used as templates for the amplicon-based Illumina COVIDSeq Test ([https://www.illumina.com/products/by-type/ivd-products/covidseq.html](https://www.illumina.com/products/by-type/ivd-products/covidseq.html)) sequencing workflow. Afterwards, the resulting sequence reads were mapped to the reference SARS-CoV-2 genome (NCBI Accession Number NC 045512.2) using Minimap2 (Li, 2018) version 2.17-r941, with the \`–x sr\` parameter for short accurate genomic reads alignment. Primer clipping and quality trimming, intrahost variant calling, removal of reads associated with mismatched primer indices, and consensus sequence assembly were all done following the suggested workflow of iVar (Grubaugh et al., 2019) version 1.2.2, using default parameters. The consensus variants were identified by comparing the resulting assemblies with the reference sequence using MUMmer (Kurtz et al., 2004) as implemented in the RATT software (Otto et al., 2011). Lastly, SARS-CoV-2 lineage classifications were assigned using Pangolin version 2.3.2 (github.com/cov-lineages/pangolin; Rambaut et al., 2020). ## PUBLIC AVAILABILITY AND ACCESSION NUMBERS The consensus sequence assemblies reported in this study are deposited in the EpiCoV™ database of the GISAID. The accession codes are shown in **Table 1**. ## Data Availability The consensus sequence assemblies reported in this study are deposited in the EpiCoV database of the GISAID. The accession codes are shown in Table 1. ## ACKNOWLEDGEMENT This project was supported by a Genomics Biosurveillance grant from the Philippine Department of Health, and a Department of Science and Technology – Philippine Council for Health Research and Development grant to BM and the University of the Philippines. The authors would also like to thank the contributing institutions to the Philippine Genomic Biosurveillance Network. The protocols used in this study were reviewed and approved by the Single Joint Research Ethics Board of the Department of Health, with approval code SJREB-2021-11, as part of a larger research program entitled “A retrospective study on the national genomic surveillance of COVID-19 transmission in the Philippines by SARS-CoV-2 genome sequencing and bioinformatics analysis”. ## Footnotes * We received a review of the paper and made minor revisions in the title, abstract, and announcement parts of the paper. We also removed Figure 1 as suggested by a reviewer and simply retained Table 1. * Received March 4, 2021. * Revision received March 8, 2021. * Accepted March 8, 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. Grabowski F., Preibisch G., Kochanczyk M., and Lipniacki T. (2021). SARS-CoV-2 Variant Under Investigation 202012/01 has more than twofold replicative advantage. medRXiV. DOI: [https://doi.org/10.1101/2020.12.28.20248906](https://doi.org/10.1101/2020.12.28.20248906). 2. Grubaugh N., Gangavarapu K., Quick J., et al. (2019). An amplicon-based sequencing framework for accurately measuring intrahost virus diversity using PrimalSeq and iVar. Genome Biology. 20(8). DOI:[https://doi.org/10.1186/s13059-018-1618-7](https://doi.org/10.1186/s13059-018-1618-7). 3. Kirby T. (2021). New variant of SARS-COV-2 in UK causes surge of COVID-19. Lancet Respiratory Medicine. DOI:[https://doi.org/10.1016/S2213-2600(21)00005-9](https://doi.org/10.1016/S2213-2600(21)00005-9). 4. Kurtz S., Phillippy A., Delcher A., Smoot M., Shumway M., Antonescu C., and Salzberg S. (2004). Versatile and open software for comparing large genomes. Genome Biology. 5:R12. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/gb-2004-5-2-r12&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14759262&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F08%2F2021.03.04.21252557.1.atom) 5. Leung K., Shum M., Leung G., Lam T., and Wu J. (2021). Early transmissibility assessment of the N501Y mutant strains of SARS-CoV-2 in the United Kingdom, October to November 2020. EuroSurveillance. 26(1):pii=2002106. DOI:[https://doi.org/10.2807/1560-7917.ES.2020.26.1.2002106](https://doi.org/10.2807/1560-7917.ES.2020.26.1.2002106). 6. Li H. (2018). Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 34:3094–3100. doi:10.1093/bioinformatics/bty191. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bty191&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29750242&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F08%2F2021.03.04.21252557.1.atom) 7. Otto T., Dillon G., Degrave W., and Berriman M. (2011). RATT: rapid annotation transfer tool. Nucleic Acids Research. 39(9):e57. DOI:[https://doi.org/10.1093/nar/gkq1268](https://doi.org/10.1093/nar/gkq1268). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkq1268&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21306991&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F08%2F2021.03.04.21252557.1.atom) 8. Rambaut, A., Holmes, E.C., O’Toole, Á. et al. (2020). A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nature Microbiology. 5:1403–140. DOI:[https://doi.org/10.1038/s41564-020-0770-5](https://doi.org/10.1038/s41564-020-0770-5). 9. Shu, Y., and McCauley, J. (2017). GISAID: Global initiative on sharing all influenza data – from vision to reality. EuroSurveillance. 22(13). DOI:10.2807/1560-7917.ES.2017.22.13.30494. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2807/1560-7917.ES.2017.22.13.30494&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28382917&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F03%2F08%2F2021.03.04.21252557.1.atom)