Next Article in Journal
Retrotransposons as Drivers of Mammalian Brain Evolution
Next Article in Special Issue
Screening, Diagnostic and Prognostic Tests for COVID-19: A Comprehensive Review
Previous Article in Journal
Implantation of the Bonebridge BCI 602 after Mastoid Obliteration with S53P4 Bioactive Glass: A Safe Method of Treating Difficult Anatomical Conditions-Preliminary Results
Previous Article in Special Issue
Self-Reported Real-World Safety and Reactogenicity of COVID-19 Vaccines: A Vaccine Recipient Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Temporal Dominance of B.1.1.7 over B.1.354 SARS-CoV-2 Variant: A Hypothesis Based on Areas of Variant Co-Circulation

by
Evangelia Georgia Kostaki
1,
Ioulia Tseti
2,
Sotirios Tsiodras
3,
George N. Pavlakis
4,
Petros P. Sfikakis
5 and
Dimitrios Paraskevis
1,*
1
Department of Hygiene, Epidemiology and Medical Statistics, Medical School, National and Kapodistrian University of Athens, 11527 Athens, Greece
2
Uni-Pharma S.A., 14564 Kifissia, Greece
3
4th Department of Internal Medicine, Attikon University Hospital, Medical School, National and Kapodistrian University of Athens, 12462 Haidari, Greece
4
Human Retrovirus Section, National Cancer Institute, Frederick, MD 21702, USA
5
1st Department of Prop. Internal Medicine, Medical School, National and Kapodistrian University of Athens, 11527 Athens, Greece
*
Author to whom correspondence should be addressed.
Life 2021, 11(5), 375; https://doi.org/10.3390/life11050375
Submission received: 5 April 2021 / Revised: 17 April 2021 / Accepted: 19 April 2021 / Published: 22 April 2021
(This article belongs to the Special Issue Ecology, Evolution and Epidemiology of Coronaviruses)

Abstract

:
Some emergent SARS-CoV-2 variants raise concerns due to their altered biological properties. For both B.1.1.7 and B.1351 variants, named as variants of concern (VOC), increased transmissibility was reported, whereas B.1.351 was more resistant to multiple monoclonal antibodies (mAbs), as well as convalescent and vaccination sera. To test this hypothesis, we examined the proportion of VOC over time across different geographic areas where the two VOC, B.1.1.7 and B.1.351, co-circulate. Our comparative analysis was based on the number of SARS-CoV-2 sequences on GISAID database. We report that B.1.1.7 dominates over B.1.351 in geographic areas where both variants co-circulate and the B.1.1.7 was the first variant introduced in the population. The only areas where B.1.351 was detected at higher proportion were South Africa and Mayotte in Africa, where this strain was associated with increased community transmission before the detection of B.1.1.7. The dominance of B.1.1.7 over B.1.351 could be important since B.1.351 was more resistant to certain mAbs, as well as heterologous convalescent and vaccination sera, thus suggesting that it may be transmitted more effectively in people with pre-existing immunity to other VOC. This scenario would lessen the effectiveness of vaccine and urge the need to update them with new strains.

1. Introduction

SARS-CoV-2 has caused a devastating pandemic with serious consequences in global health and economy. Since the first characterization of the SARS-CoV-2 genome, a large collection of sequences (>770 K until the middle of March 2021) has been submitted to the GISAID database (http://www.GISAID.org, accessed on 11 March 2021), thus facilitating the monitoring of global dispersal of the virus as well as the identification of lineages and variants termed as variants of interest (VOI) and variants of concern (VOC). According to the WHO, the definition of “SARS-CoV-2 Variant of Interest” (VOI) is an isolate phenotypically changed compared to a reference isolate or has a genome with mutations that lead to amino acid changes associated with established or suspected phenotypic implications and has been identified to cause community transmission/multiple COVID-19 cases/clusters, or has been detected in multiple countries; or is otherwise assessed to be a VOI by WHO in consultation with the WHO SARS-CoV-2 Virus Evolution Working Group [1]. Phenotypic changes include changes in the epidemiology, antigenicity, or virulence or changes that have or potentially have a negative impact on available diagnostics, vaccines, therapeutics, or public health and social measures. WHO will provide guidance on amino acid changes with established or suspected phenotypic implications, and may be informed by a database on key amino acid changes, or as reported in the scientific literature.
A VOI (as defined previously) is named a variant of concern (VOC) if, through a comparative assessment, it has been demonstrated to be associated with: (i) increase in transmissibility or detrimental change in COVID-19 epidemiology; (ii) increase in virulence or change in clinical disease presentation; or decrease in effectiveness of public health and social measures or available diagnostics, vaccines, or therapeutics; or assessed to be a VOC by WHO in consultation with the WHO SARS-CoV-2 Virus Evolution Working Group [1].
The earliest interest was focused on the variant harboring the D614G mutation in the spike protein as well as other linked mutations, was reported early on during the course of the pandemic in late January 2020 and was subsequently named C77 (Figure 1). Preliminary experimental evidence was suggestive of increased infectivity compared to the initial virus identified in China [2,3]. In December 2020, the UK authorities announced a new variant named B.1.1.7 or 501Y.V1 responsible for a virus surge across the UK (Figure 1) [4]. At the same time, a second VOC named B.1.351 or 501Y.V2 was announced to rapidly spread in different provinces of South Africa (Figure 1) [5]. The third variant (501Y.V3 or P.1) (Figure 1) was associated with an increased number of cases in Manaus, Brazil [6], a place hit hard by previous pandemic waves with approximately three quarters of its population reportedly developing immunity to SARS-CoV-2. Notably, all these variants embed several amino acid replacements in the spike including some key residues in the receptor binding site (RBD) or other sites. Specifically, all VOC include N501Y in RBD and D614G, while B.1.351 and P.1 have also the E484K in RBD region (Figure 1). B.1.1.7 has additionally two deletions in spike that were unique compared to the other variants. Additional variants with different combinations of mutations have been described such as CAL.20.C (B.1.429), B.1.525 and B.1.526 from southern California [7], Nigeria, and New York [8], respectively, or the B.1.1.7 with the addition of E484K (Figure 1).
These clades have raised concern due to their increased genetic variability and the fact that in preliminary observations they have been associated with virus surges in different geographical areas. Previous studies report that the N501Y mutation in the virus RBD may be associated with increased affinity of the spike protein with the ACE2 receptor in human cells, and therefore with higher transmissibility. In one report, affinity increased by ~3.5-fold [9], whereas, by the addition of the E484K to N501Y, a higher binding affinity (~12.7-fold increase) was documented [10]. In addition, the mean duration of B.1.1.7 associated infection was 13.3 days versus 8.2 days for non-B.1.1.7 viruses. No differences were documented for the peak viral concentrations of 8.5 log10 RNA copies/mL and 8.2 log10 RNA copies/mL for B.1.1.7 and non-B.1.1.7, respectively [11].
Effects on therapeutic efficacy of anti-viral interventions and/or immunization has been the main concern with the emerging variants. A recent study by Wang et al., suggested that B.1.1.7 was refractory to monoclonal antibodies (mAb) against the N-terminal of the spike and to a few mAbs to the RBD [12], but maintained susceptibility to convalescent plasma or vaccine induced neutralization. On the other hand, the E484K mutation (present in several emerging variants) has been reported to greatly affect the neutralization activity of monoclonal antibodies and to reduce convalescent plasma or vaccine induced immunity, suggesting that strains with this mutation can probably re-infect people who are vaccinated or have been previously infected with SARS-CoV-2 more easily than other viruses [13,14]. Similar effects have been reported for other mutations in the spike, such as the K417N [12], the K417T in B.1.351 and P.1, and the L18F mutations. The B.1.351 variant is reportedly more refractory to multiple mAbs from several individuals and also more resistant to convalescent plasma and vaccination sera due to the presence of the E484K mutation [12]. This finding was further supported by the diminished efficacy of Novavax Johnson & Johnson and ChAdOx1 vaccines in South Africa where the B.1.351 predominates [15,16,17].
Although the ability to escape pre-existing immunity provides the most reasonable explanation for the fitness advantage of the B.1.351 and P.1 viruses, it does not explain the replicative advantage of B.1.1.7. For this VOC, improved infectivity rather than immune escape appears to contribute to increased fitness and transmission dynamics as detected in the UK, Portugal, and elsewhere.
To the best of our knowledge, potential differences in the replicative advantage between VOC are largely unknown. Given the characteristics and the potential differences in the biological mechanism conferring improved fitness of the VOC, it is important to test the hypothesis of a selective advantage of any of the VOC in settings where the two variants co-circulate. Therefore, we aimed to analyze the proportion of VOC over time in different geographic areas where at least two of the most widely VOC, namely, B.1.1.7 and B.1.351, were present. Selection of these variants was based on the fact that they have been the most widely spread and the ones that co-circulate across different geographic areas.

2. Materials and Methods

The analysis was based on the proportions estimated on SARS-CoV-2 sequences available on GISAID database. Specifically, the number of available genomic sequences for B.1.1.7 and B.1.351 was estimated on 11 March 2021, by using the outbreak.info tool (https://outbreak.info, accessed on 11 March 2021). Only countries with more than 50 sequences were included in the analysis. Proportions were estimated using as reference the total number of sequences generated since the first identification of the corresponding variant in each country. The time difference in the date of the first identification of B.1.1.7 versus the B.1.351 was also estimated with positive values corresponding to earlier detection of B.1.1.7.
At the stage of manuscript revision, the number of available genomes belonging to B.1.1.7 and B.1.351 was assessed on 15 April 2021 (https://outbreak.info). Additionally, the number of B.1.1.7 and B.1.351 genomic sequences for different geographic regions was estimated for the time during 15 February and 14 March 2021 based on the data available on GISAID database (https://www.gisaid.org/, accessed on 15 April 2021).
Statistical analysis was performed by fitting a quartile (median) regression model on the difference between the proportions of B.1.1.7 and B.1.351 and the time difference between the earliest variant sequence for B.1.1.7 versus B.1.351. In addition, a multinomial logistic regression model was fitted on 202,005 observations corresponding to sequences from 30 countries published on GISAID database between 15 February and 14 March 2021. We used data only from the last month of our study’s observational period for this model to assure that there was enough time since the B.1.1.7 and B.1.351 introduction to spread in the local population of each country. Variant was defined as the outcome variable on this model, while time difference between the earliest variant sequence for B.1.1.7 versus B.1.351 was chosen as possible explanatory variable. The level of significance was set at 0.05. Both analyses were performed in Stata 13-StataCorp LP software.

3. Results

Analysis of available data showed that in the vast majority of the countries (31 out of 33), the proportion of B.1.1.7 cases was higher than the corresponding proportion of B.1.351 (Figure 2a). The only areas where B.1.351 dominated were South Africa and Mayotte, an overseas department of France located in the Indian Ocean off the coast of Southeast Africa. South Africa was the place where the B.1.351 was first identified and associated with a virus surge in early December 2020 [5]. Countries where both variants were detected at high proportions were New Zealand, Austria, Singapore, and the United Arab Emirates. New Zealand experienced a first epidemic wave between March and May 2020. Since then, virus spread has been successfully controlled with an approximate number of active cases equal to 93; the 7-day average number of cases ranges between 1 and 11 cases. Similarly, in Singapore the 7-day average number of cases was approximately between 10 and 31 cases after September 2020, suggesting that no virus surges were detected after the introduction of the variants in these countries. In the United Arab Emirates, a virus surge was detected at the beginning of the new year and since February 2021 the number of cases has been declining. The number of sequenced viruses from the United Arab Emirates was limited to 148 and 64 for B1.1.7 and B.1.351, respectively. This is only a small proportion compared to the total number of diagnosed cases in the same time period. In Austria, both variants have been detected and they have been mostly concentrated in the Tyrol state. Although the B.1.1.7 identified first in the UK, where it has been dominating at proportions higher than 75% since the beginning of the current year, it has surpassed B.1.351 across different regions where both variants co-circulate. The proportion of available genomes for the two variants was re-estimated for the time after 29 March 2021 and until 11 April 2021. A similar patten was observed about the dominance of B.1.1.7 versus B.1.251 (Table S1).
To investigate whether these differences were due to the longer presence of B.1.1.7 in remote geographic areas, we plotted the time lag since the first identification of the two strains (earliest sequence available per clade) (Figure 2b). More specifically, the time lag corresponds to the time difference between the date of the earliest variant sequence for each geographic region. The differences were more pronounced for the countries where both variants identified first, i.e., for the UK and South Africa. Except for these areas, the median time difference was 19 days (IQR: 10, 40).
Measurements on the proportion of the B.1.1.7 were available for some countries using different methodologies. Some used as a proxy the S gene target failure (SGTF) in real-time PCR assays, DNA sequencing or a combination of the two methods. The proportions shown in Table 1 were based on the number of cases diagnosed during a period of one week (week 7 until 9 of 2021). The countries were shown according to the proportion of the B.1.1.7 variant in order for the corresponding proportions of the B.1.351 to be comparable. The figures in Table 1, as expected, were different from the cumulative proportions of B.1.1.7 since the first day of the variant identification. For all countries, apart from Poland, the SGTF proportions were higher than the cumulative proportions, suggesting an increasing trend of the B.1.1.7 in the corresponding countries. Similar figures for the B.1.351 were not available.
To further investigate if the difference in the first detection of each variant was associated with the difference between the variants’ proportions, we performed statistical analysis by fitting a quartile (median) regression model to the data. Analysis revealed that for a one day increase in the time difference between the earliest variant sequence for B.1.1.7 versus B.1.351 (time lag), the predicted value of change in the difference between the proportions of B.1.1.7 and B.1.351 increases by 0.4 (95% CI: 0.1–0.8; p = 0.021).
Furthermore, the multinomial logistic regression analysis revealed that if the time lag were increased by one day, the relative risk for B.1.1.7 compared to B.1.351 would be expected to increase by a factor of 1.07 (95% CI: 1.02–1.11; p = 0.002). Particularly, a one day increase in the time lag increases the relative risk of Β.1.1.7 (compared to Β.1.351) by about 7%.

4. Discussion

Our analysis suggests that in countries where both B.1.1.7 and B.1.351 co-circulate, B.1.1.7 has dominated local transmission patterns, except for the South Africa where B.1.351 variant has been associated with community transmission across different provinces before the introduction of B.1.1.7 [5]. Statistical analysis confirmed our findings and specifically the multinomial logistic regression analysis revealed that a one day increase in the time lag increases the relative risk of Β.1.1.7 (compared to Β.1.351) by about 7%. We also found that B.1.1.7 was introduced before the South African lineage in Europe, Asia and the Americas, but in several areas the time difference since the first identification of the two strains was not pronounced. The transmission pattern of the two variants could also be explained due to the limited importation events in Europe from remote South Africa compared to the UK. On the other hand, as soon as B.1.351 has entered Europe as early as 10 December 2020 (Table 1), it could have spread in Europe through internal flights.
The selection advantage of different viruses is a dynamic process depending on the complex interplay between the pathogen and the host suggesting that different strains can dominate in different environments [18,19]. Selection of B.1.1.7 variant and dominance over the pre-existing viruses was probably due to its increased transmissibility rather than its ability to re-infect people with existing immunity [12]; however, this may change in the future when vaccination coverage increases. Selection advantage therefore cannot be determined across different populations and time points since it depends on several factors such as the levels of the existing immunity, the characteristics of the strains for which immunity was developed, climatic factors, the speed of the vaccination programs, compliance to public health measures, as well as other factors, importantly, the virus mutation rate, and its ability to recombine. In this study, we show that, under the conditions of the pandemic since the late 2020 and before we reach high vaccination levels, the B.1.351 virus does not dominate over B.1.1.7 in geographic areas where both variants co-circulate, and the B.1.1.7 was the first variant introduced in the population. However, our findings cannot exclude the risk for the B.1.351 to cause outbreaks due to its ability to escape immunity to other variants.
One major limitation of our study is that the analyzed viral genomes may not be representative of the circulated viruses in each region. Our analysis was based on sequencing data and information on the weekly proportions of the B.1.1.7 which confirmed that this variant was more prevalent than the B.1.351 across different geographic regions; thus, suggesting that B.1.1.7 dominates B.1.351 in areas where both variants co-circulate and the former variant was the first associated with community transmission. Another limitation was that the earliest variant date did not accurately reflect the day zero of the virus introduction in a particular area, but it provides a proxy of that date. Our study suggests that that the B.1.1.7 variant appears to be the only one able to dominate over the B.1.351 variant. Moreover, the minimum proportion or the minimum incidence of the B.1.1.7 variant necessary to keep the B.1.351 variant under control cannot be estimated from the current data.
The dominance of B.1.1.7 over B.1.351 could be of importance since the latter is more resistant to certain mAbs, convalescent plasma or vaccination sera [12], thus suggesting that it may be transmitted more effectively in people with pre-existing immunity. This scenario would lessen the effectiveness of vaccine and urge the need to update them with new strains. Based on the existing data, our study highlights that B.1.351 may have similar fitness or a disadvantage compared to B.1.1.7 in areas where it is associated with community transmission, therefore raising hope that a more efficient “immune-evade” VOC cannot cause a widespread dispersal over B.1.1.7. The future transmission dynamics cannot be easily forecasted, but the current message is that the transmission advantage of B.1.1.7 is probably higher than B.1.351, and this characteristic can protect us from the worldwide dispersal of a more dangerous variant for now. It remains to be seen whether a subsequent epidemic can be triggered by B.1.351 once the B.1.1.7 epidemic has subsided. In addition, the possibility of recombination generating a virus combining the better transmission of B.1.1.7 with the more severe pathologies of B.1.351 has to be considered and may shape the direction of virus evolution. Furthermore, this may create additional opportunities for new virus transmissions.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life11050375/s1, Table S1: Information related to the number of variant sequences, the earliest variant sequence, the last detection date, and the corresponding proportions for B.1.1.7 and B.1.351 variants of concern (VOC).

Author Contributions

Conceptualization, G.N.P., P.P.S. and D.P.; methodology, E.G.K. and D.P.; formal analysis, E.G.K., I.T., S.T., G.N.P., P.P.S. and D.P.; writing—original draft preparation, E.G.K. and D.P.; writing—review and editing, E.G.K., I.T., S.T., G.N.P., P.P.S. and D.P.; visualization, E.G.K.; supervision, D.P.; funding acquisition, P.P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Kleon Tsetis Foundation, in the context of promoting academic research initiatives within its COVID-19 response framework.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available data were downloaded from www.gisaid.org (accessed on 11 March 2021).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. World Health Organization. Weekly Epidemiological Update—23 February 2021 Special Edition: Proposed Working Definitions of SARS-CoV-2 Variants of Interestand Variants of Concern; WHO: Geneva, Switzerland, 2021. [Google Scholar]
  2. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B.; et al. Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus. Cell 2020, 182, 812–827.e19. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, L.; Jackson, C.B.; Mou, H.; Ojha, A.; Peng, H.; Quinlan, B.D.; Rangarajan, E.S.; Pan, A.; Vanderheiden, A.; Suthar, M.S.; et al. SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity. Nat. Commun. 2020, 11, 6013. [Google Scholar] [CrossRef] [PubMed]
  4. Rambaut, A.; Loman, N.; Pybus, O.; Barclay, W.; Barrett, J.; Carabelli, A.; Connor, T.; Peacock, T.; Robertson, D.L.; Volz, E. Preliminary Genomic Characterisation of an Emergent SARS-CoV-2 Lineage in the UK Defined by a Novel Set of Spike Mutations. Available online: 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 (accessed on 2 April 2021).
  5. Tegally, H.; Wilkinson, E.; Giovanetti, M.; Iranzadeh, A.; Fonseca, V.; Giandhari, J.; Doolabh, D.; Pillay, S.; San, E.J.; Msomi, N.; 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 2020. [Google Scholar] [CrossRef]
  6. Faria, N.R.; Claro, I.M.; Candido, D.; Franco, L.A.M.; Andrade, P.S.; Coletti, T.M.; Silva, C.A.M.; Sales, F.C.; Manuli, E.R.; Aguiar, R.S.; et al. Genomic Characterisation of an Emergent SARS-CoV-2 Lineage in Manaus: Preliminary Findings. Available online: https://virological.org/t/genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-manaus-preliminary-findings/586 (accessed on 2 April 2021).
  7. Zhang, W.; Davis, B.D.; Chen, S.S.; Martinez, J.M.S.; Plummer, J.T.; Vail, E. Emergence of a Novel SARS-CoV-2 Variant in Southern California. JAMA 2021, 325, 1324–1326. [Google Scholar] [CrossRef] [PubMed]
  8. Lasek-Nesselquist, E.; Lapierre, P.; Schneider, E.; George, K.S.; Pata, J. The localized rise of a B.1.526 SARS-CoV-2 variant containing an E484K mutation in New York State. medRxiv 2021. [Google Scholar] [CrossRef]
  9. Starr, T.N.; Greaney, A.J.; Hilton, S.K.; Ellis, D.; Crawford, K.H.D.; Dingens, A.S.; Navarro, M.J.; Bowen, J.E.; Tortorici, M.A.; Walls, A.C.; et al. Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding. Cell 2020, 182, 1295–1310.e20. [Google Scholar] [CrossRef] [PubMed]
  10. Nelson, G.; Buzko, O.; Spilman, P.; Niazi, K.; Rabizadeh, S.; Soon-Shiong, P. Molecular dynamic simulation reveals E484K mutation enhances spike RBD-ACE2 affinity and the combination of E484K, K417N and N501Y mutations (501Y.V2 variant) induces conformational change greater than N501Y mutant alone, potentially resulting in an escape mutant. bioRxiv 2021. [Google Scholar] [CrossRef]
  11. Kissler, S.; Fauver, J.R.; Mack, C.; Tai, C.G.; Breban, M.I.; Watkins, A.E.; Samant, R.M.; Anderson, D.J.; Ho, D.D.; Grubaugh, N.D.; et al. Densely sampled viral trajectories suggest longer duration of acute infection with B.1.1.7 variant relative to non-B.1.1.7 SARS-CoV-2. Harvard University’s DASH repository 2021. medRxiv 2021. [Google Scholar] [CrossRef]
  12. Wang, P.; Nair, M.S.; Liu, L.; Iketani, S.; Luo, Y.; Guo, Y.; Wang, M.; Yu, J.; Zhang, B.; Kwong, P.D.; et al. Antibody Resistance of SARS-CoV-2 Variants B.1.351 and B.1.1.7. Nature 2021, 2021, 1–6. [Google Scholar] [CrossRef]
  13. Cele, S.; Gazy, I.; Jackson, L.; Hwa, S.-H.; Tegally, H.; Lustig, G.; Giandhari, J.; Pillay, S.; Wilkinson, E.; Naidoo, Y.; et al. Escape of SARS-CoV-2 501Y.V2 from neutralization by convalescent plasma. Nature 2021, 2021, 1–9. [Google Scholar] [CrossRef]
  14. Wibmer, C.K.; Ayres, F.; Hermanus, T.; Madzivhandila, M.; Kgagudi, P.; Oosthuysen, B.; Lambson, B.E.; de Oliveira, T.; Vermeulen, M.; van der Berg, K.; et al. SARS-CoV-2 501Y.V2 escapes neutralization by South African COVID-19 donor plasma. Nat. Med. 2021, 27, 622–625. [Google Scholar] [CrossRef] [PubMed]
  15. Wadman, M.; Cohen, J. Novavax vaccine delivers 89% efficacy against COVID-19 in U.K.—But is less potent in South Africa. Science 2021. [Google Scholar] [CrossRef]
  16. Callaway, E.; Mallapaty, S. Novavax offers first evidence that COVID vaccines protect people against variants. Nature 2021, 590, 17. [Google Scholar]
  17. Madhi, S.A.; Baillie, V.; Cutland, C.L.; Voysey, M.; Koen, A.L.; Fairlie, L.; Padayachee, S.D.; Dheda, K.; Barnabas, S.L.; Bhorat, Q.E.; et al. Efficacy of the ChAdOx1 nCoV-19 Covid-19 Vaccine against the B.1.351 Variant. N. Engl. J. Med. 2021. [Google Scholar] [CrossRef] [PubMed]
  18. Baele, G.; Suchard, M.A.; Rambaut, A.; Lemey, P. Emerging Concepts of Data Integration in Pathogen Phylodynamics. Syst. Biol. 2017, 66, e47–e65. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Cobey, S. Pathogen evolution and the immunological niche. Ann. N. Y. Acad. Sci. 2014, 1320, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Amino acid mutations and deletions in spike protein for the different variants named after variants of concern (VOC). The different domains of the spike and their length are shown in boxes at the upper part of the figure. Countries are represented by ISO Alpha-2 codes (BR: Brazil, GB: United Kingdom, NG: Nigeria, US: United States, ZA: South Africa).
Figure 1. Amino acid mutations and deletions in spike protein for the different variants named after variants of concern (VOC). The different domains of the spike and their length are shown in boxes at the upper part of the figure. Countries are represented by ISO Alpha-2 codes (BR: Brazil, GB: United Kingdom, NG: Nigeria, US: United States, ZA: South Africa).
Life 11 00375 g001
Figure 2. (a) Proportion of B.1.1.7 and B.1.351 in different countries based on the number of different genomes available in the GISAID database; (b) Time difference between the earliest variant sequence for B.1.1.7 versus the B.1.351 variants of concern (VOC).
Figure 2. (a) Proportion of B.1.1.7 and B.1.351 in different countries based on the number of different genomes available in the GISAID database; (b) Time difference between the earliest variant sequence for B.1.1.7 versus the B.1.351 variants of concern (VOC).
Life 11 00375 g002
Table 1. Information related to the number of sequences, the earliest variant sequence and the corresponding proportions for B.1.1.7 and B.1.351 variants of concern (VOC).
Table 1. Information related to the number of sequences, the earliest variant sequence and the corresponding proportions for B.1.1.7 and B.1.351 variants of concern (VOC).
VOCB.1.1.7B.1.351
CountriesProportion (%)Proportion of SGTF 1—Week 8 of 2021 (%)Earliest Variant SequenceNumber of Variant SequencesNumber of Sequences 2Proportion (%)Earliest Variant SequenceNumber of Variant SequencesNumber of Sequences 2
Ireland74.5088.617/12/2020196626390.6222/12/2020162594
United Kingdom51.99NA 320/09/2020111,140213,7610.1610/12/2020212135,667
Ghana54.03NA10/12/2020671244.2106/01/2021495
Finland47.43NA18/12/20202685651.0719/12/20206563
Italy43.67NA14/12/2020224451390.2130/01/202183788
France42.7565.813/12/2020184743203.4822/12/20201444137
Slovenia40.00NA29/12/2020461153.6409/02/2021255
Belgium39.3346.330/11/2020206452485.0220/12/20202444857
Spain33.5325–3008/11/2020132339460.3424/12/2020102930
Sweden32.3941.5 420/12/202045814141.2224/12/2020161310
New Zealand48.09NA16/12/20206313116.9629/12/202019112
Portugal30.1750.509/11/202074624730.5104/01/202191782
Germany30.4454.530/11/2020442714,5431.0321/12/202014313,817
Israel28.74~90 516/12/202043415101.9231/12/202016833
Turkey29.15NA24/12/202047916433.7822/01/2021541428
Norway28.4972.509/12/202033711833.9227/12/2020401021
Netherlands22.0064.3 512/11/2020169777151.6322/12/2020996091
Croatia22.06NA20/01/2021602722.1309/02/2021294
Australia22.52NA30/11/20201345954.1410/12/202023556
Poland20.08922/12/20201517522.1110/02/2021295
Austria31.4963.222/12/2020336106715.0923/12/20201591054
South Korea16.45NA14/12/2020905470.9726/12/20204413
Switzerland15.8040.5 509/11/2020198312,5500.6312/11/20207712,313
Denmark12.4776.509/11/2020488939,1910.0604/01/20211219,655
Singapore19.70NA08/12/2020663358.4507/02/2021671
United Arab Emirates14.19NA16/11/2020211487.8126/12/2020564
Canada4.83NA15/12/20205411170.2625/12/20202761
Luxembourg8.4065.524/12/2020323814.1716/01/2021248
United States of America3.5826.217/12/2020265274,1290.0601/01/20213862,390
Thailand6.06NA08/01/202181323.6103/02/2021383
Japan1.07NA01/12/20205955270.2619/12/202093441
Mayotte0.19NA13/01/2021151852.2407/01/2021338647
South Africa0.47NA09/01/2021121365.7408/10/202010861652
1 SGTF: S gene target failure. 2 This number corresponds to the total number of characterized sequences irrespective of their classification since the identification of the corresponding variant. The numbers are different for the two variants due to the difference in dates of earliest variant sequence. 3 NA: Not available. 4 Based on sequencing of all variants with N501Y + A570D. 5 Data available for week 7.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kostaki, E.G.; Tseti, I.; Tsiodras, S.; Pavlakis, G.N.; Sfikakis, P.P.; Paraskevis, D. Temporal Dominance of B.1.1.7 over B.1.354 SARS-CoV-2 Variant: A Hypothesis Based on Areas of Variant Co-Circulation. Life 2021, 11, 375. https://doi.org/10.3390/life11050375

AMA Style

Kostaki EG, Tseti I, Tsiodras S, Pavlakis GN, Sfikakis PP, Paraskevis D. Temporal Dominance of B.1.1.7 over B.1.354 SARS-CoV-2 Variant: A Hypothesis Based on Areas of Variant Co-Circulation. Life. 2021; 11(5):375. https://doi.org/10.3390/life11050375

Chicago/Turabian Style

Kostaki, Evangelia Georgia, Ioulia Tseti, Sotirios Tsiodras, George N. Pavlakis, Petros P. Sfikakis, and Dimitrios Paraskevis. 2021. "Temporal Dominance of B.1.1.7 over B.1.354 SARS-CoV-2 Variant: A Hypothesis Based on Areas of Variant Co-Circulation" Life 11, no. 5: 375. https://doi.org/10.3390/life11050375

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop