Detection and Phylogenetic Analysis of Infectious Bursal Disease Virus Based on both Genome Segments A and B Isolated from Backyard Poultry Punjab, Pakistan

Syeda Fakhra Waheed1, Asim Aslam1*, Muti-ur-Rehman Khan1,

Kamran Ashraf2 and Beenish Zahid3

1Department of Pathology, Faculty of Veterinary and Animal Sciences, University of Veterinary and Animal Sciences, Lahore, 54000, Pakistan

2Department of Parasitology, Faculty of Veterinary and Animal Sciences, University of Veterinary and Animal Sciences, Lahore, 54000, Pakistan

3Department of Pathobiology, KBCMA, CVAS Narowal, University of Veterinary and Animal Sciences, Lahore Pakistan

ABSTRACT

In recent years, the re-emergence of virulent strains of infectious bursal disease virus (IBDV) has resulted in substantial economic losses in Pakistan despite mass and intense vaccination regimens. Infectious bursal disease (IBD) is one of the most important immunosuppressive diseases of the poultry and has been a constraint on sustainable food security around the globe, including Pakistan. This disease damages the bursa of Fabricius (BF), which causes immunosuppression in birds. A total of 50 tissue samples from backyard chicken flocks presenting suspected symptoms were collected during February 2020 and March 2021 for genetic characterization, followed by phylogentic analysis. A total of 4 isolates were sequenced based on partial VP1 and hyper-variable region of the VP2 genes simultaneously. According to phylogenetic analysis, the study isolates genotype A3B3 were identified as predominant strains in country backyard poultry. Concerning the identity matrix analysis of VP1 representative part, study isolates shared (87-88% nt; 95-96% aa) identity with the vvIBDV and (97-98% nt; 98-99% aa) with non-vvIBDV. While VP2 revealed (99-100% nt; 99-100% aa) identity with previously reported Pakistan vvIBDV strains and (91-95% nt; 91-95% aa) with non vvIBDV. Amino acid alignment analysis of VP1 revealed that current IBDVs have three characteristic aa residues of vvIBDV (287-A, 508-K, and 511-S) and four characteristic aa residues of non-vvIBDV (146-E, 147-G, 242-D and 390-L). While VP2 gene sequence alignment revealed eight characteristic aa residues of vvIBDV (222-A, 242-I, 253-Q, 256-I, 279-D, 284-A, 299-S, and 330-S and a distinct aa 384-I). Based on phylogeny, this is the first identification of IBDV segment reassortants having segment A of A3 (very virulent) and segment B of B3 (early Australian-like) genogroups reported in Pakistan backyard poultry. Further study is required to determine the pathogenicity of the IBDV reassortant and development of new policies for IBDV intervention in the country.

Article Information

Received 30 December 2021

Revised 20 March 2022

Accepted 05 April 2022

Available online 23 May 2022

(early access)

Published 16 March 2023

Authors’ Contribution

SFW conducted the research. AA supervised the study and AA, MRK and KA designed the study. SFW wrote the manuscript. BZ helped in interpretation of results and article writing.

Key words

Infectious bursal disease, Backyard poultry, Phylogenetic analysis, Pakistan

DOI: https://dx.doi.org/10.17582/journal.pjz/20211230141208

* Corresponding author: drasimaslamch@uvas.edu.pk

0030-9923/2023/0003-1249 $ 9.00/0

Copyright 2023 by the authors. Licensee Zoological Society of Pakistan.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

INTRODUCTION

Before 1960s backyard poultry production was the sole poultry meat source in the country, with flock size of about 20-50 birds of different breeds and ages. It was a low input business without veterinary care facilities (Ahmed et al., 2021). Infectious bursal disease (IBD) is a highly contagious and immunosuppressive disease of young chickens that causes considerable economic losses to the poultry industry globally (Qin and Zheng, 2017; Hussain et al., 2020; Maqbool et al., 2020). The disease is caused by the infectious bursal disease virus (IBDV), a non-enveloped dsRNA virus, which belongs to the genus Avibirnavirus of family Birnaviridae with two distinct serotypes (Müller et al., 2003). Serotype I is pathogenic and induces pathological lesions in chickens, whereas serotype II is non-pathogenic, primarily isolated from turkeys. Pathogenic IBDV is classified into four subtypes, i.e., attenuated, classical virulent, antigenic variants, and very virulent IBDV (Hussain et al., 2020). IBDV targets immature B lymphocytes (Zahid et al., 2018) in the lymphoid organ (bursa of Fabricius), leading to age-dependent immunosuppression exposing birds to other opportunistic pathogens. Very virulent IBDV (vvIBDV) has been reported in different countries with more than 70% chicken mortality, posing a severe threat to the poultry industry (Dey et al., 2019).

IBDV has a bi-segmented genome (segment A and B) of approximately 6.0 kb. Segment A is 3.2 kb long, with two overlapping open reading frames (ORF) encoding four proteins. ORF-1 encodes longer precursor polyprotein, cleaved by autoproteolysis and yields capsid protein VP2, nucleoprotein VP3 and protease VP4. In comparison, ORF-2 encodes only a nonstructural VP5 protein, a significant determinant of viral virulence and host protective antigen (Luque et al., 2009). The VP2 protein is the primary protective antigen containing specific epitopes responsible for inducing neutralizing antibody responses (Reddy et al., 2017). Amino acid substitutions in the VP2 hypervariable region can lead to variations in antibody recognition, antigenicity, virulence, and tissue tropism (Mato et al., 2020). Segment B is 2.8kb long, and it contains only one ORF which encodes RNA-dependent RNA polymerase VP1 protein, involved in virus replication and evolution (Ye et al., 2018). Both VP1 and VP2 proteins have been used to study pathogenicity, virulence, phylogenetic relationship and molecular characterizations (Fan et al., 2019).

IBDV may undergo segment reassortment due to the genome’s bi-segmented structure. It has been proposed that vvIBDV arose through segment reassortment between segment A of circulating IBDV strains and segment B from an unknown IBDV strain (Hon et al., 2006). Natural reassortant IBDV strains have been reported frequently in recent years from various parts of the world, including reassortment between serotypes 1 and 2 (Soubies et al., 2017; Stoute et al., 2019) and reassortment between different genotypes of serotype 1 IBDVs (Mato et al., 2020).

In 1980 a vvIBDV outbreak was first recorded with high mortality in Punjab Pakistan (Khan et al., 1988) and since then continues to pose a threat by various strains all over the country (Lone et al., 2009; Zahoor et al., 2011). Despite extensive vaccination programs, IBD still prevails at the poultry farms posing significant economic losses in Pakistan, which emphasizes the characterization of field circulating IBDV strains to understand the causes of vaccination failure, hence incorporating local viruses in vaccine development (Shafqat et al., 2017; Sajid et al., 2021). Traditionally, nucleotide sequencing of the segment A coded VP2 hypervariable region is used for genotypic characterization of IBDV strains. Both segments should be characterized to detect genetic variations, influencing strain pathogenicity (Gao et al., 2018). Both genome segments contribute to the virulence of IBDV and segment reassortment plays a role in virus evolution, which is not possible to be evaluated by targeting VP2 gene only. Since VP1 and VP2 proteins are involved in virus evolution, there is still insufficient information to predict the actual dynamics of IBDV in Pakistan. Recently, some authors have characterized the chicken isolated IBDV VP1 and VP2 genes (Hussain et al., 2019, 2020), but the literature on backyard poultry is silent. Therefore, this study was designed for IBDV genotyping and phylogenetic analysis covering both segments from indigenous breed of backyard poultry. The study’s primary objective was to detect and genotype IBDV strains circulating in country backyard poultry and add a new piece of knowledge about its characterization to better implement control measures.

Materials and Methods

Ethical approval

The present work was approved by the ethical review committee of the University of Veterinary and Animal Sciences, Lahore, Pakistan (Permit Number: DR/944).

Study area

The current study was conducted at the Department of Pathology, University of Veterinary and Animal Sciences Lahore, Pakistan.

Clinical samples

A total number of 50 tissue samples (bursa of Fabricius) were collected from backyard chicken flocks for genome extraction from February 2020 to March 2021. The dead bird subjected to virus isolation had typical hemorrhagic lesions on skeletal muscles and mucosal surface of the bursa of Fabricius. The bursal tissues were ground to make homogenate in phosphate buffered saline (pH 7.2) containing penicillin and streptomycin. The homogenate was finally centrifuged at 5 000×g for 5 min at 4°C. Lastly, the upper aquous layer was collected for further detection.

