Phylodynamic and phylogeographic reconstruction of beak and feather disease virus epidemiology and its implications for the international exotic bird trade

Beak and feather disease virus (BFDV) infects domestic and wild psittacine species and is able to cause progressive beak, claw and feather malformation and necrosis. In addition to having an impact on the health and welfare of domesticated birds, BFDV represents a significant threat to wild endangered species. Understanding the epidemiology, dynamics, viral migration rate, interaction between wild and domestic animals and the effect of implemented control strategies is fundamental in controlling the spread of the disease. With this in mind, a phylodynamic and phylogeographic analysis has been performed on a database of more than 400 replication‐associated protein (Rep) gene (ORF1) sequences downloaded from Genbank including some recently generated sequences from fifteen samples collected in Namibia. The results allowed us to reconstruct the variation of viral population size and demonstrated the effect of enforced international bans on these dynamics. A good correlation was found between viral migration rate and the intensity of animal trade between regions over time. A dominant flux of viral strains was observed from wild to domestic populations, highlighting the directionality of viral transmission and the risk associated with the capturing and trade of wild birds. Nevertheless, the flow of viruses from domestic to wild species was not negligible and should be considered as a threat to biodiversity. Therefore, considering the strong relationship demonstrated in this study between animal trade and BFDV viral fluxes more effort should be made to prevent contact opportunities between wild and domestic populations from different countries in order to control disease spread.