BioMed Research International
Volume 2018 (2018), Article ID 9451547, 16 pages
https://doi.org/10.1155/2018/9451547

Research Article

Comparative Hemolymph Proteomic and Enzymatic Analyses of Two Strains of Rhipicephalus (Boophilus) microplus Ticks Resistant and Susceptible to Ixodicides

H. Aguilar-Díaz,¹ M. Esquivel-Velázquez,² R. E. Quiroz-Castañeda,³ E. Miranda-Miranda,¹ R. J. P. Conde-Baeye,⁴ M. Cobaxín-Cárdenas,³ P. Ostoa-Saloma,⁵ and R. Cossío-Bayúgar¹

¹Unidad de Artropodología, Centro Nacional de Investigación Disciplinaria en Parasitología Veterinaria, Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias, INIFAP. Carr. Fed. Cuernavaca-Cuautla No. 8534, Jiutepec, Morelos, 62550, Mexico
²Dirección de Investigación, Hospital General de México “Dr. Eduardo Liceaga”, Dr. Balmis 148, Col. Doctores, Ciudad de México, 06720, Mexico
³Unidad de Anaplasmosis, Centro Nacional de Investigación Disciplinaria en Parasitología Veterinaria, Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias, INIFAP. Carr. Fed. Cuernavaca-Cuautla No. 8534, Jiutepec, Morelos, 62550, Mexico
⁴Centro de Investigaciones Sobre Enfermedades Infecciosas, Instituto Nacional de Salud Pública, Av. Universidad 655, CP 62501 Cuernavaca, Morelos, Mexico
⁵Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México, Ciudad Universitaria, México 044510, D.F., Mexico

Correspondence should be addressed to H. Aguilar-Díaz; aguilar.hugo@inifap.gob.mx

Received 27 October 2017; Revised 9 March 2018; Accepted 8 April 2018; Published 11 June 2018

Academic Editor: Lachhman das Singla

Copyright © 2018 H. Aguilar-Díaz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The cattle tick Rhipicephalus (Boophilus) microplus is one of the most harmful ectoparasites affecting bovines worldwide. It represents a major threat to livestock industry due to the economic losses caused and diseases associated with these ticks. The most important tick control strategy has been the use of ixodicides, resulting in chemically resistant tick populations. It is necessary to understand the mechanisms that result in resistance so as to create new strategies increasing the lifespan of ixodicides or finding alternative targets to produce new acaricides. In this paper, in order to obtain an insight into the mechanisms that govern ixodicides resistance, we will compare the hemolymph proteome of two tick R. microplus strains, one susceptible (MJ) and one resistant (SA) to ixodicides, using HPLC and 2D electrophoresis. Significant differences were found in protein content between strains using HPLC. 2D electrophoresis revealed that 68 hemolymph protein spots were common between strains; however, 26 spots were unique to the susceptible strain MJ and 5 to the resistant strain SA. The most distinctive protein spots on the preparative gels were selected for further analyses. Nine protein spots were identified by mass fingerprinting,  revealing proteins that may have a role in the ixodicides resistance or susceptibility. In this paper, we present the tick hemolymph proteome revealing a set of proteins which suggest a possible role in tick detoxification.

1. Introduction

Ticks are considered obligate haematophagous ectoparasites that infest wild and domestic animals. They are responsible for significant economic losses mainly associated with mortality and morbidity of livestock animals [1].

In cattle, the worldwide distributed tick Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) represents an important threat to animal health and production [2]. This tick is especially distributed in tropical and subtropical environments where it infests bovines and is considered one of the most harmful vectors of tick-borne diseases in cattle [3]. In endemic areas with R. microplus, major economic losses are due to anemia which leads to a reduction in weight gain and milk production. Indirectly, the main risk is the transmission of pathogens such as Anaplasma and Babesia spp., among other bacterial and fungal infections [4].

Tick infestations in bovine cattle are usually treated with chemical ixodicides applied directly to the skin of animals. However, uncontrolled and excessive use of these compounds exerts a selection pressure on those ticks with ixodicides multiresistance. This is a well-known widespread phenomenon that increases the occurrence of food and environmental contamination [2, 5–7]. In general, the mechanism of resistance is usually associated with increased metabolic detoxification or target site modification, as previously described [7–10]. Nonetheless, there is still much more to elucidate regarding tick physiology and immunology. Therefore, efforts can be directed to identify novel biomolecules candidates for tick control using other approaches, such as genomics, transcriptomics, and proteomics (omics). The main contribution of genomic approaches in the study of vector-borne diseases is the design of new peptides identified by in silico analysis of transcriptomic and genomic data so as to develop novel vaccines against ticks [11, 12]. For the last ten years, “omics” approaches have provided valuable information regarding tick–host interface immunobiology, mainly by means of genomic, transcriptomic, and proteomic data of several species of ticks, regularly using larvae or engorged females [3, 13–16]. However, until the last year, the first tick I. scapularis genome was released and new possibilities of study were opened, especially those related to host–tick–pathogen interactions, parasitic processes unique to ticks, tick reproduction, and so forth [17, 18].

In R. microplus, a proteomic analysis of saliva (sialomes) showed differences in saliva protein profiles and composition between partially engorged (PE) and fully engorged female ticks. This was the first proteomic study of tick saliva that contributed to understanding the role of tick salivary modulators in immunological defense of ticks. Proteins identified in this proteomic analysis included microplusin-like proteins, lipocalins, serpins, and hemelipoprotein (HeLp) [3, 7]. Hemelipoproteins have been reported as the most abundant proteins in tick saliva and in hemolymph in several tick species, in which transcriptional profile and protein localization analyses suggested that this protein may play vital roles in tick feeding and survival [19, 20]. Other proteomic studies have also been performed in Amblyomma americanum (salivary secretions), I. scapularis (saliva protein composition), Ornithodoros moubata (tick saliva), and R. sanguineus (sialoma) [19, 21–23]. So far, our knowledge of hemolymph proteins in ticks is still scarce, and its relevance lies in the roles this fluid plays in this arthropod. Hemolymph is a fluid that bathes all tissues in ticks. It is the first source of nutrients, transports molecules and hormones, and provides protection to pathogen agents to which ticks are exposed [24]. The hemolymph, as a vehicle of transport of factors related to humoral immune response, has been explored in insects and crustaceans but not in ticks, where it may have a significant role in defense and protection not only to microbial invasions but also to environmental factors such as pesticides [25]. Likewise, the relationship of the hemolymph, immune response, and susceptibility-resistance genotype still remains unclear. Also, to date, no “omics” studies on tick hemolymph have been reported, and much less, the role of hemolymph is unknown. In this work, we present proteomic profiles of hemolymph from resistant and susceptible ticks. The analysis of the proteomic profiles provides useful information for understanding the relationship between ixodicides and the resistance/susceptibility mechanisms of ticks along with the identification of proteins which may have key roles. The goal of this study is to find those proteomic differences between the hemolymph from ticks that are ixodicides-resistant when compared to reference strains susceptible to SDS/PAGE.

2. Materials and Methods

2.1. Tick Strains

The ixodicide susceptible R. microplus tick strain Media Joya (MJ) and the multiple ixodicides-resistant tick strain San Alfonso (SA) [26] were used in this study (Table 1). The strains have been maintained in controlled infested bovines for many generations and used as reference for the tick ixodicides resistance monitoring programs of the Mexican Federal Government. The resistant strain SA was reared and maintained at Departamento de Ectoparásitos y Dípteros del Servicio Nacional de Sanidad, Inocuidad y Calidad Agroalimentaria (SENASICA-SAGARPA), and susceptible strain MJ at Centro Nacional de Investigación Disciplinaria en Parasitología Veterinaria (CENID-PAVET, INIFAP). Each reference strain was obtained by infesting a bovine with 2 × 10⁴ larvae, aged between 10 and 15 days, and then engorged tick females were collected 21 days after infestation and placed in Petri dishes in groups of 10 for immediate hemolymph extraction as reported previously by [27].

Table 1: Percentage of mortality of susceptible MJ and resistant SA tick strains to ixodicides.

Percentage of mortality (%)

Tick Strain
Organophosphorous
Pyrethroids
Amidines

Coumaphos
Diazinon
Chlorpyrifos
Cypermethrin
Deltamethrin
Flumethrin
Amitraz

Susceptible (MJ)
100
100
100
100
100
100
100

Resistant (SA)
73.24
100
100
0
0
0
10.76

2.2. Bioassays on Ixodicides Discriminant Doses

Bioassays were run as reported previously by [28]. The ixodicides were diluted in trichloroethylene (Sigma-Aldrich) at the following concentrations: coumaphos 0.2%, chlorfenvinphos 0.2%, diazinon 0.08%, chlorpyrifos 0.2%, cypermethrin 0.05%, deltamethrin 0.09%, and flumethrin 0.01%. One milliliter of each dilution was applied evenly to a cm piece of filter paper (Whatman, Sigma-Aldrich). The trichloroethylene was allowed to evaporate from the filter paper for 2 h. The filter papers were then folded in half and sealed on the slides with clips, which formed a packet into which approximately 100 larvae were placed and then the top of the packet was sealed with another clip. The packets were kept at 27°C and 92% relative humidity for 24 h. The packets were removed from incubation and opened, live and dead larvae were counted, and the data was processed as percentage of mortality for each tick group under every ixodicide concentration.

2.3. Protein and Hemolymph Extraction and Quantification

Engorged tick females were placed in Petri dishes and observed under a dissecting microscope. The hemolymph was obtained by performing an incision with sterile scalpels in the anterior segment of ticks, very carefully to not damage any internal organ. The hemolymph was mixed with an equal volume of collection buffer (98 mM NaOH, 186 mM NaCI, 1.7 mM EDTA, and 41 mM citric acid, pH 4.5), in sterile tubes (Stoepler et al. 2012) and protease inhibitor (Roche). Then, the hemolymph was frozen and thawed three times by placing the tubes in liquid nitrogen for 30 seconds, followed by incubation at 37°C for 3 minutes, in order to lyse any cell. Extracts were then centrifuged at 14000xg for 10 minutes at 4°C. The supernatant was collected, quantified by conventional Bradford method (Bio-Rad), and stored at -70°C [29]. This protocol allowed the collection of 10-20 μl of hemolymph per tick, reducing free of detritus cell, melanization rest, tissue debris, or other contaminants.

2.4. HPLC Sample Preparation

The hemolymph from resistant strain SA and susceptible strain MJ was collected as previously described and supplemented with protease inhibitors with 0.01% phenyl thiourea, 1 µg/ml leupeptin 10 µl/ml, TLCK (L-1-chloro-3-(4-tosylamide)), and 10 µl/ml and 1 µl/ml PMSF (all supplied from Roche). The resultant supernatants were centrifuged at 12000xg for 30 minutes and then adjusted to 0.1% with trifluoroacetic acid (TFA) (Sigma-Aldrich).

2.5. HPLC Fractioning

The samples were applied on a C-18 reverse phase HPLC column on an Agilent 1100 series HPLC apparatus and eluted with a 0 to 50 % acetonitrile gradient (ACN) (SIGMA) with 25 minutes of collection and then with 50 to 100% ACN gradient within 20 minutes of collection. The collected fractions were lyophilized and resuspended in 50 μl of phosphate-buffered saline 1X.

2.6. Zymograms Assays

70 ng of protein extract from each hemolymph sample was electrophoresed as previously described [30]. After electrophoresis, gels were stained for detection of catalase and phosphatase enzymatic activity as described previously [31]. 70 ng of each protein extract was assayed by the Sigma® catalase assay kit and Sigma® phosphatase colorimetric assay kit according to the manufacturer’s instructions. Total proteins were stained by Coomassie brilliant blue R250 (Sigma-Aldrich) 0.5% for reference.

2.7. Determination of Phenol Oxidase (PO) Production

Hemolymph fractions from HPLC fractioning were probed for PO activity using l-DOPA (dihydroxyphenylalanine, Sigma-Aldrich) chromogenic reaction, as described elsewhere (29). 5 µl of the HPLC fractions was mixed with 50 µl of 0.13 mg 3-4-dihydroxy-L-phenylalanine/ml distilled water. The solution was incubated for one hour at ambient temperature and the 620 nm absorption was recorded every ten minutes. The slope of the 620 nm OD increase of every fraction tested was compared with the control reaction (H₂O) and samples with consistent OD increases were graphed in Figure 1.

(a)

(b)

Figure 1: HPLC chromatography profiles of susceptible MJ (a) and resistant SA tick strains (b). At higher concentrations of ACN (acetonitrile), multiple peaks eluted in susceptible strain MJ ((a), upper), while at low concentrations of ACN, peaks were observed in resistant strain SA ((b), upper). The oxidizing activity exhibited through L-DOPA oxidation is depicted below the chromatograms ((a) and (b), lower); the graph exhibits an increase of oxidative activity approximately at 10 min retention time in the susceptible MJ when compared to the resistant SA; additionally, the resistant strain SA exhibits several peaks of oxidative activity which are not present in the susceptible strain MJ during the 20 to 30 min retention time.

2.8. Two-Dimensional Electrophoresis (2D-E) and Gel Analysis

For analytical 2D gels, 35 µg of hemolymph proteins was loaded onto broad range lineal 7 cm IPG strips (pH 3-10, GE Healthcare) in rehydration buffer (8 M urea, 2% w/v CHAPS, 0.8% w/V DTT, 1.6% v/v Bio-Lyte Broad Range Ampholytes, and 0.002% w/v bromophenol blue) and left to rehydrate for 12 hours at room temperature. Isoelectric focusing was performed with PROTEAN i12 IEF System for a total of 10,000 V-h in four steps: 125 V for 2 h, rapid; 250 V for 20 min, gradual; 4000 V for 2 h, rapid; and 4000 V for 10,000 V-h, rapid.

Strips were then equilibrated with 1% w/v DTT (Promega) and then with 2.5% w/v iodoacetamide (Sigma-Aldrich) in fresh equilibration buffer (6 M urea, 0.375 M Tris-HCl, pH 8.8 (Sigma-Aldrich), 2% w/v SDS, and 20% v/v glycerol) for 10 minutes each under constant rocking. The second-dimension electrophoresis was performed in 4-20% precast gels (Bio-Rad) in denaturing conditions. Gels were then fixed with buffer containing acetic acid 7% and methanol 10% for 30 minutes.

Protein spots were stained with Sypro Ruby (Bio-Rad) according to manufacturer’s instructions. Stained gels were scanned with Typhoon FLA 9500 imaging system (GE Healthcare) and spot analysis was performed with PDQuest software (Bio-Rad). All samples were run in triplicate. Only spots that replicated in at least two out of the three replicates were considered for analysis. Resistance and susceptibility conditions were compared and spot differences were identified.

In order to pick spots for protein identification, samples were run in duplicate. 144 µg of hemolymph proteins from resistant strain (SA) ticks and 64.88 µg from susceptible (MJ) ticks were loaded onto IPG strips (pH 3-10, lineal) and proteins were separated as previously stated. Gels were fixed and stained with Coomassie G-250. Gel spot picking was performed manually with custom-cut pipet tips over a light box.

2.9. Protein Identification Analysis

A total of 13 protein spots were selected for digestion sequence and mass spectrometry analysis at the University Proteomics Laboratory at Biotechnology Institute of the National Autonomous University of México. In gel digestion of proteins was performed with trypsin (Promega) in reaction buffer containing ammonium bicarbonate 50 mM for 18 hours at 37°C. Peptides were desalted with Zip Tip C18 (Millipore) and applied to a LC-MS system (Nanoflux Pump EASY nLC II and LTQ-Orbitrap Velos, Thermo Fischer) with nanoelectrospray ionization (ESI). Proteins were identified automatically by Proteome Discoverer 1.4 software through the search engines SEQUEST-HT. Boophilus, Rhipicephalus, and protein databases obtained from NCBI were used for  identity searches. Spot analyses were performed with InterPro, SMART (Simple Modular Architecture Research Tool), ProtoNet, Pfam, and SPRINT software.

3. Results

3.1. Bioassays on Ixodicide Discriminant Doses

Resistant strain SA was assayed by using the larval package procedure for discriminant ixodicides doses previously reported [28]; the bioassay results were used as selection criterion for ticks from where the hemolymph was obtained. Once the packets with ixodicides were opened, live and dead larvae were counted and percentages of mortality were calculated. In Table 1, percentages of mortality are shown for susceptible strain MJ and resistant strain SA. 100% mortality of susceptible strain MJ was observed with all ixodicides assayed, including organophosphorous, pyrethroids, and amidines. Multiple ixodicides-resistant strain SA showed resistance to the three pyrethroids used (cypermethrin, flumethrin, and deltamethrin) with 0% mortality; on the contrary, 100% mortality was observed with two organophosphorous compounds (diazinon and chlorpyrifos) and 73.24% with coumaphos. This strain also showed 10.76% mortality with amidines (diazinon).

3.2. HPLC Chromatography

Spectroscopy profiles of hemolymph at 280 nm showed different patterns (Figure 1); the profile of the susceptible strain MJ exhibited an increasing complexity when eluted at higher ACN percentage (Figure 1(a), upper), whereas the profile of the resistant strain SA exhibited a more complex pattern when eluted at lower ACN percentage (Figure 1(b), upper).

3.3. Zymograms Assays

Catalase and phosphates are key enzymes during xenobiotics metabolism of most metazoan organisms; these enzymes, in combination with other xenobiotic metabolizing enzymes, are required by pesticide resistant arthropods for metabolic degradation of toxic chemicals into less damaging compounds. We considered a possible differential expression of these enzymes when comparing susceptible against resistant ticks.

3.3.1. Phosphatase Activity

The protein bands pattern with phosphatase activity varied between hemolymphs from susceptible MJ and resistant SA strains as observed in Figure 2(a).

(a)

(b)

Figure 2: Phosphatase and catalase activity gels of hemolymph proteins of susceptible MJ and resistant SA tick strains. (a) In phosphatase zymogram, two protein bands in hemolymph of resistant strain SA were observed while only one protein band was observed in susceptible strain MJ. (b) In catalase zymogram, the protein band patterns of resistant strain SA and susceptible strain MJ hemolymph are different. In resistant strain SA, seven protein bands with catalase activity were observed and only three of these bands were observed in susceptible strain MJ. Arrows indicate protein band position. Larvae hemolymphs of resistant SA and susceptible MJ strains were used as a known positive control of these activities.

In hemolymph extracts from resistant strain SA, protein bands with enzymatic activity were observed (molecular weight approx. 90 and 170 kDa) and only one protein band with phosphatase activity was observed in susceptible strain MJ hemolymph (molecular weight approx. 90 kDa). A different protein band pattern was also observed between resistant strain SA and susceptible strain MJ larvae extracts.

3.3.2. Catalase Activity

Catalase activity was assessed in a zymogram of hemolymph and larvae extracts of both resistant SA and susceptible MJ strains as shown in Figure 2(b).

In hemolymph extracts of resistant strain SA, catalase activity in hemolymph was assigned to seven protein bands with enzymatic activity (molecular weight approx. 8, 17, 50, 60, 90, 110, and 170 kDa), while in susceptible strain MJ hemolymph extracts, this number of bands decreased to three (molecular weight approximately 8, 90, and 170 kDa).

3.4. Phenol Oxidase (PO) Activity

The HPLC fractions showed differential PO activity: susceptible strain MJ hemolymph showed two fractions showing activity eluting from the column around 10% ACN (Figure 1(a), lower). Resistant strain SA hemolymph fractionation eluted samples with oxidizing activity in the first fractions (0.5% ACN elution) but also at 33%, 40%, and 50% ACN (Figure 1(b), lower). These multiple peaks of oxidizing activity in the resistant strain SA hemolymph might reflect distinct proteolytic activation of the same enzyme (the hydrophilic protein peak with oxidizing activity eluted at 0.5% ACN elution for the resistant strain SA and 10% ACN for the susceptible strain MJ). Also, the more hydrophobic peaks showing oxidizing activity could correspond to a set of enzymes unique to the susceptible strain MJ.

3.5. Protein Identification Analysis

In order to further explore the differences in the protein from resistant strain SA and strain susceptible MJ R. microplus, we compared their proteomes by means of 2D electrophoresis. The proteome images were analyzed with PDQuest software to detect spots and compare the triplicates. 68 protein spots were common to both strains, 26 spots were unique of the susceptible strain MJ, and 5 were unique to the resistant strain SA. From the spots unique to each strain, we selected those clearly visible in preparative gels. The selected spots were gel-picked and sent for identification by mass fingerprinting (Figure 3). A total of 13 spots were sent for identification, from which 9 were successfully identified (Figure 4, Table 2).

Table 2: Identified proteins in hemolymph from susceptible MJ and resistant SA R. microplus tick strains.

Spot ID
Accession number
Description
Score (%)
Coverage (%)
Theoretical MW (kDa)/IP
Function
Reference

SUSMJ-5204
JAC59075.1 (GenBank)
Microplusin-like 2, partial (Rhipicephalus microplus)
495.58
31.63
10.7/6.98
Antimicrobial peptide isolated from the cattle tick R. microplus. Its copper-chelating ability is putatively responsible for its bacteriostatic effect.
[32]

SUSMJ-6205
JAC59122.1 (GenBank)
Hemelipoprotein HeLp2, partial (Rhipicephalus microplus)
29.67
5.99
30.5/6.65
Protein involved in heme transports from hemolymph to tissues, suggesting to be an essential adaptation to the loss of the heme synthesis pathway
[33]

SUSMJ-6208
JAP78009.1 (GenBank)
Hypothetical protein (Rhipicephalus appendiculatus)
7.25
20.12
18.9/6.77
Protein binding to EGF. The functional significance of EGF domains in what appears to be unrelated proteins is not yet clear. However, a common feature is that these repeats are found in the extracellular domain of membrane-bound proteins or in proteins known to be secreted
[34]

SUSMJ-6302
JAC59082.1 (GenBank) 

JAC59119.1 (GenBank)
Secreted protein 33, partial (Rhipicephalus microplus) 

Secreted protein 27, partial (Rhipicephalus microplus)
10.08 

8.31
3.83 

3.86
37.8/7.99 

37.4/7.28
Leucine-rich repeat domain, L domain-like. Leucine-rich repeats (LRR) consist of 2-45 motifs of 20-30 amino acids in length that generally fold into an arc or horseshoe shape. LRRs occur in proteins ranging from viruses to eukaryotes, and appear to provide a structural framework for the formation of protein-protein interactions. Proteins containing LRRs include tyrosine kinase receptors, cell-adhesion molecules, virulence factors, and extracellular matrix-binding glycoproteins, and are involved in a variety of biological processes, including signal transduction, cell adhesion, DNA repair, recombination, transcription, RNA processing, disease resistance, apoptosis, and immune response
[35–38]

SUSMJ-6304
JAC59079.1 (GenBank)
Secreted protein 36, partial (Rhipicephalus microplus)
23.30
16.26
21.6/7.12
Copper/zinc superoxide dismutase (SOD). superoxide dismutases catalyse the conversion of superoxide radicals to molecular oxygen.
[32]

JAC59082.1 (GenBank)
Secreted protein 33, partial (Rhipicephalus microplus)
16.45
3.83
37.8/7.99
Uncharacterized.
The component of a membrane consisting of the gene products and protein complexes having at least some part of their peptide sequence embedded in the hydrophobic region of the membrane. Reported by Tirloni 2014.
[3]

De novo
100% identity to phospholipases

100
---
Phospholipase-mediated phospholipid hydrolysis leads to the production of lipid mediators or second messengers that affect signal transduction, thus regulating a variety of physiologic and pathophysiologic processes.
[39]

SMJ-7106
 P02070 (UniProtKB/Swiss-Prot)

2QSP (PDB) 

 2EU1 (PDB)

Hemoglobin subunit beta OS=Bos Taurus GN=HBB PE=1 SV=1 (HBB_BOVIN) 

 C Chain C, Bovine hemoglobin At Ph 5.7 

D Chain D, Crystal structure determination of goat hemoglobin (Capra hircus) at 3 Angstrom resolution

1919.94

163.98

1288.43

82.76

78.72

55.17

15.9/7.59

15.0/8.44

16.0/7.30
Haemoglobins (Hb) function is transport oxygen in blood plasma, and binds oxygen in the reduced [Fe(II)] state. (Hb) tetramer of two alpha and two beta chains, although embryonic and foetal forms can substitute the alpha or beta chain for ones with higher oxygen affinity, such as gamma, delta, epsilon or zeta chains. Hb transports oxygen from lungs to other tissues in vertebrates. Hb proteins are also present in unicellular organisms where they act as enzymes or sensors.
[40–42]

RESSA-3008
JAC59001.1 (GenBank)
TIL domain-containing protein 2 (Rhipicephalus microplus)
1193.73
52.94
9.2/5.77
Different analyses carried out show that thia protein contains a domain 25 -81 TIL (trypsin inhibitory-like) This domain is found in proteinase inhibitors, as well as in many extracellular proteins. The domain typically contains ten cysteine residues that form five disulphide bonds. The cysteine residues that form the disulphide bonds are 1-7, 2-6, 3-5, 4-10 and 8-9.
[3]

ACV83329.1 (GenBank)
Serine proteinase inhibitor precursor (Rhipicephalus microplus)
872.04
67.06
9.3/6.24
Serpin family members have activity as inhibitors of serine proteinases, and serve in regulation of proteinase-mediated, such as inflammatory, and the immune response with consequent pathology in deficiency state.
[43]

ABH10604.1
(GenBank)
Neutrophil elastase inhibitor (Rhipicephalus microplus)
210.52
21.15
11.5/6.51
In ticks, these inhibitors (due to their properties) play a role in the feeding process, however the crucial role is still unclear.
[44]

RESSA-5214
JAC59087.1 (GenBank)
Secreted protein 31, partial (Rhipicephalus microplus)
17.73
2
15.6/6.52
Unknown. Reported by Tirloni 2014.
[45]

(a)

(b)

Figure 3: 2D electrophoresis gel of hemolymph proteins of susceptible MJ and resistant SA tick strains. (a) 26 proteins were identified by PDQuest software in susceptible strain MJ and 5 proteins in resistant strain SA. 28 proteins were identified as common proteins of both strains. Yellow and red squares correspond to those protein spots that were selected for sequencing. (b) 3D surface plot of 2D electrophoresis. Protein concentration differences are observed in each plot.

Figure 4: 3D visualization of the identified spots. The same area from representative gels from resistant SA and susceptible MJ tick strains was selected and visualized as 3D graphs. The identified spots are indicated with arrows in the images and protein identity is indicated to the right.

4. Discussion

The massive use of ixodicides from the control of the cattle tick is collaterally producing the right conditions for the selection of resistance. This adaptation is attributed to resistance mechanisms that ticks have to use in order to survive in new and hostile environments, which results in serious problems of chemical tick control for the cattle farmers [46]. Hemolymph is required by all arthropods as blood is for vertebrate organisms; it provides osmotic balance and a mechanism of nutrients and oxygen distribution and also acts as the liquid medium for the propagation of circulating cells and molecules including those related to the immune system and detoxification of pesticides. In this work, we explored the role of hemolymph proteins of susceptible and resistant tick strains with a proteomic approach.

Besides the immunogenic function attributed to hemolymph during pathogen invasion, it also possesses elements used to detoxify ticks when they are exposed to ixodicides. This mechanism is a complex process that includes enzymatic reactions mediated by phosphatases, catalases, and glutathione S transferases [47, 48], which modify metabolic pathways resulting in metabolic resistance. During this study, we used SDS-PAGE zymograms and found that those polypeptide bands that exhibit specific catalase activity were different when comparing zymograms from resistant strain SA and susceptible strain MJ. The strain SA is regarded as pyrethroid resistance reference strain. According to previous studies, this type of resistance is based on overexpression of oxidases and cytochrome P450 monooxygenases acting in a functional enzymatic combination [49, 50]. The catalase and phosphatase zymogram analysis showed that both strains have a differential pattern of protein bands with enzymatic activity that may be related to ixodicides resistance; whether this is an effective biochemical marker for detection in cattle tick is to be further assessed in future studies, besides a significant role of other enzymes such as glutathione S transferases, oxidases, esterases, and cytochrome P450 monooxygenases [50–52] in ixodicides resistance of cattle ticks. In fact, other oxidase enzymes such as phenol oxidase (PO) have been reported in insects (Anopheles spp., Triatoma spp.) as part of the melanization process of the protective immune response. However, the role of PO in ticks still remains unknown [53]. So, we tried to correlate the enzymatic activity involved in detoxification with resistant and susceptible phenotypes. We observed higher oxidative enzymatic activity in the hemolymph of the susceptible strain MJ, when compared to the resistant strain SA. It is plausible that cypermethrin/flumethrin/deltamethrin resistance presented in the susceptible strain SA reduced PO activity may be associated with the necessity to metabolically modify pesticides in order to acquire an increased toxic effect [50]. Also, HPLC chromatography profiles of susceptible strain MJ and resistant strain SA revealed that a higher number of peaks were eluted at higher ACN concentrations in susceptible strain MJ; on the contrary, a more complex HPLC profile was obtained at low ACN concentration in resistant strain SA. These results suggested that hemolymph of resistant strain SA has more hydrophilic components than hemolymph of susceptible strain MJ.

In order to explore hemolymph protein differential expression, we carried out 2D electrophoresis. Different spots were identified in protein maps of both strains. We observed a remarkable difference in the amount of spots, especially in resistant strain SA where this number was significantly lower than in the susceptible strain MJ. So far, evolution theories propose that organisms exposed to adverse environmental conditions have to deal with a trade-off of energy and having limited resources for growth, reproduction, and defense to cope with those conditions [54]. Although this condition may be advantageous under certain circumstances, it also can simultaneously reduce several functions related to physical and chemical limitations, as well as to genetic mechanisms and energy assignment (antagonistic pleiotropy) [55, 56]. It is likely that growth, reproduction, and defense performance are linked to a reduced capacity in ixodicides resistance. According to this, the diversity in the number of spots in protein maps (resistant and susceptible phenotypes) suggests that such high number observed is not necessarily associated with a robust resistance mechanism, since we found that the resistant strain SA has fewer spots than the susceptible strain MJ. At this point, we only identified the protein maps of ticks under two conditions; however, the amount of proteins that may directly be associated with resistance or susceptibility is not clear.

Proteomic comparison analysis of both strains revealed that 68 proteins are common to resistant SA and susceptible MJ strains, 26 spots are only present in susceptible strain MJ, and 5 spots are present in resistant strain SA. Of those spots that were exclusively present in susceptible strain MJ, we identified a protein with 100% identity to phospholipases that has not been previously reported. Another protein identified was hemelipoprotein HeLp2 that was overexpressed in susceptible strain MJ. This protein has been reported as part of the tick sialome and considered one of the more abundant proteins in salivary glands of ticks and also of the hemolymph [3]. Since expression of HeLp2 is lower in resistant strain SA, it is likely that this protein does not have a crucial role in resistance mediated by hemolymph components.

In the susceptible strain MJ, we also identified microplusin-like 2 protein, one of the main components of innate immune response [45]. We speculate that this protein is present under susceptibility conditions and not in resistance conditions, probably as a result of the trade-off that ticks would be performing when they are exposed to ixodicides. Finally, protein SUSMJ-6304 identified in the susceptible strain MJ is also relevant, since in silico analysis showed that sequence matched with Cu-Zn superoxide dismutase (SOD) in an arthropods database (UnitProt), an enzyme involved in cellular detoxification processes by removing free radicals of many species. Nevertheless, it is known that SOD is highly present in sialome of ticks, so probably the enzyme identified in the susceptible MJ and resistance SA strains would be a leakage from salivary glands to hemolymph. Previously, we measured SOD activity by enzymatic assays of hemolymph of both strains with negative results (data not shown). However, many reports have shown the role of free radicals in the inhibition of some insecticides, which may explain why the spot SUSMJ-6304 of susceptible strain MJ has a relative abundance higher than the resistance strain SA.

In this work, we report the first proteomic approach of tick hemolymph in Mexican R. microplus strains. We identified both previously and never reported proteins, as well as enzymatic activities related to ixodicides susceptibility or resistance, which may be used as future targets for diagnostic, drugs, and vaccines design to contribute to ticks control.

Future perspectives of this work will be driven to develop targets that could be used in prevention and therapeutic methods for a tick integral control.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this article.

Acknowledgments

The authors want to thank César Alejandro Arreguín-Pérez for hemolymph’s extraction assistance. Hugo Aguilar and Mayra Cobaxin acknowledge the support given by Programa de Retenciones, CONACyT, no. MOD-ORD-27 PCI-187- 11-15. R.C.B. Project PN (CONACyT) 248049, and INIFAP, CENID-PAVET (Project SIGI: 1-1.6-13512934133-p-m.2-1). The authors also thank Salvador Neri-Orantes and Francisco Martínez-Ibañez (CENAPA-SENASICA) for providing SA strain discriminant doses bioassays and Ricardo Hernández-Ávila (IIB, UNAM) for his valuable technical assistance in proteomic analysis.

References

A. Domingos, S. Antunes, L. Borges, and V. E. Rosario, “Approaches towards tick and tick-borne diseases control,” Revista da Sociedade Brasileira de Medicina Tropical, vol. 46, pp. 265–269, 2013, Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0037-86822013000300265&nrm=iso.

M. Popara, M. Villar, L. Mateos-Hernández, I. G. Fernández de Mera, and J. de la Fuente, “Proteomics approach to the study of cattle tick adaptation to white tailed deer,” BioMed Research International, vol. 2013, 8 pages, 2013.

L. Tirloni, J. Reck, R. M. S. Terra et al., “Proteomic analysis of cattle tick Rhipicephalus (Boophilus) microplus saliva: A comparison between partially and fully engorged females,” PLoS ONE, vol. 9, no. 4, Article ID e94831, 2014.

N. N. Jonsson, R. E. Bock, and W. K. Jorgensen, “Productivity and health effects of anaplasmosis and babesiosis on Bos indicus cattle and their crosses, and the effects of differing intensity of tick control in Australia,” Veterinary Parasitology, vol. 155, no. 1-2, pp. 1–9, 2008.

J. De La Fuente, M. Rodríguez, M. Redondo et al., “Field studies and cost-effectiveness analysis of vaccination with Gavac(TM) against the cattle tick Boophilus microplus,” Vaccine, vol. 16, no. 4, pp. 366–373, 1998.

L. Lovis, J. Reggi, M. Berggoetz, B. Betschart, and H. Sager, “Determination of acaricide resistance in rhipicephalus (Boophilus) microplus (Acari: Ixodidae) field populations of Argentina, South Africa, and Australia with the larval tarsal test,” Journal of Medical Entomology, vol. 50, no. 2, pp. 326–335, 2013.

L. F. Parizi, P. C. Pohl, and A. Masuda, “New approaches toward anti-Rhipicephalus (Boophilus) microplus tick vaccine,” Revista Brasileira de Parasitologia Veterinária, vol. 18, pp. 1–7, 2009, Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1984-29612009000100001&nrm=iso.

R. Rosario-Cruz, C. Almazan, R. J. Miller, D. I. Dominguez-Garcia, R. Hernandez-Ortiz, and J. De La Fuente, “Genetic basis and impact of tick acaricide resistance,” Frontiers in Bioscience, vol. 14, no. 7, pp. 2657–2665, 2009.

F. D. Guerrero, L. Lovis, and J. R. Martins, “Acaricide resistance mechanisms in Rhipicephalus (Boophilus) microplus,” Revista Brasileira de Parasitologia Veterinária, vol. 21, pp. 1–6, 2012, Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1984-29612012000100002&nrm=iso.

A. Y. Li, R. B. Davey, R. J. Miller, and J. E. George, “Detection and characterization of amitraz resistance in the Southern Cattle tick, Boophilus microplus (Acari: Ixodidae),” Journal of Medical Entomology, vol. 41, no. 2, pp. 193–200, 2004.

R. Miller, A. Estrada-Peña, C. Almazán et al., “Exploring the use of an anti-tick vaccine as a tool for the integrated eradication of the cattle fever tick, Rhipicephalus (Boophilus) annulatus,” Vaccine, vol. 30, no. 38, pp. 5682–5687, 2012.

K. B. Temeyer, A. C. Chen, R. B. Davey, F. D. Guerrero, J. M. Howell, and D. M. Kammlah, “Nuevos enfoques para el control de Rhipicephalus (Boophilus) microplus,” Revista Mexicana de Ciencias Pecuarias, vol. 3, Suppl 1, pp. 25–40, 2013.

J. M. C. Ribeiro, M. B. Labruna, B. J. Mans et al., “The sialotranscriptome of Antricola delacruzi female ticks is compatible with non-hematophagous behavior and an alternative source of food,” Insect Biochemistry and Molecular Biology, vol. 42, no. 5, pp. 332–342, 2012.

F. F. Cardoso, C. C. G. Gomes, B. P. Sollero et al., “Genomic prediction for tick resistance in braford and hereford cattle,” Journal of Animal Science, vol. 93, no. 6, pp. 2693–2705, 2015.

J. Pagel Van Zee, N. S. Geraci, F. D. Guerrero et al., “Tick genomics: The Ixodes genome project and beyond,” International Journal for Parasitology, vol. 37, no. 12, pp. 1297–1305, 2007.

M. Rodriguez-Valle, A. Lew-Tabor, C. Gondro et al., “Comparative microarray analysis of Rhipicephalus (Boophilus) microplus expression profiles of larvae pre-attachment and feeding adult female stages on Bos indicus and Bos taurus cattle,” BMC Genomics, vol. 11, no. 1, article no. 437, 2010.

A. B. Nuss, M. G. Mathew, and M. Gulia-Nuss, “Genomic insights into the Ixodes scapularis tick vector of Lyme disease,” Nature Communications, 2016.

J. de la Fuente, R. M. Waterhouse, D. E. Sonenshine et al., “Tick Genome Assembled: New Opportunities for Research on Tick-Host-Pathogen Interactions,” Frontiers in Cellular and Infection Microbiology, vol. 6, 2016.

R. D. Madden, J. R. Sauer, and J. W. Dillwith, “A Proteomics Approach to Characterizing Tick Salivary Secretions,” Experimental and Applied Acarology, vol. 28, no. 1-4, pp. 77–87, 2002.

C. M. Maya-Monteiro, S. Daffre, C. Logullo et al., “HeLp, a heme lipoprotein from the hemolymph of the cattle tick, Boophilus microplus,” The Journal of Biological Chemistry, vol. 275, no. 47, pp. 36584–36589, 2000.

V. Díaz-Martín, R. Manzano-Román, L. Valero, A. Oleaga, A. Encinas-Grandes, and R. Pérez-Sánchez, “An insight into the proteome of the saliva of the argasid tick Ornithodoros moubata reveals important differences in saliva protein composition between the sexes,” Journal of Proteomics, vol. 80, pp. 216–235, 2013.

C. J. Oliveira, E. Anatriello, I. K. de Miranda-Santos et al., “Proteome of Rhipicephalus sanguineus tick saliva induced by the secretagogues pilocarpine and dopamine,” Ticks and Tick-Borne Diseases, vol. 4, no. 6, pp. 469–477, 2013.

J. G. Valenzuela, I. M. B. Francischetti, V. M. Pham, M. K. Garfield, T. N. Mather, and J. M. C. Ribeiro, “Exploring the sialome of the tick Ixodes scapularis,” The Journal of Experimental Biology, p. 205, 2002.

E. D. Sonenshine, Biology of Ticks, Press OU, New York, NY, USA, 1991.

D. E. Sonenshine and W. L. Hynes, “Molecular characterization and related aspects of the innate immune response in ticks,” Frontiers in Bioscience, vol. 13, no. 18, pp. 7046–7063, 2008.

N. C. Soberanes, M. V. Santamaría, H. S. Fragoso, and Z. V. García, “First case reported of amitraz resistance in the cattle tick Boophilus microplus in Mexico,” Journal of Visualized Experiments, vol. 40, pp. 81–92, 2002.

T. M. Stoepler, J. C. Castillo, J. T. Lill, and I. Eleftherianos, “A Simple Protocol for Extracting Hemocytes from Wild Caterpillars,” Journal of Visualized Experiments, vol. 15, no. 69, p. 4173, 2012, http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3523423/.

B. F. Stone and K. P. Haydock, “A method for measuring the acaricide susceptibility, of the cattle tick Boophilus microplus (Can.),” Bulletin of Entomological Research, vol. 53, no. 3, pp. 563–578, 1962.

M. M. Bradford, “Rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding,” Analytical Biochemistry, vol. 72, no. 1-2, pp. 248–254, 1976.

E. Miranda-Miranda, M. H. Murillo-Sanchez, and R. Cossio-Bayugar, “Expression of a Haemonchus contortus cysteine protease in the baculovirus system,” Electronic Journal of Biotechnology, vol. 11, pp. 49–55, 2008, http://www.scielo.cl/scielo.php?script=sci_arttext&pid=S0717-34582008000200007&nrm=iso.

R. Hood, Handbook of Detection of Enzymes on Electrophoretic Gels, CRC Press, 2nd edition, 2003.

D. E. Sonenshine and K. R. Macaluso, “Microbial invasion vs. tick immune regulation,” Frontiers in Cellular and Infection Microbiology, vol. 7, no. 390, 2017, https://www.frontiersin.org/article/10.3389/fcimb.2017.00390.

C. M. Maya-Monteiro, S. Daffre, C. Logullo et al., “HeLp, a heme lipoprotein from the hemolymph of the cattle tick, Boophilus microplus,” The Journal of Biological Chemistry, vol. 275, no. 47, pp. 36584–36589, 2000.

A. Sorkin, “Internalization of the epidermal growth factor receptor: Role in signalling,” Biochemical Society Transactions, vol. 29, no. 4, pp. 480–484, 2001.

P. Enkhbayar, M. Kamiya, M. Osaki, T. Matsumoto, and N. Matsushima, “Structural Principles of Leucine-Rich Repeat (LRR) Proteins,” Proteins: Structure, Function, and Genetics, vol. 54, no. 3, pp. 394–403, 2004.

B. Kobe and A. V. Kajava, “The leucine-rich repeat as a protein recognition motif,” Current Opinion in Structural Biology, vol. 11, no. 6, pp. 725–732, 2001.

N. J. Gay, L. C. Packman, M. A. Weldon, and J. C. J. Barna, “A leucine-rich repeat peptide derived from the Drosophila Toll receptor forms extended filaments with a β-sheet structure,” FEBS Letters, vol. 291, no. 1, pp. 87–91, 1991.

J. M. Rothberg, J. R. Jacobs, C. S. Goodman, and S. Artavanis-Tsakonas, “slit: An extracellular protein necessary for development of midline glia and commissural axon pathways contains both EGF and LRR domains,” Genes & Development, vol. 4, no. 12 A, pp. 2169–2187, 1990.

F. Spadaro, S. Cecchetti, and L. Fantuzzi, “Macrophages and phospholipases at the intersection between inflammation and the pathogenesis of hiv-1 infection,” International Journal of Molecular Sciences, vol. 18, no. 7, p. 1390, 2017.

A. G. Greenburg and H. W. Kim, “Hemoglobin-based oxygen carriers,” Critical Care (London, England), vol. 14, suppl 2, no. 8, pp. S61–S64, 2004, http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3226150/.

K. Umeda, T. Heike, M. Nakata-Hizume et al., “Sequential analysis of α- and β-globin gene expression during erythropoietic differentiation from primate embryonic stem cells,” Stem Cells, vol. 24, no. 12, pp. 2627–2636, 2006.

T. Egawa and S.-R. Yeh, “Structural and functional properties of hemoglobins from unicellular organisms as revealed by resonance Raman spectroscopy,” Journal of Inorganic Biochemistry, vol. 99, no. 1, pp. 72–96, 2005.

H. Rubin, “Serine protease inhibitors (serpins): Where mechanism meets medicine,” Nature Medicine, vol. 2, no. 6, pp. 632-633, 1996.

A. S. Tanaka, R. Andreotti, A. Gomes, R. J. S. Torquato, M. U. Sampaio, and C. A. M. Sampaio, “A double headed serine proteinase inhibitor - Human plasma kallikrein and elastase inhibitor - From Boophilus microplus larvae,” International Journal of Immunopharmacology, vol. 45, no. 1-3, pp. 171–177, 1999.

L. Tirloni, J. Reck, R. M. S. Terra et al., “Proteomic analysis of cattle tick Rhipicephalus (Boophilus) microplus saliva: A comparison between partially and fully engorged females,” PLoS ONE, vol. 9, no. 4, Article ID e94831, 2014.

K. M. U. Janadaree Bandara and S. H. P. Parakrama Karunaratne, “Mechanisms of acaricide resistance in the cattle tick Rhipicephalus (Boophilus) microplus in Sri Lanka,” Pesticide Biochemistry and Physiology, vol. 139, pp. 68–72, 2017.

T. A. Miller, “Mechanisms of resistance to pyrethroid insecticides,” Parasitology Today, vol. 4, no. 7, pp. S8–S12, 1988.

R. L. Metcalf, “Insect resistance to insecticides,” Journal of Pesticide Science, vol. 26, no. 4, pp. 333–358, 1989.

H. He, A. C. Chen, R. B. Davey, and G. W. Ivie, “Molecular cloning and nucleotide sequence of a new P450 gene, CYP319A1, from the cattle tick, Boophilus microplus,” Insect Biochemistry and Molecular Biology, vol. 32, no. 3, pp. 303–309, 2002.

R. Cossío-Bayúgar, E. Miranda-Miranda, A. Ortìz-Nàjera, S. Neri-Orantes, and S. Olvera-Valencia, “Cytochrome P-450 monooxygenase gene expression supports a multifactorial origin fro acaricide resistence in Rhipicephalus microplus,” Research Journal of Parasitology, vol. 3, no. 2, pp. 59–66, 2008, https://scialert.net/abstract/?doi=jp.2008.59.66.

N. Aïzoun, R. Aïkpon, G. Padonou et al., “Mixed-function oxidases and esterases associated with permethrin, deltamethrin and bendiocarb resistance in Anopheles gambiae s.l. in the south-north transect Benin, West Africa,” Parasites & Vectors, vol. 6, no. 1, p. 223, 2013.

J. David, H. M. Ismail, A. Chandor-Proust, and M. J. Paine, “Role of cytochrome P450s in insecticide resistance: impact on the control of mosquito-borne diseases and use of insecticides on Earth,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 368, no. 1612, pp. 20120429–20120429, 2013.

K. V. Donohue, S. M. Khalil, R. D. Mitchell, D. E. Sonenshine, and R. Michael Roe, “Molecular characterization of the major hemelipoglycoprotein in ixodid ticks,” Insect Molecular Biology, vol. 17, no. 3, pp. 197–208, 2008.

D. Tilman, “Causes, consequences and ethics of biodiversity,” Nature, vol. 405, no. 6783, pp. 208–211, 2000.

S. C. Stearns, “Trade-offs in life-history evolution,” Functional Ecology, vol. 3, no. 3, pp. 259–268, 1989.

D. A. Roff and D. J. Fairbairn, “The evolution of trade-offs: where are we?” Journal of Evolutionary Biology, vol. 20, no. 2, pp. 433–447, 2007.