Multiplex PCR analysis of clusters of unexplained viral respiratory tract infection in Cambodia

Background Fevers of unknown origin constitute a substantial disease burden in Southeast Asia. In majority of the cases, the cause of acute febrile illness is not identified. Methods We used MassTag PCR, a multiplex assay platform, to test for the presence of 15 viral respiratory agents from 85 patients with unexplained respiratory illness representing six disease clusters that occurred in Cambodia between 2009 and 2012. Results We detected a virus in 37 (44%) of the cases. Human rhinovirus, the virus detected most frequently, was found in both children and adults. The viruses most frequently detected in children and adults, respectively, were respiratory syncytial virus and enterovirus 68. Sequence analysis indicated that two distinct clades of enterovirus 68 were circulating during this time period. Conclusions This is the first report of enterovirus 68 in Cambodia and contributes to the appreciation of this virus as an important respiratory pathogen.


Introduction
Acute respiratory infections (ARI) are among the leading causes of morbidity and mortality in both developed and developing nations [1][2][3]. In Cambodia, pneumonia is the third highest cause of death in children under five [4,5]. In December 2006, a passive hospital-based surveillance study was initiated to identify the causes of acute undifferentiated fever in patients seeking healthcare in Cambodia. Between 2006 and 2009, 9,997 patients from peri-urban and rural Cambodian health centers were enrolled in the study. The most frequent diagnoses included influenza, dengue, malaria and leptospirosis [6]. No diagnosis was obtained in 62% of cases [6]. In a separate analysis using samples from December 2006 to December 2008, 4,233 patients with respiratory disease were screened for influenza virus by real-time reverse-transcriptase polymerase chain reaction (rRT PCR). Of these patients, 1,151 (27.2%) were positive for influenza [7]. In a separate study, Sreng et al. analyzed 2,805 patients with influenza-like illness in Cambodia between 2006 and 2008 and showed that only 9.6% tested positive for influenza [8]. To address the high rate of undiagnosed ARI cases in Cambodia, we selected six clusters of previously unexplained respiratory illness from 2009 to 2012 for analysis by MassTag PCR, a high throughput multiplex screening platform. We identified a potential viral pathogen in 44% (37/85) of the samples, confirming the utility of multiplex diagnostic assays in clinical microbiology and providing insights into ARI in Cambodia.

Results
From 2009 to 2012, 17,363 samples were collected from 18 health clinics in six provinces of Cambodia and screened for influenza virus, dengue virus and Plasmodium spp (Figure 1). From amongst the samples that tested negative for these agents, we selected 6 clusters representing 85 patients with unexplained acute respiratory illness for further analysis. 88% (75/85) of the patients exhibited influenza-like symptoms ( Table 1). The archived specimens from the 6 clusters were screened by MassTag PCR for viral respiratory pathogens [9,10]. At least one virus was detected in 44% (37/85) samples ( Figure 2A). Infection with two viruses was detected in 8 samples. Rhinovirus (HRV) was the most frequently detected virus (n = 17), followed by respiratory syncytial virus (RSV; n = 15), and enterovirus (HEV; n = 9) ( Figure 2B).

KD02 cluster
The KD02 cluster samples were collected from 20 patients in Kandal who presented for medical attention over an eight-day period in August, 2009. Eighteen had influenza-like illness (ILI), one had a presumptive diagnosis of typhoid fever and one was diagnosed with pneumonia ( Table 1). Analysis of these samples by MassTag PCR detected a viral agent in 11 samples (Table 2). HEV was detected in six samples; one of them also contained RSV-B. To determine the HEV serotype, we amplified and sequenced a 600 nucleotide (nt) VP4-VP2 genome fragment from each sample. Sequence analysis indicated that all six HEV-positive samples contained a single strain of enterovirus 68 (EV-D68) (100% nt identity in VP4-VP2). Unlike most HEVs, EV-D68 is predominately a respiratory pathogen. Recent genetic analysis of EV-D68 indicated three main clades of the virus circulating worldwide [11]. To determine the clades circulating in Cambodia, we amplified and sequenced a 339 nt fragment of VP1 from all six samples. The Cambodian EV-D68 strain was most similar to EV-D68 strains originating from Asia between 2006 and 2010 (99% identity to strains from Japan, Philippines and China) and clustered with strains belonging to clade A, which has been implicated in recent EV-D68 outbreaks worldwide [11][12][13][14] (Figure 3). Comparison of amino acids within the VP1 relative to the EV-D68 Fermon reference strain revealed the majority of the substitutions were in the putative surface-exposed BC and DE loops ( Figure 4). All six EV-D68 positive samples originated from adults (four men, two women) ranging between 24 to 71 years   of age. Five of the six patients presented with pharyngitis, headache, cough and malaise. The exception, a 36-year-old male, presented with abdominal cramps, and joint pain, and was initially diagnosed with typhoid fever. Since both symptoms are atypical of EV-D68 infection, they may have been caused by a co-infecting agent. Five other samples from the cluster contained a HRV; amplification and sequencing of the VP4-VP2 genome fragment revealed that four of these were HRV-C and one was HRV-A (Table 2). Three HRV-positive samples originated from adults and two from children.
Testing revealed that five cases were positive for RSV-B, all in patients 12 years old and under. One sample was also co-infected with HRV-C (Table 2). A single sample contained HRV-A. Of the six adult cases, one contained EV-D68. Sequence analysis indicates this virus (CSP-204) was similar to the strain from the Kandal cluster with 98.7% and 98.6% nt identity in the VP4-VP2 and VP1, respectively (100% amino acid (aa) identity in both proteins) (Figures 3 and 4).

