Next Article in Journal / Special Issue
Biomarkers Linked with Dynamic Changes of Renal Function in Asymptomatic and Mildly Symptomatic COVID-19 Patients
Previous Article in Journal
Polygenic Panels Predicting the Susceptibility of Multiple Upper Aerodigestive Tract Cancer in Oral Cancer Patients
Previous Article in Special Issue
The Challenging Diagnosis of Pediatric Multisystem Inflammatory Syndrome Associated with Sars-Cov-2 Infection-Two Case Reports and Literature Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SARS-CoV-2-Laden Respiratory Aerosol Deposition in the Lung Alveolar-Interstitial Region Is a Potential Risk Factor for Severe Disease: A Modeling Study

Department of Biosciences, Paris Lodron University of Salzburg (PLUS), 5020 Salzburg, Austria
*
Author to whom correspondence should be addressed.
Submission received: 9 March 2021 / Revised: 10 May 2021 / Accepted: 17 May 2021 / Published: 19 May 2021
(This article belongs to the Special Issue COVID-19 Related Complications)

Abstract

:
COVID-19, predominantly a mild disease, is associated with more severe clinical manifestation upon pulmonary involvement. Virion-laden aerosols and droplets target different anatomical sites for deposition. Compared to droplets, aerosols more readily advance into the peripheral lung. We performed in silico modeling to confirm the secondary pulmonary lobules as the primary site of disease initiation. By taking different anatomical aerosol origins into consideration and reflecting aerosols from exhalation maneuvers breathing and vocalization, the physicochemical properties of generated respiratory aerosol particles were defined upon conversion to droplet nuclei by evaporation at ambient air. To provide detailed, spatially-resolved information on particle deposition in the thoracic region of the lung, a top-down refinement approach was employed. Our study presents evidence for hot spots of aerosol deposition in lung generations beyond the terminal bronchiole, with a maximum in the secondary pulmonary lobules and a high preference to the lower lobes of both lungs. In vivo, initial chest CT anomalies, the ground glass opacities, resulting from partial alveolar filling and interstitial thickening in the secondary pulmonary lobules, are likewise localized in these lung generations, with the highest frequency in both lower lobes and in the early stage of disease. Hence, our results suggest a disease initiation right there upon inhalation of virion-laden respiratory aerosols, linking the aerosol transmission route to pathogenesis associated with higher disease burden and identifying aerosol transmission as a new independent risk factor for developing a pulmonary phase with a severe outcome.

1. Introduction

Virus shedding via virion-laden respiratory aerosol particles by contagious individuals has been confirmed for the common seasonal coronavirus, influenza virus and rhinoviruses [1,2]. Similarly, exposure to SARS-CoV-2-laden aerosols by sharing crowded indoor space can be expected as a substantial transmission risk with a seasonal pattern [3]. In the context of long-distance travel, this was demonstrated by a detailed analysis of 2334 index patients and 72,093 close contacts for high-speed train passengers in China [4]. The data revealed that the main determinants for transmission included co-exposure time and spatial distance. In occupational scenarios, the European Centre for Disease Prevention and Control reported that the vast majority (95%) of 447 investigated clusters occurred completely or predominantly indoors across all professional categories [5]. These reports are aligned with observations of cluster-forming events associated with recreational indoor activities and miscellaneous social gatherings [6,7]. It has also become increasingly acknowledged that aerosol transmission may play a significant, and not only a minor, role in confined spaces [8,9,10].
It has been demonstrated that SARS-CoV-2 is viable for hours in aerosols [11]. Experimental data on aerosol exhalation during different expiratory maneuvers describe size distribution, mechanisms of generation and sites of origin within the respiratory tract [12,13,14]. Due to the very specific particokinetic properties of aerosols derived from vocalization and breathing, the majority of these ambient air respiratory aerosol particles (ARAPs) resist gravitational settling and stay airborne for extended time [15,16]. Hence, virion-laden particles may contribute to short- and long-range airborne transmission by deposition and subsequent initiation of COVID-19 in the respiratory tract.
The etiology of COVID-19, however, has been predominantly linked with airborne droplet or contact transmission, with initial virus deposition and replication starting at the mucosa of the oral-nasal cavity followed by propagation through the upper respiratory tract epithelia to the alveolar-interstitial region with vulnerable cell types [17,18]. In cases with pulmonary involvement a stepwise migration via conducting airways towards the most distal regions of the lungs has been proposed [19,20]. This concept might need revision when considering transmission by immediately inhalable virion-laden aerosols with inherent potential for disease initiation in the peripheral lung. Outside the field of viral respiratory infections and with lung diseases caused by inhalable toxic particulate matter it is accepted that the site of initial deposition is also the site of initial adverse effects [21]. Such a reverse approach on disease propagation would have implications for pharmacological treatment considerations.
The hypothesis that ARAPs have the potential to initiate disease with pulmonary involvement, associated with a higher disease burden and characterized by pneumonia, dyspnea and reduced blood oxygen saturation [22], is supported by radiology findings on chest CTs. At the onset of the pandemic, classification of specific chest CT image anomalies has become key to support differential diagnosis of COVID-19-induced lung inflammation from pneumonia of other origin and the radiologist’s performance achieved sensitivity and specificity close to 90% [23,24]. Image classification by computerized deep learning methods, in particular convolutional neuronal networks, which rely on strong correlations between lesion areas in radiologic images and clinical indicators, achieve a specificity of >99% [25]. CT assessment includes the number, type, rate of anomaly development and particularly the distribution pattern in the lung. The predominant finding in chest CTs of patients (including presymptomatic patients) [26] was the peripheral distribution of ground-glass opacity (GGO) that occurred in the early stages of lung involvement [27]. This was described as a locally-confined, hazy increase in attenuation of secondary pulmonary lobules (units of three to five terminal bronchioles contained by fibrous septa) caused by partial filling and collapse of alveoli, interstitial thickening and increased capillary blood flow [28]. The key features of GGO that were observed with high frequency were bilateral involvement, posterior part or lower lobe predilection, prominent peripheral and sub-pleural distribution, absence of peri-bronchial localization [27]. Peripheral GGOs were observed in 86% (18/21) [29] and 80% (175/219) [30] of patients in two cohorts in China. In a longitudinal study of 90 cases where temporal changes in GGOs were investigated [31], the authors observed a gradual dissemination from unilateral to bilateral lung involvement with an increasing number of lesions and afflicted lobes. At the onset of symptoms, the abnormalities were predominantly sub-pleural with a peak at day 6 to 11. Another study of 121 cases that also mapped data on GGOs, demonstrated increased bilateral involvement with time (early-phase 28%, 10/36; late-phase 88%, 22/25) [32]. The highest frequency was observed in both lower lobes (65%, 79/121; 63%, 76/121).
With mounting evidence and intensified debate on the role of aerosol transmission and involvement thereof in disease etiology and pathophysiology, it is compelling to combine in silico aerosol lung deposition modeling to supplement in vivo clinical observations and to generate a synergistic new view (Figure 1).
Literature research reveals that little is known about the potential of SARS-CoV-2-laden respiratory aerosols for direct inoculation and disease initiation in the human lung. Elucidation of the key events in the initiation of COVID-19 would be highly desirable, e.g., understanding the potential contribution of virion-laden aerosols to the tissue-delivered dose; and spatially-resolved distribution of their deposition onto lung epithelia of the bronchial, bronchiolar and alveolar region. Hence, to establish a spatially-resolved (detailed lung generation vs. ARAP size) deposition heat map, a five-step in silico modeling of virion-laden aerosols originating from vocalization combined with modeling of respiratory tract deposition was performed. Our study specifically focused on aerosols that originate from vocalization maneuvers, and were chosen because these aerosols are available for transmission via ambient air for hours and possibly generated by asymptomatic, but infected individuals. In addition, it has been proposed that, to a large extent, these aerosol size characteristics are not dependent on vocalization frequency and amplitude; thus, they are sex and age independent [37].

2. Materials and Methods

A stepwise modeling and analysis procedure was performed as depicted in Figure 2.

2.1. Respiratory Aerosol Particle Modeling

The respiratory aerosol particle (RAP) modeling is based on the aerosol size distribution studies of Morawska et al. [14] and Johnson et al. [12] for RAPs originating from human exhalation activities. Their results are in good agreement with reported particle size distributions from human activity [38]. From this comprehensive set of exhalation maneuvers that also included coughing, our investigation was limited to RAP data representing vocalization and breathing. To generate the data, the authors basically combined two experimental systems: an expiratory droplet investigation system (EDIS) that includes an aerodynamic particle sizer (APS) to measure from 0.5 µm to 20 µm; and a droplet deposition analysis (DDA) to measure particles larger than 20 µm. The measured particle size distribution was then analyzed by curve fitting. Five distinct particle populations, represented by five size distributions with corresponding mid-point diameters, were identified and assigned to distinct exhalation maneuvers and anatomical origins. These five mid-point diameters, hereafter referred to as mode 1 to mode 5, were used for the in silico particle deposition simulation in our study.
At the point of exhalation, RAPs are in evaporation equilibrium with the presumable breath cloud, i.e., 90% (+/− 7%) relative humidity (RH) and 28°C (+/− 1°) [12]. With exposure to ambient air, these exhaled respiratory aerosol particles (ERAPs) immediately begin to converge at a new evaporation equilibrium. For ERAPs at room temperature (RT), a RH <60%, and a particle size ≤20 µm, the transition to ambient air respiratory aerosol particles (ARAPs) is in the order of milliseconds to tenths of seconds [13,14]. In the literature, ARAPs ≤5 µm are frequently termed “droplet nuclei”, and RAPs beyond this size as “droplets”. Nevertheless, the concept of droplet nuclei formation is not limited to these size thresholds. According to Nicas et al. [39], and Holmgren et al. [40], irrespective of particle size, the evaporation shrinkage factor from ERAPs to ARAPs was experimentally determined at 0.5. Required for the calculation of particokinetic parameters, the physicochemical properties for ERAP and ARAP were established. Each particle was defined as spherically shaped object containing a single virion (diameter 0.12 µm, density 1.2 g/cm3) [41] in association with respiratory mucus (ERAP) or mucus remnants, when being evaporated to an equilibrium with ambient air (ARAP). The definition of one virion per particle is supported by the literature which, for example, reports 106−1011 viral RNAs per milliliter in COVID-19 patients’ sputum [42]. If the highest concentration is aerosolized in our 1.6 µm ERAPs (0.8 µm ARAPs), this results in less than 1 virion per aerosol particle. Thus, given that the adsorption of an individual virion by an aerosol is a random event, most particles are free of virions and the remaining are loaded by just one virion. This observation was long before confirmed by Couch et al. [43] in the context of adenoviruses, where the vast majority had only one virion and this broadly concurs with a review of Poon et al. [44]. As shown in Table 1, the composition of the mucus associated with ERAPs is defined [45].
Due to the negligible contribution of the virion volume to the ERAP volume, for all particle sizes (1.6, 3.6, 7.0, 11.0, 145.0 µm), the density of the ERAPs was assumed to be equivalent to the mucus density of 1.041 g/cm3. The density calculation for ARAPs was determined with an evaporation shrinkage factor of 0.5 and elimination of water content with the following Equations (1)–(5):
dARAP = EF × dERAP
VARAP = 1/6 × π × d3ARAP
m evaporated H2O = ρ H2O × (VERAP − VARAP)
mARAP = mERAP − mevaporated H2O
ρARAP = mARAP/VARAP
d, diameter; V, volume; m, mass; ρ, density; EF, evaporation shrinkage factor
Table 2 provides the calculated characteristics of ARAPs used in our in silico deposition simulation and refers to the anatomical origin during exhalation activity.

2.2. Lung Modeling

The simulation of ARAP deposition is based on the anatomical regions and airway generation model defined by International Commission on Radiological Protection (ICRP) [36] and the Yeh/Schum five-Lobe model [33]. In brief, the respiratory system is divided into four anatomical regions: the extrathoracic region, the bronchial region (BB), the bronchiolar region (bb) and the alveolar-interstitial region (AI) [47]. Together, BB, bb and AI constitute the thoracic region. The “Yeh/Schum 5-Lobe model” reflects an asymmetric human lung with five lobes and describes the airways within each lobe in a specific single-path manner. From this model, the human respiratory tract was ultimately compartmentalized into 110 regions; beginning with the common anatomical region of the trachea and expanding to the most distal alveolar region of each individual lobe of the lungs. Figure 3 depicts the anatomical segmentation of the thoracic region used in our simulation.

2.3. Deposition Simulation and Analysis

To generate the ARAP deposition heat map for each of the five particle modes, a series of calculations and models were performed as briefly described here. Firstly, an ARAP exposure dose by mass was defined and converted into an exposure dose by particle number. Secondly, Multiple Path Particle Dosimetry (MPPD) [34,35] was employed to simulate ARAP deposition per thoracic region. Thirdly, an ICRP-model was used to calculate the extrathoracic deposition and the in-/exhaled particle fractions. Finally, the combined data were used to calculate an ARAP deposition probability for a single inhaled ARAP to generate the final respiratory tract deposition heat map.

2.4. Simulation of Thoracic ARAP Deposition via in Silico Model MPPD

MPPD supports mathematical modeling of aerosol deposition in the human lung. Numerical solution of equations governing air flow, particokinetics and particle deposition is based on physical and physiological parameters of the lung and physicochemical determinants of aerosols under investigation. MPPD allows us to predict total, regional, lobar and local deposition in the lung. In our study we applied MPPD with the Yeh/Schum 5-Lobe lung morphometry model. It provides the advantage of a whole-lung model; while newer computational fluid dynamics models (CFD) are still regional, they are either limited to simulate the upper respiratory tract and airway deposition, hence, without the acinar lung generations, or deep lung simulation, they are missing aerosol deposition in the bronchial and bronchiolar region [48,49]. Yeh/Schum allows a good predictive accuracy for particle deposition in the deep lung, in good agreement and validated with in vivo experimental data or new CFD models like the 3D model of the deep lung, termed DLM [50]. Table 3 provides MPPD input parameters used:
For each distinct thoracic structure, the number of deposited particles per breath was extracted from the created MPPD detailed report (output field “Tot. Dep. Particle s (#)”). The sum of deposited particles over all thoracic structures per breath represents the overall number of deposited particles in the thoracic region per breath.

2.5. Calculation of Extrathoracic Deposition via the ICRP Model

Based on selected ARAP exposure concentration, particle mode and tidal volume used in the MPPD model, the particle number per breath was calculated for each specific ARAP mode. Extrathoracic deposited particle number per breath, and aerosol particle number not deposited, were calculated according to ICRP Equations (6) and (7) as follows:
IF = 1 − 0.5 × (1 − 1/(1 + 0.00076 × d2.8))
DF = IF × (0.0587 + 0.911/(1 + e4.77+1.485×ln d) + 0.943/(1 + e0.503−2.58×ln d))
TD = DF × TP
ETP = TD − TTD
NDP = TD − ETP − TTD
d, diameter; IF, inhalable fraction; DF, total deposited fraction; TD, total deposited particle number per breath; TP, particle number per breath; ETP, extrathoracic deposited particle number per breath; TTD, overall deposited particle number in the thoracic region per breath; NDP, not deposited particle number per breath.

2.6. Calculation of Deposition Probability for Thoracic Deposition Based on Total Deposited Particles

To establish the thoracic deposition probability, the number of deposited particles per breath for each specific ARAP mode was expressed in relation to the total deposited particle numbers per breath.

3. Results

The ARAP respiratory tract deposition was analyzed at three different levels using a top-down approach. Five ARAP modes, exhaled during vocalization and breathing, represented by five size distribution modes that originated from bronchiolar and laryngeal region and the oral cavity were considered. Firstly, the overall deposition rate, that was segregated into extrathoracic and thoracic deposition fractions (Figure 4a), was investigated. Considerable thoracic deposition was observed for mode 1 (0.8 µm), mode 2 (1.8 µm) and mode 3 (3.5 µm). Mode 1 ARAPs were particularly notable with a high inhalation-exhalation rate of 69.1%. If the particles were retained, these also had the highest thoracic deposition probability of >50%. Across all ARAP modes, the ARAPs represented by mode 2 showed the highest absolute deposition probability for the thoracic region (24.8%). Mode 3 and mode 4 (5.5 µm) showed a predominant extrathoracic deposition; while mode 5 (72.5 µm) lacked any thoracic deposition. From all the particle modes generated through vocalization, mode 1 ARAPs captured the highest proportion and represented 72.8% of all ARAPs. This was followed by mode 2 (21.0%), and ARAP modes 3 to 5 contributed to 6.2% [12,14]. Weighting all modes against the relative proportion showed that only ARAP mode 1 and 2 contribute to thoracic deposition in vivo (Figure 4b).
A more detailed approach was then performed to provide spatial distribution information on ARAP deposition. This was achieved by dissecting the information into lung lobe-specific data for the three thoracic respiratory tract regions BB, bb and AI. Irrespective of particle size, Figure 5 depicts a strikingly dominant deposition probability in the AI region of all lobes followed by bb, and then BB with the lowest deposition probability. The two lower lobes have a considerably higher burden than the middle right and the upper lobes; thus, indicating that the AI regions of the lower lobes are potential “hot spots” for ARAP deposition.
The third level in refining the analysis of ARAP deposition dissected the thoracic region into up to 25 lung generations for each lobe; beginning at common central structures and expanded into the most distal alveolar region. The two-dimensional (ARAP mode vs. lung lobe generation) deposition heat map (Figure 6a) immediately revealed deposition “hot spots” that begin at the proximal alveolar regions (PAR) of both lower lobes and intensify towards the distal alveolar regions until the penultimate generation. Taking into consideration the relative abundance of the different ARAP modes from vocalization, the weighted heat map (Figure 6b) again stressed the major contribution of ARAP mode 1 and the minor contributions of the other modes. Probability values and MPPD-derived raw data are available at http://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4736854, uploaded on 2 December 2020.

4. Discussion

Our modeling approach focuses on providing detailed data on deposition of virion-laden droplet nuclei in the thoracic region of the respiratory tract; thus, screening for anatomical sites with increased susceptibility for disease initiation. Data reveal that all ARAP modes that are associated with vocalization and originate from bronchiolar or laryngeal fluid film burst have the potential to be retained after inhalation. Conversely, the oral cavity-originating aerosols are not retained. Comparing the different sized ARAPs against relative abundance in the exhalation plume revealed that mode 1 and 2 ARAPs significantly contribute to a thoracic tissue-delivered dose; and ultimately, exclude all other modes. If mode 1 and 2 ARAPs alone are considered, then 43.9% of the deposited ARAPs appear in the thoracic region; thus, emphasizing the specific role of these ARAPs in the etiology of COVID-19 in vivo. To further elaborate the peculiarities of thoracic deposition, data were mapped to all five lung lobes. Calculating the deposition probability in the BB, bb and AI regions enabled differentiation of central from peripheral deposition. This analysis revealed the striking tendency of both lower lung lobes to be involved, and a preferential peripheral deposition in the AI region that exceeded central and BB deposition by up to three-fold. To date, this observed specific spatial pattern of ARAP deposition is aligned with the spatio-temporal distribution preferences of radiology chest CT anomalies in COVID-19 patients. Thus, at least in patients with pulmonary involvement, virion-laden ARAP inhalation seems to be linked to disease initiation. Finally, the chosen top-down modeling approach provided data refined from common central structures of all lobes to the most distal regions of 25 lung generations. The resulting deposition heat maps enabled recognition of deposition “hot spots” and the associated particle modes. The most prominent “hot spots” were observed in both lower lobes, beginning at the PAR and extending to the penultimate alveolar region. This was not different for our investigated aerosol upon rehydration to its size before evaporation started (ERAP), which would resemble the situation of regrowth after inhalation. The deposition probability values and heatmaps are available at http://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4736854, uploaded on 2 December 2020.
At this point, the correlation between spatially-resolved ARAP deposition and chest CT anomalies, predominantly peripheral COVID-19-associated GGOs in the absence of peribronchial GGOs, is evident and is in agreement with findings of GGOs as a phenomenon of secondary pulmonary lobules in early disease [28]. These specific sites in the peripheral lung, especially the PAR, are also a preferred area for alveolar and interstitial particle accumulation and inflammatory, infiltrative or tissue remodeling processes after exposure to low-toxicity, low-solubility inhalable aerosols [21]. Our results elucidate major differences in the tissue-delivered dose of virion-laden aerosols between the AI, bb, BB and extrathoracic regions. The study revealed a high probability of deposition of virion-laden aerosols at the alveolar epithelia in the AI region. Thus, we propose that these structures are particularly vulnerable to the initiation of adverse events after deposition of pathogenic near- and sub-micron-sized ARAPs. In the alveolar region, the protective barrier consists of a thin (average 200 nm) [51] fluid lining with an aqueous hypophase of low viscosity. At the cellular level, this means that each SARS-CoV-2 virion that is deposited by diffusion-driven mass transport can effectively and almost immediately interact with cellular receptors. A different scenario occurs when SARS-CoV-2 virion-laden aerosol particles are deposited on ciliated airway epithelia in the bb and BB regions. In healthy individuals, these areas are more effectively protected by a highly-viscous biphasic mucus layer with a thickness of approximately 20 µm [52]. A high clearance rate by the mucociliary escalator provides an outwardly-directed trajectory with a velocity exceeding 1 mm/min in the BB region. SARS-CoV-2 virions deposited on this mucus layer are dependent on non-directed mass transfer by diffusion towards the targeted epithelia. Physicochemical properties such as the viscosity and thickness of the mucus layer, however, markedly reduce the probability of timely epithelial contact by diffusional translocation; a prerequisite for subsequent disease initiation by viral entry. Due to the high mucus outward flow rate, a very low tissue-delivered dose can also be anticipated at the original site of deposition in the bb and BB regions. This is because the initially-deposited virion dose, that slowly penetrates the mucus layer via stochastic principles, would be distributed over a wide epithelial area resulting in a very high local dilution factor. Hence, implications in the early phases of COVID-19 pathogenesis at different epithelia should be expected. By taking into consideration the local dilution factor, the thickness and viscosity differences of the fluid linings in the AI and BB regions, the variation in disease initiation potential of a deposited virion dose could exceed a magnitude of 103, which is also supported by Thomas et al. [38]. From the perspective of an immunologist, this would aid in explaining differences in disease severity. A low number of virus replicates, at early stage, in the presence of nevertheless sufficient viral antigen exposure in the BB region would initiate a standard immune response without excessive immune activation and with a low potential of advancing into the frequently-observed severe course of COVID-19.
A minimum infectious dose for SARS-CoV-2 is not established and is expected to be highly variable between different tissues, mainly depending on the availability of the viral-docking receptors and the major differences of protective barriers and clearance mechanisms. However, it is proposed for other viral respiratory infections that a single virus can serve for disease initiation [39]. More specifically, a recent study concerning the SARS-CoV-2 infectious dose in the oral-nasal-cavity proposed an infectious dose of 300 PFU [53]. Another key publication for SARS-CoV reported 43 PFU correlated to 10% infected individuals [54]. Overall, these studies propose that a low number of virions is sufficient for disease initiation. This assumption is practically underpinned by the 2020 Skagit Valley Chorale super-spreading incident, where a single COVID-19 carrier infected 52 out of 61 choral members, during a reported 2.5 h co-exposure [55].
Aside from the aforementioned aspects, the cellular regulation of viral docking receptors ACE2 and TMPRSS2 that moderate SARS-CoV-2 entry and, thus, the rate of viral replication in different lung epithelia, has to be considered. Noteworthily, ACE2 gene expression was found to be low in airway and alveolar epithelial cells in healthy individuals [17], but ACE2 gene was identified as interferon-stimulated in type II pneumocytes [56] upon viral challenge as part of a host-tissue protective mechanism, thereby particularly enhancing SARS-CoV-2 susceptibility and viral replication in the alveolar-interstitial region after an initial virus dose. This is well aligned with studies reporting ACE2 and TMPRSS2 co-expression in type II alveolar cells [57,58], but not in the upper airways [59]. Of note, new gene ontology data confirm ACE2 co-expression with immune response genes, in particular pro-inflammatory cytokines, chemokines and interferon γ-inducible protein 16 (IFI16) [60], an innate immune sensor regulating the interferon response to viral infections [61,62] by inflammasome activation. In this context, the particular susceptibility of type II pneumocytes is of major importance. These cells are progenitor cells for type I pneumocytes, and hence they are crucial for epithelial repair. Without restoration of alveolar tissue homeostasis COVID-19 can rapidly progress in disease severity with signs of diffuse alveolar damage and loss of functionality. This may in particular contribute to the increased fatality rate of the elderly age groups in SARS-CoV-2 acute lung injury, as demonstrated in a mouse model of alveolar epithelial type II cell senescence associated with impaired type II cell function in acute lung injury [63].
Limitations of our study are primarily due to the fact that the modeling presents an idealized version of actuality. In real life, virion-laden droplets would arrive with a possibly non-Gaussian size distribution. Although exhalation-generated aerosol particles are generally independent of body size, sex and age, the same cannot be said for inhalation; particularly when pre-existing lung conditions are taken into consideration with high variability in lung morphometry. These limitations have to be accepted applying a whole lung model operating at reasonable computational costs. Validation of in silico aerosol deposition models with in vivo experiments, by applying emerging techniques [64,65], would be desirable. For a general validation it would have to be done in healthy patients, which is probably a big hurdle for ethics review boards.

5. Conclusions

In conclusion, the stepwise modeling approach presented here showed that virion-laden ARAPs are primarily deposited in the secondary pulmonary lobules in lung generations beyond the terminal bronchiole. These regions are identical with the peripheral localization of GGOs observed in chest CTs; an early marker of COVID-19 infection. This agreement suggests that virion-laden ARAPs from vocalization maneuvers are significantly associated with COVID-19 pathology in the lung. Therefore, SARS-CoV-2 transmission via shared indoor aerosols appears to be an important route of infection in the peripheral lung and, thus, constitutes an independent risk factor for severe courses of disease. Our findings imply an alignment of public health care policies (i) to rebalance the risk of contact, large droplet and aerosol transmission, (ii) to address the particularity of indoor exposure and long-range transmission and (iii) to consider the substantial risk for more severe disease manifestations. Disease initiation in the alveolar-interstitial region, instead of the oral-nasal cavity with subsequent propagation towards the peripheral lung, implies optimization of treatment recommendations and disease monitoring. Both diagnosis of peripheral lung involvement early after COVID-19 confirmation by molecular tests, still in the asymptomatic stage, or early onset of pulmonary symptoms after confirmed or suspected exposure should be interpreted as a sign of disease initiation in the peripheral lung as primary anatomical site of infection. This pathogenesis imposes an increased likelihood for progression to severe outcome and, thus, increased public health care burden.
Linking our ARAP deposition results with recently provided respiratory viral emission estimation data [66] or with real world aerosol concentrations in indoor air would allow us to refine infection risk assessment and risk controlling strategies. Such combined models could be of great relevance for epidemiologists in the context of SARS-CoV-2 mutational variants of concern, such as 20B/501Y.V1, VOC 202012/01 (or B.1.1.7 lineage) with increased viral copy numbers [67] in the oral-nasal cavity. A higher fraction of virion-laden ARAPs with increased virion copy numbers increases the probability to pass the threshold for disease initiation, first of all, at the deposition hot spots identified here.

Author Contributions

S.H.: N.H., and M.H. conceived the study; S.H. and N.H. performed modeling and analyzed the data, and all authors participated in the writing, reviewing, and editing of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the EU H2020 projects NanoRigo and NanoCommons (grant numbers 814530 and 731032, respectively).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study or openly available in Zenodo at http://0-doi-org.brum.beds.ac.uk/10.5281/zenodo.4301641.

Acknowledgments

This work was supported by the Allergy-Cancer-BioNano (ACBN) Research Center of the Paris Lodron University of Salzburg (PLUS) and in parts by the EU H2020 projects NanoRigo and NanoCommons (grant numbers 814530 and 731032, respectively).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Leung, N.H.; Chu, D.K.; Shiu, E.Y.; Chan, K.-H.; McDevitt, J.J.; Hau, B.J.; Yen, H.-L.; Li, Y.; Ip, D.K.; Peiris, J.M. Respiratory virus shedding in exhaled breath and efficacy of face masks. Nat. Med. 2020, 26, 676–680. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Milton, D.K.; Fabian, M.P.; Cowling, B.J.; Grantham, M.L.; McDevitt, J.J. Influenza Virus Aerosols in Human Exhaled Breath: Particle Size, Culturability, and Effect of Surgical Masks. PLoS Pathog. 2013, 9, e1003205. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Moriyama, M.; Hugentobler, W.J.; Iwasaki, A. Seasonality of respiratory viral infections. Annu. Rev. Virol. 2020, 7, 83–101. [Google Scholar] [CrossRef] [PubMed]
  4. Hu, M.; Lin, H.; Wang, J.; Xu, C.; Tatem, A.J.; Meng, B.; Zhang, X.; Liu, Y.; Wang, P.; Wu, G.; et al. Risk of Coronavirus Disease 2019 Transmission in Train Passengers: An Epidemiological and Modeling Study. Clin. Infect. Dis. 2021, 72, 604–610. [Google Scholar] [CrossRef]
  5. ECDC. COVID-19 Clusters and Outbreaks in Occupational Settings in the EU/EEA and the UK; European Centre for Disease Prevention and Control: Stockholm, Sweden, 2020. [Google Scholar]
  6. Correa-Martínez, C.L.; Kampmeier, S.; Kümpers, P.; Schwierzeck, V.; Hennies, M.; Hafezi, W.; Kühn, J.; Pavenstädt, H.; Ludwig, S.; Mellmann, A. A pandemic in times of global tourism: Superspreading and exportation of COVID-19 cases from a ski area in Austria. J. Clin. Microbiol. 2020, 58. [Google Scholar] [CrossRef] [Green Version]
  7. Hamner, L. High SARS-CoV-2 attack rate following exposure at a choir practice—Skagit County, Washington, March 2020. Mmwr. Morb. Mortal. Wkly. Rep. 2020, 69, 606–610. [Google Scholar] [CrossRef]
  8. Anderson, E.L.; Turnham, P.; Griffin, J.R.; Clarke, C.C. Consideration of the Aerosol Transmission for COVID-19 and Public Health. Risk Anal. 2020, 40, 902–907. [Google Scholar] [CrossRef]
  9. Alwan, N.A.; Burgess, R.A.; Ashworth, S.; Beale, R.; Bhadelia, N.; Bogaert, D.; Dowd, J.; Eckerle, I.; Goldman, L.R.; Greenhalgh, T.; et al. Scientific consensus on the COVID-19 pandemic: We need to act now. Lancet 2020, 396, e71–e72. [Google Scholar] [CrossRef]
  10. Meyerowitz, E.A.; Richterman, A.; Gandhi, R.T.; Sax, P.E. Transmission of SARS-CoV-2: A review of viral, host, and environmental factors. Ann. Intern. Med. 2020, M20-5008. [Google Scholar] [CrossRef]
  11. Van Doremalen, N.; Bushmaker, T.; Morris, D.H.; Holbrook, M.G.; Gamble, A.; Williamson, B.N.; Tamin, A.; Harcourt, J.L.; Thornburg, N.J.; Gerber, S.I.; et al. Aerosol and Surface Stability of SARS-CoV-2 as Compared with SARS-CoV-1. N. Engl. J. Med. 2020, 382, 1564–1567. [Google Scholar] [CrossRef]
  12. Johnson, G.R.; Morawska, L.; Ristovski, Z.D.; Hargreaves, M.; Mengersen, K.; Chao, C.Y.H.; Wan, M.P.; Li, Y.; Xie, X.; Katoshevski, D.; et al. Modality of human expired aerosol size distributions. J. Aerosol Sci. 2011, 42, 839–851. [Google Scholar] [CrossRef]
  13. Wei, J.; Li, Y. Airborne spread of infectious agents in the indoor environment. Am. J. Infect. Control 2016, 44, S102–S108. [Google Scholar] [CrossRef] [PubMed]
  14. Morawska, L.; Johnson, G.R.; Ristovski, Z.D.; Hargreaves, M.; Mengersen, K.; Corbett, S.; Chao, C.Y.H.; Li, Y.; Katoshevski, D. Size distribution and sites of origin of droplets expelled from the human respiratory tract during expiratory activities. J. Aerosol Sci. 2009, 40, 256–269. [Google Scholar] [CrossRef] [Green Version]
  15. Xie, X.; Li, Y.; Chwang, A.; Ho, P.; Seto, W. How far droplets can move in indoor environments-revisiting the Wells evaporation-falling curve. Indoor Air 2007, 17, 211–225. [Google Scholar] [CrossRef] [PubMed]
  16. Balachandar, S.; Zaleski, S.; Soldati, A.; Ahmadi, G.; Bourouiba, L. Host-to-host airborne transmission as a multiphase flow problem for science-based social distance guidelines. Int. J. Multiph. Flow 2020, 132, 103439. [Google Scholar] [CrossRef]
  17. Aguiar, J.A.; Tremblay, B.J.M.; Mansfield, M.J.; Woody, O.; Lobb, B.; Banerjee, A.; Chandiramohan, A.; Tiessen, N.; Cao, Q.; Dvorkin-Gheva, A.; et al. Gene expression and in situ protein profiling of candidate SARS-CoV-2 receptors in human airway epithelial cells and lung tissue. Eur. Respir. J. 2020, 56, 2001123. [Google Scholar] [CrossRef]
  18. Freni, F.; Meduri, A.; Gazia, F.; Nicastro, V.; Galletti, C.; Aragona, P.; Galletti, C.; Galletti, B.; Galletti, F. Symptomatology in head and neck district in coronavirus disease (COVID-19): A possible neuroinvasive action of SARS-CoV-2. Am. J. Otolaryngol. 2020, 41, 102612. [Google Scholar] [CrossRef] [PubMed]
  19. Mason, R.J. Thoughts on the alveolar phase of COVID-19. Am. J. Physiol. Lung Cell Mol. Physiol. 2020, 319, L115–L120. [Google Scholar] [CrossRef]
  20. Mason, R.J. Pathogenesis of COVID-19 from a cell biology perspective. Eur. Respir. J. 2020, 55, 2000607. [Google Scholar] [CrossRef] [Green Version]
  21. Donaldson, K.; Borm, P.; Oberdorster, G.; Pinkerton, K.E.; Stone, V.; Tran, C. Concordance between in vitro and in vivo dosimetry in the proinflammatory effects of low-toxicity, low-solubility particles: The key role of the proximal alveolar region. Inhal. Toxicol. 2008, 20, 53–62. [Google Scholar] [CrossRef]
  22. Wu, Z.; McGoogan, J.M. Characteristics of and important lessons from the coronavirus disease 2019 (COVID-19) outbreak in China: Summary of a report of 72 314 cases from the Chinese Center for Disease Control and Prevention. JAMA 2020, 323, 1239–1242. [Google Scholar] [CrossRef]
  23. Bai, H.X.; Wang, R.; Xiong, Z.; Hsieh, B.; Chang, K.; Halsey, K.; Tran, T.M.L.; Choi, J.W.; Wang, D.-C.; Shi, L.-B. Artificial Intelligence augmentation of radiologist performance in distinguishing COVID-19 from pneumonia of other origin on chest CT. Radiology 2020, 201491. [Google Scholar] [CrossRef] [PubMed]
  24. Lieveld, A.; Kok, B.; Schuit, F.; Azijli, K.; Heijmans, J.; van Laarhoven, A.; Assman, N.; Kootte, R.; Olgers, T.; Nanayakkara, P. Diagnosing COVID-19 pneumonia in a pandemic setting: Lung Ultrasound versus CT (LUVCT) A multi-centre, prospective, observational study. ERJ Open Res. 2020, 6, 539. [Google Scholar] [CrossRef] [PubMed]
  25. Liang, S.; Liu, H.; Gu, Y.; Guo, X.; Li, H.; Li, L.; Wu, Z.; Liu, M.; Tao, L. Fast automated detection of COVID-19 from medical images using convolutional neural networks. Commun. Biol. 2021, 4, 35. [Google Scholar] [CrossRef] [PubMed]
  26. Shi, H.; Han, X.; Jiang, N.; Cao, Y.; Alwalid, O.; Gu, J.; Fan, Y.; Zheng, C. Radiological findings from 81 patients with COVID-19 pneumonia in Wuhan, China: A descriptive study. Lancet Infect. Dis. 2020, 20, 425–434. [Google Scholar] [CrossRef]
  27. Zu, Z.Y.; Jiang, M.D.; Xu, P.P.; Chen, W.; Ni, Q.Q.; Lu, G.M.; Zhang, L.J. Coronavirus disease 2019 (COVID-19): A perspective from China. Radiology 2020, 296, 200490. [Google Scholar] [CrossRef] [Green Version]
  28. Sharma, A.; Eisen, J.E.; Shepard, J.-A.O.; Bernheim, A.; Little, B.P. Case 25-2020: A 47-Year-Old Woman with a Lung Mass. N. Engl. J. Med. 2020, 383, 665–674. [Google Scholar] [CrossRef] [PubMed]
  29. Ng, M.-Y.; Lee, E.Y.; Yang, J.; Yang, F.; Li, X.; Wang, H.; Lui, M.M.-s.; Lo, C.S.-Y.; Leung, B.; Khong, P.-L. Imaging profile of the COVID-19 infection: Radiologic findings and literature review. Radiol. Cardiothorac. Imaging 2020, 2, e200034. [Google Scholar] [CrossRef] [Green Version]
  30. Bai, H.X.; Hsieh, B.; Xiong, Z.; Halsey, K.; Choi, J.W.; Tran, T.M.L.; Pan, I.; Shi, L.B.; Wang, D.C.; Mei, J.; et al. Performance of Radiologists in Differentiating COVID-19 from Non-COVID-19 Viral Pneumonia at Chest CT. Radiology 2020, 296, E46–E54. [Google Scholar] [CrossRef] [PubMed]
  31. Wang, Y.; Dong, C.; Hu, Y.; Li, C.; Ren, Q.; Zhang, X.; Shi, H.; Zhou, M. Temporal changes of CT findings in 90 patients with COVID-19 pneumonia: A longitudinal study. Radiology 2020, 296, 200843. [Google Scholar] [CrossRef] [Green Version]
  32. Bernheim, A.; Mei, X.; Huang, M.; Yang, Y.; Fayad, Z.A.; Zhang, N.; Diao, K.; Lin, B.; Zhu, X.; Li, K. Chest CT findings in coronavirus disease-19 (COVID-19): Relationship to duration of infection. Radiology 2020, 295, 200463. [Google Scholar] [CrossRef] [Green Version]
  33. Yeh, H.-C.; Schum, G. Models of human lung airways and their application to inhaled particle deposition. Bull. Math. Biol. 1980, 42, 461–480. [Google Scholar] [CrossRef]
  34. Anjilvel, S.; Asgharian, B. A multiple-path model of particle deposition in the rat lung. Fundam. Appl. Toxicol. 1995, 28, 41–50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Miller, F.J.; Asgharian, B.; Schroeter, J.D.; Price, O. Improvements and additions to the multiple path particle dosimetry model. J. Aerosol Sci. 2016, 99, 14–26. [Google Scholar] [CrossRef] [Green Version]
  36. ICRP. ICRP Publication 66: Human Respiratory Tract Model for Radiological Protection; Ann: ICRP 24(1-3); ICRP: Stockholm, Sweden, 1994; Volume 66. [Google Scholar]
  37. Asadi, S.; Wexler, A.S.; Cappa, C.D.; Barreda, S.; Bouvier, N.M.; Ristenpart, W.D. Aerosol emission and superemission during human speech increase with voice loudness. Sci. Rep. 2019, 9, 1–10. [Google Scholar] [CrossRef] [Green Version]
  38. Thomas, R.J. Particle size and pathogenicity in the respiratory tract. Virulence 2013, 4, 847–858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Nicas, M.; Nazaroff, W.W.; Hubbard, A. Toward Understanding the Risk of Secondary Airborne Infection: Emission of Respirable Pathogens. J. Occup. Environ. Hyg. 2005, 2, 143–154. [Google Scholar] [CrossRef]
  40. Holmgren, H.; Bake, B.; Olin, A.-C.; Ljungström, E. Relation Between Humidity and Size of Exhaled Particles. J. Aerosol Med. Pulm. Drug Deliv. 2011, 24, 253–260. [Google Scholar] [CrossRef]
  41. Thurber, R.V.; Haynes, M.; Breitbart, M.; Wegley, L.; Rohwer, F. Laboratory procedures to generate viral metagenomes. Nat. Protoc. 2009, 4, 470. [Google Scholar] [CrossRef]
  42. Bar-On, Y.M.; Flamholz, A.; Phillips, R.; Milo, R. SARS-CoV-2 (COVID-19) by the numbers. eLife 2020, 9, e57309. [Google Scholar] [CrossRef]
  43. Couch, R.B.; Knight, V.; Douglas, R.G., Jr.; Black, S.H.; Hamory, B.H. The minimal infectious dose of adenovirus type 4; the case for natural transmission by viral aerosol. Trans. Am. Clin. Clim. Assoc. 1969, 80, 205–211. [Google Scholar]
  44. Poon, W.C.K.; Brown, A.T.; Direito, S.O.L.; Hodgson, D.J.M.; Le Nagard, L.; Lips, A.; MacPhee, C.E.; Marenduzzo, D.; Royer, J.R.; Silva, A.F.; et al. Soft matter science and the COVID-19 pandemic. Soft Matter. 2020, 16, 8310–8324. [Google Scholar] [CrossRef] [PubMed]
  45. Lai, S.K.; Wang, Y.-Y.; Wirtz, D.; Hanes, J. Micro- and macrorheology of mucus. Adv. Drug Deliv. Rev. 2009, 61, 86–100. [Google Scholar] [CrossRef] [Green Version]
  46. App, E.; Zayas, J.; King, M. Rheology of mucus and transepithelial potential difference: Small airways versus trachea. Eur. Respir. J. 1993, 6, 67–75. [Google Scholar]
  47. Hussain, M.; Madl, P.; Khan, A. Lung deposition predictions of airborne particles and the emergence of contemporary diseases, Part-I. Health 2011, 2, 51–59. [Google Scholar]
  48. Islam, M.S.; Paul, G.; Ong, H.X.; Young, P.M.; Gu, Y.T.; Saha, S.C. A Review of Respiratory Anatomical Development, Air Flow Characterization and Particle Deposition. Int. J. Envrion. Res. Public Health 2020, 17, 380. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Mutuku, J.K.; Hou, W.-C.; Chen, W.-H. An Overview of Experiments and Numerical Simulations on Airflow and Aerosols Deposition in Human Airways and the Role of Bioaerosol Motion in COVID-19 Transmission. Aerosol Air Qual. Res. 2020, 20, 1172–1196. [Google Scholar] [CrossRef]
  50. Koullapis, P.; Stylianou, F.; Sznitman, J.; Olsson, B.; Kassinos, S. Towards whole-lung simulations of aerosol deposition: A model of the deep lung. J. Aerosol Sci. 2020, 105541. [Google Scholar] [CrossRef]
  51. Knudsen, L.; Ochs, M. The micromechanics of lung alveoli: Structure and function of surfactant and tissue components. Histochem. Cell Biol. 2018, 150, 661–676. [Google Scholar] [CrossRef] [Green Version]
  52. Atanasova, K.R.; Reznikov, L.R. Strategies for measuring airway mucus and mucins. Respir. Res. 2019, 20, 261. [Google Scholar] [CrossRef]
  53. Basu, S. Computational characterization of inhaled droplet transport to the nasopharynx. Sci. Rep. 2021, 11, 6652. [Google Scholar] [CrossRef] [PubMed]
  54. Watanabe, T.; Bartrand, T.A.; Weir, M.H.; Omura, T.; Haas, C.N. Development of a Dose-Response Model for SARS Coronavirus. Risk Anal. 2010, 30, 1129–1138. [Google Scholar] [CrossRef] [Green Version]
  55. Miller, S.L.; Nazaroff, W.W.; Jimenez, J.L.; Boerstra, A.; Buonanno, G.; Dancer, S.J.; Kurnitski, J.; Marr, L.C.; Morawska, L.; Noakes, C. Transmission of SARS-CoV-2 by inhalation of respiratory aerosol in the Skagit Valley Chorale superspreading event. Indoor Air 2021, 31, 314–323. [Google Scholar] [CrossRef]
  56. Ziegler, C.G.K.; Allon, S.J.; Nyquist, S.K.; Mbano, I.M.; Miao, V.N.; Tzouanas, C.N.; Cao, Y.; Yousif, A.S.; Bals, J.; Hauser, B.M.; et al. SARS-CoV-2 Receptor ACE2 Is an Interferon-Stimulated Gene in Human Airway Epithelial Cells and Is Detected in Specific Cell Subsets across Tissues. Cell 2020, 181, 1016–1035.e1019. [Google Scholar] [CrossRef] [PubMed]
  57. Zou, X.; Chen, K.; Zou, J.; Han, P.; Hao, J.; Han, Z. Single-cell RNA-seq data analysis on the receptor ACE2 expression reveals the potential risk of different human organs vulnerable to 2019-nCoV infection. Front. Med. 2020, 14, 1–8. [Google Scholar] [CrossRef] [Green Version]
  58. Qi, F.; Qian, S.; Zhang, S.; Zhang, Z. Single cell RNA sequencing of 13 human tissues identify cell types and receptors of human coronaviruses. Biochem. Biophys. Res. Commun. 2020, 526, 135–140. [Google Scholar] [CrossRef] [PubMed]
  59. Hamming, I.; Timens, W.; Bulthuis, M.; Lely, A.; Navis, G.V.; van Goor, H. Tissue distribution of ACE2 protein, the functional receptor for SARS coronavirus. A first step in understanding SARS pathogenesis. J. Pathol. A J. Pathol. Soc. Great Br. Irel. 2004, 203, 631–637. [Google Scholar] [CrossRef]
  60. Wruck, W.; Adjaye, J. SARS-CoV-2 receptor ACE2 is co-expressed with genes related to transmembrane serine proteases, viral entry, immunity and cellular stress. Sci. Rep. 2020, 10, 21415. [Google Scholar] [CrossRef]
  61. Unterholzner, L.; Keating, S.E.; Baran, M.; Horan, K.A.; Jensen, S.B.; Sharma, S.; Sirois, C.M.; Jin, T.; Latz, E.; Xiao, T.S.; et al. IFI16 is an innate immune sensor for intracellular DNA. Nat. Immunol. 2010, 11, 997–1004. [Google Scholar] [CrossRef] [Green Version]
  62. Kerur, N.; Veettil, M.V.; Sharma-Walia, N.; Bottero, V.; Sadagopan, S.; Otageri, P.; Chandran, B. IFI16 acts as a nuclear pathogen sensor to induce the inflammasome in response to Kaposi Sarcoma-associated herpesvirus infection. Cell Host Microbe 2011, 9, 363–375. [Google Scholar] [CrossRef] [Green Version]
  63. Yazicioglu, T.; Mühlfeld, C.; Autilio, C.; Huang, C.-K.; Bär, C.; Dittrich-Breiholz, O.; Thum, T.; Pérez-Gil, J.; Schmiedl, A.; Brandenberger, C. Aging impairs alveolar epithelial type II cell function in acute lung injury. Am. J. Physiol. Lung Cell. Mol. Physiol. 2020, 319, L755–L769. [Google Scholar] [CrossRef]
  64. Thompson, R.B.; Finlay, W.H. Using MRI to measure aerosol deposition. J. Aerosol Med. Pulm. Drug Deliv. 2012, 25, 55–62. [Google Scholar] [CrossRef] [PubMed]
  65. Biddiscombe, M.F.; Meah, S.N.; Underwood, S.R.; Usmani, O.S. Comparing lung regions of interest in gamma scintigraphy for assessing inhaled therapeutic aerosol deposition. J. Aerosol Med. Pulm. Drug Deliv. 2011, 24, 165–173. [Google Scholar] [CrossRef] [PubMed]
  66. Buonanno, G.; Stabile, L.; Morawska, L. Estimation of airborne viral emission: Quanta emission rate of SARS-CoV-2 for infection risk assessment. Environ. Int. 2020, 141, 105794. [Google Scholar] [CrossRef] [PubMed]
  67. Kidd, M.; Richter, A.; Best, A.; Mirza, J.; Percival, B.; Mayhew, M.; Megram, O.; Ashford, F.; White, T.; Moles-Garcia, E. S-variant SARS-CoV-2 is associated with significantly higher viral loads in samples tested by ThermoFisher TaqPath RT-QPCR. medRxiv 2020. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the study concept and modeling. Virion-laden aerosol particles that are exhaled in the course of vocalization are immediately converted to ambient air respiratory aerosol particles (ARAPs) and remain airborne for hours. In silico modeling was used to determine whether these particles, if inhaled by exposed individuals, advance to the thoracic region and lead to preferred deposition in the distal regions of the lung. The spatial correlation of deposition hot spots with CT anomaly pattern (marked in red) would confirm disease initiation in the peripheral lung by these aerosols. For in silico simulation a five-step modeling and analysis procedure (indicated by grey boxes) was performed: step one—defining the physicochemical properties of exhaled respiratory aerosol particle upon generation; step two—establishing the physicochemical and particokinetic properties of these particles after conversion to ARAPs by evaporation to equilibrium with ambient air conditions; step three—establishing a detailed lung data model for ARAP deposition based on the Yeh/Schum 5-Lobe model [33]; step four—in silico simulation of ARAP deposition via Multiple Path Particle Dosimetry model [34,35] and International Commission on Radiological Protection model [36]; step five—refining resulting raw data by a top-down approach to provide detailed, spatially-resolved information on particle deposition in the thoracic region of the lung.
Figure 1. Schematic representation of the study concept and modeling. Virion-laden aerosol particles that are exhaled in the course of vocalization are immediately converted to ambient air respiratory aerosol particles (ARAPs) and remain airborne for hours. In silico modeling was used to determine whether these particles, if inhaled by exposed individuals, advance to the thoracic region and lead to preferred deposition in the distal regions of the lung. The spatial correlation of deposition hot spots with CT anomaly pattern (marked in red) would confirm disease initiation in the peripheral lung by these aerosols. For in silico simulation a five-step modeling and analysis procedure (indicated by grey boxes) was performed: step one—defining the physicochemical properties of exhaled respiratory aerosol particle upon generation; step two—establishing the physicochemical and particokinetic properties of these particles after conversion to ARAPs by evaporation to equilibrium with ambient air conditions; step three—establishing a detailed lung data model for ARAP deposition based on the Yeh/Schum 5-Lobe model [33]; step four—in silico simulation of ARAP deposition via Multiple Path Particle Dosimetry model [34,35] and International Commission on Radiological Protection model [36]; step five—refining resulting raw data by a top-down approach to provide detailed, spatially-resolved information on particle deposition in the thoracic region of the lung.
Jpm 11 00431 g001
Figure 2. Overview of the five modeling and analysis steps. Bullet points: • activity; ➢ outcome. ERAP, exhaled respiratory aerosol particle; ARAP, ambient air respiratory aerosol particle; MPPD, Multiple Path Particle Dosimetry model; ICRP, International Commission on Radiological Protection model.
Figure 2. Overview of the five modeling and analysis steps. Bullet points: • activity; ➢ outcome. ERAP, exhaled respiratory aerosol particle; ARAP, ambient air respiratory aerosol particle; MPPD, Multiple Path Particle Dosimetry model; ICRP, International Commission on Radiological Protection model.
Jpm 11 00431 g002
Figure 3. Anatomical structure of the thoracic region according to Yeh/Schum 5-Lobe model applied in the ARAP deposition simulation. Common structures are indicated on the far left. G, generation number; RU, right upper; RM, right middle; RL, right lower; LU, left upper; LL, left lower.
Figure 3. Anatomical structure of the thoracic region according to Yeh/Schum 5-Lobe model applied in the ARAP deposition simulation. Common structures are indicated on the far left. G, generation number; RU, right upper; RM, right middle; RL, right lower; LU, left upper; LL, left lower.
Jpm 11 00431 g003
Figure 4. Respiratory tract deposition of the inhalable fraction of ARAPs; originating from exhalation activity “voiced counting”, a combined process of vocalization and breathing. These maneuvers resulted in an ARAP size distribution represented by five modes (0.8, 1.8, 3.5, 5.5, 72.5 µm) that originated from bronchiolar and laryngeal fluid film burst and the oral cavity. Depending on size, different fractions of ARAPs are deposited in the thoracic or extrathoracic region of the respiratory tract, or are not retained and exhaled. Deposition data are shown (a) unweighted and (b) weighted by relative abundance in the exhalation plume (mode 1–5: 72.8%, 21.0%, 2.2%, 3.4%, 0.6% [12,14]).
Figure 4. Respiratory tract deposition of the inhalable fraction of ARAPs; originating from exhalation activity “voiced counting”, a combined process of vocalization and breathing. These maneuvers resulted in an ARAP size distribution represented by five modes (0.8, 1.8, 3.5, 5.5, 72.5 µm) that originated from bronchiolar and laryngeal fluid film burst and the oral cavity. Depending on size, different fractions of ARAPs are deposited in the thoracic or extrathoracic region of the respiratory tract, or are not retained and exhaled. Deposition data are shown (a) unweighted and (b) weighted by relative abundance in the exhalation plume (mode 1–5: 72.8%, 21.0%, 2.2%, 3.4%, 0.6% [12,14]).
Jpm 11 00431 g004
Figure 5. Deposition probability of ARAPs in the thoracic region, based on total deposited particle number (extrathoracic, thoracic). Deposition probability is shown for all five ARAP modes and is dissected further to the thoracic regions of the bronchial region (BB, green), bronchiolar region (bb, yellow) and alveolar-interstitial region (AI, red) for each lobe.
Figure 5. Deposition probability of ARAPs in the thoracic region, based on total deposited particle number (extrathoracic, thoracic). Deposition probability is shown for all five ARAP modes and is dissected further to the thoracic regions of the bronchial region (BB, green), bronchiolar region (bb, yellow) and alveolar-interstitial region (AI, red) for each lobe.
Jpm 11 00431 g005
Figure 6. Heat map for the probability of ARAP deposition in the thoracic region based on total deposited particle number (extrathoracic, thoracic). Shown are all thoracic sub-structures and all ARAP modes. As lobes have different numbers of generations, the sub-structures were aligned. Common structures are indicated on the far left. Probability is color-coded according to legend. Probability of deposition is depicted (a) unweighted, and (b) weighted by relative abundance in the exhalation plume (mode 1–5: 72.8%, 21.0%, 2.2%, 3.4%, 0.6% [12,14]).
Figure 6. Heat map for the probability of ARAP deposition in the thoracic region based on total deposited particle number (extrathoracic, thoracic). Shown are all thoracic sub-structures and all ARAP modes. As lobes have different numbers of generations, the sub-structures were aligned. Common structures are indicated on the far left. Probability is color-coded according to legend. Probability of deposition is depicted (a) unweighted, and (b) weighted by relative abundance in the exhalation plume (mode 1–5: 72.8%, 21.0%, 2.2%, 3.4%, 0.6% [12,14]).
Jpm 11 00431 g006aJpm 11 00431 g006b
Table 1. Mucus definition from composition data.
Table 1. Mucus definition from composition data.
ComponentFraction Per Weight (%)Mean Fraction Per Weight (%) (a)Density (g/cm3)
Protein (mucin)2–53.51.35
Carbohydrate (glycan)7.5–98.31.5
Lipid1–21.50.985
Ions111.409
H2O80–90 (b)85.71
Mucus composition 1.041
(a) mean fraction was used to calculate the density of ERAPs, (b) content decreases from airways in the lower to the upper respiratory tract [46].
Table 2. Ambient air respiratory aerosol particle characteristics used for in silico deposition simulation.
Table 2. Ambient air respiratory aerosol particle characteristics used for in silico deposition simulation.
ModeDiameter (µm)Anatomical Origin [13]Mass (g)Volume (µm3)Density (g/cm3)
Mode 1—breathing, vocalization0.8Bronchiolar fluid film burst3.57 × 10−130.2691.328
Mode 2—vocalization1.8Laryngeal fluid film burst4.057 × 10−123.0541.328
Mode 3—vocalization3.5Laryngeal fluid film burst2.982 × 10−1122.451.328
Mode 4—vocalization5.5Laryngeal fluid film burst1.157 × 10−1087.1151.328
Mode 5—vocalization72.5Oral cavity2.6507 × 10−7199,532.0401.328
Table 3. MPPD parameter settings for breathing scenario “light exercise”.
Table 3. MPPD parameter settings for breathing scenario “light exercise”.
Input SectionScenarioParameterValue Setting
Airway MorphometryAerosolModelYeh/Schum 5-Lobe
Inhalant PropertiesConstant ExposureFRC3300 mL
Exposure Condition URT50 mL
Density1.328 g/cm3
Aspect Ratio1.0 (=spherical)
Diameter0.8, 1.8, 3.5, 5.5, 72.5 µm (a)
Body OrientationUpright
Aerosol Concentration0.5 mg/m3 (b)
Breathing Frequency15 per minute
Tidal Volume750 mL
Inspiratory Fraction0.5
Pause Fraction0
Breathing ScenarioOronasal-Normal Augmenter
Deposition/ClearanceDeposition Only
(a) One simulation per aerosol particle mode. (b) concentration used (only relevant for calculating deposition probabilities and not intended to reflect real world exposure). FRC, functional residual capacity; URT, upper respiratory tract.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hofer, S.; Hofstätter, N.; Duschl, A.; Himly, M. SARS-CoV-2-Laden Respiratory Aerosol Deposition in the Lung Alveolar-Interstitial Region Is a Potential Risk Factor for Severe Disease: A Modeling Study. J. Pers. Med. 2021, 11, 431. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11050431

AMA Style

Hofer S, Hofstätter N, Duschl A, Himly M. SARS-CoV-2-Laden Respiratory Aerosol Deposition in the Lung Alveolar-Interstitial Region Is a Potential Risk Factor for Severe Disease: A Modeling Study. Journal of Personalized Medicine. 2021; 11(5):431. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11050431

Chicago/Turabian Style

Hofer, Sabine, Norbert Hofstätter, Albert Duschl, and Martin Himly. 2021. "SARS-CoV-2-Laden Respiratory Aerosol Deposition in the Lung Alveolar-Interstitial Region Is a Potential Risk Factor for Severe Disease: A Modeling Study" Journal of Personalized Medicine 11, no. 5: 431. https://0-doi-org.brum.beds.ac.uk/10.3390/jpm11050431

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop