key: cord-1053348-snra23sy
authors: Lee, Wan‐Ji; Chung, Yoon‐Seok; Yoon, Hee Sook; Kang, Chun; Kim, Kisoon
title: Prevalence and molecular epidemiology of human coronavirus HKU1 in patients with acute respiratory illness
date: 2012-11-14
journal: J Med Virol
DOI: 10.1002/jmv.23465
sha: 2e0f9b09930e69f7ec6d026acc64a7bd2ef3036f
doc_id: 1053348
cord_uid: snra23sy

In 2005, human coronavirus HKU1 (HCoV‐HKU1) was isolated and identified from a 71‐year‐old man with pneumonia in Hong Kong. To identify and classify genotypes of HCoV‐HKU1 in Korea, a sensitive, specific, and quantitative real‐time polymerase chain reaction (PCR) assay was developed and analyzed the sequences of HCoV‐HKU1 isolated in Korea. A total of 1,985 respiratory specimens taken from patients with acute respiratory illness were tested for HCoV‐HKU1 from January 2007 to May 2008. The major clinical symptoms associated with HCoV‐HKU1 infection were examined statistically and sequence variations of the RNA‐dependent RNA polymerase (RdRp), spike, and nucleocapsid genes were also analyzed. Fifty cases (2.5%) HCoV‐HKU1 were identified by real‐time PCR and viral loads ranged from 6.7 × 10(4) to 1.6 × 10(9) copies/ml. The clinical symptoms of HCoV‐HKU1 infection included rhinorrhea (72%), cough (64%), nasal congestion (56%), fever (32%), sputum (30%), sore throat (18%), chills (16%), postnasal discharge (14%), and tonsillar hypertrophy (10%). There was a seasonal distribution of HCoV‐HKU1 infection, peaking in winter and spring. Both genotypes A and B were detected but no recombination between them was found. This is the first report on the identification and genotyping of HCoV‐HKU1 as a causative agent of acute respiratory illness in Korea. The data suggest that at least two genotypes, A and B, of HCoV‐HKU1 with scattered silent mutations were circulating in Korea from 2007 to 2008. J. Med. Virol. 85:309–314, 2013. © 2012 Wiley Periodicals, Inc.

Human coronavirus (HCoV) 229E and HCoV-OC43 account for 5-30% of cases of acute respiratory tract infection [McIntosh et al., 1969] . HCoV-NL63 was discovered in 2004 as a minor pathogen of HCoV infections in humans [van der Hoek et al., 2004] . A new human coronavirus-HKU1 (HCoV-HKU1) has been reported from a patient with pneumonia in Hong Kong. This species was classified into family Coronaviridae, genus Betacoronavirus, and further as a lineage A [Woo et al., 2005c] .

The genome contains the genes for replicase polyproteins 1a/1b (pp1a/1b), hemagglutinin esterase (HE), spike (S), envelope (E), membrane (M), and nucleocapsid (N) in the same order as in other members of the genus Betacoronavirus [Pyrc et al., 2007; Woo et al., 2009] . HCoV is associated frequently with community acquired upper respiratory tract infections showing typical clinical symptoms including a runny nose, fever, coughing, and wheezing [Labret et al., 2005; Esper et al., 2006; Gose et al., 2011] . However, other disease manifestations including bronchiolitis and pneumonia can occur [Bosis et al., 2007; Gerna et al., 2007; Vabret et al., 2007] . The respiratory symptoms of HCoV-HKU1 infection are usually similar to infections with other HCoVs, but HCoV-HKU1 can cause death in elderly people [Ren et al., 2011] and upper and lower respiratory tract infections in children [Woo et al., 2005b; Jin et al., 2012] .

In this study, a Taq-Man 1 -based real-time polymerase chain reaction (PCR) method that targets the HCoV-HKU1 open reading frame (ORF) 1a and ORF 1b genes with high sensitivity and specificity was developed and evaluated. Sequence information for the RdRp, S, and N genes of HCoV-HKU1 was also analyzed to classify the genotype prevalence of the virus in Korea [Woo et al., 2005a; Woo et al., 2006] . This study is the first to identify and classify the genotypic predisposition of HCoV-HKU1 as a causative agent of acute respiratory illness in Korea.

Throat swab specimens were collected from patients with an acute respiratory illness. Enrollment of 1,985 patients was performed under the acute respiratory tract infection surveillance (ARI-Net) program in the Korea Centers for Disease Control and Prevention (KCDC) from January 2007 to May 2008. All clinical specimens were collected under regulation of Korea National Institute of Health institutional review boards (IRBs) and written informed consents were obtained from patients or their parent/guardian.

To design primer and probes for real-time PCR, full genome sequence of HCoV-HKU1 was retrieved from GenBank accession no. NC006577. HCoV-HKU1 ORF 1a and ORF 1b gene were analyzed using Primer Express software version 2.0 (Applied Biosystems, Foster City, CA) to allocate target sites which could be used as primers and probes. Total RNAs from 140 ml of clinical specimens were extracted using QIAamp viral RNA mini kit (Qiagen, Hilden, Germany) and directly used as templates. cDNA synthesis was performed using random primers and murine mammal leukemia virus (M-MLV) reverse transcriptase (Promega, Madison, WI) in accordance with manufacturer's instruction. Real-time PCR assays targeting the ORF 1a and ORF 1b sequences were performed using 10 ml aliquots of mixtures containing 1 ml cDNA or a standard plasmid, 5 ml Taq-Man 1 universal PCR Master Mix (Applied Biosystems) containing ROX as a passive reference dye, 900 nM of forward and reverse primers and 250 nM of probe (Table I) . Amplification and detection were performed using an ABI PRISM 7900HT sequence detection system (Applied Biosystems) under the following conditions: uracil-N-glycosylase was activated at 508C for 2 min, Taq polymerase was activated at 958C for 10 min and 40 cycles of amplification (958C for 15 sec, 608C for 1 min) were run [Choi et al., 2008; Lu et al., 2008] . Positive values were decided when the copy number exceeded 10 copies/ml and the two target genes (ORF 1a and ORF 1b) were detected simultaneously. The viral load was calculated and determined from mean threshold cycle (Ct) values of specimens and a standard curve prepared from a known concentration of the DNA encoding the RdRp gene. All quantification experiments were performed using triple repetitions.

The characteristics of clinical symptoms of HCoV-HKU1 infection were investigated retrospectively from list of clinical records which were compiled up through ARI-Net. Regression analysis was performed with 95% confidential level to identify relationship between HCoV-HKU1 infection and clinical symptoms using the SAS program (version 9.2).

The RdRp, S, and N genes of HCoV-HKU1 were amplified by RT-PCR with specific primers (Table II) . cDNA synthesis was performed as described above using M-MLV reverse transcriptase (Promega). The PCR assays were performed in a 50 ml reaction volumes containing 5 ml cDNA, 5 ml 10 Â PCRbuffer, 2 ml dNTPmix (final concentration of 400 mM of each dNTP), 1 ml SP Taq (Cosmo Gentech, Seoul, Korea), 1 ml of each primer (10 pM), and nuclease-free water to 50 ml. The PCR reaction was carried out at 948C for 3 min, 35 cycles of amplification (30 sec at 948C; 30 sec at 488C for RdRp, 588C for S, 568C for N; 30 sec at 728C), and a final extension step at 728C for 10 min. The PCR products were run on a 2% agarose gel, stained with SYBR safe DNA gel stain dye (Invitrogen, Carlsbad, CA) and visualized under UV light. 

