Emergence of a new SARS-CoV-2 variant from GR clade with a novel S glycoprotein mutation V1230L in West Bengal, India

Sequence retrieval

High coverage full genome sequences of 2000 SARS-CoV-2 strains, collected during January 2021 to March 2021 from India, were retrieved from Global GISAID database on 15 May 2021 [16]. The genome sequence of the prototype SARS-CoV-2 strain hCoV-19/Wuhan/WIV04/2019 (GISAID accession no. EPI_ISL_402124) and several clades/variants was also downloaded from GISAID database.

Screening of mutations

In this study we have performed non-synonymous mutational analysis of 25 proteins (NSP1-NSP16, S glycoprotein, NS3, E, M, NS6, NS7a, NS7b, NS8 and N) encoded by each of the 2000 SARS-CoV-2 genome sequences submitted from India. For performing mutational analysis, each of the 25 protein coding regions of circulating SARS-CoV-2 genomes as well as prototype genome (hCoV-19/Wuhan/WIV04/2019) was cut and subsequently translated to amino acid sequence by using TRANSEQ nucleotide-to-protein sequence conversion tool (EMBL-EBI, Cambridgeshire, UK). Next, each of the 25 protein sequences of 2000 SARS-CoV-2 genomes was aligned with the corresponding protein sequence of the prototype strain by using MEGA software (Version X) and observed for amino acid substitutions [17].

Phylogenetic analysis

