The hard tick Ixodes ricinus is an important disease vector whose salivary secretions mediate blood-feeding success on vertebrate hosts, including humans. Here we describe the expression profiles and downstream analysis of de novo-discovered microRNAs (miRNAs) expressed in I. ricinus salivary glands and saliva. Eleven tick-derived libraries were sequenced to produce 67,375,557 Illumina reads. De novo prediction yielded 67 bona fide miRNAs out of which 35 are currently not present in miRBase. We report for the first time the presence of microRNAs in tick saliva, obtaining furthermore molecular indicators that those might be of exosomal origin. Ten out of these microRNAs are at least 100 times more represented in saliva. For the four most expressed microRNAs from this subset, we analyzed their combinatorial effects upon their host transcriptome using a novel in silico target network approach. We show that only the inclusion of combinatorial effects reveals the functions in important pathways related to inflammation and pain sensing. A control set of highly abundant microRNAs in both saliva and salivary glands indicates no significant pathways and a far lower number of shared target genes. Therefore, the analysis of miRNAs from pure tick saliva strongly supports the hypothesis that tick saliva miRNAs can modulate vertebrate host homeostasis and represents the first direct evidence of tick miRNA-mediated regulation of vertebrate host gene expression at the tick-host interface. As such, the herein described miRNAs may support future drug discovery and development projects that will also experimentally question their predicted molecular targets in the vertebrate host.

Bioinformatics Group Department of Computer Science Interdisciplinary Center for Bioinformatics University of Leipzig D 04107 Leipzig Germany

Biology Center of the Czech Academy of Sciences 37005 Budweis Czech Republic

Computational Genomics and Bioinformatics Group Genetics Department University of Granada 18071 Granada Spain

ecSeq Bioinformatics D 04275 Leipzig Germany

Laboratorio de Bioinformática Centro de Investigación Biomédica PTS 18100 Granada Spain

Zobrazit více v PubMed

Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, Gasteiger E, Huang H, Lopez R, Magrane M, et al. 2004. UniProt: the Universal Protein knowledgebase. Nucleic Acids Res 32: D115–119. PubMed PMC

Asgari S. 2014. Role of microRNAs in arbovirus/vector interactions. Viruses 6: 3514–3534. PubMed PMC

Bao W, Kojima KK, Kohany O. 2015. Repbase Update, a database of repetitive elements in eukaryotic genomes. Mob DNA 6: 11. PubMed PMC

Barrero RA, Keeble-Gagnere G, Zhang B, Moolhuijzen P, Ikeo K, Tateno Y, Gojobori T, Guerrero FD, Lew-Tabor A, Bellgard M. 2011. Evolutionary conserved microRNAs are ubiquitously expressed compared to tick-specific miRNAs in the cattle tick Rhipicephalus (Boophilus) microplus. BMC Genomics 12: 328. PubMed PMC

Barturen G, Rueda A, Hamberg M, Alganza A, Lebron R, Kotsyfakis M, Shi BJ, Koppers-Lalic D, Hackenberg M. 2014. sRNAbench: profiling of small RNAs and its sequence variants in single or multi-species high-throughput experiments. Methods Next Generation Seq 1: 21–31.

Blair CD, Olson KE. 2015. The role of RNA interference (RNAi) in arbovirus-vector interactions. Viruses 7: 820–843. PubMed PMC

Brugger K, Boehnke D, Petney T, Dobler G, Pfeffer M, Silaghi C, Schaub GA, Pinior B, Dautel H, Kahl O, et al. 2016. A density map of the tick-borne encephalitis and Lyme Borreliosis vector Ixodes ricinus (Acari: Ixodidae) for Germany. J Med Entomol 53: 1292–1302. PubMed

Chan PP, Lowe TM. 2016. GtRNAdb 2.0: an expanded database of transfer RNA genes identified in complete and draft genomes. Nucleic Acids Res 44: D184–D189. PubMed PMC

Chmelař J, Kotál J, Kopecký J, Pedra JH, Kotsyfakis M. 2016. All for one and one for all on the tick-host battlefield. Trends Parasitol 32: 368–377. PubMed PMC

Cramaro WJ, Revets D, Hunewald OE, Sinner R, Reye AL, Muller CP. 2015. Integration of Ixodes ricinus genome sequencing with transcriptome and proteome annotation of the naive midgut. BMC Genomics 16: 871. PubMed PMC

Dickinson RG, O'Hagan JE, Schotz M, Binnington KC, Hegarty MP. 1976. Prostaglandin in the saliva of the cattle tick Boophilus microplus. Aust J Exp Biol Med Sci 54: 475–486. PubMed

Fezza F, Dillwith JW, Bisogno T, Tucker JS, Di Marzo V, Sauer JR. 2003. Endocannabinoids and related fatty acid amides, and their regulation, in the salivary glands of the lone star tick. Biochim Biophys Acta 1633: 61–67. PubMed

Franceschini A, Szklarczyk D, Frankild S, Kuhn M, Simonovic M, Roth A, Lin J, Minguez P, Bork P, von Mering C, et al. 2013. STRING v9.1: protein-protein interaction networks, with increased coverage and integration. Nucleic Acids Res 41: D808–D815. PubMed PMC

Giraldo-Calderón GI, Emrich SJ, MacCallum RM, Maslen G, Dialynas E, Topalis P, Ho N, Gesing S, the VectorBase Consortium, Madey G, et al. 2015. VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic Acids Res 43: D707–D713. PubMed PMC

Griffiths RB, Gordon RM. 1952. An apparatus which enables the process of feeding by mosquitoes to be observed in the tissues of a live rodent; together with an account of the ejection of saliva and its significance in Malaria. Ann Trop Med Parasitol 46: 311–319. PubMed

Gulia-Nuss M, Nuss AB, Meyer JM, Sonenshine DE, Roe RM, Waterhouse RM, Sattelle DB, de la Fuente J, Ribeiro JM, Megy K, et al. 2016. Genomic insights into the Ixodes scapularis tick vector of Lyme disease. Nat Commun 7: 10507. PubMed PMC

Hackenberg M, Rodríguez-Ezpeleta N, Aransay AM. 2011. miRanalyzer: an update on the detection and analysis of microRNAs in high-throughput sequencing experiments. Nucleic Acids Res 39: W132–W138. PubMed PMC

He Y, Ju C, Zhang X. 2015. Roles of small RNAs in the immune defense mechanisms of crustaceans. Mol Immunol 68: 399–403. PubMed

Honig V, Svec P, Halas P, Vavruskova Z, Tykalova H, Kilian P, Vetiskova V, Dornakova V, Sterbova J, Simonova Z, et al. 2015. Ticks and tick-borne pathogens in South Bohemia (Czech Republic)—spatial variability in Ixodes ricinus abundance, Borrelia burgdorferi and tick-borne encephalitis virus prevalence. Ticks Tick Borne Dis 6: 559–567. PubMed

Horká H, Cerná-Kýcková K, Skallová A, Kopecký J. 2009. Tick saliva affects both proliferation and distribution of Borrelia burgdorferi spirochetes in mouse organs and increases transmission of spirochetes to ticks. Int J Med Microbiol 299: 373–380. PubMed

Hudson A, Bowman L, Orr CW. 1960. Effects of absence of saliva on blood feeding by mosquitoes. Science 131: 1730–1731. PubMed

Jeon YH, Youn DH. 2015. Spinal gap junction channels in neuropathic pain. Korean J Pain 28: 231–235. PubMed PMC

John B, Enright AJ, Aravin A, Tuschl T, Sander C, Marks DS. 2004. Human MicroRNA targets. PLoS Biol 2: e363. PubMed PMC

Kazimírová M, Štibrániová I. 2013. Tick salivary compounds: their role in modulation of host defences and pathogen transmission. Front Cell Infect Microbiol 3: 43. PubMed PMC

Kertesz M, Iovino N, Unnerstall U, Gaul U, Segal E. 2007. The role of site accessibility in microRNA target recognition. Nat Genet 39: 1278–1284. PubMed

Koppers-Lalic D, Hackenberg M, Bijnsdorp IV, van Eijndhoven MA, Sadek P, Sie D, Zini N, Middeldorp JM, Ylstra B, de Menezes RX, et al. 2014. Nontemplated nucleotide additions distinguish the small RNA composition in cells from exosomes. Cell Rep 8: 1649–1658. PubMed

Kotsyfakis M, Kopáček P, Franta Z, Pedra JH, Ribeiro JM. 2015a. Deep sequencing analysis of the Ixodes ricinus haemocytome. PLoS Negl Trop Dis 9: e0003754. PubMed PMC

Kotsyfakis M, Schwarz A, Erhart J, Ribeiro JM. 2015b. Tissue- and time-dependent transcription in Ixodes ricinus salivary glands and midguts when blood feeding on the vertebrate host. Sci Rep 5: 9103. PubMed PMC

Langmead B, Trapnell C, Pop M, Salzberg SL. 2009. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol 10: R25. PubMed PMC

Laslett D, Canback B. 2004. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences. Nucleic Acids Res 32: 11–16. PubMed PMC

Lorenz R, Bernhart SH, Höner Zu Siederdissen C, Tafer H, Flamm C, Stadler PF, Hofacker IL. 2011. ViennaRNA Package 2.0. Algorithms Mol Biol 6: 26. PubMed PMC

Lucas K, Raikhel AS. 2013. Insect microRNAs: biogenesis, expression profiling and biological functions. Insect Biochem Mol Biol 43: 24–38. PubMed PMC

Luhur A, Chawla G, Sokol NS. 2013. MicroRNAs as components of systemic signaling pathways in Drosophila melanogaster. Curr Top Dev Biol 105: 97–123. PubMed

Luo J, Liu GY, Chen Z, Ren QY, Yin H, Luo JX, Wang H. 2015. Identification and characterization of microRNAs by deep-sequencing in Hyalomma anatolicum anatolicum (Acari: Ixodidae) ticks. Gene 564: 125–133. PubMed

Marwaha L, Bansal Y, Singh R, Saroj P, Bhandari R, Kuhad A. 2016. TRP channels: potential drug target for neuropathic pain. Inflammopharmacology 24: 305–317. PubMed

Min H, Yoon S. 2010. Got target? Computational methods for microRNA target prediction and their extension. Exp Mol Med 42: 233–244. PubMed PMC

Mudenda L, Pierlé SA, Turse JE, Scoles GA, Purvine SO, Nicora CD, Clauss TR, Ueti MW, Brown WC, Brayton KA. 2014. Proteomics informed by transcriptomics identifies novel secreted proteins in Dermacentor andersoni saliva. Int J Parasitol 44: 1029–1037. PubMed

Mullokandov G, Baccarini A, Ruzo A, Jayaprakash AD, Tung N, Israelow B, Evans MJ, Sachidanandam R, Brown BD. 2012. High-throughput assessment of microRNA activity and function using microRNA sensor and decoy libraries. Nat Methods 9: 840–846. PubMed PMC

Nelder MP, Russell CB, Sheehan NJ, Sander B, Moore S, Li Y, Johnson S, Patel SN, Sider D. 2016. Human pathogens associated with the blacklegged tick Ixodes scapularis: a systematic review. Parasit Vectors 9: 265. PubMed PMC

Oliveira CJ, Sá-Nunes A, Francischetti IM, Carregaro V, Anatriello E, Silva JS, Santos IK, Ribeiro JM, Ferreira BR. 2011. Deconstructing tick saliva: non-protein molecules with potent immunomodulatory properties. J Biol Chem 286: 10960–10969. PubMed PMC

Ribeiro JM. 1987. Role of saliva in blood-feeding by arthropods. Annu Rev Entomol 32: 463–478. PubMed

Ribeiro JM, Makoul GT, Levine J, Robinson DR, Spielman A. 1985. Antihemostatic, antiinflammatory, and immunosuppressive properties of the saliva of a tick, Ixodes dammini. J Exp Med 161: 332–344. PubMed PMC

Ribeiro JM, Makoul GT, Robinson DR. 1988. Ixodes dammini: evidence for salivary prostacyclin secretion. J Parasitol 74: 1068–1069. PubMed

The RNAcentral Consortium. 2015. RNAcentral: an international database of ncRNA sequences. Nucleic Acids Res 43: D123–D129. PubMed PMC

Robinson MD, McCarthy DJ, Smyth GK. 2010. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26: 139–140. PubMed PMC

Rueda A, Barturen G, Lebrón R, Gómez-Martín C, Alganza A, Oliver JL, Hackenberg M. 2015. sRNAtoolbox: an integrated collection of small RNA research tools. Nucleic Acids Res 43: W467–W473. PubMed PMC

Sá-Nunes A, Bafica A, Lucas DA, Conrads TP, Veenstra TD, Andersen JF, Mather TN, Ribeiro JM, Francischetti IM. 2007. Prostaglandin E2 is a major inhibitor of dendritic cell maturation and function in Ixodes scapularis saliva. J Immunol 179: 1497–1505. PubMed

Schwarz A, Tenzer S, Hackenberg M, Erhart J, Gerhold-Ay A, Mazur J, Kuharev J, Ribeiro JM, Kotsyfakis M. 2014. A systems level analysis reveals transcriptomic and proteomic complexity in Ixodes ricinus midgut and salivary glands during early attachment and feeding. Mol Cell Proteomics 13: 2725–2735. PubMed PMC

Shao CC, Xu MJ, Chen YZ, Tao JP, Zhu XQ. 2015. Comparative profiling of microRNAs in male and female Rhipicephalus sanguineus. Appl Biochem Biotechnol 176: 1928–1936. PubMed

Sturm M, Hackenberg M, Langenberger D, Frishman D. 2010. TargetSpy: a supervised machine learning approach for microRNA target prediction. BMC Bioinformatics 11: 292. PubMed PMC

Tirloni L, Reck J, Terra RM, Martins JR, Mulenga A, Sherman NE, Fox JW, Yates JR III, Termignoni C, Pinto AF, et al. 2014. Proteomic analysis of cattle tick Rhipicephalus (Boophilus) microplus saliva: a comparison between partially and fully engorged females. PLoS One 9: e94831. PubMed PMC

Tirloni L, Islam MS, Kim TK, Diedrich JK, Yates JR III, Pinto AF, Mulenga A, You MJ, Da Silva Vaz I Jr. 2015. Saliva from nymph and adult females of Haemaphysalis longicornis: a proteomic study. Parasit Vectors 8: 338. PubMed PMC

Tsetsarkin KA, Liu G, Kenney H, Hermance M, Thangamani S, Pletnev AG. 2016. Concurrent micro-RNA mediated silencing of tick-borne flavivirus replication in tick vector and in the brain of vertebrate host. Sci Rep 6: 33088. PubMed PMC

Vay L, Gu C, McNaughton PA. 2012. The thermo-TRP ion channel family: properties and therapeutic implications. Br J Pharmacol 165: 787–801. PubMed PMC

Vayssier-Taussat M, Kazimirova M, Hubalek Z, Hornok S, Farkas R, Cosson JF, Bonnet S, Vourch G, Gasqui P, Mihalca AD, et al. 2015. Emerging horizons for tick-borne pathogens: from the ‘one pathogen-one disease’ vision to the pathobiome paradigm. Future Microbiol 10: 2033–2043. PubMed PMC

Wang F, Gong H, Zhang H, Zhou Y, Cao J, Zhou J. 2015. Lipopolysaccharide-induced differential expression of miRNAs in male and female Rhipicephalus haemaphysaloides ticks. PLoS One 10: e0139241. PubMed PMC

Wong P, Laxton V, Srivastava S, Chan YW, Tse G. 2017. The role of gap junctions in inflammatory and neoplastic disorders. Int J Mol Med 39: 498–506. PubMed PMC

Wikel S. 2013. Ticks and tick-borne pathogens at the cutaneous interface: host defenses, tick countermeasures, and a suitable environment for pathogen establishment. Front Microbiol 4: 337. PubMed PMC

Zhou J, Zhou Y, Cao J, Zhang H, Yu Y. 2013. Distinctive microRNA profiles in the salivary glands of Haemaphysalis longicornis related to tick blood-feeding. Exp Appl Acarol 59: 339–349. PubMed

Transcriptomic analysis of the tick midgut and salivary gland responses upon repeated blood-feeding on a vertebrate host

Tick extracellular vesicles enable arthropod feeding and promote distinct outcomes of bacterial infection

Tick Salivary Compounds for Targeted Immunomodulatory Therapy