Abstract
Objective:
To examine differences in mental health–related healthcare utilization for anxiety and depression between individuals who did and did not undergo COVID-19 PCR testing in Ontario. Background: The COVID-19 pandemic has been associated with changes in mental health and healthcare utilization.
Methods:
We conducted a population-based retrospective cohort study using linked ICES data, including 6,175,114 adults (January 2020–March 2021). Exposure was PCR-positive, PCR-negative, or untested. The outcome was time to first mental health–related healthcare use for anxiety and depression, identified using validated codes. Adjusted hazard ratios (aHRs) were estimated using Cox models with propensity score matching.
Results:
Individuals who underwent testing had higher mental health–related healthcare utilization than untested individuals. This was observed in PCR-positive (aHR 6.37; 95% CI 6.25–6.50) and PCR-negative groups (aHR 5.91; 95% CI 5.87–5.95). Higher utilization occurred among younger individuals, females, and socioeconomically disadvantaged groups. Results were consistent in matched analyses.
Conclusion:
Individuals underwent testing had higher mental health service utilization; similar estimates across PCR-positive and PCR-negative groups suggest testing reflects underlying vulnerability and healthcare-seeking behavior rather than a causal effect on mental health outcomes.
Introduction
The COVID-19 pandemic and the restrictive measures implemented led to a surge in mental health issues globally [1], including anxiety and depression [2–4]. According to a World Health Organization (WHO) research brief, anxiety and depression increased by 25% in the first year of the pandemic [5]. A global study reported almost 76 million cases of anxiety and 53 million cases of major depressive disorder due to the pandemic, resulting in a 26% and 28% increase, respectively [6]. Mental health has long been recognized to be impacted by epidemics and pandemics; the mental issues caused by viral outbreaks have been dubbed a “parallel epidemic” [7]. The rise in the incidence and prevalence of anxiety and depression was reported as an outcome of the COVID-19 pandemic, and mental health professionals issued a warning of an “echo pandemic” [8] of mental health issues [9]. Concerns surrounding mental health during COVID-19 could be viewed as an epidemic instead of just a delayed outcome of the pandemic [10].
Economic and social mitigation strategies implemented in Canada and other developed countries during the COVID-19 pandemic had unintended consequences that disproportionately affected vulnerable populations. These measures exacerbated existing socioeconomic inequalities, increasing exposure to poverty and its associated mental health impacts, including anxiety and depression [11, 12]. Evidence indicates that low-income individuals, racial and ethnic minorities, and those with limited access to healthcare experienced a disproportionate burden of COVID-19 infection and adverse outcomes [13, 14]. These groups were also more likely to face economic instability, food insecurity, housing challenges, and barriers to mental health services [15]. The intersecting factors including socioeconomic status, pre-existing health conditions, access to care, living conditions, sex, and age, contribute to heightened vulnerability and widening disparities in mental health outcomes [11–13].
Testing for COVID-19 has played a crucial role in identifying and managing the spread of the virus, encompassing various dimensions that highlight its role as an exposure variable [14, 15]. Primarily, it serves as a proxy for direct exposure to the disease, with individuals seeking testing typically driven either by symptoms or close contact with confirmed cases, thus facing the risk of contracting COVID-19 disease. This inherent risk can induce anxiety and stress, as individuals struggle with the uncertainty surrounding their health status and the potential ramifications of a positive test result. The fear of a positive result and its implications on one’s health, livelihood, and relationships can significantly impact mental wellbeing. It is crucial to recognize and address the associated psychological distress that testing can induce. Anxiety, fear of the unknown, stigma, discomfort, information overload, and post-traumatic stress symptoms are among the psychological distress factors associated with COVID-19 testing [12, 16–18]. Overall, while testing may not capture all aspects of disease exposure, it serves as a practical and informative proxy variable that can provide valuable insights into the dynamics of disease transmission and its impact on populations.
Administrative health databases provide comprehensive, population-level data that enable the examination of healthcare utilization patterns during the COVID-19 pandemic. Beyond infection itself, the pandemic introduced multiple stressors including uncertainty, perceived risk of illness, social disruption, and changes in healthcare access that may influence mental health and healthcare-seeking behavior. COVID-19 testing represents a key point of interaction with the healthcare system and may therefore serve as a proxy for underlying psychological vulnerability, perceived risk, and engagement with care. Integrating COVID-19 testing data with other health information allows for a more comprehensive assessment of how these factors relate to mental health service utilization.
Therefore, the objective of this study is to examine the differences in mental health–related healthcare utilization for anxiety and depression between individuals who did and did not undergo COVID-19 PCR testing in Ontario, Canada.
Methods
Data source
This study utilized secondary data from administrative health databases accessed through the Ontario Health Database Platform (OHDP). The platform permitted access to patient-level linked population-based health administrative data held at the Institute for Clinical Evaluative Sciences (ICES). The ICES data repository consists of record-level, coded, and linkable health data sets that integrate health services data, registries, and surveys. ICES is an independent, non-profit research institute that houses routinely collected health data from Ontario’s publicly funded healthcare system.
All Ontario residents covered by the province’s universal health coverage were identified using the Registered Persons Database (RPDB), which provided age, sex, and demographic information. COVID-19 status was determined using C19INTGR, a comprehensive dataset of PCR COVID-19 diagnostic results in Ontario. Chronic conditions such as asthma, hypertension, diabetes, COPD, and congestive heart failure were identified using ICES-derived cohorts. These datasets were linked using unique encoded identifiers and analyzed at ICES (Supplementary Tables S3, S4). This study followed the Reporting of Studies Conducted Using Strengthening the Reporting of Observational studies in Epidemiology (STROBE) reporting guidelines [19]. The codes for identifying mental health visits are used for mental health system performance measurement in Ontario and elsewhere in Canada [20–24].
Study design, settings, and cohort eligibility
We constructed a retrospective cohort of people who were living in Ontario, Canada, had a valid health card, and for whom information on sex and age was available between the period of January 2020 to March 2021 (Figure 1). All those cases who were not Ontario residents or were less than 18 years of age, as well as those who had a standard COVID-19 lab-based PCR test result after the study end period, i.e., beyond 31 March 2021, were excluded. We utilized the National Ambulatory Care Reporting System (NACRS) and the Ontario Health Insurance Plan (OHIP) to identify physician, outpatient visits, and emergency department visits for mental health, and we used the Canadian Institute for Health Information’s Discharge Abstract Database (DAD) and Ontario Mental Health Reporting System (OMHRS) to capture psychiatric hospitalizations. Throughout the study period, COVID-19 testing guidelines were applicable to the general population in Ontario and were not limited to individuals considered at high risk. This helped to maintain a focus on the general population and ensure that our findings accurately reflect the broader impact of the pandemic on Ontario residents.
FIGURE 1
The exposure in this study was COVID-19 testing by standard lab-based PCR tests in Ontario between January 2020 and March 2021, irrespective of COVID-19 disease status, leading to three groups for comparison: (i) those PCR tested with a positive COVID result; (ii) those PCR tested with a negative COVID result, and (iii) those who were not PCR tested. The event of interest was the time in months until first recorded health services utilization related to anxiety and depression, using DSM-V criteria, ICD-10 codes, and OHIP diagnostic codes during the study period. Individuals who were not PCR tested were selected as the reference group to provide a population-level comparator reflecting mental health–related healthcare utilization among those who did not engage with COVID-19 testing services during the study period. This approach allows for the evaluation of differences associated with testing exposure within the broader population context. However, we acknowledge that untested individuals may differ systematically from those who were tested, particularly in terms of healthcare-seeking behavior, access to care, and underlying health concerns. To mitigate these differences, we applied propensity score matching to balance measured sociodemographic and clinical characteristics between groups. The participants were censored if there was no utilization of anxiety/depression-related health services by the end of the study follow-up period on 31 March 2021. Also, those who died during the study period were censored at the date of death, assuming the event of interest was not experienced first.
Variables
Clinical variables
COVID-19 status was based on a standard lab-based PCR test, while the prevalence of chronic conditions was identified using ICES-derived cohorts. Details of each chronic condition and algorithm used are described in the Supplementary Material (Supplementary Tables S3, S5).
Service utilization variables
We utilized an algorithm [25, 26] to identify cases of health service utilization among the Ontario population aged ≥18 years within the ICES data holdings, defined as: (i) one hospitalization in OMHRS or DAD with a primary discharge diagnosis of either anxiety, depression, and anxiety or depressive episode otherwise not specified (NOS) based on DSM-V criteria and ICD-10 codes, or (ii) at least one physician billing and/or emergency department visits in the OHIP claims database or NACRS with a diagnostic code for anxiety, depression, and anxiety or depressive episode otherwise not specified (NOS) based on DSM-IV criteria and ICD-10 codes. A complete list of the codes used to define the study cohort and the study variables is available in the Supplementary Material (Supplementary Tables S3, S5).
Sociodemographic variables
We examined sociodemographic variables potentially associated with health service utilization including age, sex, place of residence, and neighborhood-level income quintile using the Ontario Marginalization Index (ON-Marg) [27]. ON-Marg allows for a comprehensive exploration of multiple dimensions of marginalization in urban and rural areas of Ontario.
Statistical analysis
Descriptive statistics (e.g., mean, median, range, standard deviation) were calculated for the study sample (N = 6,175,114) (Figure 1). Frequencies for demographic and health service use variables were examined with cross-tabulations across both COVID-19 test-based groups using t-tests and chi-squared statistics (Table 1). Criteria for statistical significance were set at p < 0.05. We followed people with and without COVID-19 standard lab-based PCR tests from January 01, 2020, to March 31, 2021, to estimate survival time. Time was defined as the interval between the COVID-19 PCR test date and an event (i.e., first contact with health services for anxiety and depression during study period).
TABLE 1
| Demographics and morbidity information | Total sample | Individuals without PCR COVID-19 tests | PCR test with COVID negative result | PCR test with COVID positive result |
|---|---|---|---|---|
| Total n (%) | 6175114 (100) | 4274037 (69.21) | 1755519 (28.43) | 145558 (2.36) |
| Mean age yrs. (±SD) | 42.95 (15.49) | 43.84 (15.60) | 41.07 (15.12) | 39.39 (14.47) |
| Anxiety/depression related health service utilization | ||||
| No | 5511050 (89.25) | 3788384 (61.35) | 1587938 (25.72) | 134728 (2.18) |
| Yes | 664064 (10.75) | 485653 (7.86) | 167581 (2.71) | 10830 (0.18) |
| Age groups n (%) | ||||
| 18–45 years | 3569070 (57.80) | 2363642 (38.28) | 1,108,188 (17.95) | 97240 (1.57) |
| 46–65 years | 2014498 (32.62) | 1454906 (23.56) | 518139 (8.39) | 41453 (0.67) |
| 66–104 years | 591546 (9.58) | 455489 (7.38) | 129192 (2.09) | 6,865 (0.11) |
| Sex n (%) | ||||
| Male | 2985377 (48.35) | 2136749 (34.60) | 778344 (12.60) | 70884 (1.14) |
| Female | 3189737 (51.65) | 2137288 (34.61) | 977175 (15.82) | 75274 (1.22) |
| Rural n (%) | ||||
| No | 5580516 (90.37) | 3845375 (62.27) | 1594901 (25.83) | 140240 (2.27) |
| Yes | 594598 (9.63) | 428662 (6.94) | 160618 (2.60) | 5,318 (0.09) |
| Ontario marginalization index (%) | ||||
| Dependency quintile | ||||
| 1 | 1802531 (29.19) | 1230040 (19.92) | 519994 (8.42) | 52497 (0.85) |
| 2 | 1281491 (20.75) | 882985 (14.30) | 365722 (5.92) | 32784 (0.53) |
| 3 | 1078592 (17.47) | 749393 (12.14) | 305374 (4.95) | 23825 (0.39) |
| 4 | 1011615 (16.38) | 708060 (11.47) | 283600 (4.59) | 19955 (0.32) |
| 5 | 1000885 (16.21) | 703559 (11.39) | 280829 (4.55) | 16497 (0.27) |
| Deprivation quintile | ||||
| 1 | 1519683 (24.61) | 1042701 (16.89) | 451348 (7.31) | 25634 (0.42) |
| 2 | 1331132 (21.56) | 922903 (14.95) | 381542 (6.18) | 26687 (0.43) |
| 3 | 1169663 (18.94) | 810983 (13.13) | 330316 (5.35) | 28364 (0.46) |
| 4 | 1086791 (17.60) | 755009 (12.23) | 302288 (4.90) | 29494 (0.48) |
| 5 | 1067845 (17.29) | 742441 (12.02) | 290025 (4.70) | 35379 (0.57) |
| Instability quintile | ||||
| 1 | 1393327 (22.56) | 970766 (15.72) | 382819 (6.20) | 39742 (0.64) |
| 2 | 1186825 (19.22) | 828652 (13.42) | 333674 (5.40) | 24499 (0.40) |
| 3 | 1121085 (18.15) | 781645 (12.66) | 317085 (5.13) | 22355 (0.36) |
| 4 | 1076219 (17.43) | 740661 (11.99) | 312060 (5.05) | 23498 (0.38) |
| 5 | 1397658 (22.63) | 952313 (15.42) | 409881 (6.64) | 35464 (0.57) |
| Ethnicity diversity quintile | ||||
| 1 | 920772 (14.91) | 654464 (10.60) | 257398 (4.17) | 8,910 (0.14) |
| 2 | 998212 (16.17) | 691866 (11.20) | 293106 (4.75) | 13240 (0.21) |
| 3 | 1114619 (18.05) | 758660 (12.29) | 336714 (5.45) | 19245 (0.31) |
| 4 | 1349162 (21.85) | 919897 (14.90) | 398032 (6.45) | 31233 (0.51) |
| 5 | 1792349 (29.03) | 1249150 (20.23) | 470269 (7.62) | 72930 (1.18) |
| Asthma | ||||
| No | 6118287 (99.08) | 4240424 (68.67) | 1734137 (28.08) | 143726 (2.33) |
| Yes | 56827 (0.92) | 33613 (0.54) | 21382 (0.35) | 1832 (0.03) |
| COPD | ||||
| No | 6115375 (99.03) | 4232216 (68.54) | 1738483 (28.15) | 144676 (2.34) |
| Yes | 59739 (0.97) | 41821 (0.68) | 17036 (0.28) | 882 (0.01) |
| Hypertension | ||||
| No | 5841788 (94.60) | 4042477 (65.46) | 1662482 (26.92) | 136829 (2.22) |
| Yes | 333326 (5.40) | 231560 (3.75) | 93037 (1.51) | 8,729 (0.14) |
| Congestive heart failure | ||||
| No | 6156440 (99.70) | 4263386 (69.04) | 1747958 (28.31) | 145096 (2.35) |
| Yes | 18674 (0.30) | 10651 (0.17) | 7,561 (0.12) | 462 (0.01) |
| Diabetes | ||||
| No | 6028038 (97.62) | 4172160 (67.56) | 1715613 (27.78) | 140265 (2.27) |
| Yes | 147076 (2.38) | 101877 (1.65) | 39906 (0.65) | 5,293 (0.09) |
Demographic characteristics and morbidity information of the study population cohort in Ontario (Ontario, Canada, 2021).
Cox proportional hazards regression models were used to estimate the association between COVID-19 testing and time to mental health–related healthcare utilization for anxiety and depression, with results reported as hazard ratios (HRs) and 95% confidence intervals (CIs) for both the overall cohort and the propensity score–matched sample. Kaplan–Meier methods were used to estimate time-to-event functions, and differences between exposure groups were assessed using the log-rank test.
To address potential confounding by indication, propensity scores were calculated using age and sex and applied through matching with replacement, resulting in a matched sample of 664,233 individuals. This approach matched each exposed individual to the nearest comparator based on propensity score, allowing reuse of controls to improve balance. Covariates were selected a priori based on clinical relevance and prior literature, including sociodemographic characteristics, comorbidities, and measures of healthcare utilization, and were incorporated into both regression and propensity score models to enhance comparability between groups. All analyses were conducted using SAS Enterprise Guide 7.1 (SAS Institute, Cary, NC).
Ethics approval
This study was approved by the University of Waterloo research ethics board (protocol reference no. 43910), and data use was authorized under section 45 of Ontario’s Personal Health Information Protection Act. ICES is a prescribed entity under section 45 of Ontario’s Personal Health Information Protection Act, which authorizes ICES to collect personal health information, without consent, for the purpose of health system evaluation and improvement.
Results
The final study cohort was extracted as shown in Figure 1, while overall cohort characteristics at baseline are shown in Table 1. The demographic characteristics of propensity score matched sample are provided in the supplement (Supplementary Table S2).
Out of 6,175,114 individuals, there were 4,274,037 (69.21%) without a standard lab-based PCR test for COVID-19. Out of those who were tested for COVID-19 infection, there were 1,755,519 (28.43%) with negative results while 145,558 (2.36%) tested positive for COVID-19 during the study period. The mean age of the overall sample was 42.9 years with a standard deviation of 15.5. The number of females in the overall sample was higher (n = 3,189,737; 51.65%), as compared to males (n = 2,985,377; 48.35%). A higher proportion of individuals (n = 5,511,050; 89.25%) did not utilize health services for anxiety or depression-related issues as compared to those who contacted health services (n = 664,064; 10.75%) for these medical reasons. Out of those who contacted health services, the highest proportion of the individuals was 18–45 years of age (57.80%), followed by those between 46 and 65 years (32.62%) and 66–104 years (9.58%), respectively. In terms of sociodemographic status, the sample included 29.19% vs. 16.21% (Q1 vs. Q5) in dependency quintile, 24.61% vs. 17.29% (Q1 vs. Q5) in deprivation quintile, 22.56% vs. 22.63% (Q1 vs. Q5) in instability quintile, and 14.91% vs. 29.03% in ethnic diversity quintile (Table 1).
As seen in Figure 2, the time until the event (i.e., mental health service utilization), in months, was fastest in the group with negative COVID-19 test results, followed by the group with positive COVID-19 test results, both noticeably trailing those who never had standard lab-based PCR testing done for COVID-19 (Figure 2). For the negative group, the survival curve showed a gradual decrease in survival over time, whereas, for those without the COVID-19 PCR test, i.e., the control group, the curve followed the same pattern throughout. In the multivariable Cox proportional hazards model, COVID-19 testing, age, female sex, and chronic conditions were significantly associated with anxiety and depression related health service utilization (Table 2), with COVID-19 testing being the most important variable.
FIGURE 2
TABLE 2
| Variables | Crude HR (95% CI) | p-value | Adjusted HR (95% CI) | p-value |
|---|---|---|---|---|
| COVID-19 test status | ||||
| Covid positive | 3.83 (3.76–3.91) | <0.0001 | 4.14 (4.39–4.44) | <0.0001 |
| Covid negative | 4.21 (4.18–4.24) | <0.0001 | 4.41 (4.06–4.22) | <0.0001 |
| No PCR test | Ref | | Ref | |
| Age groups | ||||
| 18–45 years | 1.21 (1.20–1.22) | <0.0001 | 1.38 (1.36–1.39) | <0.0001 |
| 46–65 years | 1.14 (1.13–1.15) | <0.0001 | 1.31 (1.30–1.33) | <0.0001 |
| 65–104 years | Ref | | Ref | |
| Sex | ||||
| Female | 1.05 (1.05–1.06) | <0.0001 | 1.05 (1.04–1.05) | <0.0001 |
| Male | Ref | Ref | Ref | Ref |
| Rural | ||||
| Yes | 1.00 (0.99–1.00) | <0.0001 | 0.99 (0.98–1.00) | 0.2682 |
| No | Ref | | Ref | |
| Asthma | ||||
| Yes | 1.07 (1.05–1.10) | <0.0001 | 1.00 (0.98–1.02) | 0.7479 |
| No | Ref | | Ref | |
| COPD | ||||
| Yes | 0.94 (0.92–0.96) | <0.0001 | 0.91 (0.89–0.93) | <0.0001 |
| No | Ref | | Ref | |
| Congestive heart failure | ||||
| Yes | 0.69 (0.65–0.70) | <0.0001 | 0.43 (0.41–0.45) | <0.0001 |
| No | Ref | | Ref | |
| Hypertension | ||||
| Yes | 0.97 (0.96–0.98) | <0.0001 | 0.99 (0.98–1.00) | 0.2763 |
| No | Ref | | Ref | |
| Diabetes | ||||
| Yes | 0.95 (0.94–0.97) | <0.0001 | 0.92 (0.91–0.94) | <0.0001 |
| No | Ref | | Ref | |
| Ontario marginalization index Dependency quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 0.99 (0.98–1.00) | 0.1920 | 1.00 (0.99–1.00) | 0.9588 |
| 3 | 0.99 (0.98–1.00) | 0.3700 | 1.00 (0.99–1.00) | 0.7363 |
| 4 | 0.98 (0.98–0.99) | 0.0023 | 0.99 (0.99–1.00) | 0.7545 |
| 5 | 0.97 (0.97–0.98) | <0.0001 | 0.98 (0.97–0.99) | 0.0002 |
| Deprivation quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 0.99 (0.98–1.00) | 0.0761 | 0.99 (0.99–1.00) | 0.8050 |
| 3 | 0.99 (0.98–1.00) | 0.7505 | 1.00 (0.99–1.01) | 0.2907 |
| 4 | 0.98 (0.97–0.99) | 0.0002 | 1.00 (0.99–1.00) | 0.9674 |
| 5 | 0.97 (0.97–0.98) | <0.0001 | 1.00 (0.99–1.00) | 0.9923 |
| Instability quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 1.00 (0.99–1.00) | 0.9996 | 0.99 (0.98–0.99) | 0.0285 |
| 3 | 1.00 (0.99–1.00) | 0.7777 | 0.99 (0.98–1.00) | 0.0379 |
| 4 | 0.99 (0.99–1.00) | 0.7811 | 0.97 (0.97–0.98) | <0.0001 |
| 5 | 0.98 (0.97–0.99) | <0.0001 | 0.96 (0.955–0.97) | <0.0001 |
| Ethnic diversity quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 1.01 (1.00–1.02) | 0.0060 | 1.00 (0.99–1.00) | 0.9285 |
| 3 | 1.01 (1.00–1.01) | 0.0151 | 0.99 (0.98–1.00) | 0.1295 |
| 4 | 1.00 (0.99–1.01) | 0.5133 | 0.98 (0.97–0.99) | 0.0095 |
| 5 | 0.978 (0.97–0.98) | <0.0001 | 0.97 (0.96–0.98) | <0.0001 |
Multivariable Cox regression analysis of factors associated with anxiety/depression related health services utilization among Ontario population for overall Sample (Ontario, Canada, 2021).
In the adjusted Cox proportional hazards model for the overall sample, the individuals who had PCR tests had a higher health services utilization for anxiety and depression during the pandemic, as compared to those who had no PCR testing done for COVID-19. The hazard ratios were aHR, 4.14; (95%) CI, 4.39–4.44 and aHR, 4.41; CI, 4.06–4.22; for COVID-19 test positive and COVID-19 test negative, respectively (Table 2). A higher health services utilization for anxiety and depression was observed among females (aHR, 1.05; CI, 1.04–1.05), individuals in the 18–45 years age group (aHR 1.38; CI, 1.36–1.39) as well as individuals in the 46–65 years age group (aHR, 1.31; CI, 1.30–1.33).
In the Cox proportional hazards model in the propensity score matched sample, the individuals who tested positive for COVID-19 (aHR, 6.37; CI, 6.25–6.50) had the higher health services utilization for anxiety and depression during the pandemic, followed by those who tested negative for COVID-19 (aHR, 5.91; CI, 5.87–5.95) (Table 3). Also, a higher health services utilization was noted among young individuals, females, and those residing in rural areas. Moreover, higher services utilization was observed for socioeconomic status, with material deprivation quintile significantly associated with anxiety and depression. For the neighborhood-level income quintile, similar risks of anxiety and depression were found among those in Q3-Q5 of the material deprivation quintile (Table 3).
TABLE 3
| Variables | Crude HR (95% CI) | p-value | Adjusted HR (95% CI) | p-value |
|---|---|---|---|---|
| COVID-19 test status | ||||
| Covid positive | 6.34 (6.21–6.46) | <0.0001 | 6.37 (6.25–6.50) | <0.0001 |
| Covid negative | 5.88 (5.84–5.92) | <0.0001 | 5.91 (5.87–5.95) | <0.0001 |
| No PCR test | Ref | | Ref | |
| Age groups | ||||
| 18–45 years | 1.04 (1.03–1.05) | <0.0001 | 1.000 (0.977–0.996) | <0.0001 |
| 46–65 years | 1.00 (0.991–1.01) | 0.7066 | 0.980 (0.971–0.989) | <0.0001 |
| 65–104 years | Ref | | Ref | |
| Sex | ||||
| Female | 1.00 (0.99–1.00) | 0.0080 | 1.039 (1.034–1.044) | <0.0001 |
| Male | Ref | Ref | Ref | Ref |
| Rural | ||||
| Yes | 1.02 (1.02–1.03) | <0.0001 | 1.02 (1.01–1.03) | <0.0001 |
| No | Ref | | Ref | |
| Asthma | ||||
| Yes | 1.04 (1.02–1.06) | <0.0001 | 0.94 (0.92–0.96) | <0.0001 |
| No | Ref | | Ref | |
| COPD | ||||
| Yes | 1.00 (0.98–1.03) | 0.3877 | 0.99 (0.97–1.01) | 0.4776 |
| No | Ref | | Ref | |
| Congestive heart failure | ||||
| Yes | 1.21 (1.16–1.27) | <0.0001 | 1.09 (1.05–1.14) | <0.0001 |
| No | Ref | | Ref | |
| Hypertension | ||||
| Yes | 1.01 (1.00–1.02) | 0.0014 | 1.00 (0.99–1.01) | 0.2208 |
| No | Ref | | Ref | |
| Diabetes | ||||
| Yes | 1.00 (0.98–1.01) | 0.6820 | 0.99 (0.97–1.00) | 0.4915 |
| No | Ref | | Ref | |
| Ontario marginalization index | ||||
| Dependency quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 0.99 (0.99–1.00) | 0.6058 | 0.99 (0.99–1.00) | 0.3154 |
| 3 | 1.00 (0.99–1.00) | 0.7809 | 0.99 (0.99–1.00) | 0.5263 |
| 4 | 0.99 (0.99–1.00) | 0.6277 | 1.00 (0.99–1.00) | 0.8722 |
| 5 | 1.00 (0.99–1.00) | 0.9174 | 0.99 (0.99–1.00) | 0.8452 |
| Deprivation quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 0.99 (0.98–1.00) | 0.3369 | 1.00 (0.99–1.01) | 0.1062 |
| 3 | 1.00 (0.99–1.00) | 0.6515 | 1.01 (1.00–1.02) | <0.0001 |
| 4 | 0.99 (0.98–1.00) | 0.0557 | 1.02 (1.01–1.03) | <0.0001 |
| 5 | 0.98 (0.97–0.99) | <0.0001 | 1.03 (1.02–1.04) | <0.0001 |
| Instability quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 1.00 (0.99–1.01) | 0.5453 | 0.98 (0.98–0.99) | 0.0044 |
| 3 | 1.00 (0.99–1.01) | 0.4752 | 0.98 (0.97–0.99) | 0.0001 |
| 4 | 1.00 (0.99–1.00) | 0.7743 | 0.97 (0.96–0.97) | <0.0001 |
| 5 | 0.98 (0.97–0.99) | <0.0001 | 0.95 (0.94–0.95) | <0.0001 |
| Ethnic diversity quintile | ||||
| 1 | Ref | | Ref | |
| 2 | 1.00 (0.99–1.00) | 0.9959 | 0.99 (0.98–1.00) | 0.1657 |
| 3 | 0.99 (0.98–0.99) | 0.0127 | 0.97 (0.96–0.98) | <0.0001 |
| 4 | 0.98 (0.97–0.98) | <0.0001 | 0.97 (0.96–0.98) | <0.0001 |
| 5 | 0.96 (0.96–0.97) | <0.0001 | 0.97 (0.96–0.98) | <0.0001 |
Multivariable Cox regression analysis of factors associated with anxiety/depression related health services utilization among Ontario population after propensity scores computation (Ontario, Canada, 2021).
In the Cox proportional hazards model stratified by age groups, the health services utilization for anxiety and depression was highest among younger adults (aHR, 5.56; CI, 5.51–5.61), as compared to the middle-age group (aHR, 4.05; CI, 4.00–4.09) and older adults (aHR, 1.59; CI, 1.55–1.62), who tested negative with COVID-19 (Supplementary Table S1). Females were at higher risk as compared to males across all age groups (Table 3). Rural residence was only statistically significant among younger adults (aHR, 1.02; CI, 1.01–1.03) as compared to other age groups, though the effect is small. Among chronic conditions, asthma and hypertension were observed to be associated with anxiety and depression-related visits and hospitalization, across young and middle age groups (Supplementary Table S1).
Discussion
Main findings
This study examined differences in first recorded mental health service utilization for anxiety and depression between individuals who did and did not undergo COVID-19 PCR testing in a large population-based cohort during the pandemic. We observed higher utilization among individuals who underwent testing, reflecting underlying differences between tested and untested populations rather than a direct effect of testing itself, with adjusted hazard ratios of 6.37 (95% CI 6.25–6.50) for PCR-positive individuals and 5.91 (95% CI 5.87–5.95) for PCR-negative individuals. Higher utilization was particularly evident among younger individuals, females, and those residing in rural areas, with additional socioeconomic gradients indicating greater use among more deprived populations. The magnitude and similarity of effect estimates across PCR-positive and PCR-negative groups suggest that testing functions primarily as a proxy for underlying vulnerability and healthcare-seeking behavior rather than a causal determinant of mental health outcomes.
Comparison with existing literature
The observed patterns of mental health–related healthcare utilization may be consistent with broader reports from the pandemic; however, these findings cannot distinguish between true increases in mental health burden and differences in healthcare-seeking behavior. The increased mental health–related healthcare utilization observed in this study is consistent with evidence on the psychological impact of the COVID-19 pandemic, driven by factors such as uncertainty, social isolation, and economic disruption [21, 22]. Higher risks among younger individuals and females align with prior studies reporting disproportionate mental health impacts in these groups [28–32].
Social determinants of health play a key role in shaping these outcomes. Individuals in lower-income groups experienced greater mental health challenges, consistent with evidence linking financial stress, job loss, and poverty to anxiety and depression [33–35]. Similarly, we observed higher risks among rural residents, which may reflect barriers to care, social isolation, and limited access to mental health services [34, 36].
Chronic conditions, including asthma, hypertension, and congestive heart failure, were also associated with increased mental health service utilization. These conditions, identified as risk factors for severe COVID-19 illness, may contribute to heightened psychological distress, further exacerbated by disruptions in routine care during the pandemic [37, 38].
Overall, these findings highlight the disproportionate burden of mental health service utilization among socioeconomically and clinically vulnerable populations during the pandemic [37, 39].
What does this study add to existing knowledge?
This study characterizes patterns of mental health–related healthcare utilization among individuals who engaged with COVID-19 testing during the pandemic rather than establishing a causal relationship between testing and mental health outcomes. The findings indicate higher utilization among tested individuals; however, the magnitude of these associations is unusually large for observational mental health research and is unlikely to represent causal effects. Instead, these patterns are most plausibly explained by substantial residual confounding, including underlying differences in healthcare-seeking behavior and psychosocial vulnerability between individuals who seek testing and those who do not. The similarity of effect estimates across PCR-positive and PCR-negative groups further supports the interpretation that testing functions primarily as a proxy for underlying vulnerability rather than a direct determinant of mental health outcomes. Accordingly, these results should be interpreted as descriptive of healthcare utilization patterns during the pandemic rather than evidence of a direct psychological impact of COVID-19 testing.
Individuals undergoing testing may differ systematically in ways that predispose them to higher mental health service use, including greater baseline health anxiety, increased psychological vulnerability, and a higher propensity to engage with healthcare services [40]. These factors may contribute to residual confounding, even after adjustment for measured variables. Surveillance bias is also a likely contributor, as individuals interacting with the healthcare system through testing may be more likely to have mental health conditions identified and recorded. The similar effect estimates observed in both PCR-positive and PCR-negative groups further support the interpretation that testing itself is not the primary driver, but rather a marker of shared underlying characteristics among those who access testing [41]. Reverse causality must also be considered, as individuals with pre-existing anxiety or heightened concern about COVID-19 may have been more likely to seek testing.
Taken together, these findings suggest that COVID-19 testing serves as an indicator of a population with elevated underlying risk of mental health service utilization, rather than a direct causal determinant of mental health outcomes.
What are the key implications for public health practice or policy?
The observed disparities in mental health–related healthcare utilization among socioeconomically disadvantaged populations point to important inequities in patterns of healthcare engagement during the pandemic. However, these findings should not be interpreted as evidence of greater underlying mental health need, as they may reflect differences in healthcare access, utilization, and help-seeking behavior rather than true differences in burden. Instead, the results are best understood as identifying populations with higher levels of interaction with the healthcare system during this period.
While causal inferences cannot be drawn, these patterns may still be useful for public health planning by highlighting groups that were more frequently represented in mental health service use during the pandemic. Monitoring trends in healthcare utilization can support resource planning and help ensure that services are responsive to populations engaging with care. However, caution is required in translating these findings into policy, as observed differences may not correspond to underlying need.
Future research should prioritize disentangling the roles of psychosocial factors, healthcare-seeking behavior, and structural inequalities in shaping observed patterns of mental health service utilization [42]. The findings highlight the importance of interpreting administrative healthcare data within the broader context of social and behavioral responses to the pandemic. Strengthening mental health systems and incorporating mental health considerations into pandemic preparedness remain important goals, but should be informed by evidence that distinguishes true need from patterns of healthcare engagement.
Strengths and limitations
To our knowledge, this is one of the largest population-based studies examining the association between COVID-19 testing and mental health–related healthcare utilization in Ontario. The use of linked administrative health data within a universal healthcare system allowed for comprehensive capture of physician visits, emergency department encounters, and psychiatric hospitalizations, enhancing the generalizability of the findings. Additionally, the application of propensity score matching improved comparability between tested and untested individuals by balancing measured sociodemographic and clinical characteristics.
However, several important limitations should be considered when interpreting these findings. The outcome reflects the first recorded use of mental health–related services during the study period rather than incident diagnoses of anxiety or depression. Accordingly, the results represent patterns of healthcare utilization, which may include both new-onset and pre-existing but previously unrecorded conditions, rather than true incidence of mental illness.
Moreover, despite the use of propensity score matching, residual confounding remains likely, particularly due to unmeasured psychosocial factors. Variables such as employment status, income loss, social isolation, media exposure, baseline anxiety traits, and individual propensity for healthcare-seeking were not available in administrative data. These factors are strongly associated with both the likelihood of undergoing COVID-19 testing and the use of mental health services. Individuals experiencing greater psychological distress or socioeconomic disruption during the pandemic may have been both more likely to seek testing and more likely to access care. As a result, the observed associations may be inflated, reflecting underlying vulnerability rather than the effect of testing itself.
The similar magnitude of associations in both PCR-positive and PCR-negative groups further supports this interpretation, suggesting that testing is unlikely to be a direct causal driver of increased mental health service utilization. Instead, testing may act as a marker of heightened vulnerability, capturing individuals with greater health anxiety, perceived risk, or engagement with the healthcare system. Effect sizes of this magnitude are unlikely to be explained by measured confounders alone and likely reflect selection processes and unmeasured differences between tested and untested populations.
The administrative data captures only individuals who interact with the healthcare system, potentially excluding those with untreated or unrecognized mental health conditions. This may introduce detection (surveillance) bias, as individuals who seek testing are also more likely to have mental health concerns identified and recorded.
Additional limitations include the exclusion of individuals with missing demographic data, limited information on certain populations (e.g., immigrants and Indigenous communities), and the inability to assess severity or chronicity of mental health conditions. Finally, the study period did not extend into later phases of the pandemic, which may limit the assessment of longer-term mental health impacts.
Conclusion
Our study identifies higher mental health–related healthcare utilization among individuals who underwent COVID-19 PCR testing for anxiety and depression. These patterns reflect differences between tested and untested populations and should not be interpreted as evidence of a direct psychological impact of testing. The similar magnitude of associations observed in PCR-positive and PCR-negative individuals strongly suggests that testing functions primarily as a proxy for underlying vulnerability and healthcare-seeking behavior rather than a causal determinant of mental health outcomes. These findings pertain to patterns of healthcare utilization rather than new-onset mental illness. As such, the results are best understood as descriptive of healthcare engagement during the pandemic, highlighting the importance of interpreting administrative data within the context of underlying differences in access, utilization, and psychosocial factors. Incorporating mental health considerations into pandemic preparedness remains important but should be guided by evidence that distinguishes true need from patterns of healthcare use.
Statements
Data availability statement
The data set from this study is held securely in coded form at ICES. While data sharing agreements prohibit ICES from making the data set publicly available, access may be granted to those who meet pre-specified criteria for confidential access, available at https://www.ices.on.ca/DAS. The full data set creation plan and underlying analytic code are available from the authors upon request, understanding that the computer programs may rely upon coding templates or macros that are unique to ICES and are therefore either inaccessible or may require modification.
Ethics statement
The studies involving humans were approved by University of Waterloo research ethics board (protocol reference no. 43910). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants’ legal guardians/next of kin because Data use was authorized under section 45 of Ontario’s Personal Health Information Protection Act. ICES is a prescribed entity under section 45 of Ontario’s Personal Health Information Protection Act, which authorizes ICES to collect personal health information, without consent, for the purpose of health system evaluation and improvement.
Author contributions
KS was responsible for framework conceptualization, data acquisition, analysis, and drafting as well as revising the manuscript. JV provided statistical support, contributing to the selection and application of methodologies and interpretation of results. VG and JD offered critical feedback and guidance throughout the research process. They contributed to the interpretation of findings and revisions to enhance the manuscript’s clarity and impact. ZB conceptualized the study, provided critical feedback, supervision, and guidance throughout the research process and contributed to the interpretation of findings and revisions. All authors reviewed and approved the final manuscript. KS accepts full responsibility for the finished work and the conduct of the study, had access to the data, and the decision to publish.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This study was supported by the Ontario Health Data Platform (OHDP), a Province of Ontario initiative to support Ontario’s ongoing response to COVID-19 and its related impacts. The opinions, results and conclusions reported in this paper are those of the authors and are independent from the funding sources. No endorsement by the OHDP, its partners, or the Province of Ontario is intended or should be inferred. This study was also supported by ICES, which is funded by an annual grant from the Ontario Ministry of Health (MOH) and the Ministry of Long-Term Care (MLTC). This document used data adapted from the Statistics Canada Postal CodeOM Conversion File, which is based on data licensed from Canada Post Corporation, and/or data adapted from the Ontario Ministry of Health Postal Code Conversion File, which contains data copied under license from ©Canada Post Corporation and Statistics Canada. Parts of this material are based on data and/or information compiled and provided by Canadian Institute of Health Information (CIHI) and the Ontario Ministry of Health. The analyses, conclusions, opinions, and statements expressed herein are solely those of the authors and do not reflect those of the funding or data sources; no endorsement is intended or should be inferred. We thank the Toronto Community Health Profiles Partnership for providing access to the Ontario Marginalization Index. This research was undertaken, in part, thanks to funding from the Canada Research Chairs Program, Grant number: CRC-2022-00190. Dr. Butt holds the Canada Research Chair in Interdisciplinary Research for Pandemic Preparedness.
Acknowledgments
We would like to thank ICES for their support in data access, procurement, and data management.
Conflict of interest
The author(s) declared that they do not have any conflicts of interest.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.ssph-journal.org/articles/10.3389/ijph.2026.1609134/full#supplementary-material
Abbreviations
World Health Organization (WHO), post-traumatic stress disorder (PTSD), Institute for Clinical Evaluative Sciences (ICES), Registered Persons Database (RPDB), the Ontario Mental Health Reporting System (OMHRS), Discharge Abstract Database (DAD), emergency departments (EDs), National Ambulatory Care Reporting System (NACRS), Ontario Marginalization Index (ON-Marg).
References
1.
AhrendtDCabritaJClericiEHurleyJLeončikasTMascheriniMet alLiving, working and COVID-19 | eurofound (2020). Available online at: https://www.eurofound.europa.eu/publications/report/2020/living-working-and-covid-19 (Accessed June 10, 2023).
2.
SavitskyBFindlingYEreliAHendelT. Anxiety and coping strategies among nursing students during the covid-19 pandemic. Nurse Educ Pract (2020) 46:102809. 10.1016/J.NEPR.2020.102809
3.
GritsenkoVSkugarevskyOKonstantinovVKhamenkaNMarinovaTReznikAet alCOVID 19 fear, stress, anxiety, and substance use among Russian and Belarusian university students. Int J Ment Health Addict (2021) 19:1–7. 10.1007/S11469-020-00330-Z
4.
MamunMAChandrimaRMGriffithsMD. Mother and son suicide pact due to COVID-19-Related online learning issues in Bangladesh: an unusual case report. Int J Ment Health Addict (2020) 1:1–4. 10.1007/S11469-020-00362-5
5.
Director-general WHO. COVID-19 pandemic triggers 25 % increase in prevalence of anxiety and depression worldwide. World Health Organ (2022) 25–8. Available online at: https://www.who.int/news/item/02-03-2022-covid-19-pandemic-triggers-25-increase-in-prevalence-of-anxiety-and-depression-worldwide (Accessed May 18, 2022).
6.
SantomauroDFMantilla HerreraAMShadidJZhengPAshbaughCPigottDMet alGlobal prevalence and burden of depressive and anxiety disorders in 204 countries and territories in 2020 due to the COVID-19 pandemic. The Lancet (2021) 398:1700–12. 10.1016/S0140-6736(21)02143-7/ATTACHMENT/927FDFEF-CCD4-4655-AACF-4E7D54DFECF5/MMC1.PDF
7.
YaoHChenJHXuYF. Patients with mental health disorders in the COVID-19 epidemic. Lancet Psychiatry (2020) 7:e21. 10.1016/S2215-0366(20)30090-0
8.
An echo pandemic of mental health issues? Not if we can help it. CMHA Natl (2022) Available online at: https://cmha.ca/an-echo-pandemic-of-mental-health-issues-not-if-we-can-help-it/ (Accessed May 21, 2022).
9.
Canadian Mental Health Association. Policy brief Covid-19 and mental health: heading off an echo pandemic introduction. The Impact Of Covid-19 On Ment Health (2020) 1–9.
10.
Castaldelli-MaiaJMMarzialiMELuZMartinsSS. Investigating the effect of national government physical distancing measures on depression and anxiety during the COVID-19 pandemic through meta-analysis and meta-regression. Psychol Med (2021) 51:881–93. 10.1017/S0033291721000933
11.
LorenzLPerkoniggAMaerckerA. A socio-interpersonal approach to adjustment disorder: the example of involuntary job loss. Eur J Psychotraumatol (2018) 9:1425576. 10.1080/20008198.2018.1425576
12.
PengidSPeltzerKde Moura VillelaEFFodjoJNSSiauCSChenWSet alUsing Andersen’s model of health care utilization to assess factors associated with COVID-19 testing among adults in nine low-and middle-income countries: an online survey. BMC Health Serv Res (2022) 22:265. 10.1186/S12913-022-07661-8
13.
KimTJvon dem KnesebeckO. von dem Knesebeck O. Perceived job insecurity, unemployment and depressive symptoms: a systematic review and meta-analysis of prospective observational studies. Int Arch Occup Environ Health (2016) 89:561–73. 10.1007/S00420-015-1107-1
14.
LipsitchMSwerdlowDLFinelliL. Defining the epidemiology of Covid-19 — studies needed. New Engl J Med (2020) 382:1194–6. 10.1056/nejmp2002125
15.
GuptaRKMarksMSamuelsTHALuintelARamplingTChowdhuryHet alSystematic evaluation and external validation of 22 prognostic models among hospitalised adults with COVID-19: an observational cohort study. Eur Respir J (2020) 56:2003498. 10.1183/13993003.03498-2020
16.
Guttmann-DuckeCKlingerSZiescheROtzelbergerBIdzkoMPonocnyAet alPersonality traits and mental distress after COVID-19 testing. Prospective long-term analysis in a Viennese cohort. Front Psychiatry (2023) 14:1129794. 10.3389/FPSYT.2023.1129794
17.
LeeSACrunkEA. Fear and psychopathology during the COVID-19 crisis: neuroticism, hypochondriasis, reassurance-seeking, and coronaphobia as fear factors. Omega (2020) 85:483–96. 10.1177/0030222820949350
18.
LeeSAJobeMCMathisAAGibbonsJA. Incremental validity of coronaphobia: coronavirus anxiety explains depression, generalized anxiety, and death anxiety. J Anxiety Disord (2020) 74:102268. 10.1016/J.JANXDIS.2020.102268
19.
CevallosMEggerM. STROBE (STrengthening the reporting of observational studies in epidemiology). In: Guidelines for reporting health research: a user’s manual (2014). p. 169–79. 10.1002/9781118715598.ch17
20.
Canadian Institute for Health Information. Mental illness hospitalization (2023). Available online at: https://www.cihi.ca/en/indicators/mental-illness-hospitalization (Accessed February 22, 2023).
21.
ChiuMGuttmannAKurdyakP. Mental health and addictions system performance in Ontario: an updated scorecard, 2009-2017. Healthc Q (2020) 23:7–11. 10.12927/HCQ.2020.26340
22.
MHASEF Research Team. The mental health of children and youth in Ontario: a baseline scorecard (2015). Technical Appendix.
23.
KurdyakPGandhiSHolderLRashidMSaundersNChiuMet alIncidence of access to ambulatory mental health care prior to a psychiatric emergency department visit among adults in Ontario, 2010-2018. JAMA Netw Open (2021) 4:e215902. 10.1001/JAMANETWORKOPEN.2021.5902
24.
DoktorchikCPattenSEastwoodCPengMChenGBeckCAet alValidation of a case definition for depression in administrative data against primary chart data as a reference standard. BMC Psychiatry (2019) 19:1–8. 10.1186/S12888-018-1990-6/TABLES/3
25.
MyranDTCantorNRhodesEPuglieseMHenselJTaljaardMet alPhysician health care visits for mental health and substance use during the COVID-19 pandemic in Ontario, Canada. JAMA Netw Open (2022) 5:e2143160. 10.1001/jamanetworkopen.2021.43160
26.
GeELiYWuSCandidoEWeiX. Association of pre-existing comorbidities with mortality and disease severity among 167,500 individuals with COVID-19 in Canada: a population-based cohort study. PLoS One (2021) 16:e0258154. 10.1371/journal.pone.0258154
27.
Public Health Ontario. Ontario marginalization index (ON-Marg). Public Health Ontario (2022). Available online at: https://www.publichealthontario.ca/en/data-and-analysis/health-equity/ontario-marginalization-index (Accessed November 27, 2022).
28.
CénatJMDalexisRD. The complex trauma spectrum during the COVID-19 pandemic: a threat for children and adolescents’ physical and mental health. Psychiatry Res (2020) 293:113473. 10.1016/J.PSYCHRES.2020.113473
29.
GolbersteinEWenHMillerBF. Coronavirus disease 2019 (COVID-19) and mental health for children and adolescents. JAMA Pediatr (2020) 174:819–20. 10.1001/JAMAPEDIATRICS.2020.1456
30.
JiaoWYWangLNLiuJFangSFJiaoFYPettoello-MantovaniMet alBehavioral and emotional disorders in children during the COVID-19 epidemic. J Pediatr (2020) 221:264–6.e1. 10.1016/j.jpeds.2020.03.013
31.
BhatiaR. Editorial: effects of the COVID-19 pandemic on child and adolescent mental health. Curr Opin Psychiatry (2020) 33:568–70. 10.1097/YCO.0000000000000651
32.
Watkins-MartinKOrriMPennestriMHCastellanos-RyanNLaroseSGouinJPet alDepression and anxiety symptoms in young adults before and during the COVID-19 pandemic: evidence from a Canadian population-based cohort. Ann Gen Psychiatry (2021) 20:42. 10.1186/S12991-021-00362-2
33.
CordesJCastroMC. Spatial analysis of COVID-19 clusters and contextual factors in New York City. Spat Spatiotemporal Epidemiol (2020) 34:100355. 10.1016/j.sste.2020.100355
34.
StokesEKZambranoLDAndersonKNMarderEPRazKMEl Burai FelixSet alCoronavirus disease 2019 case surveillance - united States, January 22-May 30, 2020. MMWR Morb Mortal Wkly Rep (2020) 69:759–65. 10.15585/MMWR.MM6924E2
35.
The global economic outlook during the COVID-19 pandemic: a changed world (2026) Available online at: https://www.worldbank.org/en/news/feature/2020/06/08/the-global-economic-outlook-during-the-covid-19-pandemic-a-changed-world (Accessed April 6, 2026).
36.
MillettGAJonesATBenkeserDBaralSMercerLBeyrerCet alAssessing differential impacts of COVID-19 on Black communities. Ann Epidemiol (2020) 47:37–44. 10.1016/j.annepidem.2020.05.003
37.
LundCBreenAFlisherAJKakumaRCorrigallJJoskaJAet alPoverty and common mental disorders in low and middle income countries: a systematic review. Soc Sci Med (2010) 71:517–28. 10.1016/j.socscimed.2010.04.027
38.
MendenhallEKohrtBANorrisSANdeteiDPrabhakaranD. Non-communicable disease syndemics: poverty, depression, and diabetes among low-income populations. Lancet (2017) 389:951–63. 10.1016/S0140-6736(17)30402-6
39.
IemmiVBantjesJCoastEChannerKLeoneTMcDaidDet alSuicide and poverty in low-income and middle-income countries: a systematic review. Lancet Psychiatry (2016) 3:774–83. 10.1016/S2215-0366(16)30066-9
40.
BradyREBrazAN. Challenging interactions between patients with severe health anxiety and the healthcare system: a qualitative investigation. J Prim Care Community Health (2023) 14:21501319231214876. 10.1177/21501319231214876
41.
ShultzJMCooperJLBainganaFOquendoMAEspinelZAlthouseBMet alThe role of fear-related behaviors in the 2013-2016 West Africa ebola virus disease outbreak. Curr Psychiatry Rep (2016) 18:104. 10.1007/S11920-016-0741-Y
42.
PelletierLO’DonnellSMcRaeLGrenierJ. The burden of generalized anxiety disorder in Canada. Health Promotion Chronic Dis Prev Can (2017) 37:54–62. 10.24095/hpcdp.37.2.04
Summary
Keywords
anxiety, COVID-19, depression, health administrative data, mental health
Citation
Saqib K, Goel V, Dubin JA, VanderDoes J and Butt ZA (2026) COVID-19 testing and mental health service utilization in Ontario: a population-based cohort study. Int. J. Public Health 71:1609134. doi: 10.3389/ijph.2026.1609134
Received
27 September 2025
Revised
17 April 2026
Accepted
04 May 2026
Published
14 May 2026
Volume
71 - 2026
Edited by
Abanoub Riad, Masaryk University, Czechia
Reviewed by
Fernando Del Fiol, University of Sorocaba, Brazil
Pablo Lopez, CONICET Institute of Cognitive and Translational Neuroscience (INCYT), Argentina
Updates
Copyright
© 2026 Saqib, Goel, Dubin, VanderDoes and Butt.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Zahid Ahmad Butt, zahid.butt@uwaterloo.ca
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.