Multiple alignment comparison of the non-structural genes of three strains of equine influenza viruses (H3N8) isolated in Morocco
BMC Research Notes volume 8, Article number: 471 (2015)
Three equine influenza viruses, A/equine/Nador/1/1997(H3N8), A/equine/Essaouira/2/2004(H3N8), and A/equine/Essaouira/3/2004(H3N8), were isolated from different Equidae during local respiratory disease outbreaks in Morocco in 1997 and 2004. Their non-structural (NS) genes were amplified and sequenced.
The results show high homology of NS nucleotide sequences of A/equine/Nador/1/1997 with European strains (i.e., A/equine/newmarket/2/93 and A/equine/Grobois/1/1998) and clustered into the European lineage. However, NS gene of A/equine/Essaouira/2/2004(H3N8) and A/equine/Essaouira/3/2004(H3N8) strains indicated high homology with equine influenza strains that had circulated before 1990 (A/equine/Fontainbleu/1/1979(H3N8), which belonged to a pre-divergent phase Amino acid sequence comparison of the NS1 protein with reference strain A/equine/Miami/1963(H3N8) shows that the A/equine/Nador/1/1997(H3N8) strain has 12 substitutions at the residues D/24/N, R/44/K, S/48/I, R/67/Q, A/86/V, E/139/K, A/112/T, E/186/K, L/185/F, A/223/E, S/213/T and S/228/P. In both A/equine/Essaouira/2/2004(H3N8) and A/equine/Essaouira/3/2004(H3N8) strains, the NS1 sequences present one common mutation at the residue: S/228/P.
It seems that all of these substitutions are not produced at the key residues of the RNA-binding domain (RBD) and the effector domain (ED). Consequently, we can suppose that they will not affect the potency of inhibition of cellular defences, and the virulence of the Moroccan equine strains will be maintained.
Equine influenza virus is a highly contagious agent, capable of causing explosive outbreaks of respiratory disease among susceptible horse populations in many countries . Two distinct subtypes of influenza virus, H7N7 (prototype influenza A/equine/Prague/1/56) and H3N8 (prototype influenza A/equine/Miami/1/63), have been recognized in the horse [2, 3]. The first subtype (H7N7) is believed to be extinct and has not been isolated in the horse since 1979. The other subtype (H3N8) was isolated for the first time in 1963 and was responsible for all recent equine influenza outbreaks [4, 5]. These viruses belong to the Orthomyxoviridae family type A. The genome comprises eight segments of RNA . These segments are encapsulated by a nucleoprotein (NP), which gives a helical symmetry to each. Six segments encode one protein: hemagglutinin (HA), neuraminidase (NA), nucleoprotein (NP), and three polymerases: PA, PB1, and PB2. Each of the other segments encodes two proteins: the matrix (M1 and M2) and a non-structural protein (NS1) and NS2 (NEP: nuclear export). The genomes of several virulent strains encode an 11th protein PB1-F2 . NS1, NS2, and PB1-F2 proteins participate in viral replication and are not incorporated into the viral structure.
The presence of double-stranded RNA (dsRNA) is a signal of a viral infection, causing the activation of the cellular defence system through the synthesis of the interferon protein (IFN) . However, influenza viruses have developed a strategy to antagonize the cellular defences. The non-structural protein NS1 is regarded as the major antagonist of the immune response of the host cells . In most viruses, the NS1 protein consists of 230 amino acids (aa) and the NEP protein of 121 aa . However, a size variability of the NS1 protein was observed, especially in the human and swine viruses . In the case of H3N8 equine influenza virus, truncations have been shown as well as other influenza viruses [12, 13].
The linear structure of the NS1 protein contains two domains linked by a non-structured region composed of ten amino-acids: The first one is the RNA-binding domain (RBD), representing the N-terminal protein (amino acids 1–73), the second is an effector domain (ED), which represents C-terminal protein (amino acids 74–230) .
The RNA-binding domain has an alpha-helical structure and forms a symmetrical homo-dimer with a unique six-helical chain fold. The effector domain forms a dimer, and each monomer consists of seven beta-strands and three alpha helices . The non-structural protein NS1 contains several interaction domains: the region of binding to dsRNA (residues 1–73), the specific cleavage and polyadenylation specific factor region (CPSF) (residues 175–210), the region of the binding protein (PABP)—II (poly (A)—binding protein—II) (residues 218–225), the region of nuclear localization signal 1 (NLS 1) (residues 34–38), the region nuclear localization signal 2 (NLS 2: nuclear localization signal) (residues 211–216), the region of the nuclear export signal (NES) (residues 132–141), and the region that interacts with subunit regulatory p85-β of phosphatidylinositol 3-kinase (PI3K) (residues 89–93, 137–142 and 164–167) .
The PDZ domain binding is present at the C-terminal of NS1 among equine influenza viruses (H3N8) (sequence ESEV/EPEV), and it is known as a determinant of viral virulence.
In the present study, we report analysis of the partial nucleotide sequence of the NS gene and the NS1 protein for three Moroccan equine influenza viruses: A/Equine/Nador/1/97, A/Equine/Essaouira/2/2004, and A/Equine/Essaouira/3/2004, and their genetic comparison with influenza strains (human, equine, swine, avian) available in the GenBank database.
A/equine/Nador/1/97 was isolated in Nador from a mule after four passages on eleven -day-old embryonated specific pathogen-free chicken eggs (ECE4) as described by Kissi et al. (1998) . Both A/equine/Essaouira/2/2004 and A/equine/Essaouira/3/2004 were isolated respectively from an infected donkey and a horse during 2004 outbreaks in Essaouira. The isolates were passaged on Madin Darby Canine Kidney cell line after two passages (MDCK2) at 34 °C in an atmosphere of 5 % CO2 in Eagle’s minimum essential media supplemented with 5 % foetal calf serum.
Viral RNA extraction and amplification
Viral RNA was extracted directly from isolates using a PureLink Viral RNA/DNA Mini-Kit (Invitrogen, Van Allen Way, Carlsbad, California, USA) following the manufacturer’s recommended protocol. Complementary DNA was obtained by RT reactions, which were carried out by using a SuperScript III First-Strand Synthesis System (Life Technologies, Carlsbad, CA, USA).
PCR was performed using a Platinum PCR SuperMix High-Fidelity Kit (Invitrogen, Carlsbad, CA, USA) with cDNA obtained using primers specific for NS1F (ATGGATTCCAACACTGTGTC) and NS1R (TCAAACTTCTG(A/G)CTCAATTG) at a final concentration of 0.5 µM for primers. Primer design is detailed by Tissier (2008)  and was synthesized by the Unité d’Appui Technique à la Recherche Scientifique, Centre National de Recherche Scientifique et Technique (CNRST), Rabat, Morocco. The assay was performed on the SmartCycler instrument (Cepheid, Sunnyvale, CA, USA), using the following thermocycling protocol: incubation at 95 °C for 2 min, then 35 cycles of denaturation at 95 °C for 30 s, 52 °C for 1 min for hybridization, and 72 °C for 30 s, and a final extension for 3 min at 72 °C.
Sequencing NS Genes and phylogenetic Analysis
The amplified PCR NS products were sequenced. Briefly, the PCR products were purified using EXOSAP-IT (USB Corporation, Cleveland, OH, USA) and bi-directionally sequenced by using BigDye1 Terminator v3.1 (Applied Biosystems, Foster City, Calif. USA) on a 3130xl model sequencer (Applied Biosystems). Analysis of the electro-phoregramm was carried out with the sequencing analysis software version 5.3.1 (Applied Biosystems). We performed phylogenetic analysis of 54 influenza strains (including equine, avian and human isolates) published in GenBank database, selected using the neighbour-joining method with bootstrap analyses of 500 replicates in CLUSTAL W. The tree was visualized using MEGA5.1 software (http://megasoftware.net/) .
Determination of amino acid sequences of NS genes
Multiple alignments of the deduced amino acid sequences were used by Basic Local Alignment Search Tool (BLAST) software. Comparing NS1 protein of Moroccan isolates to the reference strain (A/Equine/Miami/1/1963), the amino acid substitutions that were involved in both RNA-binding domain (RBD) and effector domain (ED) were determined.
Results and discussion
Partial NS nucleotide sequences of the Moroccan equine influenza isolates (A/equine/Nador/1/1997, A/equine/Essaouira/2/2004, and A/equine/Essaouira/3/2004) were submitted to the GenBank database, and their accession numbers of NS genes are JX182368, JQ955611, and JQ955614, respectively.
The BLAST result shows that the nucleotide sequences of the NS gene of two strains A/equine/Essaouira/2/2004 and A/equine/Essaouira/3/2004 presents high similarity between them, and towards other old strains like A/equine/Miami/1/1963 (>99 %), A/equine/Uruguay/1/63(99 %), A/eq/LaPlata/1/88 (>99 %), and A/equine/Fontainebleau/79 (98 %). These strains have circulated before 1990. Since then, the equine influenza viruses have diverged into two antigenically distinct lineages, Eurasian and American. According to several authors, the pre-divergent strains have disappeared and were supplanted by viruses that evolved in lineages and sub-lineages .
The partial nucleotide sequence of the NS gene of A/equine/Nador/1/97 strain is genetically distant towards the pre-divergent strains. Indeed, it demonstrates high homology to the strains that have circulated in Europe during the 1990s, i.e., A/equine/Newmarket/2/1993 (99 %), A/equine/Avesta/1/1993 (99 %), A/equine/Grobois/1/98 (98 %), A/equine/Italy/1199/1992 (98 %), and A/equine/Roma/5/1991 (98 %).
Data obtained from BLAST was confirmed by the phylogenetic analysis based on 54 nucleotide sequences of the NS gene of influenza A viruses, which were isolated from various species (human, avian, swine, and equine). The NS genes of the Moroccan strains are clustered into a lineage of the equine influenza virus (H3N8) that belongs to the A allele. The two strains, A/equine/Essaouira/2/2004 and A/equine/Essaouira/3/2004, belong to the pre-divergent phase, while the A/equine/Nador/1/97 belongs to the Eurasian lineage (Fig. 1). Moreover, the sequencing analysis of the HA (accession numbers JQ955607, JQ955609, and JQ955612) of A/equine/Nador/1/1997, A/equine/Essaouira/2/2004, and A/equine/Essaouira/3/2004, respectively) revealed same data. Such results show that these three strains isolated in Morocco share the same ancestors as the equine influenza (H3N8) virus and that there was no genetic reassortment between co-circulating strains in Morocco .
The data matrix containing 54 taxa and phylogenetic analyses were deposited into TreeBase under accession url: http://purl.org/phylo/treebase/phylows/study/TB2:S18066.
Comparison analysis of amino acids alignment
Compared to the reference strain A/equine/Miami/1/1963, the sequences of amino acids of the NS1 protein of the strain A/equine/Nador/1/1997 have 12 substitutions: D/24/N, R/44/K, S/48/I, R/67/Q, A/86/V, E/139/K, A/112/T, E/186/K, L/185/F, A/223/E. S/213/T, and S/228/P. These mutations occurred in both RNA-binding and effector domains (Fig. 2).
The D/24/N substitution occurred in a changing region of the RBD (19 to 28) of the NS1 protein of influenza viruses . This mutation may not affect the structure of the NS1 protein RNA binding domain or its activity. This finding is confirmed by the virulence of the A/equine/Nador/1/1997 strain, which caused severe respiratory signs among naturally infected mules. The same substitution was also reported by Parveen and Asad (2010)  concerning the human strain A/North Dakota/04/2009 (H1N1) (ACR67137).
The three mutations (R/44/K, S/48/I, and R/67/Q) revealed on the RBD do not seem to affect virulence, because they do not concern the following residues: Thr-5, Pro-31, Asp-34, Arg-35, Arg-38, Lys-41, Gly-45, Arg-46, and Thr-49, which are closely linked to the activity of the NS1 protein. Several authors have reported that these nine amino acids are the key residues of RBD activity [9, 10, 24].
Mutations that occurred on the effector domain were reported in several sites recognized as having functional roles. In fact, mutations at residues A/86/V and A/112/T occurred in a specific region of the eIF4GI protein (81 to 113), and they participate with the protein PABPI in the initiation of mRNA viral translation [25, 26].
The binding of the protein CPSF30 (cleavage and polyadenylation specific factor) inhibits the maturation of cellular pre-mRNA. This link is provided by amino acids 144, 184, and 188 (186 region) [15, 27]. It appears that both L/185/F and E/186/K mutations may not reduce the affinity of NS1 to CPSF30 protein, and consequently do not inhibit viral replication.
Three substitutions, E/139/K, S/213/T, and A/223/E, were reported in domain activity of nucleocytoplasmic distribution in the cell of the NS1 protein. The E/139/K mutation is located at the nuclear export signal region (NES: nuclear export signal) (residues 132–141) . Furthermore, similar mutations at this residue were also observed for the NS1 of the following strains: N/139/E (A/equine/North Carolina/152429/2002, ACZ47169); and E/139/A (A/equine/Xuzhou/01/2013, AHA98359) (Fig. 2). The two mutations, S/213/T and A/223/E, occurred at nuclear localization signal (NLS2) and PABPII domains, respectively. NLS2 is defined within the region between amino acids 213 and 221. This domain is present in the NS1 proteins of most influenza A virus strains. However, several substitutions affecting it were noticed, such as: S/213/T (i.e.,: A/equine/Nador/1/1997, AFN69204 and A/equine/Aboyne/1/2005, ACH95651), and P/214/S (i.e.,: A/equine/Cheshire/1/2007, ACH95645). However, nuclear localization signal 1 (NLS1) contains the stretch of basic amino acids Asp-Arg-Leu-Arg-Arg (codons 34 to 38). This sequence is conserved in three Moroccan strains as all NS1 proteins of influenza A viruses . It appears that the substitution occurred in NLS2 do not seem to interfere with the nuclear localization of NS1 protein; because the intracellular localization is regulated by one or two nuclear localization signals (NLS) and a nuclear export signal (NES) .
Multiple comparisons of the NS1 protein of 27 equine influenza viruses show that the mutations D/24/N and E/139/K characterize the A/equine/Nador/1/1997(H3N8) strain, while the A/86/V and S/213/T mutations are shared with several European lineage strains. However, the eight substitutions (R/44/K, S/48/I, R/67/Q, A/112/T, E/186/K, L/185/F, A/223/E, and S/228/P) are present in various equine influenza viruses (H3N8) belonging to different phases of evolution (pre-divergence and the American and European lineages) (Table 1).
The NS1 protein of both strains of Essaouira presents an identical sequence of them. Therefore, compared to reference strain A/equine/Miami/1/1963, the carboxyl terminus of the NS1 proteins shows a common mutations regarding substitution S/228/P. Moreover, the PDZ-binding sequence, Glu-Pro-Glu-Val (EPEV) motif appears in both avian as well as swine, human and equine viruses.
In addition, the presence of the EPEV sequence in the NS1 protein of three equine Moroccan strains is an indicator of virulence . This sequence, present at the C-terminus of NS1 of the H5N1 virus (sequence ESEV/EPEV) and the H1N1 pandemic of 1918 (sequence KSEV) but absent among seasonal human viruses, establishes binding with the PDZ domains . The ESEV sequence is common among avian viruses. Its presence in NS1 of the A/equine/Miami/1/1963(H3N8) strain is related to its ancestral avian virus, which infected horses of North America in 1963 , while the S/228/P mutation observed in the three Moroccan strains is considered to be an adaptation of avian influenza viruses (H3N8) origin to the new equine host .
The study of NS1 sequences of three Moroccan strains shows that all mutations are not produced at the key residues of RBD (RNA-binding domain) and ED (effector domain) domains, and therefore do not appear to affect the potency of inhibition of cellular factor defences, and their virulence.
basic local alignment search tool
reverse transcriptase-polymerase chain reaction
nuclear export protein
Myers C, Wilson D. Equine Influenza Virus. Clin Tech Equine Pract. 2006;5:187–96.
Sovinova O, Tumova B, Pouska F, Nemec J. Isolation of a virus causing respiratory disease in horses. Acta Virol. 1958;2:52–61.
Waddell GH, Teigland MB, Sigel MM. A new influenza virus associated with equine respiratory disease. J Am Vet Med Assoc. 1963;143:587–90.
Timoney PJ. Equine influenza. Comp. lmmun MicrobioL infect Dis. 1996;19:205–11.
Barbic L, Madic J, Turk N, Daly J. Vaccine failure caused an outbreak of equine influenza in Croatia. Vet Microbiol. 2009;133:164–71.
Lai ACK, Rogers KM, Glaser A, Tudor L, Chambers T. Alternate circulation of recent equine-2 influenza viruses (H3N8) from two distinct lineages in the United States. Virus Res. 2004;100:159–64.
Burnouf T, Griffiths E, Padilla A, Seddik S, Stephano MA, Gutierrez JM. Assessment of the viral safety of antivenoms fractionated from equine plasma. Biologicals. 2004;32:115–28.
Jiao P, Tian G, Li Y, Deng G, Jiang Y, Liu C, Liu W, Bu Z, Kawaoka Y, and Chen H. A single-amino-acid substitution in the NS1 protein changes the pathogenicity of H5N1 avian influenza viruses in mice. J Virol. 2008;82:1146–54.
Hale BG, Randall RE, Ortin G, Jackson D. The multifunctional NS1 protein of influenza A viruses. J Gen Virol. 2008;89:2359–76.
Wang L, Suarez DL, Pantin-Jackwood M, Mibayashi M, García-Sastre A, Saif YM, Lee CW. Characterization of influenza virus variants with different sizes of the Non-structural (NS) genes and their potential as a live influenza vaccine in poultry. Vaccine. 2008;26:3580–6.
Suarez DL, Perdue ML. Multiple alignment comparison of the non-structural genes of influenza A viruses. Virus Res. 1998;54:59–69.
Quinlivan M, Zamarin D, Garcı´a-Sastre A, Cullinane A, Chambers T, Palese P. Attenuation of equine influenza viruses through truncations of the NS1 protein. J Virol. 2005;79:8431–9.
Na W, Kang B, Kim HI, Hong M, Park SJ, Jeoung HY, An DJ, Moon H, Kim JK, Song D. Isolation and genetic characterization of naturally NS-truncated H3N8 equine influenza virus in South Korea. Epidemiol Infect. 2014;142:759–66.
Dundon WG, Capua I. A Closer look at the NS1 of influenza virus. Viruses. 2009;1:1057–72.
Dongzi L, Jingfang L, Zhizhen Z. Structure and function of the NS1 protein of influenza A virus. Acta Biochim Biophys Sin. 2007;39:155–62.
Darapaneni V, Prabhaker VK, Kukol A. Large-scale analysis of influenza A virus sequences reveals potential drug target sites of non-structural proteins. J Gen Virol. 2009;90:2124–33.
Kissi B, Daoudi N, El kantour A, Id sisi yahia K, Benazzou H. First isolation in the Morocco of the equine influenza virus called A/Equi/Morocco/1/97. Esp. Vét 1998; 14–10.
Tissier O. Comparative study of sequences of the equine influenza virus isolated in French. Thesis submitted to the University Claude Bernard–Lyon I (medicine–pharmacy) 2008. http://www2.vetagro-sup.fr/bib/fondoc/th_sout/phpdrupal.php?annee=2008. Accessed 11 May 2015.
Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24:1596–9.
Martella V, Elia G, Decaro N, Di Trani L, Lorusso E, Campolo M, et al. An outbreak of equine influenza virus in vaccinated horses in Italy is due to an H3N8 strain closely related to recent North American representatives of the Florida sub-lineage. Vet Microbiol. 2007;121:56–63.
Boukharta M, Touil N, El Fahim E, Terta M, Kissi B, Loutfi C, El Harrak M, Ennaji MM. Phylogenetic and molecular characterization of the equine influenza virus A (H3N8) causing the 1997 and 2004 outbreaks in Morocco. J Equine Vet Sci. 2013;34:942–8.
Suarez DL, Perdue ML. Multiple alignment comparison of the non-structural genes of influenza A viruses. Virus Res. 1998;54:59–69.
Parveen S, Asad UK. Structural and functional analysis of NS1 and NS2 proteins of H1N1 subtype. Genom Proteom Bioinform. 2010;8:190–9.
Yin C, Khan JA, Swapna GV, et al. Conserved surface features form the double-stranded RNA binding site of non-structural protein 1 (NS1) from influenza A and B viruses. J Biol Chem. 2007;282:20584–92.
Marc D. NS1 influenza viruses: an influential protein. Virologie. 2012;16:95–106.
Burgui I, Aragon T, Ortın J, Nieto A. PABP1 and eIF4GI associate with influenza virus NS1 protein in viral mRNA translation initiation complexes. J Gen Virol. 2003;84:3263–74.
Zohari S, Gyarmati P, Ejdersund A, Berglöf U, Thorén P, Ehrenberg M, Czifra G, Belák S, Waldenström J, Olsen B, Berg M. Phylogenetic analysis of the non-structural (NS) gene of influenza A viruses isolated from mallards in Northern Europe in 2005. Virol J. 2008;5:147.
Greenspan D, Palese P, Krystal M. Two nuclear location signals in the influenza virus NS1 nonstructural protein. J Virol. 1988;62:3020–6.
Tynell J, Melén K, Julkunen I. Mutations within the conserved NS1 nuclear export signal lead to inhibition of influenza A virus replication. Virol J. 2014;11:1–12.
Tombari W, Nsiri J, Larbi I, Guerin JC, Ghram A. Genetic evolution of low pathogenecity H9N2 Avian influenza viruses in Tunisia: acquisition of new mutations. Virol J. 2011;8:1–12.
Munier S, Moisy D, Marc D, Naffakh N. Interspecies transmission, adaptation to humans and pathogenicity of animal influenza viruses. Pathol Biol (Paris). 2010;58:59–68.
Dea S, Elazhary MASY, Roy RS. The human and animals influenza virus. A review of the literature. Can Vet J. 1980;21:171–8.
Bogoyavlenskiy A, Berezin V, Prilipov A, Korotetskiy I, Zaitseva I, Kydyrmanov A, Karamedin K, Ishmukhametova N, Asanova S, Sayatov M, Zhumatov K. Phylogenetic analysis of the non-structural (NS) gene of influenza A viruses Isolated in Kazakhstan in 2002–2009. Virologica sinica. 2011;26:376–85.
The work presented here was carried out in collaboration between all authors. BM participated in the design of the study, carried out the molecular genetic studies, analysis and interpretation of data, the sequence alignment phylogenetic and drafted the manuscript. SA participated equitably in the molecular genetic studies and interpretation of data. MME, EHM, are conceived of the study and helped revising the manuscript. All authors read and approved the final manuscript.
This study was supported by the University Hassan II, Mohammedia-Casablanca, Faculty of Sciences and Techniques, Mohammedia, Morocco. The authors thank the Inspection de Service de Santé Militaire/Division Vétérinaire staff for their support. Experimental and protocol designs were performed at the Laboratory of Virology, Microbiology and Quality/ETB (Faculty of Sciences and Technics, Mohammedia, University Hassan II Mohammadia-Casablanca). Strain isolations were performed at the Society of Pharmaceutical and Veterinary Products, Virology Laboratory, Rabat, Morocco. NS gene sequencing was done at Functional Genomics Platform, Technical Support Unit for Scientific Research, CNRST, and Rabat, Morocco.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Mohamed Boukharta and Souad Azlmat contributed equally to this work
About this article
Cite this article
Boukharta, M., Azlmat, S., Elharrak, M. et al. Multiple alignment comparison of the non-structural genes of three strains of equine influenza viruses (H3N8) isolated in Morocco. BMC Res Notes 8, 471 (2015). https://doi.org/10.1186/s13104-015-1441-0