Age- and Sex-Associated Pathogenesis of Cell Culture-Passaged Kemerovo Virus in IFNAR(-/-) Mice.

Kemerovo virus (KEMV) is a tick-borne orbivirus transmitted by ticks of the genus Ixodes. Previous animal experimentation studies with orbiviruses, in particular the interferon receptor double knock-out (IFNAR(-/-)) mouse model, did not indicate bias that is related to age or sex. We endeavoured to assess the effect of serial and alternated passages of KEMV in mammalian or Ixodes cells on virus replication and potential virulence in male or female IFNAR(-/-) mice, with important age differences: younger males (4-5 months old), older males (14-15 months old), and old females (14-15 months old). After 30 serial passages in mammalian or tick cells, or alternated passages in the two cell types, older female mice which were inoculated with the resulting virus strains were the first to show clinical signs and die. Younger males behaved differently from older males whether they were inoculated with the parental strain of KEMV or with any of the cell culture-passaged strains. The groups of male and female mice inoculated with the mammalian cell culture-adapted KEMV showed the lowest viraemia. While older female and younger male mice died by day 6 post-inoculation, surprisingly, the older males survived until the end of the experiment, which lasted 10 days. RNA extracted from blood and organs of the various mice was tested by probe-based KEMV real-time RT-PCR. Ct values of the RNA extracts were comparable between older females and younger males, while the values for older males were >5 Ct units higher for the various organs, indicating lower levels of replication. It is noteworthy that the hearts of the old males were the only organs that were negative for KEMV RNA. These results suggest, for the first time, an intriguing age- and sex-related bias for an orbivirus in this animal model. Changes in the amino acid sequence of the RNA-dependent RNA polymerase of Kemerovo virus, derived from the first serial passage in Ixodes cells (KEMV Ps.IRE1), were identified in the vicinity of the active polymerase site. This finding suggests that selection of a subpopulation of KEMV with better replication fitness in tick cells occurred.

Genetic structure and Rickettsia infection rates in Ixodes ovatus and Haemaphysalis flava ticks across different altitudes.

Ixodid ticks, such as Ixodes ovatus and Haemaphysalis flava, are important vectors of tick-borne diseases in Japan, such as Japanese spotted fever caused by Rickettsia japonica. This study describes the Rickettsia infection rates influenced by the population genetic structure of I.ovatus and H. flava along an altitudinal gradient. A total of 346 adult I. ovatus and 243 H. flava were analyzed for the presence of Rickettsia by nested PCR targeting the 17kDA, gltA, rOmpA, and rOmpB genes. The population genetic structure was analyzed utilizing the mitochondrial cytochrome oxidase 1 (cox1) marker. The Rickettsia infection rates were 13.26% in I. ovatus and 6.17% in H. flava. For I. ovatus, the global FST value revealed significant genetic differentiation among the different populations, whereas H. flava showed non-significant genetic differentiation. The cox1 I. ovatus cluster dendrogram showed two cluster groups, while the haplotype network and phylogenetic tree showed three genetic groups. A significant difference was observed in Rickettsia infection rates and mean altitude per group between the two cluster groups and the three genetic groups identified within I. ovatus. No significant differences were found in the mean altitude or Rickettsia infection rates of H. flava. Our results suggest a potential correlation between the low gene flow in I. ovatus populations and the spatially heterogeneous Rickettsia infection rates observed along the altitudinal gradient. This information can be used in understanding the relationship between the tick vector, its pathogen, and environmental factors, such as altitude, and for the control of tick-borne diseases in Japan.

mRNA vaccination of rabbits alters the fecundity, but not the attachment, of adult Ixodes scapularis.

19ISP is a nucleoside-modified mRNA-lipid nanoparticle vaccine that targets 19 Ixodes scapularis proteins. We demonstrate that adult I. scapularis have impaired fecundity when allowed to engorge on 19ISP-immunized rabbits. 19ISP, therefore, has the potential to interrupt the tick reproductive cycle, without triggering some of the other effects associated with acquired tick resistance. This may lead to the development of new strategies to reduce I. scapularis populations in endemic areas.

From the fat body to the hemolymph: Profiling tick immune and storage proteins through transcriptomics and proteomics.

