COVID-19, caused by SARS-CoV-2 infection, left widespread impacts worldwide. In Brazil, immunization reduced incidence rates. However, six months later, waning neutralizing antibody titers and new immune-evading variants increased cases, resulting in recurring waves. This study evaluated hospitalized COVID-19 patients after the vaccination rollout, comparing the clinical outcomes between vaccinated and unvaccinated patients. Positive samples underwent nucleotide sequencing. A total of 218 patients were included; 202 (92 %) had vaccination data, 98 received at least one dose, and 64 completed the vaccination schedule, predominantly with CoronaVac®. Vaccinated individuals were older on average since the campaign was primarily conducted among the elderly. The Gamma variant predominated during the study period. While not statistically significant, trends indicated greater respiratory assistance needs, more extended hospital stays, and increased ICU time among unvaccinated patients. Mortality was 45 % in vaccinated and 37 % in unvaccinated groups, with no notable difference. However, patients with a complete vaccination schedule showed a higher chance of survival, though not significant (p = 0.11). The factors significantly associated with higher mortality were older patients, those requiring vasopressor drugs, and mechanical ventilation. These findings provide clinical, epidemiological, and phylogenetic insights into COVID-19 patients during vaccination implementation. They underscore the need to evaluate vaccine effectiveness against circulating variants and highlight the importance of complete vaccination schedules for improving patient outcomes.
In 2019, following its initial identification in Wuhan, China, SARS-CoV-2 rapidly evolved into a pandemic pathogen, spreading across multiple continents and triggering a profound global impact.1
In Brazil, the vaccination campaign began on a large scale, with four vaccines developed using different platforms. The nationwide SARS-CoV-2 immunization program was launched in January 2021 and was gradually rolled out, initially prioritizing older adults, healthcare professionals, and individuals at high risk of severe disease. Subsequently, it expanded to include the entire adult population. Even though antivirals are available in the public health system for outpatients and hospitalized patients who meet clinical criteria, vaccination remains the most effective prevention strategy.2,3
In Brazil, the first vaccine utilized to control the pandemic was CoronaVac® (Sinovac Life Sciences, Beijing, China), based on the inactivated viruses’ platform, followed by Comirnaty® (BNT162b2, RNAm, Pfizer-BioNTech, USA), ChAdOx1 nCoV-19 (AZD1222, Oxford/AstraZeneca, UK), and Ad26.COV2.S (Johnson &. Johnson–Janssen). As a result of the slow and gradual implementation of the immunization program and the surge of viral variants with mutations and deletions in the S protein sites that enhance viral entry, replication, and immune evasion, the vaccines have been less effective in preventing hospitalization.4
This scenario, combined with the decline in antibody titer observed six months after vaccination, raises concerns about the long-term protection of the immunized population, posing a significant public health challenge. Currently, the vaccination schedule includes multiple boosters for all ages, totaling >500 million doses administered (https://infoms.saude.gov.br/extensions/SEIDIGI_DEMAS_Vacina_C19/SEIDIGI_DEMAS_Vacina_C19.html). To this day, little is known about the long-term effectiveness of COVID-19 vaccines in reducing the severity and mortality associated with the currently circulating variants. This uncertainty persists despite the availability of antiviral therapies and multiple reinfections among patients. This study aimed to describe the epidemiological, clinical profiles, and outcomes of a hospitalized cohort of SARS-CoV-2 infected patients and assess whether vaccination was a protective factor against severe disease during the emergence of Gamma and Delta variants of concern.
Material and methodsStudy designThis retrospective cohort study was conducted with hospitalized patients diagnosed with COVID-19 at the Clinical Hospital Complex of the Federal University of Paraná (CHC-UFPR), a tertiary academic hospital in southern Brazil. Patients had been previously screened and included in the multicentric study “The effectiveness of COVID-19 vaccines in Latin America, 2021: a multicenter regional case-control study”, supported by the Pan American Health Organization (PAHO) and Oswaldo Cruz Institute (Fiocruz). The Institutional Ethics Board revised and approved this study (n° 53,891,121.0.0000.0096).
The inclusion criteria were hospitalized patients between August 2021 and February 2022 who were ≥ 18-years-old, eligible for the anti-SARS-CoV-2 vaccine, presenting cough and fever, with at least 3-days of symptoms, and had not been infected with SARS-CoV-2 in the past 90-days. Patients with no clinical data, contraindication to the vaccine, or unconfirmed SARS-CoV-2 infection were excluded.
Data and samplesThe researchers obtained retrospective clinical, epidemiological, vaccination status, and outcome data from electronic medical charts. Data collection encompassed the entirety of the patient’s hospitalization, concluding at the time of discharge or death.
Patients were categorized by vaccination status: vaccinated and non-vaccinated, complete and incomplete vaccination schedules, and timely and untimely vaccination for those fully vaccinated. Timely vaccination was defined as receiving the vaccination at least 14-days before the infection. A compliance group was used to identify those patients who had received full and timely vaccine doses.
After diagnosis, swab samples with sufficient volume and a SARS-CoV-2 qRT-PCR with Ct < 30 were sent to the Respiratory Viruses and Measles Laboratory (LVRS, Fiocruz, Rio de Janeiro, Brazil) for complete genome sequencing. The Fiocruz Genomic Network performed the analysis using the ViralFlow pipeline (Available at: https://viralflow.github.io/). Data on circulating VOCs were obtained from the Genomic Surveillance of SARS-CoV-2 in Brazil. (SESA-PR, 2021, available at: https://www.saude.pr.gov.br/Pagina/Coronavirus-COVID-19).
Phylogenetic analysisTo examine the evolution and prevalence of SARS-CoV-2 strains in admitted patients, the complete genome sequences generated in this study were compared with genome sequences from reference strains identified in Brazil. These sequences were retrieved from GISAID (https://gisaid.org/) until December 2022. All sequences were aligned using the MAFFT program (Available at: https://mafft.cbrc.jp/) and manually edited to remove artifacts such as insertions, deletions, and gaps using the Aliview program (Available at: https://github.com/NCIP/alview). The Maximum Likelihood (ML) phylogenetic trees were estimated using IQ-TREE2 (Available at: http://www.iqtree.org/), using the GTR nucleotide substitution model, and the resulting tree was edited using FigTree (V.1.4.3, available at http://tree.bio.ed.ac.uk/software/figtree/).
Data analysisData were analyzed using R Core Team Statistical Software, V.4.2.1 (R Foundation for Statistics Computing, Vienna, Austria, 2014). Relative and absolute frequencies were analyzed, and Fisher's exact, Chi-Square, Mann-Whitney, and Gray tests were applied, as appropriate. An unadjusted bivariate analysis was performed comparing fully vaccinated patients with unvaccinated or partially vaccinated patients and comparing death and discharge outcomes. Multivariate analysis to assess disease severity was performed, considering that those with bivariate analysis results had presented a p-value < 0.2 as explanatory variables. All p-values are based on two-tailed comparisons, and the significance level was set at p < 0.05.
The Kaplan-Meier method compared the case survival curves according to the vaccination schedule. The log-rank test was used to verify differences between groups. Survival curves were stratified by the presence of vaccination and by the completeness of the vaccination schedule.
ResultsA total of 218 participants were initially included. Of these, 202 (92 %) had vaccination data available, with 98 vaccinated and 107 unvaccinated patients (Fig. 1). Among the vaccinated patients, only 6 (3 %) did not have data on the vaccination scheme used. In addition, some medical records reported a complete vaccination schedule before admission, but the vaccine date was not informed, hindering the ability to evaluate the timeliness of immunization.
Table 1 compares the clinical and demographic profiles of vaccinated (45 %) and non-vaccinated (49 %) patients. For both groups, most patients were male, had completed primary school, and never smoked. The most common symptoms were dyspnea (70 %, 69 %) and cough (58 %, 51 %) for vaccinated and unvaccinated patients, respectively. Regarding the vaccinated ones, more than half (60 %) received the CoronaVac® vaccine (Sinovac Biotech, Bejin, China). The second vaccine most frequently administered beyond the analyzed sample was the ChAdOx1 nCoV-19 vaccine (27), followed by Comirnaty® (9) and Ad26. COV2.S.1 Gamma and Delta VOCs were the most frequent viral genotypes found. Unvaccinated patients were significantly younger. There was a trend towards a higher need for respiratory support, longer length of hospitalization, and increased ICU time among the unvaccinated group, although the difference was not statistically significant. There was no significant difference in the mortality rate between the vaccinated and unvaccinated groups, with rates of 45 % and 37 %, respectively.
Clinical, epidemiological, and outcome characteristics of the vaccinated and unvaccinated groups, and among the vaccinated, according to compliance with the vaccine schedule.
SD, Standard Deviation; IQR, Interquartile Range. In bold: statistically significant p-values.
Out of 98 vaccinated patients, 92 confirmed the completion of immunization. Of these, 70 % (65/92) adhered to the recommended vaccination schedule (compliance), which included one dose of Ad26.COV2.S and two doses of another vaccine (Table 1). Dyspnea was the most frequent symptom, reported by90 % of the non-compliance group, which also had a significantly higher mortality rate compared to the compliance group (68 % vs. 36 %, p = 0.0063).
A survival curve was built to examine the impact of immunization between the different groups. When analyzing the vaccinated and unvaccinated groups, the chance of survival was significantly higher among the unvaccinated ones (Fig. 2A) (p = 0.035), and patients with a complete vaccination schedule had a greater chance of survival (Fig. 2B), but the difference was not significant (p = 0.11).
Subsequently, a comparison was made between patients according to clinical outcomes, discharge or death. The results are presented in Table 2. In the adjusted analysis, the factors associated with the outcome of death were older age, the need for vasopressor drugs, and mechanical ventilation.
Correlation between clinical and epidemiological findings and outcomes.
IQR, Interquartile Range; RR, Respiratory Rate; ICU, Intensive Care Unit; NS, Non-Significant result.
aOnly values with statistically significant results are presented.
All 218 admissions analyzed in this study were ordered by month. Fig. 3 compares these admissions to the distribution of Variants of Concern (VOCs) identified at the time of the survey in Paraná.
The sequencing results were compiled and are shown as a phylogenetic tree (Fig. 4). The Gama variant was the predominant VOC among Delta, Alpha, and Zeta variants during the study period.
Phylogenetic tree of SARS-CoV-2 strains from admitted patients. Note: Phylogenetic trees based on maximum-likelihood analysis were built using complete sequences of SARS-CoV-2 strains from patients hospitalized with COVID-19 between August 2021 and February 2022. The black dots at the end of the branches represent taxa. The study’s sequences are numbered, and reference sequences are identified in the figure and are available in Supplementary Material.
The rapid development of vaccines using various platforms has been the primary strategy in fighting SARS-CoV-2. However, low- and middle-developed countries faced restricted availability of vaccines with higher efficacy.5 Our study examined the outcomes of vaccinated patients compared with non-vaccinated patients hospitalized with COVID-19 a few months after starting anti-COVID-19 vaccination and showed no significant differences between the groups. Factors including adherence to vaccination schedules and the emergence of Variant of Concern (VOC) strains may have influenced these results.
Brazil faced increasing challenges in vaccine adherence, stemming from concerns about vaccine efficacy and difficulties in identifying SARS-CoV-2 infections, particularly in older adults. These issues heightened underreporting and vaccine hesitancy, as some segments of the population were hesitant to recognize the severity of the infection and the advantages of vaccination. Consequently, these factors influenced hospital admissions and vaccination data trends.
Studies have shown that individuals vaccinated for >6-months have a critical reduction in antibody titers.6 However, demographic and individual factors can affect the immune response to the immunogenic stimulus. This, combined with extensive viral genetic variability, results in the frequent emergence of viral variants. Therefore, ongoing studies must assess the relationship between vaccination, SARS-CoV-2 infection, and disease severity.7
When comparing the vaccinated and unvaccinated groups, only the age variable remained statistically different after the adjusted analysis. Specifically, older individuals had a greater likelihood of being vaccinated. This finding reflects the Brazilian vaccination scenario, where the vaccination campaign began in January 2021 and followed a list of priority groups for vaccination. Older adults, due to the immunosenescence and a higher prevalence of comorbidities, were given priority for the anti-SARS-CoV-2 vaccine, along with healthcare professionals, who were at greater risk of infection due to more frequent exposure.8,9
Regarding sociodemographic data, some reports indicate that it is essential to identify risk or protective factors, such as male sex, dark or brown skin, obesity, and other comorbidities. However, we could not observe the influence of these factors on the patient's outcomes.10
The epidemiological profile of patients with a compliance and noncompliance vaccination schedule was similar to the vaccinated and unvaccinated patients. Age continued to be a statistically significant variable, which is in line with the Brazilian scenario at the time of data collection. Patients with incomplete vaccination showed slightly more severe symptoms, including dyspnea, fatigue, muscle weakness, and altered mental status. However, desaturation and the need for supplemental oxygen were common in both groups. Notably, CHC-UFPR, as a tertiary reference hospital, frequently manages moderate to severe cases.
Among patients with a complete vaccination schedule, discharge was the most frequent outcome, while among patients with an incomplete vaccination schedule, the most frequent outcome was death. For these groups, age, type of vaccine, VOCs identified, and outcome were the most significant variables in the bivariate analysis. However, only age remained statistically significant after adjustment.
The Gamma VOC was the most prevalent compared to Zeta, Delta, and Alpha. Reports from respiratory virus surveillance in Curitiba showed alignment between the phylogenetic tree distribution and the circulating variants identified in the city. While Gamma is often associated with greater disease severity, consistent with our findings and a previous analysis conducted by our group, we could not demonstrate an increase in disease severity linked to Gamma or other VOCs.11
The survival curve did not show a statistically significant difference between vaccinated and non-vaccinated patients. In the Brazilian context, CoronaVac was the most utilized vaccine in the study population. A subsequent multicenter vaccine effectiveness study conducted in Latin America, which included data from the Curitiba Center (Brazil), reported that all vaccines administered in Brazilian territory showed effectiveness. However, there was variability between products.9 As expected, vaccine effectiveness was greater in patients with a complete vaccination schedule. CoronaVac had an efficacy rate of 53 % in this study in patients with a complete vaccination schedule. Furthermore, the vaccine's effectiveness decreased with advancing age.12 Given that the majority of participants in this study are mostly older individuals, it is possible that the vaccine effectiveness was reduced due to a decline in immunity, reinforcing the need to administer booster doses.6,13–15
Booster doses are recommended in homologous or heterologous modalities and have positive effects in increasing immunity, including against VOCs, which have emerged after the development of vaccines that targeted the infection of the ancestral strain (Wuhan).16 Furthermore, a single dose of CoronaVac was associated with low protection in symptomatic patients and hospital admissions, highlighting the need for a complete vaccination schedule.17
To understand the higher probability of survival among unvaccinated patients, an analysis of fatal cases revealed that vaccinated patients who had not developed serological immunity before becoming infected with the COVID-19 virus were more likely to die (Suppl. data ‒ Fig. 4). Another study at CHC-UFPR analyzing the serum of healthcare professionals at the hospital to identify anti-S1 and anti-N IgG antibodies as a tool to predict neutralization found a seroconversion rate of 97 % among HCWs 40-days after receiving the first dose.6 In contrast to other vaccines in Brazil, CoronaVac showed low protection until the second dose was administered (over four weeks after the first dose) .17 Ranzini and colleagues recommended prioritizing full vaccination schedules over single doses for better protection, especially for older adults.
Additionally, Lira et al. reported similar findings, showing more severe infections in individuals who had not completed the vaccination schedule or were not vaccinated at all. However, unvaccinated patients still showed a 26 % lower chance of dying. These results could be influenced by factors such as underreporting during the pandemic, the slow progress of vaccination in the country, and vaccine hesitancy among the population, all of which may have contributed to these observations as confounding factors.18 Therefore, administering the full vaccination schedule to a single portion of the population would have a positive effect compared to administering only one dose to a more significant portion of the population in the context of the CoronaVac immunizer.15
The effectiveness of CoronaVac among healthcare professionals showed 50.7 %, 83.7 %, and 100 % efficacy in symptomatic infections, hospitalizations, and severe cases, respectively.16 Other studies showed 47 % effectiveness in preventing symptomatic infections and severity in the population ≥ 70-years-old (61 % effectiveness in preventing cases of death and 55 % for hospitalizations).13 A study examining the Latin American population showed 33 % effectiveness among partially vaccinated individuals, while among those fully vaccinated, a 53 % effectiveness was found.9
This study had some limitations due to its retrospective nature and the critical period of the pandemic, which may have resulted in the potential loss of clinical data. However, it contributes by providing clinical, epidemiological, and even phylogenetic information from COVID-19 patients hospitalized during the pandemic. Furthermore, this analysis highlights the critical importance of ongoing assessments regarding the impact of decisions on implementing preventive measures on public health, which can help with future vaccination strategies.
Since the new coronavirus was described in early 2020, there has been a major effort to develop vaccines that enhance immunity against the new agent. However, after the COVID-19 vaccines were developed, the circulating lineage corresponded to the ancestral variant, and only approximately a year later, the vaccines started to be administered to the population. Additional variants were described, which had a crucial role in sustaining the pandemic through point mutations that conferred immune response escape and high rates of infection and transmissibility.16,19
In conclusion, this study found no evidence of vaccination affecting the survival of hospitalized patients in this cohort. However, the results are specific to a tertiary care hospital, where patients generally present with more severe cases, multiple comorbidities, and higher intensive care requirements, limiting the generalizability to less severe cases. Therefore, it is crucial to highlight that these findings should not be extrapolated, especially given the clear and undeniable benefits vaccines have shown in controlling the pandemic. Continuous surveillance of circulating strains, such as through the sentinel surveillance network for influenza and other respiratory virus infections, and ongoing research on vaccinated individuals remain crucial for generating data to inform public health policies.
CRediT authorship contribution statementLaura Holtman-Ferreira: Formal analysis, Investigation, Writing – original draft. Elessandra de Souza Bitencourt: Funding acquisition, Project administration. Betina Mendez de Alcantara Gabardo: Data curation, Formal analysis, Investigation, Methodology. Susanne Edinger Pereira: Data curation, Formal analysis, Investigation, Methodology. Francine Teixeira: Data curation, Formal analysis, Investigation, Methodology. Diego da Silva Magatão: Data curation, Formal analysis, Investigation, Methodology. Vitor Loureiro Dias: Data curation, Formal analysis, Investigation, Methodology. Ricardo Petterle: Data curation. Meri Bordignon Nogueira: Funding acquisition, Investigation. Sonia Mara Raboni: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – review & editing.