Skip to main content
Download PDF

Disparate thermostability profiles and HN gene domains of field isolates of Newcastle disease virus from live bird markets and waterfowl in Uganda

Virology Journal volume 13, Article number: 103 (2016) Cite this article

Abstract

Background

Uganda poultry production is still faced with frequent outbreaks of Newcastle disease (ND) in the backyard free-range systems despite the accessibility of cross protective vaccines. Live bird markets and waterfowl has long been reported as a major source of disease spread as well as potential sources of avirulent strains that may mutate to virulent strains. ND-virus has been reported enzootic in Ugandan poultry but limited studies have been conducted to ascertain thermostability phenotypes of the Ugandan ND-virus strains and to understand how these relate to vaccine strains.

Methods

This study evaluated thermostability of 168 ND-virus field isolates recovered from live bird markets and waterfowls in Uganda compared to two live commercial vaccine strains (I2 and LaSota) by standard thermostability procedures and Hemagglutinin-Neuraminidase (HN) gene domains. The known pathotypes with thermostability profiles were compared at HN amino acid sequences.

Results

Field isolates displayed disparate heat stability and HN gene domains. Thermolabile isolates were inactivated within 15 min, while the most thermostable isolates were inactivated in 120 min. Four thermostable isolates had more than 2 log2 heamaglutinin (HA) titers during heat treatment and the infectivity of 9.8 geometric mean of log10 EID50 % in embryonated eggs. One isolate from this study exhibited a comparable thermostability and stable infectivity titers after serial passages, to that of reference commercial vaccine was recommended for immunogenicity and protection studies.

Conclusion

The occurrence of ND-virus strains in waterfowl and live bird markets with disparate thermostability and varying HN gene domains indicate circulation of different thermostable and thermolabile ND-virus pathotypes in the country.

Background

Newcastle disease (ND) is a contagious poultry infection that usually invokes trade barriers and exerting great economic threats [1]. Newcastle disease virus (ND-virus) genome contains 6 genes encoding for structural proteins positioned from 3′ to 5′: Nucleoprotein (NP), Phosphoprotein (P), Matrix (M), Fusion (F) Hemagglutinin-Neuraminidase (HN) and polymerase (L) [23]. The physico-chemical properties of these proteins have been extensively used to categorize strains. According to World Organization for Animal Health [27], for a ND-virus to be notified as pathogenic, it has to meet the criteria of Intracerebral pathogenicity Index (ICPI) of beyond 0.7. Today, field strains are categorized by virulence into five pathotypes [3, 5, 6], confirmed by multiple basic amino acid sequences in the fusion (F) protein cleavage site (FPCS) [12, 26].

Great importance is attached to the thermostability of live vaccinates because of their retained potency in areas with poor transport, storage facilities, human error and power shortages. Studies in search of new thermostable isolates assess virulence of pathotypes by determining the amino acid sequences at the F protein cleavage site (FPCS), mean death time (MDT) and Intracerebral pathogenicity indices (ICPI). In areas where there is frequent use of vaccines, selective pressure leading to new strains or failed vaccination has been recorded due to the use of phylogenetically divergent vaccine strains from the circulating local ND-viruses [24, 28]. Besides, ND-virus isolates with amino acids at FPCS typical of low-virulent type and ICPI values typical of highly virulent strains have been reported [37]. Such places are experiencing continual outbreaks of velogenic ND in backyard and vaccinated flocks [1, 4].

HN gene is known to play crucial multifunctions in infectivity [14, 42]. It has also been confirmed through reverse genetics using cDNA clones (ICS) as a crucial determinant of NDV thermostability with conferred complete protection of birds [43]. Thermostability characterization of heat-stable ND-viruses relating to HN gene will improve understanding of the molecular basis of ND-virus thermostability. This will augment rational design of NDV vaccines to solve managerial problems in scavenging rural chicken [18] or further evaluate performance of established vaccine strains.

Recovery of virulent ND strains in birds vaccinated with avirulent vaccines strains or in ND-virus endemic countries raises concerns whether these strains are derived from indigenous or vaccine strains [2]. In Uganda, available ND vaccines used are namely: LaSota strain and I2 strain with well-defined genotypes. However, the recent ND-virus isolates from Ugandan live bird markets indicate a high pathogenicity and low evolution [8] with no thermostability phenotypes. Several studies have recovered low-virulent, thermostable ND-virus in several species of birds [5, 33]. Such strains have been evaluated as vaccine candidates to protect village birds in the tropics where temperatures are high and local farmers either lack or are unable to pay for the cold-chain needed to sustain live thermolabile ND-virus vaccine usage [35].

Thermostability of hemagglutinins of ND-virus is retained through many serial passages in embryonating eggs and can distinguish one culture from another [13]. The basis of thermostability testing of ND-virus is that all strains have hemagglutinin surface proteins, which agglutinates chicken RBCs in vitro, and Neuraminidase enzyme that promotes virus release from infected cells. The activity of these surface glyco-proteins is used to detect hemagglutinating viruses in the family [27], and to follow their stability when the virus is exposed at different temperatures [1]. Since the development of a thermostable Australian V4 ND vaccine strain, several thermostability testing methods have been used to evaluate live ND-virus isolates and vaccine candidates for ability to survive under different temperatures [19, 21, 22, 36].

