Abstract
Oysters, being filter feeders, can accumulate some human pathogens such as norovirus, a highly infectious calicivirus, most common cause of acute gastroenteritis worldwide. Accumulated virus decays over a period of days to weeks, possibly rendering contaminated oysters safe again. Sensitive molecular methods have been set up for shellfish analysis but without answering the question of infectious virus detection. Using the Tulane virus (TV), a norovirus surrogate that recognizes the same ligand as human norovirus in oyster tissues, the genome and infectious virus decay rates were estimated using inverse linear regression in a Bayesian framework for genome copies. Infectivity decreased faster than genome copies but infectious viruses were detected for several days. Quantifying the decrease in viral infectivity and genome detection in oysters over such a long period may help local authorities to manage production areas implicated in shellfish-borne outbreaks, and thus protect consumers.
Introduction
Shellfish, being filter feeders, are known to accumulate different types of pathogens of human fecal pollution. As a consequence, they have been implicated in many gastroenteritis outbreaks, and norovirus (HuNoV) is the agent most frequently detected following oyster consumption (; ). Several factors contribute to these outbreaks. First, among the human enteric viruses, HuNoVs are recognized as the leading cause of epidemics and sporadic cases of gastroenteritis in all age groups (). Additionally, HuNoVs are excreted in large quantities by ill people as well as by asymptomatic individuals (). Given the foregoing, huge numbers of HuNoV particles are discharged in sewage and due to their resistance to inactivation, they are frequently detected in wastewater treatment plant effluent and in surface waters (; ). Second, depuration of shellfish, which was developed to eliminate bacteria, does not efficiently eliminate viruses, which can persist for several weeks or months in bivalve tissues (). The identification of specific ligands for human HuNoV in various oyster species may contribute to this long persistence (, ; ).
Several studies have clearly demonstrated the slow decrease of HuNoV concentrations in oyster samples using genome detection (). But how long infectious virus persists remain unclear. Until recently HuNoV was not cultivable, however, a newly described culture system based on enteroids offers now the capacity to detect infectious HuNoV (). This recent development would help to solve the question of whether a surrogate virus can be identified that correctly mimics the behavior of HuNoVs, and to validate alternative method aiming to detect capsid integrity (). Selection of an appropriate surrogate requires consideration of several issues, such as environmental persistence, physical and genetic similarities, ease of propagation, and tissue binding properties (). We previously demonstrated the potential of Tulane virus (TV), a member of the recovirus genus (ReCV) of the Caliciviridae family, to mimic HuNoV behavior in oyster tissues as it showed the same tissue distribution and a comparable decrease of genomic copies over time ().
The objectives of the present study were:
• to quantify the decrease of TV genome copies as inferred from real-time RT-PCR (rRT-PCR) and the decrease of TV infectious particles over time in oysters reared in a scientific oyster farm; and
• to evaluate the measurement error estimated from standard curve data, by using inverse linear regression in a Bayesian framework to more accurately assess the decay rate for genome and infectivity detection.
Materials and Methods
Viruses and Cell Lines
Mengovirus (MgV) strain pMC0 (provided by A. Bosch, University of Barcelona, Spain) was propagated in HeLa cells, and TV strain M033 (provided by T. Farkas, Louisiana State University at Baton Rouge, LA, United States) in LLC-MK2 cells (ATCC CCL-7; ATCC, Manassas, VA, United States), as previously described (; ). To produce high titer of TV, several large flasks were infected at low multiplicity of infection (MOI). When cytopathic effect (CPE) was complete, cultures were frozen and thawed (-20°C) three times, centrifuged, and the supernatant filtered to 0.22 μm (Minisart®, Sartorius, Germany). Viral suspensions were kept frozen at -80°C in aliquots, and titrated by TCID50 or rRT-PCR as described below.
TV Bioaccumulation in Oysters
Four bioaccumulation experiments (named A, B, C, and D) were performed in March, April, May, and early June 2016, each used 250 oysters (Crassostrea gigas) directly purchased from the same producer (few days before each experiment). Aquariums were filled with 100 L of natural seawater seeded with 11.3 ± 0.12 log10 RNA copies of TV, which corresponds to 7.14 ± 0.1 log10 TCID50. After 24 h of bioaccumulation, virus concentrations were checked in a random sample of 20 oysters.
Persistence Experiments
Following bioaccumulation, the oysters were rinsed and transported to an experimental farm with direct access to natural seawater and facilities to maintain oysters for long periods (). Oysters were placed in clean seawater in large tanks (500 L) and supplied with constantly circulating aerated and filtered natural seawater (200 L/hour/tank). Parameters such as temperature and salinity were routinely measured. Prior to release into the environment, seawater was treated according to the safety rules of the experimental farm.
Sample Processing and Virus Recovery
Oysters samples, comprised of 20 individuals each, were collected at time 0 h, 2 and 3 days, and then weekly up to 4 weeks. For the two last experiments (C and D) an additional sampling was performed at day 25. Samples were received at the laboratory in less than 24 h, except for experiment B for which the four first samples were frozen. On arrival, all oysters were weighed, immediately shucked, flesh weight recorded to calculate the allometric coefficient (). Then the digestive tissues (DTs) were dissected and distributed into 1.5 g aliquots. One DT aliquot was immediately homogenized in 2 mL of phosphate-buffered saline (PBS, pH 7.4), extracted by vortexing with an equal volume of chloroform-butanol for 30 s, and treated with Cat-Floc T (Calgon, Ellwood City, PA, United States) for 5 min on the bench before centrifugation for 15 min at 13,500 × g. The resulting suspension was precipitated with polyethylene glycol 8000/NaCl (PEG 8000) (Sigma, Saint-Quentin, France) for 1 h at 4°C and centrifuged for 20 min at 11,000 × g at 4°C ().
Purification for Infectious TV Detection
Then, PEG pellets were suspended with 10 mL of medium M199 (Gibco®, Life Technologies) supplemented with antibiotics-antimitotic (Gibco®, Life Technologies). An aliquot of 100 μL was kept for nucleic acid (NA) extraction. The mixture was centrifuged at 2862 × g for 10 min at 4°C, filtered through a 0.22 μm filter, and concentrated using an Amicon R Ultra-15 (100 K) Centrifugal Filter Device by centrifugation at 5000 × g for 20 min at 4°C, to a volume of 1.5 ± 0.2 ml. This concentrate was loaded into PD-10 Desalting columns containing Sephadex G-25 Medium, and equilibrated with M199 medium up to a volume of 2.5 mL, and centrifuged at 1000 × g for 2 min. The supernatant was recovered and used for infectivity assays, except for a 100 μL aliquot, which was kept for NA extraction to verify that purifications steps did not induce loss of viral genomic copies.
Infectivity Assay and Titration Using TCID50 Assay
Purified supernatants (1 mL) were inoculated on confluent LLC-MK2 cells in T 75 flasks with 9 mL of M199 medium at 37°C. After 1 h of contact, an aliquot of 100 μL was kept and stored frozen. After incubation for 7 days in an incubator at 37°C with gentle agitation, 100 μl was collected and compared to the 1-h control by rRT-PCR. Samples showing an increase of more than three Logs (9 CT-values) were considered to be positive.
For titration of samples, 1 mL of purified supernatant was serially diluted (10-1 to 10-8) in medium and inoculated in 10-well for each dilution onto cells in 96-well culture plates (Thermo Fisher Scientific, France) under a volume of 200 μL, and incubated for 7 days at 37°C. For each experiment two 96-well culture plates were inoculated. CPE counts were converted into TCID50 units using the Reed and Muench calculation method () and were expressed per g of DT. A negative control was included on each plate.
NA Extraction and Genome rRT-PCR Detection
For each sample, two additional extractions were performed on frozen DT aliquots with MgV added as an extraction control before homogenization. All the first steps were as described as above except that the PEG pellet was recovered in 1 mL of sterile water and NA were extracted immediately.
Nucleic acid extractions were carried out using the NucliSENS kit (bioMérieux, Lyon, France) and an automatic easyMAG extractor (bioMérieux, Lyon, France) according to the manufacturer’s instructions, with minor modifications, as described previously (). NAs were recovered in 100 μL of elution buffer (bioMérieux) and triplicate 5 μL aliquots were analyzed immediately using the UltraSense One-Step quantitative RT-PCR system (Life Technologies, France). For TV, primers TVIFf (5′-CTGGGATACCCACAACATC-3′), TVIFr (5′’-GCCAGTTAACAGCTTCAGC-3′) and probe TVIFp (FAM-TGTGTGTGCCACTGGATAGCTAGCACC-BHQ) were used to amplify a region between nucleotides 3775 and 3884 (GenBank accession number EU391643.1) (). MgV spiked into the two additional DT extracts was detected as previously described (). A standard curve was amplified for each run as previously described ().
Statistical Analysis
Genome and infectivity data were analyzed jointly in a two-level model. In both the infectivity and the genome exponential decay was assumed, using linear regression on the log-transformed numbers of virus particles or genome copies. In the genome model, measurement error was estimated from standard curve data by using inverse linear regression in a Bayesian framework (see Supplementary Material). Error estimates in the infectivity data were used to weight individual infectivity measurements. Numbers of infectious viruses were expressed as a fraction of “total” viruses, as estimated from numbers of genome copies. Hence, the decay of total virus and of the fraction of infectious virus could both be estimated. The two-level model then allowed for estimation of the overall (average) decay rates in total and infectious virus.
The variations in temperature, salinity and allometric coefficient of oysters were compared between the four experiments using ANOVA Fisher and Tukey HSD tests, and were considered significant with a p-value of 0.05.
Results
Controls and Data Robustness
The oysters used in this study were grown in the natural environment and thus, there may be physiological differences between oyster batches used for each experiment due to small variations in the environmental conditions. To investigate the similarity in physiological state of oyster batches, we measured various environmental parameters and calculated the allometric coefficient (Table 1). The mean body weight per oyster and the allometric coefficients were similar in three experiments out of four. For the last experiment performed in June, these values appeared somewhat higher. However, the difference was not statistically significant (p = 0.2712 and p = 0.1091, respectively). The body weight remained stable over the 4 weeks for all experiments (data not shown). The main difference observed between experiments was water temperature (p = 0.0001), and a lower salinity was also noted for experiment A (p = 0.0001) (Table 1).
Table 1
| Experiments | Sea water | Oyster | Quality controls | |||
|---|---|---|---|---|---|---|
| Temp | Salinity | Body weight | DT weight | allom. coeff | ext. eff | |
| A | 9.1 ± 1.6 | 30.0 ± 1.2# | 5.86 ± 0.85 | 0.47 ± 0.09 | 12.99 ± 3.30 | 36 ± 6 |
| B | 11.4 ± 2.1 | 31.5 ± 0.8 | 6.05 ± 0.73 | 0.47 ± 0.03 | 12.84 ± 2.00 | 42 ± 14 |
| C | 13.7 ± 2.2 | 32.3 ± 0.5 | 5.93 ± 0.50 | 0.47 ± 0.09 | 12.55 ± 2.37 | 36 ± 10 |
| D | 17.2 ± 1.8# | 31.7 ± 1.4 | 8.29 ± 0.81 | 0.55 ± 0.11 | 15.18 ± 2.05 | 33 ± 10 |
Physical characteristics recorded for the four experiments.
The allometric coefficient (allom. coeff) was calculated by dividing the total body weight by the weight of the digestive tissues (DT), and it is expressed here as the arithmetic mean values for all samples analyzed per experiment. The extraction efficiency (ext. eff) represents the mengovirus recovery for each extraction and is expressed as the average value in % obtained for all extracts per experiment. #: statistically different from the other experiments, (p = 0.0001).
Mengovirus, added to the two extractions made for genome detection, showed comparable extraction efficiency coefficients, varying from 33 to 42% for all four experiments (Table 1).
An additional control was added on the sample prepared for infectivity detection. Indeed to verify that no TV was lost during the purification steps, a portion of the eluate after the Sephadex column was extracted and analyzed by rRT-PCR. Thus a total of three NA extracts were analyzed for each sample, and each undiluted NA extract was analyzed by rRT-PCR in triplicate (3 × 5 μL).
rRT-PCR Data
The initial virus concentration in oysters was comparable across batches (2.68–4.32 × 106 RNA/g DT); this consistency is important for monitoring and comparing virus decay (Table 2).
Table 2
| Time (Days) | Experiment A | Experiment B | Experiment C | Experiment D | ||||
|---|---|---|---|---|---|---|---|---|
| RNA | TCID50 | RNA | TCID50 | RNA | TCID50 | RNA | TCID50 | |
| 0 | 4.32 × 106 | 7.43 × 103 | 2.68 × 106 | 4.18 × 102 | 2.85 × 106 | 1.87 × 104 | 1.99 × 106 | 1.32 × 104 |
| 2 | 9.74 × 105 | 1.48 × 103 | 8.41 × 105 | 18.7 | 2.54 × 106 | 3.72 × 103 | No sample | |
| 3 | 1.48 × 106 | 1.40 × 103 | 1.66 × 106 | 9 | 1.59 × 106 | 1.87 × 103 | 2.17 × 106 | 3.72 × 103 |
| 7 | 4.84 × 105 | 4.69 × 102 | 7.45 × 105 | 7 | 1.04 × 106 | 3.72 × 102 | 6.56 × 105 | 5.90 × 102 |
| 14 | 1.91 × 105 | 7.43 | 6.41 × 104 | Positive | 3.91 × 104 | 4.7 | 9.95 × 104 | Positive |
| 21 | 1.84 × 104 | Positive | 8.00 × 103 | Positive | 4.55 × 103 | Positive | 1.18 × 104 | |
| 25 | No sample | No sample | 6.38 × 103 | 4.46 × 103 | ||||
| 28 | 2.53 × 103 | 6.63 × 102 | 3.30 × 103 | 1.11 × 103 | ||||
Virus concentrations obtained for the four separate experiments as measured by rRT-PCR (RNA genome copies) and infectivity assay (TCID50).
RNA: Geometric mean titer calculated based on the concentrations obtained from the three separate extractions and triplicates of rRT-PCR and expressed per g of DT. TCID50: mean values of two titration plates performed per sample and expressed per g of DT (values in italics for experiment B were derived from samples that were frozen prior to the infectivity test). ‘Positive’ results indicate that some infectious viruses were present but at concentrations too low for titration.
Quantification using rRT-PCR depends on standard curves. To improve error estimation, the Bayesian inverse regression model was used to jointly analyze standard curve data and sample data. This procedure combines the effects of uncertainty in CT-values and (inverse) regression errors on standard curves, on estimates of numbers of genome copies per g of DT. All standard curves were included in the statistical analysis to account for effects of measurement errors associated with CT-values. For example, a standard curve which includes the 95% prediction interval for the CT-values, as estimated in the Bayesian model, shows that the regression error is small over the range of genome copy numbers covered by the standards, but tends to increase when extrapolating outside the range (Figure 1). Using these errors for each standard curve, the number of genome copies was calculated for each extraction, based on the three replicate CT-values (after quality control, verifying that the difference among the three values was lower than 2 CT). Resulting estimates of genome copy numbers and 95% posterior credibility intervals can be given for all samples analyzed (Figure 2). A total of 12 measurements (four experiments quantified three times) were obtained for each sampling time, except for day 25, which was sampled only for the last two experiments. At low virus concentrations, observations are censored (CT > 40) and the error estimates increase since confidence in CT-values between 40 and 45 is limited. Also shown are the overall mean decay curves for all four experiments as predicted by the Bayesian hierarchical model (mean decay and 95% credibility) (Figure 2).
FIGURE 1
FIGURE 2
Estimated decay rates were calculated for each experiment separately (Figure 3), and despite some differences observed in the water temperature and salinity, values obtained are comparable (Table 3). The mean decay of 0.120 log10 units per day (95% CI 0.10–0.15 per day) corresponds to a decay rate for detectable viruses of 3.6 log10 units in 1 month (30 days) and a half-life of 2.5 days (95% CI 2.0–3.1 days).
FIGURE 3
Table 3
| Experiment | Concentration1 | Decay rate2 | ||||
|---|---|---|---|---|---|---|
| Median | P2.5 | P97.5 | Median | P2.5 | P97.5 | |
| RT-PCR (Genome Copies) | ||||||
| A | 1.54 × 106 | 0.84 × 106 | 2.74 × 106 | 0.117 | 0.100 | 0.136 |
| B | 1.55 × 106 | 0.80 × 106 | 2.83 × 106 | 0.125 | 0.105 | 0.145 |
| C | 2.38 × 106 | 1.35 × 106 | 4.28 × 106 | 0.124 | 0.108 | 0.140 |
| D | 2.11 × 106 | 1.15 × 106 | 4.00 × 106 | 0.116 | 0.100 | 0.133 |
| Overall | 1.86 × 106 | 0.85 × 106 | 3.71 × 106 | 0.121 | 0.098 | 0.146 |
| Infectious virus (TCID50) | ||||||
| A | 6.57 × 103 | 4.77 × 103 | 9.03 × 103 | 0.200 | 0.179 | 0.221 |
| B3 | NA | NA | NA | NA | NA | NA |
| C | 1.32 × 104 | 9.50 × 103 | 1.86 × 104 | 0.240 | 0.221 | 0.260 |
| D | 1.34 × 104 | 8.34 × 103 | 2.18 × 104 | 0.196 | 0.153 | 0.237 |
| Overall | 9.86 × 103 | 2.94 × 103 | 3.40 × 104 | 0.214 | 0.147 | 0.325 |
Initial numbers and decay rates of Tulane virus (TV) as measured by rRT-PCR (genome copies) and infectivity assay (TCID50), as estimated by the two-level Bayesian regression model.
The 97.5th and 2.5th percentile values are also shown. 1Concentrations are expressed as RNAc/g of DT or TCID50/g of DT with median estimates and 95% range calculated based on the four extraction results obtained. 2Decay rates, expressed as decrease in log10 numbers of RNAc or TCID50 per day, in each of the four experiments, and estimated mean decay rate over all experiments combined (median estimates and 95% range). 3TCID50 data from experiment B discarded as freezing these samples prior to analysis clearly resulted in a reduction in infectious virus. NA, nucleic acid.
Detection of Infectious TV
Genomic concentrations of inoculum prepared for infectivity showed a decrease of about 50%. Within the error range of the quantification by rRT-PCR, and thus all these inoculum were used for infectivity tests. Infectious TV was detected for up to 21 days in three experiments, and for 14 days in experiment D (Table 2). Experiments A and C showed similar decreases in infectious virus, with quantification reaching the limit of detection after 2 weeks, and after 1 week for experiment D. For this last experiment, infectious TV was detected only for 14 days presumably because of a slight increase in water temperature. The fastest decrease in infectious TV for experiment B was probably caused by freezing the samples at -20°C before analysis, and therefore, the results were not included in the data analysis.
When added to the Bayesian regression model, the infectivity data allowed the fraction of infectious virus to be estimated (compared to the total number of genome copies). Assuming exponential decay of the infectious fraction, the estimated initial numbers and decay rates are shown in the lower part of Table 3 and Figure 3. The decay rate of infectious virus is estimated at 0.214 log10 units per day (95% CI 0.147–0.325 per day), which corresponds to a 6.42 log10 unit decrease in 1 month (30 days). The corresponding half-life is 1.4 days (95% CI 0.9–2.0 days).
Discussion
To infect their human host, HuNoVs required to bind to carbohydrates of the histo-blood group antigen (HBGA) family (; ). Binding of HuNoV to similar ligands was demonstrated in oyster digestive tissues (), impacting bioaccumulation efficiency, with seasonal differences shown for some strains (, ). In addition to HuNoVs, within the Caliciviridae family, members of the lagovirus and recovirus genera also recognize HBGAs of their respective hosts (, ; ). Yet, among these viruses only some ReCVs such as TV are cultivable to high titers in susceptible cells and present the same stability as HuNoVs in seawater, as demonstrated by genome detection (). Moreover, TV and Norwalk virus have similar tissue distributions and half-lives in oysters, confirming that TV is an appropriate surrogate for studying HuNoV behavior in shellfish (). Infectious TV particles detected after the heat treatment of eastern oysters (Crassostrea virginica) confirmed its higher stability compared to murine norovirus (MNV-1) (). Additionally, TV treated at 70°C lost 100% of infectivity (84% at 56°C for 2 min) (), a rate comparable to the inactivation of GII.3 and GII.4 as demonstrated on human enteroids (), and this confirmed its interest as a HuNoV surrogate ().
Infectivity studies that can be undertaken with cell culture-adapted ReCVs are therefore of utmost importance for shellfish. Detecting infectious viral particles in oysters requires methods of virus recovery that do not damage the virus capsid, and eliminate the cytotoxic impact of the oyster matrix on cell cultures. Long before the development of molecular techniques, sensitive methods were developed to detect infectious enterovirus from various environmental samples, including shellfish (). PEG was proved an effective concentration method to recover a variety of infectious human enteric viruses such as hepatitis A virus, rotavirus and poliovirus from oyster tissues (). Later on, PEG was also used for molecular detection (), and to detect HuNoV in shellfish-related outbreaks (, ). We reasoned that the use of identical initial recovery steps for both genome and infectious virus detection would enable the results to be more meaningfully compared, and thus, we decided to focus on the PEG method. After several failed attempts, some purification steps were added to reduce the toxicity to the cells, allowing the successful detection of infectious virus while minimizing the lost of particles. This PEG method was then applied to evaluate the decrease of genome copies and infectious particles in contaminated oysters kept in a scientific farm to mimic producer’s farm environment. Another approach using direct elution with no centration step was reported efficient but need higher concentrations ().
Quantification of genome copies in oyster DTs was achieved using rRT-PCR. Despite the application of controls as recommended in the ISO/TS 16216-1 method (), questions have been raised regarding the accuracy and variability in RNA quantification using rRT-PCR between laboratories, and error propagation through the use of reference materials. Analysis of quantitative PCR data is not a trivial problem: aside from uncertainties in measured CT-values, as revealed by differences between replicate measurements, the translation of CT-values into genome copies causes additional uncertainty. Indeed, the standard curve, as defined by the same (uncertain) CT data, must be inverted to translate CT-values into genome copy number estimates. Quality criteria to validate these curves, and triplicate sample extractions provide additional assurances regarding quantification (; ). In order to use all information provided by this hierarchical setup, and estimate the resulting uncertainty, appropriate statistical models considering all sources of uncertainty are needed. The Bayesian approach, nesting the sample data into the standard curve regression model, with CT-values known and genome copy numbers missing, considers the error associated to the CT-value of the sample. Running the Bayesian model then produces estimates for missing data, or rather a sample from their joint posterior distribution. Even high CT-values above the cutoff for quantification, e.g., CT-values greater than 40, can be included as censored observations, with the mere effect of inflated error estimates (Figure 2).
Using this approach, we observed a mean decay rate value of 0.121 log10 RNAc per day, corresponding to a half-life of 2.0–3.1 days, shorter but in the same range, than our previous estimate of 4.56 days (). The accuracy of the quantification method applied here and other factors such as the use of different oyster batches and variable water content (nutrients, salinity, water temperature…) may explain this difference. Regarding the temperature parameter, the same decay rate was observed for the experiment performed with water temperature around 17°C compared to the others, consistent with earlier studies showing no significant decay of HuNoV at this temperature (; ). However, it would interesting to further evaluate the impact of this parameter by repeating experiments in various condition. The model used, estimating variability associated to the measurement for low concentrations, allows us to conclude that the concentration variability observed in oyster tissues after 1 month is likely linked to inter-animal variability rather than to quantification errors, as PCR variation at this concentration was lower than the overall variation observed. Oysters through their physiological activity may eliminate the viral particles at variable individual rates. Future developments may include a more comprehensive regression model including covariates like oyster batch, water temperature or salinity.
The major result of this study is the demonstration that a ReCV, a virus closely related to the Norovirus genus, remains infectious for 3 weeks in oysters immersed in natural running seawater. The mean decay rate for infectious particles was 0.214 log10 TCID50 per day compared to 0.121 log10 RNAc per day, thus, the long-lasting question of whether infectious particles decrease at a faster rate than genome copies is partly answered. Reports on the detection of infectious viruses from shellfish are rare, making these data original. Although the infectious virus titer decreases faster than the concentration in viral genome copies, it is only of a factor two. Indeed, infectious particles were still present after 3 weeks confirming observations made following illness outbreaks in France, from which we concluded that the contamination events occurred 2 or 3 weeks before marketing and consumption of contaminated oysters (; ). Hepatitis A cases in the Netherlands were also linked to contamination events occurring several weeks before consumption (). Our study also underscores the need for further method refinements to improve the sensitivity and to increase the sample size. Such an approach will be beneficial both for genome detection as it will move the CT-values into an area where the PCR gives more reliable results, and for the detection of infectious viruses, as after 3 weeks some infectious particles were detected, albeit below the limit of quantification. Being able to detect infectious Tulane virus is a step forward toward a better understanding of viral persistence in oysters. This model could be used to validate potential treatments such as longer depuration time, but also other methods such as in situ capture using pig mucin or HBGA (; ), keeping in mind the ultimate goal to compare with infectious HuNoV using the model based on enteroids ().
Our results, based on genome detection and infectious particle detection of the Tulane virus, member of the Caliciviridae family, show the long persistence of infectious virus, in shellfish and highlight the need to date more precisely contamination events to estimate the infectious risk associated with detection of the viral genome. This strengthens the recommendation for closure of production areas for several weeks after a contamination event, such as the 4-week closure applied in France. Such data will be useful for local authorities to manage producing areas implicated in outbreaks, and thus to protect shellfish consumers.
Statements
Author contributions
DP and JS performed the experiments and analysis. PT designed the model and analyzed the data. VB was in charge of the shellfish care and related analysis. FLG designed the study and wrote the manuscript.
Funding
This study was supported by the Direction Générale de l’Alimentation (French Government DGAl, convention LNR), by EU-H2020 Grant No. 643476 COMPARE, and by IFREMER (Direction Scientifique) through financial support for DP.
Acknowledgments
We acknowledge the help of colleagues from the Experimental Platform (Bouin) for taking care of oysters. We also acknowledge the scientific and technical advice of Tibor Farkas (Louisiana State University, Baton Rouge, LA, United States) for Tulane virus detection. We are grateful to Jacques Le Pendu and Catherine McLeod for the critical reviews of the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.00716/full#supplementary-material
References
1
AraudE.DiCaprioE.MaY.LouF.GaoY.KingsleyD.et al (2016). Thermal inactivation of enteric viruses and bioaccumulation of enteric foodborne viruses in live oysters (Crassostrea virginica).Appl. Environ. Microbiol.822986–2099. 10.1128/AEM.03573-15
2
AtmarR. L.NeillF. H.RomaldeJ. L.Le GuyaderF.WoodleyC. M.MetcalfT. G.et al (1995). Detection of Norwalk virus and Hepatitis A virus in shellfish tissues with the PCR.Appl. Environ. Microbiol.613014–3018.
3
AtmarR. L.OpekumA. R.GilgerM. A.EstesM. K.CrawfordS. E.NeillF. H.et al (2014). Determination of the 50% human infectious dose for Norwalk virus.J. Infect. Dis.2091016–1022. 10.1093/infdis/jit620
4
BoxmanI.VerhoefL.VennemaH.NguiS.FriesemaI. H.WhitesideC.et al (2016). International linkage of two food-borne hepatitis A clusters through traceback of mussels, the Netherlands, 2012.Eurosurveillance2130113. 10.2807/1560-7917.ES.2016.21.3.30113
5
CromeansT.ParkG. W.CostantiniV.LeeD.WangQ.FarkasT.et al (2014). Comprehensive comparison of cultivable norovirus surrogates in response to different inactivation and disinfection treatments.Appl. Environ. Microbiol805743–5751. 10.1128/AEM.01532-14
6
De GraafM.van BeekJ.KoopmansM. P. G. (2016). Human norovirus transmission and evolution in a changing world.Nat. Rev.14421–433. 10.1038/nrmicro.2016.48
7
DoreB.KeaveneyS.FlanneryJ.Rajko-NenowP. (2010). Managment of health risks associated with oysters harvested from a norovirus contaminated area, Ireland, February-March 2010.Eurosurveillance1510–13.
8
DrouazN.SchaefferJ.FarkasT.Le PenduJ.Le GuyaderF. S. (2015). Tulane virus as a potential surrogate to mimic norovirus behavior in oysters.Appl. Environ. Microbiol.815249–5256. 10.1128/AEM.01067-15
9
EFSA (2012). Scientific Opinion on Norovirus (NoV) in oysters: methods, limits and control options.EFSA J.10:2500. 10.2903/j.efsa.2012.2500
10
EttayebiK.CrawfordS. E.MurakamiK.BroughmanJ. R.KarandikarU.TengeV. R.et al (2016). Replication of human noroviruses in stem cell-derived human enteroids.Science3531387–1393. 10.1126/science.aaf5211
11
FarkasT.LunC. W. P.LeyB. (2014). Relationship between gentypes and serotypes of genogroup 1 recoviruses: a model for human norovirus antigenic diversity.J. Gen. Virol.951469–1478. 10.1099/vir.0.064675-0
12
FarkasT.SestakK.WeiC.JiangX. (2008). Characterization of a rhesus monkey calicivirus representing a new genus of caliciviridae.J. Virol.825408–5416. 10.1128/JVI.00070-08
13
International Organization for Standardization (2013). ISO/TS 15216-1 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.Geneva: ISO.
14
KazamaS.MasagoY.TohmaK.SoumaN.ImagawaT.SuzukiA.et al (2016). Temporal dynamics of norovirus determined through monitoring of municipal wastewater by pyrosequencing and virological surveillance of gastroenteritis cases.Water Res.92244–253. 10.1016/j.watres.2015.10.024
15
KnielK. (2014). The makings of a good human norovirus surrogate.Curr. Opin. Virol.485–90. 10.1016/j.coviro.2014.01.002
16
Le GuyaderF.ParnaudeauS.SchaefferJ.BoschA.LoisyF.PommepuyM.et al (2009). Detection and quantification of noroviruses in shellfish.Appl. Environ. Microbiol.74618–624. 10.1128/AEM.01507-08
17
Le GuyaderF. S.AtmarR. L.Le PenduJ. (2012). Transmission of viruses through shellfish: when specific ligands come into play.Curr. Opin. Virol.2103–110. 10.1016/j.coviro.2011.10.029
18
Le GuyaderF. S.KrollJ.Ambert-BalayK.Ruvoen-ClouetN.DesaubliauxB.ParnaudeauS.et al (2010). Comprehensive analysis of a norovirus-associated gastroenteritis outbreak, from the environment to the consumer.J. Clin. Microbiol.48915–920. 10.1128/JCM.01664-09
19
Le GuyaderF. S.Le SauxJ.-C.Ambert-BalayK.KrolJ.SeraisO.ParnaudeauS.et al (2008). Aichi virus, norovirus, astrovirus, enterovirus and rotavirus involved in clinical cases from a French oyster-related gastroenteritis outbreak.J. Clin. Microbiol.464011–4017. 10.1128/JCM.01044-08
20
Le GuyaderF. S.LoisyF.AtmarR. L.HutsonA. M.EstesM. K.Ruvoen-ClouetN.et al (2006). Norwalk virus specific binding to oyster digestive tissues.Emerg. Infect. Dis.12931–936. 10.3201/eid1206.051519
21
Le MennecC.ParnaudeauS.RumebeM.Le SauxJ.-C.PiquetJ.-C.Le GuyaderF. S. (2017). Follow-up of norovirus contamination in an oyster production area linked to repeated outbreaks.Food Environ. Virol.954–61. 10.007/s12560-016-9260-6
22
Le PenduJ.NystromK.Ruvoen-ClouetN. (2014). Host pathogen co-evolution and glycan interactions.Curr. Opin. Virol.788–94. 10.1016/j.coviro.2014.06.001
23
LewisG. D.MetcalfT. G. (1988). Polyethylene glycol precipitation for recovery of pathogenic virues, including Hepatitis A virus and human rotavirus, from oyster, water and sediment samples.Appl. Environ. Microbiol.54k1983–1988.
24
MaaloufH.SchaefferJ.ParnaudeauS.Le PenduJ.AtmarR. L.CrawfordS. E.et al (2011). Strain-dependent norovirus bioaccumulation in oysters.Appl. Environ. Microbiol.773189–3196. 10.1128/AEM.03010-10
25
MaaloufH.ZakhourM.Le PenduJ.Le SauxJ.-C.AtmarR. L.Le GuyaderF. S. (2010). Norovirus genogroup I and II ligands in oysters: tissue distribution and seasonal variations.Appl. Environ. Microbiol.765621–5630. 10.1128/AEM.00148-10
26
MartinL. R.DukeG. M.OsorioJ. E.HallD. J.PalmenbergA. C. (1996). Mutational analysis of the Mengovirus Poly(C) tract and surrounding heteropolymeric sequences.J. Virol.702027–2031.
27
McLeodC.HayB.GrantC.GreeningG.DayD. (2009). Inactivation and elimination of human enteric viruses by Pacific oysters.J. Appl. Microbiol.1071809–1818. 10.1111/j.1365-2672.2009.04373.x
28
McLeodC.PoloD.Le SauxJ.-C.Le GuyaderF. S. (2017). Depuration and relaying: a review on potential removal of norovirus from oysters.Comp. Rev. Food Sci. Food Saf.16692–706. 10.1111/1541-4337.12271
29
MetcalfT. G.MoultonE.EckersonD. (1980). Improved method and test strategy for recovery of enteric viruses from shellfish.Appl. Environ. Microbiol.39141–152.
30
PintoR. M.CostafredaM. I.BoschA. (2009). Risk assessment in shellfish-borne outbreaks of hepatitis A.Appl. Environ. Microbiol.757350–7355. 10.1128/AEM.01177-09
31
PoloD. S. J.FournetN.Le SauxJ.-C.ParnaudeauS.McLeodC.Le GuyaderF. S. (2016). Digital PCR for quantifying norovirus in oysters implicated in outbreaks, France.Emerg. Infect. Dis.222189–2191. 10.320/eid2212.160841
32
PouillotR.Van DorenJ. M.WoodsJ.PlanteD.SmithM.GoblickG.et al (2015). Meta-analysis of the reduction of norovirus and male-specific coliphage concentrations in wastewater treatment plants.Appl. Environ. Microbiol.814669–4681. 10.1128/AEM.00509-15
33
RandazzoW.KhezriM.OllivierJ.Le GuyaderF. S.Rodriguez-DiazJ.AznarR.et al (2018). Optimization of PMAxx pretreatment to distinguish between human norovirus with intact and altered capsids in shellfish and sewage samples.Int. J. Food Microbiol.2661–7. 10.1016/j.ijfoodmicro.2017.11.011
34
ReedL.MuenchH. (1938). A simple method of estimating fifty per cent end points.Am. J. Hyg.27493–497.
35
SchwabK. J.NeillF. H.EstesM. K.MetcalfT. G.AtmarR. L. (1998). Distribution of Norwalk virus within shellfish following bioaccumulation and subsequent depuration by detection using RT-PCR.J. Food Prot.611674–1680. 10.4315/0362-028X-61.12.1674
36
TeunisP. F. M.SukhrieF. H. A.VennemaH.BogermanJ.BeersmaM. F. C.KoopmansM. P. G. (2015). Shedding of norovirus in symptomatic and asymptomatic infections.Epidemiol. Infect.1431710–1717. 10.1017/S095026881400274X
37
TianP.EngelbrektsonA. L.JiangX.ZhongW.MandrellR. E. (2007). Norovirus recognizes histo-blood group antigens on gastrointestinal cells of clams, mussels, and oysters: a possible mechanism of bioaccumulation.J. Food Prot.702140–2147. 10.4315/0362-028X-70.9.2140
38
TianP.YangD.ShanL.QianqianL.LiuD.WangD. (2017). Estimation of human norovirus infectivity from environmental water samples by in situ capture RT-qPCR method.Food Environ. Virol.1029–38. 10.1007/s12560-017-9317-1
39
WangD.XuS.YangD.YoungG. M.TianP. (2014). New in situ capture quantitative (real-time) reverse transcription-PCR method as an alternative approach for determining inactivation of Tulane virus.Appl. Environ. Microbiol.802120–2124. 10.1128/AEM.04036-13
40
YuY.CaiH.HuL.LeiR.PanY.YanS.et al (2015). Molecular epidemiology of oyster-related human noroviruses and their global genetic diversity and temporal-geographical distribution from 1983 to 2014.Appl. Environ. Microbiol.817615–7624. 10.1128/AEM.01729-15
Summary
Keywords
infectious virus, genome detection, oysters, persistence, shellfish outbreaks
Citation
Polo D, Schaeffer J, Teunis P, Buchet V and Le Guyader FS (2018) Infectivity and RNA Persistence of a Norovirus Surrogate, the Tulane Virus, in Oysters. Front. Microbiol. 9:716. doi: 10.3389/fmicb.2018.00716
Received
19 December 2017
Accepted
27 March 2018
Published
12 April 2018
Volume
9 - 2018
Edited by
Jean-Christophe Augustin, École Nationale Vétérinaire d’Alfort, France
Reviewed by
Sandra Martin-Latil, Agence Nationale de Sécurité Sanitaire de l’Alimentation, de l’Environnement et du Travail (ANSES), France; Dario De Medici, Istituto Superiore di Sanità, Italy
Updates
Copyright
© 2018 Polo, Schaeffer, Teunis, Buchet and Le Guyader.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Françoise S. Le Guyader, soizick.le.guyader@ifremer.fr; sleguyad@ifremer.fr
This article was submitted to Food Microbiology, a section of the journal Frontiers in Microbiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.