Open Access

Hepatitis A virus subgenotyping based on RT-qPCR assays

  • Coralie Coudray-Meunier1,
  • Audrey Fraisse1,
  • Camélia Mokhtari2,
  • Sandra Martin-Latil1,
  • Anne-Marie Roque-Afonso2, 3, 4 and
  • Sylvie Perelle1Email author
BMC Microbiology201414:296

https://doi.org/10.1186/s12866-014-0296-1

Received: 7 April 2014

Accepted: 13 November 2014

Published: 25 November 2014

Abstract

Background

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.

Results

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.

Conclusions

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.

Keywords

Hepatitis A virus RT-qPCR assays Genotyping

Background

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 [1]. 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 [1]-[3].

HAV infection is the leading worldwide cause of acute viral hepatitis [4],[5]. There are an annual estimated of 1.5 million cases of hepatitis A worldwide [6]. 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 [7]. 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 [8]. 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 [9]-[16]. HAV is stable in the environment and is particularly resistant to disinfectants, heating, pressure and low pH [4],[17]. 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 [8].

The development of sensitive, reliable techniques for the detection of HAV in food and water samples contributes to the safety of these products [18]. 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 [19]-[22]. 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 [23],[24]; 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 [25].

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.

Results

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).
Table 1

GenBank accession numbers for HAV isolates

Genotype

Reference strain

GenBank accession number

HAV IA

AB020564.1

EU526088; EU526089; EU131373; AB020567; X75216; EF406357; AB623053; X83302; X75214; AB020569; AB618531; AB020565; K02990; AB618529; AF485328; EF207320; HM769724; AB020564; EU251188; X75215; AB020568; HQ437707; AF512536; AB020566; AF357222

HAV IB

M14707

HQ246217; NC_001489; M14707; HV192265; FB746524; M59810; EF406361; EF406359; EF405360; DQ646426; EF406363; EF406362; AF268396; M59809; M16632; M59808; EF406358; AF314208; M20273

HAV IIA

AY644676.1

AY574059; AY644676; GU390574; GU390572; GU390577; GU390576

HAV IIB

AY644670.1

AY644670; Z77248; Z77247; Z77245; Z77244; Z77243; Z77246

HAV IIIA

AB279732.1

AB279732; FJ360735; EU011791; FJ360730; FJ360733; DQ991030; AB279734; DQ991029; FJ360732; FJ360734; AB279733; FJ360731

HAV IIIB

AB279735.1

AB279735; AB258387; AB425339; AB300205

Nucleic acid sequences were used to design primers and probes sets.

Table 2

Sequences of primers and probes in the direction 5’-3’

Genotype

Reference strain

 

Sequence

Position

HAV IA

AB020564.1

Forward primer

GCA TTT AGG TTT TTC CTC ATT

702-722

Reverse primer

TCA ACK GAC TGA ATC ATT

837-820

Probe

TCC AAA CAA GGA ATT TTC CAG AC

743-765

HAV IB

M14707

Forward primer

AAG CTT ATT GTG TAY TGT TAT

2070-2090

Reverse primer

CAG AAT CAT CTC CAA CYT

2223-2208

Probe

TTC TCC TTC TAA C GT TGC TTC CCA

2103-2125

HAV IIA

AY644676.1

Forward primer

ACY ATG ATG AGC AGA ATT G

2978-2996

Reverse primer

GCA TAT TTT AAT CTC TGC TT

3129-3110

Probe

AGA CCT GGA ATC G TC AGT GGA TGA

3004-3027

HAV IIB

AY644670.1

Forward primer

GGA GAT TTG AAA GTC ATA TTG

3047-3067

Reverse primer

TTC CTG GGC ATA CTT TAG

3135-3118

Probe

AGT CTT AAT TCT TTG TAT GGT TTTC

3078-3094

HAV IIIA

AB279732.1

Forward primer

TCC CTT GGA TTT GAC AAT

2605-2622

Reverse primer

RGT ATT RAA CCT AAC AGC

2763-2747

Probe

AAT TAT AAC TGG G GC TAC TGA TGT T

2623-2647

HAV IIIB

AB279735.1

Forward primer

AAT CCG ATG CTT CTC AAG

1560-1577

Reverse primer

GCC TTC CTG AAT GGT ATT

1841-1824

Probe

AAA ATT ACA CAC TTY ACA ACY TGG A

1589-1613

HAV 5’UTR

M14707

Forward primer*

TCA CCG CCG TTT GCC TAG

68-85

Reverse primer*

GGA GAG CCC TGG AAG AAA G

241-223

Probe*

CCT GAA CCT GCA GGA ATT AA

169-150

The specific genotype SNP is in bold. Probes are FAM-BHQ except HAV 5’UTR which is FAM-MGB. *: Costafreda et al. [19].

Figure 1

HAV genome regions targeted for genotyping. The different genomic regions used to identify each HAV genotype are represented below the HAV genome scheme.

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 [19] was in the same range as that of the subtyping RT-qPCR assays at 500 genome copies/assay.
Table 3

Characteristics of RT-qPCR standard curves

Genome copies/RT-qPCR assay

Mean Ct values +/− SD

HAV IA

HAV IB

HAV IIA

HAV IIB

HAV IIIA

HAV IIIB

HAV 5'UTR

5 x 107

/

/

/

/

26,56 ± 0,89

/

/

5 x 106

/

/

/

/

29,49 ± 0,88

/

/

5 x 105

26.28 ± 0.59

26.54 ± 0.95

26.64 ± 0.60

22.25 ± 0.45

32.97 ± 0.94

24.50 ± 0.50

27.07 ± 1.00

5 x 104

30.01 ± 0.58

29.60 ± 0115

30.26 ± 0.52

25.53 ± 0.35

36.90 ± 1.69

27.65 ± 0.48

30.76 ± 0.74

5 x 103

33.34 ± 0.66

32.66 ± 1.47

33.82 ± 0.56

28.62 ± 0.66

41.34 ± 0.69

31.58 ± 0.47

33.32 ± 0.79

5 x 102

36.28 ± 0.69

36.16 ± 1.48

36.59 ± 0.79

32.00 ± 0.67

nd

35.64 ± 1.48

37.16 ± 1.83

5 x 101

nd

nd

nd

35.79 ± 1.09

/

nd

nd

5 x 100

nd

nd

nd

nd

/

nd

nd

E

99.0%

109.0%

95.8%

100.8%

83.9%

86.5%

102.2%

R2

0.974

0.898

0.973

0.982

0.952

0.966

0.916

Parameters of RT-qPCR amplification curves obtained for HAV detection by the RT-qPCR reference method and HAV subgenotyping by RT-qPCR assays. The limit of detection (LOD) has been defined as the lowest amount of HAV detected in the three experiments and is shown in bold. nd: not detected. / : not analyzed.

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).
Table 4

Specificity of subgenotyping RT-qPCR assays

 

Set primers/probe

 

HAV IA

HAV IB

HAV IIA

HAV IIB

HAV IIIA

HAV IIIB

Sample

HAV IA

30.01 ± 0.58 (6/6)

-

-

-

-

-

HAV IB

-

29.60 ± 1.15 (6/6)

-

-

-

-

HAV IIA

-

-

30.57 ± 0.87 (6/6)

-

-

-

HAV IIB

-

-

35.02 ± 1.31 (4/6)

25.53 ± 0.35 (6/6)

-

-

HAV IIIA

-

-

-

-

36.90 ± 1.69 (6/6)

-

HAV IIIB

-

-

-

-

-

27.65 ± 0.48 (6/6)

Six subtyping RT-qPCR assays were tested with 5 x 104 genome copies/assay for all subtypes of HAV in duplicate in three different experiments. Results are expressed as means cycle threshold (Ct) values ± standard deviations (SD). The number of positive Ct values is given in parentheses.

Fecal and serum samples analysis

Human clinical fecal and serum samples were genotyped by sequencing the VP1/2A region, as described [26] and provided by the NRC. Then, they were tested with the consensus RT-qPCR assay [19] and with all the subtype-specific RT-qPCR assays separately (Tables 5 and 6).
Table 5

Stool samples analysis

Stool

Age/Sex

Transaminases IU/mL

Travel

Genotype by sequencing

HAV genome copies/g of stool

Difference of quantification between 5’UTR and subtype RT-qPCR assays (log10)

Difference of quantification between IA and IIA subtype RT-qPCR assays (log10)

5’UTR

IA

IB

IIA

IIB

IIIA

IIIB

0780627147

21-30/M

NC

none

IA

2.30 x107

1.93 x108

-

-

-

-

-

IA : −0.92

-

1280210015

51-60/F

5680

Senegal

IB

1.50 x1010

-

7.19 x1010

-

-

-

-

IB : −0.68

-

1280514230

51-60/F

7191

Benin

IB

6.85 x1010

-

6.37 x1010

-

-

-

-

IB : +0.03

-

1181216151

51-60/F

3623

none

IA

7.75 x1011

5.75 x1011

-

3.85 x106

-

-

-

IA : +0.13 ; IIA : +5.30

−5.17

078014121

51-60/F

NC

Morocco

IB

1.12 x107

-

3.73 x105

-

-

-

-

IB : +1.48

-

Each sample was tested with the reference RT-qPCR assay targeting the 5’UTR of HAV and the 6 genotype-specific RT-qPCR assays. The subtyping results were compared with those obtained with sequencing by the NRC. Concentrations are given in genome copies per gram of stool. NC = Not communicated. The difference of quantification between 5’UTR and subtype RT-qPCR assays is calculated by the formula: (log10 (genomes copies determined by reference RT-qPCR/genomes copies determined by subgenotyping RT-qPCR assays)). The difference of quantification between IA and IIA subtypes by RT-qPCR assays is calculated by the formula: (log10 (genomes copies determined by IA RT-qPCR/genomes copies determined by IIA RT-qPCR assays)).

Table 6

Serum samples analysis

Serum

Age/Sex

Transaminases IU/mL

Travel

Genotype by sequencing

HAV genome copies/μl of serum

Difference of quantification between 5’UTR and subtype assays (log10)

Difference of quantification between IA and IB subtype assays (log10)

5’UTR

IA

IB

IIA

IIB

IIIA

IIIB

1310074855

61-70/F

3935

Ethiopia

IB

1.50 x108

-

2.87 x108

-

-

-

-

IB : –0.28

-

1310016965

11-20/M

400

NC

IB

3.13 x105

-

4.18 x103

-

-

-

-

IB : +1.87

-

1311012387

1-10/F

NC

NC

IA

6.37 x104

9.26 x104

-

-

-

-

-

IA : −0.16

-

1311018436

1-10/M

NC

NC

IA

2.86 x104

3.60 x104

-

-

-

-

-

IA : −0.10

-

1311072234

31-40/F

2018

NC

IA

7.36 x104

7.31 x104

8.58 x103

-

-

-

-

IA : 0 ; IB : +0.93

–0.93

1310064958

61-70/M

2314

Guadeloupe

IA

4.73 x105

5.50 x105

-

-

-

-

-

IA : −0.07

-

1311009716

51-60/M

2339

NC

IB

1.70 x106

-

2.30 x106

-

-

-

-

IB : −0.13

-

1311071503

1-10/M

358

NC

IA

2.72 x105

2.85 x106

3.22 x104

-

-

-

-

IA : −1.02 ; IB : +0.93

−1.95

1310024892

1-10/M

NC

Morocco

IA

6.14 x106

6.47 x106

-

-

-

-

-

IA : −0.02

-

1310066012

1-10/F

1113

Morocco

IA

2.04 x103

2.48 x103

-

-

-

-

-

IA : −0.09

-

1309036458

11-20/M

3393

Cameroun

IIA

3.42 x105

-

-

4.66 x106

-

-

-

IIA : −1.13

-

1309064503

1-10/F

2470

Algeria

IA

3.27 x103

3.99 x103

-

-

-

-

-

IA : −0.09

-

1309044888

11-20/M

1352

Guinea

IB

4.10 x105

-

7.00 x105

-

-

-

-

IB : −0.23

-

1310077770

31-40/M

NC

Ethiopia

IB

1.48 x106

-

4.35 x106

-

-

-

-

IB : −0.47

-

1310044717

51-60/F

6500

NC

IB

7.16 x105

-

1.32 x106

-

-

-

-

IB : −0.26

-

1311011402

1-10/M

993

Ethiopia

IB

7.89 x105

-

1.38 x106

-

-

-

-

IB : −0.24

-

1311018712

41-50/M

NC

Ethiopia

IB

4.63 x103

-

1.22 x104

-

-

-

-

IB : −0.42

-

1310005428

21-30/M

3400

Morocco

IA

5.87 x103

1.52 x104

1.37 x102

-

-

-

-

IA : −0.41 ; IB : +1.63

–2.05

1311011353

41-50/F

1736

Ethiopie

IB

4.43 x104

-

9.28 x104

-

-

-

-

IB : −0.32

-

1310023611

21-30/F

1388

Tunisia

IA

3.84 x104

2.09 x104

4.57 x103

-

-

-

-

IA : +0.26 ; IB : +0.92

−0.66

1380219001

51-60/M

4000

Madagascar

IIIA

4.73 x105

-

-

-

-

3.10 x104

-

IIIA : +1.18

-

1311018504

61-70/F

1334

NC

IB

6.96 x104

-

2.09 x105

-

-

-

-

IB : −0.48

-

1309047363

11-20/F

2000

NC

IB

9.76 x104

-

1.30 x105

-

-

-

-

IB : −0.12

-

1310011213

1-10/M

1831

Algeria

IA

7.89 x104

6.68 x104

4.77 x103

-

-

-

-

IA : +0.07 ; IB : +1.22

−1.15

1309064723

11-20/F

1494

Morocco

IA

3.22 x104

2.76 x104

-

-

-

-

-

IA : +0.07

-

1309039311

1-10/M

NC

Guinea

IB

1.25 x105

-

5.35 x105

-

-

-

-

IB : −0.63

-

1310053557

21-30/M

269

NC

IA

9.38 x101

6.33 x102

-

-

-

-

-

IA : −0.83

-

1310078280

1-10/F

857

Tunisia

IA

6.06 x104

3.41 x104

-

-

-

-

-

IA : +0.25

-

1380311211

61-70/M

1269

Madagascar

IIIA

1.48 x104

-

-

-

-

4.25 x104

 

IIIA : −0.46

-

1311062298

11-20/F

940

Morocco

IB

2.95 x102

2.64 x102

-

-

-

-

-

IA : +0.05

-

Each sample was tested with the reference RT-qPCR assay targeting the 5’UTR of HAV and the 6 genotype-specific assays. The subtyping results were compared with those obtained with sequencing by the NRC. Concentrations are given in genome copies per μl of serum. NC = Not communicated. The difference of quantification between 5’UTR and subtype RT-qPCR assays is calculated by the formula: (log10 (genomes copies determined by reference RT-qPCR/genomes copies determined by subgenotyping RT-qPCR assays)). The difference of quantification between IA and IB subtypes by RT-qPCR assays is calculated by the formula: (log10 (genomes copies determined by IA RT-qPCR/genomes copies determined by IB RT-qPCR assays)).

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.

Discussion

Although HAV has been shown to possess a single conserved antigenic neutralization site [27] 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 [23],[28],[29]. Genotype II isolates were originally identified in France in 1979 and Sierra Leone in 1988 [1] and appear to be limited to West Africa [30]. Genotype III has been reported in many parts of the world [28] 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 [5].

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 [31].

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 [32]. 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 [33] and of IA, IB and IIIA in Korea [34]. 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 [35]-[37]. Indeed, HAV exploits all known mechanisms of genetic variation to ensure its survival, including mutation and recombination [38],[39]. HAV recombination was originally reported in cell culture [40]. Its extent in nature was appreciated only recently [35],[36],[38],[39],[41] and it appears that recombination occurs along the entire length of the genome [38].

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 [30]. 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.

Conclusions

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.

Methods

Viral isolates

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 [40]. HAV stock was produced by propagation in foetal rhesus monkey kidney (FRhK-4) cells (ATCC, CRL-1688) [42] and titrated by plaque assay [43]. 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.

Clinical samples

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 [26]. 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. [19] 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 [44] and multiple alignment was visualized with JALVIEW software (version 2.8) [45]. 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).

RT-qPCR conditions

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 [46].

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.

Authors’ contributions

CC, AF, CM performed these experiments. All authors wrote, read and approved the final manuscript.

Abbreviations

ATCC: 

American type culture collection

BLAST: 

Basic local alignment search tool

CEN: 

European committee for standardisation

cDNA: 

Complementary DNA

DNA: 

Deoxyribonucleic acid

HAV: 

Hepatitis A virus

ISO/TS: 

International organization for standardization/technical specification

LOD: 

Limit of detection

ORF: 

Open reading frame

NCBI: 

National center for biotechnology information

NRC: 

National reference centre

PBS: 

Phosphate-buffered saline

PFU: 

Plaque-forming units

RNA: 

Ribonucleic acid

RT-qPCR: 

Quantitative reverse transcriptase PCR

SNP: 

Single nucleotide polymorphism

UTR: 

Untranslated region

Declarations

Acknowledgements

This work is part of the thesis by Coralie Coudray-Meunier, a PhD student who received financial support from ANSES.

Authors’ Affiliations

(1)
Universit Paris-Est, ANSES, Food Safety Laboratory, Enteric viruses unit
(2)
AP-HP, Hpital Paul Brousse
(3)
Univ Paris-Sud
(4)
INSERM

References

  1. 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.View ArticlePubMedGoogle Scholar
  2. 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.View ArticlePubMedGoogle Scholar
  3. 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.View ArticleGoogle Scholar
  4. Koopmans M, Duizer E: Foodborne viruses: an emerging problem. Int J Food Microbiol. 2004, 90: 23-41. 10.1016/S0168-1605(03)00169-7.View ArticlePubMedGoogle Scholar
  5. 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.View ArticleGoogle Scholar
  6. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  7. Hepatitis A vaccines. Wkly Epidemiol Rec. 2000, 75: 38-44.Google Scholar
  8. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  9. 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]
  10. 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.View ArticleGoogle Scholar
  11. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  12. 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.View ArticlePubMedGoogle Scholar
  13. 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, KoreaGoogle Scholar
  14. 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-PubMedGoogle Scholar
  15. 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-PubMedGoogle Scholar
  16. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  17. Koopmans M, von Bonsdorff CH, Vinjé J, de Medici D, Monroe S: Foodborne viruses. FEMS Microbiol Rev. 2002, 26: 187-205.View ArticlePubMedGoogle Scholar
  18. 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.View ArticlePubMedGoogle Scholar
  19. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  20. 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, SwitzerlandGoogle Scholar
  21. 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, SwitzerlandGoogle Scholar
  22. 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.View ArticleGoogle Scholar
  23. 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.View ArticleGoogle Scholar
  24. 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.View ArticlePubMedGoogle Scholar
  25. 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.View ArticlePubMedGoogle Scholar
  26. 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-PubMedGoogle Scholar
  27. Stapleton JT, Lemon SM: Neutralization escape mutants define a dominant immunogenic neutralization site on hepatitis A virus. J Virol. 1987, 61: 491-498.PubMed CentralPubMedGoogle Scholar
  28. Costa-Mattioli M: Genetic variability of hepatitis A virus. J Gen Virol. 2003, 84: 3191-3201. 10.1099/vir.0.19532-0.View ArticlePubMedGoogle Scholar
  29. 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.View ArticlePubMedGoogle Scholar
  30. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  31. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  32. 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.View ArticlePubMedGoogle Scholar
  33. 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.View ArticlePubMedGoogle Scholar
  34. 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.View ArticleGoogle Scholar
  35. 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.View ArticlePubMedGoogle Scholar
  36. 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.View ArticlePubMedGoogle Scholar
  37. 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.View ArticlePubMedGoogle Scholar
  38. 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.View ArticlePubMedGoogle Scholar
  39. 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.View ArticlePubMedGoogle Scholar
  40. 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.PubMed CentralPubMedGoogle Scholar
  41. 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.View ArticleGoogle Scholar
  42. 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.View ArticlePubMedGoogle Scholar
  43. 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.View ArticlePubMedGoogle Scholar
  44. Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMed CentralView ArticlePubMedGoogle Scholar
  45. 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.PubMed CentralView ArticlePubMedGoogle Scholar
  46. 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) 6688PubMed CentralView ArticlePubMedGoogle Scholar

Copyright

© Coudray-Meunier et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.