Open access
Announcement
20 February 2020

Complete Genome Sequences of Five Human Coronavirus NL63 Strains Causing Respiratory Illness in Hospitalized Children in China

ABSTRACT

We report the complete genome sequences of five human coronavirus NL63 (HCoV-NL63) strains obtained using next-generation sequencing. The five HCoV-NL63 strains were obtained from hospitalized children with severe acute respiratory infection detected in Guangdong, China. This study provides several complete genomes of HCoV-NL63 and improves our understanding of HCoV-NL63 evolution in China.

ANNOUNCEMENT

Human coronavirus NL63 (HCoV-NL63) is a member of the family Coronaviridae, genus Alphacoronavirus, and was first discovered in 2004 (1). HCoV-NL63 is mainly associated with the common cold in children, the elderly, and immunocompromised patients (2, 3). The genome of HCoV-NL63 is about 27 kb with a conserved gene order of 5′-orf1ab-spike (S)-orf3-envelope (E)-membrane (M)-nucleocapsid (N)-poly (A). The species tropism of HCoV-NL63 is determined by spike glycoprotein. HCoV-NL63 and severe acute respiratory syndrome coronavirus (SARS-CoV) share the same cell receptor, angiotensin converting enzyme 2 (ACE-2) (4, 5), for entry into host cells, and HCoV-NL63 is recognized as a common cause of upper respiratory tract infection and has been prevalent worldwide.
Here, nasopharyngeal swab samples were collected from hospitalized children with severe acute respiratory infection in Guangzhou, China, in 2018. This study was performed in strict accordance with human subject protection guidance provided by the Research Ethics Committee of Guangzhou Medical University. The respiratory samples were filtered with 0.22-μm filters, RNA extraction was performed using a Qiagen viral RNA extraction kit, and extracted RNA was used for sequence-independent single-primer amplification (SISPA) (6, 7) as follows: a reverse transcription reaction was performed with SuperScript III reverse transcriptase using a primer containing a fixed sequence, followed by a random hexamer at the 3′ end (FR26RV, GCCGGAGCTCTGCAGATATCNNNNNN). Then, Klenow fragment polymerase (New England Biolabs) was used for DNA synthesis. Finally, PCR amplification was conducted using primers consisting of the fixed portions of the random primers (FR26, GCCGGAGCTCTGCAGATATC). Purified DNA was used for next-generation sequencing (NGS). Libraries were prepared with the Nextera XT kit (Illumina), and paired-end reads (2 × 125 bp) determined using a HiSeq 2500 instrument were used for cleaning and assembling using CLC Genomics Workbench version 11.0. Illumina sequencing yielded about 10 million reads per sample. Reads were assembled into contigs with a de novo assembly model, and the contig sequences were then extracted for subsequent analysis. Partial genome sequences of five HCoV-NL63 strains were obtained by NGS methods. Meanwhile, sets of specific primer pairs were designed and used to amplify the gap region of HCoV-NL63, which was used for genome assemblies using the SeqMan subprogram of the DNAStar software version 7.1.0 with default parameters (Table 1). Finally, five complete genome sequences of HCoV-NL63 were obtained using next-generation sequencing and Sanger sequencing methods together and were designated strains ChinaGD01 (27,531 bp), ChinaGD02 (27,516 bp), ChinaGD03 (27,516 bp), ChinaGD04 (27,532 bp), and ChinaGD05 (27,544 bp). The five HCoV-NL63 strains presented here were aligned using MAFFT version 7.158 (8) and showed 98.5 to ∼99.1% nucleotide homology with the prototype HCoV-NL63 virus (GenBank accession number NC_005831.2) as estimated using MEGA version 5.10 software (9) (Fig. 1).
TABLE 1
TABLE 1 Primers used for the genome sequencing of HCoV-NL63
PrimerSequence (5′–3′)Target (nucleotide position) or amplification methodSize (bp)
1FCCTGGCCTCTTGCTTTTTCACATGT205041,686
1RACTTCGACGGTTGAGAAACAAATAG221891,686
2FCGCGTTAAGAGTGGTTCACCAGGTG220591,623
2RCAAAGCTGCAAGCCGTCCAGTAATT236811,623
3FTTCAATTCAAGCCGATCAACAAGTT236241,600
3RGTCATCAATTAATCGAAGGAACATC252231,600
4FCGAAGAGCCTGTTGTTGGTATAGTC251531,690
4RAACACGCTTCCAACGAGGTTTCTTC268491,690
5FCCAGGGCTGATAAGCCTTCTCAGTT26800754
5RGTGTATCCATATCAAAAACAATATC27553754
6FTGAGGATGTTTGTGTTTGTTTTGAC188641,710
6RGTCAGGAACACCTAATTGTAACATA205731,710
7FTGCGTGGTTGGTTGGGTATGGATGT172341,682
7RACGCTCATACGAACCCTGAATACTA189151,682
8FATTCAGCAACTGGTTCCTTAGATGT154791,808
8RGTTATCGCCACAAACATGAGCACTT172861,808
9FCTCCCTACTATGACACAGCTGAATC140051,596
9RAGCCGCAAAGAGTCTAAGTGTATCT156001,596
10FGACCGTACAACTATTCAAAGTGTTG123981,659
10RGTTCTTTACCACTAATAGCATACTT140561,659
11FGGGCTATGGCTAATGGTTATACAAG98011,435
11RTTTGCGATATTCATGGCACGCTTCA112351,435
12FACCCTTCAGAGTGTTGCTTCATCAT110901,378
12RAGTCGAGCTGCACTAGAACCCCTTG124671,378
13FCAACCACTGTAACTAGCTTTCATGG77582,103
13RCTGCCAAAATAGAATAGCACTCAAC98602,103
14FGTCAAAAGGGTGATGCTGAAGAGGC54242,394
14RTCAACTGACCATTCTCAATGTACTT78172,394
15FTAGAGATGAATTGGGTGTTCGTGTT34242,043
15RGGTCCAACATCACCTGTAACAAATT54662,043
16FGCAGATGTTCCAGATGCTTTTCAAT16371,864
16RGCAACTGTACAAGTGTGGTACTAAT35001,864
17FCAGCAATTATGTTCTTCAGGACTTT5651,118
17RGTGTAAATGTGCGATAAACTGATTG16821,118
18FCTTAAAGAATTTTTCTATCTATAGA11,056
18RCATGCACCAACACTCCAACTCTCAG10561,056
GSP 1CGAAGAGCCTGTTGTTGGTATAGTC3′ RACEaUnknown
APGGCCACGCGTCGACTAGTACTTTTTTTTTTTTTTTTT3′ RACEUnknown
GSP 2CCAGGGCTGATAAGCCTTCTCAGTTNested PCRUnknown
AUAPGGCCACGCGTCGACTAGTACNested PCRUnknown
AAPGGCCACGCGTCGACTAGTACGGGGGGGGGG5′ RACEUnknown
GSP 1GTGTAAATGTGCGATAAACTGATTG5′ RACEUnknown
AUAPGGCCACGCGTCGACTAGTACNested PCRUnknown
GSP 2CATGCACCAACACTCCAACTCTCAGNested PCRUnknown
GSP3CCATGGCCAAAAACAACATCAAAGTNested PCRUnknown
a
RACE, rapid amplification of cDNA ends.
FIG 1
FIG 1 Phylogenetic analysis of HCoV-NL63 based on complete genomes. The complete genomes of 69 coronavirus references were obtained from the GenBank database, and multiple alignments were performed using MAFFT version 7.158 with default parameters. The phylogenetic tree was constructed using the neighbor-joining method with 1,000 bootstraps in MEGA version 5.10 with default parameters. The numbers at the nodes represent bootstrap support. Bootstrap values greater than 70% were considered statistically significant for grouping. The HCoV-NL63 strains presented in this study are marked with red circles.
Only two complete genome sequences of HCoV-NL63 associated with acute respiratory illness have been obtained and reported in China. The complete genome sequence data from our study will provide insight into the evolution and genetic diversity of HCoV-NL63 in China.

Data availability.

The complete genome sequences of the five newly identified HCoV-NL63 strains have been deposited in GenBank under the accession numbers MK334043, MK334044, MK334045, MK334046, and MK334047. The sequencing reads are available in the SRA database under BioProject accession number PRJNA601331.

ACKNOWLEDGMENTS

This research was supported by grants from the National Key Research and Development Program of China (2018YFC1200100), National Natural Science Foundation of China (NSFC 81702047, 81772191, 91842106, and 8181101118), State Key Laboratory of Respiratory Disease (SKLRD-QN-201715 and SKLRD-QN-201912), National Key Technology R&D Program (2018YFC1311900), Guangdong Science and Technology Foundation (2019B030316028), 2015 Thousand Talents Plan Award of China, China Postdoctoral Science Foundation, and the Ph.D. Start-up Fund of the Natural Science Foundation of Guangdong Province, China.

REFERENCES

1.
van der Hoek L, Pyrc K, Jebbink MF, Vermeulen-Oost W, Berkhout RJ, Wolthers KC, Wertheim-van Dillen PM, Kaandorp J, Spaargaren J, Berkhout B. 2004. Identification of a new human coronavirus. Nat Med 10:368–373.
2.
van der Hoek L, Pyrc K, Berkhout B. 2006. Human coronavirus NL63, a new respiratory virus. FEMS Microbiol Rev 30:760–773.
3.
Kaiser L, Regamey N, Roiha H, Deffernez C, Frey U. 2005. Human coronavirus NL63 associated with lower respiratory tract symptoms in early life. Pediatr Infect Dis J 24:1015–1017.
4.
Hofmann H, Pyrc K, van der Hoek L, Geier M, Berkhout B, Pöhlmann S. 2005. Human coronavirus NL63 employs the severe acute respiratory syndrome coronavirus receptor for cellular entry. Proc Natl Acad Sci U S A 102:7988–7993.
5.
Kuba K, Imai Y, Rao S, Gao H, Guo F, Guan B, Huan Y, Yang P, Zhang Y, Deng W, Bao L, Zhang B, Liu G, Wang Z, Chappell M, Liu Y, Zheng D, Leibbrandt A, Wada T, Slutsky AS, Liu D, Qin C, Jiang C, Penninger JM. 2005. A crucial role of angiotensin converting enzyme 2 (ACE2) in SARS coronavirus-induced lung injury. Nat Med 11:875–879.
6.
Reyes GR, Kim JP. 1991. Sequence-independent, single-primer amplification (SISPA) of complex DNA populations. Mol Cell Probes 5:473–481.
7.
Wang Y, Zhu N, Li Y, Lu R, Wang H, Liu G, Zou X, Xie Z, Tan W. 2016. Metagenomic analysis of viral genetic diversity in respiratory samples from children with severe acute respiratory infection in China. Clin Microbiol Infect 22:458.e451–458.e459.
8.
Katoh K, Rozewicki J, Yamada KD. 2019. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform 20:1160–1166.
9.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. 2011. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28:2731–2739.

