SARS-CoV-2 Variant of Concern 202012/01 has about twofold replicative advantage and acquires concerning mutations

Frederic Grabowski, Grzegorz Preibisch, Stanisław Giziński, Marek Kochańczyk, Tomasz Lipniacki
doi: https://doi.org/10.1101/2020.12.28.20248906

Abstract

The novel SARS-CoV-2 Variant of Concern (VOC)-202012/01 (also known as B.1.1.7), first collected in United Kingdom on September 20, 2020, is a rapidly growing lineage that in January 2021 constituted 86% of all SARS-CoV-2 genomes sequenced in England. The VOC has been detected in 40 out of 46 countries that reported at least 50 genomes in January 2021. We have estimated that the replicative advantage of the VOC is in the range 1.83–2.18 [95% CI: 1.71–2.40] with respect to the 20A.EU1 variant that dominated in England in November 2020, and in range 1.65–1.72 [95% CI: 1.46–2.04] in Wales, Scotland, Denmark, and USA. As the VOC strain will likely spread globally towards fixation, it is important to monitor its molecular evolution. We have estimated growth rates of expanding mutations acquired by the VOC lineage to find that the L18F substitution in spike has initiated a substrain of high replicative advantage in relation to the remaining VOC substrains. The L18F substitution is of significance because it has been found to compromise binding of neutralizing antibodies. Of concern are immune escape mutations acquired by the VOC: E484K, F490S, S494P (in the receptor binding motif of spike) and Q677H, Q675H (in the proximity of the polybasic cleavage site at the S1/S2 boundary). These mutants may hinder efficiency of existing vaccines and expand in response to the increasing after-infection or vaccine-induced seroprevalence.

1 Introduction

The earliest genome belonging to the novel SARS-CoV-2 Variant of Concern (VOC)-202012/01, also known as B.1.1.7 lineage, was collected on September 20, 2020, in Kent, UK (GISAID sequence accession ID: EPI_ISL_601443). The lineage, characterized by nine spike protein mutations (deletions: 69–70 HV, 145V; substitutions: N501Y, A570D, D614G, P681H, T716I, S982A, D1118H), started to spread rapidly in mid-October 2020 to constitute in January 2021 86% of all SARS-CoV-2 genomes sequenced in England [1, 2]. Spread of the VOC-202012/01 variant, hereafter referred to as the VOC, co-occurred with a rapid surge of cases in December in Kent and Greater London [3, 4].

In England, the VOC is currently replacing the recently dominant 20A.EU1 strain, characterized by A222V substitution in spike protein [5]. Both strains are independent substrains of a spike glycoprotein D614G variant that has spread in spring 2020 in England and worldwide, almost reaching fixation [6]. The 20A.EU1 strain started expanding in England in mid-August 2020 and constituted more than 65% of genomes sequenced in England in November 2020 [2]. The VOC strain in the receptor-binding motif (RBM) of spike shares mutation N501Y with the 501Y.V2 and P.1 strains that are currently rapidly spreading in South Africa and Brazil, respectively [7, 8].

Deletion 69H–70V in spike glycoprotein, which is characteristic, but not unique, to the VOC, prevents detection of the spike gene by the dPCR probe used by some laboratories of the English diagnostic system (spike gene target failure, SGTF) [9]. As the VOC has a high multiplicative potential, it has become the most prevalent Δ69–70HV strain; consequently, the proportion of the SGTF has been used as a proxy for the prevalence of the VOC genome [9]. Based on the SGTF, the Public Health England agency determined the multiplicative advantage of the VOC in NHS STP areas of England in weeks 44–48 of 2020. As an average for the considered STP areas the authors obtained the ratio of reproduction numbers equal 1.47 [95% CI: 1.34–1.59] [9]. Leung et al., based on GISAID data from the period September 22–November 16, 2020, and a competition transmission model of two viruses, estimated this ratio as 1.75 [95% CrI: 1.70–1.80] [10]. Davies et al., based on data from the COVID-19 Genomics UK (COG-UK) Consortium from October and November 2020 estimated that the VOC is 43–82% [95% CrI: 38–106%] more transmissible than preexisting variants of SARS-CoV-2 [11, 12]. Our early estimate of the ratio of reproduction numbers of the VOC to non-VOC strains, based on GISAID data in weeks 43–47 of 2020 in England was 2.24 [95% CI: 2.03–2.48] [13].

In this study, based on GISAID data available on February 12, 2021 (when this manuscript has been revised), we calculated the growth of the ratio of the VOC to the 20A.EU1 genome sequences collected in the period between week 43 and week 51 of 2020 (October 19–December 20, 2020) and estimated the replicative advantage of the VOC strain in relation to the 20A.EU1 strain. Then, using the same approach, we estimated the replicative advantage of the 20A.EU1 strain in relation to previous D614G strains and the D614G strain in relation to the D614 strains (i.e., the strains with non-mutated residue 614). In such approach, we analyze the progression of strains with increasing replication advantage: D614G → 20A.EU1 → the VOC. In contrast to the approach of Public Health England (PHE) agency [9] and Davies et al. [11], we used England-aggregated weekly data. We prefer to use such aggregated data because of the geographically non-uniform emergence of VOC substrains that, as we will show, have diverse replication advantages. Some substrains that emerge locally may have low or no replicative advantage and in the long term may be out-competed and eventually eliminated; however, such transient substrains can considerably contribute to the genome composition in specific NHS STP areas of England.

The VOC strain, as all strains, mutates continuously. Because of its significant replicative advantage, any accrued mutations gain an opportunity to spread (potentially worldwide), depending on their replicative advantage with respect to the bulk VOC strain, or higher ability to infect seroprevalent individuals. We have systematically estimated growth rates of spreading mutations acquired by the VOC lineage to find that spike L18F substitution has initiated a sub-strain of high replicative advantage. The L18F mutation is of significance because when recently analyzed in the context of the South African strain 501Y.V2 it has been found to compromise binding of neutralizing antibodies [14, 15].

