Citation: Altaf Hussain, Tiantian Wu, Hui Li, Linjin Fan, Kai Li, Li Gao, Yongqiang Wang, Yulong Gao, Changjun Liu, Hongyu Cui, Qing Pan, Yanping Zhang, Asim Aslam, Khan Muti-Ur-Rehman, Muhammad Munir, Salman Latif Butt, Xiaomei Wang, Xiaole Qi. Pathogenic Characterization and Full Length Genome Sequence of a Reassortant Infectious Bursal Disease Virus Newly Isolated in Pakistan[J]. VIROLOGICA SINICA, 2019, 34 (1): 102-105 https://doi.org/10.1007/s12250-019-00082-8
Received: 31 July, 2018; Accepted: 28 November 2018; Published: 6 February 2019
Copyright: © Wuhan Institute of Virology, CAS 2019
Data Availability: All relevant data are within the paper and its Supporting Information files.
Altaf Hussain and Tiantian Wu these authors have contributed equally to this work.
Electronic supplementary material The online version of this article (https://doi.org/10.1007/s12250-019-00082-8) contains supplementary material, which is available to authorized users
Infectious bursal disease (IBD) is one of the most important diseases of the poultry.The IBD virus (IBDV), a nonenveloped virus belonging to the Birnaviridae family with a genome consisting of two segments of double-stranded RNA (segments A and B), targets B lymphocytes of bursa of Fabricious leading to immunosuppression.In Pakistan, poultry farming is the second biggest industry and IBD is the second biggest disease threating the poultry sector. However, there is limited genome information of IBDV available in Pakistan (Shabbir et al.2016).
In March 2017, a suspected IBD outbreak was reported in 28-day-old broiler chicken flock in Punjab, which is the most important chicken-farming area of Pakistan.A total of 30% to 50% morbidity and 13% mortality were observed without previous history of the vaccination.Seventeen (n=17) bursa samples were collected for virus detection.Using RT-PCR with primers of A628U/A1540L (Supplementary Table S1), a specific fragment of IBDV with an expected length of 930 bp was amplifiedfromall 17samples.Sequencingresults showed that all 17 sequences covering the hypervariable region (HVR) of VP2 carry 100% homology.Thereafter, one positive sample was randomly selected for virus isolation and one IBDV strain named PK2 was isolated successfully as described previously (Yuwen et al.2008).Furthermore, the genome cDNA of PK2 was cloned for segment A (primers of AU/ A1542L and A1421U2/AL2) or B (BU/B1344L and B1344U/ BL)(Supplementary Table S1).Sequencing results showed that the segment A of PK2 strain contained 3, 260 nucleotides (nt), including two partially overlapping ORFs (3, 073 nt) flanked by 50-non-coding-region (NCR, 96 nt) and 30-NCR (91 nt).The large ORF (3, 039 nt) of segment A encodes a polyprotein (NH2-VP2-VP4-VP3-COOH) while the small ORF (450 nt) encodes VP5.Segment B consists of VP1 coding ORF (2, 637 nt) flanked by 50-NCR (111 nt) and 30- NCR (79 nt).The genome sequence of PK2 is submitted to GenBank with the accession number MF996499(segment A) and MF996500(segment B), which is the first full-length genome sequence of Pakistan isolated IBDV strain appeared in GenBank.
Then the alignment and phylogenetic analyses based on the nucleotide or amino acid sequences were performed using the Clustal W (version 1.8)(Thompson et al.1997) and the MEGA (version 3.1)(Kumar et al.2004), and the confidence levels were assessed using 1000 bootstrap replications.The sequences of reference IBDV including very virulent, variant, attenuated, and serotype Ⅱ strains were obtained from GenBank.Phylogenetic analysis based on the amino acid sequence showed that the polyprotein and VP5 coded by segment A was closely related to very virulent IBDV (vvIBDV) strains (Fig. 1A), while segment-B-coded VP1 branched out of the vvIBDV group and was clustered into a unique branch (Fig. 1B).In addition to serotype Ⅱ, serotype I IBDV was further divided into three wild subtypes including classical, variants, very virulent strains, and one artificially attenuated strain.It has been reported that co-evolution of genome segments is a major evolutionary feature in IBDV (Le Nouen et al.2006).However, recently, it has been reported a few strains exhibited markedly different genetic relatives for segments A and B.(1) vv-A/Att-B (segment A from very virulent strain and segment B from attenuated strain):SH95, GX-NN-L, KZC-104, IBD13HeB01, MB11/ ABT/MVC/2016, JBN2011.(2) Att-A/vv-B:ZJ2000, TL2004, and HN04.(3) vv-A/c-B (segment A from vvIBDV and segment B from classic strain):Br/03/DR, CA-S7610.(4) vv-A/var-B (segment A from vvIBDV and segment B from variant strain):02015.2.(5) vv-A/Ⅱ-B (segment A from vvIBDV and segment B from serotype Ⅱ strains):CA-K785, CA-D495, and 100056.It was suggested that the exchange of the double-stranded genomic RNA segments contributed to the generation of new reassortant viruses because of the coinfection with different types of IBDV (Jackwood et al.2011; Le Nouen et al.2006), which have not been verified in laboratory.In fact, multiple factors may contribute to the gene reassortment, including host immunity, vaccinations, selection pressure, environment, and so on.It is identified that PK2 is a genetically segment-reassortant strain whose segment A was derived from vvIBDV, while segment B was originated from a unique ancestor (vv-A/Uniq-B).
In segment A, genetic distances derived from the amino acid sequences revealed that the polyprotein of PK2 shared 98.4%-99.8% identity with those of very virulent strains, while it has only 89.5%-97.9% identity with non-vvIBDV strains including classical, variant, attenuated, and serotype Ⅱ strains.Sixteenamino acid residuesin VP2(222A, 242I, 253Q, 256I, 279D, 284A, 299S, 330S, 451L), VP4(541I, 680Y, 685 N, 715S, 751D), and VP3(981P, 1005A), which are usually conserved in vvIBDV strains, were also present in PK2 (Fig. 1C).In addition, N-terminus quadruplet (MLSL) and the conserved amino acid residues at positions 49R, 78F in VP5 of vvIBDV (Lu et al.2015) were also observed in PK2(Fig. 1C). However, the amino acid residue at position 129 is P, which is the same with attenuated IBDV.PK2 VP5 has the amino acid sequence identity of 94.5%-96.6% with those of vv strains, while it has 93.2%-95.2% identity with non-vv strains.
In segment B, the homology analysis revealed that VP1 of PK2 had identity of 97.5%-98.3% with vvIBDV, while the identity with non-vv was 97.4%-98.7%.As was shown in Fig. 1C, among fifteen characteristic amino acid residues of VP1, PK2 has eight amino acids (4V, 61I, 145S, 287A, 508K, 511S, 546P, 646S) as same as vvIBDV, while the other seven amino acids (146E, 147G, 242D, 390L, 562S, 687S, and 695K) are the same as attenuated IBDV.The amino acid triplets at positions 145/146/147 of VP1 are important virulent sites through influencing the RNA dependant RNA polymerase (RdRp).It has been reported that there were three types of the triplet, TDN, TEG, NEG, which were the characteristic of very virulent, HLJ0504-like, and attenuated IBD (Gao et al.2014).Interestingly, VP1 of PK2 has the fourth kind of triplet (SEG), as being discovered in some Uruguay isolates (Hernández et al.2015).To be mentioned, VP1 of PK2 has 98%-98.8% identity with HLJ-0504, Gx, Harbin-1, and HeB10XS02, which was named as HLJ0504-like strain (He et al.2014; Qi et al.2011).In the group, PK2 shared 95.8%-96.6% amino acid sequence identity of polyprotein and 98%-98.8% of VP1 with other strains.It was reported that HLJ0504-like strains were the popular strains in China (He et al.2014; Qi et al.2011).Recently, this type of IBDV were also observed in Nigeria (Nwagbo et al.2016) and Algeria (Abed et al.2018).Segment B was reported to play key role in IBDV evolution through enhancing the virulence of the virus synergistically with its existing genome segment A (Hon et al.2006).However, the origin of segment B is unknown.Maybe the unique B was from other reservoir hosts like other poultry or wild birds (Jackwood et al.2005).
Furthermore, the serious pathogenicity of PK2 was evaluated.The 2 week-old SPF chickens housed separately in negative pressure isolators were randomly divided into two groups.The 9 chickens in the first group were individually infected with 1.26×103 ELD50(0.2 mL) of PK2 strain via eye and intranasal routes.The 3 chickens in the second group were provided with 0.2 mL of phosphate buffer solution as the negative control.Chickens were monitored daily for clinical signs for 7 days.The clinical signs were evaluated and showed as mean symptomatic index (MSI) as described previously (Le Nouen et al.2006). PK2 caused 100% morbidity and 55.6% mortality, whereas neither clinical signs nor mortality were observed in the PBS control group (Fig. 1D, 1E).From 2 days post-infection (d p.i.), chickens in PK2-infected group began showing typical signs of acute IBD (ruffled feathers, dehydration, and prostration), and the highest MSI (2.46 ± 0.35) was observed at 4 d p.i.(Fig. 1E).In PK2 group, the mortality peak (33.3%) occurred from 3 to 4 d p.i., and there were 2 chickens died on 6 d p.i.(Fig. 1D).At 7 d p.i., all surviving chickens were weighed and euthanized.The body weight of PK group (147.38 ± 11.32 g) was significantly lower than the control (204.43 ± 5.31 g)(P < 0.05)(Fig. 1F).IBDV was detected in the infected bursa with 6.4 9 108 viral RNA copies/g. PK2 caused severe histopathological bursal lesions including follicle atrophy, connective tissue hyperplasia, macrophages infiltration, and lymphocytic necrosis (Fig. 1H) with the bursa:body-weight index (BBIX) of 0.21 ± 0.04 (Fig. 1G).Meanwhile, the spleen/body weight ratio of PK group (1.15 ± 0.14) also had significant difference from the control (1.78 ± 0.20)(Fig. 1I).It is valuable to further study the pathogenic mechanism in more detail.
In conclusion, a naturally occurring reassortant strain of IBDV, carrying segment A from a very virulent strain and segment B from a unique ancestor, was first identified in Pakistan.For the first time, the full-length genome and the virulence of Pakistan IBDV was analyzed.The genomic and biological characterization of the strain clearly highlights the nature and dynamics of the IBDV in poultry sector of Pakistan.The present study therefore provides useful data for further understanding of the nature and evolution of IBDV, which will support a better prevention and control of the disease in Pakistan.
This work was supported by the National Key Research and Development Program of China (Nos. 2016YFE0203200, 2017YFD0500704), the Major Project of National Natural Science Foundation of China (No.314 30087), the Modern Agro-industry Technology Research System (No.CARS-41-G15).
Animal and Human Rights Statement
Animal experiments in this study were approved by the Ethics Committees of Harbin Veterinary Research Institute (HVRI), Chinese Academy of Agricultural Sciences (CAAS)(Approval Number:SQ-2017-080).
- . Abed M, Soubies S, Courtillon C, Briand FX, Allée C, Amelot M, De Boisseson C, Lucas P, Blanchard Y, Belahouel A. 2018. Infectious bursal disease virus in algeria: detection of highly pathogenic reassortant viruses. Infect Genet Evol, 60: 48-57.
- . Gao L, Li K, Qi X, Gao H, Gao Y, Qin L, Wang Y, Shen N, Kong X, Wang X. 2014. Triplet amino acids located at positions 145/146/147 of the RNA polymerase of very virulent infectious bursal disease virus contribute to viral virulence. J Gen Virol, 95: 888-897.
- . He X, Xiong Z, Yang L, Guan D, Yang X, Wei P. 2014. Molecular epidemiology studies on partial sequences of both genome segments reveal that reassortant infectious bursal disease viruses were dominantly prevalent in southern china during 2000-2012. Arch Virol, 159: 3279-3292.
- . Hernández M, Tomás G, Marandino A, Iraola G, Maya L, Mattion N, Hernández D, Villegas P, Banda A, Panzera Y. 2015. Genetic characterization of south american infectious bursal disease virus reveals the existence of a distinct worldwide-spread genetic lineage. Avian Pathol, 44: 212-221.
- . Hon CC, Lam TY, Drummond A, Rambaut A, Lee YF, Yip CW, Zeng F, Lam PY, Ng PT, Leung FC. 2006. Phylogenetic analysis reveals a correlation between the expansion of very virulent infectious bursal disease virus and reassortment of its genome segment b. J Virol, 80: 8503-8509.
- . Jackwood D, Gough R, Sommer S. 2005. Nucleotide and amino acid sequence analysis of a birnavirus isolated from penguins. Vet Rec, 156: 550-552.
- . Jackwood DJ, Sommer-Wagner SE, Crossley BM, Stoute ST, Woolcock PR, Charlton BR. 2011. Identification and pathogenicity of a natural reassortant between a very virulent serotype 1 infectious bursal disease virus (IBDV) and a serotype 2 IBDV. Virology, 420: 98-105.
- . Kumar S, Tamura K, Nei M. 2004. Mega 3: integrated software for molecular evolutionary genetics analysis and sequence alignment. Brief Bioinform, 5: 150-163.
- . Le Nouen C, Rivallan G, Toquin D, Darlu P, Morin Y, Beven V, de Boisseson C, Cazaban C, Comte S, Gardin Y. 2006. Very virulent infectious bursal disease virus: reduced pathogenicity in a rare natural segment-b-reassorted isolate. J Gen Virol, 87: 209-216.
- . Lu Z, Zhang L, Wang N, Chen Y, Gao L, Wang Y, Gao H, Gao Y, Li K, Qi X. 2015. Naturally occurring reassortant infectious bursal disease virus in northern china. Virus Res, 203: 92-95.
- . Nwagbo IO, Shittu I, Nwosuh CI, Ezeifeka GO, Odibo FJ, Michel LO, Jackwood DJ. 2016. Molecular characterization of field infectious bursal disease virus isolates from nigeria. Vet World, 9: 1420.
- . Qi X, Gao L, Qin L, Deng X, Wu G, Zhang L, Yu F, Ren X, Gao Y, Gao H. 2011. Genomic sequencing and molecular characteristics of a very virulent strain of infectious bursal disease virus isolated in china. Agr Sci Technol, 12: 1946-1949.
- . Shabbir MZ, Ali M, Abbas M, Chaudhry UN, Munir M. 2016. Molecular characterization of infectious bursal disease viruses from pakistan. Arch Virol, 161: 2001-2006.
- . Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The clustal_x windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res, 25: 4876-4882.
- . Yuwen Y, Gao Y, Gao H, Qi X, Li T, Liu W, Wang X. 2008. Sequence analysis of the VP2 hypervariable region of eight very virulent infectious bursal disease virus isolates from the northeast of china. Avian Dis, 52: 284-290.