Expression analysis of glutathione S-transferases and ferritins during the embryogenesis of the tick Haemaphysalis longicornis

The diverse functions of Mu-class Glutathione S-transferase HrGSTm1 during the development of Hyalomma rufipes with a focus on the detoxification metabolism of cyhalothrin

BACKGROUND

Glutathione S-transferases (GSTs) are a superfamily of multifunctional enzymes in living organisms with metabolic and detoxification functions, which can detoxify exogenous and endogenous compounds and thereby reduce the damage caused by toxic substances to the body. Ticks are obligate blood-sucking ectoparasites that can transmit various pathogens, and the characterization of tick-derived GSTs may help improve current understanding of the molecular mechanism of tick resistance to insecticides. In this study, a novel GST gene, named HrGSTm1, was identified from Hyalomma rufipes.

METHODS

Sequence analysis was performed by using bioinformatics techniques. A prokaryotic expression system was used to obtain the recombinant expression protein rHrGSTm1. Detection of spatiotemporal expression patterns of target genes and their response to the toxicity of cyhalothrin on female H. rufipes was performed by using a quantitative PCR platform. The optimal enzymological parameters of rHrGSTm1 using glutathione as substrate were calculated. The antioxidant capacity of the recombinant protein was evaluated by DPPH• (1,1-Diphenyl-2-picrylhydrazyl radical 2,2-Diphenyl-1-(2,4,6-trinitrophenyl) hydrazyl). Knockdown of the HrGSTm1 genes through RNA interference was used to analyze their effects on the physiological parameters of ticks. The changes in HrGSTm1 messenger RNA expression patterns under cypermethrin stress were analyzed.

RESULTS

The complementary DNA sequence of HrGSTm1 contained a 672-bp open reading frame, which potentially encoded 223 amino acids. The predicted molecular weight was 25.62 kDa, and the isoelectric point 8.22. HrGSTm1 is a Mu-class GST, belonging to the cytoplasmic GSTs with no signal peptide observed. The Vmax and Km of rHrGSTm1 were 3.367 ± 0.81 uM and 2.208 ± 0.76 uM, respectively, and its activities were dependent on different temperatures and pH conditions; the scavenging rate of rHrGSTm1 to DPPH• reached 76.4% at 1.25 mg/ml. Variable expressions of HrGSTm1 were observed under various treatment periods and in different tissues, with the highest appearing in eggs (analysis of variance [ANOVA], F(2, 9) = 279.9, P < 0.0001) and Malpighian tubules (ANOVA, F(3, 12) = 290.5, P < 0.0001). After knockdown of HrGSTm1, compared with the control group, the mortality in the treatment group was increased by 16.7%, the average oviposition rate decreased by 33.9%, the average engorged body weight decreased by 287.38 mg and egg weight decreased by 127.46 mg, although only the engorged body weight was significantly different (t-test, t(44) = 2.886, P = 0.006). After exposure to three sublethal concentrations (LC05, LC10, LC50) of cyhalothrin, the expression level of HrGSTm1 in the midgut, ovary and salivary gland was upregulated, whereas in Malpighian tubules, it showed a trend of upregulation at first and then downregulation, implying different functions during the detoxification in different tissues.

CONCLUSIONS

In this study, a novel GST of the Mu-class was successfully isolated from H. rufipes and systematically subjected to bioinformatic analysis and recombination identification. The variation trend of HrGSTm1 expression level in different tissues suggests that the gene has different detoxification functions in different tissues. The potential function of this gene was analyzed to provide basic research for further investigation of its detoxification mechanism.

Elucidating structure and dynamics of glutathione S-transferase from Rhipicephalus (Boophilus) microplus

Rhipicephalus (Boophilus) microplus is tick parasite that affects the cattle industry worldwide. In R. (B.) microplus, acaricide resistance develops rapidly against many commercial acaricides. One of main resistance strategies is to enhance the metabolic detoxification mediated by R. (B.) microplus glutathione-S-transferase (RmGST). RmGST detoxifies acaricides by catalyzing the conjugation of glutathione to acaricides. Although structural and dynamic details of RmGST are expected to elucidate the biologic activity of this molecule, these data have not been available to date. Thus, Molecular Dynamics simulations were employed to study ligand-free RmGST at an atomic level. Like other m-class GSTs, the flexible m loop (m1) of RmGST was observed. M1 seems to shield the active sites from the bulk. A RmGST dimer is stabilized by the lock-and-key motif (F57 as "key") and hydrogen bonds of R82-E91 and R82-D98 at the dimer interface. Without substrates, conserved catalytic Y116 and N209 can interact with V112, G210 (for Y116) and F215 (for N209). Overall, most residues involving in RmGST function and stability are similar to other m-class GSTs. This implies similar structural stability and catalytic activity of RmGST to other GSTs. An insight obtained here will be useful for management of acaricide resistance and tick control.Communicated by Ramaswamy H. Sarma.

