Single-cell m6A mapping in vivo using picoMeRIP-seq.

Current N6-methyladenosine (m6A) mapping methods need large amounts of RNA or are limited to cultured cells. Through optimized sample recovery and signal-to-noise ratio, we developed picogram-scale m6A RNA immunoprecipitation and sequencing (picoMeRIP-seq) for studying m6A in vivo in single cells and scarce cell types using standard laboratory equipment. We benchmark m6A mapping on titrations of poly(A) RNA and embryonic stem cells and in single zebrafish zygotes, mouse oocytes and embryos.

Altered non-coding RNA expression profile in F1 progeny 1 year after parental irradiation is linked to adverse effects in zebrafish.

Gamma radiation produces DNA instability and impaired phenotype. Previously, we observed negative effects on phenotype, DNA methylation, and gene expression profiles, in offspring of zebrafish exposed to gamma radiation during gametogenesis. We hypothesize that previously observed effects are accompanied with changes in the expression profile of non-coding RNAs, inherited by next generations. Non-coding RNA expression profile was analysed in F1 offspring (5.5 h post-fertilization) by high-throughput sequencing 1 year after parental irradiation (8.7 mGy/h, 5.2 Gy total dose). Using our previous F1-γ genome-wide gene expression data (GSE98539), hundreds of mRNAs were predicted as targets of differentially expressed (DE) miRNAs, involved in pathways such as insulin receptor, NFkB and PTEN signalling, linking to apoptosis and cancer. snRNAs belonging to the five major spliceosomal snRNAs were down-regulated in the F1-γ group, Indicating transcriptional and post-transcriptional alterations. In addition, DEpiRNA clusters were associated to 9 transposable elements (TEs) (LTR, LINE, and TIR) (p = 0.0024), probable as a response to the activation of these TEs. Moreover, the expression of the lincRNAs malat-1, and several others was altered in the offspring F1, in concordance with previously observed phenotypical alterations. In conclusion, our results demonstrate diverse gamma radiation-induced alterations in the ncRNA profiles of F1 offspring observable 1 year after parental irradiation.

Epigenetic, transcriptional and phenotypic responses in Daphnia magna exposed to low-level ionizing radiation.

Ionizing radiation is known to induce oxidative stress and DNA damage as well as epigenetic effects in aquatic organisms. Epigenetic changes can be part of the adaptive responses to protect organisms from radiation-induced damage, or act as drivers of toxicity pathways leading to adverse effects. To investigate the potential roles of epigenetic mechanisms in low-dose ionizing radiation-induced stress responses, an ecologically relevant crustacean, adult Daphnia magna were chronically exposed to low and medium level external 60Co gamma radiation ranging from 0.4, 1, 4, 10, and 40 mGy/h for seven days. Biological effects at the molecular (global DNA methylation, histone modification, gene expression), cellular (reactive oxygen species formation), tissue/organ (ovary, gut and epidermal histology) and organismal (fecundity) levels were investigated using a suite of effect assessment tools. The results showed an increase in global DNA methylation associated with loci-specific alterations of histone H3K9 methylation and acetylation, and downregulation of genes involved in DNA methylation, one-carbon metabolism, antioxidant defense, DNA repair, apoptosis, calcium signaling and endocrine regulation of development and reproduction. Temporal changes of reactive oxygen species (ROS) formation were also observed with an apparent transition from ROS suppression to induction from 2 to 7 days after gamma exposure. The cumulative fecundity, however, was not significantly changed by the gamma exposure. On the basis of the new experimental evidence and existing knowledge, a hypothetical model was proposed to provide in-depth mechanistic understanding of the roles of epigenetic mechanisms in low dose ionizing radiation induced stress responses in D. magna.

Inhibition of methyltransferase activity of enhancer of zeste 2 leads to enhanced lipid accumulation and altered chromatin status in zebrafish.

