A transcriptome-based phylogenetic study of hard ticks (Ixodidae)

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.

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.

Endogenous Viral Elements in Ixodid Tick Genomes

The documentation of endogenous viral elements (EVEs; virus-derived genetic material integrated into the genome of a nonviral host) has offered insights into how arthropods respond to viral infection via RNA interference pathways. Small non-coding RNAs derived from EVE loci serve to direct RNAi pathways in limiting replication and infection from cognate viruses, thus benefiting the host's fitness and, potentially, vectorial capacity. Here we use informatic approaches to analyze nine available genome sequences of hard ticks (Acari: Ixodidae; Rhipicephalus sanguineus, R. microplus, R. annulatus, Ixodes ricinus, I. persulcatus, I. scapularis, Hyalomma asiaticum, Haemaphysalis longicornis, and Dermacentor silvarum) to identify endogenous viral elements and to illustrate the shared ancestry of all elements identified. Our results highlight a broad diversity of viral taxa as having given rise to 1234 identified EVEs in ticks, with Mononegavirales (specifically Rhabdoviridae) well-represented in this subset of hard ticks. Further investigation revealed extensive adintovirus integrations in several Ixodes species, the prevalence of Bunyavirales EVEs (notably not observed in mosquitoes), and the presence of several elements similar to known emerging human and veterinary pathogens. These results will inform subsequent work on current and past associations with tick species with regard to the viruses from which their "viral fossils" are derived and may serve as a reference for quality control of various tick-omics data that may suffer from misidentification of EVEs as viral genetic material.

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.

Development and field evaluation of PCR assays based on minimum length Bm86 cDNA fragments required for Rhipicephalus and Hyalomma tick species delineation

Introduction

Hyalomma and Rhipicephalus ticks are important genera that can transmit diseases to both animals and humans, including Crimean-Congo hemorrhagic fever, tick-borne encephalitis, and several types of spotted fever. The accurate identification of tick species is essential for the effective control and prevention of tick-borne diseases. However, traditional identification methods based on morphology can be challenging and subjective, leading to errors. The development of DNA markers has provided more precise and efficient methods for tick species identification, but the currently available markers have limitations in their discriminatory power and sensitivity. To address this need for more sensitive and specific markers, this study aimed to identify two minimum sequence fragments required for tick Hyalomma and Rhipicephalus species identification using the Bm86 cDNA marker, which has previously been shown to be in perfect agreement with the current taxonomy of hard ticks based on its complete sequence.

Methods

Based on our in silico determination that a minimum sequence of 398 bp for Rhipicephalus spp. (from 1487 to 1884) and 559 bp for Hyalomma species (from 539 to 1097) was necessary for species delineation, two distinct PCR assays were developed to apply these sequences in practice.

Results and discussion

Discrimination between species within each genus was achieved through sequence homology and phylogenetic analysis following the sequencing of the two PCR products. Subsequently, their performance was evaluated by testing them on the field-collected ticks of the Hyalomma and Rhipicephalus genera obtained from various host animals in different geographic regions of Tunisia. The use of shorter partial sequences specific to the tick genera Rhipicephalus and Hyalomma, which target the tick's RNA banks, could represent a significant advance in the field of tick species identification, providing a sensitive and discriminatory tool for interspecific and intraspecific diversity analysis.

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.

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

Hoogstraal and Kim (1985) proposed from morphology, three groups of Haemaphysalis subgenera: (i) the "structurally advanced"; (ii) the "structurally intermediate"; and (iii) the "structurally primitive" subgenera. Nuclear gene phylogenies, however, did not indicate monophyly of these morphological groups but alas, only two mitochondrial (mt) genomes from the "structurally intermediate" subgenera had been sequenced. The phylogeny of Haemaphysalis has not yet been resolved. We aimed to resolve the phylogeny of the genus Haemaphysalis, with respect to the subgenus Alloceraea. We presented 15 newly sequenced and annotated mt genomes from 15 species of ticks, five species of which have not been sequenced before, and four new 18S rRNA and 28S rRNA nuclear gene sequences. Our datasets were constructed from 10 mt protein-coding genes, cox1, and the 18S and 28S nuclear rRNA genes. We found a 132-bp insertion between tRNA-Glu (E) gene and the nad1 gene in the mt genome of Haemaphysalis (Alloceraea) inermis that resembles insertions in H. (Alloceraea) kitaokai and Rhipicephalus (Boophilus) geigyi. Our mt phylogenies had the three species of Amblyomma (Aponomma) we sequenced embedded in the main clade of Amblyomma: Am. (Aponomma) fimbriatum, Am. (Aponomma) gervaisi and Am. (Aponomma) latum. This is further support for the hypothesis that the evolution of eyes appears to have occurred in the most-recent-common-ancestor of Amblyocephalus (i.e. Amblyomminae plus Rhipicephalinae) and that eyes were subsequently lost in the most-recent-common-ancestor of the subgenus Am. (Aponomma). Either Africaniella transversale or Robertsicus elaphensis, or perhaps Af. transversale plus Ro. elaphensis, appear to be the sister-group to the rest of the metastriate Ixodida. Our cox1 phylogenies did not indicate monophyly of the "structurally primitive", "structurally intermediate" nor the "structurally advanced" groups of Haemaphysalis subgenera. Indeed, the subgenus Alloceraea may be the only monophyletic subgenus of the genus Haemaphysalis sequenced thus far. All of our mt genome and cox1 phylogenies had the subgenus Alloceraea in a clade that was separate from the rest of the Haemaphysalis ticks. If Alloceraea is indeed the sister to the rest of the Haemaphysalis subgenera this would resonate with the argument of Hoogstraal and Kim (1985), that Alloceraea was a subgenus of "primitive" Haemaphysalis. Alectorobius capensis from Japan had a higher genetic-identity to A. sawaii, which was also from Japan, than to the A. capensis from South Africa. This indicates that A. capensis from Japan may be a cryptic species with respect to the A. capensis from South Africa.

Exploring Evolutionary Relationships within Neodermata Using Putative Orthologous Groups of Proteins, with Emphasis on Peptidases

The phylogenetic relationships within Neodermata were examined based on putative orthologous groups of proteins (OGPs) from 11 species of Monogenea, Trematoda, and Cestoda. The dataset included OGPs from BUSCO and OMA. Additionally, peptidases were identified and evaluated as phylogenetic markers. Phylogenies were inferred using the maximum likelihood method. A network analysis and a hierarchical grouping analysis of the principal components (HCPC) of orthologous groups of peptidases were performed. The phylogenetic analyses showed the monopisthocotylean monogeneans as the sister-group of cestodes, and the polyopisthocotylean monogeneans as the sister-group of trematodes. However, the sister-group relationship between Monopisthocotylea and Cestoda was not statistically well supported. The network analysis and HCPC also showed a cluster formed by polyopisthocotyleans and trematodes. The present study supports the non-monophyly of Monogenea. An analysis of mutation rates indicated that secreted peptidases and inhibitors, and those with multiple copies, are under positive selection pressure, which could explain the expansion of some families such as C01, C19, I02, and S01. Whilst not definitive, our study presents another point of view in the discussion of the evolution of Neodermata, and we hope that our data drive further discussion and debate on this intriguing topic.

A survey of hard ticks associated with cave dwelling mammals in Turkey

Members of the subgenus Eschatocephalus Frauenfeld, 1853 are highly specialized cave-dwelling ectoparasites of bats. We conducted a comprehensive field-based survey on bat ticks in Turkey and provided information about the phylogenetic placement of collected species. Sampling was carried out at 26 caves from 18 provinces around Turkey between 2019 and 2021. Eighty-one tick specimens collected from the cave environment or on various cave roosting bats resulted in five species: Ixodes vespertilionis, Ixodes simplex, Ixodes ariadnae, Ixodes kaiseri, and Haemaphysalis erinacei. While I. simplex was the most frequently collected species with a rate of 56.2% mainly from bats (Miniopterus schreibersii), I. vespertilionis was the most prevalent species (65.4%) and found mainly on cave walls. The first record of I. ariadnae was provided for Turkey. Phylogenetic trees were built using mt 16S rDNA and COI markers. Our results demonstrated the presence of two distinct lineages of I. vespertilionis in Turkey; one lineage grouped with European isolates, whereas three sequences clustered separately. The phylogenetic pattern of I. simplex was consistent with previous results; this clade was clustered distantly to other bat tick species. The significance of the surprising records of H. erinacei and I. kaiseri in caves is also discussed.

Integrated Jingmenvirus Polymerase Gene in Ixodes ricinus Genome

