DNA methylation, mediated by DNA methyltransferases (Dnmts), is a typical epigenetic process that plays an important role in affecting organism acclimatization and adaptation to environmental changes. However, information about Dnmts and their associations with the cold tolerance of ticks remains meager. Hence, in the present study, the Dnmts in important vector ticks Dermacentor silvarum and Haemaphysalis longicornis were cloned and identified, and their functions in cold response were further explored. Results showed that the length of DsDnmt and DsDnmt1 in D. silvarum, and HlDnmt1 and HlDnmt in H. longicornis were 1,284, 549, 1,500, and 1,613 bp, respectively. Bioinformatics in protein analysis revealed that they were all unstable hydrophilic proteins and were mainly characterized with Dcm (DNA cytosine methyltransferase domain), Dnmt1-RFD (DNA methyltransferase replication foci domain), zf-CXXC (zinc finger-CXXC domain), and BAH (Bromo adjacent homology domain). The relative expression of these Dnmts was reduced after cold treatment for 3 days (P < 0.05), and increased with the extension of treatment. Western blot revealed that Dnmt1 decreased first and then increased significantly (P < 0.05) in both tick species, whereas other Dnmts fluctuated at varying degrees. RNA interference significantly silenced the genes Dnmts (P < 0.01), and mortality increased significantly (P < 0.05), when exposed to sub-lethal temperature, underscoring the important roles of Dnmts during the cold response of D. silvarum and H. longicornis. The above results lay the foundation for further understanding of the epigenetic regulation of DNA methylation in cold acclimatization and adaptation of ticks.

Background The tick Dermacentor silvarum Olenev (Acari: Ixodidae) is a vital vector tick species mainly distributed in the north of China and overwinters in the unfed adult stage. The knowledge of the mechanism that underlies its molecular adaptation against cold is limited. In the present study, genes of hsp70 and hsp90 cDNA, named Dshsp70 and Dshsp90, and tubulin were cloned and characterized from D. silvarum, and their functions in cold stress were further evaluated. Methods The genome of the heat shock proteins and tubulin of D. silvarum were sequenced and analyzed using bioinformatics methods. Each group of 20 ticks were injected in triplicate with Dshsp90-, Dshsp70-, and tubulin-derived dsRNA, whereas the control group was injected with GFP dsRNA. Then, the total RNA was extracted and cDNA was synthesized and subjected to RT-qPCR. After the confirmation of knockdown, the ticks were incubated for 24 h and were exposed to − 20 °C lethal temperature (LT50), and then the mortality was calculated. Results Results indicated that Dshsp70 and Dshsp90 contained an open reading frame of 345 and 2190 nucleotides that encoded 114 and 729 amino acid residues, respectively. The transcript Dshsp70 showed 90% similarity with that identified from Dermacentor variabilis, whereas Dshsp90 showed 85% similarity with that identified from Ixodes scapularis. Multiple sequence alignment indicates that the deduced amino acid sequences of D. silvarum Hsp90, Hsp70, and tubulin show very high sequence identity to their corresponding sequences in other species. Hsp90 and Hsp70 display highly conserved and signature amino acid sequences with well-conserved MEEVD motif at the C-terminal in Hsp90 and a variable C-terminal region with a V/IEEVD-motif in Hsp70 that bind to numerous co-chaperones. RNA interference revealed that the mortality of D. silvarum was significantly increased after injection of dsRNA of Dshsp70 (P = 0.0298) and tubulin (P = 0.0448), whereas no significant increases were observed after the interference of Dshsp90 (P = 0.0709). Conclusions The above results suggested that Dshsp70 and tubulin play an essential role in the low-temperature adaptation of ticks. The results of this study can contribute to the understanding of the survival and acclimatization of overwintering ticks. Graphical abstract

No abstract

The tick Haemaphysalis longicornis has two reproductive groups: a bisexual group (HLBP) and a parthenogenetic group (HLPP). The comparative molecular regulation of ovarian development in these two groups is unexplored. We conducted transcriptome sequencing and quantitative proteomics on the ovaries of HLBP and HLPP, in different feeding stages, to evaluate the molecular function of genes associated with ovarian development. The ovarian tissues of HLBP and HLPP were divided into three feeding stages (early-fed, partially-fed and engorged). A total of 87,233 genes and 2,833 proteins were annotated in the ovary of H. longicornis in the different feeding stages. The differentially expressed genes (DEGs) of functional pathway analysis indicated that Lysosome, MAPK Signaling Pathway, Phagosome, Regulation of Actin Cytoskeleton, Endocytosis, Apoptosis, Insulin Signaling Pathway, Oxidative Phosphorylation, and Sphingolipid Metabolism were most abundant in the ovary of H. longicornis in the different feeding stages. Comparing the DEGs between HLBP and HLPP revealed that the ABC Transporter, PI3K-Akt Signaling Pathway and cAMP Signaling Pathway were the most enriched and suggested that the functions of signal transduction mechanisms may have changed during ovarian development. The functions of the annotated proteome of ovarian tissues were strongly correlated with the transcriptome annotation results, and these were further validated using quantitative polymerase chain reaction (qPCR). In the HLBP, the expression of cathepsin L, secreted proteins and glycosidase proteins was significantly up-regulated during feeding stages. In the HLPP, the lysozyme, yolk proteins, heat shock protein, glutathione S transferase, myosin and ATP synthase proteins were up-regulated during feeding stages. The significant differences of the gene expression between HLBP and HLPP indicated that variations in the genetic background and molecular function might exist in the two groups. These results provide a foundation for understanding the molecular mechanism and exploring the functions of genes in the ovarian development of different reproductive groups of H. longicornis.

The tick Dermacentor silvarum mainly overwinters in northern China as an unfed adult entering diapause; however, little is known about the involvement of genes and pathways in summer diapause and overwintering adults of this species. Hence, we investigated the differential transcriptome profile of D. silvarum ticks’ diapause and cold-responsive genes and pathways during its summer diapause, overwintering, and questing/active stages. Through transcriptomic sequencing, a total of 136,740 unigenes were assembled, and 1,549, 2,196, and 2,125 differentially expressed genes (DEG) were identified in active and diapause, diapause and overwintering, and active and overwintering ticks, respectively. Between active and diapause ticks, 875 genes were up-regulated, while 674 genes were found down-regulated in diapause ticks. When compared with active ticks, overwintering ticks had 993 up-regulated genes and 1,132 down-regulated genes. Comparison between diapause and overwintering ticks revealed 745 up-regulated and 1,451 down-regulated genes in overwintering ticks. Gene Ontology analysis revealed that, among the active, diapause, and overwintering ticks, most of the DEGs in molecular function were enriched in catalytic activity and hydrolase activity. In the cellular component, most DEGs were assigned to the integral component of the cell membrane , whereas the oxidation-reduction process was the most enriched among biological processes. In addition, the lysosome pathway was the most enriched pathway identified in KEGG pathway analysis. The above results deepened our understanding of the mechanism underlying diapause and overwintering adaptation of ticks, which is important for integrative prevention and control of ticks and tick-borne diseases.