Genome-wide DNA methylomes from discrete developmental stages reveal the predominance of non-CpG methylation in Tribolium castaneum

The emerging role of epigenetic mechanisms in insect defense against pathogens

Insects resist infection by natural selection that favors the survival and reproduction of the fittest phenotypes. Although the genetic mechanisms mediating the evolution of insect resistance have been investigated, little is known about the contribution of epigenetic mechanisms. Gene expression in response to a pathogen selection pressure is regulated by different mechanisms affecting chromatin plasticity. Whether transgenerational inheritance of genome-wide epigenetic marks contributes to the heritable manifestation of insect resistance is presently debated. Here, we review the latest works on the contributions of chromatin remodeling to insect immunity and adaptation to pathogens. We highlight DNA methylation, histone acetylation, and microRNAs in mediating the transgenerational inherited transcriptional reprogramming of defense-related gene expression and the evolution of insect resistance.

A review on the DNA methyltransferase family of insects: Aspect and prospects

The DNA methyltransferase family contains a conserved set of DNA-modifying enzymatic proteins. They are responsible for epigenetic gene modulation, such as transcriptional silencing, transcription activation, and post-transcriptional modulation. Recent research has revealed that the canonical DNA methyltransferases (DNMTs) biological roles go beyond their traditional functions of establishing and maintaining DNA methylation patterns. Although a complete DNA methylation toolkit is absent in most insect orders, recent evidence indicates the de novo DNA methylation and maintenance function remain conserved. Studies using various molecular approaches provided evidence that DNMTs are multi-functional proteins. However, still in-depth studies on their biological role lack due to the least studied area in insects. Here, we review the DNA methylation toolkit of insects, focusing on recent research on various insect orders, which exhibit DNA methylation at different levels, and for which DNMTs functional studies have become available in recent years. We survey research on the potential roles of DNMTs in the regulation of gene transcription in insect species. DNMTs participate in different physiological processes by interacting with other epigenetic factors. Future studies on insect's DNMTs will benefit to understand developmental processes, responses to various stimuli, and adaptability of insects to different environmental conditions.

Complete Metamorphosis in Manduca sexta Involves Specific Changes in DNA Methylation Patterns

The transition between morphologically distinct phenotypes during complete metamorphosis in holometabolous insects is accompanied by fundamental transcriptional reprogramming. Using the tobacco hornworm (Manduca sexta), a powerful model for the analysis of insect evolution and development, we conducted a genome-wide comparative analysis of gene expression and DNA methylation in caterpillars and adults to determine whether complete metamorphosis has an epigenetic basis in this species. Bisulfite sequencing indicated a generally low level of DNA methylation with a unimodal CpGO/E distribution. Expression analysis revealed that 24 % of all known M. sexta genes (3.729) were upregulated in last-instar larvae relative to the adult moth, whereas 26 % (4.077) were downregulated. We also identified 4.946 loci and 4.960 regions showing stage-specific differential methylation. Interestingly, genes encoding histone acetyltransferases and histone deacetylases were differentially methylated in the larvae and adults, indicating there is crosstalk between different epigenetic mechanisms. The distinct sets of methylated genes in M. sexta larvae and adults suggest that complete metamorphosis involves epigenetic modifications associated with profound transcriptional reprogramming, involving approximately half of all the genes in this species.

DNMT1 Promotes Genome Methylation and Early Embryo Development in Cockroaches

The influence of DNA methylation on gene behavior and its consequent phenotypic effects appear to be very important, but the details are not well understood. Insects offer a diversity of DNA methylation modes, making them an excellent lineage for comparative analyses. However, functional studies have tended to focus on quite specialized holometabolan species, such as wasps, bees, beetles, and flies. Here, we have studied DNA methylation in the hemimetabolan insect Blattella germanica. In this cockroach, a gene involved in DNA methylation, DNA methyltransferase 1 (DNMT1), is expressed in early embryogenesis. In our experiments, RNAi of DNMT1 reduces DNA methylation and impairs blastoderm formation. Using reduced representation bisulfite sequencing and transcriptome analyses, we observed that methylated genes are associated with metabolism and are highly expressed, whereas unmethylated genes are related to signaling and show low expression. Moreover, methylated genes show greater expression change and less expression variability than unmethylated genes.

Genome-Wide Differential DNA Methylation in Reproductive, Morphological, and Visual System Differences Between Queen Bee and Worker Bee (Apis mellifera)

There are many differences in external morphology and internal physiology between the Apis mellifera queen bee and worker bee, some of which are relevant to beekeeping production. These include reproductive traits, body size, royal jelly secreting properties, and visual system development, among others. The identification of candidate genes that control the differentiation of these traits is critical for selective honeybee breeding programs. In this study, we compared the genomic methylation of queen bee and worker bee larvae at 3, 4, and 5 days of age by whole-genome bisulfite sequencing, and found that the basic characteristics of genomic methylation in queen and worker larvae were the same. There were approximately 49 million cytosines in the Apis larvae genome, of which about 90,000 were methylated. Methylated CpG sites accounted for 99% of the methylated cytosines, and methylation mainly occurred in exons. However, methylation levels of queen and worker larvae showed different trends with age: the methylation level of queen larvae varied with age in an inverted parabola, while the corresponding trend for worker larvae with resembled an exponential curve with a platform. The methylation level of queen larvae was higher than that of worker larvae at 3 days of age, lower than that of worker larvae at 4 days of age, and similar to that of worker larvae at 5 days old. The top 10 differentially methylated genes (DMGs) and 13 caste-specific methylated genes were listed, and correlations with caste determination were speculated. We additionally screened 38 DMGs between queen larvae and worker larvae involved in specific organ differentiation as well as reproduction, morphology, and vision differentiation during caste determination. These genes are potential molecular markers for selective breeding of A. mellifera to improve fecundity, royal jelly production, body size, and foraging, and represent candidate genes for investigating specialized functional segregation during the process of caste differentiation.

Functional diversification of three delta-class glutathione S-transferases involved in development and detoxification in Tribolium castaneum

Glutathione S-transferases (GSTs) are members of a multifunctional enzyme superfamily. Forty-one GSTs have been identified in Tribolium castaneum; however, none of the 41 GSTs has been functionally characterized. Here, three delta-class GSTs, TcGSTd1, TcGSTd2 and TcGSTd3, of T. castaneum were successfully cloned and expressed in Escherichia coli. All of the studied GSTs catalysed the conjugation of reduced glutathione with 1-chloro-2,4-dinitrobenzene. Insecticide treatment showed that the expression levels of TcGSTd3 and TcGSTd2 were significantly increased after exposure to phoxim and lambda-cyhalothrin, whereas TcGSTd1 was slightly upregulated only in response to phoxim. A disc diffusion assay showed that overexpression of TcGSTD3, but not TcGSTD1 or TcGSTD2, in E. coli increased resistance to paraquat-induced oxidative stress. RNA interference knockdown of TcGSTd1 caused metamorphosis deficiencies and reduced fecundity by regulating insulin/target-of-rapamycin signalling pathway-mediated ecdysteroid biosynthesis, and knockdown of TcGSTd3 led to reduced fertility and a decreased hatch rate of the offspring, probably caused by the reduced antioxidative activity in the reproductive organs. These results indicate that TcGSTd3 and TcGSTd2 may play vital roles in cellular detoxification, whereas TcGSTd1 may play essential roles in normal development of T. castaneum. These delta-class GSTs in T. castaneum have obtained different functions during the evolution.

A new neuropeptide insect parathyroid hormone iPTH in the red flour beetle Tribolium castaneum

In the postgenomics era, comparative genomics have advanced the understanding of evolutionary processes of neuropeptidergic signaling systems. The evolutionary origin of many neuropeptidergic signaling systems can be traced date back to early metazoan evolution based on the conserved sequences. Insect parathyroid hormone receptor (iPTHR) was previously described as an ortholog of vertebrate PTHR that has a well-known function in controlling bone remodeling. However, there was no sequence homologous to PTH sequence in insect genomes, leaving the iPTHR as an orphan receptor. Here, we identified the authentic ligand insect PTH (iPTH) for the iPTHR. The taxonomic distribution of iPTHR, which is lacking in Diptera and Lepidoptera, provided a lead for identifying the authentic ligand. We found that a previously described orphan ligand known as PXXXamide (where X is any amino acid) described in the cuttlefish Sepia officinalis has a similar taxonomic distribution pattern as iPTHR. Tests of this peptide, iPTH, in functional reporter assays confirmed the interaction of the ligand-receptor pair. Study of a model beetle, Tribolium castaneum, was used to investigate the function of the iPTH signaling system by RNA interference followed by RNA sequencing and phenotyping. The results suggested that the iPTH system is likely involved in the regulation of cuticle formation that culminates with a phenotype of defects in wing exoskeleton maturation at the time of adult eclosion. Moreover, RNAi of iPTHRs also led to significant reductions in egg numbers and hatching rates after parental RNAi.

Vertebrate parathyroid hormone (PTH) and its receptors have been extensively studied with respect to their function in bone remodeling and calcium metabolism. Insect parathyroid hormone receptors (iPTHRs) have been previously described as counterparts of vertebrate PTHRs, however, they are still orphan receptors for which the authentic ligands and biological functions remain unknown. We describe an insect form of parathyroid hormone (iPTH) by analyzing its interactions with iPTHRs. Identification of this new insect peptidergic system proved that the PTH system is an ancestral signaling system dating back to the evolutionary time before the divergence of protostomes and deuterostomes. We also investigated the functions of the iPTH system in a model beetle Tribolium castaneum by using RNA interference. RNA interference of iPTHR resulted in defects in wing exoskeleton maturation and fecundity. Based on the differential gene expression patterns and the phenotype induced by RNAi, we propose that the iPTH system is likely involved in the regulation of exoskeletal cuticle formation and fecundity in insects.

Latrophilin participates in insecticide susceptibility through positively regulating CSP10 and partially compensated by OBPC01 in Tribolium castaneum

Latrophilin (LPH) is an adhesion G protein-coupled receptor (aGPCR) that participates in multiple essential physiological processes. Our previous studies have shown that lph is not only indispensable for the development and reproduction of red flour beetles (Tribolium castaneum), but also for their resistance against dichlorvos or carbofuran insecticides. However, the regulatory mechanism of lph-mediated insecticide susceptibility remains unclear. Here, we revealed that knockdown of lph in beetles resulted in opposing changes in two chemoreception genes, chemosensory protein 10 (CSP10) and odorant-binding protein C01 (OBPC01), in which the expression of TcCSP10 was downregulated, whereas the expression of TcOBPC01 was upregulated. TcCSP10 and TcOBPC01 were expressed at the highest levels in early pupal and late larval stages, respectively. High levels of expression of both these genes were observed in the heads (without antennae) of adults. TcCSP10 and TcOBPC01 were significantly induced by dichlorvos or carbofuran between 12 and 72 h (hrs) after exposure, suggesting that they are likely associated with increasing the binding affinity of insecticides, leading to a decrease in sensitivity to the insecticides. Moreover, once these two genes were knocked down, the susceptibility of the beetles to dichlorvos or carbofuran was enhanced. Additionally, RNA interference (RNAi) targeting of lph followed by exposure to dichlorvos or carbofuran also caused the opposing expression levels of TcCSP10 and TcOBPC01 compared to the expression levels of wild-type larvae treated with insecticides alone. All these results indicate that lph is involved in insecticide susceptibility through positively regulating TcCSP10; and the susceptibility could also further partially compensated for through the negative regulation of TcOBPC01 when lph was knockdown in the red flour beetle. Our studies shed new light on the molecular regulatory mechanisms of lph related to insecticide susceptibility.

Latrophilin mediates insecticides susceptibility and fecundity through two carboxylesterases, esterase4 and esterase6, in Tribolium castaneum

Latrophilin (LPH) is known as an adhesion G-protein-coupled receptor which involved in multiple physiological processes in organisms. Previous studies showed that lph not only involved the susceptibility to anticholinesterase insecticides but also affected fecundity in Tribolium castaneum. However, its regulatory mechanisms in these biological processes are still not clear. Here, we identified two potential downstream carboxylesterase (cce) genes of Tclph, esterase4 and esterase6, and further characterized their interactions with Tclph. After treatment of T. castaneum larvae with carbofuran or dichlorvos insecticides, the transcript levels of Tcest4 and Tcest6 were significantly induced from 12 to 72 h. RNAi against Tcest4 or Tcest6 led to the higher mortality compared with the controls after the insecticides treatment, suggesting that these two genes play a vital role in detoxification of insecticides in T. castaneum. Furthermore, with insecticides exposure to Tclph knockdown beetles, the expression of Tcest4 was upregulated but Tcest6 was downregulated, indicating that beetles existed a compensatory response against the insecticides. Additionally, RNAi of Tcest6 resulted in 43% reductions in female egg laying and completely inhibited egg hatching, which showed the similar phenotype as that of Tclph knockdown. These results indicated that Tclph affected fecundity by positively regulating Tcest6 expression. Our findings will provide a new insight into the molecular mechanisms of Tclph involved in physiological functions in T. castaneum.

The crayfish model (Orconectes rusticus), epigenetics and drug addiction research

Fundamental signs of epigenetic effects are variations in the expression of genes or phenotypic traits among isogenic mates. Therefore, genetically identical animals are in high demand for epigenetic research. There are many genetically identical animals, including natural parthenogens and inbred laboratory lineages or clones. However, most parthenogenetic animal taxa are very small in combined epigenetic and drug addiction research. Orconectes rusticus has a unique phylogenetic position, with 2-3 years of life span, which undergoes metamorphosis that creates developmental stages with distinctly different morphologies, unique lifestyles, and broad behavioral traits, even among isogenic mates reared in the same environment offer novel inroads for epigenetics studies. Moreover, the establishment of crayfish as a novel system for drug addiction with evidence of an automated, operant self-administration and conditioned-reward, withdrawal, reinstatement of the conditioned drug-induced reward sets the stage to investigate epigenetic mechanisms of drug addiction. We discuss behavioral, pharmacological and molecular findings from laboratory studies that document a broad spectrum of molecular and, behavioral evidence including potential hypotheses that can be tested with the crayfish model for epigenetic study in drug addiction research.

Evolutionary transition between invertebrates and vertebrates via methylation reprogramming in embryogenesis

Major evolutionary transitions are enigmas, and the most notable enigma is between invertebrates and vertebrates, with numerous spectacular innovations. To search for the molecular connections involved, we asked whether global epigenetic changes may offer a clue by surveying the inheritance and reprogramming of parental DNA methylation across metazoans. We focused on gametes and early embryos, where the methylomes are known to evolve divergently between fish and mammals. Here, we find that methylome reprogramming during embryogenesis occurs neither in pre-bilaterians such as cnidarians nor in protostomes such as insects, but clearly presents in deuterostomes such as echinoderms and invertebrate chordates, and then becomes more evident in vertebrates. Functional association analysis suggests that DNA methylation reprogramming is associated with development, reproduction and adaptive immunity for vertebrates, but not for invertebrates. Interestingly, the single HOX cluster of invertebrates maintains unmethylated status in all stages examined. In contrast, the multiple HOX clusters show dramatic dynamics of DNA methylation during vertebrate embryogenesis. Notably, the methylation dynamics of HOX clusters are associated with their spatiotemporal expression in mammals. Our study reveals that DNA methylation reprogramming has evolved dramatically during animal evolution, especially after the evolutionary transitions from invertebrates to vertebrates, and then to mammals.

Zelda and the maternal-to-zygotic transition in cockroaches

In the endopterygote Drosophila melanogaster, Zelda is an activator of the zygotic genome during the maternal-to-zygotic transition (MZT). Zelda binds cis-regulatory elements (TAGteam heptamers), making chromatin accessible for gene transcription. Zelda has been studied in other endopterygotes: Apis mellifera and Tribolium castaneum, and the paraneopteran Rhodnius prolixus. We studied Zelda in the cockroach Blattella germanica, a hemimetabolan, short germ-band, and polyneopteran species. B. germanica Zelda has the complete set of functional domains, which is typical of species displaying ancestral features concerning embryogenesis. Interestingly, we found D. melanogaster TAGteam heptamers in the B. germanica genome. The canonical one, CAGGTAG, is present at a similar proportion in the genome of these two species and in the genome of other insects, suggesting that the genome admits as many CAGGTAG motifs as its length allows. Zelda-depleted embryos of B. germanica show defects involving blastoderm formation and abdomen development, and genes contributing to these processes are down-regulated. We conclude that in B. germanica, Zelda strictly activates the zygotic genome, within the MZT, a role conserved in more derived endopterygote insects. In B. germanica, zelda is expressed during MZT, whereas in D. melanogaster and T. castaneum it is expressed beyond this transition. In these species and A. mellifera, Zelda has functions even in postembryonic development. The expansion of zelda expression beyond the MZT in endopterygotes might be related with the evolutionary innovation of holometabolan metamorphosis. DATABASES: The RNA-seq datasets of B. germanica, D. melanogaster, and T. castaneum are accessible at the GEO databases GSE99785, GSE18068, GSE63770, and GSE84253. In addition, the RNA-seq library from T. castaneum adult females is available at SRA: SRX021963. The B. germanica reference genome is available as BioProject PRJNA203136.

Dnmt1 has an essential function despite the absence of CpG DNA methylation in the red flour beetle Tribolium castaneum

Epigenetic mechanisms, such as CpG DNA methylation enable phenotypic plasticity and rapid adaptation to changing environments. CpG DNA methylation is established by DNA methyltransferases (DNMTs), which are well conserved across vertebrates and invertebrates. There are insects with functional DNA methylation despite lacking a complete set of Dnmts. But at least one of the enzymes, DNMT1, appears to be required to maintain an active DNA methylation system. The red flour beetle, Tribolium castaneum, lacks Dnmt3 but possesses Dnmt1 and it has been controversial whether it has a functional DNA methylation system. Using whole genome bisulfite sequencing, we did not find any defined patterns of CpG DNA methylation in embryos. Nevertheless, we found Dnmt1 expressed throughout the entire life cycle of the beetle, with mRNA transcripts significantly more abundant in eggs and ovaries. A maternal knockdown of Dnmt1 caused a developmental arrest in offspring embryos. We show that Dnmt1 plays an essential role in T. castaneum embryos and that its downregulation leads to an early developmental arrest. This function appears to be unrelated to DNA methylation, since we did not find any evidence for this modification. This strongly suggests an alternative role of this protein.

Consistent inverse correlation between DNA methylation of the first intron and gene expression across tissues and species

BACKGROUND

DNA methylation is one of the main epigenetic mechanisms for the regulation of gene expression in eukaryotes. In the standard model, methylation in gene promoters has received the most attention since it is generally associated with transcriptional silencing. Nevertheless, recent studies in human tissues reveal that methylation of the region downstream of the transcription start site is highly informative of gene expression. Also, in some cell types and specific genes it has been found that methylation of the first intron, a gene feature typically rich in enhancers, is linked with gene expression. However, a genome-wide, tissue-independent, systematic comparative analysis of the relationship between DNA methylation in the first intron and gene expression across vertebrates has not been explored yet.

RESULTS

The most important findings of this study are: (1) using different tissues from a modern fish, we show a clear genome-wide, tissue-independent quasi-linear inverse relationship between DNA methylation of the first intron and gene expression. (2) This relationship is conserved across vertebrates, since it is also present in the genomes of a model pufferfish, a model frog and different human tissues. Among the gene features, tissues and species interrogated, the first intron's negative correlation with the gene expression was most consistent. (3) We identified more tissue-specific differentially methylated regions (tDMRs) in the first intron than in any other gene feature. These tDMRs have positive or negative correlation with gene expression, indicative of distinct mechanisms of tissue-specific regulation. (4) Lastly, we identified CpGs in transcription factor binding motifs, enriched in the first intron, the methylation of which tended to increase with the distance from the first exon-first intron boundary, with a concomitant decrease in gene expression.

CONCLUSIONS

Our integrative analysis clearly reveals the important and conserved role of the methylation level of the first intron and its inverse association with gene expression regardless of tissue and species. These findings not only contribute to our basic understanding of the epigenetic regulation of gene expression but also identify the first intron as an informative gene feature regarding the relationship between DNA methylation and gene expression where future studies should be focused.

Comparative Transcriptomics in Two Extreme Neopterans Reveals General Trends in the Evolution of Modern Insects

The success of neopteran insects, with 1 million species described, is associated with developmental innovations such as holometaboly and the evolution from short to long germband embryogenesis. To unveil the mechanisms underlining these innovations, we compared gene expression during the ontogeny of two extreme neopterans, the cockroach Blattella germanica (polyneopteran, hemimetabolan, and short germband species) and the fly Drosophila melanogaster (endopterygote, holometabolan, and long germband species). Results revealed that genes associated with metamorphosis are predominantly expressed in late nymphal stages in B. germanica and in the early-mid embryo in D. melanogaster. In B. germanica the maternal to zygotic transition (MZT) concentrates early in embryogenesis, when juvenile hormone factors are significantly expressed. In D. melanogaster, the MZT extends throughout embryogenesis, during which time juvenile hormone factors appear to be unimportant. These differences possibly reflect broad trends in the evolution of development within neopterans, related to the germband type and the metamorphosis mode.

•

Transcriptomes of cockroaches and flies show key differences along development

•

Cockroaches and flies express metamorphosis factors with distinct timings in ontogeny

•

Cockroaches methylate DNA in early embryogenesis, whereas flies do not

•

MZT is limited to the early embryo in cockroaches, but it extends until hatching in flies

Methylation profile of a satellite DNA constituting the intercalary G+C-rich heterochromatin of the cut trough shell Spisula subtruncata (Bivalvia, Mactridae)

Tandemly repeated DNAs usually constitute significant portions of eukaryotic genomes. In bivalves, however, repetitive DNAs are habitually not widespread. In our search for abundant repetitive DNAs in trough shells, we discovered a novel satellite DNA, SSUsat, which constitutes at least 1.3% of the genome of Spisula subtruncata. As foreseen by the satellite DNA library hypothesis, we confirmed that this satellite DNA is also present in two other Mactridae species, showing a highly conserved nucleotide sequence together with a dramatic diminution in the number of repeats. Predominantly located at the G + C-rich intercalary heterochromatin of S. subtruncata, SSUsat displays several DNA methylation peculiarities. The level of methylation of SSUsat is high (3.38%) in comparison with bivalve standards and triplicates the mean of the S. subtruncata genome (1.13%). Methylation affects not only the cytosines in CpG dinucleotides but also those in CHH and CHG trinucleotides, a feature common in plants but scarce and without any clear known relevance in animals. SSUsat segments enriched in methylated cytosines partly overlap those showing higher sequence conservation. The presence of a chromosome pair showing an accumulation of markedly under-methylated SSUsat monomers additionally indicates that the methylation processes that shape repetitive genome compartments are quite complex.

DNA methylation profiles in cancer diagnosis and therapeutics

Cancer initiation and proliferation is regulated by both epigenetic and genetic events with epigenetic modifications being increasingly identified as important targets for cancer research. DNA methylation catalyzed by DNA methyltransferases (DNMTs) is one of the essential epigenetic mechanisms that control cell proliferation, apoptosis, differentiation, cell cycle, and transformation in eukaryotes. Recent progress in epigenetics revealed a deeper understanding of the mechanisms of tumorigenesis and provided biomarkers for early detection, diagnosis, and prognosis in cancer patients. Although DNA methylation biomarker possesses potential contributing to precision medicine, there are still limitations to be overcome before it reaches clinical setting. Hence, the current status of DNA methylation biomarkers was reviewed and the future use in clinic was also predicted.

5-Aza-2'-deoxycytidine in the medial prefrontal cortex regulates alcohol-related behavior and Ntf3-TrkC expression in rats

Recent studies have indicated that DNA methylation plays an important role in the development of alcohol abuse. 5-Aza-2'-deoxycytidine (5-Aza-dc), an inhibitor of DNA methyltransferases, was FDA approved for myelodysplastic syndrome treatment. However, it is unclear whether 5-Aza-dc is involved in alcohol abuse. In this study, using a chronic alcohol exposure model in rats, 5-Aza-dc was injected into the medial prefrontal cortex (mPFC). Alcohol-drinking behavior and the anxiety related behavior were evaluated by two-bottle choice and open field test. We found that 5-Aza-dc injection into the mPFC significantly decreased alcohol consumption and alcohol preference in alcohol-exposure rats, corresponding to the reduced blood alcohol levels. Although 5-Aza-dc potentiated the anxiety-like behavior of alcohol-exposure rats, it had no effect on the locomotor activity. Moreover, both of the mRNA and protein levels of DNA Methyltransferase 3A (DNMT3A) and DNMT3B in the mPFC were upregulated after 35 days of alcohol exposure and this upregulation could be reversed by 5-Aza-dc treatment. Additionally, 5-Aza-dc reversed the alcohol-induced downregulation of neurotrophin-3 (Ntf3), correspondingly the expression of its receptor-TrkC was reduced. These findings identified a functional role of 5-Aza-dc in alcohol-related behavioral phenotypes and one of the potential target genes, Ntf3. We also provide novel evidence for DNA methyltransferases as potential therapeutic targets in alcohol abuse.