RK01 cluster
Between July 25 th and September 12 th of 2011 a cluster of 11 children less than five years of age from Ratanakiri received a diagnosis of pneumonia. Ten of them had shortness of breath; eight had nausea and vomiting. MassTag PCR detected RSV-A in nine cases (82%). Four RSV-A positive samples were co-infected with other viruses, identified as, coxsackie virus B4 (CV-B4), human coronavirus (hCoV-OC43), HRV-C, or EV-D68 (Table 1). Sequence analysis of a VP1 fragment amplified from this EV-D68 strain indicated that this strain was only 86% identical (92% aa identity) to the Cambodian strains from 2009, and clustered with EV-D68 sequences from clade B (Figures 3 and 4). The highest similarity was to strains reported from the Philippines in 2008 (99% nt identity in the VP1) [11,15,16].

KD02 (2010), ST02, and KT02 clusters
The remaining three clusters comprising 32 cases all originated during the dry season in Cambodia (November through May). In cluster KD02, 10 of 11 patients were initially diagnosed with pharyngitis. All 11 patients in this cluster reported headache, cough and sore throat. ADV and HRV-B represented the lone viral agents detected in this cluster. In cluster ST02 all 11 cases were initially diagnosed with malaria and complained of fever (100%), cough (73%) and chills (63.5%). HRV-B was detected in three samples, and one sample contained CoV-OC43. In cluster KT02, 70% of the cases were diagnosed with influenza; the remaining 30% of cases were diagnosed with malaria (Table 1). HRV was detected in four of these samples (3 HRV-A and 1 HRV-C); CoV-229E was present in one sample.

Discussion
Our analysis of 6 clusters of previously undiagnosed ARI from Cambodia revealed the presence of potential viral pathogens in 44% of cases. The agents most commonly detected were HRV, RSV, and HEV. Our results are similar to published data from a Cambodian children cohort study wherein the predominant agents found were HRV and RSV [17]; but differ in that we also report findings in adults wherein the predominant organisms found were HRV and EV-D68.
RSV is a frequent cause of respiratory disease in children and was the predominant agent detected in pediatric clusters in Ratanakiri (RK01; RSV-A) and Kampong Speu (KS02; RSV-B) respectively. The two clusters of RSV occurred during the rainy season suggesting that RSV infections may represent a significant cause of viral pneumonia in Cambodian children during this period [17][18][19].
EV-D68 was the primary agent detected in a cluster initially diagnosed as pharyngitis/influenza and was also detected in individual cases from two other clusters. Our study represents the initial report linking EV-D68 with respiratory disease in Cambodia. First isolated in the 1960s, EV-D68 is atypical among EVs in respiratory tract tropism [20,21]. Since the late 2000′s, there has been a global increase of respiratory disease outbreaks associated with EV-D68. EV-D68 has been reported in outbreaks of respiratory disease in Asia (China, Thailand, Japan, the Philippines), Europe (Italy, Finland, the Netherlands, the UK, France), Oceania (New Zealand), Africa (the Gambia, Senegal, and South Africa) and the United States [11][12][13]15,16,[22][23][24]. The clinical presentation  of EV-D68 infections in these outbreaks ranged from acute, mild illness to severe pneumonia, and, in rare instances, death. Retrospective studies have shown epidemic emergence of EV-D68 in some countries, leading to the hypothesis that evolution of the virus may have increased its fitness, enhancing its global expansion [11,13]. In our study we found that multiple clades of EV-D68 were circulating in Cambodia between 2009 2011, similar to reports from other countries [12,13,16]. We detected representatives of two distinct EV-D68 clades, both with variations in the putative antigenic sites relative to the EV-D68 reference strain. It has been suggested that changes within the antigenic sites of the capsid may have enabled the recent worldwide emergence of EV-D68 [25]. The majority of EV-D68 positive cases in our study originated from adults. Although EV-D68 was originally isolated from children, recent reports link EV-D68 to ARI in both children and adults, though the severity of symptoms may differ [13,26].
During the preparation of this manuscript, a large outbreak of respiratory disease in children linked to EV-D68 has occurred in the US [27]. This outbreak, which has resulted in large number of hospitalized cases, represents the first documented large-scale outbreak of disease due to EV-D68.