2 Materials and Methods

2.1 Data sources

In Figures 1–7 we use weekly aggregated data of samples collected in England or South Africa submitted to GISAID until February 12, 2021, and samples collected in Wales, Scotland, Den-mark, and the USA and submitted until February 17, 2021. Table 1 is based on GISAID data as of February 16, 2021. In supplementary Figure S1, we use COG-UK pillar 2 data as of February 15, 2021. All data used for Figures 1–7 and supplementary Figure S1 are collected in a single supplementary file, Data File S1.

View this table:
Table 1: Fraction of VOC genomes in countries that in January 2021 reported more than 50 genomes.
Figure 1.
Figure 1.

Fixation of D614G strain and competition between its substrains: 20A.EU1 and the VOC in England. The proportion of D614G (orange), 20A.EU1 (cyan), VOC (magenta) and D614 strains (black) in weeks 10–53 of 2020 has been determined based on GISAID data available on February 12, 2021 (provided in supplementary Data File S1).

Figure 2.
Figure 2. The replicative advantage of the VOC, 20A.EU1, and the D614G strain.

(A) The ratio of the VOC to 20A.EU1 sequences collected in weeks 38–51 of 2020 in England. The trend line is fitted to data points from weeks 43–51 (blue) and from weeks 43–47 (red). The weekly growth rate is 1.88 [95% CI: 1.75–2.01] for weeks 43–51 and 2.24 [95% CI: 2.03–2.48] for weeks 43–47. (B) Stabilization of fits shown in panel A with time-shifting date of the last submission. (C) The ratio of the 20A.EU1 to non-20A.EU1 non-VOC D614G genomes collected in weeks 29–51 of 2020 in England. The trend line is fitted to data points from weeks 34–45. The weekly growth rate is 1.25 [95% CI: 1.23–1.28]. (D) The ratio of the D614G to D614 genomes collected in weeks 11–27 of 2020 in England. The trend line is fitted to data points from weeks 11–24. The weekly growth rate is 1.44 [95% CI: 1.40–1.48]. Data points for selected weeks are labeled with ratios of the counts of genome sequences. Panels A, C, and D are based on GISAID data submitted till February 12, 2021 (provided in supplementary Data File S1).

Figure 3.
Figure 3. The replicative advantage of the VOC in Denmark, Scotland, Wales, and USA.

(A) The ratio of the VOC to the 20A.EU1 sequences collected in Denmark. The trend line is fitted to data points from week 50 of 2020 to week 4 of 2021. (B) The ratio of the VOC to 20A.EU1 sequences collected in Scotland. The trend line is fitted to data points from weeks 50 of 2020 to week 4 of 2021. (C) The ratio of the VOC to 20A.EU1 sequences collected in Wales. The trend line is fitted to data points from week 50 of 2020 to week 4 of 2021. (D) The ratio of the VOC to non-VOC sequences collected in USA. The trend line is fitted to data points from week 51 of 2020 to week 3 of 2021. Data points are labeled with ratios of the counts of genome sequences. The weekly growth rates and respective confidence intervals are given in each panel. All panels are based on GISAID data submitted till February 17, 2021 (provided in supplementary Data File S1).

Figure 4.
Figure 4.

The replicative advantage of the 501Y.V2 strain in South Africa. Shown is the ratio of the 501Y.V2 to the non-501Y.V2 genomes collected in weeks 41–50 of 2020 in South Africa. The trend line is fitted to data points from weeks 43–50. The weekly growth rate is 1.58 [95% CI: 1.45–1.72]. Data points for selected weeks are labeled with ratios of the counts of genome sequences. The figure is based on GISAID data submitted till February 12, 2021 (provided in supplementary Data File S1).

Figure 5.
Figure 5. Mutations in spike protein in the VOC lineage in England in relation to the first VOC genome collected in September 20, 2020 (GISAID sequence accession ID: EPI_ISL_601443):

(A) Dots denote genomes collected in a given date (horizontal axis) with a respective number of novel (amino acid-level) mutations in relation to the first VOC genome (vertical axis). The number of the genomes with a given number of novel mutations is provided as gray numbers next to a brace. (B) The average number of novel (amino-acid level) mutations in a VOC lineage sequence (weekly average).

Figure 6.
Figure 6.

The growth of VOC substrains with respect the whole VOC strain. The figure is based on spike protein sequences from genomes collected in England and submitted to GISAID till February 12, 2021. The growth of the whole VOC strain (red line) in the time span between t and t0 = February 12, 2021 is calculated as VOC(t0)/VOC(t), i.e., the ratio of the number of the VOC genomes collected till February 12, 2021 and the number of VOC genomes collected till t (date given on horizontal axis). Disks denote VOC substrains defined by particular mutations (with at least 30 genomes submitted); horizontal axis value for each circle is the date of first collection of a sequence with a given mutation; vertical axis value is the number of genomes collected till February 12, 2021, conferring a given mutation. Red, orange, and green circles denote substrains with leading mutation in respectively in RBD, NTD, and signal peptide; black circles denote mutations in other domains. Selected substrains that may have replicative advantage and all substrains with the leading mutation in RBD are annotated. Detailed data that were analyzed to prepare this figure are provided in supplementary Data File S1 (sheet ‘Mutations’).

Figure 7.
Figure 7.

The replicative advantage of the spike L18F substrain in relation to the L18 VOC strains. The ratio vL of the number of VOC genomes conferring spike L18F mutation to the number of non-mutated (L18) VOC genomes collected in the period between week 49 of 2020 and week 3 of 2021 in England. Data aggregated into weeks indicated that vL changes from 7:2887 = 0.0024 in week 50 of 2020 to 329:8049 = 0.041 in the second week of 2021, i.e., 16.9 times, meaning that vL increases 16.91/5 = 1.76-fold per week. The trend line is fitted to data points shown as filled circles. The weekly growth rate of the ratio is 1.75. The 95% CI calculated as 1.96 × standard error of the slope is [1.66–1.85]; the 95% CrI calculated assuming binomial distribution of substrain genomes is [1.59–2.07]. Data points in week 50 of 2020 and week 2 of 2021 are labeled with ratios of the counts of both types of genome sequences. The figure is based on GISAID data submitted until February 12, 2021 (provided in supplementary Data File S1).

