Next Article in Journal
Dietary Patterns and Oral Health Behaviours Associated with Caries Development from 4 to 7 Years of Age
Previous Article in Journal
Chemical-Mediated Targeted Protein Degradation in Neurodegenerative Diseases
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Relationship of the Pine Growth Promoting Pantoea eucalypti FBS135 with Type Strains P. eucalypti LMG 24197T and P. vagans 24199T

College of Forestry, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
These authors have contributed equally to this work and share first authorship.
Life 2021, 11(7), 608; https://doi.org/10.3390/life11070608
Submission received: 6 May 2021 / Revised: 10 June 2021 / Accepted: 11 June 2021 / Published: 24 June 2021
(This article belongs to the Section Plant Science)

Abstract

:
Endophytes in woody plants are much less understood. Pantoea strain FBS135 is an endophytic bacterium isolated from Pinus massoniana with the ability to promote pine growth significantly. In this study, we demonstrated that FBS135 has the astonishing ability of low nitrogen tolerance but no ability of nitrogen fixation. To exactly determine the phylogenetic status of FBS135, we sequenced the whole genomes of P. eucalypti LMG 24197T and P. vagans 24199T, type strains of two Pantoea species, which are evolutionarily closest to FBS135. P. eucalypti LMG 24197T contained a single chromosome of 4,035,995 bp (C+G, 54.6%) plus three circular plasmids while LMG 24199T comprises a single circular chromosome of 4,050,173 bp (C+G, 55.6%) and two circular plasmids. With the genomic information, FBS135 was finally identified as a P. eucalypti strain, although it showed some different physiological traits from the two type strains. Comparative genomic analyses were performed for the three strains, revealing their common molecular basis associated with plant lifecycle as well as the differences in their gene arrangements relating to nitrogen utilization.

1. Introduction

Having been studied for decades, plant endophytes generally refer to those microorganisms that can colonize intercellular spaces or cells of plants and establish a symbiotic relationship with their host plants [1]. In recent years, Hallmann et al. raised a more practical definition of endophyte: Microorganisms separated from plant tissues that have been subjected to strict surface sterilization can be regarded as endophytes so long as they do not cause disease to the host plants [2,3]. Endophytes are widespread in plants. It is assumed that “plants without plant endophytes are a rare exception” [4].
Endophytes can promote plant growth through direct or indirect actions [5,6,7,8,9]. For example, endophytes can directly promote plant growth by producing and regulating various plant hormones such as auxin, cytokinin, or ethylene. Some endophytes can synthesize 1-aminocyclopropane-1-carboxylate (ACC) deaminase, which cleave to ACC in higher plants, thus reducing endogenous ethylene level. In addition, a large group of endophytes can obtain nitrogen, phosphorus, or iron from the environment and share with their hosts. Indirect beneficial mechanisms of endophytes not only include the production of antibiotics or volatile compounds, which suppress phytopathogens, but also the induction of systemic resistance to biotic pathogens and systemic tolerance to abiotic stress. Due to the various benefits they provide to host plants, endophytes have been drawing more and more attention in recent years.
Pantoea is a class of Gram-negative bacteria, which is ubiquitous in the environment. Many Pantoea strains are common endophytes [10,11]. The benefits of Pantoea to plants are similar to other endophytes as stated above. Intriguingly, some endophytic Pantoea strains isolated from sugarcane, sweet potato, rice, or other plants seem to share the trait of biological nitrogen-fixing activity [12,13,14,15]. For instance, P. agglomerans YS19, which has the ability to significantly promote the growth of rice, was verified for its nitrogen fixation function by acetylene reduction assay and 15N isotope labeling experiments. Thus, it is regarded suitable to be developed as a biological regulation agent [14]. Another endophyte Pantoea sp. 9C, isolated from sugarcane, can not only increases amino acid contents within sugarcane shoot tissues [16], but also significantly enhance the fresh weight of seedlings, the height of shoots, and both the length and number of rice seedlings’ roots [15]. This strain also has biological nitrogen fixation ability [13].
To date, research on endophytes mainly focuses on agricultural plants; only a few reports were performed on wood plants. Pantoea strain FBS135 is an endophytic bacterium isolated from the shoot of a healthy Pinus massoniana in Nanjing, China [17]. The greenhouse experiments evidenced that FBS135 can not only promote the growth of seedlings of P. massoniana, but also enhance the tolerance of pine seedlings under stress. Our earlier phylogenetic analysis revealed that FBS135 should be a member of P. eucalypti, a species firstly identified from eucalyptus in Uruguay [18]; however, the analysis was conducted solely based on its 16S rRNA sequences. Given the complexity of Pantoea genus [19] together with some different characters between FBS135 and P. eucalypti, we cannot exclude the possibility that FBS135 may be a novel species.
Nowadays, whole-genome information has become a requirement for exact taxonomic identification of prokaryotes. To clarify the taxonomic status of FBS135, in this study, we examined some important biological features of FBS135 and sequenced the whole genomes of the type strains for two species, which are pretty close to FBS135. Comparative genomic analysis was performed to reveal their genetic relationship.

2. Materials and Methods

2.1. Bacterial Strains and Media

The Pantoea strains and Azotobacter chroococcum ACCC10218 were routinely grown in LB medium unless otherwise stated. Ashby nitrogen-free medium (1 L): 10 g Mannitol, 0.2 g KH2PO4, 0.2 g MgSO4·7H2O, 0.2 g NaCl, 0.1 g CaSO4·2H2O, 5 g CaCO3, and 15 g Agar. NFDM nitrogen-free medium (1 L): 20 g sucrose, 0.015 g FeSO4·7H2O, 0.5 g MgSO4, 3.4 g KH2PO4, 0.01 g NaCl, 0.005 g Na2MoO4·2H2O, and 12.06 g K2HPO4 [20]. Gaseous Nitrogen-15N2(Enrichment > 99 atom%) was purchased from Shanghai Engineering Research Center of Stable Isotope.

2.2. Biological Assay

Physiological and biochemical tests were performed using GEN III MicroPlate following the instructions of the producer. FBS135 was grown on agar medium and then suspended in a special “gelling” inoculating fluid (IF) at the cell density of 90–98% T. Then the cell suspension was inoculated into the GEN III MicroPlate, 100 µL per well, and the MicroPlate was incubated at 30 °C to allow the phenotypic fingerprint to form [21,22]. The results were observed in 3–6 d.

2.3. Stable Isotope Labeling

Single colonies of FBS135 and A. chroococcum ACCC10218 were individually transferred to liquid medium and cultured at 30 °C and 200 rpm/min for four days. Then 50 μL of each bacterial culture were transferred to the plates of nitrogen-free media and spread evenly. The air containing 15N gaseous nitrogen (15N2) or normal air with 14N2 was used to incubate the plates. Three plates were placed in a sealed tank filled with the air containing 10% (v/v) 15N2 washed by sulfuric acid. Another three plates were placed in a tank filled with normal air. After both tanks had been incubated at 30 °C for 4 d, bacterial lawns were scraped off the plates and transferred into 50-mL centrifuge tubes. After being freeze-dried by an FD5-series freeze dryer (GOLD-SIM), the sample particle was passed through a clean 100-mesh sieve to obtain uniform bacteria powder, then packed in silver-foil cups and combusted within the reactor of an Organic Elemental Analyzer (FLASH 2000EA, Thermo Fisher Scientific, Bremen, Germany) using its existing reaction scheme for nitrogen analysis. Delta (15)N is determined using an isotope ratio mass spectrometer (DELTA V ADVANTAGE, Thermo Fisher Scientific, Bremen, Germany) [23].

2.4. Genome Sequencing

For DNA preparation, the bacteria were grown in LB at 37 °C with vigorous shaking. Cell pellets were collected and subjected to extraction with a Trizol® agent (Thermo Fisher Scientific, Burlington, NJ, USA). Electrophoresis showed that the DNA size of both genomes was greater than 10,000 bp without smear. The value of OD 260/OD280 ranged from 1.87 to 1.89. Careful examination indicated that the DNAs prepared were of good quality for sequencing. Sequencing was performed with Illumina Hiseq plus PacBio RS II GridION or PromethION sequencers at Nextomics Biosciences Co., Ltd. (Wuhan, China). After read filtering and adapter trimming, the reads sequenced by Illumina were de novo assembled using SOAP de novo (v2.04) (SciCrunch Registry, RRID:SCR_005400) before being corrected by PacBio sequencing result. All corrected reads were further assembled with Celera Assembler 8.0. Scaffolds were verified with Illumina data, and gaps were closed with GapCloser v1.12.
Sequences of the genomes were blasted (blastn, parameter setting: evalue = 0.00001, outfmt = 6, max_target_seqs = 100,000) against the plasmid proprietary database. When the length of total alignments accounted for > 1/5 of the total sequence and the sequence was less than 1 Mb, it was judged as a plasmid sequence.

2.5. Annotation

Annotations of protein coding, tRNA, and rRNA genes were conducted with the NCBI Prokaryotic Genome Annotation Pipeline [24]. The rpsBLAST (parameter: -evalue 0.01 -seg no -outfmt 5, ver. 89) was carried out for COG annotation [25,26]; GO annotations (ver. 20180828) were performed using Blast2go; functional annotations of encoding proteins were completed with three databases, Refseq (ver. 89), Pfam (ver. 31.0), and TIGRFAMs (ver. 15.0); protein data were retrieved from TIGRFAMs, Pfam, and GO databases, and treated by Interproscan (parameters: -appl Pfam, TIGRFAM, SMART -iprlookup -goterms -tp -f TSV, ver. 5.30–69.0); encoding proteins searched from KEGG (ver. 87.0-r20180701) and Ref-seq databases went through Blastp (parameters: -evalue 1 × 10−5- outfmt 6 std qlen slen stitle -max_target_seqs 5) and BlastX.

2.6. Phylogenetic Analysis

The rpoB, gyrB, infB, and atpD gene sequences of all strains we used were downloaded from GenBank for MLSA [27]. The maximum likelihood method and the proximity method were employed for analysis. After the four genes were aligned with Clustal W, phylogenetic profiling was conducted using MEGA software (Version 6.0).
Average nucleotide identity (ANI) and average amino acid identity (AAI) matrices were calculated by the EDGAR software framework. Blast hits between orthologous genes of the core of the genomes were analyzed for their mean/median percent identity values. The species cut-off was set to be 95% for the ANI and AAI indices [28]. JSpecies WS online (http://jspecies.ribohost.com/jspeciesws/, accessed on 19 February 2020) was used to calculate ANIb (ANI based on BLAST+) and ANIm (ANI based on MUMmer) values by pairwise genome comparisons. The tetra-nucleotide signatures (TETRA) correlation indexes were also determined using the JSpeciesWS software (version 3.7.11) [29].
The genome-to-genome distance calculator (v2.1) (http://ggdc.dsmz.de/ggdc.php#, accessed on 6 May 2021) was employed for the genome-based DNA-DNA hybridization analysis. Distances were calculated by alignment tools BLAST+ (recommended) and MUMMER. Three calculative formulas were used: Formula 1 (HSP length/total length), formula 2 (identities/HSP length), and formula 3 (identities/total length).
The phylogeny distance analysis was completed with Composition Vector Tree (CVTree, v3) (https://www.github.com/ghzuo/cvtree, accessed on 6 May 2021). The genome sequences input were in (.faa) format so as to exploit basic parameters. Built-in genomes were chosen for comparison, while eight genomes from the NCBI database were selected as outgroups. The output results were in (.nwk) format and visualization was done using the software FigTree.exe (ver. 1.4.4) and iTOL (ver. 5) [30].

2.7. Collinearity Analysis

Genome-wide collinearity analysis was performed by the online tool MAUVE (Multiple Alignment of Conserved Genomic Sequence, v20150226). Results were saved in (.xmfa) format for subsequent analysis [31]. The text editor Sublime Text3 was used to view the backbone files containing both start and end information of the common and differential conservative region sequences. The generated topological guide_tree files in (.newick) format were read by FigTree.

3. Results

3.1. FBS135 Has Prominent Low Nitrogen Tolerance and Different Biochemical Characters from the Type Strains of P. vagans and P. eucalypti

Unlike the morphology on nutrient-rich LB plates (Figure 1A), FBS135 grew very well on nitrogen-free Ashby medium plates showing smooth and translucent colonies (Figure 1B), a typical character of slime layer. This phenomenon implies that FBS135 may have the ability of biological nitrogen fixation. However, we were unable to detect acetylene reduction activity with FBS135 (data not shown). A series of measurements were taken including growing FBS135 on a different medium recipe (NFDM medium) and substituting agar with agarose in plate preparation so as to avoid possible nitrogen pollution; nevertheless, FBS135 still grew vigorously after seven successive acts of streaking.
To directly determine whether FBS135 had N2 fixation activity, we performed a tracing test with 15N isotope, and using A. chroococcus, a well-known nitrogen fixer, as the positive control. The results showed that the 15N/14N ratio detected for A. chroococcus incubated in the air with 15N2 was around 10 times higher than the same strain incubated in the normal air with 14N2. By contrast, the 15N/14N ratio detected for FBS135 incubated in the 15N air showed no significant difference with that grown in the 14N normal air, both of which were below 0.5% and approximately equal to the level for A. chroococcus incubated in the normal air. These results evidenced that FBS135 was unable to fix gaseous nitrogen.
In our previous work, we identified FBS135 as a P. eucalypti strain based on its 16S RNA sequence. In the meantime, a taxonomic classification was performed using the Biolog® system to learn more about its phenotypic patterns, including 71 carbon source utilization assays and 23 chemical sensitivity assays (Supplementary Materials Figure S1). The Biolog® identification result showed that biochemical and physiological characteristics of FBS135 was closest to P. vagans followed by P. eucalypti (Figure S1). This is somewhat different from the taxonomic result according to 16S RNAs, although it was reasonable given the close relationship between P. vagans and P. eucalypti.
FBS135, P. vagans LMG 24199T, and P. eucalypti LMG 24197T presented different capability of utilization of six carbon sources including Tween 40, α-D-lactose, gentiobiose, L-histidine, L-pyroglutamic acid, and also glucose (Table 1).
To further confirm this, an experiment on the carbon utilization of FBS135 was performed independently with glucose as the positive control (Figure S2F). On the sixth day, the optical density of FBS135 in gentiobiose culture reached its peak of approximately 5.0 (O.D.600 ≈ 5.0) (Figure S2C), followed by the culture fluid using Tween 40 or α-D-Lactose as the single carbon source (O.D.600 ≈ 1.0) (Figure S2A,B). On the contrary, FBS135 failed entirely to grow when L-pyroglutamic acid or histidine was used as the single carbon source (Figure S2D,E). Taken together, the individual assays confirmed the difference between FBS135 and the two type strains concerning their availability in the presence of the five carbon substrates. Although the differences may be only strain-specific, it is important not to exclude the possibility that FBS135 may be a novel species, since it is well known that distinguishing the species within Pantoea only by 16S sequences is challenging [19].

3.2. WGS Analysis of the Two Type Strains and FBS135

To precisely pinpoint the taxonomic status of FBS135, we did a whole-genome sequencing of two type strains, P. eucalypti LMG 24197T and P. vagans LMG 24199T, which are closely related to FBS135 according to 16S RNA profiling. The statistics showed that the qualities for both of the sequencing and the assembly were satisfactory (Table S1A,B). The genome of P. eucalypti LMG 24197T comprises a single circular chromosome (4,035,995 bp, C+G, 54.6%) and three circular plasmids, named pVag1 (~529 kb), pVag2 (138 kb), and pVag3 (94 kb) (Figure 2A). The genome sequence and the three plasmids are deposited in the DDBJ/EMBL/GenBank database under the accession numbers NZ_CP045720, NZ_CP045721, NZ_CP045722, and NZ_CP04573, respectively. The genome of LMG 24199T comprises a single circular chromosome (4,050,173 bp, C+G, 55.6%) and two circular plasmids, named pEuc1 (~560 kb) and pEuc2 (180 kb) (Figure 2B). The genome sequence and the three plasmids are deposited in the DDBJ/EMBL/GenBank database under the accession numbers NZ_CP038853, NZ_CP038854, and NZ_CP038855, respectively. Genome and plasmid information of the two strains have been summarized in Table 2. Different from the two strains, FBS135 has four plasmids, one (pPant4) of which seems to be responsible and indispensable for its growth [17].
Annotation analysis showed that the genome structures of P. eucalypti LMG 24197T and P. vagans LMG 24199T are generally similar yet with differences. P. eucalypti LMG 24197T contained a total of 4424 complete CDSs, while P. vagans LMG 24199T contained 4331 complete CDSs (Figure 2). P. vagans LMG 24199T contained more 5S rRNA sequence (Table S2) and two CRISPR (Clustered Regularly Interspaced Short Palindromic Repeat Sequences)-related genes in a total size of 1705 bp (Table S3). No CRISPR-related genes were found in the genomes of LMG 24197T or FBS135. In addition, P. vagans LMG 24199T contains one more genomic island (Table 3).
The genomes of P. eucalypti LMG 24197T and P. vagans LMG 24199T were annotated against COG, GO, KEGG, Refseq, Pfam, and TIGRFAMs databases (Table S4). Furthermore, 3144 genes (71.07%) of P. eucalypti LMG 24197T and 3180 genes (73.42%) of P. vagans LMG 24199T were successfully annotated with COG (Figure S3). Both genomes comprise more than 100 genes of unknown function. Moreover, 2739 genes (61.91%) of P. eucalypti LMGq 4197T and 2776 (64.1%) genes of P. vagans LMG 24199T were assigned to the GO database. The GO functional classification of the annotated proteins of the two genomes shared a high similarity. For both strains, the largest number of genes was assigned to the group of Biological Process, within which the sub-group of Metabolic Process (MP) and Cellular Process (CP) comprises the two highest numbers of genes (Figure S4). In addition, 2812 (63.56%) genes of P. eucalypti LMG 24197T and 2842 (65.62%) genes of P. vagans LMG 24199T were assigned to different KEGG pathways. For both strains, the pathway for membrane transports contained the largest number of their annotated genes (Figure S5). In general, the number of successfully annotated genes for P. eucalypti LMG 24197T and P. vagans LMG 24199T against different databases are similar.
Of note, in both genomes, we found neither nifH, the very conservative gene for nitrogenase, nor other common genes known to be involved in nitrogen fixation. Similarly, we did not find these genes in FBS135. These results confirm the inability of the three strains to fix inorganic nitrogen, which was just as consistent with the 15N-labelling experiment result above. Our results further showed that, similar to FBS135, both P. eucalypti LMG 24197T and P. vagans LMG 24199T could thrive on nitrogen-free Ashby plates after seven successive acts of streaking, indicating that the two strains also have strong ability of low nitrogen tolerance.

3.3. Phylogentic Analysis of the Two Type Strains and FBS135 Based on Genome Information

With the genome information above, we were able to make a more exact phylogenetic analysis for FBS135 together with P. eucalypti LMG 24197T and P. vagans LMG 24199T. We firstly constructed the phylogenetic trees based on the conversional multilocus sequence analysis (MLSA) in order to compare the result with an earlier study where P. eucalypti LMG 24197T and P. vagans LMG 24199T were initially reported [18]. Four conserved genes (rpoB, gyrB, infB, and atpD) of 21 type Pantoea strains were used for the MLSA while Budvicia aquatica DSM 5075T was used as the outgroup. According to the MLSA, FBS135, P. eucalypti LMG 24197T, and P. eucalypti LMG 24198 T were clustered into one branch (Table S5), which was basically consistent with the taxonomic result based on 16S rRNA sequences (Figure S6).
ANI, TETRA, and dDDH are three popular methods that are used currently to discriminate a wide range of bacterial taxons. An ANIb value between P. vagans LMG 24199T and P. eucalypti LMG 24197T was 90.09%, confirming that they were different species. The ANIb value between FBS135 and P. eucalypti LMG 24197T was greater than 99%, while the value between FBS135 with P. vagans LMG 24199T was merely 89.94% (Figure 3). ANIm values between the strains were almost the same as the ANIb values. The TETRA value between FBS135 and P. eucalypti LMG 24197T was 0.9997, greater than the corresponding cut-off value 0.998 (Table S6). With three different formulas, the digital DNA–DNA hybridization (dDDH) values of LMG 24197T and LMG 24199T were calculated to be 84.2%, 42.4%, and 75.9%, respectively. The dDDH between FBS135 and LMG 24199T were 95.3%, 94.6%, and 96.9% respectively, all exceeding the threshold value 70%. On the contrary, the dDDH between FBS135 and LMG 24199T were 69.6%, 42.3%, and 64.2%, respectively (Table S7), all less than 70%, indicating that FBS135 and LMG 24199T were different species. With all clues above, it was concluded that FBS135 and P. eucalypti LMG 24197T belonged to the same species.
It was interesting to know that from the perspective of GC mol% alone, the various Pantoea species showed obvious differences. The GC mol% differences between P. eucalypti LMG 24197T and P. vagans LMG 24199T is nearly 0.5%. The GC mol% differences between FBS135 and LMG 24197 is nearly 0.2%, while the difference between FBS135 and P. vagans LMG 24199T is around 1% (Table S7). Therefore, the GC mol% differences again corroborated that FBS135 was phylogenetically closer to LMG 24197T than to LMG 24199T.
Finally, we analyzed the taxonomic relationships of the strains with CVTree, a software that can construct phylogenetic trees based on whole-genome sequences but using a different algorithm from the above methods [32]. The result given by CVtree also showed that FBS135 and P. eucalypti LMG 24197T were clustered into one clade (Figure S7), which is consistent with the results from the other methods above.

3.4. Comparative Analysis Reveals Genomic Differences between FBS135 and Its Relative Strains

To better understand genomic differences among FBS135 and its phylogenetic relatives, their genomes were visualized together (Figure 4). It is easy to see that the genome of FBS135 has a higher similarity to the genome of P. eucalypti LMG 24197T than to the other Pantoea strains. Compared with FBS135, there are large fragment deletions at the region of 600 kb, 1000 kb, 1300 kb, and 3600 kb in these strains including P. eucalypti LMG 24197T and P. vagans LMG 24199T. FBS135 shared no gene overlap with the two type strains at 65 regions of different size ranging from 20 bp to 40,776 bp. Five of these regions were larger than 10,000 bp with the longest one (40,776 bp, from 1,092,468 bp to 1,133,244 bp) containing a multitude of genes encoding MFS transporter, transcriptional regulator, phage regulatory protein, and a lot of hypothetical proteins of unknown functions.
Five proteins (VgrG, TssH, TssK, B9D02_RS17585, B9D02_RS17595) related to the type VI secretion system were detected in a large region of 23,860 bp (3,673,699 bp~3,697,559 bp). In another large region of 37,425 bp (1,241,124 bp~1,278,359 bp), dozens of phage-related genes were found, indicating their possible phage origin. It is noteworthy that two genes encoding host specificity protein J (B9D02_RS06210) and serine acetyltransferase (B9D02_RS06115) were found in this region. Both of the genes were identified to be involved in FBS135 ability to tolerate low nitrogen (data now shown). The serine acetyltransferase gene (B9D02_RS06115) was also present in P. eucalypti LMG 24197T and P. vagans LMG 24199T, while the host-specific protein J gene (B9D02_RS06210) was only found in FBS135. Mechanisms of the two genes involved in low nitrogen tolerance remain unknown.
We aligned the genome sequences of the three strains according to their similarity and then split them into 1541 fragments of different sizes. FBS135 was found to have 65 exclusive fragments but lack 57 fragments present in the other two genomes; LMG 24197T had 93 exclusive fragments but lacked 44 fragments present in the other two genomes; P. vagans LMG 24199T had 368 exclusive fragments but lacked 445 fragments presented in the other two genomes. P. vagans LMG 24199T contained much more DNA fragments than P. eucalypti LMG 24197T and FBS135.
The collinearity analysis demonstrated that the genome structures of FBS135 and P. eucalypti LMG 24197T were almost identical. Their major fragments were well matched without obvious gene inversion or insertion (Figure 5A). The collinearity degree between FBS135 and P. vagans LMG 24199T was also high, although some inverted or rearranged fragments were detected. The genome collinearity between the FBS135 and other strains, for example, strain P. brenneri LMG5343T or P. conspicua LMG24534T, was poor. A large number of rearrangement events such as fragment inversion, deletion, translocation, or insertion were found in their genomes (Figure 5B). It is noteworthy that there was a large volume of differences between gene arrangements of P. brenneri LMG 5343T and FBS135, although they were phylogenetically close with a 16S rDNA similarity as high as 99.2%.
According to these collinearity analysis results, a rootless phylogenetic topology tree was generated (Figure 5C), where FBS135 and P. eucalypti LMG 24197T were clustered in the same clade while LMG 24199T was in another clade. Other type strains were even distant from FBS135. Such results were consistent with their phylogenetic relationships as analyzed above.

4. Discussion

In this study, the complete genomes of P. eucalypti LMG 24197T and P. vagans LMG 24199T have been sequenced and analyzed. With the genome information, we precisely determined the taxonomic status of Pantoea sp. FBS135, an endophyte isolated from Pinus massoniana with evidenced ability of promoting pine growth and stress tolerance [17]. We also demonstrated that FBS135 can tolerate an extremely low nitrogen environment without fixing nitrogen in the air. Comparative genomic analysis was performed for the two type strains and FBS135 revealed their close relationships and differences.
Compared with a pyramid of studies on symbiotic relationship between endophytic bacteria and agricultural plants, interactions between endophyte and forest plants are scarcely investigated. The significance of endophytes on forest plants may thus be underestimated, especially considering the huge influence of forest system on global ecology [33,34,35]. Studies on endophytes are also important for the production of economic trees. For example, endophytes have been demonstrated to have a remarkable performance as biocontrol agents and as plant growth promoters for grapevines [36,37]. In another example, P. agglomerans SWg2, an endophytic bacterium isolated from healthy mulberry roots, has been proven to not only promote mulberry growth and vitality, but also protect mulberry trees from bacterial pathogens [38].
We have previously reported the beneficial effect of FBS135 on pine growth [17]; however, the phylogenetic status of FBS135 was only roughly determined on the basis of 16S sequence, which will hinder further studies on it. In recent years, precise identification of a bacterium using whole-genome sequence has been proposed to be the standard [28,39]. Unfortunately, this has not been applied to FBS135 due to the absence of whole-genome information of type strains that are phylogenetically close to FBS135. It is challenging to distinguish species within Pantoea genera merely by 16S rRNA sequences. For instance, James T. Tambong found that many Pantoea strains have been previously misidentified by using MLSA. [40] The evolutionary relationship between P. eucalypti and P. vagans are also very close. In the GenBank, there are multiple genome sequences identified as P. vagans, which has been corrected to be P. eucalypti after a complete phylogenetic analysis. Therefore, it is extremely necessary to sequence the whole genomes of LMG 24197T and LMG 24199T, which do not only provide valuable information for the two type stains themselves, but also help to solve the taxonomic status of FBS135.
Widely distributed in the world and as one of the main tree species for artificial forestation in China, pine trees play an important role in the world’s forest system. Pine can tolerate barren soil and are often cultivated in barren soil. Ecological studies have shown that the nitrogen in the region is not enough to support the biomass of pine trees [41]. Thus, the nitrogen shortage has been a mystery for a long time. In recent years, some researchers have put forward convincing evidence that nitrogen-fixing endophytes in pine trees may be a main contributor responsible for the “mysterious” nitrogen source [33,34,42,43,44], although a deep understanding of the endophytes is still lacking. Mechanisms of biological nitrogen fixation in a highly oxygen-producing environment like pine needles also remains unclear. The possibility should not be excluded that some endophytes, during the long process of life evolution, may have formed special nitrogenases, which are not sensitive to oxygen.
As a class of typical endophytic bacteria, some Pantoea strains are found to have nitrogen fixation ability [12,13,14,15,45]. As the main objective of our study, Pantoea sp. FBS135 can grow vigorously on plates of two different nitrogen-free media (Ashby and NFDM), even after seven acts of successive streaking. Both mediums are prepared with analytically pure chemicals and the main carbon sources are different. Furthermore, we used agarose instead of agar to prepare plates in order to reduce trace nitrogen impurities. In either case, FBS135 grew very well. The odds are slim that FBS135 grew with nitrogen provided by the trace impurities from somewhere, which, therefore, caused us to speculate whether FBS135 has the ability of biological nitrogen fixation. However, the acetylene reduction assays with FBS135 showed a negative result. In addition, we failed to detect the presence of nifH, the most conserved gene for nitrogenase, in the genome of FBS135. Both facts appeared in contradiction to the speculation, indicating that either the speculation is incorrect, or FBS135 has a different nitrogen fixation mechanism from the traditional concepts.
In order not to miss a possible important discovery like novel nitrogenase, we exploited an additional method with stable isotope 15N labeling to determine whether FBS135 has the ability of N2 fixation. The detection with 15N isotope is extremely sensitive and, more importantly, it is a direct measurement, in contrast to the acetylene reduction assay, which may lead to a negative result in the case of a novel nitrogenase that cannot reduce acetylene but only N2. Unfortunately, the results proved that FBS135 cannot fix gaseous nitrogen. Its growth in/on strictly formulated nitrogen-free medium could be supported by extremely trace impurities in chemicals or agar/agarose.

5. Conclusions

In summary, in this study, we report the complete genome sequence of P. eucalypti LMG 24197T and P. vagans 24199T. As two type strains, their genome information is valuable for taxonomic identification of Pantoea strains, which are widely distributed in the environment but not identified at a species level [19,40]. In addition, the pine endophyte FBS135 was precisely identified with the genome information of the two strains. We also demonstrated that FBS135 failed to fix nitrogen but has a strong ability of tolerating extremely low nitrogen nutrients. In the future, more in-depth studies will be focused on the ability of FBS135′s tolerance to low nitrogen, which can help to elucidate the adaptability of this type of bacteria to the environment and shed light on the mechanisms of their symbiosis with host plants.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/life11070608/s1, Figure S1: ‘Phenotypic Fingerprint’ of FBS135 after six-day cultivation, Figure S2: Utilization of single carbon source by FBS135 in six days, Figure S3: Classification and statistics of COG function of genome -encoded proteins, Figure S4: Classification of GO functionally encoded proteins, Figure S5: Statistics of KEGG functional classification of genome-encoded proteins, Figure S6: Phylogenetic tree constructed based on MLSA with four conserved genes, Figure S7: Phylogenetic analysis based on the whole genome sequence of FBS135; Table S1A: Sequencing data statistics, Table S1B: Assembly result statistics, Table S2: ncRNA prediction result statistics, Table S3: CRISPR prediction result, Table S4: Functional annotation statistics of genome-encoded proteins, Table S5: Strain information for phylogenetic analysis, Table S6: Key parameters of phylogenetic analysis between FBS135 and P. eucalypti LMG 24197T, Table S7. Digital DNA-DNA hybridization results of FBS135 and two type strains.

Author Contributions

C.W., Z.S. and Y.L. conducted the main experiments and analyzed the data; B.F., Z.S. and C.W. wrote the manuscript draft; Y.Z. revised the manuscript; B.F. designed the experiments and finalized the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (No. 31970097) and the Key Scientific Project for Jiangsu Provincial Universities (17KJA220001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are greatly indebted to Teresa A. Coutinho, University of Pretoria, South Africa, for her kindly proving the two type strains P. eucalypti LMG 24197T and P. vagans 24199T.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Kado, C.I. Plant pathogenic bacteria. In Prokaryotes; Balows, A., Trüper, H.G., Dworkin, M., Harder, W., Schleifer, K.H., Eds.; Springer: New York, NY, USA, 1992; pp. 660–662. [Google Scholar]
  2. Elhady, A.; Adss, S.; Hallmann, J.; Heuer, H. Rhizosphere Microbiomes Modulated by Precrops Assisted Plants in Defense Against Plant-Parasitic Nematodes. Front. Microbiol. 2018, 9, 1133. [Google Scholar] [CrossRef] [Green Version]
  3. Hallmann, J.; Quadt-Hallmann, A.; Mahaffee, W.F.; Kloepper, J.W. Bacterial endophytes in agricultural crops. Can. J. Microbiol. 1997, 43, 895–914. [Google Scholar] [CrossRef]
  4. Partidamartínez, L.P.; Heil, M. The Microbe-Free Plant: Fact or Artifact? Front. Plant Sci. 2011, 2, 100. [Google Scholar]
  5. Santoyo, G.; Moreno-Hagelsieb, G.; Orozco-Mosqueda Mdel, C.; Glick, B.R. Plant growth-promoting bacterial endophytes. Microbiol. Res. 2016, 183, 92–99. [Google Scholar] [CrossRef] [PubMed]
  6. Hardoim, P.R.; van Overbeek, L.S.; Berg, G.; Pirttila, A.M.; Compant, S.; Campisano, A.; Doring, M.; Sessitsch, A. The Hidden World within Plants: Ecological and Evolutionary Considerations for Defining Functioning of Microbial Endophytes. Microbiol. Mol. Biol. Rev. 2015, 79, 293–320. [Google Scholar] [CrossRef] [Green Version]
  7. Glick, B.R. Beneficial Plant.-Bacterial Interactions; Springer International Publishing: Midtown Manhattan, NY, USA, 2020. [Google Scholar]
  8. Koskimäki, J.J.; Pirttilä, A.M.; Ihantola, E.L.; Halonen, O.; Frank, A.C. The intracellular Scots pine shoot symbiont Methylobacterium extorquens DSM13060 aggregates around the host nucleus and encodes eukaryote-like proteins. MBio 2015, 6, e00039-15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Sun, H.; Santalahti, M.; Pumpanen, J.; Köster, K.; Berninger, F.; Raffaello, T.; Asiegbu, F.O.; Heinonsalo, J. Bacterial community structure and function shift across a northern boreal forest fire chronosequence. Sci. Rep. 2016, 6, 32411. [Google Scholar] [CrossRef]
  10. Walterson, A.M.; Stavrinides, J. Pantoea: Insights into a highly versatile and diverse genus within the Enterobacteriaceae. FEMS Microbiol. Rev. 2015, 39, 968–984. [Google Scholar] [CrossRef] [Green Version]
  11. Smits, T.H.; Rezzonico, F.; Pelludat, C.; Goesmann, A.; Frey, J.E.; Duffy, B. Genomic and phenotypic characterization of a nonpigmented variant of Pantoea vagans biocontrol strain C9-1 lacking the 530-kb megaplasmid pPag3. FEMS Microbiol. Lett. 2010, 308, 48–54. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Asis, C.A., Jr.; Adachi, K. Isolation of endophytic diazotroph Pantoea agglomerans and nondiazotroph Enterobacter asburiae from sweetpotato stem in Japan. Lett. Appl. Microbiol. 2004, 38, 19–23. [Google Scholar] [CrossRef]
  13. Loiret, F.G.; Ortega, E.; Kleiner, D.; Ortega-Rodés, P.; Rodes, R.; Dong, Z. A putative new endophytic nitrogen-fixing bacterium Pantoea sp. from sugarcane. J. Appl. Microbiol. 2004, 97, 504–511. [Google Scholar] [CrossRef]
  14. Feng, Y.; Shen, D.; Song, W. Rice endophyte Pantoea agglomerans YS19 promotes host plant growth and affects allocations of host photosynthates. J. Appl. Microbiol. 2006, 100, 938–945. [Google Scholar] [CrossRef]
  15. Salazar-Gutierrez, I.; Loiret, F.G.; Ortega-Rodes, P.; Ortega, E. Nitrogen source affects glutamine synthetase activity in Pantoea sp.; bacterium inoculation promotes rice seedlings’ growth. Afr. J. Microbiol. Res. 2011, 5, 4774–4779. [Google Scholar] [CrossRef]
  16. Loiret, F.G.; Grimm, B.; Hajirezaei, M.R.; Kleiner, D.; Ortega, E. Inoculation of sugarcane with Pantoea sp. increases amino acid contents in shoot tissues; serine, alanine, glutamine and asparagine permit concomitantly ammonium excretion and nitrogenase activity of the bacterium. J. Plant. Physiol. 2009, 166, 1152–1161. [Google Scholar] [CrossRef]
  17. Song, Z.; Lu, Y.; Liu, X.; Wei, C.; Oladipo, A.; Fan, B. Evaluation of Pantoea eucalypti FBS135 for pine (Pinus massoniana) growth promotion and its genome analysis. J. Appl. Microbiol. 2020, 129, 958–970. [Google Scholar] [CrossRef]
  18. Brady, C.L.; Venter, S.N.; Cleenwerck, I.; Engelbeen, K.; Vancanneyt, M.; Swings, J.; Coutinho, T.A. Pantoea vagans sp. nov., Pantoea eucalypti sp. nov., Pantoea deleyi sp. nov. and Pantoea anthophila sp. nov. Int. J. Syst. Evol. Microbiol. 2009, 59, 2339–2345. [Google Scholar] [CrossRef] [PubMed]
  19. Janda, J.M.; Abbott, S.L. 16S rRNA gene sequencing for bacterial identification in the diagnostic laboratory: Pluses, perils, and pitfalls. J. Clin. Microbiol. 2007, 45, 2761–2764. [Google Scholar] [CrossRef] [Green Version]
  20. Kanvinde, L.; Sastry, G.R. Agrobacterium tumefaciens Is a Diazotrophic Bacterium. Appl. Environ. Microbiol. 1990, 56, 2087–2092. [Google Scholar] [CrossRef] [Green Version]
  21. Bochner, B.R. Sleuthing out bacterial identities. Nature 1989, 339, 157–158. [Google Scholar] [CrossRef]
  22. Bochner, B. Breathprints at the microbial level. ASM News 1989, 55, 536–539. [Google Scholar]
  23. Wang, X.; Feng, L.; Zhang, F.; Ding, Z. On-line measurements of δ15N in biological fluids by a modified continuous-flow elemental analyzer with an isotope-ratio mass spectrometer. Rapid Commun. Mass Spectrom. 2008, 22, 1196–1202. [Google Scholar]
  24. Pruitt, K.D.; Tatusova, T.; Brown, G.R.; Maglott, D.R. NCBI Reference Sequences (RefSeq): Current status, new features and genome annotation policy. Nucleic Acids Res. 2012, 40, D130–D135. [Google Scholar] [CrossRef] [Green Version]
  25. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tatusov, R.L.; Galperin, M.Y.; Natale, D.A.; Koonin, E.V. The COG database: A tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000, 28, 33–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Brady, C.; Cleenwerck, I.; Venter, S.; Vancanneyt, M.; Swings, J.; Coutinho, T. Phylogeny and identification of Pantoea species associated with plants, humans and the natural environment based on multilocus sequence analysis (MLSA). Syst. Appl. Microbiol. 2008, 31, 447–460. [Google Scholar] [CrossRef] [Green Version]
  28. Richter, M.; Rossello-Mora, R. Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 2009, 106, 19126–19131. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Richter, M.; Rosselló-Móra, R.; Oliver Glöckner, F.; Peplies, J. JSpeciesWS: A web server for prokaryotic species circumscription based on pairwise genome comparison. Bioinformatics 2015, 32, 929–931. [Google Scholar] [CrossRef]
  30. Letunic, I.; Bork, P. Interactive Tree of Life (iTOL): An online tool for phylogenetic tree display and annotation. Bioinformatics 2007, 23, 127–128. [Google Scholar] [CrossRef] [Green Version]
  31. Darling, A.C.E.; Mau, B.; Blattner, F.R.; Perna, N.T.; Darling, A.C.E.; Mau, B.; Blattner, F.R.; Perna, N.T. Mauve: Multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 2004, 14, 1394–1403. [Google Scholar] [CrossRef] [Green Version]
  32. Xu, Z.; Hao, B.L. CVTree update: A newly designed phylogenetic study platform using composition vectors and whole genomes. Nucleic Acids Res. 2009, 37, W174–W178. [Google Scholar] [CrossRef] [Green Version]
  33. Pirttilä, A.M.; Frank, C. Endophytes of Forest Trees: Biology and Applications, 2nd ed.; Springer International Publishing: Midtown Manhattan, NY, USA, 2018. [Google Scholar]
  34. Moyes, A.B.; Kueppers, L.M.; Pett-Ridge, J.; Carper, D.L.; Vandehey, N.; O’Neil, J.; Frank, A.C. Evidence for foliar endophytic nitrogen fixation in a widely distributed subalpine conifer. New Phytol. 2016, 210, 657–668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Carrell, A.A.; Frank, A.C. Bacterial endophyte communities in the foliage of coast redwood and giant sequoia. Front. Microbiol. 2015, 6, 1008. [Google Scholar] [CrossRef] [PubMed]
  36. Campisano, A.; Antonielli, L.; Pancher, M.; Yousaf, S.; Pindo, M.; Pertot, I. Bacterial endophytic communities in the grapevine depend on pest management. PLoS ONE 2014, 9, e112763. [Google Scholar] [CrossRef] [PubMed]
  37. Campisano, A.; Pancher, M.; Puopolo, G.; Puddu, A.; Lòpez-Fernàndez, S.; Biagini, B.; Yousaf, S.; Pertot, I. Diversity in Endophyte Populations Reveals Functional and Taxonomic Diversity between Wild and Domesticated Grapevines. Am. J. Enol. Vitic. 2015, 66, 12. [Google Scholar] [CrossRef]
  38. Xie, J.; Shu, P.; Strobel, G.; Chen, J.; Wei, J.H.; Xiang, Z.H.; Zhou, Z.Y. Pantoea agglomerans SWg2 colonizes mulberry tissues, promotes disease protection and seedling growth. Biol. Control. 2017, 113, 9–17. [Google Scholar] [CrossRef] [Green Version]
  39. Chun, J.; Oren, A.; Ventosa, A.; Christensen, H.; Arahal, D.R.; da Costa, M.S.; Rooney, A.P.; Yi, H.; Xu, X.W.; De Meyer, S.; et al. Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int. J. Syst. Evol. Microbiol. 2018, 68, 461–466. [Google Scholar] [CrossRef]
  40. Tambong, J.T. Taxogenomics and Systematics of the Genus Pantoea. Front. Microbiol. 2019, 10, 13. [Google Scholar] [CrossRef]
  41. Sheffer, E.; Batterman, S.A.; Levin, S.A.; Hedin, L.O. Biome-scale nitrogen fixation strategies selected by climatic constraints on nitrogen cycle. Nat. Plants 2015, 1, 15182. [Google Scholar] [CrossRef]
  42. Padda, K.P.; Puri, A.; Chanway, C. Endophytic nitrogen fixation—A possible ‘hidden’ source of nitrogen for lodgepole pine trees growing at unreclaimed gravel mining sites. FEMS Microbiol. Ecol. 2019, 95, fiz172. [Google Scholar] [CrossRef]
  43. Kandel, S.L.; Herschberger, N.; Kim, S.; Doty, S. Diazotrophic Endophytes of Poplar and Willow for Growth Promotion of Rice Plants in Nitrogen-Limited Conditions. Crop. Sci. 2015, 55, 1765–1772. [Google Scholar] [CrossRef] [Green Version]
  44. Knoth, J.L.; Kim, S.H.; Ettl, G.J.; Doty, S.L. Biological nitrogen fixation and biomass accumulation within poplar clones as a result of inoculations with diazotrophic endophyte consortia. New Phytol. 2014, 201, 599–609. [Google Scholar] [CrossRef] [PubMed]
  45. Ruppel, S.; Merbach, W. Effect of ammonium and nitrate on 15N2-fixation of Azospirillum spp. and agglomerans in association with wheat plants. Microbiol Res. 1997, 152, 377–383. [Google Scholar] [CrossRef]
Figure 1. Low nitrogen tolerance ability of FBS135. (A) Morphology of FBS135 on Luria-Bertani plate. (B) Morphology of FBS135 on Ashby medium. (C) Detection of 15N content in bacterial cells of P. eucalypti FBS135 and A. chroococcum ACCC10218 incubated under different conditions. The two strains were cultured with normal air or the air containing 15N labeled N2.
Figure 1. Low nitrogen tolerance ability of FBS135. (A) Morphology of FBS135 on Luria-Bertani plate. (B) Morphology of FBS135 on Ashby medium. (C) Detection of 15N content in bacterial cells of P. eucalypti FBS135 and A. chroococcum ACCC10218 incubated under different conditions. The two strains were cultured with normal air or the air containing 15N labeled N2.
Life 11 00608 g001
Figure 2. Genome structure diagram for Pantoea eucalypti LMG 24197T and Pantoea vagans LMG 24199T. The circle diagrams from the outermost layer to the innermost layer show their sense strands, the negative sense strands, tRNAs, rRNAs, CRISPR-related genes, gene islands, GC content percentage, GC-skew, and genome sequencing depth. (A) Main genome structure diagram of P. eucalypti LMG 24197T. (B) Main genome structure diagram of P. vagans LMG 24199T m.
Figure 2. Genome structure diagram for Pantoea eucalypti LMG 24197T and Pantoea vagans LMG 24199T. The circle diagrams from the outermost layer to the innermost layer show their sense strands, the negative sense strands, tRNAs, rRNAs, CRISPR-related genes, gene islands, GC content percentage, GC-skew, and genome sequencing depth. (A) Main genome structure diagram of P. eucalypti LMG 24197T. (B) Main genome structure diagram of P. vagans LMG 24199T m.
Life 11 00608 g002
Figure 3. Heat map of ANI comparison between FBS135 and related strains.
Figure 3. Heat map of ANI comparison between FBS135 and related strains.
Life 11 00608 g003
Figure 4. Genome-wide visualization results of FBS135 and some related strains.
Figure 4. Genome-wide visualization results of FBS135 and some related strains.
Life 11 00608 g004
Figure 5. Comparative analysis revealed genomic difference between FBS135 and its relative strains. Collinearity analysis results of FBS135 with P. eucalypti LMG24197T and P. vagans LMG24199T (A) and FBS135 with P. brenneri LMG5343T and P. conspicua LMG24534T (B). Simulated topology tree for FBS135 and the related strains according to their collinear analysis (C).
Figure 5. Comparative analysis revealed genomic difference between FBS135 and its relative strains. Collinearity analysis results of FBS135 with P. eucalypti LMG24197T and P. vagans LMG24199T (A) and FBS135 with P. brenneri LMG5343T and P. conspicua LMG24534T (B). Simulated topology tree for FBS135 and the related strains according to their collinear analysis (C).
Life 11 00608 g005
Table 1. Carbon utilization difference between FBS135 and two related strains.
Table 1. Carbon utilization difference between FBS135 and two related strains.
CharacteristicP. vagans LMG 24199TP. eucalypti LMG 24197TFBS135
Tween 40++(+)
α-D-Lactose-+(+)
Gentiobiose(+)-+
L-Histidine++-
L-pyroglutamic acid-+-
Note: +, 90–100% of strains positive in 3–6 days; (+), weakly positive; -, 90–100% of strains negative.
Table 2. Genomic and plasmid information of P. eucalypti LMG 24197T and P. vagans LMG 24199T and FBS135.
Table 2. Genomic and plasmid information of P. eucalypti LMG 24197T and P. vagans LMG 24199T and FBS135.
StrainTypeLength (bp)GC Content (%)Topology
Pantoea vagans LMG24199Tgenome4,050,17355.6circular
pVag1559,69254.0circular
pVag2180,46453.1circular
Pantoea eucalypti LMG 24197Tgenome4,035,99554.6circular
pEuc1529,30352.3circular
pEuc2138,72550.9circular
pEuc394,96752.2circular
P. eucalypti FBS135genome4,023,75156.2circular
pPant1535,55352.4circular
pPant2199,37250.7circular
pPant3133,12751.7circular
pPant473,56253.3circular
Table 3. Gene structure statistics of P. eucalypti LMG 24197T and P. vagans LMG 24199T.
Table 3. Gene structure statistics of P. eucalypti LMG 24197T and P. vagans LMG 24199T.
StrainP. vagans LMG 24199TP. eucalypti LMG24197T
TypeNumberLength (bp)%
Genome
NumberLength (bp)% Genome
tRNA8062680.138062920.13
16S rRNA710,7000.22710,7020.22
23S rRNA720,3370.42720,3550.42
5S rRNA89200.0278050.02
CDS43314,147,02086.5744244,127,71586.01
CRISPR217050.04000
genomic_
island
250,8471.06334,3730.72
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wei, C.; Song, Z.; Lu, Y.; Zhao, Y.; Fan, B. Relationship of the Pine Growth Promoting Pantoea eucalypti FBS135 with Type Strains P. eucalypti LMG 24197T and P. vagans 24199T. Life 2021, 11, 608. https://doi.org/10.3390/life11070608

AMA Style

Wei C, Song Z, Lu Y, Zhao Y, Fan B. Relationship of the Pine Growth Promoting Pantoea eucalypti FBS135 with Type Strains P. eucalypti LMG 24197T and P. vagans 24199T. Life. 2021; 11(7):608. https://doi.org/10.3390/life11070608

Chicago/Turabian Style

Wei, Chunyue, Zhongwen Song, Yiming Lu, Yinjuan Zhao, and Ben Fan. 2021. "Relationship of the Pine Growth Promoting Pantoea eucalypti FBS135 with Type Strains P. eucalypti LMG 24197T and P. vagans 24199T" Life 11, no. 7: 608. https://doi.org/10.3390/life11070608

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