Progressive Increase in Virulence of Novel SARS-CoV-2 Variants in Ontario, Canada

Introduction

Novel SARS-CoV-2 variants of concern (VOC), including viral lineages carrying the N501Y and/or E484K mutations (Alpha/B.1.1.7, Beta/B.1.351 and Gamma /P.1, were first identified in Ontario, Canada in December 2020 (1). While infection with such lineages was initially uncommon, these VOCs outcompeted earlier SARS-CoV-2 lineages in Ontario, and as of late April 2021 comprised almost all new infections in Ontario, with Alpha being the most prevalent lineage(1). In April 2021, the B.1.617.2 variant, now known as “Delta” under revised World Health Organization nomenclature, emerged in the province and outcompeted earlier VOCs, representing the majority of infections in the province as of July 2021 (2, 3).

This serial replacement by emerging variants reflects progressively higher effective reproduction numbers which allow novel variants to outcompete previously dominant strains in the face of identical disease control interventions (4-6). However, VOCs are also concerning because of an apparent increase in virulence, with increased risk of hospitalization, intensive care unit admission and death, after adjustment for age and other predictive factors (7-10). Although the increased virulence of strains with the N501Y mutation relative to strains which lack this mutation has been described (7-9), only limited information is available on the virulence of infection with the Delta variant, relative to earlier N501Y-positive VOCs (i.e., Alpha, Beta, and Gamma) (10). Our objectives were (i) to use Ontario’s COVID-19 case data to evaluate the virulence of N501Y-positive variants relative to earlier SARS-CoV-2 lineages; and (ii) to evaluate the virulence of the Delta variant of SARS-CoV-2 relative to N501Y-positive VOCs.

Methods

Data Sources

We created a retrospective cohort of people in Ontario testing positive for SARS-CoV-2 and screened for variants of concern, with dates of test report between February 7 and June 22, 2021. Case information was extracted from the Ontario provincial Case and Contact Management (CCM) database as described elsewhere (11). Although VOC cases were identified in Ontario beginning in December 2020, systematic screening was not implemented until February 2021, after which time all COVID-19 PCR positive specimens with a cycle threshold (Ct) value ≤ 35 were screened for the N501Y mutation using a single nucleotide polymorphism real-time PCR assay developed at the Public Health Ontario Laboratory. On March 22, 2021, Ontario initiated universal screening for N501Y and E484K mutations using a multiplex real-time PCR assay on all specimens testing positive for SARS-CoV-2 with a Ct value ≤ 35. Whole genome sequencing (WGS) was performed on a 5% sample of specimens regardless of the presence of mutations. Initially all specimens with the N501Y and/or E484K mutation and Ct value < 30 were sequenced; however, as of June 2021 routine WGS was no longer performed on specimens without E484K, under the presumption that such specimens were of the Alpha lineage. By late April 2021, more than 90% of infections in Ontario were screen positive for N501Y(Figure 1); subsequently, N501Y mutations became less common, and N501Y-negative specimens subjected to WGS were demonstrated to be predominantly Delta variant, with Delta representing > 60% of new infections in Ontario on July 1, 2021 (2, 3). Analyses were restricted to individuals whose viral isolate had been screened for a VOC. Cases were classified as N501Y-positive (screening positive for N501Y or identified as Alpha, Beta, or Gamma by WGS), probable Delta variant (identified by WGS at any point in time or negative for N501Y and any other mutations from May 1, 2021 onward), or not VOC (all N501Y negative cases between February 7 and April 30,2021 and all non-Delta cases testing positive for a mutation other than N501Y regardless of date). Individuals for whom VOC screening information was not available (N = 31,464), or for whom screening could not be completed (N = 13,579) were excluded from the analysis. We further restricted our analysis to individuals without a record of long-term care residence, as the epidemiology and severity of SARS-CoV-2 infection in long-term care home residents have been distinct in Ontario (12, 13). The dataset was extracted on July 13, 2021, with the most recent test report date of July 12, 2021. We included cases with a test report date at least 14 days prior (June 27, 2021) to account for delays between testing and occurrence of hospitalization, intensive care unit (ICU) admission, and death. Vaccination data for cases, including dates of administration of first and second doses of approved SARS-CoV-2 vaccines (where relevant) were obtained from COVaxON, a centralized COVID-19 vaccine information system for the province. Individuals were categorized as fully vaccinated if the case episode date occurred 7 or more days after receipt of the second dose; partially vaccinated if the case episode date occurred 14 days or more after receipt of the first dose but less than 7 days after receipt of the second dose; and unvaccinated otherwise.

• Download figure
• Open in new tab

Figure 1. Trends in SARS-CoV-2 Case Occurrence and distribution of hospitalizations, ICU admission, and deaths by VOC status in Ontario, Canada.

(A) Cases of reported SARS-CoV-2 infection by test report date, presented as 7-day average for cases reported between February 7 and June 27, 2021. (B) Distribution of COVID-19 hospitalized cases, cases admitted to ICU, and deaths by VOC status, by week of test report. Cases are coloured by assigned VOC status. Prior to May 1, 2021, Delta cases were only detected by whole genome sequencing. After May 1, 2021, all screened specimens not identified as an N501Y-positive VOC or another variant were classified as probable Delta VOC infections.

Results

We included 212,332 cases in our analysis, with test report dates between February 7 and June 27, 2021. Changes over time in incidence of infections by VOC status, and in incidence of infections resulting in ICU admission, are presented in Figure 1. Among all reported cases, 77% were infections with N501Y-positive VOC; 2.7% were classified as probable Delta infections. In crude comparisons, there were significant differences between variants in probability of hospitalization and ICU; individuals with N501Y-positive variants and probable Delta infection were significantly younger and less likely to have comorbidities than individuals with non-VOC infections; N501Y-positive infections were significantly more common in the Greater Toronto and Hamilton area (GTHA, the province’s largest metropolitan area); Delta variant infections were significantly less common in Ottawa (the province’s second largest metropolitan area) (Table 1).

After adjustment for age, sex, vaccination status, comorbidities, health unit, and temporal trend, large and significant increases in the risk of hospitalization, ICU admission, and death were seen with both N501Y-positive VOC, and probable Delta infections, relative to non-VOC (Table 2). Adjusted risks of hospitalization, ICU admission and death were 52% (43-62%), 89% (67-116%) and 51% (30-74%) higher with N501Y-positive VOC than with non-VOC; respective differences between non-VOC and Delta infections were 108% (80-138%), 234% (164-331%), and 132% (47-230%). We identified a linear increase in severity with progression from non-VOC infections to N501Y-positive infections, to Delta variant infections, with relative increases in risk of hospitalization, ICU admission and death of 49% (41-58%), 86% (68-109%) and 51% (32-72%), respectively (Table 3).

Discussion

The emergence of novel SARS-CoV-2 variants of concern has slowed progress against the pandemic in three distinct ways: (i) by increasing transmissibility and the disease’s reproduction number; (ii) by increasing immune escape and diminishing vaccine effectiveness; and (iii) by increasing the virulence of SARS-CoV-2 infection. Credible reports of enhanced virulence of VOC with the N501Y mutation have emerged over the past several months (7-9), but to our knowledge this is one of the first analyses to demonstrate enhanced virulence of the Delta variant (10).

We show that in the Canadian province of Ontario, VOCs with the N501Y mutation were associated with a markedly increased risk of hospitalization, intensive care unit admission and death among infected individuals, but that the Delta variant, which is now supplanting these other VOCs in Ontario, has increased these risks even further. Individuals infected with VOCs were, on average, younger and less likely to have comorbid conditions than those infected with non-VOC, but nonetheless had higher crude risks of hospitalization and ICU admission. Once we adjusted for confounding factors such as age, vaccination status, comorbidity, and temporal trends, elevated per-infection risk, including risk of death, remained markedly higher with VOCs, and with Delta in particular. Indeed, given the relatively small number of Delta infections in our time series, it is remarkable that a clear and significant elevation of risk of even uncommon, delayed outcomes, such as death, is readily visible in our analysis.

Canada is now one of the most widely vaccinated countries in the world with respect to SARS-CoV-2, and vaccination has undoubtedly blunted the impact of emergence of these VOCs. We find a marked reduction in risk of severe disease and death among both partially and fully vaccinated individuals in our study, consistent with the findings of Nasreen and colleagues (19). Notably, our analysis is performed in a case-only dataset, so the effects we report here (approximately 80-90% protection against death), represent a substantial degree of protection conferred by vaccines even when they fail to prevent infection. Such direct protective effects may help reduce the health impacts of ongoing SARS-CoV-2 transmission in Ontario, even if herd immunity proves elusive due to the high reproduction numbers of VOC (5, 6).

Even after adjustment for vaccination status we observed downward temporal trends in case severity, with a week-on-week reduced risk of hospitalization, ICU admission and death of approximately 0.5% (0-1.2%), 1.9% (0.8-3.1%) and 3.9% (2.5% to 5.4%), respectively. We suspect this decrease may represent attenuation of illness by vaccination in individuals who would not have been considered partly or fully vaccinated using our definitions, due to inadequate time elapsed since vaccination. This decrease may also represent residual confounding by changing age distribution of cases over time, perhaps not captured by our broad 10-year age categories, due to early prioritization of older individuals for vaccination in the province of Ontario (20).

An important limitation in our analysis is the possibility of our having misclassified early Delta variant infections as non-VOC due to the absence of routine screening for characteristic Delta mutations, with likely underestimation of the prevalence of Delta prior to May 2021. We caution that such misclassification may have biased our estimates of severity downwards, and our estimates of excess risk for Delta are likely biased towards the null.

In summary, we demonstrate that in the Canadian province of Ontario, despite widespread vaccination and VOC infections occurring more frequently in younger and healthier individuals, VOCs are associated with a substantial increase in virulence, including increased risk of death. The Delta variant is more virulent than previously dominant N501Y-positive VOCs. Combined with increased transmissibility and immune escape, these VOCs represent a significant escalation in risk to public health during the SARS-CoV-2 pandemic.