Abstract
Ticks are blood-feeding arthropods that use the components of their salivary glands to counter the host's hemostatic, inflammatory, and immune responses. The tick midgut also plays a crucial role in hematophagy. It is responsible for managing blood meals (storage and digestion) and protecting against host immunity and pathogen infections. Previous transcriptomic studies revealed the complexity of tick sialomes (salivary gland transcriptomes) and mialomes (midgut transcriptomes) which encode for protease inhibitors, lipocalins (histamine-binding proteins), disintegrins, enzymes, and several other tick-specific proteins. Several studies have demonstrated that mammalian hosts acquire tick resistance against repeated tick bites. Consequently, there is an urgent need to uncover how tick sialomes and mialomes respond to resistant hosts, as they may serve to develop novel tick control strategies and applications. Here, we mimicked natural repeated tick bites in a laboratory setting and analyzed gene expression dynamics in the salivary glands and midguts of adult female ticks. Rabbits were subjected to a primary (feeding on a naive host) and a secondary infestation of the same host (we re-exposed the hosts but to other ticks). We used single salivary glands and midguts dissected from individual siblings adult pathogen-free female Ixodes ricinus to reduce genetic variability between individual ticks. The comprehensive analysis of 88 obtained RNA-seq data sets allows us to provide high-quality annotated sialomes and mialomes from individual ticks. Comparisons between fed/unfed, timepoints, and exposures yielded as many as 3000 putative differentially expressed genes (DEG). Interestingly, when classifying the exposure DEGs by means of a clustering approach we observed that the majority of these genes show increased expression at early feeding time-points in the mid-gut of re-exposed ticks. The existence of clearly defined groups of genes with highly similar responses to re-exposure suggests the existence of molecular swiches. In silico functional analysis shows that these early feeding reexposure response genes form a dense interaction network at protein level being related to virtually all aspects of gene expression regulation and glycosylation. The processed data is available through an easy-to-use database-associated webpage (https://arn.ugr.es/IxoriDB/) that can serve as a valuable resource for tick research.