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

A role of epigenetic mechanisms in regulating female reproductive responses to temperature in a pest beetle

Many species are threatened by climate change and must rapidly respond to survive in changing environments. Epigenetic modifications, such as DNA methylation, can facilitate plastic responses by regulating gene expression in response to environmental cues. Understanding epigenetic responses is therefore essential for predicting species' ability to rapidly adapt in the context of global environmental change. Here, we investigated the functional significance of different methylation-associated cellular processes on temperature-dependent life history in seed beetles, Callosobruchus maculatus Fabricius 1775 (Coleoptera: Bruchidae). We assessed changes under thermal stress in (1) DNA methyltransferase (Dnmt1 and Dnmt2) expression levels, (2) genome-wide methylation and (3) reproductive performance, with (2) and (3) following treatment with 3-aminobenzamide (3AB) and zebularine (Zeb) over two generations. These drugs are well-documented to alter DNA methylation across the tree of life. We found that Dnmt1 and Dnmt2 were expressed throughout the body in males and females, but were highly expressed in females compared with males and exhibited temperature dependence. However, whole-genome methylation did not significantly vary with temperature, and only marginally or inconclusively with drug treatment. Both 3AB and Zeb led to profound temperature-dependent shifts in female reproductive life history trade-off allocation, often increasing fitness compared with control beetles. Mismatch between magnitude of treatment effects on DNA methylation versus life history effects suggest potential of 3AB and Zeb to alter reproductive trade-offs via changes in DNA repair and recycling processes, rather than or in addition to (subtle) changes in DNA methylation. Together, our results suggest that epigenetic mechanisms relating to Dnmt expression, DNA repair and recycling pathways, and possibly DNA methylation, are strongly implicated in modulating insect life history trade-offs in response to temperature change.

DNA methylase 1 influences temperature responses and development in the invasive pest Tuta absoluta

DNA methylase 1 (Dnmt1) is an important regulatory factor associated with biochemical signals required for insect development. It responds to changes in the environment and triggers phenotypic plasticity. Meanwhile, Tuta absoluta Meyrick (Lepidoptera: Gelechiidae)-a destructive invasive pest-can rapidly invade and adapt to different habitats; however, the role of Dnmt1 in this organism has not been elucidated. Accordingly, this study investigates the mechanism(s) underlying the rapid adaptation of Tuta absoluta to temperature stress. Potential regulatory genes were screened via RNAi (RNA interference), and the DNA methylase in Tuta absoluta was cloned by RACE (Rapid amplification of cDNA ends). TaDnmt1 was identified as a potential regulatory gene via bioinformatics; its expression was evaluated in response to temperature stress and during different development stages using real-time polymerase chain reaction. Results revealed that TaDnmt1 participates in hot/cold tolerance, temperature preference and larval development. The full-length cDNA sequence of TaDnmt1 is 3765 bp and encodes a 1254 kDa protein with typical Dnmt1 node-conserved structural features and six conserved DNA-binding active motifs. Moreover, TaDnmt1 expression is significantly altered by temperature stress treatments and within different development stages. Hence, TaDnmt1 likely contributes to temperature responses and organismal development. Furthermore, after treating with double-stranded RNA and exposing Tuta absoluta to 35°C heat shock or -12°C cold shock for 1 h, the survival rate significantly decreases; the preferred temperature is 2°C lower than that of the control group. In addition, the epidermal segments become enlarged and irregularly folded while the surface dries up. This results in a significant increase in larval mortality (57%) and a decrease in pupation (49.3%) and eclosion (50.9%) rates. Hence, TaDnmt1 contributes to temperature stress responses and temperature perception, as well as organismal growth and development, via DNA methylation regulation. These findings suggest that the rapid geographic expansion of T absoluta has been closely associated with TaDnmt1-mediated temperature tolerance. This study advances the research on 'thermos Dnmt' and provides a potential target for RNAi-driven regulation of Tuta absoluta.

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.

Potential Role of DNA Methylation as a Driver of Plastic Responses to the Environment Across Cells, Organisms, and Populations

There is great interest in exploring epigenetic modifications as drivers of adaptive organismal responses to environmental change. Extending this hypothesis to populations, epigenetically driven plasticity could influence phenotypic changes across environments. The canonical model posits that epigenetic modifications alter gene regulation and subsequently impact phenotypes. We first discuss origins of epigenetic variation in nature, which may arise from genetic variation, spontaneous epimutations, epigenetic drift, or variation in epigenetic capacitors. We then review and synthesize literature addressing three facets of the aforementioned model: (i) causal effects of epigenetic modifications on phenotypic plasticity at the organismal level, (ii) divergence of epigenetic patterns in natural populations distributed across environmental gradients, and (iii) the relationship between environmentally induced epigenetic changes and gene expression at the molecular level. We focus on DNA methylation, the most extensively studied epigenetic modification. We find support for environmentally associated epigenetic structure in populations and selection on stable epigenetic variants, and that inhibition of epigenetic enzymes frequently bears causal effects on plasticity. However, there are pervasive confounding issues in the literature. Effects of chromatin-modifying enzymes on phenotype may be independent of epigenetic marks, alternatively resulting from functions and protein interactions extrinsic of epigenetics. Associations between environmentally induced changes in DNA methylation and expression are strong in plants and mammals but notably absent in invertebrates and nonmammalian vertebrates. Given these challenges, we describe emerging approaches to better investigate how epigenetic modifications affect gene regulation, phenotypic plasticity, and divergence among populations.

De novo genome assembly of the invasive mosquito species Aedes japonicus and Aedes koreicus

BACKGROUND

Recently, two invasive Aedes mosquito species, Ae. japonicus and Ae. koreicus, are circulating in several European countries posing potential health risks to humans and animals. Vector control is the main option to prevent mosquito-borne diseases, and an accurate genome sequence of these mosquitoes is essential to better understand their biology and to develop effective control strategies.

METHODS

A de novo genome assembly of Ae. japonicus (Ajap1) and Ae. koreicus (Akor1) has been produced based on a hybrid approach that combines Oxford Nanopore long-read and Illumina short-read data. Their quality was ascertained using various metrics. Masking of repetitive elements, gene prediction and functional annotation was performed.

RESULTS

Sequence analysis revealed a very high presence of repetitive DNA and, among others, thermal adaptation genes and insecticide-resistance genes. Through the RNA-seq analysis of larvae and adults of Ae. koreicus and Ae. japonicus exposed to different temperatures, we also identified genes showing a differential temperature-dependent activation.

CONCLUSIONS

The assembly of Akor1 and Ajap1 genomes constitutes the first updated collective knowledge of the genomes of both mosquito species, providing the possibility of understanding key mechanisms of their biology such as the ability to adapt to harsh climates and to develop insecticide-resistance mechanisms.

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.

Adapting to climate with limited genetic diversity: Nucleotide, DNA methylation and microbiome variation among populations of the social spider Stegodyphus dumicola

Understanding the role of genetic and nongenetic variants in modulating phenotypes is central to our knowledge of adaptive responses to local conditions and environmental change, particularly in species with such low population genetic diversity that it is likely to limit their evolutionary potential. A first step towards uncovering the molecular mechanisms underlying population-specific responses to the environment is to carry out environmental association studies. We associated climatic variation with genetic, epigenetic and microbiome variation in populations of a social spider with extremely low standing genetic diversity. We identified genetic variants that are associated strongly with environmental variation, particularly with average temperature, a pattern consistent with local adaptation. Variation in DNA methylation in many genes was strongly correlated with a wide set of climate parameters, thereby revealing a different pattern of associations than that of genetic variants, which show strong correlations to a more restricted range of climate parameters. DNA methylation levels were largely independent of cis-genetic variation and of overall genetic population structure, suggesting that DNA methylation can work as an independent mechanism. Microbiome composition also correlated with environmental variation, but most strong associations were with precipitation-related climatic factors. Our results suggest a role for both genetic and nongenetic mechanisms in shaping phenotypic responses to local environments.

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.

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.