Next Article in Journal
The Potential of Phage Therapy against the Emerging Opportunistic Pathogen Stenotrophomonas maltophilia
Next Article in Special Issue
Human Astrovirus Outbreak in a Daycare Center and Propagation among Household Contacts
Previous Article in Journal
A Bioluminescent Biosensor for Quantifying the Interaction of SARS-CoV-2 and Its Receptor ACE2 in Cells and In Vitro
Previous Article in Special Issue
Organoids to Dissect Gastrointestinal Virus–Host Interactions: What Have We Learned?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 6
10.3390/v13061056
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evidence of Circulation of Several HAV Genetic Variants and Emergence of Potential Antigenic Variants in an Endemo-Epidemic Country before Vaccine Introduction

by
Kaouther Ayouni
1,2,*,
Anissa Chouikha
1,*,
Oussema Khamessi
3,
Henda Touzi
1,
Walid Hammemi
1 and
Henda Triki
1,4
1
Laboratory of Clinical Virology, WHO Reference Laboratory for Poliomyelitis and Measles in the Eastern Mediterranean Region, Pasteur Institute of Tunis, University Tunis El Manar (UTM), Tunis 1002, Tunisia
2
Faculty of Sciences of Tunis, University of Tunis El Manar, Campus Universitaire, El Manar, Tunis 2092, Tunisia
3
Institut Pasteur de Tunis, Université de Tunis El Manar, LR11IPT08 Venins et Biomolecules Therapeutiques, Tunis 1002, Tunisia
4
Faculty of Medicine of Tunis, University of Tunis El Manar, 15 Rue Djebel Lakhdhar, La Rabta, Tunis 1007, Tunisia
*
Authors to whom correspondence should be addressed.
Viruses 2021, 13(6), 1056; https://doi.org/10.3390/v13061056
Submission received: 15 April 2021 / Revised: 20 May 2021 / Accepted: 27 May 2021 / Published: 3 June 2021
(This article belongs to the Special Issue Gastroenteritis Viruses 2021)
Browse Figures
Versions Notes

Abstract

:
Similar to several other countries in the world, the epidemiology of hepatitis A virus changed from high to intermediate endemicity level in Tunisia, which led to the occurrence of outbreaks. This study aimed to determine the genetic and antigenic variability of HAV strains circulating in Tunisia during the last few years. Genotyping using complete VP1 gene and VP1-2A junction confirmed the predominance of genotype IA, with co-circulation of several genetic and antigenic variants. Phylogenetic analysis including Tunisian and strains from other regions of the world showed the presence of at least two IA-variants within IA subgenotype. Amino-acid analysis showed several mutations in or close to epitope regions in the VP1-region. This study provides a baseline on the genetic and antigenic variability of HAV circulating strains before the introduction of vaccination into the national immunization schedule.

1. Introduction

Hepatitis A infection remains a major public health problem despite the presence of effective vaccines since 1992. The infection is distributed worldwide but occurs more frequently in less developed regions with poor hygienic and sanitary conditions, causing substantial morbidity and mortality rates [1,2]. WHO estimates that hepatitis A caused approximately 7134 deaths in 2016, which accounts for 0.5% of the mortality worldwide due to viral hepatitis [3].
Hepatitis A infection is caused by hepatitis A virus (HAV), an RNA virus classified within the Picornaviridae family and Hepatovirus genus. HAV exists in a dual phenotype, naked and quasi-enveloped virions [4]. The quasi-enveloped virions represent the immature particles that are found in the blood, whereas the naked virions in feces are those released form their quasi-envelope by the action of bile salts in the passage from the bile ducts to the gut [4,5,6]. The naked virions are mature and contain the fully processed VP1 protein [4,5,6]. HAV genome is a positive single stranded RNA of approximately 7.5 Kb which comprises a single open reading frame (ORF) and two untranslated regions on its 3′and 5′ ends. Genotyping analysis of different HAV genomic regions revealed the existence of 6 different genotypes [7,8,9,10]. Genotypes I to III, sub-divided into sub-genotype IA, IB, IC, IIA, IIB, IIIA and IIIB, infect humans, while genotypes IV to VI are of simian origin [11,12].
The severity of the disease is age dependent. Unlike young children, teenagers and adults usually develop symptomatic clinical forms, sometimes with high severity [13,14]. HAV infection is more common in developing countries, most children get infected before the age of 10 years and older children and adults are generally immune. In these countries, symptomatic disease rates are low and outbreaks are uncommon. In countries with higher economic level, the HAV infection occurs in older ages and symptomatic cases as well as outbreaks are more frequent. During the past years, an epidemiological change of hepatitis A pattern was observed in many parts of the world [15,16,17,18,19,20,21,22,23,24]. Globally, an increase of symptomatic cases was observed [25,26]. Several outbreaks have been triggered in countries whose epidemiological status has changed from high to medium or low endemicity with more severe outcomes [25,27]. Thus, the implementation of programs for HAV surveillance and the introduction of HAV vaccine stand out as the best measure to limit and prevent the spread of the virus, especially in countries with intermediate or low endemicity levels.
In Tunisia, a transition from high to intermediate endemicity level was observed [28,29,30]. Improvement of socioeconomic and hygienic conditions had contributed to the decreasing frequency of the infection among young children, especially among primary school attendants, and led to serious outbreaks and significant disease burden. The national health authorities have committed, in line with the World Health Organization’s (WHO) resolution, to eliminate hepatitis viruses by 2030 as a public health threat. In this context, novel measures to limit the transmission of the disease were undertaken, by implementing an HAV surveillance program and introducing HAV vaccine. In October 2018, HAV vaccination was introduced for 6-year-old children at school entry. Since September 2020, the national HAV vaccination schedule has combined a single dose at 12 months with a catch-up vaccination for non-vaccinated children aged 6 years.
The knowledge of the molecular epidemiology of viral strains circulating is a key component to evaluate the vaccine efficiency, to follow progression of the national program and to detect a possible emergence of vaccine-escape variants. Only few studies have investigated the molecular patterns of HAV in Tunisian population. The most recent study identified HAV strains circulating during 2008–2013 in North Central Tunisia and its migration pattern [31]. The present retrospective study aims to provide an overview of the genetic diversity of HAV strains circulating in Tunisia between 2013 and 2018, prior to the vaccine introduction into the national immunization program.

2. Materials and Methods

2.1. Samples

Out of 103 HAV IgM positive serum samples, 77 (74.8%) were included in this study. The inclusion criteria were positivity of both anti-HAV IgM and RNA detection by PCR amplification in the 5′UTR genomic regions (position: 55–678, according to reference strains HM175). The samples were collected as part of diagnostic activity of the Laboratory of Clinical Virology of Pasteur Institute of Tunis, between 2013 and 2018, and are originated from 7 districts; in the Northern, the Eastern and the central regions of the country. The samples were stored at −20 °C. No identifying patient data were used in this study.

2.2. RNA Extraction and RT-PCR

Viral RNA was extracted from 140 μL of serum using the QIAmp viral RNA mini Kit (QIAGEN, Hilden-Germany) following the manufacturer’s instructions. The RNA extracted was used as a template for the amplification of the complete VP1 gene (954 bp) and the VP1-2A junction (394bp) by nested RT-PCR. Reverse transcription and amplification were performed using previously published primers [8,31,32]. For both regions, the cDNA was synthetized using 10 µL of RNA and anti-sens primer m2: 5′-AGTCACACCTCTCCAGGAA-3′ and HAV3381: 5′-CCATYTCAAGAGTCCACACACT-3′ for the complete VP1 gene and the VP1-2A junction, respectively. The amplification of the VP1-2A junction was then carried out with 10 µL of cDNA and 0.2 µM of outer primers: HAV2870 (forward: 5′-GACAGATTCYACATTTGGATTGGT-3′) and HAV3381 (reverse: 5′-CCATYTCAAGAGTCCACACACT-3′) at 94 °C for 5 min followed by 35 cycles (30 s at 94 °C, 30 s at 58 °C, 1 min at 72 °C) and a final elongation at 72 °C for 5 min. Nested PCR was carried out with 5 µL of the product of the first PCR and 0.2 µM of inner primers HAV2896 (forward: 5′-CTATTCAGATTGCAAATTAYAAT-3′) and HAV3289 (reverse: 5′-AAYTTCATYATTTCATGCTCCT-3′), following the same thermal condition with the exception of the hybridation temperature 55 °C for the nested PCR. Total VP1 gene was amplified using outer primers HAV-2167 (forward: 5′-GTTTTGCTCCTCTTTATCATG-3′) and m2 (reverse: 5′-AGTCACACCTCTCCAGGAA-3′) and inner primers HAV m1 (forward: 5′-GCTCCTCTTTATCATGCTATG-3′) and HAV3125 (reverse: 5′-CCTGCATTCTATATGACTCT-3′), with same amplification duration for both PCR rounds: 94 °C for 5 min followed by 40 cycles (30 s at 94 °C, 30 s at 53 °C for first PCR or at 55 °C for the second PCR, 1 min at 72 °C) and a final elongation at 72 °C for 5 min. The PCR products were visualized on a 1% agarose gel.

2.3. Amplicon Purification and Sequencing

The PCR products of the VP1-2A junction and the whole VP1 were purified using the QIAquick PCR purification Kit (QIAGEN, GmnH, Hilden, Germany) following manufacturer’s instructions. Purified PCR-products were sequenced by inner primers of the nested PCR on an ABI Prism 3130-Genetic Analyser (3130-Genetic Analyser Applied Biosystems) using the Big Dye terminator ready reaction cycle sequencing Kit (Applied Biosystems).

2.4. Sequence Analyses and Genotyping

The sequence of each isolate was deduced by aligning the respective forward and reverse sequences using Chromas software version 2.6.2. Genotyping was performed using phylogenetic analysis comparing the VP1-2A and VP1 sequences with reference sequences representing the different HAV genotypes and sub-genotypes IA to IIIB. GenBank accessions numbers of references sequences used in this study were as follows: X75215, AF357222, AB020566 for genotype IA, M14707, M20273, AF314208 for genotype IB, HQ401240 for genotype IC, AY644676 for genotype IIA, AY644670 for genotype IIB, AB279733, AY644337 for genotype IIIA and AB279735, AB258387 for genotype IIIB. Different strains previously submitted in the nucleotide GenBank collection from Tunisia and other countries were selected, to study nucleotide and antigenic variability. Phylogenetic trees were constructed using MEGA software version 7.0.26 by the maximum likelihood method and the kimura-2 parameter model [33]. Topology was supported by 1000 bootstrap replicates. The new HAV sequences generated as part of the present work were submitted to GenBank database under accession numbers: MW117967 to MW118018 and MW222090 to MW222114 for the complete VP1 with a partial 2A region, and VP1-2A junction, respectively.

2.5. Molecular Modeling

The three-dimensional structure of HAV-VP1 protein (accession number: MW118004) was predicted using Modeller version 9.24 [34]. Its amino acid sequence was compared with other sequences retrieved from nrNCBI database using FASTA and BLAST. The Crystal structure of HAV (PDB code 4QPI) was used as a template, in which the VP1 protein (Chain_A) was extracted using MMTSB tool set [12]. The model corresponding to the best value of the DOPE score was selected after generating 1000 conformers [35]. All structures were visually explored using PyMOL molecular Viewer [36].

3. Results

3.1. Genetic and Phylogenetic Analysis

The 393 nucleotides in the VP1-2A junction could be amplified and sequenced from the 77 samples included in the study. Fifty-two complete sequences of the VP1 gene (954 nt long) were obtained. The phylogenetic trees in Figure 1 compare the Tunisian sequences in the entire VP1 region (Figure 1A) and the VP1-2A junction (Figure 1B) with reference sequences representative of the different HAV genotypes and subtypes. Both trees show that all isolates belonged to HAV genotype IA with similar grouping of the Tunisian sequences into 5 clusters (designed (a) to (d)); more consistent bootstrap values in the complete VP1 region were observed.
To better understand the molecular epidemiology of circulating strains in Tunisia, the 77 VP1-2A junctions detected in this work (2013–2018) were aligned with 104 previously published Tunisian strains of genotype IA collected during 2001–2013 period (Figure 2). The 104 sequences were selected among a total of 187 sequences published in Genbank. For identical sequences detected in the same year, only one representative sequence was selected for the phylogenetic analysis. Similar to Figure 1, the tree in Figure 2 shows the presence of at least five different variants that circulated in Tunisia between 2001 and 2018. A major variant (Variant 2) comprised strains that circulated in different regions of the country from 2001 to 2018. Variants 1, 3, 4 and 5 included sequences that were detected during somewhat shorter periods in 2007–2018, 2001–2015, 2012–2018 and 2014–2018, respectively.
The phylogenetic comparison of the entire VP1 of HAV Tunisian strains with others published from 21 different countries of the world is shown in Figure 3. Based on the entire HAV VP1 sequence, the phylogenetic analysis shows the presence of three main clusters I to III. All the Tunisian sequences were grouped in one cluster (Cluster I) with sequences from France (4 sequences, represented in red), Spain (6 sequences, represented in blue) and Yugoslavia (2 sequences, represented in green). Cluster II enclosed sequences almost from all the different countries whereas the cluster III grouped only one sequence from Chile, one from France and one from Spain (Figure 3).

3.2. Amino Acid Analysis of the VP1 Protein

Amino acid analysis was performed for the 52 entire VP1 sequences obtained in this work by comparing the deduced amino acid with the prototype strain of genotype IA GBM/WT (accession numbers: X75215). Twelve different amino acid replacements were detected in the complete VP1 region (Table 1, Figure 4). Five were non-conservative replacements (position: VP1-1, VP1-2, VP1-6, VP1-26, VP1-99), one was semi-conservative (VP1-166) and the remaining ones were conservative (Table 1). The number of detected replacements varied from 3 to 5 per strain. Replacements VP1-26, VP1-37 and VP1-99 were found in all isolates; the one at position VP1-26 was reported as specific of genotype IB [37]. Figure 4 shows the location of replacements in relation to previously described epitope regions [38,39,40,41]. All amino acid replacements in the VP1 protein were located in or close to epitope regions (Figure 4). Ten of the amino acid replacements were located in B/T cell epitope regions 1–17, 6–17, 2–33, 10–33, 17–33, 11–25, 13–24, 15–20, 70–81, 99–107, 99–122 and 166–178 [12,37,38,39,40]. Replacement VP1-166 was located in epitope region 166–178 and also around the immunodominant site. Two replacements located in positions 37 and 64 were close to B/T cell epitope regions 10–33, 17–33 and 70–81. Among the detected amino acid replacements, only 2 out of the 12 replacements (positions 29 and 37) were previously described in Tunisia [42,43,44], the remaining ones are reported for the first time in the present work.
In addition to the twelve amino acid replacements, a deletion of six amino acids RWFFNL (position 128 to 133) was found in one sequence (accession number: MW118004) located at the B/T cell epitope region 115–139. To study the effect of this deletion on the VP1 protein, a model of MW118004 was generated and selected according to several evaluation criteria. Figure 5 shows the three-dimensional structure of the VP1 of this isolate (Figure 5b) as compared to the template HAV-VP1 protein (PDB: 4QPI) (Figure 5a). The analysis of the three-dimensional structure, determined by homology modeling, shows that the deleted region found in the sequence of this isolate (PSTLRWFFNLF/PASTL- - - - - -F) resulted in a change in the VP1 protein structure by causing a destruction of an α-helix in the mature protein (Figure 5). This destroyed α-helix being located in an epitope region is likely to induce physicochemical modifications of the properties of this epitope, which may affect the binding affinity of the responding B cell and also of T cell, given that this deletion are located in a peptide recognized by T cell [38].

4. Discussion

This study reports an in-depth analysis of the genetic variability of HAV in Tunisia during the last few years. The sampling is not representative of hepatitis A cases in Tunisia due to the high frequency of asymptomatic forms and of symptomatic forms that are not reported. However, this study provides an overview on the molecular characterization of circulating strains. Phylogenetic analysis of HAV strains collected from different Tunisian districts during a six year period was conducted. All of the detected sequences belonged to genotype IA. The predominance of this genotype with the co-circulation of genotype IB was previously reported in Tunisia in both environmental and clinical samples [31,42,43,44,46,47]. The predominance of genotype IA was also described globally [48]. Same genotyping results were obtained based on the VP1-2A junction or the entire VP1 gene by phylogenetic analysis. Several studies showed that HAV genotyping could be performed using either the entire VP1 region, the N terminus of the VP1 region, the VP1–2A junction or the VP1–2B region [11,49,50,51,52]. Our results confirm that both regions can be used for HAV genotyping, although the reliability of bootstrap values in the phylogenetic analysis is more consistent when the entire VP1 sequences is used. Thus, VP1-2A junction remains a good solution for HAV genotyping, but for more in-depth analysis, larger genomic regions are more helpful.
To better understand the evolution and the genetic variability of circulating HAV sequences, we compared the VP1-2A junction sequences obtained in this work with others described in previous studies. Phylogenetic analyses showed that the Tunisian HAV strains detected between 2001 and 2018 were divided into at least five variants. A major variant (Variant2) included strains that circulated from 2001 to 2018 which seems have a continuous circulation in the country. The other variants were detected during shorter periods, 2007–2018 (Variant1), 2012–2018 (Variant4), 2014–2018 (Variant5) and 2001–2015 (Variant3). Variant3 was not detected after 2015. These findings confirm the co-circulation of different HAV variants within the genotype IA in multiple areas of Tunisia, which indicates that the endemicity level remains quite important. In fact, the number of genetic variants of pathogens generally decreases together with the progress of control strategies. While variant2 has continuous circulation, variant1, 3 and 5 were detected during shorter period, but it cannot be excluded that they may circulated during longer period given the lack of performant molecular surveillance during the study period.
When genotype IA Tunisian sequences, including the ones reported herein and previously published ones, were compared to HAV sequences from other countries, we found that all sequences from Tunisia grouped together in one cluster with other sequences from France, Spain and Yugoslavia. The second Cluster II enclosed strains from Asia, America and Europe whereas the cluster III grouped only sequences from Chile, France and Spain. These results suggest the presence of at least three IA variants within the IA subgenotype. However, the relatively limited number of sequences covering the total VP1 region available in GenBank, which originates from only 21 countries, may be the reason of the detection of only these three IA variants. Most of the published sequences cover a short size segment (<900 bases) of the genome which is not enough to properly determine the genetic variability of circulating strains. The availability in the future of full VP1 sequences from other countries in the world will help to better define the number of variants circulating worldwide and assess the genetic variability within genotype IA.
In the second part of the present work, we studied the impact of nucleotide mutations in the capsid protein VP1 gene on the amino-acid sequences to better understand the antigenic diversity of circulating HAV strains. The VP1 protein is a structural protein, known to contain major immunodominant epitopes of HAV [53,54]. VP1 protein analysis is important for predicting the possible antigenic escaping mutants within the circulating strains. Comparison of the deduced VP1 amino acid sequences obtained in the present study with the reference sequence of genotype IA showed several amino acid differences, mostly located in the N terminus of the VP1 region. Three to five different amino-acid mutations were identified per strain, all these mutations were located in or close to B/T cell epitope regions [12,37,38,40,41,45]. Most of the amino acid replacements were non-conservative and located in epitope regions. Thus, strains bearing such localized amino acid changes might be considered as potential antigenic variants. Eight amino acid replacements located in positions VP1-1, VP1-2, VP1-3, VP1-6, VP1-18, VP1-64, VP1-72 and VP1-99 were described for the first time in our study. Replacements VP1-26, VP1-37 and VP1-99 were found in all isolates of the present study. Replacement VP1-26 was identified previously in Argentina [49]. Replacement VP1-37 (R → K) was previously described in Tunisia and also in Argentina [36,37,38,49]. Replacement VP1-99 (I → T) is, to our knowledge reported for the first time in this study. This replacement same as replacements VP1-26 (R → K) and VP1-29 (R → K) seem to be associated to subtype IB [37]. In fact, replacement VP1-29 (R → K) was previously identified in Tunisia and in India and was associated widely to subtype IB [37]. Furthermore, recombination between sub-genotype IA and IB was previously described in Tunisia and other countries, which may explain this result [43,55].
A non-conservative replacement at position VP1-166 was previously identified in three isolates in Spain, during two outbreaks among male patients having sex with men; one was in non-vaccinated patient and the two others were in partly vaccinated patients [41,56]. These three isolates were considered as antigenic variants. In the present work, a semi-conservative amino acid replacement (V → A) at the same position was identified in one isolate detected in non-vaccinated patient, bearing four other amino acid replacements; two non-conservatives (in B/T cell epitope regions) and two conservatives (one in and the other close to B/T cell epitope regions). Thus, this isolate can be also considered as an antigenic variant.
These findings emphasize the importance of following the genetic variability of circulating strains to determinate the effects of vaccination on the molecular epidemiology of hepatitis A and on the potential selection of vaccine-escape mutants. Indeed, the emergence of potential antigenic variants, as well as vaccine-escape mutants after vaccine implementation, were previously described in many studies, especially when a man has sex with another man, and in immunocompromised patients [41,45,56,57,58,59]. In addition, it was previously reported that after vaccination amino acid replacement occurring in and around the epitopes regions was observed in vaccinated and unvaccinated persons [56]. The abundance of these amino-acid replacements was significantly higher in vaccinated patients, suggesting a positive selection of antigenic variants in some vaccinated patients [56]. Further molecular studies on HAV strains targeting the entire capsid proteins (VP1, VP2 and VP3) will be of great interest to better understand the evolution of circulating strains before and after introduction of HAV vaccine.
Interestingly, a deletion of 6 amino acids (PSTL- - - - - - F) position 128 to 133 was found in another isolate of the present study (accession N°: MW118004), located in the B/T cell epitope region (position 115–139) and resulted in the destruction of an α-helix (Figure 5) in the mature protein. Despite this deletion, the reading frame was conserved providing a non-defective virus as it was characterized by the presence of anti-HAV IgM antibody in symptomatic patient sera. Deletion mutants of VP1 were identified in other studies [44,45]. Aragones et al., 2007, studied HAV isolates evolution under a selective pressure of monoclonal antibodies and showed the occurrence of several deletion mutants, among which four contains the total or a part of the deleted region identified in the present study [45]. In Tunisia, Khelifi et al., reported another deletion in one isolate, a 38 amino acid deletion in VP1 protein (position 150–187) located in the part of the immunodominant antigenic site [44]. The occurrence of deletions seems to play a role in the adaptation to different selective pressures such as host defense mechanisms or new environmental conditions.
Despite the high degree of conservation of the capsid amino acid sequence shown by many studies, we found some degree of heterogeneity which revealed the circulation of several antigenic variants in the Tunisian HAV viral population. A previous study reported the circulation of two antigenic variants in Tunisia during 2003 [44]. One presented 38 amino acid deletions in the part of the immunodominant antigenic site of VP1 region and the second variant presented a replacement at the positionVP1-10. Our results suggest that a variety of mutants and antigenic variants circulates in Tunisia, despite the relatively low mutation rate reported for HAV, which seems to be related to the strict structural constraints of the viral capsid and a restricted codon usage [45,60].

5. Conclusions

The present study provides recent data on the molecular characteristics of circulating HAV strains in Tunisia. Genotype IA predominates, and several genetic and antigenic variants within this genotype co-circulate. Our study provides a baseline on the genetic variability of circulating strains before the introduction of the anti-HAV vaccines into the national immunization schedule. Strengthening the national hepatitis A infection surveillance system will be of great interest to monitor the impact of the vaccination program into HAV molecular epidemiology and the possible emergence of vaccine-escaping variants under selective pressures of this vaccination.

Author Contributions

Conceptualization, H.T. (Henda Triki) and K.A.; methodology, K.A.; software, O.K., K.A and A.C.; validation, H.T. (Henda Triki), A.C. and K.A.; formal analysis, K.A. and O.K.; investigation, K.A., H.T. (Henda Touzi) and W.H.; resources, H.T. (Henda Touzi); data curation, K.A. and O.K.; writing—original draft preparation, K.A. and A.C.; writing—review and editing, A.C. and H.T. (Henda Triki); visualization, K.A. and A.C.; supervision, H.T. (Henda Triki); project administration, H.T. (Henda Triki); funding acquisition, H.T. (Henda Triki) All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Ministry of Higher Education (Reasearch Laboratory LR20IPT02: “Virus, Vectors and Hosts: One Health approach and technological innovation for a better health”) and the Clinical Investigation Center (CIC).

Institutional Review Board Statement

