Rhipicephalus bursa Sialotranscriptomic Response to Blood Feeding and Babesia ovis Infection: Identification of Candidate Protective Antigens

Introduction

Ticks are widely distributed obligate hematophagous ectoparasites, which have recognized effects on host species. During blood feeding, ticks secrete varying substances into the host bloodstream acting as remarkable vectors of numerous pathogens, some of which can cause severe diseases in vertebrate hosts, including humans (Jongejan and Uilenberg, 2004; Domingos et al., 2013; Sonenshine and Michael Roe, 2014). Reflecting the progress of feeding, salivary glands (SG) increase ~25-fold in mass and protein content, as the glands are responsible for the production of complex saliva that is capable of quelling host innate and adaptive immune responses (Sauer et al., 2000; Kazimírová and Stibraniova, 2013; Kotál et al., 2015; Šimo et al., 2017). SG play an essential role in tick survival and success as parasites by modulating host haemostasis and complement systems (Sauer et al., 2000; Francischetti et al., 2009; Kazimírová and Stibraniova, 2013). In addition to being involved with osmoregulation (Kaufman, 2010), this tissue is also responsible for the production of cement, which is an adhesive substance that surrounds the mouthparts and the host skin that ensures tick attachment (Sauer et al., 2000; Francischetti et al., 2009; Kazimírová and Stibraniova, 2013; Šimo et al., 2017). SG are also pivotal in tick pathogen interactions, because pathogens need to cross the physical barrier of SG epithelium and endure the salivary biochemical environment to gain access to the next host. Remarkably, to increase their proliferation and transmission, pathogens adapted to SG in a way that exploits tick salivary molecules (Ramamoorthi et al., 2005; Kaufman, 2010). Therefore, these features make SG an exceptional target for the identification of new candidate protective antigens that are relevant to biological functions associated with tick development, fertility, feeding, and pathogen infection and transmission (Merino et al., 2013; Shahein et al., 2013).

Research that examined tick SG made the characterization of a large number of tick salivary compounds possible, but the function of several of these molecules remains unknown (Francischetti et al., 2009). The sialomes of some tick species have been described (Francischetti et al., 2008, 2011; Anatriello et al., 2010; Karim et al., 2011; Tan et al., 2015; de Castro et al., 2016; Moreira et al., 2017), and this information represents an important data source for functional studies and analyses of gene expression dynamics during tick feeding. Moreover, high-throughput technologies have also enabled researchers to study the effects of sex, physiological stages, and different tick statuses such as the presence of pathogens in tick tissues (Chmelar et al., 2016).

Rhipicephalus bursa is a multi-host tick that is mainly associated with ruminants, but it can occasionally parasitize other animals such as wild ungulates and small mammals (Walker et al., 2000; de la Fuente et al., 2004; Santos-Silva et al., 2011; Mihalca et al., 2012). R. bursa is recognized as the primary vector of Babesia ovis (Moltmann et al., 1982a), but it transmits other pathogens such as Rickettsia spp. and Anaplasma spp. (Raele et al., 2015; Dahmani et al., 2016; Ferrolho et al., 2016b), thus demonstrating its importance in animal health, particularly in livestock. B. ovis, an intraerythrocytic apicomplexan parasite, is the main etiological agent of ovine babesiosis, which is a tick-borne disease of small ruminants, and its geographical distribution overlaps with that of R. bursa (Walker et al., 2000; Ranjbar-Bahadori et al., 2012; Erster et al., 2015; Ferrolho et al., 2016a). This highly pathogenic organism is characterized by low parasitaemia, and it causes severe infections (Habela et al., 1990; Sevinc et al., 2013; Hurtado et al., 2015). B. ovis is extremely well adapted to the vector, and it survives in the tick during several successive generations (Yeruham et al., 2001) using horizontal and vertical transmission (Friedhoff, 1988). Microscopy studies in the 1980's discovered that the B. ovis cycle within the tick is similar to other Babesia spp. (Moltmann et al., 1982a,b). Briefly, Babesia penetrates the tick midgut, undergoes meiosis, and differentiates into motile ookinetes that propagate via haemolymph to reach all tick organs. B. ovis kinetes reach SG within 48 h post-infestation, and they undergo a final step of multiplication to produce sporozoites (Moltmann et al., 1982a; Antunes et al., 2017). Adult ticks are the main vector, and both females and males are implicated in the transmission of the hemoparasite. However, females present a higher threat due to transovarial transmission and extended feeding periods (Friedhoff, 1988).

The importance of the R. bursa-B. ovis system was emphasized in a disease outbreak that resulted in animal morbidity and mortality (Hurtado et al., 2015). Pathogen and vector control methods are limited to the common usage of imidocarb dipropionate (to manage animal disease) and acaricides (McHardy et al., 1986; Belloli et al., 2006; Domingos et al., 2013). Safer and effective alternatives are urgently needed, including the development of vaccines that may reduce tick infestations and block pathogen transmission (Merino et al., 2013; Liu and Bonnet, 2014; Neelakanta and Sultana, 2015). Studies of the molecular interactions associated with the tick-pathogen interface represent a bridge for the identification of antigenic targets to implement vaccination strategy. Information about the R. bursa and B. ovis interactome is scarce. Thus, in the present study, SG of R. bursa adult females were used to assess the transcriptomic response to blood feeding and B. ovis infection. Fed-infected, fed-uninfected, and unfed-uninfected female ticks were produced, SG were isolated and used for RNA extraction. RNA-seq and de novo transcriptome assembly approaches were used to construct the sialotranscriptome of fed-infected, fed-uninfected, and unfed-uninfected R. bursa specimens. These catalogs were analyzed, and four genes were selected for further functional studies, thus allowing the evaluation of encoded proteins for inclusion in anti-tick and tick-borne pathogen vaccines. These data are essential for vaccinomics pipelines, which could enhance our knowledge of the dynamic processes that occur at the tick-pathogen-host interface.

Materials and Methods

Ethics Statement

Animal experiments were conducted with the approval of the Divisão Geral de Alimentação e Veterinária (DGAV), Portugal, under Art° 49, Portaria n°1005/92 from 23rd October (permit number 0421/2013) and the Council of Ethics of the Instituto de Higiene e Medicina Tropical (IHMT). Animals were maintained and manipulated following protocols compliant with the national and European Animal Welfare legislation, in frame with DL 113/2013 and Directive 2010/63/EUbased on the principle of the Three R's, to replace, reduce, and refine the use of animals for scientific purposes.

Rhipicephalus bursa Colony

R. bursa colony was established under laboratory conditions and further maintained. For colony initiation, adult ticks were collected either in naturally infested domestic animals or by dragging/flagging the vegetation and kept in a chamber regulated at 25 ± 1°C, 70 ± 10% relative humidity and a photoperiod of 16:8 (light: dark). During oviposition, the dark period was increased to improve female egg laying. After oviposition, each female and a sample of eggs were tested by conventional PCR for pathogens detection (Babesia spp., Anaplasma spp., Ehrlichia spp.) during two generations using the protocols and primers described elsewhere (Inokuma et al., 2000; de la Fuente et al., 2003; Aktaş et al., 2005; Harrus et al., 2011). Ticks were fed on Hyla breed rabbits at Centro de Estudos de Vetores e Doenças Infeciosas, Instituto Nacional de Saúde Doutor Ricardo Jorge (CEVDI/INSA) in appropriate conditions. Ten lineages of R. bursa were selected in order to reduce interbreeding.

In Vitro Babesia ovis Cultures

In vitro B. ovis cultures were established at IHMT in biosafety level 2 facilities, following a protocol adapted from Vega et al. (1985). Briefly, cryopreserved B. ovis (Israeli strain) infected red blood cells (RBC) were used to initiate the culture. B. ovis merozoites were cultured in lamb erythrocytes maintained in 20% lamb serum-containing medium, in an atmosphere of 5% CO₂/2% O₂/93% N₂ at 37°C, as described elsewhere (Horta et al., 2014). Half of the medium was replaced daily and cultures monitored for parasitaemia by preparing thin blood smears stained with Hemacolor^® Rapid staining of blood smear (EMD Millipore, Darmstadt, Germany). Intraerythrocytic parasites were observed under a 400x original magnification of a Nikon eclipse 80i fluorescence microscope.

Salivary Glands and RNA Samples for RNA-Seq

Fed and Unfed R. bursa

Thirty adult female ticks were carefully removed from the rabbits ear 10–12 days post attachment. Equally, thirty unfed adult female ticks were also obtained. Ticks were individually rinsed in distilled water, after in 75% (v/v) ethanol, once more in water and dissected under a stereoscopic microscope at 4x magnification (Motic SMZ-171B, China) using sterile conditions in ice-cold phosphate-buffered saline (PBS). The SG were stored in RNAlater (Ambion, Austin, TX, USA) and afterwards pooled, resulting in two samples for the fed condition and other two for the unfed. Total RNA was extracted from each sample using Tri-reagent (Sigma–Aldrich, St. Louis, MO, USA). RNA quantity was estimated using the ND-1000 Spectrophotometer (NanoDrop ND1000, Thermo Fisher Scientific, Whaltman, MA, USA).

Fed-B. ovis Infected R. bursa

A batch of 60 female ticks were inoculated with B. ovis in the trochanter—coxae articulation and allowed to feed on rabbits. After drop off, SG were carefully isolated and DNA/RNA extracted has mentioned previously. Genomic DNA was used to amplify a 549 bp fragment of B. ovis 18S ribosomal DNA (18S rRNA) using primers and conditions described elsewhere (Aktaş et al., 2005). RNA from positive samples (Supplementary Figure 1) were used for the production of two RNA pools with fifteen samples each. All samples were promptly shipped in dry ice to Parque Cientifico de Madrid for sequencing. The tick infection model and vector competence was evaluated. B. ovis inoculated R. bursa were allowed to feed in a naïve lamb. The lamb was monitored every two days for babesiosis clinical symptoms and blood collected for B. ovis detection by PCR (Supplementary Figure 1) using the above mentioned conditions. After 8 days, the ticks were recovered for analysis.

RNA-Seq

RNA quality was assessed using an Agilent RNA 6000 bioanalyzer (Agilent Technologies, CA, USA). Libraries preparation was performed with “NEBNext Ultra Directional RNA Library Prep” kit (New England Biolabs, Ipswich, MA, USA) following manufacturer instructions. Briefly, prior to cDNA library construction magnetic beads with oligo (dT) were used to enrich poly (A) mRNA from 1 μg of total-RNA. Next, the purified mRNAs were disrupted into short fragments, and double-stranded cDNAs were immediately synthesized. The cDNAs were subjected to end-repair and adenilation, then connected with sequencing adapters. Suitable fragments, purified by size selection protocol with AMPure XP beads (Beckman Coulter), were selected as templates for PCR amplification. The final library sizes and qualities were evaluated electrophoretically using an Agilent High Sensitivity DNA kit (Agilent Technologies, CA, USA); the mean fragment size was 510 bp. Subsequently, the library was sequenced using a HiSeq 2500 sequencer (Illumina, CA, USA) in rapid run mode. Cluster generation was performed, followed by 2 × 100 cycle sequencing reads separated by a paired-end turnaround. Image analysis was performed using the HiSeq control software version 1.8.4. The raw fastq files were deposited in the Sequence Read Archives (SRA) of the National center for Biotechnology information (NCBI) under the accession numbers SRR4428986, SRR4428987 and SRR4428988, Biosamples SAMN05916213, SAMN05916214, and SAMN05916215, regarding the unfed-uninfected, fed-uninfected and fed-infected populations, respectively, of Bioproject PRJNA348674. The Transcriptome Shotgun Assembly (TSA) projects have been deposited at DDBJ/EMBL/GenBank under the accessions GFZD00000000, GFZJ00000000, and GFZK00000000. The versions described in this paper are the first versions, GFZD01000000, GFZJ01000000, and GFZK01000000.

Transcriptomic Data of Female R. bursa Sialome

Assembly and Analysis of Transcripts

This project comprised de novo assembly of six transcriptomes. Three conditions and two replicas per condition: F, SG from fed ticks; NFni, SG from unfed-uninfected ticks; and Fi, SG from fed-B. ovis infected ticks. Subsequently, two comparisons were performed: F vs. NFni (response to blood feeding) and F vs. Fi (response to B. ovis infection). Quality analysis of the raw reads was done with Prinseq tool (Schmieder and Edwards, 2011). Pre-processing of reads included: (a) right trimming where quality < Q30; (b) left trimming of the first base; (c) filtering out reads with Ns; (d) quality analysis of the processed data. For each of the four transcriptomes three de novo assemblies were made with three different k-values using the de novo transcriptome assembler Oases (Velvet, version: 1.2.10) (Schulz et al., 2012). The annotation of each transcript was done based on the Basic Local Alignment Search Tool (BLAST) results comparing the transcript to a database of reference proteins. The set of reference proteins was selected from UniProt database from all the organisms belonging to the taxon “Ixodidae”. In total 76, 475 proteins were used as reference proteins. A set of unigenes for each sample was obtained. The assignment of each transcript to a protein was based on BLAST similarity. Rich functional annotation for each unigene extracted from the UniProt protein in which the read clustering process has been centered for this unigene is provided. Afterwards a unigene expression quantification was performed using eXpress. To compare the transcripts from the samples, the transcripts were clustered by protein. The protein driven transcript clusters that were done using UniProt proteins, were furtherly clustered by UniRef90 proteins. The mapping from the UniProt proteins to UniRef90 was done using UniProt retrieval tool. The quantification per UniRef90 cluster was calculated adding the quantification per protein included in each UniRef90 cluster. P-value calculation of the Z-test was based on the raw counts (total exon reads per gene). Genes were considered significantly differentially expressed if the P-value was below 0.05. Functional annotation of these genes was manually done by compiling information from UniProt, RefSeq, GO, Panther, KEGG, Pfam, and NCBI databases.

Gene Ontology Assignments

Functional data for each identified protein was obtained using Blast2GO platform version 4.0.7 available at https://www.blast2go.com (Conesa et al., 2005; Götz et al., 2008). Homology to the protein sequences was searched by BLAST against Arthropoda (nr subset) [arthropoda, taxa:6656] from 30.01.2017 as well as against to InterPro protein signature databases, using InterProScan. To retrieve gene ontology (GO) terms, a mapping step was performed gathering GO annotations and evidence codes (EC). Annotation to assign functional terms was performed next. At this step, the most specific and reliable annotation was considered. Finally, to map a set of annotations to high level GO terms, GO slim option was used. GO frequency charts were constructed using the Microsoft Office 2016 Excel tool. The most up and down-regulated genes in response to feeding and infection (P < 0.1) were analyzed using the same approach.

Validation of RNA-Seq Data

A total of 18 transcripts with differential regulation and belonging to different functional classes with a potential interference in response to blood feeding and B. ovis infection, were chosen for RNA-Seq validation through qPCR using the minimum information for publication of qPCR experiments (Bustin et al., 2009). Ten individual R. bursa SG, from each condition studied, were used to extract total RNA using the GRS FullSample Purification kit, Grisp™ (Porto, Portugal), which included DNAse treatment and 60 ng/μL of each sample were used to synthesize cDNA using the iScript™ cDNA Synthesis Kit (Bio-Rad, CA, USA). qPCR reactions of 10 μL were performed in triplicate using IQ™ SYBR^® Green Supermix kit (Bio-Rad, CA, USA) in a CFX Connect™ Real-Time PCR Detection System (Bio-Rad). The cycling conditions were as follows: an initial cycle of denaturation at 95°C for 10 min; followed by 45 cycles of 95°C for 15s and temperature of each primer set for 45s. Fluorescence readings were taken at 62°C after each cycle and a dissociation curve (60–95°C) was performed. Negative controls were prepared with water. To determine the reaction efficiency standard curves were constructed with five-fold serial dilutions of cDNA from R. bursa. Reactions specificity was assured by the absence of PCR product in control reactions and by the dissociation curves (60–95°C) run at the end the cycling protocol. The average expression stability (M-value) of the reference genes, β tubulin, β actin, elongation factor, and 16S, was assessed based in geNorm algorithm (Vandesompele et al., 2002) included in the CFX Manager™ Software (Bio-Rad, CA, USA) and gene relative quantification was evaluated using the CFX Manager™ Software including the Pfaff method (Pfaffl, 2001) using the above-mentioned reference genes for normalization. Normalized Cq-values were compared between conditions by Student's t test (P < 0.05). Primers were design using Primer3 platform (http://bioinfo.ut.ee/primer3-0.4.0/) and their conditions are summarized in Supplementary Table 1. Pearson's correlation was used to compare the expression values between RNA-Seq and qPCR methods for the 18 selected genes.

RNA Interference Assays

Lamb Infection With B. ovis

A six-month old lamb bred and maintained at the Instituto Nacional de Investigação Agrária e Veterinária (INIAV) animal facility was splenectomized and, 45 days after, intravenously inoculated with 1 mL of cryopreserved B. ovis culture with 9% parasitemia. The B. ovis infection was monitored daily by blood screening. Genomic DNA was extracted from lamb blood using the NZY Blood gDNA Isolation Kit (NZYTech, Lisboa, Portugal) as per manufacturer instructions. As previously mentioned, B. ovis infection was screened using conventional PCR with primers and conditions described elsewhere (Aktaş et al., 2005). PCRs were performed in 25 μl reactions with Supreme NZYTaq 2 × Green Master Mix (NZYTech), 1 μM primers and 5 μl of template DNA. A negative control with water and a positive B. ovis (Israeli strain) control were added. The PCR was carried out with a thermal cycling profile of 95°C for 2 min, and 35 cycles of 95°C for 30 s, 62°C for 45 s and 72°C for 45 s, followed by a 72°C extension for 5 min, in a T-100^® Thermal Cycler (Bio-Rad, CA, USA). Resulting amplicons were checked on a 0.5X TBE, 1.2% (w/v) agarose gel.

Synthesis of dsRNA

Specific primers containing T7 promoter sequences (5′-TAATACGACTCACTATAGGGTACT-3′) at the 5′- end were manually designed using as template available sequences, in particular, GACK01008016 from Rhipicephalus pulchellus, GBBO01000019 from Rhipicephalus microplus, GBBR01000108 from R. microplus, and GACK01007634 from R. pulchellus and synthesized by StabVida (Lisbon, Portugal) (Supplementary Table 2). R. bursa cDNA was synthetized using the iScript cDNA synthesis (Bio-Rad) following the manufacturer instructions and further used as template to amplify fragments of interest by PCR. Amplifications of target DNA fragments were achieved using the iProof High Fidelity PCR kit (Bio-Rad) in a 50 μl of final volume reaction, including 200 mM of each primer. Cycling conditions were for 40 cycles: 30 s at 94°C, 30 s at specific annealing temperature and 30 s at 72°C with a final extension step of 7 min at 72°C (Supplementary Table 2). All PCR assays were performed in a T100 thermal cycler (Bio-Rad). Amplification results were analyzed on a 0.5x TBE, 1.2 % (w/v) agarose gel. Amplicons were purified using the NZYGelpure kit (NZYtech) and sent for Sanger sequencing at StabVida (Lisbon, Portugal). The obtained sequences were aligned and compared to reference sequences. After validation of the amplified sequences the MEGAscript RNAi Kit (Ambion, Austin, TX, USA) was used to synthesize dsRNA according to manufacturer's instructions. The resulting dsRNA was purified and analyzed by spectrometry and agarose gel.

Inoculation of dsRNA and Tick Infestation

R. bursa adult female ticks from the established colony at CEVDI/INSA were cleaned and placed ventral side up on double sticky tape, affixed to a plane wood table. Thirty female ticks per group were injected in the trochanter-coxae articulation with 69 nL of gene specific dsRNA (1 × 10¹¹ to 1 × 10¹² molecules) or unrelated dsRNA as control, using the nano-injector (Nanoject, Drummond Scientific, PA, USA). The mouse β-2-microglobulin dsRNA (dsβ2M) (GenBank: NM_009735) was used as control (Couto et al., 2017). After dsRNA injection, female ticks were held in a humidity chamber for 4 h after which they were allowed to feed on the splenectomized lamb infected with B. ovis together with 30 male ticks per feeding cell. Tick-feeding cells (450 × 400 mm) (cotton fabric) were glued to shaved skin using Pattex^® contact glue (Henkel Nederland, Nieuwegein, Netherlands) on the day before infestation. Ticks were monitored daily and allowed to feed in the infected lamb for 8 days. After this period, attached ticks were manually removed.

Analysis of Tick Biological Parameters After Gene Knockdown

Tick mortality was evaluated as the ratio of dead ticks to the total number of initial ticks. To analyze tick mortality, the Chi-square test (P > 0.05) was used with the null hypothesis that tick mortality was independent of gene knockdown. The ability to attach to the vertebrate host was also evaluated as the ratio of attached ticks and the total number of live ticks. The Chi-square test (P > 0.05) was also used in this analysis. Tick weight was determined in individual female ticks collected after feeding and further compared between ticks injected with test genes dsRNA and control dsRNA by Student's t-test with unequal variance (P > 0.05).

Gene Knockdown Assessment and Determination of B. ovis Infection by qPCR

To assess gene knockdown efficiency in tick SG ten ticks per group were randomly selected and tissues dissected and further used to extract total RNA and DNA and synthetize cDNA, as described previously. Quantity and quality of the RNA samples was estimated using the QIAxcel Advanced system (Qiagen™, Hilden, Germany). qPCR assays were performed under the conditions aforementioned. Gene expression was analyzed by the CFX Manager™ Software (Bio-Rad) as previously referred. Infection levels in tick SG were estimated using qPCR by evaluation of the levels of B. ovis 18S ribosomal DNA (18S rRNA) normalized against tick 16S rDNA, as described previously for other Babesia spp. (Antunes et al., 2012). The primers used for detection of B. ovis were the same used previously for conventional PCR. The cycling conditions are described in the Supplementary Table 1. Normalized Cq-values were compared between ticks injected with dsRNA and control ticks by Student's t-test with unequal variance (P > 0.05).

Antigenicity Prediction

Antigenicity of the selected molecules was estimated in silico using VaxiJen Server (Doytchinova and Flower, 2007) to allow antigen classification based on the physicochemical properties of proteins without resorting to sequence alignment. Complete sequences of the proteins were retrieved from UniProt in FASTA format and antigenicity estimated using the settings of parasite as target organism and threshold level 0.4.

Results

Assembly and Annotation of Female R. bursa Sialomes

R. bursa female ticks representing the three conditions were produced and used for SG dissections, which were followed by DNA and RNA extractions. RNA qualitative and quantitative analysis are summarized in Supplementary Table 3. Infection of protozoan-exposed group (Fi) was confirmed prior to experimentation, and total RNA was used in RNA-Seq analyses. Data were collected as two sets of matched 100-bp reads and quality analysis and raw read pre-processing were performed. The de novo assembly statistics are presented in Table 1.

TABLE 1

Table 1. Assembly statistics of the six examined Rhipicephalus bursa sialotranscriptomes.

A substantial increase in the number of contigs was observed in fed-uninfected samples compared to unfed-uninfected samples. The fed-B. ovis infected samples exhibited the highest number of contigs (Table 1). Each transcript was annotated based on BLAST results that compared the transcript to a database of reference proteins. The complete list of results can be accessed in Supplementary Datasheets 1, 2.

The obtained transcriptomes were analyzed using the Blast2GO tool and a public Arthropoda database (nr subset) (arthropoda, taxa: 6656; from 30.01.2017). Molecular functions (Figure 1A) and biological processes (Figure 1B) of the three transcriptomes were analyzed.

FIGURE 1

Figure 1. Radar plots of the three transcriptomes per represented molecular functions (A) and biological processes (B). The lines represent a pattern of the three transcriptomes unfed-uninfected, fed-uninfected and fed-Babesia ovis infected, allowing a visual comparison between conditions.

The molecular functions represented in the three sialotranscriptomes included DNA, RNA, protein, and ion binding properties as well as kinase, oxidoreductase, peptidase, and transmembrane transporter activities (Figure 1A). The remaining functions represented molecular functions that were present in both fed-uninfected and fed-infected catalogs, with the exception of nucleoside-triphosphatase and structural molecule activities that were exclusive to the unfed-uninfected sialotranscriptome. Ion binding was the most represented molecular function in all three datasets (Figure 1A). Biological processes such as catabolic, cellular protein modification, single-organism cellular, small molecule metabolic processes, translation, and signal transduction were also overrepresented in all sialotranscriptomes (Figure 1B). Anatomical structure development, chromosome organization, macromolecular complex assembly, response to stress, ribonucleoprotein complex biogenesis, vesicle-mediated transport, and DNA, RNA, and lipid metabolic processes were represented in the two sialotranscriptomes associated with fed-uninfected and fed-infected conditions (Figure 1B). Single-organism development is a feeding-exclusive process, while cellular component assembly, organelle organization, transport, and nucleic acid metabolic processes were exclusive to the unfed-uninfected samples.

Profile of SG Transcriptomic Dynamics in Response to Tick Feeding and B. ovis Infection

To clarify the response of R. bursa sialotranscriptomes to Babesia infection and blood feeding, an analysis that focused on the most (P < 0.1) up-regulated and down-regulated transcripts (Supplementary Figure 2) and as well as significantly differentially expressed (P < 0.05) genes (Figures 2, 3) was conducted.

FIGURE 2

Figure 2. Rhipicephalus bursa SG transcriptional response to blood meal based on Gene Ontology functional classes assignments of encoded proteins. Yellow bars represent down regulated genes, orange bars represent up regulated genes with statistical significance (P < 0.05).

FIGURE 3

Figure 3. Rhipicephalus bursa SG transcriptional response to Babesia ovis infection based on Gene Ontology functional classes assignments of encoded proteins. Gray bars represent down regulated genes and blue bars represent up regulated genes with statistical significance (P < 0.05).

In total, 7,272 and 13,819 different expressed sequence tags (ESTs) were obtained from the SG of unfed and fed ticks, respectively. From these, 5,188 were found in both conditions, 2,884 were exclusive to the unfed population, and 8,631 were only present in the SG of fed R. bursa females. The sialotranscriptome associated with the fed-uninfected condition was compared to the fed-infected one. The results of RNA-Seq analyses indicated that 13,819 ESTs were obtained from the sialotranscriptome of the fed sample, and 15,292 ESTs were obtained from the fed-infected sample. Of these, 9,722 ESTs were present in both samples. A total of 4,097 ESTs were exclusive to the fed-uninfected ticks, and 5,570 ESTs were only present in the SG of the fed R. bursa females.

Analysis of the most up-regulated and down-regulated transcripts (P < 0.1) (Supplementary Figure 2) indicated that 500 and 216 ESTs were differentially regulated upon feeding and infection, respectively. The diversity of molecular functions and biological processes was higher in response to blood feeding compared to infection conditions. Regarding molecular functions, hydrolase activity was the only Babesia infection-exclusive function, and it was completely down-regulated. The blood-feeding exclusive functions were anion, metal ion, heterocyclic, and organic cyclic compound, and protein binding activities, and these functions were only associated with up-regulated transcripts. Regarding biological processes, B. ovis infection resulted in the induction of biosynthetic processes, cellular protein metabolic processes, gene expression, macromolecular complex assembly, organelle organization, and symbiosis (encompassing mutualism through parasitism). However, infection was also associated with the down-regulation of catabolic processes, cellular component organization, lipid metabolic and single-organism cellular process, and transmembrane transport. R. bursa blood meals predominantly induced biological processes such as oxidation-reduction, organic substance biosynthetic, and cellular biosynthetic processes, and cellular amino acid metabolic process and signal transduction were down-regulated.

SG Gene Differential Expression in Response to Blood Feeding

Fifty-two genes were considered significantly differentially expressed (P < 0.05), and these were classified based on GO for biological process and molecular functions (Figure 2). Seventy-five percent of these genes were up-regulated, and metabolism was the most up-regulated functional class in response to blood feeding. Functional classes such as transport, detoxification, and cell functions were only up-regulated, while signaling was down-regulated. Transcripts from structural, RTT (replication-transcription-translation), proteolysis, and metabolism functional classes were also differentially regulated during blood meals.

SG Gene Differential Expression in Response to B. ovis Infection

Thirty-six genes were considered differentially expressed (P < 0.05) and classified by functional classes as previously described (Figure 3). Further analyses revealed that 64 and 36% of the differentially expressed genes were up-regulated and down-regulated, respectively. Metabolism was a highly represented functional class that was associated with both up- and down-regulated genes. Structural and RTT functional classes were also affected in the R. bursa sialome by Babesia infection. Proteolysis and immunity were exclusively up-regulated, while transport was down-regulated.

Validation of RNA-Seq Results

Sixteen genes identified as differentially expressed in response to infection and blood feeding in RNA-Seq were selected for data validation by qPCR analysis. From the RNA-Seq catalog derived from the comparison of fed vs. unfed populations, nine transcripts that encoded the following proteins were selected: annexin (UniProt ID: A0A023FX57), aspartic protease (UniProt ID: Q2WFX6), yolk cathepsin (UniProt ID: Q56CZ1), a putative hydroxysteroid 17-beta dehydrogenase (UniProt ID L7M196), hirudin-like (UniProt ID: F0JA28), lachesin (UniProt ID: L7M018), lipocalin 9 (UniProt ID: A0A034WWJ8), a putative scinderin-like (UniProt ID: L7MCZ6), and vitellogenin-3 (UniProt ID: A0A034WWF8) (Figure 4). Regarding the RNA-Seq data obtained from the comparison of infected and uninfected SG, eight genes encoding the following proteins were selected: a putative chondroitin sulfate synthase 1-like (UniProt ID: V5H7Q8), lachesin (UniProt ID: L7M018), laminin receptor (UniProt ID: E2J6W6), a putative glycine rich protein (UniProt ID:L7M1K6), a mucin-like protein (UniProt ID: C9W1L9), a putative ornithine decarboxylase antizyme (UniProt ID: A0A023FCB3), a secreted cement protein (UniProt ID: A0A034WWS7), and a putative yurt (UniProt ID: V5HE08) (Figure 5). A moderate positive correlation between the mRNA levels by both RNA-Seq and qPCR methods was obtained (Pearson's correlation coefficient 0.5394, P = 0.025).

FIGURE 4

Figure 4. Differentially gene expression of Rhipicephalus bursa SG in response to blood feeding evaluated by qPCR. Red bars represent SG from fed R. bursa ticks and green bars represent the SG from unfed R. bursa ticks. *P < 0.05.

FIGURE 5

Figure 5. Differentially gene expression of Rhipicephalus bursa SG in response to Babesia ovis infection evaluation by qPCR. Red bars represent the B. ovis infected SG and green bars represent the SG from uninfected R. bursa ticks. *P < 0.05.

Selection of Genes for RNA Interference Studies

Genes for RNAi functional studies were selected based on their potential role in the condition studied and fold change of expression. The gene encoding a putative vitellogenin-3 (Vg-3) was identified herein as up-regulated in response to feeding in both RNA-Seq (fold-change 17.51, P = 0.025) and qPCR (fold-change 98.05, P < 0.001) evaluations. The GO analysis assigned the encoded protein to a lipid transporter activity function (molecular function), belonging to the lipid transport biological process. Lachesin, which was also selected for functional analysis, was found to be up-regulated in the RNA-Seq analysis (fold-change = 15.14, P = 0.045) in response to blood feeding, while it was down-regulated based on the qPCR analysis (fold-change = −3.83, P < 0.001). This gene was also identified in the transcriptomic response to infection (fold-change = −0.80, P = 0.857), so its expression during B. ovis infection was also verified by qPCR (fold-change = −2.6427, P = 0.955). Lachesin belongs to the UniRef90_A0A1E1X7K6 cluster that is related to neural cell adhesion molecules. The gene designated as secreted cement encodes a component that is potentially involved in cement cone formation and tick attachment, and it was up-regulated in response to infection based on both RNA-Seq (fold-change = 15.73, P = 0.0298) and qPCR (fold-change = 4.197, P = 0.007) results. The expression of this secreted cement protein was also characterized by qPCR in response to blood feeding, indicating high up-regulation (fold-change = 47.4, P < 0.0001) in accordance with its role in the feeding process. Lastly, an uncharacterized gene designated as glycine rich that encodes a putative glycine rich protein was selected from the catalog associated with infection response, and it was up-regulated based on the results of both RNA-Seq (fold-change = 14.76, P = 0.0382) and qPCR (fold-change = 2.931, P = 0.016) analyses.

Functional Analyses of Differentially Expressed Tick Genes in Response to Feeding and B. ovis Infection

Tick Attachment, Weight and Survival Rate After RNAi

After dsRNA injection, biological parameters such as tick mortality, attachment, and weight were determined and statistically analyzed (Table 2).

TABLE 2

Table 2. Evaluation of tick mortality, attachment, and weight after dsRNA injection in Rhipicephalus bursa ticks.

RNAi assays indicated that tick survival was significantly affected in dsRNA-injected ticks, in both dsvitellogenin (7/30; Chi-square, P < 0.001) and dslachesin (9/30; Chi-square, P < 0.001) groups compared to controls (25/30), suggesting that these genes may play an important role in tick survival. The dsvitellogenin group was most affected with the highest mortality rate (76.67%). As represented in Table 2, the dscement group was the most significantly affected by RNAi (P = 0.008), as 45.8% of the ticks were not able to correctly attach to the vertebrate host to complete blood meal. The dslachesin injected population mimicked the control group's ability to attach to the host and feed. The average body weight was also measured, and it was significantly higher in the control group (133 ± 119 mg) than the Vg-3-silenced group (40 ± 19 mg); however, no statistical study was conducted because of the low number of ticks (N = 4). Lachesin knockdown did not affect tick weight (149 ± 108 mg) (P > 0.05). The knockdown of the gene encoding the cement protein significantly reduced female weight (52 ± 46 mg, P = 0.021) and only 13 ticks were able to attach to the host.

Gene Silencing Efficiency and Babesia Infection Evaluation

Under the studied conditions, dsRNA-mediated gene knockdown efficiency and its effect on B. ovis infection was assessed (Table 3).

TABLE 3

Table 3. Efficiency of gene knockdown by RNA interference and its influence on B. ovis infection levels in Rhipicephalus bursa ticks SG.

The injection of dsRNA molecules in R. bursa ticks led to a significant reduction of vitellogenin, lachesin, and secreted cement mRNA levels in SG by 92% (P = 0.040), 51% (P = 0.047), and 65% (P = 0.018), respectively. Regarding the levels of infection acquired after feeding on an experimentally B. ovis-infected lamb, the results indicated that the knockdown of lachesin significantly reduced B. ovis infection levels by 70% (P = 0.00251) in R. bursa SG (Table 3). The remaining groups exhibited increased infection levels.

Antigenicity of vitellogenin-3, lachesin, and secreted cement proteins were predicted by VaxiJen tool selecting parasite as the target organism. The three proteins showed to be probable antigens.

Discussion

Babesiosis is one of the most important diseases transmitted by ticks that affect a wide range of vertebrates, considered an emerging zoonose (Hunfeld et al., 2008; Ord and Lobo, 2015; Antunes et al., 2017). B. ovis is a potentially lethal pathogen that is normally found in small ruminants, and it is primarily transmitted by R. bursa, a tick species that is widely distributed in the Mediterranean region (Walker et al., 2000; Ferrolho et al., 2016a). Despite the importance of the R. bursa-B. ovis-vertebrate host interactome, no studies have examined these molecular relationships. Although it is recognized that transcripts and protein levels in ticks do not always correlate because of post-transcriptional and post-translational modifications (Ayllón et al., 2015; Villar et al., 2015), transcriptomic analysis is essential for a proper understanding of the molecular constituents of cells and tissues and the interactions and relationship between parasites and disease development (Li and Biggin, 2015; Rokyta et al., 2015). The integration of different omics analyses have allowed the detailed characterization of tick-pathogen molecular pathways (Ayllón et al., 2015; Cabezas-Cruz et al., 2017a,b). Herein, to elucidate the cellular mechanisms behind blood feeding and Babesia infection, three sialotranscriptomes of R. bursa females were analyzed and SG genes were selected for further characterization with RNAi to assess their potential as tick protective antigens.

Overall Characteristics of the R. bursa Sialome in Response to Blood Feeding and Babesia Infection

A strong transcriptional response was induced after tick feeding and during B. ovis infection, since a higher and more diverse number of transcripts were detected in the fed-uninfected sample, and even more diverse transcripts were detected in the fed-infected samples (Table 1 and Figure 1A) in comparison with the unfed-uninfected SG samples. This type of response was previously described in other systems (Heekin et al., 2013; Tirloni et al., 2014; Ayllón et al., 2015; Villar et al., 2015; de Castro et al., 2016, 2017; Kim et al., 2016; Perner et al., 2016; Valdés et al., 2016; Schroeder et al., 2017), thus indicating that different tick biological processes or statuses stimulate different gene expression regulatory strategies.

Functional annotation indicated that in all transcriptomes, ion binding molecular function was the most represented category, and its representation nearly doubled in response to feeding (Figure 1A). Being obligatory hematophagous ectoparasites, ticks must deal with the iron and heme resulting from blood catabolism. Ticks are known to express iron and heme binding proteins that sequester excess iron or heme, preventing cell damage for physiologically normal cells (Galay et al., 2015; Kim et al., 2016).

Structural molecule activity is the only class more represented in the unfed-uninfected SG transcriptome, while other molecular function categories such as structural constituent of ribosome or enzyme regulator activity are exclusive to the fed-uninfected and fed-infected populations that exhibit high cellular activity (Villar et al., 2014).

The most represented biological process in all sialomes was the cellular protein modification. The transcript abundance of transcripts belonging to this biological process doubled in the fed-uninfected and fed-infected SG samples in comparison to the unfed-uninfected ones. The anatomical structure development process was only represented in the fed populations, and this possibly reflected SG enlargement during feeding as the majority of acinar cells undergo marked hypertrophy in Ixodid females (Šimo et al., 2017). Furthermore, some pathogens induce cytoskeletal rearrangement by affecting the regulation of specific mRNAs (Ayllón et al., 2013, 2015; Ireton, 2013; Cotté et al., 2014; de la Fuente et al., 2017). As expected, metabolism-related processes were markedly represented in the transcriptomes of fed samples. The response to stress was only identified in the fed-uninfected and fed-infected SG samples, and this was in accordance with previous studies that indicated high regulation of such pathways in ticks and cells infected with Anaplasma spp. (Villar et al., 2010, 2014) and during feeding (oxidative stress response) (Kim et al., 2016). The unfed sialotranscriptome profile revealed the maintenance of basal cellular metabolism (Figure 1B). Lipid metabolic processes were exclusively represented in the fed-uninfected and fed-infected samples, thus correlating with higher cellular energy requirements and saliva production (Denardi et al., 2011). Being a cellular energy source, lipids in tick SG are implied in cement cone formation, thus explaining the high representation of such metabolic activity (Denardi et al., 2011). A comparable result was obtained in Ixodes ricinus and Rhipicephalus appendiculatus SG after feeding (Kotsyfakis et al., 2015; de Castro et al., 2016). Salivary lipid interacting proteins were up-regulated in I. ricinus infected with Borrelia burgdorferi (Cotté et al., 2014) suggesting that certain pathogens can manipulate vector lipid metabolism to facilitate infection and multiplication (Perera et al., 2012; Grabowski et al., 2017).

Specific R. bursa Sialome Response to Blood Feeding

Few studies have focused on the sialotranscriptomic response to tick feeding (McNally et al., 2012; Kotsyfakis et al., 2015; Yu et al., 2015; de Castro et al., 2016, 2017; Maruyama et al., 2017), but all demonstrated that transcription was highly affected in SG. Kotsyfakis et al. (2015) showed that fed I. ricinus, SG exhibit 10 times more overexpression compared to the midgut. Herein, genes that were highly differentially expressed in response to blood meals indicated up-regulation at rates of 75.0% (P < 0.05) to 83.8% (P < 0.1). GO analyses revealed that expression of secreted proteins was induced during tick feeding, including 14 lipocalins, four metalloproteases, two glycine rich proteins, and three microplusins (Supplementary Datasheet 1, Supplementary Figure 2). Such transcriptional regulation differs throughout tick feeding, thus reflecting the necessity of the tick to first attach to the host, evade and modulate host immune defenses, and maintain this status during the prolonged feeding period (McNally et al., 2012; Kotsyfakis et al., 2015; Chmelar et al., 2016; de Castro et al., 2017). Furthermore, fatty-acid related transcripts were highly represented in the up-regulated SG genes, suggesting a significant investment in carbohydrate metabolism. After tick attachment, SG differentiate and convert from an inactive to a metabolically active status with intense biosynthesis of molecules and ion transport, which increase cell energy requirements (McNally et al., 2012). The most up-regulated transcripts identified herein using RNA-Seq analyses encoded a fatty acid synthase (fold-change = 17.67), followed by vitellogenin-3 (fold-change = 17.51) and a glycine-rich cell wall structural protein (fold-change = 17.51). Two uncharacterized proteins (fold-changes = −17.66 and −16.63) and two glycine rich proteins (fold-changes = −16.63 and −15.71) encoded transcripts were highly down-regulated (Supplementary Datasheet 1). These results suggested that in the late stage of feeding, female ticks switch the regulation of specific proteins related to the production of cement cone, thus driving drop-off in accordance with previous reports (McNally et al., 2012; Kotsyfakis et al., 2015; de Castro et al., 2017).

Specific R. bursa Sialome Response to B. ovis Infection

The sialotranscriptomes of fed-infected and fed-uninfected female R. bursa were compared to characterize SG transcriptional regulation in response to pathogen infection. As all of the SG samples belonged to fed ticks, the effect of the feeding process can be annulled. Some studies aimed to understand the effects of pathogens on tick SG at transcriptomic, proteomic, and metabolomic levels (Nene et al., 2004; Zivkovic et al., 2010; Mercado-Curiel et al., 2011; McNally et al., 2012; Cotté et al., 2014; Ayllón et al., 2015; Villar et al., 2015; Valdés et al., 2016). Because of their medical importance, many of these studies were dedicated to Anaplasma spp./Borrelia spp.-Ixodes spp. interactions This is the first study that specifically focused on the Rhipicephalus SG transcriptomic response to Babesia infection. Pathogens highly adapted to the vector such as Anaplasma-R. microplus do not induce great effects on SG, while pathogens that pose a higher threat to vector fitness would lead to a greater gene modulation (Cen-Aguilar et al., 1998; Zivkovic et al., 2010; Mercado-Curiel et al., 2011; Chmelar et al., 2016; de la Fuente et al., 2016; Šimo et al., 2017). In Babesia infections, tick development tends to be impaired, but adaptive parasite tolerance has been described in R. microplus (Cen-Aguilar et al., 1998; Antunes et al., 2017). Furthermore, a small number of genes were considered differentially expressed (36 genes at P < 0.05 and 260 genes at P < 0.1), suggesting the long co-evolution of R. bursa and B. ovis. In both analyses an up-regulation of 63–64% of the genes occurred. Our results showed that during Babesia invasion, cellular metabolism tended to increase, whereas biosynthesis and protein processing were the most represented categories (Supplementary Datasheet 2, Supplementary Figure 2). This metabolism induction was previously demonstrated in other vector-pathogen systems (Mercado-Curiel et al., 2011; Heekin et al., 2012; Ayllón et al., 2015; Villar et al., 2015). The most up-regulated genes found were related to glycine-rich proteins (GRPs), including uncharacterized protein (fold-change = 17.53), glycine rich proteins (fold-change = 16.45 and 15.65), and secreted cement protein (fold-change = 15.73). Glycine rich proteins have been identified as upregulated in response to infection and cement proteins (Nene et al., 2004; Zivkovic et al., 2010). With rare exceptions, the role of such proteins during pathogen infection/dissemination have not been investigated (Trimnell et al., 2002). Lipocalins and defensins were identified as up-regulated in our dataset, showing an investment of the tick in the immune response, as expected. To validate the RNA-Seq results, qPCR was employed targeting putatively down-regulated ornithine decarboxylase antizyme, lachesin, and chondroitin sulfate synthase genes and putative upregulated laminin receptor, yurt, glycine rich, secreted cement, and mucin. Chondroitin sulfate synthase and lachesin expression trends were not confirmed, indicating up-regulation in infected SG. Chondroitin's are known to be involved in Plasmodium spp. adhesion to cells (Dinglasan et al., 2007; Couto et al., 2017), so the up-regulation of related molecules in infected tick SG suggests that Babesia spp. (considered a Plasmodium-like parasite) may use similar strategies to invade cells.

Functional Studies for the Identification of Tick Protective Antigens

Vitellogenin-3

Multiple vitellogenins (Vgs) have been described in ticks (Thompson et al., 2007; Boldbaatar et al., 2010; Khalil et al., 2011; Taheri et al., 2014; Rodriguez et al., 2016), and they are involved in detoxification and oxidative molecular processes (Galay et al., 2015). In the sialotranscriptome obtained in response to blood feeding, the translation of one of the assembled transcripts showed high similarity to R. microplus putative Vg-3 protein (UniProt ID: A0A034WWF8). An up-regulation of the expression of the correspondent gene in the SG of fed R. bursa was demonstrated by both RNA-Seq and qPCR (RNA-Seq: fold-change = 17.509, P = 0.025; and qPCR fold-change = 98.05, P < 0.001), and the results were in accordance with those of previous studies (Horigane et al., 2010; Yang et al., 2015). Vgs are thought to be absent from SG, whereas heme transport and storage are thought to be dependent of the hemelipoglyco-carrier protein (CP) (Donohue et al., 2009). In ticks, both Vg proteins and CP bind heme (Logullo et al., 2002), which is a functional component of many hemoproteins, but it is cytotoxic in larger amounts (Ferrolho et al., 2016b; Hajdusek et al., 2016). The similarities between CPs and Vgs in ticks, as well as their common evolutionary origin, greatly complicate their differentiation and function assignments (Gudderra et al., 2002; Donohue et al., 2009; Boldbaatar et al., 2010). The present study showed that R. bursa possesses a gene very similar to Vgs in SG, and it shares several molecular features with CPs. Further studies are necessary to clarify Vgs classification in ticks as well as the function and localization of Vg-3 in R. bursa species as these Vgs are expressed in a tissue-specific manner in ticks (Rodriguez et al., 2016). Vg-3 knockdown experiments resulted in increased tick mortality. No statistical analyses were performed regarding feeding behaviors, body weight and Babesia infection, because of the low number of samples; however, decreased blood-uptake and increased Babesia infection was observed. Based on the principal functions associated to this type of molecule, we can suggest that a decrease in the expression of putative Vg-3 reduces heme and lipid binding and storage (Figure 6A). A deficient heme seizure may increase cellular toxicity, thus contributing for the formation of reactive oxygen species (ROS). Also, the role of Vgs on lipid transport is compromised, and this may unbalance normal energy production. Vgs have been consistently discovered has highly immunogenic molecules in Rhipicephalus ticks (Boldbaatar et al., 2008, 2010; Smith and Kaufman, 2013; Taheri et al., 2014; Rodriguez et al., 2016), and the results of the present study stimulates future research.

FIGURE 6

Figure 6. Proposed model of putative vitellogenin-3, cement protein and lachesin functions and its impact on Rhipicephalus bursa SG during feeding and Babesia ovis infection. (A) Vitellogenin-3 described function relates to heme detoxification and lipid storage contributing for cell survival. A decrease of the expression of putative vitellogenin-3 leads to deficient heme seizure, increasing the formation of reactive oxygen species (ROS) as well as cellular toxicity. Lipid storage is also compromised leading to an unbalance in the production of energy. (B) Putative cement protein is a component of the cement cone, which facilitates the tick attachment and feed on the host. An impact in the production of cement proteins leads to an incapacity of ticks to correctly attach and subsequently feed on the host, resulting in tick death and reduced blood ingestion. (C) Lachesin is a cell surface protein that as a potential role in cell adhesion, maintaining apical-basal polarity, vesicle trafficking, cell growth and survival, as well as parasite invasion. A negative manipulation of the expression of lachesin results in an abnormal cell growth and ultimately cell apoptosis, and also a decrease of Babesia spp. infection.

Putative Secreted Cement Protein and Glycine-Rich Protein

The genes encoding putative cement protein and GRP were found to be significantly up-regulated in response to B. ovis infection in R. bursa SG, in accordance with a previous study (Nene et al., 2004). The cement cone is composed of several molecules that are embedded in a proteinaceous matrix, presenting several GRPs (Bishop et al., 2002; Trimnell et al., 2005; Maruyama et al., 2010). Different species of ticks rely on different types and amounts of GRPs in order to attach and feed on their hosts. Briefly, ticks with short mouthparts need higher amounts of GRPs than those with long mouthparts. Moreover, one-host ticks present a greater variety of these proteins than ticks that feed on several hosts (Maruyama et al., 2010). A successfully knockdown was observed in cement-silenced ticks, but no silencing was demonstrated in glycine-rich dsRNA-injected ticks, suggesting that a higher concentration may be needed to reduce the expression of this gene. The cement-silenced ticks significantly affected tick attachment, feeding, and body weight (Figure 6B). The dsglycine-rich RNA inoculated group exhibited a slight decrease in these two parameters, reflecting its potential in tick feeding capacity and attachment to the host. Curiously, in both dsRNA-injected groups, an increase of Babesia levels was detected. Previous studies concerning cement cone proteins showed that immunization with these proteins significantly affected tick attachment to the host (Trimnell et al., 2005) and it reduced pathogen transmission (Labuda et al., 2006). Therefore, these two proteins are attractive targets for vaccine development.

Lachesin

Lachesin is a cell surface protein of the immunoglobulin superfamily (Karlstrom et al., 1993; Llimargas et al., 2004) that regulates organ size by influencing cell length and cell detachments, suggesting a role in cell adhesion and connection (Llimargas et al., 2004). In ticks, the gene encoding lachesin was first identified in the genome of Ixodes scapularis (Gulia-Nuss et al., 2016) and more recently in the sialotranscriptome of Amblyomma cajennense (Garcia et al., 2014), R. pulchellus (Tan et al., 2015), and R. appendiculatus (de Castro et al., 2016). However, no studies that focus on this molecule in ticks have been performed. In the present study, an assembled transcript translated to a protein highly similar to lachesin (UniProt ID: L7M018). A highly dynamic expression profile of lachesin in response to infection and feeding was found in the present study, and this observation aligned to its presumed role on cell-adhesion led to its selection for RNAi studies. Tick inoculation with dslachesin resulted in 51% gene knockdown that led to a significantly high tick mortality. Lachesin accumulates in specific invertebrate cell junctions, and it is responsible for establishing and/or maintaining cell polarity, cell adhesion, and cell-cell interactions (Tepass et al., 2001). Apical-basal polarity is subjected to tight regulation, as it is crucial during tissue formation, including vesicle trafficking machinery, morphogenesis, and modulation of epithelial cell growth and survival (Bonazzi and Cossart, 2011). Moreover, adhesive contacts between cells and the extracellular matrix appear as important landmarks for polarity. Therefore, manipulating the expression of genes involved in this processes can induce abnormal cell growth and cell apoptosis (Tepass, 2012). In addition, the lachesin knockdown resulted in lower pathogen infection in the SG. No statistical effect was demonstrated in the other biological parameters studied. Despite the tight organization of the epithelium barrier and its interactions with cellular factors that are crucial to cell-pathogen defense, a large number of pathogens have developed strategies to target host proteins involved in cell adhesion, to colonize epithelia, invade host cells, or even disrupt host barriers to facilitate access to other tissues (Bonazzi and Cossart, 2011). Thus, our results suggest that lachesin plays an important role in tick survival and also that B. ovis may require this molecule for tissue invasion (Figure 6C). This molecule appears to be good candidate for future vaccination assays, as it demonstrates a dual-effect targeting both tick and pathogen.

Conclusions

Tick and tick-borne diseases constitute a growing burden for human and animal health, stressing the urgency in the development of new effective tools to control this global threat. Due to the important role of tick SG in tick biology and pathogen transmission, the main objective of the present study was the identification and functional characterization of R. bursa SG genes involved in tick feeding and B. ovis infection. Quantitative transcriptome analysis showed lachesin and putative vitellogenin-3 has highly upregulated in response to blood meal and the genes encoding for a putative secreted cement and GRPs highly upregulated in response to B. ovis infection. RNAi studies suggest that lachesin and putative vitellogenin-3 affect tick survival while the putative cement protein has an impact in tick attachment to the host and tick weight after feeding. Moreover, B. ovis infection levels in tick SG were reduced, subsequently to lachesin knockdown. Overall the results of the present study endorse the inclusion of these proteins in vaccination trials.

Author Contributions

SA, JdlF, MMS-S, and AD designed the study; ASS and MMS-S were responsible for tick rearing and tick inoculation methodology. SA established and maintained B. ovis cultures. SA, JC, and JF performed transcriptomic analyses; SA and FR performed qPCR assays; SA, FR, JF, and JN performed the RNA interference studies; SA, JC, JdlF, and AD performed data analysis and wrote the manuscript. All authors edited and approved the final manuscript.

Funding

This research was supported by the project PTDC/CVT-EPI/4339/2012 funded by Fundação para a Ciência e Tecnologia (FCT). SA is the recipient of a post-doctoral grant supported by FCT (SFRH/BPD/108957/2015) and JC and JF are the recipients of Ph.D. grants supported by the FCT (SFRH/BD/121946/2016, SFRH/BD/122894/2016, respectively).

Conflict of Interest Statement

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

Acknowledgments

A special thanks to Abel Oliva, Ph.D. from Instituto de Tecnologia Química e Biológica, Universidade Nova de Lisboa and Varda Shkap, Ph.D. from Kimron Veterinary Institute, Israel, for providing the Israeli strain of B. ovis. To Olga Moreira, Ph.D. and all the staff from the Estação Zootécnica from INIAV that gently cooperated in the animal experiments. To Ana Borda D'Água from CEVDI/INSA that kindly contribute to the maintenance of ticks colonies. FCT for funds to GHTM—UID/Multi/04413/2013.

Supplementary Material

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

References

Aktaş, M., Altay, K., and Dumanli, N. (2005). Development of a polymerase chain reaction method for diagnosis of Babesia ovis infection in sheep and goats. Vet. Parasitol. 133, 277–281. doi: 10.1016/j.vetpar.2005.05.057

PubMed Abstract | CrossRef Full Text | Google Scholar

Anatriello, E., Ribeiro, J. M. C., de Miranda-Santos, I. K. F., Brandao, L. G., Anderson, J. M., Valenzuela, J. G., et al. (2010). An insight into the sialotranscriptome of the brown dog tick, Rhipicephalus sanguineus. BMC Genomics 11:450. doi: 10.1186/1471-2164-11-450

CrossRef Full Text | Google Scholar

Antunes, S., Galindo, R. C., Almazan, C., Rudenko, N., Golovchenko, M., Grubhoffer, L., et al. (2012). Functional genomics studies of Rhipicephalus (Boophilus) annulatus ticks in response to infection with the cattle protozoan parasite, Babesia bigemina. Int. J. Parasitol. 42, 187–195. doi: 10.1016/j.ijpara.2011.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Antunes, S., Rosa, C., Couto, J., Ferrolho, J., and Domingos, A. (2017). Deciphering Babesia-vector interactions. Front. Cell. Infect. Microbiol. 7:429. doi: 10.3389/fcimb.2017.00429

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayllón, N., Villar, M., Busby, A. T., Kocan, K. M., Blouin, E. F., Bonzon-Kulichenko, E., et al. (2013). Anaplasma phagocytophilum inhibits apoptosis and promotes cytoskeleton rearrangement for infection of tick cells. Infect. Immun. 81, 2415–2425. doi: 10.1128/IAI.00194-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Ayllón, N., Villar, M., Galindo, R. C., Kocan, K. M., Sima, R., Lopez, J. A., et al. (2015). Systems biology of tissue-specific response to Anaplasma phagocytophilum reveals differentiated apoptosis in the tick vector Ixodes scapularis. PLoS Genet. 11:e1005120. doi: 10.1371/journal.pgen.1005120

PubMed Abstract | CrossRef Full Text | Google Scholar

Belloli, C., Lai, O. R., Ormas, P., Zizzadoro, C., Sasso, G., and Crescenzo, G. (2006). Pharmacokinetics and mammary elimination of imidocarb in sheep and goats. J. Dairy Sci. 89, 2465–2472. doi: 10.3168/jds.S0022-0302(06)72320-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Bishop, R., Lambson, B., Wells, C., Pandit, P., Osaso, J., Nkonge, C., et al. (2002). A cement protein of the tick Rhipicephalus appendiculatus, located in the secretory e cell granules of the type III salivary gland acini, induces strong antibody responses in cattle. Int. J. Parasitol. 32, 833–842. doi: 10.1016/S0020-7519(02)00027-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Boldbaatar, D., Battsetseg, B., Matsuo, T., Hatta, T., Umemiya-Shirafuji, R., Xuan, X., et al. (2008). Tick vitellogenin receptor reveals critical role in oocyte development and transovarial transmission of Babesia parasite. Biochem. Cell Biol. 86, 331–344. doi: 10.1139/O08-071

PubMed Abstract | CrossRef Full Text | Google Scholar

Boldbaatar, D., Umemiya-Shirafuji, R., Liao, M., Tanaka, T., Xuan, X. N., and Fujisaki, K. (2010). Multiple vitellogenins from the Haemaphysalis longicornis tick are crucial for ovarian development. J. Insect Physiol. 56, 1587–1598. doi: 10.1016/j.jinsphys.2010.05.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Bonazzi, M., and Cossart, P. (2011). Impenetrable barriers or entry portals? The role of cell-cell adhesion during infection. J. Cell Biol. 195, 349–358. doi: 10.1083/jcb.201106011

PubMed Abstract | CrossRef Full Text | Google Scholar

Bustin, S. A., Benes, V., Garson, J. A., Hellemans, J., Huggett, J., Kubista, M., et al. (2009). The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin. Chem. 55, 611–622. doi: 10.1373/clinchem.2008.112797

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabezas-Cruz, A., Alberdi, P., Valdés, J. J., Villar, M., and de la Fuente, J. (2017a). Anaplasma phagocytophilum infection subverts carbohydrate metabolic pathways in the tick vector, Ixodes scapularis. Front. Cell. Infect. Microbiol. 7:23. doi: 10.3389/fcimb.2017.00023

PubMed Abstract | CrossRef Full Text | Google Scholar

Cabezas-Cruz, A., Espinosa, P. J., Obregón, D. A., Alberdi, P., and de la Fuente, J. (2017b). Ixodes scapularis tick cells control Anaplasma phagocytophilum infection by increasing the synthesis of phosphoenolpyruvate from tyrosine. Front. Cell. Infect. Microbiol. 7:375. doi: 10.3389/fcimb.2017.00375

PubMed Abstract | CrossRef Full Text | Google Scholar

Cen-Aguilar, J. F., Rodriguez-Vivas, R. I., Dominguez-Alpizar, J. L., and Wagner, G. G. (1998). Studies on the effect of infection by Babesia sp. on oviposition of Boophilus microplus engorged females naturally infected in the Mexican tropics. Vet. Parasitol. 78, 253–257.

PubMed Abstract | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Conesa, A., Gotz, S., Garcia-Gomez, J. M., Terol, J., Talon, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610

PubMed Abstract | CrossRef Full Text | Google Scholar

Cotté, V., Sabatier, L., Schnell, G., Carmi-Leroy, A., Rousselle, J. C., Arsène-Ploetze, F., et al. (2014). Differential expression of Ixodes ricinus salivary gland proteins in the presence of the Borrelia burgdorferi sensu lato complex. J. Proteomics 96, 29–43. doi: 10.1016/j.jprot.2013.10.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Couto, J., Antunes, S., Pinheiro-Silva, R., do Rosário, V., de la Fuente, J., and Domingos, A. (2017). Solute carriers affect Anopheles stephensi survival and Plasmodium berghei infection in the salivary glands. Sci. Rep. 7:6141. doi: 10.1038/s41598-017-06317-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahmani, M., Davoust, B., Rousseau, F., Raoult, D., Fenollar, F., and Mediannikov, O. (2016). Natural anaplasmataceae infection in Rhipicephalus bursa ticks collected from sheep in the French Basque Country. Ticks Tick Borne Dis. 8, 18–24. doi: 10.1016/j.ttbdis.2016.09.009

PubMed Abstract | CrossRef Full Text | Google Scholar

de Castro, M. H., de Klerk, D., Pienaar, R., Latif, A. A., Rees, D. J. G., and Mans, B. J. (2016). De novo assembly and annotation of the salivary gland transcriptome of Rhipicephalus appendiculatus male and female ticks during blood feeding. Ticks Tick Borne Dis. 7, 536–548. doi: 10.1016/j.ttbdis.2016.01.014

PubMed Abstract | CrossRef Full Text | Google Scholar

de Castro, M. H., de Klerk, D., Pienaar, R., Rees, D. J. G., and Mans, B. J. (2017). Sialotranscriptomics of Rhipicephalus zambeziensis reveals intricate expression profiles of secretory proteins and suggests tight temporal transcriptional regulation during blood-feeding. Parasit. Vec. 10:384. doi: 10.1186/s13071-017-2312-4

CrossRef Full Text | Google Scholar

de la Fuente, J., Antunes, S., Bonnet, S., Cabezas-Cruz, A., Domingos, A. G., Estrada-Pena, A., et al. (2017). Tick-Pathogen interactions and vector competence: identification of molecular drivers for tick-borne diseases. Front. Cell. Infect. Microbiol. 7:114. doi: 10.3389/fcimb.2017.00114

PubMed Abstract | CrossRef Full Text | Google Scholar

de la Fuente, J., Van Den Bussche, R. A., Prado, T. M., and Kocan, K. M. (2003). Anaplasma marginale msp1α genotypes evolved under positive selection pressure but are not markers for geographic isolates. J. Clin. Microbiol. 41, 1609–1616. doi: 10.1128/JCM.41.4.1609-1616.2003

CrossRef Full Text | Google Scholar

de la Fuente, J., Vicente, J., Hofle, U., Ruiz-Fons, F., de Mera, I. G. F., Van Den Bussche, R. A., et al. (2004). Anaplasma infection in free-ranging Iberian red deer in the region of Castilla-La Mancha, Spain. Vet. Microbiol. 100, 163–173. doi: 10.1016/j.vetmic.2004.02.007

PubMed Abstract | CrossRef Full Text | Google Scholar

de la Fuente, J., Villar, M., Cabezas-Cruz, A., Estrada-Pena, A., Ayllon, N., and Alberdi, P. (2016). Tick-Host-Pathogen interactions: conflict and cooperation. PLoS Pathog. 12, 5488–5488. doi: 10.1371/journal.ppat.1005488

PubMed Abstract | CrossRef Full Text | Google Scholar

Denardi, S. E., Furquim, K. C., Bechara, G. H., and Camargo-Mathias, M. I. (2011). Salivary glands of Amblyomma cajennense (Acari: Ixodidae): a histological and an ultrastructural overview. Exp. Appl. Acarol. 54, 177–189. doi: 10.1007/s10493-010-9421-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinglasan, R. R., Alaganan, A., Ghosh, A. K., Saito, A., van Kuppevelt, T. H., and Jacobs-Lorena, M. (2007). Plasmodium falciparum ookinetes require mosquito midgut chondroitin sulfate proteoglycans for cell invasion. Proc. Natl. Acad. Sci. U.S.A. 104, 15882–15887. doi: 10.1073/pnas.0706340104

PubMed Abstract | CrossRef Full Text | Google Scholar

Domingos, A., Antunes, S., Borges, L., and Rosario, V. E. (2013). Approaches towards tick and tick-borne diseases control. Rev. Soc. Bras. Med. Trop. 46, 265–269. doi: 10.1590/0037-8682-0014-2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Donohue, K. V., Khalil, S. M. S., Sonenshine, D. E., and Roe, R. M. (2009). Heme-binding storage proteins in the Chelicerata. J. Insect Physiol. 55, 287–296. doi: 10.1016/j.jinsphys.2009.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Doytchinova, I. A., and Flower, D. R. (2007). VaxiJen: a server for prediction of protective antigens, tumour antigens and subunit vaccines. BMC Bioinformatics 8:4. doi: 10.1186/1471-2105-8-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Erster, O., Roth, A., Wolkomirsky, R., Leibovich, B., Savitzky, I., Zamir, S., et al. (2015). Molecular detection of Babesia ovis in sheep and ticks using the gene encoding B. ovis surface protein D (BoSPD). Vet. Parasitol. 214, 282–288. doi: 10.1016/j.vetpar.2015.09.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrolho, J., Antunes, S., Sanches, G. S., Couto, J., Evora, P. M., Rosa, C., et al. (2016a). Ferritin 1 silencing effect in Rhipicephalus sanguineus sensu lato (Acari: Ixodidae) during experimental infection with Ehrlichia canis. Ticks Tick Borne Dis. 8, 174–184. doi: 10.1016/j.ttbdis.2016.10.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Ferrolho, J., Antunes, S., Santos, A. S., Velez, R., Padre, L., Cabezas-Cruz, A., et al. (2016b). Detection and phylogenetic characterization of Theileria spp. and Anaplasma marginale in Rhipicephalus bursa in Portugal. Ticks Tick Borne Dis. 7, 443–448. doi: 10.1016/j.ttbdis.2016.01.004

CrossRef Full Text

Francischetti, I. M., Anderson, J. M., Manoukis, N., Pham, V. M., and Ribeiro, J. M. C. (2011). An insight into the sialotranscriptome and proteome of the coarse bontlegged tick, Hyalomma marginatum rufipes. J. Proteomics 74, 2892–2908. doi: 10.1016/j.jprot.2011.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Francischetti, I. M., Mans, B. J., Meng, Z., Gudderra, N., Veenstra, T. D., Pham, V. M., et al. (2008). An insight into the sialome of the soft tick, Ornithodorus parkeri. Insect Biochem. Mol. Biol. 38, 1–21. doi: 10.1016/j.ibmb.2007.09.009

CrossRef Full Text

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Friedhoff, K. T. (1988). “Chapter 2: Transmission of Babesia,” in Babesiosis of Domestic Animals and Man, ed M. Ristic (Boca Raton, FL: CRC Press), 23–52.

Galay, R. L., Umemiya-Shirafuji, R., Mochizuki, M., Fujisaki, K., and Tanaka, T. (2015). Iron metabolism in hard ticks (Acari: Ixodidae): the antidote to their toxic diet. Parasitol. Int. 64, 182–189. doi: 10.1016/j.parint.2014.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Götz, S., García-Gómez, J. M., Terol, J., Williams, T. D., Nagaraj, S. H., Nueda, M. J., et al. (2008). High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 36, 3420–3435. doi: 10.1093/nar/gkn176

PubMed Abstract | CrossRef Full Text | Google Scholar

Grabowski, J. M., Gulia-Nuss, M., Kuhn, R. J., and Hill, C. A. (2017). RNAi reveals proteins for metabolism and protein processing associated with Langat virus infection in Ixodes scapularis (black-legged tick) ISE6 cells. Parasit. Vec. 10:24. doi: 10.1186/s13071-016-1944-0

CrossRef Full Text | Google Scholar

Gudderra, N. P., Sonenshine, D. E., Apperson, C. S., and Roe, R. M. (2002). Tissue distribution and characterization of predominant hemolymph carrier proteins from Dermacentor variabilis and Ornithodoros parkeri. J. Insect Physiol. 48, 161–170. doi: 10.1016/S0022-1910(01)00160-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Gulia-Nuss, M., Nuss, A. B., Meyer, J. M., Sonenshine, D. E., Roe, R. M., Waterhouse, R. M., et al. (2016). Genomic insights into the Ixodes scapularis tick vector of Lyme disease. Nat. Commun. 7:10507. doi: 10.1038/ncomms10507

PubMed Abstract | CrossRef Full Text | Google Scholar

Habela, M., Reina, D., Nieto, C., and Navarrete, I. (1990). Antibody response and duration of latent infection in sheep following experimental infection with Babesia ovis. Vet. Parasitol. 35, 1–10. doi: 10.1016/0304-4017(90)90111-N

PubMed Abstract | CrossRef Full Text | Google Scholar

Hajdusek, O., Sima, R., Perner, J., Loosova, G., Harcubova, A., and Kopacek, P. (2016). Tick iron and heme metabolism - new target for an anti-tick intervention. Ticks Tick Borne Dis. 7, 565–572. doi: 10.1016/j.ttbdis.2016.01.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Harrus, S., Perlman-Avrahami, A., Mumcuoglu, K. Y., Morick, D., Eyal, O., and Baneth, G. (2011). Molecular detection of Ehrlichia canis, Anaplasma bovis, Anaplasma platys, Candidatus Midichloria mitochondrii and Babesia canis vogeli in ticks from Israel. Clin. Microbiol. Infect. 17, 459–463. doi: 10.1111/j.1469-0691.2010.03316.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Heekin, A. M., Guerrero, F. D., Bendele, K. G., Saldivar, L., Scoles, G. A., Dowd, S. E., et al. (2013). The ovarian transcriptome of the cattle tick, Rhipicephalus (Boophilus) microplus, feeding upon a bovine host infected with Babesia bovis. Parasit. Vec. 6:276. doi: 10.1186/1756-3305-6-276

PubMed Abstract | CrossRef Full Text | Google Scholar

Heekin, A. M., Guerrero, F. D., Bendele, K. G., Saldivar, L., Scoles, G. A., Gondro, C., et al. (2012). Analysis of Babesia bovis infection-induced gene expression changes in larvae from the cattle tick, Rhipicephalus (Boophilus) microplus. Parasit. Vectors 5:162. doi: 10.1186/1756-3305-5-162

PubMed Abstract | CrossRef Full Text | Google Scholar

Horigane, M., Shinoda, T., Honda, H., and Taylor, D. (2010). Characterization of a vitellogenin gene reveals two phase regulation of vitellogenesis by engorgement and mating in the soft tick Ornithodoros moubata (Acari: Argasidae). Insect Mol. Biol. 19, 501–515. doi: 10.1111/j.1365-2583.2010.01007.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Horta, S., Barreto, M. C., Pepe, A., Campos, J., and Oliva, A. (2014). Highly sensitive method for diagnosis of subclinical B. ovis infection. Ticks Tick Borne Dis. 5, 902–906. doi: 10.1016/j.ttbdis.2014.07.005

CrossRef Full Text

Hunfeld, K. P., Hildebrandt, A., and Gray, J. S. (2008). Babesiosis: recent insights into an ancient disease. Int. J. Parasitol. 38, 1219–1237. doi: 10.1016/j.ijpara.2008.03.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Hurtado, A., Barandika, J. F., Oporto, B., Minguijón, E., Povedano, I., and García-Pérez, A. L. (2015). Risks of suffering tick-borne diseases in sheep translocated to a tick infested area: a laboratory approach for the investigation of an outbreak. Ticks Tick Borne Dis. 6, 31–37. doi: 10.1016/j.ttbdis.2014.09.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Inokuma, H., Raoult, D., and Brouqui, P. (2000). Detection of Ehrlichia platys DNA in brown dog ticks (Rhipicephalus sanguineus) in Okinawa Island, Japan. J. Clin. Microbiol. 38, 4219–4221. doi: 10.1016/S0304-4017(02)00289-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ireton, K. (2013). Molecular mechanisms of cell-cell spread of intracellular bacterial pathogens. Open Biol. 3:130079. doi: 10.1098/rsob.130079

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Karlstrom, R. O., Wilder, L. P., and Bastiani, M. J. (1993). Lachesin: an immunoglobulin superfamily protein whose expression correlates with neurogenesis in grasshopper embryos. Development 118, 509–522.

PubMed Abstract | Google Scholar

Kaufman, W. R. (2010). Ticks physiological aspects with implications for pathogen transmission. Ticks Tick Borne Dis. 1, 11–22. doi: 10.1016/j.ttbdis.2009.12.001

CrossRef Full Text | Google Scholar

Kazimírová, M., and Stibraniova, I. (2013). Tick salivary compounds: their role in modulation of host defences and pathogen transmission. Front. Cell. Infect. Microbiol. 3:43. doi: 10.3389/fcimb.2013.00043

PubMed Abstract | CrossRef Full Text | Google Scholar

Khalil, S. M. S., Donohue, K. V., Thompson, D. M., Jeffers, L. A., Ananthapadmanaban, U., Sonenshine, D. E., et al. (2011). Full-length sequence, regulation and developmental studies of a second vitellogenin gene from the American dog tick, Dermacentor variabilis. J. Insect Physiol. 57, 400–408. doi: 10.1016/j.jinsphys.2010.12.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, T. K., Tirloni, L., Pinto, A. F., Moresco, J., Yates, J. R., da Silva Vaz, I., et al. (2016). Ixodes scapularis tick saliva proteins sequentially secreted every 24 h during blood feeding. PLoS Negl. Trop. Dis. 10:e0004323. doi: 10.1371/journal.pntd.0004323

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Kotsyfakis, M., Schwarz, A., Erhart, J., and Ribeiro, J. M. (2015). Tissue- and time-dependent transcription in Ixodes ricinus salivary glands and midguts when blood feeding on the vertebrate host. Sci. Rep. 5:9103. doi: 10.1038/srep09103

PubMed Abstract | CrossRef Full Text | Google Scholar

Labuda, M., Trimnell, A. R., Lickova, M., Kazimirova, M., Davies, G. M., Lissina, O., et al. (2006). An antivector vaccine protects against a lethal vector-borne pathogen. PLoS Pathog. 2:e27. doi: 10.1371/journal.ppat.0020027

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. J., and Biggin, M. D. (2015). Gene expression. Statistics requantitates the central dogma. Science 347, 1066–1067. doi: 10.1126/science.aaa8332

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X. Y., and Bonnet, S. I. (2014). Hard tick factors implicated in pathogen transmission. PLoS Negl. Trop. Dis. 8:e2566. doi: 10.1371/journal.pntd.0002566

PubMed Abstract | CrossRef Full Text | Google Scholar

Llimargas, M., Strigini, M., Katidou, M., Karagogeos, D., and Casanova, J. (2004). Lachesin is a component of a septate junction-based mechanism that controls tube size and epithelial integrity in the Drosophila tracheal system. Development 131, 181–190. doi: 10.1242/dev.00917

PubMed Abstract | CrossRef Full Text | Google Scholar

Logullo, C., Moraes, J., Dansa-Petretski, M., Vaz, I. S., Masuda, A., Sorgine, M. H. F., et al. (2002). Binding and storage of heme by vitellin from the cattle tick, Boophilus microplus. Insect Biochem. Mol. Biol. 32, 1805–1811. doi: 10.1016/S0965-1748(02)00162-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Maruyama, S. R., Anatriello, E., Anderson, J. M., Ribeiro, J. M., Brandão, L. G., Valenzuela, J. G., et al. (2010). The expression of genes coding for distinct types of glycine-rich proteins varies according to the biology of three metastriate ticks, Rhipicephalus (Boophilus) microplus, Rhipicephalus sanguineus and Amblyomma cajennense. BMC Genomics 11:363. doi: 10.1186/1471-2164-11-363

CrossRef Full Text | Google Scholar

Maruyama, S. R., Garcia, G. R., Teixeira, F. R., Brandão, L. G., Anderson, J. M., Ribeiro, J. M. C., et al. (2017). Mining a differential sialotranscriptome of Rhipicephalus microplus guides antigen discovery to formulate a vaccine that reduces tick infestations. Parasit. Vec. 10:206. doi: 10.1186/s13071-017-2136-2

CrossRef Full Text | Google Scholar

McHardy, N., Woollon, R. M., Clampitt, R. B., James, J. A., and Crawley, R. J. (1986). Efficacy, toxicity and metabolism of imidocarb dipropionate in the treatment of Babesia ovis infection in sheep. Res. Vet. Sci. 41, 14–20.

PubMed Abstract | Google Scholar

McNally, K. L., Mitzel, D. N., Anderson, J. M., Ribeiro, J. M., Valenzuela, J. G., Myers, T. G., et al. (2012). Differential salivary gland transcript expression profile in Ixodes scapularis nymphs upon feeding or flavivirus infection. Ticks Tick Borne Dis. 3, 18–26. doi: 10.1016/j.ttbdis.2011.09.003

CrossRef Full Text | Google Scholar

Mercado-Curiel, R. F., Palmer, G. H., Guerrero, F. D., and Brayton, K. A. (2011). Temporal characterisation of the organ-specific Rhipicephalus microplus transcriptional response to Anaplasma marginale infection. Int. J. Parasitol. 41, 851–860. doi: 10.1016/j.ijpara.2011.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Merino, O., Antunes, S., Mosqueda, J., Moreno-Cid, J. A., Perez de la Lastra, J. M., Rosario-Cruz, R., et al. (2013). Vaccination with proteins involved in tick-pathogen interactions reduces vector infestations and pathogen infection. Vaccine 31, 5889–5896. doi: 10.1016/j.vaccine.2013.09.037

PubMed Abstract | CrossRef Full Text | Google Scholar

Mihalca, A. D., Dumitrache, M. O., Magdas, C., Gherman, C., Domsa, C., Mircean, V., et al. (2012). Synopsis of the hard ticks (Acari: Ixodidae) of Romania with update on host associations and geographical distribution. Exp. Appl. Acarol. 58, 183–206. doi: 10.1007/s10493-012-9566-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Moltmann, U. G., Mehlhorn, H., and Friedhoff, K. T. (1982a). Electron-microscopic study on the development of Babesia ovis (piroplasmia) in the salivary glands of the vector tick Rhipicephalus bursa. Acta Trop. 39, 29–40.

PubMed Abstract | Google Scholar

Moltmann, U. G., Mehlhorn, H., and Friedhoff, K. T. (1982b). Ultrastructural study of the development of Babesia ovis (Piroplasmia) in the ovary of the vector tick Rhipicephalus bursa. J. Protozool. 29, 30–38. doi: 10.1111/j.1550-7408.1982.tb02877.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Moreira, H N, Barcelos, R. M., Vidigal, P. M., Klein, R. C., Montandon, C. E., Maciel, T. E., et al. (2017). A deep insight into the whole transcriptome of midguts, ovaries and salivary glands of the Amblyomma sculptum tick. Parasitol. Int. 66, 64–73. doi: 10.1016/j.parint.2016.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Neelakanta, G., and Sultana, H. (2015). Transmission-blocking vaccines: focus on anti-vector vaccines against tick-borne diseases. Arch. Immunol. Ther. Exp. 63, 169–179. doi: 10.1007/s00005-014-0324-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Nene, V., Lee, D., Kang'a, S, Skilton, R., Shah, T., de Villiers, E., et al. (2004). Genes transcribed in the salivary glands of female Rhipicephalus appendiculatus ticks infected with Theileria parva. Insect Biochem. Mol. Biol. 34, 1117–1128. doi: 10.1016/j.ibmb.2004.07.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ord, R. L., and Lobo, C. A. (2015). Human babesiosis: pathogens, prevalence, diagnosis and treatment. Curr. Clin. Microbiol. Rep. 2, 173–181. doi: 10.1007/s40588-015-0025-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Perera, R., Riley, C., Isaac, G., Hopf-Jannasch, A. S., Moore, R. J., Weitz, K. W., et al. (2012). Dengue virus infection perturbs lipid homeostasis in infected mosquito cells. PLoS Pathog. 8:e1002584. doi: 10.1371/journal.ppat.1002584

PubMed Abstract | CrossRef Full Text | Google Scholar

Perner, J., Provazník, J., Schrenková, J., Urbanová, V., Ribeiro, J. M., and Kopáček, P. (2016). RNA-seq analyses of the midgut from blood- and serum-fed Ixodes ricinus ticks. Sci. Rep. 6:36695. doi: 10.1038/srep36695

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfaffl, M. W. (2001). A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 29:e45. doi: 10.1093/nar/29.9.e45

PubMed Abstract | CrossRef Full Text | Google Scholar

Raele, D. A., Galante, D., Pugliese, N., De Simone, E., and Cafiero, M. A. (2015). Coxiella-like endosymbiont associated to the “Anatolian brown tick” Rhipicephalus bursa in Southern Italy. Microbes Infect. 17, 799–805. doi: 10.1016/j.micinf.2015.09.011

CrossRef Full Text | Google Scholar

Ramamoorthi, N., Narasimhan, S., Pal, U., Bao, F., Yang, X. F., Fish, D., et al. (2005). The Lyme disease agent exploits a tick protein to infect the mammalian host. Nature 436, 573–577. doi: 10.1038/nature03812

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranjbar-Bahadori, S., Eckert, B., Omidian, Z., Shirazi, N. S., and Shayan, P. (2012). Babesia ovis as the main causative agent of sheep babesiosis in Iran. Parasitol Res. 110, 1531–1536. doi: 10.1007/s00436-011-2658-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Rodriguez, P. B. R., Cruz, R. R., Garcia, D. I. D., Gutierrez, R. H., Quintanilla, R. E. L., Sahagun, D. O., et al. (2016). Identification of immunogenic proteins from ovarian tissue and recognized in larval extracts of Rhipicephalus (Boophilus) microplus, through an immunoproteomic approach. Exp. Parasitol. 170, 227–235. doi: 10.1016/j.exppara.2016.10.005

CrossRef Full Text | Google Scholar

Rokyta, D. R., Margres, M. J., and Calvin, K. (2015). Post-transcriptional mechanisms contribute little to phenotypic variation in snake venoms. G3 5, 2375–2382. doi: 10.1534/g3.115.020578

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos-Silva, M. M., Beati, L., Santos, A. S., De Sousa, R., Nuncio, M. S., Melo, P., et al. (2011). The hard-tick fauna of mainland Portugal (Acari: Ixodidae): an update on geographical distribution and known associations with hosts and pathogens. Exp. Appl. Acarol. 55, 85–121. doi: 10.1007/s10493-011-9440-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Sauer, J. R., Essenberg, R. C., and Bowman, A. S. (2000). Salivary glands in ixodid ticks: control and mechanism of secretion. J. Insect Physiol. 46, 1069–1078. doi: 10.1016/S0022-1910(99)00210-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmieder, R., and Edwards, R. (2011). Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864. doi: 10.1093/bioinformatics/btr026

PubMed Abstract | CrossRef Full Text | Google Scholar

Schroeder, C. L. C., Narra, H. P., Sahni, A., Khanipov, K., Patel, J., Fofanov, Y., et al. (2017). Transcriptional profiling of Rickettsia prowazekii coding and non-coding transcripts during in vitro host-pathogen and vector-pathogen interactions. Ticks Tick Borne Dis. 8, 827–836. doi: 10.1016/j.ttbdis.2017.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Schulz, M. H., Zerbino, D. R., Vingron, M., and Birney, E. (2012). Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 28, 1086–1092. doi: 10.1093/bioinformatics/bts094

PubMed Abstract | CrossRef Full Text | Google Scholar

Sevinc, F., Guler, L., Sevinc, M., Ekici, O. D., and Isik, N. (2013). Determination of immunoreactive proteins of Babesia ovis. Vet. Parasitol. 198, 391–395. doi: 10.1016/j.vetpar.2013.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Shahein, Y. E., Abouelella, A. M., Hussein, N. A., Hamed, R. R., El-Hakim, A. E., Abdel-Shafy, S., et al. (2013). Identification of four novel Rhipicephalus annulatus upregulated salivary gland proteins as candidate vaccines. Protein J. 32, 392–398. doi: 10.1007/s10930-013-9498-x

PubMed Abstract | CrossRef Full Text | Google Scholar

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

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, A. D., and Kaufman, W. R. (2013). Molecular characterization of the vitellogenin receptor from the tick, Amblyomma hebraeum (Acari: Ixodidae). Insect Biochem. Mol. Biol. 43, 1133–1141. doi: 10.1016/j.ibmb.2013.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Sonenshine, D. E., and Michael Roe, R. (2014). Biology of Ticks. New York, NY: Oxford University Press.

Google Scholar

Taheri, M., Nabian, S., Ranjbar, M., Fard, R. M. N., Sadeghian, A. G., and Sazmand, A. (2014). Study of vitellogenin in Boophilus annulatus tick larvae and its immunological aspects. Trop. Biomed. 31, 398–405.

PubMed Abstract | Google Scholar

Tan, A. W. L., Francischetti, I. M. B., Slovak, M., Kini, R. M., and Ribeiro, J. M. C. (2015). Sexual differences in the sialomes of the zebra tick, Rhipicephalus pulchellus. J. Proteomics 117, 120–144. doi: 10.1016/j.jprot.2014.12.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Tepass, U. (2012). The apical polarity protein network in Drosophila epithelial cells: regulation of polarity, junctions, morphogenesis, cell growth, and survival. Annu. Rev. Cell Dev. Biol. 28, 655–685. doi: 10.1146/annurev-cellbio-092910-154033

PubMed Abstract | CrossRef Full Text | Google Scholar

Tepass, U., Tanentzapf, G., Ward, R., and Fehon, R. (2001). Epithelial cell polarity and cell junctions in Drosophila. Annu. Rev. Genet. 35, 747–784. doi: 10.1146/annurev.genet.35.102401.091415

PubMed Abstract | CrossRef Full Text | Google Scholar

Thompson, D. M., Khalil, S. M. S., Jeffers, L. A., Sonenshine, D. E., Mitchell, R. D., Osgood, C. J., et al. (2007). 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. 37, 363–374. doi: 10.1016/j.ibmb.2007.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Tirloni, L., Reck, J., Terra, R. M. S., Martins, J. R., Mulenga, A., Sherman, N. E., et al. (2014). Proteomic analysis of cattle tick Rhipicephalus (Boophilus) microplus saliva: a comparison between partially and fully engorged females. PLoS ONE 9:e94831. doi: 10.1371/journal.pone.0094831

PubMed Abstract | CrossRef Full Text | Google Scholar

Trimnell, A. R., Davies, G. M., Lissina, O., Hails, R. S., and Nuttall, P. A. (2005). A cross-reactive tick cement antigen is a candidate broad-spectrum tick vaccine. Vaccine 23, 4329–4341. doi: 10.1016/j.vaccine.2005.03.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Trimnell, A. R., Hails, R. S., and Nuttall, P. A. (2002). Dual action ectoparasite vaccine targeting ‘exposed' and ‘concealed' antigens. Vaccine 20, 3560–3568. doi: 10.1016/S0264-410X(02)00334-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Valdés, J. J., Cabezas-Cruz, A., Sima, R., Butterill, P. T., and RuŽek, D, Nuttall, P. A. (2016). Substrate prediction of Ixodes ricinus salivary lipocalins differentially expressed during Borrelia afzelii infection. Sci. Rep. 6:32372. doi: 10.1038/srep32372

PubMed Abstract | CrossRef Full Text | Google Scholar

Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:research0034. doi: 10.1186/gb-2002-3-7-research0034

PubMed Abstract | CrossRef Full Text | Google Scholar

Vega, C. A., Buening, G. M., Green, T. J., and Carson, C. A. (1985). In vitro cultivation of Babesia bigemina. Am. J. Vet. Res. 46, 416–420.

Google Scholar

Villar, M., Ayllon, N., Alberdi, P., Moreno, A., Moreno, M., Tobes, R., et al. (2015). Integrated metabolomics, transcriptomics and proteomics identifies metabolic pathways affected by Anaplasma phagocytophilum infection in tick cells. Mol. Cell Proteomics 14, 3154–3172. doi: 10.1074/mcp.M115.051938

PubMed Abstract | CrossRef Full Text | Google Scholar

Villar, M., Ayllon, N., Busby, A. T., Galindo, R. C., Blouin, E. F., Kocan, K. M., et al. (2010). Expression of heat shock and other stress response proteins in tcks and cultured tick cells in response to Anaplasma spp. infection and heat shock. Int. J. Proteomics 2010:657261. doi: 10.1155/2010/657261

CrossRef Full Text | Google Scholar

Villar, M., Popara, M., Ayllon, N., Fernandez de Mera, I. G., Mateos-Hernandez, L., Galindo, R. C., et al. (2014). A systems biology approach to the characterization of stress response in Dermacentor reticulatus tick unfed larvae. PLoS ONE 9:e89564. doi: 10.1371/journal.pone.0089564

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker, J. B., Keirans, J. E., and Horak, I. G. (2000). The genus Rhipicephalus (Acari, Ixodidae): a guide to the brown ticks of the World. Trop. Anim. Health Prod. 32, 417–418. doi: 10.1017/CBO9780511661754

CrossRef Full Text | Google Scholar

Yang, X., Yu, Z., He, Y., Xu, X., Gao, Z., Wang, H., et al. (2015). Purification of vitellin and dynamics of vitellogenesis in the parthenogenetic tick Haemaphysalis longicornis (Acari: Ixodidae). Exp. Appl. Acarol. 65, 377–388. doi: 10.1007/s10493-014-9866-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yeruham, I., Hadani, A., and Galker, F. (2001). The effect of the ovine host parasitaemia on the development of Babesia ovis (Babes, 1892) in the tick Rhipicephalus bursa (Canestrini and Fanzago, 1877). Vet. Parasitol. 96, 195–202. doi: 10.1016/S0304-4017(00)00433-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, X., Gong, H., Zhou, Y., Zhang, H., Cao, J., and Zhou, J. (2015). Differential sialotranscriptomes of unfed and fed Rhipicephalus haemaphysaloides, with particular regard to differentially expressed genes of cysteine proteases. Parasit. Vec. 8:597. doi: 10.1186/s13071-015-1213-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zivkovic, Z., Esteves, E., Almazan, C., Daffre, S., Nijhof, A. M., Kocan, K. M., et al. (2010). Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus) microplus in response to infection with Anaplasma marginale. BMC Genomics 11:12. doi: 10.1186/1471-2164-11-186

CrossRef Full Text | Google Scholar