Geographic and demographic heterogeneity of SARS-CoV-2 diagnostic testing in Illinois, USA, March to December 2020

Case definition

This work defines cases as SARS-CoV-2 infections recorded in Illinois surveillance as a result of a positive diagnostic test, regardless of symptom status.

Datasets

County-level positive tests and total tests were obtained from the Illinois National Electronic Disease Surveillance System (I-NEDSS) database maintained by IDPH. Daily testing volume data included 12,746,960 total specimens and 1,131,284 positive specimens recorded from March 17, 2020, to December 31, 2020, stratified by age, county of test, and race/ethnicity. Moving averages of daily testing volume were calculated on a seven-day lagging window. Until October 14, 2020, only molecular tests (reverse transcriptase polymerase chain reaction [RT- PCR] tests) were reported in this dataset. On October 14, 2020, IDPH began reporting antigen tests in this dataset. Testing site locations were scraped from the IDPH website on April 23, June 15, and October 26, 2020.

Individual-level case data, including date of first positive specimen, patient’s home ZIP code, race, ethnicity, hospital admission status, and date of death, were obtained from I-NEDSS. Data were pulled on March 16, 2021, and included 1,098,549 cases reported to IDPH, 907,799 of which had specimen collection dates in 2020. Among all cases that had a date of death in I-NEDSS, 18,830 were designated as having died due to COVID-19 and were considered confirmed deaths. Cases classified as died from COVID-19 met at least one of the following criteria: presence of COVID-19 on the death certificate; death within 30 days of symptom onset/diagnosis or during hospitalization, unless the cause of death is clearly unrelated to COVID-19 (e.g. accident); never returned to baseline health after diagnosis; autopsy result consistent with COVID-19. Individuals with date of first positive specimen on or before December 31, 2020, whose deaths were confirmed in I-NEDSS after March 16, 2021, would not be included in this death tally.

In the case data, individuals were assigned to a region based on the county of symptom onset, and secondarily on listed ZIP of residence if county of onset was not available. To estimate under-reporting rates, a naïve (crude) CFR was calculated as cumulative deaths divided by cumulative cases. Counties were aggregated into COVID-19 Regions as defined by IDPH specifically for the COVID-19 response (9) and into super-regions as follows: COVID-19 Regions 1 and 2 into North Central super-region; 3 and 6 into Central super-region; 4 and 5 into Southern super-region; and 7-11 into Northeast (Figure 1A). County populations were obtained from the 2018 American Community Survey (ACS) (10).

Self-reported race or ethnicity were available for 657,219 (72.4%) individual cases reported to IDPH and were missing for the remainder. Cases with multiple races reported were categorized as “Other”. All individuals with Hispanic/Latino as ethnicity were categorized as Hispanic/Latino regardless of race(s). Individuals with “unknown” ethnicity or no reported ethnicity were considered non-Hispanic/Latino. For brevity, non-Hispanic Black and non-Hispanic White populations are referred to as Black and White, respectively.

In the testing dataset, race or ethnicity was recorded for 7,143,108 total specimens (56.0%) and 652,643 positive specimens (57.7%) using a single variable indicating either a non-Hispanic race or Hispanic-Latino ethnicity.

Measuring distance to nearest testing site

We computed the distance from the centroid of each census block group to its nearest testing site location on October 26, 2020, and used the estimated population of each census block group (2016 ACS via Safegraph) to create a cumulative distribution of this distance over a region’s population. Census block groups were assigned to COVID-19 region by whether a census block group’s centroid fell within the boundaries of a COVID-19 region. Distance to the nearest testing site by ZIP was measured from each ZIP’s centroid to the nearest testing site location listed on IDPH’s website as sites open to the public on April 23, June 15, and October 26, 2020. The IDPH list of testing sites is not comprehensive as some testing sites asked not to be listed, and data on the actual number of sites offering testing were not available.

Estimation of infection detection rate