2.2 Genome sequence analysis

We analyzed mutations in spike gene sequences that had at most 5% letters other than A, C, G, T. Sequences devoid of full daily collection date were excluded from the analysis. The spike gene was localized using EMBOSS stretcher and then re-aligned using EMBOSS needle [16] to the reference sequence: GenBank’s NC_045512.2 from Wuhan or GISAID’s EPI_ISL_601443 from Kent. To minimize ambiguous reporting of mutations variants, indels were left-aligned using an in-house script; insertions and deletions on consecutive residues were collated and considered a single mutation.

2.3 Monte Carlo estimation of the confidence interval for the growth rate of the ratio of strains

To estimate the 95% credible interval for k subsequent weeks a priori we assume that genomes of two compared strains in each week of the considered period follow a respective binomial distribution having a success probability p = nx/(nx + ny), where nx and ny are the numbers of two compared strain genomes (e.g., VOC and 20A.EU1). By sampling from k such binomial distributions for the considered time window of k weeks 105 times, we obtained 105 series of k simulated sequenced genome proportions. We performed fitting to each such series to obtain 105 estimates of the weekly growth rate of the ratio of the genome sequences. In all cases but one for L18F analysis (Figure 7) the 95% credible interval obtained using the a priori method has been narrower than the 95% confidence interval calculated as 1.96 × (standard error of the slope). In these cases we reported confidence interval, while in L18F case we reported the credible interval.

2.4 Estimation of the relative replication advantage of viral strains

To estimate the ratios of replication numbers of strain x and y we estimated the weekly growth of the ratio of the number of their sequenced genomes, v. Then, assuming that both strains have the same average serial interval of 6.73 days [17], we obtained the ratio of their replication numbers Embedded Image. This estimation of the relative replicative advantage of viral strains does not involve a direct calculation of their reproduction numbers.

3 Results

3.1 Evolution of SARS-CoV-2 in England

The first prevailing mutation of SARS-CoV-2 was D614G substitution in spike protein (the first GISAID reported genome, EPI_ISL_913915, was collected on January 2, 2020, in Mexico). This substitution initiated a strain that spread worldwide nearly reaching global fixation [6]. In England, the D614G strain appeared during the spring wave of epidemic in 2020, and in summer it exceeded 98% of all sequenced genomes (Figure 1). Its substrain, 20A.EU1, started expanding in England in week 31 of 2020 and reached its maximum of 68% of all sequenced genomes in week 44 (Figure 1). The VOC strain is also a substrain of the D614G lineage, independent of 20A.EU1, which started expanding in week 43 (before week 43, less than 5 VOC genomes were collected per week), and in week 51 reached 57% of all sequenced genomes. At the same time the proportion of the 20A.EU1 strain dropped to 35% and the proportion of all other genomes dropped to 8%. One can observe that in weeks 44–51 the decrease of the proportion of 20A.EU1 is about twofold (from 68% to 35%), while the decrease of the proportion of remaining strains (that is, strains other than considered 20A.EU1 and VOC) is nearly fourfold (from 31% to 8%).

This motivates us to calculate the replicative advantage of the VOC strain in relation to the 20A.EU1 strain instead of all non-VOC strains (the replicative advantage of the VOC over non-VOC strains would depend on the proportion of the 20A.EU1 strain in all non-VOC genomes).

3.2 Replicative advantage of the VOC over the 20A.EU1 strain in England

In the eight-week period of weeks 43–51, the ratio q of the VOC to the 20A.EU1 strain increased from q = 30/4030 ≈ 0.0074 to q = 4928/3051 ≈ 1.62, that is, 217 times. This implies that q increased 2171/8 ≈ 1.96-fold per week. To estimate the growth of q in a more rigorous way, we fitted the trend line using two fitting windows (Figure 2A). Fitting in weeks 43–47 gives weekly growth v = 2.24 [95% CI: 2.03–2.48], whereas fitting in weeks 43–51 gives v = 1.88 [95% CI: 1.75–2.01]. These confidence intervals are estimated as 1.96 × standard error of the slope. However, as the estimate for weeks 43–47 is based on only 5 data points and the number of the VOC genomes in the first data point is small (n = 30), we additionally estimated the 95% credible interval for weeks 43–47 a priori (assuming binomial distributions of VOC and 20A.EU1 strain genomes, see Methods). Using this auxiliary method, we estimated that the 95% CrI is 2.09–2.45, which is somewhat narrower than the 95% CI calculated from the standard error of the slope. This demonstrates that the CI calculated from the standard error of the slope was not a result of an incidental linearity of data points.

To estimate the ratio of reproduction numbers of the VOC and the 20A.EU1 strain, Embedded Image, we assumed that both strains have the same average serial interval of 6.73 days [17]. Then, Embedded Image for the fitting window in between weeks 43 and 47, and Embedded Image for fitting window in between weeks 43 and 51. The eight-week window estimate gives a smaller advantage of the VOC than the four-week window. We think that this discrepancy is caused by two factors. As shown in the report of PHE [3] in week 51 of 2020, the VOC almost reached fixation in the Greater London and Kent, being nearly absent in central England. Additionally, more stringent measures implemented in Greater London and Kent limited the growth of absolute numbers of VOC cases. The heterogeneity in ratio q is not important as long as q is low (as it was until week 47); however, when q becomes high in some subregions, growth of q in the whole region decelerates. For this reason, we have not extended the fitting window past week 51. Since we may not rule out a chance that the higher growth of q in the window of 43–47 weeks is an artifact caused by a small number of data points, we conclude that that ratio Embedded Image is in the range 1.83–2.18 [95% CI: 1.71–2.40].

We should notice that sequenced genomes are submitted to GISAID with some time delay after sample collection; however, as of February 12, the data for weeks 44–51 of 2020 appears nearly complete. As shown in Figure 2B, replicative advantage estimates stabilize with the progression of the ‘date of last submission’. We have also verified our fits using the COG-UK genome database. We performed the analysis based on pillar 2 genomes that excludes routine tests of health and care workers and other tests made for particular purposes (pillars 1, 3, and 4)1. For the COG-UK data we obtained nearly the same (± 3%) replicative advantage of the VOC strain in the 43–47 and 43–51 week windows, see Supplementary Figure S1 for fits and values.

In the same way we estimated that the ratio of Rt of the VOC strain to Rt of other strains (that are neither VOC nor A20.EU1) is in the range of 2.03–2.47 [95% CI: 1.89–2.77], where the lower bound is the estimate within weeks 43–51 and the higher bound is the estimate for weeks 43–47 of 2020. As may be expected, Embedded Image, which reflects the fact that the A20.EU1 strain has a replicative advantage over previously dominant non-VOC D614G strains. We estimated this advantage by fitting a trend line to data points from weeks 34–45 of 2020, because for this period the exponential growth of the ratio of the A20.EU1 strain to other non-VOC D614G strains is observed (Figure 2C). In this period, the ratio of the A20.EU1 to other non-VOC D614G strains grows at a rate of 1.25 [95% CI: 1.23–1.28] per week, which gives Embedded Image. Finally, we estimated that the ratio Embedded Image. Concluding, we showed that the D614G strain that spread worldwide towards fixation had replicative advantage of 1.42 in relation to D614 strains in England. Its substrain A20.EU1 had replicative advantage of 1.24 over bulk D614G, and reached the proportion of 68% of genomes in England. Currently, A20.EU1 is outcompeted by the VOC that has about two-fold replicative advantage in relation to the A20.EU1 strain.

3.3 Worldwide spread of the VOC strain

London serves as a major transportation hub and thus, unsurprisingly, among 46 countries that reported SARS-CoV-2 genomes in January 2021, 40 countries reported a VOC genome from this period. We estimated the ratio of the VOC genomes to all genomes in these countries (Table 1). We found that, in addition to England, in 10 countries the VOC genomes constituted more than half of reported genomes. The data suggest that the strain is spreading globally, even though from countries other than England less than 20% of VOC genomes were reported (as of February 16, 2021).

Using the same method as in the previous section we estimated the replicative advantage of the VOC strain in four other countries in which the number of reported VOC genomes permits such analysis. In Denmark, Scotland, and Wales, similarly to England, the A20.EU1 strain constitutes a large share of all genomes in November 2020, respectively: 39%, 64%, and 70%. We thus compared the VOC with the A20.EU1 strain in these countries. In the USA, where we have not found a dominating D614G substrain, we compare the VOC with all non-VOC genomes. In all four countries, we found periods of the exponential growth of the ratio of compared genomes (seven-week-long in Denmark, Scotland, and Wales, and five-week-long in the USA; see linear growth in logarithmic scale in Figure 3). This allowed us to estimated the replicative advantage of the VOC strain in a narrow range: from 1.696.73/7 = 1.65 for Wales to 1.766.73/7 = 1.72 for the USA with 95% confidence intervals in the range [1.46–2.04]. We should notice that in all four countries the replicative advantage of the VOC strain has been found smaller than estimated for England.

3.4 Replicative advantage of the 501Y.V2 strain in South Africa

In addition to Δ69–70HV, mutation N501Y in the RBD of spike is considered as the most important recent mutation [9]. This mutation occurred independently in the South African strain 501Y.V2, where it is accompanied by two other mutations in spike RBD: K417N and E484K [7]. Using the same method as previously we estimated the replicative advantage of the 501Y.V2 strain over other South African strains (Figure 4). In weeks 43–50 of 2020, when an exponential growth is observed, the ratio of 501Y.V2 strain to other strains grows at the weekly rate of 1.58 [95% CI: 1.45–1.72], which gives Embedded Image. The noisiness of data is associated with the small number of genomes sequenced in the entire period of weeks 43–50 (only 407 501Y.V2 sequences and 376 non-501Y.V2 sequences collected in this period were submitted to GISAID).

3.5 Emergence of mutations in VOC genomes in England

Because of its high replicative advantage, the VOC strain will likely become globally dominant, possibly reaching fixation. It is thus crucial to track mutations that arise in this strain, that could further increase its replicative advantage. We thus performed sequence analysis of the spike gene in all genomes from England submitted to GISAID till February 12, 2021 (supplementary Data File S1, sheet ‘Mutations’). There are 2232 different mutations in genomes collected in England, including 1213 different mutations in VOC genomes with 697 ‘confirmed mutations’ found in more than one submitted VOC genome. Majority of these mutations (535 out of 697) have arisen also in non-VOC strains. This suggest that the mutational space is to a large extent already explored, however the nine mutation characterizing the spike protein of VOC may increase or decrease the replicative advantages of recurrent mutations or allow for propagation of novel mutations. In Figure 5A we show accumulation of mutations is VOC spike in time; within 53,185 analyzed VOC genomes, 20% have at least one mutation, while 1.7% has two mutation or more. In relation to the first collected VOC lineage genome (GISAID accession ID: EPI_ISL_601443), the VOC genomes collected in the end of January have on average 0.3 mutation in their spike protein (Figure 5B).