Members of the jingmenviruses group have been found in arthropods and mammals on all continents except Australia and Antarctica. Two viruses of this group were isolated from patients with fever after a tick bite. Using a nested RT-PCR assay targeting a jingmenvirus polymerase gene fragment, we screened ticks collected in seven regions of Russia and found that the abundant jingmenvirus-positive were of Ixodes ricinus species, with the prevalence ranging from 19.8% to 34.3%. In all cases, DNase/RNase treatment suggested that the detected molecule was DNA and subsequent next generation sequencing (NGS) proved that the viral polymerase gene was integrated in the I. ricinus genome. The copy number of the integrated polymerase gene was quantified by qPCR relative to the ITS2 gene and estimated as 1.32 copies per cell. At least three different genetic variants of the integrated polymerase gene were found in the territory of Russia. Phylogenetic analysis of the integrated jingmenvirus polymerase gene showed the highest similarity with the sequence of the correspondent gene obtained in Serbia from I. ricinus.

Characterization of the complete mitochondrial genome of Ixodes granulatus (Ixodidae) and its phylogenetic implications

Ticks are deemed to be second only to mosquitoes as the most common vector of human infectious diseases worldwide that give rise to human and animal diseases and economic losses to livestock production. Our understanding of the phylogenetic analysis between tick lineages has been restricted by the phylogenetic markers of individual genes. Genomic data research could help advance our understanding of phylogenetic analysis and molecular evolution. Mitochondrial genomic DNA facilitated the phylogenetic analysis of eukaryotes containing ticks. In this study, we sequenced and assembled the circular complete mitogenome information of Ixodes granulatus. The 14,540-bp mitogenome consists of 37 genes, including 13 protein-coding genes (PCGs), two genes for ribosomal RNA (rRNAs), and 22 genes for transfer RNA (tRNAs), and the origin of the L-strand replication region. The directions of the coding strand and component genes in the non-Australasian Ixodes mitochondrial genome were similar to those found in most other Australasian Ixodes, except for the loss of a lengthy control region. The phylogenetic tree based on maximum likelihood (ML) and Bayesian inference (BI) computational algorithms showed that I. granulatus exhibits a close relationship with I. hexagonus and I. ricinus. To our knowledge, this is the first study exploring the complete mitogenome for the species I. granulatus. Our results provide new insights for further research on the evolution, population genetics, systematics, and molecular ecology of ticks.

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.

A remarkable assemblage of ticks from mid-Cretaceous Burmese amber

Four fossil ticks (Arachnida: Parasitiformes: Ixodida) are described from mid-Cretaceous (ca. 100 Ma) Burmese amber of Myanmar. Ixodes antiquorum sp. nov. (Ixodidae) is the first Mesozoic record of Ixodes and the oldest representative of the most species-rich extant tick genus. Its affinities appear to lie with modern Australian forms, consistent with the hypothesis that Burmese amber hosted Gondwanan faunal elements. Even more remarkable is Khimaira fossus gen. et sp. nov. which combines a body resembling that of a soft tick (Argasidae) with a basis capitulum more like that of a hard tick (Ixodidae). We refer it to Khimairidae fam. nov. as a possible transitional form between the two main families of ticks alive today. Another member of the extinct Deinocrotonidae is described as Deinocroton copia sp. nov., while the first described adult female for Cornupalpatum burmanicum is associated with a dinosaur feather barb.

Necrotic reshaping of the glioma microenvironment drives disease progression

Glioblastoma is the most common primary brain tumor and has a dismal prognosis. The development of central necrosis represents a tipping point in the evolution of these tumors that foreshadows aggressive expansion, swiftly leading to mortality. The onset of necrosis, severe hypoxia and associated radial glioma expansion correlates with dramatic tumor microenvironment (TME) alterations that accelerate tumor growth. In the past, most have concluded that hypoxia and necrosis must arise due to "cancer outgrowing its blood supply" when rapid tumor growth outpaces metabolic supply, leading to diffusion-limited hypoxia. However, growing evidence suggests that microscopic intravascular thrombosis driven by the neoplastic overexpression of pro-coagulants attenuates glioma blood supply (perfusion-limited hypoxia), leading to TME restructuring that includes breakdown of the blood-brain barrier, immunosuppressive immune cell accumulation, microvascular hyperproliferation, glioma stem cell enrichment and tumor cell migration outward. Cumulatively, these adaptations result in rapid tumor expansion, resistance to therapeutic interventions and clinical progression. To inform future translational investigations, the complex interplay among environmental cues and myriad cell types that contribute to this aggressive phenotype requires better understanding. This review focuses on contributions from intratumoral thrombosis, the effects of hypoxia and necrosis, the adaptive and innate immune responses, and the current state of targeted therapeutic interventions.

The complete mitochondrial genome of Ixodes vespertilionis (Acari: Ixodidae)

Ixodes vespertilionis is a tick parasitizing on the bodies of bats. In our study, the complete mitogenome of I. vespertilionis was determined by using Illumina sequencing technology. The mitogenome was 14,559 bp in size and was predicted to encode 37 genes including 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and one control region. The gene order of the mitogenome is identical to Argasidae and non-Australasian Prostriata. The phylogenetic analysis by the Maximum-likelihood method reveals that I. vespertilionis is phylogenetically closest to Ixodes simplex. These data provide novel reference for further studies on the population genetics and phylogenetics of ticks.

Taxonomic Sampling and Rare Genomic Changes Overcome Long-Branch Attraction in the Phylogenetic Placement of Pseudoscorpions

Long-branch attraction is a systematic artifact that results in erroneous groupings of fast-evolving taxa. The combination of short, deep internodes in tandem with long-branch attraction artifacts has produced empirically intractable parts of the Tree of Life. One such group is the arthropod subphylum Chelicerata, whose backbone phylogeny has remained unstable despite improvements in phylogenetic methods and genome-scale data sets. Pseudoscorpion placement is particularly variable across data sets and analytical frameworks, with this group either clustering with other long-branch orders or with Arachnopulmonata (scorpions and tetrapulmonates). To surmount long-branch attraction, we investigated the effect of taxonomic sampling via sequential deletion of basally branching pseudoscorpion superfamilies, as well as varying gene occupancy thresholds in supermatrices. We show that concatenated supermatrices and coalescent-based summary species tree approaches support a sister group relationship of pseudoscorpions and scorpions, when more of the basally branching taxa are sampled. Matrix completeness had demonstrably less influence on tree topology. As an external arbiter of phylogenetic placement, we leveraged the recent discovery of an ancient genome duplication in the common ancestor of Arachnopulmonata as a litmus test for competing hypotheses of pseudoscorpion relationships. We generated a high-quality developmental transcriptome and the first genome for pseudoscorpions to assess the incidence of arachnopulmonate-specific duplications (e.g., homeobox genes and miRNAs). Our results support the inclusion of pseudoscorpions in Arachnopulmonata (new definition), as the sister group of scorpions. Panscorpiones (new name) is proposed for the clade uniting Scorpiones and Pseudoscorpiones.

Sequence diversity and evolution of a group of iflaviruses associated with ticks

We studied a group of tick-associated viruses with characteristics of members of the family Iflaviridae, a family of viruses frequently found in arthropods. Our aim was to gain insight into the evolutionary dynamics of this group of viruses, which may be linked to the biology of ticks. We explored assembled RNA-Seq data sets for different species of ticks. We identified members of five different iflavirus species, four of them novel, and discovered nine new genome sequences, including variants. Five variants represented a virus species associated with Ixodes ricinus. Unexpectedly, a sequence found in the Ixodes scapularis cell line ISE6 was nearly identical to the sequences of I. ricinus variants, suggesting a contamination of this cell line by I. ricinus material. Analysing patterns of substitutions between these variants, we detected a strong excess of synonymous mutations, suggesting evolution under strong positive selection. The phylogenies of the viruses and of their tick hosts were not congruent, suggesting recurrent host changes across tick genera during their evolution. Overall, our work constitutes a step in the understanding of the interactions between this family of viruses and ticks.

Morphology of Pholeoixodes species associated with carnivores in the western Palearctic: Pictorial key based on molecularly identified Ixodes (Ph.) canisuga, I. (Ph.) hexagonus and I. (Ph.) kaiseri males, nymphs and larvae

Three Palearctic members of the subgenus Pholeoixodes, i.e., Ixodes canisuga, Ixodes hexagonus and Ixodes kaiseri are frequently collected from dogs, cats, red foxes, badgers and other carnivorous/insectivorous hosts in Europe. While a pictorial identification key has been reported for female Pholeoixodes ticks, a similar work has not been done on their male, nymphal and larval specimens. This study was initiated in order to clarify and re-examine those morphological characters of these three tick species, which can be used relatively easily to identify/distinguish them. In the case of larvae the aims included finding alternatives to chaetotaxy, which is difficult to observe and its usefulness is also affected by uncertainties in literature data. For this, 609 Pholeoixodes ticks (males, nymphs and larvae) were collected from carnivores, hedgehogs and their environment in six European countries (representing Western, Central and Southeastern Europe), followed by detailed morphological examination and/or molecular analyses to confirm the identity of their species. Based on the morphology of 84 molecularly analyzed specimens and a new identification key compiled accordingly, altogether 116 I. canisuga, 277 I. hexagonus and 216 I. kaiseri males, nymphs and larvae were identified. Ixodes kaiseri was not found in Western Europe, where I. canisuga predominated. In Central Europe, all three Pholeoixodes species were collected, the largest number of specimens represented by I. hexagonus. On the other hand, in Southeastern Europe I. kaiseri had the highest abundance. In conclusion, the morphology of internal spur on the first coxae (as the traditionally used character to distinguish I. hexagonus from other Pholeoixodes species) is trustworthy to recognize males but is less informative in the case of nymphs and larvae. The latter can be identified more properly by observing the morphology of basis capituli. In particular, nymphs and larvae of I. canisuga have anteriorly flattened basis capituli, forming a plateau that surrounds the base of the hypostome. On the other hand, nymphs and larvae of I. hexagonus and I. kaiseri lack a similar plateau, but (unlike I. canisuga) have cornuae, which are either posterolaterally or caudally directed, respectively.

Tick species from Africa by migratory birds: a 3-year study in Italy

The role of resident or migratory birds in dispersal of tick species and tick-borne pathogens is still poorly known in Italy. We report here the results of a 3-year project based on sampling ticks from migratory birds, as well as from the vegetation at three stop-over sites for migrants, namely the islands of Ventotene (Latium), Asinara (Sardinia) and Ustica (Sicily). During the spring seasons from 2017-2019, in total 2681 ticks were collected, 2344 of which were sampled from migratory birds and 337 from the vegetation. Ticks were identified by morphology or by molecular tools when necessary. In total, 16 tick species were identified among which the following were exclusively found on birds: Hyalomma rufipes (43.3%), Hy. truncatum (0.1%), Ixodes frontalis (11.8%), Ix. inopinatus (0.2%), Ix. ricinus (3%), Haemaphysalis punctata (0.08%), Hae. erinacei (0.1%), Amblyomma variegatum (0.08%) and Argas vulgaris 0.1%), whereas five species were exclusively collected from the vegetation: Rhipicephalus bursa (10.5%), Rh. turanicus (5.9%), Rh. sanguineus sensu lato (2%), Rh. pusillus (2.4%), Hae. sulcata (0.08%). Hy. marginatum (10.3%) and Ix. ventalloi (9.3%) were found both on birds and on the vegetation on the island Ustica. It is worth noting that the search for ticks on the vegetation did not detect allochthonous tick species. Although we found several interesting local species and allochthonous ticks like Hy. rufipes, Am. variegatum and Ar. vulgaris on birds, further investigations are needed to better define the possible role of migratory birds in the introduction of ticks and tick-borne diseases in Italy, above all after the evidence of imported ticks positive to Crimean Congo hemorrhagic fever (CCHF) virus in several European countries.

Perpetuation of Borreliae

With one exception (epidemic relapsing fever), borreliae are obligately maintained in nature by ticks. Although some Borrelia spp. may be vertically transmitted to subsequent generations of ticks, most require amplification by a vertebrate host because inheritance is not stable. Enzootic cycles of borreliae have been found globally; those receiving the most attention from researchers are those whose vectors have some degree of anthropophily and, thus, cause zoonoses such as Lyme disease or relapsing fever. To some extent, our views on the synecology of the borreliae has been dominated by an applied focus, viz., analyses that seek to understand the elements of human risk for borreliosis. But, the elements of borrelial perpetuation do not necessarily bear upon risk, nor do our concepts of risk provide the best structure for analyzing perpetuation. We identify the major global themes for the perpetuation of borreliae, and summarize local variations on those themes, focusing on key literature to outline the factors that serve as the basis for the distribution and abundance of borreliae.

De novo assembly and analysis of the transcriptome of the Dermacentor marginatus genes differentially expressed after blood-feeding and long-term starvation

Background

