Evidence for positive selection of hepatitis A virus antigenic variants in vaccinated men-having-sex-with men patients: Implications for immunization policies

Figures

• 

  Fig. 1. Diagram of the fragments used for the deep-sequencing analysis covering, the 5′ non-coding-region (5′NCR), the complete capsid region (VP4, VP2, VP3, and VP1), and a fragment of the 2A region.

• 

  Fig. 2. Phylogenetic analysis of HAV strains isolated from patients of the men-having-sex-with-men (MSM) 2016–2018 outbreak from Barcelona. (a) Phylogenetic tree based on a 255-bp fragment (spanning positions 2958–3230) from the VP1/2A protein. (b) Phylogenetic tree based on a sequence of 91-amino acids (spanning positions 251–300 from VP1 and 1–41 from 2A). (c) Phylogenetic tree based on a 1004-bp fragment (spanning positions 2208–3212) including the complete VP1 sequence and a short fragment of 2A. (d) Phylogenetic tree based on a sequence of 334-amino acids (spanning positions 1–300 from VP1 and 1–35 from 2A). Nucleotide dendrograms were constructed using the neighbour-joining method with distance calculation by the Kimura-2-parameter, and performing a bootstrap of 1000 replicates. Amino acid dendograms were constructed using the Neighbour-Joining method with distance calculation using the Poisson correction method. Samples from vaccinated individuals (n = 5) are underlined, samples from HIV-infected patients (n = 6) are indicated with the “†” symbol. Samples further analysed by MiSeq Illumina (n = 13) are indicated with “*”. The three main circulating strains in Europe related to…

• 

  Fig. 3. HAV viral loads. (a) Boxplot of the average virus load in sera from the non-vaccinated and vaccinated groups. (b) Boxplot of the average virus load in sera from HIV-non-infected and HIV-infected groups. Bottom and upper limits of the boxes represent the 25th and 75th percentiles, respectively. When present (n > 9), bottom and upper whiskers represent 10th and 90th percentiles, respectively. The median is represented by a solid black line and the average as a dotted red line.

• 

  Fig. 4. Informative content (IC) of each amino acid position along the capsid, and analysis of the frequency of the NSyn mutations. (a) The analysis of the informative content was performed using overlapping windows of 10 amino acids sliding every 3 amino acids. IC values (in bits) range from 0 (maximum variability) to 4·32 (minimum variability). Each point (amino acid position) corresponds to the mean IC from all the sequences of a given group. The black line corresponds to the non-vaccinated group and the red line to the vaccinated group. Values have been normalized by the frequency of each substitution and the number of sequences per group. (b) Frequency (%) of the NSyn mutations detected in the non-vaccinated patients (n = 8; black dots) and vaccinated patients (n = 5; red dots); black dots are more common due to the higher n. The horizontal green line represents the 5% frequency; points over this line represent high-frequency mutations.

• 

  Fig. 5. Location of the variants affecting residues in and around the antigenic sites of the HAV capsid. (a) Protomers from non-vaccinated patients (M1, M2, M4, M9, M10, M17, M30, and M47). (b) Protomers from vaccinated patients (M5, M6, M7, M46 and M48). Green protein corresponds to VP2, red to VP3 and blue to VP1. Protomers used to locate the replacements were extracted from the recently resolved crystal structure of HAV [18]. Epitopes (yellow residues) were identified either by mAb escape mutants [19,20] (VP1 and VP3) or mAb footprints [21] (VP2 and VP3). Grey residues represent mutations directly affecting the epitopes. Clear blue residues represent mutations located at a distance of 1–2 residues from the epitopes identified in monoclonal antibodies (mAbs) resistant mutants.

• 

  Fig. 6. Analysis of NSyn mutations located in and around the main antigenic sites in the non-vaccinated and vaccinated groups. (a) Percent of all NSyn mutations located in or around the antigenic sites distributed in the VP1 immunodominant site (dark blue), the VP1 glycophorin A binding site (red) and the VP2, VP3 putative receptor binding (clear blue). (b) Average of the abundance of haplotypes bearing the NSyn mutations from (a).