key: cord-0727388-6izmhty0 authors: Watanabe, Toru; Bartrand, Timothy A.; Weir, Mark H.; Omura, Tatsuo; Haas, Charles N. title: Development of a Dose‐Response Model for SARS Coronavirus date: 2010-05-20 journal: Risk Anal DOI: 10.1111/j.1539-6924.2010.01427.x sha: 3b4901b9a9ead003c7f0a971d23095a2bb19649a doc_id: 727388 cord_uid: 6izmhty0 In order to develop a dose‐response model for SARS coronavirus (SARS‐CoV), the pooled data sets for infection of transgenic mice susceptible to SARS‐CoV and infection of mice with murine hepatitis virus strain 1, which may be a clinically relevant model of SARS, were fit to beta‐Poisson and exponential models with the maximum likelihood method. The exponential model (k= 4.1 × l0(2)) could describe the dose‐response relationship of the pooled data sets. The beta‐Poisson model did not provide a statistically significant improvement in fit. With the exponential model, the infectivity of SARS‐CoV was calculated and compared with those of other coronaviruses. The does of SARS‐CoV corresponding to 10% and 50% responses (illness) were estimated at 43 and 280 PFU, respectively. Its estimated infectivity was comparable to that of HCoV‐229E, known as an agent of human common cold, and also similar to those of some animal coronaviruses belonging to the same genetic group. Moreover, the exponential model was applied to the analysis of the epidemiological data of SARS outbreak that occurred at an apartment complex in Hong Kong in 2003. The estimated dose of SARS‐CoV for apartment residents during the outbreak, which was back‐calculated from the reported number of cases, ranged from 16 to 160 PFU/person, depending on the floor. The exponential model developed here is the sole dose‐response model for SARS‐CoV at the present and would enable us to understand the possibility for reemergence of SARS. The causal pathogen of severe acute respiratory syndrome (SARS) is a newly isolated coronavirus (SARS-CoV) that first appeared in late 2002 in Guangdong Province, People's Republic of China. In the spring of 2003, a large outbreak of this severe pneumonia occurred in Hong Kong and rapidly spread throughout the world. Ultimately, 8,096 cases of SARS were identified in 29 countries or areas and 774 patients reportedly died. (1) The rapid transmission and high mortality rate made SARS a global threat for which no efficacious therapy was available and empirical strategies had to be used to treat the patients. Since the SARS pandemic in 2003, developments of vaccines (2) and antivirals (3) have been rapidly proceeding for treatment of SARS patients and prevention of its reemergence. In parallel, mathematical models expressing the SARS propagation from person to person have been developed to discuss countermeasures against such a disease transmission from community to national scales. (4, 5) These models successfully explained the epidemic curve by estimating model parameters of infectious rate, removal rate, and so on. However, these models did not incorporate primary exposures or transmission via environmental reservoirs and the model parameters would not be consistent and highly influenced by numerous factors in environment and human behavior. Such models with parameters estimated on a case-bycase basis are of little use in any outbreaks occurring in other situations. On the other hand, dose-response models to characterize the interaction between human and virus seems robust and applicable to assessing the risk of SARS via any possible routes of infection (e.g., aerosols containing virus particles, surface or hand contamination). Nevertheless, no dose-response model for SARS-CoV is available at the moment mainly due to unavailability of data sets challenging humans or animals with this virus. In this article, we develop a dose-response model for SARS-CoV based on two data sets for infection of transgenic mice susceptible to SARS-CoV and infection of mice with murine hepatitis virus strain 1 (MHV-1) that may be a clinically relevant model of SARS. (6) And also, as an example of model application, we analyze the epidemiological data of a SARS outbreak that occurred at an apartment complex in Hong Kong in 2003 with the developed model. Coronaviruses cause acute and chronic respiratory, enteric, and central nervous system (CNS) diseases in many species of animals including humans. Previous to the emergence of SARS, there were two prototype human coronavirues OC43 and 229E, both etiologic agents of the common cold. SARS was the first example of serious illness in humans caused by coronaviruses. (7) Since the outbreak in 2003, many researchers have worked to elucidate structures of viral genes and proteins, mechanisms of infection and replication, and its pathogenesis. Coronaviruses are divided into three groups (groups I to III) based on the genome sequences and SARS-CoV may be a member of group II as well as murine hepatitis virus (MHV), bovine coronavirus, porcine hemagglutinating en-cephalomyelitis virus (HEV), equine coronavirus, and human coronavirues OC43 and NL63, which also cause respiratory infections. Among them, MHV that infects both mice and rats often has been studied as a suitable model of human coronavirus diseases. Coronaviruses are generally restricted in their host range and viruses associated with disease in one species can be limited in their ability to replicate in other species. (8) SARS-CoV differs from this general pattern. This virus infects and replicates in mice, ferrets, hamsters, cats, and several species of nonhuman primates (cynomolgus and rhesus macaques, African green monkeys, and marmosets). Nevertheless, most attempts to reproduce completely human clinical disease and pathological findings in these animals failed. On the other hand, De Albuquerque et al. (9) demonstrated that intranasal infection of A/J mice with MHV-1 produced pulmonary pathological features of SARS. From the fact that all MHV-1-infected A/J mice developed progressive interstitial pneumonitis, including dense macrophage infiltrates, giant cells, and hyaline membranes, resulting in death of all animals, they concluded that A/J mice infected with MHV-1 would be a potentially useful small animal model of human SARS that defines its pathogenesis and suggests treatment strategies. We employed the data set (2 in Table I ) challenging A/J mice with MHV-1 as a surrogate data set for SARS-CoV. The data set was obtained by monitoring the survival of four groups of the mice (n = 5 per group) for 21 days after intranasal inoculation with MHV-1 at 5, 50, 500, and 5,000 PFU, respectively. In other approaches, transgenic (tg) mice expressing the human receptor for SARS-CoV, which are very susceptible to SARS-CoV, have been developed and used for pathogenesis studies. (10) DeDiego et al. (11) reported a data set for infection of the tg mice with recombinant SARS-CoV (1 in Table I ). In their experiment, four groups of the tg mice (n = 2 to 6 per group) were intranasally inoculated with 240, 800, 2,400, and 12,000 PFU of rSARS-CoV, respectively, and the survival was monitored 13 days. The details of data sets 1 and 2 are shown in the Appendix. The endpoint was death of mice for both data sets 1 and 2. Since neither infection nor illness of mice has been analyzed in these experiments, we assumed that all mice with the illness died. In general, mortality of animals with severe illness like SARS would be dependent on their physical strength, which indicates a potential of recovery from the illness, and availability of medical treatment. As described above, a high mortality was observed during the SARS outbreak, although the patients were treated with some empirical strategies. In the experiments using mice, higher mortality can be expected since the mice developing illness had never received any medical treatments and basically the physical strength of mice would be much lower than that of humans. This could support the above assumption and therefore we decided to use the data sets 1 and 2 to develop the dose-response model with the endpoint of human illness, that is, SARS. In order to compare the infectivity of SARS-CoV with other coronaviruses, we also collected data sets (3 to 15 in Table I ) challenging humans and animals with HCoV-229E, MHV, HEV, and infectious bronchitis virus (IBV), which had been studied for many years before the emergence of SARS. Except for one data set (15 via tracheal inoculation), all these data sets were obtained by the intranasal inoculation, which is probably the primary route of infection with SARS-CoV. Each data set in Table I was fitted to two types of dose-response model, that is, exponential and beta-Poisson models, (18) with maximum likelihood method and goodness-of-fit to the models compared based on their likelihoods. The general forms of these models are as follows: where p(d) is the risk of illness at the dose of d; and k, N 50 , and α are parameters specific for the pathogen. Parameter k in the exponential model equals the reciprocal of the probability that a single pathogen will initiate the response. Parameters N 50 and α in the beta-Poisson model are the median dose to get the response and the exponential fitting parameter, respectively. The beta-Poisson model was fitted to only four data sets 2, 3, 14, and 15 since it is not meaningful to fit this model to data sets including less than two dose points corresponding to other than 0 or 100% response. The model fitting was done via maximum likelihood estimation, using a quasi-Newton method, implemented in the package of R version 2.6.2 for Windows. (19) The sensitivity analysis of parameters of the best-fitted model was performed with 10,000 bootstrap trials by a program written on R to estimate the doses corresponding to 10% and 50% responses (ID 10 and ID 50 ). During the pandemic of SARS in 2003, Hong Kong was the hardest hit reporting area with 1,755 cases and 299 deaths in a population of 6.7 million. (1) The large community outbreak at an apartment complex named Amoy Gardens affected more than 300 among 20,000 residents in the early stage of the outbreak in Hong Kong. According to the report by Department of Health (DOH), Hong Kong SAR, the outbreak was begun with an index case who visited his family living on the 16th floor in Block E of Amoy Gardens. (20, 21) He stayed overnight and used the toilet there. DOH concluded that the sewage contamination associated with the index case (and other infected persons) excreting coronavirus that gained entry to households through the bathroom floor drain with dried U-traps was the primary cause of the SARS outbreak. The building's ventilation system made a significant contribution to the entry of virus into households. According to McKinney et al., (22) many residents had installed high-powered fans in their small bathrooms with capacities 6 to 10 times higher than the required capacity and the fans created large negative pressure and drew air from waste pipes. The WHO environmental team verified that sewer gas and aerosolized droplets, which the hydraulic action caused by flushing toilets generated, were being drawn into the bathrooms from the waste pipe system. Although DOH mentioned other possible causes of the SARS outbreak, such as personto-person transmission, vectors acting as mechanical carriers for the virus, and environmental contamination, we assumed that all cases were infected via the airborne transmission, reported as the primary cause, in our model application. DOH investigated epidemiological data of 321 confirmed SARS cases. Although the detailed data are not reported in the literature, the location of flats for 99 cases in Block E has been published in the modeling study by Li et al. (21) Since the number of residents at each flat during the outbreak has not been investigated, their group assumed that each flat housed four persons, which was the largest number of cases in any one flat, and calculated the risk of illness, equivalent to the attack rate, as the number of cases divided by the assumed number of residents. (23) By using the dose-response model developed here, we estimated the dose (d) of SARS-CoV for cases in Block E from the attack rate (p(d)). The dose estimated on the above assumption might be lower than the true value because number of residents in some flats would be below four and there were possibly unoccupied flats. We also calculated the attack rate as the number of flats having at least one case divided by the total number of flats and estimated the dose, which would be higher than the real value, from the attack rate. Table II shows the result of fitting models to the data sets. The fits of the beta-Poisson model to data sets 2, 3, 14, and 15 could not be rejected on statistical grounds, as indicated by p fit values greater than 0.05. Similarly, the exponential model provided statistically significant fits to all data sets except 14 and 15. As for data sets 2 and 3, both beta-Poisson and exponential models were acceptable. We tested the null hypothesis that the beta-Poisson model provided a statistically significant improvement in fit to these data sets rather than a more parsimonious (exponential) model by comparing the difference of variance between two models against the critical value (3.84) on the chi-square distribution with 1 degree of freedom. As the result, since the hypothesis was rejected for both data sets, the exponential model was used for the analysis of viral infectivity in the latter part. The estimated parameter k in the exponential model for rSARS-CoV (data set 1) was very close to that for MHV-1 (data set 2). In order to decide whether these data sets could be pooled or not, the following statistic was calculated and compared with the critical value (3.84) in the chi-square distribution with 1 degree of freedom: (18) where Y 1 , Y 2 , and Y p are deviances of fits of exponential model to data sets 1 and 2 and pooled data sets, respectively. As shown in Table II , we yielded ln k = 6.01 and Y p = 1.75 when the exponential model was fitted to the pooled data sets (data set p). The value of (0.17) smaller than 3.84 indicated that data sets 1 and 2 could be pooled to p. Consequently, the exponential model ( Fig. 1) with k = 4.1 × 10 2 was employed as the dose-response model for SARS due to intranasal infection. Table II . ID 10 and ID 50 of SARS-CoV were 43 (95% CI = 20 to 81 PFU) and 280 PFU (95% CI = 130 to 530 PFU), respectively. The doses (ID 10 = 2.0 TCD 50 ; ID 50 = 13 TCD 50 ) for HCoV-229E seem relatively low; however, Schmidt et al. (24) reported that the 50% endpoint assay was about 10 to 30 times less sensitive than the plaque assay with HCoV-229E. Accordingly, the estimated infectivity of SARS-CoV would be comparable to that of HCoV that causes a mild cold in humans. Similar infectivity was also observed in data sets of three-day-old mice infected with MHV-S (data set 4) and eight-week-old rats infected with HEV-67N (data set 14). These animal coronaviruses belong to the same genetic group as SARS-CoV. The infectivity of coronavirus might be related to the viral evolution, although a lot of unknown factors exist. The effect of host age was observed in data sets for MHV-S (data sets 4 to 6) and HEV-67N (data sets 9 to 14) and it was obvious that young mice and rats are more susceptible than old ones. In contrast, infection among children and adolescents was relatively uncommon (25) and the mortality of senior people was comparatively high in the SARS pandemic. (26) This age-dependency of human susceptibility to SARS-CoV is different not only from other coronaviruses but also from other human respiratory viruses. Avian influenza A (H5N1) virus is a highly virulent respiratory virus like SARS-CoV. The epidemiologic investigations show that this virus primarily infects young people and Goicoechea (27) made a comment that the receptor recognized by this virus, which is expressed in the lower respiratory tract, may be expressed in the upper airway in children, increasing the risk of infection. In case of respiratory syncytial virus (SRV), another well-known human respiratory virus, the serious forms of the disease (principally bronchiolitis and interstitial pneumonia) are found most frequently in infants under six months of age and the disease is progressively milder with increasing age. (28) The doses of SARS-CoV for residents in Block E, Amoy Gardens, estimated based on the assumed number of residents and the number of affected flats, are described in Table III . As Yu et al. (23) analyzed, the attack rate was dependent on the floor and the middle floors between levels 14 and 23 had the highest risk since the flat where the index case visited was located there. Therefore, the estimated dose of SARS-CoV (63 to 160 PFU) for residents on the middle floor was also higher than the others. The estimated dose (42 to 117 PFU) for residents on the upper floor was slightly higher than that (16 to 49 PFU) on the lower floor. Li et al. (21) explained the reason as that viral particles were probably transmitted with upflow air movements between flats in this block. He and his colleagues (23) have confirmed through air modeling that the exhaust fans propelled virus-laden air into an outside air shaft, where it was carried upward by natural air currents and into other flats on the upper floor through open windows. Based on our estimate of dose, this transmission mode appears highly plausible, although we have been unable to locate data on viral load in either fecal matter or wastewater. Although we believe that there are no additional published data sets available for model development for SARS-CoV other than those in Table I , partial information supporting our analysis has been provided by some researchers. Roberts et al. (29) developed mouse-adapted SARS-CoV through 15 serial passages in the respiratory tract of young BALB/c mice and they observed the lethality of mice dependent on the dose of the virus. Mice receiving a dose higher than 10 3.9 TCID 50 of the virus died and the 50% lethal dose (LD 50 ) was 10 4.6 TCID 50 . In contrast, ID 50 of 2.8 × 10 2 PFU estimated with the developed model is low. Since ordinary SARS-CoV cannot cause the death of mice, there is a possibility that the virulence of the virus was still low even after its adaptation to mice. Nagata et al. (30) reported that cynomolgus monkeys inoculated intranasally with SARS-CoV of 10 6 TCID 50 did not show any clinical sign and symptoms, while slight histopathological changes and virus antigen-positive cells were detected. Kuiken et al. (31) and McAuliffe et al. (32) also demonstrated the pulmonary replication of SARS-CoV without any severe symptoms in monkeys receiving intranasal inoculation of 10 6 TCID 50 viruses. The fact that as many as 10 6 viruses did not cause lethality in monkeys means that the susceptibility to this virus may be dependent on host species even within primates. Of course, this model originated from data sets for mice not for humans. As Mizgerd and Skerrett (33) pointed out, there are the differences between mice and humans relevant to pneumonia. However, since it is not realistic to obtain the data sets for human infection with SARS-CoV by challenging to volunteers due to its high mortality, this model is still a valuable method to evaluate the risk of SARS. In prior work, it has been shown that data from animal experiments provide reasonable estimates for human susceptibility in inhaled pathogens. (34, 35) It is widely considered that SARS is a respiratory illness. However, SARS patients may also exhibit gastrointestinal symptoms, splenic atrophy, and lymphadenopathy. (7) Among them, diarrhea is a very frequent finding in SARS patients (30% to 40% of patients). (7) Cheng et al. (36) reported that geometric mean titer for SARS-CoV in feces of patients ranged from 5.1 × 10 1 copies/mL over 23 days after the onset of illness (dpi) to 8.9 × 10 4 copies/mL between 12 and 14 dpi. Number of gene copies has not been related yet to the plaque-forming capability of this virus; however, there is no doubt that a huge amount of infectious viruses was excreted from patients for a certain period. According to the DOH report, at the beginning of the outbreak in Block E, a part of infectious viruses excreted from the index case must have reached other households through the waste pipe and ventilation systems. Considering this virus transportation, we can imagine that residents at the households inhaled small doses of SARS-CoV (Table III) as estimated here by model application, although there are several unknown factors such as virus dispersion and inactivation during transportation. The SARS pandemic was brought under control through a concerted global effort, and by July 5, 2003 , no further human-to-human transmission took place. However, there are still several possibilities that might lead to the reemergence of SARS in humans. (37) The potential source of SARS-CoV might come from infected animals circulating in the geographical region, as highlighted by the four community-acquired SARS cases between December 2003 and January 2004 in Guangdong, China. In these recent cases, fortunately, all patients only developed mild symptoms and secondary transmission did not occur since the animal virus had not fully adapted in humans yet. In addition to SARS-CoV, two novel CoVs (HCoV-NL63 and HCoV-HKU1) that cause respiratory illness in humans were identified after the SARS outbreak. These facts indicate the possibility that a pandemic of SARS and other coronavirus infection may be brought by rapid viral evolution. We proposed the exponential model with k = 4. 1 × 10 2 as a dose-response model for SARS coronavirus (SARS-CoV) based on the available data sets. With this model, the doses of SARS-CoV corresponding to 10% and 50% responses (illness) were estimated at 43 PFU (95% CI = 20 to 81 PFU) and 280 PFU (95% CI = 130 to 530 PFU), respectively. The estimated infectivity of SARS-CoV was comparable to those of HCoV-229E, known as an agent of human common cold, and of some animal coronaviruses (MHV-S and HEV-67N) belonging to the same genetic group as SARS-CoV. The developed model was applied to the analysis of the epidemiological data of the SARS outbreak that occurred at an apartment complex in Hong Kong in 2003. From the reported number of cases, it was revealed that the apartment residents would be exposed to a dose of SARS-CoV between 16 and 160 PFU per person, which depends on the floor, during the outbreak. Although the susceptibility to SARS-CoV seems to be host-dependent, the developed model is the sole dose-response model for SARS-CoV at the present and would help us predict the reemergence of SARS in the future. We extend our special thanks to Ms. Stephy Y. M. Chan (the University of Hong Kong) for helping us collect the information on the SARS outbreak in Hong Kong. We acknowledge the support of the Japan Society for the Promotion of Science (JSPS) to Toru Watanabe, under JSPS Postdoctoral Fellowships for Research Abroad, and that of the Center for Advancing Microbial Risk Assessment (CAMRA) to Timothy A. Bartrand Summary of probable SARS cases with onset of illness from 1 Grand challenges in human coronavirus vaccine development Current status of antiviral severe acute respiratory syndrome coronavirus research Modeling the SARS epidemic A double epidemic model for the SARS propagation Murine hepatitis virus strain 1 as a model for severe acute respiratory distress syndrome (SARS) Coronavirus pathogenesis and the emerging pathogen severe acute respiratory syndrome coronavirus Is there an ideal animal model for SARS? Murine hepatitis virus strain 1 produces a clinically relevant model of severe acute respiratory syndrome in A/J mice Lethal infection of K18-hACE2 mice infected with severe acute respiratory syndrome coronavirus Pathogenicity of severe acute respiratory coronavirus deletion mutants in hACE-2 transgenic mice Effects of a "new" human respiratory virus in volunteers Age-dependent response of mice to a mouse hepatitis virus, MHV-S Pathogenicity of mouse hepatitis virus for mice depending upon host age and route of infection Neurotropisim of swine haemagglutinating encephalomyelitis virus (coronavirus) in mice depending upon host age and route of infection Susceptibility of rats of different ages to inoculation with swine haemagglutinating encephalomyelitis virus (a coronavirus) by various routes Cloacal inoculation with the Connecticut strain of avian infectious brouchitis virus: An attempt to produce nephropathogenic virus by in vivo passage using cloacal inoculation Quantitative Microbial Risk Assessment A Language and Environment for Statistical Computing. R Foundation for Statistical Computing Outbreak of severe acute respiratory syndrome (SARS) at Amoy Gardens, Kowloon Bay, Hong Kong, main findings of the investigation Multi-zone modeling of probable SARS virus transmission by airflow between flats in Block E Environmental transmission of SARS at Amoy Gardens Evidence of airborne transmission of the severe acute respiratory syndrome virus Plaque assay and improved yield of human coronaviruses in a human rhabdomyosarcoma cell line Epidemiology of severe acute respiratory syndrome (SARS): Adults and children Initial viral load and the outcomes of SARS Human H5N1 influenza The pathogenesis of respiratory syncytial virus infection in infant ferrets A mouse-adapted SARScoronavirus causes disease and mortality in BALB/c mice Pathology and virus dispersion in cynomolgus monkeys experimentally infected with severe acute respiratory syndrome coronavirus via different inoculation routes Newly discovered coronavirus as the primary cause of severe acute respiratory syndrome Replication of SARS coronavirus administered into the respiratory tract of African Green, rhesus and cynomolgus monkeys Animal models of human pneumonia Dose-response models for inhalation of Bacillus anthracis spores: Interspecies comparisons Quantitative microbial risk assessment model for Legionnaires' disease: Animal model selection and dose-response modeling Viral shedding patterns of coronavirus in patients with probable severe acute respiratory syndrome SARS and other coronaviruses in humans and animals