To estimate the fraction of infections detected in a particular week (f_{inf det}), the expected infection fatality ratio (IFR) among cases with new positive specimens collected during that week was calculated based upon these cases’ age distribution using either (i) the exponential meta-regression performed by Levin et al. (11), which used first-wave data from multiple countries; or (ii) estimates from O’Driscoll et al. (12), which uses an ensemble model incorporating data from multiple countries to infer age-specific infection mortality rates. Infection fatality ratio is the fraction of all SARS-CoV-2 infections that result in death. The results of the meta-regression of Levin et al. (11), with associated uncertainties for each coefficient, are reproduced below:

Due to the fact that, at any given time, the age distribution of cases (i.e. detected infections) was not necessarily representative of the age distribution of all incident infections, only cases 61-70 years of age were included in this analysis. This age range was selected for its large number of cases with confirmed COVID-19 deaths (>12 every week of specimen collection except for the week of March 8^th, 2020, the first week in which a fatal case was documented in ages 61-70 years) while being less likely than older age groups (>70 years) to be associated with the less representative transmission and testing conditions in long-term care facilities. For our estimates using IFR from O’Driscoll et al. (12), IFR was uniformly sampled from 0.39 - 1.24% for ages 61-70 years, which was obtained by combining IFR of 0.46% (95% CI: 0.39 - 0.57%) for ages 60-64 years and 1.08% (95% CI: 0.92 - 1.24%) for ages 65-69 years.

First, a naive (crude) estimate of f_{inf det} was made by dividing the expected IFR by the reported CFR among that week’s cases. To account for a decreasing IFR due to improved clinical outcomes among the infected over the course of the pandemic, a second estimate was made by adjusting the expected IFR down by the relative decrease in a sigmoid curve fitted to the hospital fatality ratio (HFR) over time among people aged 61-70 years. HFR was calculated from I-NEDSS data as fraction of admitted cases that were later recorded as a death due to COVID-19. This sigmoid curve was fitted to weekly HFR with a non-linear least squares regression.

To account for unreported deaths, a third estimate was made in which the adjusted IFR was divided by the CFR, then multiplied by the estimated fraction of all COVID-19 deaths that were reported as COVID-19 deaths (f_{death det}) on the median date of death among that week’s cases. To estimate f_{death det}, we compared observed counts of COVID-19 deaths to excess deaths in select-cause mortality data (Figure S1). Select-cause mortality data provided by the National Center for Health Statistics (NCHS), including respiratory diseases and circulatory diseases among others, showed the expected weekly count of deaths by a selection of comorbid conditions of COVID-19 alongside the reported counts of deaths by these causes that occurred in the state in 2020 (13). Excess select-cause deaths are calculated as the difference between the expected weekly select-cause death curve and the actual weekly select-cause death curve. Assuming that all excess select-cause deaths were attributable to COVID-19 and that the epidemic did not appreciably reduce deaths indirectly due to other causes in the list of select causes curated by NCHS, we calculated f_{death det} by dividing the observed number of COVID-19 deaths each week (from I-NEDSS) by the excess select-cause deaths in the same week. To account for uncertainty in our estimate of f_{death det}, we then sampled 1,000 realizations of excess deaths using a Skellam distribution, which models the difference between two Poisson random variables, and recalculated f_{death det} for each realization.

• Download figure
• Open in new tab

Figure S1.

Calculation of f_deathdet from Illinois COVID-19-attributed deaths and excess select cause deaths.

These three estimates were made from the week of March 8^th to the week of December 27^th, with 1,000 bootstrapped samples taken on a weekly basis from estimates of CFR for cases in a given week, expected IFR, the prediction band of the sigmoid curve fitted to HFR, and the estimates of f_{death det} to generate a range of estimates for f_{inf det}. For cases in a given week, estimates of f_{death det} were drawn from the week of the median death date of that week’s cases. All infection detection estimates, as well as HFR and f_{death det}, were conducted at the statewide level.

Spatio-temporal variation in testing and access to testing in Illinois

As of December 31, 2020, more than 900,000 SARS-CoV-2 cases and 18,000 COVID-19 deaths were recorded in Illinois (Figure 2) (14). The first wave of COVID-19 occurred in early May in the Northeast and Southern super-regions and in mid- to late-May in the Central and North-Central super-regions. COVID-19 Regions 1 and 7-11 experienced the bulk of the first-wave cases and deaths, and Regions 11 (city of Chicago) and 10 (suburban Cook County) recorded the highest peaks in daily detected cases during this time. By August 2020, daily detected cases in Regions 2, 3, 4, 5, and 6 had surpassed their peak numbers in May. By October 2020, all COVID-19 regions were experiencing a second wave of hospitalizations and deaths.

• Download figure
• Open in new tab

Figure 2.

Epidemic trajectory of COVID-19 cases in the 11 COVID-19 regions of Illinois in 2020.

Reflecting both population density and the regional differences in initial burden of COVID-19, the majority of testing sites were located in the Northeast super-region (64.7% of Illinois testing sites on October 26, 2020), and many were in Regions 10 and 11 (42.1% of Illinois testing sites on October 26, 2020) (Figure 1A). Although the number of diagnostic testing sites in the state has nearly quadrupled since April (Figure 1B), most new testing sites since June have been established in the Northeast super-region. Over 50% of individuals in Regions 3, 4, and 6 resided more than 10 miles from the nearest Illinois testing site (Figure 1B and 1C). This distance is not necessarily reflective of the distance any given individual in an area must travel or will travel to receive a test. Many test sites restricted testing to symptomatic individuals, close contacts, or in-network patients in terms of referrals or insurance plans, but testing criteria data were not sufficiently available or reliable to assess access to unrestricted testing. Moreover, these restrictions were subject to continuous change as the availability of resources at individual sites fluctuated. Conversely, individuals in border areas could seek testing across state lines.

Although testing was limited in all COVID-19 regions during the first wave, testing volume expanded 5- to 10-fold between early May and the end of December (Figure 3). Controlling for population size, the overall testing rate was highest in Regions 10 and 11 during the early outbreak in March to June 2020. Some regions (particularly Regions 1 and 6) that were slower to increase testing in the first wave outpaced other regions in testing intensity by late October due to prioritization of the deployment of mobile teams to areas of greatest impact (meat processing plants, low income housing areas, etc.) and establishment of community drive through testing sites where none previously existed. In November and December there was a concerted increase in testing in all Regions, with Region 3 achieving the highest testing rates in the state, before a decrease around Christmas.

• Download figure
• Open in new tab

Figure 3.

Daily SARS-CoV-2 diagnostic tests administered per 1,000 population in each COVID-19 region in Illinois in 2020. Shown: 7-day moving averages. Colors indicate super-region membership of each COVID-19 region, as indicated in Figure 1A.

In Region 6, overall testing intensity was dominated by the University of Illinois at Urbana-Champaign (UIUC) due to their efforts to conduct mass testing on their entire campus population (5) (Figure 3 Region 6 with and without Champaign County, Figure 4A). However, outside of Champaign County, the remainder of Region 6 contained some of the lowest per-capita testing in Illinois, reflecting the substantial portion of Region 6 residents who resided more than 10 miles from a testing site (Figure 1C). There was considerable county-level heterogeneity in testing intensity (Figure 4A) and positive tests per capita (Figure 4B) within Regions 1-6.

• Download figure
• Open in new tab

Figure 4.

(A) County-level average daily SARS-CoV-2 diagnostic tests per 1,000 population for a representative week ending October 26, 2020. Champaign County, in COVID-19 Region 6, stands out with the highest per capita testing rate in the state. Central and Southern counties have the lowest rates. (B) County-level average daily positive tests per 1,000 population for the week ending October 26, 2020.

Changing demographics of the tested population

Prior to mid-August, testing was most intensive in the elderly population, with those aged over 80 years receiving the most tests per capita, and intensity of testing increasing with age (Figure 5A). Routine testing in long-term care facilities contributed to higher testing intensity in the elderly population (15,16). In July, pilot testing at UIUC led testing in people aged 18-22 years to exceed testing by more than twice the rate in all other age groups, even the over 80-year-old group. Testing at other university campuses would also contribute to the increased testing rate in 18-22 year-old people, and the testing rate declined in late November following the Thanksgiving holiday and winter recess. Working-age adults may have been subject to routine testing at employers’ behest. Pediatric testing, including testing in older children, remained much lower than all other age groups. This could have been due to lower prevalence of SARS-CoV-2 infection in children because of lower susceptibility (17), low rates of test-seeking among SARS-CoV-2-infected children because they are less likely to be symptomatic, barriers to accessing testing because some providers did not test pediatric patients, or simply a lack of routine testing in children.

• Download figure
• Open in new tab

Figure 5.

SARS-CoV-2 diagnosic testing in Illinois stratified by age group. (A) Seven-day moving average of daily tests per 1,000 population, by age group. (B) Share of tests by age group.

The portion of tests that were conducted on young people steadily increased between July and September 2020 (Figure 5B). Because COVID-19 is more likely to be asymptomatic or mild in younger patients (18,19), this change in the tested population would be expected to lead to lower case fatality rates in the population as a whole even without any improvements in treating COVID-19.

We compared the per capita testing rates of non-Hispanic White, Black, and Hispanic/Latino populations in Illinois (Figure 6). Testing increased between March and October for all three groups. During the first wave, per capita testing was highest in Black and Hispanic/Latino populations for most age groups, reflecting their disproportionate share of COVID-19 burden (14). After the end of the first wave in late June, testing was consistently lowest in Hispanic/Latino populations, with only minimal expansion of testing among Hispanic/Latino elders. The testing rate in the Black population saw a sharp increase around the beginning of July but did not increase further between July and October. In contrast, the testing rate in the White population increased steadily. The impact of student testing at university campuses was visible in all three demographic groups as the step-increase from mid-August to late-November and was largest in the White population. Testing rates began to decline sharply in all groups by mid-December.

• Download figure
• Open in new tab

Figure 6.

Seven-day moving average of daily SARS-CoV-2 tests per 1,000 population, by race/ethnicity and age group in Illinois in 2020.

On a per capita basis, testing intensity was lowest in Hispanic-Latino populations and highest in Black populations, although the extent of this difference appeared to vary by COVID-19 Region (Figure S2). Since the pandemic has disproportionately affected Black and Hispanic/Latino communities in Illinois (14), the higher testing intensity in Black populations does not necessarily mean that the testing was sufficient to capture burden relative to White populations.

• Download figure
• Open in new tab

Figure S2.

Testing rate per capita by race/ethnicity and COVID-19 Region. Population denominators were drawn from ACS 2018. “Native” included both “American Indian or Alaskan Native” and “Native Hawaiian or Other Pacific Islander”. Asian and Native groups were excluded from Figures 6 and 8 due to small population denominators after stratifying by age. The small denominators are reflected in the large fluctuations in the Native timeseries in this figure.

Crude assessment of under-testing with case fatality ratio

We considered whether the naive case fatality ratio (CFR), defined as the number of COVID-19 deaths divided by number of detected cases, could assess differences in SARS-CoV-2 testing in each super-region when detailed data on testing rates and hospital admissions were unavailable. Assuming that COVID-19 death ascertainment rates were uniformly high, and there was little to no geographic variation in infection fatality rate for SARS-CoV-2, we expected areas and populations with higher testing rates to also have lower case fatality rates.

In Illinois, crude CFR decreased over the course of 2020 in all super-regions and age groups (Figure 7), concurrent with the scale-up of testing. CFR was highest in the Northeast and Southern super-regions in the first two months of the epidemic, although CFR in the Central super-region increased in May and June for adults aged 41-50 years and 61-70 years. From July onward, CFR was similar in all regions.

• Download figure
• Open in new tab

Figure 7.

COVID-19 case fatality ratios (CFR), the fraction of recorded cases with a COVID-19-attributed death, by age group and super-region, for cases detected in 2020. Error bars indicate the standard error in regional CFR.

Differences in CFR could be driven by heterogeneous access to testing and care as well as regional differences in the prevalence of comorbidities, standard of care, or hospital capacity. The higher CFR in the Southern super-region prior to July reflected the consistently lower testing rates in COVID-19 Regions 4 and 5 during the first wave. However, the Northeast’s CFR prior to July was substantially elevated over the CFRs in the Central and North-Central super-regions despite the Northeast super-region’s higher intensity of testing. This discrepancy suggested that despite its higher testing intensity during the first wave, the Northeast super-region was disproportionately under-tested relative to its share of the state’s COVID-19 burden. Alternatively, the higher CFR despite higher testing in the Northeast could have been driven by insufficient targeting of tests to the most affected populations. The convergence of regional CFRs in late 2020 suggested that earlier differences in the regional CFRs might not have been driven by differences in regional prevalence of comorbidities. In April-May 2020, the Northeast super-region experienced the greatest strain on hospital capacity (Figure S3), which could have contributed to lower quality of care. However, hospital capacity overall is also highest in the Northeast, and peak inpatient census did not exceed capacity (Table S1).

• Download figure
• Open in new tab

Figure S3.

Daily COVID-19 hospital census by COVID-19 Region (subplots) and super-region (colors) in 2020. Lines show 7-day rolling averages. Data for med/surg occupancy are not available prior to May 2020. Med/surg census covers all admitted patients who are not in the intensive care unit (ICU).

Most majority-Black and majority-Hispanic/Latino ZIP codes are in the Northeast and Southern super-regions (10), where the highest overall CFRs were observed during March-April 2020. When stratified by age and race (Figure 8), CFR increased with age, with disparities becoming less pronounced for older age groups and as the epidemic progressed. CFRs remained higher for Black and Hispanic populations for all but the over-80 year-old age group through the end of 2020. Higher fatality rates among Black and Hispanic/Latino cases could be due to a combination of under-testing leading to fewer detected cases and disparities in clinical outcomes resulting in more deaths. The under-testing of Black and Hispanic/Latino populations is reflected in the higher CFRs in these populations. While elevated prevalence of comorbidities such as diabetes and hypertension increased the underlying infection fatality rate in in Black and Hispanic/Latino populations, the enormous difference in CFR in younger age groups in Mar-Jun 2020 is unlikely to be explained by comorbidities alone. For example, if diabetes prevalence at age 45 were around 5% in non-Hispanic Whites and 11% in non-Hispanic Blacks (20), and diabetes increased the risk of severe outcomes by around 60% (21), the disparity in diabetes prevalence would increase the CFR in the Black population ages 41-50 by approximately 10%. Yet in this age group, the relative risk of death given a case was almost 300% higher in the Black population compared to White in Mar-Apr 2020 [2.94 (95% CI: 1.72-5.03)]. In Nov-Dec 2020, relative risk remained high at 3.31 (2.25-4.89).

• Download figure
• Open in new tab

Figure 8.

COVID-19 case fatality ratio by age group in non-Hispanic White, Hispanic/Latino, and non-Hispanic Black populations, for cases detected in 2020. Error bars indicate the standard error in group CFR.

In a sensitivity analysis, similar results were observed when cases with “unknown” ethnicity were removed altogether, instead of being allocated to non-Hispanic racial groups (Figure S4). As presented, the latter scenario may slightly underestimate CFR for older Black and White populations.

• Download figure
• Open in new tab

Figure S4.

Sensitivity of CFR to case race/ethnicity assignment. Alternate scenarios were run in which cases with “unknown” ethnicity were: 1) Assumed to be non-Hispanic and assigned to the recorded racial group (solid line, and identical to Figure 8), 2) Assigned to the “unknown” racial group (dashed line), or 3) Assumed to be Hispanic (dotted line). For Hispanic-Latino, scenarios #1 and #2 are identical. For Black and White, scenarios #2 and #3 are identical. Assumption of non-Hispanic ethnicity slightly, but significantly, decreases CFR for older non-Hispanic White and Black populations.

The majority of SARS-CoV-2 infections were never detected

Low per capita testing rates are not necessarily problematic if there is little SARS-CoV-2 circulation and testing is highly targeted. These conditions do not describe Illinois in 2020. To estimate the extent to which testing was able to identify all incident SARS-CoV-2 infections in Illinois as a whole, we compared the expected IFR among individuals aged 61-70 years, as estimated by Levin et al. (11) and O’Driscoll et al. (12), to the CFR among the same group (Figure 9A). We restricted the analysis to this age group because detection rates are likely highly heterogeneous across age groups, and the 61-70 age group has a sizable number of weekly cases and deaths. We generated a naive estimate (Figure 9D) as well as estimates accounting for both a non-stationary IFR due to improving clinical outcomes among the infected (Figure 9B, E) and under-reporting of deaths (Figure 9C, F). This methodology only provides estimates of the detection rate for all infections and does not account for any heterogeneity in the detection of mild symptomatic and asymptomatic infections versus severely symptomatic infections. Because O’Driscoll et al. (12) estimates a lower IFR for ages 61-70 than Levin et al. (11), the estimated fraction of infections detected using IFR from O’Driscoll et al. is slightly lower than the same estimate using IFR from Levin et al.

• Download figure
• Open in new tab

Figure 9.

Estimating fraction of infections detected. (A) Naive CFR among adults aged 61-70 years by week of specimen collection alongside expected IFR among this age distribution, as estimated by Levin et al. (11) or O’Driscoll et al. (12). Shaded regions show 95% confidence intervals. (B) HFR among admitted 61-70 year-olds by week of specimen collection (solid red) with standard error of proportions (shaded red), fitted sigmoid curve (solid black) and 95% prediction interval (dashed black). (C) Fraction of all deaths reported by date of death, based upon comparison of COVID-19 mortality and excess all-cause mortality. (D) Estimated fraction of SARS-CoV-2 infections among adults aged 61-70-years that were detected by surveillance: assuming that all deaths are reported and IFR is stationary; (E) assuming all COVID-19 deaths are reported and IFR is non-stationary due to improving clinical outcomes; (F) and assuming that excess deaths are unreported COVID-19 deaths and IFR is non-stationary. Dark blue estimates use IFR estimates from Levin et al. (11) and light blue from O’Driscoll et al. (12). Shaded regions in D-F are 95% credible intervals with 1000 bootstrapped samples.

In March and April 2020, excess deaths in Illinois greatly exceeded COVID-19-attributed deaths, driving down the estimates of infection detection rates in Figure 9F. Due to lack of data on cause of death, we made the simplifying assumption in Figure 9C and 10F that all excess select-cause deaths documented by NCHS were COVID-19 related. This assumption produced a floor on the estimated detection rate of SARS-CoV-2 infections because not all excess deaths would be directly related to COVID-19. The true fraction of SARS-CoV-2 infections detected in Illinois likely lay between the estimates in Figure 9E and 9F.

In the early epidemic, prior to mid-April, we estimated that less than 10% of SARS-CoV-2 infections among adults aged 61-70 years were detected and reported to IDPH (Figure 9D-F). This low level of detection in the early stages of the epidemic is consistent with other estimates of around 10% detection rate (3,22). Despite the 3- to 4-fold scale-up in testing volume over the summer, we estimated that the detection rate among this population had yet to exceed 40% as of late December and could have been as low as under 20%.

Compared with the 61-70 year old age group (Figure 9), we expected the detection rate of SARS-CoV-2 infections in individuals over 70 years old to be higher (22). Older people are more likely to show symptoms and thus seek testing, testing intensity generally increases with age, and routine testing in long-term care facilities may additionally identify asymptomatic infections. Infections among younger age groups might have been detected at a lower rate than that of the 61-70 age group because younger adults were less likely to present symptoms. An exception was in college-age young adults: routine testing by universities could have led to higher overall detection rates than 40% in this population, although the overall rate would have masked heterogeneity across the state and in different segments of the college-age population.