Scientifica / 2012 / Article

Research Article | Open Access

Volume 2012 |Article ID 871201 | 4 pages |

Comparative Analysis of the Full-Length Genome Sequence of a Clinical Isolate of Human Parainfluenza Virus 4B

Academic Editor: M. Clementi
Received27 May 2012
Accepted02 Jul 2012
Published19 Jul 2012


We are engaged in airborne transmission and epidemiology studies of respiratory pathogens, with particular interest in human parainfluenza virus type 4 (hPIV-4) and other lesser studied viruses. In this paper, hPIV-4 was detected in primary rhesus monkey kidney (PRMK) cells that had been inoculated with nasopharyngeal swab material obtained from a child with a mild upper respiratory tract illness. Attempts to isolate the virus in pure culture were hampered by the presence of a fast-growing simian spumavirus that was a contaminant of the PRMK cells. Total RNA was extracted from the PRMK cell culture, and PCR followed by sequencing of a subgenomic section of the fusion protein gene suggested the hPIV-4 was subtype 4B. At the time of this work, two complete but dissimilar hPIV-4B genomes had been deposited by others in GenBank. To gain better insights on hPIV-4B, and to test methods that we are developing for viral forensics, the entire genomic sequence of our virus was determined from archived RNA. The hPIV-4B genomic sequence that we determined conforms to the paramyxovirus “rule of six.” Here, we compare and contrast the genetic features of the three completely sequenced hPIV-4B genomes currently present in GenBank.

Human parainfluenza viruses (hPIVs) are single-stranded, negative sense RNA viruses of the genus Rubulavirus, family Paramyxoviridae, which cause acute respiratory tract infections in children and adults. Four hPIV serotypes (hPIV 1–4) have been identified; serotype 4 is further subdivided into two antigenic subtypes: 4A and 4B [1, 2]. The epidemiology and clinical manifestations of hPIV 1–3 are well known, whereas comparatively little is known about hPIV4s, as they are difficult to isolate in cell culture and are absent from routine respiratory virus detection tests in most clinical virology laboratories [35]. Whereas hPIV4s were formerly mostly associated with mild respiratory illnesses in young people, recent studies indicate the viruses can cause more severe infections such as pneumonia in young and older patients (mentioned in [35]).

The genomic cRNAs of hPIV-4 subtypes 4A and B are a little >17.0 kbp in length. Their viral genomes encode for nucleocapsid (NP), phospho (P), nonstructural (V), matrix (M), fusion (F), haemagglutinin-neuraminidase (HN), and large (L) proteins. Prior to this work, there were two complete hPIV-4B sequences in GenBank: those of strains 68–333 [6] and SKPIV-4 [7].

Primary monkey kidney (PMK) cells are inoculated with appropriate specimens for the detection of human parainfluenza viruses in many American diagnostic microbiology laboratories. The PMK cells available to these diagnostic laboratories are usually harvested from one of various Chlorocebus or Asian macaque species and contain a mixture of kidney cell-types. Furthermore, the PMK cells can contain endogenous simian viruses that are either latent in the kidneys or cause persistent but inapparent kidney infections in their hosts. Their presence in PMK cultures generally becomes evident after the cells are maintained in culture for more than a few days. Regardless, experience has shown that the probability of detecting human parainfluenza viruses in clinical specimens through in vitro virus culture is better with PMK cells other than cell lines commonly used in diagnostic virology laboratories.

The virus analyzed in this work was from an immunocompetent two-year-old child in Chicago with a mild upper respiratory infection of two-days duration at the time of specimen collection (October 2004). At the time of specimen collection, the patient’s symptoms included runny nose, barky cough, low fever, and decreased appetite. A nasopharyngeal swab specimen from the patient was eluted in universal virus transport medium (BD, NJ, USA), and equal aliquots of the solubilized material inoculated into A549, MDCK, WI38, and rhesus PMK cells and inoculated at 35°C. The PMK cell-culture media contained antibodies against PIV5 and SV40. The cultures tested negative by direct immunofluorescence assays (DFA) at 24 and 72 hrs p.i. using a commercial kit that detects PIV-1, -2, -3, influenza A and B viruses, adenovirus, and RSV (Respiratory Panel 1 DFA kit, Millipore, Billerica, MA, USA). However, with FITC-labeled anti-PIV-4 antibody (catalog item no. 5034, Millipore), sporadic PMK (but not the other) cells were borderline positive at 24 hr and positive at 72 hr p.i., demonstrating characteristic punctuate intracytoplasmic staining. Unfortunately, large vacuoles and widespread cell deterioration were evident in about 30% of the PMK cells by 72 hrs p.i. (including the negative controls), suggesting that a contaminating virus was present in the PMK cultures. Aliquots were therefore taken from the hPIV-4B-infected PMK culture and inoculated into NCI-292, Vero, LLC-MK2, or CV-1 cells, in hopes of isolating the hPIV-4 virus in cells not susceptible to the contaminating virus. Thereafter, an RNA stabilizing solution (RNAlater, Ambion, Austin, TX, USA) was added to PMK cells, total RNA purified as described previously [8], and the RNA archived at −80°C. Attempts to isolate hPIV-4 were not successful; the contaminant, identified as a Group VI spumavirus (foamy retrovirus) (data not shown), caused extensive CPE (large vacuoles) 24 hrs after inoculation of the NCI-292, Vero, LLC-MK2, or CV-1 cells, and all the cultures were terminated.

Two-step reverse transcription PCR of the archived RNA with primers Para4-F (5′-catgggtgtcaaaggtttatc-3′) and Para4-R (5′-tgctgctgtaacttgtgcagc-3′) amplified a 376-base pair (bp) section of the HPIV-4 F gene [8]. Sequencing of the amplicon revealed the virus was probably hPIV-4B. As a complete genomic sequence of hPIV-4B was not available for comparison in 2004, and our priorities were focused on other viruses, further analyses were postponed until an opportune time was available for the development of sequencing strategies appropriate for hPIV-4B.

We revived our sequencing efforts after two independently-derived hPIV-4B sequences were deposited in GenBank. For our work, targeted hPIV-4B sequences were RT-PCR-amplified from the archived RNA using a genome walking approach. Overlapping primers described in [6, 7] and others purpose-designed by us for our tasks were used for PCR amplification and sequencing. Superscript II reverse transcriptase (Life Technologies) was used for first-strand cDNA synthesis in the presence of SUPERase-In RNase inhibitor (Ambion), and high fidelity Platinum Taq DNA polymerase (Life Technologies) was used for PCR. The 3′ and 5′ ends of the viral genome were determined from vRNA using a RACE (rapid amplification of cDNA ends) kit (RLM RACE, Ambion, Austin, TX) following the manufacturer’s instructions. Of note, efforts for determining the 3′ and 5′ end sequences of the viral genome were laborious; these tasks are simpler using viral genomes arising from purified virus particles. Sequences were directly analyzed using an Applied Biosystem 3130 DNA analyzer by using BigDye Terminator (v. 3.1) chemistry and the same oligonucleotide primers used for amplifications. The virus sequence was designated hPIV-4B 04-13 (“04-13” signifies “unusual” isolate no. 13 of year 2004).

The complete hPIV-4B 04-13 cRNA is 17,304 bp and thus conforms to the paramyxovirus “rule of six” since it is divisible by 6. A full-genome BLAST analysis reveals 98% homology with hPIV-4B strain 68–333 and 97% homology with hPIV-4B strain SKPIV-4. Key genetic features of hPIV-4B strains SKPIV-4, 68–333, and 04-13 are given in Table 1.

Virus designationSKPIV-468–33304-13

