- Data note
- Open Access
Whole genome sequencing and taxonomic profiling of two Pantoea sp. isolated from environmental samples in Israel
BMC Genomic Data volume 23, Article number: 31 (2022)
As part of a research aiming at the isolation of bacteria secreting growth inhibiting compounds, cultures of Francisella tularensis were implanted in environmental samples and monitored for inhibition zones on agar. Two antibiotic-like secreting bacteria were isolated, their genomic sequence was deciphered and taxonomic profiling analysis classified them as belonging to the Pantoea genus.
Two bacterial isolates exhibiting growth inhibition zones to F. tularensis (LVS) were analyzed using the Oxford Nanopore Technology (ONT). Preliminary de novo assembly of the reads was performed, followed by taxonomic profiling based on Multi Locus Sequence Analysis (MLSA) and implementation of the Average Nucleotide Identity (ANI) measure. The genomic sequences resulted in the identification of two different Pantoea species, denoted EnvD and EnvH.
Subsequent de novo genome assembly generated 5 and 10 contigs for EnvD and EnvH, respectively. The largest contig (4,008,183 bps and 3,740,753 bps for EnvD and EnvH, respectively), overlaps to a major extent to the chromosome of closely related Pantoea species. ANI values calculated for both isolates revealed two apparently new species of the Pantoea genus.
Our study deciphered the identity of two bacteria producing antibiotic-like compounds, and the genomic sequence revealed they represent distinct Pantoea species.
Pantoea is a genus of Gram-negative bacteria belonging to the family Erwiniaceae. These bacteria were isolated from a multitude of environmental sources (such as plants, animals, water, soil and human). They are distributed in nature, and over 25 species were already documented [1,2,3,4]. In this study we present the taxonomic profiling of two environmental isolates which were identified by their ability to inhibit growth of F. tularensis (LVS) colonies. The isolates were sequenced by long read sequencing and assembled by a de novo assembly. We calculated the ANI measure to establish the species relatedness and the MLSA method was implemented to generate a phylogenetic tree showing these variants as new Pantoea species.
In this study, over 20 environmental samples were collected from distinct different geographic locations in Israel (natural sources). The samples were prepared by wiping a 20 cm2 area of urban asphalt roads with 3 sterile cotton swabs (Copan Italia SPA) moistened aseptically in sterile phosphate-buffered saline (PBS) solution. The swabs were discarded after extraction by vortex in 5 mL PBS. After 2 min. of passive sedimentation of dust and dirt, 4.2 mL were transferred to new sterile tube. An aliquot of 0.1 mL of each environmental sample was mixed with 0.1 mL of PBS containing 108 cfu/ml of indicator bacteria (F. tularensis LVS), the mix was spread on cysteine heart agar with hemoglobin (CHA, Difco), plates were incubated for 2 days at 37 °C and monitored for colonies demonstrating inhibition zone for the indicator strain. Several such colonies were further isolated to pure cultures on CHA, and their inhibition to LVS was verified. Two of the isolates showed differential capability to grow in Brain Heart Infusion broth, and were denoted herein as EnvD (from Talpiot industrial zone, Jerusalem) and EnvH (from Romema neighborhood, Haifa). For DNA purification, EnvD and EnvH colonies were grown on CHA, suspended to high density in 1 mL PBS, mixed 1:1 with ATL buffer, heated at 100 °C for 30 min before DNA purification using the QIAamp DNA blood minikit (Qiagen). Whole genome sequencing was conducted in GenoHub facility (https://genohub.com) using the Oxford Nanopore MinION Technology. Libraries for both EnvD and EnvH were prepared using the SQK-LSK109 ligation sequencing kit (Oxford Nanopore Technologies) and sequencing was conducted using an R9 flow cell. Basecalling was performed using Guppy v 4.4.2. A total of 210,000 reads for EnvD (mean average length of 3854 bps) and 190,000 reads for EnvH (mean average length of 1881 bps) were obtained, resulting in a coverage of 155x and 89x, respectively (data file 1) .
De novo assembly was conducted by implementing Flye , designed for the assembly of long reads generated by ONT. Five contigs were obtained for EnvD, the largest harboring a length of 4,008,183 bps, the N50 value being therefore 4,008,183. The EnvH generated reads were assembled to a total of 10 contigs, the longest being 3,740,753 bps (which is therefore the N50 value). A rough and preliminary estimation of the taxonomical relatedness of the two isolates conducted by a Blast analysis  against the nucleotide database (nt, https://www.ncbi.nlm.nih.gov/nucleotide/), disclosed similarity to the Pantoea genus. Subsequently, the longest contig of each isolate was compared to all publicly available Pantoea sequences (National Center for Biotechnology Information, NCBI), revealing that the best matching hits for EnvD and EnvH are Pantoea agglomerans and Pantoea stewartii, respectively. Accordingly, the contigs generated by the assembly were aligned to their closest species, using Mauve . The alignments display a nearly complete coverage of the reference chromosome by the largest contig, both for EnvD and EnvH. In addition, two of the remaining contigs overlap plasmid regions in the reference genomes (data file 2) .
To further characterize whether the sequenced genomes can be associated with one of the already known Pantoea species, we used the ANI measure as a well-established whole genome similarity metrics [10,11,12,13,14]. The ANI value was estimated using the FastANI algorithm . Pairwise ANI values of at most 82% were obtained for each of the sequenced genomes with genomes representing known Pantoea species . According to the generally accepted cutoff value of 95% used as a boundary for species delineation [14, 16], it appears that EnvD and EnvH constitute new species within the Pantoea genus. To note, the pairwise ANI value between EnvD and EnvH is 83%, therefore representing two distinct lineages.
To assign the taxonomic profiling of the two isolates, we implemented the Multi Locus Sequence Analysis (MLSA) typing method, tailored for phylogeny analysis of Pantoea species [17, 18], using five core genes that are effective at species-level delineation of the genus Pantoea: fusA, gyrB, leuS, pyrG and rpoB . The sequence of the genes orthologous to EnvD and EnvH were extracted from the contig sequences covering the chromosome regions (contig 4 for EnvD and contig 3 for EnvH), and concatenated into a mini-gene. Alongside, the five protein-coding genes were extracted from 37 representative, reference and/or type strains of Pantoea and Tatumella species (for a complete list of species included in the set of 37 genomes, refer to ). The multiple alignment of the 37 concatenated sequences together with the concatenated sequences originating from EnvD and EnvH was constructed using the MAFFT algorithm (Multiple Alignment using fast Fourier Transform)  of the MegAlign™ Pro (©1993–2020) (DNASTAR®). Phylogenetic analysis of the aligned sequences was conducted using the PhyML tool , which estimates Maximum-Likelihood phylogenies (NGPhylogeny.fr ). Each of the new isolates, EnvD and EnvH, forms a distinct and separate branch (data file 3) . While EnvH is related to Pantoea stewartii, Pantoea ananatis and Pantoea allii species, EnvD is not only branching from a separate cluster, but is also relatively distant from species in this cluster (which includes, among others, Pantoea agglomerans, corroborating with the preliminary assignment described above for EnvD).
To conclude, two growth-inhibiting bacteria from environmental samples collected from two distinct areas in Israel were identified and assigned as belonging to the Pantoea genus. Their taxonomical profiling unveiled that these bacteria can be classified as new Pantoea species diverging from known Pantoea species described up to date. Further experimental characterization of the mechanism involved in the LVS growth inhibition of these Pantoea is now undertaken. Please see Table 1 for links to Data files 1-4 and Data sets 1-3. 1.
The de novo assembly resulted in a number of contigs; while the longest contig very nearly cover the chromosome region, some other regions of the EnvD and EnvH sequenced genomes are fragmentally covered.
By reason of practical considerations, the phylogeny relatedness of the two Pantoea sp. EnvD and EnvH is based on a limited dataset of 37 known Pantoea species, comprising of representative, reference and/or type strains of Pantoea and Tatumella species.
Availability of data and materials
Data files 1–4 described in this Data note can be freely and openly accessible on Figshare (https://figshare.com/) [5, 9, 22]. Data sets 1–3 are available on the NCBI database. The raw reads have been submitted to the NCBI Sequence Read Archive under the accession number SRP316834 for Pantoea sp. EnvD and Pantoea sp. EnvH (Data set 1) [23,24,25]. The genome assemblies of the two samples were submitted to NCBI GenBank and are available under the accession number JAGTWO000000000 for the Pantoea sp. EnvD and the accession number JAGTWN000000000 for the Pantoea sp. EnvH (Data sets 2–3) [24, 25].
Multi Locus Sequence Analysis
Average Nucleotide Identity
Oxford Nanopore Technologies
Multiple Alignment using fast Fourier Transform
Tambong JT. Taxogenomics and systematics of the genus Pantoea. Front Microbiol. 2019;10:2463.
Agarwal G, Choudhary D, Stice SP, MyersGitaitis BKRD, Venter SN, et al. Pan-genome-wide analysis of Pantoea ananatis identified genes linked to pathogenicity in onion. Front Microbiol. 2021;12:684756.
Agarwal G, Gitaitis RD, Dutta B. Pan-Genome of novel Pantoea stewartii subsp indologenes reveals genes involved in onion pathogenicity and evidence of lateral gene transfer. Microorganisms. 2021;9(8):1761.
Stice SP, Shin GY, De Armas S, Koirala S, Galvan GA, Siri MI, et al. The distribution of onion virulence gene clusters among Pantoea spp. Front Plant Sci. 2021;12:643787.
Data file 1: Sequencing and assembly metrics. Figshare: https://doi.org/10.6084/m9.figshare.15111765.v1 (2021).
Kolmogorov M, Bickhart DM, Behsaz B, Gurevich A, Rayko M, Shin SB, et al. metaFlye: scalable long-read metagenome assembly using repeat graphs. Nat Methods. 2020;17(11):1103–10.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Darling AC, Mau B, Blattner FR, Perna NT. Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004;14(7):1394–403.
Data file 2: Alignment of contigs to a reference genome. Figshare: https://doi.org/10.6084/m9.figshare.15105198.v1 (2021).
Figueras MJ, Beaz-Hidalgo R, Hossain MJ, Liles MR. Taxonomic affiliation of new genomes should be verified using average nucleotide identity and multilocus phylogenetic analysis. Genome Announc. 2014;2(6):e00927-14.
Arahal DR. Whole-genome analyses: Average Nucleotide Identity. Methods in Microbiology. 2014;41:103–22.
Jain C, Rodriguez RL, Phillippy AM, Konstantinidis KT, Aluru S. High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat Commun. 2018;9(1):5114.
Konstantinidis KT, Tiedje JM. Genomic insights that advance the species definition for prokaryotes. Proc Natl Acad Sci U S A. 2005;102(7):2567–72.
Goris J, Konstantinidis KT, Klappenbach JA, Coenye T, Vandamme P, Tiedje JM. DNA-DNA hybridization values and their relationship to whole-genome sequence similarities. Int J Syst Evol Microbiol. 2007;57(Pt 1):81–91.
Data file 3: Table of ANI scores for Pantoea sp. EnvD and Pantoea sp. EnvH. Figshare: https://doi.org/10.6084/m9.figshare.19204658.v2 (2022).
Richter M, Rossello-Mora R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009;106(45):19126–31.
Tambong JT, Xu R, Kaneza CA, Nshogozabahizi JC. An in-depth analysis of a multilocus phylogeny identifies leuS As a reliable phylogenetic marker for the genus Pantoea. Evol Bioinform Online. 2014;10:115–25.
Deletoile A, Decre D, Courant S, Passet V, Audo J, Grimont P, et al. Phylogeny and identification of Pantoea species and typing of Pantoea agglomerans strains by multilocus gene sequencing. J Clin Microbiol. 2009;47(2):300–10.
Katoh K, Standley DM. MAFFT: iterative refinement and additional methods. Methods Mol Biol. 2014;1079:131–46.
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.
Lemoine F, Correia D, Lefort V, Doppelt-Azeroual O, Mareuil F, Cohen-Boulakia S, et al. NGPhylogeny.fr: new generation phylogenetic services for non-specialists. Nucleic Acids Res. 2019;47(W1):W260–5.
Data file 4: Phylogenetics analyses. Figshare: https://doi.org/10.6084/m9.figshare.15111588.v3 (2021).
National Center for Biotechnology Information. Sequence Read Archive. https://trace.ncbi.nlm.nih.gov/Traces/sra/?study=SRP316834 (2021).
National Center for Biotechnology Information. Assembly. https://www.ncbi.nlm.nih.gov/Traces/wgs/JAGTWO01?display=contigs (2021).
National Center for Biotechnology Information. Assembly. https://www.ncbi.nlm.nih.gov/Traces/wgs/JAGTWN01?display=contigs (2021).
We thank Dr. Adi Beth-Din for fruitful discussions.
This research received no external funding.
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.
About this article
Cite this article
Guedj-Dana, Y., Cohen-Gihon, I., Israeli, O. et al. Whole genome sequencing and taxonomic profiling of two Pantoea sp. isolated from environmental samples in Israel. BMC Genom Data 23, 31 (2022). https://doi.org/10.1186/s12863-022-01049-7
- Environmental samples
- Whole genome sequencing
- Oxford nanopore
- De novo assembly