Molecular identification

The virus genome was extracted according to the protocol as described by Qi et al. (2015). Viral RNA was extracted from the bursa homogenate using QIAamp Viral RNA Kit (Qiagen, Germany) as directed by the manufacturer. Viral cDNA was synthesized using RevertAid First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, USA) (Zahid et al., 2016). The selected DNA fragment was amplified through PCR in a 25µL reaction mixture, containing 12.5µL master mix (PrimeSTAR Max DNA polymerase, Catalog Number: R045A), 1µL DNA template, 1µL each forward and reverse primers and 9.5µL PCR grade water. Previously reported VP1 gene-based primers; F(5’-AGGAGAAGCCCAATGCGT-3’) and R(5’-GTCATCAATGGACCTCTC-3’) with amplicon size 1255bp (Hussain et al., 2020); VP2 gene-based primers; F(5’-AGCCAACATCAACGACAAA-3’) and R(5’-CAAGACGGTCCCTCTCACT-3’) with an estimated amplicon size of 782bp were employed (Ali et al., 2019). The PCR was performed under the following conditions: preheating (95℃ for 5 min), with subsequent 35 cycles of denaturation (95℃ for 30 s), annealing (56℃ for 30 s) and extension (72℃ for 1 min) followed by a final extension at 72℃ for 10 min. The RT-PCR products were checked by 1% agarose gel, stained with GelRed™ Nucleic Acid Gel Stain in gel electrophoresis and visualized in a gel documentation system (Bio-Rad Laboratories, United States). Amplified PCR product was purified using GeneJET PCR Purification Kit (Thermo Fischer Scientific, USA). The purified samples were sent to the Comate Bioscience Co., Ltd., China for Sanger sequencing.

Statistical and phylogenetic analysis

The obtained oligonucleotide sequences were properly trimmed with Chromas Software (version 2.6.6). Multiple sequence alignments were performed on BioEdit Software (version 7.2.5) using the ClustalW (version 1.8) (Thompson et al., 1997). The evolutionary analysis was conducted on MEGA Software (version 6) inferred through the Neighbor-Joining method with the maximum likelihood method. The reliability of the tree was tested by performing 1000 bootstrap replicates (Tamura et al., 2013). The pairwise nucleotide (nt) and amino acid (aa) distance matrix were calculated with UGENE software.

RESULTS AND DISCUSSION

In the present study, the results of RT-PCR revealed that the positive rate of IBDV samples was 32% (16/50). Afterwards, four positive samples were selected randomly to sequence the partial VP1 and VP2 gene fragments and the nucleotide sequences were determined and submitted to NCBI GenBank database (Table I). The phylogenetic analysis revealed that serotype I was subgrouped into vvIBDV, classical, attenuated, and variant IBDV. All the vvIBDV isolates, including reference strain (NC_004178), were placed in a single clade with a bootstrap value of 100. In contrast, classical, attenuated, and variant IBDV were grouped into another monophyletic clade. In the evolutionary tree based on the 89 nucleotide sequences of VP1 retrieved from GenBank, serotype I strains revealed distinct branches and current study isolates with previously reported Pakistan strains were branched out of the vvIBDV branch and formed a cluster with HLJ0504-like strains, including HLJ0504 and Harbin-1 (Fig. 1). Based on the 72 nucleotide sequences of VP2, the phylogenetic tree was distinctly divided into two major branches, vvIBDV and non-vvIBDv. Present study isolates belonged to vvIBDV branch with already reported isolates (Pakistan reference strains) and formed one subgroup with a bootstrap value of 100 (Fig. 2). Segment A-coded VP2 gene phylogenetic analysis exhibited that our isolate had a close relationship with other Pakistan and Indian isolates. Currently identified isolates acquired both segments (A and B) from different genetic relatives and were further confirmed by sequence alignment and phylogeny. VP2 gene (segment A) of current isolates was clustered into vvIBDV branch while VP1 gene (segment B) branched out of vvIBDV and formed a unique branch with strain like HLJ0504 previously reported in China (Hussain et al., 2019; Qi et al., 2011), Algeria (Nwagbo et al., 2016) and Pakistan (Hussain et al., 2019). Strains like HLJ0504 and Harbin-1 with unique segment (segment A originated from vvIBDV) and segment B acquired from a unique ancestor. So the current study reported unique segment-reassortant IBDV from backyard poultry. The results agree with (Hussain et al., 2019, 2020), who reported Pakistan originated reassortant strains from commercial poultry. Concerning the identity matrix analysis of VP1 representative part, study isolates shared (87-88% nt; 95-96% aa) homology with the vvIBDV and (97-98% nt; 98-99% aa) with non-vvIBDV. VP1 of isolates have comparatively higher identity of (97-98% nt; 98-99% aa) with previously reported HLJ0504-like strains which are the unique reassortant retrieving segment A from vvIBDV and B from unique ancestor (He et al., 2014; Qi et al., 2011). Concerning the identity matrix analysis of VP2 representative part, present study isolates revealed 99-100% nt and 99-100% aa identity with previously reported Pakistan vvIBDV strains and 91-95% nt and 91-95% aa homology with non vvIBDV.

 

