Data collection spanned from 2020 to 2021. A total of 819 individuals in decent health volunteered for the study in 2020, including 474 discharged patients with symptomatic COVID-19, 26 with asymptomatic COVID-19, and 319 controls (Figure 1). For those with COVID-19, initial testing was performed 3 to 173 days (within 6 months, 90% interval: 45-155 days) after discharge. 5 symptomatic patients and 19 controls were excluded for reasons outlined in Figure 1. We focused our analyses on the remaining 469 symptomatic patients, subdivided into 94 severe (82 severe, 12 critical) cases and 375 nonsevere (170 mild, 205 moderate) cases 19, and 300 controls. Their demographic and lifestyle characteristics are summarized in Table 1 and underlying conditions in Table S1. Consistent with earlier reports 20, significantly fewer COVID-19 patients were regular smokers (8.1% vs. 33.7%, χ21 = 80.75, p < 0.001) or drinkers (11.5% vs. 27.7%, χ21 = 32.61, p < 0.001) as compared with controls. More patients than controls reported that they had been taking medication for hypertension (16.8% vs. 8.0%, χ21 = 12.34, p < 0.001), hyperlipidemia (2.8% vs. 0.3%, χ21 = 6.09, p = 0.014), or digestive problems (2.1% vs. 0.3%, χ21 = 4.20, p = 0.040). Most of the patients and controls completed all 9 tests (see Procedures): odor threshold, odor discrimination (Table S2), odor identification (Table S3), gustatory function, trigeminal function, Montreal Cognitive Assessment (MoCA) 21,22, go/no-go task (inhibitory control, Figure S1A) 23, task-switching task (cognitive flexibility, Figure S1B) 24, and the short form Beck Depression Inventory (BDI-SF) 25. The valid response rate for a test was 96.4% on average. 75 of the patients were retested for odor threshold, odor discrimination, and/or taste identification –– tasks that are less susceptible to practice effects or seasonal changes –– 1 to 7 months after the initial evaluation, out of which 64 were retested once and 11 were retested twice for odor threshold and discrimination. Another 68 of the patients were followed up about 1 year after the initial evaluation, roughly 1.5 years after hospital discharge. Together with 53 newly enrolled COVID-19 patients who had been discharged for about 1.5 years, they were assessed for the following measures in 2021: odor identification, trigeminal function, MoCA, BDI-SF, as well as the generalized anxiety disorder 7-item scale (GAD-7) 26.
Table 1
Characteristics of enrolled COVID-19 patients and controls in the initial evaluation in 2020.
| COVID-19 Cases (N = 469) | Controls (N = 300) |
Severe (N = 94) | Nonsevere (N = 375) |
Age, No. (%) | | | |
≤ 39 yr. | 10 (10.6%) | 68 (18.1%) | 85 (28.3%) |
40 – 49 yr. | 10 (10.6%) | 76 (20.3%) | 49 (16.3%) |
50 – 59 yr. | 38 (40.4%) | 113 (30.1%) | 87 (29.0%) |
60 – 69 yr. | 30 (31.9%) | 103 (27.5%) | 63 (21.0%) |
≥ 70 yr. | 6 (6.4%) | 15 (4.0%) | 16 (5.3%) |
Female sex, No. (%) | 55 (58.5%) | 231 (61.6%) | 187 (62.3%) |
Regular smoker, No. (%) a | 7 (7.4%) | 31 (8.3%) | 101 (33.7%) |
Regular drinker, No. (%) a | 13 (13.8%) | 41 (10.9%) | 83 (27.7%) |
Education, mean (SD), yr. a | 12.5 (3.6) | 12.4 (3.2) | 12.0 (3.3) |
Length of hospital stay, median (SD), d. | 25 (15.5) | 17 (12.3) | |
Time since discharge, No. (%) | |
≤ 7 wk. (min. 3 d.) | 24 (5.1%) |
8 – 11 wk. | 98 (20.9%) |
12 – 15 wk. | 130 (27.7%) |
16 – 19 wk. | 75 (16.0%) |
20 – 23 wk. | 133 (28.4%) |
≥ 24 wk. (max. 173 d.) | 9 (1.9%) |
a Lifestyle and demographic factors that differed between the patients and controls in one or more age groups. Their effects were statistically accounted for in subsequent analyses. |
Based on the control data, age, sex, education, smoking, and drinking significantly contributed to performances on one or more of the chemosensory and cognitive tests (Table S4). Incidences of chronic diseases other than hypertension (13.4%) were low (< 5%) in our sample and none significantly affected performances. To quantify the impacts of COVID-19, we stratified the data by age group and compared the performances of the COVID-19 patients with those of the controls of the same age group, adjusted for the effects of significant confounders including education, smoking and/or drinking whose distributions differed between the patient and control groups (Table 1). Specifically, we standardized the continuous dependent data of all participants, adjusted for the significant confounders, based on the distributions of the bootstrapped control data per age group by z-transformation, so that for any given continuous dependent variable, the average z score of the controls fell around 0, whereas that of the patients reflected their relative standing with respect to matched controls (see Statistical analyses). Data from the 26 asymptomatic individuals (10 males, mean age ± SD = 45.6 ± 13.4 years) were separately analyzed in the same manner and presented in Figure S2.
Impaired olfactory, gustatory, and nasal trigeminal functions in discharged COVID-19 patients. Out of the 469 symptomatic patients, who were tested within 6 months (mean ± SD = 105.4 ± 37.4 days) of hospital discharge, 45.2%, 50.5%, and 70.4% subjectively reported they experienced losses of smell, taste, and appetite, respectively, when they had COVID-19. Many believed they had fully recovered by the time of testing: only 29.2%, 27.7%, and 20.5% of them considered their sense of smell, sense of taste, and appetite, respectively, to be poorer than before they contracted COVID-19, and only 19.4% rated their sense of smell to be below average. Their performances on the chemosensory tests, however, revealed a different picture (Figure 2A). Comparisons of the normalized z-scores between the patients and controls, which factored out the effects of significant demographic and lifestyle confounders and allowed direct comparisons to be made across tests, showed that these discharged patients performed significantly worse than the controls on odor threshold (Δz = -0.56, t694.7 = -7.09, p < 0.001), discrimination (Δz = -0.32, t681.2= -4.07, p < 0.001), and identification (Δz = -0.36, t760.7= -3.69, p < 0.001), with odor threshold (i.e., olfactory sensitivity) being the most impaired numerically (F1.9, 861.4 = 5.40, p = 0.006). Their overall olfactory function, as indexed by the mean of the z-scores for threshold, discrimination, and identification, was 0.41 SD below that of the controls (t748.9= -7.27, p < 0.001, Cohen’s d = 0.52). They also showed less differentiation of odor valence and rated the 18 odorants used in the odor identification test (Table S3) as more similar in valence (t753 = -2.21, p = 0.027). The pattern was found in those with an odor identification accuracy above 88% as well (≥ 16/18, t483 = -2.08, p = 0.038) and hence was unlikely a mere reflection of poor differentiation of odor quality. Eight patients, all aged below 40 years, reported signs of parosmia: 5 noted things smelled similarly unpleasant –– stinky, metallic, or ammonia-like; 3 noted several specific odors smelled different from before, like coffee, cola, milk, or celery. Impairments of gustatory and trigeminal functions were likewise apparent. The discharged patients scored 0.44 SD below the controls on taste identification (t735.7 = -4.61, p < 0.001, Cohen’s d = 0.33), and were particularly poor at identifying sourness (B = -1.25, p < 0.001) and saltiness (B = -0.79, p = 0.012) (Figure 2A inset), the two tastes mediated by ion channels, presumably due to reduced sensitivities (reflected in intensity ratings) to acid (Δz = -0.25, t727 = -3.35, p = 0.001) and salt (Δz = -0.29, t723 = -3.46, p = 0.001). Their nasal trigeminal sensitivity (indexed by ethanol lateralization accuracy) fell 0.32 SD below that of the controls (t705 = -3.86, p < 0.001, Cohen’s d = 0.30). These adverse impacts of COVID-19 on olfactory function (p < 0.001), taste identification (p = 0.014), and nasal trigeminal sensitivity (p = 0.019) held true for patients who considered their sense of smell to be at least average and their sense of smell, sense of taste, and appetite to be as good as or better than before they caught COVID-19, which underscores the necessity of objective testing in evaluating chemosensory functions. Moreover, smell loss appeared to persist after viral clearance even in asymptomatic individuals, and significant impairment of olfactory function was evident in our small sample of 26 discharged asymptomatic individuals (Δz = -0.31, t319 = -2.30, p = 0.022, Figure S2).
Table 2
Contributions of disease severity and time since discharge to chemosensory functions, cognitive performances and depressive state within 6 months of hospital discharge. The regression analyses on the raw data included demographic and lifestyle factors as additional regressors, see also Table S5. Numbers in square brackets: lower and upper bounds; B: unstandardized regression coefficient.
| Severity (severe vs. nonsevere) | Days since discharge |
| B (95% CI) | p | B (95% CI) | p |
Chemosensory functions | | | | |
Odor threshold [0-20] | -0.93 (-1.85, -0.02) | 0.046 | 0.02 (0.008, 0.03) | <0.001 |
Odor discrimination [0-10] | -0.07 (-0.50, 0.36) | 0.75 | -0.001 (-0.005, 0.004) | 0.76 |
Odor identification [0-18] | -0.71 (-1.34, -0.08) | 0.028 | 0.001 (-0.006, 0.008) | 0.76 |
Taste identification [0-5] | -0.13 (-0.30, 0.04) | 0.12 | 0.003 (0.001, 0.004) | 0.006 |
Trigeminal lateralization [0-10] | 0.07 (-0.35, 0.49) | 0.74 | 0.007 (0.003, 0.01) | 0.002 |
Cognitive performances | | | | |
MoCA [0-30] | -0.28 (-0.98, 0.42) | 0.44 | -0.001 (-0.008, 0.007) | 0.89 |
Go/No-go d’ | -0.09 (-0.29, 0.11) | 0.38 | -0.001 (-0.003, 0.001) | 0.46 |
Switch cost, ms. | -4.56 (-53.83, 44.70) | 0.86 | 0.47 (-0.07, 1.00) | 0.085 |
Depressive state | | | | |
BDI-SF [0-39] | 1.25 (0.05, 2.44) | 0.041 | -0.02 (-0.03, -0.003) | 0.017 |
Impaired cognitive functions and elevated depressive symptoms in discharged COVID-19 patients. More alarmingly, impairments in several domains of cognition were evident in the discharged COVID-19 patients (Figure 2A). They were overall 0.26 SD below the controls on MoCA (t749 = -3.19, p = 0.001, Cohen’s d = 0.24), a brief test developed to detect mild cognitive impairment (MCI) 21,22, and performed significantly worse than the controls on items pertaining to alternation (cognitive flexibility) (B = -0.70, p < 0.001), language (Δz = -0.24, t677.4 = -3.05, p = 0.002), memory (Δz = -0.16, t749 = -2.07, p = 0.039), and marginally so on items pertaining to visuospatial ability (Δz = -0.15, t676.4 = -1.86, p = 0.063) (Table S4). In addition, they exhibited heightened impulsivity (decreased inhibitory control, d’, Δz = -0.17, t724 = -2.08, p = 0.038) –– despite reacting more slowly than the controls (t502.7 = -2.69, p = 0.007) –– in the go/no-go task and increased switch cost (Δz = 0.19, t685.8 = 2.09, p = 0.037) in the task-switching task, which, along with the reduced performance on the MoCA alternation item, pointed to compromised executive function. The gap in MoCA scores was wider for individuals aged 50 years and over (Δz = -0.33, t377.1 = -3.15, p = 0.002): 55.1% of the older patients fell below the recently recommended cutoff of 25 for MCI 27, as opposed to 41.9% in the age-matched controls, the difference being 13.2% (χ21 = 7.25, p = 0.007). Whereas caution has been suggested in applying a specific MoCA cutoff score and the cutoff of 25 was likely associated with a high rate of false positives in our sample 28,29, it was clear that COVID-19 was associated with impairments in alternation/cognitive flexibility (p < 0.001), memory (p = 0.011), language (p = 0.036), and visuospatial ability (p = 0.032), and substantially increased the risk for MCI in those 50 years and older. Note that the patients were otherwise in decent health and, like the controls, had no history of dementia, stroke, or neurotraumatic, neurodegenerative, or neuropsychiatric diseases (see inclusion criteria).
Depression was also common in the discharged COVID-19 patients, who on average scored 0.37 SD above the controls on BDI-SF (t751.4 = 4.36, p < 0.001; Figure 2A) – a reliable measure of depressive symptoms 25. 16.8% of them met the screening criterion (≥10) and 7.7% met the diagnostic criterion for depression (≥14) 30, as compared with 8.0% and 3.0% in the controls, χ21s = 12.34 and 7.26, ps < 0.001 and = 0.007, respectively. Depression has been linked with deficits in olfaction and cognition 31,32. Indeed, we observed significant negative correlations between the normalized z-scores for BDI-SF and those for olfactory function and MoCA in both the patients (r458 = -0.21 and r457 = -0.15, ps < 0.001 and = 0.002, respectively) and the controls (r295 = -0.21 and r294 = -0.15, ps < 0.001 and = 0.010, respectively), indicating that the aforementioned impacts of COVID-19 were partly mediated by the elevated depressive symptoms in the patients. Critically, however, after partialling out the influence of depression level (BDI-SF), the adverse effects of COVID-19 remained significant for olfactory function (overall: F1,750 = 34.21, p < 0.001; threshold: p < 0.001, discrimination: p = 0.001, identification: p = 0.007), taste identification (F1,735 = 15.80, p < 0.001), nasal trigeminal sensitivity (F1,704 = 11.78, p = 0.001), MoCA (F1,748 = 6.41, p = 0.012 for all participants; F1,451 = 5.83, p = 0.016 for those 50 years and older), task switch cost (F1,694 = 4.41, p = 0.036), and marginally significant for inhibitory control (F1,723 = 2.75, p = 0.098). In other words, pathological factors specific to COVID-19 and independent of depressive symptoms associated with COVID-19 significantly impair chemosensory and cognitive functions, and the impacts last long after viral clearance. On the other hand, with the effect of depression level statistically accounted for, the patients’ z scores for olfactory function remained significantly correlated with those for MoCA (r446 = 0.26, p < 0.001) and inhibitory control (r440 = 0.17, p < 0.001), corroborating the link between olfactory dysfunction and cognitive impairment 17. Further inspection indicated that the correlations were driven by odor discrimination and identification (combined z scores, rs = 0.26 and 0.18, respectively, ps < 0.001), which tap into higher levels of olfactory processing, rather than odor threshold (rs = 0.14 and 0.037, ps = 0.003 and 0.43, respectively).
Slow and steady recovery of basic chemosensory but not cognitive functions independent of disease severity. We wondered whether the observed adverse impacts of COVID-19 were related to disease severity and whether they would diminish over time. As an initial step, we performed multiple regressions on the raw patient data collected within 6 months of hospital discharge, firstly using severity, days since discharge, as well as age, sex, education, smoking and drinking as the regressors (Tables 2 and S5), and then introducing BDI-SF score as an additional regressor to partial out the effect of depressive state. Results showed that the patients with severe COVID-19 performed worse than those with nonsevere COVID-19 on odor threshold (β = -0.091, p = 0.046) and identification (β = -0.10, p = 0.028), and suffered from a more severe depressive state (β = 0.093, p = 0.041). The differences in smell became marginally significant after controlling for depression level (0.05 < ps < 0.09). On the other hand, basic chemosensory functions like odor threshold (β = 0.16, p < 0.001), taste identification (β = 0.13, p = 0.006) and nasal trigeminal lateralization (β = 0.15, p = 0.002) improved steadily over time; these improvements persisted after adjusting for depression level (ps < 0.009). No significant effect of disease severity or days since discharge was detected for measures of cognitive functions including MoCA score, inhibitory control in the go/no-go task, and switch cost in the task-switching task. We next compared the normalized performances of the patients, classified by disease severity and time since discharge, respectively, and those of the controls, so as to further dissect and characterize the roles of severity and time (Figure 2A and B).
Various serious central neurological manifestations including stroke and psychosis have been reported in COVID-19 patients who developed acute respiratory distress syndrome and were admitted to intensive care units 7. In our original sample of 469 symptomatic patients, only 12 (2.56%) had critical COVID-19. Although we did not detect statistically significant differences in cognitive functions between those with severe (severe and critical) and nonsevere (mild and moderate) COVID-19, this did not rule out the possibility that patients with critical disease are more prone to cognitive decline. What was striking, however, was that cognitive impairments were evident in those who had nonsevere COVID-19, no underlying health conditions, and were on average tested 106.2 days ( > 3 months, SD = 36.7 days) after hospital discharge –– After partialling out the influence of depression level, they still scored significantly lower on MoCA (Δz = -0.20, F1, 557 = 5.05, p = 0.025) and exhibited reduced inhibitory control (Δz = -0.22, F1, 536 = 5.57, p = 0.019), as well as compromised olfactory (Δz = -0.32, F1, 554 = 24.04, p < 0.001), gustatory (taste identification: Δz = -0.34, F1, 544 = 9.53, p = 0.002) and nasal trigeminal functions (lateralization: Δz = -0.24, F1, 523 = 6.75, p = 0.010), relative to the controls. The gap in MoCA score was again more pronounced for individuals aged 50 years and over (Δz = -0.32, F1, 295 = 6.21, p = 0.013).
More concerning was that cognitive functions, unlike basic chemosensory functions, showed no sign of improvement over time. The patients tested on week 20 or more (mean = 151.4 days ≈ 5 months) significantly outperformed those tested within 11 weeks of discharge (mean = 58.9 days ≈ 2 months) on odor threshold (Δz = 0.54, t261 = 3.78, p < 0.001), taste identification (Δz = 0.49, t222.5 = 2.52, p = 0.012) and nasal trigeminal lateralization (Δz = 0.31, t237 = 2.19, p = 0.030), but not on MoCA (t251 = -0.37, p = 0.71), inhibitory control (t246 = -0.059, p = 0.95), or task switch cost (t235 = 1.44, p = 0.15). Meanwhile, their performances on odor threshold (Δz = -0.34, t235.8 = -2.90, p = 0.004), discrimination (Δz = -0.30, t432 = -2.83, p = 0.005), and identification (Δz = -0.35, t190.3= -2.28, p = 0.024) still fell below those of the controls. That is, full recovery was not attained for all chemosensory functions by 5 months of discharge.
We went on and examined data collected after a longer period of recovery. In the subset of 75 patients (52.3 ± 10.3 years) with retest data for odor threshold, odor discrimination, and/or taste identification, 74.7% and 56.8% showed improvements in odor threshold and discrimination, respectively, over the initial test, and 73.7% of those who initially misidentified at least one tastant showed an improvement in taste identification. 66 were retested around 9 months (275.9 ± 23.2 days) after discharge. By then, their performances on these chemosensory tasks were overall comparable to those of the controls (odor threshold: t363 = 1.58, p = 0.12; odor discrimination: t361 = 0.44, p = 0.66; taste identification: t74.4 = -1.21, p = 0.23; Figure 2B and C), although one 46-year-old woman remained anosmic (mild case, threshold = 0, discrimination at chance) and one 62-year-old woman (severe case) failed to recognize bitter, salty, and umami, out of the five basic tastes. 11 were retested twice for odor threshold and discrimination, once 118.7 ± 23.6 days (~3.9 months) and once 272.0 ± 26.1 days (~8.9 months) after discharge. Data from their initial (80.7 ± 26.9 days or ~2.6 months after discharge) and follow-up tests (Figure S3) partially echoed the patterns seen in Figure 2B (leftmost panel) and suggested that the recovery of olfactory sensitivity was steady and relatively faster in the first few months following acute COVID-19, whereas that of olfactory discrimination was a more heterogenous process.
68 patients (58.3 ± 6.4 years) were retested for odor identification, trigeminal lateralization, cognitive function, and mood state about a year (13.8 ± 0.7 months) following the initial test –– roughly 1.5 years (17.0 ± 0.8 months) after hospital discharge (Figure 2C). 50.7% and 53.2% of them showed improvements in odor identification and nasal trigeminal sensitivity, respectively. There was however no improvement in MoCA score (t64 = -0.67, p = 0.51) or depressive symptoms (t67 = 0.008, p = 0.99). Poor performance on MoCA (t341 = -5.17, p < 0.001) and an elevated level of depression (t60.4 = 3.30, p = 0.002) were also evident in an independent group of 53 COVID-19 patients (34 females, 59.5 ± 7.4 years) enrolled about 1.5 years (17.2 ± 0.9 months) after hospital discharge. Collectively, these 121 patients, mostly nonsevere cases (81.8%) and over 50 years (94.2%), scored 0.68 SD above the controls on BDI-SF (t153.1 = 4.01, p < 0.001) and 0.54 SD above the patients tested at 5 months (20-23 weeks) after discharge (t199.9 = 2.90, p = 0.004). 24.0% of them met the screening criterion (≥10 on BDI-SF) and 14.0% met the diagnostic criterion for depression (≥14 on BDI-SF) 30, whereas 14.0% met the criterion for general anxiety disorder (≥10 on GAD-7) 26. A recent cohort study 11 noted a similar pattern –– that depression or anxiety was more frequently reported at 12 months than 6 months following acute COVID-19. There was a strong correlation between levels of depression (BDI-SF) and anxiety (GAD-7) (r121 = 0.75, p < 0.001), in line with the documented high comorbidity of depressive and anxiety disorders 33. After controlling for the effect of mood state, no statistically significant difference was detected between the patients and controls in olfactory identification (p = 0.37) or trigeminal lateralization (p = 0.63). Nonetheless, the patients still fell 0.43 SD below the controls on MoCA (F1, 406 = 13.70, p < 0.001), which pointed to a worrisome long-lasting toll of COVID-19 on cognitive abilities.