Slight reduction in SARS-CoV-2 exposure viral load due to masking results in a significant reduction in transmission with widespread implementation =================================================================================================================================================== * Ashish Goyal * Daniel B. Reeves * E. Fabian Cardozo-Ojeda * Bryan T. Mayer * Joshua T. Schiffer ## Abstract Masks are a vital tool for limiting SARS-CoV-2 spread in the population. Here we utilize a mathematical model to assess the impact of masking on transmission within individual transmission pairs and at the population level. Our model quantitatively links mask efficacy to reductions in viral load and subsequent transmission risk. Our results reinforce that the use of masks by both a potential transmitter and exposed person substantially reduces the probability of successful transmission, even if masks only lower exposure viral load by ∼50%. Slight increases in mask adherence and/or efficacy above current levels would reduce the effective reproductive number (Re) substantially below 1, particularly if implemented comprehensively in potential super-spreader environments. Our model predicts that moderately efficacious masks will lower exposure viral load 10-fold among people who get infected despite masking, potentially limiting infection severity. Because peak viral load tends to occur pre-symptomatically, we also identify that antiviral therapy targeting symptomatic individuals is unlikely to impact transmission risk. Instead, antiviral therapy would only lower Re if dosed as post-exposure prophylaxis and if given to ∼50% of newly infected people within 3 days of an exposure. These results highlight the primacy of masking relative to other biomedical interventions under consideration for limiting the extent of the COVID-19 pandemic prior to widespread implementation of a vaccine. ## Introduction Masks are a barrier method to prevent the spread of respiratory viral infections. A mask essentially serves as a filter that prevents passage of some portion of viruses from the airway of the transmitter to the airway of exposed contacts. Mask efficacy is therefore mediated by a reduction in exposure viral load (*1, 2*). If a potential transmitter as well as an exposed contact are masked, then this filtering process occurs twice potentially amplifying protection. Mask efficacy is inferior to that of other, less permeable barrier methods to prevent infections such as condoms (*1, 3, 4*). The most commonly used cloth and hospital masks do not provide a perfect facial seal and mask fabric does not block emission or inhalation of all aerosolized viral particles (*5, 6*). N95 masks may bypass these shortcomings but are in short supply and difficult to wear for long periods of time (*7, 8*). As a result of these imperfections, recommendations for mask use have varied over the course of the COVID-19 pandemic. Nevertheless, widespread mask use is recognized as a critical component of any viable public health strategy against COVID-19 (*9-11*). Recent models demonstrate that slight increases in mask utilization could be the single most important factor that prevents exponential growth in incident cases (*12-14*). Quantifying mask efficacy in real-world settings remains challenging. Elegant experimental work demonstrated the efficacy of masks in animal models (*15*). Many studies have been performed in hospital settings where mask compliance is uniform and other complementary infection prevention methods are more commonly employed than in other public gathering or work locations (*16, 17*). To the best of our knowledge, no study has captured the impact of masking on the likelihood of super-spreader events, with specific consideration of intermittent compliance. Here we develop a mathematical model capturing viral load-mediated effects of mask use on transmission probability within transmission pairs and at the population level. We use this approach to estimate the efficacy of masks in real world settings, and to characterize effects on super-spreader events as well as exposure viral loads of those who get infected despite masking. Finally, we compare the preventative impact of masking to the use of antiviral therapies given early during symptomatic infection, or when used as post-exposure prophylaxis (PEP). ## Results ### Baseline mathematical model of SARS-CoV-2 viral load dependent transmission We employed a previously developed mathematical model which links viral load shedding at the individual level with population level epidemic spread (*18, 19*), to determine the impact of masks on epidemics. Briefly, the model is built upon the assumptions that each transmitter has a specific number of exposure contacts per day and that each exposure contact has a certain probability of successful transmission based on the viral load of the transmitter. This probability is based on a transmission dose (TD) response curve. We fit the model to frequency histograms describing heterogeneity of individual R0, or the number of secondary infections attributed to each infected person (*20-24*). The individual R0 distribution for SARS-CoV-2 transmission is highly over-dispersed, meaning that most infected people do not spread infection while a minority infect a large number of people. We also fit the model to distributions of individual serial intervals, the time from the onset of symptoms in the transmitter to symptom onset in the secondarily infected person (*25*). The overall qualitative conclusions of this model were that the period of contagiousness for SARS-CoV-2 is quite short, typically less than a day, and that super-spreader events are largely attributable to high variability in the number of exposure contacts per day among infected people. ### Predicted impact of transmitter or exposure contact masking on transmission probability within transmission pairs We added masking to this model by assuming that a mask decreases the exposure viral load in a transmission pair by a value that we refer to as the combination mask efficacy (*ε**C*). This efficacy represents the proportion of viruses filtered by masks worn by both the transmitter and exposed person. If the transmitter is wearing mask with efficacy *ε**T* and the exposed person is wearing a mask with efficacy *ε**E*, then the exposure viral load *V**E* can be related to the transmitter viral load *V**T* by: *V**E* = *V**T* (*1-ε**T*)(*1-ε**E*). The combination mask efficacy is then *ε**C*=*1-*(*1-ε**T*)(*1-ε**E*), which takes on a value of zero when both parties are not wearing a mask or wearing masks that are totally ineffective **(Fig 1)**. ![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F1) Figure 1. Mathematical model of mask impact on SARS-CoV-2 exposure viral load **a**. The viral load emitted by a potential transmitter (Vt) can be filtered, resulting in lower exposure viral loads due to a single mask worn by a transmitter or exposed individual with efficacy eT or eE respectively. Dual masking lowers exposure viral load further by filtering virus twice. **b**. Dual masking may prevent super-spreader events to a greater extent than a masked transmitter or masked exposed individual. As in our prior model, the exposure viral load impacts contagiousness, which is the probability that virus is passaged to the exposed person’s airway, as well as infectiousness, the probability of cellular infection given the presence of virus in the airway. Each of these properties is associated with a dose response curve (contagiousness dose (CD) response curve and infectiousness dose (ID) response curve), the product of which is the TD response curve. In this context, we first establish a baseline probability of transmission given no use of masks on both sides of a potential transmission pair **(Fig 2a)**. The absolute & relative reductions in transmission probability with more effective masks vary as *V**T* increases. At lower viral loads (<107.5 viral RNA copies), a moderate to highly effective mask worn by either a transmitter (0.9>*ε**T*≥0.5) or an exposed contact (0.9>*ε**E*≥0.5) is sufficient to partially lower the absolute probability of transmission **(Fig 2b-c, Sup fig 1a)**. The relative reduction in transmission probability increases linearly with increasing mask efficacy at 107 viral RNA copies with an increasingly concave, curvilinear relationship at higher viral loads **(Sup fig 1c)**. ![Supplementary figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F11.medium.gif) [Supplementary figure 1.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F11) Supplementary figure 1. Impact of masking of the transmitter or exposed contact alone or both members of the transmission pair, on transmission risk given an exposure viral load. Each panel is based on simulations of 1000 transmission pairs. **a-b**. Absolute reduction in transmission risk. **c-d**. Relative reduction in transmission risk. **a**,**c**. Transmitter is masked only **b, d**. Both transmitter and exposed contact are masked ![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F2) Figure 2. Impact of masking of the transmitter alone, the exposed contact alone or both members of the transmission pair, on transmission risk given a single exposure contact. **a-d**. Each panel is based on simulations of 1000 transmission pairs. **a**. No masking, **b**. Transmitter is masked, **c**. Exposed contact is masked, **d**. Both members are masked. At higher transmitter viral load (107.5–109 viral RNA copies), a moderate to highly effective mask worn by either a transmitter (0.9>*ε**T*≥0.5) or an exposed contact (0.9>*ε**E*≥0.5) insignificantly lowers the absolute probability of transmission **(Fig 2b-c, Sup fig 1a)**. At high viral loads, the relative reduction in transmission probability increases dramatically with extremely effective masks of efficacy ≥0.9, when the mask is worn by either a transmitter or an exposed contact **(Sup fig 1c)**. ### Predicted impact of dual masking on transmission probability within transmission pairs If both the transmitter and exposed person wear masks (dual masking), then lower mask efficacies are sufficient to significantly lower transmission risk at a wider range of exposure viral loads **(Fig 2d)**. At viral loads <108 viral RNA copies, masks worn by both transmitter and exposed contact of more than moderate efficacy (*ε**T*≥0.5, *ε**E*≥0.5 resulting in *ε**C*≥0.75), is sufficient to partially lower the absolute probability of transmission **(Fig 2d, Fig 3, Sup fig 1b)**. The relative reduction in transmission probability according to mask efficacy increases more rapidly with dual **(Sup fig 1d)** compared to single **(Sup fig 1c)** masking. If both transmitter and exposed contacts wear masks with *ε**T*= 0.9 and *ε**E* = 0.9 (*ε**C*=0.99), then transmission probability is reduced to <5% for viral loads <108.5 viral RNA copies and to ∼20% for transmitter viral load of 109 viral RNA copies **(Fig 2d, 3d & Sup fig 1c)**. ![Figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F3) Figure 3. Impact of masking of both members of a transmission pair on transmission risk according to transmitter viral load. **a-d**. Each panel is based on simulations of 1000 transmission pairs. **a**. Low viral load, **b**. Moderate viral load, **c**. High viral load, **d**. Extremely high viral load. ### Predicted impact of transmitter and exposure contact masking on effective reproduction number (Re) at different levels of implementation We next explore the impact of general masking adherence rates on population level metrics of infection by simulating 3000 potential transmitters assuming heterogeneity in viral load trajectories and exposure contact networks among individuals. Reduction in the effective reproductive number (Re) depends on both the mask efficacy levels (*ε**T* and *ε**E*) and the level of adherence to masking **(Fig 4)**. If we assume 25% of people wear masks 25% of the time, then most variability in results occur due to the stochastic nature of the model **(Fig 4a)**. If we assume that 50% of people wear masks 50% of the time, then the use of masks with high efficacy of ∼0.9 results in a drop of Re from ∼1.8 to ∼1.0 **(Fig 4b)**. With 75% of people wearing masks 75% of the time, a mask efficacy of ∼0.5 allows for a reduction of Re from ∼1.8 to ∼1.0 **(Fig 4c)**. With 100% of people wearing masks 100% of the time, then a mask efficacy of ∼0.3 is sufficient to achieve Re ∼1.0, and efficacy of 0.5 in both transmitter and exposed contacts lowers Re to less than 0.6 **(Fig 4d)**. ![Figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F4.medium.gif) [Figure 4.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F4) Figure 4. Effect of mask utilization and efficacy on mean SARS-CoV-2 R0. Each heat map is based on simulations of 3000 transmitters with varying daily exposure contacts. **a**. Low mask utilization, **b**. Moderate mask utilization, **c**. High mask utilization, **d**. Extremely high mask utilization. Current estimates of daily mask use over the last 2 months vary between states in the United States between 40 and 60% ([http://www.healthdata.org/acting-data/maps-mask-use](http://www.healthdata.org/acting-data/maps-mask-use)), and Re has varied between 0.8 and 1.2 ([https://covid19-projections.com/us](https://covid19-projections.com/us)). These results suggest that panel **Fig 4c** is likely the closest to current U.S. epidemic conditions and that εT and εE likely fall roughly between 0.4 and 0.6 in a real-world setting, if mask efficacy is equivalent between transmitter and exposed contacts. While a combined efficacy *ε**C* of 0.65-0.85 can be roughly estimated from the model, the possibility of superior efficacy of masks in transmitters versus exposed, or vice versa, cannot be excluded. ### Proportions of infections attributable to masked and unmasked transmission pairs We next project the proportion of transmission events attributable to different masking profiles among transmission pairs assuming equally effective mask (*ε*) used by both transmitters and exposed contacts. In circumstances with low mask utilization (25% of people wearing masks 25% of the time), nearly all transmissions occur from an unmasked person to an unmasked person **(Fig 5a)**. A similar trend is noted for moderate mask utilization (50%), particularly as mask efficacy increases **(Fig 5b)**. ![Figure 5.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F5.medium.gif) [Figure 5.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F5) Figure 5. Effect of mask utilization and efficacy on proportion of masked transmissions contributing to total R0. Each histogram bar in each panel is based on simulations of 3000 transmitters with varying daily exposure contacts. **a**. Low mask utilization, **b**. Moderate mask utilization, **c**. High mask utilization, **d**. Extremely high mask utilization. For high (75%) and extremely high (90%) mask utilization scenarios, if mask efficacy is moderate (*ε*∼0.6) as is currently believed, then a higher proportion of transmissions occur to or from a person wearing a mask, despite the fact that the total number of transmissions dramatically decreases **(Fig 5c, d)**. We project that at current epidemic conditions in the United States, including ∼50% mask use (similar to 75% of people wearing masks 75% of the time) and mask efficacy of 0.3-0.6, that fewer than half of ongoing transmissions occur within unmasked pairs and that transmission between masked transmitters or masked exposed contacts likely contributes significantly to Re **(Fig 5c)**. ### Predicted impact of transmitter and exposure contact masking on super-spreader events We next identified that increased mask compliance and efficacy dramatically decreases the proportion of infected people who successfully transmit to another person **(Fig 6a-d)**. If 50% of people were to wear 50% effective masks half of the time, then the likelihood of an individual transmitting decreases from 30% to 20% **(Fig 6b)**. If mask compliance is increased to 75% of people 75% of the time, then the likelihood of a person transmitting decreases to ∼15% **(Fig 6c)**. ![Figure 6.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F6.medium.gif) [Figure 6.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F6) Figure 6. Effect of mask utilization and efficacy on percent of infected people who transmit to others. Each heat map is based on simulations of 3000 transmitters with varying daily exposure contacts. **a**. Low mask utilization, **b**. Moderate mask utilization, **c**. High mask utilization, **d**. Extremely high mask utilization. When masking is applied homogeneously across the population, the proportion of infectors (transmitters who infect at least one person) who pass the infection to 5 or more people decreases, as mask utilization and efficacy increase **(Sup fig 2a-d)**. Increased mask utilization and increased mask efficacy leads to an even reduction of all types of transmission events, including transmissions to small numbers (1-3) of people, or super-spreader events to >5, >10, >20 or >50 people **(Fig 7a-d)**. Improvements in mask efficacy have a larger impact as utilization of mask use increases, including against super-spreader events **(Fig 7c, d)**. Under all simulations, super-spreader events with transmission to >5 people persist and make a nearly equivalent contribution proportionally to overall Re, though their absolute impact is considerably lessened with higher mask compliance and efficacy. ![Supplementary figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F12.medium.gif) [Supplementary figure 2.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F12) Supplementary figure 2. Effect of mask utilization and efficacy on percent of infectors who supers-spread to others. Each heat map is based on simulations of 3000 transmitters with varying daily exposure contacts. **a**. Low mask utilization, **b**. Moderate mask utilization, **c**. High mask utilization, **d**. Extremely high mask utilization. ![Figure 7.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F7.medium.gif) [Figure 7.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F7) Figure 7. Effect of mask utilization and efficacy on proportion of super-spreader transmissions. Each histogram bar in each panel is based on simulations of 3000 transmitters with varying daily exposure contacts. **a**. Low mask utilization, **b**. Moderate mask utilization, **c**. High mask utilization, **d**. Extremely high mask utilization. Our results suggest that with current levels of masking in the United States **(Fig 7c**, *ε*= 0.3-0.6**)**, most of the contribution to Re still comes from super-spreader events involving >5 secondary infections. We therefore simulated masking applied to 100% of people with >10 exposure contacts per day **(Sup fig 3a-d)** and found that even modest uptake of moderately effective masks (∼0.5) in the remainder of the population appeared to be sufficient to maintain Re <1. ![Supplementary figure 3.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F13.medium.gif) [Supplementary figure 3.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F13) Supplementary figure 3. Effect of mask utilization and efficacy on mean SARS-CoV-2 R0 assuming all persons with >10 exposure contacts per day are masked. Each heat map is based on simulations of 3000 transmitters with varying daily exposure contacts. **a**. Low mask utilization, **b**. Moderate mask utilization, **c**. High mask utilization, **d**. Extremely high mask utilization. ### Predicted impact of transmitter and exposure contact masking on viral inoculum at time of infection Another theoretical benefit of masks is reduction in exposure viral load which in animal models of SARS-CoV-1 and MERS, leads to less severe infection (*26-28*). Simulating under the assumption that 75% of people wear masks 75% of the time (i.e., a situation representing current levels of masking in the United States), we identified that viral load of a transmitter required to generate secondary infections increases slightly with higher implementation of more efficacious masks, particularly with dual masking of transmitters and exposed contacts **(Fig 8a)**. Viral load exposure at the time of a successful transmission decreased according to efficacy of mask, particularly if both transmitters and exposed contacts are masked **(Fig 8b)**. With dual masking in place with efficacies of 0.6, exposure viral load decreased by ∼1 log **(Fig 8b)**. With dual masking in place with efficacies of 0.9, exposure viral load decreased by ∼2 logs **(Fig 8b)**. ![Figure 8.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F8.medium.gif) [Figure 8.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F8) Figure 8. Effect of mask efficacy on exposure viral load during successful transmissions. Each boxplot in each panel is based on simulations of 3000 transmitters with varying daily exposure contacts. Boxplots are median and interquartile range (IQR) and lines are 1.5 the IQR. **a**. Viral load of transmitters, **b**. Viral load exposure in secondarily infected persons. ### Predicted impact of antiviral therapy during early symptomatic infection on Re Treatment as prevention is a highly effective means for reducing person-to-person transmission of HIV (*29, 30*). Our models predict that initiation of potent antiviral therapy early after symptom onset is likely to have therapeutic benefit (*18*). We therefore tested whether early symptomatic therapy which rapidly eliminates shedding might also decrease secondary transmissions. Owing to the fact that symptomatic therapy would almost invariably occur after peak viral shedding **(Sup Fig 4a)**, our simulations suggest that even widespread implementation of early symptomatic therapy would not lower Re **(Fig 9a)**, percentage of infected people who transmit to at least one other person **(Fig 9c)** or percentage of infectors who transmit to at least 5 other people **(Fig 9e)**. ![Supplementary figure 4.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F14.medium.gif) [Supplementary figure 4.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F14) Supplementary figure 4. Viral load simulations on antiviral therapy. **a-b.**Solid lines trajectories are assumed with therapy, whereas dashed lines are counterfactuals with no therapy. **b**. Simulated early symptomatic therapy. **b**. Simulated PEP. ![Figure 9.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F9.medium.gif) [Figure 9.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F9) Figure 9. Projected impact of antiviral therapy on SARS-CoV-2 R0. Each line in each panel is based on simulations of 3000 transmitters with varying daily exposure contacts. Panels a, c and e are for early symptomatic therapy. Each simulation assumes 100% efficacy. Panels b, d and f are for post-exposure prophylaxis given during pre-symptomatic infection. Colored lines assume different antiviral efficacies (inset). **a, b**. Projected R0 given different amounts of uptake in the population. **c, d**. Projected percent of transmitters who infect at least one person given different amounts of uptake in the population. **e, f**. Projected percent of infectors who infect at least five people given different amounts of uptake in the population. ### Predicted impact of PEP on Re PEP is also a potential method for lowering SARS-CoV-2 transmissions. Because PEP is given in the pre-symptomatic phase and would usually fall before or near peak viral shedding **(Sup Fig 4b)**, our simulations suggest an inverse linear relationship between uptake of PEP and Re **(Fig 9b)**, percentage of infected people who transmit to at least one other person **(Fig 9d)** or percentage of infectors who transmit to at least 5 other people **(Fig 9f)**. To achieve Re <1 would require PEP efficacy of 50% and ∼75% uptake in the population, which would in turn require 75% of SARS-CoV-2 cases to be contact traced. Increases in PEP efficacy beyond 0.5 would provide minimal to no enhancement of this benefit **(Fig 9b)**. We further determine that the timing of PEP is critical. If PEP is initiated within 2 days of exposure, then the percent of people receiving effective treatment is highly predictive of Re **(Fig 10)**. However, from day 3 onwards, the impact of effective PEP diminishes **(Fig 10)**. ![Figure 10.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2020/10/05/2020.09.13.20193508/F10.medium.gif) [Figure 10.](http://medrxiv.org/content/early/2020/10/05/2020.09.13.20193508/F10) Figure 10. Projected impact of post-exposure prophylaxis on SARS-CoV-2 according to the day of implementation. Each line in is based on simulations of 3000 transmitters with varying daily exposure contacts and are colored by percent of transmitters receiving therapy. Overall, these results highlight the fact that masking is likely to have more of an impact on Re than any form of licensed antiviral therapy that emerges during the course of the pandemic. ## Discussion Relative to other barrier methods for preventing the spread of infectious diseases such as condoms, masks are imperfect (*31*). Surgical and cloth masks, which are now used commonly by members of the public, do not completely eliminate droplet and airborne emission of viral particles by a transmitter (*1*). Nor do they prevent viral exposure to airway cells among exposed contacts. Their effectiveness in real world settings is further limited by intermittent compliance and improper masking technique. As a result, masks only prevent a proportion of person-to-person transmissions. Nevertheless, our results demonstrate that in the absence of a licensed vaccine, based on moderate efficacy, low cost, high availability, and ease of use, masks are the most effective currently available biomedical intervention. If implemented widely and strategically, on top of baseline levels of physical distancing, masking would be sufficient to suppress ongoing spread of SARS-CoV-2 until widespread deployment of a vaccine is possible. More specifically, our model suggests that increased masking would lower the effective reproduction number (Re), lower the percentage of infected people who transmit the virus, decrease the total number of super-spreader events, and lower the exposure viral loads among infected people, possibly leading to less severe infections overall (*26*). Importantly, there appears to be a critical threshold of compliance. We predict massive additional benefits accrued with an increase in masking compliance from 75% of people masking 75% of the time to 90% of people masking 90% of the time. Masking also highlights the critical nature of suppressing super-spreader events. If nearly 100% compliance could be achieved among persons with 10 or more exposure contacts per day, then this would be sufficient for maintaining Re less than one. This result highlights that policies mandating the proper use of masks at all times by all persons at sites of known super-spreader events including high risk work environments, locker rooms, weddings, social gatherings, and schools should be considered. Slight increases in mask efficacy could also drive Re to much lower levels. We believe that different types of masks should be comparatively tested with the same scientific rigor applied to clinical trials of small molecular agents and vaccines. An important artifact of widespread masking is that while the total number of incident cases is expected to decrease dramatically, the proportion of transmissions in which at least one member of the transmission pair is masked will be higher. Therefore, anecdotal documentation of successful infections between masked individuals, or even super-spreader events in which many infected people were masked, should not be misinterpreted as failure of masking policy. The counterfactual, that masks limited the severity of these events, is likely to be true. Only longitudinal incidence data, along with shifts in level of masking in a given region, are appropriate for inferring the effectiveness of masks. It is unlikely that antiviral therapy will prove nearly as useful as widespread masking for preventing transmissions. We previously demonstrated that antiviral therapy given during the early symptomatic phase of infection has the potential to limit duration of shedding and infection associated inflammation and is likely to be more efficacious than therapy given later during COVID-19 to hospitalized patients (*18*). Unfortunately, early symptomatic therapy would likely occur after peak viral load in a majority of cases. Our simulations suggest that even 100% penetrance of extremely potent antiviral therapies would have a negligible impact on population level spread of the virus. Therefore, while treatment as a prevention is a vital piece of HIV public health policy (*30*), it is unlikely to impact the COVID-19 pandemic. On the other hand, treatment in the earliest pre-symptomatic phase of infection, which could only realistically occur in the setting of PEP, happens prior to peak viral load and therefore could limit secondary transmissions. However, the gains from this approach diminish with each day following exposure. In order to meaningfully impact Re, over 50% of exposed contacts would need to receive fully effective therapy within 3 days of an exposure. Given that no available agent yet achieves this level of efficacy and that identifying 50% of post-exposure contacts is unrealistic in most countries, it is clear that relative to masking, PEP will only have an adjunctive role in managing the pandemic. Potential areas of implementation are among high risk populations such as skilled nursing facility residents or cancer center patients and among populations where masking is difficult or impossible. Our prior work strongly suggests the presence of a transmission dose response curve in which exposure viral load is a key determinant of transmission risk (*32, 33*). Our current analysis is built upon this assumption. We project that masks will lead to a lower exposure viral load among newly infected people, particularly if both the transmitter and exposed individual are successfully masked. Animal models of SARS-CoV-1 and MERS (*27, 28, 34*), as well as challenge studies with influenza H1N1 in humans (*35*), all demonstrate that lower exposure dose is associated with less severe disease, and human data from the Hong Kong SARS-CoV-1 outbreak in Hong Kong suggest a similar trend (*36*). While data for SARS-CoV-2 in humans is lacking, it is notable that age-adjusted hospitalization and death rates may be decreasing since more widespread utilization of masks. A mask related reduction in exposure viral load is a plausible but unproven reason for this observation. Our work has key limitations. First, based on available data, it is impossible to know the true average efficacy of a mask worn by a transmitter or an exposed contact. Many epidemics at the state level have demonstrated a reduction in the effective reproductive number ranging from 0.2-0.5 when more widespread masking was implemented, even as physical distancing levels waxed and waned. Our model suggests that if 50-75% of people wear masks 50-75% of the time, which is roughly in accordance with state level observations of mask compliance, then a broad estimate for real world mask efficacy is ∼0.4-0.6, assuming that efficacy is equal between transmitters and exposed contact, and that masks are properly used to optimize their efficacy. If transmitter masking is more efficacious, while exposed contact masking is proportionally less efficacious, then similar results can be expected. The real-world estimate is inclusive of multiple factors including variability in mask type and masking technique. Regardless of the precise estimate, it is clear that wider implementation would yield significant reductions in spread of SARS-CoV-2 at the population level. Second, our generalized model is not region-specific for the current pandemic. The relative impact of super-spreader events, intensity of transmission and proportion of symptomatic cases may vary from region to region based on contact network structure and age demographics. Nevertheless, the general qualitative conclusions about masking are insensitive to these differences and are likely to be generalizable across the globe. Finally, our model does not include a standard SIR format and therefore does not capture other dynamic features that might alter the force of infection such as herd immunity or time-variant shifts in degree of physical distancing. Another missing feature that could be captured with an SIR modeling framework is the possibility of an assortative mixing pattern, in which individuals with lower adherence to masking might preferentially interact with others who have low adherence to masking (*12*). Such an effect could allow persistence of SARS-CoV-2 within this sub-population, even if masking is sufficient for containment in the rest of the population. In conclusion, we developed a mechanistic model to demonstrate how masks reduce exposure SARS-CoV-2 viral load and transmission probability. Widespread use of even modestly effective masks is predicted to severely limit epidemic spread and represents the key available intervention along with physical distancing, to mitigate the number of infections, and perhaps the proportion of infections that are severe, while the world awaits a widely available and effective vaccine. ## Methods ### SARS-CoV-2 within-host model We used the within-host model describing the SARS-CoV-2 infection from our previous study (*18*). This model assumes that the contact of SARS-CoV-2 (*V*) with susceptible cells (*S*) produces infected cells at rate *βVS* which then generates new virus at a per-capita rate *π*. The model also incorporates the death of infected cells mediated by (1) the innate responses (*δI**k*) and (2) the acquired immune responses ![Graphic][1] by SARS-CoV-2-specific effector cells (*E*). The magnitude of the innate immunity is dependent on the infected cell density and the exponent *k*. The nonlinearity of the acquired responses is captured by the Hill coefficient *r* that allows for rapid saturation of the killing. Finally, the parameter *ϕ* defines level of SARS-CoV-2-specific effector cells at which the killing of infected cells becomes half maximal. In the model, the rise of SARS-CoV-2-specific effector cells rise is described in a two-stage manner. The first stage defines the proliferation of the first precursor cell compartment (*M*1) at rate *ωIM*1 and differentiation into a second precursor cell compartment (*M*2) at a per capita rate *q*. Finally, second precursor cells differentiate into effector cells at the same per capita rate *q* and are cleared at rate *δ**E*. The model is expressed as a system of ordinary differential equations: ![Formula][2] The initial conditions for the model were assumed as *S*(0) = 107 cells/mL, *I*(0) = 1 cells/mL, ![Graphic][3] copies/mL, *M*1(0) = 1, *M*2(0) = 0 and *E* = 0. For simulations we sampled parameter values from a nonlinear mixed-effect model as described in (*19*), with the following fixed effects and standard deviation of the random effects (in parenthesis): Log10*β*: - 7.23 (0.2) virions-1 day-1; *δ*: 3.13 (0.02) day-1 cells-k; *k*: 0.08 (0.02); Log10(*π*): 2.59 (0.05) day-1; *m*: 3.21 (0.33) days-1cells-1; Log10(*ω*): −4.55 (0.01) days-1cells-1. We also assumed *r* = 10; *δ**E* = 1 day-1; *q* = 2.4 × 10−5 day-1 and *c* = 15 day-1. ### Dose-response model We employed our previously developed dose-response model to estimate the probability of virus entering the airway given a transmitter viral load (i.e., contagiousness) and the probability of cellular infection given a transmitter viral load, (i.e., infectiousness) *P**t*[*V*(*t*)] (response) based on viral loads *V*(*t*) (dose) (*19*). The relation between the response and the dose follows, ![Graphic][4], being *λ* the viral load that corresponds to 50% infectiousness and 50% contagiousness and *α* the Hill coefficient that controls the sharpness in the dose-response curve. We assumed that the viral load-dependent contagiousness (i.e., the probability that virus is passaged to the exposed person’s airway) is the same as infectiousness. We estimate the transmission risk as the product of the infectiousness and contagiousness (*19*). ### Transmission model and reproduction number As in our previous model (*19*), we determined the total exposed contacts of a transmitter within a time step (Δ*t*) using a gamma distribution, i.e. ![Graphic][5], where *θ* and *ρ* represent the average daily contact rate and the dispersion parameter, respectively. The true number of exposure contacts (with viral airway exposure) was then obtained by multiplying the total exposed contacts and the contagiousness of the transmitter (i.e., ![Graphic][6]). We modelled infectiousness as a Bernoulli event with mean *P**t*, yielding the number of secondary infections within a time step as ![Graphic][7]. Finally, we summed up the number of secondary infections over 30 days since the time of exposure to obtain the individual reproduction number, i.e. ![Graphic][8]. For each successful transmission, we further assumed that it takes *τ* days for the first infected cell to produce virus. In simple steps, we followed the procedure below to estimate Re, 1. Simulate viral load *V*(*t*) of a simulated infected individual using **the within-host model**. 2. For a given combination of (*λ, τ, α, θ, ρ*) 1. For each time step Δ*t* 1. Compute *P**t*[*V*(*t*); *λ, α*] 2. Draw ![Graphic][9] 3. Calculate ![Graphic][10] 2. Calculate ![Graphic][11] 3. Repeat Steps 1 and 2 to estimate *R**e* for 3,000 infected individuals. The population level *R* can then be calculating by taking the mean of 3,000 individual *R**e* values. ### Parameter values for the transmission model For simulations, we used the parameter set [*α, λ, τ, θ, ρ*] = [0.8, 107, 0.5, 4, 40]) as they most closely reproduces empirically observed individual *R* and serial interval histograms as well as mean *R* across individuals (*R* ∈ [1.4, 2.5]) and mean serial interval across individuals (SI ∈ [4.0, 4.5]) early during the pandemic (*20-22, 24, 25*). ### Modeling mask use To evaluate the impact of the use of mask on epidemics, we first assumed that a mask decreases the exposure viral load by a fraction (1 − *ϵ*), being *ϵ* the mask efficacy or the proportion of viruses filtered by the mask of transmitter or exposed individuals. If the transmitter is wearing mask with efficacy *ϵ**T* and the exposed person is wearing a mask with efficacy *ϵ**E*, then the combined mask efficacy *ϵ**C*, is given by 1 − (1 − *ϵ**T*)(1 − *ϵ**E*). Infectiousness or contagiousness reduction by the use of mask can be computed as: ![Formula][12] Similarly, the transmission risk reduction using mask can be computed as: ![Formula][13] Finally, we modeled the *compliance* of an individual wearing mask in the population as a Bernoulli event with mean *σ* and the *adherence* of wearing it at time step Δ*t* as a Bernoulli event with mean *ϑ*. Compliance is defined as whether the person ever wears a mask. Adherence is the percentage of time that a mask wearer wears a mask. ### Simulating secondary transmissions with mask use For a specific scenario with selected *σ* and *ϑ*, we followed the procedure below to estimate the population level *R*: 1. Simulate *V*(*t*) for a transmitter using the within-host model in **eq. 1**. 2. Simulate transmitter mask compliance using *Ber*(*p* = *σ*). 3. Discretize the time-space of 30 days over time steps Δ*t*. For each time step, 1. If the transmitter is *compliant* in using a mask: 1. Simulate transmitter mask *adherence* at time step Δ*t* using *Ber*(*p* = *ϑ*) 2. If transmitter is wearing a mask at time step Δ*t* 1. Draw ![Graphic][14]. 2. Determine masking adherence among exposed contacts at time step Δ*t* using *k*∼*Ber*(*p* = *σϑ*). If *k* = 1, then there is 100% adherence among exposed contacts and if *k* = 0, then there is 0% adherence among exposed contacts. 3. Determine the number of exposed contact wearing a mask (i.e., ![Graphic][15]) and the number of exposed contact not wearing a mask (i.e., ![Graphic][16]. 1. Compute ![Graphic][17] 2. Calculate ![Graphic][18] 3. Compute ![Graphic][19] 4. Calculate ![Graphic][20] 3. If the transmitter is not adhering with the use of a mask at time step Δ*t*, which is determined at step (i). 1. Draw ![Graphic][21]. 2. Determine masking adherence among exposed contacts at time step Δ*t* using *k*∼*Ber*(*p* = *σϑ*). 3. Determine the number of exposed contact wearing a mask (i.e., ![Graphic][22]) and the number of exposed contact not wearing a mask (i.e., ![Graphic][23] 1. Compute ![Graphic][24] 2. Calculate ![Graphic][25] 3. Compute ![Graphic][26] 4. Calculate ![Graphic][27] 2. If the transmitter is not compliant in wearing a mask 1. Draw ![Graphic][28]. 2. Determine masking adherence among exposed contacts at time step Δ*t* using *k*∼*Ber*(*p* = *σϑ*). 3. Determine the number of exposed contact wearing a mask (i.e., ![Graphic][29]) and the number of exposed contact not wearing a mask (i.e., ![Graphic][30] 1. Compute ![Graphic][31] 2. Calculate ![Graphic][32] 3. Compute ![Graphic][33] 4. Calculate ![Graphic][34] 3. Calculate ![Graphic][35]. 4. Calculate ![Graphic][36]. 5. Repeat Steps **1** to **4** to estimate *Q*5 for 3,000 infected individuals (transmitters). Re can then be calculating by taking the mean of 3,000 individual *Q*5 values. ### Modeling antiviral treatment We simulate the antiviral treatment by assuming that the antiviral treatment reduces the viral production (*π*) by (1 − *ϵ**treat*), where *ϵ**treat* is the efficacy of treatment. Here, *ϵ**treat* = 0 and *ϵ**treat* = 1 represent the case of completely ineffective and 100% effective treatment, respectively. For the transmission simulations we model *coverage* to treatment with as a Bernoulli event with mean *ψ*. In the presence of treatment with mean coverage *ψ* and efficacy *ϵ**treat*, we follow the procedure below to estimate the population level *R**e*: 1. Determine adherence to treatment using *Ber*(*p* = *ψ*) 2. Determine the time of start of antiviral treatment (*T**treat*) for treatment in symptomatic phase and pre-symptomatic phase by randomly drawing a number from Uniform distributions (*U*(0.5 + *I**I*, 5*_I**I*)) and (*U*(0.5,5)), respectively, where *I**I* is the incubation period of the infected individual. 3. Simulate viral load *V*(*t*) of a simulated infected individual using the within-host model in **eq. 1** with *ϵ**treat* = 0 for *t* ≤ *T**treat* and 0 < *ϵ**treat* ≤ 1 for *t* > *T**treat*. 4. For a given combination of (*λ, τ, α, θ, ρ*) 1. For each time step Δ*t* 1. Compute *P**t*[*V*(*t*); *λ, α*] 2. Draw ![Graphic][37] 3. Calculate ![Graphic][38] 2. Calculate ![Graphic][39] 5. Repeat Steps 1 and 4 to estimate *R**e* for 3,000 infected individuals. The population level *R**e* can then be calculating by taking the mean of 3,000 individual *R**e* values. ## Data Availability All modeled data and code will be available upon request. ## Footnotes * We added a couple of references * Received September 13, 2020. * Revision received October 3, 2020. * Accepted October 3, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1. M. Liang et al., Efficacy of face mask in preventing respiratory virus transmission: A systematic review and meta-analysis. Travel Med Infect Dis, 101751 (2020). 2. 2. E. P. Fischer et al., Low-cost measurement of face mask efficacy for filtering expelled droplets during speech. Science Advances 6, eabd3083 (2020). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo4OiJhZHZhbmNlcyI7czo1OiJyZXNpZCI7czoxMzoiNi8zNi9lYWJkMzA4MyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEwLzA1LzIwMjAuMDkuMTMuMjAxOTM1MDguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 3. 3. R. Chou, T. Dana, R. Jungbauer, C. Weeks, M. S. McDonagh, Masks for Prevention of Respiratory Virus Infections, Including SARS-CoV-2, in Health Care and Community Settings: A Living Rapid Review. Ann Intern Med, (2020). 4. 4. A. S. Magaret et al., Effect of Condom Use on Per-act HSV-2 Transmission Risk in HIV-1, HSV-2-discordant Couples. Clin Infect Dis 62, 456–461 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/civ908&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26578538&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 5. 5. A. Konda et al., Aerosol Filtration Efficiency of Common Fabrics Used in Respiratory Cloth Masks. ACS Nano 14, 6339–6347 (2020). 6. 6. E. E. Sickbert-Bennett et al., Filtration Efficiency of Hospital Face Mask Alternatives Available for Use During the COVID-19 Pandemic. JAMA Intern Med, (2020). 7. 7. N. H. L. Leung et al., Respiratory virus shedding in exhaled breath and efficacy of face masks. Nat Med 26, 676–680 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0843-2&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32371934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 8. 8. J. Howard et al., Face Masks Against COVID-19: An Evidence Review. Preprints 2020040203 (doi: 10.20944/preprints202004.0203.v3). (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.20944/preprints202004.0203.v3&link_type=DOI) 9. 9. R. Zhang, Y. Li, A. L. Zhang, Y. Wang, M. J. Molina, Identifying airborne transmission as the dominant route for the spread of COVID-19. Proc Natl Acad Sci U S A 117, 14857–14863 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzI2LzE0ODU3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMDUvMjAyMC4wOS4xMy4yMDE5MzUwOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 10. 10. Y. Wang et al., Reduction of secondary transmission of SARS-CoV-2 in households by face mask use, disinfection and social distancing: a cohort study in Beijing, China. BMJ Glob Health 5, (2020). 11. 11. W. Lyu, G. L. Wehby, Community Use Of Face Masks And COVID-19: Evidence From A Natural Experiment Of State Mandates In The US. Health Aff (Millwood) 39, 1419–1425 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1377/hlthaff.2020.00818&link_type=DOI) 12. 12. D. N. Fisman, A. L. Greer, A. R. Tuite, Bidirectional impact of imperfect mask use on reproduction number of COVID-19: A next generation matrix approach. Infect Dis Model 5, 405–408 (2020). 13. 13. C. J. Worby, H. H. Chang, Face mask use in the general population and optimal resource allocation during the COVID-19 pandemic. Nat Commun 11, 4049 (2020). 14. 14. R. O. J. H. Stutt, R. Retkute, M. Bradley, C. A. Gilligan, J. Colvin, A modelling framework to assess the likely effectiveness of facemasks in combination with ‘lock-down’ in managing the COVID-19 pandemic. Proc Math Phys Eng Sci 476, 20200376 (2020). 15. 15. J. F. Chan et al., Surgical mask partition reduces the risk of non-contact transmission in a golden Syrian hamster model for Coronavirus Disease 2019 (COVID-19). Clin Infect Dis, (2020). 16. 16. X. Wang, E. G. Ferro, G. Zhou, D. Hashimoto, D. L. Bhatt, Association Between Universal Masking in a Health Care System and SARS-CoV-2 Positivity Among Health Care Workers. JAMA, (2020). 17. 17. J. J. Bartoszko, M. A. M. Farooqi, W. Alhazzani, M. Loeb, Medical masks vs N95 respirators for preventing COVID-19 in healthcare workers: A systematic review and meta-analysis of randomized trials. Influenza Other Respir Viruses 14, 365–373 (2020). 18. 18. A. Goyal, E. F. Cardozo-Ojeda, J. T. Schiffer, Potency and timing of antiviral therapy as determinants of duration of SARS CoV-2 shedding and intensity of inflammatory response. medRxiv, 2020.2004.2010.20061325 (2020). 19. 19. A. Goyal, D. B. Reeves, E. F. Cardozo-Ojeda, J. T. Schiffer, B. T. Mayer, Wrong person, place and time: viral load and contact network structure predict SARS-CoV-2 transmission and super-spreading events. medRxiv, 2020.2008.2007.20169920 (2020). 20. 20. A. Endo, S. Abbott, A. J. Kucharski, S. Funk, C. f. t. M. M. o. I. D. C.-W. Group, Estimating the overdispersion in COVID-19 transmission using outbreak sizes outside China. Wellcome Open Res 5, 67 (2020). 21. 21. Q. Bi et al., Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study. Lancet Infect Dis 20, 911–919 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 22. 22. Y. Zhang, Y. Li, L. Wang, M. Li, X. Zhou, Evaluating Transmission Heterogeneity and Super-Spreading Event of COVID-19 in a Metropolis of China. Int J Environ Res Public Health 17, (2020). 23. 23. D. Miller et al., Full genome viral sequences inform patterns of SARS-CoV-2 spread into and within Israel. medRxiv, 2020.2005.2021.20104521 (2020). 24. 24. D. Adam et al., Clustering and superspreading potential of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections in Hong Kong. Europe PMC 10.21203/rs.3.rs-29548/v1, (2020). 25. 25. Z. Du et al., Serial Interval of COVID-19 among Publicly Reported Confirmed Cases. Emerg Infect Dis 26, (2020). 26. 26. M. Gandhi, C. Beyrer, E. Goosby, Masks Do More Than Protect Others During COVID-19: Reducing the Inoculum of SARS-CoV-2 to Protect the Wearer. J Gen Intern Med, (2020). 27. 27. T. Watanabe, T. A. Bartrand, M. H. Weir, T. Omura, C. N. Haas, Development of a dose-response model for SARS coronavirus. Risk Anal 30, 1129–1138 (2010). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20497390&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 28. 28. A. Totura et al., Small Particle Aerosol Exposure of African Green Monkeys to MERS-CoV as a Model for Highly Pathogenic Coronavirus Infection. Emerg Infect Dis 26, (2020). 29. 29. M. Das et al., Decreases in community viral load are accompanied by reductions in new HIV infections in San Francisco. PLoS One 5, e11068 (2010). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0011068&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20548786&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 30. 30. M. S. Cohen et al., Antiretroviral Therapy for the Prevention of HIV-1 Transmission. N Engl J Med 375, 830–839 (2016). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1600693&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27424812&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 31. 31. S. Weller, K. Davis, Condom effectiveness in reducing heterosexual HIV transmission. Cochrane Database Syst Rev, CD003255 (2002). 32. 32. B. T. Mayer et al., Estimating the Risk of Human Herpesvirus 6 and Cytomegalovirus Transmission to Ugandan Infants from Viral Shedding in Saliva by Household Contacts. Viruses 12, (2020). 33. 33. I. Boucoiran et al., Nonprimary Maternal Cytomegalovirus Infection After Viral Shedding in Infants. Pediatr Infect Dis J 37, 627–631 (2018). 34. 34. L. M. Gretebeck, K. Subbarao, Animal models for SARS and MERS coronaviruses. Curr Opin Virol 13, 123–129 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.coviro.2015.06.009&link_type=DOI) 35. 35. M. J. Memoli et al., Validation of the wild-type influenza A human challenge model H1N1pdMIST: an A(H1N1)pdm09 dose-finding investigational new drug study. Clin Infect Dis 60, 693–702 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciu924&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25416753&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) 36. 36. C. M. Chu et al., Viral load distribution in SARS outbreak. Emerg Infect Dis 11, 1882–1886 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1112.040949&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16485474&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F05%2F2020.09.13.20193508.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233768300012&link_type=ISI) [1]: /embed/inline-graphic-1.gif [2]: /embed/graphic-11.gif [3]: /embed/inline-graphic-2.gif [4]: /embed/inline-graphic-3.gif [5]: /embed/inline-graphic-4.gif [6]: /embed/inline-graphic-5.gif [7]: /embed/inline-graphic-6.gif [8]: /embed/inline-graphic-7.gif [9]: /embed/inline-graphic-8.gif [10]: /embed/inline-graphic-9.gif [11]: /embed/inline-graphic-10.gif [12]: /embed/graphic-12.gif [13]: /embed/graphic-13.gif [14]: /embed/inline-graphic-11.gif [15]: /embed/inline-graphic-12.gif [16]: /embed/inline-graphic-13.gif [17]: /embed/inline-graphic-14.gif [18]: /embed/inline-graphic-15.gif [19]: /embed/inline-graphic-16.gif [20]: /embed/inline-graphic-17.gif [21]: /embed/inline-graphic-18.gif [22]: /embed/inline-graphic-19.gif [23]: /embed/inline-graphic-20.gif [24]: /embed/inline-graphic-21.gif [25]: /embed/inline-graphic-22.gif [26]: /embed/inline-graphic-23.gif [27]: /embed/inline-graphic-24.gif [28]: /embed/inline-graphic-25.gif [29]: /embed/inline-graphic-26.gif [30]: /embed/inline-graphic-27.gif [31]: /embed/inline-graphic-28.gif [32]: /embed/inline-graphic-29.gif [33]: /embed/inline-graphic-30.gif [34]: /embed/inline-graphic-31.gif [35]: /embed/inline-graphic-32.gif [36]: /embed/inline-graphic-33.gif [37]: /embed/inline-graphic-34.gif [38]: /embed/inline-graphic-35.gif [39]: /embed/inline-graphic-36.gif