Coronaviruses (CoVs) are enveloped viruses with a single positive-stranded RNA genome (∼26–32 kb in length). They belong to the subfamily Orthocoronavirinae under the family Coronaviridae, and are classified into four genera: Alphacoronaviruses (α), Betacoronaviruses (β), Gammacoronaviruses (γ), and Deltacoronaviruses (δ).[1,2] The viral genome normally encodes four structural proteins, spike (S), envelope (E), membrane (M), and nucleocapsid (N), as well as several non-structural proteins and multiple unique accessory proteins.[1,2]
CoVs infect humans and a variety of avian and mammalian species worldwide. There are six CoVs known to infect humans, including two α-CoVs (229E and NL63) and four β-CoVs (OC43, HKU1, severe acute respiratory syndrome [SARS]-CoV, and Middle East respiratory syndrome [MERS]-CoV).[1–4] All human CoVs are zoonotic as a distinguishing characteristic. In particular, bats are regarded as a key reservoir of CoVs, and many human CoVs are believed to have originated from bats.[5,6] Since the beginning of this century, two zoonotic CoVs, SARS-CoV and MERS-CoV, have been identified to cause severe human diseases.[3,4,7] The outbreak of SARS-CoV in 2003 was responsible for 8096 cases and 774 deaths worldwide. Since its discovery in Middle Eastern countries in 2012, MERS-CoV has infected 2494 people with a current case fatality rate of 34.4%.[9,10] These outbreaks have raised public health concerns of the potential for the emergence of another novel zoonotic CoV.
Here, we report a previously unknown bat-origin CoV causing severe and fatal pneumonia in five patients from Wuhan, China. Sequence results revealed that this virus, harboring a single open reading frame gene 8 (ORF8), is phylogenetically closest to bat SARS-like CoV, but is in a separate lineage. Furthermore, the amino acid sequence of the tentative receptor-binding domain (RBD) of this new CoV resembles that of SARS-CoV, indicating that they might use the same receptor. These findings highlight the urgent need for regular surveillance of the interspecies transmission of bat-origin CoV to human populations.
This study was conducted in accordance with the Declaration of Helsinki and was approved by the National Health Commission of the People's Republic of China and Ethics Commission of the Wuhan Jinyintan Hospital (No. KY-2020-01.01). The requirement for written informed consent was waived given the context of emerging infectious diseases.
Clinical specimen and data collection
Bronchoalveolar lavage fluid (BAL) samples were collected from five patients hospitalized with pneumonia in Wuhan Jinyintan Hospital, Wuhan, Hubei province, China from December 18 to 29, 2019. Information was gathered, including clinical data, demographic characteristics, underlying medical conditions, clinical signs and symptoms, chest radiographic findings, clinical laboratory testing results, traveling history, recent animal exposure, and outcomes. The data collected for the cases were deemed by the National Health Commission of the People's Republic of China as the contents of a public health outbreak investigation.
Nucleic acids were extracted from 200 μL BAL of each sample with the Direct-zol RNA Miniprep kit (Zymo Research, Irvine, CA, USA) and Trizol LS (Thermo Fisher Scientific, Carlsbad, CA, USA) according to the manufacturer's instructions in a biosafety III laboratory. A 50-μL elution was obtained from each sample. The DNA/RNA concentrations were measured by a Qubit Fluorometer (Thermo Fisher Scientific). The sequencing library was constructed by a transposase-based methodology and sequenced on an Illumina sequencing platform (Illumina, San Diego, CA, USA). At least 25 million single-end 76-bp reads were generated for each sample on the Illumina NextSeq platform. Quality control processes included removal of low-complexity reads by bbduk (entropy = 0.7, entropy-window = 50, entropy k = 5; version: January 25, 2018), adapter trimming, low-quality reads removal, short reads removal by Trimmomatic (adapter: TruSeq3-SE.fa:2:30:6, LEADING: 3, TRAILING: 3, SLIDING WINDOW: 4:10, MINLEN: 70, version: 0.36), host removal by bmtagger (using human genome GRCh38 and yh-specific sequences as reference), and ribosomal reads removal by SortMeRNA (version: 2.1b). Taxonomic assignment of the clean reads was performed with Kraken 2 against the reference databases, including archaea, bacteria, fungi, human, plasmid, protozoa, univec, and virus sequences (software 2.0.7-beta, database version: August 2, 2019). A negative control sample was processed and sequenced in parallel for each sequencing run as a contamination control. The data were classified by simultaneous alignment to the microbial genome databases comprising viruses, bacteria, fungi, and parasites after filtering of the adapters and human-origin reads. The sequences were confirmed by Sanger sequencing with specific primers and one-step real-time polymerase chain reaction (RT-PCR) Kit (Invitrogen, Carlsbad, CA, USA).
Multiple sequence alignment was performed with the ClustalW program using MEGA software (version 7.0.14). Phylogenetic trees were constructed by means of the maximum-likelihood method with MEGA software (version 7.0.14). The full-genome viral sequences were deposited in the dataset of Global Initiative on Sharing All Influenza Data (No. EPI_ISL_402123, EPI_ISL_403928-31) and the Genome Warehouse in National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences, under Project ID PRJCA002165 that is publicly accessible at https://bigd.big.ac.cn/gwh as of January 2020.
The BAL specimens were inoculated onto Vero cells (American Type Culture Collection [ATCC], CCL-81). All cultures were observed daily for a cytopathic effect (CPE). Maintenance medium containing tosyl-phenylalanine chloromethyl-ketone enzyme at a final concentration of 1 μg/mL was replenished at day 4, and cultures were terminated 7 days after inoculation. The viral particles were negative stained with 1% solution of phosphotungstic acid (pH 7.0) and the morphology was characterized by using 120 kV TECNAI (Thermo Fisher Scientific, Hillsboro, OR, USA) electron microscopy and camera of Gatan832 (Gatan, Pleasanton, CA, USA). The culture supernatants of cells demonstrating CPE were mixed with paraformaldehyde, dried onto formvar/carbon-coated grids, and stained. Viral nucleic acids were confirmed by RT-PCR with specific primers [Supplementary Table 1, http://links.lww.com/CM9/A191].
Spot slides were prepared by applying 20 μL of the virus-infected or non-infected cell suspension onto 12-well Teflon-coated slides. The cells were fixed with 4% paraformaldehyde in 1× phosphate-buffered saline (PBS) for 30 min, washed three times with PBS, blocked, and stained with serum from a convalescent patient or serum from a healthy person for 30 min at 37°C at a dilution of 1:200. Goat anti-human immunoglobulin G conjugated with fluorescein isothiocyanate was used as the secondary antibody (Jackson Immuno Research Laboratories, Inc., West Grove, PA, USA). Nuclei and the cytoplasm were counterstained with 4′,6-diamidino-2-phenylindole and Evans blue (Sigma-Aldrich, St. Louis, MO, USA). Fluorescent images were obtained and analyzed using laser-scanning confocal microscopy (Airyscan LSM880, Zeiss, Berlin, Germany).
General information of patients
Patient 1 was a 65-year-old man who reported a high fever and cough, with little sputum production, at the onset of illness. He had a continuous fever and developed severe shortness of breath 16 days later. He was a vendor at the Huanan Seafood Market, Wuhan, Hubei Province, China. Patient 2, a 49-year-old woman, presented with high fever and dry cough. Five days later, she developed dyspnea and was admitted to the hospital. She was also a worker in the Huanan Seafood Market. Patient 3 was a 52-year-old woman who did not report any market exposure. She was admitted to hospital because of fever, cough, and ground-glass opacity in the chest computed tomography scan. Patient 4 was a 41-year-old man who also presented with high fever and dry cough at the onset of the illness. He developed acute respiratory distress syndrome 7 days later. This patient had no known history of exposure to the Huanan Seafood Market. Patient 5, a 61-year-old man, was admitted to a local hospital with a 7-day history of fever, cough, and dyspnea. He also worked in the market.
With regards to medical history, Patient 4 had hypertension, and Patient 5 had chronic liver disease and abdominal myxoma, whereas none of the other patients had a record of underlying diseases. The demographic and clinical characteristics of the five patients are summarized in Table 1.
Novel CoV identification by next-generation sequencing
The resultant clean reads accounted for 12.0% to 92.0% of the raw reads. Most of the reads could be successfully assigned. Notably, 80.3% of the reads mapped to the viral genome for sample from Patient 5 with the highest proportion of viral reads among the five samples. Nearly all of the viral reads (97%) were classified as Coronaviridae. Similarly, in the other four patients, most of the viral reads were assigned to β-CoVs. Based on de novo assembly and careful curation, a consensus sequence of this CoV was obtained.
A substantial proportion of all sequencing reads mapped to the newly reported CoV genome (BWA mem, version: 0.7.12), ranging from 71,883 (0.3% among all reads) in Patient 4 to 37,247,818 (85.5%) in Patient 5. In addition, very few reads mapped to known bacterial pathogens, including Streptococcus, Acinetobacter baumannii, and Pseudomonas [Figure 1A–E].
The reads mapping to CoVs were assembled, and their genome sequences were confirmed by Sanger sequencing. The nucleotide (nt) similarity among the obtained five whole-genome sequences was 99.8% to 99.9%. The full length of the obtained genome was 29,870 bp with a GC content of 37.99% to 38.02%. The genome organization, 5′-ORF1ab–S–E–M–N–3′, was similar to that of the most well-known bat SARS-like (SL)-CoV [Figure 2A]. In addition, unique accessory ORFs were identified that are characterized in the subgenus Sarbecovirus, encoding putative ORF3, ORF6, ORF7, and ORF8 proteins reading from the 5′-terminus to the 3′-terminus between the structural proteins [Figure 2A].
Homology assessment showed that full-length viral genome sequences have 79.0% nt identity with that of SARS-CoV Tor2 (GenBank NC_004718), 51.8% with that of MERS-CoV (GenBank NC_019843), and 87.6% to 87.7% with those of bat SL-CoV ZC45 and ZXC21 (GenBank MG772933, MG772934), isolated from Chinese horseshoe bats (Rhinolophus sinicus) [Table 2], indicating that the novel CoVs are most similar to bat SL-CoVs.
Compared with bat SL-CoV ZC45, the novel CoVs showed 75.9%, 98.6%, 93.2% to 93.4%, and 91.1% nt identities in the S, E, M, and N genes, respectively. Overall, ORF1ab showed 89.0% nt identity between the novel CoVs and bat SL-CoV ZC45. Surprisingly, RNA-dependent RNA polymerase (RdRp), which is the most highly conserved sequence among different CoVs,[1,4] only showed 86.3% to 86.5% nt identities with bat SL-CoV ZC45. According to the International Committee on Taxonomy of Viruses criteria, a new CoV species could be defined if the nt identity is less than 90% for the conserved RdRp sequence. Thus, we considered that the novel CoVs should be classified as a new species under the subgenus Sarbecovirus of the genus Betacoronavirus.
The phylogenetic trees constructed with the sequences of the RdRp, S, and N genes, and the whole genome using a maximum-likelihood model showed that all five novel CoVs were closely related to bat SL-CoVs ZXC21 and ZC45, but in a separate evolutionary lineage under the subgenus Sarbecovirus [Figure 2B–E], which is consistent with the homology assessment results.
ORF3 and intact ORF8 gene regions were present in the novel CoVs, which are the characteristic features of bat-origin CoVs.[17,18] ORF3 of the novel CoVs showed 87.8% nt and 90.9% amino acid (aa) identities with bat SL-CoV ZC45, but less than 76.8% nt and 76.0% aa identities with the other members in the subgenus Sarbecovirus. In addition, ORF8 of the novel CoVs showed 88.5% nt and 94.2% aa identity with bat SL-CoV ZC45, respectively, and less than 67.8% nt and 58.6% aa identity, respectively, with other members of Sarbecovirus. These findings further indicated that the novel CoVs are of bat origin.
The RBD in the CoV S protein determines the host range. The RBD aa sequences of the novel CoV showed several distinct features, including higher aa identities with those of SARS-CoV (73.8–74.8%) and human angiotensin-converting enzyme 2 (hACE2)-using SL-CoVs (76.4–76.9%) than those of SL-CoVs incapable of using hACE2 (61.5–64.1%). The novel CoV does not possess the deletions commonly found in the RBD of SL-CoVs incapable of using hACE2 as a receptor [Figure 2F]. In addition, the five critical aa residues interacting with hACE2 in SARS-CoV RBD (Y442, L472, N479, D480, T487) differ from the corresponding residues in the novel CoVs (L, F, Q, S, N), although these residues possess similar polarity.[20,21] These results suggested that the novel CoVs might still use hACE2 as the receptor.
CPE was observed in 30% of Vero cells inoculated with the new CoV after two passages [Figure 3A]. The cells showed a round, refractive, and syncytium appearance. The Vero cells with CPE were further examined using negative-staining electron microscopy, demonstrating characteristic CoV particles with surface projections [Figure 3B]. Immunofluorescent assays of the culture of Vero cells showing CPE with the convalescent serum from patients showed green signals in the cytoplasm, with no signals detected in wells containing control serum, indicating the presence of viral particles in the cells [Figure 3C].
Clinical features and outcomes of the patients
The clinical features and laboratory test results of the five patients are summarized in Table 2. Fever, cough, and dyspnea were the most common symptoms. The white blood cell counts varied among these patients, but the lymphocyte counts were generally low. The alanine aminotransferase and serum creatine levels were normal or only slightly increased. Bilateral ground-glass opacities and consolidation were observed on chest radiography from two representative patients, Patient 2 based on aortic arch scan [Figure 4A] and pulmonary vein scan [Figure 4B] on day 10 after symptoms onset and Patient 5 taken on day 12 [Figure 4C] and 13 [Figure 4D] after symptoms onset.
Several complications were observed in these patients. Four of the five patients (except for Patient 3) developed acute respiratory distress syndrome requiring oxygen therapy, and two patients were given extracorporeal membrane oxygenation. Two patients (Patients 1 and 5) experienced secondary infections, and Patient 5 later developed septic shock as well as acute kidney injury, and ultimately died of multi-organ failure. Patient 3 was discharged on January 8, 2020 (day 17 after symptoms onset). The other three patients were still hospitalized at the time of manuscript preparation. The treatments for these patients were shown in Table 1.
In this study, we identified a previously unknown CoV from patients suffering from severe pneumonia. The whole-genome sequences of the viruses were obtained by a next-generation sequencing approach from all five patients, demonstrating overwhelmingly dominant viral reads in the BAL samples. Among the five novel CoV genome sequences, the nt identities reached up to 99.8% to 99.9%. The viruses successfully isolated from the patients could also be effectively recognized by serum from a convalescent patient. These findings primarily indicate that the novel CoV is associated with the pneumonia that developed in these patients. However, it remains to be determined whether this novel CoV is capable of causing similar diseases in experimental animals.
Sequence homology analysis of the viral genome showed that the CoV identified in this study is distinct from any of the known human CoVs, including SARS-CoV and MERS-CoV. The most closely related known viruses are two bat SL-CoVs (GenBank accession nos. MG772933, 772934) identified in 2005 in Zhoushan, Zhejiang, China, which is geographically distant from Wuhan; however, the nt identities among the viruses are only 87.6% to 87.7%. Phylogenetic analysis showed that this virus forms a single clade. Collectively, these data indicate that this CoV should be considered a new species. The outbreak of SARS in 2003 largely boosted awareness of threats caused by emerging CoVs. Consequently, great efforts have been made to monitor novel emerging CoVs and to trace their origins so as to establish a risk assessment and alert system for preventing potential epidemics in the human population. Clarification of the coronavirome in animals, particularly in bats as a key reservoir of a wide range of CoVs, should be a priority for any task force.[23,24]
A few striking features of these novel CoVs indicated that they are of bat-origin. First, the genome sequences of the novel CoVs show high similarity with that of bat SL-CoV ZC45. Second, the phylogenetic analysis indicated that these viruses are evolutionarily close to bat SL-CoVs ZXC21 and ZC45. Third, all of these novel CoVs contain ORF3 and intact ORF8 gene regions, which are characteristic features of bat-origin CoVs.[17,20] Moreover, the aa sequences of the N-terminal domains (NTDs) of the novel CoVs were very similar to those of ZC45 and ZXC21, whereas the RBD of the novel CoV showed higher aa sequence identity to that of SARS-CoV than to those of ZC45 and ZXC21, suggesting that a recombination event might have occurred at the region between the NTD and RBD of the S gene, facilitating the interspecies transmission.
Owing to the lack of epidemic information at present, the transmission modes of the novel CoV remain obscure. It is notable that three of the five patients had a history of recent exposure to a seafood market in Wuhan. However, the origin of infection is unknown at the time of manuscript preparation. It is assumed that the zoonotic CoV jumped to humans through an intermediate host; for example, camel is suspected as the intermediate host of MERS-CoV, whereas the palm civet may have contributed to the interspecies transmission of SARS-CoV to humans.[25,26] Bat CoVs may evolve to adapt to using humans as a host during their circulation in a mammalian host, thereby enabling them to effectively infect humans. However, two of our patients did not have a history of exposure to the seafood market. Therefore, further investigation will be needed to determine the potential of multiple infection sources responsible for this uncommon outbreak.
One of the most striking and concerning features of this virus is its ability to cause severe respiratory syndrome. The disease progressed rapidly with a major presentation of lower respiratory pathology. Notably, no obvious upper respiratory tract symptoms such as a sore throat and rhinorrhea were present in these patients. Therefore, further exploration is needed on the distribution of the viral receptor in the organs to potentially account for pathogenesis development. In addition, the possibility of unrecognized mild infections or subclinical infections should be clarified, as identification of such infections is critical to control spread of the disease. Development of serological assays would be largely beneficial to detect such types of infection at the population level.
In conclusion, we identified a novel bat-borne CoV associated with a severe and fatal respiratory disease in humans. The emergence of this virus poses a potential threat to public health. Therefore, clarification of the source and transmission mode of these infections is urgently needed to prevent a potential epidemic.
We would like to thank Dr. Chen Wang, Dr. Zhen-Dong Zhao, Dr. Fei Guo (Chinese Academy of Medical Sciences & Peking Union Medical College), and Dr. Ming-Kun Li (Beijing Institute of Genomics, Chinese Academy of Sciences) for critical reading of the manuscript and helpful discussions. We thank the colleagues who contributed to sample collection and experiments.
This study was supported by grants from the Chinese Academy of Medical Sciences (CAMS) Innovation Fund for Medical Sciences (No. 2016-I2M-1-014), the National Major Science & Technology Project for Control and Prevention of Major Infectious Diseases in China (Nos. 2017ZX10103004, 2018ZX10305409, 2017ZX10204401), and the National Natural Science Foundation (No. 81930063).
Conflicts of interest
1. Woo PC, Lau SK, Huang Y, Yuen KY. Coronavirus diversity, phylogeny and interspecies jumping. Exp Biol Med (Maywood)
2009; 234:1117–1127. doi: 10.3181/0903-MR-94.
2. Perlman S, Netland J. Coronaviruses post-SARS: update on replication and pathogenesis. Nat Rev Microbiol
2009; 7:439–450. doi: 10.1038/nrmicro2147.
3. Zhong NS, Zheng BJ, Li YM, Poon LL, Xie ZH, Chan KH, et al. Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People's Republic of China, in February, 2003. Lancet
2003; 362:1353–1358. doi: 10.1016/s0140-6736(03)14630-2.
4. Zaki AM, van Boheemen S, Bestebroer TM, Osterhaus AD, Fouchier RA. Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia. N Engl J Med
2012; 367:1814–1820. doi: 10.1056/NEJMoa1211721.
5. Corman VM, Muth D, Niemeyer D, Drosten C. Hosts and sources of endemic human coronaviruses. Adv Virus Res
2018; 100:163–188. doi: 10.1016/bs.aivir.2018.01.001.
6. Brook CE, Dobson AP. Bats as ‘special’ reservoirs for emerging zoonotic pathogens. Trends Microbiol
2015; 23:172–180. doi: 10.1016/j.tim.2014.12.004.
7. de Wit E, van Doremalen N, Falzarano D, Munster VJ. SARS and MERS: recent insights into emerging coronaviruses. Nat Rev Microbiol
2016; 14:523–534. doi: 10.1038/nrmicro.2016.81.
8. Summary of Probably SARS Cases With Onset of Illness From 1 November 2002 to 31 July 2003. Geneva: World Health Organization. Available from: http://www.who.int/csr/sars/country/table2004_04_21/en/
. [Accessed January 20, 2020].
9. Kim KH, Tandi TE, Choi JW, Moon JM, Kim MS. Middle East respiratory syndrome coronavirus (MERS-CoV) outbreak in South Korea, 2015: epidemiology, characteristics and public health implications. J Hosp Infect
2017; 95:207–213. doi: 10.1016/j.jhin.2016.10.008.
10. Middle East Respiratory Syndrome Coronavirus (MERS-CoV). Geneva: World Health Organization. Available from: https://www.who.int/emergencies/mers-cov/en/
. [Accessed January 20, 2020].
11. BBDuk Guide. Berkeley: Joint Genome Institute. Available from: https://jgi.doe.gov/data-and-tools/bbtools/bb-tools-user-guide/bbduk-guide/
. [Accessed August 1, 2019]
12. Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics
2014; 30:2114–2120. doi: 10.1093/bioinformatics/btu170.
13. Wang J, Wang W, Li R, Li Y, Tian G, Goodman L, et al. The diploid genome sequence of an Asian individual. Nature
2008; 456:60–65. doi: 10.1038/nature07484.
14. Kopylova E, Noé L, Touzet H. SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics
2012; 28:3211–3217. doi: 10.1093/bioinformatics/bts611.
15. Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol
2019; 20:257doi: 10.1186/s13059-019-1891-0.
16. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics
2009; 25:1754–1760. doi: 10.1093/bioinformatics/btp324.
17. Wu Z, Yang L, Ren X, Zhang J, Yang F, Zhang S, et al. ORF8-related genetic evidence for Chinese horseshoe bats as the source of human severe acute respiratory syndrome coronavirus. J Infect Dis
2016; 213:579–583. doi: 10.1093/infdis/jiv476.
18. Hu B, Zeng LP, Yang XL, Ge XY, Zhang W, Li B, et al. Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus. Plos Pathog
2017; 13:e1006698doi: 10.1371/journal.ppat.1006698.
19. Ge XY, Li JL, Yang XL, Chmura AA, Zhu G, Epstein JH, et al. Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor. Nature
2013; 503:535–538. doi: 10.1038/nature12711.
20. Li F. Receptor recognition and cross-species infections of SARS coronavirus. Antiviral Res
2013; 100:246–254. doi: 10.1016/j.antiviral.2013.08.014.
21. Hu D, Zhu C, Ai L, He T, Wang Y, Ye F, et al. Genomic characterization and infectivity of a novel SARS-like coronavirus in Chinese bats. Emerg Microbes Infect
2018; 7:154doi: 10.1038/s41426-018-0155-5.
22. Wong ACP, Li X, Lau SKP, Woo PCY. Global epidemiology of bat coronaviruses. Viruses
2019; 11: pii: E174. doi: 10.3390/v11020174.
23. Smith I, Wang LF. Bats and their virome: an important source of emerging viruses capable of infecting humans. Curr Opin Virol
2013; 3:84–91. doi: 10.1016/j.coviro.2012.11.006.
24. Drexler JF, Corman VM, Drosten C. Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS. Antiviral Res
2014; 101:45–56. doi: 10.1016/j.antiviral.2013.10.013.
25. Guan Y. Isolation and characterization of viruses related to the SARS coronavirus from animals in Southern China. Science
2003; 302:276–278. doi: 10.1126/science.1087139.
26. Cui J, Li F, Shi ZL. Origin and evolution of pathogenic coronaviruses. Nat Rev Microbiol
2019; 17:181–192. doi: 10.1038/s41579-018-0118-9.