Targeted silencing of the Aquaporin 2 gene of Rhipicephalus (Boophilus) microplus reduces tick fitness
© Hussein et al. 2015
Received: 14 July 2015
Accepted: 24 November 2015
Published: 2 December 2015
Ticks are blood-feeding arthropods that can affect human and animal health both directly by blood-feeding and indirectly by transmitting pathogens. The cattle tick Rhipicephalus (Boophilus) microplus is one of the most economically important ectoparasites of bovines worldwide and it is responsible for the transmission of the protozoan Babesia bovis, the etiological agent of bovine babesiosis. Aquaporins (AQPs) are water channel proteins implicated in physiological mechanisms of osmoregulation. Members of the AQP family are critical for blood-feeding arthropods considering the extreme osmoregulatory changes that occur during their feeding. We investigated the pattern of expression of a newly identified AQP2 gene of R. microplus (RmAQP2) in different tick tissues and stages. We also examined in vivo the biological implications of silencing expression of RmAQP2 silencing during tick feeding on either uninfected or B. bovis-infected cattle.
In silico gene analyses were performed by multiple alignments of amino acid sequences and topology prediction. Levels of RmAQP2 transcripts in different tick tissues and stages were analyzed by reverse transcriptase quantitative PCR. Patterns of expression of RmAQP2 protein were investigated by immunoblots. Gene silencing was performed by RNA interference and in vivo functional analyses carried out by feeding ticks on either uninfected or B. bovis-infected cattle.
RmAQP2 transcripts were found in unfed larvae, engorged nymphs, and salivary glands and guts of partially engorged females; however, of all tick tissues and stages examined, RmAQP2 protein was found only in salivary glands of partially engorged females. RmAQP2 silencing significantly reduced tick fitness and completely abrogated protein expression. The effect of RmAQP2 silencing on fitness was more pronounced in females fed on a B. bovis-infected calf than in ticks fed on an uninfected calf and none of their larval progeny survived.
Collectively, considering the gene expression and tick fitness data, we conclude that RmAQP2 is critical for tick blood feeding and may be a suitable candidate target for the development of novel strategies to control R. microplus and tick-borne parasites.
KeywordsRhipicephalus (Boophilus) microplus Aquaporin Tick fitness RNA interference Gene silencing Babesia bovis
Aquaporins (AQPs) are a family of transmembrane proteins that form pores to transport water and small solutes across cellular membranes . Members of the AQP family have been identified throughout the plant and animal kingdoms . The structures of AQPs are highly conserved among species, consisting of six transmembrane domains that are connected by two intracellular loops and three extracellular loops. Two asparagine-proline-alanine (NPA) motifs are considered AQP signature motifs and are located at the protein portion that interacts to form a pore . A total of 13 AQP members have been identified so far and classed into two subsets: those that are permeated by water and those that are permeated by water plus other small molecules, such as glycerol and urea [1, 2].
Numerous members of the AQP family have been identified in arthropods in the last few years due to the availability of several arthropod genomes and genetic information, such as transcriptomes and cDNA libraries . Studies have shown that AQPs play a pivotal role in arthropods, especially in blood-feeding species, such as mosquitoes and ticks. In fact, blood-feeding arthropods have become a model for AQP research due to the critical relevance that osmoregulation plays during feeding . Tick females can take on up to 100 times their body weight in blood, and concentrate the blood meal by returning approximately 75 % of the ingested water and ions via their saliva into the host . Therefore, the osmoregulatory system of ticks is central for their parasitic life cycle and has implications for efficient feeding and subsequent production of viable offspring.
The cattle tick Rhipicephalus (Boophilus) microplus is one of the most economically important ectoparasites of bovines, as it is the main vector of the apicomplexan protozoan Babesia bovis, the etiological agent of bovine babesiosis . Adult females R. microplus acquire B. bovis merozoites by ingesting blood from an infected bovine and pass the protozoan transovarially to their larval progeny, which then can transmit B. bovis sporozoites to naïve cattle during subsequent feeding [6–8]. In most endemic areas control of bovine babesiosis transmission relies mainly on targeting populations of R. microplus with acaricide treatments and to a lesser extent with the use of live attenuated vaccines [6–10]. However, the efficacy of commercial anti-tick vaccines is inconsistent in different regions of the world and the recent development of tick populations resistant to acaricides represents a serious threat to the cattle industry [11, 12]. Additionally, the re-emergence of R. microplus in areas that had been considered to be free of this tick, such as the regions outside the permanent quarantine zone in south Texas, USA, is causing concerns about the re-establishment of active B. bovis transmission in areas currently free of bovine babesiosis. Exposure of naïve cattle in these areas to B. bovis would lead to significant mortality since no herd immunity is present in the population.
Intake of blood meals leading to full engorgement of adult females and subsequent production of viable larval offspring are critical steps in the tick life cycle. A better understanding of mechanisms involved in these processes may lead to the identification of novel targets to control ticks and tick-borne parasites. In the present study we investigated the pattern of expression of a newly identified AQP2 gene of R. microplus (RmAQP2). Our findings showed that RmAQP2 is transcribed in unfed larvae, engorged nymphs, and salivary glands and guts of partially engorged females. However, among the tick tissues and stages analyzed, RmAQP2 protein was only found in salivary glands of partially engorged females. In vivo analyses demonstrated that silencing of RmAQP2 with RNA interference significantly reduced tick fitness and this effect was even more dramatic in females fed on a calf during acute B. bovis infection than in ticks fed on an uninfected calf.
Cattle, ticks and protozoan
Primer purposes and sequences, gene region and PCR product size for Rhipicephalus microplus Cox1 and RmAQP2
Product size (bp)
Cox1 amplification external
1 to 1298
Cox1 amplification nested primers
21 to 1072
RmAQP2 dsRNA synthesis
17 to 414
RmAQP2 dsRNA synthesis
614 to 1009
RmAQP2 Real-time PCR
349 to 453
To obtain unfed adult ticks for each experiment, approximately 40,000 larvae from 2 g of eggs were placed under a cloth patch on uninfected calves. On day 13–14, engorged nymphs were manually removed and held in an incubator at 25 °C with 96 % RH to molt to adults. After 2–3 days of incubation, freshly molted unfed adult females and males were sorted out and used for evaluation of gene expression. The first silencing experiment was carried out by placing injected ticks on a calf infected with Babesia bovis. The second silencing experiment was done on an uninfected calf. For the experiment with the Babesia bovis infected animal, the calf was injected with a stabilate containing approximately 1.4 × 108 B. bovis-infected erythrocytes (T2Bo strain) . The infected calf was monitored daily for the presence of B. bovis in peripheral blood and clinical signs of babesiosis. Parasitemia of B. bovis in peripheral blood was examined by qPCR to amplify the single copy msa-1 gene as previously described . The B. bovis-infected animal presented clinical indications of acute B. bovis infection including a drop in PCV, fever and detection of parasites in peripheral blood by qPCR (data not shown).
In silico analysis
Nucleotide sequence encoding for RmAQP2 cDNA was obtained from the R. microplus gene index project [16, 17] [GenBank ID: KP406519]. Multiple alignments of amino acid sequences were generated using the Multiple Alignment Module of LaserGene of DNASTAR software® (DNASTAR, Inc.). The Simple Modular Architecture Research Tool (SMART)  and the Transmembrane Hidden Markov Model package 2 (TMHMM2)  were used to predict domains and signal peptides in the RmAQP2 protein sequences. The tools Transmembrane Protein Representation in 2 Dimensions (TMRPres2)  and Swiss-model  were used to create 2D and 3D (respectively) visual representations of the RmAQP2 protein for topology analyses.
Reverse transcriptase real-time PCR
A reverse transcriptase quantitative real-time PCR (RT-qPCR) was standardized to assess the level of expression of RmAQP2 in different tick tissues and stages. In order to make solid significant observations regarding gene expression, qPCR was designed using the pertinent requirements of the minimum information for publication of qPCR experiments . Unfed larvae, engorged nymphs, unfed male ticks, and salivary glands, guts and ovaries from partially engorged tick females were collected in RNAlater® (Ambion) and stored at −20 °C following the manufacturer’s protocol. Total RNA was extracted using the RNAqueous® Kit (Ambion) according to the manufacturer’s protocol and quantified by Qubit™ fluorometer. Two hundred nanograms of total RNA were utilized for cDNA synthesis using the Superscript® Vilo™ cDNA Synthesis Kit (Invitrogen) following the manufacturer’s protocol. Technical replicates were performed to evaluate enzymatic variations during the synthesis of cDNA in a given RNA sample. The RmAQP2 cDNA sequence was used to design qPCR primers to amplify a fragment of 105 base pairs (Table 1). Quantitative real-time PCR were performed in a CFX96™ Real-Time PCR Detection System using the SsoFast™ EvaGreen® Supermix (Bio-Rad). Cycling conditions consisted of an enzyme activation step of 95 °C for 30 s followed by 40 cycles of denaturation at 95 °C for 5 s and annealing/extension of 60 °C for 5 s. Reactions were performed in duplicate in 20 μl using 200 nM of each primer and 2 μl of a 1/20 dilution of cDNA as template. The CFX Manager™ Software (Bio-Rad) was used to analyze the RT-qPCR data. Gene expression was normalized to the total amount of RNA used to generate the cDNA, as xpreviously described . The RmAQP2 transcription level was then calculated as a relative expression using the formula: Relative expression(sample) = 2[Cq(control) − Cq(sample)], where the control is either the highest or lowest Cq value for the gene of interest, as previously described [13, 25]. Melt curve analyses showed the absence of primer dimers and nonspecific amplification. The absence of PCR product in no RNA control reactions and in no transcriptase reverse control reactions indicated the specificity of the RmAQP2 RT-qPCR (data not shown).
Monoclonal antibody production and immunobloting assays
Three synthetic peptides from the extracellular loops of RmAQP2 were manufactured by BioSynthesis, Inc. (Texas, USA) (Additional file 1: Figure S1). Peptide 1: LGSVGLAAAP (10-mer) (Amino acids 52 to 61). Peptide 2: ADALSQVDVNLAIVYGTNATAPVFSCFPAPGV (32-mer) (Amino acids 125 to 156). Peptide 3: MCGWGSAVFSFRSYNWFWV (19-mer) (Amino acids 229 to 247). Three 6-week-old BALB/c female mice were immunized subcutaneously with 50 mg of a pool of all three RmAQP2 peptides diluted in 0.1 ml of sterile PBS emulsified with an equal volume of TiterMax® Gold Adjuvant (Sigma-Aldrich). The primary immune response was boosted twice by subcutaneous immunization at 15-day intervals with the same concentration of peptide antigen plus adjuvant. Mice immune responses were monitored by ELISA using the RmAQP2 peptides as antigens. Three days prior to cell fusion, mice were boosted intravenously with the same concentration of antigen without adjuvant. Cell fusions and cloning by limiting dilution were performed by standard procedures . Supernatants from the initial fusion and from clones obtained by limiting dilution were screened by ELISA. A total of 10 hybridoma clones were obtained and the hybridoma supernatants used to assess the expression of RmAQP2 in tick tissues and stages by immunoblot. For the immunoblot, antigens were prepared from unfed larvae, engorged nymphs, and salivary glands and guts of partially engorged females (at 5 days of feeding). Tick tissues and stages were homogenized in lysis buffer-Nonidet-P40 (NP-40) (150 mM sodium chloride 1.0 % NP-40 - 50 mM Tris, pH 8.0) and protease inhibitor (1 μg/ml). Total protein was quantified by Micro BCA Protein Assay (Thermo Scientific Pierce) and 3 μg of total protein were separated on 4–20 % Mini-PROTEAN® TGX™ Precast Gels (Bio-Rad) under reducing conditions. Proteins were transferred to a nitrocellulose membrane (Whatman, Dassel, Germany) for 1 h at 100 V. The membrane was blocked with 5 % skim milk in TBS (Tris-buffered saline: 25 mmol l−1 Tris–HCl, 150 mmol l−1 NaCl, pH7.6) for 1 h at room temperature, washed three times in TBS and incubated 1 h with the appropriate hybridoma supernatant as the primary antibody. The membrane was washed three times with TBS and incubated 30 min with anti-mouse HRP conjugated secondary antibody, then washed again three times with TBS and developed using chemiluminescent HRP antibody detection reagents.
Synthesis of double stranded RNA
Synthesis of double stranded RNA (dsRNA) was performed as previously described . Specific primers were designed based on the cDNA sequence of the RmAQP2 gene described in the R. microplus Gene Index Project  (Table 1). Two fragments of approximately 400 bp were amplified, one located at the 5′ and another one located at the 3′ of the RmAQP2 gene (Additional file 1: Figure S1). PCR products were cloned into pCR™II-TOPO® (Invitrogen), sequenced and used for in vitro transcription. The MEGAscript® Transcription Kit (Ambion) was used for the dsRNA synthesis following the manufacturer’s protocol. The two RmAQP2 dsRNA molecules were checked by electrophoresis on agarose gel, quantified by spectrophotometry and kept at −20 °C until used for tick injection.
Injection of ticks with double stranded RNA and assessment of gene silencing
Freshly molted unfed females were used for the dsRNA injection. Two independent experiments were run, each with a control group of 200 ticks injected with buffer alone and an experimental group with 200 female ticks each injected with either one or two dsRNA segments identical to RmAQP2 (Additional file 1: Figure S1). For the first experiment, adult female ticks were injected with one segment of dsRNA from the 3′ end of the gene and fed on B. bovis-infected calf. For the second experiment, adult female ticks were injected with both segments of dsRNA, one from from the 5′ end and one from the 3′ end of the gene, then fed to repletion on an uninfected calf (B. bovis-free animal). Individual females were injected with 1 μl total dsRNA (approximately 1 × 1011 molecules dissolved in 0.1 mM EDTA buffer) or buffer control (0.1 mM EDTA buffer) through the coxal membrane at the base of the 4th leg on the right ventral side, as previously described . Ticks injected with the two different dsRNA segments received 0.5 × 1011 of each molecule for a total of 1 × 1011 molecules. Injections were accomplished using a 10 μl syringe with a 33 gauge needle (Hamilton, Bonaduz, Switzerland) and the microprocessor controlled UMP3 injection pump apparatus (World Precision Instruments, Berlin, Germany). After the injection, the dsRNA-injected females, plus an equal number of males, and the buffer-injected females, plus an equal number of males, were placed under separate stockinet sleeves glued to the calf to feed on either a B. bovis-infected calf (experiment 1) or an uninfected calf (experiment 2). For both experiments silencing level of RmAQP2 was investigated using RT-qPCR as described above. Five days after the dsRNA injection, 20 partially engorged female ticks from each group were collected for dissection. Tissue collection, extraction of total RNA, cDNA Synthesis, and qPCR were performed as described above.
Evaluation of tick fitness
After the dsRNA injection, individual stockinet sleeves were checked daily for the presence of engorged female ticks. Fully engorged females were collected, weighed and put in individual wells in 24-well plates at 26 °C for oviposition. At day 14 after the beginning of oviposition, egg masses laid by each individual female were weighed and put in individual vials to evaluate hatching. Hatching was evaluated at 30 days after the egg masses were weighed and hatching positive was defined as the presence of any larvae from eggs of an individual female. The larval progeny were maintained in individual vials at 26 °C for 45 days and the larval survival was determined as the presence of any live larvae among larval progeny from each individual female.
Weights of engorged females and egg masses were compared by Student’s t-test (GraphPad Instat®, version 3.06, GraphPad Software, Inc., San Diego, CA, USA). The percentages of engorged females, oviposition, hatching and larvae survival were compared by Chi-squared (GraphPad Instat®). Gene expression in different tick tissues and stages was compared by one-way ANOVA and Tukey post-hoc test.
In silico sequence analysis of the R. microplus AQP2
Differential pattern of transcription of the RmAQP2 gene among distinct tissues and live stages of the tick
Differential Expression of the RmAQP2 protein in different tick stages and tissues
Silencing of RmAQP2 in salivary glands and its effect on protein expression
Effect of RmAQP2 silencing on tick fitness
Biological effect of RmAQP2 silencing on tick fed on a calf during acute B. bovis infection (experiment 1)
Percentage of engorged females
Weight (mg) of engorged females
Egg mass (mg)
Percentage of hatching
Percentage of larvae survival
Control females fed on a B. bovis-infected calf
46.1 % (83/180)
98.7 % (82/83)
100 % (74/74)
Silenced females fed on a B. bovis-infected calf
46.1 % (83/180)
96.3 % (80/83)
Biological effect of RmAQP2 silencing on tick fed on an uninfected calf (experiment 2)
Percentage of engorged females
Weight (mg) of engorged females
Egg mass (mg)
Percentage of hatching
Percentage of larvae survival
Control females fed on an uninfected calf
37.7 % (68/180)
98.5 % (67/68)
Silenced females fed on an uninfected calf
51.1 % (92/180)
94.5 % (87/92)
In the present study we investigated the pattern of expression of a newly identified aquaporin 2 of R. microplus ticks. RmAQP2 transcript was demonstrated in unfed larvae, engorged nymphs, and salivary glands and guts of partially engorged female ticks. Interestingly, however, among the ticks tissues and stages investigated, RmAQP2 protein was found only differentially expressed in the salivary glands of partially engorged females. In vivo functional analyses showed that silencing of RmAQP2 decreased tick fitness and this effect was more pronounced in ticks fed on a calf during acute B. bovis infection. Noteworthy, none of the larval progeny from the RmAQP2-silenced females that were fed on a B. bovis-infected calf survived. Considering that R. microplus larvae transmit B. bovis, these data may have implications for the control of tick-borne parasites .
Since R. microplus is a one-host tick that infests the host as larvae and remains until the adult females have fed to repletion, we did not think it was meaningful to assess all tick life stages of the tick for aquaporin expression. The adult female takes the majority of the blood. It is that stage that is primarily responsible for acquisition of B. bovis and it is also that stage that causes the most direct damage to the host by feeding. B. bovis is passed from the adult female transovarially to the larval offspring and it is the larval stage that transmits B. bovis. Consequently, the most biologically relevant stages of the tick for acquisition and transmission of B. bovis are the adult female and the larvae, respectively. Adult females are the most biologically relevant stages with regard to damage cause by tick feeding. Therefore, our investigations focused primarily on these two stages for expression of the aquaporin gene.
Members of the AQP family have recently been identified in several arthropods including tick species, such as D. variabilis, R. sanguineus, I. ricinus, I. scapularis and I. scapularis [28, 29]. Here we described a newly identified AQP2 of R. microplus ticks. The presence of two NPA domains and the level of amino acid identity between RmAQP2 and other AQPs prompt us to argue about a potential evolutionary conserved function for these proteins. The remarkable similarity of topology among AQP members in different species, including RmAQP2 described in this study, also supports the prediction of a putative conserved function [1, 2, 4]. A cysteine located five residues upstream of the second NPA motif suggests that RmAQP2 is likely to be sensitive to mercury as demonstrated for other AQP members . This may set the rationale for the design of novel acaricides targeting AQP to control R. microplus infestation and R. microplus-borne parasites. It has been proposed that AQPs play a critical osmoregulatory function in blood-feeding arthropods during saliva production, intake and excretion of water during blood digestion, reproduction and egg development, and off-host environment stress tolerance . Salivary glands are of vital importance for tick biology, acting as osmoregulatory organs for both off-host and on-host stages of their life cycle . We postulate that AQPs expressed uniquely in the tick salivary glands may play an important role in the osmoregulatory stress during blood feeding and digestion. RmAQP2-silenced female ticks had a significant increase in body weight at the time of repletion compared to control non-silenced ticks, suggesting a level of involvement of RmAQP2 in osmoregulatory mechanisms. Additionally, the increased body weight implies that ticks accumulate more liquids in the absence of RmAQP2. Despite being heavier, the RmAQP2-silenced group produced egg masses of reduced weight and reduced hatching success when compared with the control group, suggesting a negative correlation between intake of blood and generation of viable larval offspring.
Considering that RmAQP2 protein was only found in salivary glands, it is reasonable to assume that this AQP plays an essential role in osmoregulation during blood feeding and digestion. The data showed that knockdown of RmAQP2 interfered with the osmoregulatoty system of R. microplus during feeding as demonstrated by the increased body weight in the gene-silenced engorged females. Effect on tick fitness was even more pronounced in ticks fed on a B. bovis-infected calf perhaps due to the deleterious presence of the apicomplexa parasite infection . Additionally, it is reasonable to speculate about a potential cumulative effect of RmAQP2 silencing throughout the tick life cycle considering that none of the larval progeny survived from the silenced females that fed on the calf during acute B. bovis infection. Here we also show that dsRNA injection abolished the expression of RmAQP2 protein completely in salivary gland of partially engorged females at 5 days of feeding. This result suggests a possible rapid turnover of RmAQP2 in salivary glands. Rapid turnover of AQPs has been previously demonstrated both in animal and plant tissues and shown to be more pronounced during infections or when organisms are under stress [30–32]. The extreme osmoregulatory changes in salivary glands of ticks during feeding could induce a high level of physiological turnover of AQPs, as demonstrated by the present data for RmAQP2. It was recently demonstrated that AQP1 mRNA of I. ricinus ticks is expressed only in female tissues involved in water flux, such as guts, rectal sac and especially abundant in salivary glands . Here we demonstrated a similar pattern of expression for RmAQP2 transcripts. However, RmAQP2 protein was found in salivary glands. Therefore, it is plausible to assume the presence of a mechanism to suppress translation of RmAQP2 in female guts, larvae and nymphs. Silencing of I. ricinus AQP1 affected feeding performance, but the ticks remained feeding on the host with subsequent potential for pathogen transmission. Therefore, the authors concluded that I. ricinus AQP1 is not a suitable target to control tick infestation or block transmission of tick-borne diseases . Our data showed that RmAQP2-silenced female ticks ingested more blood and consequently showed increased engorgement body weight. Hatching and larvae survival were also affected by gene silencing and most importantly, none of the larval progeny from the RmAQP2-silenced females fed on a B. bovis-infected calf survived. Differences in site and level of expression, and permeability may explain the discrepancies between the present study and the previous published data regarding I. ricinus AQP1 . Another recent study evaluated AQP1 of R. microplus as an anti-tick vaccine . The authors demonstrated that vaccination of cattle with recombinant R. microplus AQP1 effectively reduced the number of adult female ticks that fed to repletion.
Based on the amino acid sequence, RmAQP2 has a predicted molecular weight of ≈ 30 KDa. However, preliminary in silico analysis using NetOGlyc  and NetPhos  revealed the presence of several putative glycosylation and phosphorylation sites on RmAQP2 (data not shown). Post-translational modifications of AQP have been previously described [35–37] and it may be one explanation for the detection of larger protein products in the immunoblot assays. Structural predictions suggest that this protein may exist in the membrane as a tetramer and the multiple banding pattern in the blot could represent varying stages of the breakdown of this tetrameric structure; cleavage of the intracellular tails of the RmAQP2 would result in a molecular mass of ≈ 25 KDa, suggesting that the larger 50 KDa band could be a dimer of two molecules. In additional blots (not shown) we also see another band above 100 KDa, which could represent the tetramer. Another possibility to explain the difference between expected and observed molecular weight is that the RmAQP2 cDNA sequence obtained from the R. microplus gene index project may be incomplete at the 5′ end. This sequence was derived from an EST project and although it has two stop codons in the 3′ end of the sequence, a stop codon in-frame with the presumed start ATG codon is not present. Nevertheless, this possibility is unlikely considering that the topology analyses revealed a typical predicted conformation for the AQP family [3, 4, 28]. Regardless of these considerations, the key finding is that all of the bands recognized by the monoclonal antibody raised to the RmAQP2 peptides are absent in the silenced ticks, suggesting that silencing has completely abrogated protein expression. It was beyond the scope of this study to determine the type of cell(s) in the tick salivary glands expressing the RmAQP2 protein; however, the immunoblot data clearly demonstrate that salivary glands are the unique site of expression of this AQP in R. microplus.
RNA interference has been widely used to investigate gene function in arthropods. However, off-target effects have been described in some species  and cannot be entirely ruled out as the cause of the results observed in this experiment. The R. microplus genome sequence is not available and the alignment analyses are restricted to sequences listed in databases. Blast analysis of the dsRNA sequences used in this study did not reveal significant homology to any known tick sequence other than RmAQP2. Additionally, gene silencing was validated by RT-qPCR and RmAQP2 was silenced at similar levels when we used either one or two different segments of dsRNA identical to the target gene. Taken together, these aspects support the data and represent the best possible strategy to make solid scientific observations regarding the biological effect of RmAQP2 silencing in R. microplus.
Here we demonstrate the pattern of expression of a newly identified aquaporin 2 gene of R. microplus ticks. RmAQP2 transcripts were present in unfed larvae, engorged nymphs, and salivary glands and guts of partially engorged female ticks. However, among the distinct adult tissues and life stages tested, RmAQP2 protein was only found in salivary glands of partially engorged females. In vivo investigation demonstrated that silencing of RmAQP2 decreased tick fitness. The effect of RmAQP2 silencing on tick fitness was even more pronounced in females fed on a B. bovis-infected calf and none of their larval progeny survived. Therefore, we conclude that RmAQP2 is involved in osmoregulation during feeding and is a potential target for the development of novel strategies for the control of R. microplus and R. microplus-borne diseases.
We thank Kathy Mason, Ralph Horn and James Allison for excellent technical assistance. This work was supported by USDA-ARS-CRIS project number 5348-32000-034-00D. H. E. Hussein was financially supported by the Egyptian Ministry of High Education.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Borgnia M, Nielsen S, Engel A, Agre P. Cellular and molecular biology of the aquaporin water channels. Annu Rev Biochem. 1999;68:425–58.View ArticlePubMedGoogle Scholar
- Gonen T, Walz T. The structure of aquaporins. Q Rev Biophys. 2006;39:361–96.View ArticlePubMedGoogle Scholar
- Benoit JB, Hansen IA, Szuter EM, Drake LL, Burnett DL, Attardo GM. Emerging roles of aquaporins in relation to the physiology of blood-feeding arthropods. J Comp Physiol B. 2014;184(7):811–25.View ArticlePubMedGoogle Scholar
- Campbell EM, Ball A, Hoppler S. Invertebrate aquaporins: a review. J Comp Physiol B. 2008;178:935–55.View ArticlePubMedGoogle Scholar
- Kaufman WR, Phillips JE. Ion and water balance in the ixodid tick Dermacentor andersoni. I. routes of ion and water excretion. J Exp Biol. 1973;58:523–36.Google Scholar
- Friedhoff KT. Transmission of babesia. In: Ristic M, editor. Babesiosis of Domestic Animals and Man. Boca Raton, FL, USA: CRC Press; 1988. p. 23–52.Google Scholar
- Mahoney DF, Mirre GB. A note on the transmission of Babesia bovis (syn B Argentina) by the one-host tick, Boophilus microplus. Res Vet Sci. 1979;26:253–4.PubMedGoogle Scholar
- Bock R, Jackson L, de Vos A, Jorgensen W. Babesiosis of cattle. Parasitology. 2004;129(Suppl):S247–69.View ArticlePubMedGoogle Scholar
- Fragoso H, Rad PH, Ortiz M, Rodriguez M, Redondo M, Herrera L, et al. Protection against Boophilus annulatus infestations in cattle vaccinated with the B. microplus Bm86-containing vaccine Gavac. Vaccine. 1998;16:1990–2.View ArticlePubMedGoogle Scholar
- Jonsson NN, Matschoss AL, Pepper P, Green PE, Albrecht MS, Hungerford JJ, et al. Evaluation of TickGARDPLUS, a novel vaccine against Boophilus microplus, in lactating Holstein- Friesian cows. Vet Parasitol. 2000;88:275–85.View ArticlePubMedGoogle Scholar
- Miller RJ, Davey RB, George JE. First report of organophosphateresistant Boophilus microplus (Acari: Ixodidae) within the United States. J Med Entomol. 2005;42:912–7.View ArticlePubMedGoogle Scholar
- de la Fuente J, Almazán C, Canales M, de la Lastra JMP, 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–8.View ArticlePubMedGoogle Scholar
- Bastos RG, Ueti MW, Guerrero FD, Knowles DP, Scoles GA. Silencing of a putative immunophilin gene in the cattle tick Rhipicephalus (Boophilus) microplus increases the infection rate of Babesia bovis in larval progeny. Parasites Vectors. 2009;2:57.PubMed CentralView ArticlePubMedGoogle Scholar
- Goff WL, Molloy JB, Johnson WC, Suarez CE, Pino I, Rhalem A, et al. Validation of a competitive enzyme-linked immunosorbent assay for detection of antibodies against Babesia bovis. Clin Vaccine Immunol. 2006;13:1212–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Stiller D, Goff WL, Johnson LW, Knowles DP. Dermacentor variabilis and Boophilus microplus (Acari: Ixodidae): experimental vectors of Babesia equi to equids. J Med Entomol. 2002;39:667–70.View ArticlePubMedGoogle Scholar
- Guerrero FD, Andreotti R, Bandele KG, Cunha RC, Miller RJ, Yeater K, et al. Rhipicephalus (Boophilus) microplus aquaporin as an effective vaccine antigen to protect against cattle tick infestations. Parasites Vectors. 2014;7:475.PubMed CentralPubMedGoogle Scholar
- Guerrero FD, Miller RJ, Rousseau ME, Sunkara S, Quackenbush J, Lee Y, et al. BmiGI: A database of cDNAs expressed in Boophilus microplus, the tropical/southern cattle tick. Insect Biochem Mol Biol. 2005;35:585–95.View ArticlePubMedGoogle Scholar
- McCooke JK, Guerrero FD, Barrero RA, Black M, Hunter A, Bell C, et al. The mitochondrial genome of a Texas outbreak strain of the cattle tick, Rhipicephalus (Boophilus) microplus, derived from whole genome sequencing Pacific Biosciences and Illumina reads. Gene. 2015;571:135–41.View ArticlePubMedGoogle Scholar
- Goff WL, Johnson WC, Cluff CW. Babesia bovis immunity. In vitro and in vivo evidence for IL-10 regulation of IFN-gamma and iNOS. Ann N Y Acad Sci. 1998;849:161–80.View ArticlePubMedGoogle Scholar
- Letunic I, Doerks T, Bork P. SMART: recent updates, new developments and status in 2015. Nucleic Acids Res. 2015;43:Database issue D257–D260.View ArticlePubMedGoogle Scholar
- Chaturvedi N, Shanker S, Singh VK, Sinha D. Hidden markov model for the prediction of transmembrane proteins using MATLAB. Bioinformation. 2011;7:418–21.PubMed CentralView ArticlePubMedGoogle Scholar
- Spyropoulos IC, Liakopoulos TD, Bagos PG, Hamodrakas SJ. TMRPres2D: high quality visual representation of transmembrane protein models. Bioinformatics. 2004;20:3258–60.View ArticlePubMedGoogle Scholar
- Biasini M, Bienert S, Waterhouse A, Arnold K, Studer G, Schmidt T, et al. SWISS-MODEL: modelling protein tertiary and quaternary structure using evolutionary information. Nucleic Acids Res. 2014;42(W1):W252–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55(4):611–22.View ArticlePubMedGoogle Scholar
- Bastos RG, Suarez CE, Laughery JM, Johnson WC, Ueti MW, Knowles DP. Differential expression if three members of the multidomain CCp family in Babesia bigemina. Babesia bovis and Theileria equi. PLoS ONE. 2013;8:e67765.PubMed CentralView ArticlePubMedGoogle Scholar
- Yokoyama WM. Production of monoclonal antibodies. In: Coligan JE, editor. Current protocols in immunology. 1st ed. New York: Wiley Intersciences Inc; 1994. p. 2.5.1–2.5.17.Google Scholar
- Laughery JM, Knowles DP, Schneider DA, Bastos RG, McElwain TF, Suarez CE. Targeted Surface Expression of an Exogenous Antigen in Stably Transfected Babesia bovis. PLoS ONE. 2014;9:e97890.PubMed CentralView ArticlePubMedGoogle Scholar
- Ball A, Campbell EM, Jacob J, Hoppler S, Bowmana AS. Identification, functional characterization and expression patterns of a water-specific aquaporin in the brown dog tick, Rhipicephalus sanguineus. Insect Biochem Mol Biol. 2009;39:105–12.View ArticlePubMedGoogle Scholar
- Campbell EM, Burdin M, Hoppler S, Bowmana AS. Role of an aquaporin in the sheep tick Ixodes ricinus: Assessment as a potential control target. Int J Parasitol. 2010;40:15–23.View ArticlePubMedGoogle Scholar
- Fleurat-Lessard P, Michonneau P, Maeshima M, Drevon JJ, Serral R. The distribution of aquaporin subtypes (PIP1, PIP2 and g-TIP) is tissue dependent in soybean (Glycine max) root nodules. Ann Bot. 2005;96:457–60.PubMed CentralView ArticlePubMedGoogle Scholar
- Castañeyra-Ruiz L, González-Marrero I, González-Toledo JM, Castañeyra-Ruiz A, Paz-Carmona H, Castañeyra-Perdomo A, et al. Aquaporin-4 expression in the cerebrospinal fluid in congenital human hydrocephalus. Fluids Barriers CNS. 2013;10:18.PubMed CentralView ArticlePubMedGoogle Scholar
- De Bellis M, Pisani F, Mola MG, Basco D, Catalano F, Nicchia GP, et al. A novel human aquaporin-4 splice variant exhibits a dominant-negative activity: a new mechanism to regulate water permeability. Mol Biol Cell. 2014;25:470–80.PubMed CentralView ArticlePubMedGoogle Scholar
- Steentoft C, Vakhrushev SY, Joshi HJ, Kong Y, Vester-Christensen MB, Schjoldager KT, et al. Precision mapping of the human O-GalNAc glycoproteome through SimpleCell technology. EMBO J. 2013;32(10):1478–88.PubMed CentralView ArticlePubMedGoogle Scholar
- Blom N, Gammeltoft S, Brunak S. Sequence- and structure-based prediction of eukaryotic protein phosphorylation sites. J Mol Biol. 1999;294(5):1351–62.View ArticlePubMedGoogle Scholar
- Ball LE, Garland DL, Crouch RK, Schey KL. Post-translational modifications of aquaporin 0 (AQP0) in the normal human lens: Spatial and temporal occurrence. Biochemistry. 2004;43:9856–65.View ArticlePubMedGoogle Scholar
- Moeller HB, Olesen ETB, Fenton RA. Regulation of the water channel aquaporin-2 by posttranslational modification. Am J Physiol Renal Physiol. 2011;300:1062–73.View ArticleGoogle Scholar
- Tamma G, Ranieri M, Di Mise A, Centrone M, Svelto M, Valenti G. Glutathionylation of the aquaporin-2 water channel: a novel post-translational modification modulated by the oxidative stress. J Biol Chem. 2014;289:27807–13.PubMed CentralView ArticlePubMedGoogle Scholar
- Ma Y, Creanga A, Lum L, Beachy PA. Prevalence of off-target effects in Drosophila RNA interference screens. Nature. 2006;443:359–63.View ArticlePubMedGoogle Scholar