key: cord-0736264-x6736prq
authors: Soonnarong, Rapeepun; Thongpan, Ilada; Payungporn, Sunchai; Vuthitanachot, Chanpim; Vuthitanachot, Viboonsuk; Vichiwattana, Preeyaporn; Vongpunsawad, Sompong; Poovorawan, Yong
title: Molecular epidemiology and characterization of human coronavirus in Thailand, 2012–2013
date: 2016-08-26
journal: Springerplus
DOI: 10.1186/s40064-016-3101-9
sha: 9ef23d86bba041d9fc573624cdf17e94a30fb8f6
doc_id: 736264
cord_uid: x6736prq

BACKGROUND: Coronavirus causes respiratory infections in humans. To determine the prevalence of human coronavirus (HCoV) infection among patients with influenza-like illness, 5833 clinical samples from nasopharyngeal swabs and aspirates collected between January 2012 and December 2013 were examined. RESULTS: HCoV was found in 46 (0.79 %) samples. There were 19 (0.32 %) HCoV-HKU1, 19 (0.32 %) HCoV-NL63, 5 (0.09 %) HCoV-229E, and 3 (0.05 %) HCoV-OC43. None of the sample tested positive for MERS-CoV. The majority (54 %) of the HCoV-positive patients were between the ages of 0 and 5 years. HCoV was detected throughout the 2-year period and generally peaked from May to October, which coincided with the rainy season. Phylogenetic trees based on the alignment of the spike (S) gene sequences suggest an emergence of a new clade for HCoV-229E. CONCLUSIONS: The data in this study provide an insight into the prevalence of the recent circulating HCoVs in the region.

Human coronavirus (HCoV) is an enveloped positivestrand RNA virus with ~27-32 kb genome. Its genome contains five major open reading frames (Orfs) encoding the replicase polyproteins (Orf 1a and Orf 1b), spike (S), envelope (E), membrane (M), and nucleocapsid (N) proteins (Hilgenfeld and Peiris 2013) . HCoV belongs to the order Nidovirales, family Coronaviridae and the genus Coronavirus (Lai and Cavanagh 1997) . Virions are pleomorphic with diameters between 60 and 220 nm. There are four genera of coronaviruses (Alpha-, Beta-, Gammaand Deltacoronavirus) based on serologic and phylogenetic characterization. Among the five recognized human coronaviruses, HCoV strains 229E and NL63 are alphacoronaviruses, while HCoV strains OC43, HKU1, and SARS (responsible for the severe acute respiratory syndrome) are betacoronaviruses (Gaunt et al. 2010; Gonzalez et al. 2003; Drosten et al. 2003) . In addition, a novel genotype of coronavirus was identified in Saudi Arabia in 2012, later referred to as MERS-CoV (for the Middle East respiratory syndrome coronavirus), which was isolated from the sputum of a 60-year-old man in Saudi Arabia who presented acute pneumonia and renal failure with a fatal outcome (Hilgenfeld and Peiris 2013; Zaki et al. 2012) . It was subsequently found to be most closely related to the bat betacoronavirus HKU4 and HKU5 (Gonzalez et al. 2003) . Fortunately, the conserved structure and function of the polymerase made it possible to develop a PCR assay based on the polymerase gene region to accurately differentiate new coronaviruses as was done with MERS-CoV and other novel coronaviruses (Rota et al. 2003; Moës et al. 2005) .

HCoVs are recognized as one of the most frequent causes of upper respiratory tract infection in elderly adults leading to acute pneumonia and renal failure (Memish et al. 2013; Gorse et al. 2009 ). HCoV infection is thought to participate in the exacerbations of chronic obstructive pulmonary disease, congestive heart failure, and other chronic diseases necessitating emergency care and long-term hospitalization. In addition to respiratory syncytial virus, parainfluenza virus, adenovirus, and influenza virus, HCoVs are also known to cause the common cold especially among children (Moës et al. 2005; Lu et al. 2012) . Furthermore, HCoVs are associated with mild to severe upper and lower respiratory tract illness and can cause more serious respiratory diseases in children, the elderly and people with underlying conditions (McIntosh et al. 1974; van Elden et al. 2004) .

The study of coronaviruses has sometimes been difficult due to limitations in cell culture and serology. Thus, epidemiological and viral prevalence data are valuable in investigating the emergence of HCoV infection. Using reverse-transcription polymerase chain reaction (RT-PCR) and phylogenetic analysis, we characterized HCoVs identified in Thai patients with respiratory tract infection between 2012 and 2013.