Genbank accession numberEU627591AB543337JQ241176.1
cRNA genome length17,361 bp17,304 bp17,304 bp
Leader sequenc>55 nt?55 nt55 nt
First 20 nt of viral genome complimentary to terminal 20 nt?NoYesYes
Putative 3′ peptide?nt 52–75nt 52–75
ORF 1 NP551 aa551 aa551 aa
nt 155–1,810nt 101–1,756nt 101–1,756
CTAAGAT in 1st intergenic region YesYesYes
ORF 2 399 aa399 aa399 aa
nt 2,096–3,293nt 2,041–3,238nt 2,041–3,238
229 aa229 aa229 aa
nt 2,096–2,785nt 2,041–2,730nt 2,041–2,730
CAGAAGTA in 2nd intergenic regionYesYesYes
ORF 3 M382 aa382 aa382 aa
nt 3,589–4,737nt 3,531–4,670nt 3,531–4,670
AGAGCCA AATAC in 3rd intergenic region
ORF 4 F543 aa543 aa543 aa
nt 5,232–6,863nt 5,174–6,805nt 5,174–6,805
CTATTAT follows poly A in 4th intergenic regionYesYesYes
ORF 5 HN579 aa574 aa574 aa
nt 7,563–9,302nt 7,506–9,230nt 7,506–9,230
ORF 6 L2,279 aa2,279 aa2,279 aa
nt 10,025–16,864nt 9,970–16,809nt 9,970–16,809

eader sequence as defined by Komada et al. [6].

The deduced amino acid sequences of the F, H-N, L, M, NP, and P from hPIV-4B isolate 04-13 were aligned against the homologous sequences from 7 other members of the genus Rubulavirus and 2 members of the genus Avulavirus. Sequence alignments were performed using Mafft 5.8 [9] followed by minor manual adjustments in ClustalW [10]. The E-INS-I alignment strategy was used with the following parameters: scoring matrix (BLOSUM62), gap open penalty (1.53), and offset value (0). For each gene alignment, the sequence was trimmed to the first conserved amino acid at the 5′ and 3′ ends prior to analyses. To assess gene concordance, Bayesian analyses were performed independently for each gene. Phylogenetic trees were constructed using MrBayes v. 3.1.2 [11]. A mixed prior was used on amino acid models and default priors for topology and branch lengths. The Markov chain was run for a maximum of 10 million generations, with a stopping rule implemented so that the analysis would halt when the average deviation of the split frequencies was <0.001%. Four independent analyses were conducted, each with 1 cold and 3 heated chains with the default heating parameter (temperature = 0.2). Every 50 generations were sampled and the first 25% of MCMC samples discarded as burn-in.

Preliminary phylogenetic analysis revealed that there was only a single significant incongruence among individual gene trees (defined by the presence of incompatible bipartitions that received a posterior probability of >90%, resp.). The incongruence involved the M-gene analysis that supported an alternate branching pattern for the mumps virus as has been previously observed [11]. Therefore, for the final analysis, we concatenated the sequences for the 6 genes into 1 matrix. The dataset contained 4442 amino acid characters (including gaps) for 10 viral taxa. The concatenated 6-gene Bayesian analysis demonstrated with a high level of confidence that hBIV4B isolate 04-13 is most closely related to hPIV-4B isolate 68–333 with hPIV-4B isolate SKPIV4 as the sister group to the other two isolates (Figure 1). Human parainfluenza virus 4A was found to be the sister group to the hPIV-4B isolates. The hPIV4 clade was found to be the sister group to a second Rubulavirus clade composed of mumps virus, simian virus 5, simian virus 41, and human parainfluenza virus 2.

The results of our genomic level phylogenetic analysis are consistent with previous analyses of the genus Rubulavirus [6, 7, 12]. As pointed out by Yea et al. [7], the genome of SKPIV-4 does not follow the paramyxovirus “rule of six.” Theirs is not a sequencing error; paramyxovirus genomes that violate the rule are occasionally encountered ([7], and J. Lednicky, unpublished). It will be informative henceforth to determine if the presence of “aberrant length” hPIV-4B genomes (i.e., those genomes whose length is not divisible by 6) in virus isolated from sick individuals correlates with clinical presentation and also whether genomic alterations occur as a consequence of passage of hPIV-4B in the primate cells used for the detection of the viruses by diagnostic laboratories.


