Skip to main content

Advertisement

Identification of a novel β-adrenergic octopamine receptor-like gene (βAOR-like) and increased ATP-binding cassette B10 (ABCB10) expression in a Rhipicephalus microplus cell line derived from acaricide-resistant ticks

Article metrics

Abstract

Background

The cattle tick Rhipicephalus (Boophilus) microplus is an economically important parasite of livestock. Effective control of ticks using acaricides is threatened by the emergence of resistance to many existing compounds. Several continuous R. microplus cell lines have been established and provide an under-utilised resource for studies into acaricide targets and potential genetic mutations associated with resistance. As a first step to genetic studies using these resources, this study aimed to determine the presence or absence of two genes and their transcripts that have been linked with acaricide function in cattle ticks: β-adrenergic octopamine receptor (βAOR, associated with amitraz resistance) and ATP-binding cassette B10 (ABCB10, associated with macrocyclic lactone resistance) in six R. microplus cell lines, five other Rhipicephalus spp. cell lines and three cell lines representing other tick genera (Amblyomma variegatum, Ixodes ricinus and Hyalomma anatolicum).

Methods

End-point polymerase chain reaction (PCR) was used for detection of the βAOR gene and transcripts in DNA and RNA extracted from the tick cell lines, followed by capillary sequencing of amplicons. Quantitative real-time PCR (qPCR) was performed to determine the levels of expression of ABCB10.

Results

βAOR gene expression was detected in all Rhipicephalus spp. cell lines. We observed a second amplicon of approximately 220 bp for the βAOR gene in the R. microplus cell line BME/CTVM6, derived from acaricide-resistant ticks. Sequencing of this transcript variant identified a 36 bp insertion in the βAOR gene, leading to a 12-amino acid insertion (LLKTLALVTIIS) in the first transmembrane domain of the protein. In addition, nine synonymous SNPs were also discovered in R. appendiculatus, R. evertsi and R. sanguineus cell lines. Some of these SNPs appear to be unique to each species, providing potential tools for differentiating the tick species. The BME/CTVM6 cell line had significantly higher ABCB10 (P = 0.002) expression than the other R. microplus cell lines.

Conclusions

The present study has identified a new βAOR gene and demonstrated a higher ABCB10 expression level in the BME/CTVM6 cell line, indicating that tick cell lines provide a useful experimental tool for acaricide resistance studies and further elucidation of tick genetics.

Background

Rhipicephalus (Boophilus) microplus, also known as the cattle tick or southern cattle tick, is a parasite of major global economic importance. Sales-Junior et al. [1] estimated that more than 75 % of the cattle population worldwide was infested with R. microplus, making it one of the most prevalent of bovine parasites. Rhipicephalus microplus causes significant economic losses; direct losses include damage to the hide from tick attachment, anaemia and in severe cases death [2]. More importantly, cattle ticks act as vectors of the protozoan and bacterial pathogens Babesia bigemina, Babesia bovis and Anaplasma marginale, all of which cause production losses and can be fatal [3]. Rhipicephalus appendiculatus is an important tick of cattle in eastern Africa, where it transmits Theileria parva [4]. Rhipicephalus sanguineus is a globally widespread tick that is commonly found on dogs [5]. Amblyomma variegatum infests cattle throughout sub-Saharan Africa and transmits Ehrlichia ruminantium, the causal agent of heartwater [6]. Ixodes ricinus is the most widely distributed tick species in northern, western and central Europe and is of human health significance as the vector of Lyme disease (Borrelia burgdorferi) as well as being an important vector of several diseases of cattle and sheep [4, 7]. Hyalomma anatolicum infests cattle in North Africa, the Middle East and Asia and transmits Theileria annulata, of economic significance in the region [4].

The short, one-host life-cycle of R. microplus acts together with high fecundity to result in rapid selection for development of acaricide resistance. There are several mechanisms by which acaricide resistance has evolved; those most commonly noted are target site modification and metabolic resistance [810]. Target site modification has been most clearly demonstrated in the para sodium channel, the target of synthetic pyrethroids (SP) [11, 12]. Metabolic resistance occurs as a result of modified ability to detoxify and/or sequester pesticides, mediated by families of enzymes known for detoxification of xenobiotics, such as the cytochrome P450s, esterases, glutathione S-transferases and ATP-binding cassette (ABC) transporters [10]. Penetration resistance is a little-studied mechanism which might also occur in ticks [9].

Octopamine receptors have been shown from studies in insects to be the target for formamidines such as amitraz, leading to the hypothesis that mutations in octopamine receptors would confer resistance to formamidines [13]. There are three main types of octopamine receptors in insects: α-adrenergic-like, β-adrenergic-like and tyramine/octopamine (tyraminergic) [14]. Our previous work reported an association between polymorphism in the β-adrenergic octopamine receptor (βAOR) and resistance to amitraz in R. microplus ticks in Australia [15]. A non-synonymous A → T at position 181 resulted in substitution of isoleucine to phenylalanine (I61F) in the first trans-membrane domain. It is therefore expected that any resistance conferred by mutation in βAOR would be restricted to compounds that would bind to this receptor type.

ATP-binding cassette (ABC) transporters are a superfamily of transmembrane proteins, consisting of up to 8 subfamilies [16] found in all organisms, including R. microplus [17]. They transport substrates across cell membranes as a monomeric or dimeric channel, often against concentration gradients [16]. An association between the ABC transporter and ivermectin (macrocyclic lactone) resistance amongst cattle ticks was first reported by Pohl et al. [17]. These authors detected a 3–6-fold increase in gene expression of R. microplus ABC subtype B10 (RmABCB10) in two different ivermectin-resistant strains. The importance of increased ABCB10 expression through metabolic processes was subsequently demonstrated in an ivermectin resistance-induced R. microplus cell line [18] but it is uncertain if ABCB10 plays a role in resistance to other classes of acaricides.

Over 40 continuous tick cell lines derived from species of worldwide economic importance are maintained in the Tick Cell Biobank at The Pirbright Institute [19]. The establishment of tick cell lines is slow and challenging, and the cultures are composed of cells obtained from hundreds or thousands of embryos or multiple developing nymphs or adults. Further, the cultures are composed of multiple cell types and the cells can be of variable ploidy [20], with proportions of the different cell populations in a cell line varying over time. Therefore, although a small number of acaricide resistance selection studies have been successfully conducted using continuous tick cell lines [21, 22], because they are non-clonal and phenotypically and genotypically heterogeneous it cannot be assumed that any continuous tick cell line is a reliable and consistent model for examining acaricide resistance.

The purpose of the present study was to assess the suitability of six R. microplus continuous cell lines for future genomic and transcriptomic studies on resistance to acaricides. In this study we wished to confirm that two genes that have been strongly associated with resistance (one that operates via target site modification and one that operates via metabolic resistance mechanisms) are expressed by tick cells in vitro, to identify new variants of these genes if any were present, and to determine whether there is variation in their expression level. These two genes were selected for this study because (i) amitraz and macrocyclic lactones (including ivermectin) are globally the most widely used acaricides at present, (ii) we have well-established protocols for working with these genes, and (iii) they represent two of the main mechanisms of resistance-detoxification (ABCB10) and target-site insensitivity (βAOR). The first objective was to partially sequence the βAOR gene for single nucleotide polymorphisms (SNPs) that were reported previously by Corley et al. [15]. The second objective was to determine the expression level of ABCB10 in the Rhipicephalus spp., in particular the R. microplus, cell lines. The study design did not allow a robust determination of the relationship between genotype or expression level and the resistance status of the population of ticks from which the continuous cell lines were developed.

Methods

Cell lines

Fourteen tick cell lines, derived from R. microplus, three other Rhipicephalus spp., A. variegatum, H. anatolicum and I. ricinus, were provided by the Tick Cell Biobank at The Pirbright Institute. The cell lines, passage level(s) tested, species and instar of origin, geographic origin and acaricide resistance status of the parent ticks (if known) are listed in Table 1. All tick cell lines were maintained in 2 ml complete culture medium supplemented with foetal bovine serum in flat-sided culture tubes (Nunc) [23]. Prior to nucleic acid extraction, the tubes were completely filled with culture medium and transferred to the University of Glasgow at room temperature. On arrival, the entire contents of each tube of cells was centrifuged at 300 rcf for 5 min, most of the supernatant medium was discarded, leaving only 3 ml of medium to re-suspend the cells. The cells were then divided into 3 aliquots in 1.5 ml microtubes and centrifuged at 900 rcf for 5 min. The supernatants were then discarded and the cell pellets frozen at −80 °C for nucleic acid extraction.

Table 1 Tick cell lines used in the study

Extraction of nucleic acids

Total RNA was extracted from cell pellets using the Qiagen miRNeasy Mini Kit (#217004) with on-column DNase digestion following the manufacturer’s protocol for extraction from cells. Genomic DNA (gDNA) was extracted from cell pellets using the Qiagen QiaAMP Mini Kit (#51304) following the manufacturer’s protocol. The RNA and DNA concentrations were measured spectrophotometrically using a Nanodrop ND-1000 and quality was confirmed at absorbance 260/280 nm to be ≥ 1.95.

Reverse transcription

First strand complementary DNA (cDNA) was synthesised using 800 ng total RNA with TaqMan Reverse Transcription Reagents (Life Technologies #N8080234) in a 20 μl reaction volume according to the manufacturer’s instructions. This was performed in a Life Technologies SimpliAmp™ Thermal Cycler using the recommended cycling conditions. The cDNA was then diluted with 30 μl of RNase-free water, before any subsequent experiments.

Polymerase chain reaction (PCR)

Each reaction was set up using a HotStar Taq Plus DNA Polymerase kit (Qiagen #203603) in a 20 μl volume containing final concentrations of 1× buffer (with 1.5 mM MgCl2), dNTPs (0.2 mM), forward primer (0.5 μM), reverse primer (0.5 μM), Taq polymerase (5 U/μl), RNase-free H2O and 2 μl of cDNA or 20 ng of gDNA. The primers and PCR conditions used are shown in Table 2. A non-template control (NTC) was included in every PCR. End-point PCR for ABCB10 and the housekeeping gene β-actin (ActB) using cDNA of all cell lines was done to check for presence of transcript before proceeding to quantitative PCR (qPCR). Thermal cycling conditions were initial activation at 95 °C for 5 min; 35 cycles of 95 °C for 1 min, annealing temperature (Table 2) for 1 min at 72 °C; and final extension at 60 °C for 30 min.

Table 2 Primers and PCR conditions

Sequencing of βAOR

All detected amplicons from gDNA and cDNA, with the exception of gDNA from BME/CTVM5, were cleaned for sequencing using Agencourt AMPure XP (Beckman Coulter #A63880) according to the manufacturer’s instructions and reconstituted in 40 μl of RNase-free water. For cell line BME/CTVM5, only the amplicon from cDNA was sequenced because the other two genomic fragments that were amplified were not transcribed. Cleaned amplicons were checked by 1.5 % agarose gel electrophoresis before sequencing. Sequencing reactions were set up using HPLC-purified primers and a BigDye® Terminator v3.1 Cycle Sequencing Kit (Life Technologies #4337455) in a 20 μl volume containing a final concentration of 3.2 pmol of forward or reverse primer, 0.5 μl of BigDye reaction mix, 3.5 μl of 5× sequencing buffer and 10 μl of cleaned amplicon.

Sequencing reactions were cleaned using an Agencourt CleanSEQ Dye Terminator Removal Kit (#A29151) according to the manufacturer’s instructions, and trace signal was detected through capillary electrophoresis on a 3130XL Genetic Analyzer. CLC Genomics Workbench 7.5.2 was used to manually check the trace results, generate alignments and the amino acid sequence used for domain prediction and 2-dimensional representation of protein structure.

Real-time quantitative PCR (qPCR)

Real-time quantitative measurements of ABCB10 and ActB expression were determined in a 10 μl reaction volume containing 1X Brilliant III Ultra-Fast SYBR Green Q-PCR Master Mix (Agilent Technologies #600882), 0.03 μM of ROX reference dye, 0.5 μM forward and reverse primers (Table 2) and 2 μl of cDNA. A non-template control (NTC) was included in every run. Each sample was set up in technical duplicates. The qPCR was performed using the real-time system Stratagene MX3000p and cycling conditions were initial denaturation at 95 °C for 5 min; 35 cycles of 95 °C for 1 min, 60 °C for 1 min and 72 °C for 1 min; final extension at 60 °C for 30 min; and finally a melting curve for confirmation of single amplicon. The cycle threshold (Ct) values were exported for analysis in Excel. The experiment was replicated twice with all nine cell lines investigated and twice more with only R. microplus cell lines.

Statistical analysis

Expression of the ABCB10 gene was normalised to ActB expression and the mean ratio for technical duplicates calculated based on the Pfaffl method [24], using BmVIII-SCC as the calibrator cell line. Differences between the cell lines were tested using one-way analysis of variance (ANOVA) with Tukey’s post-hoc test on Minitab v17.

Results

Detection of the βAOR gene and transcripts in tick cell lines

The βAOR gene and transcripts from the six R. microplus cell lines were detected by PCR in gDNA and cDNA respectively (Fig. 1a). An amplicon of the expected size, 183 bp, was detected in gDNA and cDNA of BME/CTVM2, BME/CTVM5 (passage 67), BME/CTVM23, BME/CTVM30 and BmVIII-SCC. An alternative amplicon of ~220 bp was detected in gDNA and cDNA of BME/CTVM6 (mixture of cells from two sublines at passages 211 and 233, separate since passage 148). There were also two additional amplicons of ~245 bp and ~220 bp in gDNA of BME/CTVM5 (passage 67); these were not detected in cDNA from this cell line (Fig. 1a). To determine whether or not the alternative/additional amplicons detected in BME/CTVM5 and BME/CTVM6 cells were consistently present in these cell lines at different passage levels, the earliest available passages of BME/CTVM5 (passage 8) and BME/CTVM6 (passage 32) and from the two high-passage sublines of BME/CTVM6 (passages 221 and 243 extracted separately) were tested. In all cases the respective amplicons detected in gDNA from the first high passages tested were also present in the respective low passages and the two sublines of BME/CTVM6 (Additional file 1: Figure S1).

Fig. 1
figure1

βAOR PCR in Rhipicephalus cell lines. a Detection of βAOR in genomic DNA (gDNA) and complementary DNA (cDNA) of Rhipicephalus microplus cell lines. The expected amplicon of 183 bp was detected in gDNA and cDNA of cell lines BME/CTVM2 (2), BME/CTVM5 (5), BME/CTVM23 (23), BME/CTVM30 (30) and BmVIII-SCC (SCC). BME/CTVM6 (6) gave an alternative amplicon of ~220 bp from both gDNA and cDNA. Additional amplicons of ~245 bp and ~220 bp were also detected in gDNA of BME/CTVM5. b Detection of βAOR in gDNA and cDNA of other tick cell lines. The expected amplicon was detected in gDNA and cDNA of cell lines RAN/CTVM3 (RAN), REN/CTVM32 (REN), RSE/PILS35 (RSE) and RA243, and only in gDNA of RML-RSE (RML). No βAOR gene or transcript was detected in cell lines AVL/CTVM13 (AVL), HAE/CTVM9 (HAE) or IRE/CTVM19 (IRE). Abbreviations: M, Marker; NTC, Non-template control; −ve, reverse transcription negative control

Amplicons of the expected 183 bp size were detected in gDNA and cDNA of RAN/CTVM3, RA243, REN/CTVM32 and RSE/PILS35, but only in gDNA of RML-RSE (Fig. 1b). The βAOR gene and transcripts were not detected in cell lines AVL/CTVM13, HAE/CTVM9 and IRE/CTVM19 (Fig. 1b) using the primers in Table 2 in two attempts. Consequently, these three cell lines were excluded from further investigation.

βAOR gene variations

Sequencing and comparison with the published sequence (GenBank JN974909) of the βAOR gene of R. microplus (strain NRFS) revealed that the transcript variant in BME/CTVM6 and the 220 bp gene variant in BME/CTVM5 contained a 36 bp insertion at position 190 from the start codon (Fig. 2a), not present in any of the other tick cell lines. This resulted in a 12-aa insertion LLKTLALVTIIS (Fig. 2b), which appeared to be an additional leucine (L) followed by a duplicated LKTLALVTIIS. Using the TMHMM prediction and TMRPres2D diagram representation, this insertion was predicted to occur within the first transmembrane domain leaving an extracellular domain comprising 66-aa instead of 54-aa (Fig. 2c). A discontiguous megablast search did not find any match for this insertion in any other species.

Fig. 2
figure2

Alignment of Rhipicehalus βAOR sequences and predicted structure of the βAOR protein. a Alignment of βAOR genomic sequences in Rhipicephalus spp. cell lines with published R. microplus strain NRFS (JN974909). Variations are highlighted in pink. b Alignment of βAOR amino acid sequences of published R. microplus strain NRFS (JN974909) and BME/CTVM6. The 12-aa insertion consists of an additional L (highlighted in red), followed by what appears to be duplication of LKTLALVTIIS (highlighted in blue). Green lines indicate the first transmembrane domain. c 2-dimensional representation of the R. microplus strain NRFS βAOR protein (left diagram) and the consequence of the 12-aa insertion (highlighted in red and blue) in the BME/CTVM6 βAOR protein (right diagram). The red line in the left diagram indicates the position of insertion. Membrane-spanning domains were predicted by the TMHMM Server at the Center for Biological Sequence Analysis, Technical University of Denmark, DTU (http://www.cbs.dtu.dk/services/TMHMM/) and 2 dimensional representation by TMRPres2D [31]

Attempts to sequence the additional 245 bp amplicon from BME/CTVM5 were unsuccessful despite successful gel extraction of the product (Additional file 2: Figure S2). There were two regions of ambiguity within the amplicon between regions of homology of the sequences, suggesting that there was more than one binding site for the primers. PCR amplification of the purified 245 bp amplicon produced the same three amplicons seen in Fig. 1a and Additional file 1: Figure S1. This confirms that there was more than one binding site for the primers.

Nine synonymous SNPs were also discovered in R. appendiculatus, R. evertsi and R. sanguineus cell lines (Table 3). Three SNPs (C141T, A225G and A231C) were present in all five R. appendiculatus, R. evertsi and R. sanguineus cell lines, while C189G was present only in the R. evertsi and R. sanguineus lines. A195C and G204C were unique to R. appendiculatus cell lines while G204A, C216G, G222A and C228G were unique to R. sanguineus cell lines.

Table 3 Variations in Rhipicephalus bAOR gene

These sequences have been submitted to GenBank under the following accession numbers: KU836738 (BmVIII-SCC), KU836739 (BME/CTVM2), KU836740 (BME/CTVM5), KU836741 (BME/CTVM6), KU836742 (BME/CTVM23), KU836743 (BME/CTVM30), KU836744 (RAN/CTVM3), KU836745 (RA243), KU836746 (REN/CTVM32), KU836747 (RML-RSE) and KU836748 (RSE/PILS35).

Expression of ABCB10

After normalisation to ActB, expression of the ABCB10 gene in BME/CTVM6 was significantly higher (ANOVA F (9,10) = 8.1, P = 0.002), by approximately seven-fold, compared to expression in all the other cell lines (Fig. 3). There was no significant difference in the ABCB10 expression level between any of the other cell lines.

Fig. 3
figure3

ABCB10 expression levels in Rhipicephalus spp. cell lines. ABCB10 expression in Rhipicephalus spp. cell lines BME/CTVM2 (2), BME/CTVM5 (5), BME/CTVM6 (6), BME/CTVM23 (23), BME/CTVM30 (30), RAN/CTVM3 (RAN), REN/CTVM32 (REN), RSE/PILS35 (RSE) and RA243 using BmVIII-SCC (SCC) as the calibrator cell line. Expression levels are shown numerically above each bar in the graph. Expression of ABCB10 was significantly higher in BME/CTVM6 (P = 0.002) compared to all other cell lines. There was no significant difference in ABCB10 expression level between all other cell lines. Gene expression is represented as mean ratio of ABCB10 to ActB, calculated according to the Pfaffl method [24], compared using one-way analysis of variance (ANOVA) with Tukey’s post-hoc test

Discussion

Acaricide resistance in R. microplus has been associated with mutations in several genes including the βAOR gene [15]. We demonstrated expression of this gene in six tick cell lines derived from four strains of R. microplus with different acaricide resistance status. Corley et al. [15] identified two non-synonymous SNPs for two highly conserved adjacent isoleucine residues of the first transmembrane domain of the β-adrenergic-like octopamine receptor (βAOR) in R. microplus. The first identified SNP resulting in an amino acid substitution of I61F was found in the Ultimo tick strain from Australia, which was resistant to amitraz, synthetic pyrethroids and organophosphates. Our study used the same primers and detected the βAOR gene and transcripts in all the R. microplus cell lines that retain the highly conserved isoleucine at amino acids 61 and 62. However, we discovered an alternative βAOR gene and transcript with a 36 bp insertion in the cell line BME/CTVM6, derived from R. microplus (Paso Ancho strain) resistant to multiple acaricides. This insertion introduced 12 aa into the first transmembrane domain of the seven-transmembrane protein structure, resulting in a longer extracellular domain. The βAOR protein is the predicted target for formamidines [25]; hence, we hypothesise that the longer extracellular domain could lead to a pore site that is less accessible to formamidines and might confer resistance to them.

The BME/CTVM5 cell line, derived from the same R. microplus strain as BME/CTVM6, has potentially three different βAOR genes, only one of which was transcribed. The transcribed gene is completely homologous to the published βAOR gene JN974909 [25]. One of the other two gene variants appears to give the same fragment size as that of BME/CTVM6 but was not transcribed. Sequencing of this 220 bp amplicon yielded the same sequence as the alternative gene and transcript detected in BME/CTVM6. Multiple attempts to sequence the 245 bp amplicon that was detected only in BME/CTVM5 were unsuccessful, often producing very noisy signals. This suggests multiple primer binding sites within the 245 bp fragment. A repeat PCR performed on the gel-extracted and cleaned 245 bp PCR amplicon from the initial PCR produced the same three amplicons, confirming the existence of multiple binding sites for the primers, most likely due to an evolutionary process of gene duplication. It is notable that despite BME/CTVM5 cells having potentially three variants of the βAOR gene, present at both low and high passage levels, only one seems to have been transcribed. As the two cell lines were each derived from different pooled egg batches laid by over 50 female ticks, a possible explanation for this phenomenon is that BME/CTVM5 and BME/CTVM6 were derived from a population of ticks in which resistance was developing but in which the resistance-conferring allelic variants were not fixed at 100 % in the population. The establishment of the two cell lines would likely have been subject to normal stochastic effects on founder populations and subsequent genetic drift, leading to fixation of the initially more frequent allele over time after establishment of the culture. If this is the case, it highlights the potential value of the heterogeneous cell culture model for studies on the evolution of acaricide resistance, and draws attention to aspects of the system that must be taken into consideration in the design of experiments. We hypothesise that the βAOR gene was first duplicated in the base population. Gene duplication or multiple gene copy number has been linked to evolutionary adaptation in the parasitic trematode Fasciola hepatica [26], the mosquito Anopheles gambiae [27] and R. microplus [28, 29].

Pohl et al. [18] reported an association between increased expression of ABC transporter genes, ABCB10 in particular, and ivermectin resistance. BME/CTVM6 demonstrated approximately seven-fold higher ABCB10 expression when compared to BmVIII-SCC, a cell line derived from acaricide-susceptible ticks (Patricia Holman, personal communication). The importance of ABCB10 expression level in ivermectin resistance was demonstrated in a R. microplus [18] cell line while in an I. scapularis cell line, another ABC transporter gene, ABCB8, was upregulated [30]. BME/CTVM5 demonstrated a trend towards higher ABCB10 expression compared to BmVIII-SCC. In order to verify if these increased gene expression levels translate to protein expression levels, further work will be necessary.

It is interesting that the BME/CTVM6 cell line derived from an acaricide-resistant strain has both a different βAOR gene and a higher level of ABCB10 expression compared to the other R. microplus cell lines, whereas the BME/CTVM5 cell line, derived from the same strain, did not show these characteristics. It is also notable that both BME/CTVM5 and BME/CTVM6 were established over 30 years ago, prior to the recent heavy and widespread use of macrocyclic lactone products and the documented emergence of ivermectin resistance. It is possible that ABCB10 played a role in resistance to acaricides other than ivermectins, but it is unlikely that any mutation in βAOR would lead to resistance to macrocyclic lactones, as the macrocyclic lactones are not known to bind to βAOR. The fact that the same tick strain has returned quite different results in terms of the expression of a detoxifying gene and the genotype of a target gene clearly demonstrates the importance of initial stochastic effects at the time of the establishment of the cell culture.

Nine synonymous SNPs were found in R. appendiculatus, R. evertsi and R. sanguineus cell lines. Some of these SNPs appear to be unique to individual species. These SNPs provide a potential tool for identifying the tick species; however, they need to be verified with R. appendiculatus, R. evertsi and R. sanguineus ticks from different regions of the world and with known acaricide resistance status.

Failure to detect the βAOR gene in gDNA extracted from the Amblyomma, Hyalomma and Ixodes spp. cell lines could be due to insufficient sequence similarity between the tick genera, such that the primers designed for R. microplus did not bind to the gene in DNA or RNA extracted from these cell lines. A search for adrenergic octopamine receptor genes in Ixodes scapularis (https://www.vectorbase.org/organisms/ixodes-scapularis and http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/nuccore) found nothing, but there were numerous putative G-protein coupled octopamine receptors. Alternatively, the βAOR gene may be absent from these species or genera. It is more likely that the annotation for β-adrenergic-like octopamine receptor is incomplete.

The R. microplus cell lines have proven to be useful tools for investigating tick genetics in relation to acaricide resistance as evidenced from the findings of this study and other reports [18, 21, 22]. This study has also shown that cultures established from field populations have the potential to include and fix different variants of the same genes, such that some care is needed in interpretation of the results.

Conclusions

The present study has identified a new βAOR gene and demonstrated a higher ABCB10 expression level in the BME/CTVM6 cell line, indicating that tick cell lines provide a useful experimental tool for acaricide resistance studies and further elucidation of tick genetics.

Abbreviations

aa, amino acid; ABCB10, ATP-binding cassette B10; ANOVA, analysis of variance; bp, base-pair; cDNA, complementary deoxyribonucleic acid; gDNA, genomic deoxyribonucleic acid; NTC, non-template control; PCR, polymerase chain reaction; qPCR, real-time quantitative PCR; βAOR, β-adrenergic octopamine receptor; SNP, single nucleotide polymorphism

References

  1. 1.

    Sales-Junior PA, Guzman F, Vargas MI, Sossai S, Patarroyo VA, Gonzalez CZ, et al. Use of biodegradable PLGA microspheres as a slow release delivery system for the Boophilus microplus synthetic vaccine SBm7462. Vet Immunol Immunopathol. 2005;107:281–90.

  2. 2.

    Jonsson NN. The productivity effects of cattle tick (Boophilus microplus) infestation on cattle, with particular reference to Bos indicus cattle and their crosses. Vet Parasitol. 2006;137:1–10.

  3. 3.

    Rosario-Cruz R, Almazan C, Miller RJ, Dominguez-Garcia DI, Hernandez-Ortiz R, de la Fuente J. Genetic basis and impact of tick acaricide resistance. Front Biosci. 2009;14:2657–65.

  4. 4.

    Coetzer JAW, Tustin RC. Infectious Diseases of Livestocks, Vol 1. 2nd ed. Cape Town: Oxford University Press; 2004.

  5. 5.

    Dantas-Torres F. The brown dog tick, Rhipicephalus sanguineus (Latreille, 1806) (Acari: Ixodidae): from taxonomy to control. Vet Parasitol. 2008;152:173–85.

  6. 6.

    Barre N, Uilenberg G. Spread of parasites transported with their hosts: case study of two species of cattle tick. Rev Sci Tech. 2010;29:149–47.

  7. 7.

    Rizzoli A, Silaghi C, Obiegala A, Rudolf I, Hubalek Z, Foldvari G, et al. Ixodes ricinus and its transmitted pathogens in urban and peri-urban areas in Europe: New hazards and relevance for public health. Front Public Health. 2014;2:251.

  8. 8.

    Abbas RZ, Zaman MA, Colwell DD, Gilleard J, Iqbal Z. Acaricide resistance in cattle ticks and approaches to its management: the state of play. Vet Parasitol. 2014;203:6–20.

  9. 9.

    Guerrero FD, Lovis L, Martins JR. Acaricide resistance mechanisms in Rhipicephalus (Boophilus) microplus. Rev Bras Parasitol Vet. 2012;21:1–6.

  10. 10.

    Guerrero FD, de León AA P, Rodríguez-Vivas RI, Jonsson NN, Miller RJ, Andreotti R. Acaricide research and development, resistance, and resistance monitoring. In: Sonenshire DE, Roe RM, editors. Biology of Ticks. Volume 2. 2nd ed. New York: Oxford University Press; 2014.

  11. 11.

    Rinkevich FD, Du Y, Dong K. Diversity and convergence of sodium channel mutations involved in resistance to pyrethroids. Pestic Biochem Physiol. 2013;106:93–100.

  12. 12.

    Morgan JA, Corley SW, Jackson LA, Lew-Tabor AE, Moolhuijzen PM, Jonsson NN. Identification of a mutation in the para-sodium channel gene of the cattle tick Rhipicephalus (Boophilus) microplus associated with resistance to synthetic pyrethroid acaricides. Int J Parasitol. 2009;39:775–9.

  13. 13.

    Baxter GD, Barker SC. Isolation of a cDNA for an octopamine-like, G-protein coupled receptor from the cattle tick, Boophilus microplus. Insect Biochem Mol Biol. 1999;29:461–7.

  14. 14.

    Evans PD, Maqueira B. Insect octopamine receptors: a new classification scheme based on studies of cloned Drosophila G-protein coupled receptors. Invert Neurosci. 2005;5:111–8.

  15. 15.

    Corley SW, Jonsson NN, Piper EK, Cutulle C, Stear MJ, Seddon JM. Mutation in the RmbetaAOR gene is associated with amitraz resistance in the cattle tick Rhipicephalus microplus. Proc Natl Acad Sci U S A. 2013;110:16772–7.

  16. 16.

    Dermauw W, Van LT. The ABC gene family in arthropods: comparative genomics and role in insecticide transport and resistance. Insect Biochem Mol Biol. 2014;45:89–110.

  17. 17.

    Pohl PC, Klafke GM, Carvalho DD, Martins JR, Daffre S, da Silva Jr VI, et al. ABC transporter efflux pumps: a defense mechanism against ivermectin in Rhipicephalus (Boophilus) microplus. Int J Parasitol. 2011;41:1323–33.

  18. 18.

    Pohl PC, Carvalho DD, Daffre S, Vaz Ida Jr S, Masuda A. In vitro establishment of ivermectin-resistant Rhipicephalus microplus cell line and the contribution of ABC transporters on the resistance mechanism. Vet Parasitol. 2014;204:316–22.

  19. 19.

    Bell-Sakyi L, Attoui H. Endogenous tick viruses and modulation of tick-borne pathogen growth. Front Cell Infect Microbiol. 2013;3:25.

  20. 20.

    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:450–7.

  21. 21.

    Cossio-Bayugar R, Wagner GG, Holman PJ. In vitro generation of organophosphate resistant Boophilus microplus (Acari: Ixodidae) cell lines. J Med Entomol. 2002;39:278–84.

  22. 22.

    Cossio-Bayugar R, Barhoumi R, Burghardt RC, Wagner GG, Holman PJ. Basal cellular alterations of esterase, glutathione, glutathione S-transferase, intracellular calcium, and membrane potentials in coumaphos-resistant Boophilus microplus (Acari: Ixodidae) cell lines. Pestic Biochem Physiol. 2002;72:1–9.

  23. 23.

    Bell-Sakyi L, Palomar AM, Bradford EL, Shkap V. Propagation of the Israeli vaccine strain of Anaplasma centrale in tick cell lines. Vet Microbiol. 2015;179:270–6.

  24. 24.

    Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.

  25. 25.

    Corley SW, Piper EK, Jonsson NN. Generation of full-length cDNAs for eight putative GPCnR from the cattle tick, R microplus using a targeted degenerate PCR and sequencing strategy. PLoS One. 2012;7:e32480.

  26. 26.

    Cwiklinski K, Dalton JP, Dufresne PJ, La CJ, Williams DJ, Hodgkinson J, et al. The Fasciola hepatica genome: gene duplication and polymorphism reveals adaptation to the host environment and the capacity for rapid evolution. Genome Biol. 2015;16:71.

  27. 27.

    Assogba BS, Djogbenou LS, Milesi P, Berthomieu A, Perez J, Ayala D, et al. An ace-1 gene duplication resorbs the fitness cost associated with resistance in Anopheles gambiae, the main malaria mosquito. Sci Rep. 2015;5:14529.

  28. 28.

    Temeyer KB, Olafson PU, Pruett JH. Sequence polymorphism in acetylcholinesterase transcripts and genotyping survey of BmAchE1 in laboratory and Mexican strains of Rhipicephalus (Boophilus) microplus. J Med Entomol. 2012;49:555–62.

  29. 29.

    Temeyer KB, Tuckow AP, Brake DK, Li AY, Perez de Leon AA. Acetylcholinesterases of blood-feeding flies and ticks. Chem Biol Interact. 2013;203:319–22.

  30. 30.

    Mangia C, Vismarra A, Kramer L, Bell-Sakyi L, Porretta D, Otranto D, et al. Evaluation of the in vitro expression of ATP binding-cassette (ABC) proteins in an Ixodes ricinus cell line exposed to ivermectin. Parasit Vectors. 2016;9:215.

  31. 31.

    Spyropoulos IC, Liakopoulos TD, Bagos PG, Hamodrakas SJ. TMRPres2D: high quality visual representation of transmembrane protein models. Bioinformatics. 2004;20:3258–60.

  32. 32.

    Holman PJ. Partial characterization of a unique female diploid cell strain from the tick Boophilus microplus (Acari: Ixodidae). J Med Entomol. 1981;18:84–8.

  33. 33.

    Bell-Sakyi L. Ehrlichia ruminantium grows in cell lines from four ixodid tick genera. J Comp Pathol. 2004;130:285–93.

  34. 34.

    Alberdi MP, Nijhof AM, Jongejan F, Bell-Sakyi L. Tick cell culture isolation and growth of Rickettsia raoultii from Dutch Dermacentor reticulatus ticks. Ticks Tick Borne Dis. 2012;3:349–54.

  35. 35.

    Alberdi MP, Dalby MJ, Rodriguez-Andres J, Fazakerley JK, Kohl A, Bell-Sakyi L. Detection and identification of putative bacterial endosymbionts and endogenous viruses in tick cell lines. Ticks Tick Borne Dis. 2012;3:137–46.

  36. 36.

    Bekker CP, Bell-Sakyi L, Paxton EA, Martinez D, Bensaid A, Jongejan F. Transcriptional analysis of the major antigenic protein 1 multigene family of Cowdria ruminantium. Gene. 2002;285:193–201.

  37. 37.

    Varma MG, Pudney M, Leake CJ. The establishment of three cell lines from the tick Rhipicephalus appendiculatus (Acari: Ixodidae) and their infection with some arboviruses. J Med Entomol. 1975;11:698–706.

  38. 38.

    Yunker CE, Cory J, Meibos H. Tick tissue and cell culture: applications to research in medical and veterinary acarology and vector-borne disease. In: Griffiths DA, Bowman CE, editors. Acarology 6, vol. 2. Chichester: Ellis Horwood; 1984.

  39. 39.

    Yunker CE. Preparation and maintenancee of arthropod cell cultures: Acari, with emphasis on ticks. In: Yunker CE, editor. Arboviruses in arthropod cells in vitro, vol. 1. Boca Raton: CRC Press; 1987.

  40. 40.

    Bell-Sakyi L, Paxton EA, Munderloh UG, Sumption KJ. Morphology of Cowdria ruminantium grown in two tick cell lines. In: Kazimirova M, Labuda M, Nuttall PA, editors. Proceedings of the 3rd International Conference “Ticks and Tick-borne Pathogens: Into the 21st Century”. Bratislava: Institute of Zoology, Slovak Academy of Sciences; 2000. p. 131–7.

  41. 41.

    Bell-Sakyi L. Continuous cell lines from the tick Hyalomma anatolicum anatolicum. J Parasitol. 1991;77:1006–8.

  42. 42.

    da Silva Jr VI, Imamura S, Nakajima C, de Cardoso FC, Ferreira CA, Renard G. Molecular cloning and sequence analysis of cDNAs encoding for Boophilus microplus, Haemaphysalis longicornis and Rhipicephalus appendiculatus actins. Vet Parasitol. 2005;127:147–55.

Download references

Acknowledgements

The authors would like to thank Dr Patricia Holman of Texas A&M University for provision of the BmVIII-SCC cell line.

Funding

This project was funded through ES’s honours project and BBSRC research experience placement (REP) awarded to KC (through allocations of REPs direct to University of Glasgow, School of Veterinary Medicine 2015). LBS is supported by the United Kingdom Biotechnology and Biological Sciences Research Council’s Institute Strategic Programme Grant (BBS/E/I/00001741) to the Pirbright Institute.

Availability of data and material

All data generated or analysed during this study are included in this published article and its supplementary information files. Sequences are submitted in the GenBank database under accession numbers KU836738–KU836748.

Authors’ contributions

CKT contributed to experimental design, completed the sequencing and statistical analysis, interpreted the data and prepared the manuscript. ES and KC performed the majority of the experiments for this study. LBS provided the tick cell lines and contributed to the preparation of this manuscript. NJ conceived and designed the study, interpreted the data and contributed to the preparation of this manuscript. All authors read and approved the final version of the manuscript.

Authors’ information

CKT is research associate for principal investigator NJ; ES was an MSc Honours project student; KC is a BVMS student; and LBS is a collaborator and the curator of Tick Cell Biobank at The Pirbright Institute.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Correspondence to Nicholas N. Jonsson.

Additional files

Additional file 1: Figure S1.

βAOR gene in earlier passages of Rhipicephalus cell lines BME/CTVM5 and BME/CTVM6. Detection of βAOR gene in Rhipicephalus microplus cell lines BME/CTVM5 passage 8 (5p8) and BME/CTVM6 passages 32, 221 and 243 (6p32, 6p221 and 6p243, respectively). Amplicons of 183 bp, 220 bp and 245 bp were detected in the gDNA of BME/CTVM5 and amplicon of only 220 bp was detected in BME/CTVM6 passages. M = Marker. (TIF 2846 kb)

Additional file 2: Figure S2.

βAOR sequences of the 245 bp amplicon from BME/CTVM5 genomic DNA. There were two regions of heterogeneity from the sequencing traces of the alternative 245 bp amplicon detected in BME/CTVM5, indicating the existence of multiple primer binding sites for the sequencing reactions. (TIF 9962 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Koh-Tan, H.H.C., Strachan, E., Cooper, K. et al. Identification of a novel β-adrenergic octopamine receptor-like gene (βAOR-like) and increased ATP-binding cassette B10 (ABCB10) expression in a Rhipicephalus microplus cell line derived from acaricide-resistant ticks. Parasites Vectors 9, 425 (2016) doi:10.1186/s13071-016-1708-x

Download citation

Keywords

  • Acaricide resistance
  • Amitraz
  • Synthetic pyrethroid
  • Ivermectin
  • Dieldrin
  • Tick cell line

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate. Please note that comments may be removed without notice if they are flagged by another user or do not comply with our community guidelines.