Ticks are blood-feeding arachnids that are known to transmit various pathogenic microorganisms to their hosts. During blood feeding, ticks activate their metabolism and immune system to efficiently utilise nutrients from the host's blood and complete the feeding process. In contrast to insects, in which the fat body is known to be a central organ that controls essential metabolic processes and immune defense mechanisms, the function of the fat body in tick physiology is still relatively unexplored. To fill this gap, we sought to uncover the repertoire of genes expressed in the fat body associated with trachea (FB/Tr) by analyzing the transcriptome of individual, partially fed (previtellogenic) Ixodes ricinus females. The resulting catalog of individual mRNA sequences reveals a broad repertoire of transcripts encoding proteins involved in nutrient storage and distribution, as well as components of the tick immune system. To gain a detailed insight into the secretory products of FB/Tr specifically involved in inter-tissue transport and humoral immunity, the transcriptomic data were complemented with the proteome of soluble proteins in the hemolymph of partially fed female ticks. Among these proteins, the hemolipoglyco-carrier proteins were predominant. When comparing immune peptides and proteins from the fat body with those produced by hemocytes, we found that the fat body serves as a unique producer of certain immune components. Finally, time-resolved transcriptional regulation of selected immune transcripts from the FB/Tr was examined in response to experimental challenges with model microbes and analyzed by RT-qPCR. Overall, our data show that the fat body of ticks, similar to insects, is an important metabolic tissue that also plays a remarkable role in immune defense against invading microbes. These findings improve our understanding of tick biology and its impact on the transmission of tick-borne pathogens.

Comparative analysis of Presence-Absence gene Variations in five hard tick species: impact and functional considerations.

Tick species are vectors of harmful human and animal diseases, and their expansion is raising concerns under the global environmental changes' scenario. Ticks host and transmit bacteria, protozoa and viruses, making the understanding of host-pathogen molecular pathways critical to development of effective disease control strategies. Despite the considerable sizes and repeat contents of tick genomes, individual tick genomics is perhaps the most effective approach to reveal genotypic traits of interest. Presence-Absence gene Variations (PAVs) can contribute to individual differences within species, with dispensable genes carried by subsets of individuals possibly underpinning functional significance at individual or population-levels. We exploited 350 resequencing datasets of Dermacentor silvarum, Haemaphysalis longicornis, Ixodes persulcatus, Rhipicephalus microplus and Rhipicephalus sanguineus hard tick specimens to reveal the extension of PAV and the conservation of dispensable genes among individuals and, comparatively, between species. Overall, we traced 550-3,346 dispensable genes per species and were able to reconstruct 5.3-7 Mb of genomic regions not included in the respective reference genomes, as part of the tick pangenomes. Both dispensable genes and de novo predicted genes indicated that PAVs preferentially impacted mobile genetic elements in these tick species.

Diversity and host specificity of Borrelia burgdorferi's outer surface protein C (ospC) alleles in synanthropic mammals, with a notable ospC allele U absence from mixed infections.

Interactions among pathogen genotypes that vary in host specificity may affect overall transmission dynamics in multi-host systems. Borrelia burgdorferi, a bacterium that causes Lyme disease, is typically transmitted among wildlife by Ixodes ticks. Despite the existence of many alleles of B. burgdorferi's sensu stricto outer surface protein C (ospC) gene, most human infections are caused by a small number of ospC alleles ["human infectious alleles" (HIAs)], suggesting variation in host specificity associated with ospC. To characterize the wildlife host association of B. burgdorferi's ospC alleles, we used metagenomics to sequence ospC alleles from 68 infected individuals belonging to eight mammalian species trapped at three sites in suburban New Brunswick, New Jersey (USA). We found that multiple allele ("mixed") infections were common. HIAs were most common in mice (Peromyscus spp.) and only one HIA was detected at a site where mice were rarely captured. ospC allele U was exclusively found in chipmunks (Tamias striatus), and although a significant number of different alleles were observed in chipmunks, including HIAs, allele U never co-occurred with other alleles in mixed infections. Our results suggest that allele U may be excluding other alleles, thereby reducing the capacity of chipmunks to act as reservoirs for HIAs.

Cross-alteration of murine skin and tick microbiome concomitant with pathogen transmission after Ixodes ricinus bite.

BACKGROUND: Ticks are major vectors of diseases affecting humans such as Lyme disease or domestic animals such as anaplasmosis. Cross-alteration of the vertebrate host skin microbiome and the tick microbiome may be essential during the process of tick feeding and for the mechanism of pathogen transmission. However, it has been poorly investigated. METHODS: We used mice bitten by field-collected ticks (nymphs and adult ticks) in different experimental conditions to investigate, by 16S rRNA gene metabarcoding, the impact of blood feeding on both the mouse skin microbiome and the tick microbiome. We also investigated by PCR and 16S rRNA gene metabarcoding, the diversity of microorganisms transmitted to the host during the process of tick bite at the skin interface and the dissemination of the pathogen in host tissues (blood, heart, and spleen). RESULTS: Most of the commensal bacteria present in the skin of control mice were replaced during the blood-feeding process by bacteria originating from the ticks. The microbiome of the ticks was also impacted by the blood feeding. Several pathogens including tick-borne pathogens (Borrelia/Borreliella, Anaplasma, Neoehrlichia, Rickettsia) and opportunistic bacteria (Williamsia) were transmitted to the skin microbiome and some of them disseminated to the blood or spleen of the mice. In the different experiments of this study, skin microbiome alteration and Borrelia/Borreliella transmission were different depending on the tick stages (nymphs or adult female ticks). CONCLUSIONS: Host skin microbiome at the bite site was deeply impacted by the tick bite, to an extent which suggests a role in the tick feeding, in the pathogen transmission, and a potentially important impact on the skin physiopathology. The diversified taxonomic profiles of the tick microbiome were also modified by the blood feeding. Video Abstract.

Bulk and single-nucleus RNA sequencing highlight immune pathways induced in individuals during an Ixodes scapularis tick bite.

Ticks are hematophagous arthropods that use a complex mixture of salivary proteins to evade host defenses while taking a blood meal. Little is known about the immunological and physiological consequences of tick feeding on humans. Here, we performed the first bulk and single-nucleus RNA sequencing (snRNA-seq) of skin and blood of four persons presenting with naturally acquired, attached Ixodes scapularis ticks. Pathways and individual genes associated with innate and adaptive immunity were identified based on bulk RNA sequencing, including interleukin-17 signaling and platelet activation pathways at the site of tick attachment or in peripheral blood. snRNA-seq further revealed that the Hippo signaling, cell adhesion, and axon guidance pathways were involved in the response to an I. scapularis bite in humans. Features of the host response in these individuals also overlapped with that of laboratory guinea pigs exposed to I. scapularis and which acquired resistance to ticks. These findings offer novel insights for the development of new biomarkers for I. scapularis exposure and anti-tick vaccines for human use.

The highly improved genome of Ixodes scapularis with X and Y pseudochromosomes.

Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the ∼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.

Molecular-phylogenetic analyses of Ixodes species from South Africa suggest an African origin of bird-associated exophilic ticks (subgenus Trichotoixodes).

BACKGROUND: Among hard ticks (Acari: Ixodidae), the genus Ixodes comprises the highest number of species, which in turn are most numerous in the Afrotropical zoogeographic region. In South Africa extensive morphological studies have been performed on Ixodes species but only few reports included molecular analyses. METHODS: In this study, 58 Ixodes spp. ticks, collected from ten mammalian and eight avian host species in South Africa, were molecularly and phylogenetically analyzed. In addition, a newly collected sample of the Palearctic Ixodes trianguliceps was included in the analyses. RESULTS: Among the ticks from South Africa, 11 species were identified morphologically. The majority of ticks from mammals represented the Ixodes pilosus group with two species (n = 20), followed by ticks resembling Ixodes rubicundus (n = 18) and Ixodes alluaudi (n = 3). In addition, single specimens of Ixodes rhabdomysae, Ixodes ugandanus, Ixodes nairobiensis and Ixodes simplex were also found. Considering bird-infesting ticks, Ixodes theilerae (n = 7), Ixodes uriae (n = 4) and ticks most similar to Ixodes daveyi (provisionally named I. cf. daveyi, n = 2) were identified. Molecular analyses confirmed two species in the I. pilosus group and a new species (I. cf. rubicundus) closely related to I. rubicundus sensu stricto. Phylogenetic trees based on concatenated mitochondrial or mitochondrial and nuclear gene sequences indicated that the subgenus Afrixodes forms a monophyletic clade with bird-associated exophilic ticks (subgenus Trichotoixodes). Ixodes trianguliceps clustered separately whereas I. alluaudi with their morphologically assigned subgenus, Exopalpiger. CONCLUSIONS: Phylogenetic analyses shed new lights on the relationships of Ixodes subgenera when including multiple sequences from subgenus Afrixodes and African as well as Palearctic species of subgenera Trichotoixodes and Exopalpiger. Subgenera Afrixodes and bird-associated Trichotoixodes share common ancestry, suggesting that the latter might have also originated in Africa. Regarding the subgenus Exopalpiger, I. alluaudi is properly assigned as it clusters among different Australian Ixodes, whereas I. trianguliceps should be excluded.

Genetic modification, characterization, and co-infection of Franken Sphingomonas and Anaplasma phagocytophilum in tick cells.

Paratransgenesis through genetic manipulation of symbiotic or commensal microorganisms has been proposed as an effective and environmentally sound approach for the control of vector-borne diseases, including tick bite-related pathologies, and reducing pathogen transmission. Here, we present a protocol for Sphingomonas transformation with Anaplasma phagocytophilum major surface protein 4 and heat shock protein 70. We describe a step-by-step protocol for in vitro study of interactions between transformed Franken Sphingomonas and Ixodes scapularis ISE6 tick cells during A. phagocytophilum infection. For complete details on the use and execution of this protocol, please refer to Mazuecos et al. (2023).1.

Anaplasma phagocytophilum in European bison (Bison bonasus) and their ticks from Lithuania and Poland.

The increasing population of European bison (Bison bonasus) can contribute to the prevalence of zoonotic pathogens. The aim of the present study was to assess the presence of A. phagocytophilum infection in European bison tissues as well as ticks removed from European bison in Lithuania and Poland. A further objective of this work was to compare the detected A. phagocytophilum strains. A total of 85 tissue samples (spleen) of European bison and 560 ticks belonging to two species, Ixodes ricinus (n = 408) and Dermacentor reticulatus (n = 152) were tested. DNA of A. phagocytophilum was detected based on RT-PCR in 40% of the European bison samples, 8.8% of the I. ricinus and 5.9% of the D. reticulatus ticks. Analysis of the obtained partial 16S rRNA gene sequences of A. phagocytophilum revealed the presence of three variants with two polymorphic sites. Furthermore, phylogenetic analysis with partial msp4 gene sequences grouped A. phagocytophilum variants into three clusters. This study revealed that the groEL gene sequences of A. phagocytophilum from European bison and their ticks grouped into ecotype I and only one sequence from Lithuanian European bison belonged to ecotype II. The results of the present study indicated that European bison may play a role as a natural reservoir of A. phagocytophilum.

Specific mRNA lipid nanoparticles and acquired resistance to ticks.

Acquired resistance to ticks can develop when animals are repeatedly exposed to ticks. Recently, acquired resistance to Ixodes scapularis was induced in guinea pigs immunized with an mRNA-lipid nanoparticle vaccine (19ISP) encoding 19 I. scapularis proteins. Here, we evaluated specific mRNAs present in 19ISP to identify critical components associated with resistance to ticks. A lipid nanoparticle containing 12 mRNAs which included all the targets within 19ISP that elicited strong humoral responses in guinea pigs, was sufficient to induce robust resistance to ticks. Lipid nanoparticles containing fewer mRNAs or a single mRNA were not able to generate strong resistance to ticks. All lipid nanoparticles containing salp14 mRNA, however, were associated with increased redness at the tick bite site - which is the first manifestation of acquired resistance to ticks. This study demonstrates that more than one I. scapularis target within 19ISP is required for resistance to ticks, and that additional targets may also play a role in this process.

German Ixodes inopinatus samples may not actually represent this tick species.

Ticks are important vectors of human and animal pathogens, but many questions remain unanswered regarding their taxonomy. Molecular sequencing methods have allowed research to start understanding the evolutionary history of even closely related tick species. Ixodes inopinatus is considered a sister species and highly similar to Ixodes ricinus, an important vector of many tick-borne pathogens in Europe, but identification between these species remains ambiguous with disagreement on the geographic extent of I. inopinatus. In 2018-2019, 1583 ticks were collected from breeding great tits (Parus major) in southern Germany, of which 45 were later morphologically identified as I. inopinatus. We aimed to confirm morphological identification using molecular tools. Utilizing two genetic markers (16S rRNA, TROSPA) and whole genome sequencing of specific ticks (n = 8), we were able to determine that German samples, morphologically identified as I. inopinatus, genetically represent I. ricinus regardless of previous morphological identification, and most likely are not I. ricinus/I. inopinatus hybrids. Further, our results showed that the entire mitochondrial genome, let alone singular mitochondrial genes (i.e., 16S), is unable to distinguish between I. ricinus and I. inopinatus. Our results suggest that I. inopinatus is geographically isolated as a species (northern Africa and potentially southern Spain and Portugal) and brings into question whether I. inopinatus exists in central Europe. Our results highlight the probable existence of I. inopinatus and the power of utilizing genomic data in answering questions regarding tick taxonomy.

A longitudinal transcriptomic analysis from unfed to post-engorgement midguts of adult female Ixodes scapularis.

The hematophagy behavior has evolved independently several times within the Arthropoda phylum. Interestingly, the process of acquiring a blood meal in ticks is considerably distinct from that observed in other blood-feeding arthropods. Instead of taking seconds to minutes to complete a blood meal, an adult female Ixodes scapularis tick can remain attached to its host for numerous days. During this extended feeding period, the tick undergoes drastic morphological changes. It is well established that the tick midgut plays a pivotal role not only in blood meal digestion but also in pathogen acquisition and transmission. However, our understanding of the underlying molecular mechanisms involved in these events remains limited. To expedite tick research, we conducted a comprehensive longitudinal RNA-sequencing of the tick midgut before, during, and after feeding. By collecting ticks in different feeding stages (unfed, slow feeding, rapid feeding, and early post-detached), we obtained a comprehensive overview of the transcripts present in each stage and the dynamic transcriptional changes that occur between them. This provides valuable insights into tick physiology. Additionally, through unsupervised clustering, we identified transcripts with similar patterns and stage-specific sequences. These findings serve as a foundation for selecting targets in the development of anti-tick control strategies and facilitate a better understanding of how blood feeding and pathogen infection impact tick physiology.

Detection of Bartonella DNA in Yellow Flies, Lone Star Ticks, and a Human Patient with Concurrent Evidence of Borrelia burgdorferi Infection in Northeast Florida, USA.

Background: Bartonella species and Borrelia burgdorferi sensu lato (Bbsl) are emerging zoonotic pathogens. The vectors and frequency of infections with both pathogen groups in the southern United States is understudied. This study describes an investigation of Bartonella and Bbsl in yellow flies collected at a residence in northeast Florida, USA, that led to subsequent discoveries of both organisms in lone star ticks (Amblyomma americanum) and a human patient. Materials and Methods: DNA samples from flies, ticks, and human patient blood specimens were tested via polymerase chain reaction assays for Bartonella or Bbsl species. DNA sequences were compared to reference strains for identification and characterization. Results: An exploratory investigation of arthropod-borne pathogens in yellow flies collected at a residence in northeast Florida revealed the presence of uncharacterized Bartonella species DNA sequences similar to ones previously detected in two lone star ticks from Virginia. Subsequent testing of several lone star ticks from the area detected similar sequences of Bartonella in three ticks. Testing of stored blood samples from a resident of the site, who had experienced chronic relapsing and remitting symptoms for over a decade, identified nearly identical Bartonella DNA sequences in multiple samples collected over a 10-year period. Two lone star ticks and several samples from the same patient and time period also tested positive for Bo. burgdorferi DNA, suggesting possible long-term coinfection of the patient with both organisms. Conclusion: This investigation identified highly similar Bartonella DNA sequences in yellow flies, lone star ticks, and a human patient in northeast Florida. Similarly, Bo. burgdorferi DNA was detected in two lone star ticks and multiple specimens from the patient. Positive PCR results from archived patient blood samples documented the presence of both organisms at multiple time points over more than a decade. More studies on human patients with chronic undefined illness and on the presence of Bartonella and Bbsl in hematophagous arthropods and animal hosts in the southeastern United States are needed.

Metagenomic surveillance for bacterial tick-borne pathogens using nanopore adaptive sampling.

Technological and computational advancements in the fields of genomics and bioinformatics are providing exciting new opportunities for pathogen discovery and genomic surveillance. In particular, single-molecule nucleotide sequence data originating from Oxford Nanopore Technologies (ONT) sequencing platforms can be bioinformatically leveraged, in real-time, for enhanced biosurveillance of a vast array of zoonoses. The recently released nanopore adaptive sampling (NAS) strategy facilitates immediate mapping of individual nucleotide molecules to a given reference as each molecule is being sequenced. User-defined thresholds then allow for the retention or rejection of specific molecules, informed by the real-time reference mapping results, as they are physically passing through a given sequencing nanopore. Here, we show how NAS can be used to selectively sequence DNA of multiple bacterial tick-borne pathogens circulating in wild populations of the blacklegged tick vector, Ixodes scapularis.

[Molecular and genetic characteristics of the multicomponent flavi-like Kindia tick virus (Flaviviridae) found in ixodes ticks on the territory of the Republic of Guinea].

INTRODUCTION: Ixodes ticks are vectors for pathogens of many infectious diseases. Recently, during the study of Rhipicephalus geigyi ticks collected from livestock in the Republic of Guinea, a new multicomponent flavi-like RNA virus, called Kindia tick virus (KITV), was discovered with an unusual mechanism for the implementation of genetic information. The aim of the work is to detect and study the genetic diversity of KITV in ixodes ticks collected in the territory of the Kindia province of the Republic of Guinea. MATERIAL AND METHODS: In 2021, 324 specimens of ticks of the species Amblyomma variegatum, Rh. geigyi, Rh. annulatus, Rh. decoloratus, Rh. senegalensis were collected from cattle. The detection of viral RNA was carried out in individual samples of ticks by RT-PCR, followed by the determination of the nucleotide sequence and phylogenetic analysis. RESULTS AND DISCUSSION: KITV detection rates in ticks of the species Rh. geigyi was 12.2%, Rh. annulatus 4.4%, Rh. decoloratus 3.3%. However, the KITV genetic material has not been identified in Am. variegatum ticks, which are one of the dominant species in West Africa. For all virus isolates, a partial nucleotide sequences of each of the four viral segments (GenBank, OK345271OK345306) were determined. The phylogenetic analysis showed a high level of identity (98.599.8%) for each of the four segments of the viral genome with those previously found in the Republic of Guinea. The obtained KITV isolates are most genetically close to Mogiana tick virus, which was previously detected in South America in Rh. microplus ticks and significantly differed from other multicomponent viruses circulating in Europe and Asia, including the Russian Federation. CONCLUSION: KITV genetic material was found in three species of ixodid ticks collected from livestock in a number of prefectures of the Republic of Guinea. The infection rate in ticks was 3.312.2%. The continuation of research in this direction remains relevant.

Characterization of the complete mitochondrial genomes of five hard ticks and phylogenetic implications.

Ticks are blood-sucking ectoparasites with significant medical and veterinary importance, capable of transmitting bacteria, protozoa, fungi, and viruses that cause a variety of human and animal diseases worldwide. In the present study, we sequenced the complete mitochondrial (mt) genomes of five hard tick species and analyzed features of their gene contents and genome organizations. The complete mt genomes of Haemaphysalis verticalis, H. flava, H. longicornis, Rhipicephalus sanguineus and Hyalomma asiaticum were 14855 bp, 14689 bp, 14693 bp, 14715 bp and 14722 bp in size, respectively. Their gene contents and arrangements are the same as those of most species of metastriate Ixodida, but distinct from species of genus Ixodes. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes with two different computational algorithms (Bayesian inference and maximum likelihood) revealed the monophylies of the genera Rhipicephalus, Ixodes and Amblyomma, however, rejected the monophyly of the genus Haemaphysalis. To our knowledge, this is the first report of the complete mt genome of H. verticalis. These datasets provide useful mtDNA markers for further studies of the identification and classification of hard ticks.