|Year : 2020 | Volume
| Issue : 5 | Page : 474-478
Mapping the genomic landscape & diversity of COVID-19 based on >3950 clinical isolates of SARS-CoV-2: Likely origin & transmission dynamics of isolates sequenced in India
Hina Singh1, Jasdeep Singh1, Mohd Khubaib1, Salma Jamal1, Javaid Ahmed Sheikh2, Sunil Kohli3, Seyed Ehtesham Hasnain4, Syed Asad Rahman5
1 Institute of Molecular Medicine, New Delhi 110 062, India
2 Department of Biotechnology, School of Chemical and Life Sciences, New Delhi 110 062, India
3 Hamdard Institute of Medical Sciences & Research, Jamia Hamdard, New Delhi 110 062, India
4 Institute of Molecular Medicine, New Delhi 110 062; 4Dr. Reddy's Institute of Life Sciences, University of Hyderabad Campus, Hyderabad, Telangana 500 046, India
5 Envirozyme Biotech Pvt Ltd., Hyderabad, Telangana 500 076, India; BioInception Pvt. Ltd, Chelmsford, Essex CM1 1GU, United Kingdom
|Date of Web Publication||20-Jun-2020|
Syed Asad Rahman
Envirozyme Biotech Pvt Ltd., Hyderabad, Telangana 500 076; BioInception Pvt. Ltd, Chelmsford, Essex CM1 1GU
Source of Support: None, Conflict of Interest: None
|How to cite this article:|
Singh H, Singh J, Khubaib M, Jamal S, Sheikh JA, Kohli S, Hasnain SE, Rahman SA. Mapping the genomic landscape & diversity of COVID-19 based on >3950 clinical isolates of SARS-CoV-2: Likely origin & transmission dynamics of isolates sequenced in India. Indian J Med Res 2020;151:474-8
|How to cite this URL:|
Singh H, Singh J, Khubaib M, Jamal S, Sheikh JA, Kohli S, Hasnain SE, Rahman SA. Mapping the genomic landscape & diversity of COVID-19 based on >3950 clinical isolates of SARS-CoV-2: Likely origin & transmission dynamics of isolates sequenced in India. Indian J Med Res [serial online] 2020 [cited 2020 Jul 11];151:474-8. Available from: http://www.ijmr.org.in/text.asp?2020/151/5/474/284485
Hina Singh, Jasdeep Singh & Syed Asad Rahman contributed equally.
The COVID-19 pandemic has stalled the world and catapulted the global health systems into unprecedented chaos. More than 200 countries have been affected by this pandemic, resulting in 2.54 million cases in a short period of time and >0.17 million deaths (as of April 23, 2020), with a mere 0.7 million recoveries. The movement of COVID-19 hotspot from China to Europe, and now to the USA, has been partly due to the staggered restrictions in global travel and partly due to potent transmission through asymptomatic carriers.
India, with 21,393 cases and 681 deaths (as of April 23, 2020), had the lowest figures for any country of the comparable population (0.5 deaths per million population). International travellers or their close contacts formed the majority of initially reported cases. The delayed onset of COVID-19 in India has given it an edge, which allowed it to impose severe restrictions to contain the local spread ,,.
In our in-depth analyses of 1500+ genomes, variability among clinical isolates was shown along the timeline, leading to distinct clustering of SARS-CoV-2 across the globe (unpublished observation). It was predicted, based on the aggregation propensity of the spike protein in the Wuhan and other isolates of SARS-CoV-2, that this virus would exhibit very high transmissibility and confer survival fitness,. Genetic diversity of the virus increases with disease progression and can be utilized to model the evolution and propagation of the disease Recently, phylogenetic network analysis of 160 SARS-CoV-2 genome samples showed a parallel evolution of the virus and its evolutionary selection in their human hosts. Similar whole-genome analyses of the Indian isolates and their comparison with global isolates can provide a better understanding of dominant clades within the population and unveil targets for developing specific interventions.
In the present study, machine learning-based t-SNE analysis of global clinical isolates has been utilized to segregate the clinical isolates into clusters while accommodating the outliers,. Whole-genome analysis of 3968 global isolates obtained from GISAID (Global initiative on sharing all influenza data), including 25 SARS-CoV-2 genomes sequenced in India [next-genome sequencing (NGS) data submitted by the ICMR-National Institute of Virology, Pune, India] and presented in [Figure 1] [Supplementary Figure 1 [Additional file 1]] and [Supplementary Table 1 [Additional file 2]], was an attempt to dissect the global genome diversity and also critically evaluate the placement of Indian isolates to understand the COVID-19 pandemic in India.
|Figure 1: Whole-genome-based t-SNE clustering of 3968 clinical isolates. (A) Comparative genome-based clustering of Indian isolates (red) with Chinease isolates (blue). (B) Comapartive genome-based clustering of Indian isolates (red) and Chinese isolates (blue) with rest of the world (green). (C) Diversity in clinical isolates showing three distinct clustering using hierarchical clustering on the t-SNE clusters. |
Click here to view
The initial cases reported from India had a travel history to China, which explained its position in a Chinese cluster [Figure 2]. The travel ban from China to India, in early February 2020, has prevented the large-scale spill-over directly from China to the Indian Sub-continent. However, various isolates transmitted from other South-East Asian countries might fall in the same cluster. The overlap of Indian samples majorly with European samples [Supplementary Figure 1], Panel III) reiterated the fact that the delayed travel restriction from the European hotspot regions affected not just India but also many countries.
|Figure 2: Position of various Indian isolates with other nations. (A-E) Clustering of SARS-CoV-2 genome sequences from India (red) with other nations around the globe. Indian samples clustered with samples from different nations – China, Kuwait, Canada, USA and Spain in whole-genome-based clustering. Figures were generated using FigTree v1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/).|
Click here to view
Hierarchical-based clustering further yielded exciting outcomes on the inter-continent transmission of COVID-19. The segregation of SARS-CoV-2 genomes into three clades indicates the emergence of evolutionary diversity [Figure 1]C. The heterogeneity of these clusters, grouped along with Chinese counterparts, validates a global spill-over event originating from Wuhan,. Hierarchical cluster 2 in Supplementary [Figure 1] Panel II (coloured by the continents) indicates the introduction of SARS-CoV-2 in India from the European, other Asian and North American nations [Supplementary Figure 1]. Detailed comparative analysis of Indian isolates with respect to other countries showed its close relationship with samples from China, USA, Canada, Spain and Kuwait, suggestive of exposure to COVID-19 due to travel history from these nations [Figure 2]. However, limited genome sequences from India make it difficult to differentiate and ascertain global transmission and transmission within the country.
The conservation of an amino acid in any protein sequence denotes its functional importance, as it undergoes fewer amino acid replacements or is more likely to substitute amino acids with similar biochemical properties. The amino acid conservation is inversely proportionate to the evolutionary rate. This is a valuable gauge of the evolutionary divergence and the analogous genomic regions. Sequence similarity between the open reading frames (ORFs) of Indian isolates and the initial sample collected in Wuhan unravels conservation in five ORFs corresponding to envelope protein, membrane glycoprotein, ORF6, ORF7b and ORF10 proteins [Figure 3]A. On the contrary, a number of mutations were observed in ORF1a, ORF1b, spike protein (surface glycoprotein), ORF3a, ORF7a, ORF8 and nucleocapsid phosphoprotein [Supplementary Figure 2 [Additional file 3]] and [Supplementary Table 2 [Additional file 4]]. Mean similarity calculated for these ORFs revealed that ORF1a in the Indian isolates was less conserved (more mutated) compared to global isolates ([Figure 3]A and [Supplementary Table 3 [Additional file 5]]). In all other ORFs, a relatively higher conservation was observed among Indian isolates compared to Wuhan strain. When compared with global isolates, Indian isolates have higher entropy for changes in ORF 1a and ORF 1b [Supplementary Figure 3 [Additional file 6]]. Further, qualitative analysis of mutations in non-conserved ORFs showed that each type of amino acid had undergone mutation in the Indian isolates [Figure 3]B. These mutations could be a major contributing factor for the separation of Indian isolates into three distinct clusters. Higher sampling rate driven by NGS of the Indian isolates would help in better understanding of actual variability in SARS-CoV-2 and assist both in identifying better diagnostic markers and in developing specific interventions in terms of vaccine candidates and drug targets.
|Figure 3: Sequence similarity and mutation analysis of open reading frames. (A) Comparison of mean sequence similarity for open reading frames between Indian and global isolates with Wuhan strain. (B) Qualitative analysis on type of mutations occurring in non-conserved open reading frames (ORFs) of Indian isolates compared to Wuhan strain.|
Click here to view
Evolutionary divergence, corroborated by epidemiological data, is a valuable tool to implement appropriate measures against this pandemic. The population density of India and the presence of functionally distinct isolates in the Indian population raise concerns and warrant an urgent need for higher sampling rate for better assessment of the evolution of SARS-CoV-2 in India. The situation is further confounded by the fact that many of these Indian isolates submitted in databanks include those of Indians living in Iran, Italian tourists visiting India, and also contains samples cultured in vitro.
In conclusion, a whole-genome diversity analysis of 3968 global clinical isolates, including 25 isolates sequenced in India, of SARS-CoV-2 was done. The variations in different open reading frames (ORFs) of SARS-CoV-2, which drives the formation of distinct Indian clusters and functional heterogeneity, were highlighted. Five ORFs corresponding to envelope protein, membrane glycoprotein, ORF6, ORF7b and ORF10 were found to be highly conserved, while a number of mutations were observed in ORF1a, ORF1b, spike protein, ORF3a, ORF7a, ORF8 and nucleocapsid phosphoprotein. Generating diverse genomic datasets will provide insight into the propagation dynamics of COVID-19, leading to a better understanding of pathogenesis and evolution of SARS-CoV-2, which will eventually lead to better intervention methods.
Acknowledgment:The seventh author (SEH) acknowledges Department of Biotechnology, Government of India for funding support (BT/PR23099/NER/95/632/2017), (BT/PR23155/NER/95/634/2017). SEH is a JC Bose National Fellow, Department of Science and Technology, Government of India & Robert Koch Fellow, Robert Koch Institute, Berlin. The first author (HS) is a recipient of Women Scientist fellowship, Department of Health Research and the second (JS) & fourth (SJ) authors received Young Scientist fellowships from the Department of Health Research, Ministry of Health and Family Welfare, Government of India. The sixth author (JAS) received UGC Startup grant and the third author (MK) received Silver Jubilee Post-Doctoral fellowship from Jamia Hamdard, New Delhi. Authors acknowledge the Originating and Submitting Laboratories for their sequences and meta-data shared through GISAID on which this study is based. Authors acknowledge BioInception Pvt. Ltd, for providing their proprietary data analysis pipeline and platform.
Financial support & sponsorship: None.
Conflict of Interest: None.
| References|| |
World Health Organization. Coronavirus disease 2019 (COVID-19) Situation report-94
. Geneva : WHO ; 2020.
Kimball A, Hatfield KM, Arons M, James A, Taylor J, Spicer K, et al
. Asymptomatic and presymptomatic SARS-CoV-2 infections in residents of a long-term care skilled nursing facility - King county, Washington, March 2020. MMWR Morb Mortal Wkly Rep
Mandal S, Bhatnagar T, Arinaminpathy N, Agarwal A, Chowdhury A, Murhekar M, et al.
Prudent public health intervention strategies to control the coronavirus disease 2019 transmission in India: A mathematical model-based approach. Indian J Med Res
Gupta N, Praharaj I, Bhatnagar T, Thangaraj JWV, Giri S, Chauhan H, et al
. Severe acute respiratory illness surveillance for coronavirus disease 2019, India, 2020. Indian J Med Res
Yadav PD, Potdar VA, Choudhary ML, Nyayanit DA, Agrawal M, Jadhav SM, et al
. Full-genome sequences of the first two SARS-CoV-2 viruses from India. Indian J Med Res
Jamal S, Singh J, Sheikh JA, Singh H, Khubaib M, Kohli S, et al
. Molecular Analyses of Over Hundred Sixty Clinical Isolates of SARS-CoV-2: Insights on Likely Origin, Evolution and Spread, and Possible Intervention. Preprints. 2020. doi: 10.20944/preprints202003.0320.v1.
Sheikh JA, Singh J, Singh H, Jamal S, Khubaib M, Kohli S, et al
. Emerging genetic diversity among clinical isolates of SARS-CoV-2: Lessons for today. Infect Genet Evol
Forster P, Forster L, Renfrew C, Forster M. Phylogenetic network analysis of SARS-CoV-2 genomes. Proc Natl Acad Sci U S A
Elbe S, Buckland-Merrett G. Data, disease and diplomacy: GISAID's innovative contribution to global health. Glob Chall
van der Maaten. Accelerating t-SNE using tree-based algorithms. J Mach Learn Res
Elbe S, Buckland-Merrett G. Data, disease and diplomacy: GISAID's innovative contribution to global health. Glob Chall. 2017; 1
Andersen KG, Rambaut A, Lipkin WI, Holmes EC, Garry RF. The proximal origin of SARS-CoV-2. Nat Med
Valdar WS, Thornton JM. Protein-protein interfaces: Analysis of amino acid conservation in homodimers. Proteins
Wu TT, Kabat EA. An analysis of the sequences of the variable regions of Bence Jones proteins and myeloma light chains and their implications for antibody complementarity. J Exp Med
[Figure 1], [Figure 2], [Figure 3]