Samples used in the present study were obtained from the routine diagnostic activity of the laboratory of clinical Virology in Pasteur Institute of Tunis.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saffar, M.J.; Abedian, O.; Ajami, A.; Abedian, F.; Mirabi, A.M.; Khalilian, A.R.; Saffar, H. Age-specific seroprevalence of anti-hepatitis a antibody among 1–30 years old population of savadkuh, mazandaran, iran with literature review. Hepat. Mon. 2012, 12, 326–332. [Google Scholar] [CrossRef] [Green Version]
  2. Wasley, A.; Fiore, A.; Bell, B.P. Hepatitis A in the era of vaccination. Epidemiol. Rev. 2006, 28, 101–111. [Google Scholar] [CrossRef] [Green Version]
  3. World Health Organization. Hepatitis A. 2020 [updated 27 July 2020]. Available online: https://www.who.int/fr/news-room/fact-sheets/detail/hepatitis-a (accessed on 6 December 2020).
  4. Feng, Z.; Hensley, L.; McKnight, K.L.; Hu, F.; Madden, V.; Ping, L.; Jeong, S.H.; Walker, C.; Lanford, R.E.; Lemon, S.M. A pathogenic picornavirus acquires an envelope by hijacking cellular membranes. Nature 2013, 496, 367–371. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Feng, Z.; Lemon, S.M. Peek-a-boo: Membrane hijacking and the pathogenesis of viral hepatitis. Trends Microbiol. 2014, 22, 59–64. [Google Scholar] [CrossRef] [Green Version]
  6. Hirai-Yuki, A.; Hensley, L.; Whitmire, J.K.; Lemon, S.M. Biliary secretion of quasi-enveloped human hepatitis a virus. mBio 2016, 7. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Ching, K.Z.; Nakano, T.; Chapman, L.E.; Demby, A.; Robertson, B.H. Genetic characterization of wild-type genotype VII hepatitis A virus. J. Gen. Virol. 2002, 83 Pt 1, 53–60. [Google Scholar] [CrossRef]
  8. Costa-Mattioli, M.; Cristina, J.; Romero, H.; Perez-Bercof, R.; Casane, D.; Colina, R.; Garcia, L.; Vega, I.; Glikman, G.; Romanowsky, V.; et al. Molecular evolution of hepatitis A virus: A new classification based on the complete VP1 protein. J. Virol. 2002, 76, 9516–9525. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Lu, L.; Ching, K.Z.; de Paula, V.S.; Nakano, T.; Siegl, G.; Weitz, M.; Robertson, B.H. Characterization of the complete genomic sequence of genotype II hepatitis A virus (CF53/Berne isolate). J. Gen. Virol. 2004, 85 Pt 10, 2943–2952. [Google Scholar] [CrossRef]
  10. Robertson, B.H.; Jansen, R.W.; Khanna, B.; Totsuka, A.; Nainan, O.V.; Siegl, G.; Widell, A.; Margolis, H.S.; Isomura, S.; Ito, K.; et al. Genetic relatedness of hepatitis A virus strains recovered from different geographical regions. J. Gen. Virol. 1992, 73 Pt 6, 1365–1377. [Google Scholar] [CrossRef]
  11. Perez-Sautu, U.; Costafreda, M.I.; Lite, J.; Sala, R.; Barrabeig, I.; Bosch, A.; Pinto, R.M. Molecular epidemiology of hepatitis A virus infections in Catalonia, Spain, 2005–2009: Circulation of newly emerging strains. J. Clin. Virol. 2011, 52, 98–102. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, X.; Ren, J.; Gao, Q.; Hu, Z.; Sun, Y.; Li, X.; Rowlands, D.J.; Yin, W.; Wang, J.; Stuart, D.I.; et al. Hepatitis A virus and the origins of picornaviruses. Nature 2015, 517, 85–88. [Google Scholar] [CrossRef]
  13. Jung, Y.K.; Kim, J.H. Epidemiology and clinical features of acute hepatitis A: From the domestic perspective. Korean J. Hepatol. 2009, 15, 438–445. [Google Scholar] [CrossRef] [PubMed]
  14. Saberifiroozi, M. Hepatitis A virus infection: Is it an important hazard to public health? Hazards of HAV for public health. Hepat. Mon. 2011, 11, 235–237. [Google Scholar]
  15. Alkhalidi, J.; Alenezi, B.; Al-Mufti, S.; Hussain, E.; Askar, H.; Kemmer, N.; Neff, G.W. Seroepidemiology of hepatitis A virus in Kuwait. World J. Gastroenterol. 2009, 15, 102–105. [Google Scholar] [CrossRef] [PubMed]
  16. Almuneef, M.A.; Memish, Z.A.; Balkhy, H.H.; Qahtani, M.; Alotaibi, B.; Hajeer, A.; Qasim, L.; Al Knawy, B. Epidemiologic shift in the prevalence of Hepatitis A virus in Saudi Arabia: A case for routine Hepatitis A vaccination. Vaccine 2006, 24, 5599–5603. [Google Scholar] [CrossRef]
  17. Bener, A.; Al-Kaabi, S.; Derbala, M.; Al-Marri, A.; Rikabi, A. The epidemiology of viral hepatitis in Qatar. Saudi J. Kidney Dis. Transplant. 2009, 20, 300–306. [Google Scholar]
  18. Jacobsen, K.H. Hepatitis A virus in West Africa: Is an epidemiological transition beginning? Niger. Med. J. 2014, 55, 279–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Jacobsen, K.H.; Koopman, J.S. The effects of socioeconomic development on worldwide hepatitis A virus seroprevalence patterns. Int. J. Epidemiol. 2005, 34, 600–609. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Melgaco, J.G.; Morgado, L.N.; Santiago, M.A.; Oliveira, J.M.; Lewis-Ximenez, L.L.; Hasselmann, B.; Cruz, O.G.; Pinto, M.A.; Vitral, C.L. A single dose of inactivated hepatitis A vaccine promotes HAV-specific memory cellular response similar to that induced by a natural infection. Vaccine 2015, 33, 3813–3820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Memish, Z.A.; Knawy, B.A.; El-Saed, A. Incidence trends of viral hepatitis A, B, and C seropositivity over eight years of surveillance in Saudi Arabia. Int. J. Infect. Dis. 2010, 14, e115–e120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Sacy, R.G.; Haddad, M.; Baasiri, G.; Khoriati, A.; Gerbaka, B.J.; Abu-Elyazeed, R. Hepatitis a in Lebanon: A changing epidemiological pattern. Am. J. Trop. Med. Hyg. 2005, 73, 453–456. [Google Scholar] [CrossRef] [Green Version]
  23. Sharar, Z.A.; Rajah, J.; Parsons, H. Childhood seroprevalence of hepatitis A in the United Arab Emirates. Trop. Dr. 2008, 38, 65–66. [Google Scholar] [CrossRef] [PubMed]
  24. Tufenkeji, H. Hepatitis A shifting epidemiology in the Middle East and Africa. Vaccine 2000, 18 (Suppl. 1), S65–S67. [Google Scholar] [CrossRef]
  25. Jacobsen, K. The Global Prevalence of Hepatitis A Virus Infection and Susceptibility: A Systematic Review; World Health Organization: Geneva, Switzerland, 2009; pp. 161–187. [Google Scholar]
  26. Rein, D.; Stevens, G.; Flaxman, A.; Wittenborn, J.; Timothy, N.; Wiktor, S.; Wiersma, S. P703 The global burden of hepatitis a virus in 1990 and 2005. J. Hepatol. 2014, 60, S303. [Google Scholar] [CrossRef]
  27. Hanafiah, K.M.; Jacobsen, K.H.; Wiersma, S.T. Challenges to mapping the health risk of hepatitis A virus infection. Int. J. Health Geogr. 2011, 10, 57. [Google Scholar] [CrossRef] [Green Version]
  28. Ayouni, K.; Kharroubi, G.; Mallekh, R.; Hammami, W.; Marouani, R.; Mhamdi, M.; Salah, A.B.; Triki, H.; Bettaieb, J. Seroprevalence of hepatitis A virus infection in Central-West of Tunisia. J. Med. Virol. 2020. [Google Scholar] [CrossRef] [PubMed]
  29. Letaief, A.; Kaabia, N.; Gaha, R.; Bousaadia, A.; Lazrag, F.; Trabelsi, H.; Ghannem, H.; Jemni, L. Age-specific seroprevalence of hepatitis a among school children in central Tunisia. Am. J. Trop. Med. Hyg. 2005, 73, 40–43. [Google Scholar] [CrossRef] [Green Version]
  30. Rezig, D.; Ouneissa, R.; Mhiri, L.; Mejri, S.; Haddad-Boubaker, S.; Alaya, N.B.; Triki, H. Seroprevalences of hepatitis A and E infections in Tunisia. Pathol. Biol. (Paris) 2008, 56, 148–153. [Google Scholar] [CrossRef] [PubMed]
  31. Beji-Hamza, A.; Taffon, S.; Mhalla, S.; Lo Presti, A.; Equestre, M.; Chionne, P.; Madonna, E.; Cella, E.; Bruni, R.; Ciccozzi, M.; et al. Migration pattern of hepatitis A virus genotype IA in North-Central Tunisia. Virol. J. 2015, 12, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Said, R.; Wolfaardt, M.; Taylor, M.B. Molecular characterisation of hepatitis A virus strains from water sources in South Africa. Water Sci. Technol. 2014, 69, 923–933. [Google Scholar] [CrossRef]
  33. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Sali, A.; Blundell, T.L. Comparative protein modelling by satisfaction of spatial restraints. J. Mol. Biol. 1993, 234, 779–815. [Google Scholar] [CrossRef]
  35. Shen, M.Y.; Sali, A. Statistical potential for assessment and prediction of protein structures. Protein Sci. 2006, 15, 2507–2524. [Google Scholar] [CrossRef] [Green Version]
  36. Seeliger, D.; de Groot, B.L. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J. Comput. Aided Mol. Des. 2010, 24, 417–422. [Google Scholar] [CrossRef] [Green Version]
  37. Joshi, M.S.; Cherian, S.S.; Bhalla, S.; Chitambar, S.D. Longer duration of viremia and unique amino acid substitutions in a hepatitis A virus strain [corrected] associated with Guillain-Barre syndrome (GBS). J. Med Virol. 2010, 82, 913–919. [Google Scholar] [CrossRef] [PubMed]
  38. Ivanov, V.S.; Kulik, L.N.; Gabrielian, A.E.; Tchikin, L.D.; Kozhich, A.T.; Ivanov, V.T. Synthetic peptides in the determination of hepatitis A virus T-cell epitopes. FEBS Lett. 1994, 345, 159–161. [Google Scholar] [CrossRef] [Green Version]
  39. Kulik, L.N.; Ivanov, V.S.; Tchikin, L.D.; Ostrovsky, A.G.; Ivanov, V.T. Localization of a sequential B-epitope in the VP2 protein of hepatitis A virus. FEBS Lett. 1995, 367, 85–88. [Google Scholar] [CrossRef] [Green Version]
  40. Lemon, S.M.; Barclay, W.; Ferguson, M.; Murphy, P.; Jing, L.; Burke, K.; Wood, D.; Katrak, K.; Sangar, D.; Minor, P.D.; et al. Immunogenicity and antigenicity of chimeric picornaviruses which express hepatitis A virus (HAV) peptide sequences: Evidence for a neutralization domain near the amino terminus of VP1 of HAV. Virology 1992, 188, 285–295. [Google Scholar] [CrossRef]
  41. Perez-Sautu, U.; Costafreda, M.I.; Cayla, J.; Tortajada, C.; Lite, J.; Bosch, A.; Pinto, R.M. Hepatitis a virus vaccine escape variants and potential new serotype emergence. Emerg. Infect. Dis. 2011, 17, 734–737. [Google Scholar] [CrossRef]
  42. Gharbi-Khelifi, H.; Ben Salem Abid, N.; Sdiri, K.; Harrath, R.; Beji, A.; Bhiri, L.; Billaudel, S.; Ferre, V.; Aouni, M. Characterization of outbreak hepatitis a isolates in five Tunisian childcare centers. Braz. J. Microbiol. 2011, 42, 1204–1212. [Google Scholar] [CrossRef] [Green Version]
  43. Gharbi-Khelifi, H.; Ferre, V.; Sdiri, K.; Berthome, M.; Fki, L.; Harrath, R.; Billaudel, S.; Aouni, M. Hepatitis A in Tunisia: Phylogenetic analysis of hepatitis A virus from 2001 to 2004. J. Virol. Methods 2006, 138, 109–116. [Google Scholar] [CrossRef]
  44. Gharbi-Khelifi, H.; Sdiri, K.; Harrath, R.; Fki, L.; Hakim, H.; Berthome, M.; Billaudel, S.; Ferre, V.; Aouni, M. Genetic analysis of HAV strains in Tunisia reveals two new antigenic variants. Virus Genes 2007, 35, 155–159. [Google Scholar] [CrossRef] [PubMed]
  45. Aragones, L.; Bosch, A.; Pinto, R.M. Hepatitis A virus mutant spectra under the selective pressure of monoclonal antibodies: Codon usage constraints limit capsid variability. J. Virol. 2008, 82, 1688–1700. [Google Scholar] [CrossRef] [Green Version]
  46. Gharbi-Khelifi, H.; Sdiri, K.; Ferre, V.; Harrath, R.; Berthome, M.; Billaudel, S.; Aouni, M. A 1-year study of the epidemiology of hepatitis A virus in Tunisia. Clin. Microbiol. Infect. 2007, 13, 25–32. [Google Scholar] [CrossRef] [Green Version]
  47. Ouardani, I.; Turki, S.; Aouni, M.; Romalde, J.L. Detection and molecular characterization of hepatitis a virus from tunisian wastewater treatment plants with different secondary treatments. Appl. Environ. Microbiol. 2016, 82, 3834–3845. [Google Scholar] [CrossRef] [Green Version]
  48. Hartard, C.; Gantzer, C.; Bronowicki, J.P.; Schvoerer, E. Emerging hepatitis E virus compared with hepatitis A virus: A new sanitary challenge. Rev. Med. Virol. 2019, 29, e2078. [Google Scholar] [CrossRef] [PubMed]
  49. Bruni, R.; Taffon, S.; Equestre, M.; Chionne, P.; Madonna, E.; Rizzo, C.; Tosti, M.E.; Alfonsi, V.; Ricotta, L.; De Medici, D.; et al. Key role of sequencing to trace hepatitis a viruses circulating in italy during a large multi-country european foodborne outbreak in 2013. PLoS ONE 2016, 11, e0149642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Desbois, D.; Couturier, E.; Mackiewicz, V.; Graube, A.; Letort, M.J.; Dussaix, E.; Roque-Afonso, A.M. Epidemiology and genetic characterization of hepatitis A virus genotype IIA. J. Clin. Microbiol. 2010, 48, 3306–3315. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Legeay, O.; Caudrelier, Y.; Cordevant, C.; Rigottier-Gois, L.; Lange, M. Simplified procedure for detection of enteric pathogenic viruses in shellfish by RT-PCR. J. Virol. Methods 2000, 90, 1–14. [Google Scholar] [CrossRef]
  52. MacDonald, E.; Steens, A.; Stene-Johansen, K.; Gillesberg Lassen, S.; Midgley, S.; Lawrence, J.; Crofts, J.; Ngui, S.L.; Balogun, K.; Frank, C.; et al. Increase in hepatitis A in tourists from Denmark, England, Germany, The Netherlands, Norway and Sweden returning from Egypt, November 2012 to March 2013. Euro. Surveill. 2013, 18, 20468. [Google Scholar] [CrossRef] [PubMed]
  53. Gauna, A.; Losada, S.; Lorenzo, M.; Bermudez, H.; Toledo, M.; Perez, H.; Chacon, E.; Noya, O. Synthetic peptides for the immunodiagnosis of hepatitis A virus infection. J. Immunol. Methods 2015, 427, 1–5. [Google Scholar] [CrossRef] [PubMed]
  54. Wang, C.H.; Tschen, S.Y.; Heinricy, U.; Weber, M.; Flehmig, B. Immune response to hepatitis A virus capsid proteins after infection. J. Clin. Microbiol. 1996, 34, 707–713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Liu, W.; Zhai, J.; Liu, J.; Xie, Y. Identification of recombination between subgenotypes IA and IB of hepatitis A virus. Virus Genes 2010, 40, 222–224. [Google Scholar] [CrossRef] [PubMed]
  56. Sabria, A.; Gregori, J.; Garcia-Cehic, D.; Guix, S.; Pumarola, T.; Manzanares-Laya, S.; Cayla, J.A.; Bosch, A.; Quer, J.; Pinto, R.M. Evidence for positive selection of hepatitis A virus antigenic variants in vaccinated men-having-sex-with men patients: Implications for immunization policies. EBioMedicine 2019, 39, 348–357. [Google Scholar] [CrossRef]
  57. D’Andrea, L.; Perez-Rodriguez, F.J.; de Castellarnau, M.; Manzanares, S.; Lite, J.; Guix, S.; Bosch, A.; Pinto, R.M. Hepatitis A virus genotype distribution during a decade of universal vaccination of preadolescents. Int. J. Mol. Sci. 2015, 16, 6842–6854. [Google Scholar] [CrossRef] [Green Version]
  58. Pinto, R.M.; D’Andrea, L.; Perez-Rodriguez, F.J.; Costafreda, M.I.; Ribes, E.; Guix, S.; Bosch, A. Hepatitis A virus evolution and the potential emergence of new variants escaping the presently available vaccines. Future Microbiol. 2012, 7, 331–346. [Google Scholar] [CrossRef] [Green Version]
  59. Pinto, R.M.; Perez-Rodriguez, F.J.; D’Andrea, L.; de Castellarnau, M.; Guix, S.; Bosch, A. Hepatitis A virus codon usage: Implications for translation kinetics and capsid folding. Cold Spring Harb. Perspect. Med. 2018, 8. [Google Scholar] [CrossRef]
  60. Sanchez, G.; Bosch, A.; Pinto, R.M. Genome variability and capsid structural constraints of hepatitis a virus. J. Virol. 2003, 77, 452–459. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Viruses 13 01056 g001 550
Figure 1. Phylogenetic trees comparing nucleotide sequences of HAV strains isolated in Tunisia during 2013–2018 with the reference sequences. The tree was performed using the maximum likelihood method and the kimura-2 parameter model. Topology was supported by 1000 bootstrap replicates. The sequences reported in this study were identified by the accession numbers. The sequences of isolates reported in this study are indicated by their accession number followed by the country code and the year of isolation. Reference HAV strains are shown in italic. Clusters are designed (a–e). The same color is assigned to the same sequences grouped in the same cluster. (A). Phylogenetic tree comparing 884nt of the VP1 region of 52 strains from this study with 12 reference HAV strains. (B). Phylogenetic tree comparing 319 nt of the VP1-2A region of 77 strains from this study with 12 reference HAV strains.
Figure 1. Phylogenetic trees comparing nucleotide sequences of HAV strains isolated in Tunisia during 2013–2018 with the reference sequences. The tree was performed using the maximum likelihood method and the kimura-2 parameter model. Topology was supported by 1000 bootstrap replicates. The sequences reported in this study were identified by the accession numbers. The sequences of isolates reported in this study are indicated by their accession number followed by the country code and the year of isolation. Reference HAV strains are shown in italic. Clusters are designed (a–e). The same color is assigned to the same sequences grouped in the same cluster. (A). Phylogenetic tree comparing 884nt of the VP1 region of 52 strains from this study with 12 reference HAV strains. (B). Phylogenetic tree comparing 319 nt of the VP1-2A region of 77 strains from this study with 12 reference HAV strains.
Viruses 13 01056 g001
Viruses 13 01056 g002 550
Figure 2. Phylogenetic tree of Tunisian HAV genotype IA isolates, constructed using the maximum likelihood method and the kimura-2 parameter model based on 168nt of the VP1-2A junction. Black squares represent isolates reported in this study (accession numbers MW117967 to MW118018 and MW222090 to MW222114). The sequences of isolates are indicated by their accession number followed by district and year of isolation. Sequences of unknown districts are indicated by their accession number followed by the country code and the year of isolation. Genetic variants are represented in different colors. Topology was supported by 1000 bootstrap replicates. Reference strain HM175 was used as an “outgroup”.
Figure 2. Phylogenetic tree of Tunisian HAV genotype IA isolates, constructed using the maximum likelihood method and the kimura-2 parameter model based on 168nt of the VP1-2A junction. Black squares represent isolates reported in this study (accession numbers MW117967 to MW118018 and MW222090 to MW222114). The sequences of isolates are indicated by their accession number followed by district and year of isolation. Sequences of unknown districts are indicated by their accession number followed by the country code and the year of isolation. Genetic variants are represented in different colors. Topology was supported by 1000 bootstrap replicates. Reference strain HM175 was used as an “outgroup”.
Viruses 13 01056 g002
Viruses 13 01056 g003 550
Figure 3. Phylogenetic tree based on the analysis of the whole VP1 region belonging to genotype IA, constructed using the maximum likelihood method and the kimura-2 parameter model. The tree includes HAV Tunisian 45 sequences reported in this study and 57 previously published sequences from Tunisia and 118 from other countries. Sequences from 21 different countries (Italy, France, Argentina, Chile, Uruguay, Spain, Venezuela, Brazil, USA, Mexico, China, Russia, Japan, South Korea, Mongolia, Hungary, Thailand, Vietnam, Indonesia, Singapore and Yugoslavia) were retrieved from GenBank. The sequences of isolates are indicated by their accession number and the country code. Black triangles represent Tunisian isolates. Sequences from others countries grouped with the Tunisian isolates are represented in different colors (sequences from France in red, from Spain in blue and from Yugoslavia in green). Topology was supported by 1000 bootstrap replicates. Bootstrap values lower than 50 were not indicated. Reference strain CF53 was used as an “outgroup”.
Figure 3. Phylogenetic tree based on the analysis of the whole VP1 region belonging to genotype IA, constructed using the maximum likelihood method and the kimura-2 parameter model. The tree includes HAV Tunisian 45 sequences reported in this study and 57 previously published sequences from Tunisia and 118 from other countries. Sequences from 21 different countries (Italy, France, Argentina, Chile, Uruguay, Spain, Venezuela, Brazil, USA, Mexico, China, Russia, Japan, South Korea, Mongolia, Hungary, Thailand, Vietnam, Indonesia, Singapore and Yugoslavia) were retrieved from GenBank. The sequences of isolates are indicated by their accession number and the country code. Black triangles represent Tunisian isolates. Sequences from others countries grouped with the Tunisian isolates are represented in different colors (sequences from France in red, from Spain in blue and from Yugoslavia in green). Topology was supported by 1000 bootstrap replicates. Bootstrap values lower than 50 were not indicated. Reference strain CF53 was used as an “outgroup”.
Viruses 13 01056 g003
Viruses 13 01056 g004 550
Figure 4. Amino acid replacements localization compared to reference strain and position within the experimentally B/Tcell epitopes regions [12,37,38,40,41,45]. Replacements located in epitope region are indicated in red. Replacements located close to epitope region are indicated in blue.
Figure 4. Amino acid replacements localization compared to reference strain and position within the experimentally B/Tcell epitopes regions [12,37,38,40,41,45]. Replacements located in epitope region are indicated in red. Replacements located close to epitope region are indicated in blue.
Viruses 13 01056 g004
Viruses 13 01056 g005 550
Figure 5. Effect of the deleted region found in the Tunisian variant on the three-dimensional structure of HAV-VP1 protein. (a) HAV capsid protein (PDB-ID: 4QPI). Gray protein corresponds to VP1, blue to VP2 and purple to VP3. (b) Top view Figure 1 protein (PDB-ID: 4QPI-A). (c) The three-dimensional predicted structure of the VP1 protein of the Tunisian variant isolated in 2018 (accession number: MW118004). The homology model was obtained using modeler and as a template HAV-VP1 protein (PDB: 4QPI-A). Circled region indicates the destroyed helix (colored in red) resulted from the deletion of 6 amino acids in the VP1 protein region.
Figure 5. Effect of the deleted region found in the Tunisian variant on the three-dimensional structure of HAV-VP1 protein. (a) HAV capsid protein (PDB-ID: 4QPI). Gray protein corresponds to VP1, blue to VP2 and purple to VP3. (b) Top view Figure 1 protein (PDB-ID: 4QPI-A). (c) The three-dimensional predicted structure of the VP1 protein of the Tunisian variant isolated in 2018 (accession number: MW118004). The homology model was obtained using modeler and as a template HAV-VP1 protein (PDB: 4QPI-A). Circled region indicates the destroyed helix (colored in red) resulted from the deletion of 6 amino acids in the VP1 protein region.
Viruses 13 01056 g005
Table
Table 1. Amino acid replacements identified in Tunisian strains isolated during 2013–2018 and their locations.
Table 1. Amino acid replacements identified in Tunisian strains isolated during 2013–2018 and their locations.
Amino Acid PositionPrototype ResidueSubstituted ResidueType of Amino Acid ReplacementNo. of Isolates
VP1-1ValGlyNonconservative2
VP1-2GlyGluNonconservative1
VP1-3AspAsnConservative1
VP1-6GlyGluNonconservative1
VP1-18ValIleConservative1
VP1-26IleThrNonconservative52
VP1-29ArgLysConservative6
VP1-37ArgLysConservative52
VP1-64ValIleConservative1
VP1-72LysArgConservative1
VP1-99IleThrNonconservative52
VP1-166ValAlaSemiconservative1
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ayouni, K.; Chouikha, A.; Khamessi, O.; Touzi, H.; Hammemi, W.; Triki, H. Evidence of Circulation of Several HAV Genetic Variants and Emergence of Potential Antigenic Variants in an Endemo-Epidemic Country before Vaccine Introduction. Viruses 2021, 13, 1056. https://doi.org/10.3390/v13061056

AMA Style

Ayouni K, Chouikha A, Khamessi O, Touzi H, Hammemi W, Triki H. Evidence of Circulation of Several HAV Genetic Variants and Emergence of Potential Antigenic Variants in an Endemo-Epidemic Country before Vaccine Introduction. Viruses. 2021; 13(6):1056. https://doi.org/10.3390/v13061056

Chicago/Turabian Style

Ayouni, Kaouther, Anissa Chouikha, Oussema Khamessi, Henda Touzi, Walid Hammemi, and Henda Triki. 2021. "Evidence of Circulation of Several HAV Genetic Variants and Emergence of Potential Antigenic Variants in an Endemo-Epidemic Country before Vaccine Introduction" Viruses 13, no. 6: 1056. https://doi.org/10.3390/v13061056

APA Style

Ayouni, K., Chouikha, A., Khamessi, O., Touzi, H., Hammemi, W., & Triki, H. (2021). Evidence of Circulation of Several HAV Genetic Variants and Emergence of Potential Antigenic Variants in an Endemo-Epidemic Country before Vaccine Introduction. Viruses, 13(6), 1056. https://doi.org/10.3390/v13061056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop