Main

The coronavirus disease 19 (COVID-19) pandemic has seen the development and deployment of vaccines at an unprecedented speed and scale. Several vaccines including ChAdOx1nCoV-19 and BNT162b2 are approved for use in multiple countries and these have been shown to reduce COVID-19 infections, transmissions, hospitalizations and deaths in randomized controlled trials and real-world effectiveness studies1,2,3,4,5,6. However, the clinical trials were underpowered to detect rare adverse events3,4 that are important for ongoing risk–benefit evaluations of these vaccines and for informing post-vaccination clinical practice. Therefore, the identification of such rare adverse events is now a global scientific priority.

The increased risk of cerebral venous sinus thrombosis following the ChAdOx1nCoV-19 vaccine is an example of a rare adverse neurological event7,8. These findings have prompted the United Kingdom, several European countries and two Canadian provinces to limit the use of the ChAdOx1nCoV-19 vaccine or restrict its use pending further risk–benefit analysis in those at low risk of severe outcomes from infection. Furthermore, two cases of transverse myelitis were identified in the treatment arm during ChAdOx1nCoV-19 clinical trials3. These cases triggered two temporary study pauses, including careful regulatory review by the UK’s Medicines and Healthcare Products Regulatory Agency (MHRA). One case of transverse myelitis was considered to be possibly causally related to the vaccine, and an association of rare neuroinflammatory side-effects with ChAdOx1nCoV-19 could not be ruled out.

Since the start of large-scale vaccine programs across the world, additional case reports have linked other neurological adverse events to COVID-19 vaccination, including Guillain–Barré syndrome9,10,11. Furthermore, surveillance studies have found a possible link between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and neurological events, including Guillain–Barré syndrome and myelitis12,13. However, case reports and surveillance studies are limited by small numbers, as well as potential selection and recording biases. Therefore, detailed assessments of potential neurological adverse events associated with COVID-19 vaccines and infection are urgently needed.

The English National Immunisation (NIMS) Database of COVID-19 vaccination includes data on vaccine type, date and doses for all people vaccinated in England. We linked NIMS, at the individual patient level, to national data for mortality, hospital admissions and SARS-CoV-2 infection data to examine the associations between the first dose of ChAdOx1nCoV-19 or BNT162b2 vaccines and neurological complications: acute central nervous system (CNS) demyelinating events, encephalitis meningitis and myelitis, Guillain–Barré syndrome, Bell’s palsy, myasthenic disorders, hemorrhagic stroke and subarachnoid hemorrhage. We use the same population to investigate the associations between a positive SARS-CoV-2 test as a secondary exposure and the same neurological conditions. Incidence rate ratios (IRRs), that is, the rates of hospital admission or death from each neurological outcome in the risk period after vaccination or after a positive test relative to the baseline period, were estimated using the self-controlled case series (SCCS) methodology (Fig. 2)14,15. A relative incidence greater than 1 indicates an increased risk after vaccination or a SARS-CoV-2-positive test. We also carried out an independent replication of the risk of neurological outcomes in a national cohort from Scotland using the SCCS design.

Results

Overall, 32,552,534 people received their first dose of COVID19 vaccine (ChAdOx1nCoV-19, n = 20,417,752; BNT162b2, n = 12,134,782) in England between 1 December 2020 and 31 May 2021. In the population of vaccinated people, 2,005,280 (~6%) had a SARS-CoV-2-positive test. Of those with a positive test, 1,838,628 (~91%) had their SARS-CoV-2 test prior to vaccination. The baseline characteristics of the study population are listed in Table 1.

Table 1 Baseline demographic characteristics of people receiving either ChAdOx1nCoV-19 or BNT162b2 vaccines or testing positive to SARS-CoV-2 virus (in those vaccinated with either vaccine), in England between 1 December 2020 and 31 May 2021

The results for the investigated neurological outcomes following each exposure (ChAdOx1nCoV-19, BNT162b2 and positive SARS-CoV-2 test) are given for pre-specified time periods: 0, 1–7, 8–14, 15–21 and 22–28 days, as well as collectively for the 1–28 day period. The baseline characteristics of people who had these neurological outcomes are listed in Table 2. The IRRs for each outcome in the risk periods immediately before and after each exposure are given in Table 3 and Fig. 1.

Table 2 Demographic characteristics of patients who had the individual outcomes in the 28 days following a COVID-19 vaccine first dose or SARS-CoV-2 infection in the vaccinated population in England from 1 December 2020 to 31 May 2021
Table 3 IRR and 95% CI for individual outcomes in pre-defined risk periods immediately before and after exposure to vaccination and before and after a positive SARS-CoV-2 test result, adjusted for calendar time from 1 December 2020 to 31 May 2021
Fig. 1: IRRs and 95% CIs for neurological outcomes following ChAdOx1nCoV-19 vaccination, BNT162b2 vaccination and positive SARS-CoV-2 test.
figure 1

