Genomics of CpG methylation in developing and developed zebrafish

Remodeling of the epigenetic landscape in rainbow trout, Oncorhynchus mykiss, offspring in response to maternal choline intake

This project focused on evaluating the effects of maternal dietary choline intake on global DNA methylation profiles and related transcriptional changes in rainbow trout offspring. Three experimental diets were formulated to test different levels of choline intake: (a) 2065 ppm choline (Low Choline, 0 % supplementation), (b) 5657 ppm choline (Medium Choline, 0.6 % supplementation), and (c) 9248 ppm choline (High Choline, 1.2 % choline supplementation). Six rainbow trout families were fed experimental diets beginning 18 months post-hatch until spawning; their offspring were fed a commercial diet. Reduced representation bisulfite sequencing (RRBS) was utilized to measure genome-wide methylation in offspring immediately after hatching. When comparing to the Medium Choline offspring, differential DNA methylation occurred more in the Low Choline offspring than High Choline, especially in genic features like promoters. The differentially methylated CpGs (q ≤ 0.01) were identified evenly between CpG islands and shores in the genome, mostly found in the introns of genes. Genes such as fabp2 and leap2B associated with protein binding, fatty acid binding, DNA binding, and response to bacteria were differentially methylated and detected as differentially regulated genes by previous RNA-seq analysis. Although these findings indicate that levels of dietary choline available in broodstock diets alter offspring DNA methylation;, most differentially expressed genes were not associated with differential DNA methylation, suggesting additional mechanisms playing a role in regulating gene expression in response to maternal choline intake.

DNA methylation differences between stick insect ecotypes

Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high-resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome-wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome-wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host-plant effects. Using methylome-environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies.

MethylRAD Sequencing Technology Reveals DNA Methylation Characteristics of Apostichopus japonicus of Different Ages

The A. japonicus industry has expanded significantly, but no research has focused on determining the age of A. japonicus during farming. Correctly estimating the age of A. japonicus can provide a decision-making basis for the breeding process and data for the protection of A. japonicus aquatic germplasm resources. DNA methylation levels in the body wall of Apostichopus japonicus at 4 months, 1 year, 2 years, and 3 years old were determined using MethylRAD-Seq, and differentially methylated genes were screened. A total of 441 and 966 differentially methylated genes were detected at the CCGG and CCWGG sites, respectively. Aspartate aminotransferase, succinate semialdehyde dehydrogenase, isocitrate dehydrogenase, the histone H2AX, heat shock protein Hsp90, aminopeptidase N, cell division cycle CDC6, Ras GTPase activating protein (RasGAP), slit guidance ligand slit1, integrin-linked kinase ILK, mechanistic target of rapamycin kinase Mtor, protein kinase A Pka, and autophagy-related 3 atg3 genes may play key roles in the growth and aging process of A. japonicus. This study provides valuable information regarding age-related genes for future research, and these candidate genes can be used to create an "epigenetic clock".

Accounting for transgenerational effects of toxicant exposure in population models alters the predicted long-term population status

Acute environmental stressors such as short-term exposure to pollutants can have lasting effects on organisms, potentially impacting future generations. Parental exposure to toxicants can result in changes to the epigenome (e.g., DNA methylation) that are passed down to subsequent, unexposed generations. However, it is difficult to gauge the cumulative population-scale impacts of epigenetic effects from laboratory experiments alone. Here, we developed a size- and age-structured delay-coordinate population model to evaluate the long-term consequences of epigenetic modifications on population sustainability. The model emulated changes in growth, mortality, and fecundity in the F0, F1, and F2 generations observed in experiments in which larval Menidia beryllina were exposed to environmentally relevant concentrations of bifenthrin (Bif), ethinylestradiol (EE2), levonorgestrel (LV), or trenbolone (TB) in the parent generation (F0) and reared in clean water up to the F2 generation. Our analysis suggests potentially dramatic population-level effects of repeated, chronic exposures of early-life stage fish that are not captured by models not accounting for those effects. Simulated exposures led to substantial declines in population abundance (LV and Bif) or near-extinction (EE2 and TB) with the exact trajectory and timeline of population decline dependent on the combination of F0, F1, and F2 effects produced by each compound. Even acute one-time exposures of each compound led to declines and recovery over multiple years due to lagged epigenetic effects. These results demonstrate the potential for environmentally relevant concentrations of commonly used compounds to impact the population dynamics and sustainability of an ecologically relevant species and model organism.

Epigenetics underpins phenotypic plasticity of protandrous sex change in fish

Phenotypic plasticity is an important driver of species resilience. Often mediated by epigenetic changes, phenotypic plasticity enables individual genotypes to express variable phenotypes in response to environmental change. Barramundi (Lates calcarifer) are a protandrous (male-first) sequential hermaphrodite that exhibits plasticity in length-at-sex change between geographic regions. This plasticity is likely to be mediated by changes in DNA methylation (DNAm), a well-studied epigenetic modification. To investigate the relationships between length, sex, and DNAm in a sequential hermaphrodite, here, we compare DNAm in four conserved vertebrate sex-determining genes in male and female barramundi of differing lengths from three geographic regions of northern Australia. Barramundi first mature as male and later sex change to female upon the attainment of a larger body size; however, a general pattern of increasing female-specific DNAm markers with increasing length was not observed. Significant differences in DNAm between males and females of similar lengths suggest that female-specific DNAm arises rapidly during sex change, rather than gradually with fish growth. The findings also reveal that region-specific differences in length-at-sex change are accompanied by differences in DNAm and are consistent with variability in remotely sensed sea temperature and salinity. Together, these findings provide the first in situ evidence for epigenetically and environmentally mediated sex change in a protandrous hermaphrodite and offer significant insight into the molecular and ecological processes governing the marked and unique plasticity of sex in fish.

Developmental remodelling of non-CG methylation at satellite DNA repeats

In vertebrates, DNA methylation predominantly occurs at CG dinucleotides however, widespread non-CG methylation (mCH) has been reported in mammalian embryonic stem cells and in the brain. In mammals, mCH is found at CAC trinucleotides in the nervous system, where it is associated with transcriptional repression, and at CAG trinucleotides in embryonic stem cells, where it positively correlates with transcription. Moreover, CAC methylation appears to be a conserved feature of adult vertebrate brains. Unlike any of those methylation signatures, here we describe a novel form of mCH that occurs in the TGCT context within zebrafish mosaic satellite repeats. TGCT methylation is inherited from both male and female gametes, remodelled during mid-blastula transition, and re-established during gastrulation in all embryonic layers. Moreover, we identify DNA methyltransferase 3ba (Dnmt3ba) as the primary enzyme responsible for the deposition of this mCH mark. Finally, we observe that TGCT-methylated repeats are specifically associated with H3K9me3-marked heterochromatin suggestive of a functional interplay between these two gene-regulatory marks. Altogether, this work provides insight into a novel form of vertebrate mCH and highlights the substrate diversity of vertebrate DNA methyltransferases.

DNA Methylation Profiles Suggest Intergenerational Transfer of Maternal Effects

The view of maternal effects (nongenetic maternal environmental influence on offspring phenotype) has changed from one of distracting complications in evolutionary genetics to an important evolutionary mechanism for improving offspring fitness. Recent studies have shown that maternal effects act as an adaptive mechanism to prepare offspring for stressful environments. Although research into the magnitude of maternal effects is abundant, the molecular mechanisms of maternal influences on offspring phenotypic variation are not fully understood. Despite recent work identifying DNA methylation as a potential mechanism of nongenetic inheritance, currently proposed links between DNA methylation and parental effects are indirect and primarily involve genomic imprinting. We combined a factorial breeding design and gene-targeted sequencing methods to assess inheritance of methylation during early life stages at 14 genes involved in growth, development, metabolism, stress response, and immune function of Chinook salmon (Oncorhynchus tshawytscha). We found little evidence for additive or nonadditive genetic effects acting on methylation levels during early development; however, we detected significant maternal effects. Consistent with conventional maternal effect data, maternal effects on methylation declined through development and were replaced with nonadditive effects when offspring began exogenous feeding. We mapped methylation at individual CpG sites across the selected candidate genes to test for variation in site-specific methylation profiles and found significant maternal effects at selected CpG sites that also declined with development stage. While intergenerational inheritance of methylated DNA is controversial, we show that CpG-specific methylation may function as an underlying molecular mechanism for maternal effects, with important implications for offspring fitness.

Understanding 'Non-genetic' Inheritance: Insights from Molecular-Evolutionary Crosstalk

Understanding the evolutionary and ecological roles of 'non-genetic' inheritance (NGI) is daunting due to the complexity and diversity of epigenetic mechanisms. We draw on insights from molecular and evolutionary biology perspectives to identify three general features of 'non-genetic' inheritance systems: (i) they are functionally interdependent with, rather than separate from, DNA sequence; (ii) precise mechanisms vary phylogenetically and operationally; and (iii) epigenetic elements are probabilistic, interactive regulatory factors and not deterministic 'epialleles' with defined genomic locations and effects. We discuss each of these features and offer recommendations for future empirical and theoretical research that implements a unifying inherited gene regulation (IGR) approach to studies of 'non-genetic' inheritance.

Comparative epigenetics in animal physiology: An emerging frontier

The unprecedented access to annotated genomes now facilitates the investigation of the molecular basis of epigenetic phenomena in phenotypically diverse animals. In this critical review, we describe the roles of molecular epigenetic mechanisms in regulating mitotically and meiotically stable spatiotemporal gene expression, phenomena that provide the molecular foundation for the intra-, inter-, and trans-generational emergence of physiological phenotypes. By focusing principally on emerging comparative epigenetic roles of DNA-level and transcriptome-level epigenetic mark dynamics in the emergence of phenotypes, we highlight the relationship between evolutionary conservation and innovation of specific epigenetic pathways, and their interplay as a priority for future study. This comparative approach is expected to significantly advance our understanding of epigenetic phenomena, as animals show a diverse array of strategies to epigenetically modify physiological responses. Additionally, we review recent technological advances in the field of molecular epigenetics (single-cell epigenomics and transcriptomics and editing of epigenetic marks) in order to (1) investigate environmental and endogenous factor dependent epigenetic mark dynamics in an integrative manner; (2) functionally test the contribution of specific epigenetic marks for animal phenotypes via genome and transcript-editing tools. Finally, we describe advantages and limitations of emerging animal models, which under the Krogh principle, may be particularly useful in the advancement of comparative epigenomics and its potential translational applications in animal science, ecotoxicology, ecophysiology, climate change science and wild-life conservation, as well as organismal health.

The round goby genome provides insights into mechanisms that may facilitate biological invasions

Background

The invasive benthic round goby (Neogobius melanostomus) is the most successful temperate invasive fish and has spread in aquatic ecosystems on both sides of the Atlantic. Invasive species constitute powerful in situ experimental systems to study fast adaptation and directional selection on short ecological timescales and present promising case studies to understand factors involved the impressive ability of some species to colonize novel environments. We seize the unique opportunity presented by the round goby invasion to study genomic substrates potentially involved in colonization success.

Results

We report a highly contiguous long-read-based genome and analyze gene families that we hypothesize to relate to the ability of these fish to deal with novel environments. The analyses provide novel insights from the large evolutionary scale to the small species-specific scale. We describe expansions in specific cytochrome P450 enzymes, a remarkably diverse innate immune system, an ancient duplication in red light vision accompanied by red skin fluorescence, evolutionary patterns of epigenetic regulators, and the presence of osmoregulatory genes that may have contributed to the round goby’s capacity to invade cold and salty waters. A recurring theme across all analyzed gene families is gene expansions.

Conclusions

The expanded innate immune system of round goby may potentially contribute to its ability to colonize novel areas. Since other gene families also feature copy number expansions in the round goby, and since other Gobiidae also feature fascinating environmental adaptations and are excellent colonizers, further long-read genome approaches across the goby family may reveal whether gene copy number expansions are more generally related to the ability to conquer new habitats in Gobiidae or in fish.

Electronic supplementary material

The online version of this article (10.1186/s12915-019-0731-8) contains supplementary material, which is available to authorized users.

PM2.5-induced extensive DNA methylation changes in the heart of zebrafish embryos and the protective effect of folic acid

We previously found that folic acid (FA) attenuated cardiac defects in zebrafish embryos exposed to extractable organic matter (EOM) from PM2.5, but the underlining mechanisms remain to be elucidated. Since DNA methylation is crucial to cardiac development, we hypothesized that EOM-induced aberrant DNA methylation changes could be diminished by FA supplementation. In this study, zebrafish embryos were exposed to EOM in the absence or presence of FA. Genomic-wide DNA methylation analysis identified both DNA hypo- and hyper-methylation changes in CCGG sites in zebrafish embryos exposed to EOM, which were attenuated by FA supplementation. We identified a total of 316 genes with extensive DNA methylation changes in EOM samples but little or no DNA methylation changes in EOM plus FA samples. The genes were involved in critical cellular processes and signaling pathways important for embryo development. In addition, the EOM-decreased SAM/SAH ratio was counteracted by FA supplementation. Furthermore, FA attenuated the EOM-induced changes in the expression of genes involved in the regulation of DNA methylation and in folate biosynthesis. In conclusion, our data suggest that FA supplementation protected zebrafish embryos from the cardiac developmental toxicity of PM2.5 by alleviating EOM-induced DNA methylation changes.

DNA methylation is maintained with high fidelity in the honey bee germline and exhibits global non-functional fluctuations during somatic development

BACKGROUND

DNA methylation of active genes, also known as gene body methylation, is found in many animal and plant genomes. Despite this, the transcriptional and developmental role of such methylation remains poorly understood. Here, we explore the dynamic range of DNA methylation in honey bee, a model organism for gene body methylation.

RESULTS

Our data show that CG methylation in gene bodies globally fluctuates during honey bee development. However, these changes cause no gene expression alterations. Intriguingly, despite the global alterations, tissue-specific CG methylation patterns of complete genes or exons are rare, implying robust maintenance of genic methylation during development. Additionally, we show that CG methylation maintenance fluctuates in somatic cells, while reaching maximum fidelity in sperm cells. Finally, unlike universally present CG methylation, we discovered non-CG methylation specifically in bee heads that resembles such methylation in mammalian brain tissue.

CONCLUSIONS

Based on these results, we propose that gene body CG methylation can oscillate during development if it is kept to a level adequate to preserve function. Additionally, our data suggest that heightened non-CG methylation is a conserved regulator of animal nervous systems.

Profiling DNA methylation patterns of zebrafish liver associated with parental high dietary arachidonic acid

Diet has been shown to influence epigenetic key players, such as DNA methylation, which can regulate the gene expression potential in both parents and offspring. Diets enriched in omega-6 and deficient in omega-3 PUFAs (low dietary omega-3/omega-6 PUFA ratio), have been associated with the promotion of pathogenesis of diseases in humans and other mammals. In this study, we investigated the impact of increased dietary intake of arachidonic acid (ARA), a physiologically important omega-6 PUFA, on 2 generations of zebrafish. Parental fish were fed either a low or a high ARA diet, while the progeny of both groups were fed the low ARA diet. We screened for DNA methylation on single base-pair resolution using reduced representation bisulfite sequencing (RRBS). The DNA methylation profiling revealed significant differences between the dietary groups in both parents and offspring. The majority of differentially methylated loci associated with high dietary ARA were found in introns and intergenic regions for both generations. Common loci between the identified differentially methylated loci in F0 and F1 livers were reported. We described overlapping gene annotations of identified methylation changes with differential expression, but based on a small number of overlaps. The present study describes the diet-associated methylation profiles across genomic regions, and it demonstrates that parental high dietary ARA modulates DNA methylation patterns in zebrafish liver.

Contrasting DNA methylation responses of inbred fish lines to different rearing environments

Epigenetic mechanisms generate plastic phenotypes that can become locally adapted across environments. Disentangling genomic from epigenomic variation is challenging in sexual species due to genetic variation among individuals, but it is easier in self-fertilizing species. We analysed DNA methylation patterns of two highly inbred strains of a naturally self-fertilizing fish reared in two contrasting environments to investigate the obligatory (genotype-dependent), facilitated (partially depend on the genotype) or pure (genotype-independent) nature of the epigenetic variation. We found higher methylation differentiation between genotypes than between environments. Most methylation differences between environments common to both strains followed a pattern where the two genotypes (inbred lines) responded to the same environmental context with contrasting DNA methylation levels (facilitated epialleles). Our findings suggest that, at least in part, DNA methylation could depend on the dynamic interaction between the genotype and the environment, which could explain the plasticity of epigenetically mediated phenotypes.

DNA Methylation Patterns in the Round Goby Hypothalamus Support an On-The-Spot Decision Scenario for Territorial Behavior

The question as to how early life experiences are stored on a molecular level and affect traits later in life is highly topical in ecology, medicine, and epigenetics. In this study, we use a fish model to investigate whether DNA methylation mediates early life experiences and predetermines a territorial male reproductive phenotype. In fish, adult reproductive phenotypes frequently depend on previous life experiences and are often associated with distinct morphological traits. DNA methylation is an epigenetic mechanism which is both sensitive to environmental conditions and stably inherited across cell divisions. We therefore investigate early life predisposition in the round goby Neogobius melanostomus by growth back-calculations and then study DNA methylation by MBD-Seq in the brain region controlling vertebrate reproductive behavior, the hypothalamus. We find a link between the territorial reproductive phenotype and high growth rates in the first year of life. However, hypothalamic DNA methylation patterns reflect the current behavioral status independently of early life experiences. Together, our data suggest a non-predetermination scenario in the round goby, in which indeterminate males progress to a non-territorial status in the spawning season, and in which some males then assume a specialized territorial phenotype if current conditions are favorable.

Developmental Dioxin Exposure Alters the Methylome of Adult Male Zebrafish Gonads

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is a persistent environmental toxicant and endocrine disrupting compound with reproductive and developmental effects in humans and model organisms, including zebrafish. Our previous microarray and histological studies found defects in spermatogenesis and fertility of zebrafish in response to acute developmental TCDD exposure. These effects are apparent following exposure during reproductive development, modeling fetal basis of adult-onset disease. Some outcomes of these previous studies (reduced fertility, changes in sex ratio, transcriptomic alterations) are also transgenerational – persisting to unexposed generations – through the male germline. We hypothesized that DNA methylation could be a possible mechanism for these reproductive effects and performed whole genome bisulfite sequencing (WGBS), which identifies whole genome DNA methylation status at the base pair level, on testes of adult zebrafish exposed to TCDD (two separate hour-long exposures to 50 pg/mL TCDD at 3 and 7 weeks post fertilization). In response to TCDD exposure, multiple genes were differentially methylated; many of which are involved in reproductive processes or epigenetic modifications, suggesting a role of DNA methylation in later-life health outcomes. Additionally, several differentially methylated genes corresponded with gene expression changes identified in TCDD-exposed zebrafish testes, indicating a potential link between DNA methylation and gene expression. Ingenuity pathway analysis of WGBS and microarray data revealed genes involved in reproductive processes and development, RNA regulation, the cell cycle, and cellular morphology and development. We conclude that site-specific changes in DNA methylation of adult zebrafish testes occur in response to acute developmental TCDD exposure.

DNA methylation and gene expression alterations in zebrafish early-life stages exposed to the antibacterial agent triclosan

There is increasing evidence that toxicant exposure can alter DNA methylation profile, one of the main epigenetic mechanisms, particularly during embryogenesis when DNA methylation patterns are being established. In order to investigate the effects of the antibacterial agent Triclosan on DNA methylation and its correlation with gene expression, zebrafish embryos were exposed during 7 days post-fertilization (starting at maximum 8-cells stage) to 50 and 100 μg/l, two conditions for which increased sensitivity and acclimation have been respectively reported. Although global DNA methylation was not significantly affected, a total of 171 differentially methylated fragments were identified by Reduced Representation Bisulfite Sequencing. The majority of these fragments were found between the two exposed groups, reflecting dose-dependant specific responses. Gene ontology analysis revealed that pathways involved in TGF-β signaling were enriched in larvae exposed to 50 μg/l, while de novo pyrimidine biosynthesis functions were overrepresented in fish exposed to 100 μg/l. In addition, gene expression analysis revealed a positive correlation between mRNA levels and DNA methylation patterns in introns, together with significant alterations of the transcription of genes involved in nervous system development, transcriptional factors and histone methyltransferases. Overall this work provides evidence that Triclosan alters DNA methylation in zebrafish exposed during embryogenesis as well as related genes expression and proposes concentration specific modes of action. Further studies will investigate the possible long-term consequences of these alterations, i.e. latent defects associated with developmental exposure and transgenerational effects, and the possible implications in terms of fitness and adaptation to environmental pollutants.

Chromatin accessibility is associated with CRISPR-Cas9 efficiency in the zebrafish (Danio rerio)

CRISPR-Cas9 technology is routinely applied for targeted mutagenesis in model organisms and cell lines. Recent studies indicate that the prokaryotic CRISPR-Cas9 system is affected by eukaryotic chromatin structures. Here, we show that the likelihood of successful mutagenesis correlates with transcript levels during early development in zebrafish (Danio rerio) embryos. In an experimental setting, we found that guide RNAs differ in their onset of mutagenesis activity in vivo. Furthermore, some guide RNAs with high in vitro activity possessed poor mutagenesis activity in vivo, suggesting the presence of factors that limit the mutagenesis in vivo. Using open access datasets generated from early developmental stages of the zebrafish, and guide RNAs selected from the CRISPRz database, we provide further evidence for an association between gene expression during early development and the success of CRISPR-Cas9 mutagenesis in zebrafish embryos. In order to further inspect the effect of chromatin on CRISPR-Cas9 mutagenesis, we analysed the relationship of selected chromatin features on CRISPR-Cas9 mutagenesis efficiency using publicly available data from zebrafish embryos. We found a correlation between chromatin openness and the efficiency of CRISPR-Cas9 mutagenesis. These results indicate that CRISPR-Cas9 mutagenesis is influenced by chromatin accessibility in zebrafish embryos.

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination

Background

DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.

Results

Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.

Conclusions

Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0136-2) contains supplementary material, which is available to authorized users.

The elephant shark methylome reveals conservation of epigenetic regulation across jawed vertebrates

BACKGROUND

Methylation of CG dinucleotides constitutes a critical system of epigenetic memory in bony vertebrates, where it modulates gene expression and suppresses transposon activity. The genomes of studied vertebrates are pervasively hypermethylated, with the exception of regulatory elements such as transcription start sites (TSSs), where the presence of methylation is associated with gene silencing. This system is not found in the sparsely methylated genomes of invertebrates, and establishing how it arose during early vertebrate evolution is impeded by a paucity of epigenetic data from basal vertebrates.

METHODS

 We perform whole-genome bisulfite sequencing to generate the first genome-wide methylation profiles of a cartilaginous fish, the elephant shark Callorhinchus milii. Employing these to determine the elephant shark methylome structure and its relationship with expression, we compare this with higher vertebrates and an invertebrate chordate using published methylation and transcriptome data.  Results: Like higher vertebrates, the majority of elephant shark CG sites are highly methylated, and methylation is abundant across the genome rather than patterned in the mosaic configuration of invertebrates. This global hypermethylation includes transposable elements and the bodies of genes at all expression levels. Significantly, we document an inverse relationship between TSS methylation and expression in the elephant shark, supporting the presence of the repressive regulatory architecture shared by higher vertebrates.

CONCLUSIONS

 Our demonstration that methylation patterns in a cartilaginous fish are characteristic of higher vertebrates imply the conservation of this epigenetic modification system across jawed vertebrates separated by 465 million years of evolution. In addition, these findings position the elephant shark as a valuable model to explore the evolutionary history and function of vertebrate methylation.

Much more than you expected: The non-DHFR-mediated effects of methotrexate

BACKGROUND

For decades, methotrexate (MTX; amethopterin) has been known as an antifolate inhibitor of dihydrofolate reductase (DHFR), and it is widely used for the treatment of various malignancies and autoimmune diseases. Although the inclusion of MTX in various therapeutic regimens is based on its ability to inhibit DHFR and consequently to suppress the synthesis of pyrimidine and purine precursors, recent studies have shown that MTX is also able to target other intracellular pathways that are independent of folate metabolism.

SCOPE OF REVIEW

The main aim of this review is to summarize the most important, up-to-date findings of studies regarding the non-DHFR-mediated mechanisms of MTX action.

MAJOR CONCLUSIONS

