Human parainfluenza virus type 2 hemagglutinin-neuramindase gene: sequence and phylogenetic analysis of the Saudi strain Riyadh 105/2009
Virology Journal volume 9, Article number: 316 (2012)
Although human parainfluenza type 2 (HPIV-2) virus is an important respiratory pathogen, a little is known about strains circulating in Saudi Arabia.
Among 180 nasopharyngeal aspirates collected from suspected cases in Riyadh, only one sample (0.56%) was confirmed HPIV-2 positive by nested RT-PCR. The sample that was designated Riyadh 105/2009 was used for sequencing and phylogenetic analysis of the most variable virus gene; the haemagglutinin-neuramindase (HN). Comparison of HN gene of Riyadh 105/2009 strain and the relevant sequences available in GenBank revealed a strong relationship with Oklahoma-94-2009 strain. Phylogenetic analysis indicated four different clusters of HPIV-2 strains (G1-4). Twenty-three amino acid substitutions were recorded for Riyadh 105/2009, from which four are unique. The majority of substitutions (n=18) had changed their amino acids characteristics. By analyzing the effect of the recorded substitutions on the protein function using SIFT program, only two located at positions 360 and 571 were predicted to be deleterious.
The presented changes of Riyadh 105/2009 strain may possess potential effect on the protein structure and/or function level. This is the first report that describes partial characterization of Saudi HPIV-2 strain.
Human parainfluenza type 2 virus (HPIV-2) is one of the important viruses that infect the human respiratory tract causing croup and pneumonia. It is an enveloped non-segmented negative single-stranded RNA virus that belongs to genus Rubulavirus of the family Paramyxoviridae. The virus genome is approximately 15 kb long and encodes seven structural proteins. The nucleocapsid (N), phosphoprotein (P), and large polymerase (L) encapsidate the genomic RNA, while hemagglutinin-neuraminidase (HN), fusion (F), and matrix (M) are envelope-associated proteins.
HN protein constitutes the prominent virus fringes responsible for binding to sialoconjugate receptors (hemagglutinating activity), and for enzymatic destruction of the receptors (neuramindase activity). Moreover, the stalk region of HN was proposed to contain the site which promotes cell fusion by F protein[3, 4]. Extensive genetic and antigenic variations among HPIV-2 strains were attributed to HN gene diversity[5, 6]. Two distinct clusters were suggested by phylogenetic analysis of HPIV-2 strains on the basis of complete HN gene sequence. However, this phylogenetic tree was inadequately informative due to the low number of analyzed sequences. In this report, the complete HN sequence of first Saudi HPIV-2 strain was analyzed and compared to the corresponding sequences available in GenBank. A phylogenetic tree was constructed to determine the predictable origin of Saudi strain and the relationship between different strains circulating globally.
Nasopharyngeal aspirates of 180 children, hospitalized with acute respiratory disease in King Khalid University Hospital, Riyadh, were obtained as described before. Sample collection was approved by the Ethical Committee of King Saud University. The aspirates were taken from patients, after getting an informed consent from their parents. Samples were chilled transported in Minimal Essential Medium supplemented with 500 U penicillin and 500 μg streptomycin per ml (Sigma, St. Louis, MO) to Virology Research Laboratory, King Saud University. After pulse-vortexing and centrifugation at 3000 rpm for 10 min, 140 μl aliquots of the clarified supernatants were utilized for RNA extraction using QIAamp viral RNA extraction kit (Qiagen, Hilden, Germany) according to the manufacturer`s guidelines. Detection of HPIV-2 RNA was accomplished by amplification of 508 bp fragment of HN gene using One-step RT-PCR kit (Qiagen) and the primer set HPIV2-cDNA/HPIV2-rev (Table1). Positive reactions were confirmed by the amplification of an internal 204 bp fragment of the PCR amplicon in a nested PCR using TopTaq Master Mix Kit (Qiagen), and the primer set HPIV2-nF/HPIV2-nR (Table1). Products of both reactions were analyzed, along with 100 bp DNA ladder (Qiagen), in 1.5% agarose gel stained with ethidium bromide.
Only one sample (0.56%), collected from two-year old female child, was confirmed positive for HPIV-2. The sample, designated Riyadh 105/2009, was used for sequencing and phylogenetic analysis of the complete HN gene. Sequencing strategy involved the amplification of three overlapping fragments, that flank the entire HN gene sequence, using FideliTaq™ RT-PCR Master Mix (GE Healthcare, Buckinghamshire, UK) and sequencing primer sets (Table1). Specific RT-PCR products were purified from agarose gel using Illustra GFX PCR DNA Kit (GE Healthcare), and were sequenced on both strands (GenArt, Rosensburg, Germany). Nucleotide sequences were edited and assembled using Bioedit Biological Sequence Alignment Editor (Ibis Biosciences, Carlsbad, CA). Complete nucleotide and deduced amino acid sequences of HN gene of Riyadh 105/2009 strain were deposited in GenBank under the accession number HM460888. Multiple alignments with 17 HPIV-2 complete HN sequences available in GenBank (Table2) were performed using the Clustal W algorithm of Megalign program, Lasergene software. Divergence analysis and tree construction were conducted using the same software, with 1,000 bootstrap replicates. The effect of non-synonymous mutations on the protein function was predicted using SIFT program (http://sift.bii.a-star.edu.sg).
As expected, the results of sequence analysis revealed a single ORF of 1716 nucleotides that encodes a polypeptide of 571 amino acids. Multiple sequence comparisons between Riyadh 105/2009 strain and the different global strains showed that the homology ranges from 91.8% to 99.1% and from 93.9% to 99% on the levels of nucleotide and deduced amino acids, respectively. A close relationship between the Saudi strain and Oklahoma-94-2009 strain (GenBank Accession number JF912194) was observed with only 16 nucleotide (0.9%) and 4 amino acid (0.7%) variations. This relationship may propose a common ancestral origin of both viruses. Seven unique nucleotide substitutions (0.41%) were found in Riyadh 105/2009 strain (G53A, A54G, T143C, T148C, T1327C, T1383, and G1401). As a result, four amino acids (0.53%) at positions 18, 48, 50, and 443 were changed.
Phylogenetic analysis of the different worldwide strains of HPIV-2 indicated four different clusters (G1-4). Clusters G1 and G4 were further subdivided into subclusters a and b (Figure1). Both Riyadh 105/2009 and Oklahoma-94-2009 strains were grouped together in cluster G3. The nucleotide and deduced amino acid sequence comparisons, between Riyadh 105/2009 strain and the different strains of each group, further established this classification. For instance, the homology between Riyadh 105/2009 and G1 viruses ranged from 91.8-92.7% (nucleotides), and 93.9-94.8% (amino acids), while for G2 and G4 viruses, the homology was 93–93.2%, and 96.2-96.4% on the nucleotide level, and 95.1-95.3%, and 96–96.2% on the amino acid level, respectively. This Phylogenetic differentiation correlates well with the tree constructed by Terrier et al., in which two clusters were defined. The first cluster included G1 and G2 strains of the conducted phylogram, while the second contained strains of G4 cluster. G3 cluster strains were not included in that study due to their recent identification. It may be concluded that new variant strains are continuing to emerge, and therefore, further updating of this phylogram may be required upon identification of more strains on spatial and temporal basis.
Analysis of the deduced amino acid sequence alignment identified 23 positions of potential variation between Riyadh 105/2009 and other international strains (Table2). Nineteen of these positions were shared with the other G3 virus (Oklahoma-94-2009), while 9–10, and 0–2 positions were only common with G4, and G1-2 viruses, respectively. The majority of the recorded substitutions (n=18) changed the amino acid characteristics on the level of charge, polarity, and/or hydropathy. However, the effect of these substitutions on the antigenic properties of the mature HN protein could not be concluded, since positions of the antigenic sites are not yet determined for HPIV-2. Aligning the positions of amino acid substitution of Riyadh 105/2009 strain with their counterparts in HPIV-3 HN protein identified a single residue (Q345R) that may be located in a distinct antigenic site.
Studying the effect of amino acid change on the protein function using SIFT program, outlined that all of the substitutions are well tolerated; except the change of asparagine (N) to tyrosine (Y) at position 360, and leucine (L) to proline (P) at position 571 that were predicted to affect the protein function with a score of 0.02 and 0.00, respectively. These two residues locate in the globular head region that carries both haemagglutinating and neutamindase activities. However, the extent of their effect on the protein function would be correctly concluded when the three-dimensional structure of HN protein becomes available, and the amino acids crucial for both activities are determined.
In conclusion, the complete HN gene sequence of the first HPIV-2 Saudi strain (Riyadh 105/2009) was analyzed and compared to the corresponding sequences/strains at GenBank. Phylogenetic analysis of these strains allowed allocation of four separate clusters G1-G4. Saudi strain was more relevant to Oklahoma-94-2009 strain, which may propose a common ancestral origin. The majority of the amino acid substitutions observed in Riyadh 105/2009 strain changed the properties of their respective amino acids. However, the effect of these substitutions on the protein structure and function requires more detailed studies.
Hall CB: Respiratory syncytial virus and parainfluenza virus. N Engl J Med 2001, 344: 1917-1928. 10.1056/NEJM200106213442507
Karron RA, Collins PL: Parainfluenza viruses. In Fields Virology. Edited by: Knipe DM, Howley PM. Philadelphia: Lippincott Williams and Wilkins; 2006:1497-1526.
Tsurudome M, Ito M, Nishio M, Kawano M, Okamoto K, Kusagawa S, Komada H, Ito Y: Identification of regions on the fusion protein of human parainfluenza virus type 2 which are required for haemagglutinin-neuraminidase proteins to promote cell fusion. J Gen Virol 1998, 79: 279-289.
Tsurudome M, Nishio M, Ito M, Tanahashi S, Kawano M, Komada H, Ito Y: Effects of hemagglutinin-neuraminidase protein mutations on cell-cell fusion mediated by human parainfluenza type 2 virus. J Virol 2008, 82: 8283-8295. 10.1128/JVI.00460-08
Terrier O, Cartet G, Ferraris O, Morfin F, Thouvenot D, Hong SS, Lina B: Characterization of naturally occurring parainfluenza virus type 2 (hPIV-2) variants. J Clin Virol 2008, 43: 86-92. 10.1016/j.jcv.2008.05.007
Tsurudome M, Nishio M, Komada H, Bando H, Ito Y: Extensive antigenic diversity among human parainfluenza type 2 virus isolates and immunological relationships among paramyxoviruses revealed by monoclonal antibodies. Virology 1989, 171: 38-48. 10.1016/0042-6822(89)90508-4
Heikkinen T, Marttila J, Salmi AA, Ruuskanen O: Nasal swab versus nasopharyngeal aspirate for isolation of respiratory viruses. J Clin Microbiol 2002, 40: 4337-4339. 10.1128/JCM.40.11.4337-4339.2002
Precious B, Southern JA, Randall RE: Sequence analysis of the HN gene of parainfluenza virus type 2. J Gen Virol 1990, 71: 1163-1168. 10.1099/0022-1317-71-5-1163
van Wyke Coelingh KL, Winter CC, Jorgensen ED, Murphy BR: Antigenic and structural properties of the hemagglutinin-neuraminidase glycoprotein of human parainfluenza virus type 3: sequence analysis of variants selected with monoclonal antibodies which inhibit infectivity, hemagglutination, and neuraminidase activities. J Virol 1987, 61: 1473-1477.
Ng PC, Henikoff S: SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res 2003, 31: 3812-3814. 10.1093/nar/gkg509
Tanabayashi K, Compans RW: Functional interaction of paramyxovirus glycoproteins: identification of a domain in Sendai virus HN which promotes cell fusion. J Virol 1996, 70: 6112-6118.
This study was partially supported by Center of Excellence in Biotechnology Research, King Saud University through the research group project number CEBR NRL 2011/02.
The authors declare that they have no conflicts of interest.
FNA and HMA developed the concept, designed experiments and interpreted the results. MSS and HMA carried most of the experiments. FNA supervised the process of sample collection by trained nurses and helped with the experiments. All of the authors read and approved the manuscript.
About this article
Cite this article
Almajhdi, F.N., Alshaman, M.S. & Amer, H.M. Human parainfluenza virus type 2 hemagglutinin-neuramindase gene: sequence and phylogenetic analysis of the Saudi strain Riyadh 105/2009. Virol J 9, 316 (2012). https://doi.org/10.1186/1743-422X-9-316
- Phylogenetic tree
- Riyadh 105/2009
- Saudi Arabia
- Sequence analysis