The following nucleotide accession numbers were used (http://www.ncbi.nlm.nih.gov/genbank/): reference genes; AU884001, AY597011, DQ415899, 

for Detecting HCoV-HKU1

The sensitivity and specificity of this real-time PCR were evaluated. Tenfold serial dilutions of the standard HCoV-HKU1 plasmid DNA ranging from 1 to 1 Â 10 10 copies/ml were used as templates for comparing the sensitivity between real-time PCR and conventional PCR. The detection limit of real-time PCR assays was as low as 10 copies/ml. The sensitivity of the real-time PCR assays used in this study was 1 Â 10 5 times higher than conventional PCR (Fig. 1) . Cross reactivity was also evaluated with cDNA of influenza viruses A/B, parainfluenza viruses 1/2/3, respiratory syncytial virus, human rhinovirus and enterovirus, DNA of adenovirus, human bocavirus, and respiratory pathogenic bacteria (Haemophilus influenzae, Legionella spp., Streptococcus pneumoniae, Mycoplasma pneumoniae, and Chlamydophila pneumoniae). No non-specific amplification was detected in either assay system (data not shown).

Fifty HCoV-HKU1-positive cases were detected (2.5%) out of 1,985 clinical specimens using real-time PCR assays targeting ORF 1a and ORF 1b. The amount of RNA detected in positive specimens ranged from 6.7 Â 10 4 to 1.6 Â 10 9 copies/ml. Based on analysis of multiple alternative clinical records, rhinorrhea (72% of patients), cough (64%), and nasal congestion (56%) were the major symptoms of patients infected with HCoV-HKU1. The other basic respiratory symptoms were fever (32%), sputum discharge (30%), sore throat (18%), chills (16%), postnasal discharge (14%), and tonsillar hypertrophy (10% ; Table III ). There were no positive cases associated with severe respiratory symptoms such as pneumonia and bronchiolitis. Seasonal distribution of HCoV-HKU1 infection was found, with the main prevalence from early winter to late spring with an 8% mean positive rate (Fig. 2) .

RdRp, S, and N genes and the putative amino acid sequences of 11 cases, of which the target genes were amplified by conventional PCR, were analyzed. Eight cases were HCoV-HKU1 genotype B (73%) and three were HCoV-HKU1 genotype A (27%) when queried 1 . Comparison of the sensitivity of the real-time polymerase chain reaction (PCR) assay for the human coronavirus HKU1 (HCoV-HKU1) ORF 1b gene with conventional PCR. The RT-PCR assay detected up to 10 copies/ml of standard DNA, whereas the detection limit of conventional PCR was 1 Â 10 6 copies/ml. through a public database [Benson et al., 2005] . Phylogenetic analysis also indicated that HCoV-HKU1 isolates in this study consisted of genotypes A and B, whereas genotype C was not identified (Fig. 3) .

Since HCoV-HKU1 was first identified in a patient with pneumonia in Hong Kong, clinical manifestations and epidemiological features of infection with this virus have been reported from several countries [Labret et al., 2005; Esper et al., 2006; Gose et al., 2011] . Chung et al. [2007] also tried to identify the virus in patients with acute expiratory wheezing in Korea but evidence of HCoV-HKU1 infection was not found. These negative results arose from several limitations such as an insensitive diagnostic strategy, selection of samples out of season as well as the relatively small number of patients studied. To circumvent these limitations, a quantitative real-time PCR was developed and targeted a study group of 1,985 throat swab specimens collected from patients with an acute respiratory illness from January 2007 to May 2008.

The aim of the real-time PCR assays developed in this study was to detect both of the HCoV-HKU1 ORF 1a and ORF 1b genes. The assays provided enhanced sensitivity and specificity, which could permit quantitative analysis of HCoV-HKU1 in clinical specimens. The ORF 1b assay gave more dynamic Ct values than the ORF 1a assay. Therefore, the ORF 1b assay could be more applicable for clinical diagnosis using respiratory specimens from patients with acute respiratory tract infections caused by HCoV-HKU1.

In this study, HCoV-HKU1 was found in 2.5% of patients (50 cases) with acute respiratory illness sampled from January 2007 to May 2008 in Korea. A similar incidence was reported previously in Hong Kong There was no HCoV-HKU1 genotype C detected in Korea during the study period. Each shadowed number indicates the specimen identification number, and reference sequences are marked as GenBank accession numbers. Phylogenetic trees are shown of (A) the RdRp gene (897 neucleotide (nt)), (B) the S gene (309 nt), and (C) the N gene (425 nt). Primer sequences were excluded from phylogenetic analysis. [Woo et al., 2009] and China [Jin et al., 2010] but an apparently higher rate was noted in Chinese study because of different clinical practices. Unlikely, apparent lower incidence (0.3%, n ¼ 4,181) of HCoV-HKU1 in patients with acute respiratory tract in Hong Kong was reported by Lau et al. [2006] . Relatively high incidence of HCoV-HKU1 in Korea might have resulted from difference of target study group and/or sensitivity of detection strategy. Even though the quantified viral loads varied from 6.7 Â 10 4 to 1.6 Â 10 9 copies/ ml in the respiratory disease specimens used, there were no significant correlations between viral load and specific clinical features. Similarly, there were no statistically significant relationships between the HCoV-HKU1 infection and any particular clinical symptoms. As for seasonal prevalence, the HCoV-HKU1 infections prevailed in the winter and spring in Korea and this result was consistent with the infection rates of HCoV-OC43, 229E and with other reports [Lau et al., 2006; Domiquez et al., 2009; Gaunt et al., 2010] . One limitation of this study was that because clinical specimens were chosen from negative cases of acute respiratory illness in the ARI-Net laboratory surveillance system [Chun et al., 2009] , dual or multiple infections of HCoV-HKU1 with other respiratory viruses could not be detected.

To identify whether there was any recombinant HCoV-HKU1 isolated from Korean patients, the nucleotide sequence information of the RdRp, S, and N genes was analyzed. Similar to previous study reported by Lau et al. [2006] in Hong Kong, our phylogenetic analysis revealed that 11 cases of HCoV-HKU1 isolated in Korea comprised two genotypes: A or B. When single reciprocal nucleotide exchanges were analyzed between the RdRp, S, and N genes, no recombination could be identified. Instead, several point mutations were observed at the nucleotide level, resulting in a transition (C to T) or a transversion (C to T or T to C) in the RdRp, S, and N genes. However, these were silent mutations so that no changes occurred at the amino acid level. Interestingly, a silent continuous mutation in the N gene was identified from 6/8 isolates of HCoV-HKU1 genotype B. Thus, there might have been a unique strain of HCoV-HKU1 prevailing in Korea at that time.

This is the first report that HCoV-HKU1 was associated with acute respiratory illness in Korea. Infections with HCoV-HKU1 peaked from winter through to spring in Korea during 2007-2008. Sequencing and phylogenetic analysis revealed that the 11 HCoV-HKU1 isolates belonged to two distinct genotypes A and B based on the RdRp, S, and N genes. Unlike a previous study, no gene recombination was identified in HCoV-HKU1. Instead, a conserved silent mutation was observed in the N gene of the putative Korean strain of HCoV-HKU1 genotype B from six specimens. Intensive molecular epidemiological investigations would be required to elucidate any possible changes in viral tropism caused by such mutations.

Authors are grateful to Mr. H.D. Cheong for excellent statistical analysis of clinical symptoms.

Coronavirus HKU1 in an Italian pre-term infant with bronchiolitis

Development of real-time PCR assays for detection and quantification of human bocavirus

Establishing a surveillance network for severe lower respiratory tract infections in Korean infants and young children

Detection of viruses identified recently in children with acute wheezing

Detection of four human coronaviruses in respiratory infections in children: A oneyear study in Colorado

Coronavirus HKU1 infection in the United States

Epidemiology and clinical presentations of the four human coronaviruses 229E, HKU1, NL63, and OC43 detected over 3 years using a novel multiplex real-time PCR method

Human respiratory coronavirus HKU1 versus other coronavirus infections in Italian hospitalized patients

Coronavirus HKU1 in children, Brazil

Prevalence and clinical characteristics of human CoV-HKU1 in children with acute respiratory tract infections in China

Newly identified respiratory viruses associated with acute lower respiratory tract infections in Lanzou

Detection of the new human coronavirus HKU1: A report of 6 cases

Coronavirus HKU1 and other coronavirus infections in Hong Kong

Development and comparison of real-time and conventional RT-PCR assay for detection of human coronavirus NL63 and HKU1

Seroepidemiologic studies of coronavirus infection in adults and children

The novel human coronaviruses NL63 and HKU1

Prevalence of human coronaviruses in adults with acute respiratory tract infections in Beijing

Human (non-severe acute respiratory syndrome) coronavirus infections in hospitalized children in France

Identification of a new human coronavirus

Characterization and complete genome sequence of a novel coronavirus, coronavirus HKU1, from patients with pneumonia

Clinical and molecular epidemiological features of coronavirus HKU1-associated community-acquired pneumonia

Phylogenetic and recombination analysis of coronavirus HKU1, a novel coronavirus from patients with pneumonia

Comparative analysis of 22 coronavirus HKU1 genomes reveals a novel genotype and evidence of natural recombination in coronavirus HKU1

More and more a coronaviruses: Human coronavirus HKU1