The effectiveness of MTX is undoubtedly caused by its capability to affect various intracellular pathways at many levels. Although the most important therapeutic mechanism of MTX is strongly based on the inhibition of DHFR, many other effects of this compound have been described and new studies bring new insights into the pharmacology of MTX every year.

GENERAL SIGNIFICANCE

Identification of these new targets for MTX is especially important for a better understanding of MTX action in new protocols of combination therapy.

Gene copy silencing and DNA methylation in natural and artificially produced allopolyploid fish

Allelic silencing is an important mechanism for coping with gene dosage changes in polyploid organisms that is well known in allopolyploid plants. Only recently, it was shown in the allotriploid fish Squalius alburnoides that this process also occurs in vertebrates. However, it is still unknown whether this silencing mechanism is common to other allopolyploid fish, and which mechanisms might be responsible for allelic silencing. We addressed these questions in a comparative study between Squalius alburnoides and another allopolyploid complex, the Amazon molly (Poecilia formosa). We examined the allelic expression patterns for three target genes in four somatic tissues of natural allo-anorthoploids and laboratory-produced tri-genomic hybrids of S. alburnoides and P. formosa. Also, for both complexes, we evaluated the correlation between total DNA methylation level and the ploidy status and genomic composition of the individuals. We found that allelic silencing also occurs in other allopolyploid organisms besides the single one that was previously known. We found and discuss disparities within and between the two considered complexes concerning the pattern of allele-specific expression and DNA methylation levels. Disparities might be due to intrinsic characteristics of each genome involved in the hybridization process. Our findings also support the idea that long-term evolutionary processes have an effect on the allele expression patterns and possibly also on DNA methylation levels.

Evolutionary Consequences of DNA Methylation in a Basal Metazoan

Gene body methylation (gbM) is an ancestral and widespread feature in Eukarya, yet its adaptive value and evolutionary implications remain unresolved. The occurrence of gbM within protein-coding sequences is particularly puzzling, because methylation causes cytosine hypermutability and hence is likely to produce deleterious amino acid substitutions. We investigate this enigma using an evolutionarily basal group of Metazoa, the stony corals (order Scleractinia, class Anthozoa, phylum Cnidaria). We show that patterns of coral gbM are similar to other invertebrate species, predicting wide and active transcription and slower sequence evolution. We also find a strong correlation between gbM and codon bias, resulting from systematic replacement of CpG bearing codons. We conclude that gbM has strong effects on codon evolution and speculate that this may influence establishment of optimal codons.

Genome-Wide Analysis of Transposon and Retroviral Insertions Reveals Preferential Integrations in Regions of DNA Flexibility

DNA transposons and retroviruses are important transgenic tools for genome engineering. An important consideration affecting the choice of transgenic vector is their insertion site preferences. Previous large-scale analyses of Ds transposon integration sites in plants were done on the basis of reporter gene expression or germ-line transmission, making it difficult to discern vertebrate integration preferences. Here, we compare over 1300 Ds transposon integration sites in zebrafish with Tol2 transposon and retroviral integration sites. Genome-wide analysis shows that Ds integration sites in the presence or absence of marker selection are remarkably similar and distributed throughout the genome. No strict motif was found, but a preference for structural features in the target DNA associated with DNA flexibility (Twist, Tilt, Rise, Roll, Shift, and Slide) was observed. Remarkably, this feature is also found in transposon and retroviral integrations in maize and mouse cells. Our findings show that structural features influence the integration of heterologous DNA in genomes, and have implications for targeted genome engineering.

Zebrafish Discoveries in Cancer Epigenetics

The cancer epigenome is fundamentally different than that of normal cells. How these differences arise in and contribute to carcinogenesis is not known, and studies using model organisms such as zebrafish provide an opportunity to address these important questions. Modifications of histones and DNA comprise the complex epigenome, and these influence chromatin structure, genome stability and gene expression, all of which are fundamental to the cellular changes that cause cancer. The cancer genome atlas covers the wide spectrum of genetic changes associated with nearly every cancer type, however, this catalog is currently uni-dimensional. As the pattern of epigenetic marks and chromatin structure in cancer cells is described and overlaid on the mutational landscape, the map of the cancer genome becomes multi-dimensional and highly complex. Two major questions remain in the field: (1) how the epigenome becomes repatterned in cancer and (2) which of these changes are cancer-causing. Zebrafish provide a tractable in vivo system to monitor the epigenome during transformation and to identify epigenetic drivers of cancer. In this chapter, we review principles of cancer epigenetics and discuss recent work using zebrafish whereby epigenetic modifiers were established as cancer driver genes, thus providing novel insights into the mechanisms of epigenetic reprogramming in cancer.