The IRRs are presented for pre-defined risk periods (0, 1–7, 8–14, 15–21 and 22–28 days) after each exposure and for the pre-risk period (28 days prior to exposure) and computed using a population of n = 32,553,534 vaccinated individuals. The horizontal bold line indicates an IRR of 1.

Neurological outcomes

Acute CNS demyelinating events

A total of 1,105 people had a hospital admission with an acute CNS demyelinating event. Of these, 131 (11.9%) had a positive SARS-CoV-2 test including 121 (10.9%) who had a positive test prior to vaccination. There were fewer than five deaths overall. There was no association between the ChAdOx1nCoV-19 or BNT162b2 vaccination and admission with an acute CNS demyelinating event in any of the pre-defined 7 day risk periods. There was an increased risk of hospital admission or death for acute CNS demyelinating events on the day of a positive SARS-CoV-2 test (IRR, 19.34; 95% CI: 8.63–43.38), and 1–7 days (IRR, 3.45; 95% CI: 1.67–7.14) and 8–14 days (IRR, 2.61; 95% CI: 1.17–5.84) after a positive test. In the 1–28 days post-exposure period we did not observe any association with ChAdOx1nCoV-19 (IRR, 0.97; 95% CI: 0.78–1.22) or with BNT162b2 (IRR, 1.02; 95% CI: 0.75–1.40), but a possible association with SARS-CoV-2 infection was observed (IRR, 1.67; 95% CI: 0.93–3.00).

Encephalitis, meningitis and myelitis

Encephalitis-, meningitis- and myelitis-related admissions occurred in 1,285 people. Of these, 255 (19.8%) had a SARS-CoV-2-positive test including 213 (16.6%) who had a positive test prior to vaccination. There were 39 deaths (in six of which the SARS-CoV-2 test was positive). There was a trend towards increased risk of encephalitis, meningitis and myelitis after ChAdOx1nCoV-19 vaccination (IRR, 1.32; 95% CI: 0.99–1.76 at 8–14 days) and no association with the BNT162b2 vaccine in any of the 7 day risk periods. There was an increased risk of hospital admission or death for this outcome on the day of a SARS-CoV-2-positive test (IRR, 38.57; 95% CI: 23.41–63.56), and 1–7 days (IRR, 5.71; 95% CI: 3.49–9.32) and 8–14 days (IRR, 5.63; 95% CI: 3.64–9.15) after the positive test. In the 1–28 days post-exposure period we did not observe an association with ChAdOx1nCoV-19 (IRR, 1.07; 95% CI: 0.87–1.31) or with BNT162b2 (IRR, 1.14; 95% CI: 0.86–1.51), but an association with SARS-CoV-2 infection was identified (IRR, 2.07; 95% CI: 1.78–4.11).

Guillain–Barré syndrome

Guillain–Barré syndrome-related admissions occurred in 622 people. Of these, 110 (17.7%) had a positive SARS-CoV-2 test including 99 (15.9%) who had a positive test prior to vaccination. There were 11 deaths (SARS-CoV-2 test was not positive in any of these deaths). There was an increased risk of hospital admission or death for Guillain–Barré syndrome at 15–21 days (IRR, 2.90; 95% CI: 2.15–3.92) and 22–28 days (IRR, 2.21; 95% CI: 1.59–3.09) following the first dose of ChAdOx1nCoV-19. No association was found with the BNT162b2 vaccine. There was an increased risk of Guillain–Barré syndrome following a positive SARS-CoV-2 test (IRR, 33.37; 95% CI: 14.21–78.36 at day 0; IRR, 7.36; 95% CI: 3.57–15.18 at 1–7 days; IRR, 5.19; 95% CI: 2.31–11.65 at 8–14 days; IRR, 6.89; 95% CI: 3.37–14.09 at 15–22 days; IRR, 3.51; 95% CI: 1.44–8.57 at 22–28 days). In the 1–28 days post-vaccination period we found an increased risk of Guillain–Barré syndrome after ChAdOx1nCoV-19 vaccination (IRR, 2.04; 95% CI: 1.60–2.60) but not after BNT162b2 vaccination (IRR, 0.86; 95% CI: 0.54–1.36). Also, an increased risk of Guillain–Barré syndrome was observed in the 1–28 days period after a SARS-CoV-2-positive test (IRR, 5.25; 95% CI: 3.00–9.18).

Bell’s palsy

There were 3,249 people admitted to hospital with Bell’s palsy. Of these, 391 (12.0%) had a positive SARS-CoV-2 test including 334 (10.3%) who had a positive test prior to vaccination. There were no deaths. There was an increased risk of hospital admission for Bell’s palsy at 15–21 days after the first dose of ChAdOx1nCoV-19 (IRR, 1.29; 95% CI: 1.08–1.56). No association was found with the BNT162b2 vaccine. There was an increased risk of Bell’s palsy following a positive SARS-CoV-2 test (IRR, 33.23; 95% CI: 22.57–48.94 at day 0; IRR, 5.84; 95% CI: 4.09–8.33 at 1–7 days; IRR, 2.17; 95% CI: 1.30–3.63 at 8–14 days). In the 1–28 days post-exposure period we did not observe any association with ChAdOx1nCoV-19 (IRR, 1.07; 95% CI: 0.94–1.21), BNT162b2 (IRR, 1.06; 95% CI: 0.90–1.26) or SARS-CoV-2 infection (IRR, 1.34; 95% CI: 0.91–1.97). In the cohort of patients with Bell’s palsy, a small proportion (6%) had a suspected concurrent diagnosis of cerebral infarction. We undertook a post-hoc sensitivity analysis excluding these people but the results did not change from the main analysis, as shown in Supplementary Table 1.

Myasthenic disorder

There were 831 hospital admissions for myasthenic disorder. Of these, 137 people (16.5%) had a positive SARS-CoV-2 test including 110 (13.2%) who had a positive test prior to vaccination. There were 30 deaths (in 10 of which the SARS-CoV-2 test was positive). There was an increased risk of hospital admission or death for myasthenic disorder at 15–21 days after the first dose of ChAdOx1nCoV-19 (IRR, 1.57; 95% CI: 1.07–2.30). No association was found with the BNT162b2 vaccine. There was an increased risk of hospital admission or death with myasthenic disorder following a positive SARS-CoV-2 test (IRR, 61.32; 95% CI: 33.43–112.50 at day 0; IRR, 13.74; 95% CI: 8.07–23.41 at 1–7 days; IRR, 3.98; 95% CI: 1.86–8.52 at 8–14 days; IRR, 2.81; 95% CI: 1.16–6.79 at 22–28 days). In the 1–28 days post-exposure period we did not observe any association with ChAdOx1nCoV-19 (IRR, 1.23; 95% CI: 0.94–1.62) or BNT162b2 (IRR, 1.18; 95% CI: 0.88–1.59), but an increased risk after a SARS-CoV-2-positive test was identified (IRR, 3.01; 95% CI: 1.70–5.36).

Hemorrhagic stroke

There were 3,503 hospital admissions due to hemorrhagic stroke. Of these, 483 people (13.8%) had a positive SARS-CoV-2 test including 392 (11.2%) who had a positive test prior to vaccination. There were 803 deaths (in 83 of which the SARS-CoV-2 test was positive). No association was found with the ChAdOx1nCoV-19 vaccine in any of the pre-specified time periods. There was an increased risk of hospital admission or death from hemorrhagic stroke at 1–7 days (IRR, 1.27; 95% CI: 1.02–1.59) and 15–21 days (IRR, 1.38; 95% CI: 1.12–1.71) after the first dose of BNT162b2. There was an increased risk of hemorrhagic stroke up to 7 days after a positive SARS-CoV-2 test (IRR, 12.42; 95% CI: 7.73–19.95 at day 0; IRR, 2.01; 95% CI: 1.29–3.15 at 1–7 days). In the 1–28 days post-vaccination period we did not observe an association with ChAdOx1nCoV-19 (IRR, 0.91; 95% CI: 0.88–0.94), but a positive association with the BNT162b2 vaccine was observed (IRR, 1.24; 95% CI: 1.07–1.43). No association was found with SARS-CoV-2 infection in the 1–28 day period after a positive test (IRR, 0.85; 95% CI: 0.57–1.26).

Subarachnoid hemorrhage

There were 2,055 hospital admissions due to subarachnoid hemorrhage. Of these, 262 people (12.7%) had a positive SARS-CoV-2 test including 220 (10.7%) who had a positive test prior to vaccination. There were of 273 deaths (in 27 of which the SARS-CoV-2 test was positive). There was no association between the ChAdOx1nCoV-19 and BNT162b2 vaccines and subarachnoid hemorrhage. There was an increased risk of hospital admission or death for subarachnoid hemorrhage on the day of a SARS-CoV-2-positive test (IRR, 24.22; 95% CI: 14.50–40.45), and at 1–7 days (IRR, 4.17; 95% CI: 2.59–6.71) and 8–14 days (IRR, 2.15; 95% CI: 1.16–3.99) after the positive test. In the 1–28 days post-exposure period we did not observe any association with ChAdOx1nCoV-19 (IRR, 1.01; 95% CI: 0.86–1.18) or with BNT162b2 (IRR, 1.05; 95% CI: 0.84–1.30), but a potentially increased risk after the SARS-CoV-2 infection (IRR, 1.51; 95% CI: 0.96–2.36) was observed.

Guillain–Barré syndrome and Bell’s palsy: composite outcome

