BACKGROUND
the intensive culturing of fish in artificially high population densities facilitates disease outbreaks  <cit> .
standard protocols for traditional fish health management, controlling and preventing diseases, have focused on vaccines, operational prophylactic measures and oral treatment using therapeutic agents  <cit> .
however, in those cases where vaccines or treatments are not available, understanding the epizootiology becomes the key to prevent outbreaks and pathogen dispersal  <cit> .
using molecular tools for strain typing of pathogens, combined with biological and ecological knowledge of both the pathogen and the host, is a prerequisite in creating an epizootiological understanding which can be applied in management of diseases in both wild and cultured populations.
historically, the atlantic cod  has been an important marine resource in norway, occurring naturally both as migratory and stationary populations  <cit> .
these populations are divided into coastal and oceanic populations, based on phenotypic and genotypic traits.
the norwegian coastal cod population is considered to consist of several stationary sub populations dispersed along the norwegian coast  <cit> .
the oceanic cod in norwegian waters can be divided into two populations, the north east arctic cod and the north sea cod, both with seasonally dependent migratory behavior.
the oceanic and some of the coastal populations will, at certain times of the year, be present in at the same locations.
during the last decade cod has been a species of increasing significance and interest for the norwegian aquaculture industry  <cit> .
cod is now intensively cultured with full control of all life-stages, except for the broodfish which are predominantly of wild origin  <cit> .
the intensive production cycle of cod consists of three separate steps, where the first step starts with the fertilization of eggs from broodstock held in large land based tanks.
the majority of broodfish consists of wild caught cod from both costal and oceanic populations.
after hatching the fry is held in indoor tanks.
the second stage is on-growth, where the fingerlings are transferred to net pens in the sea and held until they are moved to production sites.
mixing of populations from different broodstock companies occurs at the on-growth sites, i.e.
before redistribution to production sites.
mixing of populations may also occur at the production sites due to limited availability of fish from the different on-growing sites.
several generations can be present at one production site.
the different operational sites are dispersed along the norwegian coast resulting in large scale movement of cultured cod between different parts of norway  <cit> .
in  <dig> a francisella species was isolated from farmed cod in norway showing clinical signs of a chronic granulomatous infection in kidney, spleen, liver and heart  <cit> .
the bacterium was initially characterized both as a species, f. piscicida   <cit> , and later in the same year as a f. philomiragia subspecies, ssp.
noatunensis   <cit> .
the f. philomiragia ssp.
noatunensis was later elevated to the rank of species, f. noatunensis ssp.
noatunensis, with f. piscicida as a heterotypic synonym  <cit> .
following the discovery, several annual outbreaks of francisellosis have been diagnosed with main foci in the western parts of norway  <cit> .
the bacterium has been detected in cultured fish from rogaland to nordland county , however, positive wild cod has only been identified south of sogn og fjordane county , indicating a southern natural reservoir of the bacteria  <cit> .
after the discovery of f. noatunensis ssp.
noatunensis, histological material collected in  <dig> from a broodfish population in hordaland county, has been stained positive with antiserum for the f. noatunensis ssp.
noatunensis gm <dig> isolate, suggesting the presence of the bacteria in norwegian waters prior to the initial discovery  <cit> .
during the last  <dig> years members of the genus francisella have emerged as a global problem for aquaculture, causing mortality among a wide range of aquatic hosts.
these francisella isolates show high genetic similarities with f. philomiragia and with each other at the 16s rrna-gene  <cit> .
the close genetic relatedness has been confirmed with sequencing of several housekeeping genes showing identical sequences for norwegian isolates  <cit> .
a chilean isolate from atlantic salmon  ua <dig>  <cit>  showed a high similarity to the norwegian cod isolates  and was therefore described as a new isolate of the f. noatunensis ssp.
noatunensis  <cit> .
fish pathogenic francisella isolates from asia display high genetic similarity to ua <dig> and ncimb 14265t/gm <dig> and constitute a separate subspecies, f. noatunensis ssp.
orientalis   <cit> .
f. asiatica  isolated from tilapia  in costa rica is identical with f. noatunensis ssp.
orientalis  with respect to phenotype and signature sequences   <cit> .
isolates of f. noatunensis ssp.
noatunensis from norway are identical when comparing 16s and housekeeping gene sequences  <cit> , making genotyping and studies of epizootics difficult.
if an epizootiological approach is to be applied for controlling and preventing dispersal of francisella spp.
in aquaculture, a tool for the identification of isolates is a necessity.
such a tool will make it possible to differentiate between wild endemic and anthropogenic dispersed strains.
this approach also requires knowledge of the production history of cod and some knowledge of natural occurring strains in wild cod in the production areas.
targeting the housekeeping genes of bacterial fish pathogens has been applied for strain identification  <cit> , but due to identical signature sequences this has not been possible to use on f. noatunensis  <cit> .
similar problems are known from epidemiological studies of human pathogenic bacteria like the mycobacterium tuberculosis, m. avium ssp.
paratuberculosis, bacillus anthracis, yersinia pestis and the related f. tularensis  <cit>  where evolutionary or stochastic events have led to a dominance of highly fit clones making strain identification difficult  <cit> .
by using genetic markers with higher mutational rates, like variable number of tandem repeat , one can construct allelic profiles in multiple locus vntr analysis  systems, making strain differentiation of clonal bacteria populations possible  <cit> .
mlva systems have already been applied in studies of f. tularensis isolates  <cit> .
this study is an attempt to provide a sequenced based mlva tool for studying the epizootics of f. noatunensis ssp.
noatunensis in norwegian cod farming.
our findings suggest that the mlva system presented in this study is suitable for strain typing of f. noatunensis isolates from norwegian cod, and the results indicate that there are only a few clades causing francisellosis outbreaks in norwegian cod culturing.
methods
isolation of francisella noatunensis ssp.
noatunensis
all isolates of francisella noatunensis ssp.
noatunensis  included in this study were obtained from atlantic cod  suffering from francisellosis, with the exception of the type strain ncimb 14265t  <cit>  and the chilean atlantic salmon  isolate ua <dig> which were obtained from the national collection of industrial, food and marine bacteria  and intervet norbio as, respectively .
the majority of isolates were collected from cultured populations of cod geographically dispersed along  <dig> km of norwegian coast line during the period from  <dig> to  <dig>  four isolates were obtained from wild-cod, from the counties of aust-agder, vest-agder, rogaland and from the borders of rogaland and hordaland .
an overview of geographical origin of isolates is presented in figure  <dig>  isolates of closely related francisella spp., f. notunensis ssp.
orientalis  and f. philomiragia , were also included in this study.
isolate foj-001/ <dig>  was kindly provided by dr. kamaishi, and isolates foi-002/ <dig>  foi-003/ <dig> and foi-004/ <dig> from indonesia were kindly provided by intervet singapore.
most f. philomiragia isolates were obtained from the culture collections deutsche sammlung von mikroorganismen und zellkulturen  and culture collection, university of göteborg , while strains  <dig> and strain  <dig> <cit>  were kindly provided by dr. berrada and dr. friis-møller, respectively.
a complete overview of isolates and year of isolation, host and location is presented in table  <dig>
the table show source and geographical location of isolation for the francisella isolates  included in the current study.
norway is country of origin unless otherwise is stated.
isolates were kindly provided by intervet norbio asa, dr. kamaishib, intervet singaporec, dr. friis-møllerd and dr. berradae
isolated from broodfishf, isolated in chilea, w = wild caught, t = type strain
culture of f. noatunensis ssp.
noatunensis
isolations of f. noatunensis ssp.
noatunensis were performed in the field or in our laboratory by streaking tissue  from infected cod on agar plates.
initial isolations were performed using blood agar plates with  <dig> % cysteine and 1% glucose as described by nylund et al.
with subsequent growth in b <dig> broth for cryopreservation  <cit> .
however, most isolations were performed using cysteine heart agar  supplemented with 5% bovine blood  as described in olsen et al.
with slight modifications.
these included the addition of  <dig> μg ml- <dig> of ampicillin, 50- <dig> μg ml- <dig> of fungizone® and additional cysteine to a final concentration of  <dig> %.
the ampicillin, fungizone®, bovine blood and cysteine were added to the chab at 60°c.
all isolates were incubated at 20°c for 2- <dig> days  or two weeks .
a sanyo mir- <dig> incubator was used for solid cultures whereas liquid culturing was performed in an infors unitron incubator at  <dig> rpm.
colonies from chab agar were dissolved in eugon broth  for cryopreservation and stored in liquid nitrogen.
genomic analysis
a preliminary genomic sequence from francisella noatunensis ssp.
noatunensis isolate gm <dig>  was analyzed for the presence of tandem repeat regions using the software tandem repeats finder  <cit> .
this program locates and displays tandem repeats in dna sequences.
loci with tandem repeats consisting of less than  <dig> nucleotides and more than four repeats were selected .
the preliminary genomic sequence of f. noatunensis ssp.
noatunensis isolate gm <dig> was generated using the  <dig> pyro-sequencing technology  <cit>  in co-operation with intervet as/asa.
genbank accession numbers refer to sequences of the f. noatunensis ssp.
noatunensis gm <dig> isolate.
repeat sequence, number of repeats and amplicon size is that of the gm <dig> isolate.
not amplified within f. philomiragia isolatesna
dna extraction
total dna from all collected strains and tissues was extracted using the dneasy® blood & tissue kit  as described by the manufacturer.
pcr and amplification of vntr loci
based on results from the tandem repeats analysis, oligonucleotide primers flanking  <dig> potential vntr loci was designed using the vector nti suite  <dig>  program package .
in addition, two pairs of flanking primers, vntr1-2f/r and vntr2-2f/r were designed for the purpose of nested pcr, as not all isolates were amplified using the initial primers.
the oligonucleotide primers had tm ranging from  <dig>  to  <dig> °c, and the annealing-temperature for each primer combination was therefore optimized using gradient pcr with the isolate gm <dig> dna as template .
the amplification was performed in a  <dig> μl reaction mixture containing  <dig> x buffer  including  <dig>  mmol- <dig> mgcl <dig>   <dig>  mm dntp ,  <dig> μmol- <dig> of each primer ,  <dig> μl dna and  <dig>  u thermal ace™dna polymerase .
amplification was carried out in a mastercycler gradient  with denaturation at 95°c for  <dig> min;  <dig> cycles with denaturation at 95°c for  <dig> s, annealing at x°c  for  <dig> s and extension at 72°c for  <dig> min, followed by prolonged extension at 72°c for  <dig> min and a short storage at 4°c.
all pcr-products were visualized using gel-electrophoresis.
in cases where the initial pcr did not amplify a target sequence, several techniques were applied including nested- , touchdown- and gradient-pcr.
in the touchdown pcr, the amplification cycles were as follows; denaturation at 95°c for  <dig> min; followed by  <dig> times of  <dig> cycles of; denaturation at 94°c for  <dig> s, annealing for  <dig> s at 66-51°c with a reduction of 3°c at each cycle, and extension at 72°c for  <dig> min.
the rounds of touchdown were immediately followed by  <dig> cycles with denaturation at 94°c for  <dig> s, annealing at 48°c for  <dig> s, extension at 72°c for  <dig> min; and a final extension at 72°c for  <dig> min before a short storage at 4°c.
the gradient pcr amplification was as follows: denaturation at 95°c for  <dig> min; denaturation at 94°c for  <dig> s,  <dig> and 4°c below and over optimized annealing temperature  for  <dig> s, extension at 72°c for  <dig> min; and a final extension at 72°c for  <dig> min before a short storage at 4°c.
flanking primers were designed based on preliminary data of the whole genome sequence of f. noatunensis ssp.
noatunensis isolate gm <dig>  the attributes of the primers in the study; primer sequence, optimized annealing temperature  and amplicon size of the gm <dig> isolate.
sequencing of vntr loci
in order to identify repeat copy number variation among the isolates in question the resulting pcr-products were sequenced for each vntr locus.
pcr products were purified using the e.z.n.a cycle-pure kit  as described by the manufacturer.
sequencing was then performed in both directions using the pcr-primers  and the abi prism bigdye terminator chemistry  according to applied biosystems .
all sequences were assembled using the vector nti suite  <dig>  program .
possible genetic location of sequenced vntr-loci was identified by genebank blastn searches performed using the complete amplicon sequence from each primer pair from either f. philomiragia or f. noatunensis ssp.
noatunensis gm <dig> isolate.
phylogenetic analysis
at each vntr locus in a single taxon, the vntr was coded as a discrete character  based upon the specific number of repeats at the region in question.
these allele profiles were used to construct a data matrix within the mesquite system for phylogenetic analysis .
the data-matrix was exported as a nexus file into paup  <dig>  .
version  <dig>  sinauer associates, sunderland, massachusetts)) for phylogenetic analysis using the neighbor-joining  distance method  <cit> .
the phylogenetic nj analysis was not bootstrapped as there were too few characters and taxa for correct estimation of sampling error.
test of vntr stability
since vntr regions are considered to be one of the fastest evolving sequences of a genome, it was important to test the stability of the selected vntr regions both in vitro and in vivo  <cit> .
mutations in vntr regions have been shown during in vitro growth of both yersinia pestis and escherichia coli  <cit> .
therefore the effect of multiple passages of the f. noatunensis ssp.
noatunensis gm <dig> isolate on vntr stability was tested as follows: a frozen first passage culture of the gm <dig> isolate was streaked for isolation on chab, one single colony was transferred in each passage for  <dig> passages.
dna isolation, pcr and sequencing were performed on a colony from passage  <dig> as described above.
pcr and sequencing were performed for all vntr loci for verification of copy numbers .
the in vivo stability of f. noatunensis ssp.
noatunensis gm <dig> isolate was tested at different temperatures with material from challenge experiments described below.
copy numbers at each vntr locus was confirmed through direct sequencing using dna extracted from kidney tissues of challenged fish as templates .
kidney tissue from one cod in each group challenged with gm <dig>  and held at  <dig> ,  <dig>  and 18°c , was sampled for extraction of bacterial dna for direct sequencing.
tank conditions for all groups were as for the challenge experiment of the f <dig> group described in nylund et al.
<cit> .
the f <dig> and f <dig> group were bath challenged in  <dig> l of sea water containing a final concentration of  <dig>  ×  <dig> bacteria/ml.
the bacterial suspension was prepared by inoculating several chab plates with the gm <dig> isolate, harvesting + mixing of bacterial cultures in pbs and subsequently adding this to sea water.
tissues were sampled from the fish  <dig> days post challenge.
the challenge experiments were approved by the norwegian animal research authorities .
in vivo stability was shown on passage  <dig>  of the gm <dig> isolate and in vitro stability on infected cod held at  <dig> ,  <dig>  and 18°c .
results are displayed with number of nucleotides in repeat marker and numbers of repeats.
two shifts were observed in gm <dig>  vntr- <dig> after  <dig> passages on chab agar at 20°c and in vntr- <dig> after passage in cod held at 18°c.
the original sequence is shown in bold, whereas the mutations are in italics.
heterogeneous repeats with single nucleotide polymorphism , noted as nucleotide number in repeat sequence  and at repeat of occurrence ; e =  <dig> snp c-t nt <dig> r <dig> r6
vntr linkage disequilibrium
linkage disequilibrium was calculated as standardized index of association , for all norwegian f. noatunensis ssp.
noatunensis isolates  within the seven vntr loci using the lian linkage analysis  <dig>  online tool  <cit> .
accession numbers
sequences of all seven vntr loci from the gm <dig> isolate of f. noatunensis ssp.
noatunensis were assigned genebank accession numbers as follows: vntr-1: gu <dig>  vntr-2: gu <dig>  vntr-3: gu <dig>  vntr-4: gu <dig>  vntr-5: gu <dig>  vntr-6: gu <dig>  vntr-7: gu385773
RESULTS
variable number of tandem repeats  features
sequencing of francisella noatunensis and f. philomiragia isolates  showed variation in seven of the  <dig> tested vntr loci .
the remaining four vntr loci were discarded as they did not get amplified from all f. noatunesis isolates or showed no allelic variation among the isolates included in this study.
the seven vntr loci used in this study had an allelic diversity ranging from two  to  <dig>  alleles in all isolates included .
the largest observed variation in allele sizes was found in vntr- <dig> and spanned from one to  <dig> repeats in f. philomiragia isolates fp080107-ii and fp <dig>  respectively.
primers for the seven informative vntr loci  provided pcr-products for all f. noatunensis ssp.
noatunensis and orientalis isolates with the exception of vntr- <dig> in isolate ua <dig> and vntr- <dig> in ehime- <dig>  vntr- <dig> and  <dig> showed most variation among the f. noatunensis ssp.
noatunensis isolates , with seven different alleles in each, while only minor variation, two different alleles, was found in vntr- <dig>  the f. noatunensis ssp.
orientalis isolates  had identical alleles at vntr- <dig> and - <dig>  while vntr- <dig> was shown to be the most polymorphic with four different alleles.
the repeats and attributes of each vntr locus in f. noatunensis can be viewed in table  <dig>  vntr- <dig>  - <dig> and - <dig> were successfully amplified from all isolates of f. philomiragia included in this study.
it was not possible to amplify vntr- <dig> from the f. philomiragia isolate, fp080107-i  <cit> , and, with the exception of the amplification of vntr- <dig> from dsm7535t, it was not possible to amplify vntr- <dig>  - <dig> and - <dig> from these isolates .
among the f. philomiragia isolates most allele variation was seen in vntr- <dig>  whereas no variations were detected in sequences from vntr- <dig> and - <dig>  it is not known if the lack of a pcr product for a given vntr locus is a result of primer mismatch or absence of the vntr locus in question.
due to the preliminary status of the genome from which the vntr loci were identified, the genomic locations of the loci remain unknown.
possible locations of the vntr loci within the genome of f. philomiragia were identified by genebank blastn searches.
the following matches were identified ; vntr- <dig> locus  matched the intergenic segment located between the genes encoding the bor lipoprotein and nicotinamide adenine dinucleotide phosphate -quinone reductase , vntr- <dig> sequence  matched the intragenic segment in the dna-directed dna polymerase gene , vntr- <dig> sequence  matched the segment intragenically located in 50s ribosomal protein l <dig> gene  ,vntr- <dig> sequence  matched the intragenically located segment between and within the genes encoding the ftn_ <dig> hypothetical protein and trigger factor protein .
no significant matches were identified for vntr- <dig>  vntr- <dig> and - <dig> were not amplified from dsm  <dig>  but a blastn search was performed using these vntr sequences obtained from f. noatunensis ssp.
noatunensis gm <dig> isolate.
vntr- <dig> displayed intragenic match for a gene encoding the ftn_ <dig> hypothetical protein  with 99% sequence-identity and e value of  <dig> × 10- <dig>  while vntr- <dig> showed no significant match to the f. philomiragia genome.
identical allelic profiles are assigned greek numerals i-iii.
allelic profiles are presented as number of nucleotides in repeat marker  and numbers of repeats .
identical profiles occur for several isolates from cod culturing, profile i , profile ii  and profile iii
isolated in chilea, isolated from broodfishf, type straint, w = wild caught, na = not amplified with primer sets
heterogeneous repeats with single nucleotide polymorphism , noted as nucleotide number in repeat sequence  and at repeat of occurrence ; a = snp g-a nt <dig> r <dig> r <dig> b = snp.
g-a nt <dig> r <dig>  c = snp g-a nt <dig> r <dig>  d = snp t-c nt <dig> from r <dig>  e =  <dig> snp c-t nt <dig> r <dig> r <dig>  f =  <dig> snp c-t nt <dig> r <dig> r <dig>  g =  <dig> snp c-t nt <dig> r <dig> r <dig>  h =  <dig> snp c-t nt <dig> r <dig> r <dig>  i = snp t-c nt <dig> r1
the blastn search was performed using the complete amplicon of the primer pairs.
all locations were identified in f. philomiragia atcc25017t .
sequences were derived from f. philomiragia dsm7535t unless otherwise stated.
vntr-sequence obtained from f. noatunensis ssp.
noatunensis gm2212a
allele profiles
variation in the vntr loci was revealed by sequencing, and these data were used to create allele profiles for all isolates.
lack of amplifiable vntr locus was interpreted as a separate "missing" character.
a complete overview of the allele profiles is presented in table  <dig>  analysis of the seven vntr loci revealed a total of  <dig> unique allele profiles among the  <dig> f. noatunensis ssp.
noatunensis isolates.
three of the profiles contained more than one isolate and these were termed profile i , profile ii  and profile iii ; all isolates included in these profiles were collected from farmed sites with the exception of one wild isolate in profile ii.
profile i included the following isolates from farmed cod: a) four isolates collected in the period of 2004- <dig> from the rogaland county, b) one isolate from sogn and fjordane county  and c) three isolates from møre and romsdal county .
the two isolates in profile ii were collected from rogaland  and nordland  in  <dig>  a total of five isolates were included in profile iii: a) three isolates, including two originating from broodfish populations, collected in hordaland county during the period 2005- <dig>  and b) two isolates, one sampled from farmed and one from a wild/possible escaped farmed cod in sogn og fjordane county in  <dig>  the remaining six cod isolates, including three isolates from wild cod, two from broodstock populations  and one from a production site, all possessed unique allele profiles.
the two isolates from wild cod in aust- and vest agder counties share allele profiles at five out of seven loci.
the allele profile of the chilean isolate from atlantic salmon, ua <dig>  had a unique profile sharing only one vntr locus with the other isolates of f. noatinensis ssp.
noatunensis.
all seven vntrs were obtained from the isolates of f. noatunensis ssp.
orientalis  with the exception of vntr- <dig> from ehime- <dig>  none of the isolates had identical allele profiles and the allele profile of the japanese isolate ehime- <dig>  from farmed three-line grunt , was distinct from the three indonesian isolates from farmed tilapia .
none of the f. philomiragia isolates  had identical allelic profiles.
phylogenetic relationship
the nj phylogenetic analysis using the seven informative vntr loci from the  <dig> francisella spp.
isolates showed a subdivision of the f. noatunensis ssp.
noatunensis isolates .
the majority of the isolates from outbreaks of francisellosis in farmed cod group as two distinct clades.
these clades were named clade i and iii, as they correlate to the allele profiles .
the phylogenetic analysis also separates the chilean ua <dig> isolate from the norwegian f. noatunensis ssp.
noatunensis isolates, and differentiated among the few isolates of f. noatunensis ssp.
orientalis, included.
stability of vntr loci
examination of the stability of the vntr loci revealed changes in repeat numbers within two vntrs of the f. noatunensis ssp.
noatunensis gm <dig> isolate .
a shift in vntr- <dig> was observed after  <dig> passages on cysteine heart agar with blood  at 20°c.
this mutation consisted of the addition of one repeat, increasing the size from  <dig> to  <dig> repeats.
there was also a change in vntr- <dig> after passage of isolate gm <dig> in cod held at 18°c.
the number of repeats was reduced from  <dig> to nine.
the vntrs did not change when the isolate gm <dig> was passed through cod held at  <dig> and 14°c.
these results were verified by a repeating round of pcr and sequencing.
the remaining vntrs were not affected in this test.
linkage disequilibrium
to avoid effects from host and geographical separation the linkages disequilibrium was calculated for f. noatunensis ssp.
noatunensis isolates  only.
the linkages disequilibrium was shown to be significant in the lian  <dig>  analysis, indicating a clonal population structure for the cod isolates.
standardized ias was calculated to  <dig>  at a significance of ppara =  <dig>  × 10- <dig>
epizootiological data of isolates
there are several factors linking the isolates within the f. noatunensis clades .
clade i consists of eight isolates where fnnr-001- <dig>  -002- <dig>  004- <dig> and 017- <dig> were obtained from cod at one production site in rogaland county during several outbreaks in the period of 2004- <dig> .
these isolates differ from an isolate  obtained from one of the broodfish company supplying the site.
on-growth sites were not involved.
three isolates fnnmr-011- <dig>  -013- <dig> and -005- <dig> of clade i were from another county , and two of these were obtained from cod with no traceable history.
the third isolate fnnmr-005- <dig> was obtained from a site supplied by the same broodfish company as the host for isolate fnnsf-012- <dig>  which also belongs to clade i. the latter isolate comes from another county, sogn og fjordane.
the host for fnnmr-005- <dig> was kept at an on-growing site before transportation to the production site in møre og romsdal.
clade iii consists of five isolates, three from hordaland  and two from sogn of fjordane county  .
the latter two were sampled in the same fjord, fnnsf018- <dig> was isolated from farmed cod at a production site, while fnn-019- <dig> was isolated from a wild caught cod.
the wild caught cod had pellets in the gut and a morphology suggesting that it originated from a production site which lost fish during a francisellosis outbreak in  <dig>  of the last three isolates from clade iii, one  was isolated from cod at a production site in hordaland, while fnnh-007- <dig> and fnnh-014- <dig> was isolated from broodfish populations located at two sites in the same county.
this broodfish company was one of the suppliers of cod to the production site in sogn and fjordane  and to one production site in hordaland .
it is not known if the cod at the two sites, sogn and fjordane and hordaland , were offspring from the broodfish populations where fnnh-007- <dig> and fnnh-014- <dig> were isolated.
discussion
variable number of tandem repeats  system
whole genome sequencing of bacteria has presented new opportunities for identification of new genetic markers for separation of isolates.
one such marker system can be found by looking at vntrs, i.e., single locus sequences with short dna repeats  <cit> .
micro satellites, a subset of vntrs, with repeat motifs of nine bp or less are often targeted due to higher mutational rate  <cit> .
however, such hyper variability, that may occur within some vntrs, would complicate determination of genetic relationships among strains using this method, and hence, its use in phylogenetics may not be ideal  <cit> .
the mechanisms behind the length variation in micro satellites is that of slipped-strand mispairing  during dna polymerase mediated dna duplication.
however, mutations involving indels of large copynumbers have been shown consistent with recombination-mediated events  <cit> .
in escherichia coli, an average mutational rate of  <dig>  × 10- <dig> was calculated over  <dig> vntrs, and the rate seems to be dependent on intrinsic properties such as numbers of repeats  <cit> .
variation in copy numbers at vntr loci of certain sizes may affect the efficiency of promoters, thus affecting the coding potential of genes dependent on the genomic locations and the indels of repeats  <cit> .
for differentiation of bacterial isolates several vntr loci are combined in a multiple locus vntr analysis , a well established tool for epidemiological studies of bacteria with highly fit clonal dominance  <cit> .
one such bacteria is francisella tularensis, the agent for tularaemia  <cit> .
this species consist of clonal geographical subspecies  <cit>  that can be separated with the help of vntr markers  <cit> .
a similar approach could also be the solution for genotyping of francisella isolates belonging to fish pathogenic francisella spp.
most mlva systems applied in epidemiology are performed by multiplexing pcr-assays where vntr size is deduced by capillary electrophoresis allowing rapid screening of large datasets  <cit> .
due to the relatively few vntr-loci included and that there were only  <dig> isolates available for analysis in this study, sequencing was selected as a method for repeat number verification.
the sequencing also detected single nucleotide polymorphisms in addition to allele size.
the future potential for use of multiplex pcr assays with dyed primers and determining allele size by capillary electrophoresis is still an option and is recommended for larger datasets.
this study presents a sequence based mlva system consisting of seven vntr loci identified in the f. noatunensis ssp.
noatunensis gm <dig> isolate.
when applied to  <dig> francisella strains including norwegian and asian isolates of f. noatunensis, the mlva provided the best resolution shown for the fish-pathogenic francisellae so far.
due to the preliminary genome status of the f. noatunensis ssp.
noatunensis gm <dig> isolate a complete vntr locus search has not been performed, thus leaving the possibility for unidentified markers yet to be discovered.
nonetheless strain typing of f. tularensis isolates has been achieved with the use of as few as two and six vntr markers  <cit> .
the exact location of the seven vntr markers in the f. noatunensis ssp.
noatunensis is not known, however, blast searches indicate both intra- and intergenic location if the positions are homologous with that of the f. philomiragia atcc25017t isolate .
the discriminatory power of the vntr loci used in this study, was not calculated as a simpson`s index of diversity.
the reason for this is that the dataset does not fulfill the criteria of a test population proposed by van belkum et al.
<cit> , as isolates are predominately from outbreaks in cultured populations of cod with potential epizootic connections.
the vntr markers had a typeability of 100% among the norwegian isolates , and with the exception of vntr- <dig> and - <dig> for the ua <dig> and ehime isolate respectively, all the vntr markers for f. noatunensis  were amplified, indicating the usability of the typing system.
stability testing showed variability within vntr- <dig> and - <dig> from f. noatunensis ssp.
noatunensis gm <dig> isolate when grown in vitro for  <dig> passages at 20°c and in vivo in cod at 18°c.
such single repeat changes are shown to occur in 80% of mutational events with equal chance of being an insertion or a deletion  <cit> .
the sampled size was insufficient to assess the potential of hyper variable sites and mutation rates in the f. noatunensis ssp.
noatunensis gm <dig> isolate.
vntr mutational rates can differ across loci  <cit> , and could be a problem in studies of phylogenetics and epizootics.
stability and introduction of variation could be a result of bacterial phase variation and thus a result of environmental factors such as nutritional substances present in the growth medium, temperature effects  <cit> , or just simply be stochastical events.
however, four isolates  from different francisellosis outbreaks at the same site during a two year period showed identical allelic profiles indicating field stability and epizootic coherence.
the nj dendrogram presented in this study should not be viewed as a phylogenetic analysis as the use of vntrs and especially highly variable vntrs create noise in the phylogenetic signal.
the stability testing in the study shows that this may be the case of the vntr- <dig> and vntr- <dig>  there are also too few isolates and characters  included in the dataset to attempt a proper analysis of genetic relationships.
f. noatunensis strains
a total of nine allele profiles were identified among the norwegian f. noatunensis ssp.
noatunensis isolates .
three of these contained more then one isolate  and were identified using allelic profiles and visualized in the nj dendrogram.
the remaining six isolates represent unique strains of which five originated from wild cod and one from farmed cod .
although the dataset of norwegian isolates is small the trend is evident with only a few isolates occurring in outbreaks of francisellosis in cod culture, and a high diversity of strains among the isolates from wild cod.
the chilean isolate  and the norwegian cod isolates have previously been reported to be highly similar based on similarities in 16s sequence  and in five housekeeping genes   <cit> .
however, according to the vntr analysis performed in this study, the ua <dig> isolate from salmon differs at six of seven markers compared to the norwegian f. noatunensis ssp.
noatunensis isolates from cod, clearly separating them into two distinct branches in the nj dendrogram.
the clear separation of these isolates most likely reflects geographic origin and the ecological differences between them .
the allele profiles of the four f. noatunensis ssp.
orientalis isolates were unique, separating them into four different genotypes.
although the sampling set is small it correctly reflects the geographical origin and, possibly, also the variation in niche.
there was no clear pattern in the clustering of the f. philomiragia isolates.
care is needed if phylogenetic relationships are to be estimated based on vntrs  <cit> , however the deeper nodes in the dendrogram  a) between the species of f. philomiragia and f. noatunensis, b) between the subspecies f. noatunensis ssp.
orientalis and f. noatunensis ssp.
noatunensis and c) among the f. noatunensis ssp.
noatunensis isolates of norway  and chile , were similar to the tree topology obtained based on analysis of 16s rrna gene and datasets from housekeeping genes  <cit> .
epidemiology of the f. noatunensis
considering the results of the present study some reflections regarding the epizootics of the f. noatunensis are possible.
however, it is necessary to acknowledge the shortcomings of the dataset as it covers about 40% of the official diagnosed outbreaks of francisellosis in norway and that it only includes  <dig> wild isolates collected in different counties.
generally, the low allelic diversity among the isolates collected from cultured cod compared to the isolates from wild cod, indicate that the spreading could be a result of human activity.
this is also supported by the fact that the production cycle of cod in norway involves transport of farmed cod along most of the norwegian coast.
similarly, human activity has also been suggested as a cause for the geographical distribution of clonal isolates of bacillus anthracis  <cit> .
although the vntr markers provide a good starting point for studies of epizootics it may still be difficult to identify sources of infection as fish from different populations are mixed during the production, and, it is very difficult to avoid interactions between cultured and wild fish species.
in the case of clade i with the isolates fnnr-001- <dig>  -002- <dig>  004- <dig> and 017- <dig>  the new generations of cod used for restocking after initial outbreak were probably infected by the initial isolate in the population at that site.
the source of the initial infection was most likely cultured cod, since isolates belonging to the same genotype have been found at production sites in other counties, while the wild type isolate  from the same county belonged to another clade.
in clade iii, the isolates from a known broodfish supplier were identical to isolates from production sites where the aforementioned broodfish company was one of the suppliers of fish.
hence, all of these cases could be a result of transmission through movement of offspring from broodfish.
however, the situation is complicated by the fact that the sites used for on-growth are unknown.
vertical transmission would greatly impact the disease management and needs to be addressed in further studies.
CONCLUSIONS
we present the first vntr analysis of fish pathogenic francisellae.
seven polymorphic microsatellites were identified in a partial genome sequence of the f. noatunensis ssp.
noatunensis gm <dig> isolate.
a sequence based mlva system of these seven vntr-loci was applied to  <dig> aquatic francisella isolates, including both f. philomiragia  and f. noatunensis isolates .
all vntrs were amplified in the f. noatunensis isolates, with the exception of vntr- <dig> in ua <dig> and vntr- <dig> in ehime- <dig>  among the norwegian f. noatunensis ssp.
noatunensis isolates , including isolates from both farmed  and wild cod, a total of nine allelic profiles were identified.
the majority of farmed isolates were divided into two allelic profiles, indicating low allelic variation in isolates from outbreaks in cod culture compared to isolates collected from wild cod.
the allelic profile of the chilean f. noatunensis ssp.
noatunensis isolate reflected the geographical and host divergence when compared to norwegian cod isolates.
all f. noatunensis ssp.
orientalis isolates  show a unique allelic profile.
the four vntrs amplified from f. philomiragia provided a unique allelic profile for all these isolates.
the results show that this mlva system should provide a good starting point for future studies of epizootics of f. noatunensis.
authors' contributions
Øjb was responsible for carrying out this study and was the main contributor in writing the manuscript.
kfo identified the vntr loci in the draft genome, contributed to project design and writing.
an coordinated the project, contributed to project design and reviewed all drafts of the manuscript all authors read, commented on and approved the final manuscript.

 the intensive culturing of fish in artificially high population densities facilitates disease outbreaks. standard protocols for traditional fish health management, controlling and preventing diseases, have focused on vaccines, operational prophylactic measures and oral treatment using therapeutic agents. however, in those cases where vaccines or treatments are not available, understanding the epizootiology becomes the key to prevent outbreaks and pathogen dispersal. using molecular tools for strain typing of pathogens, combined with biological and ecological knowledge of both the pathogen and the host, is a prerequisite in creating an epizootiological understanding which can be applied in management of diseases in both wild and cultured populations. historically, the atlantic cod has been an important marine resource in norway, occurring naturally both as migratory and stationary populations. these populations are divided into coastal and oceanic populations, based on phenotypic and genotypic traits. the norwegian coastal cod population is considered to consist of several stationary sub populations dispersed along the norwegian coast. the oceanic cod in norwegian waters can be divided into two populations, the north east arctic cod and the north sea cod, both with seasonally dependent migratory behavior. the oceanic and some of the coastal populations will, at certain times of the year, be present in at the same locations. during the last decade cod has been a species of increasing significance and interest for the norwegian aquaculture industry. cod is now intensively cultured with full control of all life-stages, except for the broodfish which are predominantly of wild origin. the intensive production cycle of cod consists of three separate steps, where the first step starts with the fertilization of eggs from broodstock held in large land based tanks. the majority of broodfish consists of wild caught cod from both costal and oceanic populations. after hatching the fry is held in indoor tanks. the second stage is on-growth, where the fingerlings are transferred to net pens in the sea and held until they are moved to production sites. mixing of populations from different broodstock companies occurs at the on-growth sites, i.e. before redistribution to production sites. mixing of populations may also occur at the production sites due to limited availability of fish from the different on-growing sites. several generations can be present at one production site. the different operational sites are dispersed along the norwegian coast resulting in large scale movement of cultured cod between different parts of norway. in a francisella species was isolated from farmed cod in norway showing clinical signs of a chronic granulomatous infection in kidney, spleen, liver and heart. the bacterium was initially characterized both as a species, f. piscicida, and later in the same year as a f. philomiragia subspecies, ssp. noatunensis. the f. philomiragia ssp. noatunensis was later elevated to the rank of species, f. noatunensis ssp. noatunensis, with f. piscicida as a heterotypic synonym. following the discovery, several annual outbreaks of francisellosis have been diagnosed with main foci in the western parts of norway. the bacterium has been detected in cultured fish from rogaland to nordland county, however, positive wild cod has only been identified south of sogn og fjordane county, indicating a southern natural reservoir of the bacteria. after the discovery of f. noatunensis ssp. noatunensis, histological material collected in from a broodfish population in hordaland county, has been stained positive with antiserum for the f. noatunensis ssp. noatunensis gm isolate, suggesting the presence of the bacteria in norwegian waters prior to the initial discovery. during the last years members of the genus francisella have emerged as a global problem for aquaculture, causing mortality among a wide range of aquatic hosts. these francisella isolates show high genetic similarities with f. philomiragia and with each other at the 16s rrna-gene. the close genetic relatedness has been confirmed with sequencing of several housekeeping genes showing identical sequences for norwegian isolates. a chilean isolate from atlantic salmon ua showed a high similarity to the norwegian cod isolates and was therefore described as a new isolate of the f. noatunensis ssp. noatunensis. fish pathogenic francisella isolates from asia display high genetic similarity to ua and ncimb 14265t/gm and constitute a separate subspecies, f. noatunensis ssp. orientalis. f. asiatica isolated from tilapia in costa rica is identical with f. noatunensis ssp. orientalis with respect to phenotype and signature sequences. isolates of f. noatunensis ssp. noatunensis from norway are identical when comparing 16s and housekeeping gene sequences, making genotyping and studies of epizootics difficult. if an epizootiological approach is to be applied for controlling and preventing dispersal of francisella spp. in aquaculture, a tool for the identification of isolates is a necessity. such a tool will make it possible to differentiate between wild endemic and anthropogenic dispersed strains. this approach also requires knowledge of the production history of cod and some knowledge of natural occurring strains in wild cod in the production areas. targeting the housekeeping genes of bacterial fish pathogens has been applied for strain identification, but due to identical signature sequences this has not been possible to use on f. noatunensis. similar problems are known from epidemiological studies of human pathogenic bacteria like the mycobacterium tuberculosis, m. avium ssp. paratuberculosis, bacillus anthracis, yersinia pestis and the related f. tularensis where evolutionary or stochastic events have led to a dominance of highly fit clones making strain identification difficult. by using genetic markers with higher mutational rates, like variable number of tandem repeat, one can construct allelic profiles in multiple locus vntr analysis systems, making strain differentiation of clonal bacteria populations possible. mlva systems have already been applied in studies of f. tularensis isolates. this study is an attempt to provide a sequenced based mlva tool for studying the epizootics of f. noatunensis ssp. noatunensis in norwegian cod farming. our findings suggest that the mlva system presented in this study is suitable for strain typing of f. noatunensis isolates from norwegian cod, and the results indicate that there are only a few clades causing francisellosis outbreaks in norwegian cod culturing. methods isolation of francisella noatunensis ssp. noatunensis all isolates of francisella noatunensis ssp. noatunensis included in this study were obtained from atlantic cod suffering from francisellosis, with the exception of the type strain ncimb 14265t and the chilean atlantic salmon isolate ua which were obtained from the national collection of industrial, food and marine bacteria and intervet norbio as, respectively. the majority of isolates were collected from cultured populations of cod geographically dispersed along km of norwegian coast line during the period from to four isolates were obtained from wild-cod, from the counties of aust-agder, vest-agder, rogaland and from the borders of rogaland and hordaland. an overview of geographical origin of isolates is presented in figure isolates of closely related francisella spp., f. notunensis ssp. orientalis and f. philomiragia, were also included in this study. isolate foj-001/ was kindly provided by dr. kamaishi, and isolates foi-002/ foi-003/ and foi-004/ from indonesia were kindly provided by intervet singapore. most f. philomiragia isolates were obtained from the culture collections deutsche sammlung von mikroorganismen und zellkulturen and culture collection, university of göteborg, while strains and strain were kindly provided by dr. berrada and dr. friis-møller, respectively. a complete overview of isolates and year of isolation, host and location is presented in table the table show source and geographical location of isolation for the francisella isolates included in the current study. norway is country of origin unless otherwise is stated. isolates were kindly provided by intervet norbio asa, dr. kamaishib, intervet singaporec, dr. friis-møllerd and dr. berradae isolated from broodfishf, isolated in chilea, w = wild caught, t = type strain culture of f. noatunensis ssp. noatunensis isolations of f. noatunensis ssp. noatunensis were performed in the field or in our laboratory by streaking tissue from infected cod on agar plates. initial isolations were performed using blood agar plates with % cysteine and 1% glucose as described by nylund et al. with subsequent growth in b broth for cryopreservation. however, most isolations were performed using cysteine heart agar supplemented with 5% bovine blood as described in olsen et al. with slight modifications. these included the addition of μg ml- of ampicillin, 50- μg ml- of fungizone® and additional cysteine to a final concentration of %. the ampicillin, fungizone®, bovine blood and cysteine were added to the chab at 60°c. all isolates were incubated at 20°c for 2- days or two weeks. a sanyo mir- incubator was used for solid cultures whereas liquid culturing was performed in an infors unitron incubator at rpm. colonies from chab agar were dissolved in eugon broth for cryopreservation and stored in liquid nitrogen. genomic analysis a preliminary genomic sequence from francisella noatunensis ssp. noatunensis isolate gm was analyzed for the presence of tandem repeat regions using the software tandem repeats finder. this program locates and displays tandem repeats in dna sequences. loci with tandem repeats consisting of less than nucleotides and more than four repeats were selected. the preliminary genomic sequence of f. noatunensis ssp. noatunensis isolate gm was generated using the pyro-sequencing technology in co-operation with intervet as/asa. genbank accession numbers refer to sequences of the f. noatunensis ssp. noatunensis gm isolate. repeat sequence, number of repeats and amplicon size is that of the gm isolate. not amplified within f. philomiragia isolatesna dna extraction total dna from all collected strains and tissues was extracted using the dneasy® blood & tissue kit as described by the manufacturer. pcr and amplification of vntr loci based on results from the tandem repeats analysis, oligonucleotide primers flanking potential vntr loci was designed using the vector nti suite program package. in addition, two pairs of flanking primers, vntr1-2f/r and vntr2-2f/r were designed for the purpose of nested pcr, as not all isolates were amplified using the initial primers. the oligonucleotide primers had tm ranging from to °c, and the annealing-temperature for each primer combination was therefore optimized using gradient pcr with the isolate gm dna as template. the amplification was performed in a μl reaction mixture containing x buffer including mmol- mgcl mm dntp, μmol- of each primer, μl dna and u thermal ace™dna polymerase. amplification was carried out in a mastercycler gradient with denaturation at 95°c for min; cycles with denaturation at 95°c for s, annealing at x°c for s and extension at 72°c for min, followed by prolonged extension at 72°c for min and a short storage at 4°c. all pcr-products were visualized using gel-electrophoresis. in cases where the initial pcr did not amplify a target sequence, several techniques were applied including nested-, touchdown- and gradient-pcr. in the touchdown pcr, the amplification cycles were as follows; denaturation at 95°c for min; followed by times of cycles of; denaturation at 94°c for s, annealing for s at 66-51°c with a reduction of 3°c at each cycle, and extension at 72°c for min. the rounds of touchdown were immediately followed by cycles with denaturation at 94°c for s, annealing at 48°c for s, extension at 72°c for min; and a final extension at 72°c for min before a short storage at 4°c. the gradient pcr amplification was as follows: denaturation at 95°c for min; denaturation at 94°c for s, and 4°c below and over optimized annealing temperature for s, extension at 72°c for min; and a final extension at 72°c for min before a short storage at 4°c. flanking primers were designed based on preliminary data of the whole genome sequence of f. noatunensis ssp. noatunensis isolate gm the attributes of the primers in the study; primer sequence, optimized annealing temperature and amplicon size of the gm isolate. sequencing of vntr loci in order to identify repeat copy number variation among the isolates in question the resulting pcr-products were sequenced for each vntr locus. pcr products were purified using the e.z.n.a cycle-pure kit as described by the manufacturer. sequencing was then performed in both directions using the pcr-primers and the abi prism bigdye terminator chemistry according to applied biosystems. all sequences were assembled using the vector nti suite program. possible genetic location of sequenced vntr-loci was identified by genebank blastn searches performed using the complete amplicon sequence from each primer pair from either f. philomiragia or f. noatunensis ssp. noatunensis gm isolate. phylogenetic analysis at each vntr locus in a single taxon, the vntr was coded as a discrete character based upon the specific number of repeats at the region in question. these allele profiles were used to construct a data matrix within the mesquite system for phylogenetic analysis. the data-matrix was exported as a nexus file into paup. version sinauer associates, sunderland, massachusetts)) for phylogenetic analysis using the neighbor-joining distance method. the phylogenetic nj analysis was not bootstrapped as there were too few characters and taxa for correct estimation of sampling error. test of vntr stability since vntr regions are considered to be one of the fastest evolving sequences of a genome, it was important to test the stability of the selected vntr regions both in vitro and in vivo. mutations in vntr regions have been shown during in vitro growth of both yersinia pestis and escherichia coli. therefore the effect of multiple passages of the f. noatunensis ssp. noatunensis gm isolate on vntr stability was tested as follows: a frozen first passage culture of the gm isolate was streaked for isolation on chab, one single colony was transferred in each passage for passages. dna isolation, pcr and sequencing were performed on a colony from passage as described above. pcr and sequencing were performed for all vntr loci for verification of copy numbers. the in vivo stability of f. noatunensis ssp. noatunensis gm isolate was tested at different temperatures with material from challenge experiments described below. copy numbers at each vntr locus was confirmed through direct sequencing using dna extracted from kidney tissues of challenged fish as templates. kidney tissue from one cod in each group challenged with gm and held at, and 18°c, was sampled for extraction of bacterial dna for direct sequencing. tank conditions for all groups were as for the challenge experiment of the f group described in nylund et al.. the f and f group were bath challenged in l of sea water containing a final concentration of × bacteria/ml. the bacterial suspension was prepared by inoculating several chab plates with the gm isolate, harvesting + mixing of bacterial cultures in pbs and subsequently adding this to sea water. tissues were sampled from the fish days post challenge. the challenge experiments were approved by the norwegian animal research authorities. in vivo stability was shown on passage of the gm isolate and in vitro stability on infected cod held at, and 18°c. results are displayed with number of nucleotides in repeat marker and numbers of repeats. two shifts were observed in gm vntr- after passages on chab agar at 20°c and in vntr- after passage in cod held at 18°c. the original sequence is shown in bold, whereas the mutations are in italics. heterogeneous repeats with single nucleotide polymorphism, noted as nucleotide number in repeat sequence and at repeat of occurrence; e = snp c-t nt r r6 vntr linkage disequilibrium linkage disequilibrium was calculated as standardized index of association, for all norwegian f. noatunensis ssp. noatunensis isolates within the seven vntr loci using the lian linkage analysis online tool. accession numbers sequences of all seven vntr loci from the gm isolate of f. noatunensis ssp. noatunensis were assigned genebank accession numbers as follows: vntr-1: gu vntr-2: gu vntr-3: gu vntr-4: gu vntr-5: gu vntr-6: gu vntr-7: gu385773  variable number of tandem repeats features sequencing of francisella noatunensis and f. philomiragia isolates showed variation in seven of the tested vntr loci. the remaining four vntr loci were discarded as they did not get amplified from all f. noatunesis isolates or showed no allelic variation among the isolates included in this study. the seven vntr loci used in this study had an allelic diversity ranging from two to alleles in all isolates included. the largest observed variation in allele sizes was found in vntr- and spanned from one to repeats in f. philomiragia isolates fp080107-ii and fp respectively. primers for the seven informative vntr loci provided pcr-products for all f. noatunensis ssp. noatunensis and orientalis isolates with the exception of vntr- in isolate ua and vntr- in ehime- vntr- and showed most variation among the f. noatunensis ssp. noatunensis isolates, with seven different alleles in each, while only minor variation, two different alleles, was found in vntr- the f. noatunensis ssp. orientalis isolates had identical alleles at vntr- and - while vntr- was shown to be the most polymorphic with four different alleles. the repeats and attributes of each vntr locus in f. noatunensis can be viewed in table vntr- - and - were successfully amplified from all isolates of f. philomiragia included in this study. it was not possible to amplify vntr- from the f. philomiragia isolate, fp080107-i, and, with the exception of the amplification of vntr- from dsm7535t, it was not possible to amplify vntr- - and - from these isolates. among the f. philomiragia isolates most allele variation was seen in vntr- whereas no variations were detected in sequences from vntr- and - it is not known if the lack of a pcr product for a given vntr locus is a result of primer mismatch or absence of the vntr locus in question. due to the preliminary status of the genome from which the vntr loci were identified, the genomic locations of the loci remain unknown. possible locations of the vntr loci within the genome of f. philomiragia were identified by genebank blastn searches. the following matches were identified; vntr- locus matched the intergenic segment located between the genes encoding the bor lipoprotein and nicotinamide adenine dinucleotide phosphate -quinone reductase, vntr- sequence matched the intragenic segment in the dna-directed dna polymerase gene, vntr- sequence matched the segment intragenically located in 50s ribosomal protein l gene,vntr- sequence matched the intragenically located segment between and within the genes encoding the ftn_ hypothetical protein and trigger factor protein. no significant matches were identified for vntr- vntr- and - were not amplified from dsm but a blastn search was performed using these vntr sequences obtained from f. noatunensis ssp. noatunensis gm isolate. vntr- displayed intragenic match for a gene encoding the ftn_ hypothetical protein with 99% sequence-identity and e value of × 10- while vntr- showed no significant match to the f. philomiragia genome. identical allelic profiles are assigned greek numerals i-iii. allelic profiles are presented as number of nucleotides in repeat marker and numbers of repeats. identical profiles occur for several isolates from cod culturing, profile i, profile ii and profile iii isolated in chilea, isolated from broodfishf, type straint, w = wild caught, na = not amplified with primer sets heterogeneous repeats with single nucleotide polymorphism, noted as nucleotide number in repeat sequence and at repeat of occurrence; a = snp g-a nt r r b = snp. g-a nt r c = snp g-a nt r d = snp t-c nt from r e = snp c-t nt r r f = snp c-t nt r r g = snp c-t nt r r h = snp c-t nt r r i = snp t-c nt r1 the blastn search was performed using the complete amplicon of the primer pairs. all locations were identified in f. philomiragia atcc25017t. sequences were derived from f. philomiragia dsm7535t unless otherwise stated. vntr-sequence obtained from f. noatunensis ssp. noatunensis gm2212a allele profiles variation in the vntr loci was revealed by sequencing, and these data were used to create allele profiles for all isolates. lack of amplifiable vntr locus was interpreted as a separate "missing" character. a complete overview of the allele profiles is presented in table analysis of the seven vntr loci revealed a total of unique allele profiles among the f. noatunensis ssp. noatunensis isolates. three of the profiles contained more than one isolate and these were termed profile i, profile ii and profile iii; all isolates included in these profiles were collected from farmed sites with the exception of one wild isolate in profile ii. profile i included the following isolates from farmed cod: a) four isolates collected in the period of 2004- from the rogaland county, b) one isolate from sogn and fjordane county and c) three isolates from møre and romsdal county. the two isolates in profile ii were collected from rogaland and nordland in a total of five isolates were included in profile iii: a) three isolates, including two originating from broodfish populations, collected in hordaland county during the period 2005- and b) two isolates, one sampled from farmed and one from a wild/possible escaped farmed cod in sogn og fjordane county in the remaining six cod isolates, including three isolates from wild cod, two from broodstock populations and one from a production site, all possessed unique allele profiles. the two isolates from wild cod in aust- and vest agder counties share allele profiles at five out of seven loci. the allele profile of the chilean isolate from atlantic salmon, ua had a unique profile sharing only one vntr locus with the other isolates of f. noatinensis ssp. noatunensis. all seven vntrs were obtained from the isolates of f. noatunensis ssp. orientalis with the exception of vntr- from ehime- none of the isolates had identical allele profiles and the allele profile of the japanese isolate ehime- from farmed three-line grunt, was distinct from the three indonesian isolates from farmed tilapia. none of the f. philomiragia isolates had identical allelic profiles. phylogenetic relationship the nj phylogenetic analysis using the seven informative vntr loci from the francisella spp. isolates showed a subdivision of the f. noatunensis ssp. noatunensis isolates. the majority of the isolates from outbreaks of francisellosis in farmed cod group as two distinct clades. these clades were named clade i and iii, as they correlate to the allele profiles. the phylogenetic analysis also separates the chilean ua isolate from the norwegian f. noatunensis ssp. noatunensis isolates, and differentiated among the few isolates of f. noatunensis ssp. orientalis, included. stability of vntr loci examination of the stability of the vntr loci revealed changes in repeat numbers within two vntrs of the f. noatunensis ssp. noatunensis gm isolate. a shift in vntr- was observed after passages on cysteine heart agar with blood at 20°c. this mutation consisted of the addition of one repeat, increasing the size from to repeats. there was also a change in vntr- after passage of isolate gm in cod held at 18°c. the number of repeats was reduced from to nine. the vntrs did not change when the isolate gm was passed through cod held at and 14°c. these results were verified by a repeating round of pcr and sequencing. the remaining vntrs were not affected in this test. linkage disequilibrium to avoid effects from host and geographical separation the linkages disequilibrium was calculated for f. noatunensis ssp. noatunensis isolates only. the linkages disequilibrium was shown to be significant in the lian analysis, indicating a clonal population structure for the cod isolates. standardized ias was calculated to at a significance of ppara = × 10- epizootiological data of isolates there are several factors linking the isolates within the f. noatunensis clades. clade i consists of eight isolates where fnnr-001- -002- 004- and 017- were obtained from cod at one production site in rogaland county during several outbreaks in the period of 2004-. these isolates differ from an isolate obtained from one of the broodfish company supplying the site. on-growth sites were not involved. three isolates fnnmr-011- -013- and -005- of clade i were from another county, and two of these were obtained from cod with no traceable history. the third isolate fnnmr-005- was obtained from a site supplied by the same broodfish company as the host for isolate fnnsf-012- which also belongs to clade i. the latter isolate comes from another county, sogn og fjordane. the host for fnnmr-005- was kept at an on-growing site before transportation to the production site in møre og romsdal. clade iii consists of five isolates, three from hordaland and two from sogn of fjordane county. the latter two were sampled in the same fjord, fnnsf018- was isolated from farmed cod at a production site, while fnn-019- was isolated from a wild caught cod. the wild caught cod had pellets in the gut and a morphology suggesting that it originated from a production site which lost fish during a francisellosis outbreak in of the last three isolates from clade iii, one was isolated from cod at a production site in hordaland, while fnnh-007- and fnnh-014- was isolated from broodfish populations located at two sites in the same county. this broodfish company was one of the suppliers of cod to the production site in sogn and fjordane and to one production site in hordaland. it is not known if the cod at the two sites, sogn and fjordane and hordaland, were offspring from the broodfish populations where fnnh-007- and fnnh-014- were isolated. discussion variable number of tandem repeats system whole genome sequencing of bacteria has presented new opportunities for identification of new genetic markers for separation of isolates. one such marker system can be found by looking at vntrs, i.e., single locus sequences with short dna repeats. micro satellites, a subset of vntrs, with repeat motifs of nine bp or less are often targeted due to higher mutational rate. however, such hyper variability, that may occur within some vntrs, would complicate determination of genetic relationships among strains using this method, and hence, its use in phylogenetics may not be ideal. the mechanisms behind the length variation in micro satellites is that of slipped-strand mispairing during dna polymerase mediated dna duplication. however, mutations involving indels of large copynumbers have been shown consistent with recombination-mediated events. in escherichia coli, an average mutational rate of × 10- was calculated over vntrs, and the rate seems to be dependent on intrinsic properties such as numbers of repeats. variation in copy numbers at vntr loci of certain sizes may affect the efficiency of promoters, thus affecting the coding potential of genes dependent on the genomic locations and the indels of repeats. for differentiation of bacterial isolates several vntr loci are combined in a multiple locus vntr analysis, a well established tool for epidemiological studies of bacteria with highly fit clonal dominance. one such bacteria is francisella tularensis, the agent for tularaemia. this species consist of clonal geographical subspecies that can be separated with the help of vntr markers. a similar approach could also be the solution for genotyping of francisella isolates belonging to fish pathogenic francisella spp. most mlva systems applied in epidemiology are performed by multiplexing pcr-assays where vntr size is deduced by capillary electrophoresis allowing rapid screening of large datasets. due to the relatively few vntr-loci included and that there were only isolates available for analysis in this study, sequencing was selected as a method for repeat number verification. the sequencing also detected single nucleotide polymorphisms in addition to allele size. the future potential for use of multiplex pcr assays with dyed primers and determining allele size by capillary electrophoresis is still an option and is recommended for larger datasets. this study presents a sequence based mlva system consisting of seven vntr loci identified in the f. noatunensis ssp. noatunensis gm isolate. when applied to francisella strains including norwegian and asian isolates of f. noatunensis, the mlva provided the best resolution shown for the fish-pathogenic francisellae so far. due to the preliminary genome status of the f. noatunensis ssp. noatunensis gm isolate a complete vntr locus search has not been performed, thus leaving the possibility for unidentified markers yet to be discovered. nonetheless strain typing of f. tularensis isolates has been achieved with the use of as few as two and six vntr markers. the exact location of the seven vntr markers in the f. noatunensis ssp. noatunensis is not known, however, blast searches indicate both intra- and intergenic location if the positions are homologous with that of the f. philomiragia atcc25017t isolate. the discriminatory power of the vntr loci used in this study, was not calculated as a simpson`s index of diversity. the reason for this is that the dataset does not fulfill the criteria of a test population proposed by van belkum et al., as isolates are predominately from outbreaks in cultured populations of cod with potential epizootic connections. the vntr markers had a typeability of 100% among the norwegian isolates, and with the exception of vntr- and - for the ua and ehime isolate respectively, all the vntr markers for f. noatunensis were amplified, indicating the usability of the typing system. stability testing showed variability within vntr- and - from f. noatunensis ssp. noatunensis gm isolate when grown in vitro for passages at 20°c and in vivo in cod at 18°c. such single repeat changes are shown to occur in 80% of mutational events with equal chance of being an insertion or a deletion. the sampled size was insufficient to assess the potential of hyper variable sites and mutation rates in the f. noatunensis ssp. noatunensis gm isolate. vntr mutational rates can differ across loci, and could be a problem in studies of phylogenetics and epizootics. stability and introduction of variation could be a result of bacterial phase variation and thus a result of environmental factors such as nutritional substances present in the growth medium, temperature effects, or just simply be stochastical events. however, four isolates from different francisellosis outbreaks at the same site during a two year period showed identical allelic profiles indicating field stability and epizootic coherence. the nj dendrogram presented in this study should not be viewed as a phylogenetic analysis as the use of vntrs and especially highly variable vntrs create noise in the phylogenetic signal. the stability testing in the study shows that this may be the case of the vntr- and vntr- there are also too few isolates and characters included in the dataset to attempt a proper analysis of genetic relationships. f. noatunensis strains a total of nine allele profiles were identified among the norwegian f. noatunensis ssp. noatunensis isolates. three of these contained more then one isolate and were identified using allelic profiles and visualized in the nj dendrogram. the remaining six isolates represent unique strains of which five originated from wild cod and one from farmed cod. although the dataset of norwegian isolates is small the trend is evident with only a few isolates occurring in outbreaks of francisellosis in cod culture, and a high diversity of strains among the isolates from wild cod. the chilean isolate and the norwegian cod isolates have previously been reported to be highly similar based on similarities in 16s sequence and in five housekeeping genes. however, according to the vntr analysis performed in this study, the ua isolate from salmon differs at six of seven markers compared to the norwegian f. noatunensis ssp. noatunensis isolates from cod, clearly separating them into two distinct branches in the nj dendrogram. the clear separation of these isolates most likely reflects geographic origin and the ecological differences between them. the allele profiles of the four f. noatunensis ssp. orientalis isolates were unique, separating them into four different genotypes. although the sampling set is small it correctly reflects the geographical origin and, possibly, also the variation in niche. there was no clear pattern in the clustering of the f. philomiragia isolates. care is needed if phylogenetic relationships are to be estimated based on vntrs, however the deeper nodes in the dendrogram a) between the species of f. philomiragia and f. noatunensis, b) between the subspecies f. noatunensis ssp. orientalis and f. noatunensis ssp. noatunensis and c) among the f. noatunensis ssp. noatunensis isolates of norway and chile, were similar to the tree topology obtained based on analysis of 16s rrna gene and datasets from housekeeping genes. epidemiology of the f. noatunensis considering the results of the present study some reflections regarding the epizootics of the f. noatunensis are possible. however, it is necessary to acknowledge the shortcomings of the dataset as it covers about 40% of the official diagnosed outbreaks of francisellosis in norway and that it only includes wild isolates collected in different counties. generally, the low allelic diversity among the isolates collected from cultured cod compared to the isolates from wild cod, indicate that the spreading could be a result of human activity. this is also supported by the fact that the production cycle of cod in norway involves transport of farmed cod along most of the norwegian coast. similarly, human activity has also been suggested as a cause for the geographical distribution of clonal isolates of bacillus anthracis. although the vntr markers provide a good starting point for studies of epizootics it may still be difficult to identify sources of infection as fish from different populations are mixed during the production, and, it is very difficult to avoid interactions between cultured and wild fish species. in the case of clade i with the isolates fnnr-001- -002- 004- and 017- the new generations of cod used for restocking after initial outbreak were probably infected by the initial isolate in the population at that site. the source of the initial infection was most likely cultured cod, since isolates belonging to the same genotype have been found at production sites in other counties, while the wild type isolate from the same county belonged to another clade. in clade iii, the isolates from a known broodfish supplier were identical to isolates from production sites where the aforementioned broodfish company was one of the suppliers of fish. hence, all of these cases could be a result of transmission through movement of offspring from broodfish. however, the situation is complicated by the fact that the sites used for on-growth are unknown. vertical transmission would greatly impact the disease management and needs to be addressed in further studies.  we present the first vntr analysis of fish pathogenic francisellae. seven polymorphic microsatellites were identified in a partial genome sequence of the f. noatunensis ssp. noatunensis gm isolate. a sequence based mlva system of these seven vntr-loci was applied to aquatic francisella isolates, including both f. philomiragia and f. noatunensis isolates. all vntrs were amplified in the f. noatunensis isolates, with the exception of vntr- in ua and vntr- in ehime- among the norwegian f. noatunensis ssp. noatunensis isolates, including isolates from both farmed and wild cod, a total of nine allelic profiles were identified. the majority of farmed isolates were divided into two allelic profiles, indicating low allelic variation in isolates from outbreaks in cod culture compared to isolates collected from wild cod. the allelic profile of the chilean f. noatunensis ssp. noatunensis isolate reflected the geographical and host divergence when compared to norwegian cod isolates. all f. noatunensis ssp. orientalis isolates show a unique allelic profile. the four vntrs amplified from f. philomiragia provided a unique allelic profile for all these isolates. the results show that this mlva system should provide a good starting point for future studies of epizootics of f. noatunensis. authors' contributions Øjb was responsible for carrying out this study and was the main contributor in writing the manuscript. kfo identified the vntr loci in the draft genome, contributed to project design and writing. an coordinated the project, contributed to project design and reviewed all drafts of the manuscript all authors read, commented on and approved the final manuscript. 
low allelic diversity in f. noatunensis isolates from outbreaks in cod culture compared to isolates wild cod, indicate that transmission of these isolates may be a result of human activity. however, the high genetic similarity within the francisella spp. since francisella noatunensis was first isolated from cultured atlantic cod in it has emerged as a global fish pathogen causing disease in both warm and cold water species. these vntr-loci were sequenced in f. noatunensis isolates collected from atlantic cod from norway, three-line grunt from japan, tilapia from indonesia and atlantic salmon from chile. a strain typing system for use in studies of f. noatunensis epizootics would be an important tool for disease management. the aim of this study is to identify possible useful vntrs in the genome of f. noatunensis. seven polymorphic vntr loci were identified in the preliminary genome sequence of f. noatunensis ssp. vntrs can be used to separate isolates belonging to both subspecies of f. noatunensis. the norwegian isolates presented in this study show both nine allelic profiles and clades, and that the majority of the farmed isolates belong in two clades only, while the allelic profiles from wild cod are unique. these loci are known as variable numbers of tandem repeat. the sequence based mlva system presented in this study should provide a good starting point for further development of a genotyping system that can be used in studies of epizootics and disease management of francisellosis.
 since francisella noatunensis was first isolated from cultured atlantic cod in it has emerged as a global fish pathogen causing disease in both warm and cold water species. outbreaks of francisellosis occur in several important cultured fish species making a correct management of this disease a matter of major importance. currently there are no vaccines or treatments available. a strain typing system for use in studies of f. noatunensis epizootics would be an important tool for disease management. however, the high genetic similarity within the francisella spp. makes strain typing difficult, but such typing of the related human pathogen francisella tullarensis has been performed successfully by targeting loci with higher genetic variation than the traditional signature sequences. these loci are known as variable numbers of tandem repeat. the aim of this study is to identify possible useful vntrs in the genome of f. noatunensis.  seven polymorphic vntr loci were identified in the preliminary genome sequence of f. noatunensis ssp. noatunensis gm isolate. these vntr-loci were sequenced in f. noatunensis isolates collected from atlantic cod from norway, three-line grunt from japan, tilapia from indonesia and atlantic salmon from chile. the norwegian isolates presented in this study show both nine allelic profiles and clades, and that the majority of the farmed isolates belong in two clades only, while the allelic profiles from wild cod are unique.  vntrs can be used to separate isolates belonging to both subspecies of f. noatunensis. low allelic diversity in f. noatunensis isolates from outbreaks in cod culture compared to isolates wild cod, indicate that transmission of these isolates may be a result of human activity. the sequence based mlva system presented in this study should provide a good starting point for further development of a genotyping system that can be used in studies of epizootics and disease management of francisellosis. 