Characterization of cytochrome b diversity in Chinese domestic horses

Previous mitochondrial DNA (mtDNA) D‐loop and microsatellite studies have shown that Chinese horses have multiple maternal origins and high genetic diversity. To better characterize maternal genetic origins and diversity of Chinese domestic horses, we conducted a comprehensive analysis of 407 complete 1140 bp sequences of the horse mitochondrially encoded cytochrome b (CYTB) gene, including 323 horses from 13 Chinese indigenous breeds and 84 reference sequences from GenBank. A total of 114 haplotypes were identified, of which 73 appeared among the 13 Chinese horse breeds. The high mitochondrially encoded cytochrome b haplotypic diversity suggests multiple maternal origins in Chinese horses.     

Genetic diversity and population structure of Schizopygopsis younghusbandi Regan in the Yarlung Tsangpo River inferred from mitochondrial DNA sequence analysis

The genetic diversity and population structure of Schizopygopsis younghusbandi from six sites in the middle reach of the Yarlung Tsangpo River were examined based on mitochondrial DNA cytochrome b (Cyt b) gene and control region (D-loop) sequences. 1141 bp of complete Cyt b sequences and 737 bp of partial D-loop sequences for 153 individuals of S. younghusbandi were obtained by using PCR amplification and sequencing. The results showed that S. younghusbandi populations had high haplotype diversity and low nucleotide diversity, and the population genetic diversity of this species was at a low level. Analysis of molecular variance revealed that most genetic variation occurred within populations, and that genetic differentiation was at low or moderate levels. The network of haplotypes based on Cyt b gene showed that there were two groups within the examined populations. Neutrality tests and mismatch distribution analysis suggested that this species might have undergone a population expansion and the expansion time was estimated to be 0.25–0.46 Ma BP. All these results would be crucial to establish scientific strategies for effective conservation and sustainable exploitation of wild resources of S. younghusbandi.     

Population structure and biogeography of Hemiphaga pigeons (Aves: Columbidae) on islands in the New Zealand region

Aim The New Zealand avifauna includes lineages that lack close relatives elsewhere and have low diversity, characteristics sometimes ascribed to long geographic isolation. However, extinction at the population and species levels could yield the same pattern. A prominent example is the ecologically important pigeon genus Hemiphaga. In this study, we examined the population structure and phylogeography of Hemiphaga across islands in the region. Location New Zealand, Chatham Islands and Norfolk Island. Methods Mitochondrial DNA was sequenced for all species of the genus Hemiphaga. Sixty‐seven individuals from mainland New Zealand (Hemiphaga novaeseelandiae novaeseelandiae), six of the Chatham Islands sister species (Hemiphaga chathamensis), and three of the extinct Norfolk Island subspecies (Hemiphaga novaeseelandiae spadicea) were included in this study. Novel D‐loop and cytochrome b primers were designed to amplify DNA from museum samples. Additionally, five other mitochondrial genes were used to examine placement of the phylogenetic root. Results Analyses of mitochondrial DNA sequences revealed three Hemiphaga clades, consistent with the allopatric populations of recognized (sub)species on oceanic islands. Of the 23 D‐loop haplotypes among 67 New Zealand pigeons (Hemiphaga n. novaeseelandiae), 19 haplotypes were singletons and one haplotype was common and widespread. Population genetic diversity was shallow within and between New Zealand populations, indicating range expansion with high inter‐population exchange. Tentative rooting of the Hemiphaga clade with cyt b data indicates exchange between mainland New Zealand and the Chatham Islands prior to colonization of Norfolk Island. We found low genetic divergence between populations on New Zealand, the Chatham Islands and Norfolk Island, but deep phylogenetic divergence from the closest living relatives of Hemiphaga. Main conclusions The data are consistent with the hypothesis of population reduction during the Pleistocene and subsequent expansion from forest refugia. Observed mobility of Hemiphaga when feeding helps explain the shallow diversity among populations on islands separated by many hundreds of kilometres of ocean. Together with comparison of distribution patterns observed among birds of the New Zealand region, these data suggest that endemicity might represent not long occupancy of an area, but descent from geologically recent colonizations. We consider the role of lineage pruning in creating the impression of old endemicity.     

Mitochondrial DNA variation in natural populations of endangered Indian Feather-Back Fish, Chitala chitala

Genetic variation at mitochondrial cytochrome b (cyt b) and D-loop region reveals the evidence of population sub-structuring in Indian populations of highly endangered primitive feather-back fish Chitala chitala. Samples collected through commercial catches from eight riverine populations from different geographical locations of India were analyzed for cyt b region (307 bp) and D-loop region (636–716 bp). The sequences of the both the mitochondrial regions revealed high haplotype diversity and low nucleotide diversity. The patterns of genetic diversity, haplotypes networks clearly indicated two distinct mitochondrial lineages and mismatch distribution strongly suggest a historical influence on the genetic structure of C. chitala populations. The baseline information on genetic variation and the evidence of population sub-structuring generated from this study would be useful for planning effective strategies for conservation and rehabilitation of this highly endangered species.     

Maternal and paternal genetic diversity of ancient sheep in Estonia from the Late Bronze Age to the post‐medieval period and comparison with other regions in Eurasia

Sheep were among the first domesticated animals to appear in Estonia in the late Neolithic and became one of the most widespread livestock species in the region from the Late Bronze Age onwards. However, the origin and historical expansion of local sheep populations in Estonia remain poorly understood. Here, we analysed fragments of the hypervariable D‐loop of mitochondrial DNA (mtDNA; 213 bp) and the Y‐chromosome SRY gene (130 bp) extracted from 31 archaeological sheep bones dated from approximately 800 BC to 1700 AD. The ancient DNA data of sheep from Estonia were compared with ancient sheep from Finland as well as a set of contemporary sheep breeds from across Eurasia in order to place them in a wider phylogeographical context. The analysis shows that: (i) 24 successfully amplified and analysed mtDNA sequences of ancient sheep cluster into two haplogroups, A and B, of which B is predominant; (ii) four of the ancient mtDNA haplotypes are novel; (iii) higher mtDNA haplotype diversity occurred during the Middle Ages as compared to other periods, a fact concordant with the historical context of expanding international trade during the Middle Ages; (iv) the proportion of rarer haplotypes declined during the expansion of sheep from the Near Eastern domestication centre to the northern European region; (v) three male samples showed the presence of the characteristic northern European haplotype, SNP G‐oY1 of the Y‐chromosome, and represent the earliest occurrence of this haplotype. Our results provide the first insight into the genetic diversity and phylogeographical background of ancient sheep in Estonia and provide basis for further studies on the temporal fluctuations of ancient sheep populations.     

Mitochondrial DNA phylogeography of least cisco Coregonus sardinella in Alaska

This study presents the first detailed analysis of the mitochondrial DNA diversity of least cisco Coregonus sardinella in Alaska using a 678 bp segment of the control region (D‐loop) of the mitochondrial genome. Findings suggest that the history of C. sardinella in Alaska differs from that of other species of Coregonus present in the state and surrounding regions. The examined populations of C. sardinella are genetically diverse across Alaska. Sixty‐eight distinct mitochondrial haplotypes were identified among 305 individuals sampled from nine locations. The haplotype minimum spanning network and phylogeny showed a modest level of geographic segregation among haplotypes, suggesting high levels of on‐going or recent connectivity among distant populations. Observed Φ_ST values and the results of homogeneity and AMOVAs indicate incipient genetic differentiation between aggregations in three broad regional groups. Sites north of the Brooks Range formed one group, sites in the Yukon and Selawik Rivers formed a second group and sites south of the Yukon drainage formed the third group. Overall, the sequence data showed that a large proportion of mtDNA genetic variation in C. sardinella is shared across Alaska, but this variation is not homogeneously distributed across all regions and for all haplotype groups.     

Genetic diversity of the Chinese goat in the littoral zone of the Yangtze River as assessed by microsatellite and mtDNA

The objective of this study was to assess the genetic diversity and population structure of goats in the Yangtze River region using microsatellite and mtDNA to better understand the current status of those goat genetic diversity and the effects of natural landscape in fashion of domestic animal genetic diversity. The genetic variability of 16 goat populations in the littoral zone of the Yangtze River was estimated using 21 autosomal microsatellites, which revealed high diversity and genetic population clustering with a dispersed geographical distribution. A phylogenetic analysis of the mitochondrial D‐loop region (482 bp) was conducted in 494 goats from the Yangtze River region. In total, 117 SNPs were reconstructed, and 173 haplotypes were identified, 94.5% of which belonged to lineages A and B. Lineages C, D, and G had lower frequencies (5.2%), and lineage F haplotypes were undetected. Several high‐frequency haplotypes were shared by different ecogeographically distributed populations, and the close phylogenetic relationships among certain low‐frequency haplotypes indicated the historical exchange of genetic material among these populations. In particular, the lineage G haplotype suggests that some west Asian goat genetic material may have been transferred to China via Muslim migration.     

Analysis of Amblyomma sculptum haplotypes in an area endemic for Brazilian spotted fever

Amblyomma sculptum (Ixodida: Ixodidae) Berlese, 1888, a member of the Amblyomma cajennense complex, is the major vector of Brazilian spotted fever (BSF) in southeastern Brazil. In this study, the genetic diversity of A. sculptum populations in the state of Rio de Janeiro (RJ), Brazil, was investigated because genetic variability in tick populations may be related to vector competence. Samples of A. sculptum from 19 municipalities in 7 regions of RJ were subjected to DNA extraction, amplification and sequencing of D‐loop, cytochrome oxidase II and 12S rDNA mitochondrial genes. These sequences were used to map the genetic diversity of this tick. Amblyomma sculptum populations are genetically diverse in RJ, especially in the South Centre and Highland regions. Few unique haplotypes were observed in all populations, and the majority of genetic variation found was among ticks within each population. Phylogenetic reconstruction reinforced the assumption that all the haplotypes identified in RJ belong to A. sculptum. However, some RJ haplotypes are closer to A. sculptum from Argentina than to A. sculptum from elsewhere in Brazil. In RJ, A. sculptum has high genetic diversity, although little genetic differentiation. Observations also indicated a high level of gene flow among the studied populations and no evidence of population structure according to region in RJ.     

Genetic Variations of 13 Indigenous Chinese Goat Breeds Based on Cytochrome b Gene Sequences

Phylogenetic relationships among and genetic variability within 13 Chinese indigenous goat breeds and Boer goat were analyzed using cytochrome b gene sequences. There were 44 variable sites found in a 642 bp sequence, and 46 Cyt b haplotypes were subsequently defined. The phylogeny analysis of haplotypes in combination with goat Cyt b sequences from GenBank shows that Chinese goats are obviously separated from wild goats and might come from Capra aegagrus. Further analysis indicated that indigenous Chinese goats might descend from at least two lineages; most of the individuals analyzed could be classified into lineage A as defined by Luikart, but five other goats were of uncertain lineage. The Tibet plateau is a possible place of origin for Chinese goats. The neighbor-joining tree based on pairwise differences among populations shows that most Tibetan goats, except the Middle Tibet type, cluster closely with North China goats, and then with South China goats. This result confirms that differences in genetic structure exist among goats in different geographic locations. Nucleotide diversity varied among populations. Tibet and North China goats had higher genetic diversity than South China goats. The fixation index (F _st=87.72%) suggested that most of the total genetic variation was due to variation within populations. In addition, the results indicate that Cyt b gene sequence information alone might not be enough for phylogeny analysis among breeds within species, as shown by fewer polymorphic sites and lower bootstrap values on the neighbor-joining tree.     

Low genetic diversity and significant structuring of the common hamster populations <Emphasis Type="Italic">Cricetus cricetus</Emphasis> in Poland revealed by the mtDNA control region sequence variation

The genetic diversity of 12 populations in the present range of the common hamster Cricetus cricetus (Linnaeus, 1758) in Poland was established. The 366 bp of the mtDNA control region was sequenced for 195 individuals. As few as seven haplotypes were found and their distribution was geographically structured. The large geographic areas were fixed or almost fixed for a single haplotype and three groups of populations, that do not share any haplotypes, have been defined. Proportions of genetic diversity attributable to variation between groups of populations, between populations within groups and within populations were 93.64, 1.92 and 4.45% (SAMOVA: p < 0.001 for all estimates), respectively. Such pattern of variation is most probably the result of historical, postglacial bottlenecks and present genetic drift after the population decline in the last few decades.     

中国狼不同地理种群遗传多样性及系统进化分析

为了解中国狼不同地理种群遗传多样性及系统发育情况,从中国境内狼的主要分布区青海、新疆、内蒙古和吉林4个地区采集样品,用分子生物学技术手段成功地获得44个个体线粒体DNA控制区第一高变区(HVRⅠ)序列和40个线粒体Cyt b部分序列。线粒体控制区HVRⅠ共检测到51个变异位点,位点变异率为8.76%;线粒体Cyt b部分序列发现31个变异位点,位点变异率为5.33%,未见插入及缺失现象,变异类型全部为碱基置换。共定义了16个线粒体HVRⅠ单倍型,其中吉林与内蒙种群存在共享单倍型,估计这两地间种群亲缘关系较近。4个地理种群中新疆种群拥有较高的遗传多样性(0.94)。中国狼种群总体平均核苷酸多态性为2.27%,与世界其他国家地区相比,中国狼种群拥有相对较高的遗传多样性。通过线粒体HVRⅠ单倍型构建的系统进化树可以看出,中国狼在进化上分为2大支,其中位于青藏高原的青海种群独立为一支,推测其可能长期作为独立种群进化。基于青海种群与新疆,内蒙种群的线粒体Cyt b遗传距离,推测中国狼2个世系可能在更新世冰川时期青藏高原受地质作用急速隆起后出现分歧,分歧时间大约在1.1 MY前。     

Genetic diversity and conservation units in wild and captive populations of endangered freshwater fishes: a case of <Emphasis Type="Italic">Hemigrammocypris rasborella</Emphasis> in Shizuoka,Japan

Population structure and genetic diversity were examined using partial mitochondrial cytochrome b gene sequences of four wild, one reintroduced, and five captive populations of the endangered cyprinid Hemigrammocypris rasborella from three river systems in the easternmost region of the species’ range in Shizuoka Prefecture, central Honshu, Japan. We detected loss of genetic diversity from portions of the wild and captive populations, as well as suspected nonindigenous haplotypes in some captive, reintroduced, and even wild populations. Given the population structure revealed, we suggest that the populations should be managed with consideration for both the endemism and viability (avoidance of inbreeding depression) of the local populations.     

Genetic reconstruction of a bullfrog invasion to elucidate vectors of introduction and secondary spread

Reconstructing historical colonization pathways of an invasive species is critical for uncovering factors that determine invasion success and for designing management strategies. The American bullfrog (Lithobates catesbeianus) is endemic to eastern North America, but now has a global distribution and is considered to be one of the worst invaders in the world. In Montana, several introduced populations have been reported, but little is known of their sources and vectors of introduction and secondary spread. We evaluated the genetic composition of introduced populations at local (Yellowstone River floodplain) and regional (Montana and Wyoming) scales in contrast to native range populations. Our objectives were to (1) estimate the number of introductions, (2) identify probable native sources, (3) evaluate genetic variation relative to sources, and (4) characterize properties of local‐ and regional‐scale spread. We sequenced 937 bp of the mitochondrial cytochrome b locus in 395 tadpoles collected along 100 km of the Yellowstone River, from three additional sites in MT and a proximate site in WY. Pairwise Φ_ST revealed high divergence among nonnative populations, suggesting at least four independent introductions into MT from diverse sources. Three cyt b haplotypes were identical to native haplotypes distributed across the Midwest and Great Lakes regions, and AMOVA confirmed the western native region as a likely source. While haplotype (H_d = 0.69) and nucleotide diversity (π = 0.005) were low in introduced bullfrogs, the levels of diversity did not differ significantly from source populations. In the Yellowstone, two identified haplotypes implied few introduction vectors and a significant relationship between genetic and river distance was found. Evidence for multiple invasions and lack of subsequent regional spread emphasizes the importance of enforcing legislation prohibiting bullfrog importation and the need for continuing public education to prevent transport of bullfrogs in MT. More broadly, this study demonstrates how genetic approaches can reveal key properties of a biological invasion to inform management strategies.     

PCR–RFLP analysis of mitochondrial DNA in tench Tinca tinca

Polymorphism was detected at ND1, ND6, D‐loop and cyt b segments of mtDNA in 105 tench (Tinca tinca L.), using the polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) technique and five composite haplotypes were identified. The diversity indices and the results of the population comparisons revealed that the identified markers provide a powerful tool for further studies on this species.     

Genetic diversity of Chinese cattle revealed by Y‐SNP and Y‐STR markers

With its vast territory and complex natural environment, China boasts rich cattle genetic resources. To gain the further insight into the genetic diversity and paternal origins of Chinese cattle, we analyzed the polymorphism of Y‐SNPs (UTY19 and ZFY10) and Y‐STRs (INRA189 and BM861) in 34 Chinese cattle breeds/populations, including 606 males representative of 24 cattle breeds/populations collected in this study as well as previously published data for 302 bulls. Combined genotypic data identified 14 Y‐chromosome haplotypes that represented three haplogroups. Y2‐104‐158 and Y2‐102‐158 were the most common taurine haplotypes detected mainly in northern and central China, whereas the indicine haplotype Y3‐88‐156 predominates in southern China. Haplotypes Y2‐108‐158, Y2‐110‐158, Y2‐112‐158 and Y3‐92‐156 were private to Chinese cattle. The population structure revealed by multidimensional scaling analysis differentiated Tibetan cattle from the other three groups of cattle. Analysis of molecular variance showed that the majority of the genetic variation was explained by the genetic differences among groups. Overall, our study indicates that Chinese cattle retain high paternal diversity (H = 0.607 ± 0.016) and probably much of the original lineages that derived from the domestication center in the Near East without strong admixture from commercial cattle carrying Y1 haplotypes.     

Do diversification models of Madagascar’s biota explain the population structure of the endemic bat Myotis goudoti (Chiroptera: Vespertilionidae)?

Aim Three mechanisms have been proposed to explain the adaptive radiations and species diversifications of Madagascar’s biota: the ecogeographical constraint, the riverine barrier and the micro‐endemism models. On the intraspecific level, each model predicts different patterns of gene flow across the island’s physical and ecological features. To evaluate these models, phylogeographical analyses were conducted on a widespread and endemic species of bat, Myotis goudoti (Vespertilionidae). Location Madagascar. Methods In order to reconstruct the phylogeographical history of M. goudoti, the mitochondrial D‐loop and the cytochrome b gene were sequenced for 195 bats from 41 localities. Phylogenetic reconstructions and a minimum spanning tree were used to infer haplotype relationships. The effect of barriers on gene flow was evaluated using analyses of molecular variance and pairwise population differentiation. Mismatch distribution and coalescence‐based estimates were conducted to infer the demographic history of M. goudoti. Results The sequenced individuals showed 159 distinct D‐loop haplotypes, most of them being unique to a single location. Populations were significantly structured (Φ_ST = 0.170, P < 0.001) across Madagascar, but only a minor part of the overall genetic variance was explained by any of the three models. Shared ancestry of lineages across most physical or ecological barriers was common, whereas the uncovered genetic differences between southern and central‐northern populations were unexpected. Main conclusions Major barriers predicted by the three biogeographical models do not explain the segregation of mitochondrial lineages of M. goudoti across Madagascar. This is not simply attributable to the high dispersal ability of this species, as populations are notably structured. The genetic contrast between southern and central‐northern populations, separated by a zone of admixture, suggests that these areas currently support populations that expanded during the Late Pleistocene. This latitudinal differentiation of populations has been observed in less vagile animals, such as geckos and lemurs, suggesting that climate fluctuations of the Pleistocene had an impact across several groups and resulted in northern and southern refugia in Madagascar.     

Unique haplotypes in ant‐attended aphids and widespread haplotypes in non‐attended aphids

Aphid species within the genus Tuberculatus Mordvilko (Hemiptera: Aphididae) exhibit a variety of interactions with ants, ranging from close associations to non‐attendance. A previous study indicated that despite wing possession, ant‐attended Tuberculatus species exhibited low dispersal rates compared with non‐attended species. This study examined if presence or absence of mutualistic interactions and habitat continuity of host plants affected intraspecific genetic diversity and genetic differentiation in mitochondrial DNA cytochrome oxidase I (COI) sequences. Sympatric ant‐attended Tuberculatus quercicola (Matsumura) (Hemiptera: Aphididae) and non‐attended Tuberculatus paiki Hille Ris Lambers (Hemiptera: Aphididae) were collected from the daimyo oak Quercus dentata Thunberg (Fagales: Fagaceae) in Japan and examined for haplotype variability. Seventeen haplotypes were identified in 568 T. quercicola individuals representing 23 populations and seven haplotypes in 425 T. paiki representing 19 populations. Haplotype diversity, which indicates the mean number of differences between all pairs of haplotypes in the sample, and nucleotide diversity were higher in T. quercicola than T. paiki. Analysis of molecular variance (AMOVA) showed higher genetic differentiation among populations within groups of T. quercicola (39.8%) than T. paiki (22.6%). The effects of attendant ant species on genetic differentiation in T. quercicola were not distinguishable from geographic factors. Despite low dispersal rates, host plant habitat continuity might facilitate widespread dispersal of a T. quercicola haplotype in Hokkaido. These results suggested that following T. quercicola colonization, gene flow among populations was limited, resulting in genetic drift within populations. However, frequent T. paiki dispersal is clearly evident by low genetic differentiation among populations within groups, resulting in lower haplotype diversity.     

Diversity and structure of fragmented populations of a threatened endemic cyprinodontid (Aphanius sophiae) inferred from genetics and otolith morphology: Implications for conservation and management

The assessment of population structure and genetic diversity is crucial for the management and conservation of threatened species. Natural and artificial barriers to dispersal (i.e., gene flow) increase populations’ differentiation and isolation by reducing genetic exchange and diversity. Freshwater ecosystems are highly fragmented because of human activities. Threatened species with small population sizes are more sensitive to habitat fragmentation effects. Here, we investigate the genetic population structure and gene flow among seven populations of Aphanius sophiae in the Kor Basin by using sequences of the complete Cyt b gene and otolith morphometry. The Cyt b gene showed low level of genetic variation, only 4.12% of the identified sites were variable, and 2.42% were parsimony informative. Overall, haplotype diversity was low to moderate and nucleotide diversity was low to extremely low. Fish populations exhibited high levels of genetic differentiation, suggesting limited gene flow among them. These differences were obtained not only among geographically distant populations, but also among neighboring localities. Genetic population structure was supported by the AMOVA analysis and by the haplotype network (only one of 21 haplotypes were shared by two localities). Otolith morphometric analysis was in agreement with genetic results, the two most distant and isolated populations were clearly separated, and genetically close populations showed less differences in morphometry. A significant pattern of isolation by distance was also detected among A. sophiae populations, with genetic distance more correlated with hydrological distance than with geographic distance. Results suggested that limited gene flow due to habitat fragmentation is an important factor contributing to genetic structuring and to the loss of genetic variation of A. sophiae populations. Aphanius sophiae population structure seems to be the result of habitat fragmentation and water pollution, but other factors such as introduced species should be considered. Given the high degree of genetic structuring, the definition of conservation groups is of particular importance for A. sophiae, which should be considered endangered according to the IUCN criteria. Conservation plans must recognize the genetic independence of populations and manage separately preventing the loss of locally adapted genotypes.     

Conservation genetics assessment and phylogenetic relationships of critically endangered Hucho bleekeri in China

Hucho bleekeri is a critically endangered salmonid fish found in the Yangtze River drainage in China. In this study, the genetic diversity of a small population (n = 43) was first assessed with partial mitochondrial DNA sequences (D‐loop region and a cytochrome b gene [CYTB] gene fragment) and 15 microsatellite markers. Low levels of nucleotide diversity (Pi) were demonstrated in the H. bleekeri population based on the two mitochondrial DNA markers. The number of haplotypes (h) and the haplotype diversity (Hd) in the D‐loop region (12 haplotypes and Hd = 0.8208) were higher than in the CYTB gene fragment (three haplotypes and Hd = 0.0941). The number of microsatellite alleles (Na) ranged from 2 to 13 in these individuals. The mean observed heterozygosity (Ho) and the expected heterozygosity (He) were 0.59719 and 0.44735, respectively. Analysis of molecular variance showed that the degree of differentiation in the population was low (F_ST = 0.04041) and the coefficient of inbreeding (F_IS) was negative, indicating no obvious evidence of inbreeding in this population. A demographic assessment suggested that this species expanded a long time ago, but has suffered great losses in recent years. A molecular phylogenetic analysis clearly indicated that H. bleekeri is not introgressed by Brachymystax lenok tsinlingensis. The baseline population genetic information supplied by this study will be vital in monitoring this highly threatened species.     

Genetic variation in the mitochondrial 16S rRNA gene of the American dog tick,Dermacentor variabilis (Acari: Ixodidae)

Genetic variation in the mitochondrial (mt) 16S ribosomal RNA (rRNA) gene was examined for the American dog tick, Dermacentor variabilis (Say, 1821). Nine different haplotypes were detected among 369 adult D. variabilis collected from four localities in Canada. There were eight variable nucleotide positions in the 404 bp sequence alignment. Individuals of haplotype 1 occurred at frequency of >75% at all localities. Five haplotypes were detected at only one of the four localities. High haplotype diversity and low nucleotide diversity, combined with significantly negative F_s values for ticks at three localities, suggest a recent population expansion. Genetic differences were found between populations at different localities, but a Mantel regression analysis revealed no association between genetic differences and geographical distances. There was also no association between tick haplotype and the prevalence of the bacterium, Rickettsia montanensis Weiss and Moulder, 1984, in D. variabilis among localities or on opposite sides of Blackstrap Lake (Saskatchewan). The 16S rDNA haplotypes from Canadian populations of D. variabilis formed a clade with those from the eastern and central U.S.A., to the exclusion of D. variabilis from geographically isolated populations in the western U.S.A. Although sample sizes for D. variabilis in the eastern U.S.A. are small, there may be genetic divergence between populations in Canada and those in the eastern U.S.A., which may have implications for studies on the pathogenic agents transmitted by D. variabilis to its hosts.