In Figure 6 we analyze the replicative advantage of 60 VOC substrains for which there were at least 30 genomes submitted to GISAID from England. Each substrain is marked by disk on the (first collection date, total occurrences) plane. The VOC substrains that from time of its first collection grow faster than the VOC strain on average are above the red line that shows the average growth of the VOC strain sequences. Unsurprisingly, the majority of substrains that exceed 30 submitted genomes fall into this category. This analysis should be taken with great caution, as the total number of substrains is high, some of may grow faster in the considered time window just by chance without having any replicative advantage. Nevertheless, one can use this approach to screen for further analysis the substrains that can potentially have a replicative advantage. The most prevalent mutation, defining a fast growing variant is the L18F substitution (1186 genomes) in the N-terminal domain (NTD) of spike protein. The second most prevalent is L5F (658 genomes) localized in the signal peptide of spike. This mutation (at a highly homoplasic position that may be a sequencing artifact [18]) was found abundant also in non-VOC genomes. The third most prevalent mutation, that is also by far the most prevalent RBM as well as in whole RBD mutation, is S494P (441 genomes). The other two fast growing “sibling” mutations in VOC — Q677H (256 genomes) and Q675H (86 genomes) — are present in the the proximity of the polybasic cleavage site (residues 682–685) at the S1/S2 boundary influencing RBD−ACE2 binding [19]. Mutations at residue Q677 (either Q677H or Q677P) were found in several independent lineages spreading over the autumn of 2020 and into the winter of 2021 in the USA [5].

3.6 Replicative advantage of the L18F substrain

The first occurrence of the spike L18F substitution has been reported in a VOC strain genome collected on December 4, 2020 (GISAID ID: EPI_ISL_720875). As of February 12, 2021, as much as 1186 spike L18F VOC genomes have been reported in England. Of note, in Autumn 2020, that is, before the VOC lineage has become the dominant strain, the L18F substitution was a ubiquitous mutation in England. Till February 12, 2021, most of the L18F non-VOC genomes in England (97.6%, 25,655 out of 26,280) were found within the 20A.EU1 strain. The fraction of spike L18F mutation in the expanding 20A.EU1 strain was slowly increasing from 35% (1332 out of 3799) in September, 43% (5658 out of 13,046) in October, till 52% (8917 out of 17,470) in November 2020, which may suggest that this mutation was beneficial for the 20A.EU1 strain.

In Figure 7 we show the exponential growth of the L18F VOC substrain in England in the five-week period of December 7, 2020–January 17, 2021, in relation to the VOC genomes non-mutated at residue 18, denoted L18. In the considered period, the ratio of the L18F to the L18 genomes increased with the fitted weekly growth rate of 1.75, which gives Embedded Image. This credible interval is calculated assuming a binomial distribution of the number of the L18F and L18 VOC genomes in each week (see Methods). The confidence interval calculated from 1.96 × standard error of the slope, [1.63–1.81], is narrower, which means that the nearly perfect co-linearity of five data points is somewhat coincidental, and the binomial distribution-based credible interval is the proper estimate. This analysis suggests a high replication advantage of the L18F VOC substrain in relation to the remaining VOC genomes, but since it is based on very incomplete data, it must be taken with caution. The finding is supported by data from Wales, UK, where the L18F VOC genomes constituted 17% (390 out of 2302) of all VOC genomes reported in January, substantially more than the number of genomes submitted to GISAID from the same period in England, 3.0% (1333 out of 43,700) till February 18, 2020. The number of genomes in Wales is however too small to perform an analysis analogous to that in Figure 7.

3.7 VOC strain mutations in spike receptor-binding domain

Of particular concern are the VOC strain mutations occurring in the receptor-binding domain (RBD, residues 333–527), especially mutations in the receptor-binding motif (RBM, residues 438–506). These mutations may potentially lead to immune escape mutants, resulting in reinfection of convalescent individuals and aggravation of the efficacy of current vaccines. Propagation of such mutations is facilitated by high replicative advantage of the VOC strain and potential selection due to the increasing number of convalescent or immunized individuals. The VOC-202012/01 strain spike RBM mutations of special concern are substitutions E484K and S494P.

E484K

A first genome has been collected on December 17, 2020 (GISAID ID: EPI_ISL_782148), and there were 30 genomes reported from the England up till February 12, 2021. The same mutation has occurred in the fast expanding South African and Brazilian (Manaus) strains that share with the VOC substitution N501Y and additionally have a mutation of residue 417: either K417N (South African strain 501Y.V2) [7] or K417T (Manaus strain P.1) [8]. It was suggested that E484K may compromise binding of class 2 neutralizing antibodies, while the A501V mutation interferes with binding of class 1 antibodies. The P.1 strain led to the surge of infections in Manaus in December 2020 despite high seroprevalence of the population (a study of blood donors indicated that 76% [95% CI: 67–98%] of the population in Manaus had been infected with SARS-CoV-2 by October 2020 [20]).

S494P

A first genome has been collected on November 12, 2020 (GISAID ID: EPI_ISL_741039), and there were 441 genomes reported from England up till February 12, 2021. In an in silico study, this substitution has been found to increase complementarity between the RBD and ACE2 [21]. This mutation has been also characterized as an escape mutation by Koenig et al. [22], who also distinguished five additional “escape” residues in the RBM: G447, Y449, L452, F490, G496, and six outside the RBM but within the RBD: Y369, S371, T376, F377, K378, R403. Among these residues, until February 12, 2021, substitution F490S (first collected on December 13, 2020, GISAID accession ID: EPI_ISL_736026) was reported in the highest number of genomes (28 genomes in England).

4 Discussion

The mutations of SARS-CoV-2 that substantially increase replicative advantage of emerging strains will likely become dominant, either locally in countries or continents, or worldwide. Substitution D614G in spike protein (the first GISAID reported genome, EPI_ISL_913915, was collected on January 2, 2020, in Mexico) initiated a strain with replicative advantage over D614 strains estimated based on data from England as 1.42 [95% CI: 1.38–1.45]. The D614G strain has spread worldwide nearly reaching fixation; it was present in more than 99% of genomes collected worldwide in January 2021. The 20A.EU1 strain, a substrain of D614G that harbors A222V mutation in spike, emerged in Spain in early summer, 2020, spread over Europe, becoming the dominating strain (more than half of sequenced genomes) in several countries (Spain, England, Scotland, Wales, Ireland and Italy) in November, 2020, but was nearly absent outside of Europe. Based on GISAID data in England we estimated its replicative advantage over other non-VOC D614G English strains as 1.24 [95% CI: 1.22–1.27]. The VOC strain started spreading in England in October 2020, outcompeting the A20.EU1 strain, and reached about 80% of genomes in England in January 2021. We have estimated that its replicative advantage over the A20.EU1 strain is in the range 1.83–2.18. The lower bound was obtained by fit in the eight-week-long period of weeks 43–51 of 2020, when the ratio of the VOC to the A20.EU1 strain genomes increased 217 times, from 0.0074 to 1.62, whereas the upper bound was obtained in the four-week-long period of weeks 43–47 of 2020. We think that the slower growth in the period of weeks 47–51 is a consequence of (1) the fact that in Kent, Greater London, and their vicinity, the VOC strain almost reached fixation and (2) the fact that in these regions more stringent measures were implemented to suppress rapid growth of cases.

We also estimated the replicative advantage of the VOC strain in relation to the 20A.EU1 strain in Denmark, Scotland. and Wales, and in relation to bulk non-VOC strains in USA. We found find values in range from 1.69 for Wales to 1.76 for USA, with 95% confidence intervals in range [1.46–2.04], that are smaller than in England. One possible explanation is that the VOC strain is able to infect seroprevalent or exposed individuals. Such ability would increase its replicative advantage in a population, in which the fraction of seroplevalent individuals is large [23]. This mechanism would also help to explain the higher replicative advantage of the VOC strain observed in weeks 43–47 of 2020, when the VOC strain was gaining prevalence in London area [3]. The strain P.1, found in Brazil, that shares with VOC RBM mutation N501Y, caused recently the second (higher) wave of deaths in Manaus despite high seroprevalence of the population [8].

In addition to double deletion Δ69–70HV, substitution N501Y in the RBD of spike is considered the most important VOC mutation [9]. This mutation occurred independently in the South African strain 501Y.V2. We estimated that 501Y.V2 has replicative advantage over other South African strains equal 1.55 [95% CI: 1.43–1.69]. The replicative advantage of 501Y.V2 strains supports the conjecture that mutation N501Y increases infectiousness of SARS-CoV-2 by increasing the affinity of spike RBD to the angiotensin-converting enzyme 2 (ACE2) [24].

Both our estimates suggests that the replicative advantage of the VOC strain is higher than early estimate, 1.47 [95% CI: 1.34–1.59] [9]. Also Davies et al. [11] estimated that the VOC strain has 43–82% [95% CrI: 38–106%] higher transmissibility. In these studies, the authors estimated the replicative advantage of the VOC strain separately for each region and then averaged over analyzed regions. We think that our approach, based on weekly and England-averaged data, gives a more accurate estimate. This is because the VOC strain (as well as other strains) evolves and some substrains have a lower or even no replicative advantage and will become extinct in the course of evolution, and finally a substrain with the highest replicative advantage will dominate. Substrains with a small replicative advantage may contribute to the VOC to non-VOC replication ratio averaged over NHS STP areas of England but only marginally influence the expansion of the VOC strain globally. By using the aggregated approach we estimate the replicative advantage of the dominating substrain(s). To convert the weekly growth of the ratio of genomes to the ratio of respective reproduction numbers we assumed that both strains have the same mean serial interval equal 6.73 days. This calculation is approximate as the serial interval is not a number but follows a hypoexponential distribution [17]. Additionally, although no current data indicate this, it may happen that the faster spread of VOC strain partially results from a shorter generation time.

The VOC strain, because of its high replicative advantage, is likely to become globally dominant, possibly reaching fixation. It will continue to evolve so it is crucial to track which mutations already present in this strain can further increase its replicative advantage. We thus performed detailed sequence analysis of the spike protein identifying, as of February 12, 2021, 1213 different mutations in VOC genomes collected in England with 697 of them found in more than one submitted VOC genome. By systematic analysis of the propagation of VOC substrains we found that substrain(s) conferring L18F substitution is/are the most abundant and rapidly growing VOC substrain(s).

Based on data collected in the five-week period of December 7, 2020–January 17, 2021, in England, we estimated the replicative advantage of this substrain in relation to the remaining VOC strains as 1.72 [95% CrI: 1.57–2.02]. As this estimate is based on a relatively short time interval it must be taken with caution. Importantly, L18F mutation has also expanded in the South African strain 501Y.V2 defined by three spike mutations K417N, E484K, N501Y (thus sharing with the VOC strain spike mutation N501Y). Among the 501Y.V2 genomes collected after December 1, 2020, the L18F substrain constitutes 41% genomes (127 out of 309), according to GISAID as of February 12, 2021). In Brazil, in strain P.1 defined by three spike mutations K417T, E484K, N501Y (differing from the South African strain 501Y.V2 by substitution K417T instead of K417N), mutation L18F has been found in 93% of genomes (69 out of 74) collected after December 1, 2020. This data suggests a replicative advantage of L18F substrains within the VOC, 501Y.V2, and P.1 strains, in, respectively, England, South Africa, and Brazil. This replicative advantage of the VOC L18F substrain must be considered with caution until the mechanism promoting faster spread of strains containing L18F substitution is elucidated. Leucine 18 lies in the N-terminal domain (NTD), that has not been typically considered as a target for neutralizing antibodies. However, there is a growing number of studies showing that the NTD is targeted by antibodies and that NTD deletion 69H–70V (characterizing the VOC strain) compromises binding of antibodies [25–27]. With respect to L18F, an in vitro study by Cele et al. shows that an African variant L18F, D80A, D215G, K417N, E484K, N501Y, D614G, A701V propagates much faster than a variant without L18F mutation in the presence of plasma antibodies collected from donors infected in the first wave of epidemic in South Africa (June–August, 2020) [28]. Correspondingly, McCallum et al. showed that L18F substitution compromises binding of neutralizing antibodies [15]. Findings by Cele et al. and McCallum et al., together with the increase of L18F variants in 501Y.V2, P.1, and VOC strains, suggests that the replicative advantage of L18F mutants can be partly associated with their ability to infect seroprevalent individuals, and thus depend on the fraction of seroprevalent individuals in given territory. In turn, growth of strains with mutations in escape residues L18 and S494 on the VOC strain suggests an increasing selection pressure resulting from the growth of the seroprevalent fraction of the population of England. This trend can be enhanced by the ongoing English vaccination program, in which the relatively large time span between the first and second dose can be a contributing factor.

