Abstract
Background High rates of virus transmission and the presence of variants of concern can affect vaccine effectiveness (VE). Both conditions occur in low-income countries, which primarily use viral vector or inactivated virus vaccine technologies. However, few VE analyses have been conducted in such countries, and most lack the power to evaluate effectiveness in subgroups, such as the elderly.
Methods The present retrospective cohort study evaluated the effectiveness of Vaxzevria and CoronaVac vaccines for COVID-19-related infection in 60,577,870 Brazilian vaccinees from January 18 to June 30, 2021.
Study outcomes included documented infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), Covid-19–related hospitalization, ICU admission and death. We estimated VE for each outcome as one minus the hazard ratio using Cox regression adjusted for age, sex, Brazilian deprivation index, and month/region of dose administration.
Results Vaccination with Vaxzevria or CoronaVac was found to be effective against SARS-CoV-2 infection and highly effective against hospitalization, ICU admission and death in individuals up to 79 years. From 80-89 years of age, 91.2 (95CI: 89.1-92.9) VE against death was seen in Vaxzevria-vaccinated individuals versus 67.3 (95CI: 63.6-70.6) for Coronvac. Above 90 years, 70.5 (95CI: 51.4-82.1) protection was conferred to Vaxzevria-vaccinated individuals versus 35.4 (95CI: 23.8-45.1) in Coronavac-vaccinated individuals
Conclusions Both vaccines demonstrated overall effectiveness against severe COVID-19 up to 80 years of age. Our results suggest that individuals aged 90 years or older may benefit from an expedited third booster dose. Ongoing evaluations, including any additional vaccines authorized, are crucial to monitoring long-term vaccine effectiveness.
Background
Several COVID-19 vaccines have proved efficacious, and many of them are being extensively used around the world.1–4 While high-income countries preferentially administer mRNA-based vaccines, lower- and middle-income countries have employed vaccines based on viral vectors or inactivated virus technologies. A timely evaluation of the effectiveness of the currently available vaccines across different regions is essential for a comprehensive understanding of vaccine impact, considering significant variations in vaccination schedules, virus transmission and the emergence of viral variants, in addition to social and cultural standards and local health system conditions.
Brazil is one of the countries most affected by the pandemic, with high rates of transmission. The Brazilian COVID-19 vaccination program initially relied on the vaccines Vaxzevria/Fiocruz (previously Oxford-AstraZeneca or ChAdOx-1), approved in 181 countries, and Sinovac’s CoronaVac/Butantan, approved in 39 countries.5 The recommended interdose interval in Brazil for Vaxzevria is 12 weeks versus 2-4 weeks for CoronaVac. The period between doses of Vaxzevria has varied in several countries.7 However, CoronaVac has been applied at distinct intervals,1, 8 making direct comparisons difficult. Additionally, several early publications on vaccine effectiveness (VE) evaluated only the initial dose or were limited to analyzing effectiveness against symptomatic infection9, 10 and hospitalization10, 11, i.e., ICU admission and death were not addressed.
Nationwide evaluations of the effectiveness of COVID-19 vaccines in Brazil offers advantages, as this country’s large population is distributed throughout several regions with considerable differences in socio-economic aspects and access to medical facilities.
Nonetheless, data collection systems are identical throughout the entire country, offering a comprehensive source of data to perform a countrywide VE evaluations. The COVID-19 vaccination campaign was initiated nationwide on January 18, 2021. By June, a large number of vaccinees had received either Vaxzevria/Fiocruz or CoronaVac/Butantan vaccines, allowing for a detailed evaluation of the effectiveness of both vaccines while considering several outcomes and stratified age ranges, making it possible to examine in detail specific age effects previously not investigated.
A significant issue regarding the VE of vaccines against COVID-19 is the degree of circulation of distinct SARS-CoV-2 variants of concern (VOC) in different regions. During the course of the present study, the Gamma variant was the most frequent in all regions of Brazil.12 Importantly, the literature contains few reports on the VE of Vaxzevria and Coronavac against the Gamma variant.1, 10, 13
The present study aimed to evaluate the effectiveness of Vaxzevria and Coronavac vaccines in 60,577,870 Brazilian vaccinee with respect to several different outcomes: COVID-19 related infection, hospitalization, ICU admission and death, between January 18 and June 30, 2021.
Methods
Study design and datasets
We conducted a retrospective cohort using individual-level information on demographic, clinical characteristics, and SARS-COV-2 laboratory tests from the Brazilian federal health registries. The Brazilian Ministry of Health Department of Informatics (DATASUS) provided unidentified datasets of the COVID-19 Vaccination Campaign dataset (SI-PNI), the Acute Respiratory Infection Suspected Cases dataset (e-SUS Notifica), and the National Epidemiological Surveillance System registry for Severe Acute Respiratory Infection/Illness (SIVEP-Gripe). A key-coded individual identification number present in the three datasets was used for a deterministic linkage and then removed from the resulting linked dataset used in our analyses. No personally identifiable data was accessed at any stage. Codebooks, scripts and public dataset version will be available at https://vigivac.fiocruz.br
SI-PNI is a data warehouse run by DATASUS with all the vaccine doses administered by health services in Brazil. From SI-PNI, we extracted information on the COVID-19 vaccine received either Sinovac CoronaVac or Vaxzevria (under the names AstraZeneca/Fiocruz or Covishield/ChAdOx1-S), and the dates of the first and second doses. Overall and age-specific Brazilian population estimates for 2021 corrected the all-cause deaths reported in 2020 overall and age were retrieved from the Brazilian Institute of Geography and Statistics.14 Open version of the SI-PNI dataset is available at opendatasus-SI PNI.
The e-SUS Notifica is a national online health surveillance information system where acute respiratory infections cases and COVID-19 suspected or confirmed cases are registered. and has been used as a data source for epidemiological research.15 Open version of e-SUS Notifica is available at opendatasus-eSUS Notifica.
SIVEP-Gripe is the national system used to register SARI-related hospitalizations and deaths caused by influenza or other respiratory viruses. The system is a registry for new respiratory infections since the H1N1 pandemic in 2009 and widely used as a source for epidemiological studies.16–18 All COVID-19 related SARI hospitalizations and deaths (independent of hospitalization) are registered in the system. Open version of the 2021 SIVEP-Gripe dataset is available at opendatasus-SIVEP
From both SIVEP-Gripe and eSUS-Notifica, we extracted information on the date of symptom onset, RT-PCR, and antigen test results for SARS-CoV-2, and from SIVEP-Gripe, we got data of hospitalization, admission to ICU, and hospitalization outcome (discharge or death).
Study population
We included all individuals who received the COVID-19 vaccine first dose between January 18th, 2021, and June 30th, 2021. The study individuals were followed retrospectively to assess infection, hospitalization, admission to ICU, and death with a laboratory-confirmed diagnosis of SARS-CoV-2 up to June 30th, 2021.
We excluded individuals (i) vaccinated with vaccines besides Vaxzevria or CoronaVac, (ii) with inconsistent vaccine records (i.e., individuals who received the second dose without the first dose, received doses from different vaccines or interval between doses less than 14 days), (iii) with confirmed COVID-19 before the date of vaccine administration, and (iv) with missing data for essential covariates (i.e., sex or age).
Exposure and outcomes
We defined vaccination status for each vaccine based on the time elapsed since the administration of a vaccine dose:
≤13 days after the first dose (the reference period)
≥14 days after the first dose and without the second dose (partially vaccinated)
≥14 days after the second dose (fully vaccinated)
We defined the period up to 13 days after the first dose as the reference period for VE estimation based on results of a Phase III randomized controlled trial8 and three test-negative studies.11, 19, 20 The time-lapsed between the date of the first dose and the development of an effective immune response is used to detect bias in test-negative case-control studies to estimate vaccine effectiveness, the theoretical frame for such use has been discussed by Hitchings et al.21 We also analyzed vaccine effectiveness for 1 to 13 days after the second dose, with the results presented in supplementary table S1).
Laboratory confirmation of COVID-19 with a positive RT-PCR or antigen test result) was required for inclusion in the analyses. The outcomes analyzed were infection, hospitalization, admission to an intensive care unit (ICU), and death by COVID-19. We considered the time between day one of the first or second vaccination up to the symptom s onset for each’ outcome. Individuals whose symptoms started on the same day of the first vaccination dose were assigned one day of follow-up time. Death was considered at any time regardless of prior hospitalization. ICU admission was considered at any time point between the admission and the discharge or death dates.
Statistical analyses
In the primary analysis, we used a Cox regression model to estimate the hazard ratio (HR) of COVID-19 infection, hospitalization, ICU admission, and death for partially and fully vaccinated individuals. The model was adjusted for age, sex, region of residence, socioeconomic status, and month of the 1st dose. We used the Brazilian Deprivation Index (Índice Brasileiro de Privação-IBP), a municipality-level measure of material deprivation, as an indicator of socioeconomic status.21 We estimated vaccine effectiveness (VE) as 1-HR, obtained from a model including all covariates, and reported as a percentage. We also reported the crude VE for each outcome. In addition, we performed a stratified analysis by age groups (<60, 60–69, 70–79, 80–89, ≥90 years) to investigate whether VE was modified by age.
To assess the robustness of our findings, we repeated the principal analysis defining as the reference period the time elapsed up to 10 days after the date of the first dose, as it is expected that the vaccines protection increases with time. Additionally, we examined the VE’ for hospitalization, ICU admission and death using clinical suspected cases besides laboratory confirmed ones.
Analyses were performed using the R statistical software (R Core Team) and its H2O package.23, 24 Descriptive statistics were presented as frequencies and percentages. We used the 95% confidence intervals (CI) of the estimated measures of association for interpreting the findings.
RESULTS
From January 18 to June 30, 2021, 61,783,842 individuals received at least one dose of one of the two COVID-19 vaccines analyzed in this study, and 60,577,870 (98.1%) met the eligibility criteria and were included in the analysis (Figure 1). The majority (63.8%, n=38,664,633 individuals) received at least one dose of Vaxzevria and the remaining (36.2%, n=21,933,237 individuals) received at least one dose of CoronaVac. The majority of our cohort comprised women (56.1%) and individuals aged 60 years or older (44.4%). Compared to individuals that received CoronaVac, individuals that received Vaxzevria were younger (29.3% vs. 70.9% of individuals aged 60 years or older), and a lower proportion had completed the full vaccine schedule (10.6% vs. 82.7%). Vaccination with CoronaVac occurred mainly from January to April 2021, while Vaxzevria was administered predominantly after March 2021 (Figure 2). Among those who received the second dose, the median time between the first and second doses was 85 days (IQR 83–90) for Vaxzevria and 27 days (IQR 21–28) for CoronaVac. Individuals who received at least one dose of Vaxzevria or CoronaVac were mostly women (54.6% vs. 58.7% respectively) and from the southeast region of the country (44.1% vs. 46.3%, respectively) (Table 1).
Table 2 shows the COVID-19 VE analysis results, including number of events and incidence rate per 1000 person-days and supplementary table S1 shows the crude and adjusted VE analysis . We observed that individuals with full vaccination schedule (i.e., ≥ 14 days after the second dose) with Vaxzevria had a 70.0% (95% CI 68.6 to 71.3) lower risk of infection, 86.8% (95% CI 85.2 to 88.2) lower risk of hospitalization, 88.1% (95% CI 85.4 to 90.3) lower risk of ICU admission, and 90.2% (95% CI 88.3 to 91.8) lower risk of death. Partial vaccination (i.e., ≥14 days after the first dose up to the second dose) with Vaxzeria was associated with a 32.7% lower risk of infection (95% CI 31.9 to 33.5) and at least 50% lower risk of hospitalization (51.7%; 95% CI 50.4 to 52.9), ICU admission (53.6%; 95% CI 51.4 to 55.6), and death (49.3%; 95% CI 47.0 to 51.5). Complete vaccination with CoronaVac was associated with a 54.2 (95% CI 53.4-55.0) lower risk of infection, 72.6% (95% CI 71.6 to 73.6) lower risk of hospitalization, 74.2% (95% CI 72.6 to 75.7) lower risk of ICU admission, and 74.0% (95% CI 72.6 to 75.3) lower risk of death. Partial vaccination with CoronaVac was associated with less than 50% of reduction in the risk of infection (16.2%; 95% CI 15.1 to 17.4), hospitalization (26.5%; 95% CI 24.6 to 28.4), ICU admission (28.1%; 95% CI 24.9 to 31.1), and death (29.4%; 95% CI 26.7 to 32.0).
When stratifying the analysis by age, complete vaccination with Vaxzevria or CoronaVac presented a similar VE within all age groups, with the exception among individuals aged 90 years or older (Table S2, Figure 3).
In the analysis using the reference period of up to 10 days after the first dose, we found VE point and interval estimates similar to those found in the primary analysis for both Vaxzeria and Coronavac vaccines (Table S3). The results using all clinical suspected and laboratory confirmed cases for the outcomes of hospitalization, ICU admission and death were qualitatively equal to those found in primary analysis (Table S4).
DISCUSSION
Here we present nationwide results on the effectiveness of vaccination with CoronaVac/Butantan and Vaxzevria/Fiocruz after the first six months of the vaccination campaign in Brazil. Analyzing data from almost 61 million individuals vaccinated with at least one dose, our results demonstrate strong evidence of 70.0% and 54.2% protection against infection after full vaccination with Vaxzevria and CoronaVac, respectively.
Vaxzevria offered approximately 90% effectiveness against hospitalization, ICU admission and death, while CoronaVac provided approximately 75% protection following full vaccination.
Our findings regarding the Coronovac/Butantan vaccine are compatible with a previous Brazilian efficacy study24, but lower than the 83.5% protection reported by a Turkish efficacy trial.8 The effectiveness determined by a cohort study in Chile was also higher than our findings for infection (66.5% vs. 54.2%) as well as hospitalization (87.5% vs. 72.6%).
Differences between the study in Chile and the present analyses of Brazilian vaccinees may be partially explained by the higher frequency of younger individuals in the Chile study (51.2% vs. 29.1% of individuals younger than 60 years old). During the vaccination campaign, Brazil experienced health system collapse in several states, which may have influenced death rates, especially between February and May, likely affecting CoronaVac estimates more markedly due to its greater availability of this vaccine in the early stages of the vaccination program. Another reason for these differences could be the increased circulation of the Gamma lineage detected in these countries, which has been estimated at 28.6% in Chile and 69.6% in Brazil during both study periods.1, 12 In plasma samples obtained from individuals fully vaccinated with CoronaVac, a reduced capacity to neutralize the Gamma variant was observed.1 Furthermore, 9.9% of the Brazilian population was fully vaccinated from January to May 2021, compared to almost 35.4% of Chile’s population. This may have contributed to lower viral transmission in Chile compared to Brazil.1
For Vaxzevria, our findings of 70.0% effectiveness against infection exceeded the levels of 66.7% effectiveness reported in a combined analysis of four clinical trials conducted in the UK, South Africa, and Brazil.7 Effectiveness against hospitalization was consistent with the 80% and 89% protection observed in studies in Scotland3 and England,11 respectively.
Additionally, our findings support the high level of protection offered by Vaxzevria despite the abundant circulation of the Gamma variant in Brazil during the period studied. Few studies have reported on the VE of Vaxzevria in populations infected by VOCs.1, 9, 10, 13, 20 Studies analyzing effectiveness against VOCs have mainly focused on protection against symptomatic infection or hospitalization.9, 10, 13 Taken together, the findings reported herein combined with data in the literature confirm a consistently high rate of protection against moderate to severe COVID-19 in real-world studies, despite abundant circulation of VOCs.
Protection was shown to vary according to age group. The VE of CoronaVac/Butantan was close to 80% against death in individuals aged up to 79 years of age. However, a reduction in effectiveness was observed after 80 years of age, with only 35.4% protection against death seen in individuals over 90. In contrast, the Vaxzevria/Fiocruz vaccine achieved close to 90% protection against death in individuals aged less than 90 years, while a VE of 70.5% was found in those older than 90 years of age. It is reasonable to attribute the observed reduction in effectiveness to immunosenescence, which is commonly associated with a higher frequency of comorbidities, and may imply higher death rates. In the context of limited vaccine availability, the precise identification of age limits at which point immune protection becomes impaired can provide valuable evidence to inform public health measures.
Considering the current scenario in Brazil, our findings demonstrate the eventual need for a vaccine booster dose in individuals aged 80 years or older who received CoronaVac, as well as for individuals over 90 years immunized with Vaxzevria.
The differences evidenced in effectiveness between Vaxzevria and CoronaVac may be related to the distinct technologies used be each of these two products, as well as how they influence immunogenicity. Both vaccines analyzed herein activate immunological mechanisms and trigger a neutralizing antibody response against viral particles. However, CoronaVac, a whole-cell inactivated vaccine, elicits a less potent cellular response than Vaxzevria, an adenoviral-vectored vaccine.25 Additionally, Vaxzevria was shown to induce a higher peak neutralizing antibody response than CoronaVac.27 Thus, the intrinsic characteristics of each formulation may serve to explain differences observed in both clinical trials and vaccine effectivity studies.1, 26, 28
A relevant strength of our study is its large sample size, due to the use of the complete dataset covering the Brazilian COVID-19 vaccination campaign from January to June 2021. This large sample allowed us to identify the age limits in which immune protection becomes impaired, especially with regard to CoronaVac. Sensitivity analyses further confirmed the robustness of our findings. However, our study is also subject to some limitations. First, as VE was estimated using observational data, our analysis is subject to data availability and, therefore, to potential confounders. Although our analyses were not controled for comorbidities, crude and adjusted VE estimates were similar. In addition, comorbidities have been identified as the causal pathway between age and COVID-19 severity. Therefore, by controlling for age, we are also indirectly controlling for comorbidities.29 Second, in contrast to many VE studies, the reference period used herein for comparison purposes was 1-13 days after vaccination. Although using early post-vaccination as a reference may underestimate VE, previous studies have used a similar approach and obtained VE results similar to those found in clinical trials.30, 31 The early post-vaccination period can also be used as a bias indicator related to differences in SARS-CoV-2 infection risk. Additionally, the effectiveness results of the present report are similar, in the pertinent age ranges, to reports on both vaccines using distinct approaches.1, 19, 20 Finally, we also performed sensitivity analysis, which demonstrated similar results when a 0-10 day reference period was applied.
Using the data available in Brazil, we estimated overall VE for each vaccine evaluated as well as by age group. Vaxzevria/Fiocruz and CoronaVac/Butantan were both shown to be highly protective against severe COVID-19 in the population aged up to 80 years, yet due to decreased VE an early booster dose may be considered for those over 80 years of age who received CoronaVac, and especially for individuals aged over 90 years regardless of which of these two vaccines were administered. Despite high population adherence, the vaccination campaign is evolving unevenly throughout Brazil, and continuous monitoring of VE in the current context may provide sound evidence to inform public health measures.
ETHICAL CONSIDERATIONS
The Brazilian National Commission in Research Ethics approved the research protocol (CONEP approval number 4.921.308). The study was considered exempt from informed consent; no human health risks were identified. All work presented here used unidentified secondary data in accordance with the Brazilian Personal Data Protection General Law (LGPD). Data was manipulated in a secure computing environment, ensuring protection against data leakage and records reidentification.
DECLARATION OF INTERESTS
VO, VB, MB, and MB-N are employees from Fiocruz, a federal public institution, which manufactures Vaxzevria in Brazil, through a full technology transfer agreement with AstraZeneca. Fiocruz allocates all its manufactured products to the Ministry of Health for the public health service (SUS) use. All other authors report no potential competing interest.
FUNDING
This study was partially supported by a donation from the “Fazer o bem faz bem” program.
DATA SHARING
We used third-party data, provided by the Brazilian Ministry of Health. Any request for access to the data shall be directed to DATASUS - Ministry of Health Brazil: https://datasus.saude.gov.br/
Data Availability
We used third-party data, provided by the Brazilian Ministry of Health. Any request for access to the data shall be directed to DATASUS - Ministry of Health Brazil: https://datasus.saude.gov.br/
ACKNOWLEDGMENTS
The authors thank DATASUS for its excellent work in providing unidentified databases. GW, MB, and MB-N are research fellows from CNPq, the Brazilian National Research Council. The authors would also like to thank Andris K. Walter for English language revision and manuscript copyediting assistance.