BACKGROUND
flavobacterium columnare is a gram-negative bacterium belonging to the family flavobacteriaceae   <cit> .
columnaris disease caused by f. columnare represents a continuous threat to the growing aquaculture industry worldwide.
it has been ranked as the second most common disease of the channel catfish  industry in the united states  <cit> .
the bacterium is capable of causing infections in both warm and cold water species of fish  <cit> , and it infects fish species around the world, including carp, channel catfish, goldfish, eel, perch, tilapia, pike perch, rainbow trout, brown trout, salmon, tiger muskellunge and walleye  <cit> .
f. columnare causes epidermal infections affecting gills, skin and fins of the fish, producing either acute or chronic infections, depending on the virulence and genetic group  of the strain, as well as on environmental  <cit>  and host-related factors  <cit> .
f. columnare has high phenotypic homogeneity, therefore strain characterization by standard biochemical tests is not appropriate  <cit> .
however, f. columnare has been divided into three genomovars  using analysis of 16s rdna by restriction-fragment length polymorphism   <cit> .
a recent study further increased the resolution of the method to identify a new genomovar   <cit> .
of these, genomovars i and ii have been reported in europe as either common european  or likely imported   <cit> .
to obtain higher resolution on genetic diversity of f. columnare, several other molecular typing approaches have been used, including single-strand conformation polymorphism   <cit> , amplified fragment length polymorphism   <cit> , pulsed-field gel electrophoresis   <cit> , automated ribosomal intergenic spacer analysis   <cit>  and internal spacer region-single strand conformation polymorphism analysis   <cit> .
although the overall discriminatory power of these methods is high, they can suffer from poor inter-laboratory interpretability, and they are not suitable for population structure studies.
moreover, these genetic markers accumulate genetic variation rapidly, which can interfere with investigating evolutionary phylogenetic relationships or global epidemiology between closely related species of bacteria  <cit>  .
in  <dig>  the complete genome of f. columnare was published  <cit> , making it possible to compare genes from individual isolates for developing multilocus sequence typing  and multilocus sequence analysis  schemes.
mlst/mlsa schemes provide portable, universal, highly discriminatory and unambiguous data .
because this method indexes variation in housekeeping genes that have a relatively slow evolutionary rate, it has been widely used to infer population genetic structure of several different bacterial groups .
the mlst method typically uses variation in four to seven housekeeping gene sequences to characterize isolates of bacterial species.
the allele-based mlst method relies on allelic profiles  and sequence type designations  to estimate relatedness among isolates and so it ignores the number of nucleotide differences between alleles.
the mlsa method, however, relies directly on nucleotide sequences; it uses concatenated sequences of fragments of housekeeping genes to determine genus-wide phylogenetic relationships.
nevertheless, it has been shown that mlsa can also provide robust resolution at the intraspecific level, especially when inadequate phylogenetic resolution prevents mlst from distinguishing phylogenetically closely related strains  <cit> .
since the first columnaris outbreak in finland in  <dig>  f. columnare has been reported as a major threat to salmonid fish farming, particularly rainbow trout   <cit> .
despite its importance as a fish pathogen, the genetic diversity and population structure of f. columnare are poorly known.
the genetic characterization of f. columnare is essential not only to develop appropriate management strategies to minimize the risk of columnaris disease in finnish fish farms but also to better understand host specificity, pathogenicity, and distributional pattern of this bacterium, which is critical for understanding the emergence of columnaris outbreaks worldwide.
this prompted us to develop the first mlst/mlsa scheme for this species in order to investigate the population structure of f. columnare strains isolated from different geographic areas in finland.
materials and methods
bacterial strains and culture conditions
from  <dig> to  <dig>  83 f. columnare strains were obtained from nine different fish species  and water samples , in eight geographic locations in finland comprising four northern and four southern locations .
it is worth noting that all finnish f. columnare isolates have been assigned to genomovar i,  <cit> .
the f. columnare type strain ncimb 2248t isolated in the usa and two reference strains jip39/ <dig> and atcc <dig>  both isolated in france, were also included in the sample collection.
the sequences for strain atcc  <dig> has been retrieved from genbank using their accession numbers.
all the finnish strains were originally isolated using standard culture methods including shieh medium  <cit> , shieh medium supplemented with tobramycin  <cit> , or ao-agar  <cit> .
pure cultures were stored frozen at −80 °c in shieh medium containing 10 % of glycerol and 10 % of fetal calf serum.
before genomic dna extraction, the strains were revived from the stocks in 3 ml of shieh medium for 24 h while shaking  at  <dig>  °c.
overnight cultures were diluted into fresh shieh medium  and allowed to regrow at  <dig>  °c on the shaker for another 24 hours.
dna extraction, pcr amplification and sequencing of housekeeping genes
bacterial cells were harvested from fresh cultures by centrifugation at  <dig> x g for 10 min.
bacterial genomic dna was extracted using the wizard genomic dna purification kit .
the candidate genes for this study  were chosen based on the previously published mlst scheme for flavobacterium psychrophilum by nicolas et al.
<cit> .
we used the corresponding sequences derived from the whole genome sequence of f. columnare atcc  <dig>   <cit>  to redesign the primers for f. columnare.
the primers were designed using primer <dig>   <cit>  and checked for specificity using the ncbi primer-blast tool .
all primers used in this study are listed in additional file  <dig>
pcr reactions were performed in a total volume of 20 μl containing 1x phusion flash pcr master mix ,  <dig>  μm forward primer,  <dig>  μm reverse primer, and 100 ng of genomic dna.
pcr reactions were performed on a biorad 1000c thermal cycler, under the following conditions: 98 °c for 30 s, followed by 30 cycles at 98 °c for 10 s, 62 °c for 20 s, and 72 °c for 15 s, and a final extension at 72 °c for 5 min.
five microliters of the pcr products were run on a  <dig>  % agarose gel to verify correct amplification.
the pcr products were purified using  <dig> u of exonuclease i and  <dig> u of fastaptm thermosensitive alkaline phosphatase  for 15 min at 37 °c, followed by enzyme inactivation for 15 min at 85 °c.
the purified pcr products were then sequenced with the same primers used in amplification using big dye terminator  cycle sequencing kit .
briefly, each 20 μl sequencing reaction mixture contained 2 μl of pcr amplicon,  <dig>  μm of either forward or reverse pcr primer,  <dig>  μl of bigdye ready reaction mix, and  <dig> x sequencing buffer.
the sequencing reaction conditions were as follows: 30 cycles of denaturing at 96 °c for 10 s, annealing at 50 °c for 5 s, and extension at 60 °c for 4 min.
the sequencing products were purified using ethanol/edta/sodium acetate precipitation.
sequencing was performed on an abi 3130xl 16-capillary automated genetic analyzer.
mlst data treatment
the raw sequences were manually inspected and corrected using sequencher  <dig> .
<dig> .
the consensus sequence for each gene was determined by alignment of the forward and reverse sequences using bioedit sequence alignment editor version  <dig> .
<dig>   <cit>  .
the consensus sequences were aligned with clustal w implemented in mega v <dig>   <cit> .
for each housekeeping gene, allele numbers were assigned to unique sequences using the non-redundant databases  program .
for each bacterial isolate, an allelic profile was determined as the combination of alleles at the six loci selected for final analyses.
finally, strains with an identical allelic profile were assigned to the same sequence type  .
phylogenetic analysis
a neighbour-joining dendrogram was constructed using individual and concatenated sequences of the six genes that were selected for final analyses.
for comparison, the f. columnare type strain ncimb 2248t, atcc <dig> and strain jip39/ <dig> were also included in the phylogenetic analysis.
f. psychrophilum jip 02/ <dig> strain from france was used as an outgroup strain.
mega v <dig>   <cit>  was used to evaluate the model for nucleotide substitutions at each protein-coding locus and to construct a phylogenetic tree.
the best model, having the lowest bayesian information criterion  value, was used to generate the neighbour-joining tree based on  <dig> replicates.
the tamura-nei model plus a gamma distribution  model was used to infer the dendrogram for the concatenated sequences.
a chi-square test was employed to evaluate whether the presence of mlsa phylogenetic clusters is explained by regional  separation.
logistic regression was used to explore trends in outbreak likelihood over time for each clade.
moreover, the distribution of mlsa clusters  across sampling origin  were determined using pearson's chi-square statistic, or, where sample sizes were small, fisher's exact test.
a p-value less than or equal to  <dig>  was defined as statistically significant.
computational analysis
analysis of dna sequence variation of the housekeeping gene loci, including the number of alleles for each locus, gc  content, the number of segregating sites , allelic diversity, and the nucleotide diversity , was carried out using dnasp genetic software v <dig> .
<dig>  <cit> .
mega v <dig>  was used to perform tajima’s d test  <cit> .
the start  <dig> package  <cit>  was used to determine the ratios of non-synonymous to synonymous substitutions  for each locus.
the range of intraspecific sequence similarity  for each gene was resolved using bioedit program  <cit> .
population genetic and recombination analyses
evidence for clonal or recombining population structures can be estimated by assessing the level of linkage between alleles at different loci.
to test the null hypothesis, i.e.
whether alleles of the six mlst loci used in the analyses are independent , the ias  values were calculated with the start <dig> program.
the test was performed first for the entire data set of  <dig> isolates and then for only eight sts to avoid biased results due to unequal sample sizes in different sts.
ias values significantly different from  <dig> indicate that a population is clonal , whereas non-significant values indicate a recombining population structure  <cit> .
the concatenated sequence data  were formatted using xmfa2struct .
bayesian approach model with structure version  <dig>   <cit>  was used to determine the levels of inter-lineage recombination and the underlying population structure present in our data.
ten independent runs were performed for each value of the number of ancestral populations  ranging from  <dig> to  <dig>  structure was run for  <dig>  markov chain monte carlo  iterations following  <dig>  burn-in iterations.
the linkage model that reconstructs ancestral populations from dna polymorphism data was used.
the structure harvester  <cit> , which implements the evanno method  <cit> , was used to identify the most probable groups  that best fit the data.
in order to estimate the mutation and recombination rates in f. columnare, we also performed recombination analysis using clonalframe v <dig>   <cit> .
three independent runs of clonalframe were performed, each consisting of  <dig>  mcmc iterations and  <dig>  burn-in iterations.
clonalframe was also used to compare independent runs by the method of gelman and rubin  <cit> .
clonalframe estimates ρ/θ, which measures the relative frequency of occurrence of recombination and mutation in the history of the lineage, and r/m which measures the relative impact of recombination and mutation in the genetic diversification of the lineage.
the values of ρ/θ and r/m for all  <dig> isolates and for each lineage were calculated by extracting the numbers of mutation events, recombination events, and substitutions introduced by recombination from the clonalframe output.
RESULTS
mlst
thirteen housekeeping genes were successfully amplified for the 83 f. columnare isolates from finland and for both the type strain ncimb 2248t and the isolate jip39/ <dig>  the housekeeping genes ftsq, glya, murg, reca, fumc, 16s rdna, and rplb showed identical sequences among all  <dig> finnish strains and were excluded from further analyses.
six remaining housekeeping genes, trpb, tuf, atpa, rpod, gyrb and dnak were selected to develop the mlst/mlsa scheme.
the range of intraspecific sequence similarity  calculated for each gene showed that trpb had the highest level of sequence polymorphism between the strains , followed by rpod .
whereas the lowest levels of inter-strain sequence variation were found for tuf , atpa , gyrb , and dnak .
the gc content for individual genes exhibited very little variation between the strains  .
the average overall gc content for the six genes in all  <dig> finnish strains was 39 %, slightly higher than the overall gc content  for the f. columnare atcc  <dig> genome  <cit> .table  <dig> characteristics of the six loci used in the analyses for the 83 f. columnare strains
gc guanine + cytosine content; s number of segregating sites; dn/ds the ratio of non-synonymous to synonymous substitutions; ∏ nucleotide diversity; tajima’s d neutrality test values
astatistically significant, and bnon-significant value
using the mlst method, the f. columnare strains were divided into  <dig> different sequence types .
the number of alleles at each locus ranged from  <dig>  to  <dig> , and the number of variable sites varied between  <dig>  and  <dig>  .
no insertions or deletions were detected within the loci.
the allelic diversity ranged from  <dig>  for dnak to  <dig>  for trpb, with an average value of  <dig>  .
the average pairwise nucleotide diversity π for all genes was  <dig>  %, with π-values for individual loci ranging from  <dig>  for dnak to  <dig>  for trpb .
tajima’s d values for loci tuf, atpa, dnak, and gyrb were nonsignificant, supporting neutrality of the alleles of these genes  .
in contrast, trpb and rpod both showed a significantly positive value for tajima's d, indicating a decrease in population size, balancing selection  <cit> , or subdivision of the population  <cit> .
phylogenetic analysis
in total,  <dig>  bp of sequence was analyzed for each strain.
the phylogenetic tree based on the concatenated sequence shows two major lineages  which can be further divided into five mlsa clusters .
both lineages are strongly supported by high bootstrap values.
the clusters within lineages associated uniformly with the arisa genotypes of the strains .
therefore, we decided to designate the mlsa clusters in accordance with the corresponding arisa genotypes  <cit> .
lineage i includes two clusters, corresponding to arisa genotypes c and e. cluster c  contains  <dig> strains isolated from six locations: northa, northb, northc, southa, southc, and southd .
cluster e  includes  <dig> strains isolated from southc.
lineage ii includes three clusters, corresponding to arisa genotypes g, a and h. cluster g  contains strains isolated from northb, northc, southb and southc.
cluster a  includes strains isolated from northb, northd, southb, and southc.
cluster h  contains three strains isolated form northb and southc.
statistical analysis showed that there was no significant regional  clustering of f. columnare in finland, but the occurrence of cluster e isolates was more frequent in more recent years .
no other trends were apparent over the study period 2003– <dig> .
cluster h strains  were only found in years  <dig> and  <dig>  which most likely results from random sampling of a very rare genotype, rather than any systematic reasons.
furthermore, no association was found between mlsa clusters and the sampling origin including water and fish species .fig.
<dig> location of sample sites in finland for the 83 f. columnare isolates used in this study
fig.
<dig> phylogenetic tree of f. columnare strains based on a concatenated 6-housekeeping gene .
strict consensus neighbour joining tree for the 83 f. columnare strains from finland and three other strains, including the type strain ncimb 2248t, atcc  <dig>  and strain jip39/ <dig>  the numbers at the nodes represent levels  of bootstrap support from  <dig> resampled datasets.
the sequence of the flavobacterium psychrophilum  strain from france was used as an out-group.
mlsa clusters  c, e, a, h, and g are marked with purple, green, orange, light blue, and pink, respectively.
branches from finnish lineage i strains are coloured red, those from lineage ii are coloured blue, and those from other countries are coloured black
observations for regional sources were compared by nonparametric chi-square test .
observations for year of isolation were tested with logistic regression where year was fitted as a continuous covariate
adue to the small number of cluster h isolates in our dataset , we did not include cluster h in logistic regression
the phylogenetic and bayesian analyses clearly indicated that f. columnare strains are genetically distinct from f. psychrophilum.
in addition, the early-branching type strain  from the usa appears to be highly divergent from the finnish f. columnare strains.
however, the reference strains atcc  <dig> and jip 39/ <dig>  both from france, were phylogenetically close to isolates from finland and clustered with the finnish strains.
population genetic and recombination analysis
results from the clonalframe analysis  showed that mutations have occurred approximately seven times more frequently than recombination .
however, recombination has had a significant effect in the evolutionary process as shown by the r/m value of approximately three .
this indicates that even though recombination has been less frequent than mutation, each recombination event has introduced almost three times more substitutions than mutation.
according to the results extracted from the clonalframe output, the role of recombination seems to be uneven across the two mlsa lineages.
the relative frequency of occurrence of recombination versus mutation  was  <dig>  for lineage i, and  <dig>  for lineage ii.
whereas the relative impact of recombination versus point mutation  was  <dig>  for lineage i and  <dig>  for lineage ii.
however, even in case of high recombination rate, analysis of all  <dig> isolates yielded an ias  value of  <dig>  , indicating persistence of identical clonally expanded sequence types over sampling years.
moreover, ias values calculated separately for lineages i and ii were also significantly different from zero .
when only one representative of each sequence type was included, the overall ias value was decreased to  <dig> , but remained significant , meaning that recombination was not sufficient to break down linkage disequilibrium.
the analysis with structure detected two subpopulations , denoted as lineages i and ii .
lineage ii forms a genetically homogeneous group.
however, lineage i represent a genetically heterogeneous group with genetic material imported from lineage ii.fig.
<dig> estimated f. columnare population structure using structure version  <dig> .
ancestral population size of  <dig>  was chosen as the best fit for the current data.
upper labels indicate the ancestral populations ; sequence types and mlsa clusters of the isolates are listed at the bottom
discussion
in this paper, we report the results of the first mlst/mlsa scheme developed to determine the genetic variability and population genetic structure of the f. columnare strains isolated from different locations in finland.
according to classification based on 16 s rdna-rflp analysis, all finnish isolates studied thus far belong to genomovar i  <cit>  .
this is likely explained by the fact that strains from genomovars ii and iii have higher temperature requirements and thus may not be able to tolerate colder temperatures in finland; it has been shown that these genomovars thrive at higher temperatures but not at 15 °c  <cit> .
therefore, because we focus here on the finnish population, isolates from genomovars ii and iii were not included in the study.
our results demonstrated that finnish f. columnare strains were separated into different sts  and phylogenetic clusters  irrespective of the geographic origin, host species or year of isolation.
previous studies have reported a similar lack of relationship between genotypes and isolation source including both host and geographical origin.
based on three genotyping methods including rdna-rflp, isr and aflp, arias et al.
<cit>  analysed 30 f. columnare isolates from four countries and found that these strains did not cluster according to host species or geographic origin.
similarly, an analysis of ten f. columnare isolates originating from cold and warm waters using rflp analysis showed no relationship between geographic origin and genomic types  <cit> .
one possible explanation for the lack of association between sts and their geographic origin  is that the transportation of fish stocks or equipment between fish farms and natural exchanges such as wild bird could contribute to the spread of f. columnare strains between regions  <cit> .
the six loci used in the analyses  were successfully amplified and sequenced for all isolates.
based on mlst statistical analyses and concatenated phylogenetic analyses, all  <dig> isolates were grouped into eight sts and two major lineages  that were further assigned into five clusters  .
overall, the isolates were commonly designated to cluster c while cluster h strains were the rarest ones.
we found that the environmental isolates of f. columnare are phylogenetically clustered with the epidemic strains, indicating that the environmental population may be the source of epidemic strains and vice versa.
non-synonymous base substitutions in gene sequences were rare  , indicating purifying selection against amino acid changes.
this verifies that the identified sequence variability is selectively neutral at the protein level, making the sequences good candidates for multilocus sequence typing and analysis.
however, an excess of polymorphisms in the trpb and rpod sequences can also indicate balancing selection, supported also by the positive tajima’s d values.
this finding is consistent with two clearly distinct lineages in the finnish sample of f. columnare.
it is also possible that the differences in trpb and rpod are caused by both purifying selection , and potential balancing selection .
the individual trees based on trpb and rpod were found to provide higher phylogenetic resolution than corresponding tuf, atpa, dnak, and gyrb gene sequences .
although cluster h strains could only be resolved by using the concatenated sequences of all six housekeeping gene sequences, the tree constructed from the trpb gene alone could identify four major clusters  resolved by mlsa analysis.
moreover, to determine whether the trpb gene used in the concatenated sequence influenced tree topology, it was compared to the tree constructed on five mlsa gene sequences .
adding the trpb sequence to the mlsa-concatenated sequence  thus had minor effect on the clustering of strains.
furthermore, we found that only one strain, b <dig>  had discordant trpb-based identification .
thereby, we propose that the trpb gene can serve as a reliable, cost-effective, and quick molecular marker to differentiate closely related f. columnare strains.
significant linkage disequilibrium between the f. columnare mlst alleles suggests a non-random distribution of alleles and a clonal population structure where diversity increases by the accumulation of point mutations.
indeed, recovering identical or similar sts from different geographic origins suggests that f. columnare strains may undergo clonal expansion.
this conclusion is also supported by the observation that the same genotypes are repeatedly isolated over the years.
on the other hand, the low ias values  suggest that recombination may have occurred among these strains, but not frequently enough to completely remove linkage disequilibrium.
finding linkage disequilibrium between the alleles is not surprising, considering that an allele must change at least  <dig> times more frequently by recombination than by point mutation to achieve random assortment within a bacterial population  <cit> .
based on both structure and clonalframe analyses, however, recombination events have occurred both within and across lineages.
considering both the relative frequency of occurrence of recombination and mutation  and the relative effect of recombination and mutation in genetic diversification , the finnish f. columnare strains seem to have a moderate recombination rate.
the observed genetic exchange and recombination within and across the two lineages is irreconcilable with a strictly clonal population structure.
therefore, we suggest that the f. columnare population structure in finland follows an epidemic population structure where there is a background level of frequent recombination with consecutive clonal expansion of one or a few fit genotypes  <cit> , similar to that of other bacterial species such as escherichia coli  <cit> , vibrio parahaemolyticus  <cit>  and f. psychrophilum  <cit> .
compared with other fish pathogens, the estimated recombination rate for the finnish f. columnare strains seems to be slightly higher than in tenacibaculum maritimum   <cit> , and is clearly lower than in the highly recombinogenic bacterium f. psychrophilum   <cit> .
however, as our data set consists of local samples, we cannot estimate the recombination rate of f. columnare at a global  level.
our results further demonstrate that lineage i displays more evidence of recombination than lineage ii .
we also show that only lineage i seems to exhibit mixed ancestries  suggesting that limited, unidirectional inter-lineage admixture has taken place .
such infrequent recombination between lineages is also consistent with our results from the phylogenetic analysis, indicating divergence of the two lineages .
a low level of inter-lineage recombination, even among sympatric strains, coupled with high levels of intra-lineage recombination suggests that natural barriers could prevent recombination across the two lineages.
since the statistical analysis indicates that the clustering of the lineages is not caused by geography, the barrier against inter-lineage gene transfer may be caused by other, yet undefined, factors.
CONCLUSIONS
our mlst/mlsa scheme data demonstrate that both recombination and clonality play a role in shaping the population structure of f. columnare in finland.
the population structure of the finnish f. columnare is probably semi-clonal which diversifies with moderate, but variable, recombination rate.
limited association of genotypes with geography or year of isolation indicates that columnaris outbreaks in finland are caused by continuous co-circulation of f. columnare strains.
furthermore, recovering identical genotypes from both fish and from the environment confirms that environmentally persistent strains are also epidemiologically important.
availability of supporting data
dna sequences of all sts have been submitted to the european genetic database embl under the following accession numbers: trpb , rpod , atpa , tuf , gyrb , dnak .
additional files
additional file  <dig> site, year of isolation, source of isolation , location of isolation, host species, sequence type , and allelic profile data for the 83 f. columnare strains from finland analyzed by mlst.
additional file  <dig> target genes and primers for the housekeeping genes of f. columnare.
the loci used for the mlst/mlsa scheme are shown in bold font.
length refers to the length of the target sequence.
* reference for 16s rdna primers  <cit> .
additional file  <dig> neighbor joining phylogenetic trees based on the individual sequences of six mlst loci .
mega v <dig>  was used to evaluate the models for nucleotide substitution for each protein-coding locus and to construct the phylogenetic trees for finnish f. columnare strains.
the f. columnare type strain ncimb 2248t isolated in the usa and two reference strains jip39/ <dig> and atcc <dig>  both isolated in france, were also included in the phylogenetic analysis.
the best model indicated by the lowest bayesian information criterion  value was used to generate the neighbor-joining tree based on  <dig> replicates.
the t <dig> model was selected for dnak, tuf, and gyrb, while t92 + g was selected for both trpb and rpod and jc model was selected for atpa.
strains corresponding to lineage i  are coloured red and strains corresponding to lineage ii  are in blue.
other sequences assigned to reference strains are in black.
few strains representative of each genotype were used for tree construction .
additional file  <dig> the fluorescence peak profiles for the arisa genotypes analysed with abi prism 3130xl genetic analyser and the genemapper v. <dig>  software .
for the  <dig> strains isolated from finland, we also determined the arisa  genotypes following the procedure described by suomalainen et al.
<cit> .
however, the previously published method was modified so that abi prism 3130xl genetic analyser is used instead of li-cor  <dig> automatic sequencer the analysis revealed that arisa genotypes associate uniformly with the clusters from the mlsa scheme.
briefly, the pcr reaction mixture  contained 1x dreamtaq buffer ,  <dig>  mm dntps ,  <dig>  μm reverse primer ,  <dig>  μm forward primer ,  <dig>  μm 6-carboxyfluorescein  labelled forward primer ,  <dig> u of dreamtaq dna polymerase  and 1 μl of template dna.
the pcr reactions were carried out using bio-rad c <dig> or s <dig> thermal cyclers .
the thermo-cycling conditions included initial denaturation at 95 °c for 2 min followed by 30 cycles of amplification  and a final extension at 72 °c for 15 min.
the pcr products were denatured with formamide and genescantm  <dig> liz size standard was added.
products were separated with an abi prism 3130xl genetic analyser, and visualized with genemapper v. <dig>  software .
based on the fluorescence peak profiles and using the strains from suomalainen et al.
<cit>  as positive controls, the 87 f. columnare strains were designated into arisa genotypes a to h  <cit> .
additional file  <dig> phylogenetic tree based on the 16 s rdna sequence data obtained from the representatives of finnish f .columnare genotypes  studied in this study and other f .columnare sequences obtained from the genbank.
the tree was constructed by a upgma clustering method with a resampling of  <dig>  bootstrap replicates and the jukes cantor model.
two strains representative of each arisa genotype/mlsa cluster studied in this study were used for tree construction .
additional file  <dig> the concatenated tree based on six housekeeping mlsa genes  and five mlsa gene sequences .
the arrow shows disagreement concerning the position of strain g <dig>  between the two trees.
abbreviations
mlsamultilocus sequence analysis
mlstmultilocus sequence typing
atccamerican type culture collection
16s rdna-rflp16s rdna by restriction-fragment length polymorphism
sscpsingle-strand conformation polymorphism
aflpamplified fragment length polymorphism
pfgepulsed-field gel electrophoresis
arisaautomated ribosomal intergenic spacer analysis
isr-sscpinternal spacer region-single strand conformation polymorphism analysis
mcmcmarkov chain monte carlo
iasstandardized index of association
dn/dsthe ratios of non-synonymous to synonymous substitutions
t93 + gtamura-nei model and a gamma distribution
bicbayesian information criterion
gcguanine + cytosine content
sthe number of segregating sites
competing interests
the authors declare no competing interests; financial or otherwise.
authors’ contributions
r.a. designed and performed experiments, analysed data and wrote the paper.
k.p.
participated in data collection, and helped to draft the manuscript.
lr.s.
participated in designing the study and helped to draft the manuscript.
n.p.
performed arisa genotyping analysis and helped to draft the manuscript.
t.k.
supervised the project and participated in designing the study, data analysis, and helped to draft the manuscript.

 flavobacterium columnare is a gram-negative bacterium belonging to the family flavobacteriaceae. columnaris disease caused by f. columnare represents a continuous threat to the growing aquaculture industry worldwide. it has been ranked as the second most common disease of the channel catfish industry in the united states. the bacterium is capable of causing infections in both warm and cold water species of fish, and it infects fish species around the world, including carp, channel catfish, goldfish, eel, perch, tilapia, pike perch, rainbow trout, brown trout, salmon, tiger muskellunge and walleye. f. columnare causes epidermal infections affecting gills, skin and fins of the fish, producing either acute or chronic infections, depending on the virulence and genetic group of the strain, as well as on environmental and host-related factors. f. columnare has high phenotypic homogeneity, therefore strain characterization by standard biochemical tests is not appropriate. however, f. columnare has been divided into three genomovars using analysis of 16s rdna by restriction-fragment length polymorphism. a recent study further increased the resolution of the method to identify a new genomovar. of these, genomovars i and ii have been reported in europe as either common european or likely imported. to obtain higher resolution on genetic diversity of f. columnare, several other molecular typing approaches have been used, including single-strand conformation polymorphism, amplified fragment length polymorphism, pulsed-field gel electrophoresis, automated ribosomal intergenic spacer analysis and internal spacer region-single strand conformation polymorphism analysis. although the overall discriminatory power of these methods is high, they can suffer from poor inter-laboratory interpretability, and they are not suitable for population structure studies. moreover, these genetic markers accumulate genetic variation rapidly, which can interfere with investigating evolutionary phylogenetic relationships or global epidemiology between closely related species of bacteria. in the complete genome of f. columnare was published, making it possible to compare genes from individual isolates for developing multilocus sequence typing and multilocus sequence analysis schemes. mlst/mlsa schemes provide portable, universal, highly discriminatory and unambiguous data. because this method indexes variation in housekeeping genes that have a relatively slow evolutionary rate, it has been widely used to infer population genetic structure of several different bacterial groups. the mlst method typically uses variation in four to seven housekeeping gene sequences to characterize isolates of bacterial species. the allele-based mlst method relies on allelic profiles and sequence type designations to estimate relatedness among isolates and so it ignores the number of nucleotide differences between alleles. the mlsa method, however, relies directly on nucleotide sequences; it uses concatenated sequences of fragments of housekeeping genes to determine genus-wide phylogenetic relationships. nevertheless, it has been shown that mlsa can also provide robust resolution at the intraspecific level, especially when inadequate phylogenetic resolution prevents mlst from distinguishing phylogenetically closely related strains. since the first columnaris outbreak in finland in f. columnare has been reported as a major threat to salmonid fish farming, particularly rainbow trout. despite its importance as a fish pathogen, the genetic diversity and population structure of f. columnare are poorly known. the genetic characterization of f. columnare is essential not only to develop appropriate management strategies to minimize the risk of columnaris disease in finnish fish farms but also to better understand host specificity, pathogenicity, and distributional pattern of this bacterium, which is critical for understanding the emergence of columnaris outbreaks worldwide. this prompted us to develop the first mlst/mlsa scheme for this species in order to investigate the population structure of f. columnare strains isolated from different geographic areas in finland. materials and methods bacterial strains and culture conditions from to 83 f. columnare strains were obtained from nine different fish species and water samples, in eight geographic locations in finland comprising four northern and four southern locations. it is worth noting that all finnish f. columnare isolates have been assigned to genomovar i,. the f. columnare type strain ncimb 2248t isolated in the usa and two reference strains jip39/ and atcc both isolated in france, were also included in the sample collection. the sequences for strain atcc has been retrieved from genbank using their accession numbers. all the finnish strains were originally isolated using standard culture methods including shieh medium, shieh medium supplemented with tobramycin, or ao-agar. pure cultures were stored frozen at −80 °c in shieh medium containing 10 % of glycerol and 10 % of fetal calf serum. before genomic dna extraction, the strains were revived from the stocks in 3 ml of shieh medium for 24 h while shaking at  °c. overnight cultures were diluted into fresh shieh medium and allowed to regrow at  °c on the shaker for another 24 hours. dna extraction, pcr amplification and sequencing of housekeeping genes bacterial cells were harvested from fresh cultures by centrifugation at x g for 10 min. bacterial genomic dna was extracted using the wizard genomic dna purification kit. the candidate genes for this study were chosen based on the previously published mlst scheme for flavobacterium psychrophilum by nicolas et al.. we used the corresponding sequences derived from the whole genome sequence of f. columnare atcc to redesign the primers for f. columnare. the primers were designed using primer and checked for specificity using the ncbi primer-blast tool. all primers used in this study are listed in additional file pcr reactions were performed in a total volume of 20 μl containing 1x phusion flash pcr master mix,  μm forward primer,  μm reverse primer, and 100 ng of genomic dna. pcr reactions were performed on a biorad 1000c thermal cycler, under the following conditions: 98 °c for 30 s, followed by 30 cycles at 98 °c for 10 s, 62 °c for 20 s, and 72 °c for 15 s, and a final extension at 72 °c for 5 min. five microliters of the pcr products were run on a  % agarose gel to verify correct amplification. the pcr products were purified using u of exonuclease i and u of fastaptm thermosensitive alkaline phosphatase for 15 min at 37 °c, followed by enzyme inactivation for 15 min at 85 °c. the purified pcr products were then sequenced with the same primers used in amplification using big dye terminator cycle sequencing kit. briefly, each 20 μl sequencing reaction mixture contained 2 μl of pcr amplicon,  μm of either forward or reverse pcr primer,  μl of bigdye ready reaction mix, and x sequencing buffer. the sequencing reaction conditions were as follows: 30 cycles of denaturing at 96 °c for 10 s, annealing at 50 °c for 5 s, and extension at 60 °c for 4 min. the sequencing products were purified using ethanol/edta/sodium acetate precipitation. sequencing was performed on an abi 3130xl 16-capillary automated genetic analyzer. mlst data treatment the raw sequences were manually inspected and corrected using sequencher.. the consensus sequence for each gene was determined by alignment of the forward and reverse sequences using bioedit sequence alignment editor version.. the consensus sequences were aligned with clustal w implemented in mega v. for each housekeeping gene, allele numbers were assigned to unique sequences using the non-redundant databases program. for each bacterial isolate, an allelic profile was determined as the combination of alleles at the six loci selected for final analyses. finally, strains with an identical allelic profile were assigned to the same sequence type. phylogenetic analysis a neighbour-joining dendrogram was constructed using individual and concatenated sequences of the six genes that were selected for final analyses. for comparison, the f. columnare type strain ncimb 2248t, atcc and strain jip39/ were also included in the phylogenetic analysis. f. psychrophilum jip 02/ strain from france was used as an outgroup strain. mega v was used to evaluate the model for nucleotide substitutions at each protein-coding locus and to construct a phylogenetic tree. the best model, having the lowest bayesian information criterion value, was used to generate the neighbour-joining tree based on replicates. the tamura-nei model plus a gamma distribution model was used to infer the dendrogram for the concatenated sequences. a chi-square test was employed to evaluate whether the presence of mlsa phylogenetic clusters is explained by regional separation. logistic regression was used to explore trends in outbreak likelihood over time for each clade. moreover, the distribution of mlsa clusters across sampling origin were determined using pearson's chi-square statistic, or, where sample sizes were small, fisher's exact test. a p-value less than or equal to was defined as statistically significant. computational analysis analysis of dna sequence variation of the housekeeping gene loci, including the number of alleles for each locus, gc content, the number of segregating sites, allelic diversity, and the nucleotide diversity, was carried out using dnasp genetic software v.. mega v was used to perform tajima’s d test. the start package was used to determine the ratios of non-synonymous to synonymous substitutions for each locus. the range of intraspecific sequence similarity for each gene was resolved using bioedit program. population genetic and recombination analyses evidence for clonal or recombining population structures can be estimated by assessing the level of linkage between alleles at different loci. to test the null hypothesis, i.e. whether alleles of the six mlst loci used in the analyses are independent, the ias values were calculated with the start program. the test was performed first for the entire data set of isolates and then for only eight sts to avoid biased results due to unequal sample sizes in different sts. ias values significantly different from indicate that a population is clonal, whereas non-significant values indicate a recombining population structure. the concatenated sequence data were formatted using xmfa2struct. bayesian approach model with structure version was used to determine the levels of inter-lineage recombination and the underlying population structure present in our data. ten independent runs were performed for each value of the number of ancestral populations ranging from to structure was run for markov chain monte carlo iterations following burn-in iterations. the linkage model that reconstructs ancestral populations from dna polymorphism data was used. the structure harvester, which implements the evanno method, was used to identify the most probable groups that best fit the data. in order to estimate the mutation and recombination rates in f. columnare, we also performed recombination analysis using clonalframe v. three independent runs of clonalframe were performed, each consisting of mcmc iterations and burn-in iterations. clonalframe was also used to compare independent runs by the method of gelman and rubin. clonalframe estimates ρ/θ, which measures the relative frequency of occurrence of recombination and mutation in the history of the lineage, and r/m which measures the relative impact of recombination and mutation in the genetic diversification of the lineage. the values of ρ/θ and r/m for all isolates and for each lineage were calculated by extracting the numbers of mutation events, recombination events, and substitutions introduced by recombination from the clonalframe output.  mlst thirteen housekeeping genes were successfully amplified for the 83 f. columnare isolates from finland and for both the type strain ncimb 2248t and the isolate jip39/ the housekeeping genes ftsq, glya, murg, reca, fumc, 16s rdna, and rplb showed identical sequences among all finnish strains and were excluded from further analyses. six remaining housekeeping genes, trpb, tuf, atpa, rpod, gyrb and dnak were selected to develop the mlst/mlsa scheme. the range of intraspecific sequence similarity calculated for each gene showed that trpb had the highest level of sequence polymorphism between the strains, followed by rpod. whereas the lowest levels of inter-strain sequence variation were found for tuf, atpa, gyrb, and dnak. the gc content for individual genes exhibited very little variation between the strains. the average overall gc content for the six genes in all finnish strains was 39 %, slightly higher than the overall gc content for the f. columnare atcc genome.table characteristics of the six loci used in the analyses for the 83 f. columnare strains gc guanine + cytosine content; s number of segregating sites; dn/ds the ratio of non-synonymous to synonymous substitutions; ∏ nucleotide diversity; tajima’s d neutrality test values astatistically significant, and bnon-significant value using the mlst method, the f. columnare strains were divided into different sequence types. the number of alleles at each locus ranged from to, and the number of variable sites varied between and. no insertions or deletions were detected within the loci. the allelic diversity ranged from for dnak to for trpb, with an average value of. the average pairwise nucleotide diversity π for all genes was  %, with π-values for individual loci ranging from for dnak to for trpb. tajima’s d values for loci tuf, atpa, dnak, and gyrb were nonsignificant, supporting neutrality of the alleles of these genes. in contrast, trpb and rpod both showed a significantly positive value for tajima's d, indicating a decrease in population size, balancing selection, or subdivision of the population. phylogenetic analysis in total,  bp of sequence was analyzed for each strain. the phylogenetic tree based on the concatenated sequence shows two major lineages which can be further divided into five mlsa clusters. both lineages are strongly supported by high bootstrap values. the clusters within lineages associated uniformly with the arisa genotypes of the strains. therefore, we decided to designate the mlsa clusters in accordance with the corresponding arisa genotypes. lineage i includes two clusters, corresponding to arisa genotypes c and e. cluster c contains strains isolated from six locations: northa, northb, northc, southa, southc, and southd. cluster e includes strains isolated from southc. lineage ii includes three clusters, corresponding to arisa genotypes g, a and h. cluster g contains strains isolated from northb, northc, southb and southc. cluster a includes strains isolated from northb, northd, southb, and southc. cluster h contains three strains isolated form northb and southc. statistical analysis showed that there was no significant regional clustering of f. columnare in finland, but the occurrence of cluster e isolates was more frequent in more recent years. no other trends were apparent over the study period 2003–. cluster h strains were only found in years and which most likely results from random sampling of a very rare genotype, rather than any systematic reasons. furthermore, no association was found between mlsa clusters and the sampling origin including water and fish species.fig. location of sample sites in finland for the 83 f. columnare isolates used in this study fig. phylogenetic tree of f. columnare strains based on a concatenated 6-housekeeping gene. strict consensus neighbour joining tree for the 83 f. columnare strains from finland and three other strains, including the type strain ncimb 2248t, atcc and strain jip39/ the numbers at the nodes represent levels of bootstrap support from resampled datasets. the sequence of the flavobacterium psychrophilum strain from france was used as an out-group. mlsa clusters c, e, a, h, and g are marked with purple, green, orange, light blue, and pink, respectively. branches from finnish lineage i strains are coloured red, those from lineage ii are coloured blue, and those from other countries are coloured black observations for regional sources were compared by nonparametric chi-square test. observations for year of isolation were tested with logistic regression where year was fitted as a continuous covariate adue to the small number of cluster h isolates in our dataset, we did not include cluster h in logistic regression the phylogenetic and bayesian analyses clearly indicated that f. columnare strains are genetically distinct from f. psychrophilum. in addition, the early-branching type strain from the usa appears to be highly divergent from the finnish f. columnare strains. however, the reference strains atcc and jip 39/ both from france, were phylogenetically close to isolates from finland and clustered with the finnish strains. population genetic and recombination analysis results from the clonalframe analysis showed that mutations have occurred approximately seven times more frequently than recombination. however, recombination has had a significant effect in the evolutionary process as shown by the r/m value of approximately three. this indicates that even though recombination has been less frequent than mutation, each recombination event has introduced almost three times more substitutions than mutation. according to the results extracted from the clonalframe output, the role of recombination seems to be uneven across the two mlsa lineages. the relative frequency of occurrence of recombination versus mutation was for lineage i, and for lineage ii. whereas the relative impact of recombination versus point mutation was for lineage i and for lineage ii. however, even in case of high recombination rate, analysis of all isolates yielded an ias value of, indicating persistence of identical clonally expanded sequence types over sampling years. moreover, ias values calculated separately for lineages i and ii were also significantly different from zero. when only one representative of each sequence type was included, the overall ias value was decreased to, but remained significant, meaning that recombination was not sufficient to break down linkage disequilibrium. the analysis with structure detected two subpopulations, denoted as lineages i and ii. lineage ii forms a genetically homogeneous group. however, lineage i represent a genetically heterogeneous group with genetic material imported from lineage ii.fig. estimated f. columnare population structure using structure version. ancestral population size of was chosen as the best fit for the current data. upper labels indicate the ancestral populations; sequence types and mlsa clusters of the isolates are listed at the bottom discussion in this paper, we report the results of the first mlst/mlsa scheme developed to determine the genetic variability and population genetic structure of the f. columnare strains isolated from different locations in finland. according to classification based on 16 s rdna-rflp analysis, all finnish isolates studied thus far belong to genomovar i. this is likely explained by the fact that strains from genomovars ii and iii have higher temperature requirements and thus may not be able to tolerate colder temperatures in finland; it has been shown that these genomovars thrive at higher temperatures but not at 15 °c. therefore, because we focus here on the finnish population, isolates from genomovars ii and iii were not included in the study. our results demonstrated that finnish f. columnare strains were separated into different sts and phylogenetic clusters irrespective of the geographic origin, host species or year of isolation. previous studies have reported a similar lack of relationship between genotypes and isolation source including both host and geographical origin. based on three genotyping methods including rdna-rflp, isr and aflp, arias et al. analysed 30 f. columnare isolates from four countries and found that these strains did not cluster according to host species or geographic origin. similarly, an analysis of ten f. columnare isolates originating from cold and warm waters using rflp analysis showed no relationship between geographic origin and genomic types. one possible explanation for the lack of association between sts and their geographic origin is that the transportation of fish stocks or equipment between fish farms and natural exchanges such as wild bird could contribute to the spread of f. columnare strains between regions. the six loci used in the analyses were successfully amplified and sequenced for all isolates. based on mlst statistical analyses and concatenated phylogenetic analyses, all isolates were grouped into eight sts and two major lineages that were further assigned into five clusters. overall, the isolates were commonly designated to cluster c while cluster h strains were the rarest ones. we found that the environmental isolates of f. columnare are phylogenetically clustered with the epidemic strains, indicating that the environmental population may be the source of epidemic strains and vice versa. non-synonymous base substitutions in gene sequences were rare, indicating purifying selection against amino acid changes. this verifies that the identified sequence variability is selectively neutral at the protein level, making the sequences good candidates for multilocus sequence typing and analysis. however, an excess of polymorphisms in the trpb and rpod sequences can also indicate balancing selection, supported also by the positive tajima’s d values. this finding is consistent with two clearly distinct lineages in the finnish sample of f. columnare. it is also possible that the differences in trpb and rpod are caused by both purifying selection, and potential balancing selection. the individual trees based on trpb and rpod were found to provide higher phylogenetic resolution than corresponding tuf, atpa, dnak, and gyrb gene sequences. although cluster h strains could only be resolved by using the concatenated sequences of all six housekeeping gene sequences, the tree constructed from the trpb gene alone could identify four major clusters resolved by mlsa analysis. moreover, to determine whether the trpb gene used in the concatenated sequence influenced tree topology, it was compared to the tree constructed on five mlsa gene sequences. adding the trpb sequence to the mlsa-concatenated sequence thus had minor effect on the clustering of strains. furthermore, we found that only one strain, b had discordant trpb-based identification. thereby, we propose that the trpb gene can serve as a reliable, cost-effective, and quick molecular marker to differentiate closely related f. columnare strains. significant linkage disequilibrium between the f. columnare mlst alleles suggests a non-random distribution of alleles and a clonal population structure where diversity increases by the accumulation of point mutations. indeed, recovering identical or similar sts from different geographic origins suggests that f. columnare strains may undergo clonal expansion. this conclusion is also supported by the observation that the same genotypes are repeatedly isolated over the years. on the other hand, the low ias values suggest that recombination may have occurred among these strains, but not frequently enough to completely remove linkage disequilibrium. finding linkage disequilibrium between the alleles is not surprising, considering that an allele must change at least times more frequently by recombination than by point mutation to achieve random assortment within a bacterial population. based on both structure and clonalframe analyses, however, recombination events have occurred both within and across lineages. considering both the relative frequency of occurrence of recombination and mutation and the relative effect of recombination and mutation in genetic diversification, the finnish f. columnare strains seem to have a moderate recombination rate. the observed genetic exchange and recombination within and across the two lineages is irreconcilable with a strictly clonal population structure. therefore, we suggest that the f. columnare population structure in finland follows an epidemic population structure where there is a background level of frequent recombination with consecutive clonal expansion of one or a few fit genotypes, similar to that of other bacterial species such as escherichia coli, vibrio parahaemolyticus and f. psychrophilum. compared with other fish pathogens, the estimated recombination rate for the finnish f. columnare strains seems to be slightly higher than in tenacibaculum maritimum, and is clearly lower than in the highly recombinogenic bacterium f. psychrophilum. however, as our data set consists of local samples, we cannot estimate the recombination rate of f. columnare at a global level. our results further demonstrate that lineage i displays more evidence of recombination than lineage ii. we also show that only lineage i seems to exhibit mixed ancestries suggesting that limited, unidirectional inter-lineage admixture has taken place. such infrequent recombination between lineages is also consistent with our results from the phylogenetic analysis, indicating divergence of the two lineages. a low level of inter-lineage recombination, even among sympatric strains, coupled with high levels of intra-lineage recombination suggests that natural barriers could prevent recombination across the two lineages. since the statistical analysis indicates that the clustering of the lineages is not caused by geography, the barrier against inter-lineage gene transfer may be caused by other, yet undefined, factors.  our mlst/mlsa scheme data demonstrate that both recombination and clonality play a role in shaping the population structure of f. columnare in finland. the population structure of the finnish f. columnare is probably semi-clonal which diversifies with moderate, but variable, recombination rate. limited association of genotypes with geography or year of isolation indicates that columnaris outbreaks in finland are caused by continuous co-circulation of f. columnare strains. furthermore, recovering identical genotypes from both fish and from the environment confirms that environmentally persistent strains are also epidemiologically important. availability of supporting data dna sequences of all sts have been submitted to the european genetic database embl under the following accession numbers: trpb, rpod, atpa, tuf, gyrb, dnak. additional files additional file site, year of isolation, source of isolation, location of isolation, host species, sequence type, and allelic profile data for the 83 f. columnare strains from finland analyzed by mlst. additional file target genes and primers for the housekeeping genes of f. columnare. the loci used for the mlst/mlsa scheme are shown in bold font. length refers to the length of the target sequence. * reference for 16s rdna primers. additional file neighbor joining phylogenetic trees based on the individual sequences of six mlst loci. mega v was used to evaluate the models for nucleotide substitution for each protein-coding locus and to construct the phylogenetic trees for finnish f. columnare strains. the f. columnare type strain ncimb 2248t isolated in the usa and two reference strains jip39/ and atcc both isolated in france, were also included in the phylogenetic analysis. the best model indicated by the lowest bayesian information criterion value was used to generate the neighbor-joining tree based on replicates. the t model was selected for dnak, tuf, and gyrb, while t92 + g was selected for both trpb and rpod and jc model was selected for atpa. strains corresponding to lineage i are coloured red and strains corresponding to lineage ii are in blue. other sequences assigned to reference strains are in black. few strains representative of each genotype were used for tree construction. additional file the fluorescence peak profiles for the arisa genotypes analysed with abi prism 3130xl genetic analyser and the genemapper v. software. for the strains isolated from finland, we also determined the arisa genotypes following the procedure described by suomalainen et al.. however, the previously published method was modified so that abi prism 3130xl genetic analyser is used instead of li-cor automatic sequencer the analysis revealed that arisa genotypes associate uniformly with the clusters from the mlsa scheme. briefly, the pcr reaction mixture contained 1x dreamtaq buffer,  mm dntps,  μm reverse primer,  μm forward primer,  μm 6-carboxyfluorescein labelled forward primer, u of dreamtaq dna polymerase and 1 μl of template dna. the pcr reactions were carried out using bio-rad c or s thermal cyclers. the thermo-cycling conditions included initial denaturation at 95 °c for 2 min followed by 30 cycles of amplification and a final extension at 72 °c for 15 min. the pcr products were denatured with formamide and genescantm liz size standard was added. products were separated with an abi prism 3130xl genetic analyser, and visualized with genemapper v. software. based on the fluorescence peak profiles and using the strains from suomalainen et al. as positive controls, the 87 f. columnare strains were designated into arisa genotypes a to h. additional file phylogenetic tree based on the 16 s rdna sequence data obtained from the representatives of finnish f.columnare genotypes studied in this study and other f.columnare sequences obtained from the genbank. the tree was constructed by a upgma clustering method with a resampling of bootstrap replicates and the jukes cantor model. two strains representative of each arisa genotype/mlsa cluster studied in this study were used for tree construction. additional file the concatenated tree based on six housekeeping mlsa genes and five mlsa gene sequences. the arrow shows disagreement concerning the position of strain g between the two trees. abbreviations mlsamultilocus sequence analysis mlstmultilocus sequence typing atccamerican type culture collection 16s rdna-rflp16s rdna by restriction-fragment length polymorphism sscpsingle-strand conformation polymorphism aflpamplified fragment length polymorphism pfgepulsed-field gel electrophoresis arisaautomated ribosomal intergenic spacer analysis isr-sscpinternal spacer region-single strand conformation polymorphism analysis mcmcmarkov chain monte carlo iasstandardized index of association dn/dsthe ratios of non-synonymous to synonymous substitutions t93 + gtamura-nei model and a gamma distribution bicbayesian information criterion gcguanine + cytosine content sthe number of segregating sites competing interests the authors declare no competing interests; financial or otherwise. authors’ contributions r.a. designed and performed experiments, analysed data and wrote the paper. k.p. participated in data collection, and helped to draft the manuscript. lr.s. participated in designing the study and helped to draft the manuscript. n.p. performed arisa genotyping analysis and helped to draft the manuscript. t.k. supervised the project and participated in designing the study, data analysis, and helped to draft the manuscript. 
columnaris disease caused by flavobacterium columnare is a serious problem in aquaculture, annually causing large economic losses around the world. our study with reproducible methodology and comparable results establishes a robust framework for the discrimination and phylogenetic analysis of f. columnare isolates, which will help to improve our understanding about geographic distribution and epidemiology of columnaris disease. partial sequencing of six housekeeping genes revealed eight sequence types and a moderate level of genetic diversity. phylogenetic analysis of the concatenated protein-encoding gene sequence data formed two lineages, which could be further divided into five clusters. all analysed f. columnare strains appeared to have a genetic origin distinct from that of another important fish pathogen form the genus flavobacterium, f. psychrophilum. we investigated the population structure and spatiotemporal changes in the genetic diversity of f. columnare population in finland by using a multilocus sequence typing and analysis based on dna sequence variation within six housekeeping genes. the results suggest that finnish f. columnare strains have an epidemic population structure followed by clonal expansion of successful genotypes. despite considerable research, the molecular epidemiology of f. columnare remains poorly understood. a total of strains of f. columnare were collected from eight different areas located across the country between and flavobacterium columnaremlst/mlsa schemeclonalframerecombination rateclonalityissue-copyright-statement© the author 2015
 columnaris disease caused by flavobacterium columnare is a serious problem in aquaculture, annually causing large economic losses around the world. despite considerable research, the molecular epidemiology of f. columnare remains poorly understood. methods we investigated the population structure and spatiotemporal changes in the genetic diversity of f. columnare population in finland by using a multilocus sequence typing and analysis based on dna sequence variation within six housekeeping genes. a total of strains of f. columnare were collected from eight different areas located across the country between and  partial sequencing of six housekeeping genes revealed eight sequence types and a moderate level of genetic diversity. phylogenetic analysis of the concatenated protein-encoding gene sequence data formed two lineages, which could be further divided into five clusters. all analysed f. columnare strains appeared to have a genetic origin distinct from that of another important fish pathogen form the genus flavobacterium, f. psychrophilum. although the value of the index of association between alleles,, supports some degree of clonality for this species in finland, recombination has introduced molecular diversity to the population almost three times more than mutation.  the results suggest that finnish f. columnare strains have an epidemic population structure followed by clonal expansion of successful genotypes. our study with reproducible methodology and comparable results establishes a robust framework for the discrimination and phylogenetic analysis of f. columnare isolates, which will help to improve our understanding about geographic distribution and epidemiology of columnaris disease. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords flavobacterium columnaremlst/mlsa schemeclonalframerecombination rateclonalityissue-copyright-statement© the author 2015 