# Immunomodulatory Effects of Amblyomma variegatum Saliva on Bovine Cells: Characterization of Cellular Responses and Identification of Molecular Determinants

V. Rodrigues, B. Fernandez, A. Vercoutere, L. Chamayou, A. Andersen, Oana Vigy, Edith Demettre, Martial Séveno, Rosalie Aprelon, K. Giraud-Girard, et al.

► To cite this version:

V. Rodrigues, B. Fernandez, A. Vercoutere, L. Chamayou, A. Andersen, et al., Immunomodulatory Effects of Amblyomma variegatum Saliva on Bovine Cells: Characterization of Cellular Responses and Identification of Molecular Determinants. Frontiers in Cellular and Infection Microbiology, 2017, 7, pp.521. 10.3389/fcimb.2017.00521. [hal-02073207](https://hal-umontpellier.fr/hal-02073207)

HAL Id: hal-02073207

https://hal.umontpellier.fr/hal-02073207

Submitted on 12 Mar 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# Immunomodulatory Effects of *Amblyomma variegatum* Saliva on Bovine Cells: Characterization of Cellular Responses and Identification of Molecular Determinants

Valérie Rodrigues<fnref n="1" />,², Bernard Fernandez<fnref n="1" />,², Arthur Vercoutere<fnref n="1" />,², Léo Chamayou<fnref n="1" />,², Alexandre Andersen<fnref n="1" />,², Oana Vigy³, Edith Demettre⁴, Martial Seveno⁴, Rosalie Aprelon<fnref n="1" />,⁵, Ken Giraud-Girard<fnref n="1" />,⁵, Frédéric Stachurski<fnref n="1" />,², Etienne Loire<fnref n="1" />,², Nathalie Vachiéry<fnref n="1" />,²,⁵ and Philippe Holzmuller<fnref n="1" />,²*

## OPEN ACCESS

**Edited by:**

Ard Menzo Nijhof,
Freie Universität Berlin, Germany

**Reviewed by:**

Sukanya Narasimhan,
Yale School of Medicine, Yale
University, United States
Ashu Sharma,
University at Buffalo, United States

*Correspondence:*

Philippe Holzmuller
philippe.holzmuller@cirad.fr

Received: 30 September 2017

Accepted: 07 December 2017

Published: 04 January 2018

**Citation:**

Rodrigues V, Fernandez B, Vercoutere A, Chamayou L, Andersen A, Vigy O, Demettre E, Seveno M, Aprelon R, Giraud-Girard K, Stachurski F, Loire E, Vachiéry N and Holzmuller P (2018) Immunomodulatory Effects of Amblyomma variegatum Saliva on Bovine Cells: Characterization of Cellular Responses and Identification of Molecular Determinants. *Front. Cell. Infect. Microbiol.* 7:521. doi: 10.3389/fcimb.2017.00521

The tropical bont tick, *Amblyomma variegatum*, is a tick species of veterinary importance and is considered as one of major pest of ruminants in Africa and in the Caribbean. It causes direct skin lesions, transmits heartwater, and reactivates bovine dermatophilosis. Tick saliva is reported to affect overall host responses through immunomodulatory and anti-inflammatory molecules, among other bioactive molecules. The general objective of this study was to better understand the role of saliva in interaction between the *Amblyomma* tick and the host using cellular biology approaches and proteomics, and to discuss its impact on disease transmission and/or activation. We first focused on the immuno-modulating effects of semi-fed *A. variegatum* female saliva on bovine peripheral blood mononuclear cells (PBMC) and monocyte-derived macrophages *in vitro*. We analyzed its immuno-suppressive properties by measuring the effect of saliva on PBMC proliferation, and observed a significant decrease in ConA-stimulated PBMC lymphoproliferation. We then studied the effect of saliva on bovine macrophages using flow cytometry to analyze the expression of MHC-II and co-stimulation molecules (CD40, CD80, and CD86) and by measuring the production of nitric oxide (NO) and pro- or anti-inflammatory cytokines. We observed a significant decrease in the expression of MHC-II, CD40, and CD80 molecules, associated with decreased levels of IL-12-p40 and TNF-α and increased level of IL-10, which could explain the saliva-induced modulation of NO. To elucidate these immunomodulatory effects, crude saliva proteins were analyzed using proteomics with an Orbitrap Elite mass spectrometer. Among the 336 proteins identified in *A. variegatum* saliva, we evidenced bioactive molecules exhibiting anti-inflammatory, immuno-modulatory, and anti-oxidant properties (e.g., serpins, phospholipases A2,

<footnote n="1">¹Centre de Coopération Internationale en Recherche Agronomique pour le Développement, UMR ASTRE "Animal, Santé, Territoire, Risques et Ecosystèmes," Montpellier, France, ²ASTRE, Université de Montpellier (I-MUSE), CIRAD, Institut National de la Recherche Agronomique, Montpellier, France, ³Institut de Génomique Fonctionnelle, Centre National de la Recherche Scientifique, Institut National de la Santé et de la Recherche Médicale, Université de Montpellier, Montpellier, France, ⁴BioCampus Montpellier, Centre National de la Recherche Scientifique, Institut National de la Santé et de la Recherche Médicale, Université de Montpellier, Montpellier, France, ⁵CIRAD, UMR ASTRE, Petit-Bourg, Guadeloupe, France</footnote>
heme lipoprotein). We also characterized an intriguing ubiquitination complex that could be involved in saliva-induced immune modulation of the host. We propose a model for the interaction between *A. variegatum* saliva and host immune cells that could have an effect during tick feeding by favoring pathogen dissemination or activation by reducing the efficiency of host immune response to the corresponding tick-borne diseases.

**Keywords: Amblyomma variegatum, tick saliva, PBMC, immuno-modulation, proteomics**

## INTRODUCTION

*Amblyomma variegatum* is among the most important and widely distributed ticks of tropical livestock, and is the subject of veterinary and public health concerns in Africa and islands in the Indian Ocean and the Caribbean (Stachurski et al., 2010; Bournez et al., 2015). The species is the natural vector of *Ehrlichia ruminantium*, the Rickettsia that causes heartwater in ruminants (Allsopp, 2010), and can also transmit human pathogens including several species of Rickettsia and viruses (Jongejan and Uilenberg, 2004). Moreover, *A. variegatum* has been proven to reactivate the development of dermatophilosis, a skin disease with major economic impacts caused by *Dermatophilus congolensis* (Martinez et al., 1992).

Generally, at the bite site, ticks are able to inhibit pain and itch, inflammation, hemostasis, and wound healing, but also to modulate the host innate and adaptive immune responses, which consequently favor transmission of infectious agents (Šimo et al., 2017). The coevolution between ticks, hosts and pathogens has led to a balance of conflicts and cooperations between the different actors, which mainly benefit the ticks and pathogens. Nevertheless, increased production of antibodies against molecular determinants of tick saliva could also increase protection against pathogen infection (Wade, 2007; Wikel, 2013; de la Fuente et al., 2016). In tick-host interactions, tick saliva is the source of biologically active molecules that target a wide spectrum of host physiological mechanisms, mainly inhibiting host defense reactions to the benefit of the feeding ticks (Kazimírová and Štibrániová, 2013; Stibrániová et al., 2013). The immediate inflammatory response to the skin injury induces rapid infiltration of leukocytes at the tick bite site, where both resident and infiltrated cells (keratinocytes, endothelial cells, mast cells, dendritic cells, macrophages, and lymphocytes) are activated by direct contact with tick saliva (Brossard and Wikel, 2004; Francischetti et al., 2009; Wikel, 2013). Omics analyses have enabled the comprehensive characterization of the molecular determinants of tick saliva, with a functional involvement of these bioactive molecules in host immune suppression (Kotál et al., 2015; Chmelar et al., 2016a,b).

The saliva of *Amblyomma* ticks has been the subject of extensive biochemical studies since the early 2000s (Karim et al., 2004, 2011; Madden et al., 2004; Mulenga et al., 2013; Garcia et al., 2014; Araujo et al., 2016), but only a few have been devoted to the sialome of *A. variegatum*. These studies revealed new families of tick-exclusive proteins (Knizetova et al., 2006; Koh et al., 2007; Ribeiro et al., 2011), some related to suppression of inflammation (Tian et al., 2016).

In the present study, using crude saliva from *A. variegatum* semi-fed females (i.e., at the peak of salivation), our aim was to (i) explore the impact of saliva *in vitro* on bovine immune cells and (ii) determine the molecular immunomodulators. The originality of our approach was to combine cellular experiments, demonstrating modulation of lymphocyte proliferation and macrophage activation, with proteomic analysis of saliva to identify and characterize the immunomodulatory proteins potentially responsible for the observed biological effects. Consequences on host response to infestation by *A. variegatum* and transmission of pathogens are discussed.

## MATERIALS AND METHODS

### A. variegatum Saliva Production

Ten *A. variegatum* female ticks were engorged on naïve creole goats, in accordance with the experimental procedure of the project in animal experimentation approved by the ethics committee Antilles-Guyane (Project Number 69 registered by Comité National de Réflexion Ethique sur l'Expérimentation Animale). After 5-7 days, the female ticks were detached at a semi-fed stage and injected in the hemocoele with a Hamilton syringe with 10 µl of a 5% (w/v in PBS) pilocarpine solution to stimulate saliva production. A capillary micropipette placed around the tick hypostome was used to collect the saliva from ticks (duration of salivation about 30 min), and fresh crude saliva was aliquoted and immediately stored at -80°C for further use. For cellular tests, three different pools of saliva were obtained from 10 *A. variegatum* semi-fed female ticks. For proteomic analyses, crude saliva, from either individual ticks' saliva or pools of saliva, was processed directly according to the standardized protocol (see below). Protein concentration was measured with a Nanodrop 1000 (ThermoScientific). The concentration of Pilocarpine in the saliva was evaluated using an HPLC-MS/MS method proposed by Ribeiro et al. (2004). Contrary to the pilocarpine concentrations previously reported for *A. americanum* (Ribeiro et al., 2004), we did not detect in our study a concentration of salivary pilocarpine higher than 5 nM, which was set as our pilocarpine control concentration.

### Purification of Bovine PBMCs and Proliferation Assay

Peripheral blood was collected in lithium heparin tubes (BD Biosciences, Le Pont de Claix, France) from three healthy female Jersiaise cattle aged between 4 and 6 years (APAFIS#14442015081310300000). PBMC were isolated by
differential centrifugation over Histopaque 1083® (Sigma Aldrich, St. Quentin, France) and washed twice in calcium and magnesium free Hanks balanced salt solution (Life Technologies, Saint Aubin, France). PBMC were suspended at 2.5 × 10⁶ cells/ml in RPMI-1640 culture medium (Life Technologies, Saint Aubin, France) supplemented with 2 mM L-glutamine (Sigma Aldrich, St. Quentin, France), 5 × 10⁻⁵ M β₂-mercaptoethanol (Sigma Aldrich, St. Quentin, France), 50 µg/ml gentamycin (Sigma Aldrich, St. Quentin, France), and 10% heat-inactivated fetal calf serum (Life Technologies, Saint Aubin, France). One hundred microliters of PBMC (2.5 × 10⁵ cells) were distributed in triplicate in the wells of a 96-well microplate and the cells were stimulated with 2.5 µg/ml of Concanavalin A (ConA, Sigma Aldrich, St. Quentin, France) either alone (positive control corresponding to 100% of stimulation) or in the presence of A. variegatum saliva at concentrations ranging from 0.3 to 37.5 µg/ml for 72 h at 37°C with 5% CO₂. PBMC were also tested with saliva alone, without ConA stimulation, and with pilocarpine at 5 nM. Cell proliferation was measured using the CyQUANT Direct Cell Proliferation Assay kit (Molecular Probes/ThermoFisher Scientific, Paisley, UK) according to the manufacturer's instructions. Briefly, 100 µl of the 2X detection reagent were added to the cells for a further 60 min of incubation at 37°C, 5% CO₂. The fluorescence of samples was then read at 480/535 nm using an Enspire Multimode Plate Reader (Perkin Elmer, Courtaboeuf, France). Fluorescence in samples treated with ConA and saliva was compared with fluorescence measured with ConA alone to calculate the relative PBMC activation percentage. Similarly, fluorescence in unstimulated PBMC treated with saliva was compared with fluorescence measured in unstimulated cells. Experiments were performed using cells from three different animals (1973, 9567, and 9906), each of the three saliva batches being tested separately. Data are presented as relative PBMC activation, means of values obtained for the three batches of saliva (for each concentration) in triplicate experiments for each of the three animals.

## Preparation of Bovine Macrophages and Stimulation Assay

Immuno-magnetic separation of CD14⁺ monocytes from PBMC was carried out using anti-human CD14 antibody (TUK4) coated beads (Miltenyi Biotec, Bergish Gladbach, Germany), according to the manufacturer's instructions. Flow cytometry analysis was carried out on the positively selected populations and confirmed that the purity of the population was >95%. Bovine monocytes were suspended in complete IMDM (Life Technologies, Saint Aubin, France) supplemented with 2 mM L-glutamine (Sigma, Saint-Quentin, France), 10% heat-inactivated fetal calf serum (Life Technologies, Saint-Aubin, France), 5 × 10⁻⁵ M β₂-mercaptoethanol (Sigma, Saint-Quentin, France), and 50 µg/ml gentamicin (Life Technologies, Saint-Aubin, France). One hundred microliters of cells at 1.10⁶ cells/ml were seeded into 96-well plates (for NO production) or 500 µl of cells were seeded into 48-well plates (for production of cytokines and analysis of co-stimulation surface markers). Half the culture

medium was replaced every 3 days by complete IMDM, supplemented with GM-CSF and INF-γ for NO induction in the 96-well plate. After 7 days of culture, more than 95% of cells had the morphological characteristics of macrophages. At day 7, 100 ng/ml of ultrapure LPS from E. coli (InVivoGen, Toulouse, France) was added for 1 h prior to saliva stimulation, alone or in the presence of A. variegatum saliva at concentrations ranging from 7.8125 to 250 µg/ml for an additional 24 h at 37°C with 5% CO₂. Pilocarpine (5 nM) was also tested. Cells were also tested without LPS pre-stimulation, with the different concentrations of saliva or with 5 nM of pilocarpine. Supernatants were collected from the 96-well plate and NO concentration was immediately determined using the Griess method. Supernatants were collected from the 48-well plate for TNF-α, IL-12, and IL-10 quantification by ELISA, according to the protocols of Hope et al. (2002) and Kwong et al. (2002, 2010). Briefly, 5 µg/ml of capture Mabs CC327, CC301, and CC318 for TNF-α, IL-12, and IL-10, respectively (Bio-Rad AbD Serotec, UK) were used as coating antibodies, cultures were left overnight at +4°C. After washing and blocking with 1 mg/ml of casein (Sigma Aldrich, Saint Quentin, France), cell culture supernatants and serial diluted standards were incubated for 1 h at room temperature. Wells were washed and a 8 µg/ml solution of biotinylated detection Mabs (biotinylated-CC328, -CC326, and -CC320, all from Bio-Rad AbD Serotec, UK) were added for 1 h, detected by Streptavidin-HRP (eBiosciences, France) and revealed with TMB liquid substrate system for ELISA (Sigma Aldrich, France). The reaction was stopped with H₂SO₄ 0.5 M and the plates were read at 450 nm. All the results were determined with standard curves obtained with serial dilutions of recombinant proteins (ref. RBOTNFA for TNF-α from Pierce, Rockford, USA: ref. 87464 from AbCam, France for IL-12 and ref. RPO379 for IL-10 from Kingfisher, USA).

Cells in the 48-wells plate were collected and labeled with mouse anti bovine monoclonal antibodies for analysis of co-stimulation surface markers: MHC II (clone J11, kindly provided by ILRI, Kenya), CD40 (clone IL-A156, AbD Serotec, UK), CD80 (clone IL-A159, AbD Serotec, UK), and CD86 (clone IL-A190, AbD Seroteck, UK). Isotypic control was used as non-specific labeling (MCA928, AbD Serotec, UK). Antibodies were revealed using goat anti mouse IgG1-RPE (STAR132PE, AbD Serotec, UK), cells were processed for flow cytometry using a FacsCanto II cytometer (BD Biosciences, USA) and data were analyzed using Diva software (BD Biosciences, USA) for FCS/SSC and FL2 fluorescence intensities. Relative expression in MFI (median of fluorescence intensities) was calculated by comparing stimulated cell (without saliva or with different saliva concentrations) MFI with unstimulated cell MFI. Experiments were performed with cells from three different animals (1973, 9567, and 9906), and tested with a pool of saliva from the three batches.

## Proteomics Analysis of A. variegatum Saliva

A. variegatum saliva proteins (10 µg) were separated on SDS-PAGE gels (12% polyacrylamide, Mini-PROTEAN® TGX™
Precast Gels, Bio-Rad, Hercules USA) and stained with Protein Staining Solution (Euromedex, Souffelweyersheim France). Gel lanes were cut into three continuous gel pieces that were treated independently, in order to isolate majority proteins and to allow identification of low quantity proteins. Gel pieces were first destained with 50 mM triethylammonium bicarbonate (TEABC) and three washes in 100% acetonitrile. After protein reduction (with 10 mM dithiothreitol in 50 mM TEABC at 56°C for 45 min) and alkylation (55 mM iodoacetamide TEABC at room temperature for 30 min) proteins were digested in-gel using trypsin (500 ng/band, Gold, Promega, Madison USA) as previously described (Thouvenot et al., 2008). Digest products were dehydrated in a vacuum centrifuge and reduced to 3 µL. The generated peptides were analyzed by nano-flowHPLC-nanoelectrospray ionization using an Orbitrap Elite mass spectrometer (Thermo Fisher Scientific, Waltham USA) coupled to an Ultimate 3000 HPLC (Thermo Fisher Scientific). Desalting and pre-concentration of samples were performed online on a Pepmap® pre-column (0.3 × 10 mm, Dionex). A gradient consisting of 0-40% B for 60 min and 80% B for 15 min (A = 0.1% formic acid, 2% acetonitrile in water; B = 0.1% formic acid in acetonitrile) at 300 nL/min was used to elute peptides from the capillary reverse-phase column (0.075 × 150 mm, Acclaim Pepmap 100® C18, Thermo Fisher Scientific). Eluted peptides were electro-sprayed online at a voltage of 1.9 kV into an Oribtrap Elite mass spectrometer. A cycle of one full-scan mass spectrum (400-2,000 m/z) at a resolution of 120,000 (at 400 m/z), followed by 20 data-dependent MS/MS spectra was repeated continuously throughout the nanoLC separation. All MS/MS spectra were recorded using normalized collision energy (33%, activation Q 0.25 and activation time 10 ms) with an isolation window of 2 m/z. Data were acquired using Xcalibur software (v 2.2). For all full scan measurements with the Orbitrap detector a lock-mass ion from ambient air (m/z 445.120024) was used as an internal calibrant as described in Olsen et al. (2005). MS data were analyzed using the MaxQuant software package (v 1.5.0.0) as described by Cox and Mann (2008). Tandem mass spectra (MS/MS) were searched for using the Andromeda search engine (Cox et al., 2011) against the UniProtKB database (release 2015_03) for the *A. variegatum* (AMBVA) taxonomy (753 entries), Amblyomma "all species" taxonomy (22,429 entries), and the UniProtKB proteome UP000001555 database for *Ixodes scapularis* (20,473 entries) using the following parameters: enzyme specificity was set as Trypsin/P, and a maximum of two missed cleavages and a mass tolerance of 0.5 Da for fragment ion were applied. A second database of known contaminants provided with the MaxQuant suite was also used. The "match between runs" option was checked. Oxidation (M) was specified as variable modification and carbamidomethyl (C) as fixed modification. Database searches were performed with a mass tolerance of 20 ppm for precursor ion for mass calibration, and with a 4.5 ppm tolerance after calibration. The maximum false peptide and protein discovery rate was set at 0.01. The MaxQuant software generates several output files that contain information about the peptides and proteins identified. The "proteinGroups.txt" file is dedicated to identified proteins: each single row collapses

all proteins that cannot be distinguished based on identified peptides into protein groups. An in-house bioinformatics tool was developed to automatically select a representative protein ID in each protein group. First, proteins with the most identified peptides were isolated in a group called "match group" [proteins from the "Protein IDs" column with the maximum number of "peptides counts (all)"]. For 35% of the remaining match groups in which more than one protein ID existed, the "leading" protein was chosen as the best annotated protein according to the number of gene ontology annotations (retrieval from UniProtKB March 20, 2015) and/or given the following species order preference: *A. variegatum* > *A. maculatum* > *A. cajennense* > *A. parvum* > *A. americanum* > *A. brasiliense* > *A. triste* > *A. geayi* > *A. scutatum* > *A. rotundatum*. The mass spectrometry data along with the identification results were deposited in the ProteomeXchange Consortium via the PRIDE (Martens et al., 2005) partner repository with the dataset identifier PXD007821.

## Bioinformatics Analysis of Proteomics Data

We used several web servers to investigate the biochemical properties of the identified proteins.

The ProtFun 2.2 server produces *ab initio* predictions of protein function from sequences. The method queries a large number of other feature prediction servers to obtain information on various post-translational and localization aspects of the protein, which are integrated in final predictions of the cellular role, enzyme class (if any), and selected gene ontology (GO) categories of the submitted sequence (Jensen et al., 2002, 2003).

Propsearch server is designed to find the putative protein family if querying a new sequence using alignment methods has failed. Disregarding the order of amino acid residues in a sequence, Propsearch uses the amino acid composition instead. In addition, other properties like molecular weight, bulky residue content, small residue content, average hydrophobicity, average charge a.s.o. and the content of selected dipeptide-groups are also calculated from the sequence. A total of 144 such properties are weighted individually and used as query vectors. The weights have been trained on a set of protein families with known structures, using a genetic algorithm. Sequences in the database are also transformed into vectors, and the Euclidian distance between the query and database sequences is calculated. Distances are rank ordered, and sequences with lowest distance are reported on top (Hobohm and Sander, 1995). We used Propsearch first to better characterize proteins identified as "Putative Uncharacterized Protein" and second to improve the potential functional information on the identified proteins. In our analytical strategy, we conserved only data for proteins with Euclidian distances between 0.0 and 1.3 (Reliability 99.9%) and between 1.3 and 7.5 (Reliability 99.6%).

## Statistical Tests

Graphs and linear modeling were performed with R version 3.4.1 (R Core Team, 2013) and the tidyverse package (v 1.1.1). Datasets
and figures along with the script are available at: https://github.com/loire/amblyoma_paper.

# RESULTS

## Inhibition of ConA-Induced PBMC Proliferation by A. variegatum Saliva

To investigate the effect of *A. variegatum* on PBMC proliferation, cells were stimulated with ConA alone (positive control for PBMC activation) or with various concentrations of saliva for 72 h, and the relative PBMC activation for each saliva concentration vs. no saliva (positive control) was calculated. As shown in **Figure 1**, the relative PBMC activation decreased in a dose-dependent manner when cells were in contact with *A. variegatum* saliva. A concentration of 1.2 µg/ml of saliva reduced the ConA-induced proliferation of PBMC by 9–13%, depending on the animal from which the cells originated (**Figure 1**). The decrease in proliferation was linear depending on the concentration of the saliva, and reached 24–36% for 18.75 µg/ml of saliva (**Figure 1**). Saliva concentrations of <1.2 µg/ml was tested and did not result in statistically significant inhibition of ConA-induced PBMC proliferation, whereas levels of inhibition for higher concentrations (≥18.75 µg/ml) were similar to those for 18.75 µg/ml (data not shown). Unstimulated PBMC incubated with saliva did not show any proliferation and pilocarpine control had no effect on ConA-induced PBMC proliferation (data not shown).

## Modification of LPS Pre-stimulated Macrophage Activation by A. variegatum Saliva

We tested the effects of various concentrations of *A. variegatum* saliva on MHC II and co-stimulation molecules (CD40, CD80, and CD86) on LPS pre-stimulated macrophages. As shown in **Figure 2**, in all conditions tested, *A. variegatum* saliva had a biological effect on the expression of co-stimulation surface markers in LPS pre-stimulated macrophages, compared to unstimulated cells (no LPS). Saliva enhanced the LPS-induced down regulation of macrophages MHC II from 23 to 62%, from 39 to 72%, or from 60 to 75% depending on both the saliva concentration and the animal (overall *p* = 0.092, specific *p* < 0.05), and decreased the LPS-induced up regulation of both CD40 and CD80 in a dose-dependent manner (**Figure 2**). The saliva had no obvious effect on LPS-induced down regulation of CD86, which was already strongly decreased by LPS alone (**Figure 2**). The three independent animals gave similar profiles with individual variability in the level of expression of the co-stimulatory surface markers but, as shown in **Figure 2**, the trend curve shows that the individual variations all followed the same pattern. The effect of *A. variegatum* saliva alone on macrophages, without pre-stimulation with LPS, was less striking, with marked individual variability between animals: MHC II and CD86 were down regulated in a dose-dependent manner, CD40 and CD80 were up regulated in two out of three animals with no dose-dependence, CD86 was down regulated

in a dose-dependent manner (**Supplementary Figure 1**). Percent expression of CD40 and CD80 induced by saliva alone was 10 times lower than that with LPS pre-stimulation (**Supplementary Figure 1**). Macrophages derived from animal n° 1973 exhibited a difference MHC II levels depending on the concentration of the saliva compared to the other animals: expression increased for 31.25 µg/ml of saliva, then decreased for 62.5 µg/ml, increased slightly for 125 µg/ml and decreased slightly for 250 µg/ml. Same kind of inverted fluctuations were observed with animal n° 9906 for CD80 and animal n° 9567 for CD86 (**Supplementary Figure 1**).

As shown in **Figure 3**, LPS pre-stimulated macrophages produced NO and production was increased by adding up to 62.5 µg/ml of *A. variegatum* saliva, but decreased at higher saliva concentrations (125 and 250 µg/ml). The production of cytokines was also affected: IL-12 and TNF-α levels decreased dramatically, from a concentration of 31.25 µg/ml of saliva, except for animal n° 1973, in which a small dose-dependent decrease in TNF-α was observed (**Figure 3**). In contrast, the IL-10 level was increased by saliva in a dose-dependent manner (**Figure 3**). The three independent animals exhibited similar NO profiles and cytokine production, with individual variability in the levels reached but, like for co-stimulation surface markers, the trend curve showed that the individual variations all followed the same pattern. The effect of saliva alone on unstimulated cells resulted in very low levels of cytokine production, with increased

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Effect of *A. variegatum* saliva on ConA-induced proliferation of bovine PBMC. Bovine PBMC were stimulated with ConA alone or with different concentrations of *A. variegatum* saliva for 72 h. PBMC responses of cells in contact with ConA and saliva were compared to PBMC in contact with ConA alone (control) to calculate “relative PBMC activation.” Data are means of triplicate experiments, for each batch of saliva for the three animals tested. Blue lines and shading show respectively the linear regression fit and the 0.95 confidence interval. The *p*-value of the model is indicated on the panel.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Effect of *A. variegatum* saliva on the expression of surface markers on LPS-stimulated bovine macrophages. Bovine blood monocyte-derived macrophages were pre-stimulated for 1 h with LPS, and then for an additional 24 h with different concentrations of *A. variegatum* saliva. Cells were collected and stained for surface markers. Expression levels (mean of fluorescence intensity, MFI) of MHC II, CD40, CD80, and CD86 markers on macrophages compared with unstimulated control cells were measured by flow cytometry. Data are means of triplicate experiments, for a pool of saliva batches for the three animals tested. Blue lines and shading show the linear regression fit and the 0.95 confidence interval, respectively. The *p*-value of the model is indicated on the panel.</figcaption></figure>

levels of IL-10 up to 125 µg/ml of saliva, increased level of IL-12 up to 31.25 µg/ml then decreased IL-12 in a dose dependent manner, slightly decreased level of TNF-α in a dose-dependent manner (**Supplementary Figure 2**). Unstimulated cells did not produce NO in the presence of saliva (data not shown).

**Molecular Characterization of
A. variegatum Saliva**

Proteomics analysis of *A. variegatum* crude saliva led to the identification of 336 proteins, with no significant qualitative and quantitative variation between individuals (PRIDE PXD007821), allowing us to use pools of saliva for both cell biology assays and molecular characterization. The protein sequence data, consisting of the *Amblyomma* taxonomy (22,429 entries) and the *I. scapularis* proteome (UP000001555, 20,473 entries), were retrieved from UniProtKB (release 2015_03) for identification purposes. Among the 336 proteins identified, only 58 (17.3%) *A. variegatum* specific proteins were identified, knowing that the UniProtKB *A. variegatum* specific/particular database contained

753 entries (UniProtKB release 2015_03). Other identified proteins were divided into 236 *Amblyomma* (70.2%) and 42 *I. scapularis* (12.5%), which is the only ixodid tick species with a reference proteome in UniProtKB. It should be noted that 58 out of the 336 (17.3%) proteins identified were described as “Uncharacterized protein.”

Analysis of the nature of the identified proteins by the ProtFun 2.2 server indicated that 239 out of 336 (71.1%) of the proteins identified in *A. variegatum* saliva were enzymes. The predicted enzymes were further classified in four out of six different classes (51 Lyases, 25 Ligases, 12 Isomerases, 2 Hydrolases with no Oxidoreductase, and no Transferase), but 149 remained unclassified. Moreover, **Figure 4A** shows the classification of the identified proteins in 12 different functional categories by the ProtFun 2.2 server, the most widely represented for *A. variegatum* saliva being Cell envelope (21.4%), Translation (9.8%), Energy metabolism (9.5%), and Amino acid biosynthesis (8.3%). Nevertheless, 38.4% of the identified proteins were not classified in functional categories (**Figure 4A**). **Figure 4B** shows
<figure><img src="image_5.png" /><figcaption>**FIGURE 3** | Effect of *A. variegatum* saliva on the production of NO and cytokines by LPS-stimulated bovine macrophages. Bovine blood monocyte-derived macrophages were pre-stimulated for 1 h with LPS, and then for an additional 24 h with different concentrations of *A. variegatum* saliva. Cell culture supernatants were collected and NO production was titrated by the Griess method, IL-10, IL-12, and TNF-α production were titrated by ELISA. Data are means of triplicate experiments, for a pool of saliva batches for the three animals tested. Blue lines and shading show the linear regression fit and the 0.95 confidence interval, respectively. The *p*-value of the model is indicated on the panel.</figcaption></figure>

the prediction scores for 14 gene ontology (GO) categories, the most frequently represented for *A. variegatum* saliva being Growth factor (16.4%) and Immune response (9.8%). Like the functional categories, 55.2% of the identified proteins were not classified in the proposed GO categories.

In order to link the molecular data with the biological data obtained in the cell experiments, we focused on the characterization of *A. variegatum* saliva proteins with immunomodulatory properties already described in the literature in other tick species. We found 89 out of 336 (26.5%) with a potential role in the modulation of the host immune response. Among the immunomodulatory proteins, **Table 1** shows a classification into different groups of salivary proteins with immunomodulatory functions (linked with the observed effect of *A. variegatum* saliva on the host immune cells) as follows: immunogenic proteins including Da-p36 immunosuppressant and Antigen-5 family member (CAP superfamily of proteins), proteases inhibitors including serpins and other serine and systeine protease inhibitors, proteases

including serine proteases/serine carboxypeptidases/cathepsin-like cysteine proteases, ferritin and hemelipoproteins, calcium binding proteins, fatty acid, and histamine binding protein, nucleic acid binding proteins, anti-oxidant enzymes, heat shock proteins.

Twelve proteins were identified as ubiquitin with different protein IDs from *Amblyomma* and *Ixodes* taxa, of which 2 (F0J9P3, F0J9K1) specifically in *A. variegatum* (**Table 2**). Using the Propsearch server, we were able to characterize 6 saliva proteins that exhibited properties of ubiquitin-activating enzyme E1 ligase, with respective Euclidian distances, giving a reliability of 99.6% (**Table 2**). By contrast, we found ubiquitin-conjugating enzyme E2 ligase with Euclidian distances giving a percent reliability below the threshold fixed in our analytical approach (from 68 to 94%), and no ubiquitin-ligating enzyme E3 ligase (**Table 2**). However, we found 8 potential de-ubiquitinating enzymes, mainly ubiquitin carboxyl-terminal hydrolases, suggesting a complete and potentially functional ubiquitination system in *A. variegatum* saliva (**Table 2**). We
A

<figure><img src="image_6.png" /></figure>

B

<figure><img src="image_7.png" /></figure>

**FIGURE 4** | Graphic representation of cellular role (A) and selected gene ontology (GO) categories (B) attributed to proteins identified in the saliva of *A. variegatum* by the ProtFun 2.2 Server (http://www.cbs.dtu.dk/services/ProtFun/). The function prediction server produces *ab initio* predictions of protein function from sequences by querying a large number of other feature prediction servers to obtain information on different post-translational and localizational aspects of the protein, which are integrated into final predictions of the cellular role, enzyme class (if any), and selected gene ontology categories of the sequence submitted.

found the missing link due to the identification of an alpha spectrin (Table 2), which has the particularity of being described as a potential chimeric E2/E3 ubiquitin conjugating/ligating enzyme (Goodman et al., 2015).

Using Propsearch server, we were able to complete the characterization of *A. variegatum* saliva proteins with immunomodulatory properties, by classifying five additional proteins associated with modulation of the host immune response as follows (Table 3): a GTP-binding protein (nucleic-acid binding proteins), a superoxide dismutase and a D-dopachrome tautomerase (anti-oxidant enzymes), a 60S

ribosomal protein L3 (ribosomal proteins), and another alpha spectrin (Ubiquitin complex).

DISCUSSION

In the present study, we demonstrated high immunomodulation of bovine immune cells, lymphocytes and macrophages by *A. variegatum* saliva. This adds a new tick species to the increasing reports on immunomodulation by tick saliva (Kazimírová and Štibrániová, 2013). Increasing data from the literature on the interactions between tick saliva and the cells of the innate and adaptive immune system of vertebrate hosts tend to demonstrate that excreted-secreted salivary factors are critical both for completing the tick developmental cycle via blood feeding and for the transmission/reactivation of pathogens (Kotál et al., 2015; Šimo et al., 2017). Here, we showed that *A. variegatum* saliva induced inhibition of ConA-induced proliferation of bovine lymphocytes in a dose-dependent manner. This immunosuppression effect has already been observed in different tick models, and with various mitogens, as reviewed by Wikel and Kazimírová (Kazimírová and Štibrániová, 2013; Wikel, 2013). For example, *Ixodes ricinus* salivary glands secreted a protein that suppressed T lymphocyte proliferation and inhibited macrophage pro-inflammatory cytokines (Kovár et al., 2001, 2002). Among *Amblyomma* ticks, *A. cajennense* saliva has also been reported to inhibit lymphocyte proliferation in mouse and horse (Castagnolli et al., 2008). Another study showed that *Dermacentor andersoni* gland extract inhibited the pro-inflammatory cytokine TNF-α and up-regulated the anti-inflammatory cytokine IL-10, with down regulation of the expression of MHC II and of co-stimulation molecules (Wikel, 2013). These effects of saliva from different tick species, associated with our similar findings on *A. variegatum* saliva, strengthen the hypothesis that tick-host interactions generally result in induced lymphocyte anergy and reduced macrophage antigen presenting capacity and the classical activation pathway. In our study, we also showed that *A. variegatum* saliva emphasized the macrophage LPS-induced down regulation of CMH II, decreased the LPS-induced up regulation of CD40 and CD80. While CD40 binds on CD40-L on T cells for activation, inhibition of the expression of this molecule might affect directly the proliferative response of T-cells, as we observed for ConA-stimulated PBMC subjected to saliva. Similarly, while CD80 favors a Th1 response, the inhibition of the expression of LPS-induced CD80 might switch toward a Th2 response of lymphocytes, switching to an anti-inflammatory response. This switch has been demonstrated to be induced by *I. ricinus* salivary gland extracts (Kovár et al., 2001, 2002), and plays a role in countering the Th1-oriented response following tick aggression, as demonstrated for *Rhipicephalus microplus* in bovines (Brake and Pérez de León, 2012). This effect could result in the global inhibition of antigen presentation capacities and pro-inflammatory activation pathway of these immune cells, even if they are activated by the pathogen at the biting site. Furthermore, we evidenced that *A. variegatum* saliva also had an effect on NO production by
<table><caption>TABLE 1 | Potential immunomodulatory proteins identified in A. variegatum saliva.</caption><thead><tr><th>UniProt protein IDs (Leading protein)</th><th>Description</th><th>Organism</th><th>Number of proteins</th><th>Peptides</th><th>Unique peptides</th><th>Sequence coverage [%]</th><th>Mol. weight [kDa]</th></tr></thead><tbody><tr><td colspan="9"><strong>IMMUNOGENIC PROTEINS</strong></td></tr><tr><td>Q75PV7</td><td>Protein related to Da-p36</td><td>A. variegatum</td><td>1</td><td>6</td><td>6</td><td>29.7</td><td>27.443</td></tr><tr><td>F0J8E1</td><td>Antigen-5 family member (Fragment)</td><td>A.variegatum</td><td>1</td><td>2</td><td>2</td><td>32.7</td><td>12.308</td></tr><tr><td colspan="9"><strong>PROTEASE INHIBITORS INCLUDING</strong></td></tr><tr><td colspan="9"><strong>Serpins</strong></td></tr><tr><td>F0J9G9</td><td>Serpin (Fragment)</td><td>A. variegatum</td><td>6</td><td>11</td><td>9</td><td>59.3</td><td>23.07</td></tr><tr><td>F0J8P8</td><td>Serpin (Fragment)</td><td>A. variegatum</td><td>1</td><td>5</td><td>5</td><td>57.1</td><td>17.759</td></tr><tr><td>F0J8R0</td><td>Serpin (Fragment)</td><td>A. variegatum</td><td>1</td><td>2</td><td>1</td><td>29.1</td><td>17.208</td></tr><tr><td>A0A023GPF9</td><td>Putative tick serpins 14</td><td>A. triste</td><td>6</td><td>6</td><td>6</td><td>14.7</td><td>44.703</td></tr><tr><td>A0A023GN51</td><td>Putative tick serpins 13</td><td>A. triste</td><td>6</td><td>5</td><td>1</td><td>12.1</td><td>42.841</td></tr><tr><td>A0A023GP80</td><td>Putative tick serpins 27</td><td>A. triste</td><td>2</td><td>3</td><td>2</td><td>7.6</td><td>43.538</td></tr><tr><td>A0A023FPT3</td><td>Putative tick serpins 7 (Fragment)</td><td>A. cajennense</td><td>1</td><td>2</td><td>2</td><td>7.5</td><td>35.74</td></tr><tr><td>A0A023GGH1</td><td>Putative tick serpins 33 (Fragment)</td><td>A. triste</td><td>1</td><td>2</td><td>2</td><td>8.5</td><td>26.402</td></tr><tr><td>A0A023FUL1</td><td>Putative tick serpins 12 (Fragment)</td><td>A. cajennense</td><td>3</td><td>2</td><td>2</td><td>7.8</td><td>46.875</td></tr><tr><td>A0A023GPF5</td><td>Putative tick serpins 13</td><td>A. triste</td><td>2</td><td>2</td><td>1</td><td>9.3</td><td>42.797</td></tr><tr><td>G3MLB1</td><td>Putative uncharacterized protein (tick serpins 27)</td><td>A. maculatum</td><td>15</td><td>3</td><td>1</td><td>9.6</td><td>43.103</td></tr><tr><td colspan="9"><strong>Other serine protease inhibitors</strong></td></tr><tr><td>F0J9Y3</td><td>Kazal-type serine protease inhibitor domain protein</td><td>A. variegatum</td><td>4</td><td>3</td><td>3</td><td>20</td><td>26.71</td></tr><tr><td>A0A023GE28</td><td>Putative serine protease inhibitor I. scapularis serine protease inhibitor (Fragment)</td><td>A. triste</td><td>2</td><td>6</td><td>6</td><td>18.5</td><td>46.98</td></tr><tr><td>A0A023GE28</td><td>Putative serine protease inhibitor I. scapularis serine protease inhibitor (Fragment)</td><td>A. triste</td><td>2</td><td>6</td><td>6</td><td>18.5</td><td>46.98</td></tr><tr><td>A0A023FQ10</td><td>Putative serine protease inhibitor I. scapularis serine protease inhibitor</td><td>A. cajennense</td><td>1</td><td>2</td><td>2</td><td>5.1</td><td>44.715</td></tr><tr><td>A0A023GPF6</td><td>Putative serine proteinase inhibitor</td><td>A. triste</td><td>2</td><td>2</td><td>1</td><td>8.4</td><td>43.504</td></tr><tr><td colspan="9"><strong>Cysteine protease inhibitors</strong></td></tr><tr><td>A0A023FR27</td><td>Putative thyropin</td><td>A. cajennense</td><td>6</td><td>2</td><td>2</td><td>12.4</td><td>22.395</td></tr><tr><td>F0JA40</td><td>Putative thyropin</td><td>A. variegatum</td><td>7</td><td>3</td><td>3</td><td>12.7</td><td>22.532</td></tr><tr><td colspan="9"><strong>PROTEASES INCLUDING</strong></td></tr><tr><td colspan="9"><strong>Serine proteases</strong></td></tr><tr><td>F0J9J5</td><td>Trypsin-like serine protease (Fragment)</td><td>A. variegatum</td><td>1</td><td>6</td><td>6</td><td>41.1</td><td>21.867</td></tr><tr><td>F0J9E7</td><td>Trypsin-like serine protease (Fragment)</td><td>A. variegatum</td><td>3</td><td>6</td><td>4</td><td>48.9</td><td>19.618</td></tr><tr><td colspan="9"><strong>Serine carboxypeptidases</strong></td></tr><tr><td>A0A023GNB7</td><td>Putative serine carboxypeptidase lysosomal cathepsin a</td><td>A. triste</td><td>1</td><td>1</td><td>1</td><td>5.2</td><td>53.097</td></tr><tr><td colspan="9"><strong>Cathepsin-like cysteine proteases</strong></td></tr><tr><td>F0JA27</td><td>Cathepsin L-like cysteine proteinase B</td><td>A. variegatum</td><td>8</td><td>10</td><td>9</td><td>36.8</td><td>37.624</td></tr><tr><td>F0J903</td><td>Cysteine proteinase cathepsin L (Fragment)</td><td>A. variegatum</td><td>2</td><td>5</td><td>5</td><td>38.6</td><td>23.017</td></tr><tr><td>G3MHB0</td><td>Putative uncharacterized protein (Fragment) (peptidase C1 family; cathepsin L-like cysteine proteinase b)</td><td>A. maculatum</td><td>3</td><td>2</td><td>1</td><td>11.5</td><td>38.867</td></tr><tr><td colspan="9"><strong>FERRITIN AND HEMELIPOPROTEINS</strong></td></tr><tr><td>G3MPJ3</td><td>Ferritin</td><td>A. maculatum</td><td>2</td><td>3</td><td>3</td><td>19.5</td><td>22.296</td></tr><tr><td>F0J8W4</td><td>Heme lipoprotein (Fragment)</td><td>A. variegatum</td><td>2</td><td>33</td><td>22</td><td>63.5</td><td>50.007</td></tr><tr><td>F0J9N1</td><td>Heme lipoprotein (Fragment)</td><td>A. variegatum</td><td>1</td><td>20</td><td>18</td><td>67.4</td><td>27.005</td></tr><tr><td>F0J9N8</td><td>Hemelipoglycoprotein (Fragment)</td><td>A. variegatum</td><td>1</td><td>14</td><td>14</td><td>48.2</td><td>24.787</td></tr><tr><td>F0J9G0</td><td>Heme lipoprotein (Fragment)</td><td>A. variegatum</td><td>1</td><td>11</td><td>8</td><td>95.4</td><td>14.828</td></tr><tr><td>A0MVX0</td><td>Heme lipoprotein</td><td>A. americanum</td><td>1</td><td>31</td><td>5</td><td>16.3</td><td>177.04</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 1 | Continued</caption><thead><tr><th>UniProt protein IDs (Leading protein)</th><th>Description</th><th>Organism</th><th>Number of proteins</th><th>Peptides</th><th>Unique peptides</th><th>Sequence coverage [%]</th><th>Mol. weight [kDa]</th></tr></thead><tbody><tr><td>B7Q406</td><td>Hemelipoglycoprotein. putative</td><td>I. scapularis</td><td>2</td><td>3</td><td>1</td><td>1.7</td><td>177.65</td></tr><tr><td>B7QIL0</td><td>Hemelipoglycoprotein. putative (Fragment)</td><td>I. scapularis</td><td>1</td><td>2</td><td>1</td><td>1.3</td><td>172.44</td></tr><tr><td colspan="8"><strong>CALCIUM BINDING PROTEINS</strong></td></tr><tr><td>F0JA25</td><td>Calreticulin</td><td>A. variegatum</td><td>1</td><td>4</td><td>1</td><td>16.3</td><td>47.318</td></tr><tr><td>Q64KA0</td><td>Calreticulin</td><td>A. brasiliense</td><td>17</td><td>4</td><td>1</td><td>16.3</td><td>47.241</td></tr><tr><td>G3MK71</td><td>Calponin</td><td>A. maculatum</td><td>6</td><td>3</td><td>3</td><td>24.7</td><td>17.614</td></tr><tr><td colspan="8"><strong>FATTY ACID AND HISTAMINE BINDING PROTEIN</strong></td></tr><tr><td>A0A023G6C8</td><td>Putative lipocal-1 1</td><td>A. triste</td><td>1</td><td>1</td><td>1</td><td>8.5</td><td>24.71</td></tr><tr><td>A0A023GBM1</td><td>Putative phospholipase a2 I. scapularis phospholipase a2</td><td>A. triste</td><td>2</td><td>4</td><td>4</td><td>11.4</td><td>44.79</td></tr><tr><td colspan="8"><strong>NUCLEIC ACID BINDING PROTEINS</strong></td></tr><tr><td>G3MSL8</td><td>Putative uncharacterized protein (Putative gtp-binding protein)</td><td>A. maculatum</td><td>10</td><td>2</td><td>2</td><td>12.4</td><td>24.689</td></tr><tr><td>G3MRA1</td><td>Putative uncharacterized protein (gtp-binding protein)</td><td>A. maculatum</td><td>4</td><td>1</td><td>1</td><td>6</td><td>22.968</td></tr><tr><td>A0A023FGS1</td><td>Putative zinc-binding protein of the histidine triad hit family</td><td>A. cajennense</td><td>3</td><td>1</td><td>1</td><td>11.1</td><td>13.53</td></tr><tr><td colspan="8"><strong>ANTI-OXYDANT ENZYMES</strong></td></tr><tr><td>A0A023FIC2</td><td>Putative glutathione s-transferase mu class Rhipicephalus annulatus glutathione s-transferase</td><td>A. cajennense</td><td>7</td><td>6</td><td>6</td><td>29.1</td><td>25.713</td></tr><tr><td>A0A023GJMJ</td><td>Putative glutathione s-transferase</td><td>A. triste</td><td>2</td><td>2</td><td>2</td><td>14.7</td><td>24.394</td></tr><tr><td>A0A023FXR2</td><td>Putative glutathione s-transferase I. scapularis glutathione s-transferase</td><td>A. parvum</td><td>5</td><td>2</td><td>2</td><td>10.4</td><td>26.241</td></tr><tr><td>A0A023FFS7</td><td>Putative glutathione s-transferase I. scapularis glutathione s-transferase</td><td>A. cajennense</td><td>1</td><td>1</td><td>1</td><td>5.5</td><td>24.303</td></tr><tr><td>A0A023FQ81</td><td>Putative glutathione s-transferase I. scapularis glutathione s-transferase (Fragment)</td><td>A. cajennense</td><td>4</td><td>1</td><td>1</td><td>8.9</td><td>23.244</td></tr><tr><td>A0A023G5I8</td><td>Putative glutathione s-transferase I. scapularis glutathione s-transferase (Fragment)</td><td>A. triste</td><td>4</td><td>1</td><td>1</td><td>8.6</td><td>27.412</td></tr><tr><td>B7PDV0</td><td>Glutathione S-transferase. putative</td><td>I. scapularis</td><td>4</td><td>1</td><td>1</td><td>12.2</td><td>8.6628</td></tr><tr><td>G3MMK0</td><td>Putative uncharacterized protein (glutathione s-transferase I. scapularis glutathione s-transferase)</td><td>A. maculatum</td><td>2</td><td>1</td><td>1</td><td>5.9</td><td>26.159</td></tr><tr><td>A0SHR2</td><td>Protein disulfide-isomerase</td><td>A. variegatum</td><td>8</td><td>3</td><td>3</td><td>7.2</td><td>54.923</td></tr><tr><td>A0A023FZW2</td><td>Protein disulfide-isomerase</td><td>A. parvum</td><td>5</td><td>4</td><td>1</td><td>9.2</td><td>57.252</td></tr><tr><td>A0A023FNM6</td><td>Putative alkyl hydroperoxide reductase thiol specific antioxidant (Fragment)</td><td>A. cajennense</td><td>2</td><td>2</td><td>2</td><td>10.7</td><td>23.542</td></tr><tr><td>A0A023FHT1</td><td>Putative thioredoxin peroxidase</td><td>A. cajennense</td><td>14</td><td>4</td><td>4</td><td>24</td><td>20.494</td></tr><tr><td>A0A023FDL5</td><td>Putative selenoprotein w 2a</td><td>A. cajennense</td><td>1</td><td>1</td><td>1</td><td>22.1</td><td>8.6959</td></tr><tr><td>F0J9H3</td><td>Peroxidasin (Fragment)</td><td>A. variegatum</td><td>1</td><td>1</td><td>1</td><td>8.1</td><td>27.07</td></tr><tr><td>A0A023FTX2</td><td>Superoxide dismutase [Cu-Zn]</td><td>A. parvum</td><td>1</td><td>1</td><td>1</td><td>9.7</td><td>15.744</td></tr><tr><td colspan="8"><strong>HEAT SHOCK PROTEINS</strong></td></tr><tr><td>A0A023GNY6</td><td>Putative heat shock protein (Fragment)</td><td>A. triste</td><td>11</td><td>16</td><td>13</td><td>35.6</td><td>73.124</td></tr><tr><td>A0A023GP15</td><td>Putative heat shock protein 90</td><td>A. triste</td><td>9</td><td>7</td><td>2</td><td>13.5</td><td>84.001</td></tr><tr><td>A0A023GLJ3</td><td>Putative heat shock protein 70</td><td>A. triste</td><td>7</td><td>5</td><td>2</td><td>10.5</td><td>73.014</td></tr><tr><td>A0A023FRH0</td><td>Putative heat shock protein</td><td>A. cajennense</td><td>1</td><td>1</td><td>1</td><td>5.7</td><td>21.73</td></tr><tr><td>G3MF27</td><td>Putative uncharacterized protein (heat shock 70 kDa protein 5)</td><td>A. maculatum</td><td>7</td><td>4</td><td>3</td><td>8</td><td>69.078</td></tr><tr><td>G3MF42</td><td>Putative uncharacterized protein (heat shock protein)</td><td>A. maculatum</td><td>2</td><td>1</td><td>1</td><td>7.8</td><td>21.722</td></tr></tbody></table>
<table><caption>TABLE 2 | Potential ubiquitination complex identified in A. variegatum saliva.</caption><thead><tr><th>UniProt protein IDs (Leading protein)</th><th>Description</th><th>Organism</th><th>Number of proteins</th><th>Peptides</th><th>Unique peptides</th><th>Sequence coverage [%]</th><th>Mol. weight [kDa]</th><th>Propsearch distance (reliability)</th><th>Propsearch putative protein family</th></tr></thead><tbody><tr><td colspan="10"><strong>UBIQUITIN</strong></td></tr><tr><td>F0J9K1</td><td>Ubiquitin (Fragment)</td><td>A. variegatum</td><td>12</td><td>2</td><td>2</td><td>29.8</td><td>9.4959</td><td>–</td><td>–</td></tr><tr><td colspan="10"><strong>E1 UBIQUITIN LIGASE</strong></td></tr><tr><td>A0A023GDS5</td><td>Putative microtubule associated complex (Fragment)</td><td>A. triste</td><td>2</td><td>12</td><td>9</td><td>14.1</td><td>117.06</td><td>5.38 (99.6%)</td><td>Ubiquitin-activating enzyme E1 (A1S9 protein).</td></tr><tr><td>A0A023G500</td><td>Putative nidogen (Fragment)</td><td>A. triste</td><td>1</td><td>1</td><td>1</td><td>1.2</td><td>124.2</td><td>6.21 (99.6%)</td><td>Ubiquitin-activating enzyme E1 homolog (D8)</td></tr><tr><td>B7QMC8</td><td>Alpha-macroglobulin. putative (Fragment)</td><td>I. scapularis</td><td>1</td><td>8</td><td>1</td><td>7</td><td>121.21</td><td>4.83–5.37 (99.6%)</td><td>Ubiquitin-activating enzyme E1</td></tr><tr><td>B7PN44</td><td>Netrin-4, putative</td><td>I. scapularis</td><td>1</td><td>2</td><td>2</td><td>5.1</td><td>41.593</td><td>7.21 (99.6%)</td><td>Ubiquitin-activating enzyme E1 X (Fragment)</td></tr><tr><td>B7PQH0</td><td>Hemolectin. putative (Fragment)</td><td>I. scapularis</td><td>1</td><td>2</td><td>2</td><td>2.4</td><td>115.88</td><td>5.66–5.67 (99.6%)</td><td>Ubiquitin-activating enzyme E1</td></tr><tr><td>A0A023GPI0</td><td>Uncharacterized protein (Fragment)</td><td>A. triste</td><td>3</td><td>28</td><td>1</td><td>32.3</td><td>131.52</td><td>4.90–5.42 (99.6%)</td><td>Ubiquitin-activating enzyme E1</td></tr><tr><td colspan="10"><strong>E2 UBIQUITIN LIGASE</strong></td></tr><tr><td>B5M727</td><td>Alpha-2-macroglobulin (Fragment)</td><td>A. americanum</td><td>2</td><td>5</td><td>2</td><td>42.6</td><td>19.026</td><td>8.69–10.05 (94%) &gt; x &gt; 80%</td><td>Ubiquitin-conjugating enzyme E2-24 kDa (EC 6.3.2.19)</td></tr><tr><td>A0A023G356</td><td>Putative secreted protein (Fragment)</td><td>A. triste</td><td>1</td><td>3</td><td>3</td><td>12.2</td><td>21.619</td><td>11.17 (80%)</td><td>Ubiquitin-conjugating enzyme E2-24 kDa (EC 6.3.2.19)</td></tr><tr><td>A0A023FDX3</td><td>Putative I. scapularis glutathione peroxidase (Fragment)</td><td>A. cajennense</td><td>1</td><td>1</td><td>1</td><td>12.8</td><td>10.918</td><td>11.13–12.43 (80%) &gt; x &gt; 68%</td><td>Ubiquitin-conjugating enzyme variant MMS2 homolog (UEV MMS2)</td></tr><tr><td>G3MNL6</td><td>Putative uncharacterized protein (alpha crystallins?)</td><td>A. maculatum</td><td>2</td><td>1</td><td>1</td><td>13.1</td><td>19.912</td><td>9.69 (87%)</td><td>Ubiquitin-conjugating enzyme E2-24 kDa (EC 6.3.2.19)</td></tr><tr><td colspan="10"><strong>E2/E3 UBIQUITIN LIGASE</strong></td></tr><tr><td>B7P1U8</td><td>Spectrin alpha chain. putative</td><td>I. scapularis</td><td>4</td><td>4</td><td>4</td><td>2</td><td>272.3</td><td>–</td><td>–</td></tr><tr><td colspan="10"><strong>DE-UBIQUITINATING ENZYME</strong></td></tr><tr><td>A0A023GME3</td><td>Putative vitellogenin-1 (Fragment)</td><td>A. triste</td><td>1</td><td>38</td><td>10</td><td>28.4</td><td>153.63</td><td>6.12–7.44 (99.6%)</td><td>Ubiquitin carboxyl-terminal hydrolase 15 (EC 3.1.2.15)</td></tr></tbody></table>

(Continued)
<table><caption>TABLE 2 | Continued</caption><thead><tr><th>UniProt protein IDs (Leading protein)</th><th>Description</th><th>Organism</th><th>Number of proteins</th><th>Peptides</th><th>Unique peptides</th><th>Sequence coverage [%]</th><th>Mol. weight [kDa]</th><th>Propsearch distance (reliability)</th><th>Propsearch putative protein family</th></tr></thead><tbody><tr><td>A0A023GMC7</td><td>Putative vitellogenin-1 (Fragment)</td><td>A. triste</td><td>2</td><td>34</td><td>6</td><td>30.6</td><td>123.4</td><td>6.42-6.80 (99.6%)</td><td>Ubiquitin carboxyl-terminal hydrolase 15 (EC 3.1.2.15)</td></tr><tr><td>A0A023GP74</td><td>Putative titin (Fragment)</td><td>A. triste</td><td>2</td><td>21</td><td>15</td><td>10.4</td><td>389.34</td><td>7.13 (99.6%)</td><td>Probable ubiquitin carboxyl-terminal hydrolase FAF (EC 3.1.2.15)</td></tr><tr><td>A0A023G0F9</td><td>Putative peroxinectin I. scapularis peroxinectin (Fragment)</td><td>A. parvum</td><td>1</td><td>16</td><td>3</td><td>32.1</td><td>87.425</td><td>5.42 (99.6%)</td><td>Ubiquitin carboxyl-terminal hydrolase 20 (EC 3.1.2.15)</td></tr><tr><td>A0A023GCH5</td><td>Putative neurixin iii-alpha (Fragment)</td><td>A. triste</td><td>1</td><td>9</td><td>9</td><td>8.4</td><td>158.55</td><td>5.72 (99.6%)</td><td>Ubiquitin carboxyl-terminal hydrolase 19 (EC 3.1.2.15)</td></tr><tr><td>A0A023GLZ9</td><td>Putative extracellular matrix glycoprotein laminin subunit beta (Fragment)</td><td>A. triste</td><td>1</td><td>9</td><td>4</td><td>15.7</td><td>92.371</td><td>5.25 (99.6%)</td><td>Ubiquitin carboxyl-terminal hydrolase 20 (EC 3.1.2.15)</td></tr><tr><td>B7P403</td><td>Type II transmembrane protein, putative</td><td>I. scapularis</td><td>1</td><td>6</td><td>6</td><td>2.6</td><td>319.6</td><td>5.79 (99.6%)</td><td>Probable ubiquitin carboxyl-terminal hydrolase FAF (EC 3.1.2.15)</td></tr><tr><td>B7PUM7</td><td>Peroxinectin, putative</td><td>I. scapularis</td><td>1</td><td>4</td><td>1</td><td>6.7</td><td>90.84</td><td>5.18 (99.6%)</td><td>Ubiquitin carboxyl-terminal hydrolase 20 (EC 3.1.2.15)</td></tr></tbody></table>
<table><caption>TABLE 3 | Propsearch-defined potential immunomodulatory proteins identified in A. variegatum saliva.</caption><thead><tr><th>UniProt protein IDs<br>(Leading protein)</th><th>Description</th><th>Organism</th><th>Number of<br>proteins</th><th>Peptides</th><th>Unique<br>peptides</th><th>Sequence<br>coverage [%]</th><th>Mol. weight<br>[kDa]</th><th>Propsearch distance<br>(reliability)</th><th>Propsearch putative<br>protein family</th></tr></thead><tbody><tr><td colspan="10"><strong>NUCLEIC-ACID BINDING PROTEINS</strong></td></tr><tr><td>F0J9Z3</td><td>RAS oncogene family member RAB1</td><td>A. variegatum</td><td>5</td><td>1</td><td>1</td><td>11</td><td>16.156</td><td>5.65-7.34 (99.6%)</td><td>GTP-binding protein</td></tr><tr><td colspan="10"><strong>ANTI-OXYDANT ENZYMES</strong></td></tr><tr><td>A0A023G0Q2</td><td>Putative his-rich 1</td><td>A. parvum</td><td>2</td><td>2</td><td>2</td><td>14.9</td><td>15.48</td><td>1.08-2.02 (99.9 &gt; x &gt; 99.6%)</td><td>Superoxide dismutase<br>[Mn-Fe] (EC 1.15.1.1)<br>(Fragment)</td></tr><tr><td>G3MKP9</td><td>Carboxypeptidase</td><td>A. maculatum</td><td>2</td><td>1</td><td>1</td><td>2.6</td><td>52.321</td><td>5.60 (99.6%)</td><td>Dopachrome<br>tautomerase precursor<br>(EC 5.3.3.12) (DT)</td></tr><tr><td colspan="10"><strong>RIBOSOMAL PROTEINS</strong></td></tr><tr><td>G3MMA6</td><td>Putative uncharacterized protein<br>(Putative alpha crystallins?)</td><td>A. maculatum</td><td>8</td><td>1</td><td>1</td><td>6</td><td>28.586</td><td>3.77-7.18 (99.6%)</td><td>60S ribosomal protein<br>L3</td></tr><tr><td colspan="10"><strong>UBIQUITIN COMPLEX</strong></td></tr><tr><td>A0A023GMF4</td><td>Putative beta-spectrin (Fragment)</td><td>A. triste</td><td>2</td><td>1</td><td>1</td><td>0.8</td><td>205.72</td><td>5.14-5.44 (99.6%)</td><td>Spectrin alpha chain</td></tr></tbody></table>

macrophages through complex modulation of this microbicidal and immunosuppressive molecule (Wink et al., 2011): it increased its production at relatively low concentrations but decreased its production at high concentrations. Macrophages use NO together with several other free radicals as a mechanism for killing pathogens and inducing an inflammatory response, but high levels of NO has been shown to have suppressive effects on lymphocyte proliferation (MacMicking et al., 1997). We observed that NO production was increased by small quantities of saliva (31.25-62.5 µg/ml), then strongly decreased by high concentrations of saliva (up to 250 µg/ml). This inhibition correlated with the decrease in TNF-α and the increase in IL-10 we observed when the global *in vitro* immune response became clearly anti-inflammatory. It is possible that high production of NO plays a role in the establishment of lymphocyte anergy, before the deactivation of macrophages by *A. variegatum* saliva avoids tissue damage in the host that would be deleterious to tick feeding. In a previous study on *A. variegatum* saliva performed in the *in vivo* mouse model, we showed that, in skin lesions mimicking tick biting, *A. variegatum* saliva impairs leukocyte infiltration (Vachieri et al., 2015). Experiments in an *in vitro* human model further demonstrated that tick saliva impaired the migration of antigen presenting cells to the draining lymph nodes and affected the recruitment of blood monocytes to the skin lesion and prevented their maturation *in situ* (Vachieri et al., 2015). These results strengthen our findings with bovine cells, and show how the saliva of *A. variegatum*, by modulating chemotaxis, the proliferative capacity of lymphocytes, the regulatory capacities of macrophages, diverts the host's immune system to its own advantage, and the consequences this may have for the transmission or reactivation of pathogens.

Some of the strong and specific effects of tick saliva in modulating cytokine production by macrophage and lymphocytes, and hence the impacts on immune response, can be analyzed at molecular scale to characterize saliva molecules thanks to the advent of proteomics. Among the proteins we evidenced in *A. variegatum* female saliva and shown to have direct immunomodulatory functions, the CAP superfamily of proteins (comprising the CRISP, Antigen-5, and pathogen-related-1 families) has been found in most tick sialomes studied to date, with molecular characteristics similar to wasp-venom proteins and annotated as Antigen-5 (Gibbs et al., 2008). Additionally, recombinant Da-p36 has been demonstrated to suppress T-lymphocyte-mitogen-driven *in vitro* proliferation of splenocytes (Alarcon-Chaidez et al., 2003), and could therefore be part of the molecular mechanism behind our observation of inhibition lymphocyte proliferation in the presence of *A. variegatum* saliva.

The most representative category of proteins with potential immunomodulatory properties we identified were protease inhibitors, mainly serpins. Serpins are abundant in ticks, and one of their functions is to modulate the host immune system (Chmelar et al., 2017). The first tick serpin described for its effect on host defense mechanisms was named *I. ricinus* immunosuppressor (Iris) and was able to inhibit T-cell and splenocyte proliferation and to modify cytokine levels derived from PBMC (Leboulle et al., 2002), as well as to bind
to monocytes/macrophages and to suppress TNF-α secretion (Prevot et al., 2009). These results are very similar to our observations, and we hypothesize the existence of an Iris-like protein among the serpins identified in *A. variegatum* saliva. Other protease inhibitors we identified are cysteine protease inhibitors named thyropins. Recently, an exogenous thyropin was shown to affect IL-12 secretion of dendritic cells (DCs) (Zavašnik-Bergant and Bergant Marušič, 2016), which leads us think that the thyropins in *A. variegatum* saliva could be involved in the decreased production of IL-12 by bovine macrophages, like DCs. In contrast to protease inhibitors, we identified proteases in the saliva of *A. variegatum*, whose carboxypeptidases could play a role in the destruction of inflammatory peptide agonists (Ribeiro et al., 2011) and cathepsin-like cysteine proteases could destroy cellular mediators of inflammation (Radulović et al., 2014).

During tick feeding, the processing of the blood meal is accompanied by a high risk of oxidative stress linked to large amounts of heme, a bi-product of hemoglobin digestion, and free iron in the host's blood (Graca-Souza et al., 2006). We identified ferritin and hemelipoproteins in *A. variegatum* saliva, which could be used to dump heme and free iron in the host to avoid oxidative stress, and at the same time benefit tick-borne pathogens that may need iron to proliferate (Rouault, 2004). We also found numerous anti-oxidant enzymes (including glutathione-S transferase (GST), protein disulfide isomerase, thioredoxin peroxidase, peroxidasin, superoxide dismutase, and selenoproteins) in the saliva of *A. variegatum*, which could play a role in reducing the toxicity of reactive oxygen and nitrogen species (such as NO) produced as part of the wound healing mechanism and anti-microbial defenses (Vider et al.,

2001; Rojkind et al., 2002). Preventing damage caused by oxidative stress could benefit both the host and the tick, but given the susceptibility of pathogens to the products of oxidative stress, anti-oxidants in tick saliva could facilitate the transmission of tick-borne pathogens. Interestingly, we found D-dopachrome tautomerase, the functional homolog of macrophage migration inhibitory factor (MIF) (Merk et al., 2011), which is a pro-inflammatory cytokine produced by ticks to regulate host responses to tick feeding (Wasala et al., 2012). D-dopachrome tautomerase/MIF could be involved in the immunological silence induced by *A. variegatum*, with no chemotactic recruitment of host macrophages. Conversely, we evidenced several proteins dedicated to anti-inflammatory response. For instance, extracellular nucleic acids are potent proinflammatory molecules, and nucleic acid binding proteins could be part of the system used by the tick to modulate host inflammation (Radulović et al., 2014). In the same way, ribosomal and heat shock proteins have been described as potent anti-inflammatory molecules (Pockley, 2003; Poddar et al., 2013).

Fatty acid and histamine binding proteins also play roles in preventing inflammation during tick feeding. Histamine is a potent pro-inflammatory molecule that is released by cellular mediators of inflammation such as mast cells and neutrophils (White and Kaliner, 1987). Like *A. americanum*, *A. variegatum* saliva contains histamine-binding protein/lipocalin to sequester histamine and stop the corresponding inflammatory response (Radulović et al., 2014). Although other authors have demonstrated secretion of prostaglandin E2 (PGE2) and prostacyclin (PGI) in tick saliva (Ribeiro et al., 1988; Bowman et al., 1995), we did not identify these proteins among our

<figure><img src="image_8.png" /></figure>
proteomics data. However, we did identify a phospholipase A2, which could stimulate prostacyclin production (Bowman et al., 1997) and induce PGE2 production by host cells (Hara et al., 1991; Miyake et al., 1994). Tick saliva PGE2 has been shown to modulate the function of macrophages, particularly cytokine profiles (Poole et al., 2013), as we also evidenced in our study. Indirect PGI and PGE2 regulation could have been developed by *A. variegatum* to manipulate host macrophages, which, as a central immune cell, could pervert the whole immune system of the host.

Multifunctional calcium binding proteins, such as calponin and calreticulin (CRT), are present in the saliva of *A. variegatum* but apart from depleting $Ca^{2+}$ to prevent activation of blood clotting (Astrup, 2009), their function is not yet known. Nevertheless, an interesting hypothesis has emerged from *A. americanum* CRT, which could skew a Th-2 host response to tick feeding (Kim et al., 2015), and which is consistent with our observation of a M2-induced profile of macrophage cytokines/metabolism by *A. variegatum* saliva.

Finally, we also discovered an original ubiquitination complex in the saliva of *A. variegatum*, both on its molecular architecture and on its extracellular localization. In eukaryote cells, ubiquitin is transferred by three enzymes: the ubiquitin-activating enzyme E1 ligase catches ubiquitin and transfers it to the ubiquitin-conjugating enzyme E2 ligase, which associates with ubiquitin-ligating enzyme E3 ligase to transfer ubiquitin to the target substrate (Tanaka et al., 1998). Intriguingly, the ubiquitination we evidenced in *A. variegatum* is simplified with spectrin playing the role of a chemic E2/E3 ubiquitin-conjugating-ligating enzyme (Goodman et al., 2015). Whereas, ubiquitination has long been well described in connection with proteasome degradation of proteins, in the last decade, an increasing literature has proposed an important role for ubiquitination in immunity (Ben-Neriah, 2002; Zinngrebe et al., 2014), particularly in regulating both innate and adaptive immune responses (Malynn and Ma, 2010; Park et al., 2014; Hu and Sun, 2016). Interestingly, extracellular ubiquitin has been implicated in lymphocyte differentiation, suppression of the immune response, and prevention of inflammation (Sujashvili, 2016), particularly by binding to the CXCR4 receptor (Majetschak, 2011; Scofield et al., 2015). Moreover, extracellular ubiquitin can bind directly to regulatory T cells; enhancing their inhibitory effect on effector T cells proliferation (Cao et al., 2014). Exogenous ubiquitin can also reduce TNF-$\alpha$ production (Majetschak et al., 2004). The last two elements are in complete agreement with our observations on inhibition of lymphocyte proliferation and reduction of TNF-$\alpha$ production by macrophages in the presence of tick saliva, and make us suggest that *A. variegatum* uses its extracellular ubiquitination complex as an effective molecular weapon for render the immune responses of the host ineffective during the blood meal.

To conclude, our data shed new light on the molecular determinants of *A. variegatum* saliva-induced immunomodulation of bovine immunity cells and allows us to propose a new synthesis scheme of the cellular and molecular events that take place at the tick-host interface (Figure 5). Further studies on the effect of *A. variegatum* saliva on bovine

PBMC and macrophages are needed to better understand the mechanisms underlying our observations, particularly the comparative effect of saliva on different lymphocyte sub-populations such as cytotoxic T cells or NK cells, but also qualitative or quantitative cytokine modulation such as INF-$\gamma$, IL-2, IL-4, or IL-17 with impacts on Th1/Th2 or Th17 balance, as shown in *I. ricinus* (Kopecký et al., 1999; Kovár et al., 2001; Leboulle et al., 2002) and in mosquitoes (Wanasen et al., 2004). Effects on macrophages presenting capacities also need to be elucidated in connection with pathogen uptake and presentation in the presence of tick saliva. These cellular investigations will need to be combined with further molecular characterization of saliva proteins to identify the key molecules in immunomodulation mechanisms. Overall, our data provide new insights into tick-ruminant interactions for *A. variegatum*, including their involvement in vectorial competence for pathogen transmission (*E. ruminantium* infection) or on pathogen reactivation (dermatophilosis), and pave the way for integrative strategies to interfere with both the immunosuppressive and infectious processes in the corresponding tick-borne diseases.

## AUTHOR CONTRIBUTIONS

VR, BF, FS, and PH: designed the work; VR, AV, LC, RA, KG-G, and PH: generated the biological samples; VR, BF, AV, LC, AA, RA, KG-G, OV, ED, MS, and PH: performed the experiments; VR, BF, AV, LC, AA, OV, ED, MS, FS, EL, NV, and PH: acquired, analyzed, and interpreted the data; VR and PH: wrote the manuscript; VR, BF, AV, LC, AA, OV, ED, MS, RA, KG-G, FS, EL, NV, and PH: revised and approved the final version of the manuscript.

## ACKNOWLEDGMENTS

This research was supported by CIRAD and INRA and a collaborative project on *A. variegatum* with USDA-Agricultural Research Service.

This study was partly conducted in the framework of the project MALIN "Surveillance, diagnosis, control and impact of infectious diseases of humans, animals and plants in tropical islands" supported by the European Regional Development Fund (ERDF) and the Regional Council of Guadeloupe.

Mass spectrometry experiments were carried out using facilities of the Functional Proteomics Platform of Montpellier.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2017.00521/full#supplementary-material

**Supplementary Figure 1** | Effect of *A. variegatum* saliva on the expression of surface molecules of bovine unstimulated macrophages. Bovine blood monocyte-derived macrophages were cultivated for 24 h with different concentrations of *A. variegatum* saliva. Cells were collected and stained for MHC II, CD40, CD80, and CD86 surface markers. Expression levels (mean of fluorescence intensity, MFI) of markers on macrophages compared to those on
control unstimulated cells were measured by flow cytometry. Blue lines and shading show the linear regression fit and the 0.95 confidence interval, respectively. The *p*-value of the model is indicated on the panel.

**Supplementary Figure 2** | Effect of *A. variegatum* saliva on the production of cytokines by bovine unstimulated macrophages. Bovine blood monocyte-derived

macrophages were stimulated for 24 h with different concentrations of *A. variegatum* saliva. Cell culture supernatants were collected and IL-10, IL-12, and TNF-α production was titrated by ELISA. (Data are results obtained for one experiment performed on three different animals; note different scales on graphs). Blue lines and shading show the linear regression fit and the 0.95 confidence interval, respectively. The *p*-value of the model is indicated on the panel.

REFERENCES

Alarcon-Chaidez, F. J., Müller-Doblies, U. U., and Wikel, S. (2003). Characterization of a recombinant immunomodulatory protein from the salivary glands of *Dermacentor andersoni*. *Parasite Immunol.* 25, 69–77. doi: 10.1046/j.1365-3024.2003.00609.x

Allsopp, B. A. (2010). Natural history of *Ehrlichia ruminantium*. *Vet. Parasitol.* **167**, 123–135. doi: 10.1016/j.vetpar.2009.09.014

Araujo, R. N., Franco, P. F., Rodrigues, H., Santos, L. C., McKay, C. S., Sanhueza, C. A., et al. (2016). Amblyomma sculptum tick saliva: α-Gal identification, antibody response and possible association with red meat allergy in Brazil. *Int. J. Parasitol.* **46**, 213–220. doi: 10.1016/j.ijpara.2015.12.005

Astrup, T. (2009). Blood clotting and related processes. *Adv. Enzymol. Relat. Areas Mol. Biol.* **10**:1. doi: 10.1002/9780407122556.ch1

Ben-Neriah, Y. (2002). Regulatory functions of ubiquitination in the immune system. *Nat. Immunol.* **3**, 20–26. doi: 10.1038/ni0102-20

Bournez, L., Cangi, N., Lancelot, R., Pleydell, D. R., Stachurski, F., Bouyer, J., et al. (2015). Parapatric distribution and sexual competition between two tick species, *Amblyomma variegatum* and *A. hebraeum* (Acari, Ixodidae), in Mozambique. *Parasit. Vectors* **8**:504. doi: 10.1186/s13071-015-1116-7

Bowman, A. S., Gengler, C. L., Surdick, M. R., Zhu, K., Essenberg, R. C., Sauer, J. R., et al. (1997). A novel phospholipase A2 activity in saliva of the lone star tick, *Amblyomma americanum* (L.). *Exp. Parasitol.* **87**, 121–132. doi: 10.1006/expr.1997.4201

Bowman, A. S., Sauer, J. R., Zhu, K., and Dillwith, J. W. (1995). Biosynthesis of salivary prostaglandins in the lone star tick, *Amblyomma americanum*. *Insect Biochem. Mol. Biol.* **25**, 735–741. doi: 10.1016/0965-1748(95)00013-L

Brake, D. K., and Pérez de León, A. A. (2012). Immunoregulation of bovine macrophages by factors in the salivary glands of *Rhipicephalus microplus*. *Parasit. Vectors* **5**:38. doi: 10.1186/1756-3305-5-38

Brossard, M., and Wikel, S. K. (2004). Tick immunobiology. *Parasitology* **129**(Suppl.), S161-S176. doi: 10.1017/S0031182004004834

Cao, Y., Li, C., Zhang, Q., Wang, Y., and Xia, R. (2014). Extracellular ubiquitin enhances the suppressive effects of regulatory T cells on effector T cell responses. *Clin. Lab.* **60**, 1983–1991. doi: 10.7754/Clin.Lab.2014.140314

Castagnolli, K. C., Ferreira, B. R., Franzin, A. M., de Castro, M. B., and Szabó, M. P. (2008). Effect of *Amblyomma cajennense* ticks on the immune response of BALB/c mice and horses. *Ann. N.Y. Acad. Sci.* **1149**, 230–234. doi: 10.1196/annals.1428.028

Chmelar, J., Kotál, J., Karim, S., Kopacek, P., Francischetti, I. M., Pedra, J. H., et al. (2016a). Sialomes and mialomes: a systems-biology view of tick tissues and tick-host interactions. *Trends Parasitol.* **32**, 242–254. doi: 10.1016/j.pt.2015.10.002

Chmelar, J., Kotál, J., Kopecký, J., Pedra, J. H., and Kotsyfakis, M. (2016b). All for one and one for all on the tick-host battlefield. *Trends Parasitol.* **32**, 368–377. doi: 10.1016/j.pt.2016.01.004

Chmelar, J., Kotál, J., Langhansová, H., Kotsyfakis, M. (2017). Protease inhibitors in tick saliva: the role of serpins and cystatins in tick-host-pathogen interaction. *Front. Cell. Infect. Microbiol.* **7**:216. doi: 10.3389/fcimb.2017.00216

Cox, J., and Mann, M. (2008). MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. *Nat. Biotechnol.* **26**, 1367–1372. doi: 10.1038/nbt.1511

Cox, J., Neuhauser, N., Michalski, A., Scheltema, R. A., Olsen, J. V., and Mann, M. (2011). Andromeda: a peptide search engine integrated into the MaxQuant environment. *J. Proteome Res.* **10**, 1794–1805. doi: 10.1021/pr101065j

de la Fuente, J., Villar, M., Cabezas-Cruz, A., Estrada-Peña, A., Ayllón, N., and Alberdi, P. (2016). Tick-host-pathogen interactions: conflict and cooperation. *PLoS Pathog.* **12:e1005488. doi: 10.1371/journal.ppat.1005488**

Francischetti, I. M., Sa-Nunes, A., Mans, B. J., Santos, I. M., and Ribeiro, J. M. (2009). The role of saliva in tick feeding. *Front. Biosci.* **14**, 2051–2088. doi: 10.2741/3363

Garcia, G. R., Gardinassi, L. G., Ribeiro, J. M., Anatriello, E., Ferreira, B. R., Moreira, H. N., et al. (2014). The sialotranscriptome of Amblyomma triste, Amblyomma parvum and *Amblyomma cajennense* ticks, uncovered by 454-based RNA-seq. *Parasit. Vectors* **7**:430. doi: 10.1186/1756-3305-7-430

Gibbs, G. M., Roelants, K., and O'Bryan, M. K. (2008). The CAP superfamily: cysteine-rich secretory proteins, antigen 5, and pathogenesis-related 1 proteins—roles in reproduction, cancer, and immune defense. *Endocr. Rev.* **29**, 865–897. doi: 10.1210/er.2008-0032

Goodman, S. R., Petrofes Chapa, R., and Zimmer, W. E. (2015). Spectrin's chimeric E2/E3 enzymatic activity. *Exp. Biol. Med.* **240**, 1039–1049. doi: 10.1177/1535370215596084

Graca-Souza, A., Maya-Monteiro, C., Paiva-Silva, G., Braz, G. R., Paes, M. C., Sorgine, M. H., et al. (2006). Adaptations against heme toxicity in blood-feeding arthropods. *Insect Biochem. Mol. Biol.* **36**, 322–335. doi: 10.1016/j.ibmb.2006.01.009

Hara, S., Kudo, I., and Inoue, K. (1991). Augmentation of prostaglandin E2 production by mammalian phospholipase A2 added exogenously. *J. Biochem.* **110**, 163–165. doi: 10.1093/oxfordjournals.jbchem.a123550

Hobohm, U., and Sander, C. (1995). A sequence property approach to searching protein databases. *J. Mol. Biol.* **251**, 390–399. doi: 10.1006/jmbi.1995.0442

Hope, J. C., Kwong, L. S., Entrican, G., Wattegedera, S., Vordermeier, H. M., Sopp, P., et al. (2002). Development of detection methods for ruminant interleukin (IL)-12. *J. Immunol. Methods* **266**, 117–126. doi: 10.1016/S0022-1759(02)00113-8

Hu, H., and Sun, S. C. (2016). Ubiquitin signaling in immune responses. *Cell. Res.* **26**, 457–483. doi: 10.1038/cr.2016.40

Jensen, L. J., Gupta, R., Blom, N., Devos, D., Tamames, J., Kesmir, C., et al. (2002). Ab initio prediction of human orphan protein function from post-translational modifications and localization features. *J. Mol. Biol.* **319**, 1257–1265. doi: 10.1016/S0022-2836(02)00379-0

Jensen, L. J., Stærfeldt, H. H., and Brunak, S. (2003). Prediction of human protein function according to Gene Ontology categories. *Bioinformatics* **19**, 635–642. doi: 10.1093/bioinformatics/btg036

Jongejan, F., and Uilenberg, G. (2004). The global importance of ticks. *Parasitology* **129**(Suppl.), S3-S14. doi: 10.1017/S0031182004005967

Karim, S., Ramakrishnan, V. G., Tucker, J. S., Essenberg, R. C., and Sauer, J. R. (2004). *Amblyomma americanum* salivary gland homolog of nSec1 is essential for saliva protein secretion. *Biochem. Biophys. Res. Commun.* **324**, 1256–1263. doi: 10.1016/j.bbrc.2004.09.189

Karim, S., Singh, P., and Ribeiro, J. M. (2011). A deep insight into the sialotranscriptome of the gulf coast tick, *Amblyomma maculatum*. *PLoS ONE* **6**:e28525. doi: 10.1371/journal.pone.0028525

Kazimirová, M., and Štibrániová, I. (2013). Tick salivary compounds: their role in modulation of host defences and pathogen transmission. *Front. Cell. Infect. Microbiol.* **3**:43. doi: 10.3389/fcimb.2013.00043

Kim, T. K., Ibelli, A. M., and Mulenga, A. (2015). *Amblyomma americanum* tick calreticulin binds C1q but does not inhibit activation of the classical complement cascade. *Ticks Tick Borne Dis.* **6**, 91–101. doi: 10.1016/j.ttbdis.2014.10.002

Knutzova, P., Vancova, I., Kocakova, P., Slovak, M., Proost, P., and Kopacek, J. (2006). New member of the protein disulfide isomerase (PDI) family identified in *Amblyomma variegatum* tick. *Insect. Biochem. Mol. Biol.* **36**, 943–953. doi: 10.1016/j.ibmb.2006.09.005

Koh, C. Y., Kazimirova, M., Trimnell, A., Takac, P., Labuda, M., Nuttall, P. A., et al. (2007). Variegin, a novel fast and tight binding thrombin
inhibitor from the tropical bont tick. *J. Biol. Chem.* 282, 29101–29113. doi: 10.1074/jbc.M705600200

Kopecký, J., Kuthejlová, M., and Pechová, J. (1999). Salivary gland extract from *Ixodes ricinus* ticks inhibits production of interferon-gamma by the upregulation of interleukin-10. *Parasite Immunol.* 21, 351–356. doi: 10.1046/j.1365-3024.1999.00233.x

Kotál, J., Langhansová, H., Lieskovská, J., Andersen, J. F., Francischetti, I. M., Chavakis, T., et al. (2015). Modulation of host immunity by tick saliva. *J. Proteomics* 128, 58–68. doi: 10.1016/j.jprot.2015.07.005

Kovár, L., Kopecký, J., and Ríhová, B. (2001). Salivary gland extract from *Ixodes ricinus* tick polarizes the cytokine profile toward Th2 and suppresses proliferation of T lymphocytes in human PBMC culture. *J. Parasitol.* 87, 1342–1348. doi: 10.1645/0022-3395(2001)087[1342:SGEFIR]2.0.CO;2

Kovár, L., Kopecký, J., and Ríhová, B. (2002). Salivary gland extract from *Ixodes ricinus* tick modulates the host immune response towards the Th2 cytokine profile. *Parasitol. Res.* 88, 1066–1072. doi: 10.1007/s00436-002-0714-4

Kwong, L. S., Hope, J. C., Thom, M. L., Sopp, P., Duggan, S., Bembridge, G. P., et al. (2002). Development of an ELISA for bovine IL-10. *Vet. Immunol. Immunopathol.* 85, 213–223. doi: 10.1016/S0165-2427(02)00007-7

Kwong, L. S., Thom, M., Sopp, P., Rocchi, M., Wattegedera, S., Entrican, G., et al. (2010). Production and characterization of two monoclonal antibodies to bovine tumor necrosis factor alpha (TNF-alpha) and their cross-reactivity with ovine TNF-alpha. *Vet. Immunol. Immunopathol.* 135, 320–324. doi: 10.1016/j.vetimm.2010.01.001

Leboulle, G., Crippa, M., Decrem, Y., Mejri, N., Brossard, M., Bollen, A., et al. (2002). Characterization of a novel salivary immunosuppressive protein from *Ixodes ricinus* ticks. *J. Biol. Chem.* 277, 10083–10089. doi: 10.1074/jbc.M111391200

MacMicking, J., Xie, Q. W., and Nathan, C. (1997). Nitric oxide and macrophage functions. *Annu. Rev. Immunol.* 15, 323–350. doi: 10.1146/annurev.immunol.15.1.323

Madden, R. D., Sauer, J. R., and Dillwith, J. W. (2004). A proteomics approach to characterizing tick salivary secretions. *Exp. Appl. Acarol.* 32, 77–87. doi: 10.1023/B:APPA.0000018316.80224.54

Majetschak, M. (2011). Extracellular ubiquitin: immune modulator and endogenous opponent of damage-associated molecular pattern molecules. *J. Leukoc. Biol.* 89, 205–219. doi: 10.1189/jlb.0510316

Majetschak, M., Cohn, S. M., Obertacke, U., and Proctor, K. G. (2004). Therapeutic potential of exogenous ubiquitin during resuscitation from severe trauma. *J. Trauma* 56, 991–1000. doi: 10.1097/01.TA.0000127770.29009.5A

Malynn, B. A., and Ma, A. (2010). Ubiquitin makes its mark on immune regulation. *Immunity* 33, 843–852. doi: 10.1016/j.immuni.2010.12.007

Martens, L., Hermjakob, H., Jones, P., Adamski, M., Taylor, C., States, D., et al. (2005). PRIDE: the proteomics identifications database. *Proteomics* 5, 3537–3545. doi: 10.1002/pmic.200401303

Martinez, D., Barré, N., Mari, B., and Vidalenc, T. (1992). "Studies of the role of *Amblyomma variegatum* in the transmission of *Dermatophilus congolensis*," in *Tick Vector Biology*, eds B. Fivaz, T. Petney, and I. Horak (Berlin; Heidelberg: Springer), 87–99.

Merk, M., Zierow, S., Leng, L., Das, R., Du, X., Schulte, W., et al. (2011). The D-dopachrome tautomerase (DDT) gene product is a cytokine and functional homolog of macrophage migration inhibitory factor (MIF). *Proc. Natl. Acad. Sci. U.S.A.* 108, E577–E585. doi: 10.1073/pnas.1102941108

Miyake, A., Yamamoto, H., Enomori, T., and Kawashima, H. (1994). Exogenous group II phospholipase A2 induces prostaglandin E2 production in mouse peritoneal macrophages. *Eur. J. Pharmacol.* 253, 155–161. doi: 10.1016/0014-2999(94)90770-6

Mulenga, A., Kim, T. K., and Ibelli, A. M. (2013). Deorphanization and target validation of cross-tick species conserved novel *Amblyomma americanum* tick saliva protein. *Int. J. Parasitol.* 43, 439–451. doi: 10.1016/j.ijpara.2012.12.012

Olsen, J. V., de Godoy, L. M., Li, G., Macek, B., Mortensen, P., Pesch, R., et al. (2005). Parts per million mass accuracy on an Orbitrap mass spectrometer via lock mass injection into a C-trap. *Mol. Cell. Proteomics* 4, 2010–2021. doi: 10.1074/mcp.T500030-MCP200

Park, Y., Jin, H. S., Aki, D., Lee, J., and Liu, Y. C. (2014). The ubiquitin system in immune regulation. *Adv. Immunol.* 124, 17–66. doi: 10.1016/B978-0-12-800147-9.00002-9

Pockley, A. G. (2003). Heat shock proteins as regulators of the immune response. *Lancet* 362, 469–476. doi: 10.1016/S0140-6736(03)14075-5

Poddar, D., Basu, A., Baldwin, W. M. III, Kondratov, R. V., Barik, S., et al. (2013). An extraribosomal function of ribosomal protein L13a in macrophages resolves inflammation. *J. Immunol.* 190, 3600–3612. doi: 10.4049/jimmunol.12 01933

Poole, N. M., Mamidanna, G., Smith, R. A., Coons, L. B., and Cole, J. A. (2013). Prostaglandin E2 in tick saliva regulates macrophage cell migration and cytokine profile. *Parasit. Vectors* 6:261. doi: 10.1186/1756-3305-6-261

Prevot, P. P., Beschin, A., Lins, L., Beaufays, J., Grosjean, A., Bruys, L., et al. (2009). Exosites mediate the anti-inflammatory effects of a multifunctional serpin from the saliva of the tick *Ixodes ricinus*. *FEBS J.* 276, 3235–3246. doi: 10.1111/j.1742-4658.2009.07038.x

Radulović, Ž. M., Kim, T. K., Porter, L. M., Sze, S. H., Lewis, L., and Mulenga, A. (2014). A 24-48 h fed *Amblyomma americanum* tick saliva immuno-proteome. *BMC Genomics* 15:518. doi: 10.1186/1471-2164-15-518

R Core Team (2013). *R: A Language and Environment for Statistical Computing*. Vienna: R Foundation for Statistical Computing. Available online at: http://www.R-project.org/

Ribeiro, J., Makoul, G., and Robinson, D. (1988). Ixodes dammini: evidence for salivary prostacyclin secretion. *J. Parasitol.* 74, 1068–1069. doi: 10.2307/3282240

Ribeiro, J. M., Anderson, J. M., Manoukis, N. C., Meng, Z., and Francischetti, I. M. (2011). A further insight into the sialome of the tropical bont tick, *Amblyomma variegatum*. *BMC Genomics* 12:136. doi: 10.1186/1471-2164-12-136

Ribeiro, J. M., Zeidner, N. S., Ledin, K., Dolan, M. C., and Mather, T. N. (2004). How much pilocarpine contaminates pilocarpine-induced tick saliva? *Med. Vet. Entomol.* 18, 20–24. doi: 10.1111/j.0269-283X.2003.0469.x

Rojkind, M., Dominguez-Rosales, J., Nieto, N., and Greenwel, P. (2002). Role of hydrogen peroxide and oxidative stress in healing responses. *Cell. Mol. Life Sci.* 59, 1872–1891. doi: 10.1007/PL00012511

Rouault, T. A. (2004). Pathogenic bacteria prefer heme. *Science* 305, 1577–1578. doi: 10.1126/science.1102975

Scofield, S. L., Amin, P., Singh, M., and Singh, K. (2015). Extracellular ubiquitin: role in myocyte apoptosis and myocardial remodeling. *Compr. Physiol.* 6, 527–560. doi: 10.1002/cphy.c150025

Šimo, L., Kazimirova, M., Richardson, J., and Bonnet, S. I. (2017). The essential role of tick salivary glands and saliva in tick feeding and pathogen transmission. *Front. Cell. Infect. Microbiol.* 7:281. doi: 10.3389/fcimb.2017. 00281

Stachurski, F., Zoungrana, S., and Konkobo, M. (2010). Moulting and survival of *Amblyomma variegatum* (Acari: Ixodidae) nymphs in quasi-natural conditions in Burkina Faso: tick predators as an important limiting factor. *Exp. Appl. Acarol.* 52, 363–376. doi: 10.1007/s10493-010-9370-z

Stibrániová, I., Lahová, M., and Bartíková, P. (2013). Immunomodulators in tick saliva and their benefits. *Acta Virol.* 57, 200–216. doi: 10.4149/av_2013_02_200

Sujashvili, R. (2016). Advantages of extracellular ubiquitin in modulation of immune responses. *Mediators Inflamm.* 2016:4190390. doi: 10.1155/2016/4190390

Tanaka, K., Suzuki, T., and Chiba, T. (1998). The ligation systems for ubiquitin and ubiquitin-like proteins. *Mol. Cells* 8, 503–512.

Thouvenot, E., Urbach, S., Dantec, C., Poncet, J., Seveno, M., Demetre, E., et al. (2008). Enhanced detection of CNS cell secretome in plasma protein-depleted cerebrospinal fluid. *J. Proteome Res.* 7, 4409–4421. doi: 10.1021/pr8003858

Tian, Y., Chen, W., Mo, G., Chen, R., Fang, M., Yedid, G., et al. (2016). An immunosuppressant peptide from the hard tick *Amblyomma variegatum*. *Toxins* 8, 5. doi: 10.3390/toxins8050133

Vachieri, N., Puech, C., Cavelier, P., Rodrigues, V., Aprelon, R., Lefrançois, T., et al. (2015). An *in vitro* model to assess the immunosuppressive effect of tick saliva on the mobilization of inflammatory monocyte-derived cells. *Vet. Res.* 46, 117. doi: 10.1186/s13567-015-0229-5

Vider, J., Lehtmaa, J., Kullisaar, T., Vihalemm, T., Zilmer, K., Kairane, C., et al. (2001). Acute immune response in respect to exercise-induced oxidative stress. *Pathophysiology* 7, 263–270. doi: 10.1016/S0928-4680(00)00057-2

Wade, M. J. (2007). The co-evolutionary genetics of ecological communities. *Nat. Rev. Genet.* 8, 185–195. doi: 10.1038/nrg2031
Wanasen, N., Nussenzveig, R. H., Champagne, D. E., Soong, L., and Higgs, S. (2004). Differential modulation of murine host immune response by salivary gland extracts from the mosquitoes *Aedes aegypti* and *Culex quinquefasciatus*. *Med. Vet. Entomol.* 18, 191–199. doi: 10.1111/j.1365-2915.2004.00498.x

Wasala, N. B., Bowen, C. J., and Jaworski, D. C. (2012). Expression and regulation of macrophage migration inhibitory factor (MIF) in feeding American dog ticks, *Dermacentor variabilis*. *Exp. Appl. Acarol.* 57, 179–187. doi: 10.1007/s10493-012-9550-0

White, M. V., and Kaliner, M. A. (1987). Neutrophils and mast cells. I. Human neutrophil-derived histamine-releasing activity. *J. Immunol.* 139, 1624–1630.

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. doi: 10.3389/fmicb.2013.00337

Wink, D. A., Hines, H. B., Cheng, R. Y., Switzer, C. H., Flores-Santana, W., Vitek, M. P., et al. (2011). Nitric oxide and redox mechanisms in the immune response. *J. Leukoc. Biol.* 89, 873–891. doi: 10.1189/jlb.1010550

Zavašnik-Bergant, T., and Bergant Marušič, M. (2016). Exogenous thyropin from p41 invariant chain diminishes cysteine protease activity and affects IL-12 Secretion during maturation of human dendritic cells. *PLoS ONE* 11:e0150815. doi: 10.1371/journal.pone.0150815

Zinngrebe, J., Montinaro, A., Peltzer, N., and Walczak, H. (2014). Ubiquitin in the immune system. *EMBO Rep.* 15, 28–45. doi: 10.1002/embr.20138025

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Rodrigues, Fernandez, Vercoutere, Chamayou, Andersen, Vigy, Demettre, Seveno, Aprelon, Giraud-Girard, Stachurski, Loire, Vachiéry and Holzmuller. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.