Folic acid protects against arsenic-mediated embryo toxicity by up-regulating the expression of Dvr1

As a nutritional factor, folic acid can prevent cardiac and neural defects during embryo development. Our previous study showed that arsenic impairs embryo development by down-regulating Dvr1/GDF1 expression in zebrafish. Here, we investigated whether folic acid could protect against arsenic-mediated embryo toxicity. We found that folic acid supplementation increases hatching and survival rates, decreases malformation rate and ameliorates abnormal cardiac and neural development of zebrafish embryos exposed to arsenite. Both real-time PCR analysis and whole in-mount hybridization showed that folic acid significantly rescued the decrease in Dvr1 expression caused by arsenite. Subsequently, our data demonstrated that arsenite significantly decreased cell viability and GDF1 mRNA and protein levels in HEK293ET cells, while folic acid reversed these effects. Folic acid attenuated the increase in subcellular reactive oxygen species (ROS) levels and oxidative adaptor p66Shc protein expression in parallel with the changes in GDF1 expression and cell viability. P66Shc knockdown significantly inhibited the production of ROS and the down-regulation of GDF1 induced by arsenite. Our data demonstrated that folic acid supplementation protected against arsenic-mediated embryo toxicity by up-regulating the expression of Dvr1/GDF1, and folic acid enhanced the expression of GDF1 by decreasing p66Shc expression and subcellular ROS levels.

DNA methylation of genes of the main components of the telomerase complex in Danio rerio

The methylation status of the genes of telomerase reverse transcriptase (tert) and telomerase RNA (terc) was determined in brain tissues of Danio rerio of different age. It is found that, regardless of the age of fish, the regulatory region of the tert gene was completely methylated, whereas the coding region remained unmethylated in all cases. The level of methylation of the region located downstream of the coding region of the terc gene changes with age. This region was analyzed in the samples of other tissues, and its methylation status was also nonuniform. The alteration of the methylation status in the 3'-untranslated region of the terc gene suggests the possibility of transcription of the antisense strand in this region.

Developmental exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin alters DNA methyltransferase (dnmt) expression in zebrafish (Danio rerio)

DNA methylation is one of the most important epigenetic modifications involved in the regulation of gene expression. The DNA methylation reaction is catalyzed by DNA methyltransferases (DNMTs). Recent studies have demonstrated that toxicants can affect normal development by altering DNA methylation patterns, but the mechanisms of action are poorly understood. Hence, we tested the hypothesis that developmental exposure to TCDD affects dnmt gene expression patterns. Zebrafish embryos were exposed to 5nM TCDD for 1h from 4 to 5h post-fertilization (hpf) and sampled at 12, 24, 48, 72, and 96 hpf to determine dnmt gene expression and DNA methylation patterns. We performed a detailed analysis of zebrafish dnmt gene expression during development and in adult tissues. Our results demonstrate that dnmt3b genes are highly expressed in early stages of development, and dnmt3a genes are more abundant in later stages. TCDD exposure upregulated dnmt1 and dnmt3b2 expression, whereas dnmt3a1, 3b1, and 3b4 are downregulated following exposure. We did not observe any TCDD-induced differences in global methylation or hydroxymethylation levels, but the promoter methylation of aryl hydrocarbon receptor (AHR) target genes was altered. In TCDD-exposed embryos, AHR repressor a (ahrra) and c-fos promoters were differentially methylated. To characterize the TCDD effects on DNMTs, we cloned the dnmt promoters with xenobiotic response elements and conducted AHR transactivation assays using a luciferase reporter system. Our results suggest that ahr2 can regulate dnmt3a1, dnmt3a2, and dnmt3b2 expression. Overall, we demonstrate that developmental exposure to TCDD alters dnmt expression and DNA methylation patterns.