Information & Contributors

Information

Published In

cover image Microbiology Resource Announcements
Microbiology Resource Announcements
Volume 9Number 820 February 2020
eLocator: 10.1128/mra.01597-19
Editor: Simon Roux, DOE Joint Genome Institute

History

Received: 17 January 2020
Accepted: 4 February 2020
Published online: 20 February 2020

Contributors

Authors

Lu Zhang
Institute of Infectious Disease, Guangzhou Eighth People’s Hospital of Guangzhou Medical University, Guangzhou, China
Mian Gan
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Zhaoyong Zhang
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Xin Li
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Wenkuan Liu
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Airu Zhu
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Jing Sun
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Fang Li
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Yanqun Wang
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Fuchun Zhang
Institute of Infectious Disease, Guangzhou Eighth People’s Hospital of Guangzhou Medical University, Guangzhou, China
Jingxian Zhao
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Rong Zhou
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China
Jincun Zhao
Institute of Infectious Disease, Guangzhou Eighth People’s Hospital of Guangzhou Medical University, Guangzhou, China
State Key Laboratory of Respiratory Disease, National Clinical Research Center for Respiratory Disease, Guangzhou Institute of Respiratory Health, The First Affiliated Hospital of Guangzhou Medical University, Guangzhou, China

Editor

Simon Roux
Editor
DOE Joint Genome Institute

Notes

Address correspondence to Rong Zhou, [email protected], or Jincun Zhao, [email protected].

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures and Media

Figures

Media

Tables

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy