key: cord-0886153-00a53z3v authors: Shah, Richa B.; Shah, Rachna D.; Retzinger, Damien G.; Retzinger, Andrew C.; Retzinger, Deborah A.; Retzinger, Gregory S. title: Competing Bioaerosols May Influence the Seasonality of Influenza-Like Illnesses, including COVID-19. The Chicago Experience date: 2021-09-16 journal: Pathogens DOI: 10.3390/pathogens10091204 sha: 2d73dd42edae51c80f40732e2fd3ca3fb5fc9617 doc_id: 886153 cord_uid: 00a53z3v Data from Chicago confirm the end of flu season coincides with the beginning of pollen season. More importantly, 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. 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 ILIs, 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 the two measured species. Just as they do for pollens, ILIs fall as mold spores rise. In contradistinction to their temporal relationship with pollens, however, ILIs remain low when mold spores are high, rising again when mold spores fall. Perusal of the various measured pollens and mold spores reveals many are echinulated, having protuberances reminiscent of those of influenza and SARS-CoV-2 virions [5] [6] [7] [8] [9] [10] [11] . In addition, although interactions of viral 'spikes' with other host proteins have been documented [12, 13] , such protuberances seem ideally suited to interacting especially with Toll-like receptors in a fashion akin to 'hook-and-loop' adhesives. Indeed, such interactions are well described for a number of pollens and mold spores, and there is no reason to believe viruses, despite their smaller size, interact differently. Implicating Toll-like receptors, Tolllike receptor 4 (TLR4) in particular, seems appropriate on phenomenological grounds as well: (1) engagement of TLR 4 can account for the inflammatory signaling characteristic of severe COVID-19 [14, 15] , (2) TLR4 is intimately involved in the inflammation elicited both by sharply seasonal respiratory viruses and by multiple species of fungi [16] [17] [18] [19] [20] , ILI data pooled from 23 large hospitals in Chicago over the period 9 January 2015 through 18 July 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. In Chicago, the designation of ILI by emergency departments is based on fever (≥100 • F) and respiratory symptoms, i.e., cough and/or sore throat, not on any specific diagnosis. The data are Supplementary Table S3 . 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 (accessed on 18 July 2020) and https://data.cityofchicago. org/browse?limitTo=datasets&sortBy=alpha&tags=covid-19 (accessed on 18 November 2020) . Those data are also included in Supplementary Table S4 . 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 [29] . As shown in Figure 1A , ILI presentations to emergency departments in Chicago are cyclical in nature, with a periodicity of~1 year. The peak incidence occurs approximately during February, with the annual nadir occurring approximately 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, ILI obs = ILI 0 e kt + C, where ILI obs is the number of observed ILIs during the growth phase of a cycle; ILI 0 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 2 and Figure 2A ). The rate of decrease from any yearly maximum also fits reasonably well a first-order process (Table 2 and Figure 2B ). Taken at face value, the data suggest: (1) a not insignificant number of ILI presentations are due to pathogens other than influenza virus, e.g., parainfluenza virus, respiratory syncytial virus, measles, mumps, etc., and (2) the annual rate of change in influenza cases is due to change in ambient concentration of influenza virus. Figure 1A ) to a first-order rate equation, the parameters of which are given in Table 2 . The inset shows the expected linearity of the same data when plotted according to the method of Kézdy [30] , in this case ILI presentationsday n vs. ILI presentationsday n+1. Figure 1A ) to a first-order rate equation, the parameters of which are given in Table 2 . The inset shows the expected linearity of the same data when plotted according to the method of Kézdy [30] , in this case ILI presentations day n vs. ILI presentations day n+1 . Figure 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 1 May 2020, COVID-19 presentations prior to that date have been excluded from analyses. As shown in Figure 3B , the fall in presentations for the period 5 May 2020 through 27 September 2020 was roughly first-order, with parameters k = 0.065 d −1 (t 1/2~1 0.6 d), COVID-19 0 = 875 presentations and C = 262 presentations. Inasmuch as masks and physical-distancing were mandated in Chicago on 1 May 2020, the rate of fall in COVID-19 presentations, i.e., the shape of the curve, after that date was influenced to some extent by those measures, as discussed elsewhere [31] . Figure 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 1 May 2020, COVID-19 presentations prior to that date have been excluded from analyses. As shown in Figure 3B , the fall in presentations for the period 5 May 2020 through 27 September 2020 was roughly first-order, with parameters k = 0.065 d −1 (t1/2 ~ 10.6 d), COVID-190 = 875 presentations and C = 262 presentations. Inasmuch as masks and physical-distancing were mandated in Chicago on 1 May 2020, the rate of fall in COVID-19 presentations, i.e., the shape of the curve, after that date was influenced to some extent by those measures, as discussed elsewhere [31] . Starting in late September 2020, COVID-19 cases in Chicago surged in first-order fashion ( Figure 3C ), the parameters of which are k = 0.053 d −1 (t1/2 ~ 13 d), COVID-190 = 96 presentations and C = 117 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. 13 September 2020. The dashed line displayed in the primary plot is the theoretical fit of the data to a first-order rate equation, the parameters of which are given in the text. For more refined analysis, see [31] . (C) Appearance of COVID-19 presentations, 14 September 2020 to 14 November 2020. The dashed line displayed in the primary plot is the theoretical fit of the data to a first-order rate equation, the parameters of which are given in the text. The insets show expected linearities of the data when plotted according to the method of Kézdy [30] , in these cases COVID-19 presentations day n vs. COVID-19 presentations day n+1 . Starting in late September 2020, COVID-19 cases in Chicago surged in first-order fashion ( Figure 3C ), the parameters of which are k = 0.053 d −1 (t 1/2~1 3 d), COVID-19 0 = 96 presentations and C = 117 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. 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. In Chicago, pollens and mold spores are monitored from approximately mid-March to approximately mid-October, the only time during which the bioaerosols are easily measurable and also the time most problematic for persons suffering from seasonal allergies. As expected, the data indicate bioaerosol expression is cyclical with a periodicity of~1 year ( Figure 1B ). The total bioaerosol count peaks during approximately mid-September and falls sharply thereafter. Data following the peaks are somewhat limited, their collection being truncated on an arbitrary end date, i.e., approximately mid-October. In the case of pollens, the seasonal distribution is bimodal, with a dominant first mode that peaks in approximately mid-May and a smaller second mode that peaks in approximately late August ( Figure 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 ( Figures 1A and 4) . The potential relevance of this is addressed in the Discussion. In the case of mold spores, which constitute the bulk of the measured bioaerosols (Figure 1B,D and Table 3) , the peak count, which occurs during approximately late September, falls precipitously by approximately mid-October, with an empiric half-life of~10 d, Table 4 . given in the text. The insets show expected linearities of the data when plotted according to the method of Kézdy [30] , in these cases COVID-19 presentationsday n vs. COVID-19 presentationsday n+1. 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. In Chicago, pollens and mold spores are monitored from approximately mid-March to approximately mid-October, the only time during which the bioaerosols are easily measurable and also the time most problematic for persons suffering from seasonal allergies. As expected, the data indicate bioaerosol expression is cyclical with a periodicity of ~1 year ( Figure 1B ). The total bioaerosol count peaks during approximately mid-September and falls sharply thereafter. Data following the peaks are somewhat limited, their collection being truncated on an arbitrary end date, i.e., approximately mid-October. In the case of pollens, the seasonal distribution is bimodal, with a dominant first mode that peaks in approximately mid-May and a smaller second mode that peaks in approximately late August ( Figure 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 ( Figures 1A and 4) . The potential relevance of this is addressed in the Discussion. In the case of mold spores, which constitute the bulk of the measured bioaerosols ( Figures 1B and 1D and Table 3) , the peak count, which occurs during approximately late September, falls precipitously by approximately mid-October, with an empiric half-life of ~10 d, Table 4 . 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 ( Figure 5A ) and COVID-19 presentations ( Figure 6A ) are low. Mold spores, on the other hand, increase continuously in first-order fashion (Table 4 and Figure 7 ), beginning just prior to or coincident with the fall in ILI ( Figure 5B ) and (Figure 6B) , and across the entirety of the summer months, when influenza and COVID-19 cases are low. 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 ( Figure 5A ) and COVID-19 presentations ( Figure 6A ) are low. Mold spores, on the other hand, increase continuously in first-order fashion (Table 4 and Figure 7 ), beginning just prior to or coincident with the fall in ILI ( Figure 5B ) and COVID-19 presentations (Figure 6B) , and across the entirety of the summer months, when influenza and COVID-19 cases are low. Figure 4 for additional details. The onset of aerosolization of mold spores also coincides with the drop in seasonal ILI presentations. Thereafter, mold spore count increases across the entirety of the summer months-during which time ILIs remain low-and falls precipitously from mid-September to mid-October, at which time ILI cases begin to rise. See text for additional details. Table 4 . The inset shows the expected linearity of the same data when plotted according to the method of Kézdy [30] , in this case mold spore countday n vs. mold spore countday n+1. 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 (Figure 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/m 3 , is the measured mold spore count; Kd, in mold spores/m 3 , 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 the theoretical model reasonably well. From the ILI data, one calculates Po ~ 50 presentations, Kd ~ 2128 mold spores/m 3 and B ~ 16 presentations; from the COVID-19 data, one calculates Po ~ 1366 presentations, Kd ~ 1668 mold spores/m 3 and B ~ 201 presentations. The most parsimonious explanation for the near equivalence of the apparent dissociation constants is a shared receptor. Table 4 . The inset shows the expected linearity of the same data when plotted according to the method of Kézdy [30] , in this case mold spore count day n vs. mold spore count day n+1 . 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 (Figure 8 ). Because the curvatures of the plots suggest true equilibria, the data of each were fit to the equation P = P o /(1 + C/K d ) + B, where P is the observed number of presentations to emergency departments; P o is the maximum number of such presentations; C, in mold spores/m 3 , is the measured mold spore count; K d , in mold spores/m 3 , 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 the theoretical model reasonably well. From the ILI data, one calculates P o~5 0 presentations, K d~2 128 mold spores/m 3 and B~16 presentations; from the COVID-19 data, one calculates P o~1 366 presentations, K d~1 668 mold spores/m 3 and B~201 presentations. The most parsimonious explanation for the near equivalence of the apparent dissociation constants is a shared receptor. 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, 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. Made blatantly obvious by these studies are the seasonalities of ILIs, pollens and mold spores. The annual periodicities of the three indicate the rotation of the earth about the sun is ultimately responsible. Special note should be made of: (1) the nearness of the onset of bioaerosol expression to the vernal equinox, i.e., when the lengths of day and night are nearly equal, (2) the nearness of the peak in pollen count (excepting Ambrosia) to the summer solstice, i.e., the longest day of the year and (3) the nearness of the peak in mold spore count to the autumnal equinox. Because light and heat from the sun are drivers of both natural and agricultural growing seasons, these dates and their relevance to the expression and dispersal of bioaerosols should come as no surprise. 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. Alternatively, it may just represent enhanced spread of respiratory viruses by, for example, school openings. Regardless, the peak in late pollen count-as if a switch-presages the major upswing in ILI ( Figure 4 ) and COVID-19 ( Figure 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. From an anthropologic perspective, the potential of mold spores and pollens to inhibit influenza-like epidemics/pandemics, including COVID-19, certainly has great relevance and significant consequence. Still, because in comparison to plants, fungi and even viruses, humans contribute only very little to the biomass on planet Earth [32] , it seems likely some larger purpose is served by interplay between the three bioaerosols. It is tempting to speculate that any antiviral effect attributable to mold spores and/or pollens is intended to benefit primarily fungi and plants [33] , i.e., the human benefit, albeit perhaps related mechanistically, is an epiphenomenon. As just one of many possibilities, mold spores and pollens might protect primarily arthropods, birds and bats, organisms intimately involved in dissemination of reproductive elements of both fungi and plants [34] [35] [36] [37] [38] [39] [40] [41] . 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. Because certain mold spores, e.g., Aspergillus, can propagate in man if left unattended by innate immune effectors, it also follows mold spores should be prioritized over pollens. Nonviral bioaerosols 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 [42, 43] . However, 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 [44, 45] . 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 [46] [47] [48] . 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 [12, 13] , there are compelling reasons to believe TLR4, which binds the SARS-CoV-2 spike protein with greater affinity than does ACE-2 [49] , is also operative: 1) TLR4 is implicated in the inflammatory response triggered by sharply seasonal respiratory viruses [16] [17] [18] , 2) TLR4 has a significant role in innate defense against multiple species of fungi [19, 20] and polymorphisms in TLR4 are associated with invasive fungal disease [50, 51] , 3) COVID-19 prognosis correlates with radiographic involvement of alveolar spaces [21, 22] , the epithelial surfaces of which are poor in ACE-2 [52, 53] but rich in TLR4 [23] , (4) inflammation of the sort associated with acute lung injury is mediated by TLR4 [14, [54] [55] [56] [57] [58] [59] [60] [61] [62] , (5) age-dependent hyper-responsiveness of TLR4 [24] , especially in the context of interactions with TLR5 [25, 26] , can account for the age-dependent severity of COVID-19 and (6) fibrino(gen) D-dimers are markedly elevated in persons with severe COVID-19 [28] . 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 [63] , and the eponymous Toll receptor controls the antifungal response of Drosophila [64] . Given these, one can imagine the engagement of TLR4 by aerosols of all sorts and microscopic/submicroscopic sizes, including, but not limited to, viruses (diam 0.01-0.30 µm), mold spores (diam 1-50 µm) and pollens (diam 10-1000 µm), 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. The large surface area of the extracellular domain of TLR4, 6000-8500 Å 2 ensures accommodation of many such protuberances which, in turn, explains the broad specificity of the receptor [65] . In addition, just as hook-and-loop adhesives can be rendered nonfunctional/dysfunctional by nonspecific adherence of extraneous materials, so too might TLR4 hooks become saturated with one ligand to the exclusion of another. As for the role of fibrin(ogen) D-domains, their overexpression as endogenous ligands may represent an attempt by the innate immune system to purge/disengage TLR4 from pathogenic aerosols of all sorts for, perhaps, restorative purpose. 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, considering that 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. That being the case, nebulized materials that exploit competition either between the various bioaerosols or between the bioaerosols and endogenous TLR4 ligands, e.g., C-terminus of the fibrinogen γ-chain [66] [67] [68] , might prove therapeutic. Finally, and notwithstanding allergic potential, the indoor cultivation-including mold-rich fertilization-of pollenating plants might be exploited to limit occurrence of sharply seasonal ILIs. Influenza Seasonality: Underlying Causes and Modeling Theories Coronavirus Occurrence and Transmission over 8 Years in the HIVE Cohort of Households in Michigan Pollen Explains Flu-like and COVID-19 Seasonality Pollen Likely Seasonal Factor in Inhibiting Flu-like Epidemics. A Dutch Study into the Inverse Relation between Pollen Counts, Hay Fever and Flu-like Incidence Molecular and Immunological Characterization of Ragweed (Ambrosia artemisiifolia L) Pollen after Exposure of the Plants to Elevated Ozone over a Whole Growing Season Scanning Electron Microscopy of Penicillium Conidia A Comparative Analysis of Pollinator Type and Pollen Ornamentation in the Araceae and the Arecaceae, Two Unrelated Families of the Monocots Identification of Fungi of the Genus Aspergillus Section nigri Using Polyphasic Taxonomy From Colony to Rodlet. A Six Meter Long Portrait of the Xerophilic Fungus Aspergillus restrictus Decorates the Hall of the Westerdijk Institute Scanning Electron Microscopy of Ascospores of Schwanniomyces A Pneumonia Outbreak Associated with a New Coronavirus of Probably Bat Origin Structural and Functional Basis of SARS-CoV-2 Entry by Using Human ACE2 COVID-19 Patients Upregulate Toll-like Receptor 4-Mediated Inflammatory Signaling that Mimics Bacterial Sepsis COVID-19 and Toll-Like Receptor 4 (TLR4): SARS-CoV-2 May Bind and Activate TLR4 to Increase ACE2 Expression, Facilitating Entry and Causing Hyperinflammation Novel Signaling Interactions between Proteinase-Activated Receptor 2 and Toll-like Receptors In vitro and In vivo Select Targeting of Intracellular Toll-Interleukin-1 Receptor Resistance Domains for Protection Against Influenza-Induced Disease Novel Strategies for Targeting Innate Immune Responses to Influenza Innate Immunity Induced by the Major Allergen Alt a 1 from the Fungus Alternaria is Dependent upon Toll-like Receptors 2/4 in Human Toll-like Receptors as Key Mediators in Innate Antifungal Immunity COVID-19 Patients: Correlation with Disease Severity and Short-term Prognosis Maximum Chest CT Score is Associated with Progression to Severe Illness in Patients with COVID-19: A Retrospective Study from Wuhan Expression of Functional Toll-like Receptor-2 and -4 on Alveolar Epithelial Cells Aging is Associated with Chronic Immune Activation and Dysregulation of Monocyte Phenotype and Function Age-Associated Elevation in TLR5 Leads to Increased Inflammatory Responses in the Elderly Mapping Tenascin-C Interaction with Toll-like Receptor 4 Reveals a New Subset of Endogenous Inflammatory Triggers Prothrombotic Disturbances of Hemostasis of Patients with Severe COVID-19: A Prospective Longitudinal Observational Study Evaluation of Parameter Uncertainties in Nonlinear Regression Using Microsoft Excel Spreadsheet Optimizing the use of the Kézdy-Mangelsdorf-Swinbourne Method for Analysis of Data Following A exp Estimate of Benefit Attributable to Wearing Masks in Chicago during the Early Days of the Pandemic Plant-Insect-Microbe Interaction: A Love Triangle between Enemies in Ecosystem Ecology and Evolution of Insect-Fungus Mutualisms Insects as Flower Visitors and Pollinators Non-Bee Insects as Visitors and Pollinators of Crops: Biology, Ecology, and Management Dispersal of Fungi Spores by Non-Specialized Flower-Visiting Birds Fur versus Feathers: Pollen Delivery by Bats and Hummingbirds and Consequences for Pollen Production Short-Distance Pollen Dispersal by Bats in an Urban Setting: Monitoring the Movement of Vertebrate Pollinator through Fluorescent Dyes Associated with Arthropods from Bat Hibernacula in Eastern Canada, with Particular Reference to Pseudogymnoascus destructans Immune System Modulation and Viral Persistence in Bats: Understanding Viral Spillover Antimicrobial and Antiviral Properties of Different Types of Propolis Risk Factors for Severe Outcomes Following 2009 Influenza A (H1N1) Infection: A Global Pooled Analysis Epidemiological Characteristics of Patients with Severe COVID-19 Infection in Wuhan, China: Evidence from a Retrospective Observational Study Risk of Severe Influenza among Adults with Chronic Medical Conditions Prevalence and Characterization of Asthma in Hospitalized and Nonhospitalized Patients with COVID-19 Prevalence of Underlying Diseases in Died Cases of COVID-19: A Systematic Review and Meta-analysis In silico Studies on the Comparative Characterization of the Interactions of SARS-CoV-2 Spike Glycoprotein with ACE-2 Receptor Homologs and Human TLRs Increased Susceptibility for Aspergillosis and Post-Transplant Immune Deficiency in Patient with Gene Variants of TLR4 after Stem Cell Transplantation Polymorphisms in Toll-like Receptor Genes and Susceptibility to Pulmonary Aspergillosis The Protein Expression Profile of ACE2 in Human Tissues ACE2 Localizes to the Respiratory Cilia and is Not Increased by ACE Inhibitors or ARBs Role of Human Toll-like Receptors in Naturally Occurring Influenza A Infections. Influenza Other Respir. Viruses Toll-like Receptor 4-Mediated Activation of p38 Mitogen-Activated Protein Kinase is a Determinant of Respiratory Virus Entry and Tropism Involvement of TLR2 and TLR4 in Inflammatory Immune Responses Induced by Fine and Coarse Ambient Air Particulate Matter Fungal Surface and Innate Immune Recognition of Filamentous Fungi Fungal Pathogens-a Sweet and Sour Treat for Toll-like Receptors Innate Response to Pollen Allergens Pollen Lipids can Play a Role in Allergic Airway Inflammation Identification of Oxidative Stress and Toll-like Receptor 4 Signaling as a Key Pathway of Acute Lung Injury The Role of TLR4 in the Pathogenesis of Indirect Acute Lung Injury Toll-like Receptors in Bony Fish: From Genomics to Function The Dorsoventral Regulatory Gene Cassette spätzle/Toll/cactus Controls the Potent Antifugal Response in Drosophila Adults Leucine-rich Repeats and Pathogen Recognition in Toll-like Receptors Crystal structure of a 30 kDa C-Terminal Fragment from the γ Chain of Human Fibrinogen Correlating Structure and Function during the Evolution of Fibrinogen-Related Domains Related Proteins in Tissue Repair: How a Unique Domain with a Common Structure Controls Diverse Aspects of Wound Healing Acknowledgments: D.A.R. and G.S.R. thank their children, Andrew, Jonah, Ruth and Damien, for motivation and inspiration. The authors declare no conflict of interest.