A
Unraveling VP2 mutations in Infectious Bursal Disease Virus (IBDV) associated with potential vaccine escape in poultry flocks in Shandong, China
BoLiu1,2
YanliBi1,3
JinwenXie1,3
HaigeSu1,4
LisaneworkEAyalew5
QihuiLuo6
LizhongMiao1,3
NaTang1,3
WubshetAshenafiKiros1,2
AbdelrahmanSaid7
WenxiuWang1,2,3✉Email
1Shandong Binzhou Institute of Animal Husbandry and Veterinary Science256600BinzhouChina
2Shandong Academician Workstation256600BinzhouShandongChina
3Binzhou Key Laboratory of Animal and Human Diseases256600BinzhouShandongChina
4Henan Province Key Laboratory for Animal Food Pathogens Surveillance, College of Veterinary MedicineHenan Agricultural University450046ZhengzhouChina
5
A
Department of Pathology and Microbiology, Atlantic Veterinary CollegeUniversity of Prince Edward Island550 University AveC1A 4P3CharlottetownPE
6Laboratory of Experimental Animal Disease Model, College of Veterinary MedicineSichuan Agricultural University611130ChengduChina
7
A
Genetics and Genome Biology ProgramThe Hospital for Sick ChildrenTorontoON
8Canada and Parasitology and Animal Diseases DepartmentNational Research CenterDokkiGizaEgypt
Bo Liu1,2†, Yanli Bi1,3†, Jinwen Xie1,3, Haige Su1,4, Lisanework E Ayalew5, Qihui Luo6, Lizhong Miao1,3, Na Tang1,3, Wubshet Ashenafi Kiros1,2, Abdelrahman Said7, Wenxiu Wang1,2,3*
1 Shandong Binzhou Institute of Animal Husbandry and Veterinary Science, Binzhou 256600, China.
2 Shandong Academician Workstation, Binzhou 256600, Shandong, China.
3 Binzhou Key Laboratory of Animal and Human Diseases, Binzhou 256600, Shandong, China.
4 Henan Province Key Laboratory for Animal Food Pathogens Surveillance, College of Veterinary Medicine, Henan Agricultural University, Zhengzhou, 450046, China.
5 Department of Pathology and Microbiology, Atlantic Veterinary College, University of Prince Edward Island, 550 University Ave, Charlottetown, PE C1A 4P3.
6 Laboratory of Experimental Animal Disease Model, College of Veterinary Medicine, Sichuan Agricultural University, Chengdu 611130, China.
7 Genetics and Genome Biology Program, The Hospital for Sick Children, Toronto, ON, Canada and Parasitology and Animal Diseases Department, National Research Center, Dokki, Giza, Egypt.
7
*Correspondence: wwx1997@126.com
Bo Liu and Yanli Bi contributed equally to this work.
Abstract
Background
Infectious bursal disease virus (IBDV) remains a major threat to the global poultry industry, causing significant economic losses. Recent surveillance data indicate that genetically diverse strains of IBDV are circulating in different regions of China.
Methods
In this study, both commercial and backyard chicken flocks from Shandong Province were investigated for evidence of IBDV circulation using a combination of serological assays and molecular detection techniques.
Results
A
Despite routine vaccination, antibody responses were highly variable, with some flocks showing reduced maternal immunity while others exhibiting active seroconversion by 35 days of age. Many birds showed an abnormally low bursa-to-body weight ratio, indicating immunosuppression. RT-PCR detected IBDV in 68 bursal samples, most of which were classified as very virulent strains (vvIBDV) belonging to genogroup A3, while a smaller subset was identified as variant strains within genogroup A2 based on phylogenetic analysis of the VP2 gene. Comparison of the VP2 protein sequence with other Chinese strains revealed variations in key hypervariable region residues critical for antibody binding. Conserved markers of very virulent strains persisted, while the novel substitutions A222V and D279H emerged, neither of which is present in commonly used vaccine strains.
Conclusions
Our findings highlight emerging IBDV variants with altered VP2 antigenicity, potentially causing vaccine escape and increasing poultry health and economic risks.
Keywords:
Infectious bursal disease virus (IBDV)
Very virulent strains (vvIBDV)
VP2 hypervariable region (VP2_HVR)
Phylogenetic analysis
Vaccine escape
Background
Infectious bursal disease virus (IBDV), the causative agent of infectious bursal disease (IBD), is a highly contagious pathogen that primarily affects young chickens between three and six weeks of age. The disease primarily targets the bursa of Fabricius and is highly immunosuppressive, increasing susceptibility of affected birds to secondary infections [1, 2]. The disease remains to be a problem to poultry health in China and undergone significant evolution since its detection in the country [2]. Although vaccination has played a critical role in controlling the disease, evolutionary pressures from both the environment and vaccination have driven the emergence and spread of novel variant strains. These variants now predominate in China and present significant challenges to current vaccination efforts [35].
IBDV belongs to the genus Avibirnavirus within the family Birnaviridae and possesses a bi-segmented double-stranded RNA genome [6]. Segment A encodes a polyprotein that undergoes co- and post-translational cleavage to yield VP2, VP3, VP4, and a small non-structural protein VP5, while segment B encodes VP1, the viral RNA-dependent RNA polymerase [6]. Among these, VP2 forms the major structural and antigenic protein on the viral capsid surface [79]. The hypervariable region of VP2 (VP2_HVR), spanning amino acids 206 to 350, is pivotal in determining the antigenic type and virulence of IBDV strains [10]. Mutations, insertions, or deletions within this region drive antigenic drift and facilitate the emergence of immune escape variants [1113].
Since its initial emergence in the early 1960s, two distinct serotypes of IBDV have been identified. Serotype-I is pathogenic to chickens, while serotype-II infects other avian species but is non-pathogenic to chickens [7]. Within serotype-I, multiple antigenic and pathogenic variants exist, including classical, variant, very virulent, attenuated, and recombinant strains [8]. Particularly, very virulent IBDV strains have caused widespread and severe outbreaks across Europe, Asia, and Africa, becoming predominant in China since the 1990s [14]. Simultaneously, novel antigenic variants distinct from strains circulating in the United States have emerged in Chinese poultry farms [15, 16]. These variants exhibit altered antigenicity, reduced susceptibility to neutralization by antibodies induced by classical or vaccine strains, and enhanced environmental persistence, complicating disease control efforts [1618].
Vaccination remains the primary tool for controlling the disease, with both live attenuated and inactivated vaccines widely implemented [19]. However, the emergence and spread of immune-escape variants and recombinant viruses has made disease control and prevention programs increasingly challenging [14, 2022]. These variants undermine the protective efficacy of current vaccines, particularly in regions where multiple IBDV genotypes co-circulate [23]. Genetic reassortment and recombination between vaccine and field strains, along with co-infections, contribute to the virus’s extensive genetic diversity and rapid evolution of IBDV in China [13, 24, 25]. These dynamics foster the periodic emergence of highly virulent or antigenically novel strains capable of evading host immunity. In the past couple of years, there have been increasing reports of very virulent and novel immune-escape variant strains of IBDV from different provinces in China [17]. Despite widespread vaccination, Shandong Province, one of China’s largest poultry-producing regions, has reported a growing number of IBD outbreaks affecting both commercial and backyard poultry flocks [14]. This rise suggests possible antigenic mismatch and the circulation of novel immune-escape variants. Moreover, backyard flocks may serve as potential reservoirs facilitating viral transmission to commercial operations. However, genetic and serological data on circulating field strains in this region remain limited, particularly regarding the VP2_HVR, a critical determinant of antigenicity and pathogenicity. To address these gaps, the present study aimed to detect and molecularly characterize IBDV strains from commercial and backyard poultry farms in Binzhou, Shandong Province. By analyzing the VP2 region and assessing serological responses alongside immunosuppression indicators such as bursa-to-body weight ratios, we aim to better understand the genetic diversity and antigenic variation of local strains and contribute to improved region-specific IBD control strategies.
Methods
Sample collection
A
A total of 306 broiler chickens were sampled between November 2019 and June 2021 from nine commercial poultry farms (Farm: A, B, C, D, E, F, GR, NY, and SG) as well as from backyard flocks in Binzhou and adjacent regions of Shandong Province, China. All commercial farms adhered to vaccination programs against IBDV and other major avian viral pathogens (Supplementary Table S1). On six farms (A-F), 20 broilers were randomly selected at 18 days of age and another 20 at 35 days of age. Body weight and bursa of Fabricius weight were recorded, and blood samples were collected prior to euthanasia. Bursa of Fabricius tissues were then harvested for nucleic acid extraction. Additional bursal samples submitted from Farms GR, NY, SG, and backyard poultry flocks were also screened for IBDV by PCR. Blood samples could not be collected from backyard flocks and from Farms GR, NY, and SG due to logistical limitations and the nature of sample submission from external sources, which precluded on-site handling of live birds.
A
A
All procedures were conducted in accordance with the ethical guidelines of the Animal Ethics Committee of the Shandong Binzhou Animal Science and Veterinary Medicine Academy (protocol number SYXK [LU] 20210008).
Enzyme-Linked Immunosorbent Assay (ELISA)
Serum samples were collected from clinically healthy broiler chickens at the six designated farms (Farm A-F). Samples were randomly obtained at 18 days of age (n = 20) and 35 days of age (n = 20) from each farm. The samples were centrifuged at 3,000 rpm for 10 minutes to separate the serum, which was then aliquoted and stored at − 20°C until further analysis. Prior to testing, serum samples were diluted 1:2 and analyzed using an indirect ELISA with the IDEXX IBD Ab Test kit (# 99-09260, IDEXX Laboratories, Inc., Westbrook, ME, USA), following the manufacturer’s instructions. Absorbance was measured using a microplate reader, and antibody titers were calculated accordingly. A cut-off value of 396 was applied to determine seropositivity.
RNA extraction
Approximately 100 mg of each tissue sample was homogenized using a sterile mortar and pestle in 1 mL of phosphate-buffered saline (PBS, pH 7.4) containing penicillin and streptomycin at a final concentration of 700 µg/mL. Homogenization was carried out on ice to minimize RNA degradation. The homogenates were then centrifuged at 8000 rpm for 5 minutes at 4°C, and 200 µL of the clarified supernatant was collected for RNA extraction. Total viral RNA was extracted using the Simply P Virus DNA/RNA Co-Extraction Kit (Cat. No. BSC67, Bioer Technology Co., Hangzhou, China) according to the manufacturer’s protocol. All procedures were performed under RNase-free conditions to maintain RNA integrity. The quality and concentration of extracted RNA were evaluated using both agarose gel electrophoresis and spectrophotometry. High-quality RNA samples were stored at − 80°C until further analysis.
Reverse transcription and PCR amplification
Complementary DNA (cDNA) synthesis was performed using random hexamer primers and the SuperScript™ Reverse Transcriptase Kit (TransGen Biotech Co., Ltd., Beijing, China). The reverse transcription reaction was performed in a 20 µL reaction volume following the manufacturer’s instructions: RNA and primers were initially denatured at 65°C for 5 minutes, followed by reverse transcription at 42°C for 50 minutes, and final inactivation at 70°C for 15 minutes. The synthesized cDNA was stored at − 20°C until further use. To amplify the HVR of the VP2 gene of IBDV, a pair of specific primers was employed: IBDP1 (forward: 5′-TCACCGTCCTCAGCTTAC-3′) and IBDP2 (reverse: 5′-TCAGGATTTGGGATCAGC-3′), targeting a 643 bp fragment. The PCR was conducted in a 25 µL reaction mixture containing 2 µL of cDNA, 0.4 µM of each primer, 200 µM dNTPs, 2.5 mM MgCl₂, 1× PCR buffer, and 1 U of Taq DNA polymerase (TransTaq, TransGen Biotech, China). Thermal cycling conditions were: 94°C for 5 minutes (initial denaturation); 35 cycles of 94°C for 30 seconds, 55°C for 30 seconds, and 72°C for 45 seconds; with a final extension at 72°C for 10 minutes. PCR products were resolved by electrophoresis on 1% agarose gels stained with ethidium bromide (0.5 µg/mL) and visualized using a UV transilluminator. A 100 bp DNA ladder was used as a molecular weight marker to confirm the expected product size.
Immunohistochemistry (IHC) assay
Based on the availability of tissues after conducting the PCR assay, the IHC assay was performed on 17 tissue samples, as listed in Supplementary Table S3. Briefly, Bursa of Fabricius samples were excised into 1–2 cm³ sections, rinsed with physiological saline, and blotted dry using filter paper. Tissues were fixed in 10% neutral buffered formalin, dehydrated through a graded ethanol series, embedded in paraffin, and sectioned at a thickness of 0.6 µm. Sections were stained with hematoxylin and eosin (H&E), cleared with xylene, and mounted with neutral resin. Histological slides were examined and scanned using a microscope (Olympus VS200, Olympus Corporation), and the resulting images were processed and analyzed with Olyvia software.
Sanger sequencing
PCR products showing a single specific band of the expected 643 bp size were excised from the agarose gel and purified using a commercial gel extraction kit (E.Z.N.A.® Gel Extraction Kit, Omega Bio-tek, USA). Purified DNA samples were quantified and submitted for Sanger sequencing at Sangon Biotech Co. Shanghai, China. Sequencing was performed in both forward and reverse directions using the same primers (IBDP1 and IBDP2) used for PCR amplification to ensure high-fidelity reads. The resulting sequence data were quality-checked, assembled, and aligned using BioEdit and MEGA software. Sequences were then used for downstream phylogenetic analysis, including alignment with reference strains and construction of phylogenetic trees to assess genetic relationships among circulating IBDV strains.
Phylogenetic Tree and Protein Sequence Analysis
The HVR of the VP2 gene was sequenced from sixty-eight randomly selected IBDV-positive bursal samples with characteristics summarized in Supplementary Table S1, including reference number, bird type, age, farm source, and vaccination status, as well as from the VIRGO7 vaccine strain (Vigoly Bio. Tech Co. Ltd., Beijing, China). Quality control, sequence assembly, and analysis were conducted using Geneious Prime software (Dotmatics). The obtained viral sequences were aligned and analyzed alongside 50-VP2 gene reference sequences retrieved from GenBank, representing various genotypes and phenotypes of IBDV reported in China and other regions globally. Multiple sequence alignment was performed using the Clustal Omega program. A phylogenetic tree was constructed, and evolutionary relationships were inferred using the maximum likelihood method with the Jukes-Cantor (JC69) genetic distance model, implemented in the PhyML program within Geneious Prime. The phylogenetic tree was visualized and annotated using the Interactive Tree of Life (iTOL) v5 online tool [26]. Predicted amino acid sequences of the VP2 proteins were aligned, and variability at each residue was analyzed to determine viral phenotypes, following previously described methods [27].
Statistical Analysis
Data were analyzed using GraphPad Prism version 10 (Dotmatics). A two-tailed unpaired t-test was conducted to compare IBDV antibody levels and bursa-to-body weight ratios at 18 and 35 days of age within each farm. p-values < 0.05 were considered statistically significant.
Results
IBDV antibody dynamics reflect infection and immunity profiles in broiler chickens
The levels of IBDV-specific antibodies were measured in broiler chickens from six farms at 18 and 35 days of age to assess pre-existing immunity or recent viral exposure (Fig. 1). At 18 days, birds from Farms A (426.45 ± 477.84) and F (644.1 ± 493.5) exhibited significantly higher antibody titers than their 35-day-old counterparts (Farm A: 303.25 ± 318.12; Farm F: 79.9 ± 92.22; p = 0.004 and 0.001, respectively), with titers at 35 days falling below the positive cut-off value of 396, suggesting waning immunity or viral clearance. Conversely, birds from Farms B, C, D, and E showed a significant rise in antibody titers between 18 and 35 days (p < 0.0001). Titers increased from below the cut-off at 18 days (Farm B: 77 ± 72.5; Farm C: 67.3 ± 78.29; Farm D: 54.55 ± 88; Farm E: 138.2 ± 215.95) to markedly elevated levels by 35 days (Farm B: 1941.45 ± 1150.12; Farm C: 2192.3 ± 826.9; Farm D: 2584.25 ± 1266.97; Farm E: 2340.35 ± 1470.9), confirming active seroconversion following recent IBDV infection.
Fig. 1
IBDV-specific antibody titers in broiler chickens from six farms: IBDV-specific antibody levels were quantified by ELISA at 18 and 35 days of age across Farms A to F. At 18 days, birds from Farms A and F exhibited significantly higher antibody titers than at 35 days, with titers declining below the positive cut-off value of 396 at 35 days, suggesting waning immunity or viral clearance. In contrast, Farms B, C, D, and E showed significant increases in antibody titers between 18 and 35 days (****p < 0.0001), reflecting active seroconversion and recent IBDV infection. Data are presented as mean ± standard deviation.
Click here to Correct
Bursal atrophy highlights IBDV-induced immunosuppression in broiler chickens
To evaluate the immunosuppressive effects of IBDV infection, bursa-to-body weight (BBW) ratio was measured in sampled birds (Fig. 2). At 18 days, birds from Farms A-E showed significantly higher BBW ratios compared to those at 35 days (p < 0.0001). Mean values declined from 0.173 ± 0.048 to 0.092 ± 0.026 in Farm A, 0.198 ± 0.065 to 0.060 ± 0.013 in Farm B, 0.166 ± 0.046 to 0.081 ± 0.064 in Farm C, 0.157 ± 0.034 to 0.055 ± 0.017 in Farm D, and 0.186 ± 0.039 to 0.062 ± 0.026 in Farm E, indicating bursal atrophy characteristic of IBDV-induced immunosuppression. In contrast, Farm F showed no significant difference in BBW ratio between 18 days (0.293 ± 0.083) and 35 days (0.234 ± 0.145; p = 0.1232), suggesting stable bursal development and absence of active IBDV-related bursal damage during the sampling period.
Fig. 2
Bursa-to-body weight ratios in broiler chickens from six farms: Bursa-to-body weight (BBW) ratios were measured at 18 and 35 days of age to assess bursal atrophy associated with IBDV infection. Significant reductions in BBW ratios were observed in Farms A, B, C, D, and E between 18 and 35 days (****p < 0.0001), consistent with immunosuppression due to bursal damage. Conversely, Farm F showed no significant change in BBW ratio (p = 0.1232), indicating stable bursal development and absence of active IBDV-related damage during the study period. Data are expressed as mean ± standard deviation.
Click here to Correct
Phylogenetic characterization of VP2 reveals high genetic diversity among vvIBDV isolates circulating in in broiler chickens
Phylogenetic analysis was conducted on sixty-eight newly generated VP2 gene sequences, along with 50 reference IBDV_VP2 sequences retrieved from GenBank, representing established genotypes and phenotypes, including recently reported Chinese isolates. The GenBank accession numbers of the VP2 sequences are provided in Supplementary Table S2. All 68 PCR-positive samples, randomly selected for sequencing, confirmed to be infected with IBDV by both PCR (data not shown) and IHC assay (Supplementary Table S3). Of the 68-sequences analyzed, 66 clustered within genogroup A3 corresponding to vvIBDV strains, further subdividing into three distinct sub-lineages (Fig. 3). The remaining two sequences grouped with genogroup A2, indicative of variant phenotypes. None of the newly sequenced isolates clustered within genogroups A1 or A4–A7. Farm-specific VP2 sequence homology ranged from 97.88–100% (Farm A), 95.72–98.91% (Farm B), 98.58–99.66% (Farm C), 99.29–100% (Farm D), 96.99–100% (Farm E), and 98.91–100% (Farm F), suggesting co-circulation of highly related vvIBDV strains at individual sites. Several isolates (notably B12, B16 from Farm B; C2 from Farm C; and 20210604-6 from Farm SG) exhibited lower sequence identity (87.37–96.85%) relative to the main A3 cluster, suggesting the presence of more divergent vvIBDV lineages. Strains from Farm GR (IBDV-GR1, IBDV-GR2–GR7), Farm NY (IBDV-NY-1 to IBDV-NY-10), and Farm SG (20210604-6, 20210604-12), as well as all four backyard chicken isolates (IBDV-ZWT-1 to -4), also grouped within genogroup A3 and shared 91.2–100% identity, indicating viral exchange between production systems. Analysis of commercial live vaccine strains revealed distinct phylogenetic clustering: the Virgo 7 intermediate hot strain vaccine (Vigoly Bio. Tech Co., Beijing, China) fell within genogroup A3, while the B87 attenuated vaccine (Shandong Dorivo Biological Engineering Co., Ltd., Shandong, China) belonged to genogroup A8. Among the 66-genogroup A3 isolates, nucleotide identity was predominantly 96.22–100%, aside from the divergent strains noted above. Comparisons with recently reported Chinese vvIBDV strains HB20 and HLJ19 showed sequence identities of 89.66–96.55% and 90.16–97.67%, respectively. Homology with the HLJ0504-194-like vvIBDV strain (GenBank GQ451330), a major cause of IBD epidemics in China, ranged from 91.05–97.94%. Virgo 7 and B87 vaccine strains exhibited sequence identities of 87.57–98.63% and 87.37–95.53%, respectively, with the field isolates from this study. The two variant strains (A2) shared 96.2% identity with one another, 90.34–97.17% identity with reference variant strains, and 87.57–92.46% with vaccine strains. Overall, these findings confirm that vvIBDV strains of genogroup A3 are widely distributed in commercial, local, and backyard poultry in Binzhou, displaying considerable genetic diversity and multiple sub-lineages, while the low prevalence of variant strains suggests limited current circulation.
Fig. 3
Phylogenetic tree of IBDV based on VP2_HVR sequences: The tree was constructed using the maximum likelihood method with 68 newly sequenced VP2_HVR genes and 50 reference sequences from GenBank. Sixty-six field isolates clustered within genogroup A3, including four from backyard chickens (IBDV-ZWT-1 to -4), and two isolates grouped under A2. The Virgo7 and B87 vaccine strains clustered within A3 and A8, respectively. Bootstrap values > 70% are shown at major nodes. The scale bar represents the number of nucleotide substitutions per site.
Click here to Correct
Phenotypic characterization of VP2_HRV based on amino acid residues
The phenotypic analysis of the VP2_HRV protein sequences revealed a high level of amino acid conservation among most IBDV isolates classified within genogroup A3, with sequence identities ranging from approximately 91.2% to 100%. However, three field isolates C2, B12, and B16 exhibited notably lower identity values, ranging from 81.3% to 92.2%, when compared to other A3 isolates. Compared to the highly virulent Chinese HLJ0504 strain, these three isolates shared only 82.9% to 89.6% sequence identity. IBDV strains found in backyard poultry showed high genetic similarity, mostly over 97%, to those in commercial poultry, indicating that the virus maybe transmitted between these different farming systems. As illustrated in Fig. 4, eight A3 isolates, including two from backyard chickens, contained conserved amino acids at critical residues within the VP2_HVR that define very virulent IBDV. Specifically, these residues included alanine at position 222, isoleucine at 242, glutamine at 253, isoleucine at 256, aspartic acid at 279, alanine at 284, isoleucine at 294, and serine at 299. This amino acid profile mirrors that of the HLJ0504 strain. Among these, three isolates exhibited unique mutations such as alanine to valine at position 222 and aspartic acid to histidine at position 279, with the valine substitution representing a novel variation not previously reported in Chinese field isolates. The remaining 57 A3 isolates retained classical very virulent IBDV features but consistently carried this valine substitution. Two isolates classified within genogroup A2 exhibited the characteristic amino acid motif of variant IBDVs, including lysine at 221, isoleucine at 252, and serine at 299. The Virgo7 vaccine strain showed a hybrid profile combining very virulent markers with attenuating residues such as threonine at 272 and asparagine at 279, whereas the B87 vaccine strain presented classical attenuation-associated residues including proline at 222, asparagine at 279, threonine at 284, and arginine at 330. Notably, none of the IBDV strains isolated from backyard poultry harbored classical or vaccine-associated amino acid residues, supporting their classification as field-derived very virulent strains. These findings collectively confirm that the majority of circulating IBDV strains in both commercial and backyard poultry retain conserved very virulent amino acid signatures within the VP2_HVR, while the presence of unique residue substitutions in some isolates highlights early viral divergence and underscores the importance of continuous molecular surveillance to ensure effective vaccine coverage and disease control.
Fig. 4
Amino acid residue substitutions in the hypervariable region of the VP2 protein: Highlighted in light blue are amino acid residues and the characteristic phenotypic features they confer to the virus. Highlighted in light grey represents our isolates and their predicted phenotypic characteristics. Unique amino acid residues observed in our isolates are shown in red.
Click here to Correct
Discussion
This study provides valuable insights into the dynamics of maternal immunity and the serological response to IBDV in poultry farms in Binzhou, Shandong Province, China. Our data reveal two distinct antibody response profiles across the farms. In Farms A and F, birds exhibited relatively high IBDV-specific antibody titers at 18 days of age, suggesting the presence of maternally derived antibodies. However, these titers declined significantly by 35 days, falling below the protective threshold. This rapid decline likely indicates suboptimal maternal immunity transfer from breeder flocks, possibly due to factors such as improper vaccine storage, low-quality vaccines, incorrect administration, or non-adherence to vaccination schedules [21, 28, 29]. Effective maternal antibody transfer is essential for early protection against IBDV, and inadequate levels are associated with increased susceptibility to infection in progeny [30]. Of note, Farm A showed a significant increase in PCR positivity at 35 days, highlighting the susceptibility of birds following the decay of maternal antibodies In contrast, Farm F showed early PCR positivity at 18 days, but both antibody titers and bursal integrity declined by 35 days, suggesting initial protection followed by waning immunity. These observations are consistent with prior studies indicating the protective role of maternal antibodies against IBDV during early life [3133]. Conversely, Farms B, C, D, and E displayed low or undetectable antibody levels at 18 days, suggesting that maternal immunity was either absent or insufficient. These findings could reflect a lack of vaccination or insufficient immunity in breeder flocks, which may have been free of prior infection or had implemented effective biosecurity protocols. However, the significant increase in antibody titers by 35 days suggests that these birds were naturally exposed to IBDV, leading to an active immune response. This natural seroconversion was corroborated by a marked reduction in bursa-to-body weight ratios, confirming IBDV-induced bursal atrophy and immunosuppression. Although live IBDV vaccines have been shown to cause some degree of bursal damage, the amino acid signatures of circulating strains in our study were distinct from those of vaccine strains, further supporting the conclusion that field strain exposure caused the observed pathology.
Molecular analysis of the VP2_HVR revealed that most isolates belonged to genogroup A3, corresponding to the vvIBDV phenotype, while a smaller proportion clustered within genogroup A2, representing variant strains. The high sequence identity within genogroup A3 suggests close genetic relatedness; however, notable intraspecies variability and the presence of distinct sub-lineages highlight ongoing viral evolution, consistent with reports from other regions in China and globally [32, 34, 35]. Importantly, our phylogenetic data showed evidence of viral exchange between commercial and backyard poultry, with several isolates sharing > 97% identity across production systems. This suggests cross-transmission between different poultry sectors in the region. The continuous evolution of field strains is likely driven by immune selection pressure from vaccination and natural infection [3, 36, 37]. The co-circulation of multiple IBDV strains in both commercial and backyard poultry, including genetically identical strains, suggests cross-transmission between production systems. This finding underscores the need for integrated control strategies and coordinated surveillance across the poultry sector to prevent the spread of IBDV [38, 39]. The VP2 protein plays a pivotal role in IBDV virulence and antigenicity, and the mutations identified in this study provide important insights into viral evolution and immune escape. Most A3 isolates retained key amino acid residues at positions 222, 253, 279, and 284, which are well-established markers of vvIBDV strains and are critical for pathogenicity and antigenic recognition [40, 41]. However, the detection of novel mutations, particularly the A222V substitution, is noteworthy. Although this mutation has not been previously reported in Chinese strains, it has been documented in Indian isolates, is associated with immune escape mechanisms, and reduced binding affinity of neutralizing antibodies, thereby potentially undermining vaccine efficacy [15, 42]. The presence of A222V in circulating strains suggests ongoing antigenic drift and the emergence of new variants that could contribute to vaccine failure. Additionally, the identification of a D279H substitution in several isolates raises questions regarding its functional impact. While substitutions at position 279, including D279N, contribute to attenuation of virulence, the impact of D279H remains uncharacterized and requires in vivo studies to determine its role in viral pathogenicity and immune evasion. [34, 43]. Mutations at position 253, including Q253H, have also been implicated in enhanced virulence and resistance to neutralization by vaccine-induced antibodies, reinforcing the notion that these changes facilitate viral adaptation and immune escape [43]. Collectively, these findings underscore the dynamic nature of IBDV evolution in the field and highlight the challenges posed to current vaccination strategies, emphasizing the need for continuous molecular surveillance and vaccine updates to maintain effective disease control.
Our study also highlights the significant divergence between circulating field strains and commercial vaccine strains such as B87 and Virgo7, both of which possess attenuation-associated residues at key VP2 positions that are absent in the field isolates. Furthermore, the consistent identification of vvIBDV residues in backyard poultry isolates, without evidence of vaccine-associated amino acid markers, confirms that these strains are field-derived and not vaccine escape variants. This divergence underscores the potential limitations of current vaccines in providing effective protection against circulating vvIBDV strains. Experimental and field studies have demonstrated that vaccines including viral vector, immune-complex, and recombinant subunit formulations often exhibit reduced protective efficacy when challenged with contemporary vvIBDV field strains [44]. The genetic and phenotypic differences identified in our study further corroborate these findings and highlight the urgent need for vaccine reformulation to address the ongoing evolution of viral populations. Given the observed antigenic divergence, vaccination strategies require updating to include antigens more relevant to regional variants. This approach would enhance concordance with circulating field strains and may improve vaccine efficacy. Furthermore, optimizing breeder vaccination protocols is essential to ensure effective maternal antibody transfer to offspring. Ongoing molecular surveillance and regular updates to vaccine formulations are critical to prevent vaccine failure and to achieve sustained control of IBDV in poultry production [45, 46].
Conclusions
Our study demonstrated the emergence and predominant circulation of unique vvIBDV variants clustered in genogroup A3 in Binzhou, Shandong, China. These viruses primarily contain amino acid substitutions at residue positions 222 and/or 279 in the HVR of the VP2 protein, which are different from the vaccine strains. The novel A222V substitution, consistently observed in most A3 isolates, may represent an emerging molecular signature in Chinese vvIBDV populations. The observed unique mutations can affect virulence and modify the confirmation of dominant epitopes in the P domain of VP2 protein changing the antigenic properties of the emerging vvIBDV. The viruses could break through existing immunity and cause disease. Strict biosecurity measures are essential to restrict the spread of emerging viruses. In addition, continuous surveillance and animal challenge experimental studies are required to evaluate the pathogenicity and antigenicity of these viruses to select for candidates for the development of antigenically representative vaccines to circulating strains.
List of abbreviations
IBDV
Infectious bursal disease virus
IBD
Infectious bursal disease
HRV
Hypervariable region
vvIBDV
Very virulent strains Infectious bursal disease virus
VP2_HVR
VP2 hypervariable region
iTOL
Interactive Tree of Life
Figures legends
Supplementary Materials:
A
Table S1
Characteristics of sampled birds, including reference number, bird type, age, farm source, and vaccination status.
A
Table S2
GenBank accession numbers of the new VP2 sequences.
A
Table S3
Immunohistochemistry analysis of the Bursa of Fabricius
Ethics approval and consent to participate:
A
The sample collection protocol was approved by the Animal Ethics Committee of the Shandong Binzhou Institute of Animal Husbandry and Veterinary Sciences (approval number: SYXK [LU] 20210008).
Consent for publication:
N/A
Availability of data and materials:
The data will be freely available to anyone upon request. The VP2 gene sequences are deposited in the GenBank and accession numbers are provided in Supplementary Table 2.
Competing Interests:
The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.
A
Funding:
This research was funded by the Shandong Provincial Double Hundred Plan (grant numbers WST2018014 and WSR2023042), the Shandong Provincial Poultry Industry and Technology System (grant number SDAIT-11-16), and the National Key Research and Development Program of China (grant number 2023YFD1802602-3).
A
Author Contribution
Conceptualization, W.W.; methodology, Y.B., J.X., N.T., B.L.; software, B.L., Y.B., H.S., J.X., WAK and L.E.A.; validation, N.T., W.W., and A.S.; formal analysis, Y.B., and B.L.; investigation, B.L., Y.B., H.S., L.M., and J.X.,; resources, W.W.; data curation, B.L., Q.L., L.M., and Y.B., writing original draft preparation, B.L. and Y.B.; writing, review and editing, L.E.A., Q.L., A.S., and W.W.; visualization, B.L., Y.B., A.S., and L.E.A.; supervision, W.W.; project administration, W.W.; funding acquisition, W.W. All authors have read and agreed to the published version of the manuscript.
Acknowledgments:
N/A
Electronic Supplementary Material
Below is the link to the electronic supplementary material
References:
1.
Sharma JM, Kim IJ, Rautenschlein S, Yeh HY. Infectious bursal disease virus of chickens: pathogenesis and immunosuppression. Dev Comp Immunol. 2000;24:223–35.
2.
Zhang T, Wang S, Liu Y, Qi X, Gao Y. Advances on adaptive immune responses affected by infectious bursal disease virus in chicken. Front Immunol. 2023;14:1330576.
3.
Zhu K, Wu Q, Leng M, Wang Z, Lin W. Phylogenetic analysis of vp2 gene of the infectious bursal disease virus in South China during 2023. Front Vet Sci. 2025;12:1575407.
4.
Lian J, Wang Z, Xu Z, Pang Y, Leng M, Tang S, Zhang X, Qin J, Chen F, Lin W. Pathogenicity and molecular characterization of infectious bursal disease virus in China. Poult Sci. 2022;101:101502.
5.
Pikula A, Lisowska A, Jasik A, Smietanka K. Identification and assessment of virulence of a natural reassortant of infectious bursal disease virus. Vet Res. 2018;49:89.
6.
Dobos P, Hill BJ, Hallett R, Kells DT, Becht H, Teninges D. Biophysical and biochemical characterization of five animal viruses with bisegmented double-stranded RNA genomes. J Virol. 1979;32:593–605.
7.
Schnitzler D, Bernstein F, Muller H, Becht H. The genetic basis for the antigenicity of the VP2 protein of the infectious bursal disease virus. J Gen Virol. 1993;74(Pt 8):1563–71.
8.
Brandt M, Yao K, Liu M, Heckert RA, Vakharia VN. Molecular determinants of virulence, cell tropism, and pathogenic phenotype of infectious bursal disease virus. J Virol. 2001;75:11974–82.
9.
van Loon A, de Haas N, Zeyda I, Mundt E. Alteration of amino acids in VP2 of very virulent infectious bursal disease virus results in tissue culture adaptation and attenuation in chickens. J Gen Virol. 2002;83:121–9.
10.
Coulibaly F, Chevalier C, Gutsche I, Pous J, Navaza J, Bressanelli S, Delmas B, Rey FA. The birnavirus crystal structure reveals structural relationships among icosahedral viruses. Cell. 2005;120:761–72.
11.
Jackwood DJ, Sommer-Wagner S. Genetic characteristics of infectious bursal disease viruses from four continents. Virology. 2007;365:369–75.
12.
Jackwood DJ, Sreedevi B, LeFever LJ, Sommer-Wagner SE. Studies on naturally occurring infectious bursal disease viruses suggest that a single amino acid substitution at position 253 in VP2 increases pathogenicity. Virology. 2008;377:110–6.
13.
Aricibasi M, Jung A, Heller ED, Rautenschlein S. Differences in genetic background influence the induction of innate and acquired immune responses in chickens depending on the virulence of the infecting infectious bursal disease virus (IBDV) strain. Vet Immunol Immunopathol. 2010;135:79–92.
14.
Zhang W, Wang X, Gao Y, Qi X. The Over-40-Years-Epidemic of Infectious Bursal Disease Virus in China. Viruses 2022, 14.
15.
Fan L, Wu T, Hussain A, Gao Y, Zeng X, Wang Y, Gao L, Li K, Wang Y, Liu C, et al. Novel variant strains of infectious bursal disease virus isolated in China. Vet Microbiol. 2019;230:212–20.
16.
Bo Z, Wang S, Xu K, Zhang C, Guo M, Cao Y, Zhang X, Wu Y. Isolation and genomic characteristics of the novel variant infectious bursal disease virus in China. Front Vet Sci. 2023;10:1314903.
17.
Huang Y, Shu G, Huang C, Han J, Li J, Chen H, Chen Z. Characterization and pathogenicity of a novel variant infectious bursal disease virus in China. Front Microbiol. 2022;13:1039259.
18.
Aliyu HB, Hamisu TM, Hair-Bejo M, Omar AR, Ideris A. Immuno-protective effect of commercial IBD vaccines against emerging novel variant infectious bursal disease virus in specific-pathogen-free chickens. J Vet Sci. 2024;25:e70.
19.
Jackwood DJ. Advances in vaccine research against economically important viral diseases of food animals: Infectious bursal disease virus. Vet Microbiol. 2017;206:121–5.
20.
Islam MR, Nooruzzaman M, Rahman T, Mumu TT, Rahman MM, Chowdhury EH, Eterradossi N, Muller H. A unified genotypic classification of infectious bursal disease virus based on both genome segments. Avian Pathol. 2021;50:190–206.
21.
Muller H, Islam MR, Raue R. Research on infectious bursal disease–the past, the present and the future. Vet Microbiol. 2003;97:153–65.
22.
Pikula A, Lisowska A, Domanska-Blicharz K. Epidemiology of Infectious Bursal Disease Virus in Poland during 2016–2022. Viruses 2023, 15.
23.
Mosad SM, Elsayed MM, Hammad EM, Hendam BM, Ali HS, Eladl AH, Saif MA. Genotype classification and pathogenicity of infectious bursal disease virus circulating in vaccinated broiler chicken farms. Vet Res Commun. 2024;48:3089–104.
24.
Lu Z, Zhang L, Wang N, Chen Y, Gao L, Wang Y, Gao H, Gao Y, Li K, Qi X, Wang X. Naturally occurring reassortant infectious bursal disease virus in northern China. Virus Res. 2015;203:92–5.
25.
Sun JH, Lu P, Yan YX, Hua XG, Jiang J, Zhao Y. Sequence and analysis of genomic segment A and B of very virulent infectious bursal disease virus isolated from China. J Vet Med B Infect Dis Vet Public Health. 2003;50:148–54.
26.
Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49:W293–6.
27.
Ali Khan RS, Habib M, Ali W, Salah Ud Din Shah M, Ashraf A, Ali Tahir Z, Helal ZH, Khan MI, Mahboob S, K AA-G, Al-Misned F. Phylogenetic analysis of Infectious Bursal Disease viruses according to newly proposed model of classification into geno-groups. J Infect Public Health. 2019;12:410–8.
28.
Berg TP. Acute infectious bursal disease in poultry: a review. Avian Pathol. 2000;29:175–94.
29.
Mutinda WU, Nyaga PN, Mbuthia PG, Bebora LC, Muchemi G. Risk factors associated with infectious bursal disease vaccination failures in broiler farms in Kenya. Trop Anim Health Prod. 2014;46:603–8.
30.
Niewiesk S. Maternal antibodies: clinical significance, mechanism of interference with immune responses, and possible vaccination strategies. Front Immunol. 2014;5:446.
31.
Jaton J, Gomez E, Lucero MS, Rizzi L, Gravisaco MJ, Pinto S, Berinstein A, Chimeno Zoth S. Evasion of maternal antibody protection by an IBDV Argentine variant. Poult Sci. 2024;103:103431.
32.
He X, Wei P, Yang X, Guan D, Wang G, Qin A. Molecular epidemiology of infectious bursal disease viruses isolated from Southern China during the years 2000–2010. Virus Genes. 2012;45:246–55.
33.
Mondal SP, Naqi SA. Maternal antibody to infectious bronchitis virus: its role in protection against infection and development of active immunity to vaccine. Vet Immunol Immunopathol. 2001;79:31–40.
34.
Abd El-Fatah AH, Ayman D, Samir M, Eid S, Elgamal M, El-Sanousi AA, Ibrahim M, AlKhazindar M, Ali MM, Afify A. Molecular characterization of circulating infectious bursal disease viruses in chickens from different Egyptian governorates during 2023. Virol J. 2024;21:312.
35.
Wang W, Shi J, Zhang Y, Qiao Y, Zuo W, Wang Z, Nong T, Hu S, Chen Y, Kong R, et al. Genetic and pathogenic characterizations of a naturally occurring reassortant and homologous recombinant strain of the classical infectious bursal disease virus re-emerging in chickens in southern China. Front Microbiol. 2023;14:1293072.
36.
Pikula A, Lisowska A, Jasik A, Perez LJ. The Novel Genetic Background of Infectious Bursal Disease Virus Strains Emerging from the Action of Positive Selection. Viruses 2021, 13.
37.
Cubas-Gaona LL, Flageul A, Courtillon C, Briand FX, Contrant M, Bougeard S, Lucas P, Quenault H, Leroux A, Keita A, et al. Genome Evolution of Two Genetically Homogeneous Infectious Bursal Disease Virus Strains During Passages in vitro and ex vivo in the Presence of a Mutagenic Nucleoside Analog. Front Microbiol. 2021;12:678563.
38.
Tahir I, Alsayeqh AF. Phytochemicals: a promising approach to control infectious bursal disease. Front Vet Sci. 2024;11:1421668.
39.
Alkie TN, Rautenschlein S. Infectious bursal disease virus in poultry: current status and future prospects. Vet Med (Auckl). 2016;7:9–18.
40.
Xiong H, Wu J, Xie Q, Li T, Wan Z, Qin A, Ye J, Shao H. Q221K mutation in VP2 drives antigenic shift of infectious bursal disease virus. Front Immunol. 2025;16:1600371.
41.
Qi X, Zhang L, Chen Y, Gao L, Wu G, Qin L, Wang Y, Ren X, Gao Y, Gao H, Wang X. Mutations of residues 249 and 256 in VP2 are involved in the replication and virulence of infectious Bursal disease virus. PLoS ONE. 2013;8:e70982.
42.
Gao H, Zhang S, Chang H, Guo Y, Li Z, Wang Y, Gao L, Li X, Cao H, Zheng SJ. Generation of a novel attenuated IBDV vaccine strain by mutation of critical amino acids in IBDV VP5. Vaccine 2024, 42:126081.
43.
Asfor AS, Reddy V, Nazki S, Urbaniec J, Brodrick AJ, Broadbent AJ. Modeling Infectious Bursal Disease Virus (IBDV) Antigenic Drift In Vitro. Viruses 2022, 15.
44.
Ramon G, Legnardi M, Cecchinato M, Cazaban C, Tucciarone CM, Fiorentini L, Gambi L, Mato T, Berto G, Koutoulis K, Franzo G. Efficacy of live attenuated, vector and immune complex infectious bursal disease virus (IBDV) vaccines in preventing field strain bursa colonization: A European multicentric study. Front Vet Sci. 2022;9:978901.
45.
Ravikumar R, Chan J, Prabakaran M. Vaccines against Major Poultry Viral Diseases: Strategies to Improve the Breadth and Protective Efficacy. Viruses 2022, 14.
46.
Bodman-Harris O, Rollier CS, Iqbal M. Approaches to Enhance the Potency of Vaccines in Chickens. Vaccines (Basel) 2024, 12.
Total words in MS: 4818
Total words in Title: 20
Total words in Abstract: 212
Total Keyword count: 5
Total Images in MS: 4
Total Tables in MS: 0
Total Reference count: 46