Base usage and dinucleotide frequency have been extensively studied in many eukaryotic organisms and bacteria, but not for viruses. In this paper, a comprehensive analysis of these aspects for infectious bursal disease virus (IBDV) was presented. The analysis of base usage indicated that all of the IBDV genes possess equivalent overall nucleotide distributions. However when the base usage at each codon positions was analysed by using cluster analysis, the VP5 open reading frame (ORF) formed a different cluster isolated from the other genes. The unusual base usage of VP5 ORF may indicate that the gene was originated by the virus "overprinting strategy", a strategy in which virus may create novel gene by utilizing the unused reading frames of its existing genes. Meanwhile, the GC content of the IBDV genes and the chicken's coding sequences was comparable; suggesting the virus imitation of the host to increase its translational efficiency. The analysis of dinucleotide frequency indicated that IBDV genome had dinucleotide bias: the frequencies of CpG and TpA were lower and the TpG was higher than the expected. Classical methylation pathway, a process where CpG converted to TpG, may explain the significant correlation between the CpG deficiency and TpG abundance. "Principal component analysis of the dinucleotide frequencies" (DF-PCA) was used to analyse the overall dinucleotide frequencies of IBDV genome. DF-PCA on the hypervariable region and polyprotein (VPX-VP4-VP3) gene showed that the very virulent IBDV (vvIBDV) was segregated from other strains; which meant vvIBDV had a unique dinucleotide pattern. In summary, the study of base usage and dinucleotide frequency had unravelled many overlooked genomic properties of the virus.
Three isolates of Infectious bursal disease virus (IBDV), designated UPM04178, UPM04190 and UPM04238, were obtained from severe outbreaks of infectious bursal disease (IBD) in Malaysia in 2004. The hypervariable region (HPVR) of VP2 gene of these isolates was sequenced. The obtained sequences were compared with those of other isolates. The highest similarity (98%) concerning both nucleotide and amino acid sequences was found to very virulent IBDV (vvIBDV) strains. Phylogenetic analysis revealed clustering of the three isolates with vvIBDV strains. Evolutionary relatedness of the three isolates to vvIBDV strains was demonstrated by three phylogenetic methods: bootstrap values of 100%, 95% and 90% for nucleotide sequences and those of 58%, 86% and 96% for amino acid sequences were obtained by the distance, maximum parsimony and maximum likehood methods, respectively. It is concluded that UPM04178, UPM04190 and UPM04238 are vvIBDV isolates of serotype 1, which originate from a common ancestor of IBDV strains present in Malaysia.
The current available molecular method to detect infectious bursal disease virus (IBDV) is by reverse transcriptase-polymerase chain reaction (RT-PCR). However, the conventional PCR is time consuming, prone to error and less sensitive. In this study, the performances of Sybr Green I real-time PCR, enzyme-linked immunosorbent assay (ELISA) and conventional agarose detection methods in detecting specific IBDV PCR products were compared. We found the real-time PCR was at least 10 times more sensitive than ELISA detection method with a detection limit of 0.25pg. The latter was also at least 10 times more sensitive than agarose gel electrophoresis detection method. The developed assay detects both very virulent and vaccine strains of IBDV but not other RNA viruses such as Newcastle disease virus and infectious bronchitis virus. Hence, Sybr Green I-based real-time PCR is a highly sensitive assay for the detection of IBDV.
The VP2 hypervariable region of P97/302 local infectious bursal disease virus (IBDV) isolate was amplified by the reverse transcriptase (RT) nested polymerase chain reaction (PCR) and cloned. This region of P97/302 local isolate was sequenced and compared with eight other reported IBDV sequences. The result showed that P97/302 IBDV was most identical to the reported very virulent IBDV strains because it has amino acid substitutions at positions 222, 256, 294, and 299, which encode alanine, isoleucine, isoleucine, and serine, respectively. This region can be digested with restriction enzymes of Taq1, Sty1, Ssp1 but not with Sac1. The P97/302 isolate was then used for the optimization of RT nested PCR enzyme-linked immunosorbent assay (ELISA). The RT nested PCR ELISA was able to detect 10(-4) dilution of the infected bursa homogenates and was 10 times more sensitive when compared with the agarose gel detection method. The RT nested PCR ELISA can detect up to 0.48 ng of the PCR product. The specificity of this nested PCR ELISA was also high (100%).
Specific-pathogen-free (SPF) chickens infected with very virulent (vv) infectious bursal disease virus (IBDV) UPM94/273 developed lower pathogenicity compared to UPM97/61. Sequence analysis indicated that UPM94/273 is an exceptional vvIBDV. In this study, a SYBR Green I based real-time reverse transcriptase reaction assay was developed to measure viral RNA in the bursae of SPF chickens infected with IBDV. Specificity of the amplified products was confirmed by melting temperature analysis. A linear relationship was observed between the amount of input viral RNA and the threshold values for IBDV-specific product over five log10 dilutions. The viral RNA level following infection with UPM94/273 was significantly higher at day 1 and 2 post-inoculation (p.i.) compared to UPM97/61 infected chickens. However, chickens infected with UPM97/61 had significantly higher numbers of bursal cells undergoing apoptosis compared to UPM94/273 infected chickens. In both groups, the number of apoptotic cells and viral RNA levels peak at day 3 p.i. This study indicates that UPM97/61 and UPM94/273 have different efficiency of replication and percentage of apoptotic cells in bursae during the acute phase of IBDV infection.
Two Infectious bursal disease virus (IBDV) isolates, NP1SSH and NP2K were obtained from a severe infectious bursal disease (IBD) outbreak in Nepal in 2002. The hypervariable (HV) region of VP2 gene (1326 bp) of the isolates was generated by RT-PCR and sequenced. The obtained nucleotide sequences were compared with those of twenty other IBDV isolates/strains. Phylogenetic analysis based on this comparison revealed that NP1SSH and NP2K clustered with very virulent (vv) IBDV strains of serotype 1. In contrast, classical, Australian classical and attenuated strains of serotype 1 and avirulent IBDV strains of serotype 2 formed a different cluster. The deduced amino acid sequences of the two isolates showed a 98.3% identity with each other and 97.1% and 98.3% identities, respectively with very virulent IBDV (vvIBDV) isolates/strains. Three amino acids substitutions at positions 300 (E-->A), 308 (I-->F) and 334 (A-->P) within the HV region were common for both the isolates. The amino acids substitutions at positions 27 (S-->T), 28 (I-->T), 31 (D-->A), 36 (H-->Y), 135 (E-->G), 223 (G-->S), 225 (V-->I), 351 (L-->I), 352 (V-->E) and 399 (I-->S) for NP1SSH and at position 438 (I-->S) for NP2K were unique and differed from other IBDV isolates/strains. NP1SSH and NP2K showed highest similarity (97.8%) with the BD399 strain from Bangladesh as compared with other vvIBDV isolates/strains. We conclude that the NP1SSH and NP2K isolates of IBDV from Nepal represent vvIBDV of serotype 1.
The complete nucleotide sequences encoding precursor polyprotein (VP2-VP3-VP4) and VP5 of a highly virulent (hv) infectious bursal disease virus (IBDV), UPM97/61 was determined. Comparison of the deduced amino acid sequences with the published ones revealed 8 common amino acid substitutions, which were found only in the hv IBDV including the UPM97/61 strain. Three of the amino acid substitutions (222 Ala, 256 Ile and 294 Ile) were used as a marker for determining hv IBDV strains. The other five substitutions (685 Asn, 715 Ser, 751 Asp, 990 Val and 1005 Ala) were also conserved in hv IBDV strains isolated in various countries. UPM97/61 strain demonstrated also 8 unique amino acid substitutions of which 3 were in VP2, 4 in VP3 and 1 in VP4. There was 1 unique amino acid substitution in VP5 at position 19 (Asp-->Gly) not found in other strains. However, all the strains have a conserved 49 Arg. The amino acid sequence of UPM97/61 strain differed by 1.09% from the Japanese (OKYM) and Hong Kong (HK46) strains, and by 1.48% from the Israeli (IBDVKS) and European (UK661) strains. Hence, UPM97/61 is more closely related to the hv strains from Asia. However, phylogenetic analysis indicated that the origin of UPM97/61 might be the same as that of other hv strains isolated from other parts of the world.
Previously we have shown that very virulent infectious bursal disease viruses (vvIBDV) that are SspI, TaqI and StyI positive (92/04, 97/61 and 94/B551) but not SspI and TaqI positive and StyI negative (94/273) cause high mortality, up to 80% in specific-pathogen-free chickens with significant damage of the bursal as well as nonbursal tissues. In this study, we sequenced the VP2 gene (1351 bp) of the 92/04, 94/273 and 94/B551 and compared them with other IBDV strains. All the isolates have the unique amino acid residues at positions 222A, 256I, 294I and 299S found in other vvIBDV strains. The deduced VP2 amino acids encoded by 92/04 is identical to the vvIBDV strains from Israel (IBDVKS), Japan (OKYM) and Europe (UK661), whereas the 94/273 and 94/B551 isolates have one to three amino acid substitutions. The 94/273 has two amino acid substitutions at positions 254 G to S and at 270 A to E that have not been reported before from vvIBDV strains. The 94/B551 also has one amino acid substitution at position 300 E to S, which is uncommon among other vvIBDV strains. However, phylogenetic analysis suggested that the isolates are very close to each other and all of them may have derived from the same origin as vvIBDV strains isolated from China, Japan and Europe. Even though antigenic index analysis of the 94/273 and 94/B551 indicated that the isolates are unique compared to other IBDV strains, their antigenic variation remain to be determined by monoclonal antibody study.
The deduced amino acid sequences of segment A and B of two very virulent Infectious bursal disease virus (vvIBDV) isolates, UPM94/273 and UPM97/61 were compared with 25 other IBDV strains. Twenty amino acid residues (8 in VP1, 5 in VP2, 2 in VP3, 4 in VP4, 1 in VP5) that were common to vvIBDV strains were detected. However, UPM94/273 is an exceptional vvIBDV with usual amino acid substitutions. The differences in the divergence of segment A and B indicated that the vvIBDV strains may have been derived from genetic reassortment of a single ancestral virus or both segments have different ability to undergo genetic variation due to their different functional constraints.
Due to the limitations in the range of antibodies recognising avian viruses, quantitative real-time PCR (RT-qPCR) is still the most widely used method to evaluate the expression of immunologically related genes in avian viruses. The objective of this study was to identify suitable reference genes for mRNA expression analysis in chicken intraepithelial lymphocyte natural killer (IEL-NK) cells after infection with very-virulent infectious bursal disease virus (vvIBDV). Fifteen potential reference genes were selected based on the references available. The coefficient of variation percentage (CV%) and average count of these 15 genes were determined by NanoString technology for control and infected samples. The M and V values for shortlisted reference genes (ACTB, GAPDH, HMBS, HPRT1, SDHA, TUBB1 and YWHAZ) were calculated using geNorm and NormFinder. GAPDH, YWHAZ and HMBS were the most stably expressed genes. The expression levels of three innate immune response related target genes, CASP8, IL22 and TLR3, agreed in the NanoString and RNA sequencing (RNA-Seq) results using one or two reference genes for normalisation (not HMBS). In conclusion, GAPDH and YWHAZ could be used as reference genes for the normalisation of chicken IEL-NK cell gene responses to infection with vvIBDV.
Studies have shown that infectious bursal disease virus (IBDV) infects lymphoid cells, mainly B cells and macrophages. This study was aimed to examine the involvement of chicken splenic-derived dendritic cells (ch-sDCs) in specific-pathogen-free chickens following inoculation with IBDV vaccine strain (D78) and a very virulent (vv) strain (UPM0081). Following IBDV infection, enriched activated ch-sDCs were collected by using the negative selection method and were examined based on morphology and immunophenotyping to confirm the isolation method for dendritic cells (DCs). The presence of IBDV on enriched activated ch-sDCs was analyzed based on the immunofluorescence antibody test (IFAT), flow cytometry, and quantitative real-time PCR (RT-qPCR) while the mRNAs of several cytokines were detected using RT-qPCR. The isolated ch-sDCs resembled typical DC morphologies found in mammals by having a veiled shape and they grew in clusters. Meanwhile, the expression of DC maturation markers, namely CD86 and MHCII, were increased at day 2 and day 3 following vvIBDV and vaccine strain inoculation, respectively, ranging from 10% to 40% compared to the control at 2.55% (P < 0.05). At day 3 postinfection, IBDV VP3 proteins colocalized with CD86 were readily detected via IFAT and flow cytometry in both vaccine and vvIBDV strains. In addition, enriched activated ch-sDCs were also detected as positive based on the VP4 gene by RT-qPCR; however, a higher viral load was detected on vvIBDV compared to the vaccine group. Infection with vaccine and vvIBDV strains induced the enriched activated ch-sDCs to produce proinflammatory cytokines and Th1-like cytokines from day 3 onward; however, the expressions were higher in the vvIBDV group (P < 0.05). These data collectively suggest that enriched activated ch-sDCs were permissive to IBDV infection and produced a strong inflammatory and Th1-like cytokine response following vvIBDV infection as compared to the vaccine strain.
Very virulent infectious bursal disease virus (vvIBDV) targets B lymphocytes in the bursa of Fabricius (BF), causing immunosuppression and increased mortality rates in young birds. There have been few studies on the host immune response following vvIBDV infection at different inoculum doses in chickens with different genetic backgrounds. In this study, we characterized the immune responses of specific-pathogen-free (SPF) chickens and Malaysian red jungle fowl following infection with vvIBDV strain UPM0081 at 103.8 and 106.8 times the 50% embryo infectious dose (EID50). The viral burden, histopathological changes, immune cell populations, and expression of immune-related genes were measured and compared between infected and uninfected bursa at specific intervals. The populations of KUL1+, CD3+CD4+ and CD3+CD8+ cells were significantly increased in both types of chickens at 3 dpi, and there was significant early depletion of IgM+ B cells at 1 dpi in the red jungle fowl. vvIBDV infection also induced differential expression of genes that are involved in Th1 and pro-inflammatory responses, with groups receiving the higher dose (106.8 EID50) showing earlier expression of IFNG, IL12B, IL15, IL6, CXCLi2, IL28B, and TLR3 at 1 dpi. Although both chicken types showed equal susceptibility to infection, the red jungle fowl were clinically healthier than the SPF chickens despite showing more depletion of IgM+ B cells and failure to induce IFNB activation. In conclusion, high-dose vvIBDV infection caused an intense early host immune response in the infected bursa, with depletion of IgM+ B cells, bursal lesions, and cytokine expression as a response to mitigate the severity of the infection.
A SYBR Green I based one-step real-time reverse transcriptase polymerase chain reaction was developed for the detection and differentiation of very virulent (vv) and classical strains of infectious bursal disease virus (IBDV). The assay showed high PCR efficiency >93% and high reproducibility with coefficient of variation less than 0.5%. When tested on characterized IBDV strains, the very virulent and classical-specific primers detected accurately only vvIBDV and classical IBDV strains, respectively. The diagnostic efficacy of the assay was also tested on 140 bursal samples from experimental infection and 37 bursal samples from cases suspected of IBD. The assay was able to detect IBDV from bursal samples collected at days 3 and 5 post-infection with the vvIBDV strain UPM94/273 and the classical IBDV strain D78. The assay was also able to detect bursal samples infected dually with D78 and UPM94/273. The melting temperature values of the amplification products from the classical and very virulent viral infection were statistically significant (P<0.05). The specificity of the assay for detecting IBDV from suspected cases was confirmed by sequence analysis of the VP2 gene. The assay showed high sensitivity since bursal samples which were negative for IBDV were confirmed by virus isolation and PCR amplification. Hence, the new assay offers an attractive method for rapid detection of strains of IBDV.
The ability of a heat-inactivated whole virus from a highly virulent infectious bursal disease virus (hvIBDV) and VP2 protein from hvIBDV expressed in E. coli provided protection against a hvIBDV challenge in specificpathogen-free (SPF) chickens. Six out of seven chickens that were injected three times with crude VP2 protein developed significant antibody titer against IBDV. However, only four out of the seven chickens survived the hvIBDV challenge. Despite showing low antibody titer profiles, all chickens immunized with the heat-inactivated whole virus also survived the challenged with hvIBDV. However, all of these chickens had bursal atrophy and mild to moderate depletion of lymphocytes. Thus, antibodies raised against IBDV VP2 protein expressed in E. coli and denatured IBDV proteins induced some degree of protection against mortality but not against bursal damage following challenge with hvIBDV.
The characteristics of the pathogenic infectious bursal disease virus (IBDV) that infected avian species other than commercial chickens were largely unknown. In this study, by using in vivo and molecular methods, we had characterized an IBDV isolate (named 94268) isolated from an infectious bursal disease (IBD) outbreak in Malaysian village chickens--the adulterated descendant of the Southeast Asian jungle fowl (Gallus bankiva) that were commonly reared in the backyard. The 94268 isolate was grouped as the very virulent IBDV (vvIBDV) strain because it caused severe lesions and a high mortality rate in village chickens (>88%) and experimentally infected specific-pathogen-free chickens (>66%). In addition, it possessed all of the vvIBDV molecular markers in its VP2 gene. Phylogenetic analysis using distance, maximum parsimony, and maximum likelihood methods revealed that 94268 was monophyletic with other vvIBDV isolates and closely related to the Malaysian vvIBDV isolates. Given that the VP2 gene of 94268 isolate was almost identical and evolutionarily closely related to other field IBDV isolates that affected the commercial chickens, we therefore concluded that IBD infections had spread across the farm boundary. IBD infection in the village chicken may represent an important part of the IBD epidemiology because these birds could harbor the vvIBDV strain and should not be overlooked in the control and prevention of the disease.