Table I. Viral strains isolated from backyard chicken (age 3-5 weeks) during Apil and May 2020 and included in current study.

Gene	Strain	GenBank accession No.
Partial VP1	PAK-SFW PAK-SFW2 PAK-SFW3 PAK-SFW4	MW541059 OL690424 OL690425 OL690426
Partial VP2	PAK-SFW PAK-SFW2 PAK-SFW3 PAK-SFW4	MW452665 OL690421 OL690422 OL690423

 

 

In VP1 representative part sequence alignment, Pakistan IBDVs has three characteristic aa residues of vvIBDV (287-A, 508-K, and 511-S) and four characteristic aa residues of non-vvIBDV (146-E, 147-G, 242-D and 390-L) (Fig. 3). There are three domains in VP1 of IBDV; N-terminal (1-167aa), the polymerase (168-658aa), and C terminal domain (659-878aa). Pakistan IBDV isolates revealed recombinant characters in the N-terminal domain and the central polymerase domain between very virulent and attenuated strains. Self guanylylation site of RNA dependent RNA polymerase responsible for Protein priming of stains found in the N terminal domain (1–67 aa). Additionally, sequence alignment of VP2 revealed eight characteristic aa residues of vvIBDV (222-A, 242-I, 253-Q, 256-I, 279-D, 284-A, 299-S, and 330-S and a distinct aa 384-I) (Fig. 4) also presented in Pakistan IBDVs previously reported by (Hussain et al., 2020). Likewise the reassortants are also prevailing in the neighboring countries including India and China, it is thought-provoking to discover the relation among them.

As far as genotype is concerned current isolates belong to A3B3 genotype (Table II) according to the new classification scheme proposed by Wang et al. (2021) and

 

 

 

binomial scheme for genotype characterization suggested by Islam et al. (2021). In Serotype 1, segment A has genogroups classical (A1), US antigenic variant (A2), very virulent (A3), early European and recent South American distinct IBDV (A4), atypical or recombinant strains (A5), atypical Italian (A6), early Australian (A7) and Australian variant (A8). While segment B has five distinct genogroups; classical (B1), very virulent (B2), early Australian (B3), Polish and Tanzanian (B4), and Nigerian (B5). Segment A of the current four isolates clustered with vvIBDV genogroup A3, additionally, segment B of current isolates clustered with early Australian like IBDVs genogroup B3. The findings suggests that current isolates are segment reassortants with very virulent segment A and early Australian like segment B. Isolates with the same genotype A3B3 were also detected in Bangladesh in 2016. Different genotypes have appeared due to reassortment between two segments of previously circulating or newly evovled genotypes. The characterization of IBDV strains based on both genome segments information will help for a better understanding of IBDV evolution and molecular epidemiology.

CONCLUSION

The naturally occurring reassortant strains of IBDV genotype A3B3, harboring segment A from a very virulent strain and segment B from one unique ancestor, were identified through molecular characterization. To our knowledge, this is the first report on characterization of segment B coding VP1 and segment A-coding VP2 genes in Pakistan backyard poultry. The current work gives insight into a better knowledge of the nature and evolution of IBDV, which will assist in disease prevention and control programs for emerging reassortant strains.

 

Table II. Classification of IBDV isolates by phenotype and genogroups.

(Traditional classification) phenotype	Genotype	Reference strain	Accession No.
			A		B
