- Short report
- Open Access
Silencing expression of the Rhipicephalus microplus vitellogenin receptor gene blocks Babesia bovis transmission and interferes with oocyte maturation
Parasites & Vectors volume 12, Article number: 7 (2019)
Rhipicephalus microplus is an efficient biological vector of Babesia bovis, a causative agent of bovine babesiosis. Babesia bovis is passed transovarially to the next generation of ticks, which then transmit the parasite to naïve animals. Due to the importance of the R. microplus ovary for tick reproduction and transmission of B. bovis, we investigated the hypothesis that silencing vitellogenin receptor gene expression in the ovary during tick feeding on B. bovis-infected cattle would affect parasite transmission to the next generation of ticks.
Silencing expression of the vitellogenin receptor in the ovary by RNA interference, resulted in reduced tick fertility. We observed reduced egg production (i.e. reduced weight of eggs), a lower rate of embryonic development, and a reduction in hatching. Analysis of individual larvae by PCR confirmed that RNAi mediated downregulation of the R. microplus vitellogenin receptor and also interfered with transovarial transmission of B. bovis. None of the larvae (0/58) from the RmVgR dsRNA-injected group were PCR-positive, whereas 12% (7/58) and 17% (10/58) of larvae from the non-injected and buffer-injected control groups, respectively, were infected with B. bovis.
The combined effects of reduced fecundity and reduced infection in surviving larvae resulting from silencing indicate that vitellogenin receptor is essential for tick reproduction and may play a vital role in B. bovis transmission.
Rhipicephalus microplus is a blood-feeding ectoparasite of great importance to animal health, causing significant economic losses to the livestock industry, both directly and by transmitting pathogenic agents including Babesia bovis [1,2,3]. Larvae of R. microplus attach to cattle and feed to repletion before molting to the nymphal stage. Nymphs feed to repletion on the same animal and molt to adult stages. After completing a final blood meal, fully engorged females drop off to lay their eggs in the environment. Oogenesis is the development of oocytes in the ovaries, while vitellogenesis is the process whereby yolk proteins, known as vitellogenin (Vg) are synthesized in extraovarian tissues, transported by hemolymph to the ovaries and deposited inside tick oocytes . In hard ticks, Vg is a high molecular mass-precursor that is synthesized in the fat bodies and midgut [5,6,7,8,9,10]. Vg is released into the hemolymph and is taken up by oocytes via a receptor mediated endocytosis . Vg accumulates inside oocytes yolk granules as vitellin, which is a critical nutrient for tick embryo development . A vitellogenin receptor (VgR) gene has been identified in a variety of tick species including Dermacentor variabilis , Haemaphysalis longicornis , Amblyomma hebraeum . R. microplus and R. appendiculatus .
The sequence of the VgR of R. microplus (RmVgR) gene contains an open reading frame of 5400 base pairs encoding a large protein of 1800 amino acids. The amino acid sequence of RmVgR has high similarity with VgR from Ixodes scapularis, D. variabilis and A. hebraeum. Silencing with RmVgR double stranded RNA (dsRNA) reduced VgR expression resulting in reduced fertility and a decrease in the number of larvae produced . In a previous report, silencing of VgR affected transovarial transmission of B. gibsoni by female H. longicornis fed on infected dogs, as indicated by the absence of parasite DNA in the egg mass laid by the group injected with dsRNA . Similar to B. gibsoni, B. bovis is transmitted transovarially to the next tick generation. It was unknown if silencing of RmVgR gene would affect transovarial transmission of B. bovis. To fill this knowledge gap, we tested the hypothesis that RNAi-mediated silencing interfered with tick fitness and affected B. bovis transovarial transmission to the next generation of ticks. This study enhances our understanding of the essential role of RmVgR in tick reproductive physiology, and offers an additional candidate target for the development of strategies to control the tick and B. bovis transmission.
Evaluation of RmVgR transcripts by RT-PCR demonstrated that transcripts of VgR were dramatically downregulated in tick ovaries derived from the RmVgR-dsRNA group. In contrast, RmVgR transcripts were detectable in both control groups (Fig. 1). Rm α-tubulin PCR was used as a positive control to demonstrate the presence of RNA transcripts in the tested samples. These results confirm that RmVgR-dsRNA effectively targeted and silenced transcription of the VgR. Tick phenotypes including tick engorgement, oviposition, and tick hatchability were also investigated in order to evaluate effects resulting from silencing expression of the vitellogenin receptor (Table 1). There was no significant difference in the adult survival rate of females after injection or in the weight of female ticks or in the number of female ticks fed to repletion between all three groups. The adult survival rate of females was 92% in non-injected group, 90% in buffer-injected group and 91% in RmVgR-dsRNA-injected group. The average weight of engorged 306.27 ± 5.21 mg for RmVgR-dsRNA-injected female ticks and 306.27 ± 6.42 mg for the buffer-injected group (t(179) = 0. 6, P = 0.9996). The average weight of engorged non-injected female ticks was 306.25 ± 6.14 mg. The average weight of engorged non-injected female ticks was 306.25 ± 6.14 mg compared to RmVgR-dsRNA injected group (t(181) = 1.7, P = 0.9986) (Table 1). Overall more than 90% of females from the three groups successfully fed to repletion (Table 1) (Fig. 2a). However, examination of ovary tissues confirmed that in the RmVgR-dsRNA group the ovaries were abnormal and there was a reduced number of ovarioles (Fig. 2b). Furthermore, females in the RmVgR-dsRNA group laid irregularly formed eggs (Fig. 2c).
More than 90% of female ticks in all three groups laid eggs. However, there were statistically significant differences in the average weight of the egg masses between tick groups (P < 0.0001) (Fig. 3) (Table 1). The average egg diameter in non-injected and buffer control groups were 0.379 mm and 0.376 mm, respectively (Table 1). In contrast, the average egg diameter in the RmVgR-dsRNA-injected group was 0.256 mm, a significant difference was found between groups (P < 0.0001) (Fig. 4a, b). We also observed a statistically significant difference in egg viability between RmVgR-dsRNA injected females as compared to non-injected and buffer control ticks. Control ticks had viabilities of 93% and 92%, respectively, whereas 39.1% of egg masses were viable from the dsRNA-injected group (P < 0.0001) (Fig. 5a, b).
There was no significant difference between groups in infection of adult female ticks with B. bovis. The overall female infection rate was around 70% in all groups demonstrating that RmVgR silencing did not affect B. bovis acquisition during tick feeding (Table 2). However, the rates of larval infection with B. bovis were 12% (7/58) and 17% (10/58) for non-injected and buffer control groups, respectively. In contrast, there were no larvae infected with B. bovis from the RmVgR-dsRNA group (P < 0.005), confirming that transovarial transmission of B. bovis to the offspring was diminished (Table 2).
In a previous report, the phenotypic changes in the ovary from ticks of the RmVgR silenced group indicated that deposition of yolk is mediated by VgR . Similar results, using a different region of VgR gene for silencing, have been reported for A. hebraeum and R. microplus where injection with VgR-dsRNA disturbed yolk deposition and reduced tick reproduction . That data showed silencing of RmVgR impaired tick fertility, supporting a crucial role for VgR in the uptake of yolk protein by the ovary and in aiding egg formation during tick reproduction. Our data from the RmVgR-dsRNA-injected group corroborates the previous report and confirms that silencing RmVgR reduces tick reproductive fitness.
In the present study, we tested the hypothesis that silencing of vitellogenin receptor gene expression in the ovary during tick feeding on B. bovis infected cattle would affect parasite transmission to the next generation of ticks. Under normal conditions B. bovis is transmitted through infection of ovary epithelial cells by kinetes, resulting in transovarial transmission to tick eggs. Within the larvae, B. bovis develops in the salivary glands and transforms into sporozoites, an infectious stage for the mammalian host.
In this study, downregulation of VgR reduced transovarial transmission of B. bovis. As expected, due to the fact that VgR is expressed mainly by ovarian cells, there was no difference in female infection rate of B. bovis as determined by the presence of kinetes in tick hemolymph. VgR knockdown had no impact on adult female tick acquisition of B. bovis from the infected animal, as demonstrated by accumulation of kinetes in the tick hemolymph. These results, together with the success of female repletion, suggest that RmVgR-dsRNA had no effect on tick midgut epithelial cells. However, although a high female tick infection rate was found in the RmVgR-dsRNA group, the parasite was not transovarially transmitted to the next tick generation. These results are similar to the previous report using H. longicornis VgR (HlVgR) , in which silencing of HlVgR transcripts by RNAi caused abnormal oocyte development and led to atypical eggs and blocked transovarian transmission of B. gibsoni .
We postulate that VgR may be a potential mechanism for B. bovis kinete attachment and invasion of tick ovary epithelial cells resulting in infection of the next generation of ticks by transovarial transmission. One possible mechanism is that B. bovis kinete attach directly to the VgR on the ovary epithelial cells; an alternative mechanism is that kinete become coated with Vg which then facilitates indirect attachment to the Vg receptor on ovary epithelial cells. Based on these, and previously reported results, we are generating RmVgR specific synthetic peptides in order to formulate an experimental vaccine which will be tested in bovines in order to test whether RmVgR antibodies in vaccinated animals may affect tick fitness and parasite transmission.
We have demonstrated that RmVgR silencing critically reduced oocyte development and egg maturation, resulting in decreased tick reproduction fitness. Moreover, it appears that RmVgR silencing directly reduced transovarial transmission of B. bovis. This conclusion was supported by the absence of B. bovis infected larvae from the RmVgR-dsRNA group. Therefore, RmVgR may be a good potential target for the development of novel strategies to control R. microplus and babesiosis by reducing tick reproduction and blocking pathogen transmission.
Synthesis of double stranded RNA
Synthesis of dsRNA was performed as previously described . Specific primers were designed based on the RmVgR gene sequence (GenBank: KR401221) described in the R. microplus Gene Index Project  (Table 3). The sequence on GenBank (KR401221) was analyzed in silico for the identification of a sequence that is optimal for performing small interfering RNA (siRNA) using the algorithm siRNA Target Finder available from Ambion (Austin, TX, USA). A fragment of 361 bp from RmVgR gene containing the highest number of putative siRNA was identified as optimal and amplified by PCR and cloned in to PCR II-Topo (Invitrogen, Foster City, CA, USA) and used as a template for dsRNA synthesis. Reverse transcription PCR (RT-PCR) was used to amplify the RmVgR fragment from ovary mRNA. In brief, tick ovaries were removed with fine-tipped forceps and homogenized in TRIzol® reagent (Invitrogen). Total RNA was extracted following the manufacturer’s protocol and treated with DNase I (Invitrogen). Total RNA yield was quantified by Nanodrop (Thermo Fisher Scientific, Waltham, MA, USA). One hundred-fifty ng samples of total RNA were utilized to synthesize cDNA using the Superscript III kit (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions, and the cDNA stored at -20 °C. The amplified fragment was at the 5' end of the RmVgR gene and approximately 361 bp in length. PCR products were cloned into pCR™II-TOPO® (Invitrogen), sequenced and used for in vitro transcription. The MEGAscript® Transcription Kit (Ambion) was used for dsRNA synthesis following the manufacturer’s protocol. The RmVgR dsRNA molecules were confirmed by electrophoresis on agarose gel, quantified by spectrophotometry and kept at -20 oC until used for tick injection.
Cattle, ticks and protozoan
A splenectomized Holstein calf (C90979) 3–4 months of age, that tested negative for B. bovis by PCR  and cELISA , was used in this study. The animal was maintained according to protocols approved by the University of Idaho Institutional Animal Care and Use Committee. Ticks from our laboratory colony, R. microplus La Minita strain  were used. To obtain unfed adult ticks for the experiment, approximately 40,000 larvae from 2 g of eggs were placed under a cloth patch to feed. On day 13, engorged nymphs were manually removed and held in an incubator at 25 oC with 96% relative humidity to molt to adults.
After 3 days of incubation, freshly molted adult ticks were sorted by sex. Three groups of 100 female ticks were used in this study, 1) non-injected control group, 2) buffer 1× phosphate-buffered saline (PBS)-injected group and 3) RmVgR-dsRNA-injected group. In brief, individual female ticks were injected through the coxal membrane at the base of leg 4 on the right ventral side, as previously described [17, 21] with 1 μl containing either ~1 × 1011 dsRNA molecules or buffer control. The injections were accomplished using a 10 μl syringe with a 33 gauge needle (Hamilton, Bonaduz, Switzerland) and a microprocessor controlled UMP3 injection pump apparatus (World Precision Instruments, Berlin, Germany). Each group was placed, with an equal number of males, under separate stockinet sleeves glued to the back of the calf to feed. One day after applying ticks, the calf was inoculated with B. bovis S74T3Bo strain  stabilate containing approximately 1.4 × 107 B. bovis-infected erythrocytes so that the feeding of adult female ticks to repletion would be synchronized with rising B. bovis parasitemia (acquisition feeding). 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 stained blood smears and PCR as previously described .
Silencing of RmVgR after dsRNA injection was investigated using RT-PCR. Twenty-four partially and fully engorged adult female ticks were manually collected per group from the B. bovis infected calf, rinsed with 70% ethanol and dissected with sterile scalpel blades. Tick ovaries were harvested and total RNA extracted. cDNA was synthesized and PCR performed using RmVgR and Rm α-tubulin primers (Table 3) as previously described . PCR amplicons were separated by electrophoresis on 2% agarose gels and visualized under UV trans-illumination. Amplicons were TA cloned into PCR 2.1-TOPO® (Thermo Fisher Scientific) and submitted for sequencing (Eurofins MWG Operon, Louisville, KY).
Detection of Babesia bovis kinetes
Fully engorged female ticks (n = 58 per group) from non-injected, buffer-injected, and RmVgR-dsRNA-injected groups were examined to determine tick infection rates as previously described , In brief, hemolymph from individual ticks was sampled on day 8 after dropping. A distal leg segment was removed, and a drop of exuded hemolymph was placed onto a glass slide and stained with Giemsa as previously described  for the presence of B. bovis kinetes. A minimum of 50 high-power fields per sample were observed by light microscopy.
Evaluation of tick fitness
Fully engorged female ticks were collected, weighed and put in individual wells in 24-well plates at 26 oC. After 14 days of incubation, female ticks began oviposition. After completion of oviposition egg masses from individual female ticks were weighed and placed in individual vials. The embryonic development rate was determined by the percentage of eggs with an embryo per female. Effects derived from the knockdown of the VgR were estimated by measuring egg diameter, size and viability using a Leica MZ12.5 microscope. Pictures were taken by SPOT Insight 11.2 Color Mosaic (Diagnostic Instruments, Inc.) and analyzed by SPOT 5.2 Software (Diagnostic Instruments, Inc.).
Detection of Babesia bovis transovarial transmission
Whole-genomic DNA was isolated from individual unfed R. microplus larva from the three groups. Larvae were collected, washed in 75% ethanol, and stored at -20 oC. Frozen individual larva were triturated in 20 μl phosphate-buffered saline (pH=7.2) and the homogenate transferred into a 1.5 ml microcentrifuge tube to which 300 μl of cell lysis buffer (Puregene, Valencia, CA, USA) and 15 μl of Proteinase K (2 mg/ml stock) was added and vortexed. The homogenate was incubated at 56 °C overnight, then 200 μl cold Protein Precipitation Solution (Puregene) was added to the lysate, vortexed and incubated on ice for 5 min. The lysate was centrifuged at 13,200× rpm for 3 min. The supernatant was poured to a clean screw-cap 1.5 ml tube with 500 μl isopropanol and mixed by inverting 50 times. After incubation at room temperature for 10 min, the genomic DNA was precipitated by centrifugation at 13,200× rpm for 5 min. Thereafter, 2 washes were performed with 300 μl of 70% ethanol at 13,200× rpm for 5 min. The pellets were air dried for 10 min and suspended with 30 μl of TE buffer and stored at -20 oC.
Nested primer sets were used to amplify the B. bovis rap1 gene (Table 3) using the previously described method . PCR products were separated by electrophoresis on 2% agarose gels and visualized with UV trans-illumination. All the PCR amplicons were cloned into PCR 2.1-TOPO® (Thermo Fisher Scientific) and submitted for sequencing (Eurofins MWG Operon, Louisville, KY).
Weight of egg masses, egg diameter, embryonic development rate and parasite infection rate were compared by Student’s t-test (GraphPad Instat®, version 3.06, GraphPad Software, Inc., San Diego, CA, USA).
Double stranded RNA
Rhipicephalus microplus vitellogenin receptor
Reverse transcription PCR
Dantas-Torres F, Chomel BB, Otranto D. Ticks and tick-borne diseases: a One Health perspective. Trends Parasitol. 2012;28:437–46.
Roberts JA. Resistance of cattle to the tick Boophilus microplus (Canestrini). II. Stages of the life cycle of the parasite against which resistance is manifest. J Parasitol. 1968;54:667–73.
Grisi L, Leite RC, JRdS M, ATMd B, Andreotti R, PHD C, et al. Reassessment of the potential economic impact of cattle parasites in Brazil. Rev Bras Parasitol Vet. 2014;23:150–6.
Bergamo Estrela A, Seixas A, de Oliveira Nunes Teixeira V, AFM P, Termignoni C. Vitellin- and hemoglobin-digesting enzymes in Rhipicephalus (Boophilus) microplus larvae and females. Comp Biochem Physiol B Biochem Mol Biol. 2010;157:326–35.
Rosell R, Coons LB. The role of the fat body, midgut and ovary in vitellogenin production and vitellogenesis in the female tick Dermacentor variabilis. Int J Parasitol. 1992;22:341–9.
Thompson DM, Khalil SMS, Jeffers LA, Sonenshine DE, Mitchell RD, Osgood CJ, Michael Roe R. Sequence and the developmental and tissue-specific regulation of the first complete vitellogenin messenger RNA from ticks responsible for heme sequestration. Insect Biochem Mol Biol. 2007;37:363–74.
Boldbaatar D, Umemiya-Shirafuji R, Liao M, Tanaka T, Xuan X, Fujisaki K. Multiple vitellogenins from the Haemaphysalis longicornis tick are crucial for ovarian development. J Insect Physiol. 2010;56:1587–98.
Khalil SMS, Donohue KV, Thompson DM, Jeffers LA, Ananthapadmanaban U, Sonenshine DE, et al. Full-length sequence, regulation and developmental studies of a second vitellogenin gene from the American dog tick, Dermacentor variabilis. J Insect Physiol. 2011;57:400–8.
Xavier MA, Tirloni L, Pinto AFM, Diedrich JK, Yates JR, Mulenga A, et al. A proteomic insight into vitellogenesis during tick ovary maturation. Sci Rep. 2018;8:4698.
Barrero RA, Guerrero FD, Black M, McCooke J, Chapman B, Schilkey F, et al. Gene-enriched draft genome of the cattle tick Rhipicephalus microplus: assembly by the hybrid Pacific Biosciences/Illumina approach enabled analysis of the highly repetitive genome. Int J Parasitol. 2017;47:569–83.
Giorgi F, Bradley JT, Nordin JH. Differential vitellin polypeptide processing in insect embryos. Micron. 1999;30:579–96.
Raikhel AS, Dhadialla TS. Accumulation of yolk proteins in insect oocytes. Annu Rev Entomol. 1992;37:217–51.
Mitchell Iii RD, Ross E, Osgood C, Sonenshine DE, Donohue KV, Khalil SM, et al. Molecular characterization, tissue-specific expression and RNAi knockdown of the first vitellogenin receptor from a tick. Insect Biochem. Mol Biol. 2007;37:375–88.
Boldbaatar D, Battsetseg B, Matsuo T, Hatta T, Umemiya-Shirafuji R, Xuan X, Fujisaki K. Tick vitellogenin receptor reveals critical role in oocyte development and transovarial transmission of Babesia parasite. Int J Biochem Cell Biol. 2008;86:331–44.
Smith AD, Reuben Kaufman W. Molecular characterization of the vitellogenin receptor from the tick, Amblyomma hebraeum (Acari: Ixodidae). Insect Biochem. Mol Biol. 2013;43:1133–41.
Seixas A, Alzugaray MF, Tirloni L, Parizi LF, Pinto AFM, NWo G, et al. Expression profile of Rhipicephalus microplus vitellogenin receptor during oogenesis. Ticks Tick Borne Dis. 2018;9:72–81.
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. Parasit Vectors. 2009;2:57.
Guerrero FD, Miller RJ, Rousseau ME, 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–95.
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.
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.
Hussein HE, Scoles GA, Ueti MW, Suarez CE, Adham FK, Guerrero FD, Bastos RG. Targeted silencing of the Aquaporin 2 gene of Rhipicephalus (Boophilus) microplus reduces tick fitness. Parasit Vectors. 2015;8:618.
Rodriguez SD, Buening GM, Green TJ, Carson CA. Cloning of Babesia bovis by in vitro cultivation. Infect Immun. 1983;42:15–8.
Ueti MW, Palmer GH, Kappmeyer LS, Scoles GA, Knowles DP. Expression of Equi Merozoite antigen 2 during development of Babesia equi in the midgut and salivary gland of the vector tick Boophilus microplus. J Clin Microbiol. 2003;41:5803–9.
Figueroa JV, Chieves LP, Johnson GS, Buening GM. Multiplex polymerase chain reaction based assay for the detection of Babesia bigemina, Babesia bovis and Anaplasma marginale DNA in bovine blood. Vet Parasitol. 1993;50:69–81.
We thank Kathy Mason, Paul Lacy, Sara Davis, Ralph Horn and James Allison for excellent technical assistance.
This work was supported by USDA-ARS project number # 2090-32000-039-00D.
Availability of data and materials
All data generated or analyzed in this study are included within the article.
This study was approved by the Institutional Animal Care and Use Protocol Committees of the University of Idaho, Moscow, Idaho in accordance with institutional guidelines based on the U.S. National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.