Hepatitis E genotype 3 genome: A comprehensive analysis of entropy, motif conservation, relevant mutations, and clade-associated polymorphisms

Abstract

Hepatitis E virus genotype 3 (HEV-3) is an EU/EEA emergent zoonosis. HEV-3 clades/subtypes have been described. Its genome contains ORF1, which encodes nonstructural proteins for virus replication, ORF2, the capsid protein, and ORF3, a multifunctional protein involved in virion pathogenesis. The study aims with respect to HEV-3 are to: (1) calculate genome entropy (excluding hypervariable region); (2) analyze the described motifs/mutations; (3) characterize clade/subtype genome polymorphisms. Seven hundred and five sequences from the GenBank database were used. The highest entropies were identified in zoonotic genotypes (HEV-3 and HEV-4) with respect to HEV-1 in X domain, RdRp, ORF2, and ORF3. There were statistically significant differences in the entropy between proteins, protease and ORF3 being the most variable and Y domain being the most conserved. Methyltransferase and Y domain motifs were completely conserved. By contrast, essential protease H581 residue and catalytic dyad exhibited amino acid changes in 1.8% and 0.4% of sequences, respectively. Several X domain amino acids were associated with clades. We found sequences with mutations in all helicase motifs except number IV. Helicase mutations related to increased virulence and/or fulminant hepatitis were frequent, the 1,110 residue being a typical HEV-3e and HEV-3f-A2 polymorphism. RdRp motifs III, V, VII also had high mutation rates. Motif III included residues that are polymorphisms of HEV-3e (F1449) and HEV-3 m (D1451). RdRp ribavirin resistance mutations were frequent, mainly 1479I (67.4, 100% in HEV-3efglmk) and 1634R/K (10.0%, almost 100% in HEV-3e). With respect to ORF2, 19/27 neutralization epitopes had mutations. The S80 residue in ORF3 presented mutations in 3.5% of cases. Amino acids in the ORF3-PSAP motif had high substitution rates, being more frequent in the first PSAP (44.8%) than in the second (1.5%). This is the first comprehensive analysis of the HEV-3 genome, aimed at improving our knowledge of the genome, and establishing the basis for future genotype-to-phenotype analysis, given that viral features associated with severity have not been explored in depth. Our results demonstrate there are important genetic differences in the studied genomes that sometimes affect significant viral structures, and constitute clade/subtype polymorphisms that may affect the clinical course or treatment efficacy.