key: cord-0767685-ba6v30jg
authors: Ren, Li-Li; Wang, Ye-Ming; Wu, Zhi-Qiang; Xiang, Zi-Chun; Guo, Li; Xu, Teng; Jiang, Yong-Zhong; Xiong, Yan; Li, Yong-Jun; Li, Xing-Wang; Li, Hui; Fan, Guo-Hui; Gu, Xiao-Ying; Xiao, Yan; Gao, Hong; Xu, Jiu-Yang; Yang, Fan; Wang, Xin-Ming; Wu, Chao; Chen, Lan; Liu, Yi-Wei; Liu, Bo; Yang, Jian; Wang, Xiao-Rui; Dong, Jie; Li, Li; Huang, Chao-Lin; Zhao, Jian-Ping; Hu, Yi; Cheng, Zhen-Shun; Liu, Lin-Lin; Qian, Zhao-Hui; Qin, Chuan; Jin, Qi; Cao, Bin; Wang, Jian-Wei
title: Identification of a novel coronavirus causing severe pneumonia in human: a descriptive study
date: 2020-02-11
journal: Chin Med J (Engl)
DOI: 10.1097/cm9.0000000000000722
sha: a4f38656b672b4a4821a5d12ebfe78dc51fc2ba0
doc_id: 767685
cord_uid: ba6v30jg

BACKGROUND: Human infections with zoonotic coronaviruses (CoVs), including severe acute respiratory syndrome (SARS)-CoV and Middle East respiratory syndrome (MERS)-CoV, have raised great public health concern globally. Here, we report a novel bat-origin CoV causing severe and fatal pneumonia in humans. METHODS: We collected clinical data and bronchoalveolar lavage (BAL) specimens from five patients with severe pneumonia from Jin Yin-tan Hospital of Wuhan, Hubei province, China. Nucleic acids of the BAL were extracted and subjected to next-generation sequencing. Virus isolation was carried out, and maximum-likelihood phylogenetic trees were constructed. RESULTS: Five patients hospitalized from December 18 to December 29, 2019 presented with fever, cough, and dyspnea accompanied by complications of acute respiratory distress syndrome. Chest radiography revealed diffuse opacities and consolidation. One of these patients died. Sequence results revealed the presence of a previously unknown β-CoV strain in all five patients, with 99.8% to 99.9% nucleotide identities among the isolates. These isolates showed 79.0% nucleotide identity with the sequence of SARS-CoV (GenBank NC_004718) and 51.8% identity with the sequence of MERS-CoV (GenBank NC_019843). The virus is phylogenetically closest to a bat SARS-like CoV (SL-ZC45, GenBank MG772933) with 87.6% to 87.7% nucleotide identity, but is in a separate clade. Moreover, these viruses have a single intact open reading frame gene 8, as a further indicator of bat-origin CoVs. However, the amino acid sequence of the tentative receptor-binding domain resembles that of SARS-CoV, indicating that these viruses might use the same receptor. CONCLUSION: A novel bat-borne CoV was identified that is associated with severe and fatal respiratory disease in humans.

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 (a), Betacoronaviruses (b), Gammacoronaviruses (g), and Deltacoronaviruses (d). [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 a-CoVs (229E and NL63) and four b-CoVs (OC43, HKU1, severe acute respiratory syndrome [SARS]-CoV, and Middle East respiratory syndrome [MERS]-CoV). [1] [2] [3] [4] All human CoVs are zoonotic as a distinguishing characteristic. [5] 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 794 deaths worldwide. [8] Since its discovery in Middle Eastern countries in 2012, MERS-CoV has infected 2260 people with a current case fatality rate of 35.5%. [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 China and Ethics Commission of the Jin Yin-tan Hospital of Wuhan (No. KY-2020-01.01). The requirement for written informed consent was waived given the context of emerging infectious diseases.

Bronchoalveolar lavage fluid (BAL) samples were collected from five patients hospitalized with pneumonia in Jin Yintan Hospital of Wuhan, 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 mL 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 instructions's in a biosafety III laboratory. A 50-mL 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), [11] adapter trimming, lowquality 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), [12] host removal by bmtagger (using human genome GRCh38 and yh-specific sequences as reference), [13] and ribosomal reads removal by Sort-MeRNA (version: 2.1b). [14] 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). [15] 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). 

Spot slides were prepared by applying 20 mL of the virusinfected 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 0 ,6 0 -diamidino-2phenylindole 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).

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 groundglass 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 .

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 b-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), [16] ranging from 71,883 (0.27% among all reads) in Patient 4 to 37,247,818 (85%) 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 0 -ORF1ab-S-E-M-N-3 0 , 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 0 -terminus to the 3 0 -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 . Surprisingly, RNA-dependent RNA polymerase (RdRp), which is the most highly conserved sequence among different CoVs, [1, 4] only 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 batorigin 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. [19] 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 ]. [20] 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 were 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 ].

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 [22] ; however, the nt identities among the viruses are only 85.7% to 86.8%. Phylogenetic analysis showed that this virus forms 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. [26] 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.

Coronavirus diversity, phylogeny and interspecies jumping

Coronaviruses post-SARS: update on replication and pathogenesis

Epidemiology and cause of severe acute respiratory syndrome (SARS) in Guangdong, People's Republic of China

Isolation of a novel coronavirus from a man with pneumonia in Saudi Arabia

Hosts and sources of endemic human coronaviruses

Bats as 'special' reservoirs for emerging zoonotic pathogens

SARS and MERS: recent insights into emerging coronaviruses

Summary of Probably SARS Cases With Onset of Illness From

Middle East respiratory syndrome coronavirus (MERS-CoV) outbreak in South Korea, 2015: epidemiology, characteristics and public health implications

Middle East Respiratory Syndrome Coronavirus (MERS-CoV) -Saudi Arabia. Geneva: World Health Organization

Trimmomatic: a flexible trimmer for Illumina sequence data

The diploid genome sequence of an Asian individual

SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data

Improved metagenomic analysis with Kraken 2

Fast and accurate short read alignment with Burrows-Wheeler transform

ORF8-related genetic evidence for Chinese horseshoe bats as the source of human severe acute respiratory syndrome coronavirus

Discovery of a rich gene pool of bat SARS-related coronaviruses provides new insights into the origin of SARS coronavirus

Isolation and characterization of a bat SARS-like coronavirus that uses the ACE2 receptor

Receptor recognition and cross-species infections of SARS coronavirus

Genomic characterization and infectivity of a novel SARS-like coronavirus in Chinese bats

Global epidemiology of bat coronaviruses

Bats and their virome: an important source of emerging viruses capable of infecting humans

Ecology, evolution and classification of bat coronaviruses in the aftermath of SARS

Isolation and characterization of viruses related to the SARS coronavirus from animals in Southern China

Origin and evolution of pathogenic coronaviruses

Identification of a novel coronavirus causing severe pneumonia in human: a descriptive study

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. 

None.