Abstract
Sequence diversity of the hepatitis D virus (HDV) may impact viral clearance, contributing to the development of chronic infection. T-Cell-induced selection pressure and viral recombination can induce diversity throughout the viral genome including coding and noncoding regions, with the former potentially impacting viral pathogenicity and the latter exerting regulatory functions. Here, we aim to assess sequence variations of the HDV genome within and across HDV genotypes. Sequences from 721 complete HDV genomes and 793 large hepatitis D antigen (L-HDAg) regions belonging to all eight genotypes and published through December 2023 were compiled. Most retrieved sequences belonged to Genotype 1, whereas for Genotype 8, the fewest sequences were available. Alignments were conducted using Clustal Omega and Multiple Alignment using Fast Fourier Transform. Phylogeny was analysed using SplitsTree4, and recombination sites were inspected using Recombination Detection Program 4. All reported sequences were aligned per genotype to retrieve consensus and reference sequences based on the highest similarity to consensus per genotype. L-HDAg alignments of the proposed reference sequences showed that not only conserved but also highly variable positions exist, which was also reflected in the epitope variability across HDV genotypes. Importantly, in silico binding prediction analysis showed that CD8(+) T-cell epitopes mapped for Genotype 1 may not bind to major histocompatibility complex class I when examining their corresponding sequence in other genotypes. Phylogenetic analysis showed evidence of recombinant genomes within each individual genotype as well as between two different HDV genotypes, enabling the identification of common recombination sites. The identification of conserved regions within the L-HDAg allows their exploitation for genotype-independent diagnostic and therapeutic strategies, while the harmonized use of the proposed reference sequences may facilitate efforts to achieve HDV control.