The diurnal salivary glands transcriptome of Dermacentor nuttalli from the first four days of blood feeding

The ixodid tick Dermacentor nuttalli is distributed from southern Siberia to North China and is a vector of many pathogens. This species can have severe impacts on animal husbandry and human health. To date, the control of D. nuttalli is limited to the use of acaricides such as organophosphorus, synthetic pyrethroids and amidine pesticides. There are no environmentally friendly or reliable prevention and control measures, and little is known regarding key antigens involved in blood feeding. Salivary glands are major tissues involved in the blood feeding and pathogen transmission of ticks. Therefore, this study focused on salivary glands tissue to identify the dominant antigens of D. nuttalli involved in tick feeding. For this, high-throughput RNA sequencing (RNA-seq) was used for analysis. The transcriptome of female D. nuttalli ticks was assembled and characterized, and differentially expressed genes (DEGs) were identified in the salivary glands of ticks that had not fed (0 h) and of ticks after 24, 48, 72 and 96 h of feeding. There were 22,802,784, 22,275,013, 26,629,453, 24,982,389, and 22,596,230 high-quality clean reads obtained from salivary glands tissues at the five different blood feeding time points. The total number of annotated unigenes was 100,347. The differences in gene expression between different time points were compared, and functional enrichment was performed. Quantitative reverse transcription PCR (RT‒qPCR) was used to validate the RNA-seq results, the results of which showed that the differences in expressed transcripts presented similar trends. Among the identified DEGs, the most numerous were those with catalytic and binding activities and those involved in diverse metabolic pathways and cellular processes. The expression patterns of homologous and family-member proteins throughout the blood feeding period exhibited significant differences, strongly suggesting that the transcriptome composition is highly dynamic and likely subjected to important variation throughout the life cycle. Studies of gene sequences in D. nuttalli will greatly increase the information on tick protective antigens, which could potentially function as effective vaccine candidates or drug targets for the development of environmentally friendly acaricides.

RNA activation in ticks

RNA activation (RNAa) is a burgeoning area of research in which double-stranded RNAs (dsRNAs) or small activating RNAs mediate the upregulation of specific genes by targeting the promoter sequence and/or AU-rich elements in the 3'- untranslated region (3'-UTR) of mRNA molecules. So far, studies on the phenomenon have been limited to mammals, plants, bacteria, Caenorhabditis elegans, and recently, Aedes aegypti. However, it is yet to be applied in other arthropods, including ticks, despite the ubiquitous presence of argonaute 2 protein, which is an indispensable requirement for the formation of RNA-induced transcriptional activation complex to enable a dsRNA-mediated gene activation. In this study, we demonstrated for the first time the possible presence of RNAa phenomenon in the tick vector, Haemaphysalis longicornis (Asian longhorned tick). We targeted the 3'-UTR of a novel endochitinase-like gene (HlemCHT) identified previously in H. longicornis eggs for dsRNA-mediated gene activation. Our results showed an increased gene expression in eggs of H. longicornis endochitinase-dsRNA-injected (dsHlemCHT) ticks on day-13 post-oviposition. Furthermore, we observed that eggs of dsHlemCHT ticks exhibited relatively early egg development and hatching, suggesting a dsRNA-mediated activation of the HlemCHT gene in the eggs. This is the first attempt to provide evidence of RNAa in ticks. Although further studies are required to elucidate the detailed mechanism by which RNAa occurs in ticks, the outcome of this study provides new opportunities for the use of RNAa as a gene overexpression tool in future studies on tick biology, to reduce the global burden of ticks and tick-borne diseases.

Proteomic analysis of extracellular vesicles from tick hemolymph and uptake of extracellular vesicles by salivary glands and ovary cells

BACKGROUND

