Skip to main content

Rhipicephalus (Boophilus) microplus aquaporin as an effective vaccine antigen to protect against cattle tick infestations

Abstract

Background

Vaccination as a control method against the cattle tick, Rhipicephalus (Boophilus) microplus has been practiced since the introduction of two products in the mid-1990s. There is a need for a vaccine that could provide effective control of R. microplus in a more consistent fashion than existing products. During our transcriptome studies of R. microplus, several gene coding regions were discovered to encode proteins with significant amino acid similarity to aquaporins.

Methods

A cDNA encoding an aquaporin from the cattle tick, Rhipicephalus microplus, was isolated from transcriptomic studies conducted on gut tissues dissected from fully engorged adult female R. microplus.

Results

Bioinformatic analysis indicates this aquaporin, designated RmAQP1, shows greatest amino acid similarity to the human aquaporin 7 family. Members of this family of water-conducting channels can also facilitate the transport of glycerol in addition to water. The efficacy of this aquaporin as an antigen against the cattle tick was explored in cattle vaccine trials conducted in Brazil. A cDNA encoding a significant portion of RmAQP1 was expressed as a recombinant protein in Pichia pastoris, purified under native conditions using a polyhistidine C-terminus tag and nickel affinity chromatography, emulsified with Montanide adjuvant, and cattle vaccinated intramuscularly. The recombinant protein provided 75% and 68% efficacy in two cattle pen trials conducted in Campo Grande, Brazil on groups of 6 one year old Holstein calves.

Conclusion

The effectiveness of this vaccine in reducing the numbers of adult female ticks shows this aquaporin antigen holds promise as an active ingredient in cattle vaccines targeted against infestations of R. microplus.

Background

The cattle tick, Rhipicephalus (Boophilus) microplus, is an obligate parasitic cattle pest that has established populations throughout the world’s tropical and subtropical regions. R. microplus is responsible for significant economic losses to cattle producers because of direct effects through blood loss and damage to hides and indirect effects through diseases it transmits such as babesiosis and anaplasmosis. For example, Grisi et al. [1] estimated Brazil’s annual losses to R. microplus parasitism approximate US$3.2 billion. Significant efforts to control this tick are undertaken in most cattle-raising countries and these efforts presently center around the use of acaricides. However, acaricide resistant populations of R. microplus have become a major problem in most of the cattle-producing countries of the world and novel cattle tick control technologies are needed to maintain efficiencies in cattle production [2]-[4].

Vaccination as a tick control method has been practiced since the introduction of two products in the mid-1990s, TickGARD [5] and Gavac©[6], that were developed using the midgut glycoprotein Bm86 as the immunoreactive antigen. TickGARD is no longer commercially available, but Gavac© continues to be used to date, primarily in North and South America. Although the results with Gavac© have been mixed, within integrated tick management systems in some geographic regions, the vaccine has proven to reduce the number of acaricidal applications per year that are required to control R. microplus at acceptable levels [7]. As an interesting sidelight to the role of Gavac© in cattle tick control, the product has been shown to provide >99% efficacy against Rhipicephalus annulatus, a second cattle tick species which is much less prevalent and invasive than R. microplus[8],[9]. Nevertheless, the need remains for a vaccine that could provide effective control of R. microplus in a more consistent fashion than Gavac©. As part of research mining the genome of R. microplus for transcripts that would produce effective antigens for a cattle tick vaccine, focused genome [10], transcriptome [11], and proteome [12],[13] studies in R. microplus have led to identification of genes and gene coding regions that encode proteins with critical functions in the tick [14],[15]).

Several of these gene coding regions were discovered to encode proteins with significant amino acid similarity to aquaporins. Aquaporins, originally called water channels, allow the regulation of water transport across the highly hydrophobic lipid bilayer of cell membranes. Members of the aquaporin family have been found in animal taxa from mammals [16] to bacteria [17] and they are very common in certain cell types, with approximately 150,000 protein copies per red blood cell [18]. The structure of the aquaporins is such that two constrictions in the protein act as filters whose selectivity for water, glycerol, urea, and other small molecules is determined by the size and charge of the constriction pore [19]. Because cattle ticks ingest large volumes of blood relative to their body size and weight, they are required to have efficient water transport mechanisms so as to concentrate the blood components for efficient digestion [20]. Thus, the tick aquaporins are critical to tick physiology and appeared a good protein to target as an anti-cattle tick vaccine candidate. The full length transcript for one of the discovered aquaporins, designated RmAQP1, was determined and a fragment of the open reading frame (ORF) was expressed and purified as a recombinant protein in Pichia pastoris. This recombinant protein was tested in cattle pen trials for efficacy as a vaccine antigen against R. microplus.

Methods

Source of tick materials

Ticks that were the source of DNA and RNA for the transcript discovery study were obtained from engorged adult female R. microplus of the f20 La Minita strain maintained at The University of Idaho Holm Research Center (Moscow, ID). The La Minita strain was originally collected in 1996 during an outbreak in Starr County, TX and propagated at the USDA-ARS Cattle Fever Tick Research Laboratory in Edinburg, TX. Tissues that were the source of RNA for gene expression study were dissected from 25 1-2 day old adult male and female ticks from the R. microplus Deutch strain f41 generation maintained at the Cattle Fever Tick Research Laboratory. The Deutch strain is an organophosphate and pyrethroid susceptible strain originating from an outbreak in 2001 in Webb County, Texas. Tick larvae used in this study to infest cattle for the cattle vaccine trials were obtained from a laboratory colony maintained at EMBRAPA Beef Cattle. The colony was established with R. microplus ticks collected from infested cattle in Campo Grande, MS, Brazil. Larvae used for infesting cattle in the vaccine trials were 18 days post-hatch. During the vaccine trials, fully engorged adult female ticks were collected upon host detachment and brought to the laboratory to allow oviposition. Egg masses were incubated in humidity chambers at 28°C and 95% relative humidity to facilitate hatching. Larvae were used for infestation at 18 days after hatching.

RNA purification, cDNA synthesis & RACE

Total RNA was isolated using the FastPrep-24 Tissue and Cell Homogenizer and Lysing Matrix D (Qbiogene, Irvine, CA, USA) as described in Saldivar et al. [21] from gut tissue dissected from 5 engorged adult female R. microplus from the La Minita strain. The total RNA was DNAse treated following manufacturer’s protocol using Turbo DNA-free kit (Ambion, Austin, TX, USA).

One microgram DNase-free total RNA was used to make 5′ and 3′ cDNA using the SMART RACE cDNA Amplification Kit (Clontech Laboratories Inc., Mountain View, CA, USA) and Superscript III Reverse Transcriptase (Life Technologies, Grand Island, NY, USA). Primers (Sigma-Aldrich, The Woodlands, TX, USA) were designed from the available aquaporin-like sequence obtained from GenBank Accession No. CV443183 [11], synthesized, and used to obtain the full length coding region sequence. The target 5′ end was obtained using 20 μL PCR reactions containing 5′ SMART RACE cDNA, Advantage® 2 PCR Enzyme System (Clontech Laboratories Inc.), Universal Primer A and gene specific primer KB-126 (Table 1) according to manufacture’s protocol using the touchdown cycling profile suggested in the SMART RACE protocol booklet. This amplification product was used for a nested PCR using Nested Universal Primer A and gene specific primer KB-126 following manufacturer’s protocol using a cycling profile of thirty cycles including a denaturing step of 94°C for 30 sec and an annealing/extension step of 72°C for 3 min 30 sec. PCR products were analyzed by electrophoresis on 2% SeaKem Gold agarose gels in 1XTBE running buffer (Lonza Rockland, Inc., Rockland, ME, USA) and post-stained using GelStar® Nucleic Acid Gel Stain (Lonza Rockland, Inc.). The detected 965 bp amplicon was excised from the agarose gel and DNA extracted and purified using the QIAquick Gel Extraction Kit (Qiagen, Valencia, CA, USA) according to manufacturer’s protocol. The DNA was concentrated using Pellet Paint Co-Precipitant (Novagen/EMD Chemicals Inc., Gibbstown, NJ, USA), polished, ligated and transformed into XL10 Gold Kan Ultracompetent Escherichia coli cells using the PCR Script Amp Cloning Kit (Stratagene/Agilent Technologies Inc., Santa Clara, CA, USA). Individual clones were screened via PCR using internal vector primers and clones producing a correct sized product were used for plasmid DNA preparations with the QIAprep Spin Miniprep Kit (Qiagen) according to manufacturer’s instructions. Plasmid DNAs were sequenced on a 3130xl Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Identities of each nucleotide in RmAQP1 were verified on both strands to produce a high quality sequence. Sequences were assembled and analyzed using MacVector with Assembler version 10.0.2 (MacVector Inc., Cary, NC, USA). Basic Local Alignment Search Tools (BLAST) programs were run using multiple BLAST programs available at http://blast.ncbi.nlm.nih.gov/Blast.cgi[22],[23]. TOPCONS ([24], http://topcons.net/index.php) was used for the prediction of transmembrane helices in proteins. Net Phos 2.0 Server ([25], http://www.cbs.dtu.dk/services/NetPhos) was used for the predictions of serine, threonine, and tyrosine phosphorylation sites.

Table 1 Primers/Probes used for RmAQP1 transcript cloning and real-time PCR

Real-time PCR gene expression study

Quantitative PCR studies were designed with the MIQE guidelines in mind [26]. Tissue dissections were performed under phosphate-buffered saline pH =7.4. The tissues dissected from the female ticks were the synganglia, salivary glands, ovaries and midgut while tissues dissected from the male ticks were the synganglia, salivary glands, testes, accessory gland and midgut. Dissected tissues were placed in RNALater (Ambion, Austin, TX, USA) according to manufacturer’s protocol. Total RNA was isolated using the ToTALLY RNA Kit (Ambion) and DNase treated using the Turbo-DNA free kit (Ambion) according to manufacturer’s protocol. The RETROscript Kit Reverse Transcription for RT-PCR (Ambion) was used to produce cDNA from each tissue using 0.1 μg of DNase-free total RNA.

TaqMan probes and primers were designed using Beacon Designer 7.0 (PREMIER BioSoft International, Palo Alto CA; Table 1) and synthesized by Sigma-Aldrich Inc. (The Woodlands, TX, USA) for RmAQP1 and the two reference genes used for normalization, R. microplus 18S rRNA gene [21] and beta-actin. Optimization PCRs were run on all three genes to determine optimal reaction conditions, PCR efficiencies, and optimal reagent concentrations. Real-time PCR reactions were carried out in clear low profile 96 well plates (BioRad, Hercules, CA, USA) with microseal film B (BioRad) using 25 μL total volume reactions, which included TaqMan Universal Master Mix No AmpErase UNG (Applied Biosystems Inc.), 250 nM TaqMan probe, tissue specific RETROscript cDNA, and 900 nM forward and reverse primers for all three genes. The cycling profile used on the BioRad CFX96 Real-Time System was 95°C for 10 min, and 50 cycles of 95°C for 15 sec, 60°C for 1 min plus plate read. All samples were run in triplicate and both no-template and no-reverse transcriptase controls were utilized to verify DNA-free status of the samples. The fluorescence emission data analysis for the relative standard curve method for quantification was done using baseline subtracted curve fit mode with CFX Manager Software v1.5 (BioRad).

Cloning into Pichia pastoris

The partial RmAQP1 ORF used for the vaccine study was amplified with the Advantage® 2 PCR Enzyme System (Clontech Laboratories Inc.) using primers KB-156 and KB-157 (Table 1). The 597 bp amplification product was purified and gel extracted as described above. The RmAQP1 DNA was prepared for ligations by restriction enzyme digestion reactions with EcoRI and NotI (Life Technologies) per manufacturer’s protocol.

The EasySelect Pichia Expression Vector (Life Technologies), pPICZ αA restriction enzyme-digested with EcoRI and NotI and purified, was ligated onto the RmAQP1 DNA using the TA Cloning Kit (Life Technologies) using the TA Cloning Kit protocol and 137 ng RmAQP1 insert, 50 ng pPICZ αA EcoRI/NotI digested vector, and 1 unit T4 DNA ligase incubated for 17 hr at 4°C. OneShot TOP10 Electrocomp cells (Life Technologies) were transformed with ligation reaction and plated on low salt LB agar (1% tryptone, 0.5% yeast extract, 0.5% sodium chloride, 1.5% agar) with 25 μg/mL ZeocinTM (Life Technologies). Resulting colonies were screened via PCR using vector primers 5′AOX1 and 3′ AOX1 and DNA isolated from positive colonies using the QIAprep Spin Miniprep Kit (Qiagen) according to manufacturer’s instructions. The sequence of both strands of putative positive clone plasmid DNA was verified by DNA sequencing, followed by analysis with MacVector with Assembler version 10.0.2.

According to the EasySelect Pichia Expression Kit protocol, a freshly prepared 80 μL aliquot of electrocompetent P. pastoris KM71H strain and 5 μg recombinant expression vector DNA linearized with SstI was used for transformations according to the manufacturer’s instructions using the Bio-Rad Gene Pulser and Pulse Controller at pulse settings of 1.5 kV, 200′ and 25μFD. Transformation mixtures were plated on YPDS (1% yeast extract, 2% peptone, 2% dextrose, 1 M sorbitol, 2% agar) plates containing 100 μg/mL ZeocinTM and incubated at 30°C for four days to allow colonies to develop.

Analysis of Pichia pastoris transformants

Direct screening of individual Pichia KM71H colonies using PCR was done by modifying the direct screening protocol from Linder et al. [27] and the EasySelect Pichia Expression Kit manual with 25 μL reactions using the 5′ and 3′ AOX1 vector primers and 0.16 μL of a 1 vol:1 vol mix of AmpliTaq DNA polymerase (5 U/μL stock; Applied Biosystems) and TaqStart antibody (1.1 μg /μl stock; Clontech). Colonies containing the expected 1,192 bp recombinant product were re-screened using a similar approach but substituting RmAQP1-specific primers.

Selected colonies were Mut phenotyped and small-scale expression experiments used to determine the optimal method and conditions for the expression of the recombinant proteins. These protocols are described in the EasySelect Pichia Expression Kit manual for 3 mL cultures grown in BMGY (1% yeast extract, 2% peptone, 100 mM potassium phosphate pH =6.0, 1.34% Yeast nitrogen base with ammonium sulfate without amino acids, 4 × 10-5% biotin, 1% glycerol) and BMMY media (BMGY but substituting 0.5% methanol for the 1% glycerol). BMMY cultures were replenished to 0.5% final methanol concentration every 24 hr. Samples were collected at various time points and centrifuged to separate the yeast cells from the culture media supernatant.

Supernatants were frozen in liquid nitrogen and stored at -80°C. Intracellular proteins were purified by a protocol similar to that described in the EasySelect Pichia Expression Kit manual. Briefly, 100 μL of breaking buffer (50 mM sodium phosphate pH7.4, 1 mM EDTA, 5% glycerol) +1X FOCUS ProteaseArrest (GBioscience, St. Louis, MO) was used per cell pellet from a 1 ml culture sample. An equal volume of 0.5 mM acid-washed glass beads was added and the sample vortexed for 30 sec and set on ice for 30 sec. A total of 8 vortex/ice cycles were used, the sample frozen at -80°C, thawed and 8 more vortex/ice cycles used before a final short centrifugation to clarify the sample. Samples were concentrated in Amicon Ultracel units (Millipore, Billerica, MA) when necessary.

Both the intracellular cell pellets and the secreted supernatant samples were analyzed by denaturing gel electrophoresis under reducing conditions using the NuPAGE® Electrophoresis System and NuPAGE® 4-12% Bis-Tris gels in the XCell SureLockTM Mini-Cell with 1X NuPAGE MOPS SDS Running Buffer (Life Technologies) according to manufacturer’s instructions. Proteins were visualized by staining with Coomassie Brilliant Blue R-250 using a modified Fairbank’s method [28]. Recombinant aquaporin was localized in the cell pellet sample with maximal expression seen after 4 days of induction growth in BMMY.

After the optimal clone and growth conditions were determined, a large scale culture of the clone producing the highest amount of recombinant aquaporin protein was grown in 25 ml BMGY media in 500 ml baffled flasks in a shaking incubator at 30°C to an OD600 = 2-6. Cells were harvested by centrifugation and resuspended in BMMY to an OD600 = 1 and returned to the incubator for 4 days to induce expression. Every 24 hr, methanol was added to a final concentration of 0.5% to maintain induction and cells were harvested 4 days post-induction. Following centrifugation, the cell pellet was frozen at -70°C until protein extraction.

Total yeast intracellular protein was extracted similarly as described above for the small-scale expression cell pellet protocol with the exception of using 50 ml Breaking Buffer with 1X Protease Arrest and 10 cycles of 30 sec vortexing followed by 30 sec on ice. The cell pellet lysates were then frozen at -70°C overnight and thawed followed by 10 vortex/ice cycles. The protein solution was clarified by centrifugation and the resulting solution concentrated using Centricon Plus-70 Centrifugal Filter Devices (Millipore).

Purification of expressed recombinant protein

Recombinant protein was purified making use of the 6X-Histidine tag supplied by the vector sequence and the ProBond Purification System (Life Technologies) using ProBondTM nickel-chelating resin under native conditions, initially according to manufacturer’s instructions. We wished to preserve the native protein structure, thus we did not use urea, SDS, or heat in the purification steps. However, the purified protein presented solubility problems upon freezing and thawing and we adapted the ProBond purification steps to utilize buffer (50 mM NaH2PO4 pH = 8.0, 300 mM NaCl, 2 mM β-mercaptoethanol, 0.4% β -D-1-thioglucopyranoside) plus 10 mM imidazole for binding, the same buffer plus 30 mM imidazole for washing, and buffer plus 300 mM imidazole for elution [29]. Eluted protein was concentrated using Amicon Ultra-15 centrifugation units (Millipore) and, following concentration, the solution was made 50% v/v glycerol and stored at -20°C. This protein solution was used to prepare the vaccine. Protein concentration was quantified by the BioRad Protein Assay Kit I with bovine plasma gamma globulin protein standards, and purity of the protein solution verified by gel electrophoresis as described above using the NuPAGE® Electrophoresis System and NuPAGE® 4-12% Bis-Tris gels in the XCell SureLockTM Mini-Cell.

Protein identity was verified by mass spectrometry analysis and Western blotting, taking advantage of the c-myc and 6X-His tag epitopes on the recombinant protein that are provided by the expression vector sequence. The WesternBreeze Chromogenic Kit and Anti-myc-HRP and Anti-His(C-term)-HRP antibodies (Life Technologies) were utilized with standard protocols provided by the supplier. The supplier-provided alkaline phosphatase-conjugated secondary antibody was utilized to enhance sensitivity. The mass spectrometry analysis was done by Protea Bioscience Group (Morgantown, WV). The recombinant protein (in 50% glycerol solution described above) was purified by 1-D acrylamide gel electrophoresis, extracted from the gel matrix, and digested with trypsin. The resulting peptides were analyzed by LC-MS/MS using an ABSciex5500 Series QTRAP for tandem MS data acquisition followed by a search for peptide matches to the expected sequence of purified antigen.

Pen trial

Controlled pen trials were conducted to evaluate the immunogenic and protective capacity of the aquaporin antigen adjuvated with Montanide ISA 61 VG (Seppic, Paris) into doses of 2 ml containing 100 μg of the recombinant protein and phosphate-buffered saline (PBS). One-year-old Holstein calves were randomly distributed into groups of six animals each. Negative controls were injected with adjuvant prepared with PBS alone. The animals were injected intramuscularly three times with two week intervals between injections. Serum samples were taken from each animal before the first immunization and weekly thereafter. Twenty-one days after the last injection the animals were challenged with 15,000 larvae of the Campo Grande tick strain. These larvae were delivered in three applications of 5,000 larvae each during one week, placed in separate vials onto the back of the animals. As engorged female ticks detached, they were collected once a day, pooled, and weighed. This sampling was initiated upon the first day detachment started and continued until tick detachment ceased, which was 19 days for Trial 1 and 16 days for Trial 2. Twenty females from each day’s collection were pooled and incubated at 29°C and 85% humidity until egg laying was complete. Eggs collected from each pool were weighed and incubated at 29°C and 85% humidity until hatching was completed to determine the hatch percentage for each pool.

Bovine serum collection and analysis

Bovine blood was sampled weekly and separated serum frozen until analyzed by ELISA. For the ELISA, sera from all animals in each group were pooled according to day of collection. Microtiter plates were coated with antigen (50 μL per well of a 1 μg antigen/ml solution in 20 mM sodium carbonate buffer, pH 9.6) and incubated overnight at 4°C. Blocking with 2% bovine serum albumin in PBS pH 7.4 containing 0.05% Tween 20 was followed by washing five times with PBS pH 7.4. The plates were incubated for 45 min at 37°C with 100 μL per well of immunized bovine serum diluted to 1:100 in PBS. After washing in PBS pH 7.4, 50 μL of rabbit anti-bovine IgG peroxidase conjugate (Sigma, St. Louis, MO, USA) diluted 1:20,000 was added and the plate incubated for 30 min at room temperature. After washing in PBS pH 7.4, 50 μL of 1.0 mM chromogenic substrate o-phenylenediamine was added and the reaction was stopped after 15 min by adding 50 μL of 0.2 M NaOH. A microplate reader was used to assess the results with absorbance set at 490 nm.

Efficacy assessment and statistics

Reductions associated with immunization relative to the unvaccinated group were determined for numbers of adult female ticks, egg production, and larval hatching. Vaccine efficacy was calculated as 100 × [1 - (NET × EWPF × H)], where NET, EWPF, and H represent the fraction of the relevant tally in the immunized group relative to that in the control group of the total number of adult female ticks, total weight of eggs per female, and % hatch of eggs, respectively.

Ethical approval

The La Minita ticks used for the transcript discovery study were reared at The University of Idaho Holm Research Center (Moscow, ID, USA) following protocols approved by the University of Idaho Institutional Animal Care and Use Committee (IACUC). The Deutch ticks used for the gene expression studies were reared at the USDA- ARS Cattle Fever Tick Research Laboratory (Edinburg, TX, USA) with protocols approved by that Laboratory’s IACUC. The cattle vaccine studies were conducted at EMBRAPA Beef Cattle (Campo Grande, MS, Brazil) under protocols approved by the EMBRAPA review board.

Results

Aquaporin-like sequences from the cattle tick

Using ESTs from BmiGI Ver 2.0 as a starting point [11], we used 5′ RACE to isolate an 1,800 bp transcript that included the entire ORF to an aquaporin-like protein (Figure 1; GenBank Accession No. KJ626366). As this is the first aquaporin from R. microplus to be described, we designated the transcript as RmAQP1. The transcript encodes a 216 aa ORF with several stop codons flanking both the N- and C-termini of the presumptive protein, increasing our confidence that we have the authentic ORF. Analysis by TOPCONS (http://topcons.net/index.php) predicted the ORF contains 6 transmembrane helical regions and two NPA motifs, all characteristics of the aquaporin family [30]. The Net Phos 2.0 Server (http://www.cbs.dtu.dk/services/NetPhos) predicted 8 serine, 3 threonine, and 1 tyrosine phosphorylation sites.

Figure 1
figure 1

Nucleotide and translated open reading frame (ORF) sequences for RmAQP1 . Nucleotide sequence numbering starts from the 5′ end of the transcript as determined by RACE. Amino acid numbering begins from the start methionine amino acid. The entire ORF is in bold text, while the portion of the ORF used as the vaccine antigen is underlined. The six predicted transmembrane helices (aa # 11-31, 43-63, 88-108, 144-164, 176-196, and 231-251) and the two NPA motifs (aa # 68-70 and 201-203) are highlighted. The predicted phosphorylation sites at 8 serine (aa # 64, 77, 207, 265, 280, 287, 295, and 297), 3 threonine (aa # 107, 123, and 223), and 1 tyrosine (aa # 229) amino acid residues are indicated in outline font. Multiple in-frame stop codons are present both upstream of the presumptive methionine start and downstream of the stop TAG codon.

We searched other cattle tick transcriptome datasets from ongoing studies in our group and we found 2 other ESTs encoding putative full-length aquaporin-like ORFs, designated RmAQP2 and RmAQP3, and 9 ESTs encoding partial aquaporin-like ORFs (Additional file 1). We produced a ClustalW multiple alignment of RmAQP1-3 with other aquaporin-like ORFs from 5 tick species (Figure 2). The regions that align with the RmAQP1 transmembrane helices 2-6 display more amino acid similarity than other aligned regions. However, the region between predicted transmembrane helices 5 and 6 has a high number of identities in the alignment (Figure 2). In fact, that region between helices 5 and 6 contains 12 of the 40 invariant amino acids that exist over the entire alignment. The phylogenetic tree (Figure 3) that was produced from this alignment showed that RmAQP1 is most similar to the aquaporin from Ixodes scapularis (GenBank Accession No. XP_002399794). RmAQP2 was most similar to the aquaporin from Dermacentor variabilis (GenBank Accession No. ABI53034), while RmAQP3 was most similar to the aquaporins from Rhipicephalus appendiculatus (GenBank Accession No. CD780384) and Rhipicephalus sanguineus (GenBank Accession No. CAR66115). In the tissue-specific gene expression study, RmAQP1 was expressed most highly in the synganglia and lowest in the gut of male ticks, while females expressed RmAQP1 most highly in the synganglia and lowest in the ovary (Table 2).

Figure 2
figure 2

Amino acid sequence alignments of RmAQP1-3 with putative aquaporins from other tick species. Alignment was by the ClustalW multiple alignment function of MacVector 12.7.5 using the Gonnet matrix with open gap penalty of 10 and extend gap penalty of 0.05. Determination of amino acid similarity was by chemical properties of amino acid side chains with DE, AGILV, NQ, FWY, RHK, ST, CM, and P comprising the groups considered as conservative substitutions. The accession numbers for the putative tick aquaporins from R. appendiculatus, R. sanguineus, I. scapularis, I. ricinus, and D. variabilis are CD780384, CAR66115, XP_002399794, CAX48964, and ABI53034, respectively. One member each from the human aquaporin families 3 (NP 004916) and 7 (NP 001161) are also included in the alignment. The RmAQP1 was used as the model for comparing other sequences, with identities indicated by colon (:) and similarities by period (.). In the summary line below the 10 aligned sequences, a colon (:) notes amino acid positions where all 10 sequences contain the identical amino acid, a period (.) indicates all 10 sequences contain identical or similar amino acids, and an asterisk (*) indicates 9 of the 10 aligned sequences have identical or similar amino acids. Gaps inserted to optimize alignments are indicated by a dash (-). The shaded portions of the RmAQP1 sequence indicate the six predicted transmembrane helical regions.

Figure 3
figure 3

Phylogenetic tree of putative tick aquaporins. A phylogenetic tree of the multiple sequence alignment of Figure 2 was determined using the MacVector 12.7.5 package. The analysis was performed using the neighbor joining method using uncorrected distance option and bootstrap mode with 1000 replicates. Tie breaking was set as systematic and gaps distributed proportionally.

Table 2 Gene expression of RmAQP1 in various tissues of R. microplus

Production of recombinant aquaporin as vaccine antigen

During the process of determining the sequence of RmAQP1, an opportunity became available for evaluating a vaccine antigen in a controlled cattle pen test. At that time, we had cloned and sequenced only about 600 bp of the transcript, encoding only 91% of what we eventually determined to be the entire ORF (Figure 1). Nevertheless, due to the time-limited nature of the cattle pen test opportunity, we could not complete the transcript cloning within the time constraints and chose to evaluate the existing aquaporin-like antigen as a recombinant protein expressed in P. pastoris. The amino acids that comprised the vaccine antigen are shown underlined in Figure 1. The sequence encoding those amino acids was cloned into the P. pastoris expression vector pPICZ αA and the resulting recombinant protein purified as described. The full amino acid sequence of the recombinant protein is shown in Additional file 1: Table S1 and has a calculated molecular weight of 33.9 kDa. The first 91 amino acids in the 317 amino acid antigen were from the pPICZ vector, as were the final 27 amino acids. The mass spectrometry analysis verified we had produced the intended protein (Figure 4A). We detected high confidence peptides FSNSTNNGLLFINTTIASIAAK and AQASVRKFPIAK from the vector-provided N-terminus and the RmAQP1-provided ORF, respectively. The Western blot analysis of the purified protein with anti-myc-HRP (Figure 4B) and anti-His(C-term)-HRP (data not shown) antibodies confirmed the presence of these moieties on the C-terminus of the antigen.

Figure 4
figure 4

Gel electrophoresis of recombinant antigen expressed in P. pastoris and analyzed by mass spectrometry and Western blot. A) Ten μg of purified protein was added to an equal volume of sample buffer (4% sodium dodecyl sulfate, 120 mM Tris pH 6.8, 0.02% bromophenol blue, 5% β-mercaptoethanol), heated for 5 minutes at 90°C, and electrophoresed on a 10% polyacrylamide gel. Following staining with Coomassie G-250 the gel was destained with water for approximately 3 d before the indicated band (arrow) was extracted and analyzed by mass spectrometry. Lane 1: protein molecular weight standards; Lane 2: empty; Lanes 3 and 4: 10 μg of purified vaccine antigen protein with calculated molecular weight of 33.9 kDa (sequence in Additional file 1: Table S1) expressed in P. pastoris; Lanes 5 and 6: 10 μg of bovine serum albumin protein standard (MW =66.4 kDa); B) Ten μg of purified protein was electrophoresed on a NuPAGE® 4-12% Bis-Tris gel and analyzed by Western blotting using standard protocols provided with the WesternBreeze Chromogenic Kit and Anti-myc-HRP antibody (Invitrogen). Lane 1: All Blue Precision Plus Protein Standards (Bio-Rad); Lane 2: Purified recombinant Aquaporin-derived vaccine antigen. The blot image was adjusted through contrast and brightness controls to enable the visualize the minor background products of approximately 60-65 kDa.

Cattle pen tests for aquaporin antigen efficacy against R. microplus

The recombinant aquaporin-derived antigen was tested with Montanide adjuvant in cattle pen tests. Cattle were vaccinated at the beginning of the test and two and four weeks after the start date. Three weeks after the final immunization, the cattle were challenged with R. microplus larvae. The results of the cattle tests are summarized in Figure 5 and detailed in Additional file 2: Table S2, while the efficacy calculation is shown in Table 3. Trial 1 was conducted from September - December 2010 and Trial 2 was conducted from March - July 2011. Bovine blood was sampled weekly from each animal and ELISA results showed that vaccination elicited a specific humoral immune response (Figure 6; Additional file 3: Table S3). The major effect of the aquaporin antigen was on the total tick count resulting from the infestation (Table 3, NET, P <0.001). The two trial groups vaccinated with the aquaporin-derived antigen produced 29% of the adult female ticks compared to the control (vaccinated with PBS + adjuvant only). The effects on egg production and egg hatch were minor (P >0.05), however the extent of the effect on female tick production resulted in an overall efficacy of 75% and 68% for the two trials (Table 3).

Figure 5
figure 5

Tick collection data from cattle stall test vaccine evaluation trials. For both Trial 1 and Trial 2, means and standard errors are reported for A) Average number of ticks collected per animal, B) Average weight of ticks collected per animal, C) Average egg mass weight per female, and D) Average % egg hatch.

Table 3 Data from cattle stall trials evaluating RmAQP1-derived protein for efficacy as anti- R. microplus vaccine antigen
Figure 6
figure 6

ELISA results from cattle vaccination trials using recombinant expressed fragment of RmAQP1. Relative readings from ELISAs are plotted against trial day number for both the cattle group vaccinated with the recombinant expressed fragment of RmAQP1 in PBS plus adjuvant (dashed line and triangles) and the control group vaccinated with PBS plus adjuvant (solid line and circles). Blood was drawn on days 1, 7, 14, 21, 31, 41, 55, 62, 69, and 114 of the test. Three arrows indicate the days of the initial vaccination and the two boosters. The day of larval infestation is noted by an arrow with asterisk.

Discussion

The efficacy of the aquaporin-derived antigen vaccine against R. microplus in our tests was substantial enough to warrant further investigation as a potential control technology against this parasite. Prior to the two pen trials described here, we had conducted a cattle pen trial using a DNA vaccine approach and an expression vector encoding the aquaporin-derived antigen described herein [31]. We obtained approximately 50% efficacy against R. microplus (data not shown) while a rBmiTI antigen had approximately 30% efficacy as reported by Andreotti et al. [32]. Additionally, during the vaccination trials reported here, we also had other antigens being evaluated for efficacy against R. microplus. For example, in Trial 1 a salivary gland antigen and in Trial 2 a Bm86-Campo Grande antigen was evaluated as a separate group in the pen tests and showed 28% and 49% efficacy, respectively (data not shown). Thus, the aquaporin-derived vaccine was shown to outperform the other vaccines in both our pen trials. The vaccine’s major impact on R. microplus was to drastically reduce the yield of adult ticks (Table 3). Effects on average detached female tick weight, average egg mass weight, and hatch were absent or minor (Additional file 2: Table S2).

An aquaporin from I. ricinus, IrAQP1 (EMBL Accession Number FN178519), was evaluated for efficacy using in vivo feeding assays following dsRNA interference [33]. In contrast to our results, the effects from the IrAQP1 tests were manifested in significant weight reduction in treated ticks, due to reduced blood ingestion. However, reductions in adult tick mortality were not seen. There are a number of reasons that might explain the differences between these aquaporin efficacy tests. There are extensive differences between these two tick species. For example, I. ricinus is a three-host tick with an extended life cycle while R. microplus parasitizes a single host with a relatively fast life cycle. Also, IrAQP1 and RmAQP1 could be members of different aquaporin families as they have different expression patterns in different tick tissues. RmAQP1 was expressed most highly in the synganglia of both males and females (Table 2), while IrAQP1 was not detected in male adult I. ricinus or the synganglia of adult female I. ricinus[33]. Ball et al. [34] characterized the RsAQP1 from R. sanguineus and found highly similar amino acid sequence and a similar tissue expression pattern as IrAQP1. They described the tick aquaporins as falling into two families based on phylogenetic analysis of the existing aquaporin sequences at the time. Our phylogenetic analysis (Figure 3) maintains the relationships between the aquaporins of R. appendiculatus, R. sanguineus, I. ricinus, and D. variabilis noted by Ball et al. [34] with two families of aquaporins noted. However, in our phylogenetic analysis the inclusion of the additional 3 aquaporins from R. microplus discovered in our studies and our use of a different I. scapularis aquaporin appears to break out an additional aquaporin family that includes RmAQP1. This is consistent with our tissue expression results, as RmAQP1 is the first reported tick aquaporin that has substantial expression in synganglion tissue. We attempted to determine a classification of the RmAQP1-3 in conjunction with the human aquaporin classifications to perhaps learn more about the aquaporins from R. microplus. Using PSI-BLAST with 4 iterations against the NCBI nr protein database entries for Homo sapiens, all 3 R. microplus aquaporins had highest sequence similarity to several HsAQP7-like transcripts, an aquaglyceroporin (data not shown). However, there was also significant similarity to other aquaporin families, including HsAQP3, HsAQP10, HsAQP9, HsAQP4. Thus, this approach did not shed much light upon the transport capabilities of the aquaporins of R. microplus.

Conclusion

We have identified 3 aquaporin-like full length ORFs from R. microplus transcriptome datasets and a large part of one of those aquaporins, RmAQP1, was discovered to be an efficacious vaccine antigen in Brazilian Holstein calves infested with larvae from the Campo Grande strain of R. microplus. Further work is underway to evaluate the general effectiveness of this vaccine in different breeds of cattle and different geographical locations.

Authors’ contributions

FDG conducted the transcriptomic experiments that discovered the AQPs, designed the vaccine antigen and expression strategies, and led the drafting of the manuscript. RA led the design and implementation of the cattle pen tests and helped draft and revise the manuscript. KGB conducted the various PCR studies, produced the recombinant antigen, and helped draft and revise the manuscript. RCC conducted the cattle pen tests and the ELISA study. RJM and KY performed the statistical analysis of the cattle pen test data and revised the manuscript. AAPL helped design the cattle pen tests and revise the manuscript. All authors approved the final version of the manuscript submitted for publication.

Additional files

References

  1. Grisi L, Leite RC, Martins JRS, Barros ATM, Andreotti R, Cancado PHD, Perez de Leon AA, Pereira JB, Villela HS: Reassessment of the potential economic impact of cattle parasites in Brazil.Rev Bras Parasitol 2014, In press.,

  2. Andreotti R, Guerrero FD, Soares MA, Barros JC, Miller RJ, Perez de Leon A:Rhipicephalus (Boophilus) microplus acaricide resistance in the Brazilian state of Mato Grosso do Sul. Rev Bras Parasitol. 2011, 20: 127-133. 10.1590/S1984-29612011000200007.

    Article  Google Scholar 

  3. Rodríguez-Vivas RI, Rivas AL, Chowell G, Fragoso SH, Rosario CR, García Z, Smith SD, Williams JJ, Schwager SJ: Spatial distribution of acaricide profiles (Boophilus microplus strains susceptible or resistant to acaricides) in southeastern Mexico. Vet Parasitol. 2007, 146: 158-169. 10.1016/j.vetpar.2007.01.016.

    Article  PubMed  Google Scholar 

  4. Guerrero FD, de Leon AA P, Rodriguez-Vivas RI, Jonsson N, Miller RJ, Andreotti R: Acaricide Research And Development, Resistance And Resistance Monitoring. Biology of Ticks Volume 2. Edited by: Sonenshine DE, Roe RM. 2014, Oxford University Press, New York, NY, 353-381. 2

    Google Scholar 

  5. Willadsen P, Bird P, Cobon GS, Hungerford J: Commercialisation of a recombinant vaccine against Boophilus microplus. Parasitol. 1995, 110: S43-S50. 10.1017/S0031182000001487.

    Article  Google Scholar 

  6. Canales M, Enriquez A, Ramos E, Cabrera D, Dandie H, Soto A, Falcon V, Rodriguez M, de la Fuente J: Large-scale production in Pichia pastoris of the recombinant vaccine Gavac against cattle tick. Vaccine. 1997, 15: 414-422. 10.1016/S0264-410X(96)00192-2.

    Article  CAS  PubMed  Google Scholar 

  7. de la Fuente J, Almazan C, de la Lastra JM P, Kocan KM, Willadsen P: A ten-year review of commercial vaccine performance for control of tick infestations on cattle. Anim Health Res Rev. 2007, 8: 23-28. 10.1017/S1466252307001193.

    Article  PubMed  Google Scholar 

  8. Miller R, Estrada-Pena A, Almazan C, Yeater KM, Messenger M, Ellis D, de Leon AA P: Exploring the use of anti-tick vaccine as tool for integrated eradication of the cattle fever tick, Rhipicephalus (Boophilus) annulatus. Vaccine. 2012, 30: 5682-5687. 10.1016/j.vaccine.2012.05.061.

    Article  CAS  PubMed  Google Scholar 

  9. Fragoso H, Hoshman Rad P, Ortiz M, Rodriguez M, Redondo M, Herrera L, de la Fuente J: Protection against Boophilus annulatus infestations in cattle vaccinated with the B. microplus Bm86-containing vaccine Gavac. Vaccine. 1998, 16: 1990-1992. 10.1016/S0264-410X(98)00116-9.

    Article  CAS  PubMed  Google Scholar 

  10. Guerrero FD, Moolhuijzen PM, Peterson DG, Bidwell S, Caler E, Appels R, Bellgard M, Nene VM, Djikeng A: Reassociation kinetics-based approach for partial genome sequencing of the cattle tick, Rhipicephalus (Boophilus) microplus. BMC Genomics. 2010, 11: 374-10.1186/1471-2164-11-374.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Guerrero FD, Miller RJ, Rousseau M-E, Sunkara S, Quackenbush J, Lee Y, Nene V: BmiGI: a database of cDNAs expressed in Boophilus microplus, the tropical/southern cattle tick. Insect Biochem Mol Biol. 2005, 35: 585-595. 10.1016/j.ibmb.2005.01.020.

    Article  CAS  PubMed  Google Scholar 

  12. Rachinsky A, Guerrero FD, Scoles GA: Differential protein expression in ovaries of uninfected and Babesia-infected southern cattle ticks, Rhipicephalus (Boophilus) microplus. Insect Biochem Mol Biol. 2007, 37: 1291-1308. 10.1016/j.ibmb.2007.08.001.

    Article  CAS  PubMed  Google Scholar 

  13. Rachinsky A, Guerrero FD, Scoles GA: Proteomic profiling of Rhipicephalus (Boophilus) microplus midgut responses to infection with Babesia bovis. Vet Parasitol. 2008, 152: 294-313. 10.1016/j.vetpar.2007.12.027.

    Article  CAS  PubMed  Google Scholar 

  14. Rodriguez-Valle M, Lew-Tabor A, Gondro C, Moolhuijzen P, Vance M, Guerrero FD, Bellgard M, Jorgensen W: 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. 2010, 11: 437-10.1186/1471-2164-11-437.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Bellgard M, Moolhuijzen PM, Guerrero FD, Schibeci D, Rodriguez-Valle M, Peterson DG, Dowd SE, Barrero R, Hunter A, Miller RJ, Lew-Tabor AE: CattleTickBase: an integrated Internet-based bioinformatics resource for Rhipicephalus (Boophilus) microplus. Int J Parasitol. 2012, 42: 161-169. 10.1016/j.ijpara.2011.11.006.

    Article  CAS  PubMed  Google Scholar 

  16. Rojek A, Praetorius J, Frokiaer J, Nielsen S, Fenton RA: A current view of the mammalian aquaglyceroporins. Ann Rev Physiol. 2008, 70: 301-327. 10.1146/annurev.physiol.70.113006.100452.

    Article  CAS  Google Scholar 

  17. Fu D, Libson A, Miercke LJ, Weitzman C, Nollert P, Krucinski J, Stroud RM: Structure of a glycerol-conducting channel and the basis for its selectivity. Science. 2000, 290: 481-486. 10.1126/science.290.5491.481.

    Article  CAS  PubMed  Google Scholar 

  18. Denker BM, Smith BL, Kuhajda FP, Agre P: Identification, purification, and partial characterization of a novel Mr 28,000 integral membrane protein from erythrocytes and renal tubules. J Biol Chem. 1988, 263: 15634-15642.

    CAS  PubMed  Google Scholar 

  19. Beitz E, Wu B, Holm LM, Schultz JE, Zeuthen T: Point mutations in the aromatic/arginine region in aquaporin 1 allow passage of urea, glycerol, ammonia, and protons. Proc Natl Acad Sci U S A. 2006, 103: 269-274. 10.1073/pnas.0507225103.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  20. Megaw MWJ: Studies on the water balance mechanism of the tick, Boophilus microplus Canestrini. Comp Biochem Physiol. 1974, 48: 115-125. 10.1016/0300-9629(74)90859-7.

    Article  CAS  Google Scholar 

  21. Saldivar L, Guerrero FD, Miller RJ, Bendele KG, Gondro C, Brayton KA: Microanalysis of acaricide-induced gene expression in the southern cattle tick, Rhipicephalus (Boophilus) microplus. Insect Mol Biol. 2008, 17: 597-606. 10.1111/j.1365-2583.2008.00831.x.

    Article  CAS  PubMed  Google Scholar 

  22. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Altschul SF, Wootton JC, Gertz EM, Agarwala R, Morgulis A, Schaffer AA, Yu YK: Protein database searches using compositionally adjusted substitution matrices. FEBS J. 2005, 272: 5101-5109. 10.1111/j.1742-4658.2005.04945.x.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Bernsel A, Viklund H, Hennerdal A, Elofsson A: TOPCONS: consensus prediction of membrane protein topology. Nuc Acids Res. 2009, 37 (Webserver issue): W465-W468. 10.1093/nar/gkp363.

    Article  CAS  Google Scholar 

  25. Blom N, Gammeltoft S, Brunak S: Sequence and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol. 1999, 294: 1351-1362. 10.1006/jmbi.1999.3310.

    Article  CAS  PubMed  Google Scholar 

  26. Bustin SA, Benes V, Garson JA, Hellemands J, Huggett J, Kubista MK, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE Guidelines: M inimum I nformation for publication of Q uantitative real-time PCR E xperiments. Clin Chem. 2009, 55: 4-10.1373/clinchem.2008.112797.

    Article  Google Scholar 

  27. Linder S, Schliwa M, Kube-Granderath E: Direct PCR screening of Pichia pastoris clones. BioTechniques. 1996, 20: 980-982.

    CAS  PubMed  Google Scholar 

  28. Wong C, Sridhara S, Bardwell JCA, Jakob U: Heating greatly speeds Coomassie Blue staining and destaining. BioTechniques. 2000, 28: 426-432.

    CAS  PubMed  Google Scholar 

  29. Karlsson M, Fotiadis D, Sjövall S, Johansson I, Hedfalk K, Engel A, Kjellbom P: Reconstitution of water channel function of an aquaporin overexpressed and purified from Pichia pastoris. FEBS Lett. 2003, 537: 68-72. 10.1016/S0014-5793(03)00082-6.

    Article  CAS  PubMed  Google Scholar 

  30. Campbell EM, Ball A, Hoppler S, Bowman AS: Invertebrate aquaporins: a review. J Comp Physiol B. 2008, 178: 935-955. 10.1007/s00360-008-0288-2.

    Article  CAS  PubMed  Google Scholar 

  31. Guerrero Jr. F, Perez de Leon AA: Vaccination of animals to elicit a protective immune response against tick infestations and tick-borne pathogen transmission. U. S. Patent 8,722,063, May 13, 2014.

  32. Andreotti R, Cunha RC, Soares MA, Guerrero FD, Leite FPL, Perez de Leon A: Protective immunity against tick infestation in cattle vaccinated with recombinant trypsin inhibitor of Rhipicephalus microplus. Vaccine. 2012, 30: 6678-6685. 10.1016/j.vaccine.2012.08.066.

    Article  CAS  PubMed  Google Scholar 

  33. Campbell EM, Burdin M, Hoppler S, Bowman AS: Role of aquaporin in the sheep tick Ixodes ricinus: assessment as a potential control target. Int J Parasitol. 2010, 40: 15-23. 10.1016/j.ijpara.2009.06.010.

    Article  CAS  PubMed  Google Scholar 

  34. Ball A, Campbell EM, Jacob J, Hoppler S, Bowman AS: Identification, functional characterization and expression patterns of a water-specific aquaporin in the brown dog tick, Rhipicephalus sanguineus. Insect Bioch Mol Biol. 2009, 39: 105-112. 10.1016/j.ibmb.2008.10.006.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Embrapa Beef Cattle, CNPq, and Fundect are gratefully acknowledged for financial support. This article reports the results of research only. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation of endorsement by the U.S. Department of Agriculture. F.D. Guerrero, K.G. Bendele and A.A. Pérez de León are funded by USDA-ARS appropriated project 6205-32000-031-00D. The U.S. Department of Agriculture is an equal opportunity provider and employer.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Felix D Guerrero.

Additional information

Competing interests

The authors declare that they have no competing interests.

Electronic supplementary material

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Guerrero, F.D., Andreotti, R., Bendele, K.G. et al. Rhipicephalus (Boophilus) microplus aquaporin as an effective vaccine antigen to protect against cattle tick infestations. Parasites Vectors 7, 475 (2014). https://doi.org/10.1186/s13071-014-0475-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-014-0475-9

Keywords