BACKGROUND: Recent studies indicate that exposure to environmental chemicals may increase susceptibility to developing metabolic diseases. This susceptibility may in part be caused by changes to the epigenetic landscape which consequently affect gene expression and lead to changes in lipid metabolism. The epigenetic modifier enhancer of zeste 2 (Ezh2) is a histone H3K27 methyltransferase implicated to play a role in lipid metabolism and adipogenesis. In this study, we used the zebrafish (Danio rerio) to investigate the role of Ezh2 on lipid metabolism and chromatin status following developmental exposure to the Ezh1/2 inhibitor PF-06726304 acetate. We used the environmental chemical tributyltin (TBT) as a positive control, as this chemical is known to act on lipid metabolism via EZH-mediated pathways in mammals. RESULTS: Zebrafish embryos (0-5 days post-fertilization, dpf) exposed to non-toxic concentrations of PF-06726304 acetate (5 µM) and TBT (1 nM) exhibited increased lipid accumulation. Changes in chromatin were analyzed by the assay for transposase-accessible chromatin sequencing (ATAC-seq) at 50% epiboly (5.5 hpf). We observed 349 altered chromatin regions, predominantly located at H3K27me3 loci and mostly more open chromatin in the exposed samples. Genes associated to these loci were linked to metabolic pathways. In addition, a selection of genes involved in lipid homeostasis, adipogenesis and genes specifically targeted by PF-06726304 acetate via altered chromatin accessibility were differentially expressed after TBT and PF-06726304 acetate exposure at 5 dpf, but not at 50% epiboly stage. One gene, cebpa, did not show a change in chromatin, but did show a change in gene expression at 5 dpf. Interestingly, underlying H3K27me3 marks were significantly decreased at this locus at 50% epiboly. CONCLUSIONS: Here, we show for the first time the applicability of ATAC-seq as a tool to investigate toxicological responses in zebrafish. Our analysis indicates that Ezh2 inhibition leads to a partial primed state of chromatin linked to metabolic pathways which results in gene expression changes later in development, leading to enhanced lipid accumulation. Although ATAC-seq seems promising, our in-depth assessment of the cebpa locus indicates that we need to consider underlying epigenetic marks as well.

Epigenetic, transcriptional and phenotypic responses in two generations of Daphnia magna exposed to the DNA methylation inhibitor 5-azacytidine.

The water flea Daphnia magna is a keystone species in freshwater ecosystems and has been widely used as a model organism in environmental ecotoxicology. This aquatic crustacean is sensitive to environmental stressors and displays considerable plasticity in adapting to changing environmental conditions. Part of this plasticity may be due to epigenetic regulation of gene expression, including changes to DNA methylation and histone modifications. Because of the generally hypomethylated genome of this species, we hypothesized that the histone code may have an essential role in the epigenetic control and that histone modifications might be an early marker for stress. This study aims to characterize the epigenetic, transcriptional and phenotypic responses and their causal linkages in directly exposed adult (F0) Daphnia and peritoneal exposed neonates (F1) after a chronic (7-day) exposure to a sublethal concentration (10 mg/l) of 5-azacytidine, a well-studied vertebrate DNA methylation inhibitor. Exposure of the F0 generation significantly reduced the cumulative fecundity, accompanied with differential expression of genes in the one-carbon-cycle metabolic pathway. In the epigenome of the F0 generation, a decrease in global DNA methylation, but no significant changes on H3K4me3 or H3K27me3, were observed. In the F1 offspring generation, changes in gene expression, a significant reduction in global DNA methylation and changes in histone modifications were identified. The results indicate that exposure during adulthood may result in more pronounced effects on early development in the offspring generation, though interpretation of the data should be carefully done since both the exposure regime and developmental period is different in the two generations examined. The obtained results improve our understanding of crustacean epigenetics and the tools developed may promote use of epigenetic markers in hazard assessment of environmental stressors.

Gamma radiation induces locus specific changes to histone modification enrichment in zebrafish and Atlantic salmon.