In summary, we have shown that the new VOC strain has about twofold replicative advantage over the 20A.EU1 strain of SARS-CoV-2, that was dominating in England in November 2020. The strain has already spread across the world and will likely spread further towards fixation. It was present in 40 out of 46 countries that reported at least 50 viral genomes in January, 2021. Spread of the faster-replicating VOC-202012/01 strain may hinder the efforts to contain the COVID-19 epidemics prior to mass vaccinations. As the global spread of the VOC strain is very likely, it is important to monitor mutations of this strain, with particular attention to mutations interfering with immune response including the fast spreading NTD mutation L18F, and RBM mutations E484K, F490S, and S494P that may decrease the efficacy of currently available vaccines.

Data Availability

All data used in this study are referenced and gathered in the supplementary data set S1. GISAID: www.gisaid.org; UK.GOV: https://coronavirus.data.gov.uk/details/cases

Supplementary information

  • Figure S1 (appended at the end of this preprint): The replicative advantage of the VOC strain estimated based on COG-UK pillar 2 data.

  • Data File S1 (provided separately): All datasets that were extracted from public data sources and analyzed in this study.

Author contributions

Conceptualization, Frederic Grabowski, Marek Kochańczyk and Tomasz Lipniacki; Data curation, Frederic Grabowski, Grzegorz Preibisch, Stanisław Giziński and Marek Kochanczyk; Investigation, Frederic Grabowski and Marek Kochańczyk; Software, Grzegorz Preibisch and Stanisław Giziński; Supervision, Marek Kochańczyk and Tomasz Lipniacki; Visualization, Frederic Grabowski, Grzegorz Preibisch, Stanisław Giziński and Marek Kochanczyk; Writing – original draft, Tomasz Lipniacki; Writing – review & editing, Marek Kochańczyk.

Funding

This study was supported by the Norwegian Financial Mechanism GRIEG-1 grant 2019/34/H/ NZ6/00699 (operated by the National Science Centre Poland). The funding agency had no role in study design, data collection, and analysis, decision to publish, or preparation of the manuscript.

Data availability

Publicly available datasets were analyzed in this study. Genome data can be found at gisaid.org (after registration) and cogconsortium.uk, and COVID-19 cases data at ecdc.europa.eu/en/publications-data and coronavirus.data.gov.uk/details/cases. We also provide retrieved and preprocessed data in a single supplementary data file.

Conflicts of interest

The authors declare no conflict of interest.

Supplementary information

Figure S1.
Figure S1. The replicative advantage of VOC over 20A.EU1 in England, based on COG.

Ratio of VOC to 20A.EU1 genomes collected from pillar 2 tests in weeks 38-51 in England. The trend line is fitted to data points from weeks 43–51 (green) and from weeks 43–47 (blue). The weekly growth rate is 1.92 [95% CI: 1.80–2.04] for weeks 43–51 and 2.20 [95% CI: 1.98–2.43] for weeks 43-47. The estimated R ratio is RtVOC / Rt20A.EU1 = 1.87 [95% CI: 1.76-1.98] for weeks 43-51 and 2.13 [95% CI: 1.93-2.35] for weeks 43-47.

Supplementary Data Set S1 (a multi-sheet Excel file) is provied separately.

Acknowledgments

We are very grateful to the GISAID Initiative and all its data contributors, i.e., the Authors from the originating laboratories responsible for obtaining the specimens and the submitting laboratories where genetic sequence data were generated and shared via the GISAID Initiative, on which this research is based.

Footnotes