The aim of the current study was to determine the thermostability of ND-virus isolates recovered from live bird markets (LBM) and waterfowls and relate this phenotype to HN gene sequence.

Materials and methods

ND-virus isolates

A total of 168 ND-virus isolates were used in the current study. These field ND-viruses were isolated in Uganda during 2011 from LBM comprised of local chickens, turkeys and ducks. The isolates from waterfowl were recovered from freshly voided fecal droppings of migratory and resident waterfowls in selected landing sites. Sixty seven percent (112/168) of the viruses of LBM were isolated from cloacal (C) and/oro-pharyngeal (P) swabs, while thirty three percent (56/168) of isolates from waterfowls were recovered from fresh fecal (F) matter. Initially, viruses were passaged twice in 10–day-old embryonated chicken eggs (ECE) obtained from ND-free flocks and clean isolates obtained in chicken embryonated fibroblasts (CEF). The allantoic fluid was aseptically harvested, virus presence confirmed by HI test according to standard procedures and analyzed further or stored at −20 °C as virus stocks. Biological characterization for chicken ND isolates were described elsewhere [8]. Two live commercial vaccines (I2 and LaSota) were used as controls throughout the study.

The presence of the ND-virus was confirmed by PCR and serologically with positive sera in HI test [8, 27]; using 0.5 % (v/v) washed ND-virus–free chicken RBCs according to standard procedures. Serum was generated in in-house rabbits immunized with ND-virus - LaSota vaccine strain (Laborotorios Hipra, S.A, Spain).

Assessment of thermostability

HA and infectivity titers for all the stock viruses were determined by standard microtiter plate hemagglutionation assay [27] and the viability of selected isolates by end point infectivity assay [30]. Nine paired sealed vials containing 0.5 mL aliquots of each ND-virus isolate were thawed once. One pair was left on ice while the other eight were incubated on a water bath (Polypro Bath® CA, USA) kept at a constant temperature (56 ± 0.5 °C). At regular interval of 15, 30, 45, 60, 75, 90, 105 and 120 min, a pair of vials was removed and chilled quickly on ice to stop heat inactivation. All the aliquots were assayed for HA activity by standard methods [27, 30]. The isolates that showed HA titers at 56 °C after 2 h were further analyzed for regression analysis of HA activity using rate constants (k) and among these, avirulent isolates were passaged five times without heating intervals. Additionally, to determine the effect of heat on the viability or multiplication of the selected viruses, the vials incubated at 56 °C for 0, 30, and 60 min were assayed for virus viability and compared to the thermostable I2 strain. Briefly, 0.3 mL of the virus aliquot from each mentioned treatment was injected into three sets of eggs and incubated for 4 days before allantoic fluid was harvested, checked for the mean log2 HA titer and assayed for infectivity.

Stability of infectivity

The infectivity stability of the selected avirulent viruses was evaluated and compared to vaccine strains (I2 and LaSota) in 10-day-old embryonated eggs. The infectivity evaluation of the isolates involved five passages. Briefly, tenfold serial dilutions were made in normal saline solution (Phosphate Buffered Saline, PBS) and five dilutions were selected. One hundred microliters (0.1 mL) of each of these dilutions: 10-1, 10-6, 10-7, 10-8 and 10-9 was inoculated in duplicate into five eggs, incubated for 4 days, and harvested for HA and median embryo infectivity dose at fifty percent (log10EID50) assays. This was carried out for five passages where both HA and infectivity assays were carried out for each passage. Additionally, the heamagglutionation units were measured following passages, to evaluate the replication competence of the isolates.

Viral RNA extraction, RT-PCR and DNA sequencing

Viral RNA was extracted from all the HI positive allantoic fluid samples by using the QIAamp Viral RNA mini kit (Qiagen, USA) according to the manufacturer’s instructions and sequenced to obtain full HN sequences as described in [8].

HN gene sequence analysis

To predict the HN relationship between the field ND-virus isolates and vaccine strains, multiple sequence alignment of HN gene of representative thermostable ND-virus isolates was done. Using the following accession numbers previously deposited in the Gene Bank for isolates; NDV/UG/MU/007- HG937536, NDV/UG/MU/010-HG937538, NDV/UG/MU/022- HG937542, NDV/UG/MU/039- HG937548, NDV/UG/MU/059-HG937553, NDV/UG/MU/098- HG937564 and NDV/UG/MU/111- HG937566. Together with other representative thermostable vaccines: TS09C- JX110635, I2 – AY935499, I2 progenitor – AY935500, NDV4-C - JX443519 and thermolabile: Lasota –JF950510, LaSota C5 – KC844235, and D58 strain- EU305607 and others were blast searched from NCBI GeneBank using the Bio-Edit software (North Carolina State University, USA) to generate sequence analysis data. The parameters used for the sequences analysis were multiple alignment (Clustal W), sequence identity plotter and sequence matrix at both amino acid and nucleotide levels. The B- cell epitope of HN gene was predicted in a computer algorithm using web server based on software www.bioinfo.tsinghua.edu.cn/epitope/EPMLR. The transmembrane domains were predicted using Dense Alignment Surface (DAS) transmembrane domain prediction server– http:www.sbc.su.se/-miklos/DAS/maindashtml [11].

Statistical analyzes

Thermostability was analyzed by calculating the slope of the regression curve of HA activity. When appropriate, the data of exponential decline in activity of HA were presented as rate constants (k, expressed as mins-1), defined by first-order inactivation equation below as described [41].

$$ 2.303\ log10\ \left(V/Vo\right)=-kt $$

Isolates whose late time points, i.e., 105 or 120 min did not contain HA were excluded from analysis. The survival analysis for the isolates were titrated by standard median EID50 % according to standard protocols [27]. All statistical analyzes were performed in GraphPad Prism, version 5.

Results

Thermostability of field ND-virus isolates at 56 °C

The 168 field isolates showed HA titers greater than 16 log2 in embryonated eggs during the initial culture and were confirmed as ND-viruses by HI test and PCR. Regarding the source of poultry isolates, Northern region had the highest number of isolates compared to other regions. ND-viruses occurring in different species of poultry existing in the four regions of Uganda (Fig. 1). All the isolates used showed disparate thermostability at 56 °C ranging from 15 to 120 min (Table 1). Of these study isolates, 13.7 % (23/168) retained more than 2 log2 of their initial HA titers after 120 min and were regarded as thermostable while 86.3 % (145/168) without residual HA titer after the above time were considered thermolabile isolates. Of the 23 thermostable isolates, 13 were from LBM while 10 were from waterfowls (Table 1). Based on the heat inactivation rate constant (k) of HA activity, the thermostable strains were included in Table 2. When compared by rate constant of HA activity, 26.1 % (6/23) of the 23 thermostable ND-virus isolates were inactivated at a lower rate (higher heamagglutinin activity) than the reference I2 vaccine strain. Also all these isolates were from waterfowls.

Fig. 1
figure 1

Distribution of ND-virus isolates by region and species of birds. The numbers above bars indicate isolates in different regions and bird species

Full size image
Table 1 Thermostability profiles with numbers of NDV isolates by HA activity determined at 56 °C incubation
Full size table
Table 2 Hemagglutination titers and heat inactivation rate constant (k) of thermostable ND-viruses incubated at 56 °C
Full size table

By plotting the logarithmically transformed (log10) HA of the selected avirulent isolates against time, the inactivation rates were further compared (Fig. 2). Except for NDV133/UG/MU/2011 and LaSota that yielded a monophasic curve (Fig. 2a and d), other selected thermostable isolates produced biphasic curves (Fig. 2b-d). For biphasic curves, from zero to 90 time points, a plateau that was observed was followed by a linear decline. Overall, isolate NDV133/UG/MU/2011 had better survival rate compared to NDV173/UG/MU/2011, NDV177/UG/MU/2011 and NDV178/UG/MU/2011 as avirulent thermostable isolates.

Fig. 2
figure 2

In vitro HA stability test of indicated thermostable field ND-virus at different time intervals. Incubation at 56 °C revealed loss of titer for isolates after 120 min (a–d). Thermostable I2 (e) and thermolabile LaSota (f) reference strains were included as controls. Values were averages of two independent tests (Mean ± SD, n = 2)

Full size image

Stability of infectivity upon passage

Four clean thermostable isolates were passaged five times and their infectivity shown in Fig. 3. Infectivity of the isolates decreased along the five serial passages. At the end of 3rd passage, the infectivity of isolate NDV177/UG/MU/2011 rapidly declined to a geometric median log10 EID50 %/mL titer of 1.25, while that of NDV178/UG/MU/2011 was at 4.55 at the same passage. The infectivity titers of NDV133/UG/MU/2011 and NDV173/UG/MU/2011 isolates remained fairly close to that of I2 vaccine strain at geometric median log10 EID50 % of 7.4. A greater change in infectivity of the isolates was at the fourth and fifth passages. Comparatively, one isolate NDV133/UG/MU/2011 retained more than 6.5 log10EID50 %/mL of its infectivity titer compared to I2 and LaSota of 8.05 log10EID50 %/mL and 7.5 log10EID50 %/mL respectively at the 5th passage. The infectivity titers of isolates following a less than one hour time heat exposure were compared in Table 3. The ND-virus isolates NDV133/UG/MU/2011, NDV173/UG/MU/2011 and ND178/UG/MU/2011 retained more than 50 % of their log10EID50 % infectivity within one hour of exposure at 56 °C. LaSota (thermolabile) vaccine strain had no detectable infectivity after 30 min of heat exposure while I2 (thermostable) strain retained its infectivity.

Fig. 3
figure 3

Stability of infectivity of indicated ND-virus strains by geometric mean log10EID50% following serial passages. A greater change in the titers of field strains occurred at 4th and 5th passages compared to I2 reference vaccine. LaSota, a thermolabile vaccine was unheated initially but included to compare its viability with that of thermostable isolates. Values indicated are means of three independent tests (Mean ± SD, n = 3)

Full size image
Table 3 HA and infectivity titers of avirulent, thermostable ND-viruses at 56 °C for 0, 30, 60 min
Full size table

HN gene sequences

The HN gene of Ugandan isolates contains coding sequence (cds) comprising of 1713 nucleotides coding for 571 amino acids with the stalk and globular head major areas comprising of 1-143 amino acids and 125-571 amino acids respectively. The transmembrane domains of the selected thermostable were variable with three to five predictable domains at positions 24–47, 25–45, 96–97, 101–107 and 557–563. 13 cysteine (C) residues were conserved for all our isolates at position 123, 172, 186, 196, 238, 247, 251, 344, 455, 461, 465, 531 and 542. Six glycosylation sites were also conserved across all isolates at positions G1-119, G2-341, G3-433, G4-481, G5-508 and G6-538. The details is provided in Table 4. The sequence identity matrix for HN gene of the selected thermostable isolates and other strains of ND-virus are provided in Table 5 and the predicted B-cell epitopes for HN genes are given in Table 6

Table 4 Predicted transmembrane amino acid sequences analysis for HN gene of thermostable ND-viruses and vaccine strains
Full size table
Table 5 Sequence identity matrix for CDS of HN gene of thermostable ND-viruses and vaccine strains
Full size table
Table 6 B-cell epitopes of HN gene of representative thermostable field ND-viruses
Full size table

Discussions

Thermostability and infectivity of isolates

Several thermostable strains like I2, HR-V4 and V4 have been isolated and characterized [7, 16, 34]. This specific study profiled the thermostability of the Ugandan ND-virus isolates and compared them at HN protein to differentiate local strains from existing vaccines strains. While thermostability is considered to be a natural property of ND-virus isolates used to rapidly characterize field ND-virus isolates in epizootiological studies in the absence of any other marker [13], Ugandan field ND-virus isolates showed disparate thermostability profiles which is not related to the virulence of isolates, site of isolation in the bird species or geographical locality of isolation. This confirms the findings from other studies [20, 33]. Although poultry samples were collected from places of different climatic conditions of their natural sources, thermostability profiles of isolates did not reflect these conditions. The average thermostability of isolates from climatically hotter North and North-East regions of Uganda was 45 min, while that from cooler central and western regions was 60 min. This observation is in agreement with other studies supporting the fact that heat resistance in ND-virus is not a climatically determined trait [15].

The natural spread of clonal derivatives of vaccines strains in a locality following massive vaccination has been observed [3]. In the present study, thermostable isolates were recovered from unvaccinated flocks. This suggests some degree of natural area spread which explains disparate thermostability profiles of both waterfowl and live bird market bird isolates. This is further confirmed by presence of new genotypes our previous study of genotypic characterization of LBM isolates [8]. Taken together, there is a spread of field viruses with varied pathogenicity, thermostability and antigenicity in tropical developing countries, which require a carefully designed control measures.

A prolonged exposure of ND-virus isolates to a high temperature of 56 °C has been preferred to short time exposure at lower temperature during heat selection [17, 41]. This is to ensure selection of a more thermostable ND-virus from a heterogenous population containing a mixture of heat-stable and heat labile strains [19, 22]. We recovered 13.7 % (23/168) of field isolates (LBM and waterfowl) with residual HA titers after 2 h exposure at 56 °C. Of the 23 isolates, 26.2 % (6 out of 23 isolates) had greater heamagglutinin titers greater than reference I2 strain. These results are in conformity with other studies especially [19], who reported 38 % of the field ND-virus isolates with a greater heamagglutinin than reference thermostable strains. By regression analysis, we noted thermostable and thermolabile phenotypes of strains that was not influenced by site of isolation. Isolates from oro-pharyngeal (P) and cloaca (C) of the same bird had varied thermostability profiles. This confirms the idea that thermostability is not truly a genotypic character but just a trait that can readily be increased by selection involving heat shock [35].

Like in the Australian experience of enhancing thermostability of V4 vaccine strain from a subpopulation of known heat-stable virus also adopted by [16], we subjected clean ND-virus isolates to the same temperature (56 °C), followed HA titers until 2 h and selected four most thermostable isolates. Previous studies have attempted to isolate avirulent, thermostable ND-viruses from the feces of wild migratory birds for vaccine candidates [32, 33, 40]. Here, the present study identified thermostable ND-viruses from waterfowls and this could provide a promising line of research in thermostable vaccine development in Uganda and providing a base to develop an affordable thermostable ND vaccine.

DNA sequence analysis of the HN gene

HN is a type II homotetrameric glycoprotein with a monomer length of 577 amino acids for most NDV strains [25]. The heamagglutinin, Neuraminidase and thermostability are one of the multifunctional functions of HN protein that plays key roles in the steps of the NDV life cycle [14]. Mutations within the HN gene has been reported to contribute to thermostability phenotype and immunogenicity of ND-virus [39, 43]. Deletions of amino acid, R 403 were thought to influence thermostability but attempts by several studies to find alterations at such a position was futile in thermostable strains. However, evaluation of chimeric ND-viruses confirmed that thermostability phenotype of ND-virus was dependent upon the origin of HN segment and not any other ND-virus genome [43]. Still no study has offered a clear mechanism or sequence based analysis of HN gene for thermostability. By analyzing the sequences of field thermostable and thermolabile viruses, we noted no deletions at the suggested site but observed increased proportion of charged amino acid residues at the expense of uncharged polar amino acid among the thermostable isolates. This could confer rigidity and stability by minimizing deamination and backbone cleavage.

The change to a higher proportion of isoleucine (I), leucine (L), valine (V) and arginine (R) at the stalk and globular regions of NDV HN ectodomain of thermostable isolates offers an explanation to the thermostability since these aliphatic amino acid residues contribute to hydrophobic interaction, a force that maintains the conformational internal protein integrity. Besides, it has been observed that substitution of amino acids at the conserved stalk spike especially the leucine zipper motif affect the Neuraminidase activity of the globular domain where HN molecule interacts with sialic acid binding and F protein to initiate fusion. Protein stabilities have been predicted by sequence feature-based predictions upon amino acid substitutions using models [38]. Our observed amino acids substitutions at V24I, I43S, V45L, T48A, R62A, T73L, V266A, R269S, G293K, S310D, S315P, I404V, I477V and N440S in our thermostable isolates all decreased protein stability with prediction confidence of >80 % which in our case partly explains thermostability phenotype shown in Fig. 4.

Fig. 4
figure 4

Multiple amino acid sequence alignment of HN proteins. Alignment was performed by dividing amino acids into stabilizing and destabilizing groups. Black bold fonts marked by asterik (*) denote point variation in identity between thermostable and thermolabile viruses

Full size image

Several amino acid residues have been reported as functional amino acid residues, F220, S222, L224, E401, R416, and Y526 in HN protein. 13 cysteine amino acid residues at position 123, 172, 186, 196, 238, 247, 251, 344, 455, 461, 465, 531 and 542. Six potential glycosylation sites at G1-119 G2 -341, G3-433, G4-481, G5508 and G6-538, together with the different length of HN cds [29, 31]. The salient features of HN protein of our strains include: 571 amino acid size; conserved amino acids: R 174, I 175, R 197, D 198, K 236, R 416, R 498, Y 526 and E 547; potential glycosylation sites being G1-119, G2-341, G3-433, G4-481, G5-508 and G6-538, which is similar to genotype II vaccine strains Beaudette C, B1 and Ulster. Replacement of glycosylation site G5-508 in other ND-virus by tryptophan and not serine (G5-N508Y) has been used to pathotype ND-virus as avirulent [29]. All our selected thermostable isolates had this glycosylation site confirming their high pathogenicity reported earlier [8]. Further, a conserved amino acid E347 present in all thermostable and thermolabile isolates, a feature common to all vaccine strains of genotype II ND- viruses was observed. The percent similarity in relation to amino acids between isolates and thermostable I2 and TS09C is 80.6 % and that of thermolabile LaSota is 80.5 % suggesting our isolates share some attributes of these vaccine strains shown in Table 5. It could not be establish why isolates obtained from unvaccinated birds share this similarity more so with I2 vaccine that has not been widely used in the country at the time of isolation.

The most variable portion of NDV HN protein is present in the N-terminal 78 amino acids that included the transmembrane domain at position 24–47, 25–45, 557–563 using online bioinformatics webserver DAS-transmembrane prediction server [11]. The significant cut off value for the transmembrane domain considered was 2.2 at 25–45 length common to all isolates and vaccinates. However, the variable predicted number of transmembrane domains for thermostable isolates were similar to other thermostable vaccine isolates I2, V4 more than those of thermolabile ones. The sialic acid binding sites were all conserved across isolates irrespective of their thermostability phenotypes. This finding is predictable since this site plays a key attachment role of the virus to cells. The presence of cysteine amino acid residue at position 123 has been reported to be essential for intramolecular disulphide bonds that stabilize the oligomeric HN structure gene [10]. Our indicated isolates had conserved cysteine amino acid at position 123 implying absence of chymotrypsin cleavage site and these could form disulphide-linked dimers increasing the hydrophobic properties of the entire HN molecule responsible for thermostability after exposure to 56 °C. The R416 essential for receptor binding, and Neuraminidase (NA) and haemagglutination (HA) activity was likewise conserved even after exposure to 56o C, as were the receptor binding site involving E401, R416 and Y526 [9]. The predicted amino acids present in the transmembrane domain of our isolates 25FRIAVLSLIIMILVISVSILVY46 which is different from less virulent thermostable I2 or TS09C (Table 6). This variation could be due the virulence difference of these two categories.

The B-cell epitopes of thermostable isolates were 26 except for one isolate NDV/UG/MU/059 having 25 epitopes with four and 24 epitopes in the stalk and globular regions respectively (Table 6). This finding is consistent with earlier finding of the immunodorminant epitope concentrated in the globular region [10]. Earlier studies have confirmed that immune system responds to HN protein rather than the F gene and confirmed by use of neutralizing monoclonal antibodies with overlapping site at the HN gene. This variations of our isolates by HN gene suggests the existence of thermostable strains different from vaccine strains as already confirmed by our earlier findings using the HN phylogenetic tree analysis constructed by neighbor joining (NJ) algorithm with bootsrap along with distance.

Conclusion

From the current study, it can be concluded that field ND-virus (virulent or avirulent) circulating in LBM and waterfowl in Uganda display disparate thermostability profiles. Waterfowl provides a better source of avirulent, thermostable isolates with more than 4 logarithmic orders of HA within one hour of incubation at 56 °C and without loss of two logarithmic orders of infectivity titer within 30 min at serial passages. Chimeric NDV molecules has yielded new knowledge towards the understanding of NDV thermostability these require obtaining metastable HN molecules and performing mechanistic studies on how thermostability phenotype can be predicted by analysis of HN gene and such changes in amino acids of HN gene might contribute to the stability of the ND-virus strain at high temperature.

Accession numbers

The complete HN gene sequences of isolates analyzed in this study are deposited in GenBank with accession numbers; Isolates [HG937536, HG937538, HG937542, HG937548, HG937553, HG937564 and HG937566], vaccinates [JX110635, AY935499, AY935500, JX443519, JF950510, KC844235 and EU305607].

Abbreviations

EID50 % :

Median Embryo Infectivity Dose

HA:

Hemagglutination test.

ND-virus:

Newcastle disease virus

PBS:

Phosphate buffered saline

RBCs:

Red blood cells

