Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Pepper Mild Mottle Virus, a Plant Virus Associated with Specific Immune Responses, Fever, Abdominal Pains, and Pruritus in Humans

  • Philippe Colson,

    Affiliations Unité de Recherche sur les Maladies Infectieuses et Tropicales Émergentes (URMITE), Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR) 6236 – Institut de Recherche pour le Développement (IRD) 3R198, Facultés de Médecine et de Pharmacie, Université de la Méditerranée, Marseille, France, Pôle des Maladies Infectieuses et Tropicales Clinique et Biologique, Fédération de Bactériologie-Hygiène-Virologie, Centre Hospitalo-Universitaire Timone, Marseille, France

  • Hervé Richet,

    Affiliation Unité de Recherche sur les Maladies Infectieuses et Tropicales Émergentes (URMITE), Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR) 6236 – Institut de Recherche pour le Développement (IRD) 3R198, Facultés de Médecine et de Pharmacie, Université de la Méditerranée, Marseille, France

  • Christelle Desnues,

    Affiliation Unité de Recherche sur les Maladies Infectieuses et Tropicales Émergentes (URMITE), Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR) 6236 – Institut de Recherche pour le Développement (IRD) 3R198, Facultés de Médecine et de Pharmacie, Université de la Méditerranée, Marseille, France

  • Fanny Balique,

    Affiliations Unité de Recherche sur les Maladies Infectieuses et Tropicales Émergentes (URMITE), Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR) 6236 – Institut de Recherche pour le Développement (IRD) 3R198, Facultés de Médecine et de Pharmacie, Université de la Méditerranée, Marseille, France, Institut National de la Recherche Agronomique (INRA), Unité de Recherche (UR) 407, Unité de Pathologie Végétale, Montfavet, France

  • Valérie Moal,

    Affiliation Centre de Néphrologie et Transplantation Rénale, Centre Hospitalo-Universitaire Conception, Marseille, France

  • Jean-Jacques Grob,

    Affiliation Service de Dermatologie, Centre Hospitalo-Universitaire Sainte-Marguerite, Marseille, France

  • Philippe Berbis,

    Affiliation Service de Dermatologie, Centre Hospitalo-Universitaire Nord, Marseille, France

  • Hervé Lecoq,

    Affiliation Institut National de la Recherche Agronomique (INRA), Unité de Recherche (UR) 407, Unité de Pathologie Végétale, Montfavet, France

  • Jean-Robert Harlé,

    Affiliation Service de Médecine Interne, Centre Hospitalo-Universitaire Conception, Marseille, France

  • Yvon Berland,

    Affiliation Centre de Néphrologie et Transplantation Rénale, Centre Hospitalo-Universitaire Conception, Marseille, France

  • Didier Raoult

    didier.raoult@gmail.com

    Affiliations Unité de Recherche sur les Maladies Infectieuses et Tropicales Émergentes (URMITE), Centre National de la Recherche Scientifique (CNRS) Unité Mixte de Recherche (UMR) 6236 – Institut de Recherche pour le Développement (IRD) 3R198, Facultés de Médecine et de Pharmacie, Université de la Méditerranée, Marseille, France, Pôle des Maladies Infectieuses et Tropicales Clinique et Biologique, Fédération de Bactériologie-Hygiène-Virologie, Centre Hospitalo-Universitaire Timone, Marseille, France

Expression of Concern

This article [1] has been identified as one of a series of submissions for which we have concerns about the reported research ethics approval information and the article’s adherence to PLOS research ethics policies.

PLOS will be investigating these concerns in accordance with COPE guidance and journal policies. Meanwhile, the PLOS ONE Editors issue this Expression of Concern.

The journal has been informed that third author CD is deceased.

13 Dec 2022: The PLOS ONE Editors (2022) Expression of Concern: Pepper Mild Mottle Virus, a Plant Virus Associated with Specific Immune Responses, Fever, Abdominal Pains, and Pruritus in Humans. PLOS ONE 17(12): e0279046. https://doi.org/10.1371/journal.pone.0279046 View expression of concern

Abstract

Background

Recently, metagenomic studies have identified viable Pepper mild mottle virus (PMMoV), a plant virus, in the stool of healthy subjects. However, its source and role as pathogen have not been determined.

Methods and Findings

21 commercialized food products containing peppers, 357 stool samples from 304 adults and 208 stool samples from 137 children were tested for PMMoV using real-time PCR, sequencing, and electron microscopy. Anti-PMMoV IgM antibody testing was concurrently performed. A case-control study tested the association of biological and clinical symptoms with the presence of PMMoV in the stool. Twelve (57%) food products were positive for PMMoV RNA sequencing. Stool samples from twenty-two (7.2%) adults and one child (0.7%) were positive for PMMoV by real-time PCR. Positive cases were significantly more likely to have been sampled in Dermatology Units (p<10−6), to be seropositive for anti-PMMoV IgM antibodies (p = 0.026) and to be patients who exhibited fever, abdominal pains, and pruritus (p = 0.045, 0.038 and 0.046, respectively).

Conclusions

Our study identified a local source of PMMoV and linked the presence of PMMoV RNA in stool with a specific immune response and clinical symptoms. Although clinical symptoms may be imputable to another cofactor, including spicy food, our data suggest the possibility of a direct or indirect pathogenic role of plant viruses in humans.

Introduction

Little is known regarding the viral flora present in the human gastrointestinal tract [1][5]. Recently, viral communities have been examined in non-diarrheic human stool in a systematic and unbiased fashion using metagenomics [6][7]. Interestingly, the vast majority of the recovered viral sequences corresponded to uncultured pathogenic plant RNA viruses [7]. The most abundant viral sequences isolated were from the Pepper mild mottle virus (PMMoV), with up to 109 virions per gram of dry weight fecal matter. PMMoV is a non-enveloped, rod-shaped, single-stranded positive sense RNA virus classified in the genus Tobamovirus, which includes viruses extremely resistant to physical and chemical agents [8][9]. It is one of the major pathogens of Capsicum spp (chili peppers). Complementary data from Zhang et al.'s study have shown that PMMoV could be detected in non-diarrheic stool from 12 out of 18 individuals living in San Diego, USA or in Singapore, suggesting it might be geographically widespread, and in 3 out of 22 fresh and processed pepper samples. Moreover, the fecal PMMoV was viable and could infect host plants.

In the present study, we determined if PMMoV could be detected in commercial food products in our geographical area or the stool of patients and if these viruses were associated with biological and clinical symptoms.

Materials and Methods

Food products and human stool collection

Twenty-eight commercial food products were purchased in various grocery stores and supermarkets in Marseille (France) and stored at room temperature until processing (Table 1). They included four fresh peppers, thirteen manufactured products containing peppers, four products containing spice and seven products not known to contain pepper. For solid food products, 1 ml of the product was diluted in 1 ml of 1× PBS. Additionally, 565 stool samples were obtained between March and September 2008 for microbial diagnosis independent of the present study and approximately 1 ml collected randomly from each stool was stored in 1 ml of 1× PBS at +4°C. Of these samples, 357 were from adults (>18 years of age) and 208 were from children (<18 years of age). This study has been approved by our institutional ethics committee (Comité d'Ethique de l'IFR48, Université de la Méditerranée; review no. 09-001), that waived the need for consent.

thumbnail
Table 1. PMMoV RNA sequences recovered from food products.

https://doi.org/10.1371/journal.pone.0010041.t001

PMMoV RNA extraction

Viral RNA was extracted from 200 µl of the fecal supernatant or the solid food product suspensions following centrifugation at 10,000 g for 10 min or directly from the liquid food products by using the MagNA Pure LC RNA Isolation Kit (Roche Diagnostics, Meylan, France). The prepared RNA was directly analyzed or stored at −80°C until processing.

Detection of PMMoV RNA by real-time PCR

PMMoV RNA was detected using a real-time reverse transcription (RT)-PCR assay with SuperScript III Platinum One-Step Quantitative RT-PCR System (Invitrogen Life Technologies, Carlsbad, Calif., USA) on an Mx3000P thermocycler (Stratagene, La Jolla, CA 92037 USA). Each 25 µl reaction contained 10 µl of the extracted viral RNA, 12.5 µl 2× reaction mix, 0.5 µl of SuperScript III/RT Platinum Taq mix, and 0.75 µl (10 pmol/µl) of each primer and 0.5 µl (10 pmol/µl) of probe described by Zhang et al [7]. Positive samples were further tested by RT-PCR amplification and sequencing, as well as electron microscopy.

Infection of host plants with PMMoV RNA-positive food products

The viability and infectivity of PMMoV recovered from two PMMoV RNA-positive powder or dried food-products (no. 8, 17) and from a PMMoV RNA-positive Tabasco sample were tested by mechanical inoculation to host plants. Powder or dried food-products were crushed then mixed with an inoculation solution containing Na2HPO4 0.03M and 0.2% (wt/vol) Na-diethyldithiocarbamate. Twenty ml of Tabasco sauce were mixed with 80 ml of 0.5 M phosphate buffer (Na2HPO4-KH2PO4, pH 7.2), then filtered, and 8% (vol/vol) of n-butanol was added. After agitation for 15 min, centrifugation was performed at 10,000 g for 30 min at 10°C, and the supernatant was collected and filtered. Then, 4% (wt/vol) Polyethyleneglycol 6000 and 4% (wt/vol) NaCl were added prior to stirring for 15 min then centrifugation at 10,000 g and 10°C for 15 min. Pellets were resuspended in 1 ml of 0.01 M phosphate buffer. Samples were mixed with 600 mesh Carborundum and 200 µl were inoculated mechanically to three leaves of two plants of Nicotiana tabacum var Xanthi NN, a hypersensitive host for PMMoV. Negative controls were the same numbers of plants either not inoculated or inoculated with the buffer only. All plants were kept for two weeks in the same unit of an insect-proof greenhouse. Local lesions typical of PMMoV infection were excised and crushed into 0.6 ml of inoculation buffer, and 350 µl of this solution were tested for the presence of PMMoV by RT-PCR using primers published by Hamada et al. [10], by electron microscopy, and by inoculation to one plant of Capsicum annuum var Yolo Wonder or Nicotiana benthamiana, two susceptible hosts for PMMoV.

PMMoV RNA sequencing and phylogenetic analysis

PMMoV sequences corresponding to two different regions of the PMMoV genome were obtained by RT-PCR amplification using the SuperScript One-Step RT-PCR System (Invitrogen Life Technologies) and primers that had been designed with the SVARAP program [11]. The first targeted region is located near the 5′ end of the genome (primers used are PMMoV5pFwd: 5′-ATATCTGATGATGCAAGTTC-3′, 348 (the primer location was defined in reference to sequence GenBank accession nb NC_003630); PMMoV5p_Rev2: 5′-TAAACCTCTCTATTAGAGGC-3′, 928), and the second region corresponds to a fragment of the capsid gene (primers are PMMoVcaps_Fwd: 5′-CGTTAGGyAATCAGTTTCAA-3′, 5776; PMMoVcaps_Rev2: 5′-CGAACTAACTCATTCATGA-3′, 6070). For RT-PCR amplification, five µl of extracted viral RNA were added to 20 µl reaction solution containing 0.5 µl (10 pmol/µl) of primers, 12.5 µl of 2× reaction buffer, and 0.5 µl of RT/Platinum Taq mix. RT-PCR reaction was conducted under the following conditions: 45°C for 30 min for RT, 94°C for 2 min for initial denaturation, followed by 40 cycles of denaturation at 94°C for 30 s, annealing at 56°C for 45 s, elongation at 72°C for 2 min, then final elongation step at 72°C for 5 min. PCR products were purified with Sephadex G-50 Superfine on MAHVN 4550 plates (Millipore, Molsheim, France) then sequenced by using the amplification PCR primers and the Big Dye Terminator cycle sequencing kit version 1.1 on the ABI Prism 3130 genetic analyzer (Applied Biosystems, Branchburg, NJ, USA). The phylogenetic analysis was performed using Mega v.4.1 software (http://www.megasoftware.net/).

PMMoV purification from Tabasco sauce

Tabasco sauce (50 ml) was centrifuged at low speed (2,000 rpm) for 20 min. The PMMoV-containing supernatant (30 ml) was then treated with 1 M NaCl (30 ml) at 4°C overnight and filtered using a 5 µm Millex-SV filter (Millipore). The filtrate was loaded onto a caesium chloride (CsCl) density gradient consisting of 1 ml layers of 1.7, 1.5 and 1.30 g/ml CsCl in 1× PBS. After ultracentrifugation (22,000 g for 2 h at 4°C), the 1.30–1.50 g/ml fraction containing the purified PMMoV (ρ = 1.31) was collected. This fraction was washed on a Microcon 100 centrifugal filter (Millipore) with 1× PBS.

Electron microscopy study of PMMoV

A 40 µl aliquot of Tabasco sauce or fecal supernatant was deposited on a nickel 400 mesh Formvar carbon grid (EMS, Fort Washington, PA). After 15 min at 37°C, staining was performed with 1% ammonium molybdate for 10 sec. The stained grids were examined with a transmission electron microscope (Philips Morgagni 268D at 80Kv) after they were dried for 30 min at room temperature. The sensitivity of detection is estimated to be 106–108 virus particles/ml [12]. In addition, immunogold staining was performed, prior to electron microscopy, on PMMoV purified from PMMoV RNA-positive Tabasco sauce and stool, as follows: the electron microscopic grids were fixed with PBS and 1% glutaraldehyde for 2 min. After rinsing with PBS for 2 min, the grids were incubated twice for 5 min with 50 mM NH4Cl in PBS, preincubated twice for 2 min with Solution 1 (1% bovine serum albumine (BSA), 1% normal goat serum, 0.2% tween 20 in PBS) and incubated for 90 min with a polyclonal anti-PMMoV mouse antibody at 1∶1 dilution in Solution 1. These antibodies had been produced in two 6-week-old female Balb/C mice inoculated intraperitoneally with 1.3 µg of purified PMMoV using 400 µg aluminum hydroxide and 10 µg non-methylated oligonucleotides containing CpG motifs as an adjuvant, as previously described [13]. Blood collection was performed at day 38 following two booster doses administered at day 14 and day 28. After rinsing the stained grids with 0.1% BSA in PBS for 2 min, the grids were preincubated in Solution 2 (0.01% fish gelatin in PBS), transferred to a drop of gold-conjugated goat anti-mouse antibodies (diluted 1∶40 in Solution 2) for 60 min, washed six times with Solution 2 and twice with PBS for 10 min each. The grids were fixed in 2% glutaraldehyde for 15 min, then washed with PBS four times for 10 min and stained with 1% ammonium molybdate for 10 sec. Two negative controls were concurrently tested using antibodies derived from a non-immunized mouse or rotavirus instead of PMMoV.

Anti-PMMoV ELISA

ELISA Nunc Maxisorp plates were coated with 100 µl PMMoV solution per well (0.6 µg/ml purified from a PMMoV RNA-positive Tabasco sauce) overnight at 4°C. All wells were machine-washed thrice with 300µl 0.1% Tween in PBS and blocked with 200 µl of PBS containing 5% dried skimmed milk for 30 min at room temperature with gentle agitation. The wells were washed with 0.1% Tween in PBS, and 100 µl of serum diluted in 3% milk /0.1% Tween in PBS was added. The plates were incubated for 1 h at room temperature and washed as described above. After washing, 100 µl of horseradish peroxidase–conjugated goat anti-human total Ig or IgM (Beckman Coulter, Marseille, France) diluted 1∶5000 in 3% milk /0.1%Tween in PBS was added. The plates were then incubated for 1 hr at room temperature and washed as described above. The plates were developed by adding 200 µl of ortho-phenylenediamine per well and the absorbance was read at 492 nm.

Analysis of risk factors and biological and clinical signs associated with the presence of PMMoV RNA in stool or anti-PMMoV antibodies in serum

A case-control study was performed to assess the potential risk factors and the biological and clinical signs associated with the presence of PMMoV in the stool. Eighteen case-patients were compared to thirty-one control-patients with PMMoV RNA–negative stool. When possible, two control-patients from the same clinical units as the cases were selected for each case. The case- and control-patients were matched by age and sex. The variables analyzed included demographic characteristics, underlying diseases and clinical and biological signs. A Microsoft Excel spread sheet was used to enter the data. To identify risk factors for PMMoV detection, univariate analysis was performed with the Epi Info software version 3.3.2 (CDC, Atlanta, GA, USA). Proportions were compared by using the Yates Chi-square corrected test or the Fisher's exact test while continuous variables were compared by using either the ANOVA test or the Mann-Whitney/Wilcoxon two-sample test when the data were not normally distributed. For multivariate analysis, only variables with a p value ≤0.1 in univariate analysis were kept in the final model. Binary logistic regression analysis was performed using the SPSS software version 10 (Chicago, IL, USA).

Results

Detection of PMMoV RNA by real-time PCR

Food products.

PMMoV RNA sequences were recovered from twelve (57%) of the twenty-one pepper- or spice-containing food products (Table 1). Tabasco sauce contained the highest viral load, estimated to be nearly 107 PMMoV RNA copies/ml based on a cycle threshold of 22 by real-time PCR or the presence of more than one viral particle/field by electron microscopy (Figure 1a).

thumbnail
Figure 1. Electron microscopy of Pepper mild mottle virus (PMMoV)-RNA positive samples.

Footnote: The samples analyzed by electron microscopy were (a) Tabasco sauce; (b) PMMoV purified from Tabasco sauce. The samples analyzed following immunogold staining were: (c) PMMoV purified from Tabasco sauce; (d) a PMMoV-RNA positive patient's stool sample; (e) PMMoV purified from Tabasco sauce stained with antibodies from a non immunized mouse (negative control no.1) and (f) a rotavirus-positive sample (negative control no.2).

https://doi.org/10.1371/journal.pone.0010041.g001

Stool from children.

PMMoV RNA was detected in the stool from only 1 child out of 137 (0.7%) using real-time PCR (Table 2). The presence of PMMoV RNA was confirmed by sequencing and viral particles compatible with PMMoV were observed using electron microscopy. The patient was a 5-year-old boy living in Marseille who was admitted to our institution because of abdominal pains, fever, asthenia and hyperventilation.

thumbnail
Table 2. Characteristics of individuals positive and negative for real-time PCR PMMoV RNA detection in stools.

https://doi.org/10.1371/journal.pone.0010041.t002

Stool from adults.

The stool from 22 of 304 (7.2%) adults was positive for PMMoV by real-time PCR (Table 2). This proportion was significantly higher than the incidence in children (p = 0.0043). PMMoV RNA sequence could be recovered from all but two patients, whereas PMMoV-like viral particles were observed, by electron microscopy, in the stool from all but one patient, who was found positive for PMMoV RNA sequencing (Figure 1d) . The proportion of PMMoV RNA–positive patients was significantly higher in Dermatology Units than in other clinical units (7/20 (35.0%) versus 15/284 (5.3%); p = 7×10−7) (Table 2). Clinical and biological data were obtained from 18 of the 22 PMMoV-positive patients and compared to the characteristics of 31 patients with PMMoV-negative stool. By univariate analysis, abdominal pains, diverticulosis or diverticulitis, and fever were significantly associated with PMMoV RNA-positive stool (Table 3). No biological sign was significantly associated with the presence of PMMoV RNA. By multivariate analysis, using a model including age, gender and the variables significantly associated with the case-patients in the univariate analysis, fever (p = 0.0447 (relative risk (RR), 5.290 (95% confidence limits for RR (CI95), 1.040–26.911)), abdominal pains (p = 0.0376 (7.069 (1.119–44.660)), and pruritus (p = 0.0460 (6.131 (1.033–36.402)) were independently associated with the presence of PMMoV RNA. In addition, an interaction was observed between abdominal pains and diverticulitis. Overall, the model fitted the data well (Hosmer-Lemeshow test equal to 0.94). Three of four patients with concurrent abdominal pains and fever and the two patients who had abdominal pains and pruritus were PMMoV RNA-positive.

thumbnail
Table 3. Univariate analysis for risk factors and biological and clinical features in association with the detection by real-time PCR of PMMoV RNA from stools.

https://doi.org/10.1371/journal.pone.0010041.t003

Infection of host plants with PMMoV RNA-positive food products

All N. tabacum cultivar Xanthi NN plants inoculated with each of the three PMMoV RNA-positive food products developed local lesions typical of PMMoV infection within 5–7 days post-inoculation (Figures 2a–h). These lesions were observed on the three inoculated leaves for each plant, whereas control plants remained devoid of symptoms during an observation period of two weeks post-inoculation. PMMoV RNA was detected using RT-PCR from local lesions, PMMoV-like viral particles were observed by electron microscopy, and systemic mosaic symptoms typical of PMMoV were observed within 2–3 weeks post-inoculation on C. annuum cultivar Yolo Wonder or N. benthamiana.

thumbnail
Figure 2. Local lesions typical of PMMoV infection on Nicotiana tabacum cultivar Xanthi NN plants following inoculation of processed PMMoV RNA-positive food products.

Footnote: Leaves of Nicotiana tabacum cultivar Xanthi NN plants (a, e) non inoculated (negative control); (b, f) inoculated with buffer only (negative control); (c, and detail: d) inoculated with a processed PMMoV RNA-positive Tabasco sample; (g, and detail: h) inoculated with a PMMoV RNA-positive food product (no. 17; Table 1).

https://doi.org/10.1371/journal.pone.0010041.g002

PMMoV RNA sequencing and phylogenetic analysis

PMMoV sequences corresponding to the 5′ region of the PMMoV genome were obtained from four Tabasco sauce samples and thirteen patients, whereas sequences corresponding to a fragment of the capsid protein gene were obtained from eleven food products, including three Tabasco sauce samples, and eighteen patients (Figures 3a and 3b, respectively; GenBank accession numbers GQ427021–GQ427066). The BLAST analysis of the nucleotide sequences from the 5′ region produced only ten relevant hits, after which the sequences no longer corresponded to PMMoV (http://blast.ncbi.nlm.nih.gov/Blast.cgi). The sequences recovered from Tabasco sauce tended to cluster separately from those recovered from patients (Figure 3a). The mean nucleotide identity was 96.5±2.2% (range, 90.9–100%) for all PMMoV sequences, 96.4±2.4% (91.0–100%) for sequences recovered in the present study, 97.8±1.9% (93.9–100%) for sequences obtained from GenBank and 67.1±1.0% (64.1–68.6) with a Tobacco mosaic virus sequence. For the sequences of the capsid protein gene, the mean nucleotide identity was 97.1±1.8% (90.8–100%) for all PMMoV sequences, 97.2±1.9% (91.3–100%) for sequences recovered in the present study, 97.5±2.4% (90.8–100%) for sequences obtained from GenBank and 76.2±0.9% (73.4–77.7%) with Tropical soda apple mosaic virus sequence. One cluster (cluster II-1) contained four sequences from the present study (Figure 3b), and 100% nucleotide identity was observed between the sequences from two patients that received care in a Dermatology Unit and from a spicy powder.

thumbnail
Figure 3. Phylogenetic comparison of PMMoV nucleotide sequences recovered in the present study and corresponding (a) to a region located near the 5′ end of the genome or (b) to a region in the capsid gene to PMMoV sequences selected from GenBank.

Footnote: The phylogenetic trees were constructed using the neighbor-joining method based on the 5′ region (a) of the PMMoV genome or a region in the capsid gene (b). The PMMoV sequences recovered from patients' stool in the present study are in the boldface, white font and are indicated by a black square. Their name is labeled as follows: Mars_Laboratory identification number_Gender-Age_Clinical unit where the patient received care. The PMMoV sequences recovered from food products in the present study are in a boldface and are indicated by a gray square. Their name is labeled as follows: Mars_Name of the food product_Country where food product was manufactured. The remaining PMMoV sequences were obtained from GenBank. Their name is labeled as follows: GenBank Accession no._Coutry of origin_Year of submission to GenBank. Bootstrap values are indicated when greater than 50% as a percentage obtained from 1,000 resamplings of the data. The scale bar indicates the number of nucleotide substitutions per site. The sequence of Tobacco mild green mosaic virus (GenBank Accession No. NC_001556) and Tropical soda apple mosaic virus (GenBank Accession No. AY956381) were used as outgroups.

https://doi.org/10.1371/journal.pone.0010041.g003

Anti-PMMoV Total Ig and IgM ELISA

The level of serum anti-PMMoV total and IgM antibodies was examined by ELISA in a subset of 29 adult patients, 13 of which had PMMoV-positive stool, as detected by PCR (Figures 4a and 4b). For the analysis of the IgM anti-PMMoV level, the optical density (OD) for the substrate and antibody controls was 0.079 and 0.099, respectively. The mean OD was significantly higher in patients with PMMoV RNA than in the controls (0.472 versus 0.325, p = 0.013 (ANOVA); Figure 4b). Using the receiver operating characteristic (ROC) curve test with a sensitivity of 100% and a specificity of 62%, a value of 0.320 was considered positive for the detection of PMMoV-specific IgM (AUC = 0.808; 95% confidence interval: 0.645–0.971; p = 0.005). The proportion of serum samples with an OD greater than 0.320 was significantly higher for patients with PMMoV RNA-positive stool than for the controls (p = 0.0004). With this threshold, all 13 serum from patients with PCR PMMoV-positive stool were seropositive, whereas 10/16 (63%) sera from patients with PCR PMMoV-negative stool were seronegative. Thus, the positive and negative predictive values for having PCR-positive stool were 68% and 100%, respectively. For the analysis of the total anti-PMMoV level, the optical density (OD) for the substrate and antibody controls was 0.051 and 0.071, respectively. The mean OD was significantly higher in patients with PMMoV RNA than in the controls (0.663 versus 0.444, p = 0.028; Figure 4a). A cut-off value of 0.565 (estimated using the ROC curve test with a sensitivity of 69% and a specificity of 81%) was used (AUC = 0.724; 95% confidence interval: 0.529–0.918; p = 0.041). The proportion of serum samples for which the optical density (OD) was greater than 0.565 was significantly higher for patients with PMMoV-positive stool than for those with PMMoV-negative stool (p = 0.0070). Using this threshold, 9 (69%) of the 13 sera from patients with PCR PMMoV-positive stool were seropositive, whereas 13 (82%) of the 16 sera from patients with PCR PMMoV-negative stool were seronegative. Additionally, serological patterns indicating past infection (i.e., a positive total antibody result, but a negative IgM result) and negative total antibody and IgM testing were observed in two and eight PMMoV-negative patients, respectively, but in none of the PMMoV-positive patients (p = 0.0004). Finally, in the case-control study (Table 3), positive testing for anti-PMMoV IgM antibody was significantly more frequent in case-patients than in control-patients.

thumbnail
Figure 4. Comparison of the serum anti-PMMoV antibody level (a, total antibody; b, IgM) between individuals with PMMoV RNA-positive and PMMoV RNA-negative stool.

The dashed lines indicate the optical density cut-off values that were chosen for positive PMMoV serology using the receiver operating characteristic (ROC) curve test.

https://doi.org/10.1371/journal.pone.0010041.g004

Discussion

We detected PMMoV RNA by real-time PCR in stool from 22 of 304 (7.2%) adult patients and 1 of 208 children who received care in public hospitals. The present study therefore confirms the recent finding by Zhang et al. who identified PMMoV as the major RNA virus in human stool in a metagenomic study [7]. However, we found a much lower prevalence of PMMoV RNA-positive stool by real-time PCR among adult patients than Zhang et al., who detected the viral genome in 12 of 18 fecal samples (66.7%) using RT-PCR. This difference might be explained either by a lower sensitivity for real-time PCR or the lower exposure to PMMoV in our patients. We recovered PMMoV RNA sequences from 57% of food products containing pepper or spice and purchased from grocery stores in Marseille (southeastern France), indicating that these products represent a common source of viruses through ingestion. In addition, in food products and patient's stool, we could visualize by electron microscopy complete viruses whose shape and size are consistent with PMMoV and that reacted with PMMoV-specific antibodies in the immunogold staining experiments. Moreover, 46 new PMMoV sequences corresponding to two regions of the viral genome were recovered from food products or patients' stool. Four phylogenetic clusters were identified, and one of them (II-1) included sequences from both patients and a putative food source. Furthermore, three PMMoV RNA-positive food products were found to contain virus particles infectious to host plants.

We also identified statistically significant differences in the occurrence of fever, abdominal pains, and pruritus and the detection of specific immune responses to PMMoV in the case-control study. We, therefore, believe that we provide the first evidence that plant viruses may cause disease in humans. Fever and abdominal pains, which were significantly more common in PMMoV-positive patients (39% for both clinical features) than in patients with PMMoV-negative stool (13% and 7%, respectively), were also present in the only PMMoV RNA-positive child. These symptoms might not be caused by PMMoV but by consumption of spicy food. Otherwise, they may be fortuitously related to an associated unknown infectious agent or to an allergic or toxic cofactor. Nevertheless, the presence of a specific immune response to PMMoV reinforces the hypothesis that this plant virus was the cause of the symptoms in these patients. Indeed, anti-PMMoV IgM antibodies were detected in all PMMoV-positive patients, significantly more frequently than in the controls. The concurrent detection of PMMoV RNA and antibodies to PMMoV indicates that this virus is not a neutral component of the human gut flora. On the contrary, our findings suggest that PMMoV may infect humans and cause clinical symptoms. We believe that these preliminary data may prompt further studies on the pathogenic role of PMMoV in humans. Finally, the higher prevalence of PMMoV RNA among patients who received care in Dermatology Units deserves further attention, as it was high (35%).

To date, plant viruses and vertebrate viruses are believed to exist in two distinct and non overlapping biological niches, and they are not known to share common hosts in the plant or animal kingdoms [14][15]. Therefore, to date, plant viruses have not been described as pathogens for vertebrates or humans or known even to infect them [16]. This paradigm is underscored by the growing body of data about the use of plant virus–based vaccines, which are considered as safe tools that might be useful in treating chronic viral infections and cancer in the future [17][20]. However, the boundaries between plants and animals might not be so hermetical for plant viruses, and, to our knowledge, no experimental model has been described that would allow us to rule out a pathogenic role for plant viruses in vertebrates. The recent work of Zhang et al. challenged this dogma, as it was neutral, technology-driven research that was not trying to confirm a former hypothesis [7], [21]. That work identified plant viruses as the most common viral inhabitants of the human gut flora, PMMoV being the most prevalent. In addition, a recent study demonstrated that PMMoV is widespread in wastewater from the United States of America (USA), at concentrations ranging between 8.0×105 and 2.2×107 copies/ml [22]. Beyond these findings, our results, together with previous data, suggest that the relationship between plant viruses and animals may deserve further re-evaluation using modern microbiologic tools (Table 4). The oral administration of Cowpea severe mosaic virus, another plant virus, has been previously shown to induce a durable and systemic immune response in mice without requiring the co-administration of an adjuvant [23], and similar results have been observed in mice using Alfalfa mosaic virus and chimeric plant virus particles [24][25]. In addition, some plant viruses can replicate in insect tissues, as illustrated by Tomato spotted wilt virus and Maize rayado fino virus [26][28]. Furthermore, a reduction in the lifespan and fecundity of insects (Bemisia tabaci) during infection with Tomato yellow leaf curl virus has been previously observed, indicating that this plant virus had some features of an insect pathogen [29]. It has been also suggested, based on viral sequence analysis, that during evolution, plant-infecting viruses (previously known as plant circoviruses and recently reclassified in the genus Nanovirus) switched hosts to infect a vertebrate and then recombined with a vertebrate-infecting virus [16]. Moreover, it has been shown recently that Tomato spotted wilt virus could replicate in two human cell lines, HeLa and diploid fibroblasts [30]. This virus is one of the most important plant pathogens worlwide, and it causes mild infection on its main insect vector, Frankliniella occidentalis, in which it induces a strong immune response [26]. Finally, the possible involvement of plant viruses in human diseases was questioned several decades ago. Indeed, culturable Tobacco mosaic virus could be recovered from sputum and thoracentesis fluids obtained from cigarette smokers with a history of pulmonary disease, including lung cancer [31][33]. Taken together, these experimental data, coupled with the previous results reported by Zhang et al. and our findings, favor the hypothesis of a pathogenic role for PMMoV. It should be highlighted that PMMoV, like other tobamoviruses, are known to be extraordinarily resistant to heat and desiccation [34][35]. Moreover, although this deserves to be investigated in further studies, the viability and infectivity to host plants were evidenced for PMMoV recovered from food products in the present study and from stools in Zhang et al.'s study [7]. Additionally, we visualized complete viral particles in stool by electron microscopy.

thumbnail
Table 4. Summary of postulates and findings about interactions of PMMoV with the host immune system, its replication in non-plant tissues, and its potential link with clinical signs.

https://doi.org/10.1371/journal.pone.0010041.t004

In summary, we confirmed that humans may carry a high PMMoV load, likely acquired from food products, and we documented that PMMoV might not only be a common inhabitant of the human gut but may also interact with the human immune system and cause clinical symptoms. These results should prompt further studies to re-evaluate whether or not plant viruses, including PMMoV, may have a pathogenic role in humans.

Acknowledgments

We are thankful to Claude Nappez, Bernard Campagna, Audrey Borg, Audrey Ferretti, Mamadou Kaba, and Catherine Rys for their technical assistance.

Author Contributions

Conceived and designed the experiments: PC HL DR. Performed the experiments: CD FB. Analyzed the data: PC HR FB HL DR. Contributed reagents/materials/analysis tools: CD FB VM JJG PB HL JRH YB DR. Wrote the paper: PC HR CD HL DR.

References

  1. 1. Finkbeiner SR, Allred AF, Tarr PI, Klein EJ, Kirkwood CD, et al. (2008) Metagenomic analysis of human diarrhea: viral detection and discovery. PLoS Pathog 4(2): e1000011.
  2. 2. Denno DM, Klein EJ, Young VB, Fox JG, Wang D, et al. (2007) Explaining unexplained diarrhea and associating risks and infections. Anim Health Res Rev 8(1): 69–80.
  3. 3. Klein EJ, Boster DR, Stapp JR, Wells JG, Qin X, et al. (2006) Diarrhea etiology in a Children's Hospital Emergency Department: a prospective cohort study. Clin Infect Dis 43: 807–13.
  4. 4. Musher DM, Musher BL (2004) Contagious acute gastrointestinal infections. N Engl J Med 351(23): 2417–27.
  5. 5. Chiki-Brachet R, Bon F, Toubiana L, Pothier P, Nicolas JC, et al. (2002) Virus diversity in a winter epidemic of acute diarrhea in France. J Clin Microbiol 40(11): 4266–72.
  6. 6. Breitbart M, Hewson I, Felts B, Mahaffy JM, Nulton J, et al. (2003) Metagenomic analyses of an uncultured viral community from human faeces. J Bacteriol 185(20): 6220–3.
  7. 7. Zhang T, Breitbart M, Lee WH, Run JQ, Wei CL, et al. (2006) RNA viral community in human faeces: prevalence of plant pathogenic viruses. PLoS Biol 4(1): e3.
  8. 8. Anonymous. ICTVdb Management. 00.071.0.01.007 (2006) Pepper mild mottle virus. In: Büchen-Osmond C, editor. ICTVdb – The Universal Virus Database, version 4. New York, USA: Columbia University.
  9. 9. Wetter C, Conti M, Altschuh D, Tabillion R, van Regenmortel MHV (1984) Pepper Mild Mottle Virus, a Tobamovirus infecting pepper cultivars in Sicily. Phytopathology 74(4): 405–10.
  10. 10. Hamada H, Tomita R, Iwadate Y, Kobayashi K, Munemura I, Takeuchi S, et al. (2007) Cooperative effect of two amino acid mutations in the coat protein of Pepper mild mottle virus overcomes L3-mediated resistance in Capsicum plants. Virus Genes 34: 205–214.
  11. 11. Colson P, Tamalet C, Raoult D (2006) SVARAP and aSVARAP: simple tools for quantitative analysis of nucleotide and amino acid variability and primer selection for clinical microbiology. BMC Microbiol 6: 21.
  12. 12. Curry A, Appleton H, Dowsett B (2006) Application of transmission electron microscopy to the clinical study of viral and bacterial infections: present and future. Micron 37(2): 91–106.
  13. 13. Near KA, Stowers AW, Jankovic D, Kaslow DC (2002) Improved immunogenicity and efficacy of the recombinant 19-kilodalton merozoite surface protein 1 by the addition of oligodeoxynucleotide and aluminium hydroxyde gel in a murine malaria vaccine model. Infect Immun 70: 692–701.
  14. 14. Fauquet CM, Mayo MA, Maniloff J, Desselberger U, Ball LA, editors. (2005) Virus Taxonomy, VIIIth Report of the ICTV. London, UK: Elsevier/Academic. 1259 p.
  15. 15. Jones RA (2009) Plant virus emergence and evolution: Origins, new encounter scenarios, factors driving emergence, effects of changing world conditions, and prospects for control. Virus Res 141(2): 113–130.
  16. 16. Gibbs MJ, Weiller GF (1999) Evidence that a plant virus switched hosts to infect a vertebrate and then recombined with a vertebrate-infecting virus. Proc Natl Acad Sci U S A 96(14): 8022–7.
  17. 17. Yusibov V, Rabindran S, Commandeur U, Twyman RM, Fischer R (2006) The potential of plant virus vectors for vaccine production. Drugs R D 7(4): 203–17.
  18. 18. Nochi T, Takagi H, Yuki Y, Yang L, Masumura T, et al. (2007) Rice-based mucosal vaccine as a global strategy for cold-chain- and needle-free vaccination. Proc Natl Acad Sci U S A104(26): 10986–91.
  19. 19. Pascual DW (2007) Vaccines are for dinner. Proc Natl Acad Sci U S A 104(26): 10757–8.
  20. 20. Lacasse P, Denis J, Lapointe R, Leclerc D, Lamarre A (2008) Novel plant virus-based vaccine induces protective cytotoxic T-lymphocyte-mediated antiviral immunity through dendritic cell maturation. J Virol 82(2): 785–94.
  21. 21. Raoult D (2010) Technology-driven research will dominate hypothesis-driven research: the future of microbiology. Future Microbiol 5(2): 135–7.
  22. 22. Rosario K, Symonds EM, Sinigalliano C, Stewart J, Breitbart M (2009) Pepper mild mottle virus as an Indicator of Fecal Pollution. Appl Environ Microbiol [Epub ahead of print].
  23. 23. Florindo MI, de Aragao ME, da Silva AC, Otoch ML, Melo DF, et al. (2002) Immune response induced in mice oral immunization with cowpea severe mosaic virus. Braz J Med Biol Res 35(7): 827–35.
  24. 24. Modelska A, Dietzschold B, Sleysh N, Fu ZF, Steplewski K, et al. (1998) Immunization against rabies with plant-derived antigen. Proc Natl Acad Sci U S A 95(5): 2481–5.
  25. 25. Brennan FR, Bellaby T, Helliwell SM, Jones TD, Kamstrup S, et al. (1999) Chimeric plant virus particles administered nasally or orally induce systemic and mucosal immune responses in mice. J Virol 73(2): 930–8.
  26. 26. Medeiros RB, Resende Rde O, de Avila AC (2004) The plant virus Tomato Spotted Wilt Tospovirus activates the immune system of its main insect vector, Frankliniella occidentalis. J Virol 78(10): 4976–82.
  27. 27. Wijkamp I, van Lent J, Kormelink R, Goldbach R, Peters D (1993) Multiplication of tomato spotted wilt virus in its insect vector, Frankliniella occidentalis. J Gen Virol 74(Pt 3): 341–9.
  28. 28. Kitajima EW, Gámez R (1983) Electron microscopy of maize rayado fino virus in the internal organs of its leafhopper vector. Intervirology 19(3): 129–34.
  29. 29. Rubinstein G, Czosnek H (1997) Long-term association of tomato yellow leaf curl virus with its whitefly vector Bemisia tabaci: effect on the insect transmission capacity, longevity and fecundity. J Gen Virol 78(Pt10): 2683–9.
  30. 30. de Medeiros RB, Figueiredo J, de O Resende R, De Avila AC (2005) Expression of a viral polymerase-bound host factor turns human cell lines permissive to a plant- and insect-infecting virus. Proc Natl Acad Sci U S A 102(4): 1175–80.
  31. 31. Bothwell PW (1960) Lung cancer and tobacco mosaic virus. Lancet 1: 657.
  32. 32. Katsilambros L (1960) Tobacco mosaic virus and lung cancer. Lancet 2: 934.
  33. 33. LeClair RA (1967) Recovery of culturable tobacco mosaic virus from sputum and thoracentesis fluids obtained from cigarette smokers with a history of pulmonary disease. Am Rev Respir Dis 95(3): 510–1.
  34. 34. Harrison BD, Wilson TM (1999) Milestones in the research on tobacco mosaic virus. Philos Trans R Soc Lond B Biol Sci 354(1383): 521–9.
  35. 35. Rogers SO, Starmer WT, Castello JD (2004) Recycling of pathogenic microbes through survival in ice. Med Hypotheses 63(5): 773–7.
  36. 36. Hwang DJ, Roberts IM, Wilson TM (1994) Expression of tobacco mosaic virus coat protein and assembly of pseudovirus particles in Escherichia coli. Proc Natl Acad Sci U S A 91(19): 9067–71.