Given that facial weakness can be part of the spectrum of Guillain–Barré syndrome, we investigated the risk of hospital admission or death from co-occurrence of Bell’s palsy and Guillain–Barré syndrome in the same individuals before and after ChAdOx1nCoV-19 and BNT162b2 vaccination in the study period. There were 34 patients admitted to hospital for both Bell’s palsy and Guillain–Barré syndrome. Of these patients, 27 were admitted on the same day for both outcomes, while seven were first admitted for Bell’s palsy and subsequently for Guillain–Barré syndrome. We found an increased risk of co-occurrence of Bell’s palsy and Guillain–Barré syndrome in the 15–21 days (IRR, 28.86; 95% CI: 5.45–152.74) and 22–28 days (IRR, 13.35; 95% CI: 2.48–72.92) after ChAdOx1nCoV-19 vaccination (Supplementary Table 2). In the 1–28 days post-vaccination period we found a positive association with the ChAdOx1nCoV-19 vaccine (IRR, 12.66; 95% CI: 2.70–59.21) but not with the BNT162b2 vaccine (Supplementary Table 2). The IRRs associated with SARS-CoV-2 infection were not estimated due to the small number of events in those with a positive SARS-CoV-2 test.

National Scottish data

In Scotland 1,982,678 people received the first dose of ChAdOx1nCoV-19 vaccine (~10% of those in England) and 1,077,626 people received the first dose of BNT162b2 vaccine (~10% of those in England) between 1 December 2020 and 31 May 2021. Of those vaccinated, 117,554 (~3%) had a SARS-CoV-2-positive test, either prior to or after vaccination. The age distribution of vaccinated patients and the percentage of vaccinated patients having any prior conditions were similar in the Scottish and English populations (Supplementary Table 3a).

The association between ChAdOx1nCoV-19 and Guillain–Barré syndrome for the 1–28 days post-vaccination period was replicated in Scotland (IRR, 2.31; 95% CI: 1.02–5.24), as shown in Supplementary Table 3b. The timing of this increased risk was consistent with the results seen in England: IRR, 2.15 (95% CI: 0.63–7.31) for 8–14 days, 4.79 (95% CI: 1.72–13.38) for 15–21 days and 3.60 (95% CI: 1.11–11.65) for 22–28 days (Supplementary Table 3c). No association between BNT162b2 and hemorrhagic stroke was seen in Scotland (IRR, 0.65; 95% CI: 0.35–1.20) in the 1–28 days period after vaccination (Supplementary Table 3b). The findings for SARS-CoV-2 test positivity were similar to the English data, with an increased risk of neurological outcomes following a positive test, although the estimates were generally imprecise due to the smaller population and lower rates of infection.

Subgroup analyses by age group and sex

We tested for interactions between vaccination and age group (>50 years and ≤50 years) and sex for the outcomes for which we observed an increased risk after vaccination. Supplementary Table 4a lists the results of the interaction tests. When interaction terms were statistically significant, the IRR for subgroup analyses by age group or sex have been reported in Supplementary Table 4b. There was a significant interaction (P = 0.016) between age group and risk of myasthenic disorder associated with the ChAdOx1nCoV-19 vaccine (IRR, 2.44; 95% CI: 1.07–5.55 at 15–21 days) in people aged ≤50 years. The increased risk of hemorrhagic stroke associated with the BNT162b2 vaccine was significantly higher in female participants (P = 0.007): the IRR for hemorrhagic stroke in female participants was 1.44 (95% CI: 1.05–1.96) at 1–7 days and 1.84 (95% CI: 1.40–2.42) at 15–21 days after BNT162b2 vaccination compared with 1.13 (95% CI: 0.82–1.57) at 1–7 days and 0.98 (95% CI: 0.70–1.38) at 15–21 days in male participants.

Excess events due to exposures

We estimated the number of exposures needed for one excess event and the excess number of events per 10 million exposed for each outcome (Supplementary Table 5). For example, with ChAdOx1nCoV-19 there were 38 excess Guillain–Barré syndrome events per 10 million people vaccinated in the 1–28 days period after vaccination. For BNT1262b2, there were 60 extra cases of hemorrhagic stroke per 10 million people vaccinated in the 1–28 days period after vaccination. For SARS-CoV-2 infection, in the 1–28 days period after a positive test there were an estimated 123 extra events of encephalitis meningitis and myelitis and 145 of Guillain–Barré syndrome per 10 million people with a positive test.

Associations with negative and positive control outcomes

We examined the associations of exposures with celiac disease as a negative control outcome and with anaphylaxis as a positive control outcome. We found no increased risk of celiac disease (negative control) across the pre-specified time periods for the vaccine exposures, but a decreased risk on the day of vaccination. Anaphylaxis (positive control) showed the expected increased risk in the 0–7 days after the first dose for both vaccines, but not in later periods (Supplementary Table 6).

Sensitivity analyses

We conducted a number of sensitivity analyses to ascertain the robustness of our results. Overall, the main findings were not sensitive to censoring due to death (with the exception of hemorrhagic stroke), removal of patients who had a second dose of the vaccine or adjustment for potential delays in recording. Results for these sensitivity analyses are shown in Supplementary Table 7a and Extended Data Figs. 1 and 2. We also extended the exposure risk period to 5 weeks, to test the suitability of our 28 day risk period and to avoid missing events due to lag in diagnosis. Results were similar to the main analysis (Supplementary Table 7b). An increased risk of Guillain–Barré syndrome was observed in the 29–35 days after ChAdOx1nCoV-19 vaccination (IRR, 1.55; 95% CI: 1.03–2.34).

We have further restricted the analyses to patients who had a positive SARS-CoV-2 test before vaccination. The results did not change from the main analysis for the SARS-CoV-2 test exposure (Supplementary Table 7c), suggesting that there is a higher risk of neurological complications after a positive SARS-CoV-2 test for all outcomes. The number of people infected after vaccination was too small to carry out any additional meaningful analyses. We also restricted the analyses to those who did not have a positive SARS-CoV-2 test. The results did not change from the main analysis for both vaccine exposures (Supplementary Table 7d).

Discussion

This large population-based study of more than 32 million people investigated the neurological adverse events associated with the ChAdOx1nCoV-19 and BNT162b2 vaccines as well as SARS-CoV-2 infection. We identified several key findings that are timely and of high importance to the public, health policy makers and clinicians across the world. First, we found an increased risk of hospital admission for Guillain–Barré syndrome (15–21 days and 22–28 days), Bell’s palsy (15–21 days) and myasthenic disorders (15–21 days) in those who received the ChAdOx1nCoV-19 vaccine. Second, an increased risk of hospital admission for hemorrhagic stroke (1–7 days and 15–21 days) was observed in those who received the BNT162b2 vaccine. Last, we identified a much greater increase in the risk of neurological outcomes following a positive SARS-CoV-2 test, such as acute CNS demyelinating events, encephalitis meningitis and myelitis, Guillain–Barré syndrome, Bell’s palsy, myasthenic disorders, hemorrhagic stroke and subarachnoid hemorrhage. Apart from hemorrhagic stroke, the results were robust to sensitivity analyses that assessed the assumption that outcomes did not influence subsequent exposures, providing credence to our results.

In those who received the ChAdOx1nCoV-19 vaccine, we estimated 38 excess cases of Guillain–Barré syndrome per 10 million exposed in the 1–28 days risk period. However, we did not observe an increased risk of Guillain–Barré syndrome in those who received the BNT162b2 vaccine. We replicated these findings using an independent national cohort of more than 3 million people from Scotland, which provides strong support for the association between ChAdOx1nCoV-19 vaccine and Guillain–Barré syndrome. The study design and sensitivity analyses addressed some of the key issues related to this association, including confounding through fixed covariates, seasonality and variable exposures to predisposing infections during the pre- and post-vaccination periods. It remains unclear why the ChAdOx1nCoV-19 vaccine appears to contribute to the pathogenesis of Guillain–Barré syndrome while BNT162b2 vaccine does not. Further studies are needed to assess whether antibodies against adenovirus vector-based ChAdOx1nCoV-19 can cross-react with components of the peripheral nerves. Furthermore, we found that Guillain–Barré syndrome and Bell’s palsy co-occur in those who received the ChAdOx1nCoV-19 vaccine. Clinically, this is likely to represent a variant of Guillain–Barré syndrome, in line with the emerging case reports of Guillain–Barré syndrome variants with facial weakness after the ChAdOx1nCoV-19 vaccine16,17.

In those who received the BNT162b2 vaccine, we estimated 60 excess cases of hemorrhagic strokes per 10 million exposed in the 28 days after vaccination. We did not observe an increased risk of hemorrhagic stroke in those who received the ChAdOx1nCoV-19 vaccine. However, the magnitude of this association between BNT162b2 and hemorrhagic stroke was reduced in sensitivity analyses that accounted for fatal events. Intriguingly, the increased risk of hemorrhagic stroke associated with BNT162b2 was substantially higher in female participants than in male participants, suggesting that sex may further modify this risk. The mechanisms underlying this disparity in risk between the two vaccines is not clear but several reports have suggested an increased risk of immune thrombocytopenic purpura in individuals receiving mRNA vaccines18,19, which in turn could contribute to major bleeding events. It is worth highlighting that the association between the BNT162b2 vaccine and hemorrhagic stroke was not replicated using Scottish data in the context of a low number of events, and thus further validation of this finding is warranted.