Extracellular vesicles (EVs) are a heterogeneous group of cell-derived membranous structures that are important mediators of intercellular communication. Arthropods transport nutrients, signaling molecules, waste and immune factors to all areas of the body via the hemolymph. Little is known about tick hemolymph EVs.

METHODS

Hemolymph was collected from partially fed Rhipicephalus haemaphysaloides and Hyalomma asiaticum ticks by making an incision with a sterile scalpel in the middle (between the femur and metatarsus) of the first pair of legs, which is known as leg amputation. EVs were isolated from hemolymph by differential centrifugation and characterized by transmission electron microscopy (TEM) and nanoparticle tracking analysis (NTA). Proteins extracted from the hemolymph EVs were analyzed by 4D label-free proteomics. The EVs were also examined by western blot and immuno-electron microscopy analysis. Intracellular incorporation of PHK26-labeled EVs was tested by adding labeled EVs to tick salivary glands and ovaries, followed by fluorescence microscopy.

RESULTS

In this study, 149 and 273 proteins were identified by 4D label-free proteomics in R. haemaphysaloides and H. asiaticum hemolymph EVs, respectively. TEM and NTA revealed that the sizes of the hemolymph EVs from R. haemaphysaloides and H. asiaticum were 133 and 138 nm, respectively. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment analyses of identified proteins revealed pathways related to binding, catalytic and transporter activity, translation, transport and catabolism, signal transduction and cellular community. The key EV marker proteins RhCD9, RhTSG101, Rh14-3-3 and RhGAPDH were identified using proteomics and western blot. The presence of RhFerritin-2 in tick hemolymph EVs was confirmed by western blot and immuno-electron microscopy. We demonstrated that PKH26-labeled hemolymph EVs are internalized by tick salivary glands and ovary cells in vitro.

CONCLUSIONS

The results suggest that tick EVs are secreted into, and circulated by, the hemolymph. EVs may play roles in the regulation of tick development, metabolism and reproduction.

Localization of secreted ferritin (FER2) in the embryos of the tick Haemaphysalis longicornis

Despite the absence of a blood meal, embryogenesis involves many processes that require nutrients and other essential elements, including iron. Due to the lack of an external source of these nutrients, these requirements are acquired maternally. Because of the toxic nature of iron, they are transferred through iron transport molecules such as secreted ferritin (FER2). Here we tried to follow the trail of the FER2 through indirect immunofluorescence, and we observed an apparent shift of FER2 from the germ layer at the early part of development to the appendages during the late stage of embryogenesis. FER2 is also found in the middle part of the legs of the embryo. The apparent movement not only sheds light on iron processing events during embryogenesis but also indirectly guides organogenesis in the tick.

Observing How Glutathione and S-Hexyl Glutathione Bind to Glutathione S-Transferase from Rhipicephalus (Boophilus) microplus

Rhipicephalus (Boophilus) microplus is one of the most widespread ticks causing a massive loss to livestock production. The long-term use of acaracides rapidly develops acaracide resistance. In R. microplus, enhancing the metabolic activity of glutathione S-transferase (RmGST) is one of the mechanisms underlying acaracide resistance. RmGST catalyzes the conjugation of glutathione (GSH) to insecticides causing an easy-to-excrete conjugate. The active RmGST dimer contains two active sites (hydrophobic co-substrate binding site (H-site) and GSH binding site (G-site)) in each monomer. To preserve the insecticide efficacy, s-hexyl glutathione (GTX), a GST inhibitor, has been used as a synergist. To date, no molecular information on the RmGST-GSH/GTX complex is available. The insight is important for developing a novel RmGST inhibitor. Therefore, in this work, molecular dynamics simulations (MD) were performed to explore the binding of GTX and GSH to RmGST. GSH binds tighter and sits rigidly inside the G-site, while flexible GTX occupies both active sites. In GSH, the backbone mainly interacts with W8, R43, W46, K50, N59, L60, Q72, and S73, while its thiol group directs to Y7. In contrast, the aliphatic hexyl of GTX protrudes into the H-site and allows a flexible peptide core to form various interactions. Such high GTX flexibility and the protrusion of its hexyl moiety to the H-site suggest the dual role of GTX in preventing the conjugation reaction and the binding of acaracide. This insight can provide a better understanding of an important insecticide-resistance mechanism, which may in turn facilitate the development of novel approaches to tick control.

Ambivalent Roles of Oxidative Stress in Triangular Relationships among Arthropod Vectors, Pathogens and Hosts

Blood-feeding arthropods, particularly ticks and mosquitoes are considered the most important vectors of arthropod-borne diseases affecting humans and animals. While feeding on blood meals, arthropods are exposed to high levels of reactive oxygen species (ROS) since heme and other blood components can induce oxidative stress. Different ROS have important roles in interactions among the pathogens, vectors, and hosts. ROS influence various metabolic processes of the arthropods and some have detrimental effects. In this review, we investigate the various roles of ROS in these arthropods, including their innate immunity and the homeostasis of their microbiomes, that is, how ROS are utilized to maintain the balance between the natural microbiota and potential pathogens. We elucidate the mechanism of how ROS are utilized to fight off invading pathogens and how the arthropod-borne pathogens use the arthropods’ antioxidant mechanism to defend against these ROS attacks and their possible impact on their vector potentials or their ability to acquire and transmit pathogens. In addition, we describe the possible roles of ROS in chemical insecticide/acaricide activity and/or in the development of resistance. Overall, this underscores the importance of the antioxidant system as a potential target for the control of arthropod and arthropod-borne pathogens.

Virome and Blood Meal-Associated Host Responses in Ixodes persulcatus Naturally Fed on Patients

The long-lasting co-evolution of ticks with pathogens results in mutual adaptation. Blood-feeding is one of the critical physiological behaviors that have been associated with the tick microbiome; however, most knowledge was gained through the study of laboratory-reared ticks. Here we detached Ixodes persulcatus ticks at different stages of blood-feeding from human patients and performed high-throughput transcriptomic analysis on them to identify their virome and genes differentially expressed between flat and fully fed ticks. We also traced bloodmeal sources of those ticks and identified bats and three other potential mammalian hosts, highlighting the public health significance. We found Jingmen tick virus and 13 putative new viruses belonging to 11 viral families, three of which even exhibited high genetic divergence from viruses previously reported in the same tick species from the same geographic region. Furthermore, differential expression analysis suggested a downregulation of antioxidant genes in the fully fed I. persulcatus ticks, which might be related to bloodmeal-related redox homeostasis. Our work highlights the significance of active surveillance of tick viromes and suggests a role of reactive oxygen species (ROS) in modulating changes in the microbiome during blood-feeding.

Cloning, expression, and function of ferritins in the tick Haemaphysalis flava

The full-length cDNA of two ferritins of Haemaphysalis flava were cloned after which recombinant Hf-FER1 and Hf-FER2 were expressed and their function was analyzed. In addition, RNA interference (RNAi) based on the injection of Hf-fer1 or Hf-fer2 dsRNA into fully engorged female ticks was performed. The cDNA encoding Hf-FER1 is 834 bp in length. It contains an iron-responsive element in the 5' untranslated region and encodes 174 amino acid residues. The full-length cDNA of Hf-FER2 contains 696 bp and encodes 199 amino acids, including a putative signal peptide sequence. Hf-FER1 and Hf-FER2 both have the ferroxidase iron center and the ferrihydrite nucleation center. The evolutionary relationship of Hf-FER1 and Hf-FER2 was established, and the predicted quaternary structures were assembled as typical spherical shells composed of 24 subunits which was demonstrated by nature PAGE. Real-time PCR showed that Hf-fer1 and Hf-fer2 were expressed in all developmental stages, with the highest expression in fully engorged females. The expression of Hf-fer1 and Hf-fer2 were relatively high in unfed larvae. Hf-fer1 was expressed in all tissues and was especially abundant in the salivary glands of fully engorged females. In contrast, the highest levels of Hf-fer2 were found in the midgut of fully engorged females, and no expression was found in the salivary glands of this life stage. Both recombinant Hf-FER1 and Hf-FER2 had iron-binding capabilities. Silencing of both Hf-fer1 and Hf-fer2 affected fecundity. Compared to the control, the percentage of ticks that laid eggs in the Hf-fer1 and Hf-fer2 RNAi groups was 73.3% and 66.7%, respectively. The silenced ticks that laid eggs had lower egg weight to body weight ratios, and the eggs had abnormal morphologies. The hatchability of eggs with normal morphology in the Hf-fer1 and Hf-fer2 silenced groups was 47.8% and 22.8%, respectively, which was significantly different from the control group (P < 0.005). These findings indicate that Hf-FER1 and Hf-FER2 play important roles in the iron storage of H. flava.