International Journal of Microbiology

Volume 2014 (2014), Article ID 397143, 10 pages

http://dx.doi.org/10.1155/2014/397143

Research Article

Exploring Diversity among Norwegian Borrelia Strains Originating from Ixodes ricinus Ticks

Ann-Kristin Tveten

Faculty of Life Sciences, Aalesund University College, 6025 Ålesund, Norway

Received 7 March 2014; Accepted 13 August 2014; Published 27 August 2014

Academic Editor: Giuseppe Comi

Copyright © 2014 Ann-Kristin Tveten. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Characterisation of Borrelia strains from Ixodes ricinus ticks is important in the epidemiological surveillance of vector-borne pathogens. Multilocus sequences analysis (MLSA) is a molecular genotyping tool with high discriminatory power that has been applied in evolutionary studies and for the characterisation of Borrelia genospecies. MLSA was used to study genetic variations in Borrelia strains isolated from I. ricinus ticks collected from the woodlands in Skodje. The results demonstrate that the 50 Borrelia strains were separated into 36 sequence types (STs) that were not previously represented in the MLST database. A distance matrix neighbour-joining tree (bootstrapped 500 iterations) showed four deeply branched clusters, and each deeply branched cluster represented one Borrelia genospecies. The mean pairwise genetic differences confirm the genospecies clustering. The combination of alleles separates the Borrelia strains from northwest Norway from the strains in the MLST database, thus identifying new STs. Although a highly divergent B. afzelii population could be expected, the heterogeneity among the B. garinii strains is more unusual. The present study indicates that the circulation of strains between migrating birds and stationary birds in this coastal region may play a role in the evolution of B. garinii strains.

1. Introduction

Ticks are considered important carriers of pathogenic microorganisms in the northern hemisphere [1], and Ixodes ricinus ticks are associated with a diverse microbiota [2]. Borrelia spirochetes are one of the pathogens transmitted to humans by I. ricinus ticks [3], and Borrelia burgdorferi sensu lato (s.l.) represents a bacterial species complex that comprises several Borrelia genospecies associated with Lyme borreliosis (LB) [4]. Borrelia burgdorferi sensu stricto (s.s.), Borrelia afzelii, and Borrelia garinii are causative agents of human LB in Norway [5], and the prevalence of LB is monitored by the Norwegian Surveillance System for Communicable Diseases (MSIS). A fourth Borrelia genospecies, Borrelia valaisiana, was discovered in Norwegian ticks in 2010 [5]. These four genospecies have been detected in I. ricinus ticks in northwest Norway [6].

During the recent years, several genotyping methods have been established for analysis of Borrelia to characterise strains and determine the phylogenetic relationships between strains [7–11]. Multilocus sequence typing (MLST) is a molecular typing tool that is used to characterise pathogenic microorganisms [12], and MLST can be used as a tool for the epidemiological surveillance and tracking of infectious diseases [13]. MLST, as defined by Urwin and Maiden, is a genotyping tool based on the sequences of housekeeping genes that evolve slowly. The sequences used for MLST are approximately 400–500 bp in size and are located throughout the genome to avoid bias [13]. Multilocus sequence analysis (MLSA) applies MLST to characterise closely related species using a distance-based procedure for phylogenetic characterisation, and it includes pairwise genetic similarities. MLST/MLSA have been applied in population studies and for analysis of Borrelia species in geographic areas and for evolutionary studies and characterisation of new Borrelia species. While some MLST/MLSA schemes combine housekeeping genes with hypervariable regions and noncoding loci [9–11], the MLST scheme published in 2008 by Margos et al. is based on housekeeping genes that fulfill the strict criteria defined by Urwin and Maiden [7, 13]. The Borrelia MLST scheme is based on amplification, sequencing, and bioinformatic analysis of internal fragments of eight housekeeping genes (clpA, clpX, nifS, pepX, pyrG, rplB, recG, and uvrA) [7].

The Borrelia MLST scheme is available through the MLST network (http://www.mlst.net/). The MLST network enables easy access to sequences from all over the world [14]. Recently, studies of phylogenetic relationship of Borrelia genospecies using MLSA have resulted in the definition of two new Borrelia species, Borrelia bavariensis sp. [15] and Borrelia kurtenbachii [8]. These new findings demonstrate the high discriminatory power that MLST/MLSA schemes provide. Rudenko et al. utilised MLST to describe recombination at the nifS locus among B. burgdorferi s.s. and Borrelia americana strains. Their findings indicated that, to a degree, strain diversity is influenced by the host [16]. The aim of this study is to utilise MLSA to characterise Borrelia strains isolated from I. ricinus ticks collected from the fauna in northwest Norway. Multilocus sequence analysis of Norwegian strains could provide knowledge about the strain diversity among Borrelia strains found in this region and describe their evolution compared to European strains.

2. Materials and Methods

2.1. Samples

DNA from questing I. ricinus ticks collected in 2012 (531 nymphs/71 adults) and 2013 (539 nymphs/64 adults) (Figure 1) was isolated, and the DNA samples were analysed by qPCR to detect the presence of Borrelia genospecies as previously described [6]. A total of 86 samples (79 nymphs/7 adults) were positive for Borrelia in 2012, and 89 samples (83 nymphs/8 adults) were positive for Borrelia in 2013. Multilocus sequence analysis was performed using the MLST scheme developed by Margos et al. A total of 50 samples were amplified across all eight housekeeping genes (clpA, clpX, nifS, pepX, pyrG, recG, rplB, and uvrA) [7] and were submitted for sequencing. The 50 strains used in this study are listed in Table 1.

Table 1: Allele and sequence types (STs).

Source
Species
clpA
clpX
nifS
pepX
pyrG
recG
rplB
uvrA
ST

I. ricinus
B. afzelii
109
24
23
85
22
29
23
29
481

I. ricinus
B. afzelii
36
24
96
31
119
29
23
28
483

I. ricinus
B. afzelii
36
24
96
31
22
30
23
28
484

I. ricinus
B. garinii
46
30
33
42
29
88
79
36
485

I. ricinus
B. garinii
47
30
33
42
111
88
79
36
486

I. ricinus
B. garinii
47
30
33
42
111
88
79
36
486

I. ricinus
B. garinii
46
30
33
38
32
76
32
36
487

I. ricinus
B. garinii
46
30
29
38
32
76
32
36
488

I. ricinus
B. afzelii
37
24
23
35
96
30
99
142
489

I. ricinus
B. afzelii
51
82
23
86
85
27
28
29
490

I. ricinus
B. afzelii
51
82
23
86
85
27
29
29
491

I. ricinus
B. afzelii
172
24
23
85
85
27
29
29
492

I. ricinus
B. afzelii
37
135
24
87
22
29
23
28
493

I. ricinus
B. garinii
45
30
32
96
29
38
30
36
494

I. ricinus
B. burgdorferi s.s.
14
8
14
10
10
137
9
129
495

I. ricinus
B. afzelii
36
24
23
35
26
33
23
28
496

I. ricinus
B. afzelii
36
24
23
32
22
33
23
28
497

I. ricinus
B. garinii
46
30
29
96
32
76
79
33
498

I. ricinus
B. garinii
95
74
34
96
32
88
79
85
499

I. ricinus
B. afzelii
36
82
23
32
119
27
23
28
500

I. ricinus
B. afzelii
36
82
23
87
22
27
23
28
501

I. ricinus
B. afzelii
109
24
23
35
22
98
23
29
502

I. ricinus
B. afzelii
109
24
23
85
22
98
23
29
503

I. ricinus
B. afzelii
36
24
24
87
22
29
29
30
504

I. ricinus
B. afzelii
36
24
23
31
22
29
29
30
505

I. ricinus
B. afzelii
36
24
23
31
20
29
29
30
506

I. ricinus
B. afzelii
36
24
23
32
20
29
29
28
507

I. ricinus
B. afzelii
36
24
96
32
92
33
29
28
508

I. ricinus
B. afzelii
37
135
24
87
22
29
23
28
493

I. ricinus
B. afzelii
37
24
23
35
96
30
99
142
489

I. ricinus
B. afzelii
36
24
96
31
119
29
23
28
483

I. ricinus
B. afzelii
37
135
24
87
22
29
23
28
509

I. ricinus
B. afzelii
109
24
23
85
22
29
23
29
481

I. ricinus
B. afzelii
36
24
23
31
22
29
29
30
505

I. ricinus
B. garinii
46
74
29
38
35
40
32
33
510

I. ricinus
B. garinii
46
74
29
38
35
40
32
33
510

I. ricinus
B. garinii
45
30
32
96
29
38
30
36
494

I. ricinus
B. garinii
95
74
34
96
34
76
79
85
511

I. ricinus
B. garinii
45
30
32
96
29
38
30
36
494

I. ricinus
B. afzelii
36
24
23
35
85
29
28
29
512

I. ricinus
B. afzelii
51
24
23
86
85
30
99
142
513

I. ricinus
B. afzelii
36
24
23
35
85
29
28
29
512

I. ricinus
B. afzelii
51
82
23
86
119
29
28
29
514

I. ricinus
B. garinii
95
74
34
96
32
76
79
85
515

I. ricinus
B. garinii
46
30
29
96
32
76
32
33
516

I. ricinus
B. garinii
46
30
29
96
32
76
32
33
516

I. ricinus
B. valaisiana
96
76
36
98
84
44
78
86
517

I. ricinus
B. valaisiana
50
76
36
98
38
44
78
86
518

I. ricinus
B. valaisiana
96
76
36
98
84
44
78
86
517

I. ricinus
B. valaisiana
96
76
36
98
84
44
78
86
517

Figure 1: Map of the collection site. The map indicates the location site in northwest Norway and the tracking route used. The altitude of each collection site and the distance between the sites are indicated in the bottom part of the map.

2.2. Amplification and Sequencing

The amplification of each housekeeping gene was carried out in a 15 µL reaction volume containing 7.5 µL Mastermix 2.0x SYBR GREEN (Applied Biosystems), 0.75 µL of each primer (10 µM), 3.00 µL ddH₂O, and 3.0 µL template DNA. Amplification of the eight housekeeping genes (clpA, clpX, nifS, pepX, pyrG, recG, rplB, and uvrA) was performed with primers and nested touchdown PCR conditions for clpA, nifS, pepX, pyrG, recG, and uvrA as previously described [7, 8]. Amplification of clpX and rplB was performed with primers described by Margos et al. combined with a nested touchdown PCR composed of 1 cycle of denaturation (10 min, 95°C), followed by 9 cycles of denaturation (30 sec, 95°C), touchdown annealing (30 sec, 58°C–50°C), and extension (1 min, 72°C) and 30 cycles of denaturation (30 sec, 95°C), annealing (30 sec, 50°C), and extension (1 min, 72°C) with a final cycle of extension (10 min, 72°C). Gel electrophoresis was used to confirm the size and approximate concentration of the PCR products for each amplified gene. Gel electrophoresis was performed with a 2% Tris-acetate EDTA buffered (TBE) gel followed by staining with GelRed (Affymetrix, Santa Clara, USA).

The PCR products were purified in a two-step process by ExoSAP-IT (Affymetrix, Santa Clara, USA) according to the manufacturers’ protocol. The sequencing reaction was carried out in a 10 µL reaction volume containing 1 µL Big Dye Terminator 3.1 (Applied Biosystems, Carlsbad, USA), 1 µL 5x Sequencing Buffer (Applied Biosystems, Carlsbad, USA), 0.32 µL primer (10 µM), 4.68 µL ddH₂O, and 3 µL template DNA. Amplification was performed in a 2720 thermal cycler (Applied Biosystems, Carlsbad, USA) with 1 cycle of denaturation (6 min, 96°C), followed by 25 cycles of denaturation (10 sec, 96°C), annealing (5 sec, 50°C), and extension (4 min, 60°C). Finally, 10 µL ddH₂O was added to each well upon completion of the sequencing reaction, and the samples were stored at –20°C until sequencing. Sequencing was performed by Eurofins MWG Operon (Germany).

2.3. Alignment and Multilocus Sequence Analysis

The sequences for each housekeeping gene were aligned by the MUSCLE algorithm and concatenated in MEGA 5.1 [17]. The sequences were analysed with BLAST to confirm the correct Borrelia genospecies [18]. Each locus and the concatenated sequence of all eight loci were analysed. The MEGA 5.1 software was used to draw a distance matrix neighbour-joining tree (bootstrapped 500 iterations) with the concatenated DNA sequences. Sequences from B. duttonii were added from the BLAST database [18] to root the neighbour-joining tree. B. duttonii is a Borrelia species that causes relapsing fever [19]. Sequences from four reference species from the MLST database were included in the analysis, B. burgdorferi s.s. (B31), B. afzelii (VS461), B. valaisiana (70865B), and B. garinii (20047). To compare allele sequences and sequences types (STs) and to identify identical sequences, the nonredundant database (NRDB) written by Warren Gish at Washington University was used (http://pubmlst.org/analysis/). All known alleles and STs were retrieved from the Borrelia MLST database for comparison (http://borrelia.mlst.net/). STs that were not found in the MLST database were submitted to the MLST database for proper numbering. Statistical analysis was performed with MEGA 5.1. The mean pairwise genetic difference (within and between groups) was calculated using the Kimura 2-parameter model. Average nonsynonymous substitutions and synonymous substitutions (dN/dS) were calculated using the modified Nei-Gojobori method (Jukes-Cantor). The number of polymorphic sites (PS) and the nucleotide diversity (π) were calculated using Tajima’s test of neutrality. The index of association, both standard and classic, was calculated using START2 [20].

To compare the relationship between Norwegian Borrelia strains and European Borrelia strains a goeBURST minimum spanning tree was generated using PHYLOViZ 1.0 software [21].

3. Results

A total of 50 strains were sequenced across all eight loci for the present study. A distance matrix neighbour-joining tree (bootstrapped 500 iterations) (Figure 2) demonstrated that the 50 strains were separated into 36 sequence types (STs) (Table 1) that were not previously found in the MLST database. The neighbour-joining tree formed four deeply branched clusters, and, between the four clusters, there were no shared alleles. All strains cluster together with their respective reference species and contain combinations of alleles that are new entries into the MLST database. The mean pairwise genetic difference confirms the genospecies clustering (Figure 3) based on the genospecies threshold (0.017) as defined by Margos et al. [15]. Table 2 summarises the characteristics of each of the eight housekeeping genes included in the multilocus sequencing analysis of Borrelia strains.

Table 2: Characteristics of each housekeeping gene.

Species
Gene
Size (bp)
Avg. G+C content %
PS (%)a

d/d

B. burgdorferi s.s.
clpA
579
52.2
—d
—
—

B. afzelii
0,69
0,001759
0,001810

B. garinii
1,73
0,005484
0,005205

B. valaisiana
0,69
0,003454
0,003768

B. burgdorferi s.s.
clpX
624
56.8
—
—
—

B. afzelii
0,48
0,001255
0,00122

B. garinii
6,08
0,02191
0,001835

B. valaisiana
—e
—
—

B. burgdorferi s.s.
nifS
564
55.8
—
—
—

B. afzelii
0,36
0,00087
0,00090

B. garinii
1,06
0,0054
0,00434

B. valaisiana
—
—
—

B. burgdorferi s.s.
pepX
570
55.6
—
—
—

B. afzelii
0,88
0,00328
0,00223

B. garinii
1,23
0,00484
0,00501

B. valaisiana
—
—
—

B. burgdorferi s.s.
pyrG
603
51.0
—
—
—

B. afzelii
0,83
0,00240
0,00211

B. garinii
1,49
0,00805
0,00609

B. valaisiana
1,32
0,0073
0,00637

B. burgdorferi s.s.
recG
651
53.0
—
—
—

B. afzelii
0,62
0,00209
0,00157

B. garinii
1,54
0,00545
0,00463

B. valaisiana
—
—
—

B. burgdorferi s.s.
rplB
624
56.9
—
—
—

B. afzelii
0,61
0,00261
0,00143

B. garinii
1,28
0,00669
0,00386

B. valaisiana
—
—
—

B. burgdorferi s.s.
uvrA
570
54.4
—
—
—

B. afzelii
0,88
0,00234
0,00223

B. garinii
1,05
0,00335
0,003172

B. valaisiana
—
—
—

Percent polymorphic sites (PS (%)).

bNucleotide diversity ().

cAverage nonsynonymous substitutions versus synonymous substitutions (d/d).

dIt was not possible to perform statistics with only one isolate in this group.

eIt was not possible to perform statistics with only one allele variant.

Figure 2: Description of clades based on the neighbour-joining tree (bootstrapped 500 iterations) calculated based on concatenated DNA sequences showing the relationship between B. burgdorferi s.l. strains.

Figure 3: Mean pairwise difference within and between Borrelia genospecies. The line indicates species threshold (0.017) as defined by Margos et al. [15].

The comparative goeBURST analysis (level 8 grouping) displays the relationship between the Norwegian strains and European strains (Figure 4). The goeBURST analysis shows that the Norwegian Borrelia population is the most closely related to the Latvian and French populations. The dataset contains 175 unique STs and information for approximately 246 strains. A goeBURST analysis was performed at the SLV (single locus variant) level by dividing the dataset into 32 groups with 2 or more STs and 64 singletons. The largest clonal complex 0 (CC0) has 13 unique STs and represents 33 strains. CC0 represents strains identified as B. valaisiana and contains 6 SLVs. All of the SLVs share the pepX, pyrG, rplB, and uvrA alleles. In CC1 (9 STs, 22 strains), the predicted ancestral genotype is ST 86, with 8 SLVs representing 10 strains; thus this is the highest ranking ancestral genotype. All SLVs have the same clpA, clpX, and rplB alleles. A goeBURST analysis conducted at the DLV (double locus variant) level divided the dataset into 54 CCs consisting of 26 groups with 2 or more STs and 28 singletons. A goeBURST analysis conducted at the TLV (three loci variant) level divided the dataset into 33 CCs consisting of 21 groups with 2 or more STs and 11 singletons. A level 4 grouping was required for B. valaisiana to be represented by one CC in the goeBURST analysis, and a level 5 grouping was required for B. afzelii to be represented by one CC. A level 7 grouping was required for B. garinii and B. burgdorferi s.s. to be represented by one CC.

Figure 4: goeBURST. The goeBURST is coloured based on the geographic origin of each ST. Multiple colours indicated multiple origins.

The Norwegian STs are characterised as both SLVs and DLVs. The Norwegian strains are spread throughout the goeBURST tree and are clustered together with their respective genospecies from other European countries. Four branches are formed from the primary founder. These four branches contain STs related to the four genospecies, B. garinii, B. afzelii, B. valaisiana, and B. burgdorferi s.s. The highest ranking ancestral genotype, ST 86, is a B. garinii strain from Latvia.

4. Discussion

Northwest Norway is characterised by geographical features such as islands, several deep fjords, and high mountains. The landscape shows signs of overgrowth and little maintenance of the vegetation, thus providing suitable growth conditions for diverse wildlife. More than 4000 different Animalia species have been observed in this region [22]. The region has mild winters and humid summers, creating good conditions for Ixodes ricinus ticks. Borrelia spirochetes are maintained in nature by the circulation between the vector and the host, and the wide range of available hosts may influence the diversity among strains of Borrelia [23]. Borrelia spirochetes have an obligate parasitic lifestyle, and their circulation depends on an interaction between reservoir hosts and vectors [24]. Birds are reservoir hosts for B. garinii while B. afzelii are usually associated with rodents [15, 25]. Evolution and a wide range of determinants that affect reservoir hosts and vectors are likely to influence the diversity among strains differently [23, 24, 26]. Reservoir hosts and vectors are submitted to a number of factors such as climate, nutrition, and pollution that have different impacts depending on the geographic locations [23, 24]. In addition to the influence that reservoir hosts provide, the complex microbial communities associated with I. ricinus ticks also subject Borrelia spirochetes to microbial competition and natural selection. These microbial communities are composed of bacteria, endosymbionts, viruses, and protozoa from the tick’s natural habitat and hosts, and, therefore, they vary between different geographical locations [2, 27, 28].

Multilocus sequence analysis of the Borrelia strains found in I. ricinus ticks from this region shows signs of high diversity. The majority of the strains are B. afzelii and B. garinii, the most prevalent genospecies within this region [6]. Fifty Borrelia strains were separated into 36 STs, and the genetic relationship between the strains was visualised in a neighbour-joining tree. The first deeply branched cluster contains strains identified as B. afzelii, and 29 strains were separated into 22 STs. These 29 B. afzelii strains form a highly structured population consisting of STs that have not been previously described in the MLST database. In 2011, Vollmer et al. demonstrated that there is limited movement of B. afzelii strains between geographic locations and that only two B. afzelii STs have been identified in more than one geographic location [22, 29]. B. afzelii strains are described as heterogenic, and a comparative goeBURST analysis of strains from Scotland and strains from England showed that there is a limited movement of strains between the north and south in the United Kingdom [24, 29]. The dispersal of Borrelia spirochetes is linked to the movement of their hosts. B. afzelii is associated with small rodent hosts with a limited movement pattern, thus creating a limited dispersal of B. afzelii strains [30]. All strains included in the present study were collected from four sites within one geographically defined area (Figure 1). Some of the B. afzelii strains were identified in more than one sample, but they were always within the collection site.

The second deeply branched cluster contains four B. valaisiana strains divided into 2 STs that represent new entries into the MLST database. The third deeply branched cluster contains 16 B. garinii strains divided into 11 STs that are new entries into the MLST database. Both B. valaisiana and B. garinii are bird-associated genospecies, and previous studies have found that these two Borrelia genospecies have sequence types that have been identified in different geographic locations throughout Europe [29]. Vollmer et al. reported that B. valaisiana strains are not completely homogenised within Europe while B. garinii shows signs of being a mixed population across central Europe [29]. Migrating birds play an important role in the dispersal of B. garinii strains between geographic locations. Kjelland et al. studied the prevalence of Borrelia infection in 31 species of migrating birds and 6 species of resident birds in southern Norway. Borrelia infection was found in all bird species [31]. Studies also suggest that B. garinii strains have the ability to circulate both in terrestrial cycles and in seabirds with their associated tick Ixodes uriae [32–34]. A study by Comstedt et al. from 2011 determined that B. garinii was found in Norwegian populations of seabirds and I. uriae [35]. Sequencing of the intergenic spacer (IGS) from B. garinii strains has indicated that the seabird-associated B. garinii strains are more diverse than the terrestrial B. garinii strains [36]. All strains included in the present study are collected from a geographic location with a high number of observed bird species. In Møre and Romsdal county, 290 bird species have been observed in 2013 by ornithologist [37]. These birds represent a variety of different reservoir hosts for B. garinii, and the nesting on the island provides a possibility of transferring B. garinii into fauna, inducing circulation between migrating and stationary hosts. The high number of migrating birds along the coastal area of northwest Norway enables introduction of B. garinii strains from multiple different geographic origins and multiple reservoir hosts and the process of circulation between reservoir hosts has been suggested to play a role in the genetic evolution of Borrelia genospecies [29]. Some B. garinii strains from the present study are DLVs and TLVs of the strains from central Europe, while others have levels 4 and 5 connections to the central European STs. Within the Norwegian B. garinii strains, some are SLVs and DLVs, and some are more distantly related. This may indicate that the different strains have been introduced into the fauna by migrating birds at different times and that some strains might have been maintained by stationary hosts for some time. Calculation of the nucleotide diversity (π value) reveals that the π values were at least two times greater for B. garinii than for B. afzelii. Assuming a similar generation time between B. garinii and B. afzelii, the different reservoir hosts might have influenced the evolutionary rate differently between genospecies [26]. The limited geographic dispersal of B. afzelii strains has allowed them to evolve differently between geographic locations, but the nucleotide diversity within a geographic location is relatively low. B. garinii associated with migratory birds is subjected to a higher degree of ecological determinants, and the nucleotide diversity may reflect that strains have been introduced to the fauna at different times.

The fourth deep branched cluster contains one B. burgdorferi s.s. isolate and its associated reference species. The B. burgdorferi s.s. strain was isolated from an adult female tick. The low prevalence of this genospecies [6], in combination with certain challenges associated with amplification across all eight loci, caused this strain to be the only B. burgdorferi s.s. included in this study.

The number of polymorphic sites (PS) and the nucleotide diversity (π) were calculated to further describe the Borrelia strains. A sequence analysis was conducted to identify point mutations found in multiple strains. These point mutations have a minor contribution to the genetic diversity between concatenated DNA sequences from northwest Norway compared to other worldwide strains. Most of the point mutations are silent or missense mutations with little or no effect on the genetic function. The ratio between nonsynonymous substitutions and synonymous substitutions (dN/dS) < 1 indicates that the number of synonymous substitutions is higher than the number of nonsynonymous substitutions. This indicates that the genes are not under selective pressure and are suitable for MLSA analysis. The index of association describes the influence of recombination in a population and was calculated in two ways: using the classical method defined by Smith et al. and using the standard method defined by Haubold and Hudson [38, 39]. The would be significant for a clonal population. The calculated values for this population were and , and significant linkage disequilibrium was detected. Estimation of the index of association indicated that the population is clonal and that recombination has not played an important role in the evolution of this population.

Multilocus sequencing analysis has provided knowledge about the genetic variations between Borrelia genospecies originating from I. ricinus ticks from northwest Norway. The combinations of alleles separate these strains from other strains in the MLST database. This is the first multilocus sequence analysis of Borrelia strains identified from Norwegian I. ricinus ticks, and the 50 Borrelia strains from Norway represent a large contribution to the MLST database. This enabled comparative studies between Norwegian and European strains, thus creating a better understanding of the strain development of Borrelia genospecies in different geographic locations in Europe. The high level of divergence in the B. afzelii population was expected, but the divergence within the B. garinii population differs from previous MLST/MLSA studies. Migrating birds are constantly introducing new strains to the fauna along the coastline in northwest Norway, and circulation between migrating and stationary hosts would influence the strain diversity. Some of the B. garinii strains show signs of being newly introduced to the fauna while others have most likely been maintained in stationary hosts for a longer period of time.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

The author would like to thank Synnøve T. Broum for the help in collecting live ticks.

References

J. S. Gray, “The biology of Ixodes ticks, with special reference to Ixodes ricinus,” Society for Zoonotic Ecology and Epidemiology, http://www.zooeco.org/

A. K. Tveten and K. K. Sjåstad, “Identification of bacteria infecting ixodes ricinus ticks by 16S rDNA amplification and denaturing gradient gel electrophoresis,” Vector-Borne and Zoonotic Diseases, vol. 11, no. 10, pp. 1329–1334, 2011.

U. R. Hengge, A. Tannapfel, S. K. Tyring, R. Erbel, G. Arendt, and T. Ruzicka, “Lyme borreliosis,” The Lancet Infectious Diseases, vol. 3, no. 8, pp. 489–500, 2003.

M. E. Aguero-Rosenfeld, G. Q. Wang, I. Schwartz, and G. P. Wormser, “Diagnosis of lyme borreliosis,” Clinical Microbiology Reviews, vol. 18, no. 3, pp. 484–509, 2005.

V. Kjelland, S. Stuen, T. Skarpaas, and A. Slettan, “Prevalence and genotypes of Borrelia burgdorferi sensu lato infection in Ixodes ricinus ticks in southern Norway,” Scandinavian Journal of Infectious Diseases, vol. 42, no. 8, pp. 579–585, 2010.

A. K. Tveten, “Prevalence of Borrelia burgdorferi sensu stricto, Borrelia afzelii, Borrelia garinii, and Borrelia valaisiana in Ixodes ricinus ticks from the northwest of Norway,” Scandinavian Journal of Infectious Diseases, vol. 45, no. 9, pp. 681–687, 2013.

G. Margos, A. G. Gatewood, D. M. Aanensen et al., “MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi,” Proceedings of the National Academy of Sciences of the United States of America, vol. 105, no. 25, pp. 8730–8735, 2008.

G. Margos, A. Hojgaard, R. S. Lane et al., “Multilocus sequence analysis of Borrelia bissettii strains from North America reveals a new Borrelia species, Borrelia kurtenbachii,” Ticks and Tick-Borne Diseases, vol. 1, no. 4, pp. 151–158, 2010.

D. Richter, D. Postic, N. Sertour, I. Livey, F. R. Matuschka, and G. Baranton, “Delineation of Borrelia burgdorferi sensu lato species by multilocus sequence analysis and confirmation of the delineation of Borrelia spielmanii sp. nov,” International Journal of Systematic and Evolutionary Microbiology, vol. 56, no. 4, pp. 873–881, 2006.

D. Postic, M. Garnier, and G. Baranton, “Multilocus sequence analysis of atypical Borrelia burgdorferi sensu lato isolates—description of Borrelia californiensis sp. nov., and genomospecies 1 and 2,” International Journal of Medical Microbiology, vol. 297, no. 4, pp. 263–271, 2007.

W.-G. Qiu, J. F. Bruno, W. D. McCaig et al., “Wide distribution of a high-virulence Borrelia burgdorferi clone in Europe and North America,” Emerging Infectious Diseases, vol. 14, no. 7, pp. 1097–1104, 2008.

M. C. J. Maiden, J. A. Bygraves, E. Feil et al., “Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms,” Proceedings of the National Academy of Sciences of the United States of America, vol. 95, no. 6, pp. 3140–3145, 1998.

R. Urwin and M. C. J. Maiden, “Multi-locus sequence typing: a tool for global epidemiology,” Trends in Microbiology, vol. 11, no. 10, pp. 479–487, 2003.

D. M. Aanensen and B. G. Spratt, “The multilocus sequence typing network: mlst.net,” Nucleic Acids Research, vol. 33, no. 2, pp. W728–W733, 2005.

G. Margos, S. A. Vollmer, M. Cornet et al., “A new Borrelia species defined by multilocus sequence analysis of housekeeping genes,” Applied and Environmental Microbiology, vol. 75, no. 16, pp. 5410–5416, 2009.

N. Rudenko, M. Golovchenko, N. M. Belfiore, L. Grubhoffer, and J. H. Oliver Jr., “Divergence of Borrelia burgdorferi sensu lato spirochetes could be driven by the host: diversity of Borrelia strains isolated from ticks feeding on a single bird,” Parasit Vectors, vol. 7, article 4, 2014.

K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011.

S. F. Altschul, T. L. Madden, A. A. Schäffer et al., “Gapped BLAST and PSI-BLAST: a new generation of protein database search programs,” Nucleic Acids Research, vol. 25, no. 17, pp. 3389–3402, 1997.

M. Fukunaga, Y. Ushijima, L. Y. Aoki, and A. Talbert, “Detection of Borrelia duttonii, a tick-borne relapsing fever agent in central Tanzania, within ticks by flagellin gene-based nested polymerase chain reaction,” Vector-Borne and Zoonotic Diseases, vol. 1, no. 4, pp. 331–338, 2001.

K. A. Jolley, E. J. Feil, M.-S. Chan, and M. C. J. Maiden, “Sequence type analysis and recombinational tests (START),” Bioinformatics, vol. 17, no. 12, pp. 1230–1231, 2002.

A. P. Francisco, C. Vaz, P. T. Monteiro, J. Melo-Cristino, M. Ramirez, and J. A. Carriço, “PHYLOViZ: phylogenetic inference and data visualization for sequence based typing methods,” BMC Bioinformatics, vol. 13, no. 1, article 87, 2012.

Norway TNBICatGBIF, Species Map Service, The Norwegian Biodiversity Information Centre and the Global Biodiversity Information Facility Norway, 2007.

K. Kurtenbach, K. Hanincová, J. I. Tsao, G. Margos, D. Fish, and N. H. Ogden, “Fundamental processes in the evolutionary ecology of Lyme borreliosis,” Nature Reviews Microbiology, vol. 4, no. 9, pp. 660–669, 2006.

G. Margos, S. A. Vollmer, N. H. Ogden, and D. Fish, “Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato,” Infection, Genetics and Evolution, vol. 11, no. 7, pp. 1545–1563, 2011.

K. Kurtenbach, S. de Michelis, S. Etti et al., “Host association of Borrelia burgdorferi sensu lato: the key role of host complement,” Trends in Microbiology, vol. 10, no. 2, pp. 74–79, 2002.

S. A. Vollmer, E. J. Feil, C.-Y. Chu et al., “Spatial spread and demographic expansion of Lyme borreliosis spirochaetes in Eurasia,” Infection, Genetics and Evolution, vol. 14, no. 1, pp. 147–155, 2013.

L. van Overbeek, F. Gassner, C. L. van der Plas, P. Kastelein, U. Nunes-Da Rocha, and W. Takken, “Diversity of Ixodes ricinus tick-associated bacterial communities from different forests,” FEMS Microbiology Ecology, vol. 66, no. 1, pp. 72–84, 2008.

C. Schabereiter-Gurtner, W. Lubitz, and S. Rölleke, “Application of broad-range 16S rRNA PCR amplification and DGGE fingerprinting for detection of tick-infecting bacteria,” Journal of Microbiological Methods, vol. 52, no. 2, pp. 251–260, 2003.

S. A. Vollmer, A. Bormane, R. E. Dinnis et al., “Host migration impacts on the phylogeography of Lyme Borreliosis spirochaete species in Europe,” Environmental Microbiology, vol. 13, no. 1, pp. 184–192, 2011.

K. Kurtenbach, S. de Michelis, H.-S. Sewell et al., “The key roles of selection and migration in the ecology of Lyme borreliosis,” International Journal of Medical Microbiology, vol. 291, no. 33, pp. 152–154, 2002.

V. Kjelland, S. Stuen, T. Skarpaas, and A. Slettan, “Borrelia burgdorferi sensu lato in Ixodes ricinus ticks collected from migratory birds in Southern Norway,” Acta Veterinaria Scandinavica, vol. 52, article 59, 2010.

C. Larsson, P. Comstedt, B. Olsen, and S. Bergström, “First record of lyme disease Borrelia in the arctic,” Vector-Borne and Zoonotic Diseases, vol. 7, no. 3, pp. 453–456, 2007.

B. Olsen, T. G. T. Jaenson, L. Noppa, J. Bunikis, and S. Bergstrom, “A Lyme borreliosis cycle in seabirds and Ixodes uriae ticks,” Nature, vol. 362, no. 6418, pp. 340–342, 1993.

J. Bunikis, B. Olsén, V. Fingerle, J. Bonnedahl, B. Wilske, and S. Bergström, “Molecular polymorphism of the Lyme disease agent Borrelia garinii in Northern Europe is influenced by a novel enzootic Borrelia focus in the North Atlantic,” Journal of Clinical Microbiology, vol. 34, no. 2, pp. 364–368, 1996.

P. Comstedt, T. Jakobsson, and S. Bergstrom, “Global ecology and epidemiology of Borrelia garinii spirochetes,” Infection Ecology & Epidemiology, vol. 1, article 9545, 2011.

P. Comstedt, L. Asokliene, I. Eliasson et al., “Complex population structure of lyme borreliosis group spirochete Borrelia garinii in subarctic eurasia,” PLoS ONE, vol. 4, no. 6, Article ID e5841, 2009.

Norsk ornitologi forening aMoR, Artsobservasjoner (Speices observation), 2013

J. M. Smith, N. H. Smith, M. O'Rourke, and B. G. Spratt, “How clonal are bacteria?” Proceedings of the National Academy of Sciences of the United States of America, vol. 90, no. 10, pp. 4384–4388, 1993.

B. Haubold and R. R. Hudson, “LIAN 3.0: detecting linkage disequilibrium in multilocus data,” Bioinformatics, vol. 16, no. 9, pp. 847–848, 2000.