BACKGROUND
the advent of next generation sequencing technologies with their low cost, high throughput and rapid speed has strongly influenced modern taxonomy.
many taxonomists therefore turn towards integrating genome sequence analysis for formal proposals of new bacterial taxa or reclassification of old taxa .
early classification methods were based solely on phenotypic characteristics, but the introduction of dna-dna hybridization   <cit>  and 16s ribosomal rna  gene phylogeny  <cit>  in taxonomy changed the species definition .
however, due to the high level of conservation of the 16s rrna gene  <cit>  and the labour intensive and error prone nature of ddh experiments, there has been a constant demand for advanced methods at complete genome level.
average nucleotide identity  , multi locus sequence analysis of conserved genes  <cit>  and phylogenetic analysis of core genomes  <cit>  are some of the suggested methods.
brachyspira species are ggram-negative, anaerobic but oxygen-tolerant spirochaetes that colonize the intestines of many mammals and birds  <cit> .
it forms a clade within phylum spirochaetes that has undergone repeated taxonomic changes over the past few decades.
initially the isolates colonizing pigs were allocated to genus treponema  <cit> , however they were later transferred to genus serpula and then serpulina  <cit> .
in  <dig>  the genera serpulina and brachyspira were unified leading to the reclassification of serpulina species  <cit> .
genus brachyspira now encompasses seven valid species: b. aalborgi, b. alvinipulli, b. hyodysenteriae, b. intermedia, b. innocens, b. murdochii and b. pilosicoli with b. aalborgi as a type species  <cit> .
these species vary in terms of their host range and enteropathogenicity.
brachyspira hyodysenteriae and b. pilosicoli are well known porcine pathogens  <cit> .
furthermore, b. pilosicoli can cause gastrointestinal disease and egg-production losses in chickens  <cit> , and it can also colonize the human intestine leading to human intestinal spirochaetosis  <cit> .
however, most cases in humans are associated with b. aalborgi, at least in industrialized countries  <cit> .
moreover, b. alvinipulli and b. intermedia cause diarrhoea in chickens  <cit> .
however, b. intermedia, is in general considered as a harmless commensal in pigs similar to b. innocens and b. murdochii  <cit> .
there are several proposed species without current standing in nomenclature, including “b.
suanatina”  <cit> .
some members of the genus brachyspira have low interspecies 16s rrna gene variability that does not always support species delineations, especially among the indole test-positive species b. hyodysenteriae, b, intermedia and b. suanatina that require additional methods for species confirmation  <cit> .
the proposed new species b. suanatina was first isolated from pigs with dysentery-like disease in sweden and denmark, and the enteropathogenicity of strain an4859/ <dig> was confirmed by experimental challenge in pigs  <cit> .
the phenotypic characterization based on the intensity of β-haemolysis and enzyme activities showed that these isolates were indistinguishable from b. hyodysenteriae, but detection by a species-specific pcr was unsuccessful  <cit> .
based on 16s rrna and nadh oxidase  gene similarities, these brachyspira isolates from pigs were similar to some isolates from mallards  <cit>  and were distinct from all currently recognized brachyspira species, suggesting that they were members of a putatively unrecognized species.
available data suggest that b. suanatina most likely is a rare genotype in both pigs and mallards, and it has so far only been isolated in northern europe.
identification of b. suanatina currently relies on nox gene sequencing and/or nox-rflp  <cit> .
the objective of the present study was to identify the taxonomical position of b. suanatina based on phenotypic characteristics and complete genome analyses.
we have used a “taxono-genomics” approach that is based on the polyphasic combination of phenotypic and genomic characteristics  <cit> .
methods
bacterial cultures and growth conditions
type and reference strains b. hyodysenteriae b78t , b204r , b. intermedia pws/at , b. pilosicoli p43/6/78t , b. innocens b256t , b. murdochii 56-150t , b. alvinipulli c1t  and b. aalborgi 513at  were originally obtained from the american type culture collection  and stored in liquid nitrogen in the brachyspira strain collection of the national veterinary institute .
brachyspira suanatina strains an4859/ <dig>  an1681:1/ <dig>  an2384/ <dig> and dk12570- <dig> from pigs, and an3949:2/ <dig> and an1418:2/ <dig> from mallards were isolated in sweden and denmark as previously described  <cit>  and stored in liquid nitrogen.
the bacteria were retrieved from storage and were grown anaerobically  on fastidious anaerobe agar   plates with 10 % equine blood at 42 °c for 72 h and subcultured once for 48 h before use.
cells for electron microscopy were harvested from broth cultures grown anaerobically on a shaker at 37 °c to mid log phase in brain heart infusion culture broth   with 10 % foetal calf serum   .
faa plate cultures and broth culture were used to test the capability of all b. suanatina strains to grow at  <dig>   <dig>   <dig> and 46 °c.
for the generation of growth curves, broth cultures of b. suanatina strains an4859/ <dig> and an3949:2/ <dig> and b. hyodysenteriae b78t were prepared in triplicates in 10 ml bhis at  <dig> cells/ml and grown at 37 °c.
population densities were determined every 4 h between  <dig> and 70 h of incubation by colony forming unit  counts obtained from serial dilutions on faa plates in duplicates and by spectrophotometry .
comparisons were made against an in-house standard curve for b. hyodysenteriae b78t.
purity was confirmed throughout the studies by phase contrast microscopy.
biochemical phenotype and substrate utilization
the api-zym test  was applied to all b. suanatina strains with b. hyodysenteriae b78t and pseudomonas aeruginosa  as control strains.
bacterial cells were harvested from faa plates incubated at 42 °c for 3 days, suspended in 2 ml sterile ddh2o to a density of 5– <dig> mcfarland and the api-zym test was performed as described by the manufacturer.
the test strips were incubated for 4 h in aerobic atmosphere.
brachyspira suanatina strains an4859/ <dig>  an3949:2/ <dig> and b. hyodysenteriae strain b78t were used for growth substrate utilization studies.
all six b. suanatina strains were tested for growth in d-mannitol because of diverging results obtained with strains an4859/ <dig> and an3949:2/ <dig>  potential substrates were added at  <dig>  wt/v in heart infusion culture broth   and autoclaved or sterile filtered.
the fcs concentration in the hi was selected at 7 %  from test runs at  <dig>   <dig> and 10 % fcs as described previously  <cit> .
triplicate flasks with 10 ml hi and 7 % fcs culture broth  were used for each isolate and substrate, and were inoculated with  <dig> cells from a stock solution made from 48 h old faa cultures.
cultures were incubated anaerobically at 37 °c on a shaking platform.
glucose was used as positive control and his without added substrate as negative control.
growth was determined spectrophotometrically  at  <dig>   <dig> and 9 days post inoculation for all cultures, and the purity, cell density and the condition of the cells were examined by phase contrast microscopy.
compounds producing ≤ <dig>  od <dig> were not considered to be growth substrates.
these cultures contained degenerated and/or non-motile cells as determined by phase contrast microscopy.
cultures producing ≥ <dig>  od <dig> values and moderate to profuse growth of motile spirochaetes on at least one occasion were considered to be growth substrates.
the following compounds were investigated: d-glucose, d-fructose, d-fucose, d-ribose, d-mannose, sucrose, lactate, d-xylose, d-trehalose, d-galactose, d-melobiose, d-cellobiose, l-arabinose, l-fucose, l-rhamnose, d-raffinose, d-galactosamine, n-acetyl-d-glucosamine, n-acetyl-d-galactosamine, mannitol, pyruvate, esculin, porcine gastric mucin, 50x mem amino acids, pectin, hyaluronic acid, glutathione, poly-galacturonic acid, cellulose, glycerol, succinate and glycogen.
electron microscopy
strain an4859/ <dig> was grown in bhis to mid-log phase, centrifuged at  <dig>  × g for 3 min and washed twice in 1 ml  <dig>  m sodium phosphate buffer at ph  <dig> .
the pellet was resuspended in  <dig>  ml each of sterile ddh2o and negatively stained with  <dig>  ml 2 % phosphotungistic acid at ph  <dig> .
the cells were dropped on grids and examined.
for scanning electron microscopy the cells were harvested and washed as described above, dropped on grids and coated with gold-palladium with standard procedures.
the samples were examined with a jsm  <dig> scanning electron microscope .
ten cells or more were examined to determine cell dimensions, the location of attachment sites of periplasmatic flagella and the number of flagella per cell.
dna-dna hybridization
the ddh of b. suanatina strain an4859/ <dig> was performed with all the seven recognized brachyspira species.
cells were grown to mid-log phase in bhis, checked for purity and centrifuged at 4 °c at  <dig>  g followed by washing in pbs twice before mixing with 1:1 v/v iso-propanol .
dna isolation and ddh were performed by dsmz gmbh deutsche sammlung von mikroorganismen und zellkulturen gmbh dna .
dna was extracted using a french pressure cell  and was purified by chromatography on hydroxyapatite as described  <cit> .
dna-dna hybridization was performed as previously described  <cit>  using a model cary  <dig> bio uv/vis-spectrophotometer equipped with a peltier-thermostatted 6x <dig> multicell changer and a temperature controller with in situ temperature probe .
dna isolation and complete genome sequencing
genomic dna for sequencing of strain an4859/ <dig> was obtained from 60 ml mid-log bhis cultures using the dneasy blood and tissue kit  as described by the manufacturer.
the genome sequence was obtained using a combination of ion torrent and illumina miseq sequencing techniques performed at uppsala genome centre .
single end reads were generated with the read size 300 bp on the ion torrent platform and paired end reads were generated using the illumina misseq platform producing reads of 300 bp.
genomic sequences of available strains of currently valid brachyspira species and two strains of the proposed species “b.
hampsonii” were obtained from genbank .
strain names and accession numbers are provided in additional file 1: table s <dig>  draft genome scaffolds for b. aalborgi were downloaded from the metahit website .
genome assembly and average nucleotide identity
single end reads generated by the ion torrent platform were assembled using the roche gs assembler.
the resulting contigs were oriented using the b. hyodysenteriae wa <dig> genome as a reference and joining was performed manually.
reads obtained from illumina miseq paired end sequencing were assembled using mira version  <dig>  <cit> .
before performing the assembly, the reads were trimmed and filtered based on quality, using prinseq  <cit>  and further reduction in the coverage was performed using a custom perl script.
scaffolds were obtained using sspace  <cit> .
the scaffolds were aligned to ion torrent assembly using mauve  <cit>  to find possible overlaps.
overlapping scaffolds were manually joined in consed  <cit> .
a putative plasmid sequence was predicted by blasting the b. hyodysenteriae wa <dig> plasmid’s  sequence  <cit>  against the b. suanatina contigs.
the average nucleotide identity  between b. suanatina and type strains of all the seven valid brachyspira species was calculated using jspecies v <dig> .
<dig>  <cit>  using blast  and mummer  algorithms with default parameters.
genome annotation and core genome analysis
b. suanatina scaffolds were annotated using the gendb platform  <cit> .
predictions made by the gendb pipeline were also manually reviewed and curated.
prophages were predicted using phast  <cit>  and by manual curation in the annotated genome.
clusters of orthologous group  categories were assigned to genes by performing blast search against the cog database  using the threshold of > =30 % amino acid identity and e-value < =  <dig> .
in order to obtain annotations based on the latest versions of search databases, all the genomes used in the analysis were also annotated with gendb.
comparative analysis of the b. suanatina genome with other brachyspira genomes including pan and core genome calculations was performed using the edgar software framework for prokaryotic genomes  <cit> .
whole genome alignment of b. suanatina with the complete genomes of b. hyodysenteriae wa <dig> and b. intermedia pws/at was performed using the artemis comparison tool .
phylogenetic analysis based on core genomes of all brachyspira species including b. suanatina was also performed as a part of the edgar comparative analysis.
for phylogenetic analysis, multiple sequence alignments of amino acid sequences of the core genome were obtained in muscle  <cit> , unaligned blocks were removed using gblocks  <cit>  and a phylogenetic tree of concatenated alignment was generated in phylip using the neighbour joining method.
bootstrap analysis was performed on the tree using  <dig>  replicates.
phylogenetic analysis of housekeeping genes
for phylogenetic analysis of housekeeping genes,  <dig> protein coding genes that are universally conserved among bacteria  <cit>  were selected .
nucleotide sequences of these genes for b. hyodysenteriae b78t, b. hyodysenteriae wa <dig>  b. intermedia pws/at, b. murdochii 56-150t, b. pilosicoli 95/ <dig>  b. pilosicoli b <dig>  b. pilosicoli wesb, b. pilosicoli p43/6/78t, treponema denticola at and borrelia burgdorferi b31t were obtained from the ncbi gene database .
nucleotide sequences for b. aalborgi 513at, b. alvinipulli c1t, b. innocens b256t, “b.
hampsonii”  <dig>  “b.
hampsonii”  <dig> and b. suanatina an4859/ <dig> were predicted by performing homology searches against their gene sequences with b. hyodysenteriae strain wa1’s genes as reference.
sequences were aligned using muscle  <cit>  and conserved blocks were obtained using gblocks  <cit> .
genes with < 50 % alignment conservation were removed leaving  <dig> genes.
phylogenetic inference using the concatenated alignment of those  <dig> genes was drawn using the maximum likelihood method   <cit>  with  <dig> bootstrap replicates.
the general time reversible   <cit>  model with discrete gamma distribution allowing some proportion of invariable sites  as predicted by jmodeltest  <cit>  was used to model the evolutionary rates among sites.
phylogenetic analysis was performed in mega <dig>  an evolutionary analysis tool  <cit> .
RESULTS
morphology and growth characteristics
cells of strain an4859/ <dig> were ultrastructurally similar in appearance but slightly shorter and thinner and possessed fewer flagella than those of b. hyodysenteriae b78t and b. intermedia pws/at .
the cells were shaped as spirochaetes with tapered ends and periplasmic flagella that were inserted on a line subterminally at both cell ends.table  <dig> morphologic comparison of indole-positive species of genus brachyspira
a, b
b. hyodysenteriae
b. intermedia
b. suanatina
aresults for bb.
hyodysenteriae and bb.
intermedia from previous studies  <cit>
bresults from present study
strains an4859/ <dig> and an3949:2/ <dig> that were cultured anaerobically in bhis broth at 37 °c entered the log-phase at approximately 31 h of incubation and reached maximum population densities at 36–38 h post inoculation.
maximum cfu/ml was 2 ×  <dig> cells/ml at 36 h. at mid log-phase, the population doubling time was approximately 2 h. both strains grew equally well at  <dig> and 42 °c, but growth was inhibited at  <dig> and 46 °c.
biochemical tests
the enzymatic profiles are shown in table  <dig>  all six b. suanatina isolates had a similar enzymatic profile to b. hyodysenteriae strain b78t, but reactions were generally weaker.table  <dig> enzymatic profile of b. hyodysenteriae type strain  and six strains of b. suanatina
brachyspira panel
a
b78t
a
brachyspira panel according to fellström et al.,  <dig>  results according to råsbäck et al.,  <dig>  haem = intensity of β-haemolysis; is = indole spot test; hipp = hippurate cleavage capacity; α-galactosidase; β-glucosidase; s = strong reaction; + = positive reaction, − = negative reaction
bapi-zym test:  <dig> control;  <dig> alkaline phosphatase;  <dig> c <dig> esterase;  <dig> c <dig> esterase lipase;  <dig> c <dig> lipase;  <dig> leucine arylamidase;  <dig> valine arylamidase;  <dig> cystine arylamidase;  <dig> trypsin;  <dig> chymotrypsin;  <dig> acid phosphatase;  <dig> naphtolphosphohydrolase;  <dig> α-galactosidase;  <dig> β-galactosidase;  <dig> β-glucoronidase;  <dig> α-glucosidase;  <dig> β-glucosidase; 18 n-acetyl-β-glucosaminidase;  <dig> α-mannosidase;  <dig> α-fucosidase; values 0– <dig> assigned according to the manufacturer’s instruction 0 = negative reaction, 5 = maximum intensity
growth substrate utilization
the maximum growth  of strains b. hyodysenteriae b78t and b. suanatina strains an4859/ <dig> and an3949:2/ <dig> in his but without added growth substrate was  <dig> – <dig> .
the following compounds were utilized by all three strains as growth substrates: d-glucose, d-fructose, d-galactose, d-mannose, sucrose, d-trehalose, lactose, n-acetyl-d-glucosamine, pyruvate .
three out of six tested b. suanatina-strains  metabolized d-mannitol, but b. hyodysenteriae b78t did not.
substrates that were not metabolized by the three tested strains included: d-cellobiose, d-fucose, d-melobiose, d-ribose, d-raffinose, d-xylose, l-arabinose, l-fucose, l-rhamnose, n-acetyl-d-galactosamine, soluble starch, galactosamine, na-acetate, lactate, esculin, succinate, pectin, 50x mem amino acids, cellulose, glutathione, glycerol, glycogen, poly-galacturonic acid, hyaluronic acid and hog gastric mucin.
the control strain b78t showed results consistent with an earlier study  <cit> .
dna-dna hybridization
the percent dna-dna similarities of b. suanatina strain an4859/ <dig> to other brachyspira species type strains are given in table  <dig> table  <dig> percent dna-dna similarity of b. suanatina an4859/ <dig> to type strains  and percent anib/anim values
brachyspira sp.
b. hyodysenteriae
b. intermedia
b. pilosicoli
b. innocens
b. murdochii
b. alvinipulli
b. aalborgi
ain  <dig> x ssc at 59 °c
assembly statistics and average nucleotide identity
the final assembly of b. suanatina an4859/ <dig> consisted of  <dig> scaffolds comprising  <dig> ,009 bp with a gc content of 27 %.
three large scaffolds of sizes  <dig> ,936 bp,  <dig>  bp and  <dig>  bp covered 98 % of the assembled genome.
the remaining scaffolds ranged in size from 450 bp to  <dig>  bp.
one of the scaffolds contained a putative plasmid sequence of  <dig>  bp sharing 88 % identity over 51 % of its length with the b. hyodysenteriae strain wa <dig> plasmid  sequence, see below.
the assembled genome of b. suanatina strain an4859/ <dig> was larger than that of the available genomes of b. hyodysenteriae strains b78t and wa <dig>  but smaller than the genome of b. intermedia strain pws/at.
the ani of b. suanatina to type strains of other brachyspira species was lower than 95 %, with the highest value obtained between b. suanatina and b. intermedia .
general genomic features and putative plasmid
general genomic features of b. suanatina strain an4859/ <dig>  b. hyodysenteriae  and b. intermedia  are depicted in table  <dig>  the draft genome assembly of b. suanatina an4859/ <dig> consisted of  <dig>  open reading frames , including  <dig> orfs  belonging to a putative plasmid sequence.
by comparing the plasmid genes of b. suanatina with that of pbhwa <dig>  we found  <dig> cdss that were conserved in both strains.
these include four genes encoding for epimerases, six genes encoding for transferases, and one gene each encoding for radical sterile alpha motif  domain containing protein, plasmid partition protein, putative replicative dna helicase, dna primase-like protein, integrase, lipopolysaccharide biosynthesis protein, dehydratase, reductase and a hypothetical protein.
genes that were present in b. suanatina but were not found in the plasmid sequence of b. hyodysenteriae were one glycosyl transferase, one peptide abc transporter substrate-binding protein, one radical sam domain-containing protein, five hypothetical proteins, one conserved hypothetical protein and two of unknown function.table  <dig> genome assembly statisticsa, b: summary of general genomic features of b. suanatina strain an4859/ <dig> and its comparison with genomes of b. hyodysenteriae strains b78t and wa <dig> and b. intermedia strain pws/at
b. suanatina an4859/03
b. hyodysenteriae b78
t
b. hyodysenteriae wa1
b. intermedia pws/a
t
adata from references  <cit>  and from the present study
bbb.
hyodysenteriae b78t data taken from ncbi wgs assembly gcf_ <dig> with accession nz_arsy00000000
prophages and mobile genetic elements
the b. suanatina genome contained a 15 kb region similar to the core region of the prophage-like gene transfer agent  named vsh- <dig>  of b. hyodysenteriae .
by comparing this region with the gta region of b. hyodysenteriae strain wa <dig>   <cit> , we found that all the late function genes were present in the same order in both strains.
of the twelve orfs of unknown function described in b. hyodysenteriae wa1’s gta region, orfa, orfb, orfe, orff, orfg, orf <dig>  orf <dig> and orf <dig> were conserved in b. suanatina, and orfc, orf <dig>  orf <dig> and orf <dig> were not present in the b. suanatina genome.
instead six other orfs of unknown function were found in the corresponding region in b. suanatina.
two of them were also conserved in the vsh- <dig> regions of b. hyodysenteriae b78t  and b. intermedia pws/at .
one of them was homologous to a gene encoding for sulfatase  outside the vsh- <dig> region in b. intermedia pws/at and three of them were unique to the b. suanatina genome.
besides these six orfs, three additional orfs were present between hvp <dig> and hvp <dig>  two of them were also present in the vsh- <dig> region of b. intermedia encoding for hypothetical proteins and one of them was unique to b. suanatina.
another additional orf was present between hvp <dig> and lys gene in the b. suanatina genome and was homologous to hvp <dig> gene in vsh- <dig> region of b. intermedia pws/at.
a putative bacteriophage region, bsp was predicted in the b. suanatina genome, which has not been described in genomes of b. hyodysenteriae strain wa <dig> or b. intermedia strain pws/at.
analysis revealed that bsp was present approximately 33 kb upstream to the vsh- <dig> region.
the bsp region spanned over a  <dig>  kb sequence length and consisted of  <dig> cdss  including one vsh- <dig> related gene.
the bsp region was not present in the complete genome of b. hyodysenteriae strain wa <dig> genome, but a homologous bacteriophage region was found in scaffold number  <dig> of the b. hyodysenteriae b78t assembly  and consisted of  <dig> cdss .
a bacteriophage region almost identical to bsp was found in the b. pilosicoli p43/6/78t genome  spanning over 39 kb sequence and it contained  <dig> cdss  including one vsh- <dig> related gene.
the organization of genes in bsp and its comparison with identical regions in b. hyodysenteriae b78t and b. pilosicoli p43/6/78t is shown in fig.
<dig>  further analysis revealed that a significant portion of the bsp region was conserved and had been previously reported as pp <dig>  pp <dig> and pp <dig>  in b. pilosicoli 95/ <dig>  b. pilosicoli strains b <dig> and b. pilosicoli wesb respectively  <cit> , and pm <dig>  pm <dig> and pm <dig> in b. murdochii strain 56-150t  <cit> .
some of the bsp genes were also found in genomes of b. innocens b256t and b. alvinipulli c1t scattered in the genome either as single genes or in small clusters of phage related genes.
these include arqi01_ <dig>  arqi01_ <dig>  arqi01_ <dig>  arqi01_ <dig>  arqi01_ <dig>  arqi01_ <dig>  arqi01_347-arqi01_ <dig>  arqi01_462-arqi01_ <dig>  arqi01_469-arqi01_ <dig>  arqi01_ <dig>  arqi01_ <dig> in the b. innocens genome  and jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_  <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig> – jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig>  jadf01_ <dig> in the b. alvinipulli genome .
only two genes of bsp shared homology with b. intermedia strain pws/at genes.
these were: bint_ <dig> encoding hypothetical protein of b. intermedia phage and bint_ <dig> encoding putative integrase in b. intermedia.
a part from these bacteriophage regions, some phage related genes and mobile genetic elements were also found in the b. suanatina genome including putative integrases, terminases, tail protein, portal protein that were not found in b. hyodysenteriae genome except for the two reported phage related genes bhwa_ <dig> and bhwa_ <dig>  bacteriophage regions of b. intermedia pws/at  <cit>  were not found in the b. suanatina genome except one gene  that encodes for a hypothetical membrane spanning protein, and is also conserved in b. hyodysenteriae strains wa <dig> and b78t.fig.
<dig> bacteriophage region comparison: comparison of the organization of the bacteriophage region bsp in b. suanatina an4859/ <dig> with the homologous bacteriophage regions in b. pilosicoli p43/6/78t and b. hyodysenteriae b78t genomes
gene content and synteny conservation
the pan genome of b. suanatina an4859/ <dig>  b. hyodysenteriae b78t, b. hyodysenteriae wa <dig> and b. intermedia pws/at consisted of  <dig>  genes.
among them  <dig>  genes were present in all four strains,  <dig>  genes were shared between b. hyodysenteriae wa <dig> and b. suanatina,  <dig>  genes were shared between b. hyodysenteriae b78t and b. suanatina, and  <dig>  genes were shared between b. intermedia and b. suanatina .
the remaining genes had their highest similarities with other brachyspira species.
there were  <dig> genes in b. suanatina that were not present in any of the other brachyspira species, most of them consisted of hypothetical proteins and some had matches in clostridium spp.
and bacillus spp.fig.
<dig> global gene conservation and comparison of functional classification of genes: the figure shows the number of genes and their classification in different cog categories from b. suanatina an4859/ <dig>  b. hyodysenteriae b78t, b. hyodysenteriae wa <dig> and b. intermedia pws/at  venn diagram  the bar chart shows the distribution of genes in different functional categories.
genes were assigned to different categories using blastp against the cog database
by comparing the functional distribution of genes in different cog categories, the four strains b. suanatina an4859/ <dig>  b. hyodysenteriae b78t, b. hyodysenteriae wa <dig> and b. intermedia pws/at contained similar numbers of genes for particular functions  .
in contrast to gene content, gene synteny was highly conserved between b. hyodysenteriae and b. suanatina with only two major genomic rearrangement events that were found in b. hyodysenteriae wa <dig> and b. suanatina.
some insertions were also found in the b. suanatina genome, when compared with b. hyodysenteriae wa <dig>  comprising mainly of genes encoding for unknown proteins, bacteriophages and mobile genetic elements.
in b. intermedia and b. suanatina, there was a low level of conservation of synteny as compared to the gene content.
the alignment of genomic sequences for two isolates of b. suanatina and b. intermedia pws/at using act identified extensive genomic rearrangements, including at least eight major inversions .
phylogenetic inference based on the core genome and housekeeping genes
in order to infer the evolutionary relationship of b. suanatina with other brachyspira species, phylogenetic analysis was performed on the amino acid sequences of core genomes and nucleotide sequences of housekeeping genes.
in total  <dig>  genes were found conserved in all brachyspira species and their phylogenetic tree exhibited a close relationship between b. suanatina and b. intermedia whereas b. hyodysenteriae wa <dig> and b. hyodysenteriae b78t formed a distinct clade.
the phylogenetic tree is shown in fig.
<dig>  phylogenetic analysis of  <dig> housekeeping genes also showed a highly similar tree topology to the core genome analysis and confirmed a clear distinction between b. hyodysenteriae and b. suanatina.
the confidence level for the tree topology using  <dig> bootstrap replicates ranged from 60 % to 100 % .fig.
<dig> phylogenetic tree of core genome: radial unrooted phylogenetic representation based on concatenated amino acid sequences of  <dig>  core genes of available genomes of genus brachyspira and b. suanatina strain an4859/ <dig>  the clade containing b. hyodysenteriae, b. intermedia and b. suanatina strains are shown in bold.
the tree was constructed using the neighbour joining method.
numbers above the branches represent support values  obtained from  <dig> bootstrap replicates.
the tree is drawn to scale with branch lengths measured in the number of substitutions per site
fig.
<dig> phylogenetic tree of housekeeping genes: the evolutionary history of  <dig> housekeeping genes of a selection of strains representing phylum spirochaetes was inferred by the maximum likelihood method  <cit>  based on the general time reversible model  <cit> .
a discrete gamma distribution was used to model evolutionary rate differences among sites .
the model allowed some sites to be in variable .
the clade containing b. hyodysenteriae, b. intermedia and b. suanatina strains are shown in bold.
treponema denticola at and borrelia burgdorferi b31t were used as outgroup.
numbers above branches represent support values obtained from  <dig> bootstrap replicates.
the tree is drawn to scale with branch lengths measured in the number of substitutions per site.
evolutionary analysis was conducted in mega <dig>  <cit>
discussion
using the threshold value of 70 % dna-dna similarity for the delineation of bacterial species, as recommended by the ad hoc committee  <cit> , b. suanatina strain an4859/ <dig> does not belong to any of the currently recognized seven brachyspira spp.
the highest observed % dna-dna similarities  were to the type strains of b. hyodysenteriae and b. intermedia.
our finding was further confirmed by the ani values, which are considered to correspond with the % ddh values  <cit> .
since the degree of robustness of anim is higher compared to anib for genomes having >60 % ddh values  <cit>  we calculated the ani values using both algorithms and the values obtained from them are in agreement to each other.
furthermore, the evolutionary history inferred from analysis of the 16s rrna gene  <cit> , housekeeping genes and the core genome  suggests a common ancestor of the three related strains b78t, an4859/ <dig>  and pws/at, which each represent indole test-positive phenotypes of b. hyodysenteriae, b. suanatina and b. intermedia, respectively.
high bootstrap values supported the phylogram based on housekeeping genes and the assertion that b. suanatina and b. intermedia make a separate clade that is distinct from the clade formed by b. hyodysenteriae strains.
interestingly, the phenotypes of b. hyodysenteriae, b. suanatina and b. intermedia were very similar, with the strong β-haemolysis of the two former being strongly positive vs. the weak β-haemolysis in b. intermedia.
notably, in this study, no consistent diagnostic phenotypic feature was found to differentiate b. suanatina from b. hyodysenteriae.
the strains of b. suanatina in this study showed varying capabilities of using d-mannitol as growth substrate, whereas strain b78t could not.
even though this phenotype could be further investigated it would not qualify as a reliable diagnostic feature.
interestingly, among the three strains, the type strain of b. intermedia pws/at seems to have undergone the highest number of changes since its divergence from the common ancestor based on the housekeeping genes phylogram.
when comparing the conservation of gene synteny of b. suanatina with b. hyodysenteriae wa <dig> and b. intermedia pws/at, very low conservation was found between b. suanatina and b. intermedia as compared to b. suanatina and b. hyodysenteriae wa <dig>  these results also point to the high level of genomic rearrangements that have taken place in the genome of b. intermedia after these species diverged.
previous reports based on multi locus sequence typing have suggested that b. intermedia is a genetically diverse species that may possibly consist of several species  <cit> .
the same applies when strongly β-haemolytic indole-test positive isolates, to which b. hyodysenteriae and b. suanatina and other isolates of unknown species affiliation, are investigated  <cit> .
clearly, this group of brachyspira isolates needs more detailed genomic and phenotypic study to elucidate evolutionary and taxonomic relationships.
general genomic features including the distribution of genes in similar functional classes in b. suanatina, b. hyodysenteriae and b. intermedia support the fact that they are highly similar species.
however the genome size of b. hyodysenteriae is smaller as compared to b. suanatina and b. intermedia.
speculatively, the reason for the smaller genome size could be that b. hyodysenteriae has adapted to the rich environment of porcine intestinal tissues through reductive evolution.
reductive evolution is a process of reduction in genome size of a host associated bacterium by the loss of genes rendered non-essential  <cit> .
the presence of intact bacteriophage region bsp in b. suanatina and its conservation in other brachyspira genomes including b. hyodysenteriae b78t while absence in b. hyodysenteriae wa <dig> and b. intermedia pws/at also suggests the loss of this prophage in two later species due to the reductive evolution.
events of horizontal gene transfer  are also evident in the genome of b. suanatina when examining the genes that had no homology in the genomes of brachyspira species with other bacteria.
these include the presence of some genes that had their highest similarities in clostridium spp.
and bacillus spp.
both these species are also inhabitants of the large intestine where the possible exchange of genetic material with b. suanatina could have taken place.
events of reductive evolution and hgt in the genomes of brachyspira species have been described earlier  <cit> .
the unusual phage-like gene transfer agent vsh- <dig> produced by b. hyodysenteriae is likely involved in natural gene transfer and recombination within different strains of the species  <cit> .
it spans over a  <dig>  kb region of the genome of b. hyodysenteriae b204r and contains all the late function genes  <cit> .
the region containing genes homologous to vsh- <dig> has also been identified in the genomes of different strains of b. hyodysenteriae, b. intermedia, b. pilosicoli and b. murdochii  <cit> .
a region similar to vsh- <dig> was also found in the genome of b. suanatina in this study.
although it is not clear if this gta in b. suanatina strain an4859/ <dig> is functional, the conservation of the late function genes and the presence of unknown genes in the gta region suggest functional capabilities.
CONCLUSIONS
this study demonstrates the importance of integrating genomic data into the classification of different bacterial species.
based on genomic characteristics we conclude that b. suanatina does not belong to any of the currently validated seven brachyspira species.
the use of genomic data from more than one strain of b. suanatina species of the strain will facilitate future comparative genomic studies.
description of brachyspira suanatina sp.
nov.
brachyspira suanatina su.a.na.ti'na.
l. masc n. suus swine/pig, l. fem.
adj.
anatina from/of a duck, n.l.
fem.
adj.
suanatina of pig and duck, referring to the host animals  from which this spirochete was isolated.
cells are oxygen-tolerant but anaerobic and gram-negative of a typical spirochaetal ultrastructure.
cell dimensions are  <dig>  ±  <dig>  μm in length and  <dig>  ±  <dig>  μm in diameter.
the number of periplasmic flagella are 14– <dig>  per cell inserted subterminally at cells ends.
bacterial growth occurs at  <dig> and 42 °c but not at  <dig> and 46 °c.
when grown anaerobically in bhis at 37 °c cultures reach log-phase at approximately 31 h of inoculation and reach maximum population densities of approximately 2 × 108 cfu/ml at 36–38 h post inoculation.
at mid log-phase, the population doubling time is approximately 2 h. cells grow anaerobically on fastidious anaerobe agar with swarming colonies and a strong β-haemolysis on trypticase soy agar.
gives positive reactions in tests for tryptophanase activity, alkaline phosphatase, c <dig> esterase, c <dig> esterase lipase, acid phosphatase, naphtolphosphohydrolase, β-galactosidase.
α-glucosidase, β-glucosidase.
the following substrates can be utilized as single carbon, nitrogen and energy sources for growth: d-glucose, d-fructose, d-galactose, d-mannose, sucrose, d-trehalose, lactose, n-acetyl-d-glucosamine and pyruvate.
d-mannitol sustain growth of some isolates.
growth was not detected with d-cellobiose, d-fucose, d-melobiose, d-ribose, d-raffinose, d-xylose, l-arabinose, l-fucose, l-rhamnose, n-acetyl-d-galactosamine, soluble starch, galactosamine, na-acetate, lactate, esculin, succinate, pectin, 50x mem amino acids, cellulose, glutathione, glycerol, glycogen, poly-galacturonic acid, hyaluronic acid and hog gastric mucin.
the draft genomic dna g + c content of the type strain is 27 % and the genome comprises  <dig> ,009 bp, and  <dig>  orfs.
dna-dna hybridization of the type strain against all type strains of recognized brachyspira spp.
and anib/anim values  support the taxonomic classification of b. suanatina as a sp.
nov. the evolutionary history of 16s rrna and housekeeping genes, and the core genome shows a common ancestor of b. hyodysenteriae, b. intermedia and b. suanatina.
a bacteriophage region and a phage-like gta similar to vsh- <dig> in b. hyodysenteriae and a putative plasmid sequence related to b. hyodysenteriae strain wa <dig> plasmid  are present in the genome of b. suanatina strain an4859/ <dig>  type strain is enteropathogenic in pigs with swine dysentery-like symptoms when using the type strain, and mild diarrhoea was associated with challenge with an isolate of mallard origin  <cit> .
the type strain, an4859/03t  was obtained from a swedish pig with swine dysentery-like disease.
accession numbers
draft genome sequences of b. suanatina an4859/ <dig> have been deposited in the european nucleotide archive  database with scaffold accession numbers: cvlb01000001-cvlb <dig> under the study accession number: prjeb <dig>
additional files
additional file 1: table s <dig>  strain names and accession numbers of available genome used for the study.
additional file 2: table s <dig>  functional distribution of genes in different cog categories of strains b. suanatina an4859/ <dig>  b. hyodysenteriae b78t, b. hyodysenteriae wa <dig> and b. intermedia pws/at.
additional file 3: figure s <dig>  pairwise genome alignment of b. hyodysenteriae wa <dig>  b. intermedia pws/at and b. suanatina an4859/ <dig>  the genomes were aligned using artemis comparison tool .
alignment was performed according to the sequence start of b. hyodysenteriae wa <dig> and visualised in act with a cut-off set to blast scores >  <dig>  red bars indicate the similar regions in the same orientation, whereas blue bars show inverted regions.
abbreviations
actartemis comparison tool
aniaverage nucleotide identity
bhibrain heart infusion culture broth
bhisbrain heart infusion culture broth with calf serum
cdscoding sequence
cfucolony forming unit
cogclusters of orthologous group
ddhdna-dna hybridization
faafastidious anaerobe agar
fcsfoetal calf serum
gtagene transfer agent
hiheart infusion culture broth
hisheart infusion culture broth with foetal calf serum
hgthorizontal gene transfer
mlmaximum likelihood
noxnadh oxidase gene
orfopen reading frame
competing interests
the authors declare that they have no competing interests.
authors’ contributions
mm, ebr and dsj designed the experiments.
mm, sz, dsj and tr conducted experiments and data analyses.
mm and dsj drafted the manuscript.
all authors contributed to the final manuscript.
all authors read and approved the final manuscript.

 the advent of next generation sequencing technologies with their low cost, high throughput and rapid speed has strongly influenced modern taxonomy. many taxonomists therefore turn towards integrating genome sequence analysis for formal proposals of new bacterial taxa or reclassification of old taxa. early classification methods were based solely on phenotypic characteristics, but the introduction of dna-dna hybridization and 16s ribosomal rna gene phylogeny in taxonomy changed the species definition. however, due to the high level of conservation of the 16s rrna gene and the labour intensive and error prone nature of ddh experiments, there has been a constant demand for advanced methods at complete genome level. average nucleotide identity, multi locus sequence analysis of conserved genes and phylogenetic analysis of core genomes are some of the suggested methods. brachyspira species are ggram-negative, anaerobic but oxygen-tolerant spirochaetes that colonize the intestines of many mammals and birds. it forms a clade within phylum spirochaetes that has undergone repeated taxonomic changes over the past few decades. initially the isolates colonizing pigs were allocated to genus treponema, however they were later transferred to genus serpula and then serpulina. in the genera serpulina and brachyspira were unified leading to the reclassification of serpulina species. genus brachyspira now encompasses seven valid species: b. aalborgi, b. alvinipulli, b. hyodysenteriae, b. intermedia, b. innocens, b. murdochii and b. pilosicoli with b. aalborgi as a type species. these species vary in terms of their host range and enteropathogenicity. brachyspira hyodysenteriae and b. pilosicoli are well known porcine pathogens. furthermore, b. pilosicoli can cause gastrointestinal disease and egg-production losses in chickens, and it can also colonize the human intestine leading to human intestinal spirochaetosis. however, most cases in humans are associated with b. aalborgi, at least in industrialized countries. moreover, b. alvinipulli and b. intermedia cause diarrhoea in chickens. however, b. intermedia, is in general considered as a harmless commensal in pigs similar to b. innocens and b. murdochii. there are several proposed species without current standing in nomenclature, including “b. suanatina”. some members of the genus brachyspira have low interspecies 16s rrna gene variability that does not always support species delineations, especially among the indole test-positive species b. hyodysenteriae, b, intermedia and b. suanatina that require additional methods for species confirmation. the proposed new species b. suanatina was first isolated from pigs with dysentery-like disease in sweden and denmark, and the enteropathogenicity of strain an4859/ was confirmed by experimental challenge in pigs. the phenotypic characterization based on the intensity of β-haemolysis and enzyme activities showed that these isolates were indistinguishable from b. hyodysenteriae, but detection by a species-specific pcr was unsuccessful. based on 16s rrna and nadh oxidase gene similarities, these brachyspira isolates from pigs were similar to some isolates from mallards and were distinct from all currently recognized brachyspira species, suggesting that they were members of a putatively unrecognized species. available data suggest that b. suanatina most likely is a rare genotype in both pigs and mallards, and it has so far only been isolated in northern europe. identification of b. suanatina currently relies on nox gene sequencing and/or nox-rflp. the objective of the present study was to identify the taxonomical position of b. suanatina based on phenotypic characteristics and complete genome analyses. we have used a “taxono-genomics” approach that is based on the polyphasic combination of phenotypic and genomic characteristics. methods bacterial cultures and growth conditions type and reference strains b. hyodysenteriae b78t, b204r, b. intermedia pws/at, b. pilosicoli p43/6/78t, b. innocens b256t, b. murdochii 56-150t, b. alvinipulli c1t and b. aalborgi 513at were originally obtained from the american type culture collection and stored in liquid nitrogen in the brachyspira strain collection of the national veterinary institute. brachyspira suanatina strains an4859/ an1681:1/ an2384/ and dk12570- from pigs, and an3949:2/ and an1418:2/ from mallards were isolated in sweden and denmark as previously described and stored in liquid nitrogen. the bacteria were retrieved from storage and were grown anaerobically on fastidious anaerobe agar plates with 10 % equine blood at 42 °c for 72 h and subcultured once for 48 h before use. cells for electron microscopy were harvested from broth cultures grown anaerobically on a shaker at 37 °c to mid log phase in brain heart infusion culture broth with 10 % foetal calf serum. faa plate cultures and broth culture were used to test the capability of all b. suanatina strains to grow at and 46 °c. for the generation of growth curves, broth cultures of b. suanatina strains an4859/ and an3949:2/ and b. hyodysenteriae b78t were prepared in triplicates in 10 ml bhis at cells/ml and grown at 37 °c. population densities were determined every 4 h between and 70 h of incubation by colony forming unit counts obtained from serial dilutions on faa plates in duplicates and by spectrophotometry. comparisons were made against an in-house standard curve for b. hyodysenteriae b78t. purity was confirmed throughout the studies by phase contrast microscopy. biochemical phenotype and substrate utilization the api-zym test was applied to all b. suanatina strains with b. hyodysenteriae b78t and pseudomonas aeruginosa as control strains. bacterial cells were harvested from faa plates incubated at 42 °c for 3 days, suspended in 2 ml sterile ddh2o to a density of 5– mcfarland and the api-zym test was performed as described by the manufacturer. the test strips were incubated for 4 h in aerobic atmosphere. brachyspira suanatina strains an4859/ an3949:2/ and b. hyodysenteriae strain b78t were used for growth substrate utilization studies. all six b. suanatina strains were tested for growth in d-mannitol because of diverging results obtained with strains an4859/ and an3949:2/ potential substrates were added at wt/v in heart infusion culture broth and autoclaved or sterile filtered. the fcs concentration in the hi was selected at 7 % from test runs at and 10 % fcs as described previously. triplicate flasks with 10 ml hi and 7 % fcs culture broth were used for each isolate and substrate, and were inoculated with cells from a stock solution made from 48 h old faa cultures. cultures were incubated anaerobically at 37 °c on a shaking platform. glucose was used as positive control and his without added substrate as negative control. growth was determined spectrophotometrically at and 9 days post inoculation for all cultures, and the purity, cell density and the condition of the cells were examined by phase contrast microscopy. compounds producing ≤ od were not considered to be growth substrates. these cultures contained degenerated and/or non-motile cells as determined by phase contrast microscopy. cultures producing ≥ od values and moderate to profuse growth of motile spirochaetes on at least one occasion were considered to be growth substrates. the following compounds were investigated: d-glucose, d-fructose, d-fucose, d-ribose, d-mannose, sucrose, lactate, d-xylose, d-trehalose, d-galactose, d-melobiose, d-cellobiose, l-arabinose, l-fucose, l-rhamnose, d-raffinose, d-galactosamine, n-acetyl-d-glucosamine, n-acetyl-d-galactosamine, mannitol, pyruvate, esculin, porcine gastric mucin, 50x mem amino acids, pectin, hyaluronic acid, glutathione, poly-galacturonic acid, cellulose, glycerol, succinate and glycogen. electron microscopy strain an4859/ was grown in bhis to mid-log phase, centrifuged at  × g for 3 min and washed twice in 1 ml  m sodium phosphate buffer at ph . the pellet was resuspended in  ml each of sterile ddh2o and negatively stained with  ml 2 % phosphotungistic acid at ph . the cells were dropped on grids and examined. for scanning electron microscopy the cells were harvested and washed as described above, dropped on grids and coated with gold-palladium with standard procedures. the samples were examined with a jsm scanning electron microscope. ten cells or more were examined to determine cell dimensions, the location of attachment sites of periplasmatic flagella and the number of flagella per cell. dna-dna hybridization the ddh of b. suanatina strain an4859/ was performed with all the seven recognized brachyspira species. cells were grown to mid-log phase in bhis, checked for purity and centrifuged at 4 °c at  g followed by washing in pbs twice before mixing with 1:1 v/v iso-propanol. dna isolation and ddh were performed by dsmz gmbh deutsche sammlung von mikroorganismen und zellkulturen gmbh dna. dna was extracted using a french pressure cell and was purified by chromatography on hydroxyapatite as described. dna-dna hybridization was performed as previously described using a model cary bio uv/vis-spectrophotometer equipped with a peltier-thermostatted 6x multicell changer and a temperature controller with in situ temperature probe. dna isolation and complete genome sequencing genomic dna for sequencing of strain an4859/ was obtained from 60 ml mid-log bhis cultures using the dneasy blood and tissue kit as described by the manufacturer. the genome sequence was obtained using a combination of ion torrent and illumina miseq sequencing techniques performed at uppsala genome centre. single end reads were generated with the read size 300 bp on the ion torrent platform and paired end reads were generated using the illumina misseq platform producing reads of 300 bp. genomic sequences of available strains of currently valid brachyspira species and two strains of the proposed species “b. hampsonii” were obtained from genbank. strain names and accession numbers are provided in additional file 1: table s draft genome scaffolds for b. aalborgi were downloaded from the metahit website. genome assembly and average nucleotide identity single end reads generated by the ion torrent platform were assembled using the roche gs assembler. the resulting contigs were oriented using the b. hyodysenteriae wa genome as a reference and joining was performed manually. reads obtained from illumina miseq paired end sequencing were assembled using mira version. before performing the assembly, the reads were trimmed and filtered based on quality, using prinseq and further reduction in the coverage was performed using a custom perl script. scaffolds were obtained using sspace. the scaffolds were aligned to ion torrent assembly using mauve to find possible overlaps. overlapping scaffolds were manually joined in consed. a putative plasmid sequence was predicted by blasting the b. hyodysenteriae wa plasmid’s sequence against the b. suanatina contigs. the average nucleotide identity between b. suanatina and type strains of all the seven valid brachyspira species was calculated using jspecies v. using blast and mummer algorithms with default parameters. genome annotation and core genome analysis b. suanatina scaffolds were annotated using the gendb platform. predictions made by the gendb pipeline were also manually reviewed and curated. prophages were predicted using phast and by manual curation in the annotated genome. clusters of orthologous group categories were assigned to genes by performing blast search against the cog database using the threshold of > =30 % amino acid identity and e-value < =. in order to obtain annotations based on the latest versions of search databases, all the genomes used in the analysis were also annotated with gendb. comparative analysis of the b. suanatina genome with other brachyspira genomes including pan and core genome calculations was performed using the edgar software framework for prokaryotic genomes. whole genome alignment of b. suanatina with the complete genomes of b. hyodysenteriae wa and b. intermedia pws/at was performed using the artemis comparison tool. phylogenetic analysis based on core genomes of all brachyspira species including b. suanatina was also performed as a part of the edgar comparative analysis. for phylogenetic analysis, multiple sequence alignments of amino acid sequences of the core genome were obtained in muscle, unaligned blocks were removed using gblocks and a phylogenetic tree of concatenated alignment was generated in phylip using the neighbour joining method. bootstrap analysis was performed on the tree using replicates. phylogenetic analysis of housekeeping genes for phylogenetic analysis of housekeeping genes, protein coding genes that are universally conserved among bacteria were selected. nucleotide sequences of these genes for b. hyodysenteriae b78t, b. hyodysenteriae wa b. intermedia pws/at, b. murdochii 56-150t, b. pilosicoli 95/ b. pilosicoli b b. pilosicoli wesb, b. pilosicoli p43/6/78t, treponema denticola at and borrelia burgdorferi b31t were obtained from the ncbi gene database. nucleotide sequences for b. aalborgi 513at, b. alvinipulli c1t, b. innocens b256t, “b. hampsonii” “b. hampsonii” and b. suanatina an4859/ were predicted by performing homology searches against their gene sequences with b. hyodysenteriae strain wa1’s genes as reference. sequences were aligned using muscle and conserved blocks were obtained using gblocks. genes with < 50 % alignment conservation were removed leaving genes. phylogenetic inference using the concatenated alignment of those genes was drawn using the maximum likelihood method with bootstrap replicates. the general time reversible model with discrete gamma distribution allowing some proportion of invariable sites as predicted by jmodeltest was used to model the evolutionary rates among sites. phylogenetic analysis was performed in mega an evolutionary analysis tool.  morphology and growth characteristics cells of strain an4859/ were ultrastructurally similar in appearance but slightly shorter and thinner and possessed fewer flagella than those of b. hyodysenteriae b78t and b. intermedia pws/at. the cells were shaped as spirochaetes with tapered ends and periplasmic flagella that were inserted on a line subterminally at both cell ends.table morphologic comparison of indole-positive species of genus brachyspira a, b b. hyodysenteriae b. intermedia b. suanatina aresults for bb. hyodysenteriae and bb. intermedia from previous studies bresults from present study strains an4859/ and an3949:2/ that were cultured anaerobically in bhis broth at 37 °c entered the log-phase at approximately 31 h of incubation and reached maximum population densities at 36–38 h post inoculation. maximum cfu/ml was 2 ×  cells/ml at 36 h. at mid log-phase, the population doubling time was approximately 2 h. both strains grew equally well at and 42 °c, but growth was inhibited at and 46 °c. biochemical tests the enzymatic profiles are shown in table  all six b. suanatina isolates had a similar enzymatic profile to b. hyodysenteriae strain b78t, but reactions were generally weaker.table enzymatic profile of b. hyodysenteriae type strain and six strains of b. suanatina brachyspira panel a b78t a brachyspira panel according to fellström et al., results according to råsbäck et al., haem = intensity of β-haemolysis; is = indole spot test; hipp = hippurate cleavage capacity; α-galactosidase; β-glucosidase; s = strong reaction; + = positive reaction, − = negative reaction bapi-zym test: control; alkaline phosphatase; c esterase; c esterase lipase; c lipase; leucine arylamidase; valine arylamidase; cystine arylamidase; trypsin; chymotrypsin; acid phosphatase; naphtolphosphohydrolase; α-galactosidase; β-galactosidase; β-glucoronidase; α-glucosidase; β-glucosidase; 18 n-acetyl-β-glucosaminidase; α-mannosidase; α-fucosidase; values 0– assigned according to the manufacturer’s instruction 0 = negative reaction, 5 = maximum intensity growth substrate utilization the maximum growth of strains b. hyodysenteriae b78t and b. suanatina strains an4859/ and an3949:2/ in his but without added growth substrate was –. the following compounds were utilized by all three strains as growth substrates: d-glucose, d-fructose, d-galactose, d-mannose, sucrose, d-trehalose, lactose, n-acetyl-d-glucosamine, pyruvate. three out of six tested b. suanatina-strains metabolized d-mannitol, but b. hyodysenteriae b78t did not. substrates that were not metabolized by the three tested strains included: d-cellobiose, d-fucose, d-melobiose, d-ribose, d-raffinose, d-xylose, l-arabinose, l-fucose, l-rhamnose, n-acetyl-d-galactosamine, soluble starch, galactosamine, na-acetate, lactate, esculin, succinate, pectin, 50x mem amino acids, cellulose, glutathione, glycerol, glycogen, poly-galacturonic acid, hyaluronic acid and hog gastric mucin. the control strain b78t showed results consistent with an earlier study. dna-dna hybridization the percent dna-dna similarities of b. suanatina strain an4859/ to other brachyspira species type strains are given in table  table percent dna-dna similarity of b. suanatina an4859/ to type strains and percent anib/anim values brachyspira sp. b. hyodysenteriae b. intermedia b. pilosicoli b. innocens b. murdochii b. alvinipulli b. aalborgi ain x ssc at 59 °c assembly statistics and average nucleotide identity the final assembly of b. suanatina an4859/ consisted of scaffolds comprising,009 bp with a gc content of 27 %. three large scaffolds of sizes,936 bp,  bp and  bp covered 98 % of the assembled genome. the remaining scaffolds ranged in size from 450 bp to  bp. one of the scaffolds contained a putative plasmid sequence of  bp sharing 88 % identity over 51 % of its length with the b. hyodysenteriae strain wa plasmid sequence, see below. the assembled genome of b. suanatina strain an4859/ was larger than that of the available genomes of b. hyodysenteriae strains b78t and wa but smaller than the genome of b. intermedia strain pws/at. the ani of b. suanatina to type strains of other brachyspira species was lower than 95 %, with the highest value obtained between b. suanatina and b. intermedia. general genomic features and putative plasmid general genomic features of b. suanatina strain an4859/ b. hyodysenteriae and b. intermedia are depicted in table  the draft genome assembly of b. suanatina an4859/ consisted of open reading frames, including orfs belonging to a putative plasmid sequence. by comparing the plasmid genes of b. suanatina with that of pbhwa we found cdss that were conserved in both strains. these include four genes encoding for epimerases, six genes encoding for transferases, and one gene each encoding for radical sterile alpha motif domain containing protein, plasmid partition protein, putative replicative dna helicase, dna primase-like protein, integrase, lipopolysaccharide biosynthesis protein, dehydratase, reductase and a hypothetical protein. genes that were present in b. suanatina but were not found in the plasmid sequence of b. hyodysenteriae were one glycosyl transferase, one peptide abc transporter substrate-binding protein, one radical sam domain-containing protein, five hypothetical proteins, one conserved hypothetical protein and two of unknown function.table genome assembly statisticsa, b: summary of general genomic features of b. suanatina strain an4859/ and its comparison with genomes of b. hyodysenteriae strains b78t and wa and b. intermedia strain pws/at b. suanatina an4859/03 b. hyodysenteriae b78 t b. hyodysenteriae wa1 b. intermedia pws/a t adata from references and from the present study bbb. hyodysenteriae b78t data taken from ncbi wgs assembly gcf_ with accession nz_arsy00000000 prophages and mobile genetic elements the b. suanatina genome contained a 15 kb region similar to the core region of the prophage-like gene transfer agent named vsh- of b. hyodysenteriae. by comparing this region with the gta region of b. hyodysenteriae strain wa, we found that all the late function genes were present in the same order in both strains. of the twelve orfs of unknown function described in b. hyodysenteriae wa1’s gta region, orfa, orfb, orfe, orff, orfg, orf orf and orf were conserved in b. suanatina, and orfc, orf orf and orf were not present in the b. suanatina genome. instead six other orfs of unknown function were found in the corresponding region in b. suanatina. two of them were also conserved in the vsh- regions of b. hyodysenteriae b78t and b. intermedia pws/at. one of them was homologous to a gene encoding for sulfatase outside the vsh- region in b. intermedia pws/at and three of them were unique to the b. suanatina genome. besides these six orfs, three additional orfs were present between hvp and hvp two of them were also present in the vsh- region of b. intermedia encoding for hypothetical proteins and one of them was unique to b. suanatina. another additional orf was present between hvp and lys gene in the b. suanatina genome and was homologous to hvp gene in vsh- region of b. intermedia pws/at. a putative bacteriophage region, bsp was predicted in the b. suanatina genome, which has not been described in genomes of b. hyodysenteriae strain wa or b. intermedia strain pws/at. analysis revealed that bsp was present approximately 33 kb upstream to the vsh- region. the bsp region spanned over a  kb sequence length and consisted of cdss including one vsh- related gene. the bsp region was not present in the complete genome of b. hyodysenteriae strain wa genome, but a homologous bacteriophage region was found in scaffold number of the b. hyodysenteriae b78t assembly and consisted of cdss. a bacteriophage region almost identical to bsp was found in the b. pilosicoli p43/6/78t genome spanning over 39 kb sequence and it contained cdss including one vsh- related gene. the organization of genes in bsp and its comparison with identical regions in b. hyodysenteriae b78t and b. pilosicoli p43/6/78t is shown in fig.  further analysis revealed that a significant portion of the bsp region was conserved and had been previously reported as pp pp and pp in b. pilosicoli 95/ b. pilosicoli strains b and b. pilosicoli wesb respectively, and pm pm and pm in b. murdochii strain 56-150t. some of the bsp genes were also found in genomes of b. innocens b256t and b. alvinipulli c1t scattered in the genome either as single genes or in small clusters of phage related genes. these include arqi01_ arqi01_ arqi01_ arqi01_ arqi01_ arqi01_ arqi01_347-arqi01_ arqi01_462-arqi01_ arqi01_469-arqi01_ arqi01_ arqi01_ in the b. innocens genome and jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ – jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ jadf01_ in the b. alvinipulli genome. only two genes of bsp shared homology with b. intermedia strain pws/at genes. these were: bint_ encoding hypothetical protein of b. intermedia phage and bint_ encoding putative integrase in b. intermedia. a part from these bacteriophage regions, some phage related genes and mobile genetic elements were also found in the b. suanatina genome including putative integrases, terminases, tail protein, portal protein that were not found in b. hyodysenteriae genome except for the two reported phage related genes bhwa_ and bhwa_ bacteriophage regions of b. intermedia pws/at were not found in the b. suanatina genome except one gene that encodes for a hypothetical membrane spanning protein, and is also conserved in b. hyodysenteriae strains wa and b78t.fig. bacteriophage region comparison: comparison of the organization of the bacteriophage region bsp in b. suanatina an4859/ with the homologous bacteriophage regions in b. pilosicoli p43/6/78t and b. hyodysenteriae b78t genomes gene content and synteny conservation the pan genome of b. suanatina an4859/ b. hyodysenteriae b78t, b. hyodysenteriae wa and b. intermedia pws/at consisted of genes. among them genes were present in all four strains, genes were shared between b. hyodysenteriae wa and b. suanatina, genes were shared between b. hyodysenteriae b78t and b. suanatina, and genes were shared between b. intermedia and b. suanatina. the remaining genes had their highest similarities with other brachyspira species. there were genes in b. suanatina that were not present in any of the other brachyspira species, most of them consisted of hypothetical proteins and some had matches in clostridium spp. and bacillus spp.fig. global gene conservation and comparison of functional classification of genes: the figure shows the number of genes and their classification in different cog categories from b. suanatina an4859/ b. hyodysenteriae b78t, b. hyodysenteriae wa and b. intermedia pws/at venn diagram the bar chart shows the distribution of genes in different functional categories. genes were assigned to different categories using blastp against the cog database by comparing the functional distribution of genes in different cog categories, the four strains b. suanatina an4859/ b. hyodysenteriae b78t, b. hyodysenteriae wa and b. intermedia pws/at contained similar numbers of genes for particular functions. in contrast to gene content, gene synteny was highly conserved between b. hyodysenteriae and b. suanatina with only two major genomic rearrangement events that were found in b. hyodysenteriae wa and b. suanatina. some insertions were also found in the b. suanatina genome, when compared with b. hyodysenteriae wa comprising mainly of genes encoding for unknown proteins, bacteriophages and mobile genetic elements. in b. intermedia and b. suanatina, there was a low level of conservation of synteny as compared to the gene content. the alignment of genomic sequences for two isolates of b. suanatina and b. intermedia pws/at using act identified extensive genomic rearrangements, including at least eight major inversions. phylogenetic inference based on the core genome and housekeeping genes in order to infer the evolutionary relationship of b. suanatina with other brachyspira species, phylogenetic analysis was performed on the amino acid sequences of core genomes and nucleotide sequences of housekeeping genes. in total genes were found conserved in all brachyspira species and their phylogenetic tree exhibited a close relationship between b. suanatina and b. intermedia whereas b. hyodysenteriae wa and b. hyodysenteriae b78t formed a distinct clade. the phylogenetic tree is shown in fig.  phylogenetic analysis of housekeeping genes also showed a highly similar tree topology to the core genome analysis and confirmed a clear distinction between b. hyodysenteriae and b. suanatina. the confidence level for the tree topology using bootstrap replicates ranged from 60 % to 100 %.fig. phylogenetic tree of core genome: radial unrooted phylogenetic representation based on concatenated amino acid sequences of core genes of available genomes of genus brachyspira and b. suanatina strain an4859/ the clade containing b. hyodysenteriae, b. intermedia and b. suanatina strains are shown in bold. the tree was constructed using the neighbour joining method. numbers above the branches represent support values obtained from bootstrap replicates. the tree is drawn to scale with branch lengths measured in the number of substitutions per site fig. phylogenetic tree of housekeeping genes: the evolutionary history of housekeeping genes of a selection of strains representing phylum spirochaetes was inferred by the maximum likelihood method based on the general time reversible model. a discrete gamma distribution was used to model evolutionary rate differences among sites. the model allowed some sites to be in variable. the clade containing b. hyodysenteriae, b. intermedia and b. suanatina strains are shown in bold. treponema denticola at and borrelia burgdorferi b31t were used as outgroup. numbers above branches represent support values obtained from bootstrap replicates. the tree is drawn to scale with branch lengths measured in the number of substitutions per site. evolutionary analysis was conducted in mega discussion using the threshold value of 70 % dna-dna similarity for the delineation of bacterial species, as recommended by the ad hoc committee, b. suanatina strain an4859/ does not belong to any of the currently recognized seven brachyspira spp. the highest observed % dna-dna similarities were to the type strains of b. hyodysenteriae and b. intermedia. our finding was further confirmed by the ani values, which are considered to correspond with the % ddh values. since the degree of robustness of anim is higher compared to anib for genomes having >60 % ddh values we calculated the ani values using both algorithms and the values obtained from them are in agreement to each other. furthermore, the evolutionary history inferred from analysis of the 16s rrna gene, housekeeping genes and the core genome suggests a common ancestor of the three related strains b78t, an4859/ and pws/at, which each represent indole test-positive phenotypes of b. hyodysenteriae, b. suanatina and b. intermedia, respectively. high bootstrap values supported the phylogram based on housekeeping genes and the assertion that b. suanatina and b. intermedia make a separate clade that is distinct from the clade formed by b. hyodysenteriae strains. interestingly, the phenotypes of b. hyodysenteriae, b. suanatina and b. intermedia were very similar, with the strong β-haemolysis of the two former being strongly positive vs. the weak β-haemolysis in b. intermedia. notably, in this study, no consistent diagnostic phenotypic feature was found to differentiate b. suanatina from b. hyodysenteriae. the strains of b. suanatina in this study showed varying capabilities of using d-mannitol as growth substrate, whereas strain b78t could not. even though this phenotype could be further investigated it would not qualify as a reliable diagnostic feature. interestingly, among the three strains, the type strain of b. intermedia pws/at seems to have undergone the highest number of changes since its divergence from the common ancestor based on the housekeeping genes phylogram. when comparing the conservation of gene synteny of b. suanatina with b. hyodysenteriae wa and b. intermedia pws/at, very low conservation was found between b. suanatina and b. intermedia as compared to b. suanatina and b. hyodysenteriae wa these results also point to the high level of genomic rearrangements that have taken place in the genome of b. intermedia after these species diverged. previous reports based on multi locus sequence typing have suggested that b. intermedia is a genetically diverse species that may possibly consist of several species. the same applies when strongly β-haemolytic indole-test positive isolates, to which b. hyodysenteriae and b. suanatina and other isolates of unknown species affiliation, are investigated. clearly, this group of brachyspira isolates needs more detailed genomic and phenotypic study to elucidate evolutionary and taxonomic relationships. general genomic features including the distribution of genes in similar functional classes in b. suanatina, b. hyodysenteriae and b. intermedia support the fact that they are highly similar species. however the genome size of b. hyodysenteriae is smaller as compared to b. suanatina and b. intermedia. speculatively, the reason for the smaller genome size could be that b. hyodysenteriae has adapted to the rich environment of porcine intestinal tissues through reductive evolution. reductive evolution is a process of reduction in genome size of a host associated bacterium by the loss of genes rendered non-essential. the presence of intact bacteriophage region bsp in b. suanatina and its conservation in other brachyspira genomes including b. hyodysenteriae b78t while absence in b. hyodysenteriae wa and b. intermedia pws/at also suggests the loss of this prophage in two later species due to the reductive evolution. events of horizontal gene transfer are also evident in the genome of b. suanatina when examining the genes that had no homology in the genomes of brachyspira species with other bacteria. these include the presence of some genes that had their highest similarities in clostridium spp. and bacillus spp. both these species are also inhabitants of the large intestine where the possible exchange of genetic material with b. suanatina could have taken place. events of reductive evolution and hgt in the genomes of brachyspira species have been described earlier. the unusual phage-like gene transfer agent vsh- produced by b. hyodysenteriae is likely involved in natural gene transfer and recombination within different strains of the species. it spans over a  kb region of the genome of b. hyodysenteriae b204r and contains all the late function genes. the region containing genes homologous to vsh- has also been identified in the genomes of different strains of b. hyodysenteriae, b. intermedia, b. pilosicoli and b. murdochii. a region similar to vsh- was also found in the genome of b. suanatina in this study. although it is not clear if this gta in b. suanatina strain an4859/ is functional, the conservation of the late function genes and the presence of unknown genes in the gta region suggest functional capabilities.  this study demonstrates the importance of integrating genomic data into the classification of different bacterial species. based on genomic characteristics we conclude that b. suanatina does not belong to any of the currently validated seven brachyspira species. the use of genomic data from more than one strain of b. suanatina species of the strain will facilitate future comparative genomic studies. description of brachyspira suanatina sp. nov. brachyspira suanatina su.a.na.ti'na. l. masc n. suus swine/pig, l. fem. adj. anatina from/of a duck, n.l. fem. adj. suanatina of pig and duck, referring to the host animals from which this spirochete was isolated. cells are oxygen-tolerant but anaerobic and gram-negative of a typical spirochaetal ultrastructure. cell dimensions are  ±   μm in length and  ±   μm in diameter. the number of periplasmic flagella are 14– per cell inserted subterminally at cells ends. bacterial growth occurs at and 42 °c but not at and 46 °c. when grown anaerobically in bhis at 37 °c cultures reach log-phase at approximately 31 h of inoculation and reach maximum population densities of approximately 2 × 108 cfu/ml at 36–38 h post inoculation. at mid log-phase, the population doubling time is approximately 2 h. cells grow anaerobically on fastidious anaerobe agar with swarming colonies and a strong β-haemolysis on trypticase soy agar. gives positive reactions in tests for tryptophanase activity, alkaline phosphatase, c esterase, c esterase lipase, acid phosphatase, naphtolphosphohydrolase, β-galactosidase. α-glucosidase, β-glucosidase. the following substrates can be utilized as single carbon, nitrogen and energy sources for growth: d-glucose, d-fructose, d-galactose, d-mannose, sucrose, d-trehalose, lactose, n-acetyl-d-glucosamine and pyruvate. d-mannitol sustain growth of some isolates. growth was not detected with d-cellobiose, d-fucose, d-melobiose, d-ribose, d-raffinose, d-xylose, l-arabinose, l-fucose, l-rhamnose, n-acetyl-d-galactosamine, soluble starch, galactosamine, na-acetate, lactate, esculin, succinate, pectin, 50x mem amino acids, cellulose, glutathione, glycerol, glycogen, poly-galacturonic acid, hyaluronic acid and hog gastric mucin. the draft genomic dna g + c content of the type strain is 27 % and the genome comprises,009 bp, and orfs. dna-dna hybridization of the type strain against all type strains of recognized brachyspira spp. and anib/anim values support the taxonomic classification of b. suanatina as a sp. nov. the evolutionary history of 16s rrna and housekeeping genes, and the core genome shows a common ancestor of b. hyodysenteriae, b. intermedia and b. suanatina. a bacteriophage region and a phage-like gta similar to vsh- in b. hyodysenteriae and a putative plasmid sequence related to b. hyodysenteriae strain wa plasmid are present in the genome of b. suanatina strain an4859/ type strain is enteropathogenic in pigs with swine dysentery-like symptoms when using the type strain, and mild diarrhoea was associated with challenge with an isolate of mallard origin. the type strain, an4859/03t was obtained from a swedish pig with swine dysentery-like disease. accession numbers draft genome sequences of b. suanatina an4859/ have been deposited in the european nucleotide archive database with scaffold accession numbers: cvlb01000001-cvlb under the study accession number: prjeb additional files additional file 1: table s strain names and accession numbers of available genome used for the study. additional file 2: table s functional distribution of genes in different cog categories of strains b. suanatina an4859/ b. hyodysenteriae b78t, b. hyodysenteriae wa and b. intermedia pws/at. additional file 3: figure s pairwise genome alignment of b. hyodysenteriae wa b. intermedia pws/at and b. suanatina an4859/ the genomes were aligned using artemis comparison tool. alignment was performed according to the sequence start of b. hyodysenteriae wa and visualised in act with a cut-off set to blast scores >  red bars indicate the similar regions in the same orientation, whereas blue bars show inverted regions. abbreviations actartemis comparison tool aniaverage nucleotide identity bhibrain heart infusion culture broth bhisbrain heart infusion culture broth with calf serum cdscoding sequence cfucolony forming unit cogclusters of orthologous group ddhdna-dna hybridization faafastidious anaerobe agar fcsfoetal calf serum gtagene transfer agent hiheart infusion culture broth hisheart infusion culture broth with foetal calf serum hgthorizontal gene transfer mlmaximum likelihood noxnadh oxidase gene orfopen reading frame competing interests the authors declare that they have no competing interests. authors’ contributions mm, ebr and dsj designed the experiments. mm, sz, dsj and tr conducted experiments and data analyses. mm and dsj drafted the manuscript. all authors contributed to the final manuscript. all authors read and approved the final manuscript. 
the draft genome sequence of b. suanatina strain an4859/ was determined and compared with the available genomes of all valid species of brachyspira. a bacteriophage region and a putative plasmid sequence were also found in the genome of b. suanatina strain an4859/ by comparing the genomes of the three closely related species b. intermedia, b. hyodysenteriae and b. suanatina similar profiles of general genomic features and distribution of genes in different functional categories were obtained. the results of our study suggest that despite being similar to b. hyodysenteriae phenotypically, b. suanatina should be regarded as a separate species based on its genetic characteristics. phylogenetic analysis performed using housekeeping genes and core genomes of all valid brachyspira sp. according to morphological traits, growth characteristics and enzymatic profiles, b. suanatina was similar to the type strain of b. hyodysenteriae, but using the recommended threshold value of 70 % similarity by ddh it did not belong to any of the recognized brachyspira species. hampsonii” revealed that b. suanatina and b. intermedia formed a clade distinct from b. hyodysenteriae. however, the genome size of b. hyodysenteriae was smallest among the species, suggesting the possibility of reductive evolution in the divergence of this species. for this new species we propose the name b. suanatina for which the type strain is an4859/03t. in the current study characterization of b. suanatina was performed on the basis of cell morphology, growth characteristics, enzyme profiles, dna-dna hybridization and whole genome comparisons. nov.core genomedna-dna hybridizationenteropathogenhousekeeping genes phylogenyspirochaetemallardphenotypepigwhole genome sequenceissue-copyright-statement© the author 2015 the genus brachyspira currently encompasses seven valid species that colonize the intestines of mammals and birds. in a previous study a group of strongly haemolytic isolates from pigs and mallards was provisionally described as a new species within genus brachyspira, “b. based on characteristics presented in this report we propose that strains an4859/ an1681:1/ an2384/ and dk12570- from pigs in sweden and denmark, and strains an3949:2/ and an1418:2/ isolated from mallards in sweden, represent a unique species within genus brachyspira.
 the genus brachyspira currently encompasses seven valid species that colonize the intestines of mammals and birds. in a previous study a group of strongly haemolytic isolates from pigs and mallards was provisionally described as a new species within genus brachyspira, “b. suanatina”, and enteropathogenic properties were demonstrated in a porcine challenge model. methods in the current study characterization of b. suanatina was performed on the basis of cell morphology, growth characteristics, enzyme profiles, dna-dna hybridization and whole genome comparisons. the draft genome sequence of b. suanatina strain an4859/ was determined and compared with the available genomes of all valid species of brachyspira.  according to morphological traits, growth characteristics and enzymatic profiles, b. suanatina was similar to the type strain of b. hyodysenteriae, but using the recommended threshold value of 70 % similarity by ddh it did not belong to any of the recognized brachyspira species. this was further supported by average nucleotide identity values. phylogenetic analysis performed using housekeeping genes and core genomes of all valid brachyspira sp. and “b. hampsonii” revealed that b. suanatina and b. intermedia formed a clade distinct from b. hyodysenteriae. by comparing the genomes of the three closely related species b. intermedia, b. hyodysenteriae and b. suanatina similar profiles of general genomic features and distribution of genes in different functional categories were obtained. however, the genome size of b. hyodysenteriae was smallest among the species, suggesting the possibility of reductive evolution in the divergence of this species. a bacteriophage region and a putative plasmid sequence were also found in the genome of b. suanatina strain an4859/  the results of our study suggest that despite being similar to b. hyodysenteriae phenotypically, b. suanatina should be regarded as a separate species based on its genetic characteristics. based on characteristics presented in this report we propose that strains an4859/ an1681:1/ an2384/ and dk12570- from pigs in sweden and denmark, and strains an3949:2/ and an1418:2/ isolated from mallards in sweden, represent a unique species within genus brachyspira. for this new species we propose the name b. suanatina for which the type strain is an4859/03t. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords brachyspira suanatina sp. nov.core genomedna-dna hybridizationenteropathogenhousekeeping genes phylogenyspirochaetemallardphenotypepigwhole genome sequenceissue-copyright-statement© the author 2015 