In line with previous surveillance studies14,15, we identified that all investigated neurological conditions were linked to SARS-CoV-2 infection itself. Specifically, we noted an excess of inflammatory disorders including encephalitis meningitis and myelitis (123 excess cases per 10 million exposed), myasthenic disorders (163 excess cases per 10 million exposed) and Guillain–Barré syndrome (145 excess cases per 10 million exposed) in the 1–28 days period after a positive test. Our post-hoc analysis restricting to those who had a positive test prior to vaccination did not change the results. Unfortunately, the number of people with a SARS-CoV-2-positive test after vaccination was too small to evaluate the risk of neurological complications in this group. The highest IRRs for these neurological outcomes were typically seen on day 0 (the day of the positive test). We included day 0 on its own as a risk period because hospital admission can trigger SARS-CoV-2 testing; however, this could overestimate or underestimate the overall risks associated with infection. The number of people infected in Scotland (3%) was much smaller than in England (~6%), reflecting differences in infection rates and containment measures in the two countries.

Findings from this study have clear clinical and public health implications. Crucially, we found that the risk of neurological complications from infection was substantially higher than the risk of adverse events from vaccinations in our population (for example, 145 excess cases versus 38 excess cases of Guillain–Barré syndrome per 10 million exposed in those who had a positive SARS-CoV-2 test and ChAdOx1nCoV-19 vaccine, respectively). This will need ongoing analysis and monitoring as younger people are vaccinated. Clearly, neurological complications from vaccination and infection reported in this study are rare. However, these rare complications can cause lifelong disability requiring long-term care. Collectively, these results provide timely and valuable information that can help to inform clinical decision making, including facilitating earlier diagnosis, as well as resource allocation for health-care provision. This may be particularly relevant for intensive care unit resource planning, given the potential need for prolonged admission to intensive care units for a proportion of patients with Guillain–Barré syndrome20.

This study had several strengths. First, the United Kingdom was an ideal place to carry out this study given that both the ChAdOx1nCoV-19 and BNT162b2 vaccines have been rolled out at speed and scale. Second, this was a population-based study of prospectively recorded medical data and it avoided the recall and selection biases linked to case reports. Third, the large sample size provided sufficient power to investigate these rare neurological outcomes that could not be assessed through clinical trials. Fourth, the SCCS study design (Fig. 2) removes potential confounding from fixed characteristics, and the breakdown of the study period into weekly blocks accounted for temporal confounding. Last, the key finding that the ChAdOx1nCoV-19 vaccine was associated with an increased risk of subsequent Guillain–Barré syndrome was independently replicated using Scottish national data.

Fig. 2: Schematic presentation of the SCCS study design.
figure 2

Each patient is followed from the index date to the study end date and is censored if death occurred or if they had a second dose of vaccine. Risk intervals (0, 1–7, 8–14, 15–21 and 22–28 days after exposure), the pre-risk interval (the 28 days prior to exposure) and the baseline periods (from study start to 29 days before exposure, and from 29 days after exposure to study end) are shown.

There were several limitations to this study. We examined the risks associated only with the first vaccine dose because the data on outcomes following second doses were limited at the time of this study. Furthermore, it was not possible to distinguish between Guillain–Barré syndrome variants including Miller Fisher syndrome in view of the coding used in the Hospital Episode Statistics database. Also, given that we investigated only hospital admissions and mortality, a proportion of patients with milder neurological disease may not have been included and the overall burden of neurological adverse events from vaccination and infection could be underestimated. Some hospital admissions for myasthenic disorders could represent exacerbations of previously well-controlled myasthenic illness given that we excluded only those with a prior hospital contact for these neurological conditions in the preceding 2 years. The Scottish replication cohort was smaller than the study cohort, which meant that estimates were less precise and we did not have sufficient statistical power to replicate findings with smaller effect sizes.

In summary, this population-based study identifies and quantifies several rare neurological adverse events that are specific to the ChAdOx1nCoV-19 and BNT162b2 vaccines, as well as SARS-CoV-2 infection. We believe that these findings are likely to be of relevance to other countries using these vaccines and it would be useful to replicate these results in similarly large datasets internationally. Clinicians should be aware of these rare complications, and the findings from this study will be paramount to policy makers in risk–benefit evaluations and health-care resource allocation. Importantly, the risks of adverse neurological events following SARS-CoV-2 infection are much greater than those associated with vaccinations, highlighting the benefits of ongoing vaccination programs.

Methods

Ethics approval

National Health Service Research Ethics Committee (NHS REC) approval was obtained from East Midlands-Derby Research Ethics Committee [reference 04/03/2021].

Data

We used the NIMS Database of COVID-19 vaccination to identify vaccine exposure. This includes vaccine type, date and doses for all people vaccinated in England. National SARS-CoV-2 infection data were obtained from Public Health England (PHE).

This study used International Classification of Diseases 10 (ICD-10) codes to define neurological outcomes of interest (Supplementary Table 8). The diagnoses were extracted from Hospital Episode Statistics (HES), which is a database containing details of all admissions to National Health System (NHS) hospitals in England

We linked NIMS, at the individual level, to national data for mortality (Office for National Statistics, ONS), hospital admissions (HES) and SARS-CoV-2 infection data (PHE). This information is now summarized in Extended Data Fig. 3.

