Detection of a divergent Parainfluenza 4 virus in an adult patient with influenza like illness using next-generation sequencing
© Bialasiewicz et al.; licensee BioMed Central Ltd. 2014
Received: 4 February 2014
Accepted: 11 April 2014
Published: 19 May 2014
Human Parainfluenza viruses are a common cause of both upper and lower respiratory tract infections, particularly in children. Of the four Parainfluenza virus serotypes, Parainfluenza 4 is least well characterised from both the clinical, epidemiological and genetic perspectives.
Flocked nose or throat swabs from a previous study investigating viral prevalence in community-based adults suffering from influenza like illness were used as the basis for this study. Samples in which no virus was detected using a 16 viral respiratory pathogen real-time PCR panel were barcoded and pyrosequenced using the Roche 454 GS FLX Titanium chemistry. The sequences were analysed using the VirusHunter bioinformatic pipeline. Sanger sequencing was used to complete the detected Parainfluenza 4 coding region.
A variant Parainfluenza 4 subtype b strain (QLD-01) was discovered in an otherwise healthy adult who presented with influenza like illness. Strain QLD-01 shared genomic similarities with both a and b subtypes. The extent of divergence of this genome from the 5 available whole Parainfluenza 4 genomes impacted the predicted binding efficiencies of the majority of published Parainfluenza 4 PCR assays.
These findings further support a possible role for Parainfluenza 4 in the aetiology of adult respiratory disease within the community setting, and highlight the caution needed to be used in designing PCR assays from limited sequence information or in using proprietary commercial PCR assays.
KeywordsParainfluenza 4 Community infection Respiratory tract infection Influenza like illness Next generation sequencing Virus discovery PCR False negative results Adult
Human parainfluenza viruses are a common cause of both upper and lower respiratory tract infections, particularly in children [1–5]. Four serotypes are known, but most epidemiological and clinical research has been focused on parainfluenza serotypes 1–3. This has been primarily due to the poor growth characteristics in cell culture of parainfluenza 4 (PIV4), the lack of commercial diagnostic reagents, and historical exclusion from routine diagnostic testing . Two antigenically distinct PIV4 subtypes, PIV4a and PIV4b, exist . Functionally and epidemiologically, little is known about the two PIV4 subtypes, however both are capable of co-circulating within the same population . With the advent of reverse transcription PCR (RT-PCR), it has become easier to screen for an expanded range of RNA viruses, leading to a re-examination of PIV4’s epidemiology and role in human disease [2, 3, 5, 8]. The majority of research into PIV4 has focused on children within the hospital setting, however little information is available on the role of PIV4 in disease within the broader community.
Despite the resurgent interest in PIV4, a dearth of publically available PIV4 sequences, and in particular, whole genomes, still exists. The lack of appreciable sequence information hampers the design and evaluation of sensitive research and diagnostic RT-PCR assays, since these tests are reliant on oligomer homology to the target sequence. Thus it is imperative to increase the number of publically available sequences for clinically relevant pathogens where little information is currently available, particularly when variant genomes are observed.
Powerful new techniques such as next-generation sequencing have been applied to clinical samples over the last five years with the aims of discovering novel pathogens. Numerous viruses and variant strains have been identified using this approach, including a divergent PIV4 subtype a isolate late in 2013 . Unlike insensitive traditional virological methods and highly specific RT-PCR, next-generation sequencing methods have the advantage of being able to sequence total or targeted DNA and RNA from samples in an unbiased way, without a priori knowledge of the possible viral agent(s) present, thus making them the ideal tool for novel and divergent viral genome discovery.
In this study, we used a combination of RT-PCR and next-generation sequencing to identify and characterise the full coding sequence of a novel PIV4 variant from an adult participating in a community-based cohort study of respiratory illness. Furthermore we compared the primer sequences of existing RT-PCR assays to the genome of this variant strain.
Samples used for this study were derived from a previous study  investigating viral prevalence in a sub-population of participants in a community-based, randomised control trial assessing influenza vaccine effectiveness. Briefly, flocked nose or throat swabs were collected from otherwise healthy adults aged from 18–64 years who presented with influenza-like illness (ILI). ILI was defined as cough, sore throat, runny nose or nasal congestion and at least one systemic symptom (fever greater/equal than 37.8°C, feverishness, chills or myalgia). In total, 643 samples were screened for adenovirus, human metapneumovirus, parainfluenza viruses 1, 2 & 3, respiratory syncytial virus, influenza A and B, picornaviruses, bocavirus, coronaviruses (OC43, 229E, NL63 and HKU1) and WU and KI polyomaviruses using real-time PCR . 299 study samples from which no viral pathogen was detected were used as templates for novel virus discovery. Written consent was obtained from all study participants. The original and current studies were approved by the Royal Children’s Hospital Human Research Ethics Committee (Melbourne) and the Human Research Protection office of Washington University, respectively.
Viral discovery pipeline
Total nucleic acid was extracted from each sample, subjected to sample-specific barcoded random-priming cDNA synthesis and then PCR amplified using barcode-specific primers. Standard library construction and 454 GS FLX Titanium pyrosequencing was performed as previously described . The sequences were analysed using the VirusHunter bioinformatic pipeline . In brief, high quality reads with similarity to viruses at the nucleotide level or amino acid level were identified using BlastN and BlastX, respectively.
Genome sequencing & assembly
Individual 454 reads were assembled and mapped against PIV4 subtype a and PIV4 subtype b reference genomes M-25 (AB543336) and SKPIV4 (EU627591), respectively. Walking primers for amplifying and sequencing the remainder of the genome were designed based on the assembled contigs and reference genomes (see Additional file 1). Sanger sequencing was performed on the overlapping amplified cDNA bidirectionally. Contig assembly and genome characterisation was performed using CLC Bio Genomics Workbench 6.5 software. (CLC Bio, Denmark) Phylogenetic analysis was performed using MEGA 5.2 software (http://www.megasoftware.net/) . Recombination event analyses were performed using the Recombination Analyses Tool software (https://github.com/ethering/RAT) , with widow sizes of 1719, 800, and 400 applied to a whole genome alignment using QLD-01 as the reference.
A sample collected in Melbourne (Australia) from a 49 year old male produced 32 reads (see Additional file 2) with highest similarity to PIV4 (83.7-98.3%). Apart from the presentation with ILI meeting the study case definition, no other clinical information was recorded in the subject’s symptom diary. The reads mapped to six singletons and five contigs of 533, 607, 700, 759, and 904 nt in length, and showed highest similarity to SKPIV4, against which the genome walking primers were designed. Final assembly of the sequenced amplicons yielded a near-complete genome of 17,160 nt in length, and included the entire coding region. The sequence was submitted to GenBank as isolate QLD-01 (KF908238). During the study period, one other PIV4 detection was observed within the study population.
Strain QLD-01 and the recently described divergent PIV4a isolate DK(459) shared several structural features identified by Alquezar-Planas et al.; an extended C-terminal end and a 13 residue stretch within the globular head of the predicted HN protein, as well as a 57 nt insertion at the 3′ leader non-coding sequence. These features were not unique to the two isolates, and in the case of the HN protein, were more common than not across both genotypes. Conversely, QLD-01 contained the most divergent residue sequence (91.4-95.0%) of the PIV4b genotype within the variable C terminal end (residues 411–451) of the Nucleocapsid protein.
PCR target conservation
PIV4 PCR oligomer mismatches to isolate QLD-01
Improvements in sequencing and detection technologies over the past 15 years have led to increasing detection rates of existing, neglected, and unknown pathogens. The existence of PIV4 has been known since 1959; however only recently has PIV4 been more appreciated as a respiratory pathogen in its own right through the use of modern molecular methods [1–5, 8]. Modern molecular diagnostics heavily rely on PCR-based techniques. Because of their high specificity, these same methods are susceptible to decreased sensitivity or even false negative results when confronted with even minor changes in target sequences. This limitation is particularly relevant to clinically important pathogens for which little sequence data are available to guide PCR assay development and ongoing evaluation.
The variant PIV4 isolate described in this study is the most divergent of the five available whole genome sequences. Thus it is not surprising that four of the six evaluated published PCR assays contained potentially deleterious mismatches with the isolate. In particular, all four assays contain mismatches at or near their primers’ 3′ ends, which are especially sensitive to incorrect base pairing and would potentially lead to decreased primer binding efficiency, and in conjunction with the other primer mismatches, false negative results. These mismatches illustrate the difficulty in designing sensitive PCR assays based on very limited sequence information. Additionally, the use of commercial assays for which primer sequences are not readily available, such as those used in recent PIV4 epidemiological studies [1, 8], should be used with the understanding that there is no capacity to evaluate the assays’ target sequence conservation as new data on emerging variant viral strains becomes available.
Isolate QLD-01 was found in the upper respiratory tract of an adult with ILI but no other known conditions. It was the sole virus detected, despite extensive screening for other known and unknown respiratory pathogens, thus supporting a potential role as the aetiological agent of the subject’s ILI symptoms. Recent studies have reported PIV4 infections associated with both lower and upper respiratory tract symptoms within the hospital setting [1–5, 8]. This study provides further evidence of PIV4’s possible involvement in upper respiratory tract infections in otherwise healthy adults within the community setting.
In regard to its genomic structure, QLD-01’s highest similarity was to isolate SKPIV4. However its overall genomic and NP C-terminal end divergence separate it from SKPIV4 and other PIV4b genomes. In other parainfluenza viruses, the NP C terminus binds the protein-associated viral RNA to the RNA polymerase , thus QLD-01’s variant C-terminal end may alter the isolate’s viral RNA synthesis kinetics. Isolate QLD-01’s shared genomic features with both its closest PIV4b homologue and the divergent PIV4a raises the possibility of recombination events occurring between the viral subtypes, however no conclusive evidence was observed to support this hypothesis.
In this study, the utility of a combined RT-PCR and next-generation sequencing approach to identifying novel viral pathogen was demonstrated with the discovery of a variant strain of PIV4. The whole coding region of the variant strain was sequenced and showed that the majority of publically available PIV4 PCR assays contained mismatches when aligned to this variant, which may lead to decreased sensitivity and false negative results, thereby underestimating the prevalence of PIV4.
This work was supported in part by: NIH grant U54 AI057160 to the Midwest Regional Center of Excellence for Biodefense and Emerging Infectious Disease Research, ARC Linkage Grant LP0989464, and QCMRI Program grant 50006. SB is supported by the QCMRI Early Career Research Fellowship 50056.
- Frost HM, Robinson CC, Dominguez SR: Epidemiology and Clinical Presentation of Parainfluenza Type 4 in Children: A 3-Year Comparative Study to Parainfluenza Types 1–3. J Infect Dis. 2013, 209: 695-702.View ArticlePubMedPubMed CentralGoogle Scholar
- Ren L, Gonzalez R, Xie Z, Xiong Z, Liu C, Xiang Z, Xiao Y, Li Y, Zhou H, Li J, Yang Q, Zhang J, Chen L, Wang W, Vernet G, Paranhos-Baccalà G, Shen K, Wang J: Human parainfluenza virus type 4 infection in Chinese children with lower respiratory tract infections: a comparison study. J Clin Virol. 2011, 51: 209-212. 10.1016/j.jcv.2011.05.001.View ArticlePubMedGoogle Scholar
- Aguilar JC, Pérez-Breña MP, García ML, Cruz N, Erdman DD, Echevarría JE: Detection and identification of human parainfluenza viruses 1, 2, 3, and 4 in clinical samples of pediatric patients by multiplex reverse transcription-PCR. J Clin Microbiol. 2000, 38: 1191-1195.PubMedPubMed CentralGoogle Scholar
- Chiu CY, Rouskin S, Koshy A, Urisman A, Fischer K, Yagi S, Schnurr D, Eckburg PB, Tompkins LS, Blackburn BG, Merker JD, Patterson BK, Ganem D, DeRisi JL: Microarray detection of human parainfluenzavirus 4 infection associated with respiratory failure in an immunocompetent adult. Clin Infect Dis. 2006, 43: e71-e76. 10.1086/507896.View ArticlePubMedGoogle Scholar
- Lau SKP, To W, Tse PWT, Chan AKH, Woo PCY, Tsoi H, Leung AFY, Li KSM, Chan PKS, Lim WWL, Yung RWH, Chan K, Yuen K, Icrobiol JCLINM: Human Parainfluenza Virus 4 Outbreak and the Role of Diagnostic Tests. J Clin Microbiol. 2005, 43: 4515-4521. 10.1128/JCM.43.9.4515-4521.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Canchola J, Vargosko AJ, Kim HW, Parrott RH, Christmas E, Jeffries B, Chanock RB: Antigenic variation among newly isolated strains of parainfluenza type 4 virus. Am J Hyg. 1964, 79: 357-364.PubMedGoogle Scholar
- Lau SKP, Li KSM, Chau K, So L, Lee RA, Lau Y, Chan K, Lim WWL, Woo PCY, Yuen K: Clinical and molecular epidemiology of human parainfluenza virus 4 infections in hong kong: subtype 4B as common as subtype 4A. J Clin Microbiol. 2009, 47: 1549-1552. 10.1128/JCM.00047-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu W-K, Liu Q, Chen D-H, Liang H-X, Chen X-K, Huang W-B, Qin S, Yang Z-F, Zhou R: Epidemiology and clinical presentation of the four human parainfluenza virus types. BMC Infect Dis. 2013, 13: 28-10.1186/1471-2334-13-28.View ArticlePubMedPubMed CentralGoogle Scholar
- Alquezar-Planas DE, Mourier T, Bruhn CA, Hansen AJ, Vitcetz SN, Mørk S, Gorodkin J, Nielsen HA, Guo Y, Sethuraman A, Paxinos EE, Shan T, Delwart EL, Nielsen LP: Discovery of a divergent HPIV4 from respiratory secretions using second and third generation metagenomic sequencing. Sci Rep. 2013, 3: 2468-View ArticlePubMedPubMed CentralGoogle Scholar
- Howard PF, McCaw JM, Richmond PC, Nissen M, Sloots T, Lambert SB, Lai M, Greenberg M, Nolan T, McVernon J: Virus detection and its association with symptoms during influenza-like illness in a sample of healthy adults enrolled in a randomised controlled vaccine trial. Influenza Other Respi Viruses. 2013, 7: 330-339. 10.1111/j.1750-2659.2012.00395.x.View ArticleGoogle Scholar
- Holtz LR, Finkbeiner SR, Zhao G, Kirkwood CD, Girones R, Pipas JM, Wang D: Klassevirus 1, a previously undescribed member of the family Picornaviridae, is globally widespread. Virol J. 2009, 6: 86-10.1186/1743-422X-6-86.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao G, Krishnamurthy S, Cai Z, Popov VL, da Rosa AP T, Guzman H, Cao S, Virgin HW, Tesh RB, Wang D: Identification of Novel Viruses Using VirusHunter – an Automated Data Analysis Pipeline. PLoS One. 2013, 8: e78470-10.1371/journal.pone.0078470.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.View ArticlePubMedPubMed CentralGoogle Scholar
- Etherington GJ, Dicks J, Roberts IN: Recombination Analysis Tool (RAT): a program for the high-throughput detection of recombination. Bioinformatics. 2005, 21: 278-281. 10.1093/bioinformatics/bth500.View ArticlePubMedGoogle Scholar
- Wang CYT, Arden KE, Greer R, Sloots TP, Mackay IM: A novel duplex real-time PCR for HPIV-4 detects co-circulation of both viral subtypes among ill children during 2008. J Clin Virol. 2012, 54: 83-85. 10.1016/j.jcv.2012.01.013.View ArticlePubMedGoogle Scholar
- Reijans M, Dingemans G, Klaassen CH, Meis JF, Keijdener J, Mulders B, Eadie K, van Leeuwen W, van Belkum A, Horrevorts AM, Simons G: RespiFinder: a new multiparameter test to differentially identify fifteen respiratory viruses. J Clin Microbiol. 2008, 46: 1232-1240. 10.1128/JCM.02294-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Templeton KE, Scheltinga SA, Beersma MFC, Kroes ACM, Claas ECJ: Rapid and sensitive method using multiplex real-time PCR for diagnosis of infections by influenza a and influenza B viruses, respiratory syncytial virus, and parainfluenza viruses 1, 2, 3, and 4. J Clin Microbiol. 2004, 42: 1564-1569. 10.1128/JCM.42.4.1564-1569.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Van de Pol AC, van Loon AM, Wolfs TFW, Jansen NJG, Nijhuis M, Breteler EK, Schuurman R, Rossen JW: Increased detection of respiratory syncytial virus, influenza viruses, parainfluenza viruses, and adenoviruses with real-time PCR in samples from patients with respiratory symptoms. J Clin Microbiol. 2007, 45: 2260-2262. 10.1128/JCM.00848-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang S, Chen L, Zhang G, Yan Q, Yang X, Ding B, Tang Q, Sun S, Hu Z, Chen M: An amino acid of human parainfluenza virus type 3 nucleoprotein is critical for template function and cytoplasmic inclusion body formation. J Virol. 2013, 87: 12457-12470. 10.1128/JVI.01565-13.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/14/275/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.