Cold tolerance and biochemical response of unfed Dermacentor silvarum ticks to low temperature

Molecular characterization and induced changes of histone acetyltransferases in the tick Haemaphysalis longicornis in response to cold stress

BACKGROUND

Epigenetic modifications of histones play important roles in the response of eukaryotic organisms to environmental stress. However, many histone acetyltransferases (HATs), which are responsible for histone acetylation, and their roles in mediating the tick response to cold stress have yet to be identified. In the present study, HATs were molecularly characterized and their associations with the cold response of the tick Haemaphysalis longicornis explored.

METHODS

HATs were characterized by using polymerase chain reaction (PCR) based on published genome sequences, followed by multiple bioinformatic analyses. The differential expression of genes in H. longicornis under different cold treatment conditions was evaluated using reverse transcription quantitative PCR (RT-qPCR). RNA interference was used to explore the association of HATs with the cold response of H. longicornis.

RESULTS

Two HAT genes were identified in H. longicornis (Hl), a GCN5-related N-acetyltransferase (henceforth HlGNAT) and a type B histone acetyltransferase (henceforth HlHAT-B), which are respectively 960 base pairs (bp) and 1239 bp in length. Bioinformatics analysis revealed that HlGNAT and HlHAT-B are unstable hydrophilic proteins characterized by the presence of the acetyltransferase 16 domain and Hat1_N domain, respectively. RT-qPCR revealed that the expression of HlGNAT and HlHAT-B decreased after 3 days of cold treatment, but gradually increased with a longer period of cold treatment. The mortality rate following knockdown of HlGNAT or HlHAT-B by RNA interference, which was confirmed by RT-qPCR, significantly increased (P < 0.05) when H. longicornis was treated at the lowest lethal temperature (- 14 °C) for 2 h.

CONCLUSIONS

The findings demonstrate that HATs may play a crucial role in the cold response of H. longicornis. Thus further research is warranted to explore the mechanisms underlying the epigenetic regulation of the cold response in ticks.

Molecular characterization and modulated expression of histone acetyltransferases during cold response of the tick Dermacentor silvarum (Acari: Ixodidae)

BACKGROUND

Histone acetylation is involved in the regulation of stress responses in multiple organisms. Dermacentor silvarum is an important vector tick species widely distributed in China, and low temperature is a crucial factor restricting the development of its population. However, knowledge of the histone acetyltransferases and epigenetic mechanisms underlying cold-stress responses in this tick species is limited.

METHODS

Histone acetyltransferase genes were characterized in D. silvarum, and their relative expressions were determined using qPCR during cold stress. The association and modulation of histone acetyltransferase genes were further explored using RNA interference, and both the H3K9 acetylation level and relative expression of KAT5 protein were evaluated using western blotting.

RESULTS

Three histone acetyltransferase genes were identified and named as DsCREBBP, DsKAT6B, and DsKAT5. Bioinformatics analysis showed that they were unstable hydrophilic proteins, characterized by the conserved structures of CBP (ZnF_TAZ), PHA03247 super family, Creb_binding, and MYST(PLN00104) super family. Fluorescence quantitative PCR showed that the expression of DsCREBBP, DsKAT6B, and DsKAT5 increased after 3 days of cold treatment, with subsequent gradual decreases, and was lowest on day 9. Western blotting showed that both the H3K9 acetylation level and relative expression of KAT5 in D. silvarum increased after treatment at - 4, 4, and 8 °C for 3 and 6 days, whereas they decreased significantly after a 9-day treatment. RNA interference induced significant gene silencing, and the mortality rate of D. silvarum significantly increased at the respective semi-lethal temperatures.

CONCLUSION

These results imply that histone acetyltransferases play an important role in tick adaptation to low temperatures and lay a foundation for further understanding of the epigenetic regulation of histone acetylation in cold-stressed ticks. Further research is needed to elucidate the mechanisms underlying histone acetylation during cold stress in ticks.

Protein regulation mechanism of cold tolerance in Haemaphysalis longicornis

Ticks are external parasitic arthropods that can transmit a variety of pathogens by sucking blood. Low-temperature tolerance is essential for ticks to survive during the cold winter. Exploring the protein regulation mechanism of low-temperature tolerance of Haemaphysalis longicornis could help to explain how ticks survive in winter. In this study, the quantitative proteomics of several tissues of H. longicornis exposed to low temperature were studied by data independent acquisition technology. Totals of 3 699, 3 422, and 1 958 proteins were identified in the salivary gland, midgut, and ovary, respectively. The proteins involved in energy metabolism, cell signal transduction, protein synthesis and repair, and cytoskeleton synthesis changed under low-temperature stress. The comprehensive analysis of the protein regulation of multiple tissues of female ticks exposed to low temperature showed that maintaining cell homeostasis, maintaining cell viability, and enhancing cell tolerance were the most important means for ticks to maintain vital signs under low temperature. The expression of proteins involved in and regulating the above cell activities was the key to the survival of ticks under low temperatures. Through the analysis of a large amount of data, we found that the expression levels of arylamine N-acetyltransferase, inositol polyphosphate multikinase, and dual-specificity phosphatase were up-regulated under low temperature. We speculated that they might have important significance in low-temperature tolerance. Then, we performed RNA interference on the mRNA of these 3 proteins, and the results showed that the ability of female ticks to tolerate low temperatures decreased significantly.

DNA Methylation Variation Is a Possible Mechanism in the Response of Haemaphysalis longicornis to Low-Temperature Stress

Ticks are notorious ectoparasites and transmit the greatest variety of pathogens than any other arthropods. Cold tolerance is a key determinant of tick abundance and distribution. While studies have shown that DNA methylation is one of the important epigenetic regulations found across many species and plays a significant role in their response to low-temperature stress, its role in the response of ticks to low-temperature stress remains unexplored. Herein, we explored the DNA methylation profile of the tick, Haemaphysalis longicornis, exposed to low-temperature stress (4 °C) using whole-genome bisulfite sequencing (WGBS). We found that approximately 0.95% and 0.94% of the genomic C sites were methylated in the control and low-temperature groups, respectively. Moreover, the methylation level under the CG context was about 3.86% and 3.85% in the control and low-temperature groups, respectively. In addition, a total of 6087 differentially methylated regions (DMRs) were identified between the low-temperature and control groups, including 3288 hypermethylated and 2799 hypomethylated DMRs. Further, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially methylated genes revealed that most of the DMGs were significantly enriched in binding and RNA transport pathways. Taken together, this research confirmed, for the first time, the whole genome DNA methylation profile of H. longicornis and provided new insights into the DNA methylation changes relating to low-temperature stress in H. longicornis, as well as provided a foundation for future studies on the epigenetic mechanism underlying the responses of ticks to abiotic stress.

miR-2a and miR-279 are functionally associated with cold tolerance in Dermacentor silvarum (Acari: Ixodidae)

Ticks are obligate blood-sucking ectoparasites that can attack mammals, birds, reptiles as well as amphibians. Dermacentor silvarum, an important vector of various pathogenic bacteria, viruses, and protozoans, is widely distributed in China. MicroRNAs (miRNAs) are ~22 nucleotide non-coding small RNA molecules, involved in the regulation of various physiological and cellular processes. Previous studies demonstrated the vital roles of miRNAs during the reproduction and development of ticks, whereas, the regulatory/functional roles of microRNAs during the cold response of ticks remain unexplored. Here, we identified and functionally explored D. silvarum miRNAs involved in cold response to gain further understanding of the molecular regulatory mechanisms underlying cold stress in ticks. The microRNA libraries of D. silvarum were established via high-throughput sequencing after exposure to different cold treatments. A total of 147 miRNAs, including 44 known miRNAs and 103 new miRNAs, were identified. The verification of six highly differentially expressed miRNAs (miR-2a, miR-5305, miR-7, miR-279, miR-993, and novel-3) via RT-qPCR were consistent with the high-throughput sequence results. miR-2a peaked by day 6 and miR-279 expression was lowest by day 3 after cold treatment. The potential target genes of miR-2a and miR-279 were the glycogen phosphorylase (GPase) gene and serine gene, respectively. After injecting D. silvarum ticks with miR-2a and miR-279 antagonists, their respective target genes were up-regulated and vice-versa after injection with the agonists. These results indicated that these two miRNAs and their target genes may be involved in the cold response of D. silvarum ticks.

Metabolomic and transcriptomic responses of ticks during recovery from cold shock reveal mechanisms of survival

Ticks are blood-feeding ectoparasites but spend most of their life off-host where they may have to tolerate low winter temperatures. Rapid cold-hardening (RCH) is a process commonly used by arthropods, including ticks, to improve survival of acute low temperature exposure. However, little is known about the underlying mechanisms in ticks associated with RCH, cold shock, and recovery from these stresses. In the present study, we investigated the extent to which RCH influences gene expression and metabolism during recovery from cold stress in Dermacentor variabilis, the American dog tick, using a combined transcriptomics and metabolomics approach. Following recovery from RCH, 1,860 genes were differentially expressed in ticks, whereas only 99 genes responded during recovery to direct cold shock. Recovery from RCH resulted in an upregulation of various pathways associated with ion binding, transport, metabolism, and cellular structures seen in the response of other arthropods to cold. The accumulation of various metabolites, including several amino acids and betaine, corresponded to transcriptional shifts in the pathways associated with these molecules, suggesting congruent metabolome and transcriptome changes. Ticks receiving exogenous betaine and valine demonstrated enhanced cold tolerance, suggesting cryoprotective effects of these metabolites. Overall, many of the responses during recovery from cold shock in ticks were similar to those observed in other arthropods, but several adjustments may be distinct from other currently examined taxa.

Functional implication of heat shock protein 70/90 and tubulin in cold stress of Dermacentor silvarum

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

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05056-y.

DNA Methyltransferases Contribute to Cold Tolerance in Ticks Dermacentor silvarum and Haemaphysalis longicornis (Acari: Ixodidae)

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.

Cold hardening improves larval tick questing under low temperatures at the expense of longevity

Questing in ticks is essential for locating a host, and this behavioral response can occur at regionally specific low temperatures for most tick species. Little is known about the dynamics between tick questing behavior and temperature in ticks, specifically how this may impact other aspects of tick biology. Here, we examine whether cold hardening increases questing in three larval tick species (Ixodes uriae, Dermacentor variabilis, and Amblyomma americanum) at low temperatures and whether cold hardening impacts longevity. Rapid cold hardening and prolonged cold acclimation benefitted ticks by decreasing the temperature of chill coma onset, and increased survival, activity, and questing in ticks at low temperatures. Oxygen consumption increased at low temperatures following acclimation in larvae, suggesting this process has a distinct metabolic expense. This increased metabolism associated with hardening led to a substantial reduction in larval longevity as nutrient reserves are limited and cannot be replenished until a host is located. These studies suggest that tick larvae, and likely other developmental stages, require a delicate balance between the need for questing at low temperatures and survival until the first blood meal.

Distribution of Dermacentor silvarum and Associated Pathogens: Meta-Analysis of Global Published Data and a Field Survey in China

Dermacentor silvarum is an obligate blood sucking arthropod and transmits various pathogens to humans and domestic animals. Recently several new viruses were detected in D. silvarum as an emerging disease threat. In this study, we aimed to analyze its geographical distribution and associated pathogens. Data were collected from multiple sources, including a field survey, reference book, and literature review. We searched various electronic databases with the terms “Dermacentor silvarum” OR “D. silvarum” for studies published since 1963 and the positive rates for Dermacentor silvarum-associated pathogens were estimated by meta-analysis. D. silvarum was found only in four countries in Eurasia, ranging from 22° N to 57° N latitude. At least 20 human pathogens were associated with D. silvarum, including five species of spotted fever group rickettsiae, three species in the family of Anaplasmataceae, three genospecies in the complex Borrelia burgdorferi sensu lato, Francisella tularensis, Babesia venatorum, Coxiella buenetii, Borrelia miyamotoi, and five species of virus. Among them, Rickettsia raoultii was widely detected in D. silvarum, showing the highest pooled positive rate (25.15%; 95% CI 13.31–39.27). Our work presents the most comprehensive data and analysis (to our knowledge) for the geographical distribution of D. silvarum and associated pathogens, revealing an emerging threat to public health and stocking farming. Continued surveillance and further investigations should be enhanced.

Projecting the potential distribution of ticks in China under climate and land use change

Ticks are known as vectors of several pathogens causing various human and animal diseases including Lyme borreliosis, tick-borne encephalitis, and Crimean-Congo hemorrhagic fever. While China is known to have more than 100 tick species well distributed over the country, our knowledge on the likely distribution of ticks in the future remains very limited, which hinders the prevention and control of the risk of tick-borne diseases. In this study, we selected four representative tick species which have different regional distribution foci in mainland China. i.e., Dermacentor marginatus, Dermacentor silvarum, Haemaphysalis longicornis and Ixodes granulatus. We used the MaxEnt model to identify the key environmental factors of tick occurrence and map their potential distributions in 2050 under four combined climate and socioeconomic scenarios (i.e., SSP1-RCP2.6, SSP2-RCP4.5, SSP3-RCP7.0 and SSP5-RCP8.5). We found that the extent of the urban fabric, cropland and forest, temperature annual range and precipitation of the driest month were the main determinants of the potential distributions of the four tick species. Under the combined scenarios, with climate warming, the potential distributions of ticks shifted to further north in China. Due to a decrease in the extent of forest, the distribution probability of ticks declined in central and southern China. In contrast with previous findings on an estimated amplification of tick distribution probability under the extreme emission scenario (RCP8.5), our studies projected an overall reduction in the distribution probability under RCP8.5, owing to an expected effect of land use. Our results could provide new data to help identify the emerging risk areas, with amplifying suitability for tick occurrence, for the prevention and control of tick-borne zoonoses in mainland China. Future directions are suggested towards improved quantity and quality of the tick occurrence database, comprehensiveness of factors and integration of different modelling approaches, and capability to model pathogen spillover at the human-tick interface.

Supercooling Points of Adult Dermacentor variabilis (Acari: Ixodidae) From a Population Near the Northern Distribution Limit

The northern distributional limit of Dermacentor variabilis Say, the American dog tick, is expanding in Saskatchewan and Manitoba (western Canada). The ability of D. variabilis to continue to expand its range northwards will depend upon the ability of individuals within populations at the species distributional edge to withstand very low temperatures during winter. One component of cold hardiness is the supercooling point (SCP), the temperature below 0°C at which an individual freezes. In this study, the SCP was determined for 94 questing D. variabilis adults (44 females and 50 males) from an established population near Blackstrap Provincial Park in Saskatchewan. SCP values ranged from -18.2 to -6.7°C, with a median of -13.3°C. This suggests that host-seeking D. variabilis adults differ in their ability to survive exposure to subzero temperatures, for at least a short period of time, without freezing. The distribution of SCPs was bimodal, but there was no significant difference in SCP values between female and male ticks, and no relationship between SCP and tick body weight. It remains to be determined what factors contribute to the variation in SCP values among questing D. variabilis adults.

Vectors of disease at the northern distribution limit of the genus Dermacentor in Eurasia: D. reticulatus and D. silvarum

The two ixodid tick species Dermacentor reticulatus (Fabricius) and Dermacentor silvarum Olenev occur at the northern distribution limit of the genus Dermacentor in Eurasia, within the belt of [Formula: see text] latitude. Whilst the distribution area of D. reticulatus extends from the Atlantic coast of Portugal to Western Siberia, that of D. silvarum extends from Western Siberia to the Pacific coast. In Western Siberia, the distribution areas of the two Dermacentor species overlap. Although the two tick species are important vectors of disease, detailed information concerning the entire distribution area, climate adaptation, and proven vector competence is still missing. A dataset was compiled, resulting in 2188 georeferenced D. reticulatus and 522 D. silvarum locations. Up-to-date maps depicting the geographical distribution and climate adaptation of the two Dermacentor species are presented. To investigate the climate adaptation of the two tick species, the georeferenced locations were superimposed on a high-resolution map of the Köppen-Geiger climate classification. The frequency distribution of D. reticulatus under different climates shows two major peaks related to the following climates: warm temperate with precipitation all year round (57%) and boreal with precipitation all year round (40%). The frequency distribution of D. silvarum shows also two major peaks related to boreal climates with precipitation all year round (30%) and boreal winter dry climates (60%). Dermacentor silvarum seems to be rather flexible concerning summer temperatures, which can range from cool to hot. In climates with cool summers D. reticulatus does not occur, it prefers warm and to a lesser extent hot summers. Lists are given in this paper for cases of proven vector competence for various agents of both Dermacentor species. For the first time, the entire distribution areas of D. reticulatus and D. silvarum were mapped using georeferenced data. Their climate adaptations were quantified by Köppen profiles.

Ticks (Acari: Ixodoidea) in China: Geographical distribution, host diversity, and specificity

Ticks are obligate blood-sucking ectoparasites, which not only directly damage through bites but also transmit many pathogens. China has a high diversity of tick species, 125 species have been reported, including 111 hard tick and 14 soft tick species. Many of the ticks are important vectors of pathogens, resulting in zoonoses. The dynamics of ticks are affected by both the host and habitat environment. However, systematic studies on the geographical distribution, host diversity, and specificity of ticks are limited in China. To achieve this goal, the relevant available data were summarized and analyzed in this study. Ticks are distributed in all parts of China and Xinjiang has the most records of ticks. The distribution of ticks in adjacent areas is similar, indicating that the habitat environment affects their distribution. Most ticks are widely distributed, whereas some species are endemic to their distributed regions. Ticks are parasitic on mammals, birds, and reptiles, of which mammals are the main host species. Overall, most ticks parasitize different hosts, only a few ticks have strict host specificity, such as ticks that are specifically parasitic on reptiles and bats. In addition, environmental changes and control efforts also influence the dynamics of ticks. These results can better reveal tick biological traits and are valuable for tick control.