Abstract
Duck hepatitis A virus (DHAV) is one of key pathogens for duck viral hepatitis, especially in Asian duck industry. Currently, two main genotypes (DHAV-1 and -3) exist. To explore insightfully the evolutionary character, we assessed the available 141 full-length genome sequences of DHAV isolated in 1986-2020 globally and divided DHAV-1 and DHAV-3 into further seven (DHAV-1 a-g) and five (DHAV-3 a-e) sub-clades, respectively. Phylogenetic and phylogeographic network analyses indicated great genetic diversity of DHAV identified in China, where the DHAV-1 cluster and DHAV-3 cluster were linked by virus strain HDHV1-BJ (GenBank ID: FJ157172.1) and Du_CH_LSD_090612 (GenBank ID: JF828995.1) via a long mutational branch and intermediate strains. Several strains previously identified as DHAV-1 according to the partial gene sequences were actually clustered within DHAV-3 in full-length genome-based analysis. Furthermore, we identified 32 recombination events across virus genome with the recombination hotspot at the 5' end and upstream of the capsid coding region. The highest variability of DHAV polyprotein was shown at the upstream region of the N terminus P-loop region, e.g., amino acids 672-716, followed by the aa 334-359 in the Capsid encoding region. The results presented here provides a robust insight into the genetic exchange patterns of DHAV genomes during the past decades, which may be used to map the evolutionary history and facilitate preventive measures of DHAVs.