Classic strains	A1B1	F52/70 France	HG974565	HG974566
		IM USA	AY029166	AY029166
Variant strains	A2aB1	Variant E USA	AF133904	AF133905
	A2bB1	9109 USA	AY462027	AY459321
	A2cB1	GLS	AY368653	AY368654
Novel variant strains	A2dB1	SHG19 China	MH879045	MH879092
		SHG120 China	MH879063	MH879110
		SHG350 China	MH879081	MH879129
		SD19-9901 China	MT087547	MT087554
		SD19-9903 China	MT087548	MT087555
		SD19-9904 China	MT087549	MT087556
		JS19-13202 China	MT087550	MT087557
		JS19-13203 China	MT087551	MT087558
		JS19-13804 China	MT087552	MT087559
		JS19-13902 China	MT087553	MT087560
Very virulent strains	A3B2	UK661 France	NC_004178	NC_004179
		OKYM Japan	D49706	D49707
		89163 France	HG974563	HG974564
Attenuated strains	A8B1	Cu-1 Germany	X16107	AF362775
		CT France	AJ310185	AJ310186
		Gt China	DQ403248	DQ403249
Reassortant strains	A3B3	Harbin-1	EF517528	EF517529
		Gx	AY444873	AY705393
		PAK-SFW	MW452665	MW541059
		PAK-SFW2	OL690421	OL690424
		PAK-SFW3	OL690422	OL690425
		PAK-SFW4	OL690423	OL690426

 

Statement of conflict of interest

The authors have declared no conflict of interest.

References

Ahmed, T., Ameer, H.A. and Javed, S., 2021. Pakistan’s backyard poultry farming initiative: impact analysis from a public health perspective. Trop. Anim. Hlth. Prod., 53: 1-12. https://doi.org/10.1007/s11250-021-02659-6

Ali, K.R.S., Habib, M., Ali, W., Salah, D.S.M., Ashraf, A., Ali, T.Z., Helal, Z.H., Khan, M.I., Mahboob, S., Al-Ghanim, A.K. and Al-Misned, F., 2019. Phylogenetic analysis of infectious bursal disease viruses according to newly proposed model of classification into geno-groups. J. Infect. Publ. Hlth., 12: 410–418. https://doi.org/10.1016/j.jiph.2018.12.012

Beenish, Z., Asim, A., Yasin T., Tahir, Y. and Tariq, B., 2016. Conventioal and molecular detection of infectious Bursal disease virus in broiler chicken. Pakistan J. Zool., 48: 601-603.

Beenish, Z., Aslam, A., Anjum, R., Imran, S., Irshad, I., Ara, C. and Sadia, H., 2018. Study on tissue tropism and molecular detection of VP2 gene of infectious bursal disease virus in experimentally infected broiler tissues. Pak. Vet. J., 38: 291-295. https://doi.org/10.29261/pakvetj/2018.054

Dey, S., Pathak, D.C., Ramamurthy, N., Maity, H.K. and Chellappa, M.M., 2019. Infectious bursal disease virus in chickens: prevalence, impact, and management strategies. Vet. Med. Res. Rep., 10: 97. https://doi.org/10.2147/VMRR.S185159

Fan, L., Wu, T., Hussain, A., Gao, Y., Zeng, X., Wang, Yulong, Gao, L., Li, K., Wang, Yongqiang, Liu, C., Cui, H., Pan, Q., Zhang, Y., Liu, Y., He, H., Wang, X. and Qi, X., 2019. Novel variant strains of infectious bursal disease virus isolated in China. Vet. Microbiol., 230: 212–220. https://doi.org/10.1016/j.vetmic.2019.01.023

Gao, L., Li, K., Qi, X., Gao, Y., Wang, Y., Gao, H. and Wang, X., 2018. N-terminal domain of the RNA polymerase of very virulent infectious bursal disease virus contributes to viral replication and virulence. Science China Life Sci., 61: 1127-1129. https://doi.org/10.1007/s11427-017-9297-8

He, X., Xiong, Z., Yang, L., Guan, D., Yang, X. and 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. https://doi.org/10.1007/s00705-014-2195-z

Hon, C.C., Lam, T.Y., Drummond, A., Rambaut, A., Lee, Y.F., Yip, C.W., Zeng, F., Lam, P.Y., Ng, P.T. and Leung, F.C., 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. https://doi.org/10.1128/JVI.00585-06

Hussain, A., Wu, T., Li, H., Fan, L., Li, K., Gao, L., Wang, Y., Gao, Y., Liu, C., Cui, H., Pan, Q., Zhang, Y., Aslam, A., Khan, M.R., Munir, M., Butt, S.L., Wang, X. and Qi, X., 2019. Pathogenic characterization and full length genome sequence of a reassortant infectious bursal disease virus newly isolated in Pakistan. Virol. Sin., 34: 102–105. https://doi.org/10.1007/s12250-019-00082-8

Hussain, A., Wu, T.T., Fan, L.J., Wang, Y.L., Muhammad, F.K., Jiang, N., Gao, L., Li, K., Gao, Y.L., Liu, C.J., Cui, H.Y., Pan, Q., Zhang, Y.P., Aslam, A., Khan, M.R., Arshad, M.I., Abdullah, H.M., Wang, X.M. and Qi, X.L., 2020. The circulation of unique reassortment strains of infectious bursal disease virus in Pakistan. J. Integr. Agric., 19: 1867–1875. https://doi.org/10.1016/S2095-3119(20)63183-5

Islam, M.R., Nooruzzaman, M., Rahman, T., Mumu, T.T., Rahman, M.M., Chowdhury, E.H., Eterradossi, N. and Muller, H., 2021. A unified genotypic classification of infectious bursal disease virus based on both genome segments. Avian Pathol., 50: 190-206. https://doi.org/10.1080/03079457.2021.1873245

Khan, K., Shah, S. and Afzal, M., 1988. Observations on Gumboro disease (infectious bursal disease) in Pakistan. Rev. Sci. Tech. (France). https://doi.org/10.20506/rst.7.3.366

Lone, N.A., Rehmani, S.F., Kazmi, S.U., Muzaffar, R., Khan, T.A., Khan, A., Khan, S.A. and Ahmed, A., 2009. Molecular characterization of Pakistani field isolates of infectious bursal disease virus. Avian Dis., 53: 306-309. https://doi.org/10.1637/8325-042108-Reg.1

Luque, D., Saugar, I., Rejas, M.T., Carrascosa, J.L., Rodriguez, J.F. and Caston, J.R., 2009. Infectious bursal disease virus: ribonucleoprotein complexes of a double-stranded RNA virus. J. mol. Biol., 386: 901. https://doi.org/10.1016/j.jmb.2008.11.029

Maqbool, A., Safdar, M., Zia, N., Mahmood, A., Usman, M., Sharif, M., Khan, A. and Umar, S., 2020. Detection and phylogeny of infectious bursal disease virus (IBDV) during field outbreaks in broilers. Pakistan J. Zool., 52: 659–667. https://doi.org/10.17582/journal.pjz/20190210070218

Mato, T., Tatar-Kis, T., Felfoldi, B., Jansson, D.S., Homonnay, Z., Banyai, K. and Palya, V., 2020. Occurrence and spread of a reassortant very virulent genotype of infectious bursal disease virus with altered VP2 amino acid profile and pathogenicity in some European countries. Vet. Microbiol., 245: 108663. https://doi.org/10.1016/j.vetmic.2020.108663

Muller, H., Islam, M.R. and Raue, R., 2003. Research on infectious bursal disease the past, the present and the future. Vet. Microbiol., 97: 153–165. https://doi.org/10.1016/j.vetmic.2003.08.005

Nwagbo I.O., Shittu I., Nwosuh C.I., Ezeifeka G.O., Odibo F.J., Michel L.O. and Jackwood D.J., 2016. Molecular characterization of field infectious bursal disease virus isolates from Nigeria. Vet. World, 9: 1420. https://doi.org/10.14202/vetworld.2016.1420-1428

Qi, X., Gao, L., Qin, L., Deng, X., Wu, G., Zhang, L., Yu, F., Ren, X., Gao, Y. and Gao, H., 2011. Genomic sequencing and molecular characteristics of a very virulent strain of infectious bursal disease virus isolated in China. Hunan Agric. Sci. Technol., 12: 1946-1949.

Qi, X., Qin, L., Gao, Y., Gao, H., Li, Y., Gao, L., Lu, Z., Wang, N., Chen, Y., Zhang, L. and Li, K., 2015. Genetic analysis of the VP2 hypervariable region of thirty-six infectious bursal disease virus isolates in China during 2009-2012. J. Agric. Sci. Technol., 16: 1565–1569.

Qin, Y. and Zheng, S.J., 2017. Infectious bursal disease virus-host interactions: Multifunctional viral proteins that perform multiple and differing jobs. Int. J. mol. Sci., 18: 161. https://doi.org/10.3390/ijms18010161

Reddy, P.H., Johnson, A.M.A., Kumar, J.K., Naveen, T. and Devi, M.C., 2017. Heterologous expression of infectious bursal disease virus VP2 gene in Chlorella pyrenoidosa as a model system for molecular farming. Pl. Cell Tissue Organ Cult., 131: 119–126. https://doi.org/10.1007/s11240-017-1268-6

Sajid, S., Rahman, S. ur, Gilani, M.M., Sindhu, Z.D., Ali, M. Ben, Hedfi, A., Almalki, M. and Mahmood, S., 2021. Molecular characterization and demographic study on infectious bursal disease virus in Faisalabad district. PLoS One, 16: e0254605. https://doi.org/10.1371/journal.pone.0254605

Shafqat, R., Khan, A., Habib, M., Shah, M.S.U.D., Ali, W., Hussain, Z. and Tahir, Z.A., 2017. Molecular characterization of infectious bursal disease virus from commercial poultry in Pakistan. Matrix Sci. Med., 1: 01–06. https://doi.org/10.26480/msm.02.2017.01.06

Soubies, S.M., Courtillon, C., Briand, F.X., QueguinerLeroux, M., Courtois, D., Amelot, M., Grousson, K., Morillon, P., Herin, J.B. and Eterradossi, N., 2017. Identification of a European interserotypic reassortant strain of infectious bursal disease virus. Avian Pathol., 46: 19–27. https://doi.org/10.1080/03079457.2016.1200010

Stoute, S.T., Jackwood, D.J., Crossley, B.M., Michel, L.O. and Blakey, J.R., 2019. Molecular epidemiology of endemic and very virulent infectious bursal disease virus genogroups in backyard chickens in California, 2009-2017. J. Vet. Diagn. Invest., 31: 371–377. https://doi.org/10.1177/1040638719842193

Tamura, K., Stecher, G., Peterson, D., Filipski, A. and Kumar, S., 2013. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol., 30: 2725–2729. https://doi.org/10.1093/molbev/mst197

Thompson, JD., Gibson, TJ., Plewniak, F., Jeanmougin, F. and Higgins, DG., 1997. The clustal_x windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucl. Acids Res., 25: 4876–4882. https://doi.org/10.1093/nar/25.24.4876

Wang, Y.L., Fan, L.J., Jiang, N., Li, G., Kai, L., Gao, Y.L., Liu, C.J., Cui, H.Y., Qing, P. and Zhang, Y.P., 2021. An improved scheme for infectious bursal disease virus genotype classification based on both genome-segments A and B. J. Integr. Agric., 20: 1372-1381. https://doi.org/10.1016/S2095-3119(20)63424-4

Ye, C., Wang, Y., Zhang, E., Han, X., Yu, Z. and Liu, H., 2018. VP1 and VP3 are required and sufficient for translation initiation of uncapped infectious bursal disease virus genomic double-stranded RNA. J. Virol., 92: e01345-17. https://doi.org/10.1128/JVI.01345-17

Zahid, B., Aslam, A., Tipu, Y., Yaqub, T. and Butt, T., 2016. Conventional and molecular detection of infectious bursal disease virus in broiler chicken. Pakistan J. Zool., 48: 601-603.

Zahid, B., Aslam, A., Anjum, R., Imran, S., Irshad, I., Ara, C. and Sadia, H., 2018. Study on tissue tropism and molecular detection of VP2 gene of infectious bursal disease virus in experimentally infected broiler tissues. Pak. Vet. J., 38: 291-295. https://doi.org/10.29261/pakvetj/2018.054

Zahoor, M., Abubakar, M., Naim, S., Khan, Q. and Arshed, M., 2011. Molecular typing of field isolates from two outbreaks of infectious bursal disease virus from Pakistan. Vet. World, 4: 297. https://doi.org/10.5455/vetworld.4.297