This study was approved by the Faculty of Medicine of Chulalongkorn University (IRB 388/56). A total of 5833 clinical samples from patients with influenza-like illness were obtained between January 2012 and December 2013 for routine testing of respiratory viruses as part of an epidemiological surveillance (Auksornkitti et al. 2014; Prachayangprecha et al. 2013; Sriwanna et al. 2013) . From these, 5196 nasopharyngeal swabs (NPS) were collected from patients with upper respiratory tract infections who sought treatment at the Bangpakok 9 International Hospital and Chum Phae Hospital (Khon Kaen, Thailand). Additionally, 637 nasopharyngeal aspirates (NPA) were collected from patients with lower respiratory tract infections who were admitted to King Chulalongkorn Memorial Hospital and Chon Buri Hospital. All samples were stored in the viral transport medium (prepared according to The World Health Organization guideline). Patient identifiers were removed and all samples were anonymous, although data on demographics, symptoms, history of illness, results of clinical examination and laboratory investigations were retained for analysis.

Viral DNA/RNA was extracted using HiYield Viral Nucleic Acid Extraction Kit (RBC Bioscience, Taipei, Taiwan, ROC). Complementary DNA was synthesized with random hexameric primers and ImProm-II (Promega, Madison, WI, USA) according to the manufacturer's instructions. HCoV was identified using the semi-nested RT-PCR to amplify the RNA-dependent RNA polymerase (RdRp) gene of HCoV and the N gene of MERS-CoV as previously described (Corman et al. 2012; Kon et al. 2012) . For controls, a panel of viral nucleic acid extracted from samples previously tested positive for influenza A virus [subtype H1N1 (pandemic 2009), H3N2], influenza B virus, respiratory syncytial virus (RSV), human parainfluenza virus (HPIV) and human adenovirus was used to determine the specificity of the HCoV semi-nested RT-PCR. PCR was performed using the PerfectTaq TM Plus Master Mix kit (5 PRIME, Hamburg, Germany) and primers (Table 1) in 25 µL reaction volume under the following conditions: initial denaturation at 94 °C for 3 min, 40 cycles of 94 °C for 30 s, 55-58 °C for 30 s, 72 °C for 40 s, and a final extension at 72 °C for 7 min. Expected amplicons for RdRp and N were ~450 and ~260 bp, respectively. All the specimens were positive for GAPDH, which served as an internal control.

To identify the types of HCoVs (HCoV-229E, HCoV-NL63, HCoV-OC43, and HCoV-HKU1), the S gene was genotyped using nested RT-PCR. Amplification condition was the same as above except that the annealing temperature was 50 °C. Amplicons were agarose gel-purified using the Expin Gel SV kit (GeneAll, Seoul, Korea) and DNA sequencing was performed by First BASE Laboratories (Seri Kembangan, Selangor, Malaysia).

PCR products from the amplification of HCoV-229E, HCoV-NL63, HCoV-OC43, HCoV-HKU1 and MERS-CoV were cloned into pGEM-T Easy Vector System (Promega, CA, USA) according to the manufacturer's instructions. RNA transcripts from these five coronaviruses were used as standards to validate assay sensitivity. The sensitivity of the semi-nested RT-PCR assay was established for each coronavirus by testing transcripts of known concentrations from serial dilution.

Nucleotide sequences of the RdRp and S genes were edited using Chromas Lite (version 2.1.1) and compared to the HCoV reference strains available in GenBank using the Basic Local Alignment Search Tool (BLAST) program (www.ncbi.nlm.nih.gov/Blast.cig). Multiple alignments of the nucleotide sequences utilized Clustal W in BioEdit (version 7.0.9). Phylogenetic trees were constructed using the neighbor-joining method and the Kimura two-parameter distance model (MEGA version 6.06) and evaluated by 1000 bootstrap pseudo-replicates.

The nucleotide sequences of the RdRp gene (accession numbers KJ866056-KJ866101) and the S gene (accession numbers KJ866102-KJ866147) identified in this study were deposited in the GenBank database.

Serial dilutions of coronavirus RNA transcripts (HCoV-229E, HCoV-NL63, HCoV-OC43, HCoV-HKU1, and MERS-CoV) were tested using semi-nested RT-PCR. RNA transcripts for all virus types were detectable at ≤10 copies per reaction, or <400 copies/mL of sample. To evaluate detection specificity and possible cross-reactivity, nucleic acid from 18 different types of respiratory viruses were tested. All assay results were negative, and no false positive was observed.

Among the 5833 samples analyzed, 637 (10.9 %) were positive for influenza A virus, 206 (3.5 %) were positive for influenza B virus, 201 (3.4 %) were positive for respiratory syncytial virus A, 91 (1.6 %) were positive for respiratory syncytial virus B, and 78 (1.3 %) were positive for adenovirus (Table 2 ). More importantly, 46 samples (0.79 %) tested positive for HCoV. Co-infection with other respiratory viruses was not observed in the HCoV-positive samples. All samples were negative when tested for MERS-CoV N gene.

Amongst the HCoV-positive samples, 56.5 % were from men and 43.5 % were from women (gender ratio 1.3:1) (Table 3) . HCoV infection was detected in all age groups, and the mean age of HCoV-infected patients was 21.37 years (min. = 4 months, max. = 93 years, mean = 27.04 years). The percentage of HCoV infection per year was 0.81 % (23/2838) in 2012 and 0.77 % (23/2995) in 2013. There were no seasonal peaks associated with HCoV and no positive samples were identified in the typically dry months of April and November during both years (Fig. 1) . The majority of HCoV (25/46 or 54 %) were detected mainly in young children between the ages of 0-5 years and 15 % in the elderly aged over 60 years (Fig. 2) . Approximately 80 % of the positive samples (37/46) were isolated from the NPS samples, while the rest (9/46) were from the NPA samples. Furthermore, 46 % (17/37) of all positive NPS and 89 % (8/9) of all positive NPA samples belonged to the 0-5 year age group.

To further differentiate the 46 HCoV-positive samples, the S gene was sequenced. Analysis showed that these HCoV strains belonged to one of the four HCoV species (Table 2; Fig. 3 ). In all, 19 (0.32 %) were positive for HCoV-HKU1, 19 (0.32 %) were positive for HCoV-NL63, 5 patients (0.09 %) were positive for HCoV-229E, and 3 (0.05 %) were positive for HCoV-OC43. Relative to all HCoV-positive samples, therefore, the predominant genotypes were 41.3 % for both HCoV-HKU1 and HCoV-NL63 (19/46), followed by 11 % HCoV-229E (5/46) and 6.5 % HCoV-OC43 (3/46). Interestingly, HCoV-NL63 and HCoV-HKU1 appeared sporadically during the study period and were detected mainly in March 2012 and July 2013, respectively. 

Of the 46 HCoV-positive samples, 38 (82.6 %) were derived from patients with acute infections involving the upper respiratory tract such as the nose, sinuses, pharynx or larynx (Table 4 ). The remaining 8 (17.4 %) were from patients with lower respiratory tract infections. The most common clinical manifestations were fever and rhinorrhea, although several patients experienced tachypnea and hypoxemia.

To assess the relationship among the strains identified, partial nucleotide sequences of the S gene from HCoVpositive samples were subjected to phylogenetic analysis (Fig. 4) . This region was examined as it represents the antigenic determinant of the virus and correlates with strain evolution. HCoV-HKU1 strains identified in this study formed a single cluster (~99.6 % sequence identity) within clade B together with the strain identified on mainland China. Meanwhile, HCoV-NL63 strains clustered in two different groups within the same clade. Their sequences were genetically closest to the isolates from The Netherlands and U.S.A. Three HCoV-OC43 strains grouped together with previous viruses characterized in Japan in clade C. Finally, although the 5 HCoV-229E identified in this study appeared to cluster into a new group, they demonstrated nearest genetic resemblance to strains from Japan and China.

Due to the awareness of MERS-CoV and other novel human coronavirus infections worldwide, this study aimed to characterize possible circulation and prevalence of HCoV in Thailand in a 2-year period. The HCoV prevalence of 0.79 % found in this study was lower than that from previous reports, which ranged from 2. (1) the subjects included in this study were adult patients whereas other studies examined mainly hospitalized children (Talbot et al. 2009) Although a previous study in southern Thailand found that HCoV-OC43 was the predominant strain (Suwannakarn et al. 2014 ), very few HCoV-OC43 infections were detected in this study. Yearly fluctuation in the predominance of HCoV strains in Thailand has been documented (Dare et al. 2007 ). The majority of HCoVs detected in this study was HCoV-NL63 and HCoV-HKU1. HCoV-NL63 was found mainly in 2012 and HCoV-HKU1 in 2013, while HCoV-OC43 and HCoV-229E were sporadic and less frequent.

HCoV infections elicited common symptoms such as fever, cough, headache, vomiting, muscle pain and sore throat. Upper respiratory tract infection such as rhinorrhea, sputum, and lower respiratory tract infection such as tachypnea and abnormal breath sounds were frequently observed in our patients similar to past observations (Ren et al. 2011; Dominguez et al. 2009 ). However, clinical presentations were slightly different among each subtype of HCoV. In this study, vomiting was more common in HCoV-229E, HCoV-NL63 and HCoV-HKU1-positive patients but not in those infected by HCoV-OC43. Lower respiratory tract infections were associated with individuals with HCoV-OC43, HCoV-NL63, and HCoV-HKU1, but not in those with HCoV-229E. Tachypnea, hypoxemia, abnormal breathing sound and signs of respiratory distress also occurred less often for HCoV-229E. These clinical features were very similar to those observed with other respiratory viruses, such as RSV, HPIV, and human metapneumovirus.

The simultaneous screening for common respiratory viruses allowed us to investigate possible co-infection of HCoV with other viral pathogens. Although parainfluenza viruses, enteroviruses, and human metapneumovirus were not tested, co-infection of HCoV with RSV was never observed (Theamboonlers et al. 2007 ). It may be that HCoV and other respiratory viruses do not share the same seasonal circulation. The possibility of viral coinfection with bacteria, however, cannot be ruled out. Nevertheless, confirmation of the mechanism of co-infection will require good clinical sampling and further longitudinal studies conducted over several years.

Previous characterization of the genetic variation and evolution of HCoV have utilized the whole viral genome, S, RdRp, or N gene. Here, the phylogenetic analysis of HCoV was performed based on the sequences of partial spike gene because the S gene may be more discriminating than the RdRp gene (Lau et al. 2006 (Lau et al. , 2011 . While phylogenetic grouping of HCoV-HKU1, HCoV-NL63, and HCoV-OC43 were characterized by clades or clusters, HCoV-229E strains were typically subdivided by the year of emergence. HCoV-229E strains collectively appeared to form a new cluster similar to but distinct from previous strains in the phylogenetic tree.

URTI upper respiratory tract infection, LRTI lower respiratory tract infection (pneumonia, bronchiolitis, and acute infection of the pulmonary parenchyma) a Percent of HCoV-positive individuals with such symptoms b Age-related definition of tachypnea by The World Health Organization for individuals <2 months (>60 breaths/min), 2-12 months (>50 breaths/min), 1-5 years (>40 breaths/min) and ≥5 years (>20 breaths/min) c Oxygen saturation (SpO2) level below 95 % d Including retractions (subcostal, intercostal, suprasternal), nasal flaring, and grunting e Crepitation, rhonchi, and wheezing

HCoV-229E HCoV-OC43 HCoV-NL63 HCoV-HKU1 , and HCoV-229E (d) with MEGA 6.06 software using the distance method and the neighbor-joining algorithm with Kimura 2 parameters. Phylogenetic clusters are indicated to the right of the trees. CU designation is given to each of the 46 clinical isolates examined in this study (black dots). The node clusters were supported by bootstrap values >70 %

Characterization of human coronavirus etiology in Chinese adults with acute upper respiratory tract infection by real-time RT-PCR assays

Coronavirus infection in acute lower respiratory tract disease of infants

Family cluster of Middle East respiratory syndrome coronavirus infections

A novel pancoronavirus RT-PCR assay: frequent detection of human coronavirus NL63 in children hospitalized with respiratory tract infections in Belgium

Detection and typing by molecular techniques of respiratory viruses in children hospitalized for acute respiratory infection in

Epidemiology of seasonal influenza in Bangkok between 2009 and 2012

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

Characterization of a novel coronavirus associated with severe acute respiratory syndrome

Molecular characterization of human adenovirus infection in Thailand

Prevalence and genetic characterization of human coronaviruses in southern Thailand from

Coronavirus infection and hospitalizations for acute respiratory illness in young children

Human coronavirus infection among children with acute lower respiratory tract infection in Thailand

Frequent detection of human coronaviruses in clinical specimens from patients with respiratory tract infection by use of a novel real-time reverse-transcriptase polymerase chain reaction

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

In summary, this study surveyed the prevalence and clinical presentations of different HCoV infections in Thai patients from 2012 to 2013. Although the numbers of HCoV-positive samples were low, we detected four species of HCoV and phylogenetically characterized their diversity. These viruses appeared to continue to cause infections globally, thus accurate and timely diagnosis will be essential. The assay described in this study can be used to assist in the rapid and accurate detection of emerging HCoV infection.

HCoV: human coronavirus; NPA: nasopharyngeal aspirate; NPS: nasopharyngeal swab; RSV: respiratory syncytial virus; RT-PCR: reverse-transcription polymerase chain reaction.Authors' contributions RS, IT and PV performed the experiments. SP and YP designed the research study. CV and VV collected the samples. RS, SV, and YP wrote the manuscript. All authors read and approved the final manuscript.