Ionizing radiation is a recognized genotoxic agent, however, little is known about the role of the functional form of DNA in these processes. Post translational modifications on histone proteins control the organization of chromatin and hence control transcriptional responses that ultimately affect the phenotype. The purpose of this study was to investigate effects on chromatin caused by ionizing radiation in fish. Direct exposure of zebrafish (Danio rerio) embryos to gamma radiation (10.9 mGy/h for 3h) induced hyper-enrichment of H3K4me3 at the genes hnf4a, gmnn and vegfab. A similar relative hyper-enrichment was seen at the hnf4a loci of irradiated Atlantic salmon (Salmo salar) embryos (30 mGy/h for 10 days). At the selected genes in ovaries of adult zebrafish irradiated during gametogenesis (8.7 and 53 mGy/h for 27 days), a reduced enrichment of H3K4me3 was observed, which was correlated with reduced levels of histone H3 was observed. F1 embryos of the exposed parents showed hyper-methylation of H3K4me3, H3K9me3 and H3K27me3 on the same three loci, while these differences were almost negligible in F2 embryos. Our results from three selected loci suggest that ionizing radiation can affect chromatin structure and organization, and that these changes can be detected in F1 offspring, but not in subsequent generations.

Ionizing radiation induces transgenerational effects of DNA methylation in zebrafish.

Ionizing radiation is known to cause DNA damage, yet the mechanisms underlying potential transgenerational effects of exposure have been scarcely studied. Previously, we observed effects in offspring of zebrafish exposed to gamma radiation during gametogenesis. Here, we hypothesize that these effects are accompanied by changes of DNA methylation possibly inherited by subsequent generations. We assessed DNA methylation in F1 embryos (5.5 hours post fertilization) with whole genome bisulfite sequencing following parental exposure to 8.7 mGy/h for 27 days and found 5658 differentially methylated regions (DMRs). DMRs were predominantly located at known regulatory regions, such as gene promoters and enhancers. Pathway analysis indicated the involvement of DMRs related to similar pathways found with gene expression analysis, such as development, apoptosis and cancers, which could be linked to previous observed developmental defects and genomic instability in the offspring. Follow up of 19 F1 DMRs in F2 and F3 embryos revealed persistent effects up to the F3 generation at 5 regions. These results indicate that ionizing radiation related effects in offspring can be linked to DNA methylation changes that partly can persist over generations. Monitoring DNA methylation could serve as a biomarker to provide an indication of ancestral exposures to ionizing radiation.

Parental exposure to gamma radiation causes progressively altered transcriptomes linked to adverse effects in zebrafish offspring.

Ionizing radiation causes a variety of effects, including DNA damage associated to cancers. However, the effects in progeny from irradiated parents is not well documented. Using zebrafish as a model, we previously found that parental exposure to ionizing radiation is associated with effects in offspring, such as increased hatching rates, deformities, increased DNA damage and reactive oxygen species. Here, we assessed short (one month) and long term effects (one year) on gene expression in embryonic offspring (5.5 h post fertilization) from zebrafish exposed during gametogenesis to gamma radiation (8.7 or 53 mGy/h for 27 days, total dose 5.2 or 31 Gy) using mRNA sequencing. One month after exposure, a global change in gene expression was observed in offspring from the 53 mGy/h group, followed by embryonic death at late gastrula, whereas offspring from the 8.7 mGy/h group was unaffected. Interestingly, one year after exposure newly derived embryos from the 8.7 mGy/h group exhibited 2390 (67.7% downregulated) differentially expressed genes. Overlaps in differentially expressed genes and enriched biological pathways were evident between the 53 mGy/h group one month and 8.7 mGy/h one year after exposure, but were oppositely regulated. Pathways could be linked to effects in adults and offspring, such as DNA damage (via Atm signaling) and reproduction (via Gnrh signaling). Comparison with gene expression analysis in directly exposed embryos indicate transferrin a and cytochrome P450 2x6 as possible biomarkers for radiation response in zebrafish. Our results indicate latent effects following ionizing radiation exposure from the lower dose in parents that can be transmitted to offspring and warrants monitoring effects over subsequent generations.

Epigenetic marking of the zebrafish developmental program.

A characteristic of anamniote development is a relatively long period of embryonic cell divisions in the absence of on-going transcription. In zebrafish, this period lasts for 10 cell cycles, or ∼3-h postfertilization, after which zygotic genome activation (ZGA) takes place during the midblastula transition. How the embryo establishes transcriptional competence and how ZGA is spatially and temporally regulated have not been examined until recently. We review here recent data on the transitions in DNA methylation and posttranslational histone modifications occurring during early zebrafish development, as the embryo acquires transcriptional competence and initiates its own gene expression program. We also address models accounting for the origin of epigenetic states detected in early embryos. From these observations, a concept of epigenetic prepatterning of the embryonic gene expression program prior to the onset of ZGA is emerging. The recent data collectively start shedding light on how ZGA may be programmed and regulated.

