Seventy-eight entire mitochondrial genomes and nuclear rRNA genes provide insight into the phylogeny of the hard ticks, particularly the Haemaphysalis species, Africaniella transversale and Robertsicus elaphensis

Insights from entire mitochondrial genome sequences into the phylogeny of ticks of the genera Haemaphysalis and Archaeocroton with the elevation of the subgenus Alloceraea Schulze, 1919 back to the status of a genus

We used entire mitochondrial (mt) genome sequences (14.5-15 kbp) to resolve the phylogeny of the four main lineages of the Haematobothrion ticks: Alloceraea, Archaeocroton, Bothriocroton and Haemaphysalis. In our phylogenetic trees, Alloceraea was the sister to Archaeocroton sphenodonti, a tick of an archetypal reptile, the tuatara, from New Zealand, to the exclusion of the rest of the species of Haemaphysalis. The mt genomes of all four of the Alloceraea species that have been sequenced so far had a substantial insert, 132-312 bp, between the tRNA-Glu (E) gene and the nad1 gene in their mt genomes. This insert was not found in any of the other eight subgenera of Haemaphysalis. The mt genomes of 13 species of Haemaphysalis from NCBI GenBank were added to the most recent data set on Haemaphysalis and its close relatives to help resolve the phylogeny of Haemaphysalis, including five new subgenera of Haemaphysalis not previously considered by other authors: Allophysalis (structurally primitive), Aboimisalis (structurally primitive), Herpetobia (structurally intermediate), Ornithophysalis (structurally advanced) and Segalia (structurally advanced). We elevated Alloceraea Schulze, 1919 to the status of genus because Alloceraea Schulze, 1919 is phylogenetically distinct from the other subgenera of Haemaphysalis. Moreover, we propose that the subgenus Allophysalis is the sister to the rest of the Haemaphysalis (14 subgenera) and that the 'structurally primitive' subgenera Hoogstraal and Kim comprise early diverging lineages. Our matrices of the pairwise genetic difference (percent) of mt genomes and partial 16S rRNA sequences indicated that the mt genome sequence of Al. kitaokai (gb# OM368280) may not be Al. kitaokai Hoogstraal, 1969 but rather another species of Alloceraea. In a similar way, the mt genome sequence of H. (Herpetobia) nepalensis Hoogstraal, 1962 (gb# NC_064124) was only 2% genetically different to that of H. (Allophysalis) tibetensis Hoogstraal, 1965 (gb# OM368293): this indicates to us that they are the same species. Alloceraea cretacea may be better placed in a genus other than Alloceraea Schulze, 1919. Reptiles may have been the host to the most recent common ancestor of Archaeocroton and Alloceraea.

New insights into the systematics of the afrotropical Amblyomma marmoreum complex (Acari: Ixodidae) and the genome of a novel Rickettsia africae strain using morphological and metagenomic approaches

The Amblyomma marmoreum complex includes afrotropical species, such as Amblyomma sparsum, a three-host tick that parasitizes reptiles, birds, and mammals, and is a recognized vector of Ehrlichia ruminantium. However, the lack of morphological, genetic and ecological data on A. sparsum has caused considerable confusion in its identification. In this study, we used microscopy and metagenomic approaches to analyze A. sparsum ticks collected from a puff adder snake (Bitis arietans) in southwest Senegal (an endemic rickettsioses area) in order to supplement previous morphological descriptions, provide novel genomic data for the A. marmoreum complex, and describe the genome of a novel spotted fever group Rickettsia strain. Based on stereoscope and scanning electron microscopy (SEM) morphological evaluations, we provide high-quality images and new insights about punctation and enameling in the adult male of A. sparsum to facilitate identification for future studies. The metagenomic approach allowed us assembly the complete mitochondrial genome of A. sparsum, as well as the nearly entire chromosome and complete plasmid sequences of a novel Rickettsia africae strain. Phylogenomic analyses demonstrated a close relationship between A. sparsum and Amblyomma nuttalli for the first time and confirmed the position of A. sparsum within the A. marmoreum complex. Our results provide new insights into the systematics of A. sparsum and A. marmoreum complex, as well as the genetic diversity of R. africae in the Afrotropical region. Future studies should consider the possibility that A. sparsum may be a vector for R. africae.

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

The complete mitochondrial genome of cattle tick clade C reveals the genetic relationship within Rhipicephalus microplus complex

Cattle ticks (Rhipicephalus microplus) are important economic ectoparasites causing direct and indirect damage to cattle and leading to severe economic losses in cattle husbandry. It is common knowledge that R. microplus is a species complex including five clades; however, the relationships within the R. microplus complex remain unresolved. In the present study, we assembled the complete mitochondrial genome of clade C by next-generation sequencing and proved its correctness based on long PCR amplification. It was 15,004 bp in length and consisted of 13 protein genes, 22 transfer genes, and two ribosomal genes located in the two strains. There were two copies of the repeat region (pseudo-nad1 and tRNA-Glu). Data revealed that cox1, cox2, and cox3 genes were conserved within R. microplus with small genetic differences. Ka/Ks ratios suggested that 12 protein genes (excluding nad6) may be neutral selection. The genetic and phylogenetic analyses indicated that clade C was greatly close to clade B. Findings in the current study provided more data for the identification and differentiation of the R. microplus complex and made up for the lack of information about R. microplus clade C.

New insights into the molecular phylogeny, biogeographical history, and diversification of Amblyomma ticks (Acari: Ixodidae) based on mitogenomes and nuclear sequences

BACKGROUND

Amblyomma is the third most diversified genus of Ixodidae that is distributed across the Indomalayan, Afrotropical, Australasian (IAA), Nearctic and Neotropical biogeographic ecoregions, reaching in the Neotropic its highest diversity. There have been hints in previously published phylogenetic trees from mitochondrial genome, nuclear rRNA, from combinations of both and morphology that the Australasian Amblyomma or the Australasian Amblyomma plus the Amblyomma species from the southern cone of South America, might be sister-group to the Amblyomma of the rest of the world. However, a stable phylogenetic framework of Amblyomma for a better understanding of the biogeographic patterns underpinning its diversification is lacking.

METHODS

We used genomic techniques to sequence complete and nearly complete mitochondrial genomes -ca. 15 kbp- as well as the nuclear ribosomal cluster -ca. 8 kbp- for 17 Amblyomma ticks in order to study the phylogeny and biogeographic pattern of the genus Amblyomma, with particular emphasis on the Neotropical region. The new genomic information generated here together with genomic information available on 43 ticks (22 other Amblyomma species and 21 other hard ticks-as outgroup-) were used to perform probabilistic methods of phylogenetic and biogeographic inferences and time-tree estimation using biogeographic dates.

RESULTS

In the present paper, we present the strongest evidence yet that Australasian Amblyomma may indeed be the sister-group to the Amblyomma of the rest of the world (species that occur mainly in the Neotropical and Afrotropical zoogeographic regions). Our results showed that all Amblyomma subgenera (Cernyomma, Anastosiella, Xiphiastor, Adenopleura, Aponomma and Dermiomma) are not monophyletic, except for Walkeriana and Amblyomma. Likewise, our best biogeographic scenario supports the origin of Amblyomma and its posterior diversification in the southern hemisphere at 47.8 and 36.8 Mya, respectively. This diversification could be associated with the end of the connection of Australasia and Neotropical ecoregions by the Antarctic land bridge. Also, the biogeographic analyses let us see the colonization patterns of some neotropical Amblyomma species to the Nearctic.

CONCLUSIONS

We found strong evidence that the main theater of diversification of Amblyomma was the southern hemisphere, potentially driven by the Antarctic Bridge's intermittent connection in the late Eocene. In addition, the subgeneric classification of Amblyomma lacks evolutionary support. Future studies using denser taxonomic sampling may lead to new findings on the phylogenetic relationships and biogeographic history of Amblyomma genus.

Characterization of the complete mitochondrial genome and phylogenetic analyses of Haemaphysalis tibetensis Hoogstraal, 1965 (Acari: Ixodidae)

Ticks are specialized ectoparasites that feed on blood, causing physical harm to the host and facilitating pathogen transmission. The genus Haemaphysalis contains vectors for numerous infectious agents. These agents cause various diseases in humans and animals. Mitochondrial genome sequences serve as reliable molecular markers, forming a crucial basis for evolutionary analyses, studying species origins, and exploring molecular phylogeny. We extracted mitochondrial genome from the enriched mitochondria of Haemaphysalis tibetensis and obtained a 14,714-bp sequence. The mitochondrial genome consists of 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNAs (tRNAs), and two control regions. The nucleotide composition of H. tibetensis mitochondrial genome was 38.38 % for A, 9.61 % for G, 39.32 % for T, and 12.69 % for C. The A + T content of H. tibetensis mitochondrial genome was 77.7 %, significantly higher than the G + C content. The repeat units of H. tibetensis exhibited two identical repeat units of 33 bp in length, positioned downstream of nad1 and rrnL genes. Furthermore, phylogenetic analyses based on the 13 PCGs indicated that Haemaphysalis tibetensis (subgenus Allophysalis) formed a monophyletic clade with Haemaphysalis nepalensis (subgenus Herpetobia) and Haemaphysalis danieli (subgenus Allophysalis). Although the species Haemaphysalis inermis, Haemaphysalis kitaokai, Haemaphysalis kolonini, and Haemaphysalis colasbelcouri belong to the subgenus Alloceraea, which were morphologically primitive hemaphysalines just like H. tibetensis, these four tick species cannot form a single clade with H. tibetensis. In this study, the whole mitochondrial genome sequence of H. tibetensis from Tibet was obtained, which enriched the mitochondrial genome data of ticks and provided genetic markers to study the population heredity and molecular evolution of the genus Haemaphysalis.

First report of Anaplasma spp., Ehrlichia spp., and Rickettsia spp. in Amblyomma gervaisi ticks infesting monitor lizards (Varanus begalensis) of Pakistan

Ticks pose significant health risks to both wildlife and humans due to their role as vectors for various pathogens. In this study, we investigated tick infestation patterns, tick-associated pathogens, and genetic relationships within the tick species Amblyomma gervaisi, focusing on its prevalence in monitor lizards (Varanus bengalensis) across different districts in Pakistan. We examined 85 monitor lizards and identified an overall mean intensity of 19.59 ticks per infested lizard and an overall mean abundance of 11.98 ticks per examined lizard. All collected ticks (n = 1019) were morphologically identified as A. gervaisi, including 387 males, 258 females, 353 nymphs, and 21 larvae. The highest tick prevalence was observed in the Buner district, followed by Torghar and Shangla, with the lowest prevalence in Chitral. Lizard captures primarily occurred from May to October, correlating with the period of higher tick infestations. Molecular analysis was conducted on tick DNA, revealing genetic similarities among A. gervaisi ticks based on 16S rDNA and ITS2 sequences. Notably, we found the absence of A. gervaisi ITS2 sequences in the NCBI GenBank, highlighting a gap in existing genetic data. Moreover, our study identified the presence of pathogenic microorganisms, including Ehrlichia sp., Candidatus Ehrlichia dumleri, Anaplasma sp., Francisella sp., Rickettsia sp., and Coxiella sp., in these ticks. BLAST analysis revealed significant similarities between these pathogenic sequences and known strains, emphasizing the potential role of these ticks as vectors for zoonotic diseases. Phylogenetic analyses based on nuclear ITS2 and mitochondrial 16S rDNA genes illustrated the genetic relationships of A. gervaisi ticks from Pakistan with other Amblyomma species, providing insights into their evolutionary history. These findings contribute to our understanding of tick infestation patterns, and tick-borne pathogens in monitor lizards, which has implications for wildlife health, zoonotic disease transmission, and future conservation efforts. Further research in this area is crucial for a comprehensive assessment of the risks associated with tick-borne diseases in both wildlife and humans.

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.

A new subgenus, Australixodes n. subgen. (Acari: Ixodidae), for the kiwi tick, Ixodes anatis Chilton, 1904, and validation of the subgenus Coxixodes Schulze, 1941 with a phylogeny of 16 of the 22 subgenera of Ixodes Latreille, 1795 from entire mitochondrial genome sequences

A new subgenus, Australixodes n. subgen., is described for the kiwi tick, Ixodes anatis Chilton, 1904. The subgenus Coxixodes Schulze, 1941 is validated for the platypus tick, Ixodes (Coxixodes) ornithorhynchi Lucas, 1846, sister group of the subgenus Endopalpiger Schulze, 1935. A phylogeny from mitochondrial genomes of 16 of the 22 subgenera of Ixodes (32 spp.) indicates, for the first time, the relationships of the subgenera of Ixodes Latreille, 1795, the largest genus of ticks.

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.

Paradigms in tick evolution

The study of tick evolution may be classified into disciplines such as taxonomy and systematics, biogeography, evolution and development (evo-devo), ecology, and hematophagy. These disciplines overlap and impact each other to various extents. Advances in one field may lead to paradigm shifts in our understanding of tick evolution not apparent to other fields. The current study considers paradigm shifts that occurred, are in the process, or may occur in future for the disciplines that study tick evolution. Some disciplines have undergone significant changes, while others may still be developing their own paradigms. Integration of these various disciplines is essential to come to a holistic view of tick evolution; however, maturation of paradigms may be necessary before this vision can be attained.

The “southeastern Europe” lineage of the brown dog tick Rhipicephalus sanguineus (sensu lato) identified as Rhipicephalus rutilus Koch, 1844: Comparison with holotype and generation of mitogenome reference from Israel

The brown dog tick Rhipicephalus sanguineus (sensu lato) in the southeastern Mediterranean region and the Middle East is difficult to identify due to the presence of multiple mitochondrial DNA haplogroup lineages. The purpose of this study was to clarify the identity of the "southeastern Europe" lineage of this tick species complex. Our research shows that female ticks of the "southeastern Europe" lineage correspond to the morphology of R. rutilus Koch, 1844 as found in type-material at the Museum für Naturkunde Berlin in Germany. We characterised the complete mitogenomes of R. rutilus, R. turanicus Pomerantsev, 1940 and Rhipicephalus sanguineus (Latreille, 1806) in order to improve our understanding of the phylogenetic relationships among species within the R. sanguineus (sensu lato) complex. The material associated with the morphology of R. rutilus was previously labelled as the "southeastern Europe" lineage and found in Israel and Egypt, including Lower Egypt and the Nile Delta, where the original type-material was collected. Based on the morphology, genetic identity, and geographical distribution of the species, we conclude that the name R. rutilus is correctly linked to the "southeastern Europe" lineage of R. sanguineus (sensu lato).