Genetic variation in the mitochondrial 16S ribosomal RNA gene of Ixodes scapularis (Acari: Ixodidae)
Parasites & Vectors volume 7, Article number: 530 (2014)
Ixodes scapularis is a vector of several human pathogens in the United States, and there is geographical variation in the relative number of persons infected with these pathogens. Geographically isolated populations of I. scapularis have established or are in the process of establishing in southern Canada. Knowledge of the genetic variation within and among these populations may provide insight into their geographical origins in the United States and the potential risk of exposure of Canadians to the different pathogens carried by I. scapularis.
Part of the mitochondrial (mt) 16S ribosomal (r) RNA gene was amplified by PCR from 582 ticks collected from southern Canada, and Minnesota and Rhode Island in the United States. Sequence variation was examined in relation to the predicted secondary structure of the gene. Genetic diversity among populations was also determined.
DNA sequence analyses revealed 52 haplotypes. Most mutational alterations in DNA sequence occurred at unpaired sites or represented partial compensatory base pair changes that maintained the stability of the secondary structure. Significant genetic variation was detected within and among populations in different geographical regions. A greater proportion of the haplotypes of I. scapularis from the Canadian Prairie Provinces were found in the Midwest of the United States than in other regions, whereas more of the haplotypes of I. scapularis from the Canadian Central and Atlantic Provinces occurred in the Northeast of the United States. Nonetheless, 58% of I. scapularis were of a haplotype that occurs in the Midwest and Northeast of the United States; thus, their geographical origins could not be determined.
There is considerable genetic variation in the mt 16S rRNA gene of I. scapularis. There is some evidence to support the hypothesis that some lineages of I. scapularis in the Atlantic and Central Provinces of Canada may be derived from colonizing individuals originating in the Northeast of the United States, whereas those in the Prairie Provinces may be derived from individuals originating in the Midwest of the United States. However, additional genetic markers are needed to test hypotheses concerning the geographical origins of I. scapularis in Canada.
Ticks are important vectors of human and animal pathogens ,. The incidences of tick-borne diseases are increasing -, due in part, to the expansion of the distribution of some tick species into new geographical areas . Over the past 10 years, the distribution of the blacklegged tick (Ixode s scapularis) in North America has continued to expand in the Upper Midwest of the United States ,, and in southern Canada -. This range expansion is important from a public health perspective because I. scapularis is the principal vector of Borrelia burgdorferi sensu stricto, the causative agent of Lyme disease in North America . Blacklegged ticks are also important vectors of the etiologic agents of human granulocytic anaplasmosis (Anaplasma phagocytophilum) , human and rodent babesiosis (Babesia microti) , and tick-borne encephalitis (Powassan virus) . The relative occurrences of these diseases, and the prevalences of the different strains of pathogens, vary throughout the distributional range of I. scapularis,-. Understanding the evolutionary ecology of these vector-borne diseases requires detailed knowledge of the biology, ecology, and population genetics of the vector and the pathogens it carries .
The distribution of I. scapularis in the United States is divided into three geographically isolated foci: the Northeast, Midwest, and the South -. There are also geographically isolated populations of I. scapularis in several provinces in southern Canada. The first of these Canadian populations established at the Long Point peninsula, which includes the Long Point Provincial Park and adjacent National Wildlife area (Ontario) was described in the early 1970s , whereas the next two populations, at Point Peele National Park and Rondeau Provincial Park (Ontario), did not establish until the 1990s ,. Other populations have now become established or are in the process of establishing in Ontario, Nova Scotia, New Brunswick, Manitoba and Quebec ,,,,. It has been proposed that migratory passerines are transporting large numbers of I. scapularis larvae and nymphs into Canada from the United States each year during their spring migration ,.
Studies have shown there are two major lineages or clades of I. scapularis in the United States -. Individuals of the Southern clade have only been reported from North Carolina, South Carolina, Georgia, Oklahoma, Texas, Arkansas, and Florida, whereas those of the American clade occur primarily in the Northeast and Midwest, but also occur in some southern states ,-. Blacklegged ticks in southern Canada were also shown to belong to the American clade based on analyses of the DNA sequences of part of Domains IV and V of the mitochondrial (mt) 16S ribosomal RNA (rRNA) gene . The absence of individuals of the Southern clade, combined with significant differences in genetic structure among six established populations, suggested that I. scapularis populations in southern Canada were founded by colonizing individuals that originated from different populations in the Northeast and Midwest of the United States . In addition, eight of the 19 haplotypes detected among I. scapularis in southern Canada had not been reported previously in studies conducted in the United States. These eight haplotypes also represented 27% of all ticks characterized in southern Canada . This suggested that the extent of the genetic diversity in I. scapularis, based on the DNA sequences of Domains IV and V of the mt 16S rRNA gene, was more extensive than previously thought. Therefore, in the present study, we assessed the extent of the variation in the DNA sequences of the mt 16S rRNA gene for I. scapularis, both within and among populations in southern Canada, and the Midwest and Northeast of the United States. We also examined the phylogeographical relationships of I. scapularis to determine the possible geographical origins of the different tick populations in southern Canada.
A total of 582 I. scapularis were collected between 2000 and 2011 (Table 1 and Additional file 1: Table S1). Of these, 70 were adventitious (i.e., ticks not associated with known resident populations) and were collected from hosts or the environment in different Canadian provinces, while 512 ticks were collected by drag sampling  at nine localities in Canada and the United States, each of which represented an established population of I. scapularis (Figure 1). For some data analyses, the collection localities of all ticks were grouped into one of two geographical regions. The “western” region comprised the Canadian Prairie Provinces (Alberta, Saskatchewan and Manitoba) and Minnesota in the Midwest of the United States. The “eastern” region included the Central Provinces (Ontario and Quebec) and Atlantic Provinces (Newfoundland, Nova Scotia, Prince Edward Island and New Brunswick) of Canada, and Rhode Island in the Northeast of the United States.
DNA extraction and PCR amplification
The total genomic (g) DNA of each tick was extracted using the QIAamp DNA Mini Kit™ or the DNeasy Blood & Tissue Kit™ (Qiagen), as described previously ,. A region (~400 bp) spanning Domains IV and V of the mt 16S rRNA gene was amplified by PCR from the total gDNA of each tick using the primers 16S − 1 (5’-CTGCTCAATGATTTTTTAAATTGCTGTGG -3’) and 16S + 1 (5’-CCGGTCTGAACTCAGATCAAGT-3’) . PCRs were performed in reaction mixtures (25 μl or 50 μl) containing PCR buffer with KCl, 1.75 mM MgCl2, 200 μM of each dNTP, 25 ρmol of each primer, 0.5-1.25 U of Taq polymerase, and 1–2 μl of gDNA template. The cycling conditions used were 96°C for 5 min, then 30 cycles of 94°C for 30 s, 52°C for 30 s and 72°C for 30 s, and a final extension at 72°C for 5 min. All PCR products were examined on 1.5% agarose-TBE gels to verify that each amplicon represented a single band of ~450 bp.
Single-strand conformation polymorphism (SSCP) and DNA sequencing
All amplicons were subjected to single-strand conformation polymorphism (SSCP) analysis, a mutation-scanning technique that is highly effective in displaying genetic variation among amplicons (150–450 bp) that differ in DNA sequence by one or more nucleotides . The SSCP methodology used followed that described previously , except that 0–4.5 μl of DNase-free water and 5 μl of loading buffer (Gel Tracking Dye™, Promega) were added to 0.5-5 μl of each amplicon. Where possible, multiple amplicons of each SSCP profile type were subjected to automated sequencing using primers 16S − 1 and 16S + 1 in separate reactions. Amplicons were purified prior to sequencing. This was achieved by adding 1 μl of a mixture containing 3 U of exonuclease I, 0.15 U of shrimp alkaline phosphatase, and 0.7 μl of 1X PCR buffer to 10 μl of each amplicon, and incubating the samples at 37°C for 15 min. Subsequently, increasing the temperature to 80°C for 15 min inactivated the enzymes. The sequences of the haplotypes have been deposited in GenBank™ under the accession numbers HG916768-HG916804. The numerical system for haplotype designation used herein follows that of Krakowetz et al. . BLAST searches (GenBank) were performed on the sequence data obtained to determine if the haplotypes detected in the present study were identical to those in other studies, but where different haplotype designations are used (see Additional file 2: Table S2).
Sequence alignment, secondary structure and phylogenetic analyses
Sequences were aligned manually, but then modified according to the predicted secondary structure of Domains IV and V of the mt 16S rRNA gene that was constructed for I. scapularis based on the models of other organisms . Phylogenetic analyses using the neighbor-joining (NJ) method were carried out using PAUP v4.0b10 . The DNA sequences of the mt 16S rRNA gene of I. pacificus (GenBank accession no. AF309008) and several haplotypes of the Southern and American clades of I. scapularis (see Supplemental Table S2 for accession nos.) , were included in the analyses. A bootstrap analysis (1000 replicates) was used to determine the relative support for groups in the NJ tree.
Population genetics and phylogeographical analyses
The haplotype (h) and nucleotide (π) diversities of I. scapularis within established populations were determined using Arlequin . Only data from the established populations of I. scapularis were included in the analyses. Tests for selective neutrality, Tajima’s D and Fu’s F S , were also performed using Arlequin. Under the neutral model, D and F S values should be approximately zero. Significantly negative D (p <0.05) and F S (p <0.02) values are indicative of populations undergoing expansion, whereas significantly positive values are characteristic of populations undergoing bottlenecks . Arlequin was also used to conduct a Chakraborty’s test , which determines if there were significantly more haplotypes in a population than expected under neutrality, and to calculate a measure (pairwise F ST ) of genetic differentiation between each pair of populations. The significance of departures of F ST values from zero was tested using 1000 permutations. A hierarchical analysis of molecular variance (AMOVA) was conducted using Arlequin to determine if there was genetic structuring within and among populations in different geographical regions. For this analysis, the nine established populations of I. scapularis were divided into four groups based on the province (Canada) or state (United States) in which they were located. A Mantel test was also conducted using Arlequin (1000 permutations) to determine if there was a correlation between genetic (F ST ) and geographical (km) distances among populations.
Rarefaction curves were generated using EstimateS  to estimate the total number of haplotypes that can be expected in a sample (i.e., based on the asymptote of the curve), and the extent to which the majority of haplotypes have been sampled. These analyses were performed for the populations in the western and eastern geographical regions, and the pooled population data of I. scapularis. EstimateS (1000 runs) was also used to determine Chao 2 values, estimators of the expected number of haplotypes in a sample .
A minimum spanning network tree depicting the relationships of the haplotypes was produced using TCS version 1.21 . This analysis also included other haplotypes of the American clade from Canada , the Northeast (e.g., Pennsylvania, Connecticut, New York, New Jersey, Massachusetts, Maryland and Rhode Island ,), and Midwest (e.g., Illinois and Wisconsin ) of the United States (see Additional file 2: Table S2).
Fifty-two different SSCP profiles were detected among the 582 amplicons (Table 2). Amplicons with the same SSCP profile had identical DNA sequences, while those that differed in SSCP profile also differed in DNA sequence by one or more nucleotides. The DNA sequences of the 52 haplotypes varied in length from 404–407 bp and differed from one another by 1–5 bp when aligned over 408 nucleotide positions (Table 2). Genetic variation among haplotypes was detected at 42 (10.3%) positions in the sequence alignment. These mutational differences consisted of 25 transitions, 12 transversions, four indels, and one multiple nucleotide change. There was approximately a 2:1 ratio of purine:pyrimidine transitional changes (18 and seven, respectively).
Sequence variation in relation to the secondary structure
Twenty-six (61.9%) of the 42 mutational changes occurred at unpaired sites (e.g., end loops and internal loops) in the predicted secondary structure of the mt 16S rRNA gene, while another nine mutational changes represented partial-compensatory base-pair changes that maintained the secondary structure (Figure 2). A majority of the mutational changes occurred within the hypervariable region (alignment positions 154 to 279; see Figure 2). This region comprised 126 (30.9%) of the nucleotides in the 3’ terminal end of the mt 16S rRNA gene, but contained 25 (59.5%) of the variable nucleotide positions among haplotypes. Thus, the proportion of variable positions in the hypervariable region was 19.8%. In contrast, there was a significantly (χ21 = 17.99, P <0.001) lower proportion of variable positions (6.0%; 17 of 282 positions) in the 5’ and 3’ regions flanking the hypervariable region, which represented 69.1% of the 408 total nucleotides at the 3’ terminal end of the gene. These two flanking regions contained 17 (40.5%) of the 42 variable nucleotide positions among sequence types.
The tree produced from the NJ analysis of the sequence data (Additional file 3: Figure S1) separated the 16S haplotypes of I. scapularis into two major clades; the Southern clade and American clade. This was supported by the results of the bootstrap analysis. However, there was little statistical support for the different groups within the American clade. All 582 I. scapularis characterized in the present study belonged to the American clade.
Population genetic analyses
The haplotype diversities and nucleotide diversities of I. scapularis within the nine established populations ranged from 0.5994 to 0.7856, and 0.00191 to 0.00357, respectively (Table 3). For six populations, the 7–14 haplotypes detected within each population was not significantly different from the expected number based on Chakraborty’s test, whereas significantly more haplotypes were detected than expected in all three populations from Minnesota (Table 3). The results of the Tajima’s test for neutrality indicated that there were significant negative departures from zero for the three tick populations in Minnesota and the population in Pembina Valley Provincial Park, Manitoba (Table 3). Similarly, the F S statistics of Fu revealed that the F S values for five populations (three in Minnesota and two in Manitoba) differed significantly from zero. A single haplotype (Is-1) was found in all nine I. scapularis populations and comprised 34-63% of the ticks within each population (Figure 1). The second most common haplotype differed among populations and represented between 7-32% of the ticks within each population. Comparison of the F ST values revealed significant differences between most pairs of populations except between some populations in Manitoba and Minnesota (Table 4). The results of the AMOVA test (Table 5) also indicated strong genetic structure among populations both within and among different geographical regions. Most of the variance (94.3%) occurred within populations. Nonetheless, there were many shared haplotypes among populations in different geographical regions (Figure 1). The results of the Mantel Test (Figure 3), which compared pair-wise F ST values as a factor of geographical distance between tick populations, showed that there was a significant association between genetic and geographical distances (b = 0.000058, r2 = 0.299, P = 0.002).
A total of 30 haplotypes were detected in five populations in the western geographical region, while 22 haplotypes were detected in the four populations in the eastern geographical region. The Chao 2 estimates of haplotype richness were higher for the western populations than in the eastern populations (41 and 34 haplotypes, respectively; Figure 4). However, the results of the rarefaction analyses showed that, although the curves for both the western and eastern populations of I. scapularis did not converge on an asymptote, there was no significant difference in haplotype diversity between populations of the two geographical areas, as there was overlap in the 95% confidence intervals of the two curves (Figure 4). Similarly, when the data for all populations were pooled, the rarefaction curve (see Additional file 4: Figure S2) did not reach an asymptote or approach the Chao 2 estimate of the haplotype diversity (i.e., 82 haplotypes). Therefore, only 45 (55%) of the expected total number of haplotypes were detected in the nine established populations of I. scapularis.
Fifteen (29%) of the 52 haplotypes detected among the adventitious ticks and individuals from established populations in the present study were each represented by a single tick, while four haplotypes were detected in western and eastern Canada, and in the Midwest and Northeast of the United States (Figure 5 and Additional file 5: Table S3). Thirty-two haplotypes were detected among ticks from the United States; however, only four (13%) were present in populations in the Midwest and Northeast. Ticks collected from western Canada had the highest similarity, based on the proportion of shared haplotypes, with the ticks from the Midwest of the United States (13 of 33; 39%) rather than those in eastern Canada (7 of 36; 19%) or the Northeast of the United States (5 of 27; 19%). The proportion of shared haplotypes between ticks from eastern Canada and the Northeast of the United States (7 of 26; 27%) was greater than that between ticks from the Midwest of the United States (6 of 41; 15%) or western Canada (19%).
The minimum spanning network tree depicting the relationships among haplotypes of the American clade is shown in Figure 6. Most (82%) haplotypes were only found in either the western or eastern geographical region (27 and 24 haplotypes, respectively). Haplotype Is-1, the most common haplotype in both geographical regions (see Additional file 5: Table S3), represented the central haplotype of the star-shaped network tree. All other haplotypes differed from the central haplotype by 1–6 bp. Thirty-seven haplotypes differed in sequence from the central haplotype by a single nucleotide mutation. These included five of the seven most common haplotypes detected in the present study (Is-2, Is-4, Is-13, Is-15, and Is-63; Additional file 5: Table S3). Each of these haplotypes represented a link (secondary node) from the central haplotype to other haplotypes (Figure 6). Six of the eight most common haplotypes were detected in I. scapularis populations in both the western and eastern geographical regions. Another common haplotype present in eastern populations (Is-7) represented a tertiary node in the haplotype network to six other haplotypes that were also only detected in eastern tick populations.
Fifty-two haplotypes of the mt 16S rRNA gene were detected among the 582 I. scapularis individuals collected from southern Canada, and the Midwest (Minnesota) and Northeast (Rhode Island) of the United States. All 52 haplotypes belonged to the American clade, as defined by Qiu et al. . The lack of I. scapularis individuals of the Southern clade in northern United States and southern Canada is consistent with the findings of other studies -,. The total number of haplotypes detected in the present study was greater than the 7–29 haplotypes detected in other studies of the American clade ,-), and included 30 haplotypes not reported previously. The difference in number of the haplotypes detected among studies may be a consequence of differences in sample sizes. For example, the number of haplotypes among I. scapularis individuals collected from Camp Riley in Morrison County, Minnesota (14 haplotypes), Trustom Pond and Hazard Island in Washington County, Rhode Island (11 haplotypes), and Long Point Provincial Park, Ontario (13 haplotypes) were two to three times greater than that reported previously in these areas (3, 4 and 6 haplotypes, respectively) ,,; however, at least twice as many ticks were sampled from each of these localities in the present study. Nonetheless, the results of the present study indicate that genetic diversity in I. scapularis from northern United States and southern Canada, based on DNA sequences of the mt 16S rRNA gene, is considerably greater than previously thought. Furthermore, there may be a large number of undetected haplotypes within the sampled areas because the Chao 2 estimate of the total number of expected haplotypes (i.e., 82) was greater than the 52 haplotypes detected. This is likely given that 15 (29%) of the haplotypes can be considered as rare because only one tick of each of these haplotypes was detected in the present study.
Nucleotide diversities within established populations of I. scapularis were low (0.002-0.004), while haplotype diversities were relatively high (0.60-0.79) compared to those of some other species of Ixodes in North America -. For example, only 1–4 haplotypes of the mt 16S rRNA gene have been detected within populations of Ixodes angustus, Ixodes kingi and Ixodes sculptus- compared to the 7–14 haplotypes among individuals in populations of I. scapularis. Biological differences among these tick species are likely explanations for the differences in haplotype number. For instance, I. angustus, I. kingi and I. sculptus parasitize primarily small mammals (e.g., mice, shrews, voles, ground squirrels and/or pocket gophers -, which provide limited dispersal distances for ticks. In contrast, I. scapularis parasitizes a wider diversity of animals including passerine birds ,, hosts that are known to carry I. scapularis larvae and nymphs over large distances ,. Therefore, increased dispersal distance provides a greater opportunity for transfer of ticks representing different maternal lineages (haplotypes) from one population to another.
Variation in the DNA sequences among I. scapularis individuals was also examined in relation to the predicted secondary structure of Domains IV and V of the mt 16S rRNA gene. Over half (57%) of the nucleotide variation occurred within the “hypervariable” region of Domain V. Mutational alterations in DNA sequence, both within and among species of arthropod, have been shown previously to be more frequent within this part of the gene than in the flanking regions ,-. This suggests that there are fewer structural constraints for mutational changes within the hypervariable region than in other parts of Domain V . DNA sequence variation in I. scapularis at unpaired sites (62%) was also greater than that at sites involved in base pairing (38%) in the secondary structure. This difference may also be associated with reduced structural constraints for mutational alterations at nucleotide positions not involved in base pairing. In addition, substitutions at nucleotide positions involved in base pairing were more often partial compensatory base changes, rather than non-compensatory changes, hence maintaining the integrity of the secondary structure of the mt 16S rRNA gene. This is in agreement with the general patterns of mutational change in rRNA genes .
The number of nucleotide differences in DNA sequence among the 52 haplotypes ranged from 1–5 bp; however, this difference increased to 9 bp when an additional 10 haplotypes of the American clade -, were added to the analyses. Seven haplotypes (Is-1, Is-2, Is-4, Is-7, Is-13, Is-15 and Is-63) comprised 78% of the ticks collected in the present study. Of these, haplotype Is-1 was the most common (49%), which was consistent with the findings of other studies on I. scapularis conducted in southern Canada , and the Midwest and Northeast of the United States ,. Haplotype Is-1 also represented the central haplotype of the minimum spanning network. Some of the other common haplotypes represented secondary or tertiary nodes in the network. This suggests that most haplotypes were derived from the central haplotype or one of the secondary haplotypes, as a consequence of a single mutational change in the DNA sequence of the mt 16S rRNA gene. Furthermore, the star-shaped pattern of the network tree for I. scapularis is indicative of rapidly expanding populations ,, which is in agreement with the relatively recent establishment of populations of I. scapularis in the Midwest of the United States , and southern Canada ,,. Furthermore, for the populations in Manitoba and Minnesota, there were significant negative departures from zero for both the Tajima’s D and Fu’s Fs tests suggesting population expansion at these localities.
Range expansion and the establishment of geographically isolated populations of I. scapularis into southern Canada have been attributed to the transportation of large numbers of larvae and nymphs from the United States by migratory passerine birds ,. However, the geographical origins of these populations are unknown. Krakowetz et al.  reported differences in the genetic structure of I. scapularis between a population in southeastern Manitoba and several populations in southern Ontario and Nova Scotia. Subsequently, Mechai et al.  reported differences in the frequency of haplotypes of the mt cytochrome C oxidase subunit 1 gene (cox1) among blacklegged ticks from different geographical regions in Canada. It was proposed that I. scapularis populations in different geographical regions of southern Canada may be derived from populations in different regions of the United States and associated with the different routes (flyways) taken by passerines during their spring migration . Scott et al.  also suggested that there was an association between the presence of I. scapularis and other introduced tick species in different regions of Canada and the different flyways of migratory passerines. Thus, blacklegged ticks introduced into the Prairie Provinces of Canada (Manitoba, Saskatchewan and Alberta) may have originated from endemic populations in the Midwest of the United States (Minnesota and Wisconsin), while those introduced into the Central Provinces of Canada (Ontario and Quebec) may have originated from endemic populations primarily in the Northeast of the United States (Connecticut, Pennsylvania, New York, Massachusetts, Rhode Island, Maine and New Hampshire), but also from resident populations in parts of the Midwest ,. Northeastern United States was also suggested as the most likely origin of blacklegged ticks introduced into the Atlantic Provinces (Nova Scotia, New Brunswick, Newfoundland and Prince Edward Island) ,. The latter is supported by the results of the present study. For example, 38% of ticks from an established population in Lunenburg (Nova Scotia) were haplotype Is-12, a haplotype that had not been previously reported from other regions of Canada or in the United States ; however, this haplotype was detected in Rhode Island (present study). Furthermore, nine (75%) of the 12 haplotypes found in the Atlantic Provinces were also detected in Rhode Island, seven of which have also been found in other parts of the Northeast of the United States .
The results of the AMOVA test also revealed statistically significant genetic structuring of I. scapularis populations both within and among different geographical regions; however, the presence of several shared haplotypes among populations supports the hypothesis of gene flow among populations. Nonetheless, 26 (79%) of the 33 haplotypes found in the western region (Prairie Provinces and Minnesota) were not found in the eastern region (Central and Atlantic Provinces, and Rhode Island), while 19 of the 26 (73%) haplotypes in the east were not found in Minnesota or the Prairie Provinces. In addition, only 8 (32%) of the 25 haplotypes found in the three populations in Minnesota have been reported previously from the Northeast of the United States ,, which includes four haplotypes present in Rhode Island. There was also a significant positive correlation between the geographical (km) distances among populations and the magnitude of genetic differences (F ST values). Statistical analyses of the F ST data also showed significant differences in the genetic structure between some populations of I. scapularis. For example, there was a significant difference in the population genetic structure of I. scapularis from Hazard Island and Trustom Pond, two localities in South Kingstown (Rhode Island) separated by a distance of only 5 km. The reason why the genetic structure of these two tick populations differs is unclear. Although fewer haplotypes were detected at Trustom Pond, six of the seven haplotypes in this population were also detected in the population at Hazard Island. In addition, the results of the Chakraborty’s test revealed no significant difference in the number of observed and expected haplotypes for either population. In contrast, there were no significant differences between the populations in Itasca State Park (Minnesota), Pembina Valley Provincial Park and Stanley Trail (Manitoba). There were also no differences between the populations in Itasca State Park and St. Croix State Park, or between those in Camp Ripley and St. Croix State Park (Minnesota). A comparison of the haplotypes among adventitious ticks and those in established populations found in the Prairie Provinces with those in other geographical areas (Figure 6) revealed a greater similarity to I. scapularis in Minnesota than to those in the Central and Atlantic Provinces of Canada or Rhode Island based on the proportion of shared haplotypes. These results provide some support for the hypothesis that the I. scapularis populations in southern Manitoba are derived from established populations in the Midwest of the United States. However, a large proportion (58%) of the I. scapularis individuals collected in all three regions of southern Canada (Prairie, Central and Atlantic Provinces) were of a haplotype (Is-1, Is-2, Is-4 and Is-6) that also occurs in both the Midwest and Northeast of the United States ,. Therefore, other genetic markers (e.g., 12S rRNA gene  or cox1; ), in addition to the mt 16S rRNA gene, are needed to determine the geographical origin of I. scapularis introductions into southern Canada.
In conclusion, genetic variation within I. scapularis was greater than previously demonstrated based on the DNA sequence analyses of the mt 16S rRNA gene. Furthermore, a large number of rare haplotypes may still remain undetected. The results also indicated significant differences in genetic diversity both within and among populations from different geographical regions. There was also a significant positive relationship between the genetic differences between populations and the geographical distances that separated them. There was some evidence to support the hypothesis that I. scapularis in the Prairie Provinces of Canada are derived from individuals introduced from the Midwest of the United States, while those in the Atlantic and Central Provinces are derived from individuals that originated in the Northeast of the United States. However, the geographical origins of a large proportion of I. scapularis found in the different areas of southern Canada could not be inferred because they were of a haplotype that occurs in both the Midwest and Northeast of the United States. Therefore, additional studies are needed to explore other genetic markers that may be useful for understanding the trajectories of spread of I. scapularis and its pathogens on a finer scale.
St. Croix State P
- h :
Itasca State Park
Long Point Provincial Park
Polymerase chain reaction
Point Pelee National Park
Pembina Valley Provincial Park
Single-strand conformation polymorphism
Piesman J, Eisen L: Prevention of tick-borne diseases. Ann Rev Entomol. 2008, 53: 323-343. 10.1146/annurev.ento.53.103106.093429.
Parola P, Raoult D: Ticks and tickborne bacterial diseases in humans: an emerging infectious threat. Clin Infect Dis. 2001, 32: 897-928. 10.1086/319347.
Dahlgren FS, Mandel EJ, Krebs JW, Massung RF, McQuiston JH: Increasing incidence of Ehrlichia chaffeensis and Anaplasma phagocytophilum in the United States, 2000–2007. Am J Trop Med Hyg. 2011, 85: 124-131. 10.4269/ajtmh.2011.10-0613.
Jobe DA, Nelson JA, Adam MD, Martin SA: Lyme disease in urban areas, Chicago. Emerg Infect Dis. 2007, 13: 1799-1800. 10.3201/eid1311.070801.
Hamer S, Tsao JI, Walker ED, Hickling GJ: Invasion of the Lyme disease vector Ixodes scapularis: implications for Borrelia burgdorferi endemicity. EcoHealth. 2010, 7: 47-63. 10.1007/s10393-010-0287-0.
Ogden NH, Trudel L, Artsob H, Barker IK, Beauchamp G, Charron DF, Drebot MA, Galloway TD, O’Handley R, Thompson RA, Lindsay LR:Ixodes scapularis ticks collected by passive surveillance in Canada: analysis of geographic distribution and infection with Lyme borreliosis agent Borrelia burgdorferi. J Med Entomol. 2006, 43: 600-609. 10.1603/0022-2585(2006)43[600:ISTCBP]2.0.CO;2.
Ogden NH, Lindsay LR, Hanincova K, Barker IK, Bigras-Poulin M, Charron DF, Heagy A, Francis CM, O’Callaghan CJ, Schwartz I, Thompson RA: Role of migratory birds in introduction and range expansion of Ixodes scapularis ticks and of Borrelia burgdorferi and Anaplasma phagocytophilum in Canada. Appl Environ Microbiol. 2008, 74: 1780-1790. 10.1128/AEM.01982-07.
Ogden NH, Lindsay LR, Morshed MG, Sockett PN, Artsob H: The rising challenge of Lyme borreliosis in Canada. Can Commun Dis Rep. 2008, 34: 1-19.
Bouchard C, Beauchamp G, Nguon S, Trudel L, Milord F, Lindsay LR, Bélanger D, Ogden NH: Associations between Ixodes scapularis ticks and small mammal hosts in a newly endemic zone in southeastern Canada: implications for Borrelia burgdorferi transmission. Ticks Tick-borne Dis. 2011, 2: 183-190. 10.1016/j.ttbdis.2011.03.005.
Koffi JK, Leighton PA, Pelcat Y, Trudel L, Lindsay LR, Milord F, Ogden NH: Passive surveillance for I. scapularis ticks: enhanced analysis for early detection of emerging Lyme disease risk. J Med Entomol. 2012, 49: 400-409. 10.1603/ME11210.
Burgdorfer W, Barbour AG, Hayes SF, Benach JL, Grunwaldt E, Davis JP: Lyme disease-a tick-borne spirochetosis?. Science. 1982, 216: 1317-1319. 10.1126/science.7043737.
Chen SM, Dumler JS, Bakken JS, Walker DH: Identification of a granulocytotropic Ehrlichia species as the etiologic agent of human disease. J Clin Microbiol. 1994, 32: 589-595.
Spielman A: Human babesiosis on Nantucket Island - transmission by nymphal Ixodes ticks. Am J Trop Med Hyg. 1976, 25: 784-787.
Costero A, Grayson MA: Experimental transmission of Powassan virus (Flaviviridae) by Ixodes scapularis ticks (Acari: Ixodidae). Am J Trop Med Hyg. 1996, 55: 536-546.
Steiner FE, Pinger RR, Vann CN, Grindle N, Civitello D, Clay K, Fuqua C: Infection and co-infection rates of Anaplasma phagocytophilum variants, Babesia spp., Borrelia burgdorferi, and the rickettsial endosymbiont in Ixodes scapularis (Acari: Ixodidae) from sites in Indiana, Maine, Pennsylvania, and Wisconsin. J Med Entomol. 2008, 45: 289-297. 10.1603/0022-2585(2008)45[289:IACROA]2.0.CO;2.
Gatewood AG, Liebman KA, Vourc’h G, Bunikis J, Hamer SA, Cortinas R, Melton F, Cislo P, Kitron U, Tsao J, Barbour AG, Fish D, Diuk-Wasser MA: Climate and tick seasonality are predictors of Borrelia burgdorferi genotype distribution. Appl Environ Microbiol. 2009, 75: 2476-2483. 10.1128/AEM.02633-08.
Brisson D, Vandermause MF, Meece JK, Reed KD, Dykhuizen DE: Evolution of northeastern and midwestern Borrelia burgdorferi, United States. Emerg Infect Dis. 2010, 16: 911-917. 10.3201/eid1606.090329.
Humphrey PT, Caporale DA, Brisson D: Uncoordinated phylogeography of Borrelia burgdorferi and its tick vector, Ixodes scapularis. Evolution. 2010, 64: 2653-2663. 10.1111/j.1558-5646.2010.01001.x.
Kurtenbach K, Hanincová K, Tsao JI, Margos G, Fish D, Ogden NH: Fundamental processes in the evolutionary ecology of Lyme borreliosis. Nat Rev Microbiol. 2006, 4: 660-669. 10.1038/nrmicro1475.
Keirans JE, Hutcheson HJ, Durden LA, Klompen JSH:Ixodes (Ixodes) scapularis (Acari: Ixodidae): redescription of all active stages, distribution, hosts, geographical variation, and medical and veterinary importance. J Med Entomol. 1996, 33: 297-318.
Dennis DT, Nekomoto TS, Victor JC, Paul WS, Piesman J: Reported distribution of Ixodes scapularis and Ixodes pacificus (Acari: Ixodidae) in the United States. J Med Entomol. 1998, 35: 629-638.
Diuk-Wasser MA, Hoen AG, Cislo P, Brinkerhoff R, Hamer SA, Rowland M, Corinas R, Vourc’h G, Melton F, Hickling GJ, Tsao JI, Bunikis J, Barbour AG, Kitro U, Piesman J, Fish D: Human risk of infection with Borrelia burgdorferi, the Lyme Disease agent in eastern United Sates. Am J Trop Med. 2012, 86: 320-327. 10.4269/ajtmh.2012.11-0395.
Hoen AG, Margos G, Bent SJ, Diuk-Wasser MA, Barbout A, Kurtenbach K, Fish D: Phylogeography of Borrelia burgdorferi in the eastern United States reflects multiple independent Lyme disease emergence events. Proc Natl Acad Sci. 2009, 106: 15013-15018. 10.1073/pnas.0903810106.
Watson TG, Anderson RC:Ixodes scapularis Say on white-tailed deer (Odocoileus virginianus) from Long Point, Ontario. J Wildl Dis. 1976, 12: 66-71. 10.7589/0090-3558-12.1.66.
Barker IK, Lindsay LR: Lyme borreliosis in Ontario: determining the risks. Can Med Assoc J. 2000, 162: 1573-1574.
Morshed MG, Scott JD, Fernando K, Mann RB, Durden LA: Lyme disease spirochete Borrelia burgdorferi endemic at epicentre in Rondeau Provincial Park, Ontario. J Med Entomol. 2003, 40: 91-94. 10.1603/0022-2585-40.1.91.
Ogden NH, Bigras-Poulin M, O’Callaghan CJ, Barker IK, Lindsay LR, Maarouf A, Smoyer-Tomic KE, Waltner-Toews D, Charron D: A dynamic population model to investigate effects of climate on geographic range and seasonality of the tick Ixodes scapularis. Int J Parasitol. 2005, 35: 375-389. 10.1016/j.ijpara.2004.12.013.
Ogden NH, Lindsay LR, Morshed M, Sockett PN, Artsob H: The emergence of Lyme disease in Canada. Can Med Assoc J. 2009, 180: 1221-1224. 10.1503/cmaj.080148.
Scott JD, Anderson JF, Durden LA: Widespread dispersal of Borrelia burgdorferi-infected ticks collected from songbirds across Canada. J Parasitol. 2012, 98: 49-59. 10.1645/GE-2874.1.
Rich SM, Caporale DA, Telford SR, Kocher TD, Hartl DL, Spielman A: Distribution of the Ixodes ricinus-like ticks of eastern North America. Proc Natl Acad Sci U S A. 1995, 92: 6284-6288. 10.1073/pnas.92.14.6284.
Norris DE, Klompen JSH, Keirans JE, Black WC: Population genetics of Ixodes scapularis (Acari: Ixodidae) based on mitochondrial 16S and 12S genes. J Med Entomol. 1996, 33: 78-89.
Qiu W-G, Dykhuizen DE, Acosta MS, Luft BJ: Geographic uniformity of the Lyme disease spirochete (Borrelia burgdorferi) and its shared history with tick vector (Ixodes scapularis) in the northeastern United States. Genetics. 2002, 160: 833-849.
Trout RT, Steelman CD, Szalanski AL: Population genetics and phylogeography of Ixodes scapularis from canines and deer in Arkansas. Southwest Entomol. 2009, 34: 273-287. 10.3958/059.034.0308.
Krakowetz CN, Lindsay LR, Chilton NB: Genetic diversity in Ixodes scapularis (Acari: Ixodidae) from six established populations in Canada. Ticks Tick-borne Dis. 2011, 2: 143-150. 10.1016/j.ttbdis.2011.05.003.
Lindsay LR, Mathison SW, Barker IK, McEwen SA, Surgeoner GA: Abundance of Ixodes scapularis (Acari: Ixodidae) larvae and nymphs in relation to host density and habitat on Long Point, Ontario. J Med Entomol. 1999, 36: 243-254.
Dergousoff SJ, Chilton NB: Differentiation of three species of ixodid tick, Dermacentor andersoni, D. variabilis and D. albipictus, by PCR-based approaches using markers in ribosomal DNA. Mol Cell Probes. 2007, 21: 343-348. 10.1016/j.mcp.2007.04.003.
Gasser RB, Hu M, Chilton NB, Campbell BE, Jex AJ, Otranto D, Cafarchia C, Beveridge I, Zhu X: Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation. Nat Protoc. 2006, 1: 3121-3128. 10.1038/nprot.2006.485.
Gutell RR, Gray MW, Schnare MN: A compilation of large subunit (23S and 23S-like) ribosomal RNA structures: 1993. Nucleic Acids Res. 1993, 21: 3055-3074. 10.1093/nar/21.13.3055.
Swofford DL: PAUP*: phylogenetic analysis using parsimony (and other methods) 4.0 Beta. 4th edition. Sunderland: Sinauer Associates; 2002.
Excoffier L, Lischer HEL: Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.
Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphisms. Genetics. 1989, 123: 585-595.
Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.
Aris-Brosou S, Excoffier L: The impact of population expansion and mutation rate heterogeneity on DNA sequence polymorphism. Mol Biol Evol. 1996, 13: 494-504. 10.1093/oxfordjournals.molbev.a025610.
Chakraborty R: Mitochondrial DNA polymorphism reveals hidden heterogeneity within some Asian populations. Am J Hum Genet. 1990, 47: 87-94.
Colwell RK: EstimateS: statistical estimation of species richness and shared species from samples. 2013
Colwell RK, Coddington JA: Estimating terrestrial biodiversity through extrapolation. Phil Trans R Soc Lond B Biol Sci. 1994, 345: 101-118. 10.1098/rstb.1994.0091.
Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1659. 10.1046/j.1365-294x.2000.01020.x.
Anstead CA, Chilton NB: Comparison of the DNA sequences and secondary structure of the mitochondrial 16S rRNA gene of Ixodes kingi, Ixodes sculptus and Ixodes angustus. Mol Cell Probes. 2014, 28: 155-162. 10.1016/j.mcp.2014.02.001.
Anstead CA, Chilton NB: Ticks feeding on northern pocket gophers (Thomomys talpoides) in central Saskatchewan and the unexpected detection of Ixodes scapularis larvae. J Vector Ecol. 2011, 36: 355-360. 10.1111/j.1948-7134.2011.00176.x.
Anstead CA, Hwang YT, Chilton NB: Ticks (Acari: Ixodidae) on small mammals in Kootenay National Park, British Columbia, Canada. J Med Entomol. 2013, 50: 1208-1214. 10.1603/ME13067.
Anstead CA, Wallace SB, Chilton NB: Mutation scanning-based identification of larval and nymphal ticks (Acari: Ixodidae) from Richardson’s ground squirrels (Spermophilus richardsonii). Mol Cell Probes. 2014, 28: 6-9. 10.1016/j.mcp.2013.09.004.
Allan SA: Ticks (Class Arachnida: Order Acarina). Parasitic Diseases of Wild Mammals. Edited by: Samuel WM, Pybus MJ, Kocan AA. 2001, Iowa State University Press, Iowa, 72-106. 10.1002/9780470377000.ch4.
Smith SD, Bond JE: An analysis of the secondary structure of the mitochondrial large subunit rRNA gene (16S) in spiders and its implication for phylogenetic reconstruction. J Arachnol. 2003, 31: 44-54. 10.1636/0161-8202(2003)031[0044:AAOTSS]2.0.CO;2.
Misof B, Anderson CL, Buckley TR, Erpenbeck D, Rickert A, Misof K: An empirical analysis of mt 16S rRNA covarion-like evolution in insects: site-specific rate variation is clustered and frequently detected. J Mol Evol. 2002, 55: 460-469. 10.1007/s00239-002-2341-1.
Machado EG, Dennebouy N, Suarez MO, Mounolou J-C, Monnerot M: Mitochondrial 16S-rRNA gene of two species of shrimps: sequence variability and secondary structure. Crustaceana. 1993, 65: 279-286. 10.1163/156854093X00711.
Whitfield J, Cameron SA: Hierarchical analysis of variation in the mitochondrial 16S rRNA gene among Hymenoptera. Mol Biol Evol. 1998, 15: 1728-1742. 10.1093/oxfordjournals.molbev.a025899.
Buckley TR, Simon C, Flook PK, Misof B: Secondary structure and conserved motifs of the frequently sequenced Domains IV and V of the insect mitochondrial large subunit rRNA gene. Insect Mol Biol. 2000, 9: 565-580. 10.1046/j.1365-2583.2000.00220.x.
Sorokina SY, Mugue NS, Andrianov BV, Mitrofanov VG: Variation of 3’-terminal fragment of 16S rRNA gene in closely related species of Drosophila virilis group. Russ J Genet. 2005, 41: 853-858. 10.1007/s11177-005-0172-z.
Hillis DM, Dixon MT: Molecular evolution and phylogenetic inference. Quart Rev Biol. 1991, 66: 411-452. 10.1086/417338.
Lee X, Hardy K, Johnson DH, Paskewitz SM: Hunter-killed deer surveillance to assess changes in the prevalence and distribution of Ixodes scapularis (Acari: Ixodidae) in Wisconsin. J Med Entomol. 2013, 50: 632-639. 10.1603/ME12234.
Mechai S, Feil EJ, Gariepy TD, Gregory TR, Lindsay LR, Millien V, Ogden NH: Investigation of the population structure of the tick vector of Lyme disease Ixodes scapularis (Acari: Ixodidae) in Canada using mitochondrial cytochrome C oxidase subunit I gene sequences. J Med Entomol. 2013, 50: 560-570. 10.1603/ME12178.
We thank Thomas Mather (Rhode Island), David Nietzel and Melissa Kemperman (Minnesota) for kindly providing ticks from the United States, and Shaun Dergousoff for technical advice. This research was supported by grants (to NBC) from the Natural Sciences and Engineering Research Council of Canada (NSERC) and the Canadian Foundation for Innovation. Scholarship support (to CNK) by NSERC (CGSM & Vanier Canada Graduate Scholarships) and the Provincial Government of Saskatchewan in combination with the University of Saskatchewan (Saskatchewan Innovation and Opportunity Scholarships) are gratefully acknowledged.
The authors declare that they have no competing interests.
NBC conceived the project. LRL organized the collection of samples. CNK carried out laboratory work. CNK and NBC performed the data analyses. All authors interpreted the data, wrote the manuscript, and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1.: The number of nymphal or adult male and female I. scapularis collected in different years from 11 geographical regions in North America. (DOCX 72 KB)
Additional file 2: Table S2.: Comparison of the haplotype designations of I. scapularis for the mt 16S rRNA gene used in the minimum spanning network tree (see Figure 6) of the present study in relation to those (e.g., haplotype, specimen number or GenBank accession no.) used in previous studies. (DOCX 85 KB)
Additional file 3: Figure S1.: Neighbor-joining tree depicting the relationships of the 52 mt 16S rDNA haplotypes of I. scapularis detected in the present study. Also included are an additional 10 haplotypes from other studies of the American , and Southern clades ,. Numbers above branches indicate the bootstrap values (>70%). The scale bar represents the inferred substitutions per nucleotide site. Haplotypes identical to those of haplotypes A-M of Qiu et al.  are indicated by an *. (PPTX 73 KB)
Additional file 4: Figure S2.: Rarefaction curve (with 95% confidence intervals) of the haplotype diversity for the I. scapularis populations sampled in the present study. The number next to the curve indicates the total estimated number of haplotypes using the non-parametric Chao 2 estimator. See list of abbreviations for the complete names of localities of the tick populations. (PDF 28 KB)
Additional file 5: Table S3.: The number of I. scapularis of the different mt 16S rRNA gene haplotypes collected from different geographical regions. (DOCX 79 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Krakowetz, C.N., Lindsay, L.R. & Chilton, N.B. Genetic variation in the mitochondrial 16S ribosomal RNA gene of Ixodes scapularis (Acari: Ixodidae). Parasites Vectors 7, 530 (2014). https://doi.org/10.1186/s13071-014-0530-6
- Ixodes scapularis
- Genetic variation
- Mitochondrial 16S rRNA gene
- Secondary structure
- Population genetics