References

  1. Alders R, Costa R, Dias P, Fringe R, Fumo A, Lobo Q, Mata B, Silva A, Young M. Investigations into the Control of Newcastle Disease in Village Chickens in Mozambique: Review of work done in Mozambique. In: Report on the ACIAR/INIVE Newcastle Disease Control Project Coordination Meeting, Tofo, Inhambane Province, Appendix. 2001.

    Google Scholar 

  2. Alexander D. Newcastle disease manual of standards for diagnostic tests and vaccines. Paris: Office International des Epizooties; 1996.

    Google Scholar 

  3. Alexander DJ. Newcastle disease and avian paramyxovirus infections. In: Calnek BW, Barnes HJ, Beard CW, McDougald LR, editors. Diseases of poultry. 10th ed. Ames, Iowa.: Iowa State University Press; 1997. p. 541–69.

    Google Scholar 

  4. Alexander DJ. Newcastle disease virus and other avian paramyxoviruses. In: Swayne DE, Glisson JR, Jackwood MW, Pearson JE, Reed WM, editors. A Laboratory Manual for the Isolation and Identification of Avian Pathogens. 5th ed. Athens, GA: American Association of Avian Pathologists; 2008. p. 156–63.

    Google Scholar 

  5. Alexander D, Mackenzie J, Russell P. Two types of Newcastle disease viruses isolated from feral birds in Western Australia detected by monoclonal antibodies. Aust Vet J. 1986;63:365–7.

    Article  CAS  PubMed  Google Scholar 

  6. Alexander D, Russell P, Parsons G, Elzein EA, Ballouh A, Cernik K, Engstrom B, Fevereiro, Fleury, Guittet M. Antigenic and biological characterisation of avian paramyxovirus type I isolates from pigeons‐an international collaborative study. Avian Pathol. 1985;14:365–76.

  7. Bensink Z, Spradbrow P. Newcastle disease virus strain I 2–a prospective thermostable vaccine for use in developing countries. Vet Microbiol. 1999;68:131–9.

    Article  CAS  PubMed  Google Scholar 

  8. Byarugaba DK, Mugimba KK, Omony JB, Okitwi M, Wanyana A, Otim MO, Kirunda H, Nakavuma JL, Teillaud A, Paul MC. High pathogenicity and low genetic evolution of avian paramyxovirus type I (Newcastle disease virus) isolated from live bird markets in Uganda. Virol J. 2014;11:1.

    Article  Google Scholar 

  9. Connaris H, Takimoto T, Russell R, Crennell S, Moustafa I, Portner A, Taylor G. Probing the sialic acid binding site of the hemagglutinin-neuraminidase of Newcastle disease virus: identification of key amino acids involved in cell binding, catalysis, and fusion. J Virol. 2002;76:1816–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Crennell S, Takimoto T, Portner A, Taylor G. Crystal structure of the multifunctional paramyxovirus hemagglutinin-neuraminidase. Nat Struct Mol Biol. 2000;7:1068–74.

    Article  CAS  Google Scholar 

  11. Cserzö M, Wallin E, Simon I, von Heijne G, Elofsson A. Prediction of transmembrane alpha-helices in prokaryotic membrane proteins: the dense alignment surface method. Protein Eng. 1997;10:673–6.

    Article  PubMed  Google Scholar 

  12. Glickman RL, Syddall RJ, Iorio RM, Sheehan JP, Bratt MA. Quantitative basic residue requirements in the cleavage-activation site of the fusion glycoprotein as a determinant of virulence for Newcastle disease virus. J Virol. 1988;62:354–6.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Hanson R, Spalatin J: Thermostability of the hemagglutinin of Newcastle disease virus as a strain marker in epizootiologic studies. Avian diseases 1978;22:659–65.

  14. Huang Z, Elankumaran S, Panda A, Samal S. Recombinant Newcastle disease virus as a vaccine vector. Poult Sci. 2003;82:899–906.

    Article  CAS  PubMed  Google Scholar 

  15. Ibu O, Okoye J, Baba S, Soyinka S, Chah K, Antiabong J, Eze D, Salihu E, Oladele S. Thermostability Profile of Newcastle Disease Viruses Isolated from Wild Birds in Central Nigeria and the Selection of a Thermostable Clone from the Sub-Population. Int J Poult Sci. 2010;9:791–4.

    Article  Google Scholar 

  16. Ideris A, Ibrahim AL, Spradbrow P. Vaccination of chickens against Newcastle disease with a food pellet vaccine. Avian Pathol. 1990;19:371–84.

    Article  CAS  PubMed  Google Scholar 

  17. Jae-Hong K, Haan-Woo S, Il-Hwan K, Eun-Kyoung L, Kang-Seuk C, Daniel Jack K. Characterization and comparison of the pathogenicity of viscerotropic velogenic Newcastle disease virus isolates in Korea. Korean J Vet Res. 2012;52:213–21.

    Google Scholar 

  18. Kinde H, Utterback W, Takeshita K, McFarland M. Survival of exotic Newcastle disease virus in commercial poultry environment following removal of infected chickens. Avian Dis. 2004;48:669–74.

    Article  PubMed  Google Scholar 

  19. King D. Selection of thermostable Newcastle disease virus progeny from reference and vaccine strains. Avian Dis. 2001;45:512–6.

    Article  CAS  PubMed  Google Scholar 

  20. King DJ, Seal BS: Biological and molecular characterization of Newcastle disease virus (NDV) field isolates with comparisons to reference NDV strains. Avian diseases 1998;42: 507–16.

  21. Lomniczi B. Thermostability of Newcastle disease virus strains of different virulence. Arch Virol. 1975;47:249–55.

    Article  CAS  PubMed  Google Scholar 

  22. Martin PA. The epidemiology of Newcastle disease virus in Village Chicken. In: Spradbrow PB, editor. Newcastle disease in village chickens, Control with Thermostable oral vaccines proceedings; International workshop held in Kwala Lumpar, Malaysia 1991. Canbera: ACIAR; 1992. p. 40–5.

  23. Mayo M. Virus taxonomy-Houston 2002. Arch Virol. 2002;147:1071–6.

    Article  CAS  PubMed  Google Scholar 

  24. Miller PJ, Decanini EL, Afonso CL. Newcastle disease: evolution of genotypes and the related diagnostic challenges. Infect Genet Evol. 2010;10:26–35.

    Article  CAS  PubMed  Google Scholar 

  25. Mirza AM, Deng R, Iorio RM. Site-directed mutagenesis of a conserved hexapeptide in the paramyxovirus hemagglutinin ± neuraminidase glycoprotein : effects on antigenic structure and function. J Virol. 1994;68:5093–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Nagai Y, Klenk H-D, Rott R. Proteolytic cleavage of the viral glycoproteins and its significance for the virulence of Newcastle disease virus. Virology. 1976;72:494–508.

    Article  CAS  PubMed  Google Scholar 

  27. OIE. Newcastle disease: Manual of standard for diagnostic tests and vaccines. Paris: Office international des Epizootics; 2012. p. 221–3.

    Google Scholar 

  28. Perozo F, Marcano R, Afonso CL. Biological and phylogenetic characterization of a genotype VII Newcastle disease virus from Venezuela: efficacy of field vaccination. J Clin Microbiol. 2012;50:1204–8.

    Article  PubMed  PubMed Central  Google Scholar 

  29. Ponnusamy P, Kirubaharan JJ, Chandramohan A. Sequence analysis of HN and F genes of a less virulent Newcastle disease virus isolated from unvaccinated village chicken. Int J Poult Sci. 2009;8:985–94.

    Article  CAS  Google Scholar 

  30. Reed LJ, Muench H. A simple method of estimating fifty per cent endpoints. Am J Epidemiol. 1938;27:493–7.

    Google Scholar 

  31. Römer-Oberdörfer A, Werner O, Veits J, Mebatsion T, Mettenleiter TC. Contribution of the length of the HN protein and the sequence of the F protein cleavage site to Newcastle disease virus pathogenicity. J Gen Virol. 2003;84:3121–9.

    Article  PubMed  Google Scholar 

  32. Rosenberger J, Klopp S, Krauss W. Characterization of Newcastle disease viruses isolated from migratory waterfowl in the Atlantic flyway. Avian diseases 1975;19:142–49.

  33. Shim J-B, So H-H, Won H-K, Mo I-P. Characterization of avian paramyxovirus type 1 from migratory wild birds in chickens. Avian Pathol. 2011;40:565–72.

    Article  PubMed  Google Scholar 

  34. Simmons G. The isolation of Newcastle disease virus in Queensland. Aust Vet J. 1967;43:29–30.

    Article  CAS  PubMed  Google Scholar 

  35. Spradbrow P. Newcastle disease in village chickens. Poult Sci Rev. 1993;5:57–96.

    Google Scholar 

  36. Spradbrow P, Bensink Z, Grimes S. Small-scale production and testing of Newcastle disease vaccine. In: Manual produced for an ACIAR-funded workshop held at the Poultry Reference Laboratory, Onderstepoort. 1995. p. 4–9.

    Google Scholar 

  37. Tan L-T, Xu H-Y, Wang Y-L, Qin Z-M, Sun L, Liu W-J, Cui Z-z et. Molecular characterization of three new virulent Newcastle disease virus variants isolated in China. J Clin Microbiol. 2008;46:750–3.

  38. Teng S, Srivastava AK, Wang L. Sequence feature-based prediction of protein stability changes upon amino acid substitutions. BMC Genomics. 2010;11:1.

    Article  Google Scholar 

  39. Uthrakumar A, Vijayarani K, Kumanan K, Thiagarajan V. Molecular Characterization of Thermostable Newcastle disease virus Isolated from Pigeon. Indian J Virol. 2013;24:42–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Vickers M, Hanson R: Characterization of isolates of Newcastle disease virus from migratory birds and turkeys. Avian diseases 1982;26:127–33.

  41. Wambura P, Meers J, Spradbrow P. Thermostability profile of Newcastle disease virus (strain I–2) following serial passages without heat selection. Tropl Anim Health Prod. 2006;38:527–31.

    Article  CAS  Google Scholar 

  42. Wehmann E, Herczeg J, Tanyi J, Nagy É, Lomniczi B. Lentogenic field isolates of Newcastle disease virus isolated in Canada and Hungary are identical with the vaccine type used in the region. Avian Pathol. 1999;28:6–12.

    Article  CAS  PubMed  Google Scholar 

  43. Wen G, Hu X, Zhao K, Wang H, Zhang Z, Zhang T, Yang J, Luo Q, Zhang R, Pan Z. Molecular basis for the thermostability of Newcastle disease virus. Sci Rep. 2016;6:22492.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The World Bank Millennium Science Initiative Project (Grant #MSI/03/32/2010) to DKB through the Uganda National Council of Science and Technology supported collection of samples, isolation and sequencing of the isolates. We recognize the Department of Animal Health, INRA, France and the French Ministry of Agriculture for sequencing isolates by the “IDEX-nouveaux entrants” program. We are grateful to the members of Emerging Infectious Disease Program (EIDP) laboratory for housing the laboratory work. We thank Prof. P.N. Nyaga for reading through the manuscript and providing useful comments.

Funding

Collection of samples, analysis of data and part sequencing of isolates was funded by the World Bank Millennium Science Initiative (MSI) Project (Grant # MSI/03/32/2010) to DKB through the Uganda National Council of Science and Technology.

Availability of data and materials

This is not applicable since this is part of a large project pending settling the intellectual property issues.

Authors’ contributions

JBO performed thermostability studies, analyzed data and sequences and drafted the manuscript; JBO, AW and KKM performed the virus isolation in eggs. DKB, JBO, MOO, JLN and HK designed this study; DKB, MOO planned and coordinated the study. All authors’ have read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

This study was part of a bigger study approved by Uganda National Council for Science and Technology under reference number HS776. This study was reviewed and approved by the College of Veterinary Medicine research and Higher Degrees Committee with regard to scientific content and ethics.

Author information

Authors and Affiliations

  1. College of Veterinary Medicine, Makerere University, P.O. Box 7062, Kampala, Uganda

    John Bosco Omony, Agnes Wanyana, Kizito K. Mugimba, Jessica Lukanga Nakavuma & Denis Karuhize Byarugaba

  2. Uganda Industrial Research Institute (UIRI), P.O. Box 7086, Kampala, Uganda

    John Bosco Omony

  3. Uganda National Council of Science & Technology (UNCST), P.O. Box 6884, Kampala, Uganda

    Maxwell Otim-Onapa

  4. Mbarara Zonal Agricultural Research and Development Institute (MBAZARDI), National Agricultural Research Organisation (NARO), P.O. Box 389,, Mbarara, Uganda

    Halid Kirunda

Authors
  1. John Bosco Omony
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Agnes Wanyana
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Kizito K. Mugimba
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Halid Kirunda
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Jessica Lukanga Nakavuma
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Maxwell Otim-Onapa
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Denis Karuhize Byarugaba
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Denis Karuhize Byarugaba.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Omony, J.B., Wanyana, A., Mugimba, K.K. et al. Disparate thermostability profiles and HN gene domains of field isolates of Newcastle disease virus from live bird markets and waterfowl in Uganda. Virol J 13, 103 (2016). https://doi.org/10.1186/s12985-016-0560-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-016-0560-0

Keywords

Virology Journal

ISSN: 1743-422X