Epigenetic complexity during the zebrafish mid-blastula transition.

The zebrafish developmental transcription program is determined by temporal post-translational histone modifications established in a step-wise and combinatorial manner on specific promoters around the time of zygotic genome activation (ZGA). Here, we characterize this increasing epigenetic complexity before, during and after ZGA. H3K4me3/H3K27me3 co-enrichment prevails over H3K4me3/H3K9me3 at the time of ZGA. Whereas most H3K4me3-marked promoters are devoid of transcriptionally repressive H3K9me3 or H3K27me3, the latter marks rarely occur in absence of H3K4me3. On co-enriched genomic regions, H3K4me3 and H3K27me3 can overlap regardless of H3K9me3 enrichment, but H3K4me3 and H3K9me3 are mutually exclusive. H3K4me3 and H3K9me3 may however overlap only when H3K27me3 also marks the overlapping domain, suggesting that H3K27me3 may modulate chromatin states. On metagenes, H3K27me3 enrichment correlates with local alteration in H3K4me3 density, and co-enrichment in H3K9me3 is linked to alterations in both H3K27me3 and H3K4me3 profiles. This suggests physical proximity of these marks and supports a view of existence of bi- or tri-valent chromatin domains. Thus enrichment in trimethylated H3K9 or H3K27 is associated with local remodeling of chromatin manifested by changes in H3K4me3 density. We propose that metagenes can provide information on the multivalency of chromatin sates.

Prepatterning of developmental gene expression by modified histones before zygotic genome activation.

A hallmark of anamniote vertebrate development is a window of embryonic transcription-independent cell divisions before onset of zygotic genome activation (ZGA). Chromatin determinants of ZGA are unexplored; however, marking of developmental genes by modified histones in sperm suggests a predictive role of histone marks for ZGA. In zebrafish, pre-ZGA development for ten cell cycles provides an opportunity to examine whether genomic enrichment in modified histones is present before initiation of transcription. By profiling histone H3 trimethylation on all zebrafish promoters before and after ZGA, we demonstrate here an epigenetic prepatterning of developmental gene expression. This involves pre-ZGA marking of transcriptionally inactive genes involved in homeostatic and developmental regulation by permissive H3K4me3 with or without repressive H3K9me3 or H3K27me3. Our data suggest that histone modifications are instructive for the developmental gene expression program.

Tiling histone H3 lysine 4 and 27 methylation in zebrafish using high-density microarrays.

BACKGROUND: Uncovering epigenetic states by chromatin immunoprecipitation and microarray hybridization (ChIP-chip) has significantly contributed to the understanding of gene regulation at the genome-scale level. Many studies have been carried out in mice and humans; however limited high-resolution information exists to date for non-mammalian vertebrate species. PRINCIPAL FINDINGS: We report a 2.1-million feature high-resolution Nimblegen tiling microarray for ChIP-chip interrogations of epigenetic states in zebrafish (Danio rerio). The array covers 251 megabases of the genome at 92 base-pair resolution. It includes ∼15 kb of upstream regulatory sequences encompassing all RefSeq promoters, and over 5 kb in the 5' end of coding regions. We identify with high reproducibility, in a fibroblast cell line, promoters enriched in H3K4me3, H3K27me3 or co-enriched in both modifications. ChIP-qPCR and sequential ChIP experiments validate the ChIP-chip data and support the co-enrichment of trimethylated H3K4 and H3K27 on a subset of genes. H3K4me3- and/or H3K27me3-enriched genes are associated with distinct transcriptional status and are linked to distinct functional categories. CONCLUSIONS: We have designed and validated for the scientific community a comprehensive high-resolution tiling microarray for investigations of epigenetic states in zebrafish, a widely used developmental and disease model organism.

Chromatin states of developmentally-regulated genes revealed by DNA and histone methylation patterns in zebrafish embryos.

Embryo development proceeds from a cascade of gene activation and repression events controlled by epigenetic modifications of DNA and histones. Little is known about epigenetic states in the developing zebrafish, despite its importance as a model organism. We report here DNA methylation and histone modification profiles of promoters of developmentally-regulated genes (pou5f1, sox2, sox3, klf4, nnr, otx1b, nes, vasa), as well as tert and bactin2, in zebrafish embryos at the mid-late blastula transition, shortly after embryonic genome activation. We identify four classes of promoters based on the following profiles: (i) those enriched in marks of active genes (H3K9ac, H4ac, H3K4me3) without transcriptionally repressing H3K9me3 or H3K27me3; (ii) those enriched in H3K9ac, H4ac and H3K27me3, without H3K9me3; one such gene was klf4, shown by in situ hybridization to be mosaically expressed, likely accounting for the detection of both activating and repressive marks on its promoter; (iii) those enriched in H3K4me3 and H3K27me3 without acetylation; and (iv) those enriched in all histone modifications examined. Culture of embryo-derived cells under differentiation conditions leads to H3K9 and H4 deacetylation and H3K9 and H3K27 trimethylation on genes that are inactivated, yielding an epigenetic profile similar to those of fibroblasts or muscle. All promoters however retain H3K4me3, indicating an uncoupling of H3K4me3 occupancy and gene expression. All non-CpG island developmentally-regulated promoters are DNA unmethylated in embryos, but hypermethylated in fibroblasts. Our results suggest that differentially expressed embryonic genes are regulated by various patterns of histone modifications on unmethylated DNA, which create a developmentally permissive chromatin state.

Fish'n ChIPs: chromatin immunoprecipitation in the zebrafish embryo.

Chromatin immunoprecipitation (ChIP) is arguably the assay of choice to determine the genomic localization of DNA- or chromatin-binding proteins, including post-translationally modified histones, in cells. The increasing importance of the zebrafish, Danio rerio, as a model organism in functional genomics has recently sparked investigations of ChIP-based genome-scale mapping of modified histones on promoters, and studies on the role of specific transcription factors in developmental processes. ChIP assays used in these studies are cumbersome and conventionally require relatively large number of embryos. To simplify the procedure and to be able to apply the ChIP assay to reduced number of embryos, we re-evaluated the protocol for preparation of embryonic chromatin destined to ChIP. We found that manual homogenization of embryos rather than protease treatment to remove the chorion enhances ChIP efficiency and quickens the assay. We also incorporated key steps from a recently published ChIP assay for small cell numbers. We report here a protocol for immunoprecipitation of modified histones from mid-term blastula zebrafish embryos.

Histone H3 modifications associated with differentiation and long-term culture of mesenchymal adipose stem cells.

Long-term culture of mesenchymal stem cells leads to a loss of differentiation capacity, the molecular mechanism of which remains not understood. We show here that expansion of adipose stem cells (ASCs) to late passage (replicative senescence) is associated with promoter-specific and global changes in epigenetic histone modifications. In undifferentiated ASCs, inactive adipogenic and myogenic promoters are enriched in a repressive combination of trimethylated H3K4 (H3K4m3) and H3K27m3 in the absence of H3K9m3, a heterochromatin mark. Sequential chromatin immunoprecipitation assays indicate that H3K4m3 and H3K27m3 co-occupy a fraction of nucleosomes on some but not all lineage-specific promoters examined. However in cultured primary keratinocytes, adipogenic and myogenic promoters are enriched in trimethylated H3K4, K27, and K9, illustrating two distinct epigenetic states of inactive promoters related to potential for activation. H3K4m3 and H3K27m3 stably mark promoters during long-term ASC culture indicating that loss of differentiation capacity is not due to alterations in these histone modifications on these loci. Adipogenic differentiation in early passage results in H3K27 demethylation and H3K9 acetylation specifically on adipogenic promoters. On induction of differentiation in late passage, however, transcriptional upregulation is impaired, H3K27 trimethylation is maintained and H3K9 acetylation is inhibited on promoters. In addition, the polycomb proteins Ezh2 and Bmi1 are targeted to promoters. This correlates with global cellular Ezh2 increase and H3K9 deacetylation. Promoter targeting by Ezh2 and Bmi1 in late passage ASCs suggests the establishment of a polycomb-mediated epigenetic program aiming at repressing transcription.