Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus)microplus in response to infection with Anaplasma marginale
BMC Genomics volume 11, Article number: 186 (2010)
Bovine anaplasmosis, caused by the rickettsial tick-borne pathogen Anaplasma marginale (Rickettsiales: Anaplasmataceae), is vectored by Rhipicephalus (Boophilus)microplus in many tropical and subtropical regions of the world. A. marginale undergoes a complex developmental cycle in ticks which results in infection of salivary glands from where the pathogen is transmitted to cattle. In previous studies, we reported modification of gene expression in Dermacentor variabilis and cultured Ixodes scapularis tick cells in response to infection with A. marginale. In these studies, we extended these findings by use of a functional genomics approach to identify genes differentially expressed in R. microplus male salivary glands in response to A. marginale infection. Additionally, a R. microplus-derived cell line, BME26, was used for the first time to also study tick cell gene expression in response to A. marginale infection.
Suppression subtractive hybridization libraries were constructed from infected and uninfected ticks and used to identify genes differentially expressed in male R. microplus salivary glands infected with A. marginale. A total of 279 ESTs were identified as candidate differentially expressed genes. Of these, five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. The impact of selected tick genes on A. marginale infections in tick salivary glands and BME26 cells was characterized by RNA interference. Silencing of the gene encoding for putative flagelliform silk protein (100Silk) resulted in reduced A. marginale infection in both tick salivary glands and cultured BME26 cells, while silencing of the gene encoding for subolesin (4D8) significantly reduced infection only in cultured BME26 cells. The knockdown of the gene encoding for putative metallothionein (93 Meth), significantly up-regulated in infected cultured BME26 cells, resulted in higher A. marginale infection levels in tick cells.
Characterization of differential gene expression in salivary glands of R. microplus in response to A. marginale infection expands our understanding of the molecular mechanisms at the tick-pathogen interface. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. These tick genes found to be functionally relevant for tick-pathogen interactions will likely be candidates for development of vaccines designed for control of both ticks and tick-borne pathogens.
Bovine anaplasmosis, caused by the obligate intracellular rickettsial pathogen, Anaplasma marginale (Rickettsiales: Anaplasmataceae), is characterized in cattle by anemia, fever, weight loss and reduced milk production . Transmission of A. marginal e occurs mechanically by biting flies and blood-contaminated fomites, while ticks are biological vectors . Approximately 20 tick species have been incriminated worldwide as vectors of A. marginale. Of these tick species, the southern cattle tick, Rhipicephalus (Boophilus)microplus, found in tropical and subtropical regions of the world, is considered to be the most economically important ectoparasite of livestock . R. microplus vectors several pathogens and A. marginale is among the most important, causing notable economic loss in milk and beef production .
The developmental cycle of A. marginale was described in Dermacentor ticks, and male ticks were shown to become persistently infected with A. marginale and able to transmit infection repeatedly when transferred among cattle . The A. marginal e development, multiplication in the tick and transmission to the vertebrate host are coordinated with tick feeding. Within Dermacentor ticks, A. marginale undergoes a complex developmental cycle that begins with the infection of the gut cells. After a second tick feeding, many other tick tissues become infected, including the salivary glands from where the A. marginale is transmitted to cattle . Although the developmental cycle of A. marginale has not been described for Rhipicephalus (Boophilus) spp., the developmental cycle in R. microplus is most likely similar and males may also play an important role in pathogen transmission .
Molecular interactions at the tick-pathogen interface ensure survival and development of both the pathogen and vector. While recent studies on several pathogens have demonstrated that tick gene expression is modified in response to pathogen infection [8–10], information on the function of the differentially expressed genes is limited . RNA interference (RNAi) has been shown to be a useful tool for the characterization of the function of genes involved in tick-host-pathogen interactions and the transmission of tick-borne pathogens and for screening for tick protective antigens . Recently, genes differentially expressed in cultured IDE8 tick cells in response to A. marginale infection were identified and their impact on pathogen infection in D. variabilis ticks was characterized by RNAi during the pathogen developmental cycle .
Tick cell lines have been used successfully to study vector-pathogen interactions . However, these studies were conducted in the IDE8 and ISE6 tick cell lines derived from Ixodes scapularis embryos which is not a natural vector of A. marginale. Recently, a Brazilian isolate of A. marginale was propagated successfully in the BME26 cell line derived from R. microplus which provided the opportunity to study the A. marginale-tick interface in the cells cultured from a natural tick vector.
The objective of this study was to identify R. microplus genes differentially expressed in male salivary glands in response to infection with A. marginale by using suppression subtractive hybridization (SSH) approach and to characterize the function of those genes by RNAi. SSH enables identification of low-abundant rare transcripts through the comparison of two cDNA populations by selective amplification of the genes expressed in one population but not in the other [14, 15]. The results of these SSH studies were validated by real-time RT-PCR in R.microplus ticks and cultured BME26 tick cells for selected genes. Finally, functional analyses were conducted on selected genes by RNAi in both R. microplus male ticks and cultured BME26 cells to determine the putative role of these genes in A. marginale-tick interactions.
Identification of differentially expressed genes in R. microplus male salivary glands
Two SSH libraries, forward and reverse, were constructed to identify genes up- and down-regulated in R. microplus male salivary glands in response to A. marginale infection. From each library, 288 randomly selected clones were identified and sequenced. After removing vector sequences and eliminating EST clones with poor sequence quality, datasets of 128 and 151 ESTs from forward and reverse subtracted libraries were obtained, respectively, and used for bioinformatics analyses. Clustering and assembly of ESTs from forward subtracted library (up-regulated in infected cells) resulted in 43 unique transcripts of which 10 were derived from two or more ESTs (consensus sequences) and 33 were derived from a single EST (singletons). Assembly of the ESTs in the reverse subtracted library (down-regulated in infected cells) yielded 56 unique sequences (24 consensus sequences and 32 singletons). Automated annotation was then used to search public domain protein databases for putative functions (Additional file 1: Table S1). Gene ontology assignments demonstrated that up-regulated genes encoded for proteins with putative functions of binding (21%), structural molecules (11%), catalytic/enzymatic activity (6%), DNA/RNA metabolism (4%), and 58% had no known function (Figure 1A). Putative functions assigned to down-regulated genes included binding (20%), structural molecules (20%), catalytic/enzymatic activity (7%), transport (5%) and 48% had no known function (Figure 1B).
Differential gene expression in A. marginale-infected R. microplus male salivary glands and cultured BME26 cells
Fourteen candidate genes with putative functions in tick-pathogen interactions were selected for validation of SSH results by real-time RT-PCR (Additional file 2: Table S2). Real-time RT-PCR analyses were done on the same pooled RNA samples used for SSH. Statistically significant differences in expression were obtained for 5 genes (Figure 2). Similar to the SSH results, genes encoding for putative female-specific histamine-binding protein (22Hbp), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz), and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were significantly down-regulated in infected tick salivary glands. The gene encoding for the putative von Willebrand factor (94Will), identified to be up-regulated by SSH, was shown by real-time RT-PCR to be significantly down-regulated in the infected tick salivary glands. For the other 9 genes, mRNA levels were not significantly different between infected and uninfected ticks. Subolesin (4D8), used as a positive control, was down-regulated in A. marginale-infected tick salivary glands.
To evaluate the use of cultured BME26 cells for studying A. marginale-tick interactions, the same primers were used for real-time RT-PCR analysis of uninfected and A. marginale-infected BME26 tick cells. Twelve of the 14 selected genes were amplified from BME26 cultured cells. Gene expression profiles were studied for each gene at 6, 24 and 72 hours post-infection (hpi) and compared with uninfected cells collected at the same time points (Figure 3). Genes encoding for putative vacuolar H+-ATPase (36vATP) and ribosomal protein S29 (59Rib) were significantly up-regulated at 6 hpi, while putative Kunitz-like protease inhibitor precursor (108Kunz), metallothionein (93 Meth) and von Willebrand factor (94Will) were significantly up-regulated 24 and 72 hpi. The mRNA levels for the rest of the genes were not statistically different between infected and uninfected BME26 cells. The subolesin (4D8) control was significantly downregulated by 24 hpi.
All the tick sequence-derived primers were tested against bovine RNA by RT-PCR. Amplicons were not obtained for any of the primer pair tested (data not shown).
Functional roles of genes differentially expressed in R. microplus ticks and cultured BME26 cells in response to infection with A. marginale
The five genes corroborated by real-time RT-PCR to be differentially expressed in A. marginale infected tick salivary glands, female-specific histamine-binding protein (22Hbp), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz), proline-rich protein BstNI subfamily 3 precursor (7BstNI3) and von Willebrand factor (94Will) (Figure 2), were selected for functional analyses in ticks. The effect of gene knockdown on A. marginale infection and multiplication in R. microplus male tick salivary glands was evaluated by RNAi. The mRNA levels after RNAi were reduced for putative von Willebrand factor (94Will), flagelliform silk protein (100Silk) and subolesin (4D8), while silencing of the other genes did not result in statistically significant differences from the control (Table 1).
The effect of RNAi of selected genes on male B. microplus mortality was determined. Tick mortality was significantly higher in groups injected with dsRNA for von Willebrand factor (94Will), flagelliform silk protein (100Silk) and subolesin when compared with the unrelated dsRNA-injected control ticks (Table 1). Despite the fact that individual variation in gene expression affected the statistical significance of results, silencing of the genes encoding for putative von Willebrand factor (94Will), flagelliform silk protein (100Silk) and subolesin (4D8) resulted in 63%, 79% and 60% decrease in A. marginale infection levels in R. microplus male salivary glands, respectively (Table 1). Characterization of the effect of gene knockdown on A. marginale infections in cultured BME26 tick cells was attempted for all the genes which were shown to be differentially expressed in tick cells and/or tick salivary glands. However, mRNA levels were reduced only for putative ribosomal protein S29 (59Rib), metallothionein (93 Meth), flagelliform silk protein (100Silk) and subolesin (4D8) genes (Table 2). Of these genes, two genes encoding for putative flagelliform silk protein (100Silk) and subolesin resulted in 12% and 17% reduction of A. marginale infection levels, respectively, when compared with controls (Table 2). The knockdown of the gene encoding for putative metallothionein (93 Meth) resulted in higher A. marginale infection levels in tick cells (Table 2).
In the present study we identified R. microplus male salivary gland genes differentially expressed in response to A. marginale infection by use of SSH and real-time RT-PCR. Development and multiplication of A. marginale in salivary gland cells involves molecular interactions between pathogen- and tick-derived molecules. Salivary gland, the tissue of interest in this study, is a critical site in the developmental cycle from where the pathogen is transmitted to cattle. Recently, tick salivary gland proteins were shown to play a role in the infection and transmission of Borrelia burgdorferi[16, 17], A. phagocytophilum and A. marginale. A. marginale membrane surface proteins involved in tick salivary gland colonization have been identified and partially characterized [20, 21]. Understanding the molecular mechanisms of A. marginale-tick interactions for R. microplus, one of the most important vectors of A. marginale worldwide, is fundamental toward development of novel control measures .
Some of the genes identified by SSH, including those genes encoding for putative tick cement proteins, female specific histamine binding protein, IgG binding protein C, salivary gland-associated protein 64P, flagelliform silk protein and von Willebrand factor, were identified previously in different tick species and appear to be involved in tick feeding or pathogen infection [10, 23–25]. However, most of the differentially expressed genes identified in this study have not been shown to be associated with tick-pathogen interaction previously. Some cellular functions affected by A. marginale infection of R. microplus, such as cell structure and enzymatic processes, were reported previously in infected tick IDE8 cultured cells . The discrepancy observed for some studied genes between SSH and real-time RT-PCR results may reflect differences between both methods for identifying differentially expressed genes or the presence of multiple sequences targeted during RT-PCR reactions that affect the results of mRNA quantification for some genes.
In a recent study, genes differentially expressed in cultured IDE8 tick cells in response to A. marginale infection were identified and functional studies conducted in D. variabilis suggested that these genes may play different roles during pathogen infection, development and trafficking from midguts to salivary glands . Some of the genes identified by de la Fuente et al.  such as gluthathione S-tranferase, selenoprotein M and ferritin were also shown to be differentially expressed in R. microplus salivary glands in response to A. marginale infection. However, these genes were absent from the current EST dataset which could be due to differences in the system used for EST discovery (cultured IDE8 tick cells versus R. microplus salivary glands) and/or other factors such as tick species and/or A. marginale strain and infection levels.
While tick cell lines have been used successfully in A. marginale functional genomics studies , this is the first report of the use of the BME26 tick cell line derived originally from a natural vector of A. marginale for functional studies of tick-pathogen interactions. Since these studies were conducted on ticks and tick cells of the same species, most of the genes identified in tick salivary glands were also amplified from cultured BME26 tick cells. However, expression profiles of selected genes observed in cultured BME26 cells were not identical to that found in tick salivary glands. For example, the expression of the putative von Willebrand factor (94Will) was down-regulated in tick salivary glands but up-regulated in cultured BME26 tick cells infected with A. marginale. These differences may have resulted from tissue-specific regulation of gene expression or because we only observed early stages of infection in the cultured BME26 tick cells (6-72 hpi). As reported previously , results of studies using cultured tick cells must be validated in naturally infected ticks. Interestingly, expression of putative vacuolar H+-ATPase (36vATP) was significantly up-regulated in A. marginale-infected cultured BME26 cells, as reported for previous gene expression studies of cultured IDE8 cells in response to A. marginale infection .
RNAi was used in this study to assign the effect of selected gene knockdown on A. marginale infection and multiplication in ticks. Although statistically significant for flagelliform silk protein (100Silk) only, results of RNAi experiments suggested that putative von Willebrand factor (94Will), flagelliform silk protein (100Silk) and subolesin could play a role in pathogen infection of R. microplus male salivary glands. RNAi experiments in cultured BME26 tick cells provided further evidence that flagelliform silk protein (100Silk) and subolesin may play a role in A. marginale infection and/or multiplication in tick cells and suggested that metallothionein (93 Meth) may be involved tick defense against pathogen infection.
The flagelliform silk protein was identified previously in tick and orb weaving spider salivary glands but its function was not linked to pathogen infection [26–28]. Mulenga et al.  demonstrated that the flagelliform silk protein may be involved in tick attachment. In previous studies of I. ricinus after B. burgdorferi infection, the von Willebrand factor was isolated from tick salivary glands and shown to be up-regulated but its possible role in infection was not studied . A von Willebrand factor-like motif is present in the major hemelipoglycoprotein found in ixodid ticks and this protein has been shown to play a role as a heme-sequestering factor during tick feeding . Therefore, silencing of these genes may affect tick feeding, mortality and development of A. marginale in salivary glands. However, as shown previously for subolesin , gene expression studies in cultured BME26 tick cells have provided evidence that that the flagelliform silk protein may play a role in the infection of ticks with A. marginale.
The results for gene expression and silencing of subolesin in R. microplus male salivary glands and cultured BME26 cells infected with A. marginale reported herein are in agreement with previous studies in which subolesin knockdown reduced A. marginale infection in D. variabilis and cultured IDE8 cells [11, 19]. Subolesin, discovered as a tick protective antigen in I. scapularis, has been shown to be conserved in many tick species [31, 32]. Subolesin was shown by both RNAi gene knockdown and immunization trials using the recombinant protein to protect vertebrate hosts against tick infestations, reduce tick survival and reproduction, and cause degeneration of gut, salivary gland, reproductive tissues and embryos [31–37]. Targeting of subolesin by RNAi or vaccination also decreased tick vector capacity for A. marginale and A. phagocytophilum. In addition, subolesin was shown to function in the control of gene expression in ticks [38, 39] and to be differentially expressed in Anaplasma-infected ticks and cultures tick cells [11, 40]. However, subolesin expression in R. microplus tick salivary glands and cultured BME26 cells was different to previous reports showing up-regulation in A. marginale-infected D. variabilis salivary glands and cultured IDE8 cells . These differences could be due to tick species-specific differences in gene regulation or to other factors such as pathogen strain and infection levels. Nonetheless, these results expanded our knowledge on the role of subolesin in tick-Anaplasma interactions.
Metallothioneins are a family of low molecular weight proteins with a high affinity for divalent metals that function in cell detoxification, apoptosis, stress response and immunity [41–43]. Metallothioneins control the cellular zinc ion levels, which are known to be important in the immune system, and their expression has been associated with protective response against pathogens [44–48]. The results suggested a role for tick metallothioneins in defense against bacterial infections. Interestingly, selenoproteins that regulate the levels of another important trace mineral in the organism were suggested to participate in the cellular response to limit A. marginale infection in tick cells .
Although dsRNA sequences used in this study do not contain any significant overlap with other known R. microplus genes, the possibility of off-target gene silencing effects cannot be excluded due to the limited amount of sequence data available. However, RNAi seems to be very sequence-specific in ticks with little off-target effects . Availability of the complete R. microplus genome sequence data will facilitate screening for potential off-target effects. These can subsequently be minimized by avoiding the use of dsRNAs or siRNAs containing sequences which are present in multiple genes.
In our study R. microplus male salivary gland genes differentially expressed in response to A. marginale infection were identified by using SSH approach. Recently a R. microplus microarray (NimbleGen) has been developed and used for the analysis of acaricide- inducible genes in R. microplus. Microarray chip hybridization could be an alternative approach for identifying R. microplus differentially expressed genes in response to A. marginale infection.
In this study, A. marginale infection of R. microplus was shown to modulate in male salivary glands the expression of genes encoding for putative proteins involved in binding, catalytic/enzymatic activity, transport, DNA/RNA metabolism and structural molecules. Five genes encoding for putative histamine-binding protein (22Hbp), von Willebrand factor (94Will), flagelliform silk protein (100Silk), Kunitz-like protease inhibitor precursor (108Kunz) and proline-rich protein BstNI subfamily 3 precursor (7BstNI3) were confirmed by real-time RT-PCR to be down-regulated in tick salivary glands infected with A. marginale. Functional studies suggested that differentially expressed genes encoding for subolesin, putative von Willebrand factor and flagelliform silk protein could play a role in A. marginale infection and multiplication in ticks. Additionally, for the first time RNAi in cultured BME26 tick cells was used to study A. marginale-tick molecular interactions and suggested that subolesin and flagelliform silk protein may be required by A. marginale for infection and multiplication in these cultured cells. Collectively these data are important toward understanding the molecular mechanisms involved in R. microplus-A. marginale interactions.
A suppression subtractive hybridization (SSH) approach was used to identify genes differentially expressed in R. microplus male salivary gland genes in response to A. marginale infection. Sequences obtained by SSH were used to search for homology/identity to nucleotide and protein databases. Real-time RT-PCR was used to validate differential expression of selected genes in uninfected and A. marginale- infected R. microplus salivary glands. Differential expression patterns of selected genes were also studied in cultured BME26 cells at 6, 24 and 72 hpi by real-time RT-PCR. The genes that proved to have significantly different mRNA levels between uninfected and A. marginale-infected ticks or cultured BME26 tick cells were then selected for functional studies. RNAi was used to characterize the function of selected genes in A. marginale infection in vivo in R. microplus male ticks and in vitro in cultured BME26 tick cells.
Rhipicephalus microplus ticks
The R. microplus ticks used for construction of the SSH libraries originated from Mozambique and were provided by ClinVet International (Pty), Bloemfontein, South Africa. The R. microplus (Susceptible, CENAPA, Mexico strain) ticks used for the RNAi experiments were obtained from laboratory colonies maintained at the University of Tamaulipas, Mexico. Originally, these tick strains were collected from infested cattle in Tapalpa, Jalisco, Mexico. The ticks were maintained on cattle at the tick rearing facilities at the Utrecht Center for Tick-Borne Diseases, Utrecht University or the University of Tamaulipas. Larvae were kept off-host in an incubator at 20°C with 95% relative humidity and 12 hours light: 12 hours dark photoperiod. Cattle were cared for in both The Netherlands and Mexico in accordance with standards specified in the Guide for Care and Use of Laboratory Animals.
Tick cell line
The tick cell line BME26 was derived originally from embryos of R. microplus following the protocol established by Kurrti et al. The cells were maintained in L-15B300 medium  supplemented with 5% heat-inactivated FBS (Gibco/Invitrogen, NY, USA), 10% TPB (Difco, Detroid, MI, USA), 0,1% bovine lipoprotein (ICN, Irvine, CA), 100 units ml-1 penicillin (Gibco/Invitrogen) and 100 μg ml-1 streptomycin (Gibco/Invitrogen) at pH 7.2. Cultures were grown at 34°C in 25 cm2 plastic flasks (Nalge Nunc Int. Rochester, NY) containing 5 ml of the complete medium, which was changed weekly. Monolayers were subcultured when they reached a density of approximately 107 cells/ml and approximately 8 × 105 cells/ml were transferred to the new flask.
A. marginale isolates
The A. marginale isolate used for infection of R. microplus ticks and construction of SSH libraries was obtained from an infected calf in Texas, USA, in 1977. The isolate was subsequently passaged in splenectomized calves and blood samples were collected at the peak parasitemia (40%), prepared as a stabilates with 10% DMSO and stored in 2 ml aliquots in liquid nitrogen. A Brazilian strain of A. marginale with an inclusion appendage (UFMG1) , which was recently established and propagated in IDE8 tick cells, was used to infect BME26 cells .
Uninfected and A. marginale-infected R. microplus male ticks for SSH libraries construction
To obtain A. marginale-infected R. microplus male ticks, eight month-old non-splenectomized, tick-naive Holstein-Friesian calf (No.4280) was infected intravenously with the Texas isolate of A. marginale. R. microplus larvae were then fed on the calf with ascending parasitemia. Approximately 200 partially fed male ticks were collected after 21 days of feeding and the presence of A. marginale infection in salivary glands was confirmed in 20 individually dissected tick salivary glands by msp4 PCR . Salivary glands from the remaining ticks were dissected in cold PBS, pooled and immediately stored in TriReagent (Sigma-Aldrich, Zwijndrecht, The Netherlands) at -80°C. Uninfected ticks were fed in a similar way on an uninfected calf and the male tick salivary glands were dissected and stored as described previously. Cattle and tick feeding studies were conducted in accordance with approval of the Animal Experiments Committee (DEC) of the Faculty of Veterinary Medicine, Utrecht University (DEC No. 0604.0801).
Uninfected and A. marginale-infected cultured BME26 tick cells
The tick cell line BME26 was cultured as described above. Approximately 4 × 106 cells from the passage 72 were plated in 24-well plates (Nunc) and maintained in Anaplasma medium  at 34°C for 3 days. The cells were infected with the 15 day-old culture of A. marginale in BME26 cells. Five milliliters of this suspension were transferred to a plastic tube. The tube was immersed in liquid nitrogen for 5 min for cell disruption and rickettsia releasing, followed by thawing in a water bath at 37°C. Approximately 500 μl of infected cell suspension was inoculated into each well onto uninfected BME26 monolayer. The plate was maintained at 34°C and cells harvested at 6, 24 and 72 hpi from 4 wells for each time point. Uninfected control cells (n = 4 wells) were cultured in the same way but with the addition of 500 μl of Anaplasma medium only and the cells were harvested at the same time points. Cells were centrifuged at 800 g for 10 min and RNA/DNA was isolated using Trizol reagent (Gibco/Invitrogen), as recommended by the manufacturer. The infection of the tick cells by A. marginale was corroborated by msp4 PCR .
Total RNA was isolated from uninfected and A. marginale-infected tick salivary glands using Tri reagent (Sigma-Aldrich) following the manufacturer's protocol. RNA quality was checked by gel electrophoreses to confirm integrity of RNA preparations. Pools of 2 μg of total RNA were made from uninfected and A. marginale-infected salivary glands. The cDNA synthesis was done using the Super SMART PCR cDNA synthesis kit (Clontech-Takara, Saint-Germain-en-Laye, France), a method for producing high quality cDNA from a low amount of starting material. The cDNA was then directly used for PCR select subtraction (Clontech-Takara) based on a technique previously described by Diatchenko et al. [14, 15] and SSH libraries were constructed according to manufacturer's instructions. The double stranded cDNA from both groups (infected and uninfected salivary glands) was Rsa I digested. Part of the digested cDNA was ligated with Adapter 1 and part with the Adapter 2R, and the rest was saved for use as a driver in preparation for hybridization. The forward subtracted library was made by hybridizing adapter ligated cDNA from A. marginale- infected tick salivary glands as the tester in the presence of uninfected tick salivary gland CDNA as the driver. This forward reaction library was designed to produce clones that are overexpressed or up-regulated in infected salivary glands. The reverse library was made in the same way but in this case the adapter ligated cDNA from uninfected tick salivary glands was used as the tester and infected salivary gland cDNA as the driver. The reverse reaction library was designed to produce clones underexpressed or down-regulated in infected salivary glands. In either case the driver cDNA was added in excess during each hybridization to remove common cDNAs by hybrid selection and leaving over expressed and novel tester cDNA to be recovered and cloned. Differentially expressed cDNAs were PCR amplified with Advantage PCR polymerase mix (Clontech-Takara), cloned using pGEM Easy T/A cloning kit (Promega, Madison, WI, USA), transformed into JM109 and plated on LB with ampicillin, X-gal and IPTG. Approximately 300 clones were randomly picked up from each plated library in 96- well plates containing LB medium supplemented with ampicillin and grown overnight. Plasmids were purified using Wizard SV 96 Plasmid DNA purification system (Promega). Plasmid inserts were PCR amplified and PCR products were fully sequenced in an accredited service laboratory (BaseClear, Leiden, The Netherlands) using vector specific primers.
Sequence analysis and database search
Partial sequences were obtained for 279 out of 576 randomly selected sequenced SSH library clones (288 from each forward- and reverse-subtracted libraries). The cDNA Annotation System software (CAS; Bioinformatics and Scientific IT Program (BSIP), Office of Technology Information Systems (OTIS), National Institute of Allergy and Infectious Diseases (NIAID), Bethesda, MD, USA) http://exon.niaid.nih.gov was used for automated sequence clean up, contig assembly, Blast analysis  against multiple sequence databases (nonredundant sequence database and databases of tick-specific sequences. (http://www.vectorbase.org/index.php and BmGI2; http://compbio.dfci.harvard.edu/), identification and locating of signal peptide cleavage sites and gene ontology (GO) assignments. Genes for further analyses were annotated by manual curation.
Nucleotide sequence accession numbers
The gene sequences reported in this paper are deposited in the GenBank Data Library under accession numbers: GO496166-G0496262.
Real-time RT- PCR analysis
The same RNA samples prepared before for SSH from uninfected and A. marginale- infected tick salivary glands were used for real-time RT-PCR. Total RNA was extracted as described above from uninfected and A. marginale-infected cultured BME26 tick cells from quadruplicate cultures at 6, 24 and 72 hpi and used for real-time RT-PCR. Oligonucleotide primers were synthesized based on the sequences determined for candidate differentially expressed genes (Table 3) and used in 25 μl RT-PCR reactions performed using the iScript one step RT-PCR kit with SYBR green (Bio-Rad, Hercules, CA, USA) and a Bio-Rad iQ5 thermal cycler following the manufacturer's protocol. The mRNA levels were normalized against tick β-actin (Genbank accession number AY255624) using the comparative Ct method . The mRNA levels were compared between infected and uninfected tick salivary glands and tick cells by Student's t test (p = 0.05). Total RNA was isolated from bovine blood and RT-PCR reactions were performed using the oligonucleotide primers and conditions described in Table 3. PCR products were electrophoretically separated using 1.5% agarose gel stained with ethidium bromide for visualization.
RNA interference in ticks
Oligonucleotide primers containing T7 promoter sequences at the 5'-end were synthesized for in vitro transcription of dsRNA using the Acess RT-PCR system (Promega, Madison, WI, USA) and the Megascript RNAi kit (Ambion, Austin, TX, USA) following manufacturer's instructions. The dsRNA was purified and quantified by spectophotometry. Newly molted uninfected R. microplus males, were injected with approximately 0.3 μl of dsRNA (5 × 1010 molecules/μl) in the lower right quadrant of the ventral surface of the exoskeleton of the tick . Ticks (35/group) were injected using a Hamilton syringe with a 1 inch, 33 gauge needle. Control ticks were injected with R. microplus subolesin dsRNA (positive control) or the unrelated GIII dsRNA (negative control). The GIII sequence was identified in R. microplus and did not affect tick feeding, mortality and oviposition after RNAi (unpublished results). Ticks were held in a humidity chamber for 3-5 hours and mortality was recorded before the living ticks were allowed to feed in seven separate patches (five test genes and two controls), each one for a different group, glued on the back of a calf naturally infected with A. marginale in Tamaulipas, Mexico (approximately 4% rickettsemia during tick feeding). Ten females were placed in each patch simultaneously with injected males. Unattached ticks were removed 2 days after infestation and attached ticks were allowed to feed for 7 days. All the males were collected and salivary glands were dissected from individual ticks from each group. The salivary glands were stored in RNAlater (Qiagen) until used for DNA and RNA extraction to determine A. marginale infection levels by quantitative msp4 PCR  and to confirm gene knockdown by real-time RT-PCR. The mRNA levels of the target gene and the A. marginale infection were compared between test and control ticks by Student's t-test (p = 0.05). Tick mortality was evaluated as the ratio of the dead male ticks 7 days after dsRNA injection to the total number of attached male ticks feeding on the animal and was compared between test and control groups by χ2-test (α = 0.025).
RNA interference in cultured BME26 tick cells
Approximately 5 × 105 BME26 cells/well were placed in 24-well plates (Nunc). The dsRNA prepared in the way described above for each of the 8 target genes was added to the culture wells (n = 4 for each treatment). Subolesin dsRNA was used as a positive control and control wells received elution buffer only. Each of the treated wells received 10 μl of dsRNA (5 × 1010 molecules per microliter) and 190 μl of L15B Anaplasma medium and was incubated for 24 h. After this period, additional 300 μl of medium were added to each well. After 48 h cells were infected as described above. Three days post infection the cells were harvested from the plate, centrifuged at 800 × g for 10 min and used to extract RNA and DNA with TriReagent (Sigma). A. marginale infection levels were determined by msp4 PCR . Gene expression silencing was confirmed by real-time RT-PCR using sequence-specific primers (Table 3) as described above.
Kocan KM, de la Fuente J, Guglielmone AA, Melendez RD: Antigens and alternatives for control of Anaplasma marginale infection in cattle. Clin Microbiol Rev. 2003, 16 (4): 698-712. 10.1128/CMR.16.4.698-712.2003.
Kocan KM, de la Fuente J, Blouin EF, Garcia-Garcia JC: Anaplasma marginale (Rickettsiales: Anaplasmataceae): recent advances in defining host-pathogen adaptations of a tick-borne rickettsia. Parasitology. 2004, 129 (Suppl): S285-300. 10.1017/S0031182003004700.
Guerrero FD, Nene VM, George JE, Barker SC, Willadsen P: Sequencing a new target genome: the Boophilus microplus (Acari: Ixodidae) genome project. J Med Entomol. 2006, 43 (1): 9-16. 10.1603/0022-2585(2006)043[0009:SANTGT]2.0.CO;2.
de Castro JJ: Sustainable tick and tickborne disease control in livestock improvement in developing countries. Vet Parasitol. 1997, 71 (2-3): 77-97. 10.1016/S0304-4017(97)00033-2.
Kocan KM, Goff WL, Stiller D, Claypool PL, Edwards W, Ewing SA, Hair JA, Barron SJ: Persistence of Anaplasma marginale (Rickettsiales: Anaplasmataceae) in male Dermacentor andersoni (Acari: Ixodidae) transferred successively from infected to susceptible calves. J Med Entomol. 1992, 29 (4): 657-668.
Kocan KM, Stiller D, Goff WL, Claypool PL, Edwards W, Ewing SA, McGuire TC, Hair JA, Barron SJ: Development of Anaplasma marginale in male Dermacentor andersoni transferred from parasitemic to susceptible cattle. Am J Vet Res. 1992, 53 (4): 499-507.
Chevillon C, Koffi BB, Barre N, Durand P, Arnathau C, de Meeus T: Direct and indirect inferences on parasite mating and gene transmission patterns. Pangamy in the cattle tick Rhipicephalus (Boophilus) microplus. Infect Genet Evol. 2007, 7 (2): 298-304. 10.1016/j.meegid.2006.11.007.
Macaluso KR, Mulenga A, Simser JA, Azad AF: Differential expression of genes in uninfected and rickettsia-infected Dermacentor variabilis ticks as assessed by differential-display PCR. Infect Immun. 2003, 71 (11): 6165-6170. 10.1128/IAI.71.11.6165-6170.2003.
Nene V, Lee D, Kang'a S, Skilton R, Shah T, de Villiers E, Mwaura S, Taylor D, Quackenbush J, Bishop R: Genes transcribed in the salivary glands of female Rhipicephalus appendiculatus ticks infected with Theileria parva. Insect Biochem Mol Biol. 2004, 34 (10): 1117-1128. 10.1016/j.ibmb.2004.07.002.
Rudenko N, Golovchenko M, Edwards MJ, Grubhoffer L: Differential expression of Ixodes ricinus tick genes induced by blood feeding or Borrelia burgdorferi infection. J Med Entomol. 2005, 42 (1): 36-41. 10.1603/0022-2585(2005)042[0036:DEOIRT]2.0.CO;2.
de la Fuente J, Kocan KM, Almazan C, Blouin EF: RNA interference for the study and genetic manipulation of ticks. Trends Parasitol. 2007, 23 (9): 427-433. 10.1016/j.pt.2007.07.002.
Bell-Sakyi L, Zweygarth E, Blouin EF, Gould EA, Jongejan F: Tick cell lines: tools for tick and tick-borne disease research. Trends Parasitol. 2007, 23 (9): 450-457. 10.1016/j.pt.2007.07.009.
Esteves E, Bastos CV, Zivkovic Z, de La Fuente J, Kocan K, Blouin E, Ribeiro MF, Passos LM, Daffre S: Propagation of a Brazilian isolate of Anaplasma marginale with appendage in a tick cell line (BME26) derived from Rhipicephalus (Boophilus) microplus. Vet Parasitol. 2008, 161 (1-2): 150-3. 10.1016/j.vetpar.2008.12.006.
Diatchenko L, Lukyanov S, Lau YF, Siebert PD: Suppression subtractive hybridization: a versatile method for identifying differentially expressed genes. Methods Enzymol. 1999, 303: 349-380. full_text.
Diatchenko L, Lau YF, Campbell AP, Chenchik A, Moqadam F, Huang B, Lukyanov S, Lukyanov K, Gurskaya N, Sverdlov ED: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Natl Acad Sci USA. 1996, 93 (12): 6025-6030. 10.1073/pnas.93.12.6025.
Pal U, Li X, Wang T, Montgomery RR, Ramamoorthi N, Desilva AM, Bao F, Yang X, Pypaert M, Pradhan D: TROSPA, an Ixodes scapularis receptor for Borrelia burgdorferi. Cell. 2004, 119 (4): 457-468. 10.1016/j.cell.2004.10.027.
Ramamoorthi N, Narasimhan S, Pal U, Bao F, Yang XF, Fish D, Anguita J, Norgard MV, Kantor FS, Anderson JF: The Lyme disease agent exploits a tick protein to infect the mammalian host. Nature. 2005, 436 (7050): 573-577. 10.1038/nature03812.
Sukumaran B, Narasimhan S, Anderson JF, DePonte K, Marcantonio N, Krishnan MN, Fish D, Telford SR, Kantor FS, Fikrig E: An Ixodes scapularis protein required for survival of Anaplasma phagocytophilum in tick salivary glands. J Exp Med. 2006, 203 (6): 1507-1517. 10.1084/jem.20060208.
de la Fuente J, Almazan C, Blouin EF, Naranjo V, Kocan KM: Reduction of tick infections with Anaplasma marginale and A. phagocytophilum by targeting the tick protective antigen subolesin. Parasitol Res. 2006, 100 (1): 85-91. 10.1007/s00436-006-0244-6.
de la Fuente J, Lew A, Lutz H, Meli ML, Hofmann-Lehmann R, Shkap V, Molad T, Mangold AJ, Almazan C, Naranjo V: Genetic diversity of anaplasma species major surface proteins and implications for anaplasmosis serodiagnosis and vaccine development. Anim Health Res Rev. 2005, 6 (1): 75-89. 10.1079/AHR2005104.
Brayton KA, Kappmeyer LS, Herndon DR, Dark MJ, Tibbals DL, Palmer GH, McGuire TC, Knowles DP: Complete genome sequencing of Anaplasma marginale reveals that the surface is skewed to two superfamilies of outer membrane proteins. Proc Natl Acad Sci USA. 2005, 102 (3): 844-849. 10.1073/pnas.0406656102.
de la Fuente J, Kocan KM, Almazan C, Blouin EF: Targeting the tick-pathogen interface for novel control strategies. Front Biosci. 2008, 13: 6947-6956. 10.2741/3201.
Bishop R, Lambson B, Wells C, Pandit P, Osaso J, Nkonge C, Morzaria S, Musoke A, Nene V: 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. 2002, 32 (7): 833-842. 10.1016/S0020-7519(02)00027-9.
Paesen GC, Adams PL, Harlos K, Nuttall PA, Stuart DI: Tick histamine-binding proteins: isolation, cloning, and three-dimensional structure. Mol Cell. 1999, 3 (5): 661-671. 10.1016/S1097-2765(00)80359-7.
Labuda M, Trimnell AR, Lickova M, Kazimirova M, Davies GM, Lissina O, Hails RS, Nuttall PA: An antivector vaccine protects against a lethal vector-borne pathogen. PLoS Pathog. 2006, 2 (4): e27-10.1371/journal.ppat.0020027.
Hayashi CY, Lewis RV: Molecular architecture and evolution of a modular spider silk protein gene. Science. 2000, 287 (5457): 1477-1479. 10.1126/science.287.5457.1477.
Alarcon-Chaidez FJ, Sun J, Wikel SK: Transcriptome analysis of the salivary glands of Dermacentor andersoni Stiles (Acari: Ixodidae). Insect Biochem Mol Biol. 2007, 37 (1): 48-71. 10.1016/j.ibmb.2006.10.002.
Santos IK, Valenzuela JG, Ribeiro JM, de Castro M, Costa JN, Costa AM, da Silva ER, Neto OB, Rocha C, Daffre S: Gene discovery in Boophilus microplus, the cattle tick: the transcriptomes of ovaries, salivary glands, and hemocytes. Ann N Y Acad Sci. 2004, 1026: 242-246. 10.1196/annals.1307.037.
Mulenga A, Blandon M, Khumthong R: The molecular basis of the Amblyomma americanum tick attachment phase. Exp Appl Acarol. 2007, 41 (4): 267-287. 10.1007/s10493-007-9064-3.
Donohue KV, Khalil SM, Mitchell RD, Sonenshine DE, Roe RM: Molecular characterization of the major hemelipoglycoprotein in ixodid ticks. Insect Mol Biol. 2008, 17 (3): 197-208. 10.1111/j.1365-2583.2008.00794.x.
Almazan C, Kocan KM, Bergman DK, Garcia-Garcia JC, Blouin EF, de la Fuente J: Identification of protective antigens for the control of Ixodes scapularis infestations using cDNA expression library immunization. Vaccine. 2003, 21 (13-14): 1492-1501. 10.1016/S0264-410X(02)00683-7.
de la Fuente J, Almazan C, Blas-Machado U, Naranjo V, Mangold AJ, Blouin EF, Gortazar C, Kocan KM: The tick protective antigen, 4D8, is a conserved protein involved in modulation of tick blood ingestion and reproduction. Vaccine. 2006, 24 (19): 4082-4095. 10.1016/j.vaccine.2006.02.046.
Almazan C, Blas-Machado U, Kocan KM, Yoshioka JH, Blouin EF, Mangold AJ, de la Fuente J: Characterization of three Ixodes scapularis cDNAs protective against tick infestations. Vaccine. 2005, 23 (35): 4403-4416. 10.1016/j.vaccine.2005.04.012.
Almazan C, Kocan KM, Blouin EF, de la Fuente J: Vaccination with recombinant tick antigens for the control of Ixodes scapularis adult infestations. Vaccine. 2005, 23 (46-47): 5294-5298. 10.1016/j.vaccine.2005.08.004.
de la Fuente J, Almazan C, Naranjo V, Blouin EF, Meyer JM, Kocan KM: Autocidal control of ticks by silencing of a single gene by RNA interference. Biochem Biophys Res Commun. 2006, 344 (1): 332-338. 10.1016/j.bbrc.2006.03.109.
Nijhof AM, Taoufik A, de la Fuente J, Kocan KM, de Vries E, Jongejan F: Gene silencing of the tick protective antigens, Bm86, Bm91 and subolesin, in the one-host tick Boophilus microplus by RNA interference. Int J Parasitol. 2007, 37 (6): 653-662. 10.1016/j.ijpara.2006.11.005.
Kocan KM, Manzano-Roman R, de la Fuente J: Transovarial silencing of the subolesin gene in three-host ixodid tick species after injection of replete females with subolesin dsRNA. Parasitol Res. 2007, 100 (6): 1411-1415. 10.1007/s00436-007-0483-1.
de la Fuente J, Maritz-Olivier C, Naranjo V, Ayoubi P, Nijhof AM, Almazan C, Canales M, Perez de la Lastra JM, Galindo RC, Blouin EF: Evidence of the role of tick subolesin in gene expression. BMC Genomics. 2008, 9 (1): 372-10.1186/1471-2164-9-372.
Galindo RC, Doncel-Perez E, Zivkovic Z, Naranjo V, Gortazar C, Mangold AJ, Martin-Hernando MP, Kocan KM, de la Fuente J: Tick subolesin is an ortholog of the akirins described in insects and vertebrates. Dev Comp Immunol. 2009, 33 (4): 612-617. 10.1016/j.dci.2008.11.002.
de la Fuente J, Blouin EF, Manzano-Roman R, Naranjo V, Almazan C, Perez de la Lastra JM, Zivkovic Z, Massung RF, Jongejan F, Kocan KM: Differential expression of the tick protective antigen subolesin in anaplasma marginale- and A. phagocytophilum-infected host cells. Ann N Y Acad Sci. 2008, 1149: 27-35. 10.1196/annals.1428.056.
Borghesi LA, Lynes MA: Stress proteins as agents of immunological change: some lessons from metallothionein. Cell Stress Chaperones. 1996, 1 (2): 99-108. 10.1379/1466-1268(1996)001<0099:SPAAOI>2.3.CO;2.
Dallinger R: Metallothionein research in terrestrial invertebrates: synopsis and perspectives. Comp Biochem Physiol C Pharmacol Toxicol Endocrinol. 1996, 113 (2): 125-133. 10.1016/0742-8413(95)02078-0.
Dutsch-Wicherek M, Sikora J, Tomaszewska R: The possible biological role of metallothionein in apoptosis. Front Biosci. 2008, 13: 4029-4038. 10.2741/2991.
Mannick EE, Schurr JR, Zapata A, Lentz JJ, Gastanaduy M, Cote RL, Delgado A, Correa P, Correa H: Gene expression in gastric biopsies from patients infected with Helicobacter pylori. Scand J Gastroenterol. 2004, 39 (12): 1192-1200. 10.1080/00365520410003588.
Wiese L, Kurtzhals JA, Penkowa M: Neuronal apoptosis, metallothionein expression and proinflammatory responses during cerebral malaria in mice. Exp Neurol. 2006, 200 (1): 216-226.
Baerwald MR, Welsh AB, Hedrick RP, May B: Discovery of genes implicated in whirling disease infection and resistance in rainbow trout using genome-wide expression profiling. BMC Genomics. 2008, 9: 37-10.1186/1471-2164-9-37.
Hirano T, Murakami M, Fukada T, Nishida K, Yamasaki S, Suzuki T: Roles of zinc and zinc signaling in immunity: zinc as an intracellular signaling molecule. Adv Immunol. 2008, 97: 149-176. full_text.
Uzarowska A, Dionisio G, Sarholz B, Piepho HP, Xu M, Ingvardsen CR, Wenzel G, Lubberstedt T: Validation of candidate genes putatively associated with resistance to SCMV and MDMV in maize (Zea mays L.) by expression profiling. BMC Plant Biol. 2009, 9: 15-10.1186/1471-2229-9-15.
Saldivar L, Guerrero FD, Miller RJ, Bendele KG, Gondro C, Brayton KA: Microarray analysis of acaricide-inducible gene expression in the southern cattle tick, Rhipicephalus (Boophilus) microplus. Insect Mol Biol. 2008, 17 (6): 597-606. 10.1111/j.1365-2583.2008.00831.x.
Kurtti TJ, Munderloh UG, Ahlstrand GG, Johnson RC: Borrelia burgdorferi in tick cell culture: growth and cellular adherence. J Med Entomol. 1988, 25 (4): 256-261.
Munderloh UG, Jauron SD, Fingerle V, Leitritz L, Hayes SF, Hautman JM, Nelson CM, Huberty BW, Kurtti TJ, Ahlstrand GG: Invasion and intracellular development of the human granulocytic ehrlichiosis agent in tick cell culture. J Clin Microbiol. 1999, 37 (8): 2518-2524.
Ribeiro MF, Passos LM, Guimaraes AM: Ultrastructure of Anaplasma marginale with an inclusion appendage, isolated in Minas Gerais State, Brazil. Vet Parasitol. 1997, 70 (4): 271-277. 10.1016/S0304-4017(97)00004-6.
de la Fuente J, Garcia-Garcia JC, Blouin EF, McEwen BR, Clawson D, Kocan KM: Major surface protein 1a effects tick infection and transmission of Anaplasma marginale. Int J Parasitol. 2001, 31 (14): 1705-1714. 10.1016/S0020-7519(01)00287-9.
Munderloh UG, Blouin EF, Kocan KM, Ge NL, Edwards WL, Kurtti TJ: Establishment of the tick (Acari:Ixodidae)-borne cattle pathogen Anaplasma marginale (Rickettsiales:Anaplasmataceae) in tick cell culture. J Med Entomol. 1996, 33 (4): 656-664.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Schefe JH, Lehmann KE, Buschmann IR, Unger T, Funke-Kaiser H: Quantitative real-time RT-PCR data analysis: current concepts and the novel "gene expression's CT difference" formula. J Mol Med. 2006, 84 (11): 901-910. 10.1007/s00109-006-0097-6.
This research was supported by grants from the European Community, INCO-DEV program (project no. 003713), entitled 'Epidemiology and new generation vaccines for Ehrlichia and Anaplasma infections of ruminants', the Ministerio de Ciencia e Innovación, Spain (project BFU2008-01244/BMC), the CSIC intramural project 200830I249 to JF, FAPESP (Fundação de Amparo a Pesquisa do Estado de São Paulo, Brazil), and CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico, Brazil) and was facilitated through the Integrated Consortium on Ticks and Tick-borne Diseases (ICTTD-3), financed by the International Cooperation Program of the European Union, coordination action project No. 510561. We are grateful to Biankynet Gonzalez for the technical assistance.
ZZ prepared SSH libraries, participated in sequence analyses, did the molecular work and drafted the manuscript. EE did the cell culture work and helped with drafting the manuscript. CA helped with the tick RNAi experiments and revision of the manuscript. AN helped with tick infection experiment and revision of the manuscript. SD and KK helped with drafting and revision of the manuscript. FJ and JF participated in design and coordination of the study and helped with drafting and revision of the manuscript. JF contributed to acquisition of molecular biology data. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Genes identified by SSH as differentially expressed in A. marginale-infected R. microplus male salivary glands.(DOC 189 KB)
Additional file 2: Differentially expressed genes selected based on their putative role during A. marginale infection to validate SSH results by real-time RT-PCR.(DOC 40 KB)
About this article
Cite this article
Zivkovic, Z., Esteves, E., Almazán, C. et al. Differential expression of genes in salivary glands of male Rhipicephalus (Boophilus)microplus in response to infection with Anaplasma marginale. BMC Genomics 11, 186 (2010). https://doi.org/10.1186/1471-2164-11-186