Conclusions
Our findings demonstrate the utility of multiplex diagnostic assays for differential diagnosis of respiratory tract infections and underscore the increasing importance of EV-68 in respiratory disease.

Sample collection
Beginning in December 2006, hospital-based surveillance of acute febrile illness was conducted in six provinces of eastern Cambodia: Kandal, Kampong Speu, Kratie, Ratanakiri, Stung Treng and Svay Rieng (Figure 1) [6]. Over a six-year period subjects were enrolled at 18 health clinics within these provinces. Ten were located in Kandal and Kampong Speu, approximately 40 to 75 Km from the capital Phnom Penh, respectively, whereas the others were located in more remote areas with difficult access during rainy seasons.
Between 2009 and 2012, 17,363 patients were enrolled and tested by PCR for influenza virus and dengue virus (upper respiratory swab (URS) and serum, respectively) and by Giemsa stain of whole blood for Plasmodium spp. Inclusion/Exclusion criteria, testing algorithm, specimen processing and testing were as described by Kasper et al. However the catchment area of the study changed as per Figure 1 [6]. Samples were collected under the auspices of study protocols #NAMRU2.2005.0004 and #037 NECHR approved by the U.S. Naval Medical Research Center Institutional Review Board and The Cambodian Ethical Committee for Health and Research, respectively. All research was in compliance with applicable federal regulations governing the protection of human subjects. All enrolled patients provided written consent. For patients <18 years of age, consent was obtained from a parent or legal guardian.

Cluster definition
Patients enrolled at the same site within a two-week span with similar clinical symptoms were included in individual clusters. We selected five such case clusters for further analysis. Two clusters originated from Kandal, and one each from Kampong Speu, Steung Treng, and Kratie provinces ( Figure 1). Additionally, we selected a cluster of children under five years of age who were diagnosed at the Ratanakiri (RK01) field site with pneumonia. Patients in this cluster were enrolled within a 50-day period. Together the six clusters represented 85 patients.

MassTag PCR
Total nucleic acid was extracted from 250 μl of URS using the EasyMag automated extraction platform (Bio-M? rieux). Total nucleic acid (TNA) was eluted in 40 μl. cDNA was generated using 10 μl of TNA and Super-Script II Reverse Transcriptase kit (Life Technologies). All samples were screened by MassTag PCR using respiratory panels composed of 15 known respiratory viruses [9,10]. The panels included primer sets for detection of 14 RNA viruses (influenza A and B, human respiratory syncytial virus A and B (RSV-A, RSV-B), human parainfluenza virus 1, 2, 3, 4A and 4B, human metapneumovirus, human coronavirus (hCoV) OC43 and 229E, human enterovirus (HEV)/rhinovirus (HRV)), and 1 DNA virus (adenovirus (ADV). All samples positive by Mas-sTag PCR were confirmed by single agent PCR and dideoxy sequencing of the resulting amplification products.

Data analysis
Mass Tag PCR data were analyzed using non-parametric statistics. The confidence score was calculated using fold inter-quartile range distance of the positive control or the test sample from the 95th percentile cut-point of the negative control distribution for a tag. The score is similar to parametric z-score but is distribution free.

PCR and dideoxy sequencing
The list of primers can be found in Table 3. For each reaction, quantified positive control standards, samples, and non-template controls were used. In HEV-positive samples the serotype was identified using a consensus nested PCR assay that amplifies a 600-nucleotide region within the VP4-2 gene. The cycling conditions used for MassTag PCR were also used in individual PCR assays with the exception of the HEV VP4/2 nested assay (annealing temperature of 55? C).
All EV-D68 sequences were deposited in GenBank under accession numbers KJ556320 KJ556335. All base pair coordinates are provided relative to the original 1962 EV-D68 isolate, the Fermon strain (GenBank accession no. AY426531). A Maximum-likelihood phylogenetic tree of EV-D68 based on a 339 nucleotide fragment of the VP1 gene was constructed using Mega 5.2 software.  Submit your next manuscript to BioMed Central and take full advantage of: