Sputnik V Effectiveness against Hospitalization with COVID-19 during Omicron Dominance

Abstract

Mass vaccination campaigns against COVID-19 affected more than 90% of the population in most developed countries. The new epidemiologic wave of COVID-19 has been ongoing since the end of 2021. It is caused by a virus variant B.1.1.529, also known as “Omicron” and its descendants. The effectiveness of major vaccines against Omicron is not known. The purpose of this study is to evaluate the efficacy of the Sputnik V vaccine. The main goal is to assess its protection against hospitalization in the period of Omicron dominance. We conducted our study based on a large clinical center in Moscow (Russia) where 1112 patients were included. We used the case-population method to perform the calculations. The data we obtained indicate that the Omicron variant causes at least 90% of infections in the studied cohort. The effectiveness of protection against hospitalization with COVID-19 in our study was 85.9% (95% CI 83.0–88.0%) for those who received more than one dose. It was 87.6% (95% CI 85.4–89.5%) and 97.0% (95% CI 95.9–97.8%) for those who received more than two or three doses. The effectiveness in cases of more severe forms was higher than for less severe ones. Thus, present study indicates the high protective efficacy of vaccination against hospitalization with COVID-19 in case of Omicron lineage.

1. Introduction

The emergence of several vaccines with high efficacy for COVID-19 prevention was the greatest success of the first pandemic year, including Sputnik V (Gamaleya Research Centre) [1], mRNA-1273 (Moderna) [2], BNT162b2 (Pfizer) [3], and ChAdOx1 nCoV-19 Vaccine AZD1222 (Oxford, UK, AstraZeneca) [4]. The efficacy of adenoviral vector and mRNA platforms was rather high. It exceeded 90% concerning protection against the disease and 95% concerning protection against hospitalization, according to randomized double-blind clinical studies. This prompted the launch of extensive vaccination programs around the world [5]. According to the field studies, the achievement of high vaccination levels in several developed [6] and developing [7,8] countries has demonstrated the high effectiveness of such vaccines. What is more, it has significantly reduced the incidence and need for the hospitalization of patients with COVID-19 for some time [9].

Unfortunately, new variants of concern (VOC) began to appear (including Alpha B.1.1.7, Beta B.1.351, Gamma P.1, Delta B.1.1617.2, and Omicron B.1.1.529) [10]. All VOC had mutations in the receptor-binding domain (RBD) of the S protein, as well as many mutations in other parts of the genome. This led to increased transmissivity and the ability to partially escape from the virus-neutralizing activity of antibodies gained following vaccination and previous disease.

Omicron appeared to be the champion in the number of new mutations [11]. In total, this variant has more than 50 mutations, 17 of which are in the RBD region. The wide distribution of the Omicron variant has led to the largest wave of COVID-19 morbidity in the world. The increase in cases has been recorded since November 2021 [12]. In April 2022, this wave is decreasing. To note, the largest wave of Omicron incidence occurred amid a high percentage of vaccinations in developed countries. Booster immunization (re-vaccination) had already been widely used in many countries during the emergence of Omicron. Re-vaccination occurred 5–6 months following the initial course of two injections. The recommendation for booster immunization was developed during the period of the Delta variant distribution [13,14].

The wide distribution of the Omicron variant amid a significant number of vaccinations and re-vaccinations was observed. It raises the question of the need to study the efficacy of COVID-19 vaccines against this variant. This requires a better understanding of the extent to which the epidemiological efficacy decreases concerning primary infection and hospitalization for vaccinated people during the spread of the Omicron variant.

The Sputnik V vaccine action (Gam-COVID-Vac) implies the use of adenoviral vectors that are not capable of replication in the human body. Yet, these vectors can deliver the SARS-CoV-2 S glycoprotein gene into human cells. Heterologous prime-boost immunization used for the Sputnik V vaccine requires two adenoviruses of different serotypes: Ad26 and Ad5. This makes it possible to recommend a three-week interval between the first and second vaccine doses [15]. The efficacy of the Sputnik V vaccine based on the results of a randomized placebo-controlled trial amounted to 91.6% (95% CI 85.6–95.2%) [1]. The widespread use of the vaccine in Russia, Hungary, Argentina, and Bahrain made it possible to estimate the epidemiological effectiveness of the vaccine in relation to the initial variants of the virus, as well as the Delta variant. Thus, the effectiveness of the first dose of Sputnik V in Argentina was 78.6% (95% CI 74.8–81.7%), three weeks after vaccination. It was estimated during the first five weeks, during the spread of the initial variants of the virus [16]. As for laboratory-confirmed infection and hospitalization in the 60–79-year-old group, the efficacy was 87.6% (95% CI 80.3–92.2%). The effectiveness of first-dose vaccination was estimated in the age group over 60 years during Delta variant dominance. The effectiveness reached 51.98% (95% CI 35.61–64.19%) during the first three months following immunization. At the same time, in the group up to 60 years old, effectiveness was 75.28% (95% CI 69.24–80.13%). Thus, two shots among patients over 60 years old during the period of Delta variant domination is recommended [17].

The national observational study in Hungary HUN-VE (Hungarian Vaccine Efficacy) showed the effectiveness of double immunization with Sputnik V at the level of 85.7% during the period of Alpha variant dominance. It presented 85.7% (95% CI 84.3–86.9%) and 97.5% (95% CI 95.6–98.6%) in protection against infection and death, respectively. These numbers are based on observations during a few months following immunization [18]. The epidemiological efficacy of two doses of Sputnik V was shown in Bahrain for laboratory-confirmed infection and hospitalization [19]. Beta B.1.351, Delta B.1.617.2, and Kappa B.1.617.1 were common during the observation period. According to the data obtained in Russia, in relation to a retrospective cohort study conducted during Delta variant, the effectiveness of Sputnik V was rated at 80% [20].

In a cohort of people living with HIV (PLWH), Sputnik V effectiveness varied. It depended not only on the dominant SARS-CoV-2 variant, but also on the level of CD4+ cells [21]. The effectiveness of Sputnik V among PLWH with CD4 ≥ 350 cells/mL was 81.17% (95% CI 49.13–93.03%) against the original strain. Concerning the Delta variant, it was 65.35% (95% CI 52.61–74.66%). However, the epidemiological efficiency of protection against a severe course of the disease in PLWH with CD4 levels ≥ 350 cells/mL remained at the level of 93.05%.

The effectiveness of Sputnik V against the Omicron variant requires an additional study. The aim of our work is to evaluate the dependence of Sputnik V effectiveness from the number of doses received. We also evaluate the degree of protection in relation to conditions of varying severity.

2. Materials and Methods

2.1. Research Design and Ethical Issues

The purpose of the present study was to evaluate the effectiveness of the vaccines Sputnik Light, Sputnik V, and additional booster doses during the period of Omicron variant dominance. For this purpose, the random sampling of 1112 patients was conducted at the large City Clinical Hospital No. 67 named after L.A. Vorokhobov of the Moscow City Department (Moscow, Russia). It took place during the period of increasing Omicron incidence, from 11 January to 21 February 2022. Moreover, subsequent monitoring was performed. The only inclusion criterion was a positive PCR test for COVID-19. There was no exclusion criterion. At the time of admission, all hospitalized patients included in the study had nasopharyngeal secretion samples analyzed for the presence of SARS-CoV-2 virus. Their blood serum was also analyzed for the presence of IgG antibodies in NC and RBD antigens. The detection of SARS-CoV-2 RNA was performed by real-time RT-PCR and whole-genome sequencing to differentiate the virus’s genetic line. The severity of the disease was determined by WHO recommendations.

Written informed consent was obtained from all subjects in accordance with the order of the Ministry of Health of the Russian Federation of 21 July 2015 N 474n. All samples were de-identified prior to receipt by the research team. The study was submitted to the Local Ethics Committee of the Gamaleya Center. The committee concluded that the study did not use identifiable biological samples and did not provide any confidential information (decision #14 29 September 2021.).

2.2. Sample Collection and RT-PCR

Nasopharyngeal swabs were collected in transport medium for viruses (HEM, catalog number G00155, Russia). Total RNA was extracted from patients’ swabs using a kit for isolation RNA from animal cells/bacteria, smear/scraping of epithelial cells, and viruses with spin-column catalog number RU-250 (Biolabmix, Novosibirsk, Russia). Extracted RNA was tested for SARS-CoV-2 by using a one-step «SARS-CoV-2 FRT» commercial kit with catalog number EA-128 (from N.F. Gamaleya NRCEM, Moscow, Russia). The concentrations of viral samples were determined using in-house plasmid controls with known copy numbers. Specimens with Ct values less than 30 were selected for whole-genome sequencing. In addition, samples were tested for the identification of SARS-CoV-2 Omicron and Delta genetic lines with «AmpliTest SARS-CoV-2 VOC v.3» catalog number K-V017-FRT (from the Centre for Strategic Planning of FMBA, Moscow, Russia).

2.3. SARS-CoV-2 Sequencing

2.3.1. Library Preparation and Sequencing

Whole-genome amplification of SARS-CoV-2 virus genome was performed using ARTIC primer V4 with reverse transcription polymerase chain reaction using BioMaster RT-PCR—Premium (Biolabmix, catalog number RM05-200, Russia). DNA libraries were constructed using NEBNext Fast DNA Fragmentation and Library Prep Set for Ion Torrent (New England Biolabs, Ipswich, MA, USA), according to the manufacturer’s instructions. DNA sequencing was performed using Ion 540 Chip and Ion S5XL System (Thermo Fisher Scientific, Waltham, MA, USA).

2.3.2. NGS Data Analysis

The analysis of NGS data was performed, as previously described [22]. Briefly, the quality control of raw data was performed using FastQC [23]. The ARTIC primers were trimmed using cutadapt v3.1 [24]. The trimming of reads by a quality filter was performed using vsearch v2.17.0, and reads smaller than 50 nt were discarded [25]. Trimmed reads were aligned to SARS-CoV-2 Wuhan-Hu-1 (Accession ID GenBank MN908947.3) using BWA-MEM v0.7.17-r1188 [26]. Variant calling and generating consensus sequences were performed using FreeBayes v1.3.5 [27], bcftools v1.12 [28], and bedtools v2.30.0 [29]. Regions with less than 10-fold coverage were masked. Lineages were assigned with Pangolin v3.1.17 using pangoLearn module 2021-12-06. All sequencing data are available online (EPI_ISL_11864996-EPI_ISL_11865125, EPI_ISL_11872910).

2.4. Immunological Analysis

IgG antibodies against the nucleocapsid protein (NC) and receptor-binding domain (RBD) were measured using an enzyme immunoassay. A commercial test system «SARS-CoV-2-RBD-ELISA-Gamaleya» (from N.F. Gamaleya NRCEM, Moscow, Russia) was used by the manufacturer’s instructions to detect IgG in the receptor-binding domain of surface glycoprotein S (spike).

To interpret the results, the cut-off value of the optical density OD (cut-off) was calculated using the following formula:

cut-off = OD_K− + 3 × SD_K−

where:

• –. SD_K−—mean squared deviation of the optical density of the negative control.

The test serum sample was evaluated as positive if its OD value (A450 mod.) was 2.5 times higher than the cut-off value.

IgG detection of the SARS-CoV-2 (NC) nucleocapsid was performed using an “in-house” ELISA test system. NC was cloned, expressed in E. coli, and purified in house. For coating into 96-well ELISA plates (Costar 2592, Corning, NY, USA), 100 μL of a solution of 1 μg/mL of recombinant protein was added to a phosphate-salt buffer (PBS pH 7.4) and incubated overnight at +4 °C. The next day, the plates were blocked for 2 h at room temperature with a S002X buffer containing 0.5% casein (Xema, Moscow, Russia).

The samples were pre-diluted 10 times with an ELISA buffer S011 (Xema, Moscow, Russia). Convalescent sera samples of patients with COVID-19 confirmed by RT-PCR were used as a positive control, and archived sera samples collected before 2020 were used as a negative control.

Before using the ELISA plate, the blocking buffer S002X was removed. The diluted samples were added into the wells of the ELISA plate and incubated for 1 h at 37 °C. Following incubation, the plate was washed 3 times with PBS containing 0.1% Twin-20 and incubated with 100 μL of HRP-conjugated anti-human IgG (Novex A18823; USA) for 1 h at 37 °C. After incubation, the plate was washed 6 times and the reaction was visualized using a chromogen-substrate solution R055 (Xema, Russia). After 10 min, the reaction was stopped with a 10% HCl and the OD of the substrate mi xture was measured on a spectrophotometer at a wavelength of 450 nm (Multiscan FC, Thermo Scientific, MA, USA).

To interpret the results, the cut-off value of the OD was calculated using the following equation:

Cut off = X + 0.2

where:

• X is the mean OD value of negative control samples 1 and 2.

For each sample, the positivity coefficient (PC) was calculated using the following formula:

PC = OD sample/cut-off

• PC < 0.9—negative, PC > 1.1—positive.

As an additional immunological analysis, the Mindray CL-900i^® anti-SARS-CoV-2 IgG (Shenzhen Mindray Bio-Medical Electronics Co., Shenzhen, China) test was used according to the manufacturer’s instructions.

2.5. Statistical Analysis

The algorithm for initial data processing and cleaning is shown in Figure S1.

Since the quantitative parameters in most groups varied from the normal distribution (the Shapiro–Wilk and Kolmogorov–Smirnov tests p < 0.05), nonparametric parameters were used to describe the median measures of distributions (Me (Q1–Q3)). To compare the parameters in the groups, The Mann–Whitney test or the Kruskal–Wallis test with subsequent posterior comparisons with Bonferroni correction were used.

To compare the qualitative parameters, Pearson’s chi-squared test and Fisher’s exact test were used (explanations are given in the text and captions of the figures). When comparing monopole tables, posteriori comparisons with Benjamin–Hochberg multiplicity correction were used.

The epidemiological effectiveness of vaccination was evaluated by the formula VE = (1 − RR) × 100%, where RR is the risk ratio of becoming sick with various severity in the study groups. Since it was not possible to obtain data on people who were not ill, we used the ratios of vaccinated and ill people in Moscow from open data sources on 26 January 2022 [30,31]. To simplify access to the statistical data from these sources, we converted the Russian language references data to English, and submitted the translation as Supplementary Materials Tables S10 and S11. The 95% confidence interval was calculated according to the previously described method [32].

For statistical analysis, we used R with RStudio (Boston, MA, USA) and SPSS Statistics v.26 (IBM, Armonk, NY, USA).

3. Results

We performed the random sampling of patients at the large clinical center in Moscow to assess the effectiveness of the Sputnik Light and Sputnik V vaccines in protection from hospitalization during the prevalence of Omicron. This also included people who received booster doses. Patients were recruited from 11 January to 21 February during the period of increasing COVID-19 incidence in the context of the emerging Omicron variant in Moscow. The only inclusion criterion was a positive COVID-19 PCR test upon admission. For all patients included in the study, nasopharyngeal swabs and blood serum samples were obtained at the time of admission. The severity of COVID-19 was determined on the WHO scale at the time of the analysis and subsequent severity monitoring. PCR analysis to confirm COVID-19 infection was performed. Subsequently, PCR analysis to detect Omicron was also conducted. Then, some samples were sequenced using NGS. These samples confirmed the diagnosis and determined the SARS-CoV-2 variant. Levels of IgG antibodies of NC and RBD were determined for sera samples using ELISA. An impersonal database was formed based on the results of patient recruitment. It included demographic information, vaccination status, laboratory parameters, and the severity of the disease.

The created database with the demographic data of patients was used to analyze the clinical and laboratory data of the patients. The initial sample was 1112 patients. Yet, due to the lack of some data, subsamples were formed for each analysis (Figure S1; Table S1). In total, 240 patients with mild severity, 420—moderate severity, 60—severe, and 78 in critical condition in the intensive care unit (ICU) at the time of performing the analysis were included in the study. The severity groups were homogeneous by gender (p = 0.849, Pearson’s chi-squared test), with a slight predominance of women (95% CI 53.8–58.8%) (Table 1).

Mild, severe, and critical groups had no statistically significant differences in age. The median age was 72–74 years. However, the group of moderate course patients included patients that were slightly younger. The median age was 69 years. At the same time, there were more vaccinated patients in this group in comparison to others (n = 119, 29.7%). The minimum proportion of vaccinated patients was in the group with patients in a critical condition at the time of sampling (n = 11, 14.7%). The differences in the proportions of vaccinated patients among the moderate and critical groups turned out to be statistically significant (p = 0.044, chi-squared test following the Benjamin–Hochberg multiplicity correction). Patients who had received other preventive drugs were excluded from the proportion calculations of vaccinated patients.

The age of vaccinated and unvaccinated participants in the study did not vary statistically (p = 0.058, Mann–Whitney test). There were more women in the group of unvaccinated patients (p = 0.008, Pearson’s chi-squared test) (Table S2). At the same time, a more detailed analysis of the vaccinated group shows that the age of those with two vaccine doses was significantly lower than in the unvaccinated groups or those who received 3–4 doses (Table S3). At the same time, the number of men in the groups increased, depending on the number of doses received. It reached its maximum in the group with 3–4 doses.

From the sample, we obtained information about re-infection cases that was available for 547 people. Only 9 people from these 547 (1.6%) cases were re-infected. Of these, five were in the group of unvaccinated individuals. One patient had previously received Sputnik Light. Three had received Sputnik V. Considering such a low proportion of patients who previously recovered from COVID-19, a special account of the earlier disease was not conducted in further calculations.

3.1. Dominant Variants of the Virus and Their Features

For the patients’ enrollment, a PCR analysis was performed using a test system that detected all SARS-CoV-2 variants. To determine the dominant virus variant, in addition to the standard PCR analysis, a test system differentiating between Delta and Omicron was used. Some of the samples with a sufficient viral load (CT below 30) were sequenced. Then, a whole-genome analysis was performed. A differentiating PCR system was used for 620 samples. The severity of the diseases was known in 506 cases from 620, at the time of the analysis. In 506 samples, with the result of a specific PCR analysis, the Delta variant was found in 41 patients, which was 8.1%. At the same time, in the group of patients in a critical condition, the proportion of patients infected with the Delta virus variant was 22.2% (n = 10). Table 2 shows a tendency of the increasing proportion of Delta variant infection cases that corresponds to an increase in the severity of the disease. Statistically significant differences were achieved when comparing mild and moderate cases with critical ones.

The viral load for the Omicron variant tended to decrease in the range from mild to severe forms of the disease. It increased slightly for patients in a critical condition, regardless of the vaccination status. The median viral load values were 5.51 × 10⁵, 8.71 × 10⁴, 2.42 × 10³, and 9.63 × 10⁵ copies/mL in the case of mild, moderate, severe, and critical forms, respectively. Statistically significant differences were found between mild and moderate (p = 0.010, with Bonferroni multiplicity correction), light and severe (p = 0.001, with Bonferroni multiplicity correction), and between severe and critical (p = 0.019, with Bonferroni multiplicity correction) forms. The results are presented in Table 3 and Figure 1.

A higher viral load was noted for the Delta variant. However, statistically significant differences were found only for moderate cases. The viral load of the Delta variant in this case was two orders of magnitude higher than for the Omicron variant (the median viral load for the Delta variant was 1.37 × 10⁶, p = 0.010, Mann–Whitney test). The Delta variant is characterized by a similar decrease in viral load. Such a decrease corresponds to an increase in severity from mild to severe and an increase in the quantity of critical cases, although, statistically, these differences were not significant (p = 0.184, Kruskal–Wallis test). This may be due to the insufficient sample size (Table 3). A significant decrease in viral load in vaccinated patients, compared to unvaccinated patients, was not observed. This applied to those infected with the Delta variant (p = 0.111, Mann–Whitney test) and those infected with the Omicron strain (p = 0.528, Mann–Whitney test) (Table 4).

We performed whole-genome sequencing of 131 samples to confirm the results of the PCR analysis for samples with a sufficient viral load. According to the obtained data, Omicron was 79.4% in the structure of all variants (n = 104) (Figure 2). Line BA.1.1 was the most represented among the subsidiary lines of Omicron, which occurred in 69.2% of cases (n = 72). As for PCR, there was a tendency of an increasing proportion of the Delta variant corresponding to the severity of the course (Figure 3, Table 5). Yet, some discrepancy in numbers between the results of the PCR and sequencing are obviously related to the differences in the viral loads of Omicron and Delta. This results in an increase in the latter in genomic sequencing results.

3.2. Features of the Immune Response Depending on the Vaccination Status

The ELISA method was used to study the level of IgG antibodies to NC and RBD in different groups, depending on the vaccination status and severity of COVID-19 (Table 4). In unvaccinated patients, there was a tendency of increasing the number of antibodies in both the NC protein and the RBD domain of the SARS-CoV-2 S protein with an increase in the severity of the course of the disease (p < 0.001 and p = 0.007, respectively). This trend was more pronounced for NC than RBD. In vaccinated patients, there was also a tendency for IgG levels to increase in accordance with the increasing severity of the disease from mild to severe with a further decline in the case of a critical course (Figure 2).

There was a decrease in the number of antibodies, even with a larger number of vaccine components administered to the group of vaccinated patients in a critical condition. This decrease may be due to the time that elapsed after the administration of the final vaccine component (Table S6). Thus, there is a tendency to increase the median time between the last component and the date of hospitalization. The median time escalated from 132 days (4.4 months) for mild cases to 183 (6.0 months) for critical cases. In the case of a critical course of the disease, at least 75% of patients received the final dose of the vaccine more than 120 days (3.9 months) before hospitalization. Statistically, these differences were not significant p > 0.05 (p = 0.519, Kruskal–Wallis test).

We did not observe any differences in the level of antibodies for patients infected with different variants of the virus (Table S7). The described patterns regarding the increase in the level of antibodies depending on the severity of COVID-19 and the number of vaccine doses administered were confirmed using an alternative immunochemiluminescent analysis of the amount of IgG antibodies to SARS-CoV-2 (Table S8).

3.3. Effectiveness of Vaccination with Sputnik Drugs

To assess the effectiveness of the Sputnik vaccine, we used the case–population method that had been recommended by the WHO for situations in which the proportion of vaccinated does not significantly change over time in the population [33]. Data for the number of unvaccinated, vaccinated, and re-vaccinated were used from the open registry.

Epidemiological efficacy was evaluated by the formula VE = (1 − RR) × 100%, where RR is the risk ratio of getting sick with various severity in the study groups. The proportion of vaccinated and non-immune individuals in the population of Moscow was calculated to analyse the epidemiological efficacy of vaccinations (Table 5, Figure 3).

Thus, even though the effectiveness of vaccinations against hospitalization decreased in comparison to the time when other strains were predominant [1,16,17,18,19,20,21], according to our data, it remained at the level of 85.9% (95% CI 83.0–88.0%) for vaccinated patients with at least one component (Table 5). For those who were vaccinated with at least two doses, the epidemiological efficacy rose to 87.6% (95% CI 85.4–89.5%), whereas for those vaccinated with 3 or 4 components (re-vaccination with Sputnik Light or Sputnik V after Sputnik V), it was 97.0% (95% CI 95.9–97.8%).

A lower number of severe course cases among vaccinated people may be due to the higher frequency of the Delta virus variant in this group. We considered that the effectiveness of vaccination against the Delta variant may be greater. Therefore, we conducted a separate analysis of effectiveness. Samples with the Delta virus variant were excluded from the calculation (Table S9). The analysis shows the low impact of the Delta variant residual proportion on the calculation of effectiveness.

4. Discussion

According to the international Sputnik vaccine supplier RDIF, the vaccine is registered for use in 71 countries around the world [34]. In the Russian Federation, Argentina, Hungary, Bahrain, and San Marino, the Sputnik vaccine is one of the main vaccines used for COVID-19 prophylaxis. The results of the clinical trials, as well as their effectiveness estimation in field studies, laced Sputnik among the most effective drugs [18,19]. Previously, a very moderate decrease in viral neutralizing activity was shown in relation to the main variants that are of concern. These variants included Alpha, Beta, Gamma, and Delta [35]. Experiments on animal models also did not reveal a significant efficiency drop [36]. This reflects the stability of the protective parameters regarding the main variants of concern.

Meanwhile, the decrease in virus neutralizing activity in matters of the Omicron variant turned out to be significantly more pronounced. Patients, vaccinated with Sputnik V after 6–12 months, had 12 times less virus neutralizing capacity of sera against Omicron compared to original variant B.1.1.1 [16]. At the same time, those who received the vaccine from 0 to 3 months and 3 to 6 months had a virus neutralizing titer detected in 80% and 68.75% of cases, respectively [37]. In turn, the presence of previous infection and booster doses among vaccinated patients increased the proportion of people with a virus neutralizing titer against Omicron to 100% [16,37]. At the same time, it remained unclear to what extent virus neutralization correlates with field effectiveness. Thus, field experiments needed to validate in vitro data obtained using virus-neutralization analysis.

In the case of the current analysis, the virus genetic study showed that the evaluation of the effectiveness mainly relates to the Omicron variant during the period of study (from 11 January to 21 February). More than 90% of infection cases in patients included in the study were caused by the Omicron variant. The composition analysis of the Omicron line showed that 37% was the BA.1 variant, 68%—BA.1.1, and 4%—BA.2, according to the whole-genome sequencing of SARS-CoV-2. The proportion of the Delta variant increased with the severity of the COVID-19 course. While the proportion of the Delta variant among patients with a mild course of the disease was 4.9%, it reached 22.2% (4.5 times more often) among critical-condition patients. This corresponds with the data for the Omicron variant severity. In this data, it is several times milder than the Delta variant [38].

The data on the assessment of IgG antibodies against RBD show that, in all cases with different disease severities, the values in the group of patients vaccinated with two doses were significantly higher than those not vaccinated. Among the vaccinated patients, there was a trend towards an increase in the level of antibodies in RBD corresponding with the increasing severity of the disease, from mild to severe, with a further decline in critical cases (Figure 2). In the case of mild and moderate patients, an increase in the level of antibodies in RBD was observed, depending on the number of components that were administered. This indicates a correlation between the number of vaccine shots and the level of IgG RBD virus-neutralizing antibodies. In this regard, it is expected that an increase in the epidemiological efficacy of the vaccine against the Omicron variant corresponds to an increase in the number of vaccine doses received. Indeed, the epidemiological efficacies regarding protection against hospitalization of patients with COVID-19 in our study were 85.9% (95% CI 83.0–88.0%), 87.6% (95% CI 85.4–89.5%), and 97.0% (95% CI 95.9–97.8%) for those who received more than one, more than two, and more than three injections, respectively. As for the protection from a critical condition requiring ICU, the epidemiological efficacies were 93.2% (95% CI 87.1–96.4%), 94.5% (95% CI 88.9–97.2%), and 99.4% (95% CI 95.6–99.9%) for those who received more than one, more than two, and more than three doses, respectively. The increase in the effectiveness of COVID-19 vaccines in terms of the prevention of more severe forms is well known and characteristic not only of Sputnik, but also of other vaccines, including Pfizer-BioNTech, Moderna, and ChAdOx1 nCoV-19 (AstraZeneca) [18,19]. In the current study, this trend continued for the Omicron variant. This brings hope that a certain level of intense immunity, which is the result of the widespread use of various vaccines based on the original strain’s Spike protein gene, will maintain its effectiveness against new variants of SARS-CoV-2. The emergence of such variants is expected [39]. An alternative to vaccination programs is the isolation and quarantine strategy [40], which is a considerably big-budget plan [41]. Our study indicates that protection against hospitalization and severe forms of COVID-19 remain consistently high for the Sputnik Light and Sputnik V vaccines.

Our research is not without its limitations. In particular, the study is based on observations made in a single medical center. According to Moscow statistics, 10–15% of patients were hospitalized with COVID-19 in the City Clinical Hospital No. 67. Population data for the number of vaccines were obtained from an open register on the date of analysis. Sputnik Light and Sputnik V are the main vaccines used for the prevention of COVID-19 in Russia. The number of doses for other COVID-19 vaccines is more than an order of magnitude lower. However, the register data used did not allow us to include vaccines and booster vaccinations only for Sputnik. In this study, it was also not possible to account for herd immunity and potential undiagnosed cases of both symptomatic and asymptomatic forms of COVID-19, as well as cases of hospitalization with undiagnosed COVID-19. We did not use special methods to optimize our sample. This could result in a slight distortion of the effectiveness calculations. It was also not possible to determine the direct risks of hospitalization for vaccinated and non-vaccinated patients with COVID-19, only the ratio of risks.

5. Conclusions

The study shows that vaccination with Sputnik V and Sputnik Light has a high effectiveness for protection against hospitalization. The reduction in the severity of COVID-19 regarding the Omicron variant was also observed. The greatest effectiveness was evident in protection against a critical course of the disease requiring patients to be admitted to the intensive care unit. An increase in the level of antibodies correlated with an increase in the number of doses, ranging from 1 to 4. As a result, there was a lower chance of hospitalization. This should not preclude the need to create new vaccines aimed at SARS-CoV-2 transmission reduction in the population. This suggests that the effectiveness of original vaccines based on the Wuhan S antigen against future VOCs, in terms of protection against hospitalization with Omicron, can be quite high, at least within the recommended vaccination and re-vaccination regimens.

Author Contributions

Conceptualization, A.S.S., V.A.G., D.Y.L. and A.L.G. (Alexander L. Gintsburg); methodology, D.A.O., N.A.K., A.A.P., D.A.K., D.V.V., S.V.S. and O.E.O.; software, D.A.O. and A.A.P.; validation, D.A.O.; formal analysis, V.A.G., D.A.O. and A.A.P.; investigation, S.V.S., O.E.O., E.N.K., A.V.P., V.G.K., N.V.S., A.L.G. (Alexey L. Gorbunov), E.A.M., A.V.S., L.L.P., V.A.G., D.A.O., N.A.K., A.A.P., D.A.K., D.V.V., E.V.S., D.D.K., S.R.K., I.V.G., A.S.O., E.D.B., R.R.A., A.A.B., L.A.V., I.P.L., A.P.T. and E.A.T.; resources, A.S.S., V.A.G., A.P.T., D.Y.L. and A.L.G. (Alexander L. Gintsburg); data curation, D.A.O., N.A.K. and A.A.P.; writing—original draft preparation, V.A.G. and D.A.O.; writing—review and editing, V.A.G., D.A.O., N.A.K., A.A.P., D.A.K. and D.V.V.; visualization, D.A.O., D.D.K. and A.A.P.; supervision, A.S.S. and V.A.G.; project administration, A.S.S. and V.A.G.; funding acquisition, A.S.S., V.A.G., D.Y.L. and A.L.G. (Alexander L. Gintsburg). All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Written informed consent was obtained from all subjects in accordance with the order of the Ministry of Health of the Russian Federation of 21 July 2015 N 474n. All samples were de-identified prior to receipt by the research team. The study was submitted to the Local Ethics Committee of the Gamaleya Center. The committee concluded that the study did not use identifiable biological samples and did not provide any confidential information (decision #14 29 September 2021).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All sequence data are available online (EPI_ISL_11864996-EPI_ISL_11865125, EPI_ISL_11872910). Inquiries about access to the original clinical data should be directed to the corresponding authors.

Acknowledgments

We would like to thank the study participants.

Conflicts of Interest

D.Y.L. and A.L.G. (Alexander L. Gintsburg) report patents for a Sputnik V immunobiological expression vector, for a pharmaceutical agent, and for its method of use to prevent COVID-19. The funders had no role in the design and conduct of the study; in the collection, management, analysis, and interpretation of the data; in the preparation, review, or approval of the manuscript; or in the decision to submit the manuscript for publication. All other authors declare no competing interests.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Tables

Figure 1. Viral load (copies/mL) depending on the severity of the disease and the strain. (A) Omicron strain; (B) Delta strain.

Figure 2. Levels of IgG antibodies in the RBD domain and NC protein depending on severity. Two groups are presented—unvaccinated and vaccinated with two components of Sputnik V. The Y-axis indicates the number of antibodies (PC), and shows the median, upper, and lower quantiles. The unvaccinated samples are characterized by a wide variation in the level of IgG antibodies in the RBD domain of the S protein. (A)—shows the level of antibodies IgG RBD in unvaccinated patients, (B)—the level of IgG NC in unvaccinated patients, (C)—the level of IgG RBD in vaccinated patient and (D)—the level of IgG NC in vaccinated patients.

Figure 3. Effectiveness of vaccination in relation to the varying severity of the disease among hospitalized patients, considering the number of vaccine doses received. (A) Protection from any cases of hospitalization, (B) protection from moderate-severe and more severe cases, (C) protection from severe and critical cases, and (D) protection from critical cases (VE 95% CI).

Table 1

Demographic analysis of patients (gender, age, and severity are known; n = 798).

	Mild, n = 240	Moderate, n = 420	Severe, n = 60	Critical, n = 78	p Value
GenderWomen, n, %Men, n, %	136 (56.7%)104 (43.3%)	247 (58.8%)173 (41.2%)	35 (58.3%)25 (41.7%)	42 (53.8%)36 (46.2%)	0.849 (Pearson’s chi-squared test)
Age Me(Q1–Q3)	74 (64–83)	69 (59–76)	72 (63.5–82)	73 (66–82)	<0.001 *p24 = 0.019 *p21 < 0.001 *
Vaccinated (starting with one component)	Taken into account 23667 (28.4%)	Taken into account 401119 (29.7%)	Taken into account 5913 (22.0%)	Taken into account 7511 (14.7%)	0.043 * (Pearson’s chi-squared test)p24 = 0.044 *
Not considered due to not being vaccinated with Sputnik V	4	19	1	3

*—The differences are statistically significant (p < 0.05).

Table 2

Genetics of the virus depending on the severity of COVID-19.

	Mild, n = 142	Moderate, n = 284	Severe, n = 35	Critical, n = 45	p
Delta variant	7 (4.9%)	20 (7.0%)	4 (11.4%)	10 (22.2%)	0.004 * (Fisher’s exact test)p14 = 0.002 *p24 = 0.003
Omicron variant	135 (95.1%)	264 (93.0%)	31 (88.6%)	35 (77.8%)

*—The differences are statistically significant (p < 0.05).

Table 3

Viral load depending on severity (n = 506).

	Mild, n = 142	Moderate, n = 284	Severe, n = 35	Critical, n = 45	p Value
Viral load of the Omicron variant
PCR +	127 (94.1%)	247 (93.6%)	31 (100.0%)	32 (91.4%)
Viral load, CtMe(Q1–Q3)	26.39 (24.23–31.42)	29.23 (25.78–33.06)	34.53 (27.57–36.76)	26.66 (24.00–32.67)	<0.001 * (the Kruskal–Wallis test)p12 = 0.009 *p13 < 0.001 *p43 = 0.020 *
Viral load, copies/mLMe(Q1–Q3)	5.51 × 105 (1.25 × 104–3.03 × 106)	8.71 × 104 (4.64 × 103–7.43 × 105)	2.42 × 103 (2.43 × 102–4.80 × 105)	9.63 × 105(9.78 × 103–3.73 × 106 )	<0.001 * (the Kruskal–Wallis test)p12 = 0.010 *p13 = 0.001 *p43 = 0.019 *
Viral load of the Delta variant
PCR+	7 (100.0%)	19 (95.0%)	4 (100.0%)	10 (100.0%)
Viral load, CtMe(Q1–Q3)	25.16 (21.02–25.54)	26.22 (24.70–28.63)	30.22 (27.58–31.60)	26.13 (25.23–27.85)	0.063 (the Kruskal–Wallis test)
Viral load, copies/mLMe(Q1–Q3)	1.79 × 106 (1.37 × 106–3.92 × 106)	1.37 × 106 (1.31 × 105–2.00 × 106)	4.71 × 104 (1.33 × 104–9.64 × 105)	1.38 × 106(9.93 × 104–2.09 × 106)	0.184 (the Kruskal–Wallis test)
p (differences between various strains)	0.094 (The Mann–Whitney test)	0.010 *(The Mann–Whitney test)	0.352 (The Mann–Whitney test)	0.805 (The Mann–Whitney test)

*—The differences are statistically significant (p < 0.05).

Table 4

Features of the immune response in immunized groups depending on severity. Only unvaccinated and vaccinated patients with at least one component were studied (n = 749). Immunological analysis of the coefficient of positivity (PC).

IgG NC	Mild	Moderate	Severe	Critical	p
UnvaccinatedMe(Q1–Q3)	n = 1680.86 (0.24–5.97)	n = 2772.16 (0.38–8.77)	n = 464.61 (1.05–11.09)	n = 657.22 (1.33–12.18)	<0.001 * (the Kruskal–Wallis test)p12 = 0.017 *p13 = 0.001 *p14 < 0.001 *p24 = 0.002 *
Vaccinated with 1 componentMe(Q1–Q3)	n = 50.63 (0.19–0.75)	n = 172.23 (0.20–12.57)	n = 1--	n = 25.03 (3.29–6.76)	0.593 (the Kruskal–Wallis test)
Vaccinated with 2 componentsMe(Q1–Q3)	n = 400.59 (0.21–3.05)	n = 830.68 (0.22–3.51)	n = 92.42 (1.62–6.94)	n = 81.36 (0.46–13.11)	0.306 (the Kruskal–Wallis test)
Vaccinated with 3–4 componentsMe(Q1–Q3)	n = 220.45 (0.23–4.72)	n = 190.37 (0.16–1.79)	n = 30.13 (0.11–0.98)	n = 1	0.357 (the Kruskal–Wallis test)
p	0.682	0.008 *p02 = 0.044 *	0.154	0.788
IgG RBD	Mild	Moderate	Severe	Critical	p
UnvaccinatedMe(Q1–Q3)	n = 1681.17 (0.64–5.9)	n = 2771.08 (0.69–3.22)	n = 462.66 (1.03–17.79)	n = 651.89 (0.80–11.93)	0.007 * (the Kruskal–Wallis test)p13 =0.008 *p23 = 0.003 *p24 = 0.041 *
Vaccinated with 1 componentMe(Q1–Q3)	n = 532.19 (3.18–72.02)	n = 173.10 (1.68–24.90)	n = 1--	n = 269.49 (56.29–82.69)	0.090 (the Kruskal–Wallis test)
Vaccinated with 2 componentsMe(Q1–Q3)	n = 4029.07 (6.19–57.00)	n = 8344.86 (5.78–83.11)	n = 993.00 (7.06–95.83)	n = 821.41 (4.44–57.70)	0.265 (the Kruskal–Wallis test)
Vaccinated with 3–4 componentsMe(Q1–Q3)	n = 2238.54 (14.56–59.07)	n = 1953.01 (36.39–79.26)	n = 344.52 (28.66–45.49)	n = 1	0.357 (the Kruskal–Wallis test)
p	<0.001 *p01 = 0.046 *p02 < 0.001 *p03 < 0.001 *	<0.001 *p02 < 0.001 *p03 < 0.001 *	0.003 *p02 = 0.012 *	0.018 *posteriori comparisons revealed no differences

*—The differences are statistically significant (p < 0.05).

Table 5

Analysis of the epidemiological efficacy of vaccination in protection against the severe course of COVID-19 in hospitalized patients.

Quantity of 18+ People in Moscow 10,600,000	Quantity of Vaccinated6,000,000	Quantity of Previously ill 2,232,836	Quantity of Non-Immune2,367,164
	Patients with at least 1 component of vaccine, nVE % (95% CI)	Patients with at least 2 components of vaccine,nVE % (95% CI)	Patients with at least 3–4 components of vaccine,n VE % (95% CI)	Infected and not vaccinated, n
Patients with mild and more severe forms of the disease	21085.9%(95% CI 83.0–88.0%)	18587.6%(95% CI 85.4–89.5%)	4597.0%(95% CI 95.9–97.8%)	561
Patients with moderate and more severe forms of the disease	14385.6%(95% CI 82.6–88.1%)	12387.6%(95% CI 84.8–89.9%)	2397.7%(95% CI 96.5–98.5%)	392
Patients with severe and critical forms of the disease	2491.4%(95% CI 86.6–94.5%)	2192.5%(95% CI 88.0–95.3%)	498.6%(95% CI 96.1–99.5%)	110
Patients in critical condition	1193.2%(95% CI 87.1–96.4%)	994.5%(95% CI 88.9–97.2%)	199.4%(95% CI 95.6–99.9%)	64

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/vaccines10060938/s1. Figure S1: Algorithm of data pre-processing; Figure S2: The general distribution of proportions of genetic lines according to the results of SARS-CoV-2 genome sequencing; Figure S3: Distribution of proportions of genetic lines by severity, according to the results of SARS-CoV-2 genome sequencing; Table S1: Patients, included in the study; Table S2: Demographic characteristics according to the vaccination status; Table S3: Demographic characteristics of hospitalized patients according to the number of vaccine doses received; Table S4: Viral load of all severity groups (including patients with unidentified course of the disease) depending on the strain and vaccination status; Table S5: Distribution of SARS-CoV-2 genetic variants by patients’ severity; Table S6: The time elapsed since the injection of the last Gam-Covid-Vac component before hospitalization with COVID-19; Table S7: Viral load and immune response of various virus variants depending on the vaccination status (those vaccinated with drugs other than Gam-Covid-Vac were not included, excluded-4); Table S8: Features of the immune response in immunized groups depending on the severity. Only those who were not vaccinated and those vaccinated with at least one “Gam-Covid-Vac” component were studied (n = 749). Immunological analysis for IgG antibodies (Mindray); Table S9: Analysis of the epidemiological efficacy of vaccination in protection against the severe course of COVID-19 in hospitalized patients; Table S10: The number of people vaccinated against COVID-19 in Moscow; Table S11: Coronavirus statistics in Moscow.

References

1. Logunov, D.Y.; Dolzhikova, I.V.; Shcheblyakov, D.V.; Tukhvatulin, A.I.; Zubkova, O.V.; Dzharullaeva, A.S.; Kovyrshina, A.V.; Lubenets, N.L.; Grousova, D.M.; Erokhova, A.S. et al. Safety and Efficacy of an rAd26 and rAd5 Vector-Based Heterologous Prime-Boost COVID-19 Vaccine: An Interim Analysis of a Randomised Controlled Phase 3 Trial in Russia. Lancet; 2021; 397, pp. 671-681. [DOI: https://dx.doi.org/10.1016/S0140-6736(21)00234-8]

2. Baden, L.R.; El Sahly, H.M.; Essink, B.; Kotloff, K.; Frey, S.; Novak, R.; Diemert, D.; Spector, S.A.; Rouphael, N.; Creech, C.B. et al. Efficacy and Safety of the mRNA-1273 SARS-CoV-2 Vaccine. N. Engl. J. Med.; 2021; 384, pp. 403-416. [DOI: https://dx.doi.org/10.1056/NEJMoa2035389] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33378609]

3. Polack, F.P.; Thomas, S.J.; Kitchin, N.; Absalon, J.; Gurtman, A.; Lockhart, S.; Perez, J.L.; Pérez Marc, G.; Moreira, E.D.; Zerbini, C. et al. Safety and Efficacy of the BNT162b2 mRNA COVID-19 Vaccine. N. Engl. J. Med.; 2020; 383, pp. 2603-2615. [DOI: https://dx.doi.org/10.1056/NEJMoa2034577] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33301246]

4. Voysey, M.; Clemens, S.A.C.; Madhi, S.A.; Weckx, L.Y.; Folegatti, P.M.; Aley, P.K.; Angus, B.; Baillie, V.L.; Barnabas, S.L.; Bhorat, Q.E. et al. Safety and Efficacy of the ChAdOx1 nCoV-19 Vaccine (AZD1222) against SARS-CoV-2: An Interim Analysis of Four Randomised Controlled Trials in Brazil, South Africa, and the UK. Lancet; 2021; 397, pp. 99-111. [DOI: https://dx.doi.org/10.1016/S0140-6736(20)32661-1]

5. Butt, A.A.; Omer, S.B.; Yan, P.; Shaikh, O.S.; Mayr, F.B. SARS-CoV-2 vaccine effectiveness in a high-risk national population in a real-world setting. Ann. Intern. Med.; 2021; 174, pp. 1404-1408. [DOI: https://dx.doi.org/10.7326/M21-1577]

6. Pawlowski, C.; Lenehan, P.; Puranik, A.; Agarwal, V.; Venkatakrishnan, A.J.; Niesen, M.J.M.; Venkatakrishnan, A.J.; Niesen, M.J.M.; O’Horo, J.C.; Virk, A. et al. FDA-authorized mRNA COVID-19 vaccines are effective per real-world evidence synthesized across a multi-state health system. Med; 2021; 2, pp. 979-992.E8. [DOI: https://dx.doi.org/10.1016/j.medj.2021.06.007]

7. Jara, A.; Undurraga, E.A.; González, C.; Paredes, F.; Fontecilla, T.; Jara, G.; Pizarro, A.; Acevedo, J.; Leo, K.; Leon, F. et al. Effectiveness of an inactivated SARS-CoV-2 vaccine in Chile. N. Engl. J. Med.; 2021; 385, pp. 875-884. [DOI: https://dx.doi.org/10.1056/NEJMoa2107715]

8. Al Kaabi, N.; Zhang, Y.; Xia, S.; Yang, Y.; Al Qahtani, M.M.; Abdulrazzaq, N.; Al Nusair, M.; Hassany, M.; Jawad, J.S.; Abdalla, J. et al. Effect of 2 inactivated SARS-CoV-2 vaccines on symptomatic COVID-19 infection in adults: A randomized clinical trial. JAMA; 2021; 326, e218565. [DOI: https://dx.doi.org/10.1001/jama.2021.8565]

9. Mahase, E. COVID-19: Moderna applies for US and EU approval as vaccine trial reports 94.1% efficacy. BMJ; 2020; 371, m4709. [DOI: https://dx.doi.org/10.1136/bmj.m4709]

10. Mohapatra, R.K.; Kuppili, S.; Kumar Suvvari, T.; Kandi, V.; Behera, A.; Verma, S.; Kudrat-E-Zahan Biswal, S.K.; Al-Noor, T.H.; El-Ajaily, M.M.; Sarangi, A.K. et al. SARS-CoV-2 and its variants of concern including Omicron: A never ending pandemic. Chem. Biol. Drug Des.; 2022; 99, pp. 769-788. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35184391][DOI: https://dx.doi.org/10.1111/cbdd.14035]

11. Jung, C.; Kmiec, D.; Koepke, L.; Zech, F.; Jacob, T.; Sparrer, K.M.J.; Kirchhoff, F. Omicron: What makes the latest SARS-CoV-2 variant of concern so concerning?. J. Virol.; 2022; 96, jvi0207721. [DOI: https://dx.doi.org/10.1128/jvi.02077-21] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35225672]

12. Dong, E.; Du, H.; Gardner, L. An Interactive Web-Based Dashboard to Track COVID-19 in Real Time. Lancet Infect. Dis.; 2020; 20, pp. 533-534. Available online: https://coronavirus.jhu.edu/map.html (accessed on 24 April 2022). [DOI: https://dx.doi.org/10.1016/S1473-3099(20)30120-1]

13. Mahase, E. Covid-19: Booster dose reduces infections and severe illness in over 60s, Israeli study reports. BMJ; 2021; 374, n2297. [DOI: https://dx.doi.org/10.1136/bmj.n2297] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34535470]

14. Mahase, E. COVID-19: Antibody boost after third dose varies greatly by vaccine, study finds. BMJ; 2021; 375, n3011. [DOI: https://dx.doi.org/10.1136/bmj.n3011] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34862175]

15. Logunov, D.Y.; Dolzhikova, I.V.; Zubkova, O.V.; Tukhvatulin, A.I.; Shcheblyakov, D.V.; Dzharullaeva, A.S.; Grousova, D.M.; Erokhova, A.S.; Kovyrshina, A.V.; Botikov, A.G. Safety and immunogenicity of an rAd26 and rAd5 vector-based heterologous prime-boost COVID-19 vaccine in two formulations: Two open, non-randomised phase 1/2 studies from Russia. Lancet; 2020; 396, pp. 887-897. [DOI: https://dx.doi.org/10.1016/S0140-6736(20)31866-3]

16. González, S.; Olszevicki, S.; Salazar, M.; Calabria, A.; Regairaz, L.; Marín, L.; Campos, P.; Varela, T.; Martínez, V.V.G.; Ceriani, L. et al. Efficacy of the first component of Gam-COVID-Vac (Sputnik V) on reduction of SARS-CoV-2 confirmed infections, hospitalisations and mortality in patients aged 60-79: A retrospective cohort study in Argentina. EClinicalMedicine; 2021; 40, 101126. [DOI: https://dx.doi.org/10.1016/j.eclinm.2021.101126]

17. Dolzhikova, I.V.; Gushchin, V.A.; Shcheblyakov, D.V.; Tsybin, A.N.; Shchetinin, A.M.; Pochtovyi, A.A.; Komissarov, A.B.; Kleymenov, D.A.; Kuznetsova, N.A.; Tukhvatulin, A.I. et al. One-shot immunization with Sputnik Light (the first component of Sputnik V vaccine) is effective against SARS-CoV-2 Delta variant: Efficacy data on the use of the vaccine in civil circulation in Moscow. medRxiv; 2021; [DOI: https://dx.doi.org/10.1101/2021.10.08.21264715]

18. Vokó, Z.; Kiss, Z.; Surján, G.; Surján, O.; Barcza, Z.; Pályi, B.; Formanek-Balku, E.; Molnár, G.A.; Herczeg, R.; Gyenesei, A. et al. Nationwide effectiveness of five SARS-CoV-2 vaccines in Hungary-the HUN-VE study. Clin. Microbiol. Infect.; 2022; 28, pp. 398-404. [DOI: https://dx.doi.org/10.1016/j.cmi.2021.11.011]

19. AlQahtani, M.; Bhattacharyya, S.; Alawadi, A.; Al Mahmeed, H.; Al Sayed, J.; Justman, J.; El-Sadr, W.M.; Hidary, J.; Mukherjee, S. Morbidity and mortality from COVID-19 post-vaccination breakthrough infections in association with vaccines and the emergence of variants in Bahrain. Res. Sq.; 2021; preprint [DOI: https://dx.doi.org/10.21203/rs.3.rs-828021/v1]

20. Barchuk, A.; Cherkashin, M.; Bulina, A.; Berezina, N.; Rakova, T.; Kuplevatskaya, D.; Stanevich, O.; Skougarevskiy, D.; Okhotin, A. Vaccine Efficacy against Referral to Hospital and Severe Lung Injury Associated with COVID- 19: A Population-Based Case- Control Study in St. Petersburg, Russia. medRxiv; 2022; [DOI: https://dx.doi.org/10.1101/2021.07.31.21261428]

21. Gushchin, V.A.; Tsyganova, E.V.; Ogarkova, D.A.; Adgamov, R.R.; Shcheblyakov, D.V.; Glukhoedova, N.V.; Zhilenkova, A.S.; Kolotii, A.G.; Zaitsev, R.D.; Logunov, D.Y. et al. Sputnik V protection from COVID-19 in people living with HIV under antiretroviral therapy. EClinicalMedicine; 2022; 46, 101360. [DOI: https://dx.doi.org/10.1016/j.eclinm.2022.101360] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35340627]

22. Klink, G.V.; Safina, K.; Nabieva, E.; Shvyrev, N.; Garushyants, S.; Alekseeva, E.; Komissarov, A.B.; Danilenko, D.M.; Pochtovyi, A.A.; Divisenko, E.V. et al. The rise and spread of the SARS-CoV-2 AY.122 lineage in Russia. Virus Evol.; 2022; 8, veac017. [DOI: https://dx.doi.org/10.1093/ve/veac017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35371558]

23. FastQC: A Quality Control Analysis Tool for High throughput Sequencing Data. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 7 April 2022).

24. Martin, M. Cutadapt Removes Adapter Sequences from High-throughput Sequencing Reads. EMBnet. J.; 2011; 17, pp. 10-12. [DOI: https://dx.doi.org/10.14806/ej.17.1.200]

25. Rognes, T.; Flouri, T.; Nichols, B.; Quince, C.; Mahé, F. VSEARCH: A versatile open source tool for metagenomics. PeerJ; 2016; 2016, pp. 1-22. [DOI: https://dx.doi.org/10.7717/peerj.2584]

26. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics; 2009; 25, pp. 1754-1760. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp324] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19451168]

27. Garrison, E.; Marth, G. Haplotype-based variant detection from short-read sequencing. arXiv; 2012; [DOI: https://dx.doi.org/10.48550/arXiv.1207.3907] arXiv: 1207.3907

28. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R. The Sequence Alignment/Map format and SAMtools. Bioinformatics; 2009; 25, pp. 2078-2079. [DOI: https://dx.doi.org/10.1093/bioinformatics/btp352]

29. Quinlan, A.R.; Hall, I.M. BEDTools: A flexible suite of utilities for comparing genomic features. Bioinformatics; 2010; 26, pp. 841-842. [DOI: https://dx.doi.org/10.1093/bioinformatics/btq033]

30. The Number of People Vaccinated against COVID-19 in Moscow. Available online: https://gogov.ru/covid-v-stats/msk#data (accessed on 24 April 2022).

31. Coronavirus Statistics in Moscow. Available online: https://gogov.ru/covid-19/msk#data (accessed on 24 April 2022).

32. Hightower, A.; Orenstein, W.; Martin, S. Recommendations for the use of Taylor series confidence intervals for estimates of vaccine efficacy. Bull. World Health Organ.; 1988; 66, pp. 99-105.

33. Evaluation of COVID-19 Vaccine Effectiveness. WHO Interim Guidance 17 March 2021; Available online: https://www.who.int/publications/i/item/WHO-2019-nCoV-vaccine_effectiveness-measurement-2021.1 (accessed on 24 April 2022).

34. Available online: https://sputnikvaccine.com/rus/about-vaccine (accessed on 24 April 2022).

35. Gushchin, V.A.; Dolzhikova, I.V.; Shchetinin, A.M.; Odintsova, A.S.; Siniavin, A.E.; Nikiforova, M.A.; Pochtovyi, A.A.; Shidlovskaya, E.V.; Kuznetsova, N.A.; Burgasova, O.A. et al. Neutralizing Activity of Sera from Sputnik V-Vaccinated People against Variants of Concern (VOC: B.1.1.7, B.1.351, P.1, B.1.617.2, B.1.617.3) and Moscow Endemic SARS-CoV-2 Variants. Vaccines; 2021; 9, 779.Available online: https://github.com/artic-network/artic-ncov2019/tree/master/primer_schemes/nCoV-2019 (accessed on 24 April 2022). [DOI: https://dx.doi.org/10.3390/vaccines9070779] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34358195]

36. Grousova, D.M.; Dolzhikova, I.V.; Zorkov, I.D.; Ilyukhina, A.A.; Kovyrshina, A.V.; Kuzina, A.V.; Botikov, A.G.; Tukhvatulin, A.I.; Shcheblyakov, D.V.; Izhaeva, F.M. et al. Sputnik V vaccine protects human ACE2-transgenic mice from lethal challenge with SARS-CoV-2 Variants of Concern (VOC: B.1.1.7, B.1.351, B.1.1.28/P.1, B.1.617.2). 2022; in press

37. Lapa, D.; Grousova, D.M.; Matusali, G.; Meschi, S.; Colavita, F.; Bettini, A.; Gramigna, G.; Francalancia, M.; Garbuglia, A.R.; Girardi, E. et al. Retention of Neutralizing Response against SARS-CoV-2 Omicron Variant in Sputnik V-Vaccinated Individuals. Vaccines; 2022; 10, 817. [DOI: https://dx.doi.org/10.3390/vaccines10050817] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35632574]

38. Wolter, N.; Jassat, W.; Walaza, S.; Welch, R.; Moultrie, H.; Groome, M.; Amoako, D.G.; Everatt, J.; Jinal, N.B.; Scheepers, C. et al. Early assessment of the clinical severity of the SARS-CoV-2 Omicron variant in South Africa. medRxiv; 2021; [DOI: https://dx.doi.org/10.1101/2021.12.21.21268116]

39. Parra-Lucares, A.; Segura, P.; Rojas, V.; Pumarino, C.; Saint-Pierre, G.; Toro, L. Emergence of SARS-CoV-2 Variants in the World: How Could This Happen?. Life; 2022; 12, 194. [DOI: https://dx.doi.org/10.3390/life12020194] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35207482]

40. Ali, A.; ur Rahman, M.; Arfan, M.; Shah, Z.; Kumam, P.; Deebani, W. Investigation of time-fractional SIQR Covid-19 mathematical model with fractal-fractional Mittage-Leffler kernel. AEJ; 2022; 61, pp. 7771-7779. [DOI: https://dx.doi.org/10.1016/j.aej.2022.01.030]

41. Yakusheva, O.; van den Broek-Altenburg, E.; Brekke, G.; Atherly, A. Lives saved and lost in the first six month of the US COVID-19 pandemic: A retrospective cost-benefit analysis. PLoS ONE; 2022; 17, e0261759. [DOI: https://dx.doi.org/10.1371/journal.pone.0261759]