Shuai Xu
Yanxia Shi
Baoju Li
Lei Li- Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China
Many Rahnella strains have been widely described as plant growth-promoting rhizobacteria with the potential to benefit plant growth and protect plants from pathogens. R. aceris ZF458 is a beneficial plant bacterium isolated from swamp soil with the potential for biocontrol. Strain ZF458 has shown broad-spectrum antagonistic activities against a variety of plant pathogens and exhibited a dramatic effect on controlling Agrobacterium tumefaciens in sunflowers. The R. aceris ZF458 genome sequence contained a 4,861,340-bp circular chromosome and two plasmids, with an average G + C content of 52.20%. Phylogenetic analysis demonstrated that R. aceris ZF458 was closely related to R. aceris SAP-19. Genome annotation and comparative genomics identified the conservation and specificity of large numbers of genes associated with nitrogen fixation, plant growth hormone production, organic acid biosynthesis and pyrroloquinoline quinone production that specific to benefiting plants in strain ZF458. In addition, numerous conserved genes associated with environmental adaption, including the bacterial secretion system, selenium metabolism, two-component system, flagella biosynthesis, chemotaxis, and acid resistance, were also identified in the ZF458 genome. Overall, this was the first study to systematically analyze the genes linked with plant growth promotion and environmental adaption in R. aceris. The aim of this study was to derive genomic information that would provide an in-depth insight of the mechanisms of plant growth-promoting rhizobacteria, and could be further exploited to improve the application of R. aceris ZF458 in the agriculture field.
Introduction
Rahnella was a Gram-negative, rod-shaped, facultatively anaerobic bacterium belonging to the family Yersiniaceae (Adeolu et al., 2016), and was widely distributed in a variety of environments, including soil, phyllosphere, water, seeds, food, some clinical samples, and even in American mastodon remains (Berge et al., 1991; Rhodes et al., 1998). Rahnella was first proposed in 1979, and normally consisted of six species including R. aquatilis (which was described as the type species of Rahnella), R. variigena, R. inusitata, R. bruchi, R. woolbedingensis, and R. victoriana by multilocus sequence analysis (Brady et al., 2014). Recently, a new species named R. aceris was proposed in 2020, which shared a high similarity with the R. aquatilis (Lee et al., 2020). To date, only one R. aceris strain SAP-19 has been reported for general physiological and biochemical characteristics, shotgun genome information and phylogenetic analyses (Lee et al., 2020). Since there have been very few studies on R. aceris, its genetic background, taxonomic classification and potential biological function were still unclear.
Decades of research have proven that Rahnella spp. could be applied to preventing and controlling plant disease or promoting plant growth for their beneficial properties. Many Rahnella strains have been described to protect plants from a wide range of phytopathogenic organisms in different ways. R. aquatilis JZ-GX1 secreted volatile organic compounds (VOCs) that destroyed the mycelial growth of Colletotrichum gloeosporioides, and thus restrained the infection and expansion of anthracnose disease in leaves (Kong et al., 2020a). R. aquatilis Ra36 could exploit fungal chemotropism to efficiently colonize the roots of host plants, resulting in efficient protection from vascular wilt disease in tomatoes by interfering with the alkaline-triggered infection of Fusarium oxysporum via gluconic acid secretion (Palmieri et al., 2020). R. aquatilis HX2 was proven to decrease the disease incidence of crown gall disease caused by Agrobacterium vitis in grapevines by producing antibacterial substances (Chen et al., 2009). Besides, R. aquatilis Ra39 exhibited a significant effect on controlling Erwinia amylovora in apple, and the control efficiency was close to the efficacy of streptomycin when the strain was mixed with the Acibenzolar-S-methyl (ASM) (Abo-Elyousr et al., 2010). Application of R. aquatilis Ra to bean plants increased the content of phenolic compounds and the activity of peroxidase (PO) enzyme than untreated plants, and resulted in a marked disease suppression against Xanthomonas axonopodis pv. phaseoli under greenhouse and field conditions (Sallam, 2011). R. aquatilis 17 and 55 exhibited a marked siderophore production and a significant effect on controlling the bacterial spot of cucumber (el-Hendawy et al., 2003). In addition, many Rahnella strains have been demonstrated to benefit plant growth through various mechanisms, such as nitrogen fixation, phytohormone production, phosphate solubilization, biosynthesis of organic acids and pyrroloquinoline quinone. For example, R. aquatilis ZF7 significantly promoted the weight of aboveground parts and the root length in Chinese cabbage due to its high Indole-3-acetic acid (IAA) biosynthetic capacity (Yuan et al., 2020). R. aquatilis HX2 exhibited a significant promotion of the growth of corn, due to its ability to solubilize mineral phosphate, produce pyrroloquinoline quinine and IAA (Guo et al., 2009). R. aquatilis JZ-GX1 exhibited a significant greening effect on Cinnamomum camphora by producing organic acids (Kong et al., 2020b), and was also reported as a phytate-degrading rhizobacteria (PDRB) that could improve the growth of poplar and Masson pine (Li et al., 2013). To date, the majority of studies on plant growth promotion in Rahnella spp. focused on R. aquatilis strains. Given the close relationship between R. aceris and R. aquatilis, the R. aceris strains may possess similar mechanisms of biocontrol and plant growth promotion.
Previous studies have demonstrated that Rahnella strains were widely distributed and adapted to diverse ecological environments, which might be related to their resistance to acids, salts, selenium, antibiotics and heavy metals. For example, R. aquatilis HX2 exhibited a tolerance to high salt, strong acids, antibiotics, and stress tolerance (Li et al., 2019, 2021). In addition, R. aquatilis HX2 could grow in high selenium concentrations and could reduce selenate and selenite to BioSeNPs (Zhu et al., 2018). R. aquatilis JZ-GX1 could stimulate the production of exopolysaccharides, and resulting in high salt tolerance to 0–9% NaCl (Li et al., 2021). R. aquatilis ZF7 harbored 24 genes associated with the biosynthesis and metabolism of β-lactamase, and exhibited tolerance to a variety of antibiotics including ampicillin, carbenicillin and vancomycin (Yuan et al., 2020). Rahnella sp. SMO-1 was reported to encode extended-spectrum β-lactamase which could efficiently hydrolyze penicillins and cefotaxime, and the strain was also showed resistance to amoxicillin, ticarcillin, cefalotin and cefuroxime (Lartigue et al., 2013). In addition, Rahnella sp. JN6 showed a high resistance to Cd, Pb, and Zn and could improve the efficiency of remediation in heavy metal contaminated soils (He et al., 2013). Previous studies have indicated that a range of regulatory mechanisms including secretory system, two-component regulatory system (TCS) and acid-resistance genes played important roles in the environmental adaption of bacteria (Liu and Ochman, 2007). For instance, multiple secretion systems were found in Yersinia enterocolitica and were essential for virulence, the biosynthesis of antibiotics and the life activities (Heermann and Fuchs, 2008). Photorhabdus luminescens harbored 18 TCSs which had been found to be involved in metabolite utilization and adaptation to various stress factors (Heermann and Fuchs, 2008). Nevertheless, relatively few studies have focused on adapting to environmental variations in Rahnella strains. Therefore, to better promote the application prospects of Rahnella spp., studying the environmental adaption of Rahnella strains, especially at the molecular level was very necessary.
In this study, R. aceris ZF458 with broad spectrum of antagonistic activities against many plant pathogens was isolated from swamp soil. Phylogenetic analysis, ANI and DDH analysis were accomplished to definite the taxonomic position of ZF458 and the relationship with R. aceris SAP-19 and other R. aquatilis strains. The whole genome of ZF458 was sequenced, annotated and compared with the genomes of other typical and widely reported Rahnella strains HX2, ZF7, Y9602 and ATCC 33071T. Comparative genome analysis was used to reveal various genes involved in biocontrol factors and environmental adaption. These data would provide an in-depth insight into the mechanisms of plant growth promotion, and improve the bio-control application of R. aceris ZF458 in the future.
Materials and Methods
Strain Isolation, Antagonistic Assays and Biocontrol Assays
Stain ZF458 was isolated from swamp soil in Sichuan, China, in August, 2019, according to a standard 10-fold dilution plating assay as described by previous study (Wang B. et al., 2013). The antagonistic activities of R. aceris ZF458 against plant pathogens were performed by plate tests, and the control effect of ZF458 against Agrobacterium tumefaciens was tested on sunflowers. The growth curve of strain ZF458 was measured according to the following steps: strain ZF458 was incubated overnight at 28°C, and the suspension was diluted to an OD600 = 0.6 with sterile Nutrient Broth (NB) medium. Subsequently, 1 μL of the suspension was added to 999 μL NB culture via vortexing. Then, 200 μL of the diluent was dripped to a 96 well plate. The microtiter plates were kept in a stable shaking speed (180 rpm⋅min–1), and the OD600 of ZF458 culture was measured at 4 h intervals over 40 h. In addition, plate assays were conducted to assess siderophore and phosphatase activities based on reported approaches (Milagres et al., 1999; Luis et al., 2014). If siderophores were produced, the chrome azurol S (CAS) agar would change from blue to orange, while if phosphatases were produced, the Pikovskaya’s (PVK) agar plate would change from white to transparent.
Growth Conditions, Microscopic Analysis and Genomic DNA Extraction
Strain ZF458 was incubated in NB medium at 28°C with shaking for 24 h. The morphology of strain ZF458 was observed by scanning electron microscope (SEM) Hitachi-S3400N and transmission electron microscope (TEM) Hitachi-7700. The genomic DNA of strain ZF458 was extracted from cultured cells (OD600 = 0.8) by TIANamp Bacteria DNA kit (Tiangen Biotech (Beijing) Co., Ltd.).
Genome Sequencing and Annotation
The whole genome of R. aceris ZF458 was sequenced in Allwegene Technologies Corporation, China. The pacific Biosciences (PacBio) RS II platform was used for the whole-genome sequencing of strain ZF458, and a 10-kb SMRT Bell was used for template library construction. The sequence reads were assembled de novo using SMRT Link v.5.1.0.1 CGView was used to generate the circular genome visualization (Stothard and Wishart, 2005). GeneMarkS (version 4.17) software was used to predict the coding genes of the sequenced genome.2 Transfer RNA (tRNA), ribosomal RNA (rRNA) and small nuclear RNAs (snRNAs) were predicted using tRNAscan-SE version 1.3.1 (Lowe and Eddy, 1997), rRNAmmer version 1.2 (Lagesen et al., 2007) and cmsearch version 1.1 (Nawrocki et al., 2009), respectively. Functional gene annotation was carried out through multiple general databases, including Non-Redundant protein databases (NR) (Li et al., 2002), Gene Ontology database (GO) (Ashburner et al., 2000), Kyoto Encyclopedia of Genes and Genomes (KEGG) (Kanehisa et al., 2006), Cluster of Orthologous Groups of proteins (COG) (Galperin et al., 2015), Transporter Classification Database (TCDB) (Saier et al., 2009), Pfam,3 Swiss-Prot4 (Amos and Rolf, 2000), Carbohydrate-Active enZTmes Database (CAZy) (Cantarel et al., 2009). Furthermore, signal peptides and transmembrane structure were predicted using SignalP 4.1 and TMHMM 2.0c (Petersen et al., 2011). In addition, prophages were predicted by using phiSpy 2.3 (You et al., 2011), and CRISPR (Clustered Regularly Interspaced Short Palindromic Repeat Sequences) was identified via CRISPRdigger 1.0 (Ibtissem et al., 2007).
Phylogenetic Analysis and Comparative Genomic Analysis
Multilocus gene sequence analysis (MLSA) based on the four housekeeping genes (16S rRNA, gyrB, atpD, and rpoB) was used to determine the taxonomic position of R. aceris ZF458. The housekeeping gene sequences of diverse strains (Supplementary Table 1) were aligned using MUSCLE and the phylogenetic tree was generated using the maximum likelihood method by MEGA 6.0 (Tamura et al., 2013). Average nucleotide identities (ANI) and in silico DNA-DNA hybridization (DDH) analysis were performed using the OrthoANIu algorithm and the Genome-to-Genome Distance Calculator (GGDC), respectively (Meier-Kolthoff et al., 2013). Based on the phylogenetic analysis, four closely related Rahnella strains including R. aquatilis ZF7 (CP032296.1), R. aquatilis HX2 (CP003403.1), Rahnella sp. Y9602 (CP002505.1) and R. aquatilis ATCC 33071T (CP003244.1) were selected for comparative genomic analysis. Moreover, pairwise alignment of genomes was generated using the Mauve 2.3.1 comparison software, and Venn diagrams were conducted using R package (Richter and Rossello-Mora, 2009).
Functional Gene Analysis Linked With Plant Growth Promotion and Environmental Adaption
Functional genes beneficial to plants such as IAA production, nitrogen fixation, phosphate solubilization, organic acid and pyrroloquinoline quinone (PQQ) biosynthesis, were searched for in the NCBI database. In addition, functional genes involved in environmental adaption including bacteria secretion system, flagella biosynthesis and chemotaxis, two-component system, quorum sensing, SeNPs production, acid-resistance were detected using the NCBI database. All the functional genes retrieved from the ZF458 genome were compared at the amino acid level with four closely related Rahnella strains ZF7, HX2, Y9602 and ATCC 33071T by KEGG database.
Results
Antagonistic Activity and Biological Characteristics of Strain ZF458
Biocontrol Traits of Strain ZF458
Strain ZF458 which was isolated from swamp soil exhibited broad, obvious antagonistic activities against diverse plant-pathogenic fungi and bacteria, including Stemphylium solani, Corynespora cassiicola, Ascochyta citrulline, Colletotrichum sp., Phytophthora capsici, Fusarium oxysporum, Pectobacterium brasiliense, Pseudomonas amygdali pv. lachrymans, Pseudomonas syringae pv. tomato, Ralstonia solanacearum, Xanthomonas campestris pv. campestris, Clavibacter michiganensis subsp. michiganensis, Agrobacterium tumefaciens, and Acidovorax citrulli (Figure 1). Importantly, strain ZF458 presented a dramatic effect in controlling Agrobacterium tumefaciens on sunflowers (Supplementary Figure 1), with a control efficiency of 80.20%. Moreover, orange halos were appeared around colonies of strain ZF458 on CAS agar plates, revealing that strain ZF458 produced siderophores (Supplementary Figure 2A). Similarly, transparent halos were observed around ZF458 colonies on PVK agar plates, indicating that strain ZF458 produced phosphatases (Supplementary Figure 2B).
Figure 1. Antagonistic assays of R. aceris ZF458 against various plant fungal and bacterial pathogens. (A,a) Botrytis cinerea, (B,b) Stemphylium solani, (C,c) Rhizoctonia solani, (D,d) Corynespora cassiicola, (E,e) Ascochyta citrulline, (F,f) Colletotrichum sp., (G,g) Fusarium oxysporum, (H,h) Phytophthora capsici, (I,i) Pectobacterium brasiliense, (J,j) Pseudomonas amygdali pv. lachrymans, (K,k) Pseudomonas syringae pv. tomato, (L,l) Ralstonia solanacearum, (M,m) Xanthomonas campestris pv. campestris, (N,n) Clavibacter michiganensis subsp. michiganensis, (O,o) Agrobacterium tumefaciens, (P,p) Acidovorax citrulli. The capital letters represent plate containing strain ZF458, the small letters represent control.
Organism Information
Strain ZF458 was determined to be a Gram-negative, facultatively anaerobic bacterium belonging to the Yersiniaceae family. The strain produced white to light yellow-colored, circular, translucent, convex with margin colonies and reached 0.1–5.0 mm in a diameter on 1-day culture on NA plate at 28°C (Supplementary Figure 3A). Strain ZF458 displayed shot rod-shaped cells with lengths of 1.5–3.0 μm and diameters of 0.7–1.3 μm (Supplementary Figures 3B,C). The growth curve showed that strain ZF458 was in exponential growth phase between 4–20 h after inoculation, and attained the stationary phase at 20 h of incubation (Supplementary Figure 4A). Furthermore, the pH of strain ZF458 culture declined rapidly from approximately 7.0 to 3.2 within 12 h of incubation, and remained stable after 12 h (Supplementary Figure 4B).
Genome Structure and Genome Comparison Between R. aceris ZF458 and Other Completely Sequenced Rahnella Strains
General Genomic Features of R. aceris ZF458
The whole genome of R. aceris ZF458 was 5.60 Mb with an overall G + C content of 52.20% (Supplementary Table 2). The entire genome comprised a circular 4,861,340-bp chromosome (CP067057.1) and two plasmids with sizes of 188,256-bp (CP067058.1) and 553,387-bp (CP067059.1). The general genome structure and functions of strain ZF458 was represented by graphical circular genome map (Figure 2). In total 5,248 predicted genes were identified in the genome, including 4,988 protein coding genes, 109 RNA genes and 88 pseudogenes. According to KEGG database, a total of 4,954 genes of ZF458 were annotated to 40 different pathways which were associated with metabolism, cellular processes, environmental information processing, genetic information processing and organismal systems. By comparing the Swiss-Prot, Pfam, TCDB, SinalIP and CAZy databases, 3,316 (63.19%), 3,533 (67.32%), 939 (17.89%), 428 (8.16%), and 210 (4.00%) of the ORFs were annotated into different groups, respectively. Besides, the ZF458 genome contained 342 secreted proteins, 10 genomic islands and 7 prophage regions (Supplementary Table 2).
Figure 2. Graphical circular maps of the R. aceris ZF458 chromosome and plasmids pPlas1 and pPlas2 genome generated using the CGview server. From outside to center, ring 1 showed the genome sequence position, rings 2 and 5 showed protein-coding genes oriented in the forward (colored by COG categories) and reverse (colored by COG categories) directions, respectively. Rings 3 and 4 showed coding genes in the forward (blue) and reverse (green) directions, respectively. Ring 6 showed the G + C% content and the inner most ring showed the GC skews, where blue indicated positive values and yellow indicated negative values.
Moreover, a total of 4,237 functional genes of ZF458 were further annotated using COG database. Most of the genes were related to amino acid transport and metabolism (10.57%), followed by carbohydrate transport and metabolism (9.86%), transcription (8.68%), inorganic ion transport and metabolism (6.66%), cell wall/membrane biogenesis (6.21%) and translation, ribosomal structure and biogenesis (5.17%) (Table 1). However, 83 genes were not annotated by COG database, whose specific functions need to be further verified.
Table 1. Number of genes associated with general COG functional categories.
Phylogenetic Tree
To understand the genetic relationships of R. aceris ZF458 with other related strains especially the Rahnella strains, a phylogenetic tree was established based on 16S rRNA and three other housekeeping genes (gyrB, atpD and rpoB). As expected, seventeen Rahnella strains were clustered into one major clade, and three Rouxiella strains, two Yersinia strains, two Serratia strains were in other clades (Supplementary Figure 5). The phylogenetic analysis indicated that 17 Rahnella strains including ZF458 were in seven Rahnella MLSA groups. Based on the observed genetic distance relationships, strain ZF458 was closely clustered together with R. aceris SAP-19, R. aquatilis ZF7, R. aquatilis HX2, Rahnella sp. Y9602 and the type strain R. aquatilis ATCC 33071T, successively followed by R. victoriana strains, R. variigena strains, R. bruchi strains, R. woolbedingensis strains, and R. inusitata strains.
Average Nucleotide Identities and DNA-DNA Hybridization Analysis
Average nucleotide identities and DDH were widely used to certify the genetic distance of bacteria at the genomic level, and strains revealing ANI values ≥96% and DDH values ≥70% were typically considered as the same species (Zhang et al., 2016). In the study, ANI and DDH calculations among Rahnella strains were performed. The results showed that the ANI values between ZF458 and SAP-19, ZF7, HX2, Y9602, KM05, KM12, MEM40 were both >98%, and the DDH values among them were both >80% (Supplementary Figure 6). These findings demonstrated that the eight R. aceris strains or R. aquatilis strains (ZF458, SAP-19, ZF7, HX2, Y9602, KM05, KM12, and MEM40) were closely clustered with one another and occupied a close genetic distance relationship. However, the ANI and DDH values between strains ZF458 and ATCC 33071T were 92.75 and 49.4% respectively, although the two strains were clustered closely in the phylogenetic tree (Supplementary Figure 5). Smaller ANI and DDH values were calculated when other Rahnella strains (except the R. aceris and R. aquatilis species) were used as reference genomes.
Genomic Features of Different Rahnella Strains
In the study, four Rahnella strains that have completed whole genome sequencing including R. aquatilis ZF7, R. aquatilis HX2, Rahnella sp. Y9602, and R. aquatilis ATCC 33071T were selected for comparative genomic analysis (Table 2). In comparison, the entire genome size of the five Rahnella strains ranged from 5.45 to 5.66 Mb, the G + C content ranged from 52.08 to 52.36%, and the predicted coding genes ranged from 4,804 to 5,065. The results illustrated that the genome size of strain ZF458 was larger than that of ZF7 and ATCC 33071T, but smaller than that of HX2. Furthermore, the genomes of strains ZF458 and Y9602 contained one circular chromosome and two plasmids, while HX2 and ATCC 33071T contained one circular chromosome and three plasmids, and the sequenced strain ZF7 harbored a circular chromosome and only one plasmid (Table 2).
Table 2. Genomic features of Rahnella aceris ZF458 and other Rahnella spp.
Comparison of R. aceris ZF458 With the Four Typical Rahnella Strains
In addition, the entire genome of ZF458, SAP-19, ZF7, HX2, Y9602 and ATCC 33071T were compared using Mauve to assess the homology among these Rahnella strains. The results showed that more local collinear blocks (LCB) inversion was obviously observed between strains ZF458 and ATCC 33071T than that among strains ZF458, ZF7, HX2 and Y9602 (Figure 3A). This indicated that the ZF458 genome was highly syntenic with R. aquatilis ZF7, R. aquatilis HX2 and Rahnella sp. Y9602. At the species level, less local collinear blocks inversion was found between strains ZF458 and SAP-19, but compared with the ZF458 genome, large local collinear blocks (LCB) translocation occurred in the genome of SAP-19 (Figure 3B).
Figure 3. Global alignment of the genome sequence of completely sequenced R. aceris ZF458 against five other Rahnella genome sequences. (A) Mauve progressive alignment of the ZF458, ZF7, HX2, Y9602, and ATCC 33071 genomes. (B) Mauve progressive alignment of the ZF458 and SAP-19 genomes. ZF458 genome was used as the reference genome. Boxes with the same color indicate syntenic regions. Boxes below the horizontal strain line indicate inverted regions. Rearrangements are indicated with the colored lines. The scale is in nucleotides. (C) Venn diagram showing the numbers of shared and unique Clusters of Orthologous Genes among ZF458, ZF7, HX2, Y9602, and ATCC 33071 genomes. (D) Venn diagram showing the numbers of shared and unique Clusters of Orthologous Genes between the ZF458 and SAP-19 genomes.
As shown in Figure 3C, 3,530 conserved genes were shared by strains ZF458, ZF7, HX2, Y9602 and ATCC 33071T. ZF458 shared 3,893 genes with ZF7, 3,884 genes with HX2, and 3,928 genes with counterparts in the Y9602 genome, while only 3,647 orthologous genes existed in strains ZF458 and ATCC 33071. Interestingly, no unique genes were found in the genome of R. aquatilis strains ZF7, HX2, Y9602 and ATCC 33071T, while one unique gene was found in ZF458 genome (Figure 3C). In addition, 3,795 orthologous genes were found between the ZF458 and SAP-19 genomes, and 356 unique genes existed in the ZF458 genome were not retrieved in the genome of SAP-19 (Figure 3D).
Synergistic Interactions With the Host Plant: Benefiting Plants
Nitrogen Fixation and Metabolism
To compare the key genes associated with the plant growth promotion in R. aceris ZF458, four strains including ZF7, HX2, Y9602 and ATCC 33071T were selected. The results showed that a relatively whole nif gene cluster with 21 genes related to nitrogen fixation was retrieved in the genomes of strains ZF458, HX2 and ATCC 33071T (Table 3). All the genes were assigned to Mo nitrogenase structural genes, FeMo-co synthesis, maturation of nitrogenase, electron transport, regulation of nif gene expression and three unknown functional genes, and shared high similarities among the three strains. However, only one gene nifJ encoding pyruvate ferredoxin/flavodoxin oxidoreductase was present in the ZF7 genome, and indeed, the nif gene cluster which existed in strains ZF458, HX2 and ATCC 33071T was found in the plasmid genome, except the gene nifJ, which was retrieved from the chromosome genome. In addition, 7 highly homologous nitrogen regulation related genes (glnABDGO, ntrB, and nac) were found in the genomes of R. aceris ZF458 and four other Rahnella strains (Table 3).
Table 3. Homolog analysis of nitrogen fixing gene cluster in R. aceris ZF458 and other Rahnella strains.
Production of the Plant Growth Hormone Indole-3-Acetic Acid
The comparative genome analysis indicated that eleven key genes responsible for IAA production were found in the ZF458 genome, which contained eight genes related to tryptophan operon (trpG, trpE, trpD, trpCF, trpB, trpA, trpS, and trpR), mtr gene encoding tryptophan permease, ipdC gene encoding indolepyruvate decarboxylase, and acdS gene encoding 1-aminocyclopropane-1-carboxylate (ACC) deaminase. Moreover, all of these genes were found in the genomes of ZF458, ZF7, HX2, Y9602 and ATCC 33071T with sequence identities exceeding 93% at the amino acid level (Supplementary Table 3).
Phosphate Solubilization
According to the comparative genomic analysis, 26 genes involved in phosphate solubilization were found in the ZF458 genome. Among these genes, six genes (aphA, phoN, iap, phoA, phoB, and phoR) were related to organic phosphate acquisition, and twenty genes were associated with carbon-phosphorus lyase, which contained the whole carbon-phosphorus lyase operon (phnNMLKJIHGFAEDC) and phosphate transporter (pstSCAB) (Supplementary Table 4). All these genes retrieved from the ZF458 genome shared high similarities (greater than 94%) with those of four other Rahnella strains ZF7, HX2, Y9602 and ATCC 33071, except the gene aphA, which was absent in the ATCC 33071 genome (Supplementary Table 4).
Biosynthesis of Organic Acid
Furthermore, 44 genes linked with organic acid biosynthesis were found in the genome of strains ZF458, ZF7, HX2, and Y9602, while three genes eda, idn, and ldh were absent in the genome of ATCC 33071 (Table 4). The 44 genes were assigned for diverse pathways, including 3 genes for Glycolytic pathway (EMP), 24 genes for Tricarboxylic Acid Cycle pathway (TCA), 4 genes for Entner-Doudoroff pathway (ED), 3 genes for direct oxidation pathway, and 10 genes for another undefined pathway. The identities of these genes between ZF458 and ZF7, HX2, Y9602 were higher than those between ZF458 and ATCC 33071T strains, especially for the gene gntK, which only showed 52% sequence identities between ZF458 and ATCC 33071T (Table 4). The results demonstrated that these strains may play important roles in inorganic phosphorus dissolution.
Table 4. Homolog analysis of organic acid biosynthesis genes in R. aceris ZF458 and other Rahnella strains.
Production of Plant Promotion Hormone Pyrroloquinoline Quinone
According to the comparative analysis, the R. aceris ZF458 plasmid genome contained a highly conserved pyrroloquinoline quinone gene cluster (pqqABCDEF) and upstream promoter orfX, covering 6,631 bp with 7 ORFs (Supplementary Figure 7), and the gcd gene encoding PQQ-dependent dehydrogenase. The pqq gene cluster of ZF458 shared a high homology with those in other four Rahnella strains (>95% identity), except that the sequence identity of pqqF between ZF458 and ATCC 33071T was only 80% (Supplementary Table 5).
Survival and Adapting to Environmental Variations
Bacterial Secretion System
Type II Secretion System
In order to explore the different secretion systems, five released complete genome sequences of Rahnella strains namely ZF458, ZF7, HX2, Y9602, ATCC 33071T were compared. The genome of ZF458 contained a variety of secretion systems, which were closely associated with the life activities of bacteria. The comparative analysis showed that the R. aceris ZF458 chromosome contained a highly conserved T2SS (type II secretion system) gene cluster (gspACDEFGHIJKLM and pilD) (Figure 4D), covering 12,797 bp with 13 ORFs. All 13 genes existed in strains ZF458, HX2, Y9602 and ATCC 33071T, while only four genes (gspE, gspL, gspM, and pilD) were found in ZF7 genome. The gsp gene cluster from the ZF458 genome shared a high consistency with strains HX2 and Y9602 (exceeding 97%) at the amino acid level, whereas lower sequence identities for genes gspA, gspC, gspI, gspJ, gspL, and gspM (under 80%) were found between ZF458 and ATCC 33071T (Supplementary Table 6).
Figure 4. Comparison of the bacterial secretion system gene clusters of R. aceris ZF458 against four previously fully sequenced Rahnella genomes. (A) Type I secretion system. (B) Type III secretion system. (C) Type IV secretion system. (D) Type II secretion system. (E) Sec (secretion) system. (F) Conjugal transfer region of type IV secretion. (G) Type VI secretion system. The same color represented genes with the same or similar biological function. Arrows denoted putative transcriptional units. The length of blocks represented the size of genes (1 cm = 1000 bp).
Type IV Secretion System
In addition, nine key genes related to T4SS (virB2,3-4,5,6,8,9,10,11) were found in strains ZF7, HX2, and ATCC 33071T, but not in ZF458 and Y9602 (Figure 4C), and all the genes were retrieved from the chromosome genome. However, another gene mltE associated with T4SS was found in the five Rahnella strains and shared a 99% similarity at the amino acid level (Supplementary Table 6). Furthermore, the tra and trb gene clusters with 19 ORFs associated with the conjugal transfer region of the type IV secretion system were discovered in the plasmid genomes of ZF458 and HX2 (Figure 4F). Most of these genes were highly conserved with sequence identities exceeding 90%, except that traG and traI with sequence identities 87 and 89%, respectively (Supplementary Table 6). Interestingly, none of these genes existed in the genome of strains ZF7, Y9602 and ATCC 33071T.
Type VI Secretion System
In this study, 27 genes related to T6SS were found in the genome of R. aceris ZF458, and 13 genes were regarded as core genes (Figure 4G). The 13 core T6SS genes were highly conserved in the five Rahnella strains with a consistency of amino acid sequences over 90%, although the genes were in different encoding orders. However, 14 core T6SS genes existed in the genome of strains ZF7, HX2, Y9602, and ATCC 33071T, for the core gene tssL absent in the ZF458 genome (Supplementary Table 6). Moreover, other function genes involved in inner membrane proteins (ImpL and ImpK), extracellular structural components (vgrG and hcp) and regulatory proteins or structural proteins (VasL and ImpF) were present in the ZF458 genome, and shared a high similarity with those of the other 4 Rahnella strains at the amino acid level (Supplementary Table 6).
Other Secretion Systems
In addition, three genes related to T1SS (raxA, raxB, and tolC) (Figure 4A), three genes associated with T3SS (hrpA, hrpB, and radD) (Figure 4B), ten genes associated with Sec system (Figure 4E), four genes connected with twin arginine targeting (tatB, tatC, and tatE) and two genes related to SRP components (ffh and ftsY) were all present in the genomes of ZF458, ZF7, HX2, Y9602 and ATCC 33071T, most of which were highly conserved (Supplementary Table 6).
Two-Component System
Moreover, 19 TCSs (two-component system) were found in the genomes of strains ZF458, ZF7, HX2, Y9602 and ATCC 33071T, most of which were highly conserved (Supplementary Table 7). According to the homology, the topological characteristics of sensor histidine kinase (HK) and response regulator (RR) (Lavín et al., 2007), the 19 TCSs were divided into diverse subfamilies. For insurance, the motility TCS (CheA/Y) existed in ZF458 with a high identity (more than 98%) at the amino acid level among of the five Rahnella strains. The nitrogen regulation TCS (glnL/G) was searched for in ZF458 and shared a high homology with that in ZF7, HX2, Y9602 and ATCC 33071T (>99% identity). Phosphate regulation TCS (phoR/B) was observed in the genomes of strains ZF458, ZF7, HX2, Y9602 and ATCC 33071T, with sequence identities exceeding 98% (Supplementary Table 7). All these genes might play important roles in R. aceris ZF458 in terms of physiological metabolism and environmental adaptation.
Selenium Metabolism
In addition, many genes associated with selenium metabolism were found in the genome of R. aceris ZF458, which were similar to strains ZF7, HX2, Y9602 and ATCC 33071T (identities >90%) (Table 5). Among these genes, 8 genes including cysA, cysW, cysT, cysP, cysJ, cysI, tsgA, and hfq were associated with selenium transportation and regulation, 19 genes containing nar gene cluster and tat gene cluster were associated with selenate reduction, 15 genes were related to selenite reduction, and 8 genes were involved in selenium association (Table 5).
Table 5. Genes related to selenium metabolism in R. aceris ZF458 and other Rahnella strains.
Acid Resistance
Rahnella strains adapted to diverse ecological environments, including water, soil and plants, which may be attributed to their stress tolerance traits. In this study, 28 genes associated with acid resistance were present in ZF458, and were grouped into different systems according to the enzyme activity. All the genes shared high similarities (greater than 90%) with those of other Rahnella strains, except for the gene ompF (encoding porin), which showed a lower homology with ZF7, HX2, Y9602 and ATCC 33071T (identity <90%) (Supplementary Table 8).
Motility, Chemotaxis and Quorum Sensing
A group of tightly gene clusters (flh, flg, and fli) related to flagella biosynthesis encoding 41 flagella-associated proteins (FlhDC, MotAB, FlhBAE, FlgA-N, FliC-T, FliA, and FliZ) were found in the genome of R. aceris ZF458, and had high amino acid similarities among the five Rahnella strains, except for the gene fliC, and fliD, which displayed low similarities between ZF458 and ATCC 33071T (49% identities) (Supplementary Table 9). In addition, 8 genes (cheA, cheW, mcp, tar, cheR, cheB, cheY, and cheZ) associated with chemotaxis were found in the genome of ZF458, with sequence similarities exceeding 97% between ZF458 and other four Rahnella strains (Supplementary Table 9). 12 key genes related to the autoinducer-2 (AI-2)- dependent signaling systems were detected in the five Rahnella strains, with sequences identities exceeding 94% between ZF458 and four other strains ZF7, HX2, Y9602 and ATCC 33071 (Supplementary Table 10).
Discussion
Rahnella was widely distributed in a variety of environments, and had been reported as PGPR that benefited plant growth (Sallam, 2011; Li et al., 2014). However, their specific mechanisms of plant-growth promotion and environmental adaption at the molecular level especially for R. aceris were still unclear. In this study, R. aceris ZF458 which was isolated from swamp soil was displayed. Strain ZF458 produced siderophores and phosphatases, which were key factors to promote plant growth (Li et al., 2013). ZF458 presented broad spectrum antagonistic activities against various plant pathogens and exhibited a significant effect in controlling Agrobacterium tumefaciens on sunflowers. All of these characteristics indicated that ZF458 might be a plant growth-promoting rhizobacteria. The complete genome of R. aceris ZF458 was sequenced and compared with other Rahnella strains to better understand the molecular mechanisms for plant-growth promotion and environmental adaption.
The phylogenetic trees showed that ZF458 clustered together with SAP-19, and the two strains both belonged to the R. aceris species, forming a tight cluster with R. aquatilis, which was consistent with a previous study (Lee et al., 2020). A previous study also indicated that R. aceris SAP-19 shared an ANI value of 92.7% and a DDH value under 48.6% with the genome sequence of R. aquatilis (Lee et al., 2020), however, our study showed that higher ANI and DDH values (>96and >70%, respectively) were calculated between R. aceris ZF458 and R. aquatilis strains. The results also indicated that strain ZF458 was closely related to the R. aquatilis strains, although it was classified as a new species of Rahnella genus by the phylogenetic analysis. According to the phylogenetic tree, strain ATCC 33071T was in a farther clade. Interestingly, the analysis of ANI and DDH revealed the same result, as the ANI and DDH values between ZF458 and ATCC 33071T were 92.75 and 49.4%. The results also indicated that the ANI and DDH values between ATCC 33071T and other R. aquatilis strains (including HX2, ZF7, Y9602) were below 95 and 70%, respectively, and this was consistent with the previous report showing that the ANI value between R. aquatilis strains ZF7 and ATCC 33071T was only 92.9% (Yuan et al., 2020). Previous studies indicated that strain ZF7, HX2 and Y9602 were both isolated from soil (Guo et al., 2012; Martinez et al., 2012a), while ATCC 33071T was isolated from drinking water (Martinez et al., 2012b), so, it may suggest that the genetic evolutionary distance could be associated with the habitat and adaptation. All five Rahnella strains harbored one or more plasmids, which could be related to genetic evolution, based on the previous report that 19% of strains in the genus Rahnella were plasmid-containing and had highly homologous regions for the same species (Rozhon et al., 2010).
In the study, an obvious gene transfer between ZF458 and other R. aquatilis strains indicated that strain ZF458 was another species not belonging to R. aquatilis. Furthermore, the nucleotide level similarities between ZF458 and ZF7, HX2, Y9602 were markedly higher than the similarity between ZF458 and ATCC 33071T, demonstrating the far phylogenetic distance between ZF458 and ATCC 33071T and supporting the phylogenetic result that ZF458 and ATCC 33071T were in different subclades. Low similarity regions were found between ZF458 and SAP-19, although the two strains both belonged to the R. aceris species, for the reason of strain SAP-19 was just whole genome shotgun sequence, not completed genome sequence. Based on the comparative analysis, only one unique gene existed in the ZF458 genome, while no unique genes existed in ZF7, HX2, Y9602 and ATCC 33071T, which confirmed the phylogenetic analysis described above that ZF458 belonged to R. aceris but was closely clustered with R. aquatilis.
Nitrogen was an important element for plant growth since plants only absorbed reduced forms of nitrogen, such as ammonia and nitrates. Biological nitrogen fixation (BNF), a microbiological process which converted atmospheric nitrogen into a plant usable form, offered an alternative to N fertilizers (Brewin and Legocki, 1992). This capability presented an opportunity to improve the nitrogen source utilization and crop yields, through the introduction of nitrogen fixing bacteria into crops, or the nitrogenase enzyme responsible for nitrogen fixation (Oldroyd and Dixon, 2014). Recently, a limited number of archaea and bacteria were found to fix nitrogen, such as Proteobacteria, Firmicutes, Cyanobacteria, and Actinobacteria (Santos et al., 2012). Rahnella spp. were also reported to possess the capacity for nitrogen fixation (Berge et al., 1991). For, insurance, R. aquatilis HX2 was reported as a PGPR and could promote the growth of corn for the ability to fix nitrogen (Guo et al., 2012). Previous studies indicated that the nifHDK gene cluster encoding Mo-nitrogenase played important role in fixing nitrogen, and nifBENXV was important for the synthesis and maturation of the FeMo cofactor (Wang L. et al., 2013). In this study, 21 nif genes were detected in the genome of ZF458, which predicted that ZF458 had a strong nitrogen fixation ability.
Auxin IAA was a primary plant hormone synthesized by plant growth-promoting bacteria that had a profound influence on benefiting plant (Bal et al., 2013). Currently, Rahnella spp. isolated from different environments were reported to have a strong capacity to produce IAA. Such as R. aquatilis ZF7 exhibited a high IAA biosynthesis capacity of 30.86 μg mL–1 and had a significant growth promoting effect on Chinese cabbage (Yuan et al., 2020). R. aquatilis HX2 could biosynthesize IAA and significantly promote maize growth (Guo et al., 2012). The endophyte Rahnella strains isolated from sweet potato plants, showed the ability to produce IAA and revealed obvious growth promotion in hybrid poplar (Khan and Doty, 2009). In this study, 11 genes associated with IAA biosynthesis were found in the ZF458 genome, including 9 genes (trp gene cluster) linked with tryptophan, ipdC gene encoding indole-3-pyruvate decarboxylase, and acdS gene encoding ACC deaminase, which revealed high similarities to R. aquatilis strains HX2 and ZF7. Previous study showed that acdS gene played a major role in PGPR activities. For example, the disruption of acdS gene in R. aquatilis HX2 reduced IAA production and decreased the growth promotion activity on corn (Peng et al., 2019). Previous studies demonstrated that the IPyA pathway was the main IAA synthetic pathway in a broad range of bacteria (Eastman et al., 2014), and the ipdC gene played an important role in regulating the IPyA pathway (Spaepen et al., 2007). So, the IPyA pathway may play an important role in the biosynthesis of IAA for R. aceris ZF458.
Phosphorus was considered as an important macronutrient for the growth and development of plants. Microorganisms were reported to play an important role in dissolving the insoluble phosphates in the rhizosphere (Passariello et al., 2006). Currently, phosphate solubilization of microorganisms was usually divided into two aspects, inorganic phosphorus solubilization and organic phosphorus solubilization, and organic phosphate in soil was degraded mainly through enzymes such as phosphatase, phytase and carbon-phosphorus lyase (White and Metcalf, 2007). Previous study showed that the carbon-phosphorus lyase system which encoded by phn gene cluster played important roles in organic phosphate dissolution (Jochimsen et al., 2011). In this study, the entire phn gene cluster was observed in the ZF458 genome, which was highly conserved in Rahnella strains, thus indicating that strain ZF458 may possess the strong ability to dissolve organic phosphorus. In addition, mineral phosphate solubilization could be improved by the production of organic acids, including gluconic, formic, citric, acetic, lactic, and acetic (Vassilev et al., 2006). According to the comparative genomic analysis, 44 genes related to organic acid biosynthesis were retrieved in the ZF458 genome, which were highly similar to other Rahnella strains. Thus, the organic acids generated by these genes played an important role in inorganic phosphorus solubilization for R. aceris ZF458. Many Rahnella strains have been reported to have a strong phosphate-solubilizing activity or organic acid synthesis ability. For example, R. aquatilis JZ-GX1 was able to secrete organic acids in an iron-deficient environment, which facilitated the production of phosphatases and siderophores, and thus promoted plant growth (Kong et al., 2020b). Rahnella sp. W25 isolated from the crop rhizosphere of calcareous soil, showed the maximum phosphate-solubilizing capability on tricalcium phosphate, aluminum phosphate, and ferric phosphate of 385.5, 110.4, and 216.6 mg⋅L–1 respectively (Qiao et al., 2013). In addition, the phosphorous solubilized by Rahnella strains W25 was significantly negatively correlated with the pH of the culture medium (Qiao et al., 2013). The aphA gene encoded a molecular class B bacterial phosphatase was reported to exhibit a better activity at acidic pH values (Passariello et al., 2006). Importantly, the pH value of the liquid culture of R. aceris ZF458 reduced to 3.24 from the original value of 7.0 with the growth of culture time, which might indicate that the AphA (protein id, WP_134705824.1) enzyme generated strong activity when solubilizing organic phosphate.
Gluconic acid was reported as a crucial antibacterial substance produced by R. aquatilis, and pyrroloquinoline quinone (PQQ) was necessary for the production of gluconic acid (Choi et al., 2008). Previous studies had demonstrated that the GDH (glucose dehydrogenase)-PQQ holoenzyme was important for the biosynthesis of an antimicrobial substance (Matsushita et al., 2002; Werra et al., 2009), and PQQ was reported as a plant growth promotion factor, which in addition to gluconic acid production, has also been associated with its antioxidant properties (Choi et al., 2008). Bacterial genes associated with PQQ biosynthesis have been identified in various species were clustered in pqqABCDEF operons (Guo et al., 2009). Previous studies revealed that pqqA gene encoded a small peptide that contained tyrosine and glutamate and served as the precursor and rate-determining step for PQQ production (Goosen et al., 1992). Also, pqqB was reported as a carrier that facilitated the secretion of PQQ across the plasma-membrane into the periplasm in K. pneumoniae (Velterop et al., 1995). For example, the disruption of the pqqA or pqqB gene in R. aquatilis HX2 reduced the plant growth promotion activity, and eliminated the ability to produce antibacterial substances thus decreasing the control effect against grapevine crown gall (Li et al., 2014). In this study, the whole pqq operons were observed in the ZF458 genome, and shared a high similarity with those in HX2, which indicated that R. aceris ZF458 could produce PQQ normally, although accurate biochemical functions remain to be confirmed.
PGPR exerted a series of regulatory mechanisms to adapt to the various complicated environments, among them, the bacterial secretion system, two-component system were well studied regulatory mechanisms associated with environmental adaption. Previous studies had indicated that large numbers of bacterial genomes (up to 17% in Proteobacteria genomes) encoded proteins in the general secretory pathway (GSP) (Bendtsen et al., 2005). Bacteria exported different kinds of substrates, including small molecules, extracellular enzymes, effector proteins and DNA via several sophisticated secretion systems, which were important for bacteria to adapt and survive the diverse ecological environments (Gerlach and Hensel, 2007).
Type II secretion system depended on the sec machinery and secreted folded proteins including various hydrolyzing enzymes such as the pseudolysin from the periplasm across the extracellular environment, which was important for bacterial survival in host or environment (Bendtsen et al., 2005). Many T2SS components had been well characterized, such as pectinases and cellulases associated with pathogenicity were secreted by T2SS in Pectobacterium (Filloux, 2004). T2SS was well-conserved and composed of 12–15 components, which were called Gsp proteins (Cianciotto, 2005). GspD formed a multimeric pore for transportation of secreted proteins and interacted with GspC, which played an important role in V. cholerae (Korotkov et al., 2011). The GspC, GspF, GspL, and GspM constituted the inner membrane platform and were associated with ATPase secretion in the cytoplasm (Johnson et al., 2006). In this study, ZF458 harbored 12 gsp genes which were highly conserved in the Rahnella strains, and the Gsps might be associated with the protein secretion and transportation, but the role of T2SS in R. aceris ZF458 needs to be determined in the future.
T4SS was characterized by the ability to regulate the translocation of single-stranded DNA or complex proteins, and existed in both Gram-negative and Gram-positive bacteria (Alvarez-Martinez and Christie, 2009). The T-DNA transfer system of A. tumefaciens was the prototypical type A T4SS, and together with the T4SS in E. coli was the best-characterized T4SS, which consisted of 12 proteins including VirB1-VirB11 and VirD4 (Christie et al., 2014). T4SS was recognized as the most ubiquitous secretion system in nature, for the ability to conjugate which was a common bacteria trait. However, no vir genes was found in the ZF458 genome, except for the conjugal transfer region (including tra and trb gene clusters) which existed in the plasmid genome. Previous studies revealed that T4SS evolved from bacterial conjugation machinery according to the sequence similarities (Christie and Cascales, 2003). Such as the octopine-type Ti-plasmids which were recognized as the conjugal transfer system in Agrobacterium tumefaciens was closely related to Ti plasmid vir genes (Alt-Mörbe et al., 1996). In addition, the tra genes related to conjugal transfer system were reported to be involved in the synthesis and possible retraction of the sex pilus (Winans and Walker, 1985), while the vir genes were also proven to be required for pilus biogenesis (Trokter et al., 2014). Thus, the gene clusters in the conjugal transfer region of ZF458 might possess the same function as T4SS, which needs to be further verified.
T6SS was reported to transfer effector proteins into eukaryotic and prokaryotic cells and had a crucial role in bacterial competition (Ho et al., 2014). T6SS was composed of 13 conserved core genes and few accessory genes which were necessary for function though the genes varied in composition among species (Chang et al., 2014). ImpL, an integral polytopic inner membrane protein, interacted with the essential T6SS component, and ImpK was proved to be necessary for its function of mediating the Hcp secretion in A. tumefaciens (Ma et al., 2009). TssL, TssM, and TssJ composed a membrane complex in the T6SS, and played an important role in protein secretion and the assembly of cell surface appendages in E. coli (Felisberto-Rodrigues et al., 2011). Notably, the comparative genomic analysis indicated the variation in the gene composition between Rahnella strains, but contained the 13 conserved core genes, which might reveal that Rahnella perform the normal function of T6SS. To date, T6SS had been proved to be associated with the symbiotic interactions with hosts in many bacteria, but more studies need to explore the function of T6SS in the Rahnella strains.
Two-component regulatory system was an essential mechanism for signal transduction, generally consisting of a sensor histidine kinase (HK) and a response regulator (RR) (Lavín et al., 2007). Previous studies showed that GacS/GacA system controlled the biosynthesis of secondary metabolites and extracellular enzymes that were linked with the ecological fitness and tolerance to stress (Heeb and Haas, 2001). PhoB/R TCS played a pivotal role in regulating the inorganic phosphate regulatory system (Yamada et al., 1989). The CheA/CheW TCS regulated chemotaxis and RcsC/RcsB were involved in the motility of bacteria (Majdalani and Gottesman, 2005). Overall, TCS played an essential role in the life activity of bacteria, for the diverse functional TCSs that could adapt to various environments. Recently, it was reported that TCS had a close relationship with biocontrol factors in bacteria, however, the functions of TCSs in ZF458 need to be elucidated.
Selenium (Se) was an important trace element for humans and animals and existed in four states with chemical forms of selenide, elemental selenium, selenite and selenate. Microorganism was reported to play a key role in the biogeochemical cycle of Se in the ecological environment (Skalickova et al., 2017). Many bacteria have been reported to reduce selenate and/or selenite to elemental Se nanoparticles (SeNPs) (Stolz et al., 2006). SeNPs could be used for medicine, therapeutics and environmental remediation due to a variety of advantages including biocompatibility and low toxicity (Wang et al., 2007). Recently, the SeNPs synthesized by bacteria had drawn extensive attention due to their great potential for the bioremediation of polluted environments. R. aquatilis HX2 was reported to show a strong tolerance to Se and reduce selenate and selenite to BioSeNPs (Zhu et al., 2018). R. aquatilis ZF7 could produce SeNPs from Na2SeO3, and possessed the most production of SeNPs at 0.85 mM (Yuan et al., 2020). Previous studies proved that the small RNA chaperone Hfq was a key regulator for the reduction of selenium nanoparticles in R. aquatilis HX2 (Xu et al., 2020). In this study, many genes linked with selenium metabolism were found in the ZF458 genome, and shared a high homology with those in the genomes of R. aquatilis strains ZF7 and HX2, and the results indicated that R. aceris ZF458 might synthesize SeNPs and had a widely applications in nano-agriculture.
Acid resistance genes were essential for successful colonization in acidic environments however, most studies on acid resistance genes had been conducted in E. coli. For instance, the alternative stress sigma factor (RpoS), was reported to be associated with the oxidative AR system, which allowed E. coli O157:H7 to survive in acidic environments (Masuda and Church, 2010). In addition, the glutamate and arginine decarboxylase systems in E. coli were also revealed as part of its acid resistance mechanism (Kim et al., 2016). Based on the comparative genomic analysis, many acid resistance genes were observed in the ZF458 genome, but their function still needs to be confirmed.
Conclusion
Rahnella aceris ZF458, which was isolated from the swamp soil, exhibited broad-spectrum antagonistic activities against 14 diverse plant pathogens, and displayed a significant effect in controlling Agrobacterium tumefaciens on sunflowers. A whole genomic analysis and a phylogenetic analysis showed that strain ZF458 belonged to R. aceris, and was closely related to R. aquatilis. A comparative genomic analysis indicated that R. aceris ZF458 harbored many genes related to nitrogen fixation, IAA production, organic phosphate dissolution, organic acid biosynthesis and PQQ production, which had been proved to be beneficial to plant growth. In addition, large numbers of genes associated with environmental adaption, such as the secretion system, two-component system, selenium metabolism and acid-resistance were found in the ZF458 genome and shared a high homology with Rahnella strains ZF7, HX2, Y9602 and ATCC 33071T. As far as we know, this was the first study to systematically analyze the genes associated with plant growth promotion and environmental adaption in R. aceris. Overall, these features of R. aceris ZF458 provided new insight into the comprehensive metabolic pathways (Figure 5), and would assist in studying plant growth promotion and biocontrol application in R. aceris.
Figure 5. Metabolic pathway map of R. aceris ZF458 related to plant growth promotion and environmental adaption.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA796123.
Author Contributions
SX, LL, and BL conceived and designed the experiments. SX, YZ, and YP performed the experiments and analyzed the data. SX and LL wrote the manuscript. LL, XX, YS, and AC revised the manuscript. All authors have read and approved of the final version of the manuscript.
Funding
This work was performed with the support of the China Agriculture Research System of MOF and MARA (CARS-25); the Science and Technology Innovation Program if the Chinese Academy of Agricultural Sciences (CAAS-ASTIP-IVFCAAS); and the Key Laboratory of Horticultural Crops Genetic Improvement, Ministry of Agriculture in China (IVF2020).
Conflict 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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to acknowledge Yang Cao (Allwegene Technologies Corporation) for the genome data analysis. Finally, we thank the editor and reviewers for their kind suggestions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.850084/full#supplementary-material
Supplementary Figure 1 | Control efficiency of R. aceris ZF458 against Agrobacterium tumefaciens on the stem of sunflowers. (A) Agrobacterium tumefaciens ACCC 19185, (B) Agrobacterium tumefaciens ACCC 19185 and R. aceris ZF458 (V:V = 1:1), (C) ZF458, (D) control of sterilized water.
Supplementary Figure 2 | Production of extracellular enzymes and siderophores by R. aceris ZF458. Plate assays for the production of siderophores (A) and phosphatases (B).
Supplementary Figure 3 | General characteristics of R. aceris ZF458. (A) Image of ZF458 colony morphology. (B) Image of ZF458 cells using transmission electron microscopy (Hitachi 7700, Japan). (C) Image of ZF458 cells using scanning electron microscope (Hitachi S3400N, Japan).
Supplementary Figure 4 | Growth curve aND pH curve of R. aceris ZF458 for different culture times. (A) Optical density of ZF458 at different culture times, (B) pH of ZF458 at different culture times.
Supplementary Figure 5 | Phylogenetic tree highlighting the relative positions of R. aceris ZF458 among other Rahnella strains. The phylogenetic tree was constructed based on four housekeeping genes (16S rRNA, gyrB, atpD, and rpoB) according to the aligned gene sequences using maximum likelihoods derived from MEGA 6.0 software. Bootstrap values (1,000 replicates) were shown at the branch points. The scale bar indicates 0.05 nucleotide substitutions per nucleotide position.
Supplementary Figure 6 | Percentage of the average nucleotide identities (ANI) and in silico DNA-DNA hybridization (DDH) among the selected Rahnella strains. ANI values were computed for a pairwise genome comparison using the OrthoANIu algorithm. The percentage of ANI was shown on the top right. DDH values were calculated by using the Genome-to-Genome Distance Calculator (GGDC). The percentage of DDH was shown on the bottom left.
Supplementary Figure 7 | Comparisons of Pyrroloquinoline Quinone genes of R. aceris ZF458 against four other previously fully sequenced Rahnella genomes. The same color represented genes with the same biological function.
Footnotes
- ^ https://www.pacb.com/support/software-downloads/
- ^ http://topaz.gatech.edu/GeneMark/
- ^ http://pfam.xfam.org/
- ^ http://www.ebi.ac.uk/uniprot/
References
Abo-Elyousr, A. M. K., Sallam, M., Hassan, M., and Zeller, W. (2010). Effect of Acibenzolar-S-methyl and Rahnella aquatilis (Ra39) on Chitinase and β-1, 3-glucanase activities and disease resistance of apple plants. Plant Pathol. J. 26, 63–69.
Adeolu, M., Alnajar, S., Naushad, S., and Gupta, S. (2016). Genome-based phylogeny and taxonomy of the ‘Enterobacteriales’: proposal for Enterobacterales ord. nov. divided into the families Enterobacteriaceae, Erwiniaceae fam. nov., Pectobacteriaceae fam. nov., Yersiniaceae fam. nov., Hafniaceae fam. nov., Morganellaceae fam. nov., and Budviciaceae fam. nov. Int. J. Syst. Evol. Microbiol. 66, 5575–5599. doi: 10.1099/ijsem.0.001485
Alt-Mörbe, J., Stryker, J. L., Fuqua, C., Li, P. L., Farrand, S. K., and Winans, S. C. (1996). The conjugal transfer system of Agrobacterium tumefaciens octopine-type Ti plasmids is closely related to the transfer system of an IncP plasmid and distantly related to Ti plasmid vir genes. J. Bacteriol. 178, 4248–4257. doi: 10.1128/jb.178.14.4248-4257
Alvarez-Martinez, C. E., and Christie, P. J. (2009). Biological diversity of prokaryotic type IV Secretion. Microbiol. Mol. Biol. Rev. 73, 775–808. doi: 10.1371/10.1128/MMBR.00023-09
Amos, B., and Rolf, A. (2000). The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucleic Acids Res. 28, 45–48. doi: 10.1093/nar/28.1.45
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. Gene Ontol. Consort. Nat. Genet. 25, 25–29.. doi: 10.1038/75556
Bal, H. B., Das, S., Dangar, T. K., and Adhya, T. K. (2013). ACC deaminase and IAA producing growth promoting bacteria from the rhizosphere soil of tropical rice plants. J. Basic Microbiol. 53, 972–984.. doi: 10.1002/jobm.201200445
Bendtsen, J. D., Binnewies, T. T., Hallin, P. F., Sicheritz-Pontén, T., and Ussery, D. W. (2005). Genome update: prediction of secreted proteins in 225 bacterial proteomes. Microbiology 151, 1725–1727. doi: 10.1099/mic.0.28029-0
Berge, O., Heulin, T., Achouak, W., Richard, C., Bally, R., and Balandreau, J. (1991). Rahnella aquatilis, a nitrogen-fixing enteric bacterium associated with the rhizosphere of wheat and maize. Can. J. Microbiol. 37, 195–203. doi: 10.1139/m91-030
Brady, C., Hunter, G., Kirk, S., Arnold, D., and Denman, S. (2014). Rahnella rasiliensana sp. nov., Rahnella bruchi sp. nov., Rahnella woolbedingensis sp. nov., classification of Rahnella genomospecies 2 and 3 as Rahnella variigena sp. nov. and Rahnella inusitata sp. nov., respectively and emended description of the genus Ra. Syst. Appl. Microbiol. 37, 545–552. doi: 10.1016/j.syapm.2014.09.001
Brewin, N. J., and Legocki, A. B. (1992). Biological nitrogen fixation for sustainable agriculture: A perspective. Plant Soil. 141, 1–11.
Cantarel, B. L., Coutinho, P. M., Corinne, R., Thomas, B., Vincent, L., and Bernard, H. (2009). The Carbohydrate-Active EnZymes database (CAZy): an expert resource for Glycogenomics. Nucleic Acids Res. 37, D233–D238. doi: 10.1093/nar/gkn663
Chang, J. H., Desveaux, D., and Creason, A. L. (2014). The ABCs and 123s of Bacterial Secretion Systems in Plant Pathogenesis. Annu. Rev. Phytopathol. 52, 317–345. doi: 10.1146/annurev-phyto-011014-015624
Chen, F., Li, J. Y., Guo, Y. B., Hui, J., and Wang, H. M. (2009). Biological control of grapevine crown gall: partial characterisation of an antibacterial substance produced by Rahnella aquatilis strain HX2. Eur. J. Plant Pathol. 124, 427–437. doi: 10.1007/s10658-009-9429-z
Choi, O., Kim, J., Kim, J. G., Jeong, Y., and Moon, J. S. (2008). Pyrroloquinoline quinone is a plant growth promotion factor produced by Pseudomonas fluorescens B16. Plant Physiol. 146, 657–668. doi: 10.1104/pp.107.112748
Christie, P. J., and Cascales, E. (2003). The versatile bacterial type IV secretion systems. Nat. Rev. Microbiol. 1, 137–149. doi: 10.1038/nrmicro753
Christie, P. J., Whitaker, N., and González-Rivera, C. (2014). Mechanism and structure of the bacterial type IV secretion systems. Biochim. Biophys. Acta. 1843, 1578–1591. doi: 10.1016/j.bbamcr.2013.12.019
Cianciotto, N. P. (2005). Type II secretion: a protein secretion system for all seasons. Trends Microbiol. 13, 581–588. doi: 10.1016/j.tim.2005.09.005
Eastman, A. W., Heinrichs, D. E., and Yuan, Z. C. (2014). Comparative and genetic analysis of the four sequenced Paenibacillus polymyxa genomes reveals a diverse metabolism and conservation of genes relevant to plant-growth promotion and competitiveness. BMC Genomics 15:851. doi: 10.1186/1471-2164-15-851
el-Hendawy, H. H., Osman, M. E., and Sorour, N. M. (2003). Characterization of two antagonistic strains of Rahnella aquatilis isolated from soil in Egypt. Folia Microbiol. 48, 799–804. doi: 10.1007/BF02931517
Felisberto-Rodrigues, C., Durand, E., Aschtgen, M. S., Blangy, S., Ortiz-Lombardia, M., Douzi, B., et al. (2011). Towards a structural comprehension of bacterial type VI Secretion Systems: Characterization of the TssJ-TssM Complex of an Escherichia coli Pathovar. PLoS Pathog. 7:e1002386. doi: 10.1371/journal.ppat.1002386
Filloux, A. (2004). The underlying mechanisms of type II protein secretion. Biochim. Biophys. Acta. 1694, 163–179. doi: 10.1016/j.bbamcr.2004.05.003
Galperin, M. Y., Makarova, K. S., Wolf, Y. I., and Koonin, E. V. (2015). Expanded microbial genome coverage and improved protein family annotation in the COG database. Nucleic Acids Res. 43, D261–D269. doi: 10.1093/nar/gku1223
Gerlach, R. G., and Hensel, M. (2007). Protein secretion systems and adhesins: The molecular armory of Gram-negative pathogens. Int. J. Med. Microbiol. 297, 401–415. doi: 10.1016/j.ijmm.2007.03.017
Goosen, N., Huinen, R. G., and van de Putte, P. (1992). A 24-amino-acid polypeptide is essential for the biosynthesis of the coenzyme pyrrolo-quinoline-quinone. J. Bacteriol. 174, 1426–1427. doi: 10.1128/jb.174.4.1426-1427.1992
Guo, Y. B., Jiao, Z., Li, L., Wu, D., Crowley, D. E., Wang, Y., et al. (2012). Draft Genome Sequence of Rahnella aquatilis Strain HX2, a plant growth-promoting rhizobacterium isolated from vineyard soil in beijing. China J. Bacteriol. 194, 6646–6647. doi: 10.1128/JB.01769-12
Guo, Y. B., Li, J. Y., Li, L., Chen, F., Wu, W. L., Wang, J. H., et al. (2009). Mutations That Disrupt Either the pqq or the gdh Gene of Rahnella aquatilis abolish the production of an antibacterial substance and result in reduced biological control of grapevine crown gall. Appl. Environ. Microbiol. 75, 6792–6803. doi: 10.1128/AEM.00902-09
He, H., Ye, Z., Yang, D., Yan, J., Xiao, L., Zhong, T., et al. (2013). Characterization of endophytic Rahnella sp. JN6 from Polygonum pubescens and its potential in promoting growth and Cd, Pb, Zn uptake by Brassica napus. Chemosphere 90, 1960–1965. doi: 10.1016/j.chemosphere.2012.10.057
Heeb, S., and Haas, D. (2001). Regulatory roles of the GacS/GacA two-component system in plant-associated and other gram-negative bacteria. Mol. Plant Microbe. Interact. 14, 1351–1363. doi: 10.1094/MPMI.2001.14.12.1351
Heermann, R., and Fuchs, T. M. (2008). Comparative analysis of the Photorhabdus luminescens and the Yersinia enterocolitica genomes: uncovering candidate genes involved in insect pathogenicity. BMC Genomics 9:40. doi: 10.1186/1471-2164-9-40
Ho, B. T., Dong, T. G., and Mekalanos, J. J. (2014). A view to a kill: the bacterial type VI secretion system. Cell Host Microbe. 15, 9–21. doi: 10.1016/j.chom.2013.11.008
Ibtissem, G., Gilles, V., and Christine, P. (2007). CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats. Nucleic Acids Res. 35, W52–W57 doi: 10.1093/nar/gkm360
Jochimsen, B., Lolle, S., McSorley, F. R., Nabi, M., Stougaard, J., Zechel, D. L., et al. (2011). Five phosphonate operon gene products as components of a multi-subunit complex of the carbon-phosphorus lyase pathway. Proc. Natl. Acad. Sci. USA 108, 11393–11398. doi: 10.1073/pnas.1104922108
Johnson, T. L., Abendroth, J., Hol, W., and Sandkvist, M. (2006). Type II secretion: From structure to function. FEMS Microbiol. Lett. 255, 175–186.
Kanehisa, M., Goto, S., Hattori, M., Aoki-Kinoshita, K. F., and Hirakawa, M. (2006). From genomics to chemical genomics: New developments in KEGG. Nucleic Acids Res. 34, D354–D357. doi: 10.1093/nar/gkj102
Khan, Z., and Doty, S. L. (2009). Characterization of bacterial endophytes of sweet potato plants. Plant Soil 322, 197–207.
Kim, G. H., Fratamico, P., Breidt, F., and Oh, D. H. (2016). Survival and expression of acid resistance genes in Shiga toxin-producing Escherichia coli acid adapted in pineapple juice and exposed to synthetic gastric fluid. J. Appl. Microbiol. 121, 1416–1426. doi: 10.1111/jam.13223
Kong, W. L., Rui, L., Ni, H., and Wu, X. Q. (2020a). Antifungal Effects of Volatile Organic Compounds Produced by Rahnella aquatilis JZ-GX1 Against Colletotrichum gloeosporioides in Liriodendron chinense × tulipifera. Front. Microbiol. 11:1114. doi: 10.3389/fmicb.2020.01114
Kong, W. L., Wu, X. Q., and Zhao, Y. J. (2020b). Effects of Rahnella aquatilis JZ-GX1 on Treat Chlorosis Induced by Iron Deficiency in Cinnamomum camphora. J. Plant Growth Regual. 39, 877–887
Korotkov, K. V., Johnson, T. L., Jobling, M. G., Pruneda, J., Pardon, E., Héroux, A., et al. (2011). Structural and functional studies on the interaction of GspC and GspD in the type II secretion system. PLoS Pathog. 7:e1002228. doi: 10.1371/journal.ppat.1002228
Lagesen, K., Hallin, P., Rødland, E. A., Staerfeldt, H. H., Rognes, T., and Ussery, D. W. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108. doi: 10.1093/nar/gkm160
Lartigue, M. F., Nordmann, P., Edelstein, M. V., Cuzon, G., Brisse, S., and Poirel, L. (2013). Characterization of an extended-spectrum class A β-lactamase from a novel enterobacterial species taxonomically related to Rahnella spp./Ewingella spp. J. Antimicrob. Chemother. 68, 1733–1736. doi: 10.1093/jac/dkt122
Lavín, J., Kiil, K., Resano, O., Ussery, D. W., and Oguiza, J. A. (2007). Comparative genomic analysis of two-component regulatory proteins in Pseudomonas syringae. BMC Genomics 8:397–390. doi: 10.1186/1471-2164-8-397
Lee, S. D., Jeon, D., Kim, I. S., Choe, H., and Kim, J. S. (2020). Rahnella aceris sp. nov., isolated from sap drawn from Acer pictum. Arch Microbiol. 202:9. doi: 10.1007/s00203-020-01961-5
Li, G. E., Wu, X. Q., Ye, J. R., Hou, L., and Zhou, A. D. (2013). Isolation and identification of phytate-degrading rhizobacteria with activity of improving growth of poplar and Masson pine. World J. Microbiol. Biotechnol. 29, 2181–2193. doi: 10.1007/s11274-013-1384-3
Li, L., Jiao, Z., Hale, L., Wu, W. L., and Guo, Y. B. (2014). Disruption of Gene pqqA or pqqB Reduces plant growth promotion activity and biocontrol of crown gall disease by Rahnella aquatilis HX2. PLoS One 9:e115010. doi: 10.1371/journal.pone.0115010
Li, L., Li, J., Peng, J., Wu, W., and Guo, Y. (2019). Identification of atpD as an optimal reference gene to explore antibiotic resistance and stress tolerance in Rahnella aquatilis. J. Appl. Microbiol. 126, 1096–1107. doi: 10.1111/jam.14215
Li, P. S., Kong, W. L., and Wu, X. Q. (2021). Salt tolerance mechanism of the rhizosphere bacterium JZ-GX1 and its effects on tomato seed germination and seedling growth. Front. Microbiol. 12:657238. doi: 10.3389/fmicb.2021.657238
Li, W., Lukasz, J., and Adam, G. (2002). Tolerating some redundancy significantly speeds up clustering of large protein databases. Bioinformatics 18, 77–82. doi: 10.1093/bioinformatics/18.1.77
Liu, R., and Ochman, H. (2007). Stepwise formation of the bacterial flagellar system. Proc. Natl. Acad. Sci. U.S.A. 104, 7116–7121. doi: 10.1073/pnas.0700266104
Lowe, T. M., and Eddy, S. R. (1997). tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955–964. doi: 10.1093/nar/25.5.955
Luis, S. D. P., Prieto, C., and Marcela, C. F. (2014). Screening phosphate solubilizing actinobacteria isolated from the rhizosphere of wild plants from the Eastern Cordillera of the Colombian Andes. Afr. J. Microbiol. Res. 8, 734–742. doi: 10.5897/AJMR2013.5940
Ma, L. S., Lin, J. S., and Lai, E. M. (2009). An IcmF Family Protein, ImpLM, is an integral inner membrane protein interacting with ImpKL, and its walker a motif is required for type VI secretion system-mediated Hcp Secretion in Agrobacterium tumefaciens. J. Bacteriol. 191, 4316–4329. doi: 10.1128/JB.00029-09
Majdalani, N., and Gottesman, S. (2005). The Rcs phosphorelay: a complex signal transduction system. Annu. Rev. Microbiol. 59, 379–405. doi: 10.1146/annurev.micro.59.050405.101230
Martinez, R. J., Bruce, D., De Tter, C., Goodwin, L. A., Han, J., Han, C. S., et al. (2012a). Complete Genome Sequence of Rahnella sp. Strain Y9602, a gammaproteobacterium isolate from metal- and radionuclide-contaminated soil. J. Bacteriol. 194, 2113–2114. doi: 10.1128/JB.00095-12
Martinez, R. J., Bruce, D., Detter, C., Goodwin, L. A., Han, J., Han, C. S., et al. (2012b). Complete Genome Sequence of Rahnella aquatilis CIP 78.65. J. Bacteriol. 194, 3020–3021. doi: 10.1128/JB.00380-12
Masuda, N., and Church, G. M. (2010). Regulatory network of acid resistance genes in Escherichia coli. Mol. Microbiol. 48, 699–712. doi: 10.1046/j.1365-2958.2003.03477.x
Matsushita, K., Toyama, H., Yamada, M., and Adachi, O. J. (2002). Quinoproteins: structure, function, and biotechnological applications. Plant Physiol. 58, 13–22. doi: 10.1007/s00253-001-0851-1
Meier-Kolthoff, J. P., Auch, A. F., Klenk, H. P., and Gker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 14:60. doi: 10.1186/1471-2105-14-60
Milagres, A. M., Machuca, A., and Napoleão, D. (1999). Detection of siderophore production from several fungi and bacteria by a modification of chrome azurol S (CAS) agar plate assay. J. Microbiol. Methods. 37, 1–6. doi: 10.1016/s0167-7012(99)00028-7
Nawrocki, E. P., Kolbe, D. L., and Eddy, S. R. (2009). Infernal 1.0: inference of RNA alignments. Bioinformatics 25, 1335–1337. doi: 10.1093/bioinformatics/btp157
Oldroyd, G. E., and Dixon, R. (2014). Biotechnological solutions to the nitrogen problem. Curr. Opin. Biotechnol. 26, 19–24.
Palmieri, D., Vitale, S., Lima, G., Di Pietro, A., and Turra, D. (2020). A bacterial endophyte exploits chemotropism of a fungal pathogen for plant colonization. Nat. Commun. 11, 52–64. doi: 10.1038/s41467-020-18994-5
Passariello, C., Forleo, C., Micheli, V., Schippa, S., Leone, R., Mangani, S., et al. (2006). Biochemical characterization of the class B acid phosphatase (AphA) of Escherichia coli MG1655. Biochim. Biophys. Acta 1764, 13–19. doi: 10.1016/j.bbapap.2005.08.028
Peng, J., Wu, D., Liang, Y., Li, L., and Guo, Y. (2019). Disruption of acdS gene reduces plant growth promotion activity and maize saline stress resistance by Rahnella aquatilis HX2. J. Basic Microbiol. 59, 402–411. doi: 10.1002/jobm.201800510
Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786. doi: 10.1038/nmeth.1701
Qiao, Z. W., Hong, J. P., Xie, Y. H., and Xuan, L. L. (2013). Screening, identification and phosphate-solubilizing characteristics of Rahnella sp. phosphate-solubilizing bacteria in calcareous soil. Ying Yong Sheng Tai Xue Bao 24, 2294–2300.
Rhodes, A. N., Urbance, J. W., Youga, H. Corlew-Newman, H., Reddy, C. A., Klug, M. J., et al. (1998). Identification of bacterial isolates obtained from intestinal contents associated with 12,000-year-old mastodon remains. Appl. Environ. Microbiol. 64, 651–658. doi: 10.1089/oli.1.1998.8.67
Richter, M., and Rossello-Mora, R. (2009). Shifting the genomic gold standard for the prokaryotic species definition. Proc. Natl. Acad. Sci. USA 106, 19126–19131. doi: 10.1073/pnas.0906412106
Rozhon, W., Petutschnig, E., Khan, M., Summers, D. K., and Poppenberger, B. (2010). Frequency and diversity of small cryptic plasmids in the genus Rahnella. BMC Microbiol. 10:56. doi: 10.1186/1471-2180-10-56
Saier, M. H., Reddy, V. S., Tamang, D. G., and Research, V. A. (2009). The Transporter Classification Database: recent advances. Nucleic Acids Res. 37, D274–D278. doi: 10.1093/nar/gkn862
Sallam, N. M. (2011). Biological control of common blight of bean (Phaseolus vulgaris) caused by Xanthomonas axonopodis pv. phaseoli by using the bacterium Rahnella aquatilis. Arch. Phytopathol. Plant Protec. 44, 1966–1975.
Santos, P. D., Fang, Z., Mason, S. W., Setubal, J. O. C., and Dixon, R. (2012). Distribution of nitrogen fixation and nitrogenase-like sequences amongst microbial genomes. BMC Genomics 13:162. doi: 10.1186/1471-2164-13-162
Skalickova, S., Milosavljevic, V., Cihalova, K., Horky, P., Richtera, L., and Adam, V. (2017). Selenium nanoparticles as a nutritional supplement. Nutrition 33, 83–90. doi: 10.1016/j.nut.2016.05.001
Spaepen, S., Vanderleyden, J., and Remans, R. (2007). Indole-3-acetic acid in microbial and microorganism-plant signaling. FEMS Microbiol. Rev. 31, 425–448. doi: 10.1111/j.1574-6976.2007.00072.x
Stolz, J. F., Basu, P., Santini, J. M., and Oremland, R. S. (2006). Arsenic and selenium in microbial metabolism. Annu. Rev. Microbiol. 60, 107–130. doi: 10.1146/annurev.micro.60.080805.142053
Stothard, P., and Wishart, D. C. (2005). Circular genome visualization and exploration using CGView. Bioinformatics 21, 537–539. doi: 10.1093/bioinformatics/bti054
Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol. 30, 2725–2729.
Trokter, M., Felisberto-Rodrigues, C., Christie, P. J., and Waksman, G. (2014). Recent advances in the structural and molecular biology of type IV secretion systems. Curr. Opin. Struct. Biol. 27, 16–23. doi: 10.1016/j.sbi.2014.02.006
Vassilev, N., Vassileva, M., and Nikolaeva, I. (2006). Simultaneous P-solubilizing and biocontrol activity of microorganisms: potentials and future trends. Appl. Microbiol. Biotechnol. 71, 137–144. doi: 10.1007/s00253-006-0380-z
Velterop, J. S., Sellink, E., Meulenberg, J. J., David, S., Bulder, I., and Postma, P. W. (1995). Synthesis of pyrroloquinoline quinone in vivo and in vitro and detection of an intermediate in the biosynthetic pathway. J. Bacteriol. 177, 5088–5098. doi: 10.1128/jb.177.17.5088-5098.1995
Wang, B., Yuan, J., Jian, Z., Zhang, M., Li, R., Ruan, Y., et al. (2013). Effects of novel bioorganic fertilizer produced by Bacillus amyloliquefaciens W19 on antagonism of Fusarium wilt of banana. Biol. Fert Soils. 49, 435–446. doi: 10.1007/s00374-012-0739-5
Wang, H., Zhang, J., and Yu, H. (2007). Elemental selenium at nano size possesses lower toxicity without compromising the fundamental effect on selenoenzymes: comparison with selenomethionine in mice. Free Radic. Biol. Med. 42, 1524–1533. doi: 10.1016/j.freeradbiomed.2007.02.013
Wang, L., Zhang, L., Liu, Z., Zhao, D., Liu, X., Zhang, B., et al. (2013). A Minimal Nitrogen Fixation Gene Cluster from Paenibacillus sp. WLY78 Enables Expression of Active Nitrogenase in Escherichia coli. PLoS Genet. 9:e1003865. doi: 10.1371/journal.pgen.1003865
Werra, P. D., Péchy-Tarr, M., Keel, C., and Maurhofer, M. (2009). Role of Gluconic Acid Production in the Regulation of Biocontrol Traits of Pseudomonas fluorescens CHA0. Appl. Environ. Microbiol. 75, 4162–4174. doi: 10.1128/AEM.00295-09
White, A. K., and Metcalf, W. W. (2007). Microbial metabolism of reduced phosphorus compounds. Annu. Rev. Microbiol. 61, 379–400. doi: 10.1146/annurev.micro.61.080706.093357
Winans, S. C., and Walker, G. C. (1985). Conjugal transfer system of the IncN plasmID pKM101. J. Bacteriol. 161, 402–410. doi: 10.1128/jb.161.1.402-410.1985
Xu, Q., Song, Y., Lin, Z., Bauelos, G., and Guo, Y. B. (2020). The small RNA chaperone Hfq is a critical regulator for bacterial biosynthesis of selenium nanoparticles and motility in Rahnella aquatilis. Appl. Microbiol. Biotechnol. 104, 1721–1735. doi: 10.1007/s00253-019-10231-4
Yamada, M., Makino, K., Amemura, M., Shinagawa, H., and Nakata, A. (1989). Regulation of the phosphate regulon of Escherichia coli: analysis of mutant phoB and phoR genes causing different phenotypes. J. Bacteriol. 171, 5601–5606. doi: 10.1128/jb.171.10.5601-5606.1989
You, Z., Liang, Y., Lynch, K. H., Dennis, J. J., and Wishart, D. S. (2011). PHAST: A Fast Phage Search Tool. Nucleic Acids Res. 39, W347–W352. doi: 10.1093/nar/gkr485
Yuan, L., Li, L., Zheng, F., Shi, Y., Xie, X., Chai, A., et al. (2020). The complete genome sequence of Rahnella aquatilis ZF7 reveals potential beneficial properties and stress tolerance capabilities. Arch Microbiol. 202, 483–499. doi: 10.1007/s00203-019-01758-1
Zhang, Y. C., Fan, Q. R., and Loria, R. (2016). A re-evaluation of the taxonomy of phytopathogenic genera Dickeya and Pectobacterium using whole-genome sequencing data. Syst. Appl. Microbiol. 39, 252–259. doi: 10.1016/j.syapm.2016.04.001
Keywords: Rahnella aceris ZF458, comparative genomic analysis, biological control, plant growth promotion, environmental adaptation
Citation: Xu S, Zhao Y, Peng Y, Shi Y, Xie X, Chai A, Li B and Li L (2022) Comparative Genomics Assisted Functional Characterization of Rahnella aceris ZF458 as a Novel Plant Growth Promoting Rhizobacterium. Front. Microbiol. 13:850084. doi: 10.3389/fmicb.2022.850084
Received: 07 January 2022; Accepted: 24 February 2022;
Published: 04 April 2022.
Edited by:
Martin Filion, Agriculture and Agri-Food Canada (AAFC), CanadaReviewed by:
Yunzeng Zhang, Yangzhou University, ChinaVíctor González, National Autonomous University of Mexico, Mexico
Copyright © 2022 Xu, Zhao, Peng, Shi, Xie, Chai, Li and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lei Li, lilei01@caas.cn; Baoju Li, libaojuivf@163.com