- Research article
- Open Access
Hepatitis A virus subgenotyping based on RT-qPCR assays
BMC Microbiologyvolume 14, Article number: 296 (2014)
The hepatitis A virus (HAV) is the most frequent cause of viral hepatitis worldwide and is recognized as one of the most widespread foodborne pathogens. HAV genotypes and subtypes differ in their geographic distribution and the incidence of HAV infection varies considerably among countries, and is particularly high in areas with poor sanitation and hygiene. Phylogenetic analyses are traditionally used in clinical microbiology for tracing the geographic origin of HAV strains. In food microbiology, this approach is complicated by the low contamination levels of food samples. To date, real-time reverse-transcription PCR has been one of the most promising detection methods due to its sensitivity, specificity and ability to deliver quantitative data in food samples, but it does not provide HAV subtyping information.
Six subtype-specific RT-qPCR assays were developed for human HAV. The limit of detection of HAV was 50 genome copies/assay for subtype IIB, 500 genome copies assay for IA, IB, IIA and IIIB and 5000 genome copies/assay for IIIA. The specificity of the assays was evaluated by testing reference isolates and in vitro HAV RNA transcripts. No significant cross reactivity was observed. Subtyping results concordant with sequencing analysis were obtained from 34/35 clinical samples. Co-infection with a minor strain of a different subtype was suggested in 5 cases and a recombinant event in one case.
These RT-qPCR assays may be particularly useful for accurately tracing HAV in low-level contaminated samples such as food matrices but also to allow co-infection identification in human samples.
Hepatitis A virus (HAV) is a small, non-enveloped hepatotropic virus classified into the Hepatovirus genus within the Picornaviridae family. Its genome consists of an approximately 7.5 kilobase positive single-strand RNA comprising a 5’ untranslated region (5’UTR), a single open reading frame (ORF) that encodes both structural and non-structural proteins, and a 3’ UTR with a short poly(A) tail. There is only one serotype of HAV. Genomic characterization of HAV has been carried out mainly by sequencing of strains from different geographic regions of the world. Firstly, using a short fragment of the VP1/2A junction region, strains were classified in to seven genotypes on the basis of >15% nucleotide variation between isolates, and the subgenotypes with >7.5% to <15% nucleotide variation . Then, the complete genomic data indicated that genotypes II and VII should be considered a single genotype, based upon the complete VP1 sequence [2; 3]. So, by sequencing of the VP1/2A junction and the VP1 gene, three genotypes (I, II, III) divided in two subtypes (A and B) have been described for humans and three others (IV, V, VI) for primates -.
HAV infection is the leading worldwide cause of acute viral hepatitis ,. There are an annual estimated of 1.5 million cases of hepatitis A worldwide . Optimal use of vaccination can significantly reduce the hepatitis A disease burden and the World Health Organization position on hepatitis A vaccines depend on the level of endemicity in countries. In highly endemic countries, large-scale vaccination programmes are not recommended. In countries of intermediate endemicity, large-scale childhood vaccination may be considered as a supplement to health education and improved sanitation. And in regions of low endemicity, vaccination against hepatitis A is indicated for individuals with increased risk of contracting the infection such as travelers to areas of intermediate or high endemicity . HAV’s geographical distribution is dependent on socioeconomic development and sanitation levels. In areas with high and very high endemicity (Africa, Middle East, India, Central and South America), where infections are mostly asymptomatic and epidemics are rare, 50% seroprevalence is reached between the ages of 5 and 14 . In areas with moderate endemicity (Eastern Europe and south-eastern Asia), 50% seroprevalence is reached between the ages of 14 and 34 and epidemics can occur within the general population. In areas with low endemicity (North America, Western Europe and Australia), most of the population is still susceptible to HAV, particularly people over 50 years old, and the risk of fulminant hepatitis is higher.
HAV is transmitted mainly by the fecal-oral route, either by person-to-person contact or by ingestion of contaminated water and food, particularly shellfish, soft fruits and raw vegetables -. HAV is stable in the environment and is particularly resistant to disinfectants, heating, pressure and low pH ,. Contamination may occur during growth in the field as well as during processing, storage, distribution or final preparation. In developed countries, low incidence and low vaccine coverage have led to a high proportion of susceptible individuals, which creates a potential for expanded hepatitis A outbreaks when contaminated products are widely distributed .
The development of sensitive, reliable techniques for the detection of HAV in food and water samples contributes to the safety of these products . However, detection of HAV on the basis of its infectivity is complicated by the absence of a reliable cell culture method and the low contamination levels of food samples. HAV detection is currently based on nucleic acid testing methods. The International Organization for Standardization/Technical specification (ISO/TS) 15216 standard was published in the first half of 2013 and will be published as ISO standard methods after validation. These protocols target the 5’UTR which shows the lowest diversity across HAV genotypes -. Currently, HAV genotyping relies on amplification, sequencing and phylogenetic analysis of a portion of the viral genome. However, these techniques are time-consuming and may lack sensitivity, particularly with food samples, where the level of contamination by enteric viruses is often very low. Alternative approaches for HAV genotyping in complex samples (food, environmental) may help to better manage the risk. Indeed, although genotypes I and III are the most frequently reported worldwide, HAV genotypes and HAV strains differ in their geographic distribution ,; strain genotyping can thus give clues to understanding food contamination routes. Currently, very few studies describe alternatives to sequencing for HAV genotyping. In recent years single-nucleotide polymorphism (SNP) genotyping has become an area of intense investigation and a valuable tool for diagnosing various pathologies. Various methods for SNP detection have been reported including real-time PCR performed with primers and a probe spanning the SNP site .
The aim of this study was to develop a new approach for the subgenotyping of human HAV based on six simplex SNP genotyping RT-qPCR assays and to apply this approach to human clinical samples.
Design of HAV subtype RT-qPCR assays
The HAV subtype RT-qPCR assays were designed to give subtype-specific amplification on the basis of SNP differentiating the targeted subtype from the others. In other words, at the SNP position, the same nucleotide was found for all the subtypes except for the subtype of targeted HAV. Consequently, different regions of HAV genome were chosen given their subtype specificity and the absence of major nonspecific homologies on BLAST analysis. Moreover, degenerated bases were used to detect genetic variation within a given subtype (Table 1; Table 2; Figure 1).
Sensitivity of subtype-specific RT-qPCR assays
The sensitivity of the simplex subtype–specific RT-qPCR assays was evaluated with serial 10-fold dilutions of in vitro transcribed RNA for HAV IIA, IIB, IIIA and IIIB and genomic RNA for HAV IA and IB. From 5 × 105 to 5 genome copies/assay for IA, IB, IIA, IIB, and IIIB and from 5 × 107 to 5 × 102 genome copies/assay for IIIA were tested. As shown on Table 3, mean RT-qPCR efficiency, derived from the slope parameters, ranged from 83.9% for IIIA to 109% for IB. R2 values were ≥0.898. The limit of detection (LOD) obtained for IA, IB, IIA and IIIB was 500 genome copies/assay, whereas the LOD of IIB was 50 genome copies/assay and the LOD of IIIA was 5000 genome copies/assay. The LOD of the consensus RT-qPCR assay  was in the same range as that of the subtyping RT-qPCR assays at 500 genome copies/assay.
Specificity of subtype-specific RT-qPCR assays
The specificity of the RT-qPCR assays was assessed by testing HAV RNA of each subtype at a concentration of 5 × 104 genome copies/assay with all the subtype-specific RT-qPCR assays. As shown on Table 4, detection of the specific target was observed for all assays with Ct values comprised between 25 and 37, consistent with assay sensitivity. All but one were entirely specific. The IIA-specific assay occasionally allowed amplification of the IIB target. However, this non-specific IIB amplification was not observed when as much as 5 × 104 genome copies/assay of IIB RNA was tested in the presence of a low concentration of the specific IIA target (50 genome copies/assay) (data not shown).
Fecal and serum samples analysis
Human clinical fecal and serum samples were genotyped by sequencing the VP1/2A region, as described  and provided by the NRC. Then, they were tested with the consensus RT-qPCR assay  and with all the subtype-specific RT-qPCR assays separately (Tables 5 and 6).
Four of the five stool samples and 24 out of the 30 sera were detected by a single subtype-specific assay that provided a subtype result consistent with VP1/2A sequencing. The consensus and specific RT-qPCR assays gave similar results with differences of quantification that did not exceed 1.9 log10 (genome copies/μL or genome copies/g).
In stool sample 1181216151 provided as a IA subtype by the NRC, subtype-specific assays detected both IA and IIA RNA, with a IIA concentration 5.2 log10 lower than that of the IA subtype (Table 5). Similarly, in the 5 sera provided as the IA subtype by the NRC, subtype-specific assays detected both IA and IB RNA, with IB concentrations 0.7 to 2 log10 lower than IA. A single discrepant result was observed for serum sample 1311062298 provided as a IB subtype by VP1/ 2A region sequencing and identified as a IA subtype by the subtype-specific RT-qPCR assays (Table 6).
In conclusion, the subgenotyping RT-qPCR assays allowed detecting 100% (35/35) of the clinical samples for the presence of HAV. In total, 80.0% (28 samples) of the clinical samples were found to correlate with the genotyping by sequencing the VP1/2A region. Furthermore, positivity for more than one genotype identified by sequencing appeared in 17.1% (6 samples) of the clinical samples and a subtype discrepancy in 2.9% (1 sample) of the clinical samples.
Although HAV has been shown to possess a single conserved antigenic neutralization site  leading to a single serotype, HAV strains isolated from different parts of the world have been classified into six genotypes (I to VI), of which genotype I, II, and III can infect humans. Genotype I is the most prevalent worldwide, and subtype IA is more common than IB. The other human genotypes are infrequent. In areas of low endemicity such as the United States and Western Europe, IA dominates but all genotypes and subtypes have been reported ,,. Genotype II isolates were originally identified in France in 1979 and Sierra Leone in 1988  and appear to be limited to West Africa . Genotype III has been reported in many parts of the world  but is prevalent in the Indian subcontinent. An increase in genotype IIIA infections has been reported in Korea, Russia, Estonia and in Japan. Moreover, IIIA and IIIB co-circulate broadly with IA and IB strains .
Phylogenetic analysis is useful to trace back the geographical origin of a given strain and for tracking transmissions of HAV. Accurate typing of HAV from food samples could thus be helpful for transmission investigations. However, HAV typing from food samples by a classical sequencing approach is often impaired by the low contamination levels, and does not give access to potential contamination by several strains. Indeed, implicated items (such as seafood, fruits and salads) in foodborne outbreaks can harbour a heterogeneous HAV population that reflects the diversity of the viral strains circulating at the geographic location of item contamination .
Two commercial quantitative HAV RT-qPCR assays have been described. The detection limit was 2 TCID50 /mL for the Roche kit and 5 TCID50/mL for the Artus kit. Both kits have been found suitable for detection and quantification of HAV but only the Roche kit allowed the differentiation between genotype IA and IB after melting curve analysis . The present study introduces six RT-qPCR-based assays for specific molecular genotyping of hepatitis A virus. To our knowledge, this is the first time that HAV subtyping has been achieved by specific qPCR probes. This subtype identification method can be implemented in diagnostic and research laboratories, avoiding post-PCR analysis and avoiding the problem of low viral loads in food samples.
All subtype assays were found suitable for quantification measurement for comparison with the data obtained with the reference RT-qPCR assay (detecting all genotypes). The minimal variations (around 1log10) observed for the quantification were potentially due to the differences in amplification efficiencies and calibration curves used. Most of the samples were correctly identified with regard to the genotype provided by VP1/2A sequencing. In 6 samples (1 stool and 5 sera), the specific RT-qPCR assay identified a major IA strain, the same one determined by VP1/2A sequencing, also in addition to a second subtype, present in a lower concentration.
The conventional genotyping used as a reference assay is a “golden standard assay”. The design of HAV subgenotyping RT-qPCR assays was based on SNP in the probe associated with degenerated bases in the primers to enhance the specificity. Nevertheless, cross-reactivity could be only definitively excluded with the entire genome sequencing for the tested samples. However, co-circulation of the subgenotypes IA and IIIA has been reported in India  and of IA, IB and IIIA in Korea . Co-circulation of the subgenotypes IA and IB in South Africa, South America, Europe and the US and the existence of recombination events between subgenotypes have also been observed -. Indeed, HAV exploits all known mechanisms of genetic variation to ensure its survival, including mutation and recombination ,. HAV recombination was originally reported in cell culture . Its extent in nature was appreciated only recently ,,,, and it appears that recombination occurs along the entire length of the genome .
The present finding from the stool sample of a patient who had not traveled abroad may be due to a co-infection by IA and IIA subtypes. Indeed, HAV IA is the dominant strain in France but IIA strains have been isolated among French travelers returning from Africa as well as from autochthonous cases . A co-infection rather than an event of recombination is suggested because of the huge difference in the concentration of the subtypes. Regarding these two signals, although non-specific amplification due to a very high viral load cannot be excluded, it should be noted that no IIA amplification was detected from any of the 14 HAV IA serum samples.
The discovery of a major IA signal, combined with a 10- to 100-fold lower IB signal in 5 sera from patients having traveled abroad (at least for three of them) may suggest either an event of recombination or, more likely, a co-infection. For these samples, the genome copy numbers determined by the 5’-UTR assay was not the sum of those determined by subgenotyping RT-qPCR assays together which can be explained by the lack of accurate quantification or by cross reactivity. As conventional Sanger sequencing does not allow accurate identification of multiple species within a sample, the hypotheses could be investigated by cloning and sequencing or by next generation sequencing.
A single sample from a patient contaminated in Morocco provided a discrepant result by specific RT-qPCR and sequencing; this sample may correspond to an IA/IB recombinant in the P1 region of the HAV genome since IA-specific amplification targets the VP4 region (nt 702 to 820) and sequencing targets the VP1/2A region (nt 2870 to 3381). The sequencing of this sample was attempted but has been unsuccessful may be because of the low viral load.
It was concluded that the RT-qPCR assays developed in this study are suitable tools for quantification of HAV and subtype identification. They need to be validated by testing a larger number of clinical, environmental and food samples. Conventional genotyping used as a reference assay is a “golden standard assay”, and the RT-qPCR assays described here could be recommended as an additional test to the conventional genotyping and for use in cases of failure of the conventional typing method. They may be particularly useful for accurately tracing HAV in samples with low-level contamination such as food matrices, but also can provide easy identification of a co-infection in human samples.
The genotype IB HM175/18f strain, clone B (VR-1402) was obtained from the American Type Culture Collection (ATCC). This clone replicates rapidly and has cytopathic effects in cell culture . HAV stock was produced by propagation in foetal rhesus monkey kidney (FRhK-4) cells (ATCC, CRL-1688)  and titrated by plaque assay . Results were expressed in plaque-forming units/mL (PFU/mL) and HAV stock contained 107 PFU/mL. Aliquots of 100 μL were kept frozen at −80°C for later use.
Ethics statement: Hepatitis A virus infection is a notifiable disease in France. The current system of mandatory reporting was approved by the Commission Nationale de l’Informatique et des Libertés (deliberation n° 02–082, November 19 2002). Patients receive oral and written information on the finality of the notification and on the modalities of information recording. This information is available on line on the web site of the Institut de Veille Sanitaire (IVS) at http://www.invs.sante.fr/content/download/6498/42945/version/2/file/fiche_info_patient.pdf.
All clinical and biological parameters are treated anonymously. The virological surveillance of strain diversity is performed on stored samples obtained for hepatitis A diagnosis (no need for any additional blood draw). Diagnostic laboratories are asked to contribute to HAV strains surveillance by sending samples to the National Reference Centre (NRC) for HAV. All data and samples are anonymously collected and analyzed. The study was conducted in accordance with the ethics principles of the Declaration of Helsinki.
HAV genotyping from stools and serum samples collected by the French NRC for Hepatitis A was determined by sequencing of the VP1/2A junction region as previously described . Stool samples were suspended in 10 mM Phosphate Buffered Saline (PBS), pH 7.4, to obtain a final 10% suspension (w/v), vortexed and centrifuged at 3000 g for 30 min at 4°C. Aliquots of 100 μL supernatant were kept frozen at –80°C for later use. Serum samples were kept frozen at –80°C until later use.
Viral RNA extraction
Aliquots of frozen fecal samples or viral stocks were supplemented with NucliSens® easyMAG™ lysis buffer (BioMérieux, Marcy l’Etoile, France) up to 3 mL and subjected to the NucliSens® easyMAG™ platform (Biomérieux) for total nucleic acid extraction by the “off board Specific A protocol” according to manufacturer’s instructions. Nucleic acids were finally eluted in 70 μL of elution buffer and stored at −80°C.
Two hundred μL of frozen sera samples were subjected to the NucliSens® easyMAG™ platform (Biomérieux) for total nucleic acid extraction by the “Specific B protocol” according to manufacturer’s instructions. Nucleic acids were then eluted in 50 μL of elution buffer and stored at -80°C.
HAV RNA in vitro transcripts
The cDNA corresponding to nucleotides 39–518 (5’UTR) of the IB genomic sequence (M59808.1) was cloned into the pGEM-T Easy vector (Promega, Charbonnières-les-Bains, France) and propagated in E. coli One Shot® TOP10F’ (Life technologies, Saint Aubin, France). High quality DNA plasmid containing HAV regions (p-HAV5) was purified using the Qiagen Plasmid midi kit (Qiagen, Courtaboeuf, France) according to the manufacturer’s protocol.
HAV cDNA of genotypes IIA, IIB, IIIA or IIIB corresponding respectively to the 588–3183, 587–3183, 618–3210 and 618–3210 positions of the genomic sequence (AY644676.1, AY644670.1, AB279732.1, AB279735.1) were cloned into the pBluescriptIISK + vector by Genecust (Dudelange, Luxembourg). All recombinant plasmids were purified by Genecust and used to produce RNA transcripts. HAV IIA, IIB, and IIIB DNA plasmids (0,5 μg) were digested with Hind III (Life technologies) and HAV 5’UTR and HAV IIIA DNA plasmids were digested with Spe I (Life technologies). Digested plasmids were transcribed by using the MEGAscript® kit (Life technologies) according to the manufacturer’s protocol. Synthesized RNA was treated twice with Turbo™ DNase (Life technologies) according to the manufacturer’s protocol in order to remove the DNA template following transcription, and purified by using the MEGAclear kit (Life technologies) according to manufacturer’s instructions. The synthesized RNA was confirmed with RT-qPCR and quantified by measuring absorbance at 260/280 nm with a Nanodrop ND-100 (Thermoscientific, France) and the free software available on the “http://endmemo.com/bio/dnacopynum.php” website. RNA stocks were diluted to contain 109 copies/μL and aliquoted and stored at - 80°C.
Titers of the clarified fecal suspensions, serum samples and HM175/18f supernatants were obtained by RT-qPCR targeting the 5’UTR (see below), using a standard curve derived from ten-fold dilutions of the 5’UTR transcript RNA from p-HAV5. Titer was expressed in genome copies.
Primers and probes
The RT-qPCR assay targeting the 5’UTR described by Costafreda et al.  was used to detect all HAV genotypes. This RT-qPCR assay is referred to as “consensus RT-qPCR”. Primers and probe sets were designed by using Beacon Designer software (Bio-Rad, Marnes-la-Coquette, France) to give subtype-specific amplification on the basis of single nucleotide polymorphisms. To identify genotype-specific conserved regions of HAV, complete sequences available from GenBank (NCBI) were aligned (Table 1) with MUSCLE software  and multiple alignment was visualized with JALVIEW software (version 2.8) . Hydrolysis probes were labeled at the 5’ end with 6-carboxyfluorescein (FAM) and at the 3’ end with black hole quencher 1 (BHQ1) (Table 2 and Figure 1).
Primers and probes were purchased from Life Technologies or Eurofins MWG Operon (Les Ulis, France).
Quantitative one-step RT-PCR for detection of HAV was carried out on a CFX96™ real-time PCR detection system from Bio-Rad. Reactions were performed in a 15 μL reaction mixture containing 1X of RNA Ultrasense™ master mix and 0.63 μL of RNA Ultrasense™ enzyme mix, which are components of RNA Ultrasense™ One-Step Quantitative RT-PCR System (Life technologies), 2 U RNAse inhibitor (Life technologies), 5 μg of bovine serum albumin (Life Technologies), 500 nM of forward primer, 900 nM of reverse primer, 250 nM of probe, and 5 μL of sample. A negative control containing all the reagents except the RNA template was included in each set of reactions. The one-step RT-qPCR program involved 60 min reverse transcription of RNA at 55°C, followed by a 5 min denaturation step at 95°C, and finally 45 cycles of 15 s at 95°C, 1 min at 56°C and 1 min at 65°C. Fluorescence was automatically recorded by the instrument at the end of the elongation steps (1 minute at 65°C) for each amplification cycle. All samples were characterized by a corresponding cycle threshold (Ct) value. Negative samples gave no Ct value. For each specific RT-qPCR assay, a standard curve was generated using 10-fold dilutions of titered RNA corresponding to each subtype. For the consensus RT-qPCR assay, a standard curve was generated using 10-fold dilutions of titered RNA transcripts from p-HAV5. The slopes (S) of the regression lines were used to calculate the amplification efficiency (E) of the RT-qPCR reactions, according to the formula E =10|-1/s| -1 .
Assay performance assessment
Genotype IB HAV RNA obtained from HM175/18f, genotype IA RNA obtained from a fecal sample (stool number 128061099) and genotype IIA, IIIA, IIB, IIIB RNA transcripts were used to determine the sensitivity and the specificity of the subgenotyping RT-qPCR assays. All samples were analyzed in duplicate in three different experiments resulting from 6 Ct values.
CC, AF, CM performed these experiments. All authors wrote, read and approved the final manuscript.
American type culture collection
Basic local alignment search tool
European committee for standardisation
Hepatitis A virus
International organization for standardization/technical specification
Limit of detection
Open reading frame
National center for biotechnology information
National reference centre
Quantitative reverse transcriptase PCR
Single nucleotide polymorphism
Robertson BH, Jansen RW, Khanna B, Totsuka A, Nainan OV, Siegl G, Widell A, Margolis HS, Isomura S, Ito K: Genetic relatedness of hepatitis A virus strains recovered from different geographical regions. J Gen Virol. 1992, 73: 1365-1377. 10.1099/0022-1317-73-6-1365.
Lu L, Ching KZ, de Paula VS, Nakano T, Siegl G, Weitz M, Robertson BH: Characterization of the complete genomic sequence of genotype II hepatitis A virus (CF53/Berne isolate). J Gen Virol. 2004, 85: 2943-2952. 10.1099/vir.0.80304-0.
Costa-Mattioli M, Cristina J, Romero H, Perez-Bercoff R, Casane D, Colina R, Garcia L, Vega I, Glikman G, Romanowsky V, Castello A, Nicand E, Bassin M, Billaudel S, Ferre V: Molecular evolution of hepatitis A virus: a new classification based on the complete VP1 protein. J Virol. 2002, 76: 9515-9525. 10.1128/JVI.76.18.9516-9525.2002.
Koopmans M, Duizer E: Foodborne viruses: an emerging problem. Int J Food Microbiol. 2004, 90: 23-41. 10.1016/S0168-1605(03)00169-7.
Vaughan G, Goncalves Rossi LM, Forbi JC, de Paula VS, Purdy MA, Xia G, Khudyakov YE: Hepatitis A virus: host interactions, molecular epidemiology and evolution. Infect Genet Evol. 2014, 21C: 227-243. 10.1016/j.meegid.2013.10.023.
Franco E, Meleleo C, Serino L, Sorbara D, Zaratti L: Hepatitis A: epidemiology and prevention in developing countries. World J Hepatol. 2012, 4: 68-73. 10.4254/wjh.v4.i3.68.
Hepatitis A vaccines. Wkly Epidemiol Rec. 2000, 75: 38-44.
Mohd Hanafiah K, Jacobsen KH, Wiersma ST: Challenges to mapping the health risk of hepatitis A virus infection. Int J Health Geogr. 2011, 10: 57-10.1186/1476-072X-10-57.
Anonymous: Hepatitis A outbreak in Australia [Internet]. National Travel Health Network and Centre; 2009 November 13 [cited 2010 March 25]. Available from:., [http://www.nathnac.org/pro/clinical_updates/hepatitisaoutbreakaustralia_131109healthprofessionals.htm]
Beuchat LR: Vectors and conditions for preharvest contamination of fruits and vegetables capable of causing enteric diseases. Brit Food J. 2006, 108: 38-53. 10.1108/00070700610637625.
Gallot C, Grout L, Roque-Afonso A-M, Couturier E, Carrillo-Santisteve P, Pouey J, Letort M-J, Hoppe S, Capdepon P, Saint-Martin S, De Valk H, Vaillant V: Hepatitis A associated with semidried tomatoes, France, 2010. Emerg Infect Dis. 2011, 17: 566-567. 10.3201/eid1703.101479.
Hernández F, Monge R, Jiménez C, Taylor L: Rotavirus and hepatitis A virus in market lettuce (Latuca sativa) in Costa Rica. Int J Food Microbiol. 1997, 37: 221-223. 10.1016/S0168-1605(97)00058-5.
Lacey S, Tracy L, Hammond R, Revill P, Donnan E, Lalor K, Bowden S: Hepatitis A Virus Outbreak: Molecular Epidemiology Indicates a Common Source. 2009, The 5th TEPHINET Southeast Asia and Western Pacific Bi-Regional Scientific Conference, Seoul, Korea
Petrignani M, Verhoef L, van Hunen R, Swaan C, van Steenbergen J, Boxman I, Ober HJ, Vennema H, Koopmans M: A possible foodborne outbreak of hepatitis A in the Netherlands, January–February 2010. Euro Surveill. 2010, 15 (11): 19512-
Petrignani M, Harms M, Verhoef L, van Hunen R, Swaan C, van Steenbergen J, Boxman I, Perani Sala R, Ober H, Vennema H, Koopmans M, van Pelt W: Update: a food-borne outbreak of hepatitis A in the Netherlands related to semi-dried tomatoes in oil, January–February 2010. Euro Surveill. 2010, 15 (20): 19572-
Rosenblum LS, Mirkin IR, Allen DT, Safford S, Hadler SC: A multifocal outbreak of hepatitis A traced to commercially distributed lettuce. Am J Public Health. 1990, 80: 1075-1079. 10.2105/AJPH.80.9.1075.
Koopmans M, von Bonsdorff CH, Vinjé J, de Medici D, Monroe S: Foodborne viruses. FEMS Microbiol Rev. 2002, 26: 187-205.
Sánchez G, Bosch A, Pintó RM: Hepatitis A virus detection in food: current and future prospects. Lett Appl Microbiol. 2007, 45: 1-5. 10.1111/j.1472-765X.2007.02140.x.
Costafreda MI, Bosch A, Pinto RM: Development, evaluation and standardization of a real time TaqMan reverse transcription-PCR assay for quantification of hepatitis A virus in clinical and shellfish samples. Appl Environ Microbiol. 2006, 72: 3846-3855. 10.1128/AEM.02660-05.
Microbiology of Food and Animal Feed - Horizontal Method for Determination of Hepatitis A Virus and Norovirus in Food Using Real-Time RT-PCR - Part 1: Method for Quantification. 2013, International Organization for Standardization, Geneva, Switzerland
Microbiology of Food and Animal Feed - Horizontal Method for Determination of Hepatitis A Virus and Norovirus in Food Using Real-Time RT-PCR - Part 2: Method for Qualitative Detection. 2013, International Organization for Standardization, Geneva, Switzerland
Lees D: International standardization of a method for detection of human pathogenic viruses in molluscan shellfish. Food Environ Virol. 2010, 2: 146-155. 10.1007/s12560-010-9042-5.
Nainan OV, Xia G, Vaughan G, Margolis HS: Diagnosis of hepatitis a virus infection: a molecular approach. Clin Microbiol. 2006, 19: 63-79. 10.1128/CMR.19.1.63-79.2006.
Tjon G, Xiridou M, Coutinho R, Bruisten S: Different transmission patterns of hepatitis A virus for two main risk groups as evidenced by molecular cluster analysis. J Med Virol. 2007, 79: 488-494. 10.1002/jmv.20821.
Woodward J: Bi-allelic SNP genotyping using the TaqMan(®) assay. Methods Mol Biol. 2014, 1145: 67-74. 10.1007/978-1-4939-0446-4_6.
Schwarz NG, Revillion M, Roque-Afonso AM, Dussaix E, Giraud M, Liberpre C, Couturier E, Delarocque Astagneau E: A food-borne outbreak of hepatitis A virus (HAV) infection in a secondary school in Upper Normandy, France, in November 2006. Euro Surveill. 2008, 13 (22): 18885-
Stapleton JT, Lemon SM: Neutralization escape mutants define a dominant immunogenic neutralization site on hepatitis A virus. J Virol. 1987, 61: 491-498.
Costa-Mattioli M: Genetic variability of hepatitis A virus. J Gen Virol. 2003, 84: 3191-3201. 10.1099/vir.0.19532-0.
Tjon GM, Wijkmans CJ, Coutinho RA, Koek AG, van den Hoek JA, Leenders AC, Schneeberger PM, Bruisten SM: Molecular epidemiology of hepatitis A in Noord-Brabant, The Netherlands. J Clin Virol. 2005, 32: 128-136. 10.1016/j.jcv.2004.03.008.
Desbois D, Couturier E, Mackiewicz V, Graube A, Letort MJ, Dussaix E, Roque-Afonso AM: Epidemiology and genetic characterization of hepatitis A virus genotype IIA. J Clin Microbiol. 2010, 48: 3306-3315. 10.1128/JCM.00667-10.
Vaughan G, Xia G, Forbi JC, Purdy MA, Rossi LM, Spradling PR, Khudyakov YE: Genetic relatedness among hepatitis A virus strains associated with food-borne outbreaks. PLoS One. 2013, 8: e74546-10.1371/journal.pone.0074546.
Sánchez G, Populaire S, Butot S, Putallaz T, Joosten H: Detection and differentiation of human hepatitis A strains by commercial quantitative real-time RT-PCR tests. J Virol Methods. 2006, 132: 160-165. 10.1016/j.jviromet.2005.10.010.
Hussain Z, Das BC, Husain SA, Asim M, Chattopadhyay S, Malik A, Poovorawan Y, Theamboonlers A, Kar P: Hepatitis A viral genotypes and clinical relevance: clinical and molecular characterization of hepatitis A virus isolates from northern India. Hepatol Res. 2005, 32: 16-24. 10.1016/j.hepres.2005.01.014.
Song HU, Hwang SG, Kwon CI, Lee JE, Ko KH, Hong SP, Park PW, Rim KS: Molecular epidemiology of hepatitis A virus in the South-East area of Gyeonggi-do in Korea. Yonsei Me J. 2009, 50: 512-516. 10.3349/ymj.2009.50.4.512.
Aguirre S, Malirat V, Scodeller E, Mattion N: First full-length genomic sequence of a hepatitis A virus isolated in Argentina shows recombination between subgenotypes IA and IB. Virus Res. 2011, 155: 316-324. 10.1016/j.virusres.2010.10.030.
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. 10.1007/s11262-009-0431-5.
Villar LM, Morais LM, Aloise R, Melo MM, Calado IA, Lampe E, Gaspar AM: Co-circulation of genotypes IA and IB of hepatitis A virus in Northeast Brazil. Braz J Med Biol Res. 2006, 39: 873-881. 10.1590/S0100-879X2006000700004.
Belalov IS, Isaeva OV, Lukashev AN: Recombination in hepatitis A virus: evidence for reproductive isolation of genotypes. J Gen Virol. 2011, 92: 860-872. 10.1099/vir.0.024786-0.
Cristina J, Costa-Mattioli M: Genetic variability and molecular evolution of hepatitis A virus. Virus Res. 2007, 127: 151-157. 10.1016/j.virusres.2007.01.005.
Lemon SM, Murphy PC, Shields PA, Ping LH, Feinstone SM, Cromeans T, Jansen RW: Antigenic and genetic variation in cytopathic hepatitis A virus variants arising during persistent infection: evidence for genetic recombination. J Virol. 1991, 65: 2056-2065.
Costa-Mattioli M, Ferré V, Casane D, Perez-Bercoff R, Coste-Burel M, Imbert-Marcille BM, Andre EC, Bressollette-Bodin C, Billaudel S, Cristina J: Evidence of recombination in natural populations of hepatitis A virus. Virology. 2003, 20: 51-59. 10.1016/S0042-6822(03)00109-0.
Cromeans T, Sobsey MD, Fields HA: Development of a plaque assay for a cytopathic, rapidly replicating isolate of hepatitis A virus. J Med Virol. 1987, 22: 45-56. 10.1002/jmv.1890220107.
Dubois E, Hennechart C, Deboosère N, Merle G, Legeay O, Burger C, Le Calvé M, Lombard B, Ferré V, Traoré O: Intra-laboratory validation of a concentration method adapted for the enumeration of infectious F-specific RNA coliphage, enterovirus, and hepatitis A virus from inoculated leaves of salad vegetables. Int J Food Microbiol. 2006, 108: 164-171. 10.1016/j.ijfoodmicro.2005.11.007.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.
Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview version 2: a multiple sequence alignment and analysis workbench. Bioinformatics. 2009, 25: 1189-1191. 10.1093/bioinformatics/btp033.
Tichopad A, Dilger M, Schwarz G, Plaffl MW: Standardized determination of real-time PCR efficiency from a single reaction set-up. Nucleic Acids Res. 2003, 31 (20): e122-10.1093/nar/gng122. Erratum in: Nucleic Acids Res 31(22) 6688
This work is part of the thesis by Coralie Coudray-Meunier, a PhD student who received financial support from ANSES.
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.