Pandemic viruses do not always show decreased virulence. Instead, it appears there is an evolutionary equilibrium between transmissibility and virulence. It can be termed this concept “intermittent virulence”. There is a high possibility that new virulent variants will arise in the near future, although it is improbable that SARS-CoV-2’s virulence will be the same as was seen during the alpha or delta waves, due to the fact that the human population has reached a sufficient level of herd immunity through natural infection or due to the vaccination programs. The most recent global mortality data raised the question of whether this pandemic is really over. Furthermore, it is uncertain when the endemic phase will begin. Darwin’s words: “the survival of the fittest” are still valid, and the virus will continue killing nonvaccinated old people, vaccinated old people, and those with comorbidities. We have underestimated the SARS-CoV-2 mastery of immune escape and have not yet seen the full adaptive potential this virus can develop through natural selection.
1. Introduction
The current pandemic has generated a unique opportunity to study SARS-CoV-2 in real time. Several institutions have implemented online monitoring resources to track the progress of the ongoing pandemic. Furthermore, specialized genomic sequencing technologies have contributed to our comprehension of its biology and evolution [
1]. The evolution of viral adaptability is the result of selective pressure acting across different mechanisms [
2], which include adaptation and evolvability. The process of adaptation is how a population organically evolves new or altered features that enable members in the population to more effectively address problems imposed by their prevailing milieu. Evolvability, which is the potential to evolve through natural selection, is a distinct but related concept. A population may be comparatively advantageous with respect to other populations in its ability to adapt to environmental variations [
3].
For example, purifying selection removes virus strains with unfavorable characteristics (e.g., disadvantageous conformational changes in some proteins) within an infected organism, whereas positive selection promotes strains with an evolutionary advantage (e.g., immune evasion) [
3,
4,
5]. Natural selection can have a significant impact on virus-binding proteins, which control cell entry and are the main objective for immune recognition [
6]. Genetic changes in the SARS-CoV-2 spike protein may affect the viral efficiency by promoting immune evasion or increasing ACE2 receptor binding. Importantly, variant lineages that appeared at the end of 2020 have numerous spike protein mutations in common, including K417N (discovered in the beta and gamma variants); E484K; N501Y; and D614G (found in the Alpha, Beta, and Gamma variants) [
5,
7,
8,
9]. In vitro, both D614G and N501Y cause a structural rearrangement of the spike polypeptide chains, which enhances the probability of attaching to ACE2, improves the infectivity of a cell, and uncovers the cleavage region [
8,
9]. K417N reduces monoclonal antibody neutralization while moderately increasing the ACE2-binding affinity. E484K binds to ACE2 with a higher affinity and has a lower monoclonal antibody neutralizing efficacy [
10].
The evolution of a pathogen should be towards a lesser virulence, according to classical evolutionary theory, because killing the host has a detrimental influence on the pathogen’s spread [
11]. Several works published from 1980 to 1990 explained how a pathogen’s evolutionary route to intermediate virulence—rather than zero virulence—might be generated by a positive relationship between disease transmission and virulence [
12,
13,
14,
15,
16]. Some researchers expect that, with the Omicron surge producing a more contagious but less severe disease, the end of COVID-19 is near. Even though the threat of a “super killer” virus is unfounded, the widespread belief that a virus will change to become more lethal during an outbreak exemplifies this phenomenon [
17].
Nevertheless, other experts believe that Omicron’s lessened pathogenicity is coincidental and that continued fast antigenic evolution will certainly yield new variants that will evade immunity and become more virulent. Furthermore, they also believe that one of the most enduring misconceptions about pathogen evolution is that viruses will evolve to become less pathogenic to preserve the life of their hosts [
18]. Viruses mutate to maximize their propagation, which can potentially be associated with increased virulence, such as when large viral loads promote transmission while also increasing morbidity [
18]. If this is the case, pathogens may change to become more virulent. If higher pathogenicity appears later in the disease, upon the canonical transmission window (as it does in influenza virus, hepatitis C virus, HIV, SARS-CoV-2, and many other viruses), it has a minimal impact on viral fitness, and natural selection will not act against it. Predicting virulence evolution is complicated, and Omicron’s decreased virulence is not a strong forecaster regarding the appearance of new variants [
18].
To date, it is very difficult to predict whether Omicron will define the end of the pandemic or new, more virulent variants will emerge and the pandemic will continue. Therefore, in this work, we will compare the main genetic mutations in the Delta and Omicron (BA.1, BA.2. BA.4, and BA.5) variants to unravel their relevance in transmission and virulence, since the latest findings [
19,
20] revealed an enhanced Omicron BA.2, BA.4, and BA.5 pathogenicity in animal models that seems to confirm the predictions for a higher virulence [
18]. However, animal models may not reflect the real situation in humans, where previous exposure to the virus (via either the natural disease or vaccination) confers an effective protection level and could explain why these new Omicron variants are not causing high mortality compared, for example, to the Delta variant.
2. The Delta Variant
The B.1.617.2 (Delta) variant was reported for the first time in India in December 2020 and then in several other nations across the globe. Several variants in the Delta variant spike protein have been recently found (L452R, tbl478K, D614G, tbl19R, 157–158, P681R, and D950N) [
21]. Due to the greater severity of the clinical symptoms elicited by the Delta variant, it quickly outperformed the other variants of concern (VOCs, i.e., variants for which there is evidence of diagnostic detection failures or higher disease severity, as evidenced by increased hospitalizations or deaths or an increase in the transmissibility or a significant reduction in the neutralization by antibodies generated during a previous infection or vaccination or the reduced effectiveness of treatments or vaccines) [
22,
23]. Delta seems to have enhanced its ability to propagate in human cells as a result. Some hypotheses have been presented to clarify its increased infectiousness, including receptor-binding domain (RBD) mutations that enhance its interaction with the receptor [
24], a P681R mutation near the S1-S2 zone that leads to more efficacious furin cleavage [
24,
25], and modifications in its RNA polymerase that increase the viral multiplication. In comparison to the Wuhan strain, Delta was also more virulent and highly fusogenic, according to in vitro tests [
26]. The Delta variant was distinguished by the P681R mutation in the spike protein. Such a change occurred near the SARS-CoV-2 S protein’s furin cleavage site (FCS), and as a consequence, a glycosylation site was lost, leaving the furin cleavage site uncovered, thus promoting syncytia formation and a greater pathogenicity [
27].
Researchers have also discovered two features that were unique to the Delta variant and could explain its infectiousness. Firstly, when the cell membrane has the highest level of Delta spike protein expression, these cells show a greater capacity to fuse with targeted cells that generate lower amounts of ACE2 compared with cells from the other variants [
28]. As the degree of ACE2 expression rises, the distinctions between the variations become less pronounced. Second, pseudoviruses with the Delta spike design penetrate cells expressing the ACE2 receptor faster than other variants. These findings show that the Delta spike protein mutated to improve the merging process for penetrating cells that produce low amounts of ACE2 receptors. This advancement may help to explain why Delta quickly propagated after contact and infect more host cells, shortening the incubation stage and increasing the viral load during the disease [
28].
SARS-CoV-2 infectiousness and virulence can be increased by genetic variations in the spike protein genetic code, which impact its configuration, stability, and functionality [
7,
29]. SARS-CoV-2 infectiousness was discovered to be increased by the S-protein D614G mutation [
30,
31]. Moreover, the structural assessment revealed that changing aspartic acid (D) to glycine (G) at location 614 of the S-protein modified its structure, facilitating with cleaving the furin site [
32,
33,
34]. Patients with the D614G mutation G614 had reduced RT-PCR cycle thresholds, indicating a greater viral load in the upper respiratory system [
7]. A study comparing death rates across nations found a link between the G614 variant and greater mortality [
35]. Compared to the other SARS-CoV-2 strains, the Delta variant has been linked to a 120% enhanced probability of hospitalization, a 287% enhanced probability of being admitted to an intensive care unit, and a 137% enhanced probability of mortality [
36].
3. The Omicron Variant
In November 2021, Omicron was first discovered in South Africa [
37]. Then, the BA.1 lineage of Omicron expanded quickly over the whole world, outcompeting the other variants, such as Delta. Another Omicron variant, the BA.2 lineage, was discovered in numerous nations, including Denmark and the United Kingdom, as of February 2022 [
38]. BA.2 rapidly surpassed BA.1, showing that BA.2 had a higher transmission rate than BA.1. [
19]. Omicron had the same P681H mutation as the Alpha and Delta strains, but it also harbors a new and unique glycosite, threonine (Thr376), that has only been found in Omicron. This was placed close to a proline amino acid that controls O-glycosylation [
39]. When compared with the wild-type or Delta, researchers noticed a notable rise in the use of Core 2-type O-glycans for the Omicron variant, which is compatible with O-glycan insertion [
39]. This discovery is groundbreaking, as it was formerly reported that, in vitro, the impairment of GALNT1 glycosyltransferase function (which regulates the insertion of O-glycans in the vicinity of the furin cleavage area) was caused by a mutation of proline 681 [
27].
This means that Omicron’s new Thr376 mutation recovered the ability to insert O-glycans that cover the furin cleavage site. In conclusion, a mutation in p681 in the Delta variant resulted in glycans loss (
Figure 1A), uncovering the furin cleavage site and making this variant more pathogenic due to enhanced syncytia formation [
27]. On the contrary, a mutation in Thr376 in the Omicron variant induced the addition of O-glycans (
Figure 1B), thus covering the furin cleavage site [
39]. As a result, this mutation has increased immune evasion while decreasing pathogenicity. This phenomenon has also been documented by research on the Nipah virus, where numerous N-glycans on the fusion protein were removed, resulting in hyperfusogenic phenotypes but also showing an enhanced susceptibility to antibody neutralization [
40]. Although glycans are known to perform an important role in immune evasion, they are also emerging as essential determinants of virulence [
41,
42].
In a recent study, it was found that coculturing spike-expressing cells with HEK293-ACE2/TMPRSS2 cells dramatically increased the amount of syncytia formation caused by BA.2 spike compared to BA.1 spike but less than Delta (
Figure 2) [
19]. Since the efficiency of S1/S2 cleavage has been related to cell fusion induced by SARS-CoV-2’s spike protein [
26,
44], it was proposed that BA.2 spike is cleaved in a more effective way than BA.1 spike. Nevertheless, a Western blot test revealed that BA.2 spike is cleaved less efficiently than BA.1 spike, implying that BA.2 spike produces more syncytia than BA.1 spike by using a S1/S2 cleavage-independent mechanism. To find out whether BA.2 S uses TMPRSS2, 293-ACE2 cells with or without TMPRSS2 expression were used in a cell-based fusion assay [
19].
The findings demonstrated that BA.2’s relatively greater fusogenic potential depends on TMPRSS2’s expression in the target cells [
19].
It is worth noting that BA.1 makes inefficient use of TMPRSS2 during infection [
45]. Thus, the new mutations found in BA.2 apparently restored its capacity to use TMPRSS2, which has resulted in a higher fusogenicity and pathogenic potential compared with BA.1. [
19].
Analyses of these genomic differences between the Delta and Omicron variants allow us to understand the relevance of glycans in virulence. Notably, Delta showed a greater fusogenicity due to the p681 mutation that eliminated an O-glycan molecule from the furin cleavage site [
27], and also, the D614G mutation increased syncytia creation and the viral load through enhanced furin-induced spike cleavage [
46]. Interestingly, a similar mutation occurred with the influenza virus. When the virulence of the H3N2 strains in mice was compared, it was discovered that viruses isolated after 1980 had high glycan numbers and caused mild disease in mice. An N-linked glycan from the hemagglutinin receptor of the influenza virus was lost due to a mutation in the gene codifying for such a receptor in the Beijing/89 strain, which was linked to enhanced virulence in mice [
47]. The virulence of the 2009 H1N1 virus was also driven by a sudden glycan loss near the receptor-binding site [
48]. A subsequent study in mice confirmed this discovery, providing evidence that a comparable mutation increased the virulence of the 1918 H1N1 pandemic virus [
49]. The 1918 H1N1 virus contained fewer glycosylation sites on the hemagglutinin receptor than seasonal influenza viruses with lower virulence. Utilizing site-directed mutagenesis, it was discovered that, by incorporating two extra glycosylation sites (asparagine Asn71 and Asn286) in one flank of the hemagglutinin receptor, a highly virulent 1918 HA chimeric virus was considerably attenuated in mice, and removing the glycosylation sites enhanced the virulence in vivo [
49].
Similarly, the Omicron variant BA.1 developed a mutation in Thr376, which is located near the FCS, and consisted of the addition of O-glycan molecules, thus covering the FCS (
Figure 1B). That interference was traduced in a lower fusogenicity and pathogenicity [
40,
45]. This is also consistent with the clinical symptomatology being milder, as a significantly reduced Omicron reproduction velocity was recently detected in lung epithelial cells [
50,
51]. In vitro experiments showed that the Omicron variant BA.1 was less likely to disseminate by cell merging than the other variants, adding to the proof that the virus on its own caused less severe illness [
45,
50].
All these data may clarify why Omicron-BA.1-infected individuals suffered fewer serious complications [
52,
53]. Regarding the pathogenicity of BA.2 in humans, the risk of hospitalization in Germany after a BA.1 or BA.2 Omicron variant infection was approximately 80% lower than after a Delta variant infection, especially in people under 35 years old [
54]. Both BA.1 and BA.2 showed a similar impact on hospitalization or intensive care unit (ICU) admission, implying that, despite evidence of the greater transmissibility of BA.2 [
19], this variant is pathogenically equivalent to BA.1 [
54]. Confirming these results, no clinical differences among individuals infected with BA.1 and BA.2 were detected in Denmark [
55] and South Africa [
56]. In a study carried out in France, the first 207 Omicron BA.2 cases were recorded, and it was discovered that severe Omicron BA.2 infections were exclusively seen in individuals over the age of 80. Three patients (1.5%) died at the ages of 80, 97 and 99, and two of them had diabetes. Two people (ages 80 and 97) received the COVID-19 vaccine (three doses). Individuals who passed away from Omicron BA.2 infection were much older than those who passed away from Omicron BA.1 infection [
57].
The latest research findings showed that a genetic mutation in the E protein from Omicron BA.1 and BA.2 also contributed to its reduced pathogenic potential. The SARS-CoV-2 envelope protein (2-E) creates a homopentameric cation channel that is essential for virulence. SARS-CoV-2 does not travel through the traditional biosynthetic secretory channel; instead, it enters lysosomes and exits via lysosomal deacidification (alkalization). The lysosomes are where coronavirus envelope (E) protein channels are found. The E protein is therefore in charge of the calcium outflow to achieve an alkaline pH, since, otherwise, the virus would be destroyed by the lysosome’s acidic pH. While mutated T9I channels had less of an impact on the luminal pH, the expression of wild-type channels significantly alkalized the lysosomal pH. The viral load was decreased as a result of decreased lysosomal deacidification. T9I’s cytotoxic potential and cytokine production were both roughly 150-fold lower than those of the wild-type, which could further diminish the virulence of Omicron variants [
58].
4. Discussion
Our proposal of an “intermittent virulence” reconciles two contrary hypotheses. The first claims that the Omicron variants (BA.1 and BA.2) were going to define the end of the current pandemic, whereas the other suggests that new virulent variants will arise in the future. We propose that, in order to survive, viruses may develop what we have termed “intermittent virulence”. This concept implies that, when the virus has become more contagious but less virulent, natural selection could create new variants with a higher pathogenic capacity; otherwise, the dominant variant could disappear or become endemic due to the community having reached herd immunity through natural immunity or due to global vaccination programs.
It is possible that the pandemic phase of COVID-19 is about to end but the virus is not going to disappear and an endemic phase will begin, where intermittent virulence could guarantee the survival of this pathogen. Recent findings showed that it is not a rule that viruses become less virulent over time [
64]. A particularly aggressive variant of the human immunodeficiency virus (HIV) has been discovered in the Netherlands, where it has been circulating for some years. More than 100 individuals with HIV-1 subtype B infection experienced a two-fold decline in CD4+ cell numbers than predicted. These people were already at risk of contracting AIDS within two to three years of being diagnosed. This virus lineage, which appears to have emerged
de novo around the 1990s, has undergone considerable changes in its genes, changing almost 300 amino acids, making it difficult to determine the reason for the increased virulence [
64].
Since it is known that genetic recombination can result in the appearance of new, extremely pathogenic variants, the most concerning evolutionary aspect of SARS-CoV-2 is that it could experiment with widespread recombinations [
65,
66]. Simultaneous infections with various subtypes of the same virus might cause this outcome. The activity of the Nsp14 protein regulates this mechanism in SARS-CoV-2 [
67]. In February 2022, for example, evidence of Deltacron XD recombinant SARS-CoV-2 transmission and circulation in Northwest France was reported. Following virological and epidemiological studies, 17 cases of this recombinant SARS-CoV-2 were validated by genotyping or inferred due to epidemiological relationships, indicating an extensive propagation incident and transmission of this virus but not showing evidence of severe clinical symptoms [
68]. Another contemporary research discovered a Delta variant sub-lineage spreading throughout the United States, particularly in Colorado (CO), Texas (TX), and Wyoming (WY). This sub-lineage is characterized by a spike protein mutation at position 112 that has been found in low-prevalence lineages around the world, as well as in circulating U.S. isolates since the end of April 2021. Two unique mutations are found in this Delta S:S112L sub-lineage group: ORF1b:V2354F, which corresponds to nonstructural protein NSP15 at position 303 (NSP15:V303F), and a premature stop codon (Q94 *) that truncates ORF7a [
69]. Unfortunately, the clinical characteristics of the people infected with this sub-lineage were not examined in this investigation.
To predict the epidemic’s future, researchers must accurately investigate how the virus’ tropism evolves. A respiratory tropism is prevalent at this time, but it is known from other coronaviruses that this can evolve very quickly [
70]. The avian coronavirus spike protein, for example, had three amino acid modifications that enabled the virus to attach to kidney cells [
71]. Coronaviruses also are neurotropic in mice [
72] and SARS-CoV-2 also shows neurotropism [
73,
74,
75]. The changing structure of the spike protein is significantly related to natural selection in cell ingress and fusion. SARS-CoV-2 was able to swap hosts and adapt to the human receptor ACE2 [
76] after the acquisition of a furin cleavage site. At least in vitro, the virus has encountered another receptor for entry, specifically the CD147 receptor, a protein present in several tissues, which include epithelial and neural cells [
77]. This is significant, because using a variety of receptors allows for switching between different cell types and, thus, different entry gateways [
70]. Notably, SARS-CoV-2 can infect nonpermissive cells (which do not harbor ACE2 receptors) by using an innovative intracytoplasmic connection mechanism that could serve as an alternative viral transmission pathway, independent of the canonical extra-cytoplasmic ACE2-binding mechanism [
78].
An increased virulence produced by new variants could guarantee the survival of this pathogen, thus confirming the validity of the term “intermittent virulence”. However, it is improbable that the virulence of SARS-CoV-2 would be the same as was seen during the alpha or delta waves, due to the fact that the human population has reached a sufficient level of herd immunity through natural infection or due to the employment of vaccination programs. The most recent global mortality data (
Figure 3) makes us question whether this pandemic is really over, and it also makes uncertain when the endemic phase will begin. Darwin’s words “the survival of the fittest” are still correct. The virus will continue killing nonvaccinated old people, vaccinated old people, and those with comorbidities. For example, it is widely recognized that patients with COVID-19 may have a worse prognosis if their blood glucose levels are high, increasing the risk of multiple organ failure, shock, and the need for intensive care unit (ICU) placement [
79]. People with type 2 diabetes mellitus exhibited a greater incidence, severity of symptoms, and mortality after COVID-19 infection [
80]. According to a recent investigation, glycemic control significantly affects how well the patients would respond immunologically to the SARS-CoV-2 messenger ribonucleic acid (mRNA) vaccine [
81]. The immune response is impaired by hyperglycemia during immunization: a recent study found that 21 days after the first dose of the vaccine, neutralizing antibody titers and CD4 cytokine responses involving type 1 helper T cells were lower in type 2 diabetic patients with glycosylated hemoglobin (HbA1c) levels > 7% than in diabetics with HbA1c levels ≤ 7% [
81]. This aspect is very important, because according to the International Diabetes Atlas, it was estimated that, in 2021, there were 537 million diabetics in the world. Three out of four adults with diabetes live in low- and middle-income countries, and almost one in two (240 million) adults living with diabetes are undiagnosed [
82]. An undiagnosed person will surely have elevated blood glucose levels, which will prevent a satisfactory immune response when vaccinated, as previously reported [
81].
A recent study on the efficacy of the COVID-19 vaccination and the protection decline over time [
83] reported that 8 months after receiving two doses of the COVID-19 vaccine, the immunological function was inferior in those who were older and had preexisting conditions. These findings suggest that booster doses for older people and people with known inadequate or declining vaccine-elicited immunogenicity should be prioritized because they also have the highest likelihood of experiencing severe COVID-19 symptoms if infected [
83].
While this work was in peer review, a new preprint was posted on Biorxiv [
84], which showed that Omicron BA.4 and BA.5 replication is linked to the lower activation of the epithelial innate immune responses. These subvariants have improved transmission and potentially reduce immune protection from severe disease by combining the evolution of antibody escape with the increased antagonism of interferon signaling. The study also discovered the increased expression of the innate immune antagonist proteins Orf6 and N, which are comparable to Alpha, implying that human adaptation processes are similar [
84].
As a final comment, one should remember that the in vitro antibody neutralization tests performed to evaluate the efficacy of vaccines or monoclonal antibodies do not mirror the real situation within the host, where SARS-CoV-2 can hijack infected cells and release decoy targets to avoid being neutralized by vaccine-derived antibodies, or can infect non-permissive cells by using tunneling nanotubes or through the cell-to-cell infection, thus avoiding immune surveillance and causing syncytia-induced lymphopenia [
85]. We have underestimated this master of immune escape, and have not yet seen the full adaptive potential this virus can develop through natural selection.
This entry is adapted from the peer-reviewed paper 10.3390/covid2080080