Identification of a suitable nonhuman primate (NHP) model of COVID-19 remains challenging. Here, we characterized severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in three NHP species: Old World monkeys Macaca mulatta (M. mulatta) and Macaca fascicularis (M. fascicularis) and New World monkey Callithrix jacchus (C. jacchus). Infected M. mulatta and M. fascicularis showed abnormal chest radiographs, an increased body temperature and a decreased body weight. Viral genomes were detected in swab and blood samples from all animals. Viral load was detected in the pulmonary tissues of M. mulatta and M. fascicularis but not C. jacchus. Furthermore, among the three animal species, M. mulatta showed the strongest response to SARS-CoV-2, including increased inflammatory cytokine expression and pathological changes in the pulmonary tissues. Collectively, these data revealed the different susceptibilities of Old World and New World monkeys to SARS-CoV-2 and identified M. mulatta as the most suitable for modeling COVID-19.
Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection was first reported in December 2019. SARS-CoV-2 has become a global pandemic, leading to about 10 million confirmed cases and ~500 thousand deaths. COVID-19 is thought to be the second most widespread infectious disease since the 1918 flu1 and has promoted the global collaboration of clinicians and scientists to investigate the diagnosis, prevention and treatment of COVID-192 and to understand genomic and structural information regarding SARS-CoV-2.3 However, some challenging and critical questions, including the roles of cytokine storm,4 SARS-CoV-2 tropism,5 and virusaEUR”host interaction,6 remain unanswered, and antiviral drugs and vaccines have not been developed,7 which requires an animal model that can recapitulate the important features of COVID-19 in patients.8 In fact, several animal models of COVID-19 were recently reported; these models include hACE2-transgenic mouse models,5,9 a non-transgenic mouse model,10 a golden hamster model,11 a ferret model12 and nonhuman primate (NHP) models.13,14,15 However, in addition to the different genetic backgrounds of the experimental animals, the use of different viral strains, dosages and routes of challenge in these models resulted in different and even controversial conclusions. Therefore, in this study, three species of NHPs were challenged with the same strain of SARS-CoV-2 to identify the suitability of NHP models of COVID-19.
To determine the suitability of NHP species for modeling COVID-19, we experimentally inoculated three species of monkeys belonging to two families (Old World and New World monkeys) with SARS-CoV-2 (Fig. 1a and Supplementary Table 1) according to the methods described for a reported SARS model,16 followed by comprehensive comparisons of the clinical signs of COVID-19, viral shedding and replication, and host responses to viral infection.
First, a rise in body temperature, one of the clinical signs of SARS-CoV-2 infection, was assessed. Increases in body temperature of more than 0.5aEUR?A?C were observed in 11 out of 14 Macaca mulatta (HHH-1aEUR”HHH-14), 4 out of 6 Macaca fascicularis (SXH-1aEUR”SXH-6) and 3 out of 6 Callithrix jacchus (RH-1aEUR”RH-6) (Fig. 1b). Notably, M. fascicularis SXH-6 showed a continuously high body temperature over the 21-day experimental period. The loss of body weight (BW) was noticed at 10 days post inoculation (dpi), when 11 out of 13 M. mulatta showed a decrease in BW by 5.88aEUR”28.57%. Five out of six M. fascicularis showed a decrease in BW by 2.17aEUR”10.51% (Supplementary Table 2), suggesting more severe loss of BW in M. mulatta than in M. fascicularis.
Next, the effects of SARS-CoV-2 infection on the lung were evaluated on the basis of chest radiographs using a four-pattern approach.17 We observed lung abnormalities in the form of nodules, masses, and interstitial patterns in both M. mulatta and M. fascicularis at 12 dpi. In particular, aged M. mulatta showed more interstitial patterning than young or adult M. mulatta (Fig. 1c). Importantly, overall comparison of the chest radiographs showed that the pulmonary abnormalities were more severe in M. mulatta than in M. fascicularis, suggesting that M. mulatta are more susceptible to SARS-CoV-2 infection than M. fascicularis. However, no other respiratory symptoms were observed in any of the three monkey species during this study.
Furthermore, we studied the kinetics of infection by monitoring viral shedding from the respiratory and alimentary tracks and replication in tissues at the indicated time points after infection (Fig. 1a). Nasal swab samples from all three species of monkeys showed peaks in the levels of viral genomic RNA (~aEUR?5.0aEUR?A–aEUR?108 copies/I 1/4 l) at 2 dpi, followed by two smaller peaks from 6 to 8 dpi. Viral RNA was still detectable in some of the swab samples from Old World monkeys taken at 14 dpi. In contrast, we found lower viral RNA levels in the swab samples from C. jacchus. In the peripheral blood, viral RNA levels were very low but detectable in the majority of monkeys at 2 dpi and disappeared from blood samples taken from nearly all monkeys at 10 dpi (Fig. 2a). Moreover, we evaluated the SARS-CoV-2 viral genome in various tissues of the infected monkeys. Viral RNA was detected in the following ten tissues from two M. mulatta (HHH-13 and HHH-14) necropsied at the early stage of infection: the lung, trachea, bronchia, spleen, stomach, rectum, bladder, uterus, hilar lymph node, and mesenteric lymph node. However, at the late stage of infection, viral RNA was found in only the bronchus, laryngopharynx, and spleen of M. mulatta and M. fascicularis (Fig. 2b). Viral RNA and spike protein levels in respiratory tissue samples were confirmed by RNA hybridization with the RNAscope platform (Fig. 2c and Supplementary Table 3) and immunofluorescent assays (Supplementary Fig. 1). We did not find viral RNA in any of the 16aEUR”18 tissues collected from the RH2 and RH5 strains of C. jacchus. These results suggest that SARS-CoV-2 could replicate in multiple tissues of model M. mulatta and M. fascicularis and that in addition to the respiratory route, the fecalaEUR”oral route may be involved in viral transmission.
Finally, host responses to SARS-CoV-2 infection were evaluated by assessments of the production of virus-specific antibodies and secretion of inflammatory cytokines, hematology, and histopathology. We found that SARS-CoV-2 infection induced the production of virus-specific antibodies in the majority of monkeys, with these antibodies becoming detectable at 4 dpi (Fig. 3a, b). Overall, antibody levels in old and adult M. mulatta were higher than those in young M. mulatta. Consistent with the results of virus monitoring, we did not detect virus-specific antibodies in serum samples from infected C. jacchus. Neutralization assays revealed that virus-specific antibodies from most of the animals had neutralizing activity against SARS-CoV-2, with the highest viral titer (1:256) observed in M. fascicularis at 21 dpi (Supplementary Table 4).
Inflammatory responses to infection were determined via assays to detect multiple cytokines in serum samples. Eight cytokines (G-CSF, IL-1A, IL-8, IL-15, IL-18, MCP-1, MIP-1B, and sCD40-L) were detected in most of the infected animals. Notably, M. mulatta HHH-11 and M. fascicularis SXH-6 showed high expression of the inflammatory cytokines IL-10, IL-1A, IL-8, IL-15, and MCP-1 (Fig. 3c, d). Interestingly, the nasal swabs, anal swabs and fecal samples from these two animals exhibited prolonged viral shedding (Fig. 2a), indicating that these inflammatory cytokines were related to SARS-CoV-2 infection. High levels of IL-15, IL-18, and MCP-1 were observed in M. mulatta HHH-3, which showed no antibody induction after SARS-CoV-2 infection (Fig. 3c).
Cellular responses to SARS-CoV-2 infection were analyzed by flow cytometry. Three major populations of immune cells (T cells, B cells, and monocytes) in red blood cell (RBC)-depleted peripheral blood were measured. Regardless of the age of M. mulatta or sex of M. fascicularis, the frequencies of CD4+ T cells, CD8+ T cells, and monocytes, but not B cells, peaked at 2 or 4 dpi and then gradually declined (Supplementary Fig. 2a, b). The lowest frequencies of these three types of lymphocytes were observed at 10 or 12 dpi, suggesting that SARS-CoV-2 infection caused transient lymphopenia in M. mulatta and M. fascicularis.
More specifically, to determine the pathological responses to SARS-CoV-2 infection, we necropsied eight animals at the early and late stages of infection according to viral load and chest radiograph. Severe gross lesions were observed mainly on the lungs, spleen and lymph nodes dissected from M. mulatta and M. fascicularis but not C. jacchus. The main gross lesions included massive pulmonary punctate hemorrhage, swollen hilar and mediastinal lymph nodes, and swollen mesenteric lymph nodes (Supplementary Fig. 3). Histopathological examination showed thickening of the pulmonary septum and the infiltration of inflammatory cells, diffuse hemorrhage and necrosis in the lung. Mild lesions were also observed in the liver, spleen, and kidney (Fig. 4a). The detailed histopathological changes were described in Fig. 4a. In summary, M. mulatta showed pathological changes in all six examined tissues, whereas M. fascicularis had four tissues (lung, trachea, bronchus, and kidney) with lesions and C. jacchus had pathological changes in three examined tissues (lung, liver, and spleen). Furthermore, typical ultrastructural lesions were observed in the lung (increased numbers of type II pneumocytes, RBCs and macrophages in the alveolar airspace), mesenteric lymph node (intact structure with a large number of lymphocytes, secretory granules in the cytoplasm of cells) and spleen (increased number of granules in the macrophages) (Fig. 4b), which was consistent with the results of microscopic inspection. These results revealed that SARS-CoV-2 infection induced damage to multiple other tissues in addition to the pulmonary organs.
SARS-CoV-2, which remains enigmatic, has rapidly spread worldwide, leading to high economic losses and threats to public health. Within 6 months after SARS-CoV-2 was first reported, five kinds of animals were developed into models in which to study COVID-19: the mouse, ferret, tree shrew, golden hamster, and a NHP species. Murine models of COVID-19 were developed by genetic manipulation of human ACE2 through transgene expression,9 CRISPR-cas9 editing,5 and the transient overexpression of ACE2,18 which makes mice susceptible to SARS-CoV-2 infection. However, due to the many differences in their genetics, anatomy and physiology, these murine models of COVID-19 are not suitable for studies of COVID-19 pathogenesis. The ferret and golden hamster are naturally susceptible to SARS-CoV-2 infection.11,12 However, they are also genetically different from humans and are not ideal animal models for studies of vaccines and pathogenesis. The tree shrew is relatively genetically close to primates and can act as an alternative experimental animal to NHPs.19 Unfortunately, our recent study showed that the tree shrew is not susceptible to SARS-CoV-2 infection and may instead be an asymptomatic carrier.20
NHPs are widely recognized as ideal animal models for emerging and re-emerging diseases, including SARS16,21,22 and Middle East respiratory syndrome.23 To date, there have been four reports of NHP (Old World monkeys) models of COVID-1914,15,24,25 that recapitulate different aspects of COVID-19, but these models have also raised some debate, probably due to differences in animal species (M. mulatta and M. fascicularis), virus strain, challenge route, and methods of evaluation. Using comparative analysis, our study revealed that M. mulatta is a more susceptible species to SARS-CoV-2 infection than M. fascicularis and C. jacchus.
Host ACE2 is one of the most important receptors for SARS-CoV-2.26 At the molecular level, amino acid alignment of the ACE2 fragments that interact with the receptor-binding domain of the SARS-CoV-2 spike glycoprotein revealed the ACE2 fragments from M. mulatta, M. fascicularis, and humans share the same sequence, while those from C. jacchus and humans differ by four amino acids (Supplementary Fig. 4). These four amino acids are critically involved in the binding of ACE2 to SARS-CoV-2.27 Y41 and Q42 in hACE2, which bind the S protein via the formation of hydrogen bonds, are substituted by H and E, respectively, in C. jacchus ACE2.28 These findings may partially explain the lower susceptibility of C. jacchus compared to the other species examined to SARS-CoV-2 infection. The sequences of this critical domain in the ACE2 receptors from mouse and ferret differ from that in the hACE2 receptor by eight and seven amino acids, respectively. However, this domain in tree shrew ACE2 differs from that in human ACE2 by ten amino acids. Surprisingly, the golden hamster and human ACE2 receptors differ by only two amino acids. In fact, differences in the number of these amino acids parallel the susceptibility to SARS-CoV-2 infection.
The NHP model in this study recapitulates several clinical features of COVID-19. First, the increased body temperature and abnormal chest radiograph of our model simulate two typical clinical signs in COVID-19 patients, fever and pneumonia, respectively.29,30,31 Importantly, gross lesions of the lungs after necropsy strikingly matched the abnormal spots on chest radiographs taken right before necropsy.
Detection of viral genomic RNA is the most important evidence for the diagnosis of COVID-19 due to SARS-CoV-2 infection. Genomic SARS-CoV-2 RNA is detectable in samples of bronchoalveolar lavage fluid, sputum, nasal swabs, feces, and blood/serum.32 Specifically, nasal swabs had higher levels of viral RNA at the early stage of infection for ~2 weeks.33 In the present study, the dynamics of viral shedding in the NHP models were similar to those reported in COVID-19 patients. In addition, we noticed a high level of viral shedding via feces that persisted for a long time in four model monkeys (Fig. 2a), suggesting that the fecalaEUR”oral route shows great potential risk in the transmission of SARS-CoV-2.34,35 The lung is widely considered to be the target organ of SARS-CoV-2 replication since a large viral copy number was detected in bronchoalveolar lavage fluid from COVID-19 patients.32 At the early stage of infection, we detected high levels of SARS-CoV-2 viral genomic RNA in the lung and nine other tissues (Fig. 2b, c), including the stomach, rectum, bladder, and uterus, which suggests that SARS-CoV-2 may be transmitted through the fecalaEUR”oral, urea, and maternalaEUR”fetal routes.36,37 Importantly, we also detected viral RNA in the spleen and blood of our model, which is consistent with the clinical findings that SARS-CoV-2 is present in patientsaEUR(TM) spleens38 and blood and correlated with the severity of COVID-19.39
In terms of host responses to SARS-CoV-2 infection, our NHP model partially recapitulates COVID-19 through at least three features. Typical histopathology of pneumonia is observed in COVID-19 patients in additional to mild inflammation in liver.40 In addition, secondary lymphoid organs could be directly infected by SARS-CoV-2.38 Consistently, our study revealed that after SARS-CoV-2 infection, M. mulatta and M. fascicularis showed severe histopathological changes in lung, including pneumonia, and inflammation in the spleen and liver (Fig. 4a). The COVID-19-associated cytokine storm is reported to be an important pathogenic factor and a potential promising therapeutic target in COVID-19 patients.4,41 Although not consistent with those in COVID-19 patients, the cytokines induced in our model could be used as markers for the evaluation of anti-inflammatory drugs.42 Elevated serum IL-6 is clinically correlated with the severity of COVID-19 and regarded as an important biomarker.43 In a reported NHP model of COVID-19, the level of IL-6 peaked at 1 dpi and decreased to 0 at 3 dpi.14 In this study, we did not detect IL-6, possibly due to our use of a lower challenge dose of virus, which may also have led to the different cytokines profiles and differences in other host responses in the reported NHP models.14 By flow cytometric analysis of peripheral blood from M. mulatta and M. fascicularis, we noticed changes in the abundance of T cells, B cells, and monocytes after viral infection (Supplementary Fig. 2), indicating the induction of host immune responses and transient lymphopenia by SARS-CoV-2 infection.40
The NHP models of COVID-19 play critical roles in the research and development of SARS-CoV-2 vaccines. So far, among four NHPs modes of COVID-19 reported, three models are in rhesus monkey (M. mulatta) and one in cynomolgus monkey (M. fascicularis). Rhesus monkey models are widely used in evaluation of SARS-CoV-2 vaccine. In fact, four SARS-CoV-2 vaccines are all evaluated in rhesus monkeys, including two inactivated vaccines,44,45 one adenoviral vectored vaccine46 and one DNA vaccine.47 Our results in this study further confirmed that rhesus monkey is a good animal model for evaluation of vaccines against SARS-CoV-2.
In conclusion, we have established a NHP model of COVID-19 through comparative analysis of the clinical symptoms, viral replication and tissue tropism, and host responses to SARS-CoV-2 infection among three species from two families of NHPs in one system. Ultimately, this study identified M. mulatta as suitable for the study of SARS-CoV-2 infection, as it most closely recapitulated human-like conditions. The use of this model is promising for the development of COVID-19 drugs and vaccines and other basic studies.
Materials and methods
All animal procedures were approved by the Institutional Animal Care and Use Committee of the Institute of Medical Biology, Chinese Academy of Medical Science (ethics number: DWSP202002 001), and performed in strict accordance with the guidelines for the National Care and Use of Animals approved by the National Animal Research Authority (PR China) and the Experimental Animal Ethics Committee of the Institute in the ABSL-4 facility of the National Kunming High-level Biosafety Primate Research Center, Yunnan, China.
Three species of monkeys from two primate families (Old World and New World monkeys) used for this study were bred and provided by the Kunming Primate Center of the Chinese Academy of Medical Sciences (Laboratory animal production license# SCXK(Dian)-K2015-0004). The monkeys used in our study were divided into M. mulatta group (14 rhesus monkeys, child 4, adults 6, and old 4), M. fascicularis group (6 adult monkeys, male 3, and female 3), and C. jacchus (6 adult monkeys, male 3, and female 3). All animals were screened for specific pathogens, including SARS-CoV-2, before the experiment (Supplementary Table 1 and Fig. 1a). Throughout this study, animals were fed daily with commercial monkey chow, treat and fruit, and monitored twice a day by the trained personnel. Referencing a NHP model of SARS,16 we inoculated adult and aged Old world monkeys with a total of 4.75aEUR?ml of 106aEUR?pfu/ml SARS-CoV-2 intratracheally (4aEUR?ml), intranasally (0.5aEUR?ml) and on the conjunctiva (0.25aEUR?ml), and a half-dose of the virus was administered to young monkeys. New World monkeys were intranasally challenged with 1.0aEUR?ml of the same virus due to their small body size. After viral inoculation, animals were checked daily for clinical signs. Then, we anesthetized the animals with ketamine (6aEUR?mg/kg) and performed the following experimental procedures: body temperature and weight monitoring, sample collection, chest radiography, and necropsies at the indicated stages (Fig. 1a).
A chest X-ray image of each anesthetized monkey was taken at 55aEUR”75aEUR?V and 8aEUR”12.5aEUR?mA every other day using a MobileCooper mobile digital medical X-ray photography system (Browiner, China). Data were independently evaluated in a double-blind manner by two radiologists via a four-pattern approach (analyses of consolidation, interstitial areas, nodules or masses, and atelectasis).
Blood was collected from the inguinal veins of anesthetized animals for flow cytometric analysis conducted on a CytoFLEX instrument (Beckman, USA) using fluorescent-labeled antibodies against the cellular surface markers CD45, CD3, CD4, CD8, CD14, and CD20.
Real-time RT-PCR was used to quantify the viral RNA genome using TaqMan Fast Virus 1-Step Master Mix (Thermo Fisher Scientific, USA) and purified SARS-CoV-2 RNA as a standard curve on a CFX384 Touch Real-Time PCR Detection System (Bio-Rad, USA). Primers and the following probes specific for the NP gene were synthesized according to sequences reported by the China CDC: Target-2-F:5aEUR2-GGGGAACTTCTCCTGCTAGAAT-3aEUR2, Target-2-R:5aEUR2-CAGACATTTTGCTCTCAAGCTG-3aEUR2, Target-2-P:5aEUR2-FAM-TTGCTGCTGCTTGACAGATT-TAMRA-3aEUR2. One-step RT-PCR was conducted under the following conditions: 25aEUR?A?C for 2aEUR?min, 50aEUR?A?C for 15aEUR?min, 95aEUR?A?C for 2aEUR?min, and 40 cycles at 95aEUR?A?C for 5aEUR?s and 60aEUR?A?C for 31aEUR?s.
Formalin-fixed paraffin-embedded (FFPE) sections (5aEUR?I 1/4 m) were prepared for hematoxylin and eosin (H&E) staining. The slides were evaluated in a double-blind manner by a pathologist. To prepare ultrastructural sections, tissues were embedded in epoxy resin and cut into 60aEUR”70aEUR?nm sections. After staining with 2% uranium acetate, ultrathin sections were examined under an H-7650 Hitachi transmission electron microscope (Japan).
A novel nucleic acid hybridization technique, RNAscope, was performed to detect viral RNA in the FFPE sections using a SARS-CoV-2-specific probe (Advanced Cell Diagnostics, USA) designed for NHP tissues. This probe specifically targets the S gene of the SARS-CoV-2 genome. The experiment was conducted according to the procedures provided by the manufacturer.
SARS-CoV-2-specific antibodies in serum samples were evaluated via a commercial SARS-CoV-2 Antibody Assay Kit (ELISA) (AnyGo, China) based on the protocol provided by the company.
Serum samples were inactivated at 56aEUR?A?C for 30aEUR?min and diluted 1:8, 1:16, 1:32, 1:64, 1:128, and 1:256 with phosphate-buffered solution. A mixture of SARS-CoV-2 and diluted serum was incubated at 37aEUR?A?C for 1aEUR?h and transferred to Vero E6 cells in a 96-well plate. Cells were continuously incubated at 37aEUR?A?C for 72aEUR?h. Cytopathogenic effects were recorded to evaluate the neutralizing activities of serum samples.
To assay 23 inflammatory cytokines, the MILLIPLEX MAP NonHuman Primate Cytokine Magnetic Bead PanelaEUR”Immunology Multiplex Assay (Millipore, USA) was performed according to the manufactureraEUR(TM)s protocol on a Bio-Plex machine (Bio-Rad, USA). The inflammatory cytokines in this panel are IL-1I2, IL-4, IL-5, IL-6, IL-8/CXCL8, G-CSF, GM-CSF, IFN-I3, IL-1ra, IL-2, IL-10, IL-12p40, IL-13, IL-15, IL-17A/CTLA8, MCP-1/CCL2, MIP-1I2/CCL4, MIP-1I+-/CCL3, sCD40-L, TGF-I+-, TNF-I+-, VEGF-A, and IL-18.
The data sets used for the current study are available from the corresponding author upon reasonable request.
Morens, D. M., Daszak, P. & Taubenberger, J. K. Escaping PandoraaEUR(TM)s boxaEUR”another Novel Coronavirus. N. Engl. J. Med. 382, 1293aEUR”1295 (2020).
Wang, D. et al. Clinical characteristics of 138 hospitalized patients with 2019 Novel Coronavirus-infected pneumonia in Wuhan, China. JAMA 323, 1061aEUR”1069 (2020).
Wrapp, D. et al. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 367, 1260aEUR”1263 (2020).
Vaninov, N. In the eye of the COVID-19 cytokine storm. Nat. Rev. Immunol. 20, 277 (2020).
Sun, S.-H. et al. A mouse model of SARS-CoV-2 infection and pathogenesis. Cell Host Microbe. https://doi.org/10.1016/j.chom.2020.05.020 (2020).
To, K. K. et al. Temporal profiles of viral load in posterior oropharyngeal saliva samples and serum antibody responses during infection by SARS-CoV-2: an observational cohort study. Lancet Infect. Dis. 20, P565aEUR”P574 (2020).
Thanh Le, T. et al. The COVID-19 vaccine development landscape. Nat. Rev. Drug Discov. 19, 305aEUR”306 (2020).
Callaway, E. Monkeys and mice enlisted to fight coronavirus. Nature 579, 183 (2020).
Bao, L. et al. The pathogenicity of SARS-CoV-2 in hACE2 transgenic mice. Nature. https://doi.org/10.1038/s41586-020-2312-y (2020).
Hassan, A. O. et al. A SARS-CoV-2 infection model in mice demonstrates protection by neutralizing antibodies. Cell. https://doi.org/10.1016/j.cell.2020.06.011 (2020).
Chan, J. F. et al. Simulation of the clinical and pathological manifestations of Coronavirus disease 2019 (COVID-19) in golden Syrian hamster model: implications for disease pathogenesis and transmissibility. Clin. Infect. Dis. https://doi.org/10.1093/cid/ciaa325 (2020).
Kim, Y. I. et al. Infection and rapid transmission of SARS-CoV-2 in ferrets. Cell Host Microbe 27, 704aEUR”709 (2020).
Yu, P. et al. Age-related rhesus macaque models of COVID-19. Anim. Model Exp. Med. 3, 93aEUR”97 (2020).
Munster, V. J. et al. Respiratory disease and virus shedding in rhesus macaques inoculated with SARS-CoV-2. bioRxiv. https://doi.org/10.1101/2020.03.21.001628 (2020).
Chao, S. et al. Infection with Novel Coronavirus (SARS-CoV- 2) causes pneumonia in the Rhesus macaques. Cell Res. https://doi.org/10.21203/rs.2.25200/v1 (2020).
Kuiken, T. et al. Newly discovered coronavirus as the primary cause of severe acute respiratory syndrome. Lancet 362, 263aEUR”270 (2003).
Weyant, E. Chest X-ray resources. J. Electron. Resour. Med. Libr. 11, 159aEUR”164 (2014).
Sun, J. et al. Generation of a broadly useful model for COVID-19 pathogenesis vaccination, and treatment. Cell 182, 1aEUR”10 (2020).
Fan, Y. et al. Chromosomal level assembly and population sequencing of the Chinese tree shrew genome. Zool. Res. 40, 506aEUR”521 (2019).
Zhao, Y. et al. Susceptibility of tree shrew to SARS-CoV-2 infection. bioRxiv. https://doi.org/10.1101/2020.04.30.029736 (2020).
Rowe, T. et al. Macaque model for severe acute respiratory syndrome. J. Virol. 78, 11401aEUR”11404 (2004).
Qin, C. et al. An animal model of SARS produced by infection of Macaca mulatta with SARS Coronavirus. J. Pathol. 206, 251aEUR”259 (2005).
Van Doremalen, N. & Munster, V. J. Animal models of middle east respiratory syndrome Coronavirus infection. Antivir. Res. 122, 28aEUR”38 (2015).
Rockx, B. et al. Comparative pathogenesis Of COVID-19, MERS and SARS in a non-human primate model. Science 368, 1012aEUR”1015 (2020).
Pan, L. et al. Clinical characteristics of COVID-19 patients with digestive symptoms in Hubei, China: a descriptive, cross-sectional, multicenter study. Am. J. Gastroenterol. 115, 766aEUR”773 (2020).
Hoffmann, M. et al. SARS-CoV-2 cell entry depends on ACE2 and TMPRSS2 and is blocked by a clinically proven protease inhibitor. Cell 181, 271aEUR”280 (2020).
Wang, Q. et al. Structural and functional basis of SARS-CoV-2 entry by using human ACE2. Cell 181, 894aEUR”904 (2020).
Liu, Y. et al. Functional and genetic analysis of viral receptor ACE2 orthologs reveals broad potential host range of SARS-CoV-2. bioRxiv https://doi.org/10.1101/2020.04.22.046565 (2020).
Huang, C. et al. Clinical features of patients infected with 2019 Novel Coronavirus in Wuhan, China. Lancet 395, 497aEUR”506 (2020).
Xu, X. et al. Imaging and clinical features of patients with 2019 novel Coronavirus SARS-CoV-2. Eur. J. Nucl. Med. Mol. Imaging 47, 1275aEUR”1280 (2020).
C, Q. et al. Dysregulation of immune response in patients with COVID-19 in Wuhan, China. Clin. Infect. Dis. https://doi.org/10.1093/cid/ciaa248 (2020).
Wang, W. et al. Detection of SARS-CoV-2 in different types of clinical specimens. JAMA 323, 1843aEUR”1844 (2020).
Zou, L. et al. SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N. Engl. J. Med. 382, 1177aEUR”1179 (2020).
Xiao, F. et al. Evidence for gastrointestinal infection of SARS-CoV-2. Gastroenterology 158, 1831aEUR”1833 (2020).
Xu, Y. et al. Characteristics of pediatric SARS-CoV-2 infection and potential evidence for persistent fecal viral shedding. Nat. Med. 26, 502aEUR”505 (2020).
Dong, L. et al. Possible vertical transmission of SARS-CoV-2 from an infected mother to her newborn. JAMA 323, 1846aEUR”1848 (2020).
Gu, J., Han, B. & Wang, J. COVID-19: gastrointestinal manifestations and potential fecal-oral transmission. Gastroenterology 158, 1518aEUR”1519 (2020).
Chen, Y. et al. The novel severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2) directly decimates human spleens and lymph nodes. medRxiv. https://doi.org/10.1101/2020.03.27.20045427 (2020).
Liao, M. et al. Single-cell landscape of bronchoalveolar immune cells in patients with COVID-19. Nat. Med. 26, 842aEUR”844 (2020).
Xu, Z. et al. Pathological findings of COVID-19 associated with acute respiratory distress syndrome. Lancet Respir. Med. 8, 420aEUR”422 (2020).
DaEUR(TM)Elia, R. V., Harrison, K., Oyston, P. C., Lukaszewski, R. A. & Clark, G. C. Targeting the aEURoeCytokine stormaEUR? for therapeutic benefit. Clin. Vaccin. Immunol. 20, 319aEUR”327 (2013).
Mehta, P. et al. COVID-19: consider cytokine storm syndromes and immunosuppression. Lancet 395, 1033aEUR”1034 (2020).
Zhao, M. Cytokine storm and immunomodulatory therapy in COVID-19: role of chloroquine and anti-IL-6 monoclonal antibodies. Int. J. Antimicrob. Agents 55, 105982 (2020).
Gao, Q. et al. Development of an inactivated vaccine candidate for SARS-CoV-2. Science 369, 77aEUR”81 (2020).
Wang, H. et al. Development of an inactivated vaccine candidate, BBIBP-CorV, with potent protection against SARS-CoV-2. Cell. https://doi.org/10.1016/j.cell.2020.06.008 (2020).
van Doremalen, N. et al. ChAdOx1 nCoV-19 vaccination prevents SARS-CoV-2 pneumonia in rhesus macaques. bioRxiv. https://doi.org/10.1101/2020.05.13.093195 (2020).
Yu, J. et al. DNA vaccine protection against SARS-CoV-2 in rhesus macaques. Science, eabc6284, https://doi.org/10.1126/science.abc6284 (2020).
The authors would like to thank all staff at the National Kunming High-level Biosafety Primate Research Center for providing ABSL-4-related services and Jacky Sun for critical comments on our paper. This study was supported by the National Research and Development Project of China (2020YFC0841100, 2020YFC0846400, and 2020YFA0707600), CAMS Innovation Fund for Medical Sciences (2016-I2M-2-006 and 2020-I2M-CoV19-012), the Talents Project of Yunnan Province (2017HB068), the Special Funds for High-level Health and Family Planning Technical Personnel training of Yunnan Province (D-201653).
The authors declare no competing interests.