Figure 1. Main changes in spike protein reported in articles analyzed.
● Protomer of the spike protein;
● RBD;
● RBM;
● amino acid substitutions,
● ACE2 protein, (PDB structure
[31][32], PyMOL v.4.6).
The changes in SARS-CoV-2 are distributed in various sites in its sequence-like spike protein, N protein, RNA-dependent RNA polymerase (RdRp), NSP3, NSP4 and other open reading frames (ORFs) (Table 1).
Table 1. Changes in SARS-CoV-2 sequences reported in the studies. ¤ Changes in nucleotide sequences.
2.2. SARS-CoV-2 Variants and Clinical Outcomes
Prior to the reporting of Variants of Interest (VOI) and Variants of Concern (VOC), changes in the SARS-CoV-2 sequence that could have an impact on clinical outcomes had been determined [6]. The D614G variant in the spike protein was initially considered to be related to a higher rate of hospitalizations and moderate to severe clinical outcomes [6][17]; however, analyses in different cohorts showed no relationship with disease severity; this change increases the adaptability of the virus in human populations, without necessarily causing more severe disease [19][37]. The same scenario was visualized for the N439K variant in the spike protein, which was also not found to have a direct effect on clinical outcomes, compared to the original virus. However, it was reported that this substitution had emerged in different clades independently and that it increased affinity for ACE2 and resistance against various neutralizing monoclonal antibodies [25].
A study determined that polygenic mutations in SARS-CoV-2 had different outcomes. For mild disease, the following amino acid changes were detected: L84S, G196V in ORF8 and ORF3a, respectively, as well as L37F substitutions in NSP6, F308Y in NSP4 and S197L in the N protein. When analyzing sequences of hospitalized patients, 15 changes distributed in seven genes were found: three in the spike protein, two in RdRp, two in ORF3a, five in N protein, one in ORF6 and two in NSP3; while in fatal outcomes, L71F changes were found in NSP7 and S253P in ORF3a
[18].
In a study where associations between different mutations and clinical outcomes were analyzed, Zekri et al.
[38] found in a sample of 50 patients that the V6 deletion in the spike protein was associated with an increased risk and duration of fever and nasal congestion, while the L3606-Nsp6 deletion was associated with an increased presence of cough and conjunctival congestion.
When variants with changes in P504L, as well as Y541C in NSP13 were analyzed, an association was found between these with infection and mortality rates, without correlation with other studies
[35]. Likewise, the N501Y variant in the spike protein was found to have an increase, without statistical significance, of 18% in terms of risk of fatal outcome
[28].
In silico studies have allowed for a proposal that there are mutation signatures responsible for promoting mild and severe outcomes, in which 20 mutations could be used to separate both groups. These are distributed in the gene encoding the spike protein, as well as in other viral proteins and in untranslated regions (UTRs).
[29] This has allowed for development of models to predict the degree of severity by adjusting the age of patients and analyzing their viral sequences (
https://covidoutcome.com/, accessible from 27 December 2021).
It has been proposed that mutations in ORF1a, ORF1b and in genes encoding N protein were related to a high prevalence of asymptomatic scenarios. However, when D614G, Q57H (ORF3a) and S194L (N protein) changes were present, they were associated with mild and severe outcomes. Likewise, a single nucleotide change (nt14408) in RdRp was associated with severe cases of the disease
[30].
Regarding prolonged viral RNA shedding, which can be up to 100 days in patients with severe disease, one study reported that viral shedding time decreases when A1,430G or C12,473T mutations are present and increases when G227A is present (
p < 0.05). Likewise, mutations in G227A, C7,392T, C15,324T, and C25,626T were mostly represented in severe disease cases
[42].
The analysis of SARS-CoV-2 variants and their impact on clinical outcomes must be seen from an integral perspective; thus, the different levels of structural organization that make up the variants must be evaluated. In this context, it was determined that three structural changes at the RNA and protein levels, specifically A26194T (T268S) and C25611A (synonymous mutation) in the ORF3a region and C28854T (S194L) in the N protein were associated with an increase in severe cases and fatal outcomes (
p < 0.05)
[39].
Methylation at the m6 A loci of the spike protein has been identified in patients debuting with gastrointestinal symptoms, which could provide underlying mechanisms for its change in virulence and transmission capacity during outbreaks and affect the outcome for serious and severe disease
[16].
2.3. Rise and Spread of Variants of Concern
Chronologically, the reported VOCs in the studies analyzed were:
Beta (B.1.351): it was first documented in May 2020, in addition to the D614G substitution, this variant presents other changes such as E484K and N501Y that confer the capacity of immune escape by effect of previous infection or vaccination; the increase in its transmission has been estimated at around 50% compared to the Wuhan variant
[24].
Alpha (B.1.1.7): identified in September 2020, presents a 70% increase in transmissibility, consequence of key changes, specifically in the RBM (N501Y) and near the furin cleavage site (P681H), which could increase the affinity for ACE2 and have an impact on infection and transmission, respectively;
[24] this could have contributed to the rapid dispersion and dominance of this variant in the world before the arrival of the Delta variant (B.1.617.2).
[11][43][44].
Delta: identified in October 2020, it has become the most common variant globally, its main changes are D614G, E484Q and L452R, it has been reported that this variant has biological and clinical implications such as increased risk of hospitalization, longer duration of virus release by infected persons, low Ct values in PCR, greater affinity to the ACE2 receptor, mechanisms of escape to the effect of antibodies and transmissibility increased by 50%
[43][44].
Gamma (P.1): first documented in November 2020, highlighting the presence of three changes that confer affinity for the ACE2 receptor, these are K417T, E484K and N501Y which contribute to its increased transmissibility estimated at 40% in relation to the first variants
[24].
When independently analyzing the clinical outcomes associated with VOC, it was identified that the Alpha and Delta variants affect individuals with similar demographic and comorbidity characteristics, while patients infected with the Gamma variant are older people, mainly between 45 and 64 years old, with a higher probability of presenting cough and anosmia, compared to the other variants
[23][24].
2.4. Other Variants Related with Clinical Outcomes
The dynamics of the SARS-CoV-2 variants analyzed throughout the pandemic has been complex. In France, after the first outbreak there were new variants that had an epidemiological impact; in the comparative study by Fournier et al.
[22] it was determined that the Marseille-4 variant had 13 changes, one of which (S477N) was associated with hypoxemia (
p < 0.05). This variant could be associated with changes in the affinity for ACE2 and decrease the sensitivity of the virus to neutralizing antibodies. In this same context, a cohort study conducted in France determined that lineages B.1.177 and B.1.160, Marseille-2 and Marseille-4, respectively, during the second phase of the pandemic, were associated with more severe clinical outcomes and consequently higher mortality and hospitalization rates
[23], however in this study the association between variants and disease severity was not clear.
Conversely, the B.1.243 lineage was found to be significantly associated with a high degree of disease severity and fatal outcomes. This lineage shows several substitutions in NSP12:P323L, N:S194L as well as D614G and P681H changes in the spike protein
[24].
The B.1.616 lineage whose differences from the original SARS-CoV-2 are centered on nine changes and one deletion in the spike protein (H66D, G142V, Y144del, D215G, V483A, D614G, H655Y, G669S, Q949R, N1187D), as well as changes in other regions, was associated with a high 28-day fatality rate when compared to VOC and other unknown lineages (
p < 0.05)
[26].
Conversely, when analyzing the degree of disease severity with SARS-CoV-2 variants, Al Khatib et al.
[21] identified changes in specific regions of the B.1 and B.2 lineages associated with severe symptoms; patients who developed worse clinical scenarios had greater variability in the SARS-CoV-2 analyzed sequences (
p value 0.001).
When different clades were analyzed with respect to their clinical outcomes, it was determined that the L/V clades (variant of the ORF3a coding protein NS3-G251) were associated with more severe outcomes as they had more pronounced systemic inflammation with higher concentrations of proinflammatory cytokines, chemokines and growth factors compared to the G, S and O clades
[41]. Conversely, when outcomes were analyzed with respect to infection by the G and S/L clades, it was observed that, regardless of clade, the results were similar in terms of rate of hospitalizations and death
[33]. One study reported that clade V was statistically related to increased mortality in uni- and multivariate analyses compared to other variants
[36].
It has been reported that the M1V variant has lower rates of dyspnea, rhinitis and hospitalizations, which has been related to its infection in younger age groups, while the M4V variant infects mainly older adults and has a higher probability of producing fever, lower frequency of cough, rhinitis and olfactory and gustatory disorders, as well as a higher rate of hospitalization associated with hypoxemia. It has also been noted that the M4V variant confers some immunological escape and has been the responsible for cases of reinfection
[22][45].
3. Conclusions
The most identified SARS-CoV-2 variants presented changes in the spike protein, N protein, RdRp, NSP3, as well as in different ORFs sequences. In most of the articles, possible associations between SARS-CoV-2 variants and clinical outcomes were found. However, only eight articles reported significant associations adjusting for age, sex, comorbidities, and other variables. There are multiple factors, such as age and pre-existing diseases, involved in the course of COVID-19 disease, that have been determinant in the degree of severity. Nevertheless, the association between variants and clinical outcomes has not been fully explored at present; more research is required to establish possible associations between SARS-CoV-2 variants and illness behavior.