Clinical and Virological Outcome of Monoclonal Antibody Therapies Across SARS-CoV-2 Variants in 245 Immunocompromised Patients: A Multicenter Prospective Cohort Study.

BACKGROUND: Immunocompromised patients (ICPs) have an increased risk for a severe and prolonged COVID-19. SARS-CoV-2 monoclonal antibodies (mAbs) were extensively used in these patients, but data from randomized trials that focus on ICPs are lacking. We evaluated the clinical and virological outcome of COVID-19 in ICPs treated with mAbs across SARS-CoV-2 variants. METHODS: In this multicenter prospective cohort study, we enrolled B-cell- and/or T-cell-deficient patients treated with casirivimab/imdevimab, sotrovimab, or tixagevimab/cilgavimab. SARS-CoV-2 RNA was quantified and sequenced weekly, and time to viral clearance, viral genome mutations, hospitalization, and death rates were registered. RESULTS: Two hundred and forty five patients infected with the Delta (50%) or Omicron BA.1, 2, or 5 (50%) variant were enrolled. Sixty-seven percent were vaccinated; 78 treated as outpatients, of whom 2 required hospital admission, but both survived. Of the 159 patients hospitalized at time of treatment, 43 (27%) required mechanical ventilation or died. The median time to viral clearance was 14 days (interquartile range, 7-22); however, it took >30 days in 15%. Resistance-associated spike mutations emerged in 9 patients in whom the median time to viral clearance was 63 days (95% confidence interval, 57-69; P < .001). Spike mutations were observed in 1 of 42 (2.4%) patients after treatment with 2 active mAbs, in 5 of 34 (14.7%) treated with actual monotherapy (sotrovimab), and 3 of 20 (12%) treated with functional monotherapy (ie, tixagevimab/cilgavimab against tixagevimab-resistant variant). CONCLUSIONS: Despite treatment with mAbs, morbidity and mortality of COVID-19 in ICPs remained substantial. Combination antiviral therapy should be further explored and may be preferred in severely ICPs.

SARS-CoV-2 journey: from alpha variant to omicron and its sub-variants.

The COVID-19 pandemic has affected hundreds of millions of individuals and caused more than six million deaths. The prolonged pandemic duration and the continual inter-individual transmissibility have contributed to the emergence of a wide variety of SARS-CoV-2 variants. Genomic surveillance and phylogenetic studies have shown that substantial mutations in crucial supersites of spike glycoprotein modulate the binding affinity of the evolved SARS-COV-2 lineages to ACE2 receptors and modify the binding of spike protein with neutralizing antibodies. The immunological spike mutations have been associated with differential transmissibility, infectivity, and therapeutic efficacy of the vaccines and the immunological therapies among the new variants. This review highlights the diverse genetic mutations assimilated in various SARS-CoV-2 variants. The implications of the acquired mutations related to viral transmission, infectivity, and COVID-19 severity are discussed. This review also addresses the effectiveness of human neutralizing antibodies induced by SARS-CoV-2 infection or immunization and the therapeutic antibodies against the ascended variants.

SARS-CoV-2 Recombinants: Genomic Comparison between XBF and Its Parental Lineages.

Recombination events are very common and represent one of the primary drivers of RNA virus evolution. The XBF SARS-CoV-2 lineage is one of the most recently generated recombinants during the COVID-19 pandemic. It is a recombinant of BA.5.2.3 and BA.2.75.3, both descendants of lineages that caused many concerns (BA.5 and BA.2.75, respectively). Here, we performed a genomic survey focused on comparing the recombinant XBF with its parental lineages to provide a comprehensive assessment of the evolutionary potential, epidemiological trajectory, and potential risks. Genetic analyses indicated that although XBF initially showed the typical expansion depicted by a steep curve, causing several concerns, currently there is no indication of significant expansion potential or a contagion rate surpassing that of other currently active or previously prevalent lineages. BSP indicated that the peak has been reached around 19 October 2022 and then the genetic variability suffered slight oscillations until early 5 March 2023 when the population size reduced for the last time starting its last plateau that is still lasting. Structural analyses confirmed its reduced potential, also indicating that properties of NTDs and RBDs of XBF and its parental lineages present no significant difference. Of course, cautionary measures must still be taken and genome-based monitoring remains the best tool for detecting any important changes in viral genome composition.

Molecular determinants associated with temporal succession of SARS-CoV-2 variants in Uttar Pradesh, India.

The emergence and rapid evolution of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) caused a global crisis that required a detailed characterization of the dynamics of mutational pattern of the viral genome for comprehending its epidemiology, pathogenesis and containment. We investigated the molecular evolution of the SASR-CoV-2 genome during the first, second and third waves of COVID-19 in Uttar Pradesh, India. Nanopore sequencing of the SARS-CoV-2 genome was undertaken in 544 confirmed cases of COVID-19, which included vaccinated and unvaccinated individuals. In the first wave (unvaccinated population), the 20A clade (56.32%) was superior that was replaced by 21A Delta in the second wave, which was more often seen in vaccinated individuals in comparison to unvaccinated (75.84% versus 16.17%, respectively). Subsequently, 21A delta got outcompeted by Omicron (71.8%), especially the 21L variant, in the third wave. We noticed that Q677H appeared in 20A Alpha and stayed up to Delta, D614G appeared in 20A Alpha and stayed in Delta and Omicron variants (got fixed), and several other mutations appeared in Delta and stayed in Omicron. A cross-sectional analysis of the vaccinated and unvaccinated individuals during the second wave revealed signature combinations of E156G, F157Del, L452R, T478K, D614G mutations in the Spike protein that might have facilitated vaccination breach in India. Interestingly, some of these mutation combinations were carried forward from Delta to Omicron. In silico protein docking showed that Omicron had a higher binding affinity with the host ACE2 receptor, resulting in enhanced infectivity of Omicron over the Delta variant. This work has identified the combinations of key mutations causing vaccination breach in India and provided insights into the change of [virus's] binding affinity with evolution, resulting in more virulence in Delta and more infectivity in Omicron variants of SARS-CoV-2. Our findings will help in understanding the COVID-19 disease biology and guide further surveillance of the SARS-CoV-2 genome to facilitate the development of vaccines with better efficacies.

Using In Silico Bioinformatics Algorithms for the Accurate Prediction of the Impact of Spike Protein Mutations on the Pathogenicity, Stability, and Functionality of the SARS-CoV-2 Virus and Analysis of Potential Therapeutic Targets.

Coronavirus disease 2019 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). We have used bioinformatics to investigate seventeen mutations in the spike protein of SARS-CoV-2, as this mediates infection of human cells and is the target of most vaccine strategies and antibody-based therapies. Two mutations, H146Y and S221W, were identified as being most pathogenic. Mutations at positions D614G, A829T, and P1263L might also have deleterious effects on protein function. We hypothesized that candidate small molecules may be repurposed to combat viral infection. We investigated changes in binding energies of the ligands and the mutant proteins by assessing molecular docking. For an understanding of cellular function and organization, protein-protein interactions are also critical. Protein-protein docking for naïve and mutated structures of SARS-CoV-2 S protein was evaluated for their binding energy with the angiotensin-converting enzyme 2 (ACE2). These interactions might limit the binding of the SARS-CoV-2 spike protein to the ACE2 receptor or may have a deleterious effect on protein function that may limit infection. These results may have important implications for the transmission of SARS-CoV-2, its pathogenesis, and the potential for drug repurposing and immune therapies.

Antibody response and intra-host viral evolution after plasma therapy in COVID-19 patients pre-exposed or not to B-cell-depleting agents.

Administration of plasma therapy may contribute to viral control and survival of COVID-19 patients receiving B-cell-depleting agents that impair humoral immunity. However, little is known on the impact of anti-CD20 pre-exposition on the kinetics of SARS-CoV-2-specific antibodies. Here, we evaluated the relationship between anti-spike immunoglobulin G (IgG) kinetics and the clinical status or intra-host viral evolution after plasma therapy in 36 eligible hospitalized COVID-19 patients, pre-exposed or not to B-cell-depleting treatments. The majority of anti-CD20 pre-exposed patients (14/17) showed progressive declines of anti-spike IgG titres following plasma therapy, contrasting with the 4/19 patients who had not received B-cell-depleting agents (p = 0.0006). Patients with antibody decay also depicted prolonged clinical symptoms according to the World Health Organization (WHO) severity classification (p = 0.0267) and SARS-CoV-2 viral loads (p = 0.0032) before complete virus clearance. Moreover, they had higher mutation rates than patients able to mount an endogenous humoral response (p = 0.015), including three patients with one to four spike mutations, potentially associated with immune escape. No relevant differences were observed between patients treated with plasma from convalescent and/or mRNA-vaccinated donors. Our study emphasizes the need for an individualized clinical care and follow-up in the management of COVID-19 patients with B-cell lymphopenia.

Genomic Surveillance of SARS-CoV-2 in the Southern Province of Zambia: Detection and Characterization of Alpha, Beta, Delta, and Omicron Variants of Concern.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) have significantly impacted the global epidemiology of the pandemic. From December 2020 to April 2022, we conducted genomic surveillance of SARS-CoV-2 in the Southern Province of Zambia, a region that shares international borders with Botswana, Namibia, and Zimbabwe and is a major tourist destination. Genetic analysis of 40 SARS-CoV-2 whole genomes revealed the circulation of Alpha (B.1.1.7), Beta (B.1.351), Delta (AY.116), and multiple Omicron subvariants with the BA.1 subvariant being predominant. Whereas Beta, Delta, and Omicron variants were associated with the second, third, and fourth pandemic waves, respectively, the Alpha variant was not associated with any wave in the country. Phylogenetic analysis showed evidence of local transmission and possible multiple introductions of SARS-CoV-2 VOCs in Zambia from different European and African countries. Across the 40 genomes analysed, a total of 292 mutations were observed, including 182 missense mutations, 66 synonymous mutations, 23 deletions, 9 insertions, 1 stop codon, and 11 mutations in the non-coding region. This study stresses the need for the continued monitoring of SARS-CoV-2 circulation in Zambia, particularly in strategically positioned regions such as the Southern Province which could be at increased risk of introduction of novel VOCs.

Unique mutations in SARS-CoV-2 Omicron subvariants' non-spike proteins: Potential impacts on viral pathogenesis and host immune evasion.

SARS-CoV-2 is the causative agent behind the ongoing COVID-19 pandemic. This virus is a cumulative outcome of mutations, leading to frequent emergence of new variants and their subvariants. Some of them are a matter of high concern, while others are variants of interest for studying the mutational effect. The major five variants of concern (VOCs) are Alpha (B.1.1.7), Beta (B.1.315), Gamma (P.1), Delta (B.1.617.2), and Omicron (B.1.1.529.*/BA.*). Omicron itself has >100 subvariants at present, among which BA.1 (21K), BA.2 (21L), BA.4 (22A), BA.5 (22B), and BA.2.12.1 (22C) are the dominant ones. Undoubtedly, these variants and sometimes their progeny subvariants have significant differences in their spike region that impart them the unique properties they harbor. But alongside, the mutations in their non-spike regions could also be responsible elements behind their characteristics, such as replication time, virulence, survival, host immune evasion, and such. There exists a probability that these mutations of non-spike proteins may also impart epistatic effects that are yet to be brought to light. The focus of this review encompasses the non-spike mutations of Omicron, especially in its widely circulating subvariants (BA.1, BA.2, BA.4, BA.5, and BA.2.12.1). The mutations such as in NSP3, NSP6, NSP13, M protein, ORF7b, and ORF9b are mentioned few of all, which might have led to the varying properties, including growth advantages, higher transmission rate, lower infectivity, and most importantly better host immune evasion through natural killer cell inactivation, autophagosome-lysosome fusion prevention, host protein synthesis disruption, and so on. This aspect of Omicron subvariants has not yet been explored. Further study of alteration of expression or interaction profile of these non-spike mutations bearing proteins, if present, can add a great deal of knowledge to the current understanding of the viral properties and thus effective prevention strategies.

A herd immunity approach to the COVID-19 pandemic?

