Tick mitochondrial genomes: structural characteristics and phylogenetic implications

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.

The evolutionary characterization of Gamasida based on mitochondrial genes codon usage pattern

Mites belonging to the suborder Gamasida are species-rich and habitat-diverse, with a worldwide distribution. To adapt to the environment and obtain better living conditions, all species of the suborder Gamasida have been undergoing constant evolution. The complete mitochondrial genome (mitogenome) is an invaluable molecular marker for studying the origin of species, genetic differentiation between closely related species, and between intraspecific groups. In some species of the suborder Gamasida, mitochondrial tRNA genes are truncated and carried unstable genetic information. This study presents a comparative analysis of codon usage pattern and preference of 13 protein-coding genes of 24 species in 17 genera and 10 families of the suborder Gamasida. Results showed that have an obvious AT preference (0.664-0.829) for codon usage in the suborder Gamasida. Most of the optimal and high-frequency codons also end in A/T. The degree of natural selection varies between the same protein-coding genes of different gamasid mites or among different protein-coding genes within the same gamasid mites. Base and codon usage pattern and preference are very similar between the same species and genus, namely the closer species, the more similar their bases and codons usage patterns and preference are. T bases and C bases were the preference bases for codon usage of 24 species in the suborder Gamasida. Evolution of the suborder Gamasida was dominated by natural selection (64.1%). This study provides the first comprehensive analysis of codon usage in the suborder Gamasida, which will greatly improve our understanding of codon usage patterns and preference, genetics, and evolution of the suborder Gamasida. It will help to evaluate the degree of molecular adaptation in the suborder Gamasida and to further explore evolutionary features of the suborder Gamasida.

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.

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.

Epidemiological survey of tick-borne pathogens in dogs in Anhui Province, China

Tick-borne diseases have continued to increase worldwide in both developing and many developed countries due to the widespread of different tick species and tick's adaptability to different climatic weather. In order to investigate the prevalence of the tick-borne pathogens, EDTA-anticoagulated whole blood samples were aseptically collected from 765 pet dogs in twenty veterinary clinics located in sixteen prefecture-level cities in Anhui Province, China, and the samples were examined and analyzed for tick-borne pathogens using both microscopy and PCR. Our result analysis revealed 17(2.22%) positive samples to Babesia spp and 4(0.52%) positive samples to Hepatozoon spp, of which case of co-infection was recorded in Lu'An and Chuzhou. The BLAST analysis results of the 18S rRNA gene revealed that the dogs were infected with Babesia gibsoni and Hepatozoon canis. All samples were negative for Anaplasma spp., Ehrlichia spp., and Rickettsia spp. This is the first molecular report of B. gibsoni and H. canis in dogs in Anhui, China.

Complete mitochondrial genome of Ctenophthalmus quadratus and Stenischia humilis in China provides insights into fleas phylogeny

Fleas (Order Siphonaptera) are common blood-feeding ectoparasites, which have important economic significance. Limited mitochondrial genome information has impeded the study of flea biology, population genetics and phylogenetics. The Ctenophthalmus quadratus and Stenischia humilis complete mt genomes are described in this study. The samples were collected from Jianchuan, Yunnan plague foci, China. The mt genomes of C. quadratus and S. humilis were 15,938 bp and 15,617 bp, respectively. The gene arrangement of mt genome was consistent with that of other fleas, which include 22 tRNA genes, 13 protein-coding genes, and two rRNA genes, with a total of 37 genes. The relationship between C. quadratus and S. humilis in fleas was inferred by phylogenetic analysis of mt genome sequence datasets. Phylogenetic analyzes showed that the C. quadratus and S. humilis belonged to different species in the same family, and were closely related to Hystrichopsylla weida qinlingensis in the same family; and revealed that the family Hystrichopsyllidae is paraphyletic, supporting the monophyly of the order Siphonaptera. This study decodes the complete mt genomes of the C. quadratus and S. humilis for the first time. The results demonstrate that the C. quadratus and S. humilis are distinct species, and fleas are monophyletic. Analysis of mt genome provides novel molecular data for further studying the phylogeny and evolution of fleas.

The complete mitochondrial genome of Eulaelaps huzhuensis (Mesostigmata: Haemogamasidae)

Some mites of the family Haemogamasidae can transmit a variety of zoonotic diseases and have important public health and safety implications. Currently, however, little attention has been paid to molecular data of Haemogamasidae species, limiting our understanding of their evolutionary and phylogenetic relationships. In this study, the complete mitochondrial genome of Eulaelaps huzhuensis was determined for the first time, and its genomic information was analyzed in detail. The mitochondrial genome of E. huzhuensis is 14,872 bp in length with 37 genes and two control regions. The base composition showed a distinct AT preference. Twelve protein-coding genes have a typical ATN as the start codon, and three protein-coding genes have incomplete stop codons. During the folding of tRNA genes, a total of 30 mismatches occurred, and three tRNA genes had an atypical cloverleaf secondary structure. The order of the E. huzhuensis mitochondrial genome arrangement is a new type of rearrangement in Mesostigmata. The phylogenetic analysis confirmed that the family Haemogamasidae is a monophyletic branch and does not belong to a subfamily of the Laelapidae. Our results lay the foundation for subsequent studies on the phylogeny and evolutionary history of the family Haemogamasidae.

Canine piroplasmids: Molecular detection and laboratory characterization in dogs from Brasilia, Brazil, with the first molecular evidence of dog exposure to a novel opossum-associated Babesia sp

Canine piroplasmid infections can be caused by Babesia spp., Theileria spp. and Rangelia vitalii. In Brazil, canine babesiosis caused by Babesia vogeli is endemic and reported throughout the country. On the other hand, Rangeliosis caused by R. vitalii has only been described so far in the South and Southeast regions. Despite that, studies analyzing the laboratory and molecular characterization of these hemoprotozoa are still scarce. To investigate the occurrence, the laboratory features, the molecular characterization, and the diversity of piroplasmids from Midwestern Brazil, a survey was performed using blood samples obtained from 276 domestic dogs from Brasília, Federal District, Midwestern Brazil. A broad-range quantitative PCR (qPCR) targeting the mitochondrial large subunit ribosomal DNA (LSU4) was used to detect piroplasmid DNA. The overall molecular occurrence of piroplasmids was 11.2% (31/276), with 9.7% (27/276) of the sequences identified as Babesia vogeli (98-100% identity to B. vogeli isolate from the USA). Based on a partial 18S rRNA sequence pairwise alignment (-250 bp), 1.4% (4/276) of the sequences showed only 76.8% identity with B. vogeli but 100% identity with opossum-associated Babesia sp. (MW290046-53). These findings suggest the exposure of dogs from Brazil to a recently described Babesia sp. isolated from white-eared opossum. None of the analyzed dogs was positive for Theileria spp. or R. vitalii. Subsequently, all positive sequences were submitted to three additional PCR assays based on the 18S rRNA, cox-1, and cytb genes, aiming at performing a haplotype network analysis. Haplotype network using cox-1 sequences showed the presence of six different haplotypes of B. vogeli; one of them was shared with isolates from Brazil, the USA, and India. When including animals co-infected with other vector-borne diseases, piroplasmid-positive dogs had 2.3 times higher chance of having thrombocytopenia than the negative ones. The molecular results demonstrated that the compared Babesia vogeli sequences showed a low variability as well as evidence of exposure to a putative novel opossum-associated Babesia sp. in dogs from Midwestern Brazil.

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.

Phylogenetic relationships of the Amblyomma cajennense complex (Acari: Ixodidae) at mitogenomic resolution

The genus Amblyomma is the third most diverse in the number of species within the Ixodidae, with practically half of its species distributed in the Americas, though there are also species occurring in Africa, Asia, and Australia. Within the genus, there are several species complexes with veterinary and public health importance. The Amblyomma cajennense complex, in the Americas, is represented by six species with a wide distribution, from Texas to northern Argentina. We combined two sequencing techniques to generate complete mitogenomes of species belonging to the Amblyomma cajennense complex: genome skimming and long-range PCRs sequencing methods. Thus, we generated seven new mitochondrial genomes for all species of the Amblyomma cajennense complex, except for Amblyomma interandinum. Genetic distances between the mitogenomes corroborate the clear differentiation between the five species of the Amblyomma cajennense complex. The phylogenetic relationships of these species had previously been evaluated by combining partial nuclear and mitochondrial genes and here these relationships are corroborated with a more robust framework of data, which demonstrates that the conjunction of mitochondrial and nuclear partial genes can resolve close relationships when entire genes or genomes are unavailable. The gene order, structure, composition, and length are stable across these mitogenomes, and they share the general characteristics of Metastriata. Future studies should increase the number of available mitogenomes for this genus, especially for those species from the Indo-Pacific region and Africa, by means of a better understanding of their relationships and evolutionary process.

Mitogenomic and Phylogenetic Analysis of the Entomopathogenic Fungus Ophiocordyceps lanpingensis and Comparative Analysis with Other Ophiocordyceps Species

