Next Article in Journal
Effects of Below-Ground Microbial Biostimulant Trichoderma harzianum on Diseases, Insect Community, and Plant Performance in Cucurbita pepo L. under Open Field Conditions
Next Article in Special Issue
Investigating the Biocontrol Potential of the Natural Microbiota of the Apple Blossom
Previous Article in Journal
Enhanced Electron Uptake and Methane Production by Corrosive Methanogens during Electromethanogenesis
Previous Article in Special Issue
A Case Study Demonstrates That the Litter of the Rare Species Cinnamomum migao Composed of Different Tissues Can Affect the Chemical Properties and Microbial Community Diversity in Topsoil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparison of the Peel-Associated Epiphytic Bacteria of Anthocyanin-Rich “Sun Black” and Wild-Type Tomatoes under Organic and Conventional Farming

by
Susanna Gorrasi
1,*,
Marcella Pasqualetti
1,2,
Barbara Muñoz-Palazon
1,3,
Giorgia Novello
4,
Andrea Mazzucato
5,
Enio Campiglia
5 and
Massimiliano Fenice
1,6,*
1
Department of Ecological and Biological Sciences (DEB), University of Tuscia, Largo Università snc, 01100 Viterbo, Italy
2
Laboratory of Ecology of Marine Fungi, CoNISMa, Department of Ecological and Biological Sciences, University of Tuscia, Largo Università snc, 01100 Viterbo, Italy
3
Institute of Water Research, University of Granada, 18071 Granada, Spain
4
Department of Science, Technology and Innovation (DISIT), Università del Piemonte Orientale, Viale Teresa Michel, 11, 15121 Alessandria, Italy
5
Department of Agricultural and Forest Sciences (DAFNE), University of Tuscia, Via San Camillo de Lellis snc, 01100 Viterbo, Italy
6
Laboratory of Applied Marine Microbiology, CoNISMa, University of Tuscia, Largo Università snc, 01100 Viterbo, Italy
*
Authors to whom correspondence should be addressed.
Microorganisms 2022, 10(11), 2240; https://doi.org/10.3390/microorganisms10112240
Submission received: 26 October 2022 / Revised: 9 November 2022 / Accepted: 10 November 2022 / Published: 12 November 2022

Abstract

:
Tomatoes are among the most consumed vegetables worldwide and represent a source of health-beneficial substances. Our study represents the first investigating the peel-associated epiphytic bacteria of red and purple (anthocyanin-rich) tomatoes subjected to organic and conventional farming systems. Proteobacteria was the dominant phylum (relative abundances 79–91%) in all experimental conditions. Enterobacteriaceae represented a large fraction (39.3–47.5%) of the communities, with Buttiauxella and Atlantibacter as the most represented genera. The core microbiota was composed of 59 operational taxonomic units (OTUs), including the majority of the most abundant ones. The occurrence of the most abundant OTUs differed among the experimental conditions. OTU 1 (Buttiauxella), OTU 2 (Enterobacteriales), and OTU 6 (Bacillales) were higher in red and purple tomatoes grown under organic farming. OTU 5 (Acinetobacter) had the highest abundance in red tomatoes subjected to organic farming. OTU 3 (Atlantibacter) was among the major OTUs in red tomatoes under both farming conditions. OTU 7 (Clavibacter) and OTU 8 (Enterobacteriaceae) had abundances ≥1% only in red tomatoes grown under conventional farming. PCA and clustering analysis highlighted a high similarity between the bacterial communities of red and purple tomatoes grown under organic farming. Furthermore, the bacterial communities of purple tomatoes grown under organic farming showed the lowest diversity and evenness. This work paves the way to understand the role of nutritional superior tomato genotypes, combined with organic farming, to modulate the presence of beneficial/harmful bacteria and supply healthier foods within a sustainable agriculture.

1. Introduction

The tomato (Solanum lycopersicum L.) represents an economically important vegetable crop, with a world production of over 186 million tons [1]. It is one of the most consumed vegetables in the world and a common food in the Mediterranean diet, and it is produced and consumed both as fresh and processed products (e.g., paste, peeled tomatoes, diced products, and sauces) [2].
Tomatoes are a dietary source of health-beneficial substances, such as fibers, vitamins, minerals, and antioxidants (e.g., phenolic compounds, β-carotene, and in particular lycopene) [3,4]. Due to its nutritional importance and to an increasing consumers’ interest in food production and quality, high attention has been paid both to sustainable agricultural practices and to breeding strategies aiming at improving the shelf-life, the quality traits, and the nutraceutical value of tomato (e.g., high content in phytochemicals and bio-active compounds) [5,6].
In recent decades, there has been a worldwide awareness increase about the need of sustainable agricultural practices, especially focused to increase yield while also addressing environmental conservation and human health. Sustainable cropping systems are intended to preserve soil health and improve its chemical-physical characteristics, reduce synthetic fertilizer and pesticide inputs, and increase product quality [7].
The techniques used to increase tomato yield, paying attention to environmental sustainability and human health, include the hairy vetch mulch system, the use of plant-beneficial microorganisms to improve the fruit quality and reduce fertilization, the development of new cultivars with enhanced nutritional qualities and tolerance/resistance to abiotic and biotic stresses [8,9,10].
The development of nutritionally superior genotypes, in particular the selection of cultivars having high contents of antioxidants (such as phenolic compounds), is becoming a significant research trend. Among them, anthocyanins are secondary metabolites of the flavonoid family, playing a significant role in plant tolerance to biotic/abiotic stress and in the prevention and treatment of various important human diseases [11,12].
Several studies based on both conventional breeding and genetic engineering strategies have allowed obtaining anthocyanin-rich tomato lines [11]. These dark-skinned colored tomatoes (so-called “purple” or “black” tomatoes) result from mutations in carotenoid/flavonoid biosynthetic pathways [13]. Compared to the wild type, tomatoes enriched in anthocyanins have shown higher antioxidant activity and increased resistance to some pathogens, resulting in a significant extension of the shelf life [14].
A 20-year breeding activity carried out at the University of Tuscia (Viterbo, Italy) allowed to select and fix the anthocyanin-rich tomato line called “Sun Black” (SB). It was obtained by combining the mutations of Anthocyanin fruit (Aft) and atroviolaceum (atv) alleles, affecting the anthocyanin biosynthesis, in the same line [15]. The SB tomato features an accumulation of anthocyanins in the peel and other bioactive compounds (e.g., gentistic acid, rutin, and vitamin C) in the whole fruit [13].
The studies carried out on purple tomatoes have been mainly focused on the characterization of their nutritional/bioactive compounds and health beneficial effects and, to a lesser extent, on their longer shelf life and resistance to some pathogens [13,14,15,16].
Phyllospheric microbiota refers to the microorganisms found on the above-ground part of the plants (e.g., stems, leaves, flowers, and fruits), which may be epiphytic and endophytic, and can be both beneficial or harmful to the host plant [17,18]. Actually, investigating the epiphytic microbial communities of fruits and vegetables is important to improve food quality and safety, being the consumers directly exposed to these microorganisms. This is particularly true for widely consumed fresh products. In fact, fresh fruits and vegetables have been recognized as reservoirs of foodborne pathogens and spoilage microorganisms [19,20]. On the other hand, the microbiota naturally present on these foodstuffs harbors also beneficial representatives, which are involved in plant fitness and produce a number of industrially important compounds (e.g., antibiotics, enzymes, organic acids, and vitamins) [21,22,23]. Among them, the bioactive compound producers (having antagonistic effects on plant pathogens) can be employed as biological control agents [18,24,25]. Although the plant episphere, besides being a reservoir of possible harmful microorganisms, has been proven to be an important source of beneficial microorganisms (and their bioactive compounds), epiphytic microbial communities are the least investigated among the plant microbiota.
To the best of our knowledge, in general, no study has compared the microbiota of red and purple tomatoes, particularly in relation to peel-associated epiphytic bacteria. Furthermore, the bacterial diversity harbored by two genotypes subjected to different farming practices has never been studied. Since microorganisms are known to be involved in ecosystem functioning, plant physiology and response to abiotic/biotic stresses, as well as in phytopathology and infectious disease outbreaks, a focus on the two tomato phyllospheric microbiota would be of scientific and economic interest.
This work aimed to characterize the peel-associated epiphytic bacteria of SB and the near-isogenic wild type (WT) tomatoes grown under two farming systems (organic and conventional), shedding light on possible assemblage differences in their bacterial communities. This would provide useful information on the importance of combining nutritional superior tomato genotypes with the right farming management, to provide healthier foods in a sustainable agriculture framework.

2. Materials and Methods

2.1. Plant Materials and Growth Conditions

The research was carried out on the two tomato lines (SB and WT) grown in open field at the Experimental Farm “Nello Lupori” (42°25′35.71818″ N; 12°4′49.34524″ E) of the University of Tuscia (Viterbo, Italy) under organic and conventional farming (referred here as OF and CF, respectively). The two agrotechnical protocols, differing in the precession crop and soil management, were applied as previously reported [26,27]. The experimental conditions (genotype × farming condition) considered in this study were: T1, WT/OF; T2, WT/CF; T3, SB/OF; T4, SB/CF.
The two lines were previously selected by Mazzucato et al. [15] and, except for the different fruit colors (SB and WT are purple- and red-colored, respectively), they are characterized by the same phenotypic traits.
For both lines, seeds were germinated in Petri dishes, then plantlets were transplanted in twin rows (planting distances: 160 cm between twins, 40 cm between rows, and 50 cm between plants within the row) at the 4–5th true leaf stage. Ten plants per line were arranged in twin rows using a randomized block design with three replicates. Plants were grown on tutors with the standard agronomic practices for fresh market tomatoes. Lateral shoots were weekly removed, and plants were left to open pollination.
The soil in the experimental area was a Typic Xerofluvent (volcanic) with the following characteristics in the 0–30 cm layer: 668 g kg−1 of dry soil sand, 194 g kg−1 of dry soil silt, 138 g kg−1 of dry soil clay, bulk density 1.40 g cm−3, pH 6.43 (water, 1:2.5), organic matter 17.5 g kg−1 of dry soil, total nitrogen 1.01 g kg−1 of dry soil, and available phosphate 13.92 g kg−1 of dry soil.

2.2. Sample Collection and Processing

For each experimental condition (genotype × farming), a total of 16 healthy and representative fruits (composite sample) were aseptically sampled (they were collected wearing gloves and put into sterile plastic bags) from the three crop replicates. To obtain the composite sample, four fruits were collected from four plants located in different crop replicates. The fruits were taken from various positions on the plant, avoiding those showing the skin covered by soil debris [28]. The samples were immediately transported to the nearby laboratory to be processed.
The harvested fruits were softly pre-washed to remove debris and microorganisms that may have been incidentally deposited on the fruit surface, as reported by Janisiewicz et al. [29] with some adjustments: each sample was transferred to a sterile plastic bag with 1000 mL of sterile deionized water and gently hand massaged for 1 min. The pre-washing water was discarded and the recovery of the skin-associated microorganisms was performed as reported by Fernández-Suárez et al. [30] with slight modifications. The washing solution (1 L of sterile water containing 0.1% Tween 20) was added to the bag that was closed and placed into another bag. The double-bagged fruits were gently shaken for 45 min. The above washings were carried out 25 °C. The washing solutions were vacuum-filtered on sterile 0.22 µm nitrocellulose membranes (Millipore, Burlington, MA, USA), which were kept frozen until DNA extraction.

2.3. DNA Extraction, 16S rRNA Gene Amplicon Libraries and Sequencing

Total DNA was extracted from filters using the ZR Fungal/Bacterial DNA MiniPrep kit (Zymo Research Corp., Irvine, CA, USA) following the manufacturer’s instructions. A two-step PCR was used to prepare the multiplexed amplicon libraries. To profile the peel-associated bacterial diversity, the V5–V6 hypervariable regions of the 16S rRNA gene were amplified using the primer set 783F/1046R [31,32]. The strategy and the protocols used to prepare the libraries have already been described in detail by Gorrasi et al. [33]. The amplicon libraries were sequenced at Nuova Genetica Italiana by Illumina MiSeq (Illumina, San Diego, CA, USA) using a 2 × 250 bp paired-end protocol.
The sequence data have been submitted to the Sequence Read Archive (SRA) database under the accession number PRJNA798869.

2.4. Sequence Processing and Data Analysis

Sequence processing and data analyses were performed as previously reported [34,35].
Raw reads were demultiplexed according to the indices and internal barcodes and processed using the mothur software version 1.44.3. Forward and reverse reads were merged into contigs with a quality score of mismatched base = 0. The reads were subjected to quality trimming to remove sequences with more than eight homopolymers and ambiguous bases. The suspected chimeras were identified using the VSEARCH algorithm and then removed [35]. Unique sequences were determined and then mapped against the original dataset to calculate the abundance data. The high-quality reads were clustered into operational taxonomic units (OTUs) using a 97% similarity cut-off. Singletons (sequences appearing only once in the entire dataset) were removed.
Taxonomic profiling was assessed using the RDP classifier [36], applying a 50% confidence cut-off for the assignation [37].

2.5. Statistical Methods

Alpha diversity was described by the Shannon diversity and Pielou evenness indices. The Shannon index ranges from 0 to +∞ and increases as the community richness and evenness increase [38]. Pielou evenness index ranges from 0 to 1, and decreasing values indicate a decrease in evenness and a prevalence of few dominant species [39]. Alpha diversity indices were calculated using Species Diversity & Richness IV (SDR-IV) software (PISCES Conservation Ltd., Lymington, UK), and index comparison among the experimental conditions was assessed with the Solow test (based on 10,000-iterated randomization test) [40,41]. Indices were calculated on samples rarefied to 5299 randomly chosen sequences, corresponding to the minimum number of sequences at a sample [33].
Beta diversity was analyzed by principal component analysis (PCA) and clustering analysis, based on square root transformed data (to avoid rare OTU overweighting). The PCA analysis was performed using the CANOCO v. 5.1 software package (Microcomputer Power, Ithaca, NY, USA). The clustering analysis (complete linkage) was performed based on Bray–Curtis dissimilarity data, applying the similarity profile (SIMPROF) test [42] to individuate distinct bacterial assemblages among the experimental conditions (statistically different at p < 0.05; 1000 permutations and 999 simulations). The clustering analysis was run using the software package PRIMER-E v.6.1.18 (Plymouth, UK).
The shared OTUs among the experimental conditions were visualized by a 4-way Venn diagram, generated using the InteractiVenn tool (http://www.interactivenn.net/index2.html, accessed on 1 October 2022) [43].

3. Results

3.1. Composition of the Epiphytic Peel-Associated Bacteria

Amplicon sequencing generated a total of 69,948 high-quality bacterial sequences, resulting in 10,287 OTUs defined by the clustering process.
Applying a 50% cut-off [37], OTUs were assigned to 10 phyla, 21 classes, 46 orders, 84 families, and 161 genera. Unassigned OTUs ranged in-between 0.6–5.1%, 0.8–6.4%, 4.3–14.7%, 14.8–38.7%, and 44.0–58.8% at phylum, class, order, family, and genus level, respectively.
Overall, the major taxa (having relative abundance (Ra) ≥ 1% in at least one sample) in the peel-associated community were Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes (Figure 1). Among them, Proteobacteria was the dominant phylum (79–91%) in all experimental conditions, with the highest abundance in T4 (SB/CF). Firmicutes (Ras 2.1–9.7%) was found with the highest abundance in T3 (SB/OF). Both Actinobacteria and Bacteroidetes showed the highest abundances (11.2% and 6% for Actinobacteria and Bacteroidetes, respectively) in T2 (WT/CF), but they represented minor taxa in some experimental conditions. Actinobacteria showed Ras < 1% in SB tomatoes grown under organic farming (T3), whereas Bacteroidetes showed Ras < 1% in SB tomatoes grown both under organic and conventional farming (T3 and T4).
Overall, the studied peel-associated bacterial flora mostly belonged to the family Enterobacteriaceae (39.3–47.5%), being Buttiauxella and Atlantibacter the most represented genera (with 11,005 and 7156 sequences over the whole dataset, respectively).
However, Enterobacteriaceae, Pseudomonadaceae, Moraxellaceae, Microbacteriaceae, Erwiniaceae, Yersiniaceae, Rhodobacteraceae, Alcaligenaceae, Cytophagaceae, Weeksellaceae, Streptococcaceae, and Oceanospirillaceae were recognized as major (Ra ≥ 1% in at least one sample) families (Figure 2).
Among the 161 genera found across all samples, Buttiauxella, Atlantibacter, Acinetobacter, Siccibacter, Clavibacter, Pseudomonas, Phytobacter, Pantoea, Chania, Cedecea, Azomonas, and Lactococcus were recorded as major taxa (Figure 3).
Buttiauxella was most abundant in T1, T3, and T4 (Ras in-between 13.2–27%), with the highest abundances on SB tomato peel; it was least abundant (2.4%) in T2 (WT/CF).
Atlantibacter was mainly found in WT tomatoes, having the highest abundance in those grown under CF (23.3%); in SB grown under CF, it was recorded with the lowest abundance (0.7%), representing a rare taxon.
Acinetobacter (Ras 1.8–7.1%) showed the lowest abundances in T1 and T4 (WT/OF and SB/CF, respectively).
Siccibacter, Clavibacter, Pseudomonas, and Phytobacter, overall recorded with abundances lower than 2.5% and representing sometimes rare taxa, showed a peak in one experimental condition: Siccibacter in T3 (9.7%), Clavibacter in T2 (10.7%), Pseudomonas and Phytobacter in T4 (9.3% and 6.4%, respectively).
Pantoea showed higher abundances in tomatoes grown under OF (2.5% and 4.2%, in T1 and T3, respectively).
Chania, Cedecea, Azomonas, and Lactococcus were recognized as major genera only in one experimental condition: Chania in T4 (3.9%), Cedecea in T2 (2.3%), Azomonas in T4 (1.3%), and Lactococcus in T3 (1.1%).
Among the 10,287 defined OTUs, the most abundant ones (having a number of mapping sequences > 1000 on the whole dataset) were OTU 1 (Buttiauxella), OTU 2 (order Enterobacteriales), OTU 3 (Atlantibacter), OTU 5 (Acinetobacter), OTU 6 (order Bacillales), OTU 7 (Clavibacter), and OTU 8 (family Enterobacteriaceae). Overall, OTU 1, OTU 2, and OTU 3 were the most represented (Figure 4).
OTU 1 was the sole showing Ra ≥ 1% in all experimental conditions, and it was recorded with the highest abundances (7.9% and 19.6% in T1 and T3, respectively) in both tomato genotypes grown under OF.
OTU 2 showed Ra < 1% in both tomato genotypes grown under CF, and it showed a peak (14.4%) in T3 (SB/OF).
Different from OTU 1 and OTU 2, showing different abundances according to the farming condition, OTU 3 had higher abundances in WT than SB tomatoes. It was detected with Ra < 1% in SB tomatoes, and it showed a peak (12.9%) in T2 (WT/CF).
Also OTU 6 revealed different abundances according to the farming condition, although with less marked differences than those observed for OTU 1 and OTU 2. It was detected with higher Ras in tomatoes grown under OF (1.1% and 5.6% in T1 and T3, respectively), whereas it was a rare OTU in those grown under CF.
OTU 5 did not reveal abundance differences according to the tomato genotype or the farming condition; it was a rare OTU in T1 and T4 and had Ra ≥ 1% in T2 (2.2%) and T3 (4.7%).
Finally, OTU 7 and OTU 8 had Ra ≥ 1% only in WT tomatoes grown under CF (OTU 7:5.5%; OTU 8:4.1%).

3.2. Unique and Shared Bacterial OTUs

The comparison of the peel-associated bacterial communities revealed a very low number of ubiquitous OTUs (common to all experimental conditions) and a high number of specific members (OTUs exclusively found in one experimental condition).
The four experimental conditions shared 59 out of 10,287 OTUs, representing the core microbiota (Figure 5). Except for OTU 7, which was shared only by three experimental conditions (T1, T2, and T3), the core microbiota included all the other most abundant OTUs.
The core microbiota represented 22.6%, 24.5%, 46.9%, and 5.1% of the peel-associated bacterial communities in T1, T2, T3, and T4, respectively. It is worth noting that it accounted for a similar proportion in WT tomatoes grown under OF and CF, whereas it represented a very different moiety in communities found in SB, depending on the farming condition. Only one OTU belonged to Firmicutes (order Bacillales), whereas all the remaining OTUs belonged to Proteobacteria (Figure 6). Furthermore, only representatives of the following eight genera were found in the core microbiota: Buttiauxella, Atlantibacter, Acinetobacter, Siccibacter, Pseudomonas, Phytobacter, Pantoea, and Cedecea.
Differently from the core microbiota, including a low number of OTUs, a high number of exclusive OTUs were found in each experimental condition. The number of exclusive OTUs was 636, 3503, 2394, and 1244 for T1, T2, T3, and T4, respectively (Figure 5 Diagramma Venn OTU).
Among the 2082 OTUs found in T1, 31% were exclusive to this experimental condition. These exclusive OTUs (having Ras in-between 0.04–0.26%) accounted for 30.6% of the peel epiphytic bacterial community in WT tomatoes grown under OF. These OTUs belonged to the phyla Proteobacteria, Actinobacteria, Bacteroidetes, Firmicutes, Cyanobacteria, Verrucomicrobia, and Rhodothermaeota (Table S1). Among these OTUs, 126 were classified at genus level; they belonged to 42 genera, with Atlantibacter and Metakosakonia as the most represented genera (including 26 and 13 OTUs, respectively) (Table S1).
In T2, 66% of the recorded OTUs were exclusive. These OTUs (with Ras in-between 0.008–0.31%) collectively accounted for 45.7% of the bacteria found on WT tomatoes grown under CF. Moreover, they were representatives of the phyla Proteobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Cyanobacteria, Tenericutes, Verrucomicrobia, Campilobacterota, and Deinococcus-Thermus (Table S2). The OTUs classified at the genus level were representatives of 112 genera; in particular, Atlantibacter and Clavibacter (411 and 291 OTUs, respectively) were the most retrieved taxa (Table S2).
Fifty-seven percent of the OTUs found in T3 were exclusive. These OTUs, showing Ras in the range 0.006–0.14%, accounted for 27.1% of the bacterial communities of SB tomatoes grown under OF. At the phylum level, these OTUs belonged to Proteobacteria, Firmicutes, Actinobacteria, Rhodothermaeota, Bacteroidetes, and Cyanobacteria (Table S3). In addition, the OTUs exclusively found in T3 were representatives of 22 genera; among them, Siccibacter, Buttiauxella, and Pantoea (344, 199, and 102 OTUs, respectively) were the most found (Table S3).
In T4, 74% of the recovered OTUs were exclusive. These OTUs (with Ras in-between 0.03–1.66%; among them, only three OTUs had Ra ≥ 1%) accounted for 71.6% of the peel epiphytic bacteria of SB tomatoes grown under CF. At the phylum level, these OTUs were ascribed to Proteobacteria, Actinobacteria, Firmicutes, Bacteroidetes, Cyanobacteria, Tenericutes, and Verrucomicrobia (Table S4). In addition, the OTUs classified at the genus level were ascribed to 63 genera; Pseudomonas, Buttiauxella, Phytobacter, Chania, and Atlantibacter (80, 71, 65, 55, and 51 OTUs, respectively) were the most abundant (Table S4).

3.3. Alpha- and Beta-Diversity Analyses

Alpha diversity analysis was assessed by the Shannon diversity and Pielou evenness indices. The Shannon index was 6.734, 6.821, 5.515, and 6.673 for T1, T2, T3, and T4, respectively. It significantly differed between SB tomatoes grown under OF and CF, whereas no differences were observed between WT tomatoes in relation to the farming condition (based on Solow test). Considering the farming conditions, Shannon index values significantly differed between WT and SB grown under OF, but not between WT and SB grown under CF (Solow test). The Pielou index was 0.881, 0.795, 0.661, and 0.899 for T1, T2, T3, and T4, respectively. According to the Solow test, significant differences in community evenness were revealed among all experimental conditions except for T1 (WT/OF) and T4 (SB/CF). Overall, alpha diversity indicated that the T3 (SB/OF) peel-associated bacterial communities showed the lowest diversity and evenness (lowest Shannon and Pielou index values), and therefore they were characterized by the presence of some dominant OTUs.
PCA and clustering analyses were performed to evaluate the effect of the tested conditions (genotype × farming condition) on the epiphytic bacterial communities found on the tomato peel.
Four plot visualizations were provided for the same PCA analysis to better show the OTU occurrence variation according to the four experimental conditions (Figure 7).
Overall, the PCA (explaining 75.87% of the total variance) showed that T1 and T3 grouped together, indicating that the bacterial communities of both tomato genotypes grown under OF were similar. Conversely, T2 and T4 were separated from each other and from the other conditions.
As an overview, it is possible to observe a cloud of OTUs at the axis intersection, indicating a cluster of OTUs showing no notable differences (in distribution and abundance) among the various experimental conditions. By contrast, some OTUs, well separated from the cluster, appear to be much more oriented toward a certain condition. In particular, OTU 1, OTU 2, OTU 5, and OTU 6 were close to T1 and T3 (Figure 7a,c). As reported above, OTU 1, OTU 2, and OTU 6 showed different abundances according to the farming condition, having higher abundances in both tomato genotypes grown under OF; OTU 5 was detected with the highest abundance in T3. OTU 3, OTU 7, and OTU 8 were closer to T2, showing their highest abundances in WT tomatoes grown under CF (Figure 7b). OTU 13, OTU 14, OTU 19, OTU 20, OTU 24, OTU 25, and OTU 26 were ordered along the T4 vector (Figure 7d); these OTUs were detected with Ras in-between 1–3% only in SB grown under CF, whereas they were rare OTUs in all other conditions.
The clustering analysis (Figure 8) confirmed that the bacterial communities found in T1 and T3 did not show significant differences (p < 0.05, SIMPROF test). In addition, the hierarchical clustering dendrogram showed that T4 communities were the most dissimilar.

4. Discussion

The selection of tomato cultivars with a high content of antioxidants is one of the research trends in the production of nutritionally superior genotypes. In this context, no investigation has been carried out to compare the microbiota of red and purple tomatoes, and to get a more extended overview of the possible bacterial assemblage differences of the two genotypes subjected to different farming practices.
The current study sought to provide the first information on this topic, comparing the peel-associated epiphytic bacteria of the WT (red) and SB (purple) tomatoes grown under organic and conventional farming systems.
The taxonomic profiling showed that Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes were the major taxa across all experimental conditions, with Proteobacteria as the dominant phylum. These taxa are ubiquitous in the soil and are widely found on the surface of fruits and vegetables, including tomatoes [44,45,46].
Enterobacteriaceae, Pseudomonadaceae, Moraxellaceae, Microbacteriaceae, Erwiniaceae, Yersiniaceae, Rhodobacteraceae, Alcaligenaceae, Cytophagaceae, Weeksellaceae, Streptococcaceae, and Oceanospirillaceae were detected as major families. Most of these families or some of their genera are commonly associated with edible plants and are recorded in both epiphytic and endophytic bacterial communities [44,47,48].
It is worth noting that in this work most of the peel-associated bacteria found across all the experimental conditions belonged to the family Enterobacteriaceae. These results are consistent with those reported in the literature. Enterobacteriaceae is among the most represented families in both epiphytic and endophytic plant microbiota [44,49]. Many members are natural inhabitants of several food plants, and a high abundance of this taxon on the surface of various fruits and vegetables (including tomatoes) has been reported [28,44,46,50]. In particular, a high Enterobacteriaceae prevalence is a feature of the bacterial communities of fruits and vegetables grown close to the soil (e.g., strawberries, lettuce, peppers, and tomatoes) [28,44].
Among the 161 annotated genera, Buttiauxella and Atlantibacter were the most abundant (with 11,005 and 7156 sequences over the whole dataset, respectively), and they showed an opposite pattern of prevalence according to the tomato genotype. Buttiauxella had the highest abundances in SB tomatoes, in particular in those grown under CF (Ra = 27%); Atlantibacter was mainly found in WT tomatoes and showed the highest abundance (23.3%) in those grown under CF.
Buttiauxella species have been isolated especially from the intestines of snails and slugs, but have been also retrieved in water and soil [51,52]. Various Buttiauxella members have been found in the plant rhizosphere, where they can play important roles as phosphate solubilizers and plant growth promoters [53,54,55]. Some members have been found in plant phyllosphere [56], including that of tomato [57,58].
Atlantibacter has been proposed as a novel genus to include the two species A. hermannii and A. subterranea (formerly Escherichia hermannii and Salmonella subterranea) [59]. A. hermannii members have clinical significance, being mostly found in human wounds, sputum, and stool samples [60], whereas A. subterranea strains were isolated from a river and heavy metal contaminated sediments [61,62]. Scarce information is available on Atlantibacter members found in association with plants and their possible roles. However, some of them were revealed among plant microbiotas (e.g., in tobacco leaves, maize aerial root mucilage, and wheat bran) [63,64,65], and a representative of this genus (A. hermannii DDE1) featuring potential plant growth-promoting ability was found in pumpkin roots [66].
The other major genera that we detected on the tomato peels were Acinetobacter, Siccibacter, Clavibacter, Pseudomonas, Phytobacter, Pantoea, Chania, Cedecea, Azomonas, and Lactococcus. These taxa have members that are common plant and/or soil residents [67,68,69,70,71,72,73,74,75], mostly featuring beneficial activities toward the plants [74,75,76,77]. However, some of these genera include also detrimental bacteria [78,79].
Acinetobacter, Pseudomonas, and Pantoea are commonly found in high abundance in the rhizosphere and phyllosphere (including carposphere) of various plants, including tomato [28,46,55,58,65,80]. Various strains of these taxa have shown beneficial interactions with various plants (including tomato), owing to their abilities as plant growth promoters, bio-fertilizers, and biocontrol agents [22,74,77,81]. As for Pseudomonas and Pantoea, although encompassing a great number of beneficial members (even those having antagonistic activity toward tomato phytopathogens [81,82,83,84]), they include also plant pathogens. Notable examples are Pseudomonas syringae and Pantoea ananatis, which cause diseases in various economically important crops, including the bacterial speck and graywall of tomato [78,85].
The current study evidenced that Pantoea and Pseudomonas showed a different distribution according to the farming condition, with an opposite pattern of prevalence: Pantoea had the highest abundances (2.5–4.2%) in tomatoes grown under OF, whereas Pseudomonas was a minor taxon (Ras in-between 0.02–0.7%) in OF tomatoes and a major genus in CF tomatoes, showing a peak of abundance (9.3%) in SB/CF. However, it is not possible to discriminate possible beneficial/detrimental representatives among the OTUs ascribed to these genera, since the length of the sequenced fragment does not allow to gain deeper insight. Differently, Acinetobacter occurrence on tomato peel seemed not to be related to the tomato genotype or farming system, being present in all the experimental conditions as a major taxon and showing the highest abundances in T2 (WT/CF) and T3 (SB/OF).
Another taxon deserving discussion is Clavibacter. It is generally recognized as a genus of phytopathogens of agricultural significance, including those involved in the bacterial canker of tomato [86]. Nevertheless, the occurrence of tomato non-pathogenic Clavibacter strains has been reported [71].
Our results showed that Clavibacter had a different distribution according to the tomato genotype, being a rare taxon (Ra = 0.1%) in SB and a major genus in WT tomatoes, with an abundance peak (10.7%) in those grown under CF. However, since the tomatoes collected in this study did not show bacterial canker lesions, it is possible to assume that the detected peel-associated Clavibacter OTUs were non-pathogenic representatives.
In the current investigation, 10,287 OTUs were recognized across the whole dataset. Our results evidenced that a high number of exclusive OTUs were found in each experimental condition, but they showed a very low abundance. Differently, the tomato peel-associated core microbiota consisted of a much lower OTU number (59 OTUs, representing ~0.6% of the total number of OTUs). It is worth noting that the core microbiota accounted for a similar proportion in WT tomatoes grown under OF and CF (22.6% and 24.5%, respectively), whereas it represented a very different moiety in SB communities, depending on the farming system. Among the communities related to the various experimental conditions, the core microbiota represented the highest proportion in those retrieved on SB grown under OF (46.9%) and the lowest fraction (accounting for 5.1% of the community only) in those of SB grown under CF.
Overall, the core microbiota OTUs were representative of eight genera only: Buttiauxella, Atlantibacter, Acinetobacter, Siccibacter, Pseudomonas, Phytobacter, Pantoea, and Cedecea. As discussed above, these taxa include various beneficial members.
Looking deeper at the shared OTUs, it is possible to note that among them there were almost all the most abundant OTUs (the exception was OTU 7, being absent in T4), suggesting a possible key role in the structuring and preservation of tomato peel communities.
However, the most abundant OTUs occurrence was different among the investigated conditions. The PCA analysis evidenced that OTU 1 (Buttiauxella), OTU 2 (order Enterobacteriales), OTU 5 (Acinetobacter), and OTU 6 (order Bacillales) were related to T1 and T3 (WT/OF and SB/OF, respectively) (Figure 7a,c). Actually, OTU 1, OTU 2, and OTU 6 showed different abundances according to the farming condition, having higher abundances in tomatoes (both genotypes) grown under OF (Figure 4). Instead, OTU 5 was detected with the highest abundance in WT/OF condition (T3) (Figure 4). OTU 3 (Atlantibacter), OTU 7 (Clavibacter), and OTU 8 (family Enterobacteriaceae) were related to T2 (Figure 7b), showing their highest abundances in WT tomatoes grown under CF (Figure 4). In particular, OTU 3 was a rare OTU in SB tomatoes and among the major OTUs in WT tomatoes, showing an abundance peak (12.9%) in WT/CF condition. OTU 7 and OTU 8 represented rare OTUs in almost all experimental conditions, except in WT tomatoes grown under CF.
Alpha diversity analysis indicated that the farming system affected the peel-associated bacterial community diversity of SB tomatoes, being significantly lower in SB grown under OF than those grown under CF. Conversely, no significant differences in the community diversity were observed in WT tomatoes according to the farming condition. Considering the same farming condition, significant differences in bacterial diversity between the two different genotypes were observed only for the organic crop management. The bacterial communities of WT/OF tomatoes showed higher diversity than those of SB/OF. Both genotype and farming conditions seemed to affect bacterial evenness. Significant differences were revealed in the community evenness of both WT and SB tomatoes according to the farming condition. The bacterial communities of WT grown under OF showed a higher evenness than those grown under CF. On the contrary, the bacterial communities of SB grown under OF showed a lower evenness than those grown under CF. Considering the same farming condition, significant differences in bacterial evenness between the two different genotypes were observed for both crop managements. In OF, the bacterial communities of SB tomatoes showed a lower evenness than those of WT. By contrast, in CF, the bacterial communities of SB tomatoes were characterized by a higher evenness than those of WT.
Overall, the PCA analysis showed that T1 and T3 grouped together, indicating that the bacterial communities of both tomato genotypes grown under OF were similar. On the contrary, T2 and T4 were well separated, indicating that the bacterial communities found on WT and SB tomatoes grown under CF differed from each other and from those of the other experimental conditions. This was confirmed by the clustering analysis (Figure 8), revealing that the bacterial communities of T1 and T3 did not show significant differences (p < 0.05, SIMPROF test).

5. Conclusions

To the best of our knowledge, this work was the first investigation regarding the peel-associated epiphytic bacteria of red and purple (anthocyanin-rich) tomatoes subjected to different farming systems. This study provided the first insights into similarities/differences among the bacterial assemblages revealed in both genotypes grown under organic and conventional farming. Further investigations would be useful to give a more extended view of the tomatoes-bacteria associations, to provide a comparative analysis of the bacterial assemblages in other compartments of the tomato plant (in particular in pulp, surveying also the possible presence of bacteria with implications for human health). In addition, it would be interesting to investigate if new bacterial strains with potential beneficial effects for the plant (to be also employed in sustainable agriculture) are present and/or enriched in the purple tomatoes. Moreover, future developments of the study could involve additional factors, such as different soils and/or locations for the farming to address possible effects of biogeography on the microbiota.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/microorganisms10112240/s1, Table S1: OTUs exclusively found in T1 (WT/OF); Table S2: OTUs exclusively found in T2 (WT/CF); Table S3: OTUs exclusively found in T3 (SB/OF); Table S4: OTUs exclusively found in T4 (SB/CF).

Author Contributions

Conceptualization, S.G. and M.F.; methodology, S.G., M.P., B.M.-P. and M.F.; software, S.G., M.P. and B.M.-P.; validation, S.G., M.P., B.M.-P., G.N., A.M., E.C. and M.F.; formal analysis, S.G., M.P., B.M.-P., G.N. and M.F.; investigation, S.G., M.P., B.M.-P. and M.F.; resources, S.G., M.P., B.M.-P., A.M., E.C. and M.F.; data curation, S.G. and M.F.; writing—original draft preparation, S.G., M.P., B.M.-P., G.N., A.M., E.C. and M.F.; writing—review and editing, S.G., M.P., B.M.-P., G.N., A.M., E.C. and M.F.; visualization, S.G., M.P., B.M.-P., G.N., A.M., E.C. and M.F.; supervision, S.G., M.P., B.M.-P., G.N., A.M., E.C. and M.F.; project administration, A.M., E.C. and M.F.; funding acquisition, A.M., E.C. and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by Filas project MIGLIORA from Regione Lazio, Italy.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequence data are available at the Sequence Read Archive (SRA; https://www.ncbi.nlm.nih.gov/sra, accessed on 1 October 2022) repository; study accession number: PRJNA798869.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. FAOSTAT, Crops and Livestock Products, Tomatoes. Available online: http://www.fao.org/faostat/en/#data/QCL/visualize (accessed on 1 October 2022).
  2. Kandel, D.R.; Marconi, T.G.; Badillo-Vargas, I.E.; Enciso, J.; Zapata, S.D.; Lazcano, C.A.; Crosby, K.; Avila, C.A. Yield and fruit quality of high-tunnel tomato cultivars produced during the off-season in South Texas. Sci. Hortic. 2020, 272, 109582. [Google Scholar] [CrossRef]
  3. Borghesi, E.; González-Miret, M.L.; Escudero-Gilete, M.L.; Malorgio, F.; Heredia, F.J.; Meléndez-Martínez, A.J. Effects of salinity stress on carotenoids, anthocyanins, and color of diverse tomato genotypes. J. Agric. Food Chem. 2011, 59, 11676–11682. [Google Scholar] [CrossRef] [PubMed]
  4. Caseiro, M.; Ascenso, A.; Costa, A.; Creagh-Flynn, J.; Johnson, M.; Simões, S. Lycopene in human health. LWT 2020, 127, 109323. [Google Scholar] [CrossRef]
  5. Ronga, D.; Gallingani, T.; Zaccardelli, M.; Perrone, D.; Francia, E.; Milc, J.; Pecchioni, N. Carbon footprint and energetic analysis of tomato production in the organic vs. the conventional cropping systems in Southern Italy. J. Clean. Prod. 2019, 220, 836–845. [Google Scholar] [CrossRef] [Green Version]
  6. Natalini, A.; Acciarri, N.; Cardi, T. Breeding for nutritional and organoleptic quality in vegetable crops: The case of tomato and cauliflower. Agriculture 2021, 11, 606. [Google Scholar] [CrossRef]
  7. Roberts, D.P.; Mattoo, A.K. Sustainable crop production systems and human nutrition. Front. Sustain. Food Syst. 2019, 3, 72. [Google Scholar] [CrossRef]
  8. Bona, E.; Cantamessa, S.; Massa, N.; Manassero, P.; Marsano, F.; Copetta, A.; Lingua, G.; D’Agostino, G.; Gamalero, E.; Berta, G. Arbuscular mycorrhizal fungi and plant growth-promoting pseudomonads improve yield, quality and nutritional value of tomato: A field study. Mycorrhiza 2017, 27, 1–11. [Google Scholar] [CrossRef]
  9. Bona, E.; Todeschini, V.; Cantamessa, S.; Cesaro, P.; Copetta, A.; Lingua, G.; Gamalero, E.; Berta, G.; Massa, N. Combined bacterial and mycorrhizal inocula improve tomato quality at reduced fertilization. Sci. Hortic. 2018, 234, 160–165. [Google Scholar] [CrossRef]
  10. Roberts, D.; Mattoo, A. Sustainable agriculture—Enhancing environmental benefits, food nutritional quality and building crop resilience to abiotic and biotic stresses. Agriculture 2018, 8, 8. [Google Scholar] [CrossRef] [Green Version]
  11. Gonzali, S.; Perata, P. Anthocyanins from Purple Tomatoes as Novel Antioxidants to Promote Human Health. Antioxidants 2020, 9, 1017. [Google Scholar] [CrossRef]
  12. Luo, D.; Xiong, C.; Lin, A.; Zhang, C.; Sun, W.; Zhang, J.; Yang, C.; Lu, Y.; Li, H.; Ye, Z.; et al. SlBBX20 interacts with the COP9 signalosome subunit SlCSN5-2 to regulate anthocyanin biosynthesis by activating SlDFR expression in tomato. Hortic. Res. 2021, 8, 163. [Google Scholar] [CrossRef]
  13. Blando, F.; Berland, H.; Maiorano, G.; Durante, M.; Mazzucato, A.; Picarella, M.E.; Nicoletti, I.; Gerardi, C.; Mita, G.; Andersen, Ø.M. Nutraceutical Characterization of Anthocyanin-Rich Fruits Produced by “Sun Black” Tomato Line. Front. Nutr. 2019, 6, 133. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Zhang, Y.; Butelli, E.; De Stefano, R.; Schoonbeek, H.J.; Magusin, A.; Pagliarani, C.; Wellner, N.; Hill, L.; Orzaez, D.; Granell, A.; et al. Anthocyanins double the shelf life of tomatoes by delaying overripening and reducing susceptibility to gray mold. Curr. Biol. 2013, 23, 1094–1100. [Google Scholar] [CrossRef] [Green Version]
  15. Mazzucato, A.; Willems, D.; Bernini, R.; Picarella, M.E.; Santangelo, E.; Ruiu, F.; Tilesi, F.; Soressi, G.P. Novel phenotypes related to the breeding of purple-fruited tomatoes and effect of peel extracts on human cancer cell proliferation. Plant Phys. Biochem. 2013, 72, 125–133. [Google Scholar] [CrossRef]
  16. Tilesi, F.; Lombardi, A.; Mazzucato, A. Scientometric and Methodological Analysis of the Recent Literature on the Health-Related Effects of Tomato and Tomato Products. Foods 2021, 10, 1905. [Google Scholar] [CrossRef]
  17. Gupta, R.; Patil, R. Phyllospheric microbes: Diversity, functions, interaction, and applications in agriculture. In Current Trends in Microbial Biotechnology for Sustainable Agriculture; Yadav, A.N., Singh, J., Singh, C., Yadav, N., Eds.; Springer: Singapore, 2021; pp. 301–323. [Google Scholar] [CrossRef]
  18. Mamphogoro, T.P.; Babalola, O.O.; Aiyegoro, O.A. Exploitation of epiphytic bacterial antagonists for the management of post-harvest diseases of sweet pepper and other fresh produce—A viable option. Biocontrol Sci. Technol. 2020, 30, 741–761. [Google Scholar] [CrossRef]
  19. Berger, C.N.; Sodha, S.V.; Shaw, R.K.; Griffin, P.M.; Pink, D.; Hand, P.; Frankel, G. Fresh fruit and vegetables as vehicles for the transmission of human pathogens. Environ. Microbiol. 2010, 12, 2385–2397. [Google Scholar] [CrossRef]
  20. de São José, J.F.B.; Faria-Silva, L.; Lepaus, B.M. Decontamination of Vegetables. In Microbial Decontamination of Food; Shah, M.A., Mir, S.A., Eds.; Springer: Singapore, 2022; pp. 71–92. [Google Scholar] [CrossRef]
  21. Vandenkoornhuyse, P.; Quaiser, A.; Le Van Duhamel, M.; Dufresne, A. The importance of the microbiome of the plant holobiont. New Phytol. 2015, 206, 1196–1206. [Google Scholar] [CrossRef]
  22. Abadi, V.A.J.M.; Sepehri, M.; Rahmani, H.A.; Zarei, M.; Ronaghi, A.; Taghavi, S.M.; Shamshiripour, M. Role of dominant phyllosphere bacteria with plant growth–promoting characteristics on growth and nutrition of maize (Zea mays L.). J. Soil Sci. Plant Nutr. 2020, 20, 2348–2363. [Google Scholar] [CrossRef]
  23. Newman, D.J.; Cragg, G.M. Plant endophytes and epiphytes: Burgeoning sources of known and “unknown” cytotoxic and antibiotic agents? Planta Med. 2020, 86, 891–905. [Google Scholar] [CrossRef] [Green Version]
  24. Al-Maawali, S.S.; Al-Sadi, A.M.; Alsheriqi, S.A.K.; Al-Sabahi, J.N.; Velazhahan, R. The potential of antagonistic yeasts and bacteria from tomato phyllosphere and fructoplane in the control of Alternaria fruit rot of tomato. All Life 2021, 14, 34–48. [Google Scholar] [CrossRef]
  25. Hammami, R.; Oueslati, M.; Smiri, M.; Nefzi, S.; Ruissi, M.; Comitini, F.; Romanazzi, G.; Cacciola, S.O.; Sadfi Zouaoui, N. Epiphytic Yeasts and Bacteria as Candidate Biocontrol Agents of Green and Blue Molds of Citrus Fruits. J. Fungi 2022, 8, 818. [Google Scholar] [CrossRef] [PubMed]
  26. Campiglia, E.; Mancinelli, R.; Radicetti, E. Influence of no-tillage and organic mulching on tomato (Solanum lycopersicum L.) production and nitrogen use in the mediterranean environment of central Italy. Sci. Hortic. 2011, 130, 588–598. [Google Scholar] [CrossRef]
  27. Campiglia, E.; Radicetti, E.; Mancinelli, R. Cover crops and mulches influence weed management and weed flora composition in strip-tilled tomato (Solanum lycopersicum). Weed Res. 2015, 55, 416–425. [Google Scholar] [CrossRef] [Green Version]
  28. Ottesen, A.R.; Peña, A.G.; White, J.R.; Pettengill, J.B.; Li, C.; Allard, S.; Rideout, S.; Allard, M.; Hill, T.; Evans, P.; et al. Baseline survey of the anatomical microbial ecology of an important food plant: Solanum lycopersicum (tomato). BMC Microbiol. 2013, 13, 114. [Google Scholar] [CrossRef] [Green Version]
  29. Janisiewicz, W.J.; Ii, W.M.J.; Peter, K.A.; Kurtzman, C.P.; Buyer, J. Yeasts associated with plums and their potential for controlling brown rot after harvest. Yeast 2014, 31, 207–218. [Google Scholar] [CrossRef]
  30. Fernández-Suárez, R.; Ramírez-Villatoro, G.; Díaz-Ruiz, G.; Eslava, C.; Calderón, M.; Navarro-Ocaña, A.; Trejo-Márquez, A.; Wacher, C. Effect of postharvest UV-C treatment on the bacterial diversity of Ataulfo mangoes by PCR-DGGE, survival of E. coli and antimicrobial activity. Front. Microbiol. 2013, 4, 134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Huber, J.A.; Welch, D.B.M.; Morrison, H.G.; Huse, S.M.; Neal, P.R.; Butterfield, D.A.; Sogin, M.L. Microbial Population Structures in the Deep Marine Biosphere. Science 2007, 318, 97–100. [Google Scholar] [CrossRef]
  32. Wang, Y.; Qian, P.Y. Conservative Fragments in Bacterial 16S rRNA Genes and Primer Design for 16S Ribosomal DNA Amplicons in Metagenomic Studies. PLoS ONE 2009, 4, e7401. [Google Scholar] [CrossRef] [Green Version]
  33. Gorrasi, S.; Franzetti, A.; Ambrosini, R.; Pittino, F.; Pasqualetti, M.; Fenice, M. Spatio-temporal variation of the bacterial communities along a salinity gradient within a thalassohaline environment (Saline di Tarquinia salterns, Italy). Molecules 2021, 26, 1338. [Google Scholar] [CrossRef]
  34. Gorrasi, S.; Pasqualetti, M.; Braconcini, M.; Muñoz-Palazon, B.; Fenice, M. Could Pontimonas Harbour Halophilic Members Able to Withstand Very Broad Salinity Variations? Microorganisms 2022, 10, 790. [Google Scholar] [CrossRef] [PubMed]
  35. Muñoz-Palazon, B.; Rosa-Masegosa, A.; Vilchez-Vargas, R.; Link, A.; Gorrasi, S.; Gonzalez-Lopez, J.; Gonzalez-Martinez, A. Biological Removal Processes in Aerobic Granular Sludge for Treating Synthetic Hospital Wastewater: Effect of Temperature. J. Water Process Eng. 2022, 47, 102691. [Google Scholar] [CrossRef]
  36. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Claesson, M.J.; O’Sullivan, O.; Wang, Q.; Nikkilä, J.; Marchesi, J.R.; Smidt, H.; de Vos, W.M.; Ross, R.P.; O’Toole, P.W. Comparative analysis of pyrosequencing and a phylogenetic microarray for exploring microbial community structures in the human distal intestine. PLoS ONE 2009, 4, e6669. [Google Scholar] [CrossRef] [Green Version]
  38. Shannon, C.E. A Mathematical Theory of Communication. Bell. Sys. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  39. Pielou, E.C. The measurement of diversity in different types of biological collections. J. Theor. Biol. 1966, 13, 131–144. [Google Scholar] [CrossRef]
  40. Solow, A.R. A simple test for change in community structure. J. Anim. Ecol. 1993, 2, 191–193. [Google Scholar] [CrossRef]
  41. Seaby, R.M.; Henderson, P.A. Species Diversity and Richness, Version 4; Pisces Conservation Ltd.: Lymington, UK, 2006. [Google Scholar]
  42. Clarke, K.R.; Somerfield, P.J.; Gorley, R.N. Testing of null hypotheses in exploratory community analyses: Similarity profiles and biota-environment linkage. J. Exp. Mar. Biol. Ecol. 2008, 366, 56–69. [Google Scholar] [CrossRef]
  43. Heberle, H.; Meirelles, G.V.; da Silva, F.R.; Telles, G.P.; Minghim, R. InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinform. 2015, 16, 169. [Google Scholar] [CrossRef] [Green Version]
  44. Leff, J.W.; Fierer, N. Bacterial Communities Associated with the Surfaces of Fresh Fruits and Vegetables. PLoS ONE 2013, 8, e59310. [Google Scholar] [CrossRef] [Green Version]
  45. Constancias, F.; Saby, N.P.A.; Terrat, S.; Dequiedt, S.; Horrigue, W.; Nowak, V.; Guillemin, J.-P.; Biju-Duval, L.; Prévost-Bouré, N.C.; Ranjard, L. Contrasting Spatial Patterns and Ecological Attributes of Soil Bacterial and Archaeal Taxa across a Landscape. MicrobiologyOpen 2015, 4, 518–531. [Google Scholar] [CrossRef]
  46. Dong, C.-J.; Wang, L.-L.; Li, Q.; Shang, Q.-M. Bacterial communities in the rhizosphere, phyllosphere and endosphere of tomato plants. PLoS ONE 2019, 14, e022384. [Google Scholar] [CrossRef] [PubMed]
  47. Bill, M.; Chidamba, L.; Gokul, J.K.; Korsten, L. Mango endophyte and epiphyte microbiome composition during fruit development and post-harvest stages. Horticulturae 2021, 7, 495. [Google Scholar] [CrossRef]
  48. Kumar, A.; Zhimo, Y.; Biasi, A.; Salim, S.; Feygenberg, O.; Wisniewski, M.; Droby, S. Endophytic Microbiome in the Carposphere and Its Importance in Fruit Physiology and Pathology. In Postharvest Pathology: Plant Pathology in the 21st Century; Spadaro, D., Droby, S., Gullino, M.L., Eds.; Springer: Cham, Switzerland, 2021; Volume 11, pp. 73–88. [Google Scholar] [CrossRef]
  49. Berg, G.; Erlacher, A.; Grube, M. The Edible Plant Microbiome: Importance and Health Issues. In Principles of Plant-Microbe Interactions: Microbes for Sustainable Agriculture; Lugtenberg, B., Ed.; Springer International Publishing: Cham, Switzerland, 2015; pp. 419–426. [Google Scholar] [CrossRef]
  50. Ottesen, A.; Ramachandran, P.; Reed, E.; Gu, G.; Gorham, S.; Ducharme, D.; Newell, M.; Rideout, S.; Turini, T.; Hill, T.; et al. Metagenome tracking biogeographic agroecology: Phytobiota of tomatoes from virginia, maryland, north carolina and california. Food Microbiol. 2019, 79, 132–136. [Google Scholar] [CrossRef]
  51. Kämpfer, P. Buttiauxella. In Bergey’s Manual of Systematics of Archaea and Bacteria; Whitman, W.E.B., Ed.; John Wiley & Sons: New York, NY, USA, 2015; pp. 1–16. [Google Scholar] [CrossRef]
  52. Gorrasi, S.; Pasqualetti, M.; Franzetti, A.; Gonzalez-Martinez, A.; Gonzalez-Lopez, J.; Muñoz-Palazon, B.; Fenice, M. Persistence of Enterobacteriaceae drawn into a marine saltern (Saline di Tarquinia, Italy) from the adjacent coastal zone. Water 2021, 13, 1443. [Google Scholar] [CrossRef]
  53. Wu, K.; Luo, J.; Li, J.; An, Q.; Yang, X.; Liang, Y.; Li, T. Endophytic bacterium Buttiauxella sp. SaSR13 improves plant growth and cadmium accumulation of hyperaccumulator Sedum alfredii. Environ. Sci. Pollut. Res. 2018, 25, 21844–21854. [Google Scholar] [CrossRef]
  54. Almasia, R.; Henríquez, M.; Levican, A.; Poblete-Morales, M. Genome sequence of a potentially new Buttiauxella species, Strain B2, isolated from rhizosphere of olivillo trees (Aextoxicon punctatum). Microbiol. Resour. Announc. 2020, 9, e01351-19. [Google Scholar] [CrossRef] [Green Version]
  55. Xin, Y.-Y.; Rahman, A.; Li, H.-X.; Xu, T.; Ding, G.-C.; Li, J. Modification of total and phosphorus mineralizing bacterial communities associated with Zea mays L. through plant development and fertilization regimes. J. Integr. Agric. 2021, 20, 3026–3038. [Google Scholar] [CrossRef]
  56. Sahu, K.P.; Patel, A.; Kumar, M.; Sheoran, N.; Mehta, S.; Reddy, B.; Eke, P.; Prabhakaran, N.; Kumar, A. Integrated metabarcoding and culturomic-based microbiome profiling of rice phyllosphere reveal diverse and functional bacterial communities for blast disease suppression. Front. Microbiol. 2021, 12, 780458. [Google Scholar] [CrossRef] [PubMed]
  57. Ottesen, A.R.; Gorham, S.; Pettengill, J.B.; Rideout, S.; Evans, P.; Brown, E. The impact of systemic and copper pesticide applications on the phyllosphere microflora of tomatoes. J. Sci. Food Agric. 2015, 95, 1116–1125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Mechan Llontop, M.E.; Tian, L.; Sharma, P.; Heflin, L.; Bernal-Galeano, V.; Haak, D.C.; Clarke, C.R.; Vinatzer, B.A. Experimental evidence pointing to rain as a reservoir of tomato phyllosphere microbiota. Phytobiomes J. 2021, 5, 382–399. [Google Scholar] [CrossRef]
  59. Hata, H.; Natori, T.; Mizuno, T.; Kanazawa, I.; Eldesouky, I.; Hayashi, M.; Miyata, M.; Fukunaga, H.; Ohji, S.; Hosoyama, A.; et al. Phylogenetics of family Enterobacteriaceae and proposal to reclassify Escherichia hermannii and Salmonella subterranea as Atlantibacter hermannii and Atlantibacter subterranea gen. Nov., Comb. Nov. Microbiol. Immunol. 2016, 60, 303–311. [Google Scholar] [CrossRef]
  60. Janda, J.M.; Abbott, S.L. The Changing Face of the Family Enterobacteriaceae (Order: “Enterobacterales”): New Members, Taxonomic Issues, Geographic Expansion, and New Diseases and Disease Syndromes. Clin. Microbiol. Rev. 2021, 34, e00174-20. [Google Scholar] [CrossRef]
  61. Shelobolina, E.S.; Sullivan, S.A.; O’Neill, K.R.; Nevin, K.P.; Lovley, D.R. Isolation, characterization, and U (VI)-reducing potential of a facultatively anaerobic, acid-resistant bacterium from low-pH, nitrate- and U (VI)-contaminated subsurface sediment and description of Salmonella subterranea sp. nov. Appl. Environ. Microbiol. 2004, 70, 2959–2965. [Google Scholar] [CrossRef] [Green Version]
  62. Watanapokasin, R.Y.; Boonyakamol, A.; Sukseree, S.; Krajarng, A.; Sophonnithiprasert, T.; Kanso, S.; Imai, T. Hydrogen production and anaerobic decolorization of wastewater containing Reactive Blue 4 by a bacterial consortium of Salmonella subterranea and Paenibacillus polymyxa. Biodegradation 2009, 20, 411–418. [Google Scholar] [CrossRef]
  63. Higdon, S.M.; Pozzo, T.; Kong, N.; Huang, B.C.; Yang, M.L.; Jeannotte, R.; Brown, C.T.; Bennett, A.B.; Weimer, B.C. Genomic Characterization of a Diazotrophic Microbiota Associated with Maize Aerial Root Mucilage. PLoS ONE 2020, 15, e0239677. [Google Scholar] [CrossRef]
  64. Besaury, L.; Rémond, C. Culturable and metagenomic approaches of wheat bran and wheat straw phyllosphere’s highlight new lignocellulolytic microorganisms. Lett. Appl. Microbiol. 2022, 74, 840–850. [Google Scholar] [CrossRef]
  65. Hu, H.; Liu, Y.; Huang, Y.; Zhang, Z.; Tang, H. The Leaf Microbiome of Tobacco Plants across Eight Chinese Provinces. Microorganisms 2022, 10, 450. [Google Scholar] [CrossRef]
  66. PLaBAse-PLant-Associated BActeria Web Resource, Atlantibacter. Available online: https://plabase.informatik.uni-tuebingen.de/pb/plaba_db_q.php (accessed on 1 October 2022).
  67. Pati, B.R.; Chandra, A.K. Diazotrophic bacterial population and other associated organisms on the phyllosphere of some crop plants. Zentralbl. Mikrobiol. 1993, 148, 392–402. [Google Scholar] [CrossRef]
  68. Pande, B.N.; Krysinska-Traczyk, E.; Prazmo, Z.; Skórska, C.; Sitkowska, J.; Dutkiewicz, J. Occupational biohazards in agricultural dusts from India. Ann. Agric. Environ. Med. 2000, 7, 133–139. [Google Scholar] [PubMed]
  69. Stephan, R.; Van Trappen, S.; Cleenwerck, I.; Vancanneyt, M.; De Vos, P.; Lehner, A. Enterobacter turicensis sp. nov. and Enterobacter helveticus sp. nov., isolated from fruit powder. Int. J. Syst. Evol. Microbiol. 2007, 57, 820–826. [Google Scholar] [CrossRef] [Green Version]
  70. Zhang, G.X.; Peng, G.X.; Wang, E.T.; Yan, H.; Yuan, Q.H.; Zhang, W.; Lou, X.; Wu, H.; Tan, Z.Y. Diverse endophytic nitrogen-fixing bacteria isolated from wild rice Oryza rufipogon and description of Phytobacter diazotrophicus gen. nov. sp. nov. Arch. Microbiol. 2007, 189, 431–439. [Google Scholar] [CrossRef] [PubMed]
  71. Osdaghi, E.; Rahimi, T.; Taghavi, S.M.; Ansari, M.; Zarei, S.; Portier, P.; Briand, M.; Jacques, M.-A. Comparative genomics and phylogenetic analyses suggest several novel species within the genus Clavibacter, including nonpathogenic tomato-associated strains. Appl. Environ. Microbiol. 2020, 86, e02873-19. [Google Scholar] [CrossRef] [Green Version]
  72. Laroute, V.; Tormo, H.; Couderc, C.; Mercier-Bonin, M.; Le Bourgeois, P.; Cocaign-Bousquet, M.; Daveran-Mingot, M.-L. From Genome to Phenotype: An Integrative Approach to Evaluate the Biodiversity of Lactococcus lactis. Microorganisms 2017, 5, 27. [Google Scholar] [CrossRef]
  73. Latt, Z.K.; San, S.Y.; Kyaw, E.P.; Lynn, T.M.; Nwe, M.T.; Mon, W.W.; Aye, K.N. Using Cellulolytic Nitrogen Fixing Bacterium, Azomonas agilis for effective degradation of agricultural residues. Open Microbiol. J. 2018, 12, 154–162. [Google Scholar] [CrossRef]
  74. Compant, S.; Samad, A.; Faist, H.; Sessitsch, A. A review on the plant microbiome: Ecology, functions, and emerging trends in microbial application. J. Adv. Res. 2019, 19, 29–37. [Google Scholar] [CrossRef]
  75. Madhaiyan, M.; Saravanan, V.S.; Blom, J.; Smits, T.H.; Rezzonico, F.; Kim, S.J.; Weon, H.Y.; Kwon, S.W.; Whitman, W.B.; Ji, L. Phytobacter palmae sp. nov., a novel endophytic, N2 fixing, plant growth promoting Gammaproteobacterium isolated from oil palm (Elaeis guineensis Jacq.). Int. J. Syst. Evol. Microbiol. 2020, 70, 841–848. [Google Scholar] [CrossRef]
  76. Trias, R.; Bañeras, L.; Montesinos, E.; Badosa, E. Lactic acid bacteria from fresh fruit and vegetables as biocontrol agents of phytopathogenic bacteria and fungi. Int. J. Microbiol. 2008, 11, 231–236. [Google Scholar] [CrossRef] [Green Version]
  77. Trotel-Aziz, P.; Couderchet, M.; Biagianti, S.; Aziz, A. Characterization of new bacterial biocontrol agents Acinetobacter, Bacillus, Pantoea and Pseudomonas spp. mediating grapevine resistance against Botrytis cinerea. Environ. Exp. Bot. 2008, 64, 21–32. [Google Scholar] [CrossRef]
  78. Mansfield, J.; Genin, S.; Magori, S.; Citovsky, V.; Sriariyanum, M.; Ronald, P.; Dow, M.; Verdier, V.; Beer, S.V.; Machado, M.A.; et al. Top 10 plant pathogenic bacteria in molecular plant pathology. Mol. Plant Pathol. 2012, 13, 614–629. [Google Scholar] [CrossRef] [Green Version]
  79. Thapa, S.P.; Davis, E.W.; Lyu, Q.; Weisberg, A.J.; Stevens, D.M.; Clarke, C.R.; Coaker, G.; Chang, J.H. The Evolution, Ecology, and Mechanisms of Infection by Gram-Positive, Plant-Associated Bacteria. Annu. Rev. Phytopathol. 2019, 57, 341–365. [Google Scholar] [CrossRef] [PubMed]
  80. Telias, A.; White, J.R.; Pahl, D.M.; Ottesen, A.R.; Walsh, C.S. Bacterial community diversity and variation in spray water sources and the tomato fruit surface. BMC Microbiol. 2011, 11, 81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Enya, J.; Shinohara, H.; Yoshida, S.; Tsukiboshi, T.; Negishi, H.; Suyama, K.; Tsushima, S. Culturable leaf-associated bacteria on tomato plants and their potential as biological control agents. Microb. Ecol. 2007, 53, 524–536. [Google Scholar] [CrossRef]
  82. Mohammed, A.F.; Oloyede, A.R.; Odeseye, A.O. Biological control of bacterial wilt of tomato caused by Ralstonia solanacearum using Pseudomonas species isolated from the rhizosphere of tomato plants. Arch. Phytopathol. Plant Prot. 2020, 53, 1–16. [Google Scholar] [CrossRef]
  83. Abo-Elyousr, K.A.M.; Hassan, S.A. Biological control of Ralstonia solanacearum (Smith), the causal pathogen of bacterial wilt disease by using Pantoea spp. Egypt. J. Biol. Pest. Control 2021, 31, 113. [Google Scholar] [CrossRef]
  84. Ghadamgahi, F.; Tarighi, S.; Taheri, P.; Saripella, G.V.; Anzalone, A.; Kalyandurg, P.B.; Catara, V.; Ortiz, R.; Vetukuri, R.R. Plant Growth-Promoting Activity of Pseudomonas aeruginosa FG106 and Its Ability to Act as a Biocontrol Agent against Potato, Tomato and Taro Pathogens. Biology 2022, 11, 140. [Google Scholar] [CrossRef]
  85. Ashmawy, N.A.; Shoeib, A.A.; Youssef, H.F.B.; Mahmoud, S.M. Pathological, biochemical and molecular characterization of the seed-borne bacteria “Pantoea spp., Xanthomonas spp. and Pseudomonas spp.” from solanaceous plants in Egypt. J. Microbiol. Biotech. Food Sci. 2020, 10, 289–295. [Google Scholar] [CrossRef]
  86. Tian, Q.; Chuan, J.; Sun, X.; Zhou, A.; Wang, L.; Zou, J.; Zhao, W.; Li, X. Description of Clavibacter zhangzhiyongii sp. nov., a phytopathogenic actinobacterium isolated from barley seeds, causing leaf brown spot and decline. Int. J. Syst. Evol. Microbiol. 2021, 71, 004786. [Google Scholar] [CrossRef]
Figure 1. Bacterial composition at the phylum level. Distribution of the major phyla (Ra ≥ 1% in at least one sample); phyla with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 1. Bacterial composition at the phylum level. Distribution of the major phyla (Ra ≥ 1% in at least one sample); phyla with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g001
Figure 2. Bacterial composition at the family level. Distribution of the major families (Ra ≥ 1% in at least one sample); families with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 2. Bacterial composition at the family level. Distribution of the major families (Ra ≥ 1% in at least one sample); families with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g002
Figure 3. Bacterial composition at the genus level. Distribution of the major genera (Ra ≥ 1% in at least one sample); genera with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 3. Bacterial composition at the genus level. Distribution of the major genera (Ra ≥ 1% in at least one sample); genera with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g003
Figure 4. Bacterial composition at the OTU level. Distribution of the major OTUs (Ra ≥ 1% in at least one sample); OTUs with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 4. Bacterial composition at the OTU level. Distribution of the major OTUs (Ra ≥ 1% in at least one sample); OTUs with Ra < 1% were gathered in “Others”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g004
Figure 5. Venn diagram depicting the number of exclusive and shared OTUs among the investigated experimental conditions. The total number of OTUs found in each experimental condition is reported in square brackets. Both the number of exclusive/shared OTUs and the corresponding percentage (in round brackets) on the total of the OTUs are reported. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 5. Venn diagram depicting the number of exclusive and shared OTUs among the investigated experimental conditions. The total number of OTUs found in each experimental condition is reported in square brackets. Both the number of exclusive/shared OTUs and the corresponding percentage (in round brackets) on the total of the OTUs are reported. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g005
Figure 6. Distribution of the core microbiota OTUs across the studied experimental conditions. For each OTU, the taxonomic annotation is provided. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 6. Distribution of the core microbiota OTUs across the studied experimental conditions. For each OTU, the taxonomic annotation is provided. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g006
Figure 7. Principal component analysis (PCA) based on the bacterial OTU found across the investigated experimental conditions. Four plot visualizations of the same PCA analysis are provided to better show the OTU occurrence in the four experimental conditions; OTU abundance variation in T1 (a), T2 (b), T3 (c), and T4 (d). Circle size is proportional to the OTU relative abundance in the samples; OTU absence in an experimental condition is indicated by “+”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Figure 7. Principal component analysis (PCA) based on the bacterial OTU found across the investigated experimental conditions. Four plot visualizations of the same PCA analysis are provided to better show the OTU occurrence in the four experimental conditions; OTU abundance variation in T1 (a), T2 (b), T3 (c), and T4 (d). Circle size is proportional to the OTU relative abundance in the samples; OTU absence in an experimental condition is indicated by “+”. T1: WT/OF, T2: WT/CF, T3: SB/OF, T4: SB/CF.
Microorganisms 10 02240 g007
Figure 8. Cluster analysis dendrogram showing similarities/dissimilarities among the experimental conditions based on their OTU composition. The clustering analysis (complete linkage) was performed on Bray–Curtis dissimilarity data, applying SIMPROF routine to test significance. Solid black lines indicate samples that are significantly different (p < 0.05), whereas dashed red lines indicate groups of samples that are not significantly different (p > 0.05), according to SIMPROF test.
Figure 8. Cluster analysis dendrogram showing similarities/dissimilarities among the experimental conditions based on their OTU composition. The clustering analysis (complete linkage) was performed on Bray–Curtis dissimilarity data, applying SIMPROF routine to test significance. Solid black lines indicate samples that are significantly different (p < 0.05), whereas dashed red lines indicate groups of samples that are not significantly different (p > 0.05), according to SIMPROF test.
Microorganisms 10 02240 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gorrasi, S.; Pasqualetti, M.; Muñoz-Palazon, B.; Novello, G.; Mazzucato, A.; Campiglia, E.; Fenice, M. Comparison of the Peel-Associated Epiphytic Bacteria of Anthocyanin-Rich “Sun Black” and Wild-Type Tomatoes under Organic and Conventional Farming. Microorganisms 2022, 10, 2240. https://doi.org/10.3390/microorganisms10112240

AMA Style

Gorrasi S, Pasqualetti M, Muñoz-Palazon B, Novello G, Mazzucato A, Campiglia E, Fenice M. Comparison of the Peel-Associated Epiphytic Bacteria of Anthocyanin-Rich “Sun Black” and Wild-Type Tomatoes under Organic and Conventional Farming. Microorganisms. 2022; 10(11):2240. https://doi.org/10.3390/microorganisms10112240

Chicago/Turabian Style

Gorrasi, Susanna, Marcella Pasqualetti, Barbara Muñoz-Palazon, Giorgia Novello, Andrea Mazzucato, Enio Campiglia, and Massimiliano Fenice. 2022. "Comparison of the Peel-Associated Epiphytic Bacteria of Anthocyanin-Rich “Sun Black” and Wild-Type Tomatoes under Organic and Conventional Farming" Microorganisms 10, no. 11: 2240. https://doi.org/10.3390/microorganisms10112240

APA Style

Gorrasi, S., Pasqualetti, M., Muñoz-Palazon, B., Novello, G., Mazzucato, A., Campiglia, E., & Fenice, M. (2022). Comparison of the Peel-Associated Epiphytic Bacteria of Anthocyanin-Rich “Sun Black” and Wild-Type Tomatoes under Organic and Conventional Farming. Microorganisms, 10(11), 2240. https://doi.org/10.3390/microorganisms10112240

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop