- Open Access
Glutathione S-transferases play a role in the detoxification of flumethrin and chlorpyrifos in Haemaphysalis longicornis
Parasites & Vectorsvolume 11, Article number: 460 (2018)
Haemaphysalis longicornis is a tick of importance to health, as it serves as a vector of several pathogens, including Theileria orientalis, Babesia ovata, Rickettsia japonica and the severe fever with thrombocytopenia syndrome virus (SFTSV). Presently, the major method of control for this tick is the use of chemical acaricides. The glutathione S-transferase (GST) system is one mechanism through which the tick metabolizes these acaricides. Two GSTs from H. longicornis (HlGST and HlGST2) have been previously identified.
Enzyme kinetic studies were performed to determine the interaction of acaricides with recombinant H. longicornis GSTs. Recombinant HlGST activity was inhibited by flumethrin and cypermethrin, while recombinant HlGST2 activity was inhibited by chlorpyrifos and cypermethrin. Using real-time RT-PCR, the upregulation of the HlGST gene was observed upon exposure to sublethal doses of flumethrin, while the HlGST2 gene was upregulated when exposed to sublethal doses of chlorpyrifos. Sex and strain dependencies in the induction of GST gene expression by flumethrin were also observed. Knockdown of the HlGST gene resulted in the increased susceptibility of larvae and adult male ticks to sublethal doses of flumethrin and the susceptibility of larvae against sublethal doses of chlorpyrifos was increased upon knockdown of HlGST2.
HlGST could be vital for the metabolism of flumethrin in larvae and adult male ticks, while HlGST2 is important in the detoxification of chlorpyrifos in larval ticks.
The hard tick Haemaphysalis longicornis is a blood-sucking arthropod widely distributed in East Asia and Australia. They are known vectors of Theileria orientalis, Babesia ovata and Rickettsia japonica, as well as the severe fever with thrombocytopenia syndrome virus (SFTSV) [1, 2]. Until now, tick control has relied mainly upon the application of acaricides such as amitraz, synthetic pyrethroids and organophosphates [3, 4]. The continuous use of acaricides has some ill effects, such as the development of resistance against these acaricides . One factor that has contributed to the development of acaricide resistance is the improper application of acaricides, with a particular emphasis on sublethal doses .
Amitraz is a formamidine acaricide that has been around since the 1960s. It acts as an agonist of the octopaminergic receptors of arthropods, leading to the stimulation of monoamine oxidases and G proteins. This stimulation leads to the synthesis of cAMP and cGMP, and could alter the behavior of the arthropods . On the other hand, organophosphates are acetylcholinesterase inhibitors. Their action results in acetylcholine continuous stimulation, causing hyperactivity and, eventually, the death of the arthropods. Organophosphates possess a broad spectrum of activity against insects and acarians . Pyrethroids are esters capable of opening Na+ channels, resulting in the depolarization of nerve cell membranes. The effect on arthropods involves two phases. The first phase is the reversible “knockdown” effect, wherein arthropods cease all movements and act as if they are dead. This is caused by pyrethroids acting on the cerebral ganglia. Arthropods may still wake-up, and then go to the second phase, wherein the pyrethroids may affect peripheral nerves. This could result in brief, rapid and inconsistent movements of arthropods, perhaps leading eventually to death. The pyrethroids permethrin, cypermethrin and deltamethrin are both acaricides and insecticides, while flumethrin is mainly an acaricide .
Until now, the exact mechanisms of acaricide metabolism by the ticks are not fully understood. Three metabolic pathways are believed to play roles in this detoxification process: carboxylesterases, monooxygenases and glutathione S-transferases . Glutathione S-transferases (GSTs) are multifunctional enzymes that are responsible for the metabolism and detoxification of both xenobiotic and physiological substances. Metabolism involves the catalysis of thiol additions of the reduced glutathione to organic compounds through their electrophilic centers. The formation of more water-soluble conjugates would facilitate their elimination . The transport of molecules is facilitated by an ATP-requiring active transport system through the multidrug resistance-related protein (MRP) .
We have already identified two GST molecules of H. longicornis . GSTs act on a wide variety of substrates, and each GST isoenzyme may function very differently from the others, wherein not all GSTs are involved in the detoxification of acaricides [10, 11]. Therefore, it is important to determine the role of the GSTs of ticks in acaricide metabolism. Targeting specific GSTs that impede the ability of the arthropod to survive acaricide application could be included in the development of a tick control plan . In this study, we determined the possibility of the interaction of recombinant GSTs and several acaricides. We also observed the ability of sublethal doses of acaricides to induce GST gene and protein expression. Finally, we were able to establish the significance of GSTs in the metabolism of sublethal doses of acaricides through RNA interference (RNAi) experiments.
Ticks and experimental animals
The parthenogenetic Okayama strain and the bisexual Oita strain of H. longicornis were used in the experiments throughout this study. In Japan, no evidence of resistance against acaricides was reported in this tick species. Ticks were maintained by feeding on the ears of Japanese white rabbits (KBT Oriental, Saga, Japan) for several generations at the Laboratory of Infectious Diseases, Joint Faculty of Veterinary Medicine, Kagoshima University, Kagoshima, Japan . Experimental animals were kept at 25 °C and 40% relative humidity, with a constant supply of water and commercial feeds. The ticks, on the other hand were kept in glass tubes sealed with cotton plug and maintained at 15 °C and 80–85% relative humidity in an incubator. The ticks were maintained for 2 to 3 months after hatching or molting before use. The care and use of experimental animals in this study were approved by the Animal Care and Use Committee of Kagoshima University (Approval number VM15055).
Organophosphate acaricides (ethion, coumaphos, chlorpyrifos and diazinon), pyrethroids (cypermethrin and flumethrin), an avermectin (ivermectin), and a formamidine compound (amitraz) were evaluated for their interaction with the GSTs of H. longicornis. Ethion, chlorpyrifos, diazinon, cypermethrin and amitraz were purchased from Wako Pure Chemical Industries, Ltd. (Osaka, Japan). Coumaphos, ivermectin and flumethrin were purchased from Sigma-Aldrich (St. Louis, MO, USA).
Enzyme activity assay
The inhibition activity of recombinant GSTs was measured according to the methods of Habig and da Silva [3, 14], using 1-chloro-2,4-dinitrobenzene (CDNB) (Sigma-Aldrich) as a substrate. The recombinant GSTs used in this study were expressed as described previously . Two-hundred microliters of the reaction mixture consisting of varying concentrations of CDNB (0.125, 0.25, 0.5, 1 and 2 mM) dissolved in methanol, 5 mM glutathione, 0.1 mM of acaricide dissolved in methanol, and 120 μM recombinant GST in 100 mM Tris-HCl (pH 7.5) or without recombinant GST for the blank was tested in a 96-well plate. The methanol concentration was maintained at 5%. Equine liver GST and recombinant H. longicornis 2-cys-peroxiredoxin  were used as the positive and negative controls, respectively. Absorbance (A340nm) was measured each minute in an SH-9000 microplate reader (Corona Electric, Ibaraki, Japan) at 25 °C for 5 min. The extinction coefficient of 9.6 mM-1cm-1, corrected for the 96-well microplate light path, was used. Each assay was done in triplicate, and the results were expressed as the mean of three separate experiments. Kinetic constants Km and Vmax were calculated from the double-reciprocal plot of 1/v versus the 1/[S] or Lineweaver-Burk plot in which Vmax = 1/y-intercept of the regression line and Km = Vmax × slope of the regression line.
Determination of acaricide sublethal dose
The following acaricides were dissolved in methanol at 4 dilutions: flumethrin (0.4 μM, 4 μM, 40 μM and 400 μM), chlorpyrifos (0.01 mM, 0.1 mM, 1 mM and 10 mM), and amitraz (0.01 mM, 0.1 mM, 1 mM and 10 mM). These concentrations were based on previous studies [6, 16].
For exposure studies, the methods of Duscher  were used with some modifications. Briefly, 0.5 ml of each dilution was spotted onto a 10 × 5 cm piece of filter paper in scattered dots and dried under the fume hood for at least 2 h. Each group of 10 parthenogenetic females, nymphs and larvae, as well as bisexual male and female ticks, was placed in the acaricide-impregnated filter paper and exposed for 48 h. Mortality was checked after 48 h. For the parthenogenetic larvae and nymphs, further 10-fold dilutions were made to determine the sublethal dose. The maximum sublethal dose in this experiment is the highest dose that has either failed to cause any mortality or caused just a single mortality out of all the ticks tested .
GST gene and protein expression analysis of parthenogenetic female ticks exposed to flumethrin, chlorpyrifos and amitraz
Parthenogenetic female ticks were exposed to different sublethal concentrations of flumethrin (0, 0.4, 4 and 40 μM), chlorpyrifos (0, 0.01, 0.1 and 1 mM), and amitraz (0, 0.01, 0.1 and 1 mM). Total RNA was extracted from whole tick samples by homogenizing using an automill (Tokken, Chiba, Japan) and were added to TRI Reagent® (Sigma-Aldrich). RNA extraction was performed following the manufacturer’s protocol. Subsequently, single-stranded cDNA was prepared by reverse transcription using the ReverTra Ace® cDNA Synthesis Kit (Toyobo, Osaka, Japan), following the manufacturer’s protocol. Transcription analysis of HlGST and HlGST2 genes was performed through real-time RT-PCR using THUNDERBIRD™ SYBR® qPCR Mix (Toyobo) with an Applied Biosystems 7300 Real-Time PCR System using HlGST and HlGST2 real-time gene-specific primers (Table 1). Standard curves were made from fourfold serial dilutions of the cDNA of adult parthenogenetic ticks fed for 3 days. The PCR cycle profile used is as follows: 95 °C for 10 min, 40 cycles of a denaturation step at 95 °C for 15 s, and an annealing/extension step at 60 °C for 60 s. The data were analyzed with Applied Biosystems 7300 system SDS software. Actin, tubulin, P0 and L23 genes were evaluated for standardization at the first step of real-time RT-PCR. P0 genes were selected as an internal control for the ticks.
The protein was also extracted from the abovementioned parthenogenetic female ticks exposed to flumethrin, chlorpyrifos and amitraz. Whole tick samples were also homogenized using an automill (Tokken), and then suspended in phosphate-buffered saline (PBS) treated with Complete Mini Protease Inhibitor Cocktail Tablets (Roche, Mannheim, Germany). After sonication and recovery of the supernatant, tick proteins were separated with 12% SDS-polyacrylamide gel electrophoresis (SDS-PAGE) and transferred to a polyvinylidene difluoride (PVDF) membrane (Millipore, Billerica, MA, USA). The membrane was blocked overnight with 3% skim milk in PBS with 0.05% Tween 20, and it was then incubated with a primary antibody using mouse anti-GST sera (1:1000 dilution) for 1 h. β-tubulin was used as a control . After incubation with horseradish peroxidase-conjugated goat anti-mouse IgG (1:50,000 dilution; DakoCytomation, Glostrup, Denmark) for 1 h, the signal was detected using Clarity™ Western ECL Substrate (Bio-Rad Laboratories, Hercules, CA, USA) and analyzed using FluorChem FC2 software (Alpha Innotech, San Leandro, CA, USA).
RNA interference using double-stranded RNA (dsRNA) was performed to determine the effect of GST on survival upon acaricide exposure. The PCR primers used for the synthesis of dsRNA are listed in Table 1. The HlGST and HlGST2 fragments were amplified by PCR from plasmid clones using oligonucleotides, including HlGST T7 forward with HlGST RNAi reverse and HlGST T7 reverse with HlGST RNAi forward primers, and HlGST2 T7 forward with HlGST2 RNAi reverse and HlGST2 T7 reverse with HlGST2 RNAi forward primers, to attach the T7 promoter recognition sites on both forward and reverse ends. Enhanced green fluorescent protein (EGFP) gene fragments were amplified from pEGFP through PCR using oligonucleotides containing EGFP T7 forward and EGFP T7 reverse primers as well. PCR products were purified using a GENECLEAN II Kit (MP Biomedicals, Ilkrich, France). The T7 RiboMAX Express RNAi System (Promega, Madison, WI, USA) was used to synthesize dsRNA by in vitro transcription. The successful construction of dsRNA was confirmed by running 1 μl of the dsRNA products in 1.5% agarose gel in a Tris-acetate-EDTA (TAE) buffer. HlGST, HlGST2, or HlGST 1/2 dsRNA (0.5 μl) dissolved in high purity water was injected to the hemocoel of unfed adult ticks through the fourth coxa at 1μg/tick concentration . A total of 35 ticks per group were injected with dsRNA. The control group was injected with EGFP dsRNA. After injection, the ticks were held for 24 h in a 25 °C incubator to check for mortality resulting from injury during injection. The ticks were then kept in vials sealed with cotton plug, placed in a glass chamber, and maintained at 25 °C and 80–85% relative humidity in an incubator for another 72 h. For larvae and nymphs, the dsRNA immersion method described by Galay et al.  was performed. Briefly, a total of 35 larvae or nymphs were immersed in 40 μl of HlGST, HlGST2, or HlGST 1/2 dsRNA dissolved to a concentration of 0.5 μg/μl in high purity water for 12 h. After 12 h, the dsRNA solution was removed, and the ticks were checked for mortality resulting from immersion. Ticks were also kept in vials sealed with cotton plug, placed in a glass chamber, and maintained at 25 °C and 80–85% relative humidity in an incubator for another 72 h. Total RNA was extracted from 5 ticks of each developmental stage and sex, and their cDNA was synthesized. The cDNA was subjected to RT-PCR with a Hot Start Pol system (Jena Bioscience, Jena, Germany) using GST-specific primers, such as HlGST RT forward and HlGST RT reverse primers, and HlGST2 RT forward and HlGST2 RT reverse primers (Table 1), following the manufacturer’s instructions. The PCR cycle profile was as follows: 94 °C for 8 min, 30 cycles of a denaturation step at 94 °C for 30 s, an annealing step at 68 °C for 60 s, and an extension step at 72 °C for 60 s. The PCR products were run in 1.5% TAE agarose gel and stained with ethidium bromide in TAE buffer. Actin was used as a loading control. The absence of bands corresponding to HlGST and HlGST2 genes in their corresponding GST knockdown group demonstrates that silencing was successful (Additional file 1: Figure S1).
After confirmation of the knockdown, ticks were then exposed to sublethal doses of acaricides using the method described above. Mortality was checked after exposure. Ticks lying on their backs that could not turn over were considered dead.
Welch’s t-test was used to analyze data from the enzymatic inhibition, real-time RT-PCR of ticks, and tick survival studies. A significant difference is defined as P < 0.05. All experiments were done at least twice.
Interaction of recombinant GSTs with acaricides
Enzyme kinetic analysis in the presence or absence of acaricides was used to determine their ability to inhibit the activity of recombinant GST to catalyze the conjugation of CDNB to glutathione (GSH). The effect of acaricide interaction was determined by the change in the kinetic constants Vmax and Km, in accordance with the procedure of Mathews and van Holde . An inhibition that causes an increase in the Km without a change in the Vmax is a competitive type of inhibition. An inhibition in which the Km is not affected by decreased Vmax is a noncompetitive type of inhibition. An inhibition in which the Vmax and Km are decreased is an uncompetitive type of inhibition. In this experiment, flumethrin and cypermethrin showed uncompetitive inhibition of recombinant HlGST. On the other hand, chlorpyrifos and cypermethrin showed noncompetitive inhibition of recombinant HlGST2 (Table 2). Other acaricides, with the exception of coumaphos, did not show significant changes in the Vmax and Km. Coumaphos significantly decreased the Km of recombinant HlGST, indicating the apparent activation of the enzyme. These results demonstrated that the interactions of the recombinant GSTs with acaricides depend on the variety of acaricide.
Sublethal dose of flumethrin, chlorpyrifos and amitraz on different stages and strains of H. longicornis
Based on the results of enzyme kinetic analysis, flumethrin and chlorpyrifos were selected to analyze the importance of GSTs in their metabolism. Amitraz was also selected as a representative of the formamidine group. For adult parthenogenetic Okayama and bisexual Oita strains of ticks, the dilution of 40 μM proved to be the highest sublethal dose of flumethrin (Table 3), while 1 mM was the highest sublethal dose of chlorpyrifos (Table 4) and amitraz (Table 5). For nymphs, the dilution of 4 μM proved to be the highest sublethal dose of flumethrin (Table 3), while 100 μM was the highest sublethal dose of chlorpyrifos (Table 4) and amitraz (Table 5). For larvae, the dilution of 4 nM proved to be the highest sublethal dose of flumethrin (Table 3), while 100 nM was the highest sublethal dose of chlorpyrifos (Table 4) and amitraz (Table 5).
Effect of flumethrin, chlorpyrifos and amitraz on the gene and protein expression of GSTs of parthenogenetic ticks
The effect of flumethrin, chlorpyrifos and amitraz on the mRNA expression level of parthenogenetic female ticks were also tested (Fig. 1a). Exposure to 4 μM and 40 μM flumethrin resulted in the overexpression of HlGST. HlGST2 genes at 0.4 μM and 40 μM flumethrin exposure did not show any significant increase in mRNA expression, but they showed a significant increase at 4 μM flumethrin exposure (t(2) = 4.67, P = 0.043). Chlorpyrifos exposure did not result in any significant change in HlGST gene expression. On the other hand, the HlGST2 gene was overexpressed when ticks were exposed to 1 mM chlorpyrifos. Although amitraz exposure did not cause overexpression in either GST, HlGST showed a significant decrease in expression at a concentration of 0.01 mM (t(2) = 7.43, P = 0.018), while it increased significantly at a 0.1 mM concentration (t(2) = 6.92, P = 0.020). HlGST2 genes did not show any significant change when ticks were exposed to amitraz. These results demonstrated that flumethrin and chlorpyrifos trigger the overexpression of HlGST and HlGST2 genes, respectively. This might then indicate their possible role in the metabolism of these acaricides.
GST proteins of adult parthenogenetic H. longicornis ticks exposed to sublethal doses of flumethrin, chlorpyrifos, and amitraz were detected by Western blot analysis (Fig. 1b). The expression of HlGST protein was induced with exposure to 0.4 μM, 4 μM and 40 μM flumethrin. The protein expression of HlGST2 increased in a dose-dependent manner. Chlorpyrifos and amitraz at sublethal doses (0.01, 0.1 and 1 mM) did not cause any significant change in the expression of GST proteins. These results demonstrate that acaricides utilize GST proteins differently.
To be able to determine whether larval GSTs are induced in the same manner as the parthenogenetic female, parthenogenetic larvae exposed to sublethal doses of flumethrin (0, 0.4 and 4 nM) and chlorpyrifos (0, 10 and 100 nM) were checked for their gene and protein expression. Similar upregulation of HlGST genes and proteins was observed upon larval exposure to sublethal doses of flumethrin (Additional file 2: Figure S2). Unlike in adults, both the HlGST2 gene and protein were upregulated upon larval exposure to a sublethal dose of chlorpyrifos (Additional file 2: Figure S2). This indicates that the utilization of GSTs could vary between tick stages.
Effect of GST knockdown on different stages of parthenogenetic ticks upon exposure to flumethrin and chlorpyrifos
To further establish the importance of GST in acaricide detoxification, GST knockdown experiments were performed, and GST knockdown ticks were then exposed to different sublethal doses of flumethrin and chlorpyrifos. Nymphs and adults showed no significant increase in mortality in GST knockdowns as compared to EGFP knockdown groups (Figs. 2 and 3). The knockdown of HlGST in larvae resulted in the death of almost all larvae tested in 4 nM flumethrin (Fig. 2a). Additionally, a significant decrease in survival was also observed when both HlGST and HlGST2 were knocked down and larvae were exposed to 0.4 nM flumethrin as compared to HlGST knockdown alone (t(2) = 5.66, P = 0.030). On the other hand, the knockdown of HlGST2 and both HlGST and HlGST2 caused a significant increase in the mortality of larvae exposed to 10 μM (HlGST2: t(2) = 5.66, P = 0.030; HlGST and HlGST2: t(2) = 5.66, P = 0.030) and 100 μM chlorpyrifos (HlGST2: t(2) = 7.07, P = 0.019; HlGST and HlGST2: t(2) = 7.07, P = 0.019) (Fig. 3a). These results showed that HlGST is vital for the survival of larval ticks against sublethal doses of flumethrin, while HlGST2 is important in larval tick survival against chlorpyrifos.
Effect of GST knockdown on different sexes of ticks upon exposure to flumethrin and chlorpyrifos
In mammals, as well as insects, sexual differences in GST expression have been observed [20,21,22]; therefore, gene knockdown experiments with subsequent exposure to acaricides were also performed with male and female H. longicornis ticks to determine whether sex is a factor in tick survival against flumethrin and chlorpyrifos. Results have shown that the knockdown of HlGST and HlGST2 with subsequent exposure to sublethal doses of flumethrin leads to a significant increase in the mortality of male ticks (4 μM: t(2) = 4.92, P = 0.039; 40 μM: t(2) = 5.66, P = 0.030) (Fig. 2d). To check whether males and females have the same induction response, real-time RT-PCR and Western blot analysis were performed to check the expression levels of HlGST and HlGST2. Interestingly, although both GST genes are induced upon exposure to flumethrin, protein expressions were constant in male ticks (Fig. 4). On the other hand, female ticks showed no induction in either gene or protein expression. These results showed that HlGST is vital for the survival of male ticks against sublethal doses of flumethrin. These results also demonstrated that different strains of ticks have different induction responses to acaricides.
The predominant tick-control measure is the use of acaricides. Ticks make use of several mechanisms to metabolize these compounds. Therefore, several factors could be considered in the development of new tick-control strategies, such as the type of acaricide and its schedule of application. It is then important to understand the mechanism through which ticks metabolize these substances. Interference with these mechanisms would make the tick more prone to an acaricide and, eventually, lead to a more efficient tick control method [6, 23].
In this study, the role of two kinds of GSTs in the detoxification of several acaricides was investigated. We previously identified and expressed two mu-class GSTs from the hard tick H. longicornis , on which we used expressed recombinant GSTs to perform enzyme kinetic analysis in the presence of acaricides. The inhibition of GST activity by acaricides or insecticides has been observed previously [3, 24]. Differences in the type of inhibitions demonstrated that each GST isoenzyme has a specific profile of interaction with chemical compounds, even though the same chemical compounds would have different interactions with GST isoenzymes .
In cases of uncompetitive inhibition, inhibitors such as flumethrin and cypermethrin (Table 2) would only bind to an enzyme-substrate complex. This binding could also result in an irreversible interaction that may inactivate the enzyme . This inhibition also could be a result of the combination of GST and acaricide rather than conjugation of the acaricide with a reduced glutathione. This kind of binding of GSTs with pyrethroids was also observed in insects such as Tenebrio sp. and Aedes sp. Pyrethroids are believed to bind to the enzyme’s active site but did not yield a conjugated product. This suggests that the enzyme acts as a binding protein rather than a conjugating protein [25, 26]. This type of combination reaction was also observed in other detoxification enzymes, such as esterase. Tetranychus cinnabarinus esterase (TCE2) protein combined with abamectin rather than hydrolyzing it. It is believed that this binding decreased or delayed the noxious compound from reaching its target site; therefore, it is still considered an important mechanism in the metabolism of abamectin .
Noncompetitive inhibition occurs when an inhibitor, such as cypermethrin or chlorpyrifos (Table 2), binds to a non-substrate binding site. The presence of a non-substrate binding site was also observed in the pi class of GSTs . In this class of GSTs, the presence of low-affinity and high-affinity binding sites of the enzyme for bilirubin has been observed. Also, the inhibition of GST activity by bromosulfophthalein was attributed to the non-substrate binding site . Moreover, the noncompetitive inhibition of flumethrin was also observed in the recombinant Rhipicephalus (Boophilus) microplus GST .
In the present study, recombinant HlGST was activated by coumaphos through the lowering of the Km (Table 2). This kind of activation has also been observed in the recombinant R. microplus GST, wherein its activity was activated by coumaphos, the biological significance of which is still uncertain and could be the subject of a future study . Since the in vitro reaction showed a potential role of GSTs in interacting with acaricides, several types of acaricides that have interacted with GSTs were testedor their effect on the GST gene and protein expression levels of ticks exposed to acaricides.
Dose-dependent gene expression that led to overexpression was observed on the HlGST gene when adult ticks were exposed to sublethal doses of flumethrin (Fig. 1a). Although pyrethroids are not supposed to be substrates of GSTs, our enzyme kinetic analysis has shown the ability of recombinant HlGST to bind with flumethrin. This binding or sequestering mechanism of GSTs could give passive protection by either decreasing the level of free pyrethroids or facilitating the binding of other enzymes with it [26, 28]. This binding could be facilitated by alkyl or aryl hydrogen groups in the pyrethroids that could interact with GSTs. Also, the role of GST in detoxifying lipid peroxidation products caused by pyrethroids should be considered in determining the role of GST in the metabolism of an acaricide . The overexpression of GST genes of H. longicornis upon exposure to flumethrin was also observed by Hatta et al. . The overexpression of mu-class GST genes was also observed in mites after exposure to pyrethroids [28, 30]. Although the larva has a higher basal GST expression as compared with the adult , the ability of flumethrin to induce HlGST gene and protein expression remains the same. Interestingly, when the bisexual strain of H. longicornis was used, a different expression pattern was observed (Fig. 4a). No induction of GST genes was observed in the female bisexual strain as compared with the parthenogenetic H. longicornis. In the case of bisexual males, the induction of gene expression was observed for both HlGST and HlGST2 genes. This different expression of GSTs among arthropod strains was also observed in Anopheles gambiae, in which two strains (G3 and ZANDS) demonstrated different expressions in response to DDT. It is believed that GST expression is greatly influenced by environmental factors. Environmental factors, especially those that could confer resistance to a pesticide, include temperature, the type and dose of pesticide to which the arthropod has been exposed, and even the solvent quality used to dilute the pesticide . The continuous effects of these environmental factors have already resulted in changes at the genetic level [11, 32]. In the case of sex, cDNA analysis showed no difference between parthenogenetic, male and female ticks (Additional file 3: Figure S3). It should be noted that male and female ticks have demonstrated different gene expression patterns upon exposure to sublethal doses of flumethrin (Fig. 4a). Although the exact mechanism of the difference in male and female GST gene induction still remains unknown, differences in GST activity have been observed previously in mosquitoes . In the locust, different GST gene expressions were also observed between spermaries and ovaries . Male and female ticks could have different GST gene composition levels. In this way, when ticks are exposed to flumethrin, a GST with a higher activity against flumethrin would result in increased transcription. Therefore, if there is another unknown GST that has higher activity than the identified GSTs, the unknown GST could have proliferated instead of the identified ones. Therefore, if there is an unknown GST gene present in the female and absent in the male, it would result in different transcriptions in the known GST genes. On the other hand, when mosquitoes were exposed to pyrethroids, the overexpression of multiple cytochrome P450 genes was observed aside from the upregulated GST genes [33, 34].
The HlGST2 gene showed overexpression at the highest sublethal dose of chlorpyrifos (Fig. 1a). However, it is well known that organophosphates are metabolized by cytochrome P450 monooxygenases and hydrolases. Specifically, chlorpyrifos is activated to chlorpyrifos oxon by cytochrome P450 enzymes before being deethylated or dearylated. The metabolism of chlorpyrifos could yield metabolites that could be subjected to GSH conjugation by GSTs . Therefore, GSTs could have an indirect role in organophosphate detoxification. In R. (B.) microplus, the overexpression of GST was observed when ticks were exposed to coumaphos . On the other hand, chlorpyrifos has shown the ability to overexpress the GST gene of the migratory locust, cotton leaf worm and rice plant hopper [36,37,38]. Interestingly, HlGST2 genes appear to be maintained at lower concentrations of chlorpyrifos (Fig. 1a). The noninduction of HlGST2 at low doses of chlorpyrifos could be because of the organism’s ability to specifically produce the appropriate GST genes depending on its needs. Different responses of GSTs to chlorpyrifos were also observed in migratory locusts . On the other hand, amitraz did not show any effect on the GST gene expression level. Since amitraz is considered to be an agonist in the octopaminergic system of arthropods , this mimicry could have resulted in its non-recognition as a foreign or xenobiotic compound and, therefore, not a metabolite of the GSTs.
In accordance with the GST gene expression level, the protein expression level of HlGST2 also increased depending on the concentration of flumethrin (Fig. 1b). This also held true in larvae of ticks exposed to sublethal doses of chlorpyrifos. When ticks were exposed to a sublethal dose of flumethrin, HlGST protein initially increased but appeared to maintain its expression when the dosage of flumethrin was increased, even though gene expression increased dramatically (Fig. 1b). On the other hand, no difference was observed in the protein expression level regardless of the concentration when adult ticks were exposed to sublethal doses of chlorpyrifos (Fig. 1b). It has always been believed that GST proteins are transcriptionally regulated . Based on our results and previous studies, GST not only functions as a conjugate with GSH but also as a binding protein, wherein they “sacrifice” themselves by binding covalently with reactive compounds . The fate of the GST conjugates is yet to be demonstrated. They could possibly be released from the cytoplasm, as some studies have shown the ability of GST to be secreted or to move across the plasma membrane through the facilitation of MRP [22, 41]. This could be why non-substrates that bind with GST, such as flumethrin with HlGST or chlorpyrifos with HlGST2, did not have a drastic increase in protein expression, as GST proteins are being released from the cell after being bound to flumethrin or chlorpyrifos.
To further understand the role of GSTs in flumethrin and chlorpyrifos metabolism, RNAi was performed on GST genes, and ticks at different stages were exposed to sublethal doses of the acaricides. No significant differences in mortality were observed in the knockdown groups of adult female and nymph ticks exposed to sublethal doses of flumethrin and chlorpyrifos (Figs. 2 and 3). Other GST isoenzymes could have possibly compensated for the silenced GST . Several tick species have demonstrated the presence of multiple isoenzymes of GSTs. Thirty-five genes of GSTs, of which 14 belong to the mu-class GST were shown in an in silico analysis of the Ixodes scapularis gene database . Multiple GSTs have also been found in Dermacentor variabilis and Rhipicephalus (Boophilus) annulatus [3, 44]. Interestingly, the HlGST knockdown in male ticks resulted in increased susceptibility to flumethrin (Fig. 2d). This could mean that HlGST could be the main GST in the male tick’s detoxification mechanism against flumethrin. However, further testing and studies need to be conducted on the different detoxification mechanisms between male and female ticks. Notably, the larval stage of ticks also showed increased susceptibility to the effects of sublethal doses of flumethrin and chlorpyrifos. The importance of the GST system at the early stage of development was also observed in insects such as Tenebrio and Anopheles, as well as the red mite Panonychus, wherein increased expression and/or activity was observed at the younger stage of development as compared with the adult [20, 28, 45, 46]. Our previous results have also shown relatively higher GST protein expression in the unfed larval stages of ticks as compared to the unfed nymph and adult H. longicornis . Since the younger stages do not have a well-developed system of protection, such as integuments, it is possible that the GST system is vital for protecting larvae by detoxification.
Flumethrin may not be readily conjugated to glutathione by GSTs; however, GST could bind to flumethrin to decrease its ability to reach the tick nerve cell membrane. This, in turn, could result in a reduced toxic effect of flumethrin. The knockdown of HlGST in larvae and adult male ticks also could have resulted in increased intracellular flumethrin, eventually leading to cellular toxemia. The knockdown of GST genes in Rhipicephalus sanguineus also leads to increased mortality upon exposure to a sublethal dose of pyrethroids . In mites, exposure to pyrethroids after GST inhibition also resulted in increased susceptibility to acaricides . Higher mortality was also observed when both GSTs were knocked down, as compared to the knockdown of HlGST alone, when larvae were exposed to a sublethal dose of flumethrin. This could possibly show the ability of HlGST2 to compensate, to a certain degree, for the loss of HlGST. On the other hand, the metabolism of chlorpyrifos leads to the production of toxic metabolites . These metabolites could have increased when HlGST2 is knocked down in larvae. The abundance of these metabolites could have resulted in intracellular toxicity and, eventually, the death of ticks. The same increase in susceptibility to chlorpyrifos upon GST knockdown has also been observed in migratory locusts .
In conclusion, GSTs have been known to be involved in the detoxification of xenobiotic compounds . HlGST could play an important role in the detoxification of pyrethroids such as flumethrin, wherein the inhibition of enzymatic activity was observed. HlGST gene expression was also induced by sublethal doses of flumethrin, and the knockdown of HlGST resulted in the increased susceptibility of larvae and male ticks to flumethrin. On the other hand, HlGST2 plays an important function in the metabolism of organophosphates, such as chlorpyrifos. Chlorpyrifos inhibited the GST activity of recombinant HlGST2. HlGST2 gene induction upon exposure to sublethal doses of chlorpyrifos was also observed. More importantly, larval susceptibility to chlorpyrifos significantly increased upon HlGST2 knockdown. The data also showed that GSTs have a more important role in larval survival as compared to ticks at other stages upon exposure to sublethal doses of acaricides. We have also demonstrated that environmental factors as well as GST pool composition could play a role in the ability of acaricides to induce GST gene expression as observed in male and female ticks. These results have shown the importance of specific GSTs in the acaricide detoxification mechanism and could be considered by tick scientists in the development of new tick control strategies.
Enhanced green fluorescent protein
GST of H. longicornis
Sodium dodecyl sulphate polyacrylamide gel electrophoresis
Severe fever with thrombocytopenia virus
T. cinnabarinus esterase
Sivakumar T, Igarashi I, Yokoyama N. Babesia ovata: Taxonomy, phylogeny and epidemiology. Vet Parasitol. 2016;229:99–106.
Sonenshine DE, Roe RM. Biology of Ticks, vol. 2. 2nd ed. New York: Oxford University Press; 2014.
da Silva VI, Torino Lermen T, Michelon A, Sanchez Ferreira CA, de Joaquim Freitas DR, Termignoni C, et al. Effect of acaricides on the activity of a Boophilus microplus glutathione S-transferase. Vet Parasitol. 2004;119:237–45.
Ghosh S, Kumar R, Nagar G, Kumar S, Sharma AK, Srivastava A, et al. Survey of acaricides resistance status of Rhipicephalus (Boophilus) microplus collected from selected places of Bihar, an eastern state of India. Ticks Tick Borne Dis. 2015;6:668–75.
Lu H, Ren Q, Li Y, Liu J, Niu Q, Yin H, et al. The efficacies of 5 insecticides against hard ticks Hyalomma asiaticum, Haemaphysalis longicornis and Rhipicephalus sanguineus. Exp Parasitol. 2015;157:44–7.
Duscher GG, Galindo RC, Tichy A, Hummel K, Kocan KM, de la Fuente J. Glutathione S-transferase affects permethrin detoxification in the brown dog tick, Rhipicephalus sanguineus. Ticks Tick Borne Dis. 2014;5:225–33.
Beugnet F, Franc M. Insecticide and acaricide molecules and/or combinations to prevent pet infestation by ectoparasites. Trends Parasitol. 2012;28:267–79.
Wilce MCJ, Parker MW. Structure and function of glutathione S-transferases. Biochim Biophys Acta. 1994;1205:1–18.
Hernandez EP, Kusakisako K, Talactac MR, Galay RL, Hatta T, Matsuo T, et al. Characterization and expression analysis of a newly identified glutathione S-transferase of the hard tick Haemaphysalis longicornis during blood-feeding. Parasit Vectors. 2018;11:91.
Qin G, Jia M, Liu T, Zhang X, Guo Y, Zhu KY, et al. Characterization and functional analysis of four glutathione S-transferases from the migratory locust, Locusta migratoria. PLoS One. 2013;8:e58410.
Prapanthadara L, Promtet N, Koottathep S, Somboon P, Ketterman AJ. Isoenzymes of glutathione S-transferase from the mosquito Anopheles dirus species B: the purification, partial characterization and interaction with various insecticides. Insect Biochem Mol Biol. 2000;30:395–403.
Bartley K, Wright HW, Bull RS, Huntley JF, Nisbet AJ. Characterisation of Dermanyssus gallinae glutathione S-transferases and their potential as acaricide detoxification proteins. Parasit Vectors. 2015;8:350.
Fujisaki K. Development of acquired resistance precipitating antibody in rabbits experimentally infested with females of Haemaphysalis longicornis (Ixodoidea: Ixodidae). Natl Inst Anim Health Q (Tokyo). 1978;18:27–38.
Habig WH, Pabst MJ, Jakoby WB. Glutathione S-transferases. The first enzymatic step in mercapturic acid formation. J Biol Chem. 1974;249:7130–9.
Kusakisako K, Masatani T, Miyata T, Galay RL, Maeda H, Talactac MR, et al. Functional analysis of recombinant 2-Cys peroxiredoxin from the hard tick Haemaphysalis longicornis. Insect Mol Biol. 2016;25:16–23.
Hatta T, Matsubayashi M, Tsuji N. Microarray analysis of acaricide-inducible gene expression in the hard tick, Haemaphysalis longicornis. Bull Natl Inst Anim Health. 2014;120:69–72.
Umemiya-Shirafuji R, Tanaka T, Boldbaatar D, Tanaka T, Fujisaki K. Akt is an essential player in regulating cell/organ growth at the adult stage in the hard tick Haemaphysalis longicornis. Insect Biochem Mol Biol. 2012;42:164–73.
Galay RL, Hernandez EP, Talactac MR, Maeda H, Kusakisako K, Umemiya-Shirafuji R, et al. Induction of gene silencing in Haemaphysalis longicornis ticks through immersion in double-stranded RNA. Ticks Tick Borne Dis. 2016;7:813–6.
Matthews CK, van Holde KE. Biochemistry. Redwood City, California: The Benjamin/Cummings Publishing Company; 1990.
Tripathy A, Kar SK. Feeding stage, species, body part and sex-specific activity of glutathione S-transferase in mosquito. Trop Biomed. 2015;32:65–75.
Qin G, Jia M, Liu T, Xuan T, Yan Zhu K, Guo Y, et al. Identification and characterisation of ten glutathione S-transferase genes from oriental migratory locust, Locusta migratoria manilensis (Meyen). Pest Manag Sci. 2011;67:697–704.
Sherratt PJ, Hayes JD. Glutathione S-transferases. In: Chichester IC, editor. Enzyme Systems that Metabolise Drugs and Other Xenobiotics. West Sussex: John Wiley & Sons, Ltd; 2001.
Avinash B, Venu R, Alpha Raj M, Srinivasa Rao K, Srilatha C, Prasad TNVKV. In vitro evaluation of acaricidal activity of novel green silver nanoparticles against deltamethrin resistance Rhipicephalus (Boophilus) microplus. Vet Parasitol. 2017;237:130–6.
Hemingway J. The molecular basis of two contrasting metabolic mechanisms of insecticide resistance. Insect Biochem Mol Biol. 2000;30:1009–15.
Grant DF, Matsumura F. Glutathione S-transferase 1 and 2 in susceptible and insecticide resistant Aedes aegypti. Pestic Biochem Physiol. 1989;33:132–43.
Kostaropoulos I, Papadopoulos AI, Metaxakis A, Boukouvala E, Papadopoulou-Mourkidou E. Glutathione S-transferase in the defence against pyrethroids in insects. Insect Biochem Mol Biol. 2001;31:313–9.
Shi L, Wei P, Wang X, Shen G, Zhang J, Xiao W, et al. Functional analysis of esterase TCE2 gene from Tetranychus cinnabarinus (Boisduval) involved in acaricide resistance. Sci Rep. 2016;6:18646.
Liao C-Y, Zhang K, Niu J-Z, Ding T-B, Zhong R, Xia W-K, et al. Identification and characterization of seven glutathione S-transferase genes from citrus red mite, Panonychus citri (McGregor). Int J Mol Sci. 2013;14:24255–70.
Enayati AA, Ranson H, Hemingway J. Insect glutathione transferases and insecticide resistance. Insect Mol Biol. 2005;14:3–8.
Mounsey KE, Pasay CJ, Arlian LG, Morgan MS, Holt DC, Currie BJ, et al. Increased transcription of glutathione S-transferases in acaricide exposed scabies mites. Parasit Vectors. 2010;3:43.
Hadaway AB. Some factors affecting the distribution and rate of action of insecticides. Bull World Health Organ. 1971;44:221–4.
Prapanthadara LA, Hemingway J, Ketterman AJ. Partial purification and characterization of glutathione S-transferases involved in DDT resistance from the mosquito Anopheles gambiae. Pestic Biochem Physiol. 1993;47:119–33.
Gong Y, Li T, Zhang L, Gao X, Liu N. Permethrin induction of multiple cytochrome P450 genes in insecticide resistant mosquitoes, Culex quinquefasciatus. Int J Biol Sci. 2013;9:863–71.
Li C, Guo X, Zhang Y, Dong Y, Xing D, Yan T, et al. Identification of genes involved in pyrethroid-, propoxur-, and dichlorvos- insecticides resistance in the mosquitoes, Culex pipiens complex (Diptera: Culicidae). Acta Trop. 2016;157:84–95.
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:597–606.
Qin G, Liu T, Guo Y, Zhang X, Ma E, Zhang J. Effects of chlorpyrifos on glutathione S-transferase in migratory locust, Locusta migratoria. Pestic Biochem Physiol. 2014;109:1–5.
Zhou W-W, Liang Q-M, Xu Y, Gurr GM, Bao Y-Y, Zhou X-P, et al. Genomic insights into the glutathione S-transferase gene family of two rice planthoppers, Nilaparvata lugens (Stål) and Sogatella furcifera (Horváth) (Hemiptera: Delphacidae). PLoS One. 2013;8:e56604.
Zhang N, Liu J, Chen S-N, Huang L-H, Feng Q-L, Zheng S-C. Expression profiles of glutathione S-transferase superfamily in Spodoptera litura tolerated to sublethal doses of chlorpyrifos. Insect Sci. 2016;23:675–87.
Vandenberghe Y, Tee L, Rogiers V, Yeoh G. Transcriptional- and post-transcriptional-dependent regulation of glutathione S-transferase expression in rat hepatocytes as a function of culture conditions. FEBS Lett. 1992;313:155–9.
Jakoby WB, Keen JH. A triple-threat in detoxification: the glutathione S-transferases. Trends Biochem Sci. 1977;2:229–31.
Mukherjee SB, Aravinda S, Gopalakrishnan B, Nagpal S, Salunke DM, Shaha C. Secretion of glutathione S-transferase isoforms in the seminiferous tubular fluid, tissue distribution and sex steroid binding by rat GSTM1. Biochem J. 1999;340:309–20.
Bhattacharjee P, Paul S, Banerjee M, Patra D, Banerjee P, Ghoshal N, et al. Functional compensation of glutathione S-transferase M1 (GSTM1) null by another GST superfamily member, GSTM2. Sci Rep. 2013;3:2704.
Reddy BPN, Prasad GBKS, Raghavendra K. In silico analysis of glutathione S-transferase supergene family revealed hitherto unreported insect specific δ- and ε-GSTs and mammalian specific μ-GSTs in Ixodes scapularis (Acari: Ixodidae). Comput Biol Chem. 2011;35:114–20.
Dreher-Lesnick SM, Mulenga A, Simser JA, Azad AF. Differential expression of two glutathione S-transferases identified from the American dog tick, Dermacentor variabilis. Insect Mol Biol. 2006;15:445–53.
Kostaropoulos I, Mantzari A, Papodopoulos AI. Alterations of some glutathione characteristics during the development of Tenebrio molitor (Insecta: Coleoptera). Insect Biochem Mol Biol. 1996;26:962–9.
Strode C, Steen K, Ortelli F, Ranson H. Differential expression of the detoxification genes in the different life stages of the malaria vector Anopheles gambiae. Insect Mol Biol. 2006;15:523–30.
This study was supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI Grant Numbers 15H05264, 16H05028 and 17K19328; the Takeda Science Foundation; and the Japanese Government’s Ministry of Education, Culture, Sports, Science and Technology Scholarship (Monbukagakusho: MEXT) for doctoral fellowship. The funding bodies had no role in the design of the study; in the collection, analysis and interpretation of data; or in writing the manuscript.
Availability of data and materials
The data supporting the conclusions in this study are included within the article and its additional files.
The care and use of experimental animals in this study were approved by the Animal Care and Use Committee of Kagoshima University (approval number VM15055).
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. RT-PCR of knockdown ticks. Total RNA was extracted from whole GST and EGFP knockdown ticks 4 days post-injection/immersion to dsRNA. cDNA was synthesized and subjected to RT-PCR. PCR products were run on 1.5% TAE agarose gel and stained with ethidium bromide. Actin was used as a loading control. (PDF 1275 kb)
Figure S2. Gene (a) and protein (b) expressions of larval ticks upon exposure to sublethal doses of flumethrin and chlorpyrifos. Ticks were exposed to sublethal doses of flumethrin and chlorpyrifos for 48 h. a Total mRNA was extracted from ticks, and cDNA was then transcribed for real-time RT-PCR. P0 primers were used as a control. The error bar represents the mean ± standard deviation. *P < 0.05: significantly different by Welch’s t-test as compared to no treatment. The dotted line indicates overexpression. Overexpression is determined if there is at least a twofold increase in expression level, as shown by Bhattacharjee et al. . b Proteins were prepared from ticks exposed to sublethal doses of acaricides. Antiserum against tubulin was used as a control for Western blot analysis. Western blot analysis results are shown as representative data of three separate experiments showing the same trend. (PDF 817 kb)
Figure S3. Multiple sequence alignments of the nucleotide sequence of HlGST (a) and HlGST2 (b) of parthenogenetic, male and female ticks. Sequencing was performed using an automated sequencer (ABI PRISM 3100, Genetic Analyzer; Applied Biosystems, Foster City, CA, USA). Multiple sequence alignments of GST genes were done using MAFFT version 7 program (www.mafft.cbrc.jp). Asterisks indicate identity. Start and stop codons are in red. (ZIP 69 kb)