Ophiocordyceps lanpingensis (O. lanpingensis) belongs to the genus Ophiocordyceps, which is often found in Yunnan Province, China. This species is pharmacologically important for the treatment of renal disorders induced by oxidative stress and an inadequate immune response. In the present study, the mitogenome of O. lanpingensis was determined to be a circular molecule 117,560 bp in length, and to have 31% G + C content and 69% A + T content. This mitogenome comprised 82% of the whole genome that codes for significant genes. The protein-coding regions of the O. lanpingensis mitogenome, containing 24 protein-coding genes, were associated with respiratory chain complexes, such as 3 ATP-synthase complex F0 subunits (atp6, atp8, and atp9), 2 complex IV subunits/cytochrome c oxidases (cox2 and cox3), 1 complex III subunit (cob), 4 electron transport complex I subunits/NADH dehydrogenase complex subunits (nad1, nad4, nad5, and nad6), 2 ribosomal RNAs (rns, rnl), and 11 hypothetical/predicted proteins, i.e., orf609, orf495, orf815, orf47, orf150, orf147, orf292, orf127, orf349, orf452, and orf100. It was noted that all genes were positioned on the same strand. Further, 13 mitochondrial genes with respiratory chain complexes, which presented maximum similarity with other fungal species of Ophiocordyceps, were investigated. O. lanpingensis was compared with previously sequenced species within Ophiocordycepitaceae. Comparative analysis indicated that O. lanpingensis was more closely related to O. sinensis, which is one of the most remarkable and expensive herbs due to its limited availability and the fact that it is difficult to culture. Therefore, O. lanpingensis is an important medicinal resource that can be effectively used for medicinal purposes. More extensive metabolomics research is recommended for O. lanpingensis.

Comparative analysis of the mitochondrial genome of Dermacentor steini from different regions in China

Ticks are a group of blood-sucking ectoparasites that play an important role in human health and livestock production development as vectors of zoonotic diseases. The phylogenetic tree of single genes cannot accurately reflect the true kinship between species. Based on the complete mitochondrial genome analysis one can help to elucidate the phylogenetic relationships among species. In this study, the complete mitochondrial genome of Dermacentor steini (isolate Longyan) was sequenced and compared with the mitochondrial genes of 3 other Chinese isolates (Nanchang, Jinhua and Yingtan). In Dermacentor steini 4 isolates had identical or similar mitochondrial genome lengths and an overall variation of 0.76% between sequences. All nucleotide compositions showed a distinct AT preference. The most common initiation and stop codons were ATG and TAA, respectively. Fewer base mismatches were found in the tRNA gene of D. steini (isolate Longyan), and the vicinity of the control region and tRNA gene was a hot rearrangement region of the genus Dermacentor. Maximum likelihood trees and Bayesian trees indicate that D. steini is most closely related to Dermacentor auratus. The results enrich the mitochondrial genomic data of species in the genus Dermacentor and provide novel insights for further studies on the phylogeographic classification and molecular evolution of ticks.

Mitochondrial Genomes of the American Dog Tick (Dermacentor variabilis) Isolated from Horses in the Midwestern United States

Here, we report two complete and three partial mitochondrial genome sequences of Dermacentor variabilis specimens collected from horses in the United States. The complete genomes are 14,837 bp long and contain 13 protein-coding genes, 2 rRNA genes, and 22 tRNA genes. The sequences have been deposited under GenBank accession numbers ON052120 to ON052124.

Complete Mitogenomes of Ticks Ixodes acutitarsus and Ixodes ovatus Parasitizing Giant Panda: Deep Insights into the Comparative Mitogenomic and Phylogenetic Relationship of Ixodidae Species

Ticks rank second in the world as vectors of disease. Tick infestation is one of the factors threatening the health and survival of giant pandas. Here, we describe the mitogenomes of Ixodes acutitarsus and Ixodes ovatus parasitizing giant pandas, and perform comparative and phylogenetic genomic analyses on the newly sequenced and other available mitogenomes of hard ticks. All six newly determined mitogenomes contain a typical gene component and share an ancient Arthropoda gene arrangement pattern. Our study suggests that I. ovatus is a species complex with high genetic divergence, indicating that different clades of I. ovatus represent distinct species. Comparative mitogenomic analyses show that the average A + T content of Ixodidae mitogenomes is 78.08%, their GC-skews are strongly negative, while AT-skews fluctuate around 0. A large number of microsatellites are detected in Ixodidae mitogenomes, and the main microsatellite motifs are mononucleotide A and trinucleotide AAT. We summarize five gene arrangement types, and identify the trnY-COX1-trnS1-COX2-trnK-ATP8-ATP6-COX3-trnG fragment is the most conserved region, whereas the region near the control region is the rearrangement hotspot in Ixodidae mitogenomes. The phylogenetic trees based on 15 genes provide a very convincing relationship (Ixodes + (Robertsicus + ((Bothriocroton + Haemaphysalis) + (Amblyomma + (Dermacentor + (Rhipicentor + (Hyalomma + Rhipicephalus))))))) with very strong supports. Remarkably, Archaeocroton sphenodonti is embedded in the Haemaphysalis clade with strong supports, resulting in paraphyly of the Haemaphysalis genus, so in-depth morphological and molecular studies are essential to determine the taxonomic status of A. sphenodonti and its closely related species. Our results provide new insights into the molecular phylogeny and evolution of hard ticks, as well as basic data for population genetics assessment and efficient surveillance and control for the giant panda-infesting ticks.