Viral contamination is one of the most urgent and important topics of environmental pollution. COVID-19 is primarily transmitted from person to person, but can also be transmitted from person to animal. Herd immunity must meet the requirements in order to fulfill the goal of mitigating and ending COVID-19. This paper shows five reasons or conditions why herd immunity is not achieved in the present policies without proposed effective strategies in this paper. Unless one of the five reasons for the herd immunity model is met, the promise of herd immunity will not be fulfilled. Many COVID-19 policies worldwide with current vaccines do not meet the requirements. Policymakers have been relying on unreliable R. The number of daily deaths instead of the number of cases is a good indicator of the pandemic which will be mainly used in this paper. Currently, even in vaccinated countries, resurgences are being observed with new variants with spike mutations and immune escape. This paper proposes an effective multipronged approach such as a pharmacological approach and a non-pharmacological approach including digital fencing. Two tools such as scorecovid and deathdaily were used for justifying the claims. Digital fencing as well as pharmacological approaches may be able to overcome the pandemic. Two tools such as scorecovid and deathdaily showed that the proposed multipronged approach will be effective for mitigating the pandemic.

From Alpha to Delta-Genetic Epidemiology of SARS-CoV-2 (hCoV-19) in Southern Poland.

In Poland, the first case of SARS-CoV-2 infection was confirmed in March 2020. Since then, many circulating virus lineages fueled rapid pandemic waves which inflicted a severe burden on the Polish healthcare system. Some of these lineages were associated with increased transmissibility and immune escape. Mutations in the viral spike protein, which is responsible for host cell recognition and serves as the primary target for neutralizing antibodies, are of particular importance. We investigated the molecular epidemiology of the SARS-CoV-2 clades circulating in Southern Poland from February 2021 to August 2021. The 921 whole-genome sequences were used for variant identification, spike mutation, and phylogenetic analyses. The Pango B.1.1.7 was the dominant variant (n = 730, 89.68%) from March 2021 to July 2021. In July 2021, the B.1.1.7 was displaced by the B.1.617.2 lineage with 66.66% in July 2021 and 92.3% in August 2021 frequencies, respectively. Moreover, our results were compared with the sequencing available on the GISAID platform for other regions of Poland, the Czech Republic, and Slovakia. The analysis showed that the dominant variant in the analyzed period was B.1.1.7 in all countries and Southern Poland (Silesia). Interestingly, B.1.1.7 was replaced by B.1.617.2 earlier in Southern Poland than in the rest of the country. Moreover, in the Czech Republic and Slovakia, AY lineages were predominant at that time, contrary to the Silesia region.

Genomic surveillance of the Lambda SARS-CoV-2 variant in a global phylogenetic context.

The massive sequencing of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and global genomic surveillance strategies allowed the detection of many variants of concern and interest. The variant of interest Lambda (C.37), which originated in South America, has been the most prevalent in Peru and Chile, but its dispersion in other continents still remains unknown. The current study aims to determine the phylogenetic relationship among C.37 isolates worldwide, focusing on spike mutations to understand the spread of Lambda in pandemics. A total of 7441 sequences identified as C.37 were downloaded from the GISAID database; local analysis was carried out to identify spike mutations and phylogenetic analysis was carried out to determine the rate of spread of the virus. Our results showed some spike mutations of Lambda that allowed us to detect small local outbreaks in different countries that occurred in the past and identify several clades that have not yet been designated. Although the lineage C.37 is not epidemiologically relevant in Europe or North America, the endemic behavior of this variant in Peru had a major impact on the second SARS-CoV-2 wave.

Mapping of SARS-CoV-2 spike protein evolution during first and second waves of COVID-19 infections in India.

Aim: The aim of this study was to investigate the SARS-CoV-2 spike protein evolution during the first and second wave of COVID-19 infections in India. Materials & Methods: Detailed mutation analysis was done in 763 samples taken from GISAID for the ten most affected Indian states between March 2020 to August 2021. Results: The study revealed 242 mutations corresponding to 207 sites. Fifty one novel mutations emerged during the assessment period, including many with higher transmissibility and immune evasion functions. Highest number of mutations per spike protein also rose from 5 (first wave) to 13 (second wave). Conclusion: The study identified mutation-rich and no mutation regions in the spike protein. The conserved spike regions can be useful for designing future diagnostics, vaccines and therapeutics.

