COVID-19 pandemic has been the first of its kind caused by coronavirus. Since its outbreak in December 2019, it has led to 207,784,507 confirmed cases and 4,370,424 deaths worldwide (data accessed on 18 August 2021) [11]. There have been ups and downs of the situation, and in 2021 a number of countries have faced fresh COVID-19 resurgence, including but not limited to the 6 countries in this study. In addition to social factors like public health policy relaxation and pandemic fatigue, emergence of vaccine-escape variants was also an important contributing factor. Our data showed that in July, delta VOC has become the predominant lineage in Africa, Asia, Europe, North America and Oceania. In South America, its proportion has been increasing (23.24%) and approached that of the predominant yet diminishing gamma VOC (56.87%), which is not surprising as delta VOC was reported to be 40–60% more transmissible than alpha VOC [12].
From country-level temporal data, we might appreciate the expansion-diminishment patterns of various variants. For instance, the mean expansion-diminishment duration of alpha VOC was about 7.5 months, and the ascending and descending phases lasted for 4.75 and 2.75 months in average, respectively. In addition, another interesting observation was that the ascending phases of new cases were usually accompanied with rise and predominance of alpha (except the USA), beta, gamma, zeta, and most recently delta variants, or right after the time points of major lineage interchange. This might be (or might not be) accompanied with a rise in monthly new deaths, which appeared to be associated with immunization coverage. For instance, proportion of delta VOC in India rose sharply from March to April, surpassing alpha VOC and kappa VOI. In the meantime, both new cases and deaths increased drastically by 26.90 and 43.66 folds in May, with only 12% of Indian population vaccinated. Similar phenomenon was observed for the under-vaccinated population in South Africa. For the USA, the UK and Israel, more than half of their populations were vaccinated when delta VOC grew drastically to predominance. Their death tolls remained stable despite surges of new cases. Besides immunization coverage, availability of medical resources might also contribute to difference in mortality [13].
As of 1 June 2021, WHO validated 7 COVID-19 vaccines for emergency use, including 2 inactivated (Sinopharm and CoronaVac by Sinovac) and 5 RNA-based vaccines (BNT162b2 by Pfizer-BioNTech, AZD1222 by AstraZeneca (SK Bioscience, Serum Institute of India and European Union), and mRNA-1273 by Moderna) [14]. Emergence of SARS-CoV-2 variants has raised concern over vaccine efficacy against growing variety of lineages. A number of research groups reported vaccine-breakthrough cases involving SARS-CoV-2 variants and nearly all validated vaccines, for instance, gamma VOC against CoronaVac and AZD1222 in Brazil [15, 16]; 2 spike variants harbouring E484K and triple mutations (T95I, del142-144, and D614G) against mRNA-1273 and BNT162b2, respectively in the USA [17]; alpha and delta VOCs against AZD1222 and BNT162b2 in England [18, 19] and Scotland [20]; as well as delta VOC infecting an Everest trekker vaccinated with mRNA-1273 in Nepal [21]. Despite a broad range of cycle threshold values (15-33.3), the authors reported that all these vaccine-breakthrough cases involved reduced disease severity.
On the other hand, a number of studies focused on neutralizing activity of post-second-dose sera. For CoronaVac, Chen and coworkers reported that the neutralizing efficiency against pseudotyped lentiviruses of beta and gamma VOCs, and iota VOI were significantly reduced by factors of 5.27, 3.92 and 4.03, respectively, whereas efficiency against D614G, alpha VOC and epsilon variant were equally effective compared with the wildtype [22]. Their data also revealed that 5–34% of sera were capable of neutralizing the former 3 variants, compared to 82% against wildtype. Another study led by Cao and coworkers revealed that beta VOC pseudovirus or the authentic virus caused a major reduction in neutralization [23]. For BNT162b2, Lustig and coworkers revealed that neutralizing titers against alpha, delta and gamma VOCs were reduced by 1.7–2.6 folds, and 10.4 folds for beta VOC [24]. For AZD1222, Madhi and coworkers performed both pseudovirus and live virus assays [25]. They observed that 46% and 85% of post-second-dose sera lacked neutralization activity against pseudovirus with triple mutations (K417N, E484K and N501Y) and beta VOC pseudovirus, respectively. For live virus, 61.54% of post-second-dose sera had undetectable neutralization activity against beta VOC, and the remainder with detectable activity were reduced by factors of 4.1–31.5. For mRNA-1273, Shen and coworkers reported that neutralization activity against epsilon variant and beta VOC pseudoviruses were 2–3 and 9–14 times lower than that against D614G pseudovirus, respectively [26].
Summing up the serological findings, delta VOC or other variants reduced neutralizing activity of vaccination-elicited antibodies at varying extents, with beta VOC being apparently more resistant. Beta VOC is characterized by 9 mutations in S1 subunit (D80A, D215G, 241del, 242del and 243del in N-terminal domain (NTD); K417N, E484K and N501Y in receptor-binding domain (RBD); as well as D614G) and a single mutation in S2 subunit (A701V near S1/S2 furin cleavage site) (Table 1). Among these signature amino acid changes, E484K and combination of K417N, E484K and N501Y might reduce the effectiveness of specific monoclonal antibody treatments [27]. K417, E484 and N501 are 3 of the 21 amino acids in RBD, the major target of neutralizing antibodies [28], which interact with human angiotensin-converting enzyme 2 (hACE2) [29]. Results of in vitro studies suggested E484K might play an indispensable role in escape of beta VOC from neutralizing antibodies [30] and loss of neutralizing activity of certain monoclonal antibodies [30, 31, 32]. Despite the difference in resistance against vaccine protection, beta and gamma VOCs share very similar RBD mutation pattern, that is, same signature mutations E484K and N501Y, except for K417 which is less impactful than E484 on antibody neutralization. This suggests possible involvement of other mutations, for instance, in NTD, which may contribute to enhanced resistance of beta VOC against neutralizing antibodies, and this gap of knowledge warrants further investigation.
The unprecedentedly dynamic and vast volume of SARS-CoV-2 genomic data from GISAID database has embodied the power of global and concerted effort on tracking the still-ongoing COVID-19 pandemic. The database is very informative and user-friendly, however, care must be taken when drawing conclusions based on the sequence data, and the following considerations infer the limitations of this study. First, as individual sequences were processed and submitted by various institutes worldwide, there might be variation in quality that we could hardly monitor. For instance, details on coverage and assembly method of some sequences were not provided by the submitters. Second, our data might not reflect the latest situation (especially for June and July 2021) due to highly variable lag time between collection and submission dates. Another factor to consider was inhomogeneous sampling. For instance, at continent level, sequences submitted by the USA comprised about 90% of North American sequences, and the situation of other countries might be under-represented. At country level, we could not ensure submitted sequences were representative of the population without sampling bias towards particular lineages, nor for statistical data such as reported death tolls in certain countries. Another limitation was that we could not distinguish between locally transmitted and imported COVID-19 cases.