Analysis of the complete mitochondrial genomes of Dermacentor albipictus suggests a species complex

Dermacentor albipictus is a one-host tick broadly distributed across North America. There are two easily recognizable color variants - ornate and inornate/brown - that have been taxonomically synonymized. Based on mt-cox1 and mt-16S data, there is also evidence for two genetic lineages which do not match the color variants. We present for the first time the complete mitochondrial genomes of the two color variants of D. albipictus including representatives of each lineage. The AT-rich genomes are 14,822 bp - 14,865 bp in length and contain 13 protein coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes, arranged in the conserved type 3 metastriate mitochondrial genome order. The overall differences were 10.66% between the mitochondrial genomes of D. albipictus ornate variant lineage 1 and lineage 2, 10.51% between lineage 1 and inornate/brown variant and 5.87% between lineage 2 and inornate/brown variant. The inornate/brown variant did not form a separate lineage and all inornate isolates were found to belong to lineage 2. Ornate variant isolates occurred in both lineage 1 and 2. The high divergence of the mitochondrial genome suggests that D. albipictus may represent a species complex. Other barcoding genes that may help capture the genetic differences between color and lineage variants include nad1, nad2, nad5, cox1 and atp8 loci. The mtDNA data generated in this study are available in GenBank (Accession numbers: OM678457 - OM678459 and ON032564 - ON032573) for future studies on tick taxonomy, phylogenetics and molecular epidemiology.

Tick-borne diseases in Egypt: A one health perspective

Background

Ticks are important arthropod vectors that transmit pathogens to humans and animals. Owing to favourable climatic and environmental conditions, along with animal importation from neighbouring countries, ticks and tick-borne diseases (TBDs) are widespread in Egyptian localities. Here, we review the current knowledge on the epidemiology of TBDs in Egypt in light of the One Health paradigm.

Methods and results

Five scientific databases, including "Web of Science", "Scopus", "PubMed", "Science Direct", and "Google Scholar", were searched for articles describing TBDs in Egypt. A total of 18 TBDs have been reported in humans and animals, including three protozoal diseases (babesiosis, theileriosis, and hepatozoonosis), 12 bacterial diseases (anaplasmosis, ehrlichiosis, Lyme borreliosis, bovine borreliosis, tick-borne relapsing fever, Mediterranean spotted fever, African tick-borne fever, lymphangitis-associated rickettsiosis, bartonellosis, tularaemia, Q fever, and aegyptianellosis), and three viral diseases (Crimean-Congo haemorrhagic fever, Alkhurma haemorrhagic fever, and Lumpy skin disease).

Conclusions

Despite the circulation of zoonotic tick-borne pathogens among livestock and tick vectors, human infections have been overlooked and are potentially limited to infer the actual communicable disease burden. Therefore, facility-based surveillance of TBDs, combined with capacity building for laboratory diagnostics in healthcare facilities, is urgently required to improve diagnosis and inform policy-making in disease prevention. Additionally, collaboration between expert researchers from various disciplines (physicians, biologists, acarologists, and veterinarians) is required to develop advanced research projects to control ticks and TBDs. Considering that domestic livestock is integral to many Egyptian households, comprehensive epidemiological studies on TBDs should assess all disease contributors, including vertebrate hosts (animals, humans, and rodents) and ticks in the same ecological region, for better assessment of disease burden. Additionally, upscaling of border inspections of imported animals is required to stop crossover movements of ticks and TBDs.

The diversity and evolutionary relationships of ticks and tick-borne bacteria collected in China

Background

Ticks (order Ixodida) are ectoparasites, vectors and reservoirs of many infectious agents affecting humans and domestic animals. However, the lack of information on tick genomic diversity leaves significant gaps in the understanding of the evolution of ticks and associated bacteria.

Results

We collected > 20,000 contemporary and historical (up to 60 years of preservation) tick samples representing a wide range of tick biodiversity across diverse geographic regions in China. Metagenomic sequencing was performed on individual ticks to obtain the complete or near-complete mitochondrial (mt) genome sequences from 46 tick species, among which mitochondrial genomes of 23 species were recovered for the first time. These new mt genomes data greatly expanded the diversity of many tick groups and revealed five cryptic species. Utilizing the same metagenomic sequence data we identified divergent and abundant bacteria in Haemaphysalis, Ixodes, Dermacentor and Carios ticks, including nine species of pathogenetic bacteria and potentially new species within the genus Borrelia. We also used these data to explore the evolutionary relationship between ticks and their associated bacteria, revealing a pattern of long-term co-divergence relationship between ticks and Rickettsia and Coxiella bacteria.

Conclusions

In sum, our study provides important new information on the genetic diversity of ticks based on an analysis of mitochondrial DNA as well as on the prevalence of tick-borne pathogens in China. It also sheds new light on the long-term evolutionary and ecological relationships between ticks and their associated bacteria.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-022-05485-3.

Complete mitochondrial genomes of Thyreophagus entomophagus and Acarus siro (Sarcoptiformes: Astigmatina) provide insight into mitogenome features, evolution, and phylogeny among Acaroidea mites

Mites from the Acaroidea (Sarcoptiformes: Astigmatina) are important pests of various stored products, posing potential threats to preserved foods. In addition, mites can cause allergic diseases. Complete mitochondrial genomes (mitogenomes) are valuable resources for different research fields, including comparative genomics, molecular evolutionary analysis, and phylogenetic inference. We sequenced and annotated the complete mitogenomes of Thyreophagus entomophagus and Acarus siro. A comparative analysis was made between mitogenomic sequences from 10 species representing nine genera within Acaroidea. The mitogenomes of T. entomophagus and A. siro contained 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region. In Acaroidea species, mitogenomes have highly conserved gene size and order, and codon usage. Among Acaroidea mites, most PCGs were found to be under purifying selection, implying that most PCGs might have evolved slowly. Our findings showed that nad4 evolved most rapidly, whereas cox1 and cox3 evolved most slowly. The evolutionary rates of Acaroidea vary considerably across families. In addition, selection analyses were also performed in 23 astigmatid mite species, and the evolutionary rate of the same genes in different superfamilies exhibited large differences. Phylogenetic results are mostly consistent with those identified by previous phylogenetic studies on astigmatid mites. The monophyly of Acaroidea was rejected, and the Suidasiidae and Lardoglyphidae appeared to deviate from the Acaroidea branch. Our research proposed a review of the current Acaroidea classification system.

Comparative mitogenomics elucidates the population genetic structure of Amblyomma testudinarium in Japan and a closely related Amblyomma species in Myanmar

Ticks are the second most important vector capable of transmitting diseases affecting the health of both humans and animals. Amblyomma testudinarium Koch 1844 (Acari: Ixodidae), is a hard tick species having a wide geographic distribution in Asia. In this study, we analyzed the composition of A. testudinarium whole mitogenomes from various geographical regions in Japan and investigated the population structure, demographic patterns, and phylogeographic relationship with other ixodid species. In addition, we characterized a potentially novel tick species closely related to A. testudinarium from Myanmar. Phylogeographic inference and evolutionary dynamics based on the 15 mitochondrial coding genes supported that A. testudinarium population in Japan is resolved into a star‐like haplogroup and suggested a distinct population structure of A. testudinarium from Amami island in Kyushu region. Correlation analysis using Mantel test statistics showed that no significant correlation was observed between the genetic and geographic distances calculated between the A. testudinarium population from different localities in Japan. Finally, demographic analyses, including mismatch analysis and Tajima’s D test, suggested a possibility of recent population expansion occurred within Japanese haplogroup after a bottleneck event. Although A. testudinarium has been considered widespread and common in East and Southeast Asia, the current study suggested that potentially several cryptic Amblyomma spp. closely related to A. testudinarium are present in Asia.

Characterization of the complete mitochondrial genome of Haemaphysalis (Alloceraea) kolonini (Ixodidae) and its phylogenetic implications

Ticks transmit diverse pathogens that cause human and animal diseases, leading to an increasing number of new challenges around the world. Genomic data research could help advance our learning of phylogenetic analysis and molecular evolution. Mitochondrial genome DNA has been helpful in illustrating the phylogenetic analysis of eukaryotes containing ticks. In this research, we sequenced and assembled the circular complete mitogenome information of Haemaphysalis kolonini. The 14,948-bp mitogenome consists of 37 genes which included 13 genes for protein-coding, two genes for ribosomal RNA, 22 genes for transfer RNA, and two control regions (D-loops). Overall, the composition and arrangement of genes were compared with Haemaphysalis ticks previously recorded in Genbank. The phylogenetic tree based on Maximum likelihood (ML) and Bayesian inference (BI) computational algorithms showed that H. kolonini has a close relationship with Haemaphysalis inermis. The complete mitogenome data provide a preferable perception to the phylogenetic relationship than the single-gene data analysis. To our knowledge, this is the first research exploring the complete mitogenome for the species H. kolonini. Our results provide new insights for further research on the evolution, population genetics, systematics, and molecular ecology of ticks.

Comparative hologenomics of two Ixodes scapularis tick populations in New Jersey

Tick-borne diseases, such as those transmitted by the blacklegged tick Ixodes scapularis, are a significant and growing public health problem in the US. There is mounting evidence that co-occurring non-pathogenic microbes can also impact tick-borne disease transmission. Shotgun metagenome sequencing enables sampling of the complete tick hologenome—the collective genomes of the tick and all of the microbial species contained therein, whether pathogenic, commensal or symbiotic. This approach simultaneously uncovers taxonomic composition and allows the detection of intraspecific genetic variation, making it a useful tool to compare spatial differences across tick populations. We evaluated this approach by comparing hologenome data from two tick samples (N = 6 ticks per location) collected at a relatively fine spatial scale, approximately 23 km apart, within a single US county. Several intriguing variants in the data between the two sites were detected, including polymorphisms in both in the tick’s own mitochondrial DNA and that of a rickettsial endosymbiont. The two samples were broadly similar in terms of the microbial species present, including multiple known tick-borne pathogens (Borrelia burgdorferi, Babesia microti, and Anaplasma phagocytophilum), filarial nematodes, and Wolbachia and Babesia species. We assembled the complete genome of the rickettsial endosymbiont (most likely Rickettsia buchneri) from both populations. Our results provide further evidence for the use of shotgun metagenome sequencing as a tool to compare tick hologenomes and differentiate tick populations across localized spatial scales.

Complete mitochondrial genome of Asian longhorned tick, Haemaphysalis longicornis, Neumann, 1901 (Acari: Ixodida: Ixodidae) identified in the United States

Haemaphysalis longicornis (Ixodida: Ixodidae), the Asian longhorned tick, which is native to temperate East Asia, has been recently detected in the northeastern region of the United States, drawing concerns about its potential impact on the US animal and public health sectors. Knowledge about the genetic features of H. longicornis found in the US is limited. Therefore, we sequenced the complete mitochondrial genome (mt-genome) from two H. longicornis ticks recently collected in the State of New York, USA, in 2020. These ticks were morphologically identified and tested for tick-borne pathogens at the Connecticut Veterinary Medical Diagnostic Laboratory (Storrs, CT). The mt-genome was 14,694 bp in length and encoded 37 genes, including 13 protein-coding genes, 22 transfer RNAs, and two ribosomal RNAs. Phylogenetic analysis showed that the mt-genome clustered with those of other H. longicornis identified in China. The mt-genome sequence was 99.7% identical to a H. longicornis mt-genome (GenBank: MK439888) collected in China. The cox1 gene haplotype in these ticks belonged to the H1 type, which is the dominant haplotype present in central NJ and Staten Island, NY. The complete mt-genome data are needed to provide insights into genetic changes and phylogenetic studies of H. longicornis ticks.

Pitfalls in Tick and Tick-Borne Pathogens Research, Some Recommendations and a Call for Data Sharing

An understanding of the relationships of ticks and tick-borne pathogens can only be achieved by integrating data from multiple studies. The publication of raw material is a necessary step for wide-area meta-analyses and study design, data collection and reporting require harmonization. This is an opinion paper, not a consensus position, and is open to debate. This work reflects our view about how data should be communicated in mainstream journals. We indicate rules that should be observed in recording weather data, to avoid serendipitous correlations between the density of ticks and climate variables and recommend the inclusion of raw data in reports. We stress the need for standardized methods to collect ticks that cannot be obtained by standard flagging. The reporting of infection rates of pathogens in ticks should avoid conclusions based on pure molecular findings in feeding ticks. Studies demonstrating the vectorial capacity of ticks should not be supported only by molecular surveys of feeding ticks. Vacuous conclusions about vectorial or reservoir status based solely on the finding of genomic material of a pathogen should be discouraged. We stress that phylogenetic studies based on random selection of sequences from GenBank are unsuitable. We firmly support the development of a dedicated server of curated sequences of ticks and pathogens as a standard for future studies.

Molecular investigation of bacterial and protozoal pathogens in ticks collected from different hosts in Turkey

Background

The emergence of tick-borne disease is increasing because of the effects of the temperature rise driven by global warming. In Turkey, 19 pathogens transmitted by ticks to humans and animals have been reported. Based on this, this study aimed to investigate tick-borne pathogens including Hepatozoon spp., Theileria spp., Babesia spp., Anaplasma spp., Borrelia spp., and Bartonella spp. in tick samples (n = 110) collected from different hosts (dogs, cats, cattle, goats, sheep, and turtles) by molecular methods.

Methods

To meet this objective, ticks were identified morphologically at the genus level by microscopy; after DNA isolation, each tick sample was identified at the species level using the molecular method. Involved pathogens were then investigated by PCR method.

Results

Seven different tick species were identified including Rhipicephalus sanguineus, R. turanicus, R. bursa, Hyalomma marginatum, H. anatolicum, H. aegyptium, and Haemaphysalis erinacei. Among the analyzed ticks, Hepatozoon spp., Theileria spp., Babesia spp., and Anaplasma spp. were detected at rates of 6.36%, 16.3%, 1.81%, and 6.36%, respectively while Borrelia spp. and Bartonella spp. were not detected. Hepatozoon spp. was detected in R. sanguineus ticks while Theileria spp., Babesia spp., and Anaplasma spp. were detected in R. turanicus and H. marginatum. According to the results of sequence analyses applied for pathogen positive samples, Hepatozoon canis, Theileria ovis, Babesia caballi, and Anaplasma ovis were identified.

