- Open Access
Enterobacteriaceae dominate the core microbiome and contribute to the resistome of arugula (Eruca sativa Mill.)
- Tomislav Cernava†1,
- Armin Erlacher†1, 3,
- Jung Soh2,
- Christoph W. Sensen2, 4,
- Martin Grube3 and
- Gabriele Berg1Email authorView ORCID ID profile
© The Author(s). 2019
- Received: 24 May 2018
- Accepted: 10 January 2019
- Published: 29 January 2019
Arugula is a traditional medicinal plant and popular leafy green today. It is mainly consumed raw in the Western cuisine and known to contain various bioactive secondary metabolites. However, arugula has been also associated with high-profile outbreaks causing severe food-borne human diseases. A multiphasic approach integrating data from metagenomics, amplicon sequencing, and arugula-derived bacterial cultures was employed to understand the specificity of the indigenous microbiome and resistome of the edible plant parts.
Our results indicate that arugula is colonized by a diverse, plant habitat-specific microbiota. The indigenous phyllosphere bacterial community was shown to be dominated by Enterobacteriaceae, which are well-equipped with various antibiotic resistances. Unexpectedly, the prevalence of specific resistance mechanisms targeting therapeutic antibiotics (fluoroquinolone, chloramphenicol, phenicol, macrolide, aminocoumarin) was only surpassed by efflux pump assignments.
Enterobacteria, being core microbiome members of arugula, have a substantial implication in the overall resistome. Detailed insights into the natural occurrence of antibiotic resistances in arugula-associated microorganisms showed that the plant is a hotspot for distinctive defense mechanisms. The specific functioning of microorganisms in this unusual ecosystem provides a unique model to study antibiotic resistances in an ecological context.
Plants are generally colonized by a vastly diverse microbiota, which contributes to the host plant’s health and productivity [40, 46]. However, the structure, abundance, and occurrence of microorganisms can be variable across different plant habitats [51, 62] and plant species/genotypes [9, 15]. All members of the Brassicaceae plant family characterized thus far have shown an extraordinary degree of specificity in their associated microbiota [6, 7, 65]. They are known for a bacteria-dominated composition of the microbiome, as they do not harbor a typical mycorrhiza [30, 54]. The effective defense mechanisms of the plant are based on a unique composition of antimicrobial secondary metabolites, including glucosinolates [42, 52]. Their resulting degradation products (isothiocyanates) have also received scientific interest, as they might be involved in cancer prevention . Arugula (Eruca sativa Mill.) is the only salad plant from this family, which is quite popular for human consumption. Commonly known as rucola (or garden rocket, respectively), it has its origin in the Mediterranean and has been cultivated at least since Roman times [43, 68]. The unique aroma is caused mainly by the aforementioned isothiocyanates and butane derivates , respectively. In addition, arugula also contains a number of other health-promoting agents . While the health-promoting properties of arugula are undisputed , almost nothing is known about the associated microbiota. Our hypothesis was that arugula would be associated with a highly specialized microbiota, due to metabolites it produces, which were expected to be highly different from salad species from the Asteraceae family, which have been well studied previously [18, 53].
Arugula has been reported to be the potential source of pathogen outbreaks, especially associated with Enterobacteriaceae [32, 50]. Only in recent years, Enterobacteriaceae have been confirmed as indigenous components of the plant microbiome in other species [12, 26, 27, 60]. These pathogenic bacteria come also into focus for their potential role in the transfer of antibiotic resistance by mobile genetic elements [28, 34, 58]. The environmental resistome has now been recognized as the origin and reservoir of antibiotic resistance genes and is considered to be dynamic and ever-expanding . With increasing resistances, the success of treatment options for infectious diseases, e.g., cases of sepsis, after major surgery or during cancer chemotherapy will be compromised [20, 22]. Wellington et al.  discussed an extensive antibiotic resistance (AR) “pollution” in our environment. In a large-scale metagenomics-based study of antibiotic resistances in the environment, Nesme et al.  depicted soil habitats as the most diverse reservoir for AR across four different environmental clusters. Lettuce was also reported as a source of multi-resistant bacteria [17, 64], but such information thus far lacks for arugula. Therefore, we were interested in the potential for antibiotic resistance in this unique ecosystem and hypothesized that the arugula resistome may be affected by its particular microbiome.
We have investigated the arugula phyllosphere in comparison with the rhizosphere and the surrounding bulk soil, using whole metagenome sequencing and functional annotation of the resulting genomic sequences. In addition, complementary amplicon sequencing to further characterize the Enterobacteriaceae populations was performed in order to complete the integrative study of the arugula microbiome.
Sampling of arugula plants
Arugula plants were grown in garden soil (hereafter referred to as bulk soil) in a suburban region of Graz (Austria; approx. 47°4′ 13′′ N 15° 28′ 19″ E). All samples used are summarized in Additional file 1: Table S1. Edible plant parts consist of elongated leafs with short stalks and are termed as phyllosphere throughout this study when referring to the microbial habitat. The microbial fraction of the Eruca sativa system was extracted separately for both the phyllo- and rhizosphere in two consecutive samplings. During the first sampling, arugula plants were collected to assess the enterobacterial microbiome by amplicon sequencing of 16S rRNA gene fragments. Plants were harvested in their final growth stage of leaf development in July. During the follow-up sampling, samples from the same garden were obtained in early November at their final growth stage when the plants had already formed flowers and seed bodies. They were grown in a warm summer with continental or hemiboreal climate, and the recorded climate profiles (Additional file 1: Figure S1), which started 2 weeks before sampling, show relatively warm temperatures, moderate but fluctuating precipitation, sunshine and atmospheric pressure levels, overall low wind levels, and average air humidity. In addition to the phyllo- and rhizosphere samples, bulk soil was sampled for metagenome sequencing. All samples were stored on ice and immediately processed after arrival at a nearby laboratory.
Amplicon sequencing and processing of enterobacterial 16S rRNA gene fragments
The microbial fractions were extracted separately for the rhizo- and phyllospere, and three independent single replicates per habitat, consisting of 15–20 leaves or roots, respectively, were collected and stored separately. For the homogenization of the samples, 5 g of plant material per replicate was physically disrupted with a sterile pestle and mortar, re-suspended in 10 ml of 0.85% NaCl, transferred in two 2 ml Eppendorf tubes, and subsequently centrifuged (16,500 g, 20 min, 4 °C). The pellet obtained was used for isolation of the total-community DNA with the FastDNA® SPIN Kit for Soil (MP Biomedicals, Solon, OH, USA). For mechanical lysis, the cells were homogenized twice in a FastPrep® FP120 Instrument (Qbiogene, BIO101, Carlsbad, CA, USA) for 30 s at a speed of 5.0 m s−1 and treated according to the manufacturer’s protocol.
Amplification and sequencing of 16S rRNA fragments
The 16S rRNA gene fragments of the Enterobacteriaceae were amplified in a dual-phase nested PCR approach using multiplex identifier (MID) tagged primers. The PCR was conducted with specific primers for enterics, using Entero-F234 and Entero-R1423 (V3-V8), according to the method described by Binh et al. . The nested PCR was carried out as described by Heuer et al.  using the primer pair F984 and R1378 modified for multiplex 454 sequencing according to the specification. The products of two independent PCR reactions per sample were pooled and purified using the Wizard SV Gel and PCR Clean-Up System (Promega, Madison, USA). Purified PCR products were pooled (200 ng each) and sequenced on a Roche GS FLX+ 454 Titanium platform (Macrogen Korea, Seoul, South Korea).
Bioinformatics sequence processing
Sequences were analyzed with the QIIME software version 1.8.0 . MID, primer, and adapter sequences were removed, and the sequences were quality- (minimal score 50) and length filtered (minimal raw fragment length: E1 ≥ 350; E2 and E3 ≥ 430), followed by a denoising step (using the denoise wrapper and denoiser script; both contained within the QIIME pipeline). Chimeras and remaining sequences of non-target, plastidal, and mitochondrial origin were removed. OTU tables were created with UCLUST  at a 100% cut-off level. Unique reads were classified with the RDP classifier  v2.5 and the NCBI RefSeq database. The datasets were rarified on a naїve Bayesian classifier  to the number of least reads within each experiment, to compute alpha and beta diversity indices. An OTU-based network correlating the OTUs at 100% cut-off level was constructed with Cytoscape v.3.1.0 . Centiscape v.2.2  was used to calculate centroid values for each node of the network and label them accordingly.
Metagenome shotgun sequencing and bioinformatic processing
Five grams of each sample (phyllo-, rhizosphere, and bulk soil) were weighed, transferred into sterile plastic bags together with 10 ml 0.85% NaCl, and mechanically processed twice for 210 s in a Bagmixer (Interscience, St. Nom, France). Samples were placed at interims of 5 min on ice. Homogenized cell suspensions were further transferred into S34 tubes, and centrifuged at 10,000 rpm for 20 min. The supernatant was discarded for each approach and the cell pellets were stored at − 20 °C. In total, we processed 24 × S34 tubes (35 ml/S34) foliage samples from a total 90 bags (cell pellets were stacked three times and the supernatant discarded between each centrifugation round in order to increase the total yield), 8 × S34 tubes (35 ml/S34) rhizosphere samples from a total of 35 bags, and 6 × S34 tubes (35 ml/S34) soil samples from a total of 35 bags, respectively. DNA was extracted using the FastDNA® SPIN Kit for Soil (MP Biomedicals, Solon, OH, USA) from six pellets per habitat (foliage, rhizosphere, soil) using 300 mg each of the crude cell pellets and processing as stated in the manufacturer’s protocol. The DNA samples were eluted in 100 μL H2Oultra pure and checked for quality and quantity using NanoDrop Photometer. Then, 2 μg of each replicate were pooled for a total 12 μg of metagenomic DNA per habitat.
Shotgun sequencing of total community DNA
Metagenomic sequencing was performed on an Illumina HiSeq2000 system (2 × 150 bp) by Eurofins MWG Operon (Ebersberg, Germany) following the Eurofins MWG Operon protocol.
Sequence quality filtering
The raw read data were uploaded to the Galaxy Main server (http://usegalaxy.org) and processed using the FASTQ Joiner, FASTQ Groomer, FASTQ Quality Trimmer, and Filter FASTQ (min. length = 75, min. quality = 20, allowed outside quality range = 2). For each of the datasets, only high-quality reads were retained for further processing, resulting in the following read numbers to process further: phyllosphere: 41,867,724 (79.8%), rhizosphere: 35,463,395 (74.7%), and soil: 27,085,866 (79.3%).
Reads were aligned to the reference genomes of Brassica oleracea (BOL; ncbi.nlm.nih.gov/genome), Brassica rapa (Brapa_1.0), and Raphanus sativus (Rs1.0), respectively, with the goal of discarding aligned reads in order to reduce plant-derived genetic bias. To align the quality-filtered reads to the reference genomes, Bowtie 2  was used, with the default alignment parameters. Additional file 1: Table S2 shows the alignment results. Those reads that aligned to at least one of the three genomes were excluded from further processing. Thus, only those reads that were not aligned to any of the three reference genomes were used as input to the assembly program. The phyllosphere and soil samples had the largest and smallest numbers of aligned reads among the three samples, respectively, which was to be expected.
Assembly of the reads into contigs was performed using the Velvet (version 1.2.10; ) de novo assembly software. Since the choice of the k-mer length is crucial, multiple assemblies with different k values were performed on each sample in order to determine the best k-mer length. In all assemblies, the expected coverage value was automatically found by Velvet and the insert length for paired-end reads was set at 350. Additional file 1: Table S3 shows the key statistics of the final assemblies, where k is the length of the k-mer used. Only those contigs that had at least 90% of the median coverage and were at least as long as twice the k value were selected for subsequent analysis.
Taxonomic and functional analyses
For each sample, taxonomy and functional analyses were performed using the MEGAN program (version 5.11.3; ). For taxonomic analyses, the blastn program was run on all the filtered contigs on the TimeLogic DeCypher boards (www.timelogic.com), against the “nt” database with 34,646,553 nucleotide sequences downloaded from NCBI on February 12, 2016. In order to create input for functional analyses, the blastx (nucleotide-to-protein) program was run on all the filtered contigs on the TimeLogic DeCypher boards, against the “nr” database with 81,622,391 protein sequences downloaded from NCBI on February 12, 2016. The blastn results from each sample were imported into MEGAN to perform the taxonomy analysis for the sample, where each contig is assigned to a bacterial lineage at the lowest rank possible with sufficient confidence. An identical number of contigs assigned to bacteria (the smallest such number among the phyllosphere, rhizosphere and bulk soil samples) were randomly extracted from each sample and used to compare subsequent bacterial taxonomic and functional analyses results among the three samples. Therefore, the taxonomy and functional analyses of the bacteria are based on the blastn and blastx results, respectively, for the 89,462 subsampled contigs from each sample. Similarly, the taxonomy and functional analyses of Enterobacteriaceae are based on the blastn and blastx results, respectively, for the 5494 subsampled contigs from the phyllosphere and rhizosphere samples, respectively. A very low number (356) of contigs from the soil sample was assigned to the Enterobacteriaceae; therefore, we were not able to include them in a meaningful comparison.
Antibiotic resistance analysis
Microbial antibiotic resistance profiles were analyzed using the Comprehensive Antibiotic Resistance Database (CARD; ). The antibiotic resistance protein sequences for metagenomic data were downloaded from the CARD website (card.mcmaster.ca) as a FASTA file and used as the target sequences for the nucleotide-to-protein Blast (blastx) analysis of the assembled contigs for each sample. The analysis was performed by using only the top blast hits, with a percentage identity of at least 50%. A single hit of a contig to an antibiotic resistance protein was counted as multiple antibiotic resistance category hits, if the protein mapped to multiple categories.
MG-RAST analysis of non-assembled reads
Abundances within the enterobacterial fraction were analyzed in a separate bioinformatics approach based on non-assembled metagenome sequences. The metagenomes (phyllosphere, rhizosphere, and bulk soil) were uploaded on the Metagenomic RAST (MG-RAST) server  and initially processed using the default parameters. The processing included the removal of artificial replicate sequences , low-quality sequences , short sequences, and sequences containing ambiguous bases. Phylogenetic annotations were obtained from the RefSeq database, with a maximum e value of 10−5, minimum identity of 60%, and minimum alignment length of 15 bp for RNA databases.
Isolation of a strain collection and antibiotic susceptibility tests
For the isolation of enterobacteria, homogenized cell suspensions from the second sampling were plated on growth media. In total, 180 aerobic bacteria were isolated from the arugula phyllosphere and rhizosphere as well as bulk soil. Different media including PDA, MacConkey, KingsB, R2A (all obtained from Carl Roth GmbH, Karlsruhe, Germany), MIS (prepared according to ), and NAII (Sifin, Berlin, Germany) allowed a broad-spectrum coverage to create a representative strain library. None of the media used for the initial isolation was supplemented with antibiotics, in order to avoid an isolation bias caused by selective conditions. A disk diffusion antimicrobial susceptibility test method  was used to test for susceptibility towards a defined set of antibiotics. Instead of MH agar, NAII agar was used for the test plates. Unique morphotypes were tested for their sensitivity against Ampicillin (10 μg), Chloramphenicol (30 μg), Erythromycin (15 μg), Gentamicin (10 μg), Penicillin G (10 Units; 6 μg), Polymyxin B (300 Units; 30 μg), Streptomycin (10 μg), and Tetracycline (30 μg).
The arugula microbiome and the community structure of the enterobacterial fraction
The enterobacterial core community
Antibiotic resistances in the arugula metagenome
Antibiotic resistance analysis with the Kirby Bauer testing system
A strain library with 180 isolates was constructed following the isolation approach with various cultivation media and samples from the two plant habitats. This strain library was subjected to the standardized testing system (Kirby Bauer; ). A high antibiotic resistance profile was observed against most tested antibiotics, which comprised Ampicillin, Chloramphenicol, Erythromycin, Gentamicin, Penicillin G, Polymyxin B, Streptomycin, and Tetracycline. Isolates derived from the rhizosphere and phyllosphere showed higher antibiotic resistance per strain when compared to the soil isolates. More than 50% of the strain library subjected to the antibiotics showed resistance at the tested concentrations. Polymyxin was the only antibiotic with resistance abundance lower than 50%. Eleven strains showed multi resistance toward all subjected antibiotics (Additional file 1: Table S4). These strains were further sequenced and BLAST-searched against the 16S rRNA and nt databases of NCBI. Eight strains belong to the group of Gammaproteobacteria (including three Enterobacteriaceae members), while two were assigned to Firmicutes and one to Bacteroidetes. The genera Pseudomonas (five isolates) and Erwinina (two isolates) were repeatedly identified, while the other resistant isolates belonged to different bacterial lineages. When the isolates were compared at species level (Additional file 1: Table S4), no overlaps in taxonomic assignments between the three isolate sources were observed.
Our integrative profiling of the microbiota in arugula provided evidence for a specific microbiome composition with an enrichment of members of the Enterobacteriaceae especially in the edible plant parts. Their representatives, e.g., Pantoea, Erwinia, Pectobacterium, Serratia, and Yersinia, are well-equipped with various antibiotic resistances and thus contribute to the overall resistome. These findings confirm our initial hypotheses related to the indigenous arugula microbiome. Unexpectedly, the prevalence of specific resistance mechanisms targeting therapeutic antibiotics in plant-associated microorganisms was only surpassed by the number of identified genes coding for efflux pumps. High frequencies for fluoroquinolone, chloramphenicol, phenicol, macrolide, and aminocoumarin resistances were detected.
Naturally occurring Enterobacteriaceae members dominate arugula’s phyllosphere
The predominance of enterics is in good agreement with previous findings targeting the lettuce (Lactuca sativa L.) microbiome [12, 26, 27, 60]. The Enterobacteriaceae are a large family of Gram-negative bacteria that are excellent degraders of organic compounds. This family includes many harmless symbionts, serious pathogens, as well as those of ambivalent character, e.g., Pantoea and Erwinia ( [23, 56]), which were among the predominant lineages in the phyllosphere metagenome, as well as in the amplicon datasets dominant in arugula. They were also part of the enterobacterial core-biome and thus likely contribute to the host’s natural biota. Interpretation of the network analyses suggests that Enterobacteriaceae with a high relative abundance are not the only opportunistic niche occupiers on the plants. According to our data, it is likely that the soil provides the major reservoir from which enterics can establish and spread within the whole plant. It has already been shown that the leaf age, plant lesions, and nitrogen content affect the bacterial communities and promote the proliferation of enterics [13, 14, 27]. We assume that the arugula phyllobiome favors bacteria that can cope well with the rapidly expanding habitat of the growing leaves. Similar to lettuce or spinach, which also contains high amounts of enterobacteria, leaves of arugula have a similar structure. This provides a different habitat compared to the phyllosphere of other plants, which endorse structurally different microbiota . One aspect that might contribute to the survival of the enterobacterial community could be the tolerance against a broad range of plant-specific metabolites  that may be involved in the regulation of the microbial community. The presence of bioactive compounds may shrink the possibilities of other microorganisms to rapidly establish in affected habitats.
Antibiotic resistances are an additional risk factor associated with opportunistic pathogens
Our data showed an accumulation of Enterobacteriaceae in the phyllosphere of arugula with approximately threefold higher abundances, when compared to the rhizosphere. Other plants where high proportions of enterobacteria were found include spinach and bean sprouts sold in grocery stores . Thus, the question arises if a high occurrence of antibiotic resistances could be an additional risk factor for immunocompromised consumers, when coupled with distinct pathogenicity mechanisms. These protective mechanisms naturally occur in microorganism and are even strengthened by various human activities aimed at preventing and treating diseases . An extensive presence of antibiotic resistances was repeatedly shown for various food production systems—mostly due to the preventive use of antibiotics in animal mass production . Our study complements these findings by targeting AR in a leafy green from non-intensive farming, i.e., home gardening. We could show that a broad spectrum of antibiotic resistances was present in the indigenous microbiota and that resistances targeting therapeutic antibiotics unexpectedly prevailed. Arugula is a representative model plant for vegetables that are primarily consumed uncooked in Western cuisine, thus potential pathogens cannot be deactivated by heating. A holistic assessment of the plant resistome showed that the phyllosphere, which constitutes an interface to human health, accounts for most of the identified AR genes. Antibiotic resistance mechanisms were shown to be highly diverse within the microbiota. Similar AR profiles are also found in intensive farming . Remarkably, the most frequent resistances in the arugula microbiome have a significant overlap with the World Health Organization’s list of the most critical antibiotic resistances . While the sole presence of genes related to antibiotic resistances is not evidential for their expression, cultivation-dependent studies provided good means to confirm such findings. In the present study, the high occurrence of resistance-conferring genes in the metagenome was reflected by the tolerance observed in a broad screening of plant-associated isolates. Moreover, members of the Enterobacteriales order were shown to contribute to the resistome profile found in the microbiome. We hypothesize that the bioactive metabolites produced by the host plant facilitate the prevalence of AR in the indigenous arugula microbiome. Co-selections were previously shown for metal and antibiotic resistances . Cross-resistances, where the same genetic determined confers resistance to different substances, play an important role in these mechanisms. In the well-studied enterobacterium Escherichia coli, exposure to salicylate induces multiple antibiotic resistance by activation of the corresponding mar operon . It has been already highlighted that plant-associated bacteria and human opportunists are often equipped with the same genetic basis that facilitates manifestation of disease . Arugula-associated enterobacteria reflect this observation particularly well. The results are based on two plant development stages and include plants from a single sampling site; therefore, they provide a first insight into this highly interesting microenvironment. In order to obtain deeper insights related to the formation, dynamics, and localization (plasmid or genomic) of the antibiotic resistances, additional studies of larger sample sets are required.
The conducted study allowed a holistic assessment of the enterobacterial population in the Eruca sativa Mill. microbiome by combining different cultivation-dependent and cultivation-independent analyses. The occurrence of various Enterobacteriaceae members in the core microbiome indicated that they are well adapted to this particular host plant and part of the natural biota. Moreover, their enrichment in above-ground plant parts is an indicator for their competitiveness in a restrictive environment, which is aggravated by the plant’s bioactive metabolites. Metagenomic analyses in combination with isolates screenings provided evidence for the occurrence of various antibiotic resistance mechanisms in the microbiome. These defense mechanisms were mainly against therapeutic antibiotics. So far, such an AR spectrum was not observed in a natural habitat. Therefore, the arugula system provides a unique model to study AR in an ecological context. The implications of the resistance carriers found in this study for human health as well as the interactions of such microorganisms with the surrounding environment and intruding pathogens remain to be elucidated.
We appreciate the help of Nora Spanberger (Graz) and Olivia Laggner (Graz) during the sample preparations and strain isolations. We also want to thank Melanie Marek (Graz) for her help during susceptibility testing of the bacterial cultures.
This project was funded by the European Funds for Regional Development (EFRE), co-supported by the regional government of Styria (Das Land Steiermark, Austria), project code A3-11.P-33/2011-6.
Availability of data and materials
The metagenomes are publicly accessible on the MG-RAST server (mg-rast.org; ) under the accession numbers mgm4551355.3 (phyllosphere), mgm4551391.3 (rhizosphere), and mgm4551574.3 (bulk soil). All amplicon sequencing datasets were submitted to ENA (ebi.ac.uk/ena) and are accessible under the project accession number PRJEB25651.
The study was designed by AE, MG, and GB. AE prepared the samples for next-generation sequencing and executed the microbiological analyses. TC, AE, and JS conducted the bioinformatics analyses. JS and CS complemented the analyses with their specific expertise. TC, AE, and GB interpreted the overall data. TC and AE visualized the results. TC designed a model of the arugula resistome. TC, AE, MG, and GB wrote the manuscript. All authors read and approved the final version of the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Allen HK, Donato J, Wang HH, Cloud-Hansen KA, Davies J, Handelsman J. Call of the wild: antibiotic resistance genes in natural environments. Nat Rev Microbiol. 2010;8:251.View ArticleGoogle Scholar
- Alqasoumi S, Al-Sohaibani M, Al-Howiriny T, Al-Yahya M, Rafatullah S. Rocket “Eruca sativa”: a salad herb with potential gastric anti-ulcer activity. World J Gastroenterol. 2009;15:1958.View ArticleGoogle Scholar
- Baker-Austin C, Wright MS, Stepanauskas R, McArthur JV. Co-selection of antibiotic and metal resistance. Trends Microbiol. 2006;14:176–82.View ArticleGoogle Scholar
- Bauer AW, Kirby WM, Sherris JC, Turck M. Antibiotic susceptibility testing by a standardized single disk method. Am J Clin Pathol. 1966;45:493.View ArticleGoogle Scholar
- Barillari J, Canistro D, Paolini M, Ferroni F, Pedulli GF, Iori R, Valgimigli L. Direct antioxidant activity of purified glucoerucin, the dietary secondary metabolite contained in rocket (Eruca sativa Mill.) seeds and sprouts. J Agri Food Chem. 2005;53:2475–82.View ArticleGoogle Scholar
- Berg G, Marten P, Ballin G. Stenotrophomonas maltophilia in the rhizosphere of oilseed rape—occurrence, characterization and interaction with phytopathogenic fungi. Microbiol Res. 1996;151:19–27.View ArticleGoogle Scholar
- Berg G, Fritze A, Roskot N, Smalla K. Evaluation of potential biocontrol rhizobacteria from different host plants of Verticillium dahliae Kleb. J Appl Microbiol. 2001;91:963–71.View ArticleGoogle Scholar
- Berg G, Eberl L, Hartmann A. The rhizosphere as a reservoir for opportunistic human pathogenic bacteria. Environ Microbiol. 2005;7:1673–85.View ArticleGoogle Scholar
- Berg G, Smalla K. Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol Ecol. 2009;68:1–13.View ArticleGoogle Scholar
- Binh CT, Heuer H, Gomes NC, Kaupenjohann M, Smalla K. Similar bacterial community structure and high abundance of sulfonamide resistance genes in field-scale manures. In: Dellaguardia CS, editor. Manure: management, uses and environmental impacts. Hauppauge: Nova Science Publishers. 2010; pp. 141–166.Google Scholar
- Blair JM, Webber MA, Baylay AJ, Ogbolu DO, Piddock LJ. Molecular mechanisms of antibiotic resistance. Nat Rev Microbiol. 2015;13:42–51.View ArticleGoogle Scholar
- Brandl MT. Fitness of human enteric pathogens on plants and implications for food safety. Annu Rev Phytopathol. 2006;44:367–92.View ArticleGoogle Scholar
- Brandl MT. Plant lesions promote the rapid multiplication of Escherichia coli O157: H7 on postharvest lettuce. Appl Environ Microbiol. 2008;74:5285–9.View ArticleGoogle Scholar
- Brandl MT, Amundson R. Leaf age as a risk factor in contamination of lettuce with Escherichia coli O157: H7 and Salmonella enterica. Appl Environ Microbiol. 2008;74:2298–306.View ArticleGoogle Scholar
- Bulgarelli D, Rott M, Schlaeppi K, van Themaat EVL, Ahmadinejad N, Assenza F, et al. Revealing structure and assembly cues for Arabidopsis root-inhabiting bacterial microbiota. Nature. 2012;488:91.View ArticleGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Gonzalez Peña A, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R. QIIME allows analysis of high-throughput community sequencing data. Nature Methods. 2010;7(5):335–6.Google Scholar
- Carvalheira A, Silva J, Teixeira P. Lettuce and fruits as a source of multidrug resistant Acinetobacter spp. Food Microbiol. 2017;64:119–25.View ArticleGoogle Scholar
- Cardinale M, Grube M, Erlacher A, Quehenberger J, Berg G. Bacterial networks and co-occurrence relationships in the lettuce root microbiota. Environ Microbiol. 2015;17:239–52.View ArticleGoogle Scholar
- Cohen SP, Levy SB, Foulds J, Rosner JL. Salicylate induction of antibiotic resistance in Escherichia coli: activation of the mar operon and a mar-independent pathway. J Bacteriol. 1993;175:7856–62.View ArticleGoogle Scholar
- Collignon PC, Conly JM, Andremont A, McEwen SA, Aidara-Kane, et al. World Health Organization ranking of antimicrobials according to their importance in human medicine: a critical step for developing risk management strategies to control antimicrobial resistance from food animal production. Clin Infect Dis. 2016;63:1087–93.View ArticleGoogle Scholar
- Cox MP, Peterson DA, Biggs PJ. SolexaQA: At-a-glance quality assessment of Illumina second-generation sequencing data. BMC Bioinformatics. 2010;11 (1).Google Scholar
- Crofts TS, Gasparrini AJ, Dantas G. Next-generation approaches to understand and combat the antibiotic resistome. Nat Rev Microbiol. 2017;15:422–34.View ArticleGoogle Scholar
- Cruz AT, Cazacu AC, Allen CH. Pantoea agglomerans, a plant pathogen causing human disease. J Clin Microbiol. 2007;45:1989–92.View ArticleGoogle Scholar
- Delmotte N, Knief C, Chaffron S, Innerebner G, Roschitzki B, Schlapbach R, von Mering C, Vorholt JA. Community proteogenomics reveals insights into the physiology of phyllosphere bacteria. PNAS. 2009;106:16428–33.View ArticleGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticleGoogle Scholar
- Erlacher A, Cardinale M, Grosch R, Grube M, Berg G. The impact of the pathogen Rhizoctonia solani and its beneficial counterpart Bacillus amyloliquefaciens on the indigenous lettuce microbiome. Front Microbiol. 2014;5:175.View ArticleGoogle Scholar
- Erlacher A, Cardinale M, Grube M, Berg G. Biotic stress shifted structure and abundance of Enterobacteriaceae in the lettuce microbiome. PLoS One. 2015;10:e0118068.View ArticleGoogle Scholar
- Fornefeld E, Schikora A, Berg G, Grosch R, Erlacher A, Kühne T, Smalla K. Humanpathogene Bakterien auf Pflanzen. J Kulturpflanzen. 2015;67:297–309.Google Scholar
- Gomez-Alvarez V, Teal TK, Schmidt TM. Systematic artifacts in metagenomes from complex microbial communities. The ISME Journal 3. 2009;(11):1314–1317.Google Scholar
- Granér G, Persson P, Meijer J, Alström S. A study on microbial diversity in different cultivars of Brassica napus in relation to its wilt pathogen, Verticillium longisporum. FEMS Microbiol Lett. 2003;224:269–76.View ArticleGoogle Scholar
- Halkier BA, Gershenzon J. Biology and biochemistry of glucosinolates. Annu Rev Plant Biol. 2006;57:303–33.View ArticleGoogle Scholar
- Herman KM, Hall AJ, Gould LH. Outbreaks attributed to fresh leafy vegetables, United States, 1973–2012. Epidemiol Infect. 2015;143:3011–21.View ArticleGoogle Scholar
- Heuer H, Krsek M, Baker P, Smalla K, Wellington EM. Analysis of actinomycete communities by specific amplification of genes encoding 16S rRNA and gel-electrophoretic separation in denaturing gradients. Appl Environ Microbiol. 1997;63:3233–41.PubMedPubMed CentralGoogle Scholar
- Heuer H, Schmitt H, Smalla K. Antibiotic resistance gene spread due to manure application on agricultural fields. Curr Opin Microbiol. 2011;14:236–43.View ArticleGoogle Scholar
- Hudzicki J. Kirby-Bauer disk diffusion susceptibility test protocol. Am Soc Microbiol. 2009; http://asmscience.org/content/education/protocol/protocol.3189. Accessed 18 May 2018
- Huson DH, Weber N. Microbial community analysis using MEGAN. Methods Enzymol. 2013;531:465–85.View ArticleGoogle Scholar
- Jirovetz L, Smith D, Buchbauer G. Aroma compound analysis of Eruca sativa (Brassicaceae) SPME headspace leaf samples using GC, GC-MS, and olfactometry. J Agr Food Chem. 2002;50:4643–6.View ArticleGoogle Scholar
- Kurze S, Bahl H, Dahl R, Berg G. Biological control of fungal strawberry diseases by Serratia plymuthica HRO-C48. Plant Dis. 2001;85:529–34.View ArticleGoogle Scholar
- Langmead B, Salzberg S. Fast gapped-read alignment with Bowtie 2. Nature Meth. 2012;9:357–9.View ArticleGoogle Scholar
- Leach JE, Triplett LR, Argueso CT, Trivedi P. Communication in the Phytobiome. Cell. 2017;169:587–96.View ArticleGoogle Scholar
- Leff JW, Fierer N. Bacterial communities associated with the surfaces of fresh fruits and vegetables. PLoS One. 2013;8:e59310.View ArticleGoogle Scholar
- Lüthy B, Matile P. The mustard oil bomb: rectified analysis of the subcellular organisation of the myrosinase system. Biochem Physiol Pflanz. 1984;179:5–12.View ArticleGoogle Scholar
- Martínez-Sánchez A, Marín A, Llorach R, Ferreres F, Gil MI. Controlled atmosphere preserves quality and phytonutrients in wild rocket (Diplotaxis tenuifolia). Postharvest Biol Technol. 2006;40:26–33.View ArticleGoogle Scholar
- Martínez JL. Antibiotics and antibiotic resistance genes in natural environments. Science. 2008;321:365–7.View ArticleGoogle Scholar
- McArthur AG, Waglechner N, Nizam F, Yan A, Azad MA, Baylay AJ, et al. The comprehensive antibiotic resistance database. Antimicrob Agents Chemother. 2013;57:3348–57.View ArticleGoogle Scholar
- Mendes R, Kruijt M, de Bruijn I, Dekkers E, van der Voort M, Schneider JH, Piceno YM, DeSantis TZ, Andersen GL, Bakker PA, Raaijmakers JM. Deciphering the rhizosphere microbiome for disease-suppressive bacteria. Science. 2011;332:1097–100.View ArticleGoogle Scholar
- Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server—a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics. 2008;9:386.View ArticleGoogle Scholar
- Mie A, Andersen HR, Gunnarsson S, Kahl J, Kesse-Guyot E, Rembiałkowska E, Quaglio G, Grandjean P. Human health implications of organic food and organic agriculture: a comprehensive review. Environ Health. 2017;16:111.View ArticleGoogle Scholar
- Nesme J, Cécillon S, TO D, Monier JM, Vogel TM, Simonet P. Large-scale metagenomic-based study of antibiotic resistance in the environment. Curr Biol. 2014;24:1096–100.View ArticleGoogle Scholar
- Nygård K, Lassen J, Vold L, Andersson Y, Fisher I, Löfdahl S, et al. Outbreak of Salmonella Thompson infections linked to imported rucola lettuce. Foodborne Pathog Dis. 2008;5:165–73.View ArticleGoogle Scholar
- Philippot L, Raaijmakers JM, Lemanceau P, van der Putten WH. Going back to the roots: the microbial ecology of the rhizosphere. Nat Rev Microbiol. 2013;11:789–99.View ArticleGoogle Scholar
- Piślewska-Bednarek M, Nakano RT, Hiruma K, Pastorczyk M, Sanchez-Vallet A, Singkaravanit-Ogawa S, Ciesiołka D, Takano Y, Molina A, Schulze-Lefert P, Bednarek P. Glutathione transferase U13 functions in pathogen-triggered glucosinolate metabolism. Plant Physiol. 2017;176:538–51.View ArticleGoogle Scholar
- Rastogi G, Coaker GL, Leveau JH. New insights into the structure and function of phyllosphere microbiota through high-throughput molecular approaches. FEMS Microbiol Lett. 2013;348:1–10.View ArticleGoogle Scholar
- Rumberger A, Marschner P. 2-Phenylethylisothiocyanate concentration and microbial community composition in the rhizosphere of canola. Soil Biol Biochem. 2003;35:445–52.View ArticleGoogle Scholar
- Scardoni G, Petterlini M, Laudanna C. Analyzing biological network parameters with CentiScaPe. Bioinformatics. 2009;25:2857–9.View ArticleGoogle Scholar
- Schuerger AC, Batzer JC. Identification and host range of an Erwinia pathogen causing stem rots on hydroponically grown plants. Plant Dis. 1993;77:472–7.View ArticleGoogle Scholar
- Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.View ArticleGoogle Scholar
- Smalla K, Jechalke S, Top EM. Plasmid detection, characterization and ecology. Microbiology Spectr. 2015; https://doi.org/10.1128/microbiolspec.PLAS-0038-2014.
- Solana M, Boschiero I, Dall’Acqua S, Bertucco A. Extraction of bioactive enriched fractions from Eruca sativa leaves by supercritical CO2 technology using different CO-solvents. J Supercrit Fluids. 2014;94:245–51.View ArticleGoogle Scholar
- Teplitski M, Warriner K, Bartz J, Schneider KR. Untangling metabolic and communication networks: interactions of enterics with phytobacteria and their implications in produce safety. Trends Microbiol. 2011;19:121–7.View ArticleGoogle Scholar
- Verginer M, Siegmund B, Cardinale M, Müller H, Choi Y, Míguez CB, et al. Monitoring the plant epiphyte Methylobacterium extorquens DSM 21961 by real-time PCR and its influence on the strawberry flavor. FEMS Microbiol Ecol. 2010;74:136–45.View ArticleGoogle Scholar
- Vorholt JA. Microbial life in the phyllosphere. Nat Rev Microbiol. 2012;10:828–40.View ArticleGoogle Scholar
- Wang QG, Garrity M, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.View ArticleGoogle Scholar
- Wang FH, Qiao M, Chen Z, Su JQ, Zhu YG. Antibiotic resistance genes in manure-amended soil and vegetables at harvest. J Hazard Mater. 2015;299:215–21.View ArticleGoogle Scholar
- Wassermann B, Rybakova D, Müller C, Berg G. Harnessing the microbiomes of Brassica vegetables for health issues. Sci Rep. 2017;7:17649.View ArticleGoogle Scholar
- Wellington EM, Boxall AB, Cross P, Feil EJ, Gaze WH, Hawkey PM, et al. The role of the natural environment in the emergence of antibiotic resistance in gram-negative bacteria. Lancet Infect Dis. 2013;13:155–65.View ArticleGoogle Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.View ArticleGoogle Scholar
- Zeven AC, de Wet JMJ. Dictionary of cultivated plants and their regions of diversity: excluding most ornamentals, forest trees and lower plants. Wageningen: The Netherlands: CAPD; 1982.Google Scholar