A phylogenetic dendrogram was constructed based on the whole genome sequences of 111 SARS-CoV-2 strains, including 38 sequences of the new variant and 73 reference sequences of different clades/variants (11 reference sequences of Indian variant B.1.617; 5 reference sequences of each of G clade, GR clade, GV clade, L clade, V clade, GRY clade/UK variant, South African variant, Brazilian variant, California variant, Nigerian variant and Indian variant B.1.618; 4 reference sequences of S clade; 3 reference sequences of GH clade), using Molecular Evolutionary Genetics Analysis (MEGA) version X [18]. Initially, genome sequences of 122 SARS-CoV-2 strains were aligned multiple alignment program MAFT version 7 [19]. The alignment file was then used to build phylogenetic tree by maximum-likelihood method using general time reversal (GTR) statistical model with 1000 bootstrap replicates [17]. We further performed the phylogenetic analysis of the 70 sequences of the new variant with the sequence analysis web app Nextclade^beta version 0.14.3 (https://clades.nextstrain.org/).

Results and Discussion

By performing the whole genome mutational analysis of 2000 SARS-CoV-2 strains collected during January 2021 to March 2021 from India, we have identified 70 SARS-CoV-2 strains having new set of 11 co-existing mutations: NSP4 has two mutations D279N and L353F; NSP8 has a single mutation of V26F; NSP12 has a single mutation of P323L; S gene is characterized by three mutations D614G, P681H and V1230L; NS3 has a single mutation of G172C; NS8 has a single mutation of V62L; and N gene is characterized by two mutations R203K and G204R (Table 1, Supplementary Table 1). Interestingly, among 70 sequences, 16 sequences also harbored E484K mutation along with 614G, P681H and V1230L in the S glycoprotein. This new variant has evolved within the GR clade (Pango lineage B.1.1 or 20B), which is characterized by four co-existing mutations D614G in S glycoprotein, P323L in NSP12, and R203K and G204R in N protein, all of which are present within this new variant. Phylogenetic analysis of 38 representative strains of this new variant along with 73 reference sequences of different clades and variants by MEGA X revealed that strains of the new variant emerged from GR clade (B.1.1) and formed a new cluster (Figure 1). In consistent with this result, the phylogenetic dendrogram produced by Nextclade^bera v.0.14.3 also showed a new cluster that encompasses the 70 sequences of the new variant within the 20B clade (GR clade) (Figure 2). We propose to name this variant as GRL or lineage B.1.1/S:V1230L or 20B/ S:V1230L due to the presence of V1230L in S glycoprotein along with GR clade specific mutations (Figure 3). Interestingly, all the 70 sequences of this new variant have been exclusively observed in the state of West Bengal out of 412 sequences (16.99%) submitted from there and distributed in 8 different districts (Kolkata, n=23; Nadia, n=4; North 24 Parganas, n=21; Howrah, n=10; Hooghly, n=4; South 24 Parganas, n=4; Bardhaman, n=2; Paschim Medinipur, n=2) (Supplementary Table 1). Among 70 strains, 34, 26 and 10 strains were collected during January, February, and March 2021 respectively. In addition to 12 co-existing mutations including E484K, 113 different mutations were found throughout the genome of these 70 SARS-CoV-2 sequences. Among these 113 mutations, 21 were found in NSP3, 16 were found in S glycoprotein, 13 were found in NSP2, 8 were found in NS3 and N, 7 were found in NS8, 6 were found in NSP6, 5 were found in both NSP5 and NSP16, 4 were found in each of NSP4, NSP12, NSP14 and NS7a, 3 were found in NSP15, 2 were found in E gene, and single mutation was found in each of NSP8, NSP9 and NS6. No mutation was observed in NSP1, NSP7, NSP10, NSP11, NSP13, M and NS7b genes (Table 2, Supplementary Table 1).

• Download figure
• Open in new tab

Figure 1:

Molecular phylogenetic analysis of new variant strains by MEGA X using maximum-likelihood method. The phylogenetic dendrogram is based on whole genome sequences of 38 representative strains of new variant and 73 reference strains of 16 clades/variants. The scale bar represents 0.00020 nucleotide substitution per site. The best-fit model used for constructing the phylogenetic dendrogram was the general time-reversal model (GTR).

• Download figure
• Open in new tab

Figure 2:

Molecular phylogenetic analysis of new variant strains by Nextclade^beta v.0.14.3. Left phylogenetic tree representing all the reference strains of different Nextstrain clades (Automatically selected by Nextclade) as well as 70 strains of the new variant. Each color depicting a different clade. Right phylogenetic tree highlighting only the 70 sequences of the new variant.

• Download figure
• Open in new tab

Figure 3:

A schematic presentation illustrating the evolution of various clades and variants with their characteristic signature mutations.

Emergence of new mutations in the S glycoprotein, which is involved in the binding of SARS-CoV-2 to the host cell receptor angiotensin-converting enzyme 2 (ACE2) and entry into the host cell, could have influence on virus binding, entry, immune invasion and antibody neutralization. The D614G, the first mutation appeared in the S glycoprotein, was first observed in Germany in January 2020 and became the dominant mutation in all the circulating strains worldwide by June, 2020 [4]. Now it is present in all the circulating strains including the newly emerged VOI (Variant of Interest) and VOC [20]. Patients infected with D614G mutant had higher nasopharyngeal viral RNA loads, indicating its role in increased infectivity [4]. Afterwards, experimental evidences confirmed that the D614G variant enhances virus replication in human lung epithelial cells and primary human airway tissues by increasing the infectivity and stability of the virions. Hamster infected with D614G variant produces higher infectious titers in nasal washes and the trachea, supporting the role of the G614G mutation in high transmissibility. It has also been demonstrated that D614G mutation confers higher susceptibility to serum neutralization [21]. Another study showed that pseudovirus particles carrying S-G614 enter ACE2 overexpressing cells more efficiently than those with S-D614. This increased entry correlates with less S1 domain shedding and higher S protein incorporation into the virion. However, S614 does not alter S protein binding to ACE2 or neutralization sensitivity to pseudoviruses [22]. P681 is located adjacent to the polybasic furin cleavage site ₆₈₁P-R-R-A-R-S₆₈₆ where furin mediated proteolytic cleavage is expected to occur between arginine (R₆₈₅) and serine (S₆₈₆). This furin cleavage site is positioned at the junction between S1 domain, essential for virus binding to ACE2 receptor, and S2 domain, necessary for the fusion between virus and the cell membrane. The furin mediated cleavage at the junction of S1/S2 is essential for virus entry into the host cells [23-24]. Therefore, any mutation at this site could influence S1/S2 cleavability by furin-like proteases and hence the infection properties of the virus. The P681H mutation has previously described in the UK variant [7]. A recent study has shown that P681H mutation increases spike protein cleavage by furin-like proteases but does not efficiently impact viral entry and cell-cell spread [25]. The E484K mutation, present within the receptor binding domain of S glycoprotein, has previously reported in Brazilian variant (P.1) and South African variant (B.1.351) [8, 10]. SARS-COV2 variants having E484K mutation is less sensitive to neutralization by serum antibodies induced by previous infection or vaccination [26]. The V1230L mutation, not observed previously in other country, is located at the transmembrane (TM) domain of S2 subunit of S glycoprotein. This transmembrane domain anchors S glycoprotein to the virus envelope. The replacement of valine with more hydrophobic amino acid leucin could tighten the association of S glycoprotein to the viral envelope and also increase S glycoprotein incorporation into the viral envelope. However, further studies are required to reveal the functional relevance of this mutation.

The new variant having D614G, P681H, V1230L and E484K mutations in the S glycoprotein is expected to have increased infectivity, high transmissibility, and lower sensitivity to neutralization antibodies. As India is already reporting high infection rates in the 2^nd wave compared to the 1^st wave in 2020, the emerging variant needs to be monitored closely. Therefore, further studies involving full genome sequencing of both severe and moderate COVID-19 cases in this region as well as antibody neutralization assay of this variant are required.