A�nibrenneria salicis gen. nov. sp. nov. isolated from Salix matsudana bark canker

L3-3HA T , a Gram-negative-staining, facultatively anaerobic, motile bacterial strain, was isolated from the symptomatic bark of Salix matsudana canker in China. 16S rRNA gene analysis revealed that the novel strain shares the highest sequence similarity with Brenneria goodwinii FRB141 T (95.5 %). In phylogenetic trees based on four housekeeping genes (gyrB, rpoB, atpD and infB) and the 16S rRNA gene sequence, the novel strain formed a separate branch from the �ve genera of the family Pectobacteriaceae (Lonsdalea, Brenneria, Dickeya, Pectobacterium and Sodalis), suggesting that the novel strain should belong to a novel species of a novel genus within the family Pectobacteriaceae. The result was also supported by phylogenomics, amino acid identity and average nucleotide identity. The major fatty acids were C 14:0 , C 16:0 , C 17:0 cyclo, and C 19:0 cyclo (cid:0) 8c. Genome analysis showed that the novel strain has a large genome (5.89 Mb) with 5,052 coding genes, including 181 virulence genes by searching the pathogen-host interactions database (PHI-base), indicating that the novel strain is a potential pathogen of plants and animals. Based on phenotypic and genotypic characteristics, the L3-3HA T strain represents a novel species of a novel genus in the Pectobacteriaceae family, for which the name A�nibrenneria salicis gen nov. sp. nov. is proposed. The strain type is L3-3HA T (= CFCC 15588 T = LMG 31209 T ).

A novel strain, L3-3HA T , was isolated from the symptomatic bark of Salix matsudana canker.The novel strain belongs to the family Pectobacteriaceae based on 16S rRNA gene sequence analysis.In the family Pectobacteriaceae, many species are associated with plant disease, such as Pectobacterium carotovorum, Brenneria salicis, and Lonsdalea populi.To explore the possible function of the novel strain, genome comparison among the novel and ve reference strains was performed, and virulence genes and carbohydrate-active enzymes of the novel strain were predicted by searching the carbohydrate-active enzyme database (CAZy) and pathogen-host interactions database (PHI-base).To con rm the taxonomy status of the novel strain, polyphasic approaches, including 16S rRNA gene sequence analysis, multilocus sequence analysis (MLSA), phylogenomics, average nucleotide identity (ANI), amino acid identity (AAI), digital DNA-DNA hybridization (dDDH) and phenotypic characterization, were performed in the present study.

Materials And Methods
Isolation and cultivation of the novel strain The novel strain was isolated from the bark samples of Salix matsudana canker as described by Li et al. (2017a).Samples were sterilized in 4 % (v/v) NaClO solution for 3 min and immersed in 75 % alcohol for 30 s.After rinsing three times with sterile water, the sterilized bark samples were cut into 2×2-mm pieces and soaked in sterile water for 30 min.The solutions were diluted with sterile water to a series of 1:10 -1 to 10 -4 and spread on tryptic soy agar plates (TSA; Difco).Single colonies were obtained after incubating at 30 ℃ for 48 h.The strain was deposited at the Laboratorium voor Microbiologie, Universiteit Gent (LMG 31209 T ) and China Forestry Culture Collection Center (CFCC 15588 T ).
Phylogenetic analysis of the 16S rRNA gene was performed according to Li et al. (2016).35 almost complete (1344 bp) 16S rRNA gene sequences were aligned and overhangs trimmed using CLUSTAL W in MEGA 7.0 software (Kumar et al. 2016).The trees were reconstructed using the neighbor-joining and maximum likelihood methods and were evaluated with 1000 replicates.
MLSA based on four housekeeping gene partial sequences (gyrB, 745 bp; rpoB, 637 bp; atpD, 642 bp; and infB, 615 bp) was performed on the novel strain and reference type strains (Brady et al. 2013).The reference strain sequences of genera Lonsdalea, Brenneria, Dickeya, Pectobacterium, Sodalis, Izhakiella, Pantoea and Mixta were downloaded from GenBank (Table S1), and the sequences of the novel strain were retrieved from its genome sequence data.Phylogenetic trees based on the maximum likelihood and neighbor-joining method were reconstructed by Mega 7.0.The bootstrap method with 1,000 replicates was performed.

Whole-genome sequencing
Genome sequencing of the novel strain was performed using Illumina Miseq sequencing strategies as described by Li et al. (2016).The library was constructed with 300-to 500-bp DNA fragmentation using the TruSeq ™ Sample Prep Kit.Sequence assembly was performed using SOAPdenovo v2.04.Coding sequences were predicted using Glimmer 3.02 software.Functional annotation of the protein sequences was performed against the nonredundant GenBank database using BLASTp.
An almost full-length 16S rRNA gene sequence obtained by conventional Sanger sequencing was compared with the sequence extracted from whole-genome assembly.Two sequences (1,409 bp and 276 bp) from the genome sequence shared 100 and 99.26 % sequence similarity with the 16S rRNA gene sequence obtained by conventional Sanger sequencing.However, the 276-bp sequence may not be a 16S rRNA gene sequence, as con rmed by two methods.First, the DNA used for genome sequencing was PCR (polymerase chain reaction) ampli ed and sequenced by conventional Sanger sequencing, and no doublet and impure peaks were found.Second, the K-mer and GC-depth gures were regular and normal in the genome sequence analysis (Fig. S3 and S4).
As recommended by Chun et al. (2018), phylogenomics, DNA-DNA hybridization and ANIb analysis were used for the taxonomy of the novel strain.Phylogenomics was performed on the novel and its closely related reference species.The genomic coding sequences were clustered using CD-HIT v4.6.1 (c 0.5 -n 3 -p 1 -T 4 -g 1 -d 0 -s 0.7 -aL 0.7 -aS 0.7) (Fu et al. 2012).Six hundred eighteen single-copy genes were selected using Perl scripts and aligned using Muscle 3.8.31(Edgar 2004).A phylogenetic ML tree was constructed based on the resulting multiple alignments using PhyML V3.3 (Guindon et al. 2010) (parameters: phyml -i temp. phy -d aa -b -4 -m LG -a e -c 4 -t e >ml.log).The level of support for internal branches was assessed using the Bayesian test in PhyML.
Average nucleotide identity (ANI) is a similarity measure between two genome sequences that can be used to replace DNA-DNA hybridization (Richter et al. 2009;Goris et al. 2007).The OrthoANI values between the L3-3HA T strain and ve reference type strains of the ve genera of the family Pectobacteriaceae and other reference strains were calculated using the ANI calculator tool of the EzBioCloud website (https://www.ezbiocloud.net/tools/ani).dDDH was conducted using the genome-to-genome distance calculator (GGDC) and the free web service at http://ggdc.dsmz.de/based on genome sequences (Meier-Kolthoff et al. 2013).Amino acid identity (AAI) analysis among the novel and reference strains of genera Lonsdalea, Brenneria, Dickeya, Pectobacterium, and Sodalis was performed as described by Konstantinidis and Tiedje (2005).
Genome comparison among the novel strain and D. chrysanthemi NCPPB 402 T , S. praecaptivus HS T , P.
To understand the potential function of the novel strain, virulence genes and carbohydrate-active enzymes were predicted by searching the Carbohydrate-Active enZYmes Database (CAZy) (Lombard et al. 2014) and Pathogen-Host Interactions database (PHI-base) (Urban et al. 2020) using DIAMOND (Buch nk et al. 2015) based on the genome amino acid sequence of the novel strain (evalue≤1e-5, identity≥50 %, coverage ≥40 %), respectively.PHI-base is a web-accessible database that catalogs experimentally veri ed pathogenicity, virulence and effector genes from fungal, oomycete and bacterial pathogens (Urban et al. 2020).CAZy is a specialist database dedicated to the display and analysis of genomic, structural and biochemical information on carbohydrate-active enzymes (CAZymes) (Lombard et al. 2014).

Morphological, physiological and biochemical analyses
Colony morphology, motility and growth conditions (temperature and NaCl tolerance) experiments were performed according to the description of Li et al. (2017b).Cell morphology was examined by transmission electron microscopy with cells grown on TSA for 24 h at 30 ℃.The motility of the novel isolate was examined by observing cells developed on wet mounts using a phase-contrast microscope (TMS-F; Nikon).
The novel strain was inoculated in tryptic soy broth (TSB; Difco) and incubated at different temperatures (4, 10, 30, 37 and 50 ℃) as described by Li et al. (2015), different pH values (4.0-11.0, at intervals of 1.0) and various tolerances to NaCl (0 %-10 % [w/v], at intervals of 1 %).As reported by Gomori (1955), the pH values were adjusted using the following buffers: citric acid/sodium citrate (pH 4.0-5.5),Na 2 HPO 4 /NaH 2 PO 4 (pH 6.0-8.0),Na 2 CO 3 /NaHCO 3 (pH 8.0-9.5) and Na 2 HPO 4 /NaOH (pH 10.0-11.0).Anaerobic growth was examined after the novel strain was incubated on TSA plates at 30 ℃ for 1 week in anaerobic jars lled with N 2 /CO 2 (80: 20, v/v) (Gerhardt et al. 1981).Catalase and oxidase activities were measured as recommended by Smibert and Krieg (1994).The oxidase activity of the novel strain was examined by determining whether a color change to purple occurred after using 1 % (w/v) tetramethyl-p-phenylenediamine.Gram staining was tested according to the description of Jenkins (2003).Other biochemical and physiological tests were performed using API 20E, API 20NE, and API 50CHB/E kits (bioMérieux) and the Biolog GN 2 MicroPlate according to the manufacturer's instructions.The strains tested in these analyses were cultured on TSA at 30 ℃ for 24 h and then were harvested during the exponential phase.

Chemotaxonomic analyses
The strains were cultured on TSA at 30 ℃ for 24 h and then were harvested during the exponential phase.
Fatty acids were analyzed using the Microbial Identi cation System, Sherlock version 6.0 (Sasser 1990).
Pathogenicity determination 4 year-old poplar trees were chosen for inoculation at the poplar nursery in Puyang City, Henan Province.The strain was cultured on nutrient agar plates at 30 °C for 24 hours, then were washed and centrifuged to pellet the cells.The collected cells were diluted in SDW to 10 8 CFU/ml.After the bark was surface sterilized with 75 % ethanol, 1 cm-long cross-wounds were cut into the bark with a sterilized knife.A total of 100 μl of bacterial suspension was delivered to the inoculation point.Negative control plants were inoculated with sterile distilled water.The results of pathogenicity tests were recorded 7, 14, and 21 days after inoculation.

Results And Discussion
Morphological and phenotypic characteristics Transmission electron microscopy showed that the cells were short rods or globose, with peritrichous agella and 0.5-1.0 × 0.8-2.5 μm in size (Fig. S5).The colonies were white, circular, and approximately 1-2 mm in diameter after 48 h of growth at 30 ℃ on TSA.The strain is a motile, facultative anaerobic bacterium with positive catalase and negative oxidase activities.The detailed differentiated biochemical and physiological characteristics of the novel strain and close genera are listed in Tables 1.

Phylogenetic analysis
The 16S rRNA gene sequencing of novel strain L3-3HA T shares the highest sequence similarity with Brenneriagoodwinii FRB141 T (95.5 %) and lower sequence similarity values (95.4 %) with other validly published species.In ML and NJ phylogenetic trees based on 16S rRNA gene sequencing (Fig. 1, Fig. S1), the novel strain formed an dependent lineage, separated from those of ve Pectobacteriaceae genera and other family reference species of Izhakiella, Pantoea and Mixta.And in the ML and NJ phylogenetic trees based on MLSA data (Fig. 2, Fig. S2), the novel strain clustered into a distinct branch from other genera in the present study, suggesting that the novel strain should represent a novel species in a novel genus within the Pectobacteriaceae family.

Genomic characteristics
The sequencing depth of coverage was 130 ´.The genome of strain L3-3HA T is 5.89 Mb in size across 47 contigs (N50 is 322317 bp), including 5,052 coding genes, 10 rRNAs, 74 tRNAs, 9 other RNAs and 197 pseudogenes.The DNA G+C content is 57.1 %.Also, the strain L3-3HA T formed a distinct branch from the ve genera within the Pectobacteriaceae family clades in the phylogenomic tree (Fig. 3), demonstrating that the novel strain represents a novel species in a novel genus within the family Pectobacteriaceae.
The novel strain showed 72.43-77.56% ANI values with the reference strains of the family Pectobacteriaceae and Erwiniaceae, and the detailed data were shown in Table 2.These values were all lower than the proposed species boundary ANI cut-off (95-96 %) (Fu et al. 2012).The dDDH values between L3-3HA T strain and the same reference strains used in ANI analysis were 13 %-22 %, lower than the 70 % cut-off for species determination (Table 2) (Wayne et al. 1987).
The AAI analysis results revealed that the ve reference type strains of the genera Lonsdalea, Brenneria, Dickeya, Pectobacterium, Sodalis shared 64 %-81 % AAI values with each other.The novel strain showed 64.2 %-74.9 % AAI values with these ve reference type strains (Table S3), lower than the AAI values among the ve reference type strains.These results supported that the novel isolate represents a novel species in a novel genus in family Pectobacteriaceae.
Genome comparison shows that the novel strain (5.89 Mb) has large genome size , comprises 5,052 coding genes, 1,907, 1,617, 990, and 1,045 more genes than the reference strains L. quercina ATCC 29281 T , B. salicis DSM 30166 T , P. carotovorum LMG 2404 T and D. chrysanthemi NCPPB 402 T .L. quercina and B. salicis are pathogens of Quercus spp.and Salix spp., respectively (Brady et al. 2012;Sakamoto et al. 1999), and they have a small genome size (3.9-4.0Mb) with a single host compared with P. carotovorum and D. chrysanthemi.The plant pathogens P. carotovorum and D. chrysanthemi, with a diverse host range (Toth et al. 2003), have a large genome (4.8 Mb in size), which is seemingly important to adopt diverse hosts for the bacterial pathogen.Thus, the novel strains may be a multihost bacterial plant pathogen.
In the present study, 181 virulence genes were found in the genome of the novel strain, including 15 "loss of pathogenicity" genes, 155 "reduced virulence" genes and 15 "increased virulence" genes.The results indicated that the novel strain is a potential pathogen of plants, although it is not a pathogen of the poplar tree.And in total, 145 carbohydrate-active enzymes were found from the genome sequence of the novel strain.Eleven enzymes were found to be associated with degradation, including 2 cellulases 8 chitinases and 1 pectin lyase, which are important enzymes for the pathogens of plants.

Chemotaxonomic characteristics
The major fatty acids detected in the novel strain were C 14:0 , C 16:0 , C 17:0 cyclo and C 19:0 cyclo 8c.The percentages of the peak area for fatty acids of the novel and reference strains are listed in Table S2.The amount of C 17:0 cyclo and C 19:0 cyclo 8c can be used to distinguish the novel strain from the ve reference strains used in the present study.

Pathogenicity test
The results of pathogenicity experiment showed that the novel strain did not show pathogenicity to poplar trees.

Taxonomic conclusion
In the trees based on the 16S rRNA gene and MLSA data, the strains of the genus Brenneria were clustered into two separate clades, indicating that the taxonomy of the Brenneria species may be polyphyletic.Although the Brenneria strains clustered together in the phylogenomic trees, the Brenneria species were still divided into two clades.Therefore, we speculated that B. alni, B. corticis, B. goodwinnii, B. nigri uens, and B. populi may represent a novel genus within the family Pectobacteriaceae.
The 16S rRNA gene and MLSA analysis revealed that the novel strain belongs to a novel species in a novel genus within the family Pectobacteriaceae, which is also supported by AAI, fatty acid and phyologenomic analyses.Moreover, many physiological and biochemical characteristics can be used to distinguish the novel strain from close genera of the family Pectobacteriaceae (detailed in Table 2).Based on those data, the L3-3HA T strain represents a novel species of a novel genus in the Pectobacteriaceae family, for which the name A nibrenneria salicis gen nov.sp.nov. is proposed.
Description of A nibrenneria gen.nov.
Gram-negative, facultative anaerobic, motile, oxidase-negative, and catalase-positive.Growth occurs between 4 ℃ and 50 ℃ and at pH 6-8 and 0 %-5 % salinity.The cells are also positive for β-galactosidase activity and negative for the activities of urease and gelatinase and the production of H 2 S and indole.Citrate is not utilized, and nitrates are not reduced to nitrites.The main fatty acids are C 14:0 , C 16:0 , C 17:0 cyclo and C 19:0 cyclo 8c.The DNA G+C content is 57.1 %, and the type species is A nibrenneria salicis.A nibrenneria salicis gen.nov.sp.nov. is a member of the family Pectobacteriaceae.
Description of A nibrenneria salicis sp.nov.
A nibrenneria salicis (sa'li.cis.L. gen.n. salicis of Salix matsudana) The cells are Gram-negative, facultative anaerobic, motile, short rod or globose with peritrichous agella.The colonies are circular, convex with complete margins.The cells show negative activities for arginine dihydrolase, lysine decarboxylase, ornithine decarboxylase, tryptophan deaminase, urease, and gelatinase and are negative for the production of H 2 S and indole.The cells do not utilize citrate and do not reduce nitrates to nitrites.The cells are positive for acid production from d-arabinose, d-ribose, d-galactose, d-glucose, d-fructose, d-mannose, methyl α-d-mannopyranoside, methyl nacetylglucosamine,amygdalin,arbutin,esculin ferric citrate,salicin,dmelibiose,gentiobiose,.

Figures Figure 1
Figures

Figure 2 Maximum
Figure 2