Large hypomethylated domains serve as strong repressive machinery for key developmental genes in vertebrates

Cell cycle length governs heterochromatin reprogramming during early development in non-mammalian vertebrates

Heterochromatin marks such as H3K9me3 undergo global erasure and re-establishment after fertilization, and the proper reprogramming of H3K9me3 is essential for early development. Despite the widely conserved dynamics of heterochromatin reprogramming in invertebrates and non-mammalian vertebrates, previous studies have shown that the underlying mechanisms may differ between species. Here, we investigate the molecular mechanism of H3K9me3 dynamics in medaka (Japanese killifish, Oryzias latipes) as a non-mammalian vertebrate model, and show that rapid cell cycle during cleavage stages causes DNA replication-dependent passive erasure of H3K9me3. We also find that cell cycle slowing, toward the mid-blastula transition, permits increasing nuclear accumulation of H3K9me3 histone methyltransferase Setdb1, leading to the onset of H3K9me3 re-accumulation. We further demonstrate that cell cycle length in early development also governs H3K9me3 reprogramming in zebrafish and Xenopus laevis. Together with the previous studies in invertebrates, we propose that a cell cycle length-dependent mechanism for both global erasure and re-accumulation of H3K9me3 is conserved among rapid-cleavage species of non-mammalian vertebrates and invertebrates such as Drosophila, C. elegans, Xenopus and teleost fish.

Influence of DNA-methylation at multiple stages of limb chondrogenesis

Precise regulation of gene expression is of utmost importance during cell fate specification. DNA methylation is a key epigenetic mechanism that plays a significant role in the regulation of cell fate by recruiting repression proteins or inhibiting the binding of transcription factors to DNA to regulate gene expression. Limb development is a well-established model for understanding cell fate decisions, and the formation of skeletal elements is coordinated through a sequence of events that control chondrogenesis spatiotemporally. It has been established that epigenetic control participates in cartilage maturation. However, further investigation is required to determine its role in the earliest stages of chondrocyte differentiation. This study investigates how the DNA methylation environment affects cell fate divergence during the early chondrogenic events. Our research has shown for the first time that inhibiting DNA methylation in interdigital tissue with 5-azacytidine results in the formation of an ectopic digit. This discovery suggested that DNA methylation dynamics could regulate the fate of cells between chondrogenesis and cell death during autopod development. Our in vitro findings indicate that DNA methylation at the early stages of chondrogenesis is integral in regulating condensation by controlling cell adhesion and proapoptotic genes. As a result, the dynamics of methylation and demethylation are crucial in governing chondrogenesis and cell death during different stages of limb chondrogenesis.

Evolutionary conservation of embryonic DNA methylome remodelling in distantly related teleost species

Methylation of cytosines in the CG context (mCG) is the most abundant DNA modification in vertebrates that plays crucial roles in cellular differentiation and identity. After fertilization, DNA methylation patterns inherited from parental gametes are remodelled into a state compatible with embryogenesis. In mammals, this is achieved through the global erasure and re-establishment of DNA methylation patterns. However, in non-mammalian vertebrates like zebrafish, no global erasure has been observed. To investigate the evolutionary conservation and divergence of DNA methylation remodelling in teleosts, we generated base resolution DNA methylome datasets of developing medaka and medaka-zebrafish hybrid embryos. In contrast to previous reports, we show that medaka display comparable DNA methylome dynamics to zebrafish with high gametic mCG levels (sperm: ∼90%; egg: ∼75%), and adoption of a paternal-like methylome during early embryogenesis, with no signs of prior DNA methylation erasure. We also demonstrate that non-canonical DNA methylation (mCH) reprogramming at TGCT tandem repeats is a conserved feature of teleost embryogenesis. Lastly, we find remarkable evolutionary conservation of DNA methylation remodelling patterns in medaka-zebrafish hybrids, indicative of compatible DNA methylation maintenance machinery in far-related teleost species. Overall, these results suggest strong evolutionary conservation of DNA methylation remodelling pathways in teleosts, which is distinct from the global DNA methylome erasure and reestablishment observed in mammals.

Teratorn and its relatives - a cross-point of distinct mobile elements, transposons and viruses

Mobile genetic elements (e.g., transposable elements and plasmids) and viruses display significant diversity with various life cycles, but how this diversity emerges remains obscure. We previously reported a novel and giant (180 kb long) mobile element, Teratorn, originally identified in the genome of medaka, Oryzias latipes. Teratorn is a composite DNA transposon created by a fusion of a piggyBac-like DNA transposon (piggyBac) and a novel herpesvirus of the Alloherpesviridae family. Genomic survey revealed that Teratorn-like herpesviruses are widely distributed among teleost genomes, the majority of which are also fused with piggyBac, suggesting that fusion with piggyBac is a trigger for the life-cycle shift of authentic herpesviruses to an intragenomic parasite. Thus, Teratorn-like herpesvirus provides a clear example of how novel mobile elements emerge, that is to say, the creation of diversity. In this review, we discuss the unique sequence and life-cycle characteristics of Teratorn, followed by the evolutionary process of piggyBac-herpesvirus fusion based on the distribution of Teratorn-like herpesviruses (relatives) among teleosts. Finally, we provide other examples of evolutionary associations between different classes of elements and propose that recombination could be a driving force generating novel mobile elements.

Incomplete erasure of histone marks during epigenetic reprogramming in medaka early development

Epigenetic modifications undergo drastic erasure and reestablishment after fertilization. This reprogramming is required for proper embryonic development and cell differentiation. In mammals, some histone modifications are not completely reprogrammed and play critical roles in later development. In contrast, in nonmammalian vertebrates, most histone modifications are thought to be more intensively erased and reestablished by the stage of zygotic genome activation (ZGA). However, histone modifications that escape reprogramming in nonmammalian vertebrates and their potential functional roles remain unknown. Here, we quantitatively and comprehensively analyzed histone modification dynamics during epigenetic reprogramming in Japanese killifish, medaka (Oryzias latipes) embryos. Our data revealed that H3K27ac, H3K27me3, and H3K9me3 escape complete reprogramming, whereas H3K4 methylation is completely erased during cleavage stage. Furthermore, we experimentally showed the functional roles of such retained modifications at early stages: (i) H3K27ac premarks promoters during the cleavage stage, and inhibition of histone acetyltransferases disrupts proper patterning of H3K4 and H3K27 methylation at CpG-dense promoters, but does not affect chromatin accessibility after ZGA; (ii) H3K9me3 is globally erased but specifically retained at telomeric regions, which is required for maintenance of genomic stability during the cleavage stage. These results expand the understanding of diversity and conservation of reprogramming in vertebrates, and unveil previously uncharacterized functions of histone modifications retained during epigenetic reprogramming.

Targeted Manipulation of Histone Modification in Medaka Embryos

Recent development of targeted manipulation of histone modification enables us to experimentally and directly test the functional relevance of histone modifications accumulated at specific genomic regions. In particular, dCas9 epigenome editing has been widely used for site-specific manipulation of epigenetic modification. Here, we describe how to apply dCas9 epigenome editing in fish (medaka, Oryzias latipes) embryos and how to analyze induced changes in histone modification.

Genome-wide methylation in the panmictic European eel (Anguilla anguilla)

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, likely reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggests interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight hox clusters, and the results suggest important roles of methylation at hox genes for adaptive processes.

Wnt11 acts on dermomyotome cells to guide epaxial myotome morphogenesis

The dorsal axial muscles, or epaxial muscles, are a fundamental structure covering the spinal cord and vertebrae, as well as mobilizing the vertebrate trunk. To date, mechanisms underlying the morphogenetic process shaping the epaxial myotome are largely unknown. To address this, we used the medaka zic1/zic4-enhancer mutant Double anal fin (Da), which exhibits ventralized dorsal trunk structures resulting in impaired epaxial myotome morphology and incomplete coverage over the neural tube. In wild type, dorsal dermomyotome (DM) cells reduce their proliferative activity after somitogenesis. Subsequently, a subset of DM cells, which does not differentiate into the myotome population, begins to form unique large protrusions extending dorsally to guide the epaxial myotome dorsally. In Da, by contrast, DM cells maintain the high proliferative activity and mainly form small protrusions. By combining RNA- and ChIP-sequencing analyses, we revealed direct targets of Zic1, which are specifically expressed in dorsal somites and involved in various aspects of development, such as cell migration, extracellular matrix organization, and cell-cell communication. Among these, we identified wnt11 as a crucial factor regulating both cell proliferation and protrusive activity of DM cells. We propose that dorsal extension of the epaxial myotome is guided by a non-myogenic subpopulation of DM cells and that wnt11 empowers the DM cells to drive the coverage of the neural tube by the epaxial myotome.

Comparative Study in Zebrafish and Medaka Unravels the Mechanisms of Tissue Regeneration

Tissue regeneration has been in the spotlight of research for its fascinating nature and potential applications in human diseases. The trait of regenerative capacity occurs diversely across species and tissue contexts, while it seems to decline over evolution. Organisms with variable regenerative capacity are usually distinct in phylogeny, anatomy, and physiology. This phenomenon hinders the feasibility of studying tissue regeneration by directly comparing regenerative with non-regenerative animals, such as zebrafish (Danio rerio) and mice (Mus musculus). Medaka (Oryzias latipes) is a fish model with a complete reference genome and shares a common ancestor with zebrafish approximately 110–200 million years ago (compared to 650 million years with mice). Medaka shares similar features with zebrafish, including size, diet, organ system, gross anatomy, and living environment. However, while zebrafish regenerate almost every organ upon experimental injury, medaka shows uneven regenerative capacity. Their common and distinct biological features make them a unique platform for reciprocal analyses to understand the mechanisms of tissue regeneration. Here we summarize current knowledge about tissue regeneration in these fish models in terms of injured tissues, repairing mechanisms, available materials, and established technologies. We further highlight the concept of inter-species and inter-organ comparisons, which may reveal mechanistic insights and hint at therapeutic strategies for human diseases.

CTCF looping is established during gastrulation in medaka embryos

Chromatin looping plays an important role in genome regulation. However, because ChIP-seq and loop-resolution Hi-C (DNA-DNA proximity ligation) are extremely challenging in mammalian early embryos, the developmental stage at which cohesin-mediated loops form remains unknown. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation (the onset of cell differentiation) and characterize transcription, protein binding, and genome architecture. We find that gastrulation is associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and a large increase in the size of contact domains. In contrast, the binding of the CTCF is fixed throughout embryogenesis. Loops form long after genome-wide transcriptional activation, and long after domain formation seen in mouse embryos. These results suggest that, although loops may play a role in differentiation, they are not required for zygotic transcription. When we repeated our experiments in zebrafish, loops did not emerge until gastrulation, that is, well after zygotic genome activation. We observe that loop positions are highly conserved in synteny blocks of medaka and zebrafish, indicating that the 3D genome architecture has been maintained for >110–200 million years of evolution.

Conserved Epigenetic Regulatory Logic Infers Genes Governing Cell Identity

Determining genes that orchestrate cell differentiation in development and disease remains a fundamental goal of cell biology. This study establishes a genome-wide metric based on the gene-repressive trimethylation of histone H3 at lysine 27 (H3K27me3) across hundreds of diverse cell types to identify genetic regulators of cell differentiation. We introduce a computational method, TRIAGE, which uses discordance between gene-repressive tendency and expression to identify genetic drivers of cell identity. We apply TRIAGE to millions of genome-wide single-cell transcriptomes, diverse omics platforms, and eukaryotic cells and tissue types. Using a wide range of data, we validate the performance of TRIAGE in identifying cell-type-specific regulatory factors across diverse species including human, mouse, boar, bird, fish, and tunicate. Using CRISPR gene editing, we use TRIAGE to experimentally validate RNF220 as a regulator of Ciona cardiopharyngeal development and SIX3 as required for differentiation of endoderm in human pluripotent stem cells. A record of this paper’s transparent peer review process is included in the Supplemental Information.

Perturbing genes controlling cell decisions have major implications in development or disease. However, identifying key regulatory genes from the thousands expressed in a cell is challenging. TRIAGE is a computational method that distills patterns of epigenetic repression across diverse cell types to infer regulatory genes using input gene expression data from any cell type. Demonstrating its utility, we combine single-cell RNA-seq and TRIAGE to identify and experimentally confirm novel regulators of heart development in evolutionarily distant species.

Dynamic transcriptional and chromatin accessibility landscape of medaka embryogenesis

Medaka (Oryzias latipes) has become an important vertebrate model widely used in genetics, developmental biology, environmental sciences, and many other fields. A high-quality genome sequence and a variety of genetic tools are available for this model organism. However, existing genome annotation is still rudimentary, as it was mainly based on computational prediction and short-read RNA-seq data. Here we report a dynamic transcriptome landscape of medaka embryogenesis profiled by long-read RNA-seq, short-read RNA-seq, and ATAC-seq. By integrating these data sets, we constructed a much-improved gene model set including about 17,000 novel isoforms and identified 1600 transcription factors, 1100 long noncoding RNAs, and 150,000 potential cis-regulatory elements as well. Time-series data sets provided another dimension of information. With the expression dynamics of genes and accessibility dynamics of cis-regulatory elements, we investigated isoform switching, as well as regulatory logic between accessible elements and genes, during embryogenesis. We built a user-friendly medaka omics data portal to present these data sets. This resource provides the first comprehensive omics data sets of medaka embryogenesis. Ultimately, we term these three assays as the minimum ENCODE toolbox and propose the use of it as the initial and essential profiling genomic assays for model organisms that have limited data available. This work will be of great value for the research community using medaka as the model organism and many others as well.

Targeted in vivo epigenome editing of H3K27me3

BACKGROUND

Epigenetic modifications have a central role in transcriptional regulation. While several studies using next-generation sequencing have revealed genome-wide associations between epigenetic modifications and transcriptional states, a direct causal relationship at specific genomic loci has not been fully demonstrated, due to a lack of technology for targeted manipulation of epigenetic modifications. Recently, epigenome editing techniques based on the CRISPR-Cas9 system have been reported to directly manipulate specific modifications at precise genomic regions. However, the number of editable modifications as well as studies applying these techniques in vivo is still limited.

RESULTS

Here, we report direct modification of the epigenome in medaka (Japanese killifish, Oryzias latipes) embryos. Specifically, we developed a method to ectopically induce the repressive histone modification, H3K27me3 in a locus-specific manner, using a fusion construct of Oryzias latipes H3K27 methyltransferase Ezh2 (olEzh2) and dCas9 (dCas9-olEzh2). Co-injection of dCas9-olEzh2 mRNA with single guide RNAs (sgRNAs) into one-cell-stage embryos induced specific H3K27me3 accumulation at the targeted loci and induced downregulation of gene expression.

CONCLUSION

In this study, we established the in vivo epigenome editing of H3K27me3 using medaka embryos. The locus-specific manipulation of the epigenome in living organisms will lead to a previously inaccessible understanding of the role of epigenetic modifications in development and disease.

In Vivo Analysis of Embryo Development and Behavioral Response of Medaka Fish under Static Magnetic Field Exposures

The static magnetic field (SMF) in human exposure has become a health risk concern, especially with respect to prolonged exposure. The International Commission on Non-Ionizing Radiation Protection (ICNIRP) has been considering cell or animal models to be adopted to estimate the possible human health impacts after such exposure. The medaka fish is a good animal model for human-related health assessment studies; this paper examines both the embryo development and behavioral responses in medaka fish in vivo to long-term SMF exposure at the mT level. SMF exposure was examined for the complete developmental period of embryos until hatched; the embryos were monitored and recorded every 24 h for different morphological abnormalities in their developmental stages. The behavioral response of adult fish was also examined by analyzing their swimming velocities and positioning as compared with that of the control group. It was observed that there were no impacts on embryo development under prolonged exposure up to about 100 mT while the swimming behavior of the adult fish under exposure was different to the control group-the swimming movement of the treated group was more static, with an average velocity of 24.6% less as observed over a 24-h duration.

The DNA Methylation Landscape of Stickleback Reveals Patterns of Sex Chromosome Evolution and Effects of Environmental Salinity

Epigenetic mechanisms such as DNA methylation are a key component of dosage compensation on sex chromosomes and have been proposed as an important source of phenotypic variation influencing plasticity and adaptive evolutionary processes, yet little is known about the role of DNA methylation in an ecological or evolutionary context in vertebrates. The threespine stickleback (Gasterosteus aculeatus) is an ecological and evolutionary model system that has been used to study mechanisms involved in the evolution of adaptive phenotypes in novel environments as well as the evolution heteromorphic sex chromosomes and dosage compensation in vertebrates. Using whole genome bisulfite sequencing, we compared genome-wide DNA methylation patterns between threespine stickleback males and females and between stickleback reared at different environmental salinities. Apparent hypermethylation of the younger evolutionary stratum of the stickleback X chromosome in females relative to males suggests a potential role of DNA methylation in the evolution of heteromorphic sex chromosomes. We also demonstrate that rearing salinity has genome-wide effects on DNA methylation levels, which has the potential to lead to the accumulation of epigenetic variation between natural populations in different environments.

Differential hypomethylation of the repetitive Tol2/Alu-rich sequences in the genome of Bodianus species (Labriformes, Labridae)

Representatives of the order Labriformes show karyotypes of extreme conservatism together with others with high chromosomal diversification. However, the cytological characterization of epigenetic modifications remains unknown for the majority of the species. In the family Labridae, the most abundant fishes on tropical reefs, the genomes of the genus Bodianus Bloch, 1790 have been characterized by the occurrence of a peculiar chromosomal region, here denominated BOD. This region is exceptionally decondensed, heterochromatic, argentophilic, GC-neutral and, in contrast to classical secondary constrictions, shows no signals of hybridization with 18S rDNA probes. In order to characterize the BOD region, the methylation pattern, the distribution of Alu and Tol2 retrotransposons and of 18S and 5S rDNA sites, respectively, were analyzed by Fluorescence In Situ Hybridization (FISH) on metaphase chromosomes of two Bodianus species, B. insularis Gomon & Lubbock, 1980 and B. pulchellus (Poey, 1860). Immunolocalization of the 5-methylcytosine revealed hypermethylated chromosomal regions, dispersed along the entire length of the chromosomes of both species, while the BOD regions exhibited a hypomethylated pattern. Hypomethylation of the BOD region is associated with the precise co-location of Tol2 and Alu elements, suggesting their active participation in the regulatory epigenetic process. This evidence underscores a probable differential methylation action during the cell cycle, as well as the role of Tol2/Alu elements in functional processes of fish genomes.

Unlinking the methylome pattern from nucleotide sequence, revealed by large-scale in vivo genome engineering and methylome editing in medaka fish

The heavily methylated vertebrate genomes are punctuated by stretches of poorly methylated DNA sequences that usually mark gene regulatory regions. It is known that the methylation state of these regions confers transcriptional control over their associated genes. Given its governance on the transcriptome, cellular functions and identity, genome-wide DNA methylation pattern is tightly regulated and evidently predefined. However, how is the methylation pattern determined in vivo remains enigmatic. Based on in silico and in vitro evidence, recent studies proposed that the regional hypomethylated state is primarily determined by local DNA sequence, e.g., high CpG density and presence of specific transcription factor binding sites. Nonetheless, the dependency of DNA methylation on nucleotide sequence has not been carefully validated in vertebrates in vivo. Herein, with the use of medaka (Oryzias latipes) as a model, the sequence dependency of DNA methylation was intensively tested in vivo. Our statistical modeling confirmed the strong statistical association between nucleotide sequence pattern and methylation state in the medaka genome. However, by manipulating the methylation state of a number of genomic sequences and reintegrating them into medaka embryos, we demonstrated that artificially conferred DNA methylation states were predominantly and robustly maintained in vivo, regardless of their sequences and endogenous states. This feature was also observed in the medaka transgene that had passed across generations. Thus, despite the observed statistical association, nucleotide sequence was unable to autonomously determine its own methylation state in medaka in vivo. Our results apparently argue against the notion of the governance on the DNA methylation by nucleotide sequence, but instead suggest the involvement of other epigenetic factors in defining and maintaining the DNA methylation landscape. Further investigation in other vertebrate models in vivo will be needed for the generalization of our observations made in medaka.

The genomes of vertebrate animals are naturally and extensively modified by methylation. The DNA methylation is essential to normal functions of cells, hence the whole animal, since it governs gene expression. Defects in the establishment and maintenance of proper methylation pattern are commonly associated with various developmental abnormalities and diseases. How exactly is the normal pattern defined in vertebrate animals is not fully understood, but recent researches with computational analyses and cultured cells suggested that DNA sequence is a primary determinant of the methylation pattern. This study encompasses the first experiments that rigorously test this notion in whole animal (medaka fish). In statistical sense, we observed the very strong correlation between DNA sequence and methylation state. However, by introducing unmethylated and artificially methylated native genomic DNA sequences into the genome, we demonstrated that the artificially conferred methylation states were robustly maintained in the animal, independent of the sequence and native state. Our results thus demonstrate that genome-wide DNA methylation pattern is not autonomously determined by the DNA sequence, which underpins the vital role of DNA methylation pattern as a core epigenetic element.

Centromere evolution and CpG methylation during vertebrate speciation

Centromeres and large-scale structural variants evolve and contribute to genome diversity during vertebrate speciation. Here, we perform de novo long-read genome assembly of three inbred medaka strains that are derived from geographically isolated subpopulations and undergo speciation. Using single-molecule real-time (SMRT) sequencing, we obtain three chromosome-mapped genomes of length ~734, ~678, and ~744Mbp with a resource of twenty-two centromeric regions of length 20-345kbp. Centromeres are positionally conserved among the three strains and even between four pairs of chromosomes that were duplicated by the teleost-specific whole-genome duplication 320-350 million years ago. The centromeres do not all evolve at a similar pace; rather, centromeric monomers in non-acrocentric chromosomes evolve significantly faster than those in acrocentric chromosomes. Using methylation sensitive SMRT reads, we uncover centromeres are mostly hypermethylated but have hypomethylated sub-regions that acquire unique sequence compositions independently. These findings reveal the potential of non-acrocentric centromere evolution to contribute to speciation.

Centromere evolution and CpG methylation during vertebrate speciation

Centromeres and large-scale structural variants evolve and contribute to genome diversity during vertebrate speciation. Here, we perform de novo long-read genome assembly of three inbred medaka strains that are derived from geographically isolated subpopulations and undergo speciation. Using single-molecule real-time (SMRT) sequencing, we obtain three chromosome-mapped genomes of length ~734, ~678, and ~744Mbp with a resource of twenty-two centromeric regions of length 20–345kbp. Centromeres are positionally conserved among the three strains and even between four pairs of chromosomes that were duplicated by the teleost-specific whole-genome duplication 320–350 million years ago. The centromeres do not all evolve at a similar pace; rather, centromeric monomers in non-acrocentric chromosomes evolve significantly faster than those in acrocentric chromosomes. Using methylation sensitive SMRT reads, we uncover centromeres are mostly hypermethylated but have hypomethylated sub-regions that acquire unique sequence compositions independently. These findings reveal the potential of non-acrocentric centromere evolution to contribute to speciation.

Centromeres and large-scale structural variants evolve and contribute to genome diversity during vertebrate speciation. Here Ichikawa et al perform de novo long-read genome assembly of three inbred medaka strains, and report long-range structure of centromeres and their methylation as well as correlation of structural variants with differential gene expression.

Hypomethylated domain-enriched DNA motifs prepattern the accessible nucleosome organization in teleosts

Background

Gene promoters in vertebrate genomes show distinct chromatin features such as stably positioned nucleosome array and DNA hypomethylation. The nucleosomes are known to have certain sequence preferences, and the prediction of nucleosome positioning from DNA sequence has been successful in some organisms such as yeast. However, at gene promoters where nucleosomes are much more stably positioned than in other regions, the sequence-based model has failed to work well, and sequence-independent mechanisms have been proposed.

Results

Using DNase I-seq in medaka embryos, we demonstrated that hypomethylated domains (HMDs) specifically possess accessible nucleosome organization with longer linkers, and we reassessed the DNA sequence preference for nucleosome positioning in these specific regions. Remarkably, we found with a supervised machine learning algorithm, k-mer SVM, that nucleosome positioning in HMDs is accurately predictable from DNA sequence alone. Specific short sequences (6-mers) that contribute to the prediction are specifically enriched in HMDs and distribute periodically with approximately 200-bp intervals which prepattern the position of accessible linkers. Surprisingly, the sequence preference of the nucleosome and linker in HMDs is opposite from that reported previously. Furthermore, the periodicity of specific motifs at hypomethylated promoters was conserved in zebrafish.

Conclusion

This study reveals strong link between nucleosome positioning and DNA sequence at vertebrate promoters, and we propose hypomethylated DNA-specific regulation of nucleosome positioning.

Electronic supplementary material

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

Complete fusion of a transposon and herpesvirus created the Teratorn mobile element in medaka fish

Mobile genetic elements (e.g., transposable elements and viruses) display significant diversity with various life cycles, but how novel elements emerge remains obscure. Here, we report a giant (180-kb long) transposon, Teratorn, originally identified in the genome of medaka, Oryzias latipes. Teratorn belongs to the piggyBac superfamily and retains the transposition activity. Remarkably, Teratorn is largely derived from a herpesvirus of the Alloherpesviridae family that could infect fish and amphibians. Genomic survey of Teratorn-like elements reveals that some of them exist as a fused form between piggyBac transposon and herpesvirus genome in teleosts, implying the generality of transposon-herpesvirus fusion. We propose that Teratorn was created by a unique fusion of DNA transposon and herpesvirus, leading to life cycle shift. Our study supports the idea that recombination is the key event in generation of novel mobile genetic elements.

Teratorn is a large mobile genetic element originally identified in the small teleost fish medaka. Here, the authors show that Teratorn is derived from the fusion of a piggyBac superfamily DNA transposon and an alloherpesvirus and that it is widely found across teleost fish.

Individual recognition and the 'face inversion effect' in medaka fish (Oryzias latipes)

Individual recognition (IR) is essential for maintaining various social interactions in a group, and face recognition is one of the most specialised cognitive abilities in IR. We used both a mating preference system and an electric shock conditioning experiment to test IR ability in medaka, and found that signals near the face are important. Medaka required more time to discriminate vertically inverted faces, but not horizontally shifted faces or inverted non-face objects. The ability may be comparable to the classic 'face inversion effect' in humans and some other mammals. Extra patterns added to the face also did not influence the IR. These findings suggest the possibility that the process of face recognition may differ from that used for other objects. The complex form of recognition may promote specific processing adaptations, although the mechanisms and neurological bases might differ in mammals and medaka. The ability to recognise other individuals is important for shaping animal societies.

Dynamics of DNA methylomes underlie oyster development

DNA methylation is a critical epigenetic regulator of development in mammals and social insects, but its significance in development outside these groups is not understood. Here we investigated the genome-wide dynamics of DNA methylation in a mollusc model, the oyster Crassostrea gigas, from the egg to the completion of organogenesis. Large-scale methylation maps reveal that the oyster genome displays a succession of methylated and non methylated regions, which persist throughout development. Differentially methylated regions (DMRs) are strongly regulated during cleavage and metamorphosis. The distribution and levels of methylated DNA within genomic features (exons, introns, promoters, repeats and transposons) show different developmental lansdscapes marked by a strong increase in the methylation of exons against introns after metamorphosis. Kinetics of methylation in gene-bodies correlate to their transcription regulation and to distinct functional gene clusters, and DMRs at cleavage and metamorphosis bear the genes functionally related to these steps, respectively. This study shows that DNA methylome dynamics underlie development through transcription regulation in the oyster, a lophotrochozoan species. To our knowledge, this is the first demonstration of such epigenetic regulation outside vertebrates and ecdysozoan models, bringing new insights into the evolution and the epigenetic regulation of developmental processes.

Elucidating the mechanisms which govern the development of multicellular animals and their evolution is a fundamental task. Epigenetic mechanisms like DNA methylation have recently emerged as critical regulators of mammalian development through the control of genes that determine the identity of cells and the transmission of parental imprints. In invertebrates however, DNA is mostly unmethylated and does not play a role in development except in the peculiar case of social insects. Therefore the significance of DNA methylation in development is thought to be restricted to vertebrates, and thereby considered a recent evolutionary acquisition, and the situation in more distant organisms is unknown. Here we investigated the dynamics of genome-wide DNA methylation patterns in a mollusc, the oyster C. gigas, throughout its development. We found that the dynamics of DNA methylation correspond to the expression dynamics of distinct functional gene clusters that control two critical development steps, cleavage and metamorphosis, and we provide insights into the underlying molecular mechanisms in a non-vertebrate species. These findings challenge the present considerations on the evolution of developmental processes and their epigenetic regulation, and open a new area of research in molecular and developmental biology in invertebrates.

Comparative Analysis of Genome and Epigenome in Closely Related Medaka Species Identifies Conserved Sequence Preferences for DNA Hypomethylated Domains

The genomes of vertebrates are globally methylated, but a small portion of genomic regions are known to be hypomethylated. Although hypomethylated domains (HMDs) have been implicated in transcriptional regulation in various ways, how a HMD is determined in a particular genomic region remains elusive. To search for DNA motifs essential for the formation of HMDs, we performed the genome-wide comparative analysis of genome and DNA methylation patterns of the two medaka inbred lines, Hd-rRII1 and HNI-II, which are derived from northern and southern subpopulations of Japan and exhibit high levels of genetic variations (SNP, ∼ 3%). We successfully mapped > 70% of HMDs in both genomes and found that the majority of those mapped HMDs are conserved between the two lines (common HMDs). Unexpectedly, the average genetic variations are similar in the common HMD and other genome regions. However, we identified short well-conserved motifs that are specifically enriched in HMDs, suggesting that they may play roles in the establishment of HMDs in the medaka genome.

Toward high-resolution population genomics using archaeological samples

The term ‘ancient DNA’ (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of ‘molecular paleontology’. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.

Medaka as a model for studying environmentally induced epigenetic transgenerational inheritance of phenotypes

Ability of environmental stressors to induce transgenerational diseases has been experimentally demonstrated in plants, worms, fish, and mammals, indicating that exposures affect not only human health but also fish and ecosystem health. Small aquarium fish have been reliable model to study genetic and epigenetic basis of development and disease. Additionally, fish can also provide better, economic opportunity to study transgenerational inheritance of adverse health and epigenetic mechanisms. Molecular mechanisms underlying germ cell development in fish are comparable to those in mammals and humans. This review will provide a short overview of long-term effects of environmental chemical contaminant exposure in various models, associated epigenetic mechanisms, and a perspective on fish as model to study environmentally induced transgenerational inheritance of altered phenotypes.

Tet1 and Tet2 Protect DNA Methylation Canyons against Hypermethylation

DNA methylation is a dynamic epigenetic modification with an important role in cell fate specification and reprogramming. The Ten eleven translocation (Tet) family of enzymes converts 5-methylcytosine to 5-hydroxymethylcytosine, which promotes passive DNA demethylation and functions as an intermediate in an active DNA demethylation process. Tet1/Tet2 double-knockout mice are characterized by developmental defects and epigenetic instability, suggesting a requirement for Tet-mediated DNA demethylation for the proper regulation of gene expression during differentiation. Here, we used whole-genome bisulfite and transcriptome sequencing to characterize the underlying mechanisms. Our results uncover the hypermethylation of DNA methylation canyons as the genomic key feature of Tet1/Tet2 double-knockout mouse embryonic fibroblasts. Canyon hypermethylation coincided with disturbed regulation of associated genes, suggesting a mechanistic explanation for the observed Tet-dependent differentiation defects. Based on these results, we propose an important regulatory role of Tet-dependent DNA demethylation for the maintenance of DNA methylation canyons, which prevents invasive DNA methylation and allows functional regulation of canyon-associated genes.

LSD1/KDM1A promotes hematopoietic commitment of hemangioblasts through downregulation of Etv2

The hemangioblast is a progenitor cell with the capacity to give rise to both hematopoietic and endothelial progenitors. Currently, the regulatory mechanisms underlying hemangioblast formation are being elucidated, whereas those controllers for the selection of hematopoietic or endothelial fates still remain a mystery. To answer these questions, we screened for zebrafish mutants that have defects in the hemangioblast expression of Gata1, which is never expressed in endothelial progenitors. One of the isolated mutants, it627, showed not only down-regulation of hematopoietic genes but also up-regulation of endothelial genes. We identified the gene responsible for the it627 mutant as the zebrafish homolog of Lys-specific demethylase 1 (LSD1/KDM1A). Surprisingly, the hematopoietic defects in lsd1(it627) embryos were rescued by the gene knockdown of the Ets variant 2 gene (etv2), an essential regulator for vasculogenesis. Our results suggest that the LSD1-dependent shutdown of Etv2 gene expression may be a significant event required for hemangioblasts to initiate hematopoietic differentiation.

Genome-wide epigenetic cross-talk between DNA methylation and H3K27me3 in zebrafish embryos

DNA methylation and histone modifications are epigenetic marks implicated in the complex regulation of vertebrate embryogenesis. The cross-talk between DNA methylation and Polycomb-dependent H3K27me3 histone mark has been reported in a number of organisms [1], [2], [3], [4], [5], [6], [7] and both marks are known to be required for proper developmental progression. Here we provide genome-wide DNA methylation (MethylCap-seq) and H3K27me3 (ChIP-seq) maps for three stages (dome, 24 hpf and 48 hpf) of zebrafish (Danio rerio) embryogenesis, as well as all analytical and methodological details associated with the generation of this dataset. We observe a strong antagonism between the two epigenetic marks present in CpG islands and their compatibility throughout the bulk of the genome, as previously reported in mammalian ESC lines (Brinkman et al., 2012). Next generation sequencing data linked to this project have been deposited in the Gene Expression Omnibus (GEO) database under accession numbers GSE35050 and GSE70847.

A linear time algorithm for detecting long genomic regions enriched with a specific combination of epigenetic states

Background

Epigenetic modifications are essential for controlling gene expression. Recent studies have shown that not only single epigenetic modifications but also combinations of multiple epigenetic modifications play vital roles in gene regulation. A striking example is the long hypomethylated regions enriched with modified H3K27me3 (called, "K27HMD" regions), which are exposed to suppress the expression of key developmental genes relevant to cellular development and differentiation during embryonic stages in vertebrates. It is thus a biologically important issue to develop an effective optimization algorithm for detecting long DNA regions (e.g., >4 kbp in size) that harbor a specific combination of epigenetic modifications (e.g., K27HMD regions). However, to date, optimization algorithms for these purposes have received little attention, and available methods are still heuristic and ad hoc.

Results

In this paper, we propose a linear time algorithm for calculating a set of non-overlapping regions that maximizes the sum of similarities between the vector of focal epigenetic states and the vectors of raw epigenetic states at DNA positions in the set of regions. The average elapsed time to process the epigenetic data of any of human chromosomes was less than 2 seconds on an Intel Xeon CPU. To demonstrate the effectiveness of the algorithm, we estimated large K27HMD regions in the medaka and human genomes using our method, ChromHMM, and a heuristic method.

Conclusions

We confirmed that the advantages of our method over those of the two other methods. Our method is flexible enough to handle other types of epigenetic combinations. The program that implements the method is called "CSMinfinder" and is made available at: http://mlab.cb.k.u-tokyo.ac.jp/~ichikawa/Segmentation/