Effects of Spike Mutations in SARS-CoV-2 Variants of Concern on Human or Animal ACE2-Mediated Virus Entry and Neutralization.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a zoonotic agent capable of infecting humans and a wide range of animal species. Over the duration of the pandemic, mutations in the SARS-CoV-2 spike (S) protein have arisen, culminating in the spread of several variants of concern (VOCs) with various degrees of altered virulence, transmissibility, and neutralizing antibody escape. In this study, we used pseudoviruses that express specific SARS-CoV-2 S protein substitutions and cell lines that express angiotensin-converting enzyme 2 (ACE2) from nine different animal species to gain insights into the effects of VOC mutations on viral entry and antibody neutralization capability. All animal ACE2 receptors tested, except mink, support viral cell entry for pseudoviruses expressing the ancestral prototype S at levels comparable to human ACE2. Most single S substitutions did not significantly change virus entry, although 614G and 484K resulted in a decreased efficiency. Conversely, combinatorial VOC substitutions in the S protein were associated with increased entry of pseudoviruses. Neutralizing titers in sera from various animal species were significantly reduced against pseudoviruses expressing the S proteins of Beta, Delta, or Omicron VOCs compared to the parental S protein. Especially, substitutions in the S protein of the Omicron variant significantly reduced the neutralizing titers of the sera. This study reveals important insights into the host range of SARS-CoV-2 and the effect of recently emergent S protein substitutions on viral entry, virus replication, and antibody-mediated viral neutralization. IMPORTANCE The ongoing coronavirus disease 2019 (COVID-19) pandemic, caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), continues to have devastating impacts on global health and socioeconomics. The recent emergence of SARS-CoV-2 variants of concern, which contain mutations that can affect the virulence, transmission, and effectiveness of licensed vaccines and therapeutic antibodies, are currently becoming the common strains circulating in humans worldwide. In addition, SARS-CoV-2 has been shown to infect a wide variety of animal species, which could result in additional mutations of the SARS-CoV-2 virus. In this study, we investigate the effect of mutations present in SARS-CoV-2 variants of concern and determine the effects of these mutations on cell entry, virulence, and antibody neutralization activity in humans and a variety of animals that might be susceptible to SARS-CoV-2 infection. This information is essential to understand the effects of important SARS-CoV-2 mutations and to inform public policy to create better strategies to control the COVID-19 pandemic.

Delta variant (B.1.617.2) of SARS-CoV-2: Mutations, impact, challenges and possible solutions.

Since commencement of COVID-19 pandemic, several SARS-CoV-2 variants have emerged amid containment efforts via vaccination. The Delta variant (B.1.617.2), discovered in October 2020, was designated as a VOC by the WHO on May 11, 2021. The enhanced transmissibility of Delta variant has been associated with critical mutations such as D614G, L452R, P681R, and T478K in the S-protein. The increased affinity of the S-protein and ACE2 has been postulated as a key reason for decreased vaccine efficacy. As per evidence, the Delta variant possesses increased transmissibility and decreased vaccine efficacy compared to other VOCs like Alpha and Beta. This has led to concerns regarding the acquisition of novel mutations in the Delta variant and outbreaks in vulnerable communities, including vaccinated people. In this mini-review of Delta variant, we have explained its evolution and characteristics, the impact of spike mutations on infectivity and immune evasion, and measures to combat future outbreaks.

SARS-CoV-2 variants - Evolution, spike protein, and vaccines.

Despite the rising natural and vaccines mediated immunity, several countries have experienced a resurgence of the Coronavirus disease of 2019 (COVID-19) due to the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants. From Alpha to Omicron, the variants of concern (VOC) have evolved several spike protein mutations that may have an impact on virus characteristics, such as transmissibility and antigenicity. In this review, we describe the evolution of SARS-CoV-2, summarize current knowledge of epidemiological and clinical features of the variants, and discuss the response strategies in terms of vaccines to reduce the burden of COVID-19.

Omicron - The new SARS-CoV-2 challenge?

SARS-CoV-2 virus has infected nearly 300 M people worldwide and has been associated with over 6 M deaths by March 2022. Since the virus emergence in December 2019 in Wuhan, several new mutations have been described. The World Health Organization has developed a working name for these emerging variants according to their impact on the worldwide population. In this context a high alert has been paid to variants of concern (VOC) due to their high infectiousness and transmissibility patterns. The most recent VOC, Omicron (B.1.1.529), has become dominant in the shortest time ever and has placed Europe under an overwhelming and unprecedented number of new cases. This variant has numerous mutations in regions that are associated with higher transmissibility, stronger viral binding, affinity and antibody escape. Moreover, the mutations and deletions present in the spike protein suggest that the SARS-CoV-2 specific attachment inhibitors may not be the best option for Omicron therapy. Omicron is the dominant variant circulating worldwide and, at the end of February 2022, it was responsible for nearly all sequences reported to GISAID. Omicron is made up of several sublineages, where the most common are BA.1 and BA.2 (or Nextstrain clade 21K and 21L, respectively). At a global level, it is possible to say that the proportion of BA.2 has been increasing relative to BA.1 and in some countries it has been replacing it at high rates. In order to better assess the Omicron effectiveness on antibody escape, spread and infectious ability it is of the highest relevance to maintain a worldwide tight surveillance. Even though this variant has been associated with a lower death rate, it is important to highlight that the number of people becoming infected is concerning and that further unpredictable mutations may emerge as the number of infected people rises.

SARS-CoV-2 variants and vaccination.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused global destruction since its emergence in late 2019. Over the past two years, the virus has continuously evolved in human host, leading to emergence of variants with changed viral transmission, disease severity, and evasion of immunity. Although vaccines have been developed for the coronavirus disease 2019 (COVID-19) at an unprecedently pace, the variants have constantly posed threats to the effectiveness of the approved vaccines. In this short communication, we review the key variants and discuss their implications in viral replication, transmission, and immune evasion.

Versatile role of ACE2-based biosensors for detection of SARS-CoV-2 variants and neutralizing antibodies.

Since the beginning of the COVID-19 pandemic, accumulating mutations have led to marked changes in the genetic sequence of SARS-CoV-2. Of these, mutations in the spike (S) protein can alter the properties of the virus, particularly transmissibility and antigenicity. However, it is difficult to detect antigenic variants of the SARS-CoV-2 S protein by immunoassay. Here, we developed an ACE2-based biosensor designed to detect both SARS-CoV-2 S1 mutations and neutralizing antibodies. In "binding mode", the biosensor works by detecting binding of the S protein to an immobilized ACE2 receptor. The ACE2-based biosensor was able to detect S1 proteins of the alpha (500 pg/mL) and beta variants (10 ng/mL), as well as wild-type S1 (10 ng/mL), of SARS-CoV-2. The biosensor distinguished wild-type SARS-CoV-2 S1 from the S1 alpha and beta variants via color differences. In addition, a slight modification to the protocol enabled the ACE2-based biosensor to operate in "blocking mode" to detect neutralizing antibodies in serum samples from COVID-19 patients. Therefore, the ACE2-based biosensor is a versatile test for detecting wild-type S1, S1 mutants, and neutralizing antibodies against SARS-CoV-2. This approach to targeting both the mechanism by which SARS-CoV-2 enters host cells and the subsequent adaptive immune response will facilitate the development of various biosensors against SARS-CoV-2.

Evidence of a SARS-CoV-2 double Spike mutation D614G/S939F potentially affecting immune response of infected subjects.

OBJECTIVES: Despite extensive efforts to monitor the diffusion of COVID-19, the actual wave of infection is worldwide characterized by the presence of emerging SARS-CoV-2 variants. The present study aims to describe the presence of yet undiscovered SARS-CoV-2 variants in Italy. METHODS: Next Generation Sequencing was performed on 16 respiratory samples from occasionally employed within the Bangladeshi community present in Ostia and Fiumicino towns. Computational strategy was used to identify all potential epitopes for reference and mutated Spike proteins. A simulation of proteasome activity and the identification of possible cleavage sites along the protein guided to a combined score involving binding affinity, peptide stability and T-cell propensity. RESULTS: Retrospective sequencing analysis revealed a double Spike D614G/S939F mutation in COVID-19 positive subjects present in Ostia while D614G mutation was evidenced in those based in Fiumicino. Unlike D614G, S939F mutation affects immune response by the slight but significant modulation of T-cell propensity and the selective enrichment of potential binding epitopes for some HLA alleles. CONCLUSION: Collectively, our findings mirror further the importance of deep sequencing of SARS-CoV-2 genome as a unique approach to monitor the appearance of specific mutations as for those herein reported for Spike protein. This might have implications on both the type of immune response triggered by the viral infection and the severity of the related illness.