Study design

The SCCS design was originally developed to examine vaccine safety14,15. The method includes only individuals who have had the outcome of interest (cases) and compares the incidence rate (in each case) of the outcome of interest during a time-limited exposed period (for example, immediately after vaccination) with rates during unexposed periods (that is, baseline). A relative incidence greater than 1 indicates an increased risk after vaccination or a positive SARS-CoV-2 test. Given that the analyses are conditional on each case, any fixed characteristics during the study period, such as sex, ethnicity or chronic conditions, are inherently controlled for. Any time-varying factors, such as seasonal variation, need to be adjusted for in the analyses.

A flowchart showing how each dataset used for the SCCS analysis of each outcome was created is presented in Extended Data Fig. 3.

Study period and population

We examined the associations between ChAdOx1nCoV-19 or BNT162b2 vaccines and neurological complications during the ongoing COVID-19 vaccination program in England (commenced on 8 December 2020). We also investigated the association between a SARS-CoV-2-positive test and the neurological diseases of interest in these vaccinated individuals. Separate analyses were carried out for each neurological outcome of interest. People were considered eligible for inclusion in each study cohort if they had received at least one vaccine dose, were at least 16 years old, and were admitted to hospital or died from the outcome of interest between 1 December 2020 and 31 May 2021 (last data update). Patients were followed up from the start (1 December 2020) to the earliest of either the end of the study period (31 May 2021), the date of second dose or when they died. Only the first outcome event during the study period was included. Patients with a hospital admission for the same outcome in the 2 years prior to the start of the study period were excluded. We excluded patients who received the mRNA-1273 (Moderna) vaccine given that the number of these patients was too small for meaningful analysis.

Outcomes

The neurological outcomes in this study were selected a priori based on a literature search for any reports of emerging neurological complications from SARS-CoV-2 infection and vaccination. These included acute CNS demyelinating events, encephalitis meningitis and myelitis (excluding encephalopathy), Guillain–Barré syndrome (encapsulating all Guillain–Barré syndrome variants including Miller Fisher syndrome), Bell’s palsy, myasthenic disorders, hemorrhagic stroke and subarachnoid hemorrhage. Other outcomes, such as peripheral neuropathies, were excluded because these conditions were unlikely to result in hospital admission (particularly within the 28 day risk period in this study). We used the ICD-10 codes to define each outcome, as listed in Supplementary Table 8. The outcomes were identified as the first hospital admission due to the event of interest or as a cause of death recorded on the death certificate with the event.

Exposures

The exposure variables were a first dose of the ChAdOx1nCoV-19 or BNT162b2 vaccines and first infection with SARS-CoV-2, defined as a positive COVID-19 polymerase chain reaction (PCR) test during the study period (either before or after vaccination). For all outcomes we defined the exposure risk intervals as the following pre-specified time periods: 0, 1–7, 8–14, 15–21 and 22–28 days after each exposure date, under the assumption that the adverse events under consideration are unlikely to be related to exposure beyond 28 days after the exposure. We assumed that patients who experienced an outcome before vaccination would probably delay vaccination until symptoms had improved. Therefore we included a pre-risk period in the analyses, defined as the 1–28 days immediately before the vaccination or SARS-CoV-2 test. The baseline period was then defined as extending from 1 December 2020 to 29 days before the exposure date, and from 29 days after the exposure date to 31 May 2021 or the censored date, if earlier. Hospital admissions for the events of interest can trigger COVID-19 testing. Such events may well be caused by SARS-CoV-2 infection, but the reverse causality involved in their detection induces bias. To reduce the bias, which could overestimate or underestimate the effect of infection, we allocated day 0 to a risk period on its own.

Seasonality and COVID-19 pandemic period

Some outcomes of interest have seasonal variation: for example, Guillain–Barré syndrome is more common in the winter months21. Moreover, hospital admissions were probably influenced by the pressure on the health systems due to COVID-19, which was not uniform during the pandemic period. To allow for these underlying seasonal effects, we split the observational study period into weeks and adjusted for week as a factor variable in the statistical models. In doing so, we estimated the relative incidence in exposure periods up to 28 days after vaccination compared with the unexposed periods, adjusted for underlying seasonality, which accounted for changes in hospital admission rates caused by the pandemic. If an exposure period is long in relation to the observation period, it might be impossible to distinguish between seasonal variation and vaccine effects. However, given that the exposure period was relatively short (28 days), this was unlikely to happen given that after 28 days a vaccinated patient returns to the baseline period, ensuring that in each interval both exposed and unexposed cases appear (Extended Data Fig. 4).

Statistical analysis

We described the characteristics of each cohort (vaccinated patients with the outcomes of interest) in terms of age and sex. The SCCS models were fitted using a conditional Poisson regression model with an offset for the length of the risk period. IRRs, the rates of hospital admission or death due to each outcome of interest in risk periods relative to baseline periods, and their 95% CI were estimated using the SCCS model adjusted for week number as time-varying covariates. Exposure terms for both vaccines and for infection with SARS-CoV-2 were included in the same model.

We conducted sensitivity analyses to assess the robustness of results to assumptions such as that the occurrence of an outcome event did not influence the probability of subsequent exposures, by (1) excluding those who died from the outcome and (2) restricting analysis to the post-vaccination period, without censoring at death due to the outcome; to assess potential reporting delays in the data by (3) restricting the study period until 1 May 2021; and to test the choice of the follow-up period by (4) restricting analysis to patients who received only one vaccine dose and (5) censoring at 12 weeks after vaccination or at the date of a second dose if earlier. To test our assumption that events after 28 days are not related to the vaccines, we did an additional analysis with a 35 day risk period after exposure.

To further validate the estimated IRR associated with a positive SARS-CoV-2 test, we also restricted the analyses to only patients who had a positive SARS-CoV-2 test before vaccination.

Stata version 17 was used for these analyses.

Absolute measures of effect

Absolute risk differences cannot be obtained using the SCCS methodology. We supplemented our estimates of IRRs with measures of the effect of each exposure in absolute terms, using a method22 developed to estimate the number of exposures needed to produce one excess adverse outcome and the excess number of events per 10 million exposed for each outcome.

Interactions with age and sex

We investigated whether the association between vaccine exposures and outcomes is sex or age dependent by including, in the analysis, interaction terms between these covariates and the exposures. We used the likelihood ratio test to assess the interaction terms. We present results in separate subgroups when an interaction was significant.

Negative and positive control

We examined the association of exposures with celiac disease as a negative control outcome23, which is assumed not to be associated with exposure to vaccination or SARS-CoV-2 infection; and with anaphylaxis as a positive control outcome given that it could occur shortly after vaccination with either vaccine24.

Assessing the SCCS assumptions

To further assess the assumptions of the SCCS and our modeling choices, we visually examined the data. We plotted a histogram of the number of occurrences of an event by time prior to or since vaccination for each outcome to assess the possibility that a hospital admission for that event might affect subsequent vaccination (Extended Data Fig. 5). We plotted a histogram of the time from the event to the actual end of observation in censored and uncensored cases to assess whether event-dependent observation periods may be a problem for the analysis (Extended Data Fig. 6).

Event-dependent exposures

Extended Data Fig. 5 shows the number of occurrences of an event by time before or since vaccination. We notice that there is a decrease in the 28 days immediately before vaccination, indicating that occurrence of an event might have reduced the likelihood of vaccination. This pattern is similar for most of the outcomes and for both vaccines. Therefore, we have added the pre-risk period of 28 days.

Event-dependent observation periods

Extended Data Fig. 6 shows the frequency of days from the event to the actual end of observation in censored and uncensored cases. A spike close to zero is apparent in the censored data histogram for hemorrhagic stroke and subarachnoid hemorrhage. This indicates the presence of event-dependent observation periods (censoring on death date due to outcome), which we tested further with sensitivity analyses 1 and 2. These additional analyses were in agreement with the main analysis, suggesting that there should be little concern for these outcomes.

Replication using national Scottish data

To assess the robustness of the findings, we carried out an independent replication of the risks of neurological outcomes in vaccine recipients using a similar SCCS design. NHS Scotland provides comprehensive health care that is free at the point of care for all residents, with a base population that includes 5.4 million residents (~95% of the population) registered with a general medical practice in Scotland. We used deterministic linkage (based on the Community Health Index) to link vaccination information (extracted from general medical practice records and the Turas Vaccination Management Tool) and hospitalization records (from the Scottish Morbidity Records 01), mortality data (from National Records of Scotland) and SARS-Cov-2 infection data (from Electronic Communication of Surveillance in Scotland).

We classified exposures and outcomes in a similar way to the main analysis. An individual was defined as exposed to either the ChAdOx1nCoV-19 or BNT162b2 vaccine on the date that they received their first dose between 8 December 2020 and 31 May 2021. We defined outcomes of interest using the same ICD-10 codes as for the main analysis. We excluded patients who experienced the outcome in the 2 years prior to the study start date and excluded patients under 16 years of age. We fitted conditional Poisson regression models with an offset for the length of the risk period and adjusted for the weekly time period as a time-varying covariate. The control (baseline) period was the period from 1 December 2020 to 29 days before vaccination, allowing for a 28 day clearance (that is, pre-risk) period, and from 29 days after the exposure date to 31 May 2021 or the censored date, if earlier. We also stratified the post-vaccination follow-up period using the same time periods as the analysis in England when the number of events allowed. We also repeated the analysis by replacing vaccination with SARS-CoV-2 infection, which was defined as a positive PCR test within the study period. Given the smaller population size of Scotland, we focused on assessing the consistency of the point estimates of associations observed in the main analyses.

Reporting Summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.