References

  1. 1.↵
    Rambaut, A., Loman, N. & Pybus, O. Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations Virological. https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563 (2021).
  2. 2.↵
    Shu, Y. & McCauley, J. GISAID: Global initiative on sharing all influenza data – from vision to reality. Eurosurveillance 22, 30494. https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2017.22.13.30494 (2017).
    OpenUrl
  3. 3.↵
    Public Health England. Investigation of novel SARS-CoV-2 variant: 202012/01. Technical Briefing 2 https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/949639/Technical_Briefing_VOC202012-2_Briefing_2_FINAL.pdf (2021).
  4. 4.↵
    Greater London Authority. Coronavirus (COVID-19) Cases London Datastore. https://data.london.gov.uk/dataset/coronavirus--covid-19--cases (2021).
  5. 5.↵
    Hodcroft, E. B. et al. Emergence and spread of a SARS-CoV-2 variant through Europe in the summer of 2020. medRxiv. http://medrxiv.org/content/early/2020/11/27/2020.10.25.20219063.abstract (2020).
  6. 6.↵
    Yurkovetskiy, L. et al. Structural and Functional Analysis of the D614G SARS-CoV-2 Spike Protein Variant. Cell 183, 739–751.e8. ISSN: 0092-8674. https://www.sciencedirect.com/science/article/pii/S0092867420312290 (2020).
    OpenUrl
  7. 7.↵
    Tegally, H. et al. Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa. medRxiv. http://medrxiv.org/content/early/2020/12/22/2020.12.21.20248640.abstract (2020).
  8. 8.↵
    Sabino, E. C. et al. Resurgence of COVID-19 in Manaus, Brazil, despite high seroprevalence. The Lancet 397, 452–455. ISSN: 0140-6736. https://dx.doi.org/10.1016/S0140-6736(21)00183-5 (2021).
    OpenUrl
  9. 9.↵
    Public Health England. Investigation of novel SARS-CoV-2 variant: 202012/01. Technical Briefing 1 https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/947048/Technical_Briefing_VOC_SH_NJL2_SH2.pdf (2021).
  10. 10.↵
    Leung, K., Shum, M. H., Leung, G. M., Lam, T. T. & Wu, J. T. Early transmissibility assessment of the N501Y mutant strains of SARS-CoV-2 in the United Kingdom, October to November 2020. Eurosurveillance 26, 2002106. https://www.eurosurveillance.org/content/10.2807/1560-7917.ES.2020.26.1.2002106 (2021).
    OpenUrl
  11. 11.↵
    Davies, N. G. et al. Estimated transmissibility and severity of novel SARS-CoV-2 Variant of Concern 202012/01 in England. medRxiv. https://www.medrxiv.org/content/early/2021/02/07/2020.12.24.20248822 (2021).
  12. 12.↵
    COVID-19 Genomics UK Consortium https://www.cogconsortium.uk (2021).
  13. 13.↵
    Grabowski, F., Preibisch, G., Kochańczyk, M. & Lipniacki, T. SARS-CoV-2 Variant Under Investigation 202012/01 has more than twofold replicative advantage. medRxiv. https://www.medrxiv.org/content/early/2021/01/04/2020.12.28.20248906 (2021).
  14. 14.↵
    Wibmer, C. K. et al. SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma. bioRxiv. https://www.biorxiv.org/content/early/2021/01/19/2021.01.18.427166 (2021).
  15. 15.↵
    McCallum, M. et al. N-terminal domain antigenic mapping reveals a site of vulnerability for SARS-CoV-2. bioRxiv. https://www.biorxiv.org/content/early/2021/01/14/2021.01.14.426475 (2021).
  16. 16.↵
    Rice, P., Longden, I. & Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends in Genetics 16, 276–277. https://dx.doi.org/10.1016/S0168-9525(00)02024-2 (2000).
    OpenUrlCrossRefPubMedWeb of Science
  17. 17.↵
    Kochańczyk, M., Grabowski, F. & Lipniacki, T. Super-spreading events initiated the exponential growth phase of COVID-19 with R0 higher than initially estimated. Royal Society Open Science 7, 200786. https://royalsocietypublishing.org/doi/10.1098/rsos.200786 (2020).
    OpenUrlCrossRef
  18. 18.↵
    De Maio, N. et al. Issues with SARS-CoV-2 sequencing data Virological. https://virological.org/t/issues-with-sars-cov-2-sequencing-data (2021).
  19. 19.↵
    Qiao, B. & Olvera de la Cruz, M. Enhanced Binding of SARS-CoV-2 Spike Protein to Receptor by Distal Polybasic Cleavage Sites. ACS Nano 14. PMID: 32806067, 10616–10623. https://doi.org/10.1021/acsnano.0c04798 (2020).
  20. 20.↵
    Buss, L. F. et al. Three-quarters attack rate of SARS-CoV-2 in the Brazilian Amazon during a largely unmitigated epidemic. Science 371, 288–292. https://dx.doi.org/10.1126/science.abe9728 (2021).
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    Chakraborty, S. Evolutionary and structural analysis elucidates mutations on SARS-CoV2 spike protein with altered human ACE2 binding affinity. Biochemical and Biophysical Research Communications 534, 374–380. http://www.sciencedirect.com/science/article/pii/S0006291X20321203 (2021).
    OpenUrl
  22. 22.↵
    Koenig, P.-A. et al. Structure-guided multivalent nanobodies block SARS-CoV-2 infection and suppress mutational escape. Science 371, eabe6230. ISSN: 0036-8075.]https://science.sciencemag.org/content/early/2021/01/11/science.abe6230 (2021).
    OpenUrlAbstract/FREE Full Text
  23. 23.↵
    Chen, S., Flegg, J. A., White, L. J. & Aguas, R. Levels of SARS-CoV-2 population exposure are considerably higher than suggested by seroprevalence surveys. medRxiv. https://www.medrxiv.org/content/early/2021/01/09/2021.01.08.21249432 (2021).
  24. 24.↵
    Fratev, F. The N501Y and K417N mutations in the spike protein of SARS-CoV-2 alter the interactions with both hACE2 and human derived antibody: A Free energy of perturbation study. bioRxiv. http://biorxiv.org/content/early/2020/12/31/2020.12.23.424283.abstract (2020).
  25. 25.↵
    Collier, D. A. et al. SARS-CoV-2 B.1.1.7 escape from mRNA vaccine-elicited neutralizing antibodies. medRxiv. https://www.medrxiv.org/content/early/2021/02/03/2021.01.19.21249840 (2021).
  26. 26.
    Kemp, S. A. et al. Neutralising antibodies in Spike mediated SARS-CoV-2 adaptation. medRxiv. https://dx.doi.org/10.1101/2020.12.05.20241927 (2020).
  27. 27.↵
    Voss, W. N. et al. Prevalent, protective, and convergent IgG recognition of SARS-CoV-2 non-RBD spike epitopes in COVID-19 convalescent plasma. bioRxiv. https://www.biorxiv.org/content/early/2020/12/21/2020.12.20.423708 (2020).
  28. 28.↵
    Cele, S. et al. Escape of SARS-CoV-2 501Y.V2 variants from neutralization by convalescent plasma. medRxiv. https://www.medrxiv.org/content/early/2021/01/26/2021.01.26.21250224 (2021).
Back to top
Posted February 21, 2021.
Download PDF

Supplementary Material

Data/Code
Email
SARS-CoV-2 Variant of Concern 202012/01 has about twofold replicative advantage and acquires concerning mutations
Frederic Grabowski, Grzegorz Preibisch, Stanisław Giziński, Marek Kochańczyk, Tomasz Lipniacki
medRxiv 2020.12.28.20248906; doi: https://doi.org/10.1101/2020.12.28.20248906
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area