Abstract
Data from Chicago confirm the end of flu season coincides with the beginning of pollen season. The end of flu season also coincides with onset of seasonal aerosolization of mold spores. Overall, the data suggest bioaerosols, especially mold spores, compete with viruses for a shared receptor, with the periodicity of influenza-like illnesses, including COVID-19, a consequence of seasonal factors that influence aerosolization of competing species.
Introduction
Influenza-like illnesses (ILIs) attributable to influenza viruses and to coronaviruses are sharply seasonal.1,2 Importantly, recent data from The Netherlands indicate there exists an inverse relationship between the seasonal incidence of influenza-like illnesses, including COVID-19, and pollen count.3,4 To discern whether such a relationship might be the case generally, pollen count in Chicago was related to ILIs reported by local emergency departments. In Chicago, as in The Netherlands, ILIs fall as total pollen count rises.
Because bioaerosols measured in Chicago include not only pollens but also mold spores, ILIs were related to counts of both. Just as they do for pollens, ILIs fall as mold spore count rises. In contradistinction to their temporal relationship with pollens, however, ILIs remain low when mold spore count is high, rising again when mold spore count falls.
Perusal of the various measured pollens and mold spores reveals many are echinulated, having protuberances reminiscent of those of influenza and SARS-CoV-2 viruses.5–11 Such protuberances seem ideally suited to interacting with Toll-like receptor 4 (TLR4) in a fashion akin to ‘hook-and-loop’ adhesives. It is proposed that the various pollens and mold spores compete with relevant respiratory viruses for TLR4, limiting viral engagement and, consequently, ILIs.
Results
Kinetics of Presentations to Emergency Departments
As shown in Fig. 1A, ILI presentations to emergency departments in Chicago are cyclical in nature, with a periodicity of ~ 1 year. The peak incidence occurs during ~ February, with the annual nadir occurring during ~ August. Over the entirety of a seasonal cycle, ILI presentations never reach zero. Although there are subtleties associated with the kinetics of ILIs for each of the individual years 2015 through 2020, the annual increase in ILI presentations for each season is characterized by a leading ‘bump,’ followed thereafter by a major rise. The relevance of the bump is addressed below, under both Results and Discussion. As for the major rise, it is approximated empirically by a first-order process, ILIobs = ILI0ekt + C, where ILIobs is the number of observed ILIs during the growth phase of a cycle; ILI0 is the number of ILI presentations at the start of a cycle; t is time, in d; k is a first-order rate constant, in d-1; and C is the number of background presentations ‘masquerading’ as ILIs, Table 1 and Fig. 2A. The rate of decrease from any yearly maximum also fits reasonably well a first-order process, Table 1 and Fig. 2B. Taken at face value, the data suggest: 1) a not insignificant number of ILI presentations are not flu, and 2) the annual rate of change in flu cases is due to change in ambient concentration of influenza virus.
Fig. 3A shows the time course of the 7-day moving average of COVID-19 presentations to emergency departments of all Chicago hospitals. Because reporting was not uniformly rigorous before May 1, 2020, COVID-19 presentations prior to that date have been excluded from analyses. As shown in Fig. 3B, the fall in presentations for the period May 5, 2020 through September 27, 2020 was roughly first-order, with parameters k = 0.065 (0.055 – 0.077) d-1 (t1/2 ~ 10.6 [9 – 13] d), COVID0 = 875 (799 – 954) presentations and C = 262 (241 – 281) presentations.
Starting in late September 2020, COVID-19 cases in Chicago surged in first-order fashion, Fig. 3C, the parameters of which are k = 0.053 (0.050 – 0.057) d-1 (t1/2 ~ 13 [12– 14] d), COVID-190 = 96 (76 – 119) presentations and C = 117 (63 – 165) presentations. Taken at face value, the data suggest: 1) some of the individuals for whom a diagnosis of COVID-19 was made did not have COVID-19, and 2) changes in the number of COVID-19 presentations in Chicago were due to changes in the ambient concentration of SARS-CoV-2.
Kinetics of Pollen and Mold Spore Counts
Pollens are fertilizing elements of flowering plants whilst mold spores are reproductive elements of fungi. In published studies,3,4 pollens alone were counted and related to ILIs. Left uncounted were mold spores, important seasonal contributors to the total bioaerosol burden. For the studies reported herein, both pollens and mold spores were counted. Those counts were then analyzed, in aggregate and individually.
Usually, pollens and mold spores in Chicago are monitored from ~ mid-March to ~ mid-October, the period most problematic for persons suffering from seasonal allergies. As expected, the data indicate bioaerosol expression is cyclical with a periodicity of ~ 1 year, Fig. 1B. The total bioaerosol count peaks during ~ mid-September and falls sharply thereafter. Unfortunately, data following the peaks are somewhat limited, their collection being truncated on an arbitrary end date, i.e., ~ mid-October.
In the case of pollens, the seasonal distribution is bimodal, with a dominant first mode that peaks in ~ mid-May and a smaller second mode that peaks in ~ late August, Fig. 1C. The pollens that constitute the second mode, here termed ‘late pollens’, are predominantly Ambrosia and the other Asteraceae. Importantly, the peak of the second mode always coincides with the leading bump in ILI presentations, Figs. 1A and 4. The potential relevance of this is addressed under Discussion. In the case of mold spores, which constitute the bulk of the measured bioaerosols, Figs. 1B and 1D, and Table 2, the peak count, which occurs during ~ late September, falls precipitously by ~ mid-October, with an empiric half-life of ~ 10 d, Table 3.
Although these data substantiate the claim of an inverse relationship between the onset of pollen season and the end of flu season, pollen count declines rapidly and is not elevated when ILI presentations, Fig. 5A, and COVID-19 presentations, Fig. 6A, are low. Mold spores, on the other hand, increase continuously in first-order fashion, Table 3 and Fig. 7, beginning just prior to or coincident with the fall in ILI, Fig. 5B, and COVID-19 presentations, Fig. 6B, and across the entirety of the summer months, when influenza and COVID-19 cases are low.
Inhibition of ILI and COVID-19 Presentations by Mold Spores
If one assumes ILI and COVID-19 presentations are consequences of the binding of relevant viruses to specific receptors, then one can treat the presentations as proxies for those receptors, for which mold spores compete. Toward that end, ILI and COVID-19 presentations were plotted as functions of total mold spore count, Fig. 8. Because the curvatures of the plots suggest true equilibria, the data of each were fit to the equation P = Po/(1+C/Kd) + B, where P is the observed number of presentations to emergency departments, Po is the maximum number of such presentations, C, in mold spores/m3, is the measured mold spore count, Kd, in mold spores/m3, is the apparent dissociation constant of the receptor – mold spore complex, and B is a constant representing presentations not influenced by mold spores. As shown in the figures, the data of each plot fit reasonably well the theoretical model. From the ILI data, one calculates Po ~ 50 (41 – 60) presentations, Kd ~ 2128 (997 – 4242) mold spores/m3 and B ~ 16 (9 – 20) presentations; from the COVID-19 data, one calculates Po ~ 1366 (994 – 2679) presentations, Kd ~ 1668 (487 – 4143) mold spores/m3.and B ~ 201 (103 – 257) presentations. The most parsimonious explanation for the near equivalence of the apparent dissociation constants is a shared receptor.
Discussion
The data presented herein are consistent with those presented earlier by others,3,4 namely, the incidence of ILIs falls as pollen count rises. Because the data of the present study derive from an urban area in North America (Chicago, IL USA: latitude 41.85003, longitude −87.65005) whilst those of the earlier study derive from North Central Europe (Helmond, The Netherlands: latitude 51.48167, longitude 5.66111), it appears the inverse relationship may be generally valid.
With special regard to the late pollens, changes in their atmospheric concentration invariably coincide with the annual leading bump in ILIs. Inasmuch as Ambrosia, the dominant species, is a major respiratory allergen, the leading bump may represent ragweed sensitivities manifesting as ILI. Regardless, the peak in late pollen count --as if a switch --presages the major upswing in ILI, Fig. 4, and COVID-19, Fig. 6C, presentations. Thus, aside from any contribution to mechanistic understanding it might provide, the peak in late pollens could be exploited when contemplating an upcoming ILI season.
Separate and distinct from pollen count, mold spore count in Chicago correlates inversely with ILIs. Indeed, given their higher atmospheric concentration as well as the duration of their seasonal expression, mold spores seem more likely than pollens to be principals in any abatement of ILIs, including COVID-19. Mold spores and pollens could abate viral activity by either direct or indirect means. By direct means, they might produce substances that limit viral propagation or they might complex with viruses, limiting viral infectivity.13,14 But if direct antiviral activity is an attribute of the bioaerosols themselves, then one would not expect, a priori, significant disparity between individual susceptibilities to severe flu or COVID-19.15,16 As indirect means, others have proposed pollens stimulate the human immune system in such a way as to either potentiate endogenous antiviral activity or elicit a protective allergic response.3 Against these proposals, asthma does not confer protection against either influenza or COVID-19.17–19
The similarity of the proposed mold spore dose dependencies for abatement of flu and COVID-19 suggests a shared receptor. Although much attention has been given to angiotensin-converting enzyme 2 (ACE-2) and its role in COVID-19,20,21 there are compelling reasons to believe TLR4, which binds the SARS-CoV-2 spike protein with greater affinity than does ACE-2,22 is also operative: 1) TLR4 is implicated in the inflammatory response triggered by sharply seasonal respiratory viruses,23–25 2) TLR4 has a significant role in innate defense against multiple species of fungi,26,27 and polymorphisms in TLR4 are associated with invasive fungal disease,28,29 3) COVID-19 prognosis correlates with radiographic involvement of alveolar spaces,30,31 the epithelial surfaces of which are poor in ACE-232,33 but rich in TLR4,34 4) inflammation of the sort associated with acute lung injury is mediated by TLR4,35–44 and 5) age-dependent hyper-responsiveness of TLR445, especially in the context of interactions with TLR5,46,47 can account for the age-dependent severity of COVID-19. That TLR4 may be involved in the processing of bioaerosols is also expected on phylogenetic grounds: the receptor has been retained by some fish that breathe air, but lost by those that do not,48 and the eponymous Toll receptor controls the antifungal response of Drosophila.49
Given these, one can imagine the engagement of TLR4 by aerosols of all sorts, including, but not limited to, viruses, pollens and mold spores, in a fashion analogous to the engagement of hook-and-loop adhesives, i.e., Velcro®. Instead of loops, however, spinous processes of the various aerosols engage TLR4 ‘hooks,’ effecting an innate immune response, the nature of which depends on the arrangement and density of the engagement. And just as hook-and-loop adhesives can be rendered nonfunctional/dys-functional by nonspecific adherence of extraneous materials, so, too, might TLR4 hooks become saturated with one ligand, e.g., mold spores and/or pollens, to the exclusion of another, e.g., a respiratory virus.
The data presented herein bring new appreciation and understanding to seasonality and suggest a remarkable interplay between bioaerosols that influence the health of man. Indeed, inasmuch as humans have co-existed with plants, fungi and viruses for some time, it stands to reason that, over the course of evolution, the respiratory system of the former would have developed means to cope with the significant recurring, i.e., annual, inhalational exposure to reproductive elements of the latter. As the environment-facing interface of the respiratory tree, epithelial cells and their entourage of innate immune effectors seem ideally positioned to provide that coping mechanism.
Materials and Methods
Collection and Counting of Pollens and Mold Spores
A volumetric spore trap (Burkard Manufacturing, Hertfordshire, England) equipped with a 24 h sampling head was used to collect pollens and mold spores. The trap was fixed ~ 70 feet above ground, on a roof in Melrose Park, Il, USA. A standard glass microscope slide coated with grease was placed in a carriage that moved at a rate of 2 mm/h past the trap orifice (14 mm x 2 mm). Air was drawn through the orifice at a rate of 10 l/min, thereby impacting airborne particles against the greased slide. Slides so exposed were stained with glycerin jelly supplemented with basic fuchsin. After applying a coverslip, a slide was evaluated microscopically for both pollens and mold spores, Table 4. A new slide was placed in the trap daily, and the carriage was re-oriented to its start position. Bioaerosol counts were made Monday through Friday, generally between mid-March and mid-October.
ILI and COVID-19 Data
ILI data pooled from 23 large hospitals in Chicago over the period January 9, 2015 through July 18, 2020 were obtained from the Chicago Department of Public Health (CDPH). The 23 hospitals were chosen because they alone of Chicago-area hospitals consistently reported ILI presentations over the entirety of the study interval. The data are included in the Supplement, Table S1. COVID-19 data from all Chicago hospitals were obtained through portals of the CDPH, https://www.chicago.gov/city/en/sites/covid19/home/covid-dashboard.html and https://data.cityofchicago.org/browse?limitTo=datasets&sortBy=alpha&tags=covid-19. Those data are also included in the Supplement, Table S2.
Data Availability
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Data Analysis
Time- or dose-dependent data were paired with the corresponding time or dose and fit to equations described in the text. The best values for the parameters of the equations, as well as their corresponding 95% confidence intervals, were then determined using the paired data and a nonlinear least squares regression method.50
Author Contributions
R.B.S. collated all the data, and contributed to preparing and writing the manuscript. R.D.S. collected the pollen and mold spore data, and contributed to writing the manuscript. D.G.R. processed and analyzed the data, and contributed to preparing and writing the manuscript. A.C.R. analyzed the data and contributed to preparing and writing the manuscript. D.A.R. conceived parts of the manuscript and contributed to its preparation and writing. G.S.R. conceived and designed the study, processed and analyzed the data, and contributed to preparing and writing the manuscript.
Competing Interests
None of the authors has competing interests.
Acknowledgment
D.A.R. and G.S.R. thank their children, Andrew, Jonah, Ruth and Damien, for inspiration.