The authors thank Dr. James F. Wellehan DVM, M.S., Ph.D., DACZM, DACVM, College of Veterinary Medicine, University of Florida, Gainesville, for advice over phylogenetic analyses of paramyxoviruses and for critical review of this paper.


  1. J. Canchola, A. J. Vargosko, H. W. Kim et al., “Antigenic variation among newly isolated strains of parainfluenza type 4 virus,” American Journal of Epidemiology, vol. 79, no. 3, pp. 357–364, 1964. View at: Google Scholar
  2. R. A. Lamb and G. D. Parks, “Paramyxoviridae: the viruses and their replication,” in Fields Virology, D. M. Knipe, P. M. Howley, D. E. Griffin et al., Eds., pp. 1449–1496, Wolters Kluwer/Lippincott Williams & Wilkins, Philadelphia, Pa, USA, 5th edition, 2007. View at: Google Scholar
  3. S. K. P. Lau, K. S. M. Li, K. Y. Chau et al., “Clinical and molecular epidemiology of Human parainfluenza virus 4 infections in Hong Kong: subtype 4B as common as subtype 4A,” Journal of Clinical Microbiology, vol. 47, no. 5, pp. 1549–1552, 2009. View at: Publisher Site | Google Scholar
  4. L. Ren, R. Gonzalez, Z. Xie et al., “Human parainfluenza virus type 4 infection in Chinese children with lower respiratory tract infections: a comparison study,” Journal of Clinical Virology, vol. 51, no. 3, pp. 209–212, 2011. View at: Publisher Site | Google Scholar
  5. M. L. Vachon, N. Dionne, E. Leblanc, D. Moisan, M. G. Bergeron, and G. Boivin, “Human parainfluenza type 4 infections, Canada,” Emerging Infectious Diseases, vol. 12, no. 11, pp. 1755–1758, 2006. View at: Google Scholar
  6. H. Komada, M. Kawano, A. Uefuji et al., “Completion of the full-length genome sequence of Human parainfluenza virus types 4A and 4B: sequence analysis of the large protein genes and gene start, intergenic and end sequences,” Archives of Virology, vol. 156, no. 1, pp. 161–166, 2011. View at: Publisher Site | Google Scholar
  7. C. Yea, R. Cheung, C. Collins, D. Adachi, J. Nishikawa, and R. Tellier, “The complete sequence of a Human parainfluenza virus 4 genome,” Viruses, vol. 1, no. 1, pp. 26–41, 2009. View at: Google Scholar
  8. T. C. Rubinas, R. B. Carey, M. C. Kampert, S. Alkan, and J. A. Lednicky, “Fatal hemorrhagic pneumonia concomitant with Chlamydia pneumoniae and parainfluenza virus 4 infection,” Archives of Pathology and Laboratory Medicine, vol. 128, no. 6, pp. 640–644, 2004. View at: Google Scholar
  9. J. P. Huelsenbeck, F. Ronquist, R. Nielsen, and J. P. Bollback, “Bayesian inference of phylogeny and its impact on evolutionary biology,” Science, vol. 294, no. 5550, pp. 2310–2314, 2001. View at: Publisher Site | Google Scholar
  10. K. Katoh, K. I. Kuma, H. Toh, and T. Miyata, “MAFFT version 5: improvement in accuracy of multiple sequence alignment,” Nucleic Acids Research, vol. 33, no. 2, pp. 511–518, 2005. View at: Publisher Site | Google Scholar
  11. A. J. McCarthy and S. J. Goodman, “Reassessing conflicting evolutionary histories of the Paramyxoviridae and the origins of respiroviruses with Bayesian multigene phylogenies,” Infection, Genetics and Evolution, vol. 10, no. 1, pp. 97–107, 2010. View at: Publisher Site | Google Scholar
  12. J. D. Thompson, D. G. Higgins, and T. J. Gibson, “CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice,” Nucleic Acids Research, vol. 22, no. 22, pp. 4673–4680, 1994. View at: Google Scholar

Copyright © 2012 John A. Lednicky et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1226 Views | 509 Downloads | 2 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder