The phylogenetic and phylogeographic landscape of the beak and feather disease virus, 1996-2022

The beak and feather disease virus (BFDV), causative agent of Psittacine beak and feather disease (PBFD), is a highly fatal and widespread virus that infects both the wild and captive Psittaciformes around the world. The BFDV genome is a ssDNA of approximately 2 kb in size, making it among the smallest known pathogenic viruses. Though, the virus is placed in Circoviridae family of the Circovirus genus, there is no classification system on clade and sub-clade level according to the International Committee on Taxonomy of Viruses and the strains are grouped on the bases of geographic locations. Thus, we provide the latest and robust phylogenetic classification of BFDVs in this study based on full-length genomic sequences, grouping all the available 454 strains detected during 1996-2022 into two distinct clades, e.g., GI and GII. The GI clade is further divided into six sub-clades (GI a-f), while GII into two sub-clades (GII a and b). In addition, the phylogeographic network identified high variability among the BFDV strains, showing several branches, where all the branches are connected to four strains, e.g., BFDV-ZA-PGM-70A(GenBank ID: HM748921.1, 2008-South Africa), BFDV-ZA-PGM-81A(GenBank ID: JX221009.1, 2008-South Africa), BFDV14(GenBank ID: GU015021.1, 2010-Thailand) and BFDV-isolate-9IT11(GenBank ID: KF723390.1, 2014-Italy). Furthermore, we identified 27 recombination events in the rep (replication-associated protein) and cap (capsid protein) coding regions using the complete genomes of BFDVs. Similarly, the amino acids variability analysis indicated that both the rep and cap regions are highly variable with values exceeding the variability coefficient estimation limit of 1.00, speculating the possible amino acids drift with the emergence of new strains. The findings provided in this study may offer the latest phylogenetic, phylogeographic and evolutionary landscape of the BFDVs.