Molecular phylodynamics of infectious bursal disease viruses

The present study was conducted to study the molecular phylodynamics of the Indian field IBDVs. A total of 13 organized commercial poultry farms and 29 village poultry flocks were recruited in the study. The broiler flocks showed 15.25-60.18% mortality, followed by 12.4% in improved native poultry varieties and 5% in indigenous birds. The 664 bp hypervariable VP2 gene fragment of Western and Central Indian vvIBDVs showed 97.14-98.79 and 94.49-96.69% identity to Pakistani and South Indian vvIBDVs, respectively. An isolate was 99.54% identical to the Ventri-Plus vaccine strain, while three IBDVs showed maximum identity with the Georgia strain. Out of 22, 19 strains showed typical vvIBDV amino acid signature, while three strains showed substitutions specific to classical IBDVs. Central Indian vvIBDVs showed conserved substitutions at N212D and E300A, which can be used as a regional marker. Phylogenetic genogrouping placed global IBDVs into seven genogroups based upon virulence and geographical distribution. Nineteen field vvIBDVs were placed in the G3 genogroup, and the other three were grouped with classical IBDVs in G1 genogroup. A nucleotide span from 584 to 1248 covering VP2 hypervariable fragment was found suitable for correct genogrouping of field IBDVs. The Bayesian evolutionary analysis showed tMRCA of the year 2009 for 8 Western Indian vvIBDVs with vvIBDV from Pakistan. Central Indian vvIBDVs were evolved in the year 1991 from BD-3 and PY12 strains of vvIBDVs from Bangladesh and Pondicherry, respectively. An isolate showed evolution in year 2010 from the Nigerian ABIC strain, while three classical strains showed tMRCA of the year 2009 with the Georgia strain as a recent common ancestor.