The ixodid tick Dermacentor marginatus is a vector of many pathogens wide spread in Eurasia. Studies of gene sequence on many tick species have greatly increased the information on tick protective antigen which might have the potential to function as effective vaccine candidates or drug targets for eco-friendly acaricide development. In the current study, RNA-seq was applied to identify D. marginatus sequences and analyze differentially expressed unigenes.

Methods

To obtain a broader picture of gene sequences and changes in expression level, RNA-seq was performed to obtain the whole-body transcriptome data of D. marginatus adult female ticks after engorgement and long-term starvation. Subsequently, the real-time quantitative PCR (RT-qPCR) was applied to validate the RNA-seq data.

Results

RNA-seq produced 30,251 unigenes, of which 32% were annotated. Gene expression was compared among groups that differed by status as newly molted, starved and engorged female adult ticks. Nearly one third of the unigenes in each group were differentially expressed compared to the other two groups, and the most numerous were genes encoding proteins involved in catalytic and binding activities and apoptosis. Selected up-regulated differentially expressed genes in each group were associated to protein, lipids, carbohydrate and chitin metabolism. Blood-feeding and long-term starvation also caused genes differentially expressed in the defense response and antioxidant response. RT-qPCR results indicated 6 differentially expressed transcripts showed similar trends in expression changes with RNA-seq results confirming that the gene expression profiles in transcriptome data is in consistent with RT-qPCR validation.

Conclusions

Obtaining the sequence information of D. marginatus and characterizing the expression pattern of the genes involved in blood-feeding and during starvation would be helpful in understanding molecular physiology of D. marginatus and provides data for anti-tick vaccine and drug development for controlling the tick.

Phylogenies from mitochondrial genomes of 120 species of ticks: Insights into the evolution of the families of ticks and of the genus Amblyomma

The evolution and phylogenetic relationships of the ticks at both the family and genus levels are contested. The genus Amblyomma and its subgenera are in a state of flux; moreover, the relationships among the three tick families are controversial due to conflicting phylogenetic support for different arrangements of the three families of living ticks. With 18 newly sequenced mitochondrial (mt) genomes of ticks included, we executed the largest mt genome phylogenetic study of ticks so far. Phylogenetic trees were inferred from one sea spider mt genome, one horseshoe crab, five mite mt genomes and 146 tick mt genomes from 120 species: 153 mt genomes in total. Sixteen phylogenetic trees were inferred from 10 datasets using both maximum likelihood and Bayesian inference methods. We describe the first novel mt gene-arrangement for the metastriate Ixodidae in Amblyomma (Africaniella) transversale. Also, three unusual partial 16S rRNA gene inserts were found in the mt genome of Haemaphysalis (Alloceraea) kitaokai: we consider the possible role of past genome translocation events in the formation of these inserts. Our phylogenies revealed evidence that: (i) the genus Amblyomma is polyphyletic with respect to Amblyomma (Africaniella) transversale; (ii) the subgenus Aponomma is apparently embedded in the genus Amblyomma; (iii) Haemaphysalis (Segalia) parva and Haemaphysalis (Alloceraea) kitaokai form a clade to the exclusion of other Haemaphysalis species; and (iv) the phylogenetic position of the family Nuttalliellidae is unstable among phylogenies from different datasets.

Quantitative Visions of Reality at the Tick-Host Interface: Biochemistry, Genomics, Proteomics, and Transcriptomics as Measures of Complete Inventories of the Tick Sialoverse

Species have definitive genomes. Even so, the transcriptional and translational products of the genome are dynamic and subject to change over time. This is especially true for the proteins secreted by ticks at the tick-host feeding interface that represent a complex system known as the sialoverse. The sialoverse represent all of the proteins derived from tick salivary glands for all tick species that may be involved in tick-host interaction and the modulation of the host's defense mechanisms. The current study contemplates the advances made over time to understand and describe the complexity present in the sialoverse. Technological advances at given periods in time allowed detection of functions, genes, and proteins enabling a deeper insight into the complexity of the sialoverse and a concomitant expansion in complexity with as yet, no end in sight. The importance of systematic classification of the sialoverse is highlighted with the realization that our coverage of transcriptome and proteome space remains incomplete, but that complete descriptions may be possible in the future. Even so, analysis and integration of the sialoverse into a comprehensive understanding of tick-host interactions may require further technological advances given the high level of expected complexity that remains to be uncovered.