Ogün ADEBALİ, Aylin BİRCAN, Defne ÇİRCİ, Burak İŞLEK, Zeynep KILINÇ, Berkay SELÇUK, Berk TURHAN

Phylogenetic analysis of SARS-CoV-2 genomes in Turkey

COVID-19 has effectively spread worldwide. As of May 2020, Turkey is among the top ten countries with the most cases. A comprehensive genomic characterization of the virus isolates in Turkey is yet to be carried out. Here, we built a phylogenetic tree with globally obtained 15,277 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes. We identified the subtypes based on the phylogenetic clustering in comparison with the previously annotated classifications. We performed a phylogenetic analysis of the first 30 SARS-CoV-2 genomes isolated and sequenced in Turkey. We suggest that the first introduction of the virus to the country is earlier than the first reported case of infection. Virus genomes isolated from Turkey are dispersed among most types in the phylogenetic tree. We find 2 of the seventeen subclusters enriched with the isolates of Turkey, which likely have spread expansively in the country. Finally, we traced virus genomes based on their phylogenetic placements. This analysis suggested multiple independent international introductions of the virus and revealed a hub for the inland transmission. We released a web application to track the global and interprovincial virus spread of the isolates from Turkey in comparison to thousands of genomes worldwide.

PDF

___

Balaban M, Moshiri N, Mai U, Jia X, Mirarab S (2019). TreeCluster: Clustering biological sequences using phylogenetic trees. PLoS One 14: e0221068.
Bhattacharyya C, Das C, Ghosh A, Singh AK, Mukherjee S et al. (2020). Global spread of SARS-CoV-2 subtype with spike protein mutation D614G is shaped by human genomic variations that regulate expression of TMPRSS2 and MX1 genes. BioRxiv. doi: 10.1101/2020.05.04.075911
Dos Ramos F, Carrasco M, Doyle T, Brierley I (2004). Programmed -1 ribosomal frameshifting in the SARS coronavirus. Biochemical Society Transactions 32 (6): 1081-1083. doi: 10.1042/BST0321081
Hadfield J, Megill C, Bell SM, Huddleston J, Potter B et al. (2018). Nextstrain: real-time tracking of pathogen evolution. Bioinformatics 34 (23): 4121-4123. doi: 10.1093/ bioinformatics/bty407
Huerta-Cepas J, Serra F, Bork P (2016). ETE 3: reconstruction, analysis, and visualization of phylogenomic data. Molecular Biology and Evolution 33 (6): 1635-1638. doi: 10.1093/molbev/ msw046
Jaimes JA, Andre NM, Chappie JS, Millet JK, Whittaker GR (2020). Phylogenetic analysis and structural modeling of SARSCoV-2 spike protein reveals an evolutionary distinct and proteolytically sensitive activation loop. Journal of Molecular Biology 432 (10): 3309-3325. doi: 10.1016/j.jmb.2020.04.009
Katoh K, Standley DM (2016). A simple method to control overalignment in the MAFFT multiple sequence alignment program. Bioinformatics 32 (13): 1933-1942. doi: 10.1093/ bioinformatics/btw108
Kelly JA, Dinman JD (2020). Structural and functional conservation of the programmed -1 ribosomal frameshift signal of SARSCoV-2. BioRxiv. doi: 10.1101/2020.03.13.991083
Li C, Yang Y, Ren L (2020). Genetic evolution analysis of 2019 novel coronavirus and coronavirus from other species. Infection, Genetics and Evolution 82: 104285. doi: 10.1016/j. meegid.2020.104285
Li Q, Guan X, Wu P, Wang X, Zhou L et al. (2020). Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia. The New England Journal of Medicine 382: 1199- 1207. doi: 10.1056/NEJMoa2001316
Li X, Giorgi EE, Marichann MH, Foley B, Xiao C et al. (2020). Emergence of SARS-CoV-2 through recombination and strong purifying selection. BioRxiv. doi: 10.1101/2020.03.20.000885
Morais Júnior IJ, Polveiro RC, Souza GM, Bortolin DI, Sassaki FT et al. (2020). The global population of SARS-CoV-2 is composed of six major subtypes. BioRxiv. doi: 10.1101/2020.04.14.040782
Rambaut A, Holmes EC, Hill V, O’Toole Á, McCrone J et al. (2020). A dynamic nomenclature proposal for SARSCoV-2 to assist genomic epidemiology. BioRxiv. doi: 10.1101/2020.04.17.046086
Sagulenko P, Puller V, Neher RA (2018). TreeTime: maximumlikelihood phylodynamic analysis. Virus Evoluation 4 (1): vex042. doi: 10.1093/ve/vex042
Shu Y, McCauley J (2017). GISAID: Global initiative on sharing all influenza data - from vision to reality. Euro Surveill 22 (13): 30494. doi: 10.2807/1560-7917.ES.2017.22.13.30494
Strimmer K, Pybus OG (2001). Exploring the demographic history of DNA sequences using the generalized skyline plot. Molecular Biology and Evoluation 18 (12): 2298-2305. doi: 10.1093/ oxfordjournals.molbev.a003776
Tang X, Wu C, Li X, Song Y, Yao X et al. (2020). On the origin and continuing evolution of SARS-CoV-2. National Science Review. doi: 10.1093/nsr/nwaa036
Van Dorp L, Richard D, Tan CC, Shaw LP, Acman M et al. (2020). No evidence for increased transmissibility from recurrent mutations in SARS-CoV-2. BioRxiv. doi: 10.1101/2020.05.21.108506
Yu G, Smith DK, Zhu H, Guan Y, Lam TT-Y (2017). GGTREE: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods in Ecology and Evolution 8: 28-36. doi: 10.1111/2041-210X.12628
Zhao Z, Sokhansanj BA, Rosen GL (2020). Characterizing geographical and temporal dynamics of novel coronavirus SARS-CoV-2 using informative subtype markers. BioRxiv. doi: 10.1101/2020.04.07.030759