J Genomics 2024; 12:26-34. doi:10.7150/jgen.91337 This volume Cite

Research Paper

Isolation and draft genome sequence of Enterobacter asburiae strain i6 amenable to genetic manipulation

Akinori Kato Corresponding address

Department of Advanced Bioscience, Graduate School of Agriculture, Kindai University, 3327-204 Nakamachi, Nara, Nara 631-8505, Japan.

Citation:
Kato A. Isolation and draft genome sequence of Enterobacter asburiae strain i6 amenable to genetic manipulation. J Genomics 2024; 12:26-34. doi:10.7150/jgen.91337. https://www.jgenomics.com/v12p0026.htm
Other styles

File import instruction

Abstract

Graphic abstract

Enterobacter asburiae is a species of Gram-negative bacteria that is found in soil, water, and sewage. E. asburiae is generally considered to be an opportunistic pathogen, but has also been reported as a plant growth-promoting bacterium (PGPB), which may have beneficial effects on plant growth and development. However, genetic analysis of E. asburiae has been limited, possibly due to its redundant enzymes that digest exogenous DNA in the cell. Here, an E. asburiae strain i6 was isolated from soil in Nara, Japan. This strain was amenable to transformation and the one-step gene inactivation method based on λ Red recombinase. The transformation efficiency of the i6 strain with the 10 kb plasmid DNA pCF430 was at least four orders of magnitude higher than that of the previously sequenced E. asburiae strain ATCC 35953, which could not be transformed with the same plasmid DNA. A draft genome sequence of the i6 strain was determined and deposited into the database, allowing several factors that may determine transformation efficiency to be perturbed and tested. Together with the amenability of the i6 strain to genetic manipulation, the information from the i6 genome will facilitate characterization and fine-tuning of the beneficial and detrimental traits of this species.

Keywords: Enterobacter asburiae, gene deletion, λ Red recombinase, plant growth-promoting bacterium (PGPB), cell-cell interactions

Introduction

Enterobacter asburiae is a Gram-negative, facultative anaerobic, oxidase-negative, non-motile, and non-pigmented rod-shaped species belonging to the genus Enterobacter, which includes species that are difficult to identify with biochemical and phylogenetic tests [1-3]. The Enterobacter cloacae complex (ECC) is a group of common opportunistic pathogens consisting of Enterobacter asburiae, Enterobacter cloacae, Enterobacter hormaechei, Enterobacter kobei, Enterobacter ludwigii, Enterobacter mori, and Enterobacter nimipressuralis. In addition, recently identified species, including Enterobacter roggenkampii, Enterobacter chengduensis, and Enterobacter bugandensis, are clustered with the ECC species [4]. While E. asburiae, as a part of ECC, is considered to be an opportunistic pathogen, not only E. asburiae but also some of Enterobacter genus have been reported as plant growth-promoting bacteria (PGPB). Examples include E. asburiae PDA 134 from date palm [5], E. cloacae from citrus and maize plants [6, 7], and E. asburiae from sweet potato [8]. Enterobacter sp. strain P23 promotes rice growth under salt stress. In addition, E. mori, E. asburiae, E. ludwigii, and E. sp. J49 have been shown to promote wheat growth under stress conditions [9]. At least some strains of E. asburiae reduce the epiphytic fitness of the human enteric pathogens E. coli O157:H7 and Salmonella on lettuce and Arabidopsis by at least 100-fold [9, 10]. While E. asburiae is a preferable candidate species for genome editing to safely further enhance the ability of PGPB and to study cell-cell interactions, genetic analysis of this important species is limited [11, 12], possibly due to its redundant enzymes that digest exogenous DNA within the cell.

The one-step gene inactivation method based on λ Red recombinase is a powerful and efficient technique used to disrupt specific genes [13] and to construct gene fusions [14] in the bacterial genome. Soon after this method was originally invented in Escherichia coli [13], it was widely used to construct deletion mutants and insertion mutants of other Gram-negative species, such as Salmonella enterica [15], Yersinia pestis [16], Klebsiella pneumoniae [17], Pantoea ananatis [18]. Basically, the technique uses knockout cassettes with short (40-60 bp) homologous arms that can be produced in a single PCR reaction, but some species or strains only accept knockout cassettes with extended arms (200~1,000 bp), which require additional PCR steps to successfully manipulate the gene with λ Red recombinase [19, 20]. To date, the one-step gene inactivation method has been unsuccessful with the E. asburiae ATCC 35953 (NBRC 109912 T) strain. In addition to genetic factors that reduce the stability of linear DNA fragments, genetic factors that reduce transformation efficiency have generally been thought to hinder genetic manipulations, such as the one-step gene inactivation method. Here, I report a strain of E. asburiae i6, isolated from soil in Nara Japan, that is amenable to genetic manipulation.

Materials and Methods

Bacterial strains, plasmids, and growth conditions

Bacterial strains and plasmids used in this study are listed in Table 1. Bacteria were grown at 37ºC in LB Broth (Lennox). Ampicillin and spectinomycin were used at 100 µg/ml, kanamycin at 50 µg/ml, chloramphenicol at 25 µg/ml, and tetracycline at 12.5 µg/ml. Primers used in this study are listed in Table 2.

PCR. PCRs for 16S and groEL were performed using primer sets 63f/1387r [21] and Hsp60-F/Hsp60-R [22] according to their protocols [21, 22]. Purified PCR products were analyzed by DNA sequencing.

Construction of chromosomal gene deletion mutants

Strain AK1602, which has a deletion of the rcsB gene and a CmR cassette, was constructed by the one-step inactivation method [13] using primers A2414 and A2415 with pKD3 as the template.

Strain AK1603, which has a deletion of the klcA1 gene and a tetA, was constructed by the one-step inactivation method [13] using primers A2418 and A2419 with the phoQ::Tn10 chromosomal DNA as the template.

 Table 1 

Bacterial strains and plasmids used in this study.

Strain or plasmidDescriptionReference or source
E. asburiae
109912 T(ATCC 35953)Wild typeNBRC
i6Wild typeThis work
AK1599i6 wzc-FRT-KmR-FRTThis work
AK1602i6 ∆rcsB::FRT-CmR-FRTThis work
AK1603i6 ∆klcA1::TcRThis work
AK1604i6 ∆klcA2::FRT-KmR-FRTThis work
AK1605i6 ∆dgeC::FRT-KmRThis work
E. coli
DH5αF- Φ80dlacZΔM15 Δ(lacZYA-argF)U169 deoR recA1 endA1 hsdR17(rK- mK+) phoA supE44 thi-1 gyrA96 relA1Laboratory stock
JM109λpirrecA1 endA1 gyrA96 thi1 hsdR17(rK- mK+) e14- (mcrA-) supE44 relA1 Δ(lac-proAB)/F' [traD36, proAB+, lacIq, lacZΔM15] λpirLaboratory stock
Salmonella enterica
MS5996phoQ::Tn10[40]
EG16468∆PpmrD-pmrD::SpR[23]
Plasmids
pKD3repR6Kγ ApR FRT CmR FRT[36]
pKD4repR6Kγ ApR FRT KmR FRT[36]
pCP20reppSC101ts ApR CmR cl857 λPR flp[41]
pCF430repRK2 oriTRK2 pBAD araC TcR[33]
pACYC177repp15A ApR KmR[35]
RSFRedTERrepRSF1010 oriVRSF1010 lacI sacB CmR γ β exo[18]
RSFRedTER-SprepRSF1010 oriVRSF1010 lacI sacB SpR γ β exoThis work
pAK1001reporiS, CmR, FRT-KmR-FRT[34]
 Table 2 

Primers used in this study.

Primer nameSequence
63fCAGGCCTAACACATGCAAGTC
1387rGGGCGGWGTGTACAAGGC
Hsp60-FGGTAGAAGAAGGCGTGGTTGC
Hsp60-RATGCATTCGGTGGTGATCATCAG
A1190CAGCATCCTTGAACAAGGACAATTAACAGTTAACAAATAAGCTGTAATGCAAGTAGCG
A1191AGGTGGGACCACCCGCGCTACTGCCGCCAGGCAAAGAATCTTTATTTGCCGACTACCTTGA
A2414CTCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGGATTTTTAGTGTAGGCTGGAGCTGCTTC
A2415CGGCTATTACGAGTACGAATATAAATCTGACAGCAAATAACATATGAATATCCTCCTTAG
A2416AAAAGCTTGCTGTAGCAAGGTAGCCTTATACATGAACAATGTGTAGGCTGGAGCTGCTTC
A2417TCCCCGCGGGAGAGGGACGGGGTGAGACACCCGTCCGGGACATATGAATATCCTCCTTAG
A2418GCGATGCCGGAACCGAGAATATGTTAACCGTGGAGGATCTGTGTAGGCTGGAGCTGCTTC
A2419TCTCCTTTAGTCGGGCGGCATCGCCGCCCGTGGTCAGTCACATATGAATATCCTCCTTAG
A2420CGGGCATTAAACAACAGCACCAGCATAAGGAAATAGTTTACTCTAATGCGCTGTTAATCACT
A2421GCACCGGGGATAATGCAGCATTCATGTTCGGCATTCCTCACTAAGCACTTGTCTCCTGTT
A2424AATTATCAAATTCACGCGTTCGGACATCTTCCCTGACGGCGGAATAGGAACTTCAAGATC
A2425GCTGCTGCATATCTGGTATACCGACGGCTGGACGCCGTCACATATGAATATCCTCCTTAG

Strain AK1599, which has an insertion of the KmR cassette behind the wzc gene, was constructed by the one-step inactivation method [13] using primers A2416 and A2417 with pKD4 as the template.

Strain AK1604, which has a deletion of the klcA2 gene and a KmR cassette, was constructed by the one-step inactivation method [13] using primers A2420 and A2421 with pKD4 as the template.

Strain AK1605, which has a deletion of the dgeC gene and a KmR cassette, was constructed by the one-step inactivation method [13] using primers A2424 and A2425 with pKD4 as the template.

Plasmid construction

Plasmid RSFRedTER-Sp for λ Bet Exo Gam expression was constructed as a spectinomycin-resistant version of RSFRedTER [18] by the one-step inactivation method [13] using primers A1190 and A1191 with EG16468 chromosomal DNA [23] as the SpR marker template.

All strains constructed using PCR reactions were analyzed by DNA sequencing to confirm that the PCR-generated DNA regions had the predicted sequences.

Transformation efficiency analysis

Bacterial cells cultured overnight in LB (Lennox) were added to 10 ml of fresh medium diluted 1:50 in test tubes and shaken at 37°C until the OD595 value reached 0.3~0.8. Cells were collected by centrifugation and the supernatant was removed. The cell pellet was washed twice with 1 ml of 10% glycerol, and then the cells were suspended in 400 µl of 10% glycerol. 1-3 µl of plasmid DNA was mixed with 100 µl of competent cells (the cell suspension) and electroporated using a 2 mm cuvette and a Gene Pulser II (Bio-Rad) at 2.5 kV, 200 Ω, and 25 µF. Immediately after the pulse, SOB medium was added to the cuvette and collected in 1.5 ml tubes. Cells were spread on LB (Lennox) plates containing kanamycin (pAK1001 and pACYC177) or tetracycline (pCF430) and incubated at 37°C overnight. Colonies were then counted.

Genome sequencing, assembly and annotation

Whole genome sequencing was performed using DNBSEQ-G400 (MGI Tech Co., Ltd.). Libraries were prepared using IMGIEasy FS DNA Library Prep Set (MGI Tech Co., Ltd.) and fragments were validated using Fragment Analyzer and dsDNA 915 Reagent Kit (Agilent Technologies). Sequences were de novo assembled using SPAdes assembler version 3.15.5 [42]. Genome annotation was performed using DFAST.

Average nucleotide identity analysis

Average nucleotide identity (ANI) analysis [24] was performed using the complete genome sequences of all Enterobacter species available in the NCBI database.

Results and Discussion

Isolation and genome sequencing of the Enterobacter asburiae i6

On July 6, 2016, an E. asburiae strain was isolated from the soil near the lawn of the campus of Kindai University Faculty of Agriculture in Nara, Japan, as a strain that somewhat affects the luminescence of a Salmonella enterica strain carrying plasmid luciferase reporter under the control of the RcsB-regulated wza gene promoter, when cross-streaked. In addition, the colony morphology of this strain is slightly mucoid (or encapsulated) on LB agar plate. 16s rRNA and groEL sequencing determined that the strain was approximately E. asburiae, and this strain was designated i6. Because the i6 strain may produce a signal molecule that is recognized by S. enterica by an unknown mechanism, it would be useful to sequence the genome of strain i6 for future analyses. Therefore, chromosomal DNA was extracted from strain i6 and the genome was sequenced on a next-generation sequencer (DNBSEQ-G400) with X290 coverage. Automated assembly resulted in the formation of 28 contigs, and the longer contigs were dissected by redundant copies of 5s rRNA, 16S rRNA, 23 rRNA, and tRNA. The assembled genome sequences and their annotations were deposited in a database (accession numbers: BTPF01000001~BTPF01000028, DRR495940, PRJDB16388, and SAMD00634955). The sequences of the ATCC 35953, FDAARGOS_892, L1, and RHBSTW-01009 strains, all of which belong to E. asburiae (sensu stricto) and not to the four provisional classes (B-E) of E. asburiae according to GTDB [25], were the closest to i6 among the complete E. asburiae genomes (Table 3). The FDAARGOS_892 strain had the highest ANI value of 98.75%, whereas the L1 strain had the longest average alignment length of 70.08%. Thus, the i6 strain is most likely to be E. asburiae (sensu stricto). Alignment of the ends of each contig sequence of i6 to these sequences suggested that the large contigs 1-7, 9, and 10 form an almost complete chromosome of i6 (Table 4), and that contig 8 is a single plasmid DNA itself. Furthermore, the gaps between these contigs were shown to correspond to repetitive 16S rRNA and 23 rRNA, tRNA-Ala, tRNA-Ile, and tRNA-Glu (Table 4).

Possible genes involved in plant growth promoting traits in the E. asburiae i6 genome

Although strain i6 was not originally isolated as PGPB, a list of candidate genes that may contribute to plant growth was compiled from the genome sequence of strain i6 (Table 5). Soil beneficial bacteria can promote plant growth by synthesizing molecules similar to plant hormones [26, 27]. Auxin like indole acetic acid (IAA) is quantitatively the most plant hormones secreted by Phytophthora rhizobacteria [28, 29], suggesting that auxin is a signaling molecule in microorganisms [30]. Similar to the genome features previously reported for Enterobacter sp. J49 [30], genes responsible for the synthesis of indole-pyruvate decarboxylase and indole-3-acetaldehyde dehydrogenase were detected in the i6 genome, and no other IAA pathway-related genes were detected. The bacterial volatiles 2,3-butanediol and acetoin are also known to induce plant growth promotion [31]. Acetolactate synthase BudB converts pyruvate to acetolactate, which is subsequently converted to acetoin by acetoin decarboxylase BudA. The (S)-acetoin-forming diacetyl reductase BudC catalyzes the conversion of acetoin to 2,3-butanediol, which is reversible. All genes for butA, butB and butC were found, as well as other genes encoding acetolactate synthases (ilvB, ilvN_2, ilvG etc.). Similar to the J49 genome [30], the gcd gene encoding GDH synthesis, responsible for the production of gluconic acid, the major organic acid in the phosphate solubilization mechanism most widely used by soil bacteria, was detected in the i6 genome, but the pqq gene cluster, required for the biosynthesis of the PQQ cofactor, was not detected. Redundant genes for siderophore-production, iron ABC transporters, and type VI secretion systems which may compete with (plant) pathogens for low iron levels and cell growth, were also found in the i6 genome (data not shown).

 Table 3 

Average nucleotide identity analysis (ANI) data calculated from the nearly complete genome of Enterobacter asburiae i6 strain and the whole genome sequence of each Enterobacter strains.

Bacterial strainsTotal length (bp)Average aligned length (bp)Average aligned length (%) respect to Enterobacter asburiae i6ANI value (%) respect to Enterobacter asburiae i6
Enterobacter asburiae i64,623,6604,623,660100.00100.00
Enterobacter asburiae (sensu stricto) FDAARGOS_8924,784,8203,045,12065.8698.75
Enterobacter asburiae (sensu stricto) ATCC 359534,804,2003,122,44667.5398.71
Enterobacter asburiae (sensu stricto) RHBSTW-010094,652,2202,968,01964.1998.67
Enterobacter asburiae (sensu stricto) L14,626,7203,240,16470.0898.63
Enterobacter asburiae (B) MNCRE144,491,0603,000,47064.8997.01
Enterobacter asburiae (B) UBA118994,111,6202,892,02162.5596.99
Enterobacter asburiae (B) 17Nkhm-UP24,709,3403,041,45365.7896.95
Enterobacter asburiae (B) UBA82644,433,9403,091,30566.8696.93
Enterobacter asburiae (B) 1808-0134,769,5203,236,90670.0196.91
Enterobacter asburiae (D) R_A5.MM4,837,8603,144,46468.0195.90
Enterobacter asburiae (C) TN1525,042,8802,899,55062.7195.22
Enterobacter roggenkampii DSM 166904,899,0603,026,51765.4693.03
Enterobacter chengduensis WCHECl-C45,183,6402,917,80463.1193.21
Lelliottia nimipressuralis 51 (Enterobacter nimipressuralis)4,862,3402,888,08262.4692.70
Enterobacter bugandensis EB-2474,717,5002,825,93159.9091.52
Enterobacter mori ACYC.E9L4,807,2602,814,94260.8890.17
Enterobacter cloacae subsp. cloacae ATCC 13047 (GCF_000025565.1)5,596,7402,783,20160.1988.73
Enterobacter ludwigii EN-1194,952,1002,777,48856.0988.54
Enterobacter hormaechei ATCC 491624,884,7802,696,35155.2087.69
Enterobacter asburiae (E) INSAq1464,456,3802,514,95554.3986.68

GTDB [25] categories for E. asburiae, sensu stricto or provisional classes (B-E), are indicated in bold.

 Table 4 

Contig sequences and gaps in E. asburiae i6 assigned to the complete genome of the closest E. asburiae strains.

Bacterial strainEnterobacter asburiae (sensu stricto) ATCC 35953
i6 contig#12475 comp9 comp1063 comp1
Aligned 5'-end position on ATCC 35953 (nt)0169175626877143150068331519535888083689839377671439734044652519
Aligned 3'-end position on ATCC 35953 (nt)1691839268307731449433310295358390936846943771799396859446478664713742
Gap to right next contig (nt)-8346375125490048995145491548104653
Gene assigned at gap right next contig23S rRNA23S rRNA23S rRNA16S rRNA16S rRNA16S rRNA16S rRNA16S rRNA
tRNA-AlatRNA-Glu16S rRNAtRNA-GlutRNA-IletRNA-IletRNA-GlutRNA-Glu
tRNA-Ile16S rRNA23S rRNAtRNA-AlatRNA-Ala23S rRNA23S rRNA
16S rRNA5S rRNA23S rRNA23S rRNA
5S rRNA
Bacterial strainEnterobacter asburiae (sensu stricto) FDAARGOS_892
i6 contig#1 comp475 comp9 comp1063 comp2 comp1 comp
Aligned 5'-end position on FDAARGOS_892 (nt)01764447999864512391875710197911106666130335919818652972903
Aligned 3'-end position on FDAARGOS_892 (nt)12628474876640226913839101464611017531298549197760229729864717539
Gap to right next contig (nt)50165122489749185145491348104263-83
Gene assigned at gap right next contig23S rRNA23S rRNA23S rRNA16S rRNA16S rRNA16S rRNA16S rRNA16S rRNA
tRNA-AlatRNA-Glu16S rRNAtRNA-GlutRNA-IletRNA-IletRNA-GlutRNA-Glu
tRNA-Ile16S rRNA23S rRNAtRNA-AlatRNA-Ala23S rRNA23S rRNA
16S rRNA5S rRNA23S rRNA23S rRNA
5S rRNA
Bacterial strainEnterobacter asburiae (sensu stricto) L1
i6 contig#3 comp2 comp1 comp475 comp9 comp1063 comp
Aligned 5'-end position on L1 (nt)048298913327443039704352743336737983956299406059541474964348070
Aligned 3'-end position on L1 (nt)478726133282730346913522308366865139513794055451414258543432614561905
Gap to right next contig (nt)4263-835013512551474920514449114809
Gene assigned at gap right next contig16S rRNA23S rRNA23S rRNA5S rRNA16S rRNA16S rRNA16S rRNA16S rRNA
tRNA-GlutRNA-AlatRNA-Glu23S rRNAtRNA-GlutRNA-IletRNA-IletRNA-Glu
23S rRNAtRNA-Ile16S rRNA16S rRNA23S rRNAtRNA-AlatRNA-Ala23S rRNA
16S rRNA5S rRNA23S rRNA23S rRNA
5S rRNA
Bacterial strainEnterobacter asburiae (sensu stricto) RHBSTW-01009
i6 contig#57 comp4 comp1236 comp10 comp95
Aligned 5'-end position on RHBSTW-01009 (nt)0252913396771889563249713434034924148822435370844408144541619
Aligned 3'-end position on RHBSTW-01009 (nt)2480123916468845432495763339922941440104348790443567045367004586750
Gap to right next contig (nt)490151255020137142634812491851444919
Gene assigned at gap right next contig16S rRNA16S rRNA16S rRNALrp/AsnC family transcriptional regulator23S rRNA23S rRNA23S rRNA5S rRNA5S rRNA
23S rRNAtRNA-GlutRNA-IleDMT family transportertRNA-GlutRNA-GlutRNA-Ala23S rRNA23S rRNA
23S rRNAtRNA-Ala16S rRNA16S rRNAtRNA-IletRNA-AlatRNA-Glu
23S rRNA16S rRNAtRNA-Ile16S rRNA
16S rRNA

comp: complementary.

 Table 5 

Genes in the genome of E. asburiae i6 strain that may contribute to plant growth promotion.

genelocus_taggene productfunction
ipdCEAI6_22530indolepyruvate decarboxylasephytohormone synthesis
iaaHEAI6_10770Indole-3-acetyl-aspartic acid hydrolaseindole acetic acid (IAA) synthesis
budAEAI6_03130acetolactate decarboxylaseacetoin synthesis
budBEAI6_03140acetolactate synthase AlsSacetoin synthesis
budCEAI6_03150(S)-acetoin forming diacetyl reductase2,3-butanediol synthesis
-EAI6_31530acetolactate synthase small subunitacetoin synthesis
-EAI6_31540acetolactate synthase 3 large subunitacetoin synthesis
ilvBEAI6_35410acetolactate synthase large subunitacetoin synthesis
ilvN_2EAI6_35420acetolactate synthase small subunitacetoin synthesis
ilvGEAI6_42550acetolactate synthase 2 catalytic subunitacetoin synthesis
ilvMEAI6_42560acetolactate synthase 2 small subunitacetoin synthesis
gcdEAI6_31050quinoprotein glucose dehydrogenasegluconic acid synthesis

Genes involved in restriction and anti-restriction in the E. asburiae i6 genome

Interestingly, KlcA1 and KlcA2, two homologues of KlcA, which have been reported to function as an anti-type I restriction modification (RM) system in other species [32], were predicted on this genome and plasmid respectively, but only one type IV R enzyme was predicted to be encoded and no type I RM system (Table 6). This is in contrast to other E. asburiae strains ATCC 35953, FDAARGOS_892, and RHBSTW-01009, which encode two type I RM systems, one type IV restriction enzyme, and KlcA1 (Table 6). L1 also encodes a type I RM system, a type IV R enzyme, and other restriction systems, but no KlcA homologs. Because it is important for successful genetic manipulation to suppress restriction enzymes and allow bacteria to introduce foreign plasmids, the transformation efficiency of several plasmid DNAs was tested by electroporation using the i6 and ATCC 35953 strains (Fig. 1).

The results showed that the ATCC 35953 strain could not be transformed by the relatively large plasmid DNAs pCF430 (10 kb) [33] and pAK1001 (8.4 kb) [34], whereas i6 formed colonies, showing transformation efficiency at least four orders of magnitude higher than that of ATCC 35953 especially with pCF430 (Fig. 1). Using a relatively small plasmid DNA, pACYC177 (3.9 kb) [35], the ATCC 35953 strain also formed colonies on LB plates containing kanamycin, and the transformation efficiency of i6 was two orders of magnitude higher than that of ATCC 35953 (Fig. 1).

 Figure 1 

Transformation efficiency of E. asburiae strains i6, i6 ∆klcA1 (AK1603), i6 ∆klcA2 (AK1604), and ATCC 35953 (NBRC 109912 T) and E. cloacae ATCC 1304 (ENC) strain with pACYC177 (dark bar), pAK1001 (gray bar), and pCF430 (white bar). Transformation efficiency (colony count/µg plasmid DNA) was calculated by counting colonies after electroporation, recovery in SOB for 1 h, and overnight growth at 37°C on LB (Lennox) plates containing kanamycin (pACYC177 and pAK1001) or tetracycline (pCF430).

J Genomics Image
 Table 6 

Highly identical (>90%) ortholog list of restriction enzymes and antirestriction proteins among E. asburiae strains.

Enterobacter asburiae (sensu stricto) strain
i6L1ATCC 35953FDAARGOS_892RHBSTW-01009
gene productlocus_tag
antirestriction protein (KlcA1)EAI6_00720NDACJ69_22495I6G49_09840HV349_19355
antirestriction protein (KlcA2)EAI6_40610NDNDNDND
type IV restriction enzymeEAI6_32920NDNDNDND
type IV restriction enzymeNDNDACJ69_09575I6G49_13575ND
PD-(D/E)XK nuclease superfamily proteinNDNDACJ69_09580I6G49_13570ND
type I restriction enzyme, R subunitNDNDACJ69_09595I6G49_13555HV349_13365
type I restriction enzyme, S subunitNDNDACJ69_09610I6G49_13540HV349_13375*
type I restriction enzyme M proteinNDNDACJ69_09615I6G49_13535HV349_13390
type II restriction enzyme M proteinNDNDACJ69_21260I6G49_08600HV349_10105/HV349_12830
type I restriction enzyme M proteinNDNDACJ69_22055I6G49_09395ND
type I restriction enzyme, S subunitNDNDACJ69_22060I6G49_09400ND
type I restriction enzyme, R subunitNDNDACJ69_22070I6G49_09410ND
type I restriction enzyme, R subunitNDNDNDNDHV349_19450
type I restriction enzyme, S subunitNDNDNDNDHV349_13375
type I restriction enzyme M proteinNDNDNDNDHV349_13380
5-methylcytosine-specific restriction enzyme BNDDI57_15815NDNDND
5-methylcytosine-specific restriction enzyme subunit McrCNDDI57_15820NDNDND
putative restriction endonuclease, HNH endonucleaseNDDI57_15875NDNDND
restriction methylaseNDDI57_15880NDNDND
Type IV restriction system proteinNDDI57_15975NDNDND
type I restriction enzyme, R subunitNDDI57_15980NDNDND
type I restriction enzyme, S subunitNDDI57_15985NDNDND
type I restriction enzyme M proteinNDDI57_15990NDNDND

ND: not detected. *Identity was less than 90%. Type I restriction enzyme, R subunit was shown in bold.

Application of the one-step gene inactivation method to the i6 strain

Unlike the ATCC 35953 strain, the i6 strain was able to accept foreign DNA with high efficiency (Fig. 1), so I attempted to apply the one-step gene inactivation method [36]. Prior to this, RSFRedTER-Sp was constructed by replacing the CmR marker of the RSFRedTER plasmid [18], which can express λ Red recombinase, with the SpR marker, so that the CmR marker could be used in addition to the KmR and TcR markers in the i6 strain. (Ap-resistant pKD46 was not used because two β-lactamase genes were predicted in the i6 genome.) Gene deletions in klcA1, klcA2, and dgeC were constructed in the i6 strain expressing λ Red recombinase from RSFRedTER-Sp using TcR, CmR, and KmR markers, respectively. Similarly, a deletion in rcsB and a KmR insertion behind wzc were made in the i6 strain using CmR and KmR markers, respectively. The recombinants of interest grew with high colony formation rates (Table 7). Colony PCR in the junction region of the recombination site confirmed that the desired recombinants were constructed with high accuracy (Table 7). The ∆klcA1 andklcA2 strains were included in the transformation efficiency analysis because KlcA was originally reported as a restriction enzyme inhibitor in other bacteria, such as E. coli and Klebsiella pneumoniae [32], and may have contributed to the high transformation efficiency of the plasmid DNA of strain i6. However, contrary to expectations, these deletions increased rather than decreased the transformation efficiency of some plasmid DNA (Fig. 1). Yet, this result may have been confounded by the fact that the type I RM system, a potential target of KlcAs, is not encoded in the i6 genome. In conclusion, this report successfully applied the one-step gene inactivation method to the newly identified E. asburiae strain i6, whose transformation efficiency was much higher than that of ATCC 35953.

Among the ECC, E. cloacae and E. hormaechei are the most frequently isolated species in clinical infections, especially in immunocompromised patients and those admitted to intensive care units [37], whereas E. asburiae had lower survival rates against serum than E. cloacae, E. hormaechei, and E. ludwigii isolates [4]. In this report, the E. asburiae i6 genome unexpectedly exhibited most, if not all, of the genetic features of PGPB previously reported in Enterobater sp. J49. This was also true for the genome closest to i6, ATCC 35953, FDAARGOS_892, L1, RHBSTW-01009 (data not shown). Thus, E. asburiae is a preferred PGPB candidate in the ECC. Even if the current form of the i6 strain is not optimal in promoting plant growth, it would be advantageous and useful as a starting material for testing and fine-tuning any beneficial and detrimental traits. This is because the i6 strain, with all its potential as a PGPB, can be successfully genetically engineered. First of all, potential risk factors such as genes responsible for antibiotic resistance, biofilm formation, and some of the type VI secretion apparatus/effectors, etc. should be deleted by the one-step gene disruption method or its derivative, scarless genome editing [38]. Such a risk negative i6 derivative strain could then be chemically mutagenized to express PGPB factors at optimal levels. Alternatively, or in addition, the ability of strain i6 to be a PGPB can be further enhanced by incorporating plasmid clones harboring other PGPB factors, even from different species, but only as genetically modified organisms. On the other hand, the amenability of this strain must also be useful to identify the relevant genes by deleting candidate genes from this genome that could produce a signal molecule that activates RcsB in Salmonella or even within the i6 strains. The recently developed conjugation-mediated versatile site-specific single-copy luciferase fusion system [39], which is broadly applicable to Gram-negative bacteria, should also be helpful in detecting intercellular interactions between the i6 strain and Salmonella, etc., and even within the i6 strains.

Nucleotide Sequence Accession Number

The draft genome sequence of E. asburiae i6 has been deposited in the DDBJ/EMBL/GenBank databases under accession numbers BTPF01000001 to BTPF01000028. The raw sequence reads were deposited in DDBJ under BioProject number PRJDB16388 and BioSample number SAMD00634955.

 Table 7 

Recombination efficiency and accuracy for deletion and insertion of targeted genes in the i6 strain

Target geneChromosome or plasmidDeletion or insertionMarkerNumber of recombinantsTargeting success rate (%)
klcA1chromosomedeletionTcR289100 (8/8)
klcA2plasmiddeletionCmR143100 (8/8)
dgeCchromosomedeletionKmR48100 (8/8)
wzcchromosomeinsertionKmR6100 (6/6)

Acknowledgements

I thank Dr. Yoshihiko Hara for providing the plasmid RSFRedTER.

Competing Interests

The authors have declared that no competing interest exists.

References

1. Paauw A, Caspers MP, Schuren FH, Leverstein-van Hall MA, Deletoile A, Montijn RC. et al. Genomic diversity within the Enterobacter cloacae complex. PLoS One. 2008;3:e3018

2. Humann JL, Wildung M, Pouchnik D, Bates AA, Drew JC, Zipperer UN. et al. Complete genome of the switchgrass endophyte Enterobacter clocace P101. Stand Genomic Sci. 2014;9:726-34

3. Brady C, Cleenwerck I, Venter S, Coutinho T, De Vos P. Taxonomic evaluation of the genus Enterobacter based on multilocus sequence analysis (MLSA): proposal to reclassify E nimipressuralis and E amnigenus into Lelliottia gen nov as Lelliottia nimipressuralis comb nov and Lelliottia amnigena comb nov, respectively, E gergoviae and E pyrinus into Pluralibacter gen nov as Pluralibacter gergoviae comb nov and Pluralibacter pyrinus comb nov, respectively, E cowanii, E radicincitans, E oryzae and E arachidis into Kosakonia gen nov as Kosakonia cowanii comb nov, Kosakonia radicincitans comb nov, Kosakonia oryzae comb nov and Kosakonia arachidis comb nov, respectively, and E turicensis, E helveticus and E pulveris into Cronobacter as Cronobacter zurichensis nom nov, Cronobacter helveticus comb nov and Cronobacter pulveris comb nov, respectively, and emended description of the genera Enterobacter and Cronobacter. Syst Appl Microbiol. 2013;36:309-19

4. Ganbold M, Seo J, Wi YM, Kwon KT, Ko KS. Species identification, antibiotic resistance, and virulence in Enterobacter cloacae complex clinical isolates from South Korea. Front Microbiol. 2023;14:1122691

5. Yaish MW. Draft Genome Sequence of Endophytic Bacterium Enterobacter asburiae PDA134, Isolated from Date Palm (Phoenix dactylifera L.) Roots. Genome Announc. 2016; 4

6. Araujo WL, Marcon J, Maccheroni W Jr, Van Elsas JD, Van Vuurde JW, Azevedo JL. Diversity of endophytic bacterial populations and their interaction with Xylella fastidiosa in citrus plants. Appl Environ Microbiol. 2002;68:4906-14

7. Hinton DM, Bacon CW. Enterobacter cloacae is an endophytic symbiont of corn. Mycopathologia. 1995;129:117-25

8. Asis CA 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

9. Zhang G, Sun Y, Sheng H, Li H, Liu X. Effects of the inoculations using bacteria producing ACC deaminase on ethylene metabolism and growth of wheat grown under different soil water contents. Plant Physiol Biochem. 2018;125:178-84

10. Li G, Hu Z, Zeng P, Zhu B, Wu L. Whole genome sequence of Enterobacter ludwigii type strain EN-119T, isolated from clinical specimens. FEMS Microbiol Lett. 2015 362

11. Bi C, Zhang X, Ingram LO, Preston JF. Genetic engineering of Enterobacter asburiae strain JDR-1 for efficient production of ethanol from hemicellulose hydrolysates. Appl Environ Microbiol. 2009;75:5743-9

12. Edoamodu CE, Nwodo UU. Marine sediment derived bacteria Enterobacter asburiae ES1 and Enterobacter sp. Kamsi produce laccase with high dephenolisation potentials. Prep Biochem Biotechnol. 2022;52:748-61

13. Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci U S A. 2000;97:6640-5

14. Ellermeier CD, Janakiraman A, Slauch JM. Construction of targeted single copy lac fusions using lambda Red and FLP-mediated site-specific recombination in bacteria. Gene. 2002;290:153-61

15. Uzzau S, Figueroa-Bossi N, Rubino S, Bossi L. Epitope tagging of chromosomal genes in Salmonella. Proc Natl Acad Sci U S A. 2001;98:15264-9

16. Winfield MD, Latifi T, Groisman EA. Transcriptional regulation of the 4-amino-4-deoxy-L-arabinose biosynthetic genes in Yersinia pestis. J Biol Chem. 2005;280:14765-72

17. Mitrophanov AY, Jewett MW, Hadley TJ, Groisman EA. Evolution and dynamics of regulatory architectures controlling polymyxin B resistance in enteric bacteria. PLoS Genet. 2008;4:e1000233

18. Katashkina JI, Hara Y, Golubeva LI, Andreeva IG, Kuvaeva TM, Mashko SV. Use of the lambda Red-recombineering method for genetic engineering of Pantoea ananatis. BMC Mol Biol. 2009;10:34

19. Derbise A, Lesic B, Dacheux D, Ghigo JM, Carniel E. A rapid and simple method for inactivating chromosomal genes in Yersinia. FEMS Immunol Med Microbiol. 2003;38:113-6

20. Lesic B, Rahme LG. Use of the lambda Red recombinase system to rapidly generate mutants in Pseudomonas aeruginosa. BMC Mol Biol. 2008;9:20

21. Marchesi JR, Sato T, Weightman AJ, Martin TA, Fry JC, Hiom SJ. et al. Design and evaluation of useful bacterium-specific PCR primers that amplify genes coding for bacterial 16S rRNA. Appl Environ Microbiol. 1998;64:795-9

22. Hoffmann H, Roggenkamp A. Population genetics of the nomenspecies Enterobacter cloacae. Appl Environ Microbiol. 2003;69:5306-18

23. Kato A, Latifi T, Groisman EA. Closing the loop: the PmrA/PmrB two-component system negatively controls expression of its posttranscriptional activator PmrD. Proc Natl Acad Sci U S A. 2003;100:4706-11

24. Richter M, Rossello-Mora R. Shifting the genomic gold standard for the prokaryotic species definition. Proc Natl Acad Sci U S A. 2009;106:19126-31

25. Parks DH, Chuvochina M, Rinke C, Mussig AJ, Chaumeil PA, Hugenholtz P. GTDB: an ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy. Nucleic Acids Res. 2022;50:D785-D94

26. Glick BR. Plant growth-promoting bacteria: mechanisms and applications. Scientifica (Cairo). 2012;2012:963401

27. Coulson TJ, Patten CL. Complete Genome Sequence of Enterobacter cloacae UW5, a Rhizobacterium Capable of High Levels of Indole-3-Acetic Acid Production. Genome Announc. 2015 3

28. Kim B, Park AR, Song CW, Song H, Kim JC. Biological Control Efficacy and Action Mechanism of Klebsiella pneumoniae JCK-2201 Producing Meso-2,3-Butanediol Against Tomato Bacterial Wilt. Front Microbiol. 2022;13:914589

29. Chandarana KA, Amaresan N. Predation pressure regulates plant growth promoting (PGP) attributes of bacterial species. J Appl Microbiol. 2023 134

30. Luduena LM, Anzuay MS, Angelini JG, McIntosh M, Becker A, Rupp O. et al. Genome sequence of the endophytic strain Enterobacter sp. J49, a potential biofertilizer for peanut and maize. Genomics. 2019;111:913-20

31. Ryu CM, Farag MA, Hu CH, Reddy MS, Wei HX, Pare PW, Kloepper JW. Bacterial volatiles promote growth in Arabidopsis. Proc Natl Acad Sci U S A. 2003;100:4927-32

32. Serfiotis-Mitsa D, Herbert AP, Roberts GA, Soares DC, White JH, Blakely GW. et al. The structure of the KlcA and ArdB proteins reveals a novel fold and antirestriction activity against Type I DNA restriction systems in vivo but not in vitro. Nucleic Acids Res. 2010;38:1723-37

33. Newman JR, Fuqua C. Broad-host-range expression vectors that carry the L-arabinose-inducible Escherichia coli araBAD promoter and the araC regulator. Gene. 1999;227:197-203

34. Kato A. In vivo cloning of large chromosomal segments into a BAC derivative by generalized transduction and recombineering in Salmonella enterica. J Gen Appl Microbiol. 2016;62:225-32

35. Schottel JL, Bibb MJ, Cohen SN. Cloning and expression in streptomyces lividans of antibiotic resistance genes derived from Escherichia coli. J Bacteriol. 1981;146:360-8

36. Datsenko KA, Wanner BL. One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products. Proc Natl Acad Sci USA. 2000;97:6640-5

37. Davin-Regli A, Lavigne JP, Pages JM. Enterobacter spp.: Update on Taxonomy, Clinical Aspects, and Emerging Antimicrobial Resistance. Clin Microbiol Rev. 2019; 32

38. Fels U, Gevaert K, Van Damme P. Bacterial Genetic Engineering by Means of Recombineering for Reverse Genetics. Front Microbiol. 2020;11:548410

39. Kato A. Development of conjugation-mediated versatile site-specific single-copy luciferase fusion system. J Gen Appl Microbiol. 2023

40. Fields PI, Groisman EA, Heffron F. A Salmonella locus that controls resistance to microbicidal proteins from phagocytic cells. Science. 1989;243:1059-62

41. Cherepanov PP, Wackernagel W. Gene disruption in Escherichia coli: TcR and KmR cassettes with the option of Flp-catalyzed excision of the antibiotic-resistance determinant. Gene. 1995;158:9-14

42. Prjibelski A. et al. Using SPAdes De Novo Assembler. Curr Protoc Bioinformatics. 2020;70:e102

Author contact

Corresponding address Corresponding author: Akinori Kato, Department of Advanced Bioscience, Graduate School of Agriculture, Kindai University, 3327-204 Nakamachi, Nara, Nara 631-8505, Japan.

Tel: +81 742 43 7213 Fax: +81 742 43 8976 E-mail: kato_akindai.ac.jp


Received 2023-10-20
Accepted 2024-1-8
Published 2024-1-20


Citation styles

APA
Kato, A. (2024). Isolation and draft genome sequence of Enterobacter asburiae strain i6 amenable to genetic manipulation. Journal of Genomics, 12, 26-34. https://doi.org/10.7150/jgen.91337.

ACS
Kato, A. Isolation and draft genome sequence of Enterobacter asburiae strain i6 amenable to genetic manipulation. J. Genomics 2024, 12, 26-34. DOI: 10.7150/jgen.91337.

NLM
Kato A. Isolation and draft genome sequence of Enterobacter asburiae strain i6 amenable to genetic manipulation. J Genomics 2024; 12:26-34. doi:10.7150/jgen.91337. https://www.jgenomics.com/v12p0026.htm

CSE
Kato A. 2024. Isolation and draft genome sequence of Enterobacter asburiae strain i6 amenable to genetic manipulation. J Genomics. 12:26-34.

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
Popup Image