Coronaviruses are a large family of diverse enveloped RNA viruses found in numerous animal species. They are well-known for their ability to cross species barriers and have been transmitted from bats or intermediate hosts to humans on several occasions. Four of the seven human coronaviruses (hCoVs) are responsible for approximately 20% of common colds (hCoV-229E, -NL63, -OC43, -HKU1). Two others (SARS-CoV-1 and MERS-CoV) cause severe and frequently lethal respiratory syndromes but have only spread to very limited extents in human populations. In contrast, SARS-CoV-2, while exhibiting intermediate pathogenicity, has a profound impact on public health due to its enormous spread. In this review, we discuss which initial features of the SARS-CoV-2 Spike protein and subsequent adaptations to the new human host may have helped this pathogen to cause the COVID-19 pandemic. Our focus is on host forces driving changes in the Spike protein and their consequences for virus infectivity, pathogenicity, immune evasion and resistance to preventive or therapeutic agents. In addition, we discuss the significance and perspectives of broad-spectrum therapeutics and vaccines.

Conservation history of D614 residue is valuable in predicting the consequences of D614G mutation in the SARS-CoV-2 spike glycoprotein (SGP). We report here that the D614 belonged to an extraordinarily conserved, densely hydrophobic eleven amino acid peptide-motif vavlyqdvnct (11-aa), in the Sarbecovirus group and the variant carrying vavlyqdvnct had appeared in Chinese samples and became predominant in several geographical hotspots in late March, 2020. Interestingly a 2009 annotation of SARS-CoV contained the same mutation.

As reported by us and others previously (1, 2), the D614G mutation appeared in the spike glycoprotein (SPG) of the SARS-CoV-2 (the pathogen behind COVID-19) at the early stages of the pandemic and then G614 containing variant of SARS-CoV-2 became the predominant strain in most human populations across the world. However, one of the most recent reports from India (3) stated the incidence of G614 to be only 26% in the Indian population. This report is contradictory to the information available through the GenBank (4) SARS-CoV-2 sequence deposits made by various laboratories from India. The above stated report currently circulating in the Indian media is likely to create a public perception that the Indian strain is less contagious and such a notion could be harmful to people’s welfare. In view of this concern we have re-evaluated, updated and recalculated the incidence of the G614 variant in the Indian population by analyzing 395 Indian SARS-CoV-2 genomic sequences available in the GenBank as of June 26, 2020. In our analysis we have categorized the samples by the month in which the samples were collected. We have used an alignment-free software tool named Compare (5, 6), and the Basic Local Alignment Search Tool (BLAST) (7) in the present analysis. We finally inspected each of the 395 sequences physically for the presence of aspartic acid (D) or glycine (G) at the 614^th position of the spike glycoprotein. We analyzed an Australian cohort in parallel for comparison. We have found that the prevalence of G614 variant in the Indian samples for the month of June 2020 is 90.6%. The trends are similar with the Australian samples.

Both SARS-CoV-2 (COVID-19) and SARS coronaviruses (CoVs) are members of the subgenus Sarbecovirus. To understand the origin of SARS-CoV-2, protein sequences from sarbecoviruses were analyzed to identify highly-specific molecular markers consisting of conserved inserts or deletions (termed CSIs) in the spike (S) and nucleocapsid (N) proteins that are specific for either particular clusters/lineages of these viruses or are commonly shared by specific lineages. Three novel CSIs in the N-terminal domain of the spike protein S1-subunit (S1-NTD) are uniquely shared by the SARS-CoV-2, BatCoV-RaTG13 and most pangolin CoVs, distinguishing this cluster of viruses (SARS-CoV-2r) from all others. In the same positions, where these CSIs are found, related CSIs are also present in two other sarbecoviruses (viz. CoVZXC21 and CoVZC45 forming CoVZC cluster), which form an out group of the SARS-CoV-2r cluster. These three CSIs are not found in the SARS-CoVs. However, both SARS and SARS-CoV-2r CoVs contain two large CSIs in the C-terminal domain of S1 (S1-CTD), which binds the human ACE-2 receptor, that are absent in the CoVZC cluster of CoVs. These results indicate that while the S1-NTD of the SARS-CoV-2r viruses possesses the sequence characteristics of the CoVZC cluster of CoVs, their S1-CTD resembles the SARS viruses. Thus, the spike protein of SARS-CoV-2r viruses has likely originated from a recombination event between the S1-NTD of the CoVZC viruses and the S1-CTD of SARS viruses. This inference is also supported by the amino acid sequence similarity of the S1-NTD and S1-CTD from SARS-CoV-2 compared to the CoVZC and SARS CoVs. We also present evidence that one of the pangolin-CoV_MP789, whose receptor-binding domain is most similar to the SARS-CoV-2, is also derived by a recent recombination between the S1-NTD of the CoVZC CoVs and the S1-CTD of a SARS-CoV-2 related virus. Several other identified CSIs are specific for others clusters of sarbecoviruses including a clade consisting of bat SARS-CoVs (BM48-31/BGR/2008 and SARS_BtKY72). Structural mappings studies show that the identified CSIs are located within surface-exposed loops and form distinct patches on the surface of the spike protein. These surface loops/patches are predicted to interact with other host components and play important role in the biology/pathology of SARS-CoV-2 virus. Lastly, the CSIs specific for the SARS-CoV-2r clade provide novel means for development of new diagnostic and therapeutic targets for these viruses.