Conclusion

Theileria ovis and Anaplasma ovis were detected for the first time to our knowledge in H. marginatum and R. turanicus collected from Turkey, respectively. Also, B. caballi was detected for the first time to our knowledge in ticks in Turkey.

Complete mitochondrial genome characterization and phylogenetic analyses of the main vector of Crimean-Congo haemorrhagic fever virus: Hyalomma marginatum Koch, 1844

The Mediterranean tick, Hyalomma marginatum, is the most important vector of Crimean-Congo haemorrhagic fever virus and several pathogens that cause animal and human diseases and economic losses to livestock production. Given the medical and veterinary importance of this tick species, we sequenced and characterized its mitochondrial genome (mitogenome) for the first time. We designed two new primer sets and combined long-range PCR with next generation sequencing to generate complete mitogenomes with deep coverage from 10 H. marginatum adults. The mitogenomes contained 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal subunits, two control regions, and three tick-box motifs. The nucleotide composition of the H. marginatum mitogenomes were A+T biased (79.76%) and exhibited negative AT- and GC- skews across most PCGs. All PCGs were initiated by ATK codons and two truncated termination codons were seen in the COX2 and COX3 genes. All tRNAs exhibited typical cloverleaf structures, except for tRNA^Cys and tRNA^Ser1. A total of 62 polymorphic sites defined ten unique haplotypes. Phylogenetic analyses based on the 13 PCGs of 56 tick species revealed that four Hyalomma species (H. marginatum, H. asiaticum, H. rufipes, and H. truncatum) formed a monophyletic clade with strong support. The results of this study provide a comprehensive resource for further studies on the systematics, population genetics, molecular epidemiology, and evolution of ticks.

The complete mitochondrial genome sequences of the cat flea Ctenocephalides felis felis (Siphonaptera: Pulicidae) support the hypothesis that C. felis isolates from China and USA were the same C. f. felis subspecies

The cat flea Ctenocephalides felis (Siphonaptera: Pulicidae) is the most important ectoparasite in cats and dogs worldwide. Over the years, there has been much dispute regarding the taxonomic and systematic status of C. felis. Mitochondrial (mt) genome sequences are useful genetic markers for the identification and differentiation of ectoparasites, but the mt genome of C. felis and its subspecies has not yet been entirely characterized. In the present study, the entire mt genome of C. f. felis from China was sequenced and compared with that of C. felis from the USA. Both contain 37 genes and a long non-coding region of >6 kbp. The molecular identity between the Chinese and American isolates was 99%, except for the non-coding region. The protein-coding genes showed differences at both the nucleotide (1.2%) and amino acid (1%) levels. Interestingly, the cox1 gene of the Chinese isolate had an unusual putative start codon (TTT). Taken together, our analyses strongly support the hypothesis that C. felis isolates from China and the USA were the same C. f. felis subspecies. The mt genome sequence of the C. f. felis China isolate presented in this study provides useful molecular markers to further address the taxonomy and systematics of C. felis.

Genetic diversity analysis of Dermacentor nuttalli within Inner Mongolia, China

Background

Ticks (Arthropoda, Ixodida), after mosquitoes, are the second most prevalent vector of infectious diseases. They are responsible for spreading a multitude of pathogens and threatening the health and welfare of animals and human beings. However, given the history of tick-borne pathogen infections in the Inner Mongolia Autonomous Region of China, surprisingly, neither the genetic diversity nor the spatial distribution of haplotypes within ticks has been studied.

Methods

We characterized the haplotype distribution of Dermacentor nuttalli in four main pastoral areas of the Inner Mongolia Autonomous Region, by sampling 109 individuals (recovered from sheep) in April–August 2019. The 16S rRNA gene, cytochrome c oxidase subunit I (COI), and the internal transcribed spacer 2 region (ITS2) were amplified and sequenced from extracted DNA.

Results

Twenty-six haplotypes were identified using 16S rRNA sequences, 57 haplotypes were identified with COI sequences, and 75 haplotypes were identified with ITS2 sequences. Among the three genes, total haplotype diversity was greater than 0.7, while total nucleotide diversity was greater than 0.06. Neutrality tests revealed a significantly negative Tajima’s D result, while Fu's Fs was not significantly positive. Fixation index values (F_ST) indicated that the degree of genetic differentiation among some sampled populations was small, while for others it was moderate. Analysis of molecular variance (AMOVA) revealed that the variation within populations was greater than that among populations. The mismatch analysis of D. nuttalli exhibited double peaks.

Conclusion

The genetic diversity of D. nuttalli populations in our region can likely adapt to different geographical environments, thereby leading to genetic diversity, and creating genetic differentiation among different populations. However, genetic differentiation is cryptic and does not form a pedigree geographical structure.

Proteases and pseudoproteases in parasitic arthropods of clinical importance

Parasitic arthropods feed on blood or skin tissue and share comparable repertoires of proteases involved in haematophagy, digestion, egg development and immunity. While proteolytically active proteases of multiple classes dominate, an increasing number of pseudoproteases have been discovered that have no proteolytic function but are pharmacologically active biomolecules, evolved to carry out alternative functions as regulatory, antihaemostatic, anti-inflammatory or immunomodulatory compounds. In this review, we provide an overview of proteases and pseudoproteases from clinically important arthropod parasites. Many of these act in central biological pathways of parasite survival and host-parasite interaction and may be potential targets for therapeutic interventions.

An Assessment of the Molecular Diversity of Ticks and Tick-Borne Microorganisms of Small Ruminants in Pakistan

This study investigated ticks and tick-borne microorganisms of small ruminants from five districts of the Federally Administered Tribal Area (FATA) of Pakistan. Morphological (n = 104) and molecular (n = 54) characterization of the ticks revealed the presence of six ixodid ticks: Rhipicephalus (Rh.) haemaphysaloides, Rh. microplus, Rh. turanicus, Haemaphysalis (Hs.) punctata, Hs. sulcata and Hyalomma anatolicum. Phylogenetic analyses of nucleotide sequence data for two mitochondrial (16S and cytochrome c oxidase 1) and one nuclear (second internal transcribed spacer) DNA regions provided strong support for the grouping of the six tick species identified in this study. Microfluidic real-time PCR, employing multiple pre-validated nuclear and mitochondrial genetic markers, detected 11 potential pathogens and endosymbionts in 72.2% of the ticks (n = 54) tested. Rickettsia (R.) massiliae was the most common pathogen found (42.6% of ticks) followed by Theileria spp. (33.3%), Anaplasma (A.) ovis and R. slovaca (25.9% each). Anaplasma centrale, A. marginale, Ehrlichia spp., R. aeschlimannii, R. conorii and endosymbionts (Francisella- and Coxiella-like) were detected at much lower rates (1.9–22.2%) in ticks. Ticks from goats (83.9%) carried significantly higher microorganisms than those from sheep (56.5%). This study demonstrates that ticks of small ruminants from the FATA are carrying multiple microorganisms of veterinary and medical health significance and provides the basis for future investigations of ticks and tick-borne diseases of animals and humans in this and neighboring regions.

Characterization of the complete mitochondrial genome of Amblyomma ovale, comparative analyses and phylogenetic considerations

In this study, we sequenced two complete mitochondrial genomes of Amblyomma ovale, a tick of public health importance. Sequencing two distinct individuals, the resulting mitochondrial genomes were 14,756 and 14,760 bp in length and maintained the same gene order previously reported in Amblyomma. These were combined with RNA-seq derived mitochondrial sequences from three additional species, Amblyomma aureolatum, Amblyomma maculatum, and Amblyomma moreliae, to carry out mitogenome comparative and evolutionary analyses against all previously published tick mitochondrial genomes. We described a derivative genome rearrangement that isolates Ixodes from the remaining Ixodidae and consists of both a reverse translocation as well as an event of Tandem Duplication Random Loss. Genetic distance analyses indicated that cox2, nd1, nd5, and 16S are good candidates for future population studies in A. ovale. The phylogenetic analyses corroborated the utility of complete mitochondrial genomes as phylogenetic markers within the group. This study further supplements the genome information available for Amblyomma and facilitates future evolutionary and population genetic studies within the genus.

Leishmania Mitochondrial Genomes: Maxicircle Structure and Heterogeneity of Minicircles

The mitochondrial DNA (mtDNA), which is present in almost all eukaryotic organisms, is a useful marker for phylogenetic studies due to its relative high conservation and its inheritance manner. In Leishmania and other trypanosomatids, the mtDNA (also referred to as kinetoplast DNA or kDNA) is composed of thousands of minicircles and a few maxicircles, catenated together into a complex network. Maxicircles are functionally similar to other eukaryotic mtDNAs, whereas minicircles are involved in RNA editing of some maxicircle-encoded transcripts. Next-generation sequencing (NGS) is increasingly used for assembling nuclear genomes and, currently, a large number of genomic sequences are available. However, most of the time, the mitochondrial genome is ignored in the genome assembly processes. The aim of this study was to develop a pipeline to assemble Leishmania minicircles and maxicircle DNA molecules, exploiting the raw data generated in the NGS projects. As a result, the maxicircle molecules and the plethora of minicircle classes for Leishmania major, Leishmania infantum and Leishmania braziliensis have been characterized. We have observed that whereas the heterogeneity of minicircle sequences existing in a single cell hampers their use for Leishmania typing and classification, maxicircles emerge as an extremely robust genetic marker for taxonomic studies within the clade of kinetoplastids.