Evolution, expansion and expression of the Kunitz/BPTI gene family associated with long-term blood feeding in Ixodes Scapularis

The effect of feeding on different hosts on the egg proteins in Haemaphysalis qinghaiensis tick

The majority of ixodid ticks display host-specificity to varying extents. Feeding on different hosts affects their development and reproduction. Consequences can be analyzed at the level of the egg, as it is the initial stage of tick development. Tick egg proteins are abundant and diverse, providing nutrients for embryonic development. However, studies on tick egg profiles are scarce. In this study, we aimed to analyze whether feeding Haemaphysalis qinghaiensis ticks on the yaks (Bos grunniens) and domestic sheep (Ovis aries) has an impact on the variety and variability of the egg proteome. Detached engorged females were used to lay eggs, which were then collected, dewaxed, and subjected to protein extraction. The extracted egg proteins were enzymatically digested using Filter-Aided Sample Preparation (FASP), and the unique peptides were separated and detected by Liquid Chromatography-tandem Mass Spectrometry (LC-MS/MS). The MS data were searched against the previously constructed whole tick transcriptome library of H. qinghaiensis, and the UniProt database for the identification of tick-derived egg proteins. The analysis revealed 49 and 53 high-confidence proteins identified in eggs collected from B. grunniens (EggBg) and O. aries (EggOa), respectively. Of these, 46 high-confidence proteins were common to both egg types, while three were unique to EggBg and seven to EggOa. All the identified proteins mainly belonged to enzymes, enzyme inhibitors, transporters, and proteins with unknown functions. The differential abundance analysis showed that nine proteins were significantly more present in EggBg, while six were significantly more present in EggOa. Overall, enzymes were the most diverse group, while vitellogenin (Vg) was the most abundant. Blood meal uptake on different hosts has a certain effect on the egg proteome composition and the abundance of some proteins, but it may also lead to compensation of protein roles.

Bulk and single-nucleus RNA sequencing highlight immune pathways induced in individuals during an Ixodes scapularis tick bite

Ticks are hematophagous arthropods that use a complex mixture of salivary proteins to evade host defenses while taking a blood meal. Little is known about the immunological and physiological consequences of tick feeding on humans. Here, we performed the first bulk and single-nucleus RNA sequencing (snRNA-seq) of skin and blood of four persons presenting with naturally acquired, attached Ixodes scapularis ticks. Pathways and individual genes associated with innate and adaptive immunity were identified based on bulk RNA sequencing, including interleukin-17 signaling and platelet activation pathways at the site of tick attachment or in peripheral blood. snRNA-seq further revealed that the Hippo signaling, cell adhesion, and axon guidance pathways were involved in the response to an I. scapularis bite in humans. Features of the host response in these individuals also overlapped with that of laboratory guinea pigs exposed to I. scapularis and which acquired resistance to ticks. These findings offer novel insights for the development of new biomarkers for I. scapularis exposure and anti-tick vaccines for human use.

Host Immune Responses to Salivary Components - A Critical Facet of Tick-Host Interactions

Tick sialome is comprised of a rich cocktail of bioactive molecules that function as a tool to disarm host immunity, assist blood-feeding, and play a vibrant role in pathogen transmission. The adaptation of the tick's blood-feeding behavior has lead to the evolution of bioactive molecules in its saliva to assist them to overwhelm hosts' defense mechanisms. During a blood meal, a tick secretes different salivary molecules including vasodilators, platelet aggregation inhibitors, anticoagulants, anti-inflammatory proteins, and inhibitors of complement activation; the salivary repertoire changes to meet various needs such as tick attachment, feeding, and modulation or impairment of the local dynamic and vigorous host responses. For instance, the tick's salivary immunomodulatory and cement proteins facilitate the tick's attachment to the host to enhance prolonged blood-feeding and to modulate the host's innate and adaptive immune responses. Recent advances implemented in the field of "omics" have substantially assisted our understanding of host immune modulation and immune inhibition against the molecular dynamics of tick salivary molecules in a crosstalk between the tick-host interface. A deep understanding of the tick salivary molecules, their substantial roles in multifactorial immunological cascades, variations in secretion, and host immune responses against these molecules is necessary to control these parasites. In this article, we reviewed updated knowledge about the molecular mechanisms underlying host responses to diverse elements in tick saliva throughout tick invasion, as well as host defense strategies. In conclusion, understanding the mechanisms involved in the complex interactions between the tick salivary components and host responses is essential to decipher the host defense mechanisms against the tick evasion strategies at tick-host interface which is promising in the development of effective anti-tick vaccines and drug therapeutics.

Tick Saliva and Salivary Glands: What Do We Know So Far on Their Role in Arthropod Blood Feeding and Pathogen Transmission

Ticks are blood-sucking arthropods that have developed myriad of strategies to get a blood meal from the vertebrate host. They first attach to the host skin, select a bite site for a blood meal, create a feeding niche at the bite site, secrete plethora of molecules in its saliva and then starts feeding. On the other side, host defenses will try to counter-attack and stop tick feeding at the bite site. In this constant battle between ticks and the host, arthropods successfully pacify the host and completes a blood meal and then replete after full engorgement. In this review, we discuss some of the known and emerging roles for arthropod components such as cement, salivary proteins, lipocalins, HSP70s, OATPs, and extracellular vesicles/exosomes in facilitating successful blood feeding from ticks. In addition, we discuss how tick-borne pathogens modulate(s) these components to infect the vertebrate host. Understanding the biology of arthropod blood feeding and molecular interactions at the tick-host interface during pathogen transmission is very important. This information would eventually lead us in the identification of candidates for the development of transmission-blocking vaccines to prevent diseases caused by medically important vector-borne pathogens.

Sea Anemone Kunitz-Type Peptides Demonstrate Neuroprotective Activity in the 6-Hydroxydopamine Induced Neurotoxicity Model

Kunitz-type peptides from venomous animals have been known to inhibit different proteinases and also to modulate ion channels and receptors, demonstrating analgesic, anti-inflammatory, anti-histamine and many other biological activities. At present, there is evidence of their neuroprotective effects. We have studied eight Kunitz-type peptides of the sea anemone Heteractis crispa to find molecules with cytoprotective activity in the 6-OHDA-induced neurotoxicity model on neuroblastoma Neuro-2a cells. It has been shown that only five peptides significantly increase the viability of neuronal cells treated with 6-OHDA. The TRPV1 channel blocker, HCRG21, has revealed the neuroprotective effect that could be indirect evidence of TRPV1 involvement in the disorders associated with neurodegeneration. The pre-incubation of Neuro-2a cells with HCRG21 followed by 6-OHDA treatment has resulted in a prominent reduction in ROS production compared the untreated cells. It is possible that the observed effect is due to the ability of the peptide act as an efficient free-radical scavenger. One more leader peptide, InhVJ, has shown a neuroprotective activity and has been studied at concentrations of 0.01–10.0 µM. The target of InhVJ is still unknown, but it was the best of all eight homologous peptides in an absolute cell viability increment on 38% of the control in the 6-OHDA-induced neurotoxicity model. The targets of the other three active peptides remain unknown.

Harm or protection? The adaptive function of tick toxins

The existence of tick toxins is an old enigma that has intrigued scientists for a long time. The adaptive value of using deadly toxins for predatory animals is obvious: they try to kill the prey in the most effective way or protect themselves from their natural enemies. Ticks, however, are blood-sucking parasites, and it seems paradoxical that they have toxins similar to spiders, scorpions and snakes. Based on published data, here we examine the potential adaptive function of different types of toxins produced by soft and hard ticks. We hypothesize that there are diverse evolutionary roles behind (a) to attack and reduce the tick-transmitted pathogens inside the vertebrate host systemically to protect the tick, (b) to paralyse the host to stop grooming, (c) to speed up host heartbeat to improve blood supply and (d) to inhibit the process of necroptosis to prevent the rejection of hard ticks. We will provide published evidence that supports the above-mentioned hypotheses, and we will give an outlook how these new scientific results might be applied in modern pharmacology and medicine.

Kunitz-Type Peptides from the Sea Anemone Heteractis crispa Demonstrate Potassium Channel Blocking and Anti-Inflammatory Activities

The Kunitz/BPTI peptide family includes unique representatives demonstrating various biological activities. Electrophysiological screening of peptides HCRG1 and HCRG2 from the sea anemone Heteractis crispa on six Kv1.x channel isoforms and insect Shaker IR channel expressed in Xenopus laevis oocytes revealed their potassium channels blocking activity. HCRG1 and HCRG2 appear to be the first Kunitz-type peptides from sea anemones blocking Kv1.3 with IC₅₀ of 40.7 and 29.7 nM, respectively. In addition, peptides mainly vary in binding affinity to the Kv1.2 channels. It was established that the single substitution, Ser5Leu, in the TRPV1 channel antagonist, HCRG21, induces weak blocking activity of Kv1.1, Kv1.2, and Kv1.3. Apparently, for the affinity and selectivity of Kunitz-fold toxins to Kv1.x isoforms, the number and distribution along their molecules of charged, hydrophobic, and polar uncharged residues, as well as the nature of the channel residue at position 379 (Tyr, Val or His) are important. Testing the compounds in a model of acute local inflammation induced by the introduction of carrageenan administration into mice paws revealed that HCRG1 at doses of 0.1–1 mg/kg reduced the volume of developing edema during 24 h, similar to the effect of the nonsteroidal anti-inflammatory drug, indomethacin, at a dose of 5 mg/kg. ELISA analysis of the animals blood showed that the peptide reduced the synthesis of TNF-α, a pro-inflammatory mediator playing a leading role in the development of edema in this model.

TickSialoFam (TSFam): A Database That Helps to Classify Tick Salivary Proteins, a Review on Tick Salivary Protein Function and Evolution, With Considerations on the Tick Sialome Switching Phenomenon

Tick saliva contains a complex mixture of peptides and non-peptides that counteract their hosts' hemostasis, immunity, and tissue-repair reactions. Recent transcriptomic studies have revealed over one thousand different transcripts coding for secreted polypeptides in a single tick species. Not only do these gene products belong to many expanded families, such as the lipocalins, metalloproteases, Antigen-5, cystatins, and apyrases, but also families that are found exclusively in ticks, such as the evasins, Isac, DAP36, and many others. Phylogenetic analysis of the deduced protein sequences indicate that the salivary genes exhibit an increased rate of evolution due to a lower evolutionary constraint and/or positive selection, allowing for a large diversity of tick salivary proteins. Thus, for each new tick species that has its salivary transcriptome sequenced and assembled, a formidable task of annotation of these transcripts awaits. Currently, as of November 2019, there are over 287 thousand coding sequences deposited at the National Center for Biotechnology Information (NCBI) that are derived from tick salivary gland mRNA. Here, from these 287 thousand sequences we identified 45,264 potential secretory proteins which possess a signal peptide and no transmembrane domains on the mature peptide. By using the psiblast tools, position-specific matrices were constructed and assembled into the TickSialoFam (TSF) database. The TSF is a rpsblastable database that can help with the annotation of tick sialotranscriptomes. The TSA database identified 136 tick salivary secreted protein families, as well as 80 families of endosomal-related products, mostly having a protein modification function. As the number of sequences increases, and new annotation details become available, new releases of the TSF database may become available.

Salmon immunological defence and interplay with the modulatory capabilities of its ectoparasite Lepeophtheirus salmonis

The salmon louse Lepeophtheirus salmonis (Lsal) is an ectoparasitic copepod that exerts immunomodulatory and physiological effects on its host Atlantic salmon. Over 30 years of research on louse biology, control, host responses and the host-parasite relationship has provided a plethora of information on the intricacies of host resistance and parasite adaptation. Atlantic salmon exhibit temporal and spatial impairment of the immune system and wound healing ability during infection. This immunosuppression may render Atlantic salmon less tolerant to stress and other confounders associated with current management strategies. Contrasting susceptibility of salmonid hosts exists, and early pro-inflammatory Th1 type responses are associated with resistance. Rapid cellular responses to larvae appear to tip the balance of the host-parasite relationship in favour of the host, preventing severe immune-physiological impacts of the more invasive adults. Immunological, transcriptomic, genomic and proteomic evidence suggests pathological impacts occur in susceptible hosts through modulation of host immunity and physiology via pharmacologically active molecules. Co-evolutionary and farming selection pressures may have incurred preference of Atlantic salmon as a host for Lsal reflected in their interactome. Here, we review host-parasite interactions at the primary attachment/feeding site, and the complex life stage-dependent molecular mechanisms employed to subvert host physiology and immune responses.

Tick-Borne Encephalitis Virus Infection Alters the Sialome of Ixodes ricinus Ticks During the Earliest Stages of Feeding

Ticks are hematophagous arthropods that transmit a number of pathogens while feeding. Among these is tick-borne encephalitis virus (TBEV), a flavivirus transmitted by Ixodes ricinus ticks in the temperate zone of Europe. The infection results in febrile illness progressing to encephalitis and meningitis with a possibility of fatality or long-term neurological sequelae. The composition of tick saliva plays an essential role in the initial virus transmission during tick feeding. Ticks secrete a diverse range of salivary proteins to modulate the host response, such as lipocalins to control the itch and inflammatory response, and both proteases and protease inhibitors to prevent blood coagulation. Here, the effect of viral infection of adult females of Ixodes ricinus was studied with the goal of determining how the virus alters the tick sialome to modulate host tissue response at the site of infection. Uninfected ticks or those infected with TBEV were fed on mice and removed and dissected one- and 3-h post-attachment. RNA from the salivary glands of these ticks, as well as from unfed ticks, was extracted and subjected to next-generation sequencing to determine the expression of key secreted proteins at each timepoint. Genes showing statistically significant up- or down-regulation between infected and control ticks were selected and compared to published literature to ascertain their function. From this, the effect of tick viral infection on the modulation of the tick-host interface was determined. Infected ticks were found to differentially express a number of uncategorized genes, proteases, Kunitz-type serine protease inhibitors, cytotoxins, and lipocalins at different timepoints. These virus-induced changes to the tick sialome may play a significant role in facilitating virus transmission during the early stages of tick feeding.

Gain and loss events in the evolution of the apolipoprotein family in vertebrata

Background

Various apolipoproteins widely distributed among vertebrata play key roles in lipid metabolism and have a direct correlation with human diseases as diagnostic markers. However, the evolutionary progress of apolipoproteins in species remains unclear. Nine human apolipoproteins and well-annotated genome data of 30 species were used to identify 210 apolipoprotein family members distributed among species from fish to humans. Our study focused on the evolution of nine exchangeable apolipoproteins (ApoA-I/II/IV/V, ApoC-I~IV and ApoE) from Chondrichthyes, Holostei, Teleostei, Amphibia, Sauria (including Aves), Prototheria, Marsupialia and Eutheria.

Results

In this study, we reported the overall distribution and the frequent gain and loss evolutionary events of apolipoprotein family members in vertebrata. Phylogenetic trees of orthologous apolipoproteins indicated evident divergence between species evolution and apolipoprotein phylogeny. Successive gain and loss events were found by evaluating the presence and absence of apolipoproteins in the context of species evolution. For example, only ApoA-I and ApoA-IV occurred in cartilaginous fish as ancient apolipoproteins. ApoA-II, ApoE, and ApoC-I/ApoC-II were found in Holostei, Coelacanthiformes, and Teleostei, respectively, but the latter three apolipoproteins were absent from Aves. ApoC-I was also absent from Cetartiodactyla. The apolipoprotein ApoC-III emerged in terrestrial animals, and ApoC-IV first arose in Eutheria. The results indicate that the order of the emergence of apolipoproteins is most likely ApoA-I/ApoA-IV, ApoE, ApoA-II, ApoC-I/ApoC-II, ApoA-V, ApoC-III, and ApoC-IV.

Conclusions

This study reveals not only the phylogeny of apolipoprotein family members in species from Chondrichthyes to Eutheria but also the occurrence and origin of new apolipoproteins. The broad perspective of gain and loss events and the evolutionary scenario of apolipoproteins across vertebrata provide a significant reference for the research of apolipoprotein function and related diseases.

The Immunomodulatory Effect of IrSPI, a Tick Salivary Gland Serine Protease Inhibitor Involved in Ixodes ricinus Tick Feeding

Ticks are the most important vectors of pathogens affecting both domestic and wild animals worldwide. Hard tick feeding is a slow process—taking up to several days—and necessitates extended control over the host response. The success of the feeding process depends upon injection of tick saliva, which not only controls host hemostasis and wound healing, but also subverts the host immune response to avoid tick rejection that creates a favorable niche for the survival and propagation of diverse tick-borne pathogens. Here, we report on the molecular and biochemical features and functions of an Ixodes ricinus serine protease inhibitor (IrSPI). We characterize IrSPI as a Kunitz elastase inhibitor that is overexpressed in several tick organs—especially salivary glands—during blood-feeding. We also demonstrated that when IrSPI is injected into the host through saliva, it had no impact on tissue factor pathway-induced coagulation, fibrinolysis, endothelial cell angiogenesis or apoptosis, but the protein exhibits immunomodulatory activity. In particular, IrSPI represses proliferation of CD4⁺ T lymphocytes and proinflammatory cytokine secretion from both splenocytes and macrophages. Our study contributes valuable knowledge to tick-host interactions and provides insights that could be further exploited to design anti-tick vaccines targeting this immunomodulator implicated in I. ricinus feeding.

Construction and Composition of the Squid Pen from Doryteuthis pealeii

The pen, or gladius, of the squid is an internalized shell. It serves as a site of attachment for important muscle groups and as a protective barrier for the visceral organs. The pen's durability and flexibility are derived from its unique composition of chitin and protein. We report the characterization of the structure, development, and composition of pens from Doryteuthis pealeii. The nanofibrils of the polysaccharide β-chitin are arranged in an aligned configuration in only specific regions of the pen. Chitin is secreted early in development, enabling us to characterize the changes in pen morphology prior to hatching. The chitin and proteins are assembled in the shell sac surrounded by fluid that has a significantly different ionic composition from squid plasma. Two groups of proteins are associated with the pen: those on its surface and those embedded within the pen. Only 20 proteins are identified as embedded within the pen. Embedded proteins are classified into six groups, including chitin associated, protease, protease inhibitors, intracellular, extracellular matrix, and those that are unknown. The pen proteins share many conserved domains with proteins from other chitinous structures. We conclude that the pen is one of the least complex, load-bearing, chitin-rich structures currently known and is amenable to further studies to elucidate natural construction mechanisms using chitin and protein.

Structural characterization of PPTI, a kunitz-type protein from the venom of Pseudocerastes persicus

The main purpose of this report is to investigate the structural property and new potential function of PPTI (Pseudocerastes Persicus Trypsin Inhibitor), a kunitz-type protein with inhibitory effect against trypsin proteolytic activity. Besides kunitz-type serine protease inhibitors, PPTI shows clear-cut similarities with dendrotoxins (DTXs), the other kunitz-type protein subfamily. The most important reason is the presence of functionally important residues of DTXs at correspondingly the same positions in PPTI. As such, we proposed the new ability of PPTI for inhibiting voltage-gated potassium channels and consequently its dual functionality. At first, we determined the solution structure of PPTI via Nuclear Magnetic Resonance (NMR) spectroscopy. Then by homology modeling, we constructed the model structure of trypsin-PPTI complex to confirm the same interaction pattern as trypsin-BPTI at complex interface. Finally, by Brownian dynamics (BD) simulations of PPTI NMR derived ensemble structure as ligand against homology model of human Kv1.1 potassium channel as receptor, we evaluated the potential DTX-like activity of PPTI. The results of our study support the proposed dual functionality of PPTI.

Tick Paralysis: Solving an Enigma

In comparison to other arachnids, ticks are major vectors of disease, but less than 8% of the known species are capable of inducing paralysis, as compared to the ~99⁻100% arachnids that belong to venomous classes. When considering the potential monophyly of venomous Arachnida, this review reflects on the implications regarding the classification of ticks as venomous animals and the possible origin of toxins. The origin of tick toxins is compared with scorpion and spider toxins and venoms based on their significance, functionality, and structure in the search to find homologous venomous characters. Phenotypic evaluation of paralysis, as caused by different ticks, demonstrated the need for expansion on existing molecular data of pure isolated tick toxins because of differences and discrepancies in available data. The use of in-vivo, in-vitro, and in-silico assays for the purification and characterization of paralysis toxins were critically considered, in view of what may be considered to be a paralysis toxin. Purified toxins should exhibit physiologically relevant activity to distinguish them from other tick-derived proteins. A reductionist approach to identify defined tick proteins will remain as paramount in the search for defined anti-paralysis vaccines.

Novel Kunitz-like Peptides Discovered in the Zoanthid Palythoa caribaeorum through Transcriptome Sequencing

Palythoa caribaeorum (class Anthozoa) is a zoanthid that together jellyfishes, hydra, and sea anemones, which are venomous and predatory, belongs to the Phyllum Cnidaria. The distinguished feature in these marine animals is the cnidocytes in the body tissues, responsible for toxin production and injection that are used majorly for prey capture and defense. With exception for other anthozoans, the toxin cocktails of zoanthids have been scarcely studied and are poorly known. Here, on the basis of the analysis of P. caribaeorum transcriptome, numerous predicted venom-featured polypeptides were identified including allergens, neurotoxins, membrane-active, and Kunitz-like peptides (PcKuz). The three predicted PcKuz isotoxins (1-3) were selected for functional studies. Through computational processing comprising structural phylogenetic analysis, molecular docking, and dynamics simulation, PcKuz3 was shown to be a potential voltage gated potassium-channel inhibitor. PcKuz3 fitted well as new functional Kunitz-type toxins with strong antilocomotor activity as in vivo assessed in zebrafish larvae, with weak inhibitory effect toward proteases, as evaluated in vitro. Notably, PcKuz3 can suppress, at low concentration, the 6-OHDA-induced neurotoxicity on the locomotive behavior of zebrafish, which indicated PcKuz3 may have a neuroprotective effect. Taken together, PcKuz3 figures as a novel neurotoxin structure, which differs from known homologous peptides expressed in sea anemone. Moreover, the novel PcKuz3 provides an insightful hint for biodrug development for prospective neurodegenerative disease treatment.

Genome-wide comparison of the protein-coding repertoire reveals fast evolution of immune-related genes in cephalochordates and Osteichthyes superclass

Amphioxus is used to investigate the origin and evolution of vertebrates. To better understand the characteristics of genome evolution from cephalochordates to Osteichthyes, we conducted a genome-wide pairwise comparison of protein-coding genes within amphioxus (a comparable group) and parallel analyses within Osteichthyes (two comparable groups). A batch of fast-evolving genes in each comparable group was identified. Of these genes, the most fast-evolving genes (top 20) were scrutinized, most of which were involved in immune system. An analysis of the fast-evolving genes showed that they were enriched into gene ontology (GO) terms and pathways primarily involved in immune-related functions. Similarly, this phenomenon was detected within Osteichthyes, and more well-known and abundant GO terms and pathways involving innate immunity were found in Osteichthyes than in cephalochordates. Next, we measured the expression responses of four genes belonging to metabolism or energy production-related pathways to lipopolysaccharide challenge in the muscle, intestine or skin of B. belcheri; three of these genes (HMGCL, CYBS and MDH2) showed innate immune responses. Additionally, some genes involved in adaptive immunity showed fast evolution in Osteichthyes, such as those involving "intestinal immune network for IgA production" or "T-cell receptor signaling pathway". In this study, the fast evolution of immune-related genes in amphioxus and Osteichthyes was determined, providing insights into the evolution of immune-related genes in chordates.

Identification and characterization of serine protease inhibitors in a parasitic wasp, Pteromalus puparum

Serine protease inhibitors (SPIs) regulate protease-mediated activities by inactivating their cognate proteinases, and are involved in multiple physiological processes. SPIs have been extensively studied in vertebrates and invertebrates; however, little SPI information is available in parasitoids. Herein, we identified 57 SPI genes in total through the genome of a parasitoid wasp, Pteromalus puparum. Gene structure analyses revealed that these SPIs contain 7 SPI domains. Depending on their mode of action, these SPIs can be categorized into serpins, canonical inhibitors and alpha-2-macroglobulins (A2Ms). For serpins and canonical inhibitors, we predicted their putative inhibitory activities to trypsin/chymotrypsin/elastase-like enzymes based on the amino acids in cleaved reactive sites. Sequence alignment and phylogenetic tree indicated that some serpins similar to known functional inhibitory serpins may participate in immune responses. Transcriptome analysis also showed some canonical SPI genes displayed distinct expression patterns in the venom gland and this was confirmed by quantitative real-time PCR (qPCR) analysis, suggesting their specific physiological functions as venom proteins in suppressing host immune responses. The study provides valuable information to clarify the functions of SPIs in digestion, development, reproduction and innate immunity.

Insights into the molecular evolution of peptidase inhibitors in arthropods

Peptidase inhibitors are key proteins involved in the control of peptidases. In arthropods, peptidase inhibitors modulate the activity of peptidases involved in endogenous physiological processes and peptidases of the organisms with which they interact. Exploring available arthropod genomic sequences is a powerful way to obtain the repertoire of peptidase inhibitors in every arthropod species and to understand the evolutionary mechanisms involved in the diversification of this kind of proteins. A genomic comparative analysis of peptidase inhibitors in species belonging to different arthropod taxonomic groups was performed. The results point out: i) species or clade-specific presence is shown for several families of peptidase inhibitors; ii) multidomain peptidase inhibitors are commonly found in many peptidase inhibitor families; iii) several families have a wide range of members in different arthropod species; iv) several peptidase inhibitor families show species-specific (or clade-specific) gene family expansions; v) functional divergence may be assumed for particular clades; vi) passive expansions may be used by natural selection to fix adaptations. In conclusion, conservation and divergence of duplicated genes and the potential recruitment as peptidase inhibitors of proteins from other families are the main mechanisms used by arthropods to fix diversity. This diversity would be associated to the control of target peptidases and, as consequence, to adapt to specific environments.

Gene Duplication and Protein Evolution in Tick-Host Interactions

Ticks modulate their hosts' defense responses by secreting a biopharmacopiea of hundreds to thousands of proteins and bioactive chemicals into the feeding site (tick-host interface). These molecules and their functions evolved over millions of years as ticks adapted to blood-feeding, tick lineages diverged, and host-shifts occurred. The evolution of new proteins with new functions is mainly dependent on gene duplication events. Central questions around this are the rates of gene duplication, when they occurred and how new functions evolve after gene duplication. The current review investigates these questions in the light of tick biology and considers the possibilities of ancient genome duplication, lineage specific expansion events, and the role that positive selection played in the evolution of tick protein function. It contrasts current views in tick biology regarding adaptive evolution with the more general view that neutral evolution may account for the majority of biological innovations observed in ticks.

Sialotranscriptomics of Rhipicephalus zambeziensis reveals intricate expression profiles of secretory proteins and suggests tight temporal transcriptional regulation during blood-feeding

BACKGROUND

Ticks secrete a diverse mixture of secretory proteins into the host to evade its immune response and facilitate blood-feeding, making secretory proteins attractive targets for the production of recombinant anti-tick vaccines. The largely neglected tick species, Rhipicephalus zambeziensis, is an efficient vector of Theileria parva in southern Africa but its available sequence information is limited. Next generation sequencing has advanced sequence availability for ticks in recent years and has assisted the characterisation of secretory proteins. This study focused on the de novo assembly and annotation of the salivary gland transcriptome of R. zambeziensis and the temporal expression of secretory protein transcripts in female and male ticks, before the onset of feeding and during early and late feeding.

RESULTS

The sialotranscriptome of R. zambeziensis yielded 23,631 transcripts from which 13,584 non-redundant proteins were predicted. Eighty-six percent of these contained a predicted start and stop codon and were estimated to be putatively full-length proteins. A fifth (2569) of the predicted proteins were annotated as putative secretory proteins and explained 52% of the expression in the transcriptome. Expression analyses revealed that 2832 transcripts were differentially expressed among feeding time points and 1209 between the tick sexes. The expression analyses further indicated that 57% of the annotated secretory protein transcripts were differentially expressed. Dynamic expression profiles of secretory protein transcripts were observed during feeding of female ticks. Whereby a number of transcripts were upregulated during early feeding, presumably for feeding site establishment and then during late feeding, 52% of these were downregulated, indicating that transcripts were required at specific feeding stages. This suggested that secretory proteins are under stringent transcriptional regulation that fine-tunes their expression in salivary glands during feeding. No open reading frames were predicted for 7947 transcripts. This class represented 17% of the differentially expressed transcripts, suggesting a potential transcriptional regulatory function of long non-coding RNA in tick blood-feeding.

CONCLUSIONS

The assembled sialotranscriptome greatly expands the sequence availability of R. zambeziensis, assists in our understanding of the transcription of secretory proteins during blood-feeding and will be a valuable resource for future vaccine candidate selection.

Protease Inhibitors in Tick Saliva: The Role of Serpins and Cystatins in Tick-host-Pathogen Interaction

The publication of the first tick sialome (salivary gland transcriptome) heralded a new era of research of tick protease inhibitors, which represent important constituents of the proteins secreted via tick saliva into the host. Three major groups of protease inhibitors are secreted into saliva: Kunitz inhibitors, serpins, and cystatins. Kunitz inhibitors are anti-hemostatic agents and tens of proteins with one or more Kunitz domains are known to block host coagulation and/or platelet aggregation. Serpins and cystatins are also anti-hemostatic effectors, but intriguingly, from the translational perspective, also act as pluripotent modulators of the host immune system. Here we focus especially on this latter aspect of protease inhibition by ticks and describe the current knowledge and data on secreted salivary serpins and cystatins and their role in tick-host-pathogen interaction triad. We also discuss the potential therapeutic use of tick protease inhibitors.

Serine Protease Inhibitors in Ticks: An Overview of Their Role in Tick Biology and Tick-Borne Pathogen Transmission

New tick and tick-borne pathogen control approaches that are both environmentally sustainable and which provide broad protection are urgently needed. Their development, however, will rely on a greater understanding of tick biology, tick-pathogen, and tick-host interactions. The recent advances in new generation technologies to study genomes, transcriptomes, and proteomes has resulted in a plethora of tick biomacromolecular studies. Among these, many enzyme inhibitors have been described, notably serine protease inhibitors (SPIs), whose importance in various tick biological processes is only just beginning to be fully appreciated. Among the multiple active substances secreted during tick feeding, SPIs have been shown to be directly involved in regulation of inflammation, blood clotting, wound healing, vasoconstriction and the modulation of host defense mechanisms. In light of these activities, several SPIs were examined and were experimentally confirmed to facilitate tick pathogen transmission. In addition, to prevent coagulation of the ingested blood meal within the tick alimentary canal, SPIs are also involved in blood digestion and nutrient extraction from the meal. The presence of SPIs in tick hemocytes and their involvement in tick innate immune defenses have also been demonstrated, as well as their implication in hemolymph coagulation and egg development. Considering the involvement of SPIs in multiple crucial aspects of tick-host-pathogen interactions, as well as in various aspects of the tick parasitic lifestyle, these molecules represent highly suitable and attractive targets for the development of effective tick control strategies. Here we review the current knowledge regarding this class of inhibitors in tick biology and tick-borne pathogen transmission, and their potential as targets for future tick control trials.

Functional diversity of secreted cestode Kunitz proteins: Inhibition of serine peptidases and blockade of cation channels

We previously reported a multigene family of monodomain Kunitz proteins from Echinococcus granulosus (EgKU-1-EgKU-8), and provided evidence that some EgKUs are secreted by larval worms to the host interface. In addition, functional studies and homology modeling suggested that, similar to monodomain Kunitz families present in animal venoms, the E. granulosus family could include peptidase inhibitors as well as channel blockers. Using enzyme kinetics and whole-cell patch-clamp, we now demonstrate that the EgKUs are indeed functionally diverse. In fact, most of them behaved as high affinity inhibitors of either chymotrypsin (EgKU-2-EgKU-3) or trypsin (EgKU-5-EgKU-8). In contrast, the close paralogs EgKU-1 and EgKU-4 blocked voltage-dependent potassium channels (Kv); and also pH-dependent sodium channels (ASICs), while showing null (EgKU-1) or marginal (EgKU-4) peptidase inhibitory activity. We also confirmed the presence of EgKUs in secretions from other parasite stages, notably from adult worms and metacestodes. Interestingly, data from genome projects reveal that at least eight additional monodomain Kunitz proteins are encoded in the genome; that particular EgKUs are up-regulated in various stages; and that analogous Kunitz families exist in other medically important cestodes, but not in trematodes. Members of this expanded family of secreted cestode proteins thus have the potential to block, through high affinity interactions, the function of host counterparts (either peptidases or cation channels) and contribute to the establishment and persistence of infection. From a more general perspective, our results confirm that multigene families of Kunitz inhibitors from parasite secretions and animal venoms display a similar functional diversity and thus, that host-parasite co-evolution may also drive the emergence of a new function associated with the Kunitz scaffold.

Identification and anticoagulant activity of a novel Kunitz-type protein HA11 from the salivary gland of the tick Hyalomma asiaticum

Kunitz/bovine pancreatic trypsin inhibitor proteins are abundant in the salivary glands of ticks and perform multiple functions in blood feeding, including inhibiting blood coagulation, regulating host blood supply and disrupting host angiogenesis. In this study, we identified a novel gene designated HA11 (Hyalomma asiaticum 11 kDa protein) from the salivary gland of the tick H. asiaticum. HA11 is encoded by a gene with an open reading frame of 306 bp that is translated into a deduced 101 amino acid 11 kDa protein that shares 27% sequence identity with a Kunitz-like protease inhibitor precursor in Amblyomma variegatum. Bioinformatic analysis confirmed HA11 as a member of the Kunitz-type family of inhibitors. Real time-PCR detected HA11 mRNA transcripts in tick larvae and nymphae stages, with levels highest in salivary gland tissue, and transcription was induced by blood feeding. HA11 anticoagulant activity was demonstrated by its ability to delay normal clotting of rabbit plasma in an activated partial thromboplastin time assay. Furthermore, RNA interference confirmed that HA11 influences H. asiaticum development and blood feeding, and the recombinant protein exerted low hemolytic activity. These results suggest HA11 is a novel Kunitz-type anticoagulant protein involved in tick blood feeding that may have potential as an anticoagulant drug or vaccine.

Strong and widespread action of site-specific positive selection in the snake venom Kunitz/BPTI protein family

S1 family of serine peptidases is the largest family of peptidases. They are specifically inhibited by the Kunitz/BPTI inhibitors. Kunitz domain is characterized by the compact 3D structure with the most important inhibitory loops for the inhibition of S1 peptidases. In the present study we analysed the action of site-specific positive selection and its impact on the structurally and functionally important parts of the snake venom Kunitz/BPTI family of proteins. By using numerous models we demonstrated the presence of large numbers of site-specific positively selected sites that can reach between 30-50% of the Kunitz domain. The mapping of the positively selected sites on the 3D model of Kunitz/BPTI inhibitors has shown that these sites are located in the inhibitory loops 1 and 2, but also in the Kunitz scaffold. Amino acid replacements have been found exclusively on the surface, and the vast majority of replacements are causing the change of the charge. The consequence of these replacements is the change in the electrostatic potential on the surface of the Kunitz/BPTI proteins that may play an important role in the precise targeting of these inhibitors into the active site of S1 family of serine peptidases.

TSGP4 from Ornithodoros moubata: molecular cloning, phylogenetic analysis and vaccine efficacy of a new member of the lipocalin clade of cysteinyl leukotriene scavengers

Recently obtained evidence indicated that an orthologue of the O. savignyi TSGP4 salivary lipocalin was present in the saliva of O. moubata. TSGP4 is known to act as a cysteinyl leukotrienes scavenger helping in the prevention of inflammation and oedema at the tick bite site. Since this function seems to be crucial for successful tick feeding, the novel O. moubata TSGP4 turned into a potential vaccine target. The purposes of the current work were: (i) to clone and characterize the O. moubata TSGP4 and, (ii) to produce it as recombinant to evaluate its protective efficacy as vaccine antigen. The results of these experiments indicated that the O. moubata TSGP4 shows high sequence and structural identity with the O. savignyi orthologue suggesting identical function in the physiology of the tick-host relationship. The mature native TSGP4 is not immunogenic when it is inoculated to host with tick saliva during feeding, but host vaccination with the recombinant protein TSGP4 in Freund's adjuvants induced strong humoral immune responses that recognized both the recombinant and native TSGP4 and protected the host with a 14.1% efficacy. So, the O. moubata TSGP4 can be considered a silent salivary antigen; however, in the light of the current results, its inclusion in the current repertory of protective antigens to be targeted by anti-tick vaccines could be controversial.

A sodium channel inhibitor ISTX-I with a novel structure provides a new hint at the evolutionary link between two toxin folds

Members of arachnida, such as spiders and scorpions, commonly produce venom with specialized venom glands, paralyzing their prey with neurotoxins that specifically target ion channels. Two well-studied motifs, the disulfide-directed hairpin (DDH) and the inhibitor cystine knot motif (ICK), are both found in scorpion and spider toxins. As arachnids, ticks inject a neurotoxin-containing cocktail from their salivary glands into the host to acquire a blood meal, but peptide toxins acting on ion channels have not been observed in ticks. Here, a new neurotoxin (ISTX-I) that acts on sodium channels was identified from the hard tick Ixodes scapularis and characterized. ISTX-I exhibits a potent inhibitory function with an IC₅₀ of 1.6 μM for sodium channel Nav1.7 but not other sodium channel subtypes. ISTX-I adopts a novel structural fold and is distinct from the canonical ICK motif. Analysis of the ISTX-I, DDH and ICK motifs reveals that the new ISTX-I motif might be an intermediate scaffold between DDH and ICK, and ISTX-I is a clue to the evolutionary link between the DDH and ICK motifs. These results provide a glimpse into the convergent evolution of neurotoxins from predatory and blood-sucking arthropods.

Extensive polymorphism of Ra86 genes in field populations of Rhipicephalus appendiculatus from Kenya

Commercial vaccines based on recombinant forms of the Bm86 tick gut antigen are used to control the southern cattle tick, Rhipicephalus microplus, a 1-host species, in Australia and Latin America. We describe herein sequence polymorphism in genes encoding Ra86 homologues of Bm86 in the brown ear tick, Rhipicephalus appendiculatus, isolated from four Kenyan field populations and one laboratory colony. Sequencing of 19 Ra86 sequences defined two alleles differentiated by indels, encoding 693 amino acids (aa) and 654 aa respectively, from the Muguga laboratory reference strain. Ra86 sequences were also determined from gut cDNA from four field populations of R. appendiculatus collected in different livestock production systems in Kenya. Analysis of approximately 20 Ra86 sequences from each of the four field sites in central and Western Kenya; Makuyu, Kiambu, Kakamega and Uasin Gishu, revealed three additional size types differentiated by 39-49 amino acid indels resulting in a total of 5 indel-defined genotypes. The 693 aa type 5 was isolated only from the laboratory tick stock; genotypes 1, 2 and 3 were identified in ticks from the four Kenyan field sites and appeared to be derivatives of the shorter RA86 genotype found in Muguga laboratory stock genotype 4. By contrast no large indels have yet been observed between R. microplus sequences from Australia, South America or Africa. Evidence that selection contributes to the observed sequence variation was provided by analysis of ratio of synonymous and non-synonymous substitutions and application of the selective neutrality and neutral evolution tests to the primary data. Phylogenetic analysis clustered sequences from all Ra86 size types and Bm86, into four major clades based on amino acid substitutions, but there was no evidence that these groupings correlated with geographical separation of R. appendiculatus populations.

All For One and One For All on the Tick-Host Battlefield

The saliva of ixodid ticks contains a mixture of bioactive molecules that target a wide spectrum of host defense mechanisms to allow ticks to feed on the vertebrate host for several days. Tick salivary proteins cluster in multigenic protein families, and individual family members display redundancy and pluripotency in their action to ameliorate or evade host immune responses. It is now clear that members of different protein families can target the same cellular or molecular pathway of the host physiological response to tick feeding. We present and discuss our hypothesis that redundancy and pluripotency evolved in tick salivary immunomodulators to evade immune recognition by the host while retaining the immunomodulatory potential of their saliva.

A shrimp pacifastin light chain-like inhibitor: molecular identification and role in the control of the prophenoloxidase system

Pacifastin is a recently classified family of serine proteinase inhibitors that play essential roles in various biological processes, including in the regulation of the melanization cascade. Here, a novel pacifastin-related gene, termed PmPacifastin-like, was identified from a reverse suppression subtractive hybridization (SSH) cDNA library created from hemocytes of the prophenoloxidase PmproPO1/2 co-silenced black tiger shrimp Penaeus monodon. The full-length sequences of PmPacifastin-like and its homologue LvPacifastin-like from the Pacific white shrimp Litopenaeus vannamei were determined. Sequence analysis revealed that both sequences contained thirteen conserved pacifastin light chain domains (PLDs), followed by two putative kunitz domains. Expression analysis demonstrated that the PmPacifastin-like transcript was expressed in all tested shrimp tissues and larval developmental stages, and its expression responded to Vibrio harveyi challenge. To gain insight into the functional roles of PmPacifastin-like protein, the in vivo RNA interference experiment was employed; the results showed that PmPacifastin-like depletion strongly increased PO activity. Interestingly, suppression of PmPacifastin-like also down-regulated the expression of the proPO-activating enzyme PmPPAE2 transcript; the PmPacifastin-like transcript was down-regulated after the PmproPO1/2 transcripts were silenced. Taken together, these results suggest that PmPacifastin-like is important in the shrimp proPO system and may play an essential role in shrimp immune defense against bacterial infection. These results also expand the knowledge of how pacifastin-related protein participates in the negative regulation of the proPO system in shrimp.

Sialomes and Mialomes: A Systems-Biology View of Tick Tissues and Tick-Host Interactions

Tick saliva facilitates tick feeding and infection of the host. Gene expression analysis of tick salivary glands and other tissues involved in host-pathogen interactions has revealed a wide range of bioactive tick proteins. Transcriptomic analysis has been a milestone in the field and has recently been enhanced by next-generation sequencing (NGS). Furthermore, the application of quantitative proteomics to ticks with unknown genomes has provided deeper insights into the molecular mechanisms underlying tick hematophagy, pathogen transmission, and tick-host-pathogen interactions. We review current knowledge on the transcriptomics and proteomics of tick tissues from a systems-biology perspective and discuss future challenges in the field.

Unraveling dual feeding associated molecular complexity of salivary glands in the mosquito Anopheles culicifacies

Mosquito salivary glands are well known to facilitate meal acquisition, however the fundamental question on how adult female salivary gland manages molecular responses during sugar versus blood meal uptake remains unanswered. To investigate these responses, we analyzed a total of 58.5 million raw reads generated from two independent RNAseq libraries of the salivary glands collected from 3–4 day-old sugar and blood fed Anopheles culicifacies mosquitoes. Comprehensive functional annotation analysis of 10,931 contigs unraveled that salivary glands may encode diverse nature of proteins in response to distinct physiological feeding status. Digital gene expression analysis and PCR validation indicated that first blood meal significantly alters the molecular architecture of the salivary glands. Comparative microscopic analysis also revealed that first blood meal uptake not only causes an alteration of at least 12–22% of morphological features of the salivary glands but also results in cellular changes e.g. apoptosis, confirming together that adult female salivary glands are specialized organs to manage meal specific responses. Unraveling the underlying mechanism of mosquito salivary gene expression, controlling dual feeding associated responses may provide a new opportunity to control vector borne diseases.

The venomous cocktail of the vampire snail Colubraria reticulata (Mollusca, Gastropoda)

Background

Hematophagy arose independently multiple times during metazoan evolution, with several lineages of vampire animals particularly diversified in invertebrates. However, the biochemistry of hematophagy has been studied in a few species of direct medical interest and is still underdeveloped in most invertebrates, as in general is the study of venom toxins. In cone snails, leeches, arthropods and snakes, the strong target specificity of venom toxins uniquely aligns them to industrial and academic pursuits (pharmacological applications, pest control etc.) and provides a biochemical tool for studying biological activities including cell signalling and immunological response. Neogastropod snails (cones, oyster drills etc.) are carnivorous and include active predators, scavengers, grazers on sessile invertebrates and hematophagous parasites; most of them use venoms to efficiently feed. It has been hypothesized that trophic innovations were the main drivers of rapid radiation of Neogastropoda in the late Cretaceous.

We present here the first molecular characterization of the alimentary secretion of a non-conoidean neogastropod, Colubraria reticulata. Colubrariids successfully feed on the blood of fishes, throughout the secretion into the host of a complex mixture of anaesthetics and anticoagulants. We used a NGS RNA-Seq approach, integrated with differential expression analyses and custom searches for putative secreted feeding-related proteins, to describe in detail the salivary and mid-oesophageal transcriptomes of this Mediterranean vampire snail, with functional and evolutionary insights on major families of bioactive molecules.

Results

A remarkably low level of overlap was observed between the gene expression in the two target tissues, which also contained a high percentage of putatively secreted proteins when compared to the whole body. At least 12 families of feeding-related proteins were identified, including: 1) anaesthetics, such as ShK Toxin-containing proteins and turripeptides (ion-channel blockers), Cysteine-rich secretory proteins (CRISPs), Adenosine Deaminase (ADA); 2) inhibitors of primary haemostasis, such as novel vWFA domain-containing proteins, the Ectonucleotide pyrophosphatase/phosphodiesterase family member 5 (ENPP5) and the wasp Antigen-5; 3) anticoagulants, such as TFPI-like multiple Kunitz-type protease inhibitors, Peptidases S1 (PS1), CAP/ShKT domain-containing proteins, Astacin metalloproteases and Astacin/ShKT domain-containing proteins; 4) additional proteins, such the Angiotensin-Converting Enzyme (ACE: vasopressive) and the cytolytic Porins.

Conclusions

Colubraria feeding physiology seems to involve inhibitors of both primary and secondary haemostasis, anaesthetics, a vasoconstrictive enzyme to reduce feeding time and tissue-degrading proteins such as Porins and Astacins. The complexity of Colubraria venomous cocktail and the divergence from the arsenal of the few neogastropods studied to date (mostly conoideans) suggest that biochemical diversification of neogastropods might be largely underestimated and worth of extensive investigation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1648-4) contains supplementary material, which is available to authorized users.

Prediction of Kunitz ion channel effectors and protease inhibitors from the Ixodes ricinus sialome

In the next generation sequencing era we are encountering hundreds of thousands of sequences from specific organisms. Such massive data must be accurately classified both functionally and structurally. Determining appropriate sequences with a specific function from next generation sequencing, however, is a daunting experimental task. A recent salivary gland transcriptome from the hard tick Ixodes ricinus, a European disease vector, has been made publicly available. Among the protein families sequenced by the salivary gland transcriptome of I. ricinus, the Kunitz-domain is one of the highly represented protein families. Thus far, recent tick transciptomes solely classify (computationally) Kunitz sequences as putative serine protease inhibitors. We present here a novel method using a machine-learning algorithm to "fish" for candidate ion-channel effectors and loss of serine protease inhibitor function within the Kunitz-domain protein family of the I. ricinus salivary gland transcriptome. The models, data and scripts used in this work are available online from http://life.bsc.es/pid/web/imoal/kunitz-classification.html.

Are ticks venomous animals?

INTRODUCTION

As an ecological adaptation venoms have evolved independently in several species of Metazoa. As haematophagous arthropods ticks are mainly considered as ectoparasites due to directly feeding on the skin of animal hosts. Ticks are of major importance since they serve as vectors for several diseases affecting humans and livestock animals. Ticks are rarely considered as venomous animals despite that tick saliva contains several protein families present in venomous taxa and that many Ixodida genera can induce paralysis and other types of toxicoses. Tick saliva was previously proposed as a special kind of venom since tick venom is used for blood feeding that counteracts host defense mechanisms. As a result, the present study provides evidence to reconsider the venomous properties of tick saliva.

RESULTS

Based on our extensive literature mining and in silico research, we demonstrate that ticks share several similarities with other venomous taxa. Many tick salivary protein families and their previously described functions are homologous to proteins found in scorpion, spider, snake, platypus and bee venoms. This infers that there is a structural and functional convergence between several molecular components in tick saliva and the venoms from other recognized venomous taxa. We also highlight the fact that the immune response against tick saliva and venoms (from recognized venomous taxa) are both dominated by an allergic immunity background. Furthermore, by comparing the major molecular components of human saliva, as an example of a non-venomous animal, with that of ticks we find evidence that ticks resemble more venomous than non-venomous animals. Finally, we introduce our considerations regarding the evolution of venoms in Arachnida.

CONCLUSIONS

Taking into account the composition of tick saliva, the venomous functions that ticks have while interacting with their hosts, and the distinguishable differences between human (non-venomous) and tick salivary proteins, we consider that ticks should be referred to as venomous ectoparasites.

IrSPI, a tick serine protease inhibitor involved in tick feeding and Bartonella henselae infection

Ixodes ricinus is the most widespread and abundant tick in Europe, frequently bites humans, and is the vector of several pathogens including those responsible for Lyme disease, Tick-Borne Encephalitis, anaplasmosis, babesiosis and bartonellosis. These tick-borne pathogens are transmitted to vertebrate hosts via tick saliva during blood feeding, and tick salivary gland (SG) factors are likely implicated in transmission. In order to identify such tick factors, we characterized the transcriptome of female I. ricinus SGs using next generation sequencing techniques, and compared transcriptomes between Bartonella henselae-infected and non-infected ticks. High-throughput sequencing of I. ricinus SG transcriptomes led to the generation of 24,539 isotigs. Among them, 829 and 517 transcripts were either significantly up- or down-regulated respectively, in response to bacterial infection. Searches based on sequence identity showed that among the differentially expressed transcripts, 161 transcripts corresponded to nine groups of previously annotated tick SG gene families, while the others corresponded to genes of unknown function. Expression patterns of five selected genes belonging to the BPTI/Kunitz family of serine protease inhibitors, the tick salivary peptide group 1 protein, the salp15 super-family, and the arthropod defensin family, were validated by qRT-PCR. IrSPI, a member of the BPTI/Kunitz family of serine protease inhibitors, showed the highest up-regulation in SGs in response to Bartonella infection. IrSPI silencing impaired tick feeding, as well as resulted in reduced bacterial load in tick SGs. This study provides a comprehensive analysis of I. ricinus SG transcriptome and contributes significant genomic information about this important disease vector. This in-depth knowledge will enable a better understanding of the molecular interactions between ticks and tick-borne pathogens, and identifies IrSPI, a candidate to study now in detail to estimate its potentialities as vaccine against the ticks and the pathogens they transmit.

Understanding the evolutionary structural variability and target specificity of tick salivary Kunitz peptides using next generation transcriptome data

BACKGROUND

Ticks are blood-sucking arthropods and a primary function of tick salivary proteins is to counteract the host's immune response. Tick salivary Kunitz-domain proteins perform multiple functions within the feeding lesion and have been classified as venoms; thereby, constituting them as one of the important elements in the arms race with the host. The two main mechanisms advocated to explain the functional heterogeneity of tick salivary Kunitz-domain proteins are gene sharing and gene duplication. Both do not, however, elucidate the evolution of the Kunitz family in ticks from a structural dynamic point of view. The Red Queen hypothesis offers a fruitful theoretical framework to give a dynamic explanation for host-parasite interactions. Using the recent salivary gland Ixodes ricinus transcriptome we analyze, for the first time, single Kunitz-domain encoding transcripts by means of computational, structural bioinformatics and phylogenetic approaches to improve our understanding of the structural evolution of this important multigenic protein family.

RESULTS

Organizing the I. ricinus single Kunitz-domain peptides based on their cysteine motif allowed us to specify a putative target and to relate this target specificity to Illumina transcript reads during tick feeding. We observe that several of these Kunitz peptide groups vary in their translated amino acid sequence, secondary structure, antigenicity, and intrinsic disorder, and that the majority of these groups are subject to a purifying (negative) selection. We finalize by describing the evolution and emergence of these Kunitz peptides. The overall interpretation of our analyses discloses a rapidly emerging Kunitz group with a distinct disulfide bond pattern from the I. ricinus salivary gland transcriptome.

CONCLUSIONS

We propose a model to explain the structural and functional evolution of tick salivary Kunitz peptides that we call target-oriented evolution. Our study reveals that combining analytical approaches (transcriptomes, computational, bioinformatics and phylogenetics) improves our understanding of the biological functions of important salivary gland mediators during tick feeding.

De novo Ixodes ricinus salivary gland transcriptome analysis using two next-generation sequencing methodologies

Tick salivary gland (SG) proteins possess powerful pharmacologic properties that facilitate tick feeding and pathogen transmission. For the first time, SG transcriptomes of Ixodes ricinus, an important disease vector for humans and animals, were analyzed using next-generation sequencing. SGs were collected from different tick life stages fed on various animal species, including cofeeding of nymphs and adults on the same host. Four cDNA samples were sequenced, discriminating tick SG transcriptomes of early- and late-feeding nymphs or adults. In total, 441,381,454 pyrosequencing reads and 67,703,183 Illumina reads were assembled into 272,220 contigs, of which 34,560 extensively annotated coding sequences are disclosed; 8686 coding sequences were submitted to GenBank. Overall, 13% of contigs were classified as secreted proteins that showed significant differences in the transcript representation among the 4 SG samples, including high numbers of sample-specific transcripts. Detailed phylogenetic reconstructions of two relatively abundant SG-secreted protein families demonstrated how this study improves our understanding of the molecular evolution of hematophagy in arthropods. Our data significantly increase the available genomic information for I. ricinus and form a solid basis for future tick genome/transcriptome assemblies and the functional analysis of effectors that mediate the feeding physiology and parasite-vector interaction of I. ricinus.

Characterization of Ixophilin, a thrombin inhibitor from the gut of Ixodes scapularis

Ixodes scapularis, the black-legged tick, vectors several human pathogens including Borrelia burgdorferi, the agent of Lyme disease in North America. Pathogen transmission to the vertebrate host occurs when infected ticks feed on the mammalian host to obtain a blood meal. Efforts to understand how the tick confronts host hemostatic mechanisms and imbibes a fluid blood meal have largely focused on the anticoagulation strategies of tick saliva. The blood meal that enters the tick gut remains in a fluid state for several days during the process of feeding, and the role of the tick gut in maintaining the blood-meal fluid is not understood. We now demonstrate that the tick gut produces a potent inhibitor of thrombin, a key enzyme in the mammalian coagulation cascade. Chromatographic fractionation of engorged tick gut proteins identified one predominant thrombin inhibitory activity associated with an approximately 18 kDa protein, henceforth referred to as Ixophilin. The ixophilin gene was preferentially transcribed in the guts of feeding nymphs. Expression began after 24 hours of feeding, coincident with the flow of host blood into the tick gut. Immunity against Ixophilin delayed tick feeding, and decreased feeding efficiency significantly. Surprisingly, immunity against Ixophilin resulted in increased Borrelia burgdorferi transmission to the host, possibly due to delayed feeding and increased transmission opportunity. These observations illuminate the potential drawbacks of targeting individual tick proteins in a functional suite. They also underscore the need to identify the “anticoagulome” of the tick gut, and to prioritize a critical subset of anticoagulants that could be targeted to efficiently thwart tick feeding, and block pathogen transmission to the vertebrate host.

Tryptogalinin is a tick Kunitz serine protease inhibitor with a unique intrinsic disorder

BACKGROUND

A salivary proteome-transcriptome project on the hard tick Ixodes scapularis revealed that Kunitz peptides are the most abundant salivary proteins. Ticks use Kunitz peptides (among other salivary proteins) to combat host defense mechanisms and to obtain a blood meal. Most of these Kunitz peptides, however, remain functionally uncharacterized, thus limiting our knowledge about their biochemical interactions.

RESULTS

We discovered an unusual cysteine motif in a Kunitz peptide. This peptide inhibits several serine proteases with high affinity and was named tryptogalinin due to its high affinity for β-tryptase. Compared with other functionally described peptides from the Acari subclass, we showed that tryptogalinin is phylogenetically related to a Kunitz peptide from Rhipicephalus appendiculatus, also reported to have a high affinity for β-tryptase. Using homology-based modeling (and other protein prediction programs) we were able to model and explain the multifaceted function of tryptogalinin. The N-terminus of the modeled tryptogalinin is detached from the rest of the peptide and exhibits intrinsic disorder allowing an increased flexibility for its high affinity with its inhibiting partners (i.e., serine proteases).

CONCLUSIONS

By incorporating experimental and computational methods our data not only describes the function of a Kunitz peptide from Ixodes scapularis, but also allows us to hypothesize about the molecular basis of this function at the atomic level.

The "Vampirome": Transcriptome and proteome analysis of the principal and accessory submaxillary glands of the vampire bat Desmodus rotundus, a vector of human rabies

Vampire bats are notorious for being the sole mammals that strictly feed on fresh blood for their survival. While their saliva has been historically associated with anticoagulants, only one antihemostatic (plasminogen activator) has been molecularly and functionally characterized. Here, RNAs from both principal and accessory submaxillary (submandibular) salivary glands of Desmodus rotundus were extracted, and ~200 million reads were sequenced by Illumina. The principal gland was enriched with plasminogen activators with fibrinolytic properties, members of lipocalin and secretoglobin families, which bind prohemostatic prostaglandins, and endonucleases, which cleave neutrophil-derived procoagulant NETs. Anticoagulant (tissue factor pathway inhibitor, TFPI), vasodilators (PACAP and C-natriuretic peptide), and metalloproteases (ADAMTS-1) were also abundantly expressed. Members of the TSG-6 (anti-inflammatory), antigen 5/CRISP, and CCL28-like (antimicrobial) protein families were also sequenced. Apyrases (which remove platelet agonist ADP), phosphatases (which degrade procoagulant polyphosphates), and sphingomyelinase were found at lower transcriptional levels. Accessory glands were enriched with antimicrobials (lysozyme, defensin, lactotransferrin) and protease inhibitors (TIL-domain, cystatin, Kazal). Mucins, heme-oxygenase, and IgG chains were present in both glands. Proteome analysis by nano LC-MS/MS confirmed that several transcripts are expressed in the glands. The database presented herein is accessible online at http://exon.niaid.nih.gov/transcriptome/D_rotundus/Supplemental-web.xlsx. These results reveal that bat saliva emerges as a novel source of modulators of vascular biology.

BIOLOGICAL SIGNIFICANCE

Vampire bat saliva emerges as a novel source of antihemostatics which modulate several aspects of vascular biology.

Evolution of the tissue factor pathway inhibitor-like Kunitz domain-containing protein family in Rhipicephalus microplus

One of the principle mechanisms utilised by ticks to obtain a blood meal is the subversion of the host's haemostatic response. This is achieved through the secretion of saliva containing anti-haemostatic proteins into the feeding lesion. Lineage-specific expansion of predicted secretory protein families have been observed in all previously studied ticks and occurred in response to adaptation to a blood-feeding environment. Of these, the predominant families are common between both hard and soft ticks. One of these families, namely the Kunitz domain-containing protein family, includes proven tissue factor pathway inhibitor-like (TFPI-like) anti-haemostatics such as ixolaris and penthalaris that play a crucial role during tick feeding. Although Kunitz-type proteins have been found in Rhipicephalus microplus, the TFPI-like Kunitz protein family has not yet been studied. We report a comprehensive search for TFPI-like Kunitz domain-containing proteins in R. microplus expressed sequence tag libraries, resulting in the identification of 42 homologues. The homologues were bioinformatically and phylogenetically studied, including the application of an intensive Bayesian Markov Chain Monte Carlo (MCMC) analysis of the individual Kunitz domain nucleotide sequences. We show that the R. microplus TFPI-like Kunitz protein family groups into two main clades that presumably underwent ancient duplication, which indicates that a whole genome duplication event occurred at least 150 million years ago. Evidence for recent and ancient gene and domain duplication events was also found. Furthermore, the divergence times of the various tick lineages estimated in this paper correspond with those presented in previous studies. The elucidation of this large protein family's evolution within R. microplus adds to current knowledge of this economically important tick.