Epigenetically distinct synaptic architecture in clonal compartments in the teleostean dorsal pallium.

The dorsal telencephalon (i.e. the pallium) exhibits high anatomical diversity across vertebrate classes. The non-mammalian dorsal pallium accommodates various compartmentalized structures among species. The developmental, functional, and evolutional diversity of the dorsal pallium remain unillustrated. Here, we analyzed the structure and epigenetic landscapes of cell lineages in the telencephalon of medaka fish (Oryzias latipes) that possesses a clearly delineated dorsal pallium (Dd2). We found that pallial anatomical regions, including Dd2, are formed by mutually exclusive clonal units, and that each pallium compartment exhibits a distinct epigenetic landscape. In particular, Dd2 possesses a unique open chromatin pattern that preferentially targets synaptic genes. Indeed, Dd2 shows a high density of synapses. Finally, we identified several transcription factors as candidate regulators. Taken together, we suggest that cell lineages are the basic components for the functional regionalization in the pallial anatomical compartments and that their changes have been the driving force for evolutionary diversity.

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.

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.

Methylome inheritance and enhancer dememorization reset an epigenetic gate safeguarding embryonic programs.

Marked epigenetic reprogramming is essential to convert terminally differentiated gametes to totipotent embryos. It remains puzzling why postfertilization global DNA reprogramming occurs in mammals but not in nonmammalian vertebrates. In zebrafish, global methylome inheritance is however accompanied by extensive enhancer "dememorization" as they become fully methylated. By depleting maternal dnmt1 using oocyte microinjection, we eliminated DNA methylation in early embryos, which died around gastrulation with severe differentiation defects. Notably, methylation deficiency leads to derepression of adult tissue­specific genes and CG-rich enhancers, which acquire ectopic transcription factor binding and, unexpectedly, histone H3 lysine 4 trimethylation (H3K4me3). By contrast, embryonic enhancers are generally CG-poor and evade DNA methylation repression. Hence, global DNA hypermethylation inheritance coupled with enhancer dememorization installs an epigenetic gate that safeguards embryonic programs and ensures temporally ordered gene expression. We propose that "enhancer dememorization" underlies and unifies distinct epigenetic reprogramming modes in early development between mammals and nonmammals.

Mobility of the forearm skeleton in the Asiatic black (Ursus thibetanus), brown (U. arctos) and polar (U. maritimus) bears.

In several primates and carnivores, pronation/supination angles of the forearm skeleton were examined, and it is thought that a larger angle is useful to acquire dexterous behaviors in feeding and/or life style, including climbing. In this study, the pronation/supination angles in Asiatic black, brown and polar bears were nondestructively examined. These specimens were classified as adult or non-adult. Three or four carcasses of each group of Asiatic black and brown bears were used for CT analysis, whereas only one adult polar bear was used. The forearms were positioned within the gantry of a CT scanner in both maximally supinated and pronated states. Extracted cross-sectional CT images of two positions were superimposed by overlapping the outlines of each ulna. The centroids of the radii were detected, and then the centroid of each radius and the midpoint of a line which connects between both ends of the surface of each radius facing the ulna, were connected by lines to measure the angle of rotation as an index of pronation/supination. In adult brown and polar bears, the angles were smaller as compared with the other groups (Asiatic black and non-adult brown bears). Asiatic black and non-adult brown bears can climb trees, whereas adult brown bears and polar bears cannot. This suggests that the pronation/supination angle is related to arboreal activity in Ursidae.

Chromatin architecture transitions from zebrafish sperm through early embryogenesis.

Chromatin architecture mapping in 3D formats has increased our understanding of how regulatory sequences and gene expression are connected and regulated in a genome. The 3D chromatin genome shows extensive remodeling during embryonic development, and although the cleavage-stage embryos of most species lack structure before zygotic genome activation (pre-ZGA), zebrafish has been reported to have structure. Here, we aimed to determine the chromosomal architecture in paternal/sperm zebrafish gamete cells to discern whether it either resembles or informs early pre-ZGA zebrafish embryo chromatin architecture. First, we assessed the higher-order architecture through advanced low-cell in situ Hi-C. The structure of zebrafish sperm, packaged by histones, lacks topological associated domains and instead displays "hinge-like" domains of ∼150 kb that repeat every 1-2 Mbs, suggesting a condensed repeating structure resembling mitotic chromosomes. The pre-ZGA embryos lacked chromosomal structure, in contrast to prior work, and only developed structure post-ZGA. During post-ZGA, we find chromatin architecture beginning to form at small contact domains of a median length of ∼90 kb. These small contact domains are established at enhancers, including super-enhancers, and chemical inhibition of Ep300a (p300) and Crebbpa (CBP) activity, lowering histone H3K27ac, but not transcription inhibition, diminishes these contacts. Together, this study reveals hinge-like domains in histone-packaged zebrafish sperm chromatin and determines that the initial formation of high-order chromatin architecture in zebrafish embryos occurs after ZGA primarily at enhancers bearing high H3K27ac.

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.

Recent evolution of a TET-controlled and DPPA3/STELLA-driven pathway of passive DNA demethylation in mammals.

Genome-wide DNA demethylation is a unique feature of mammalian development and naïve pluripotent stem cells. Here, we describe a recently evolved pathway in which global hypomethylation is achieved by the coupling of active and passive demethylation. TET activity is required, albeit indirectly, for global demethylation, which mostly occurs at sites devoid of TET binding. Instead, TET-mediated active demethylation is locus-specific and necessary for activating a subset of genes, including the naïve pluripotency and germline marker Dppa3 (Stella, Pgc7). DPPA3 in turn drives large-scale passive demethylation by directly binding and displacing UHRF1 from chromatin, thereby inhibiting maintenance DNA methylation. Although unique to mammals, we show that DPPA3 alone is capable of inducing global DNA demethylation in non-mammalian species (Xenopus and medaka) despite their evolutionary divergence from mammals more than 300 million years ago. Our findings suggest that the evolution of Dppa3 facilitated the emergence of global DNA demethylation in mammals.

A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona.

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

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.

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.

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.

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.

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.

Assessing Cell-to-Cell DNA Methylation Variability on Individual Long Reads.

Understanding cell-to-cell variability in cytosine methylation is essential for understanding cellular perturbation and its molecular machinery. However, conventional methylation studies have focused on the differences in the average levels between cell types while overlooking methylation heterogeneity within cell types. Little information has been uncovered using recent single-cell methods because of either technical limitations or the great labor required to process many single cells. Here, we report the highly efficient detection of cell-to-cell DNA methylation variability in liver tissue, based on comparing the methylation status of adjacent CpG sites on long sequencing reads. This method provides abundant methylation linkage information and enables genome-wide estimation of cell-to-cell variability. We observed repressed methylation variability in hypomethylated regions compared with the variability in hypomethylated regions across the genome, which we confirmed using public human sperm data. A gradual change in methylation status at the boundaries of hypomethylated regions was observed for the first time. This approach allows the concise, comprehensive assessment of cell-to-cell DNA methylation variability.

Associations between nucleosome phasing, sequence asymmetry, and tissue-specific expression in a set of inbred Medaka species.

BACKGROUND: Transcription start sites (TSSs) with pronounced and phased nucleosome arrays downstream and nucleosome-depleted regions upstream of TSSs are observed in various species. RESULTS: We have characterized sequence variation and expression properties of this set of TSSs (which we call "Nucleocyclic TSSs") using germline and somatic cells of three medaka (Oryzias latipes) inbred isolates from different locations. We found nucleocyclic TSSs in medaka to be associated with higher gene expression and characterized by a clear boundary in sequence composition with potentially-nucleosome-destabilizing A/T-enrichment upstream (p < 10(-60)) and nucleosome- accommodating C/G-enrichment downstream (p < 10(-40)) that was highly conserved from an ancestor. A substantial genetic distance between the strains facilitated the in-depth analysis of patterns of fixed mutations, revealing a localization-specific equilibrium between the rates of distinct mutation categories that would serve to maintain the conserved sequence anisotropy around TSSs. Downstream of nucleocyclic TSSs, C to T, T to C, and other mutation rates on the sense strand increased around first nucleosome dyads and decreased around first linkers, which contrasted with genomewide mutational patterns around nucleosomes (p < 5 %). C to T rates are higher than G to A rates around nucleosome associated with germline nucleocyclic TSS sites (p < 5 %), potentially due to the asymmetric effect of transcription-coupled repair. CONCLUSIONS: Our results demonstrate an atypical evolutionary process surrounding nucleocyclic TSSs.

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

DNA methylation is a fundamental epigenetic modification in vertebrate genomes and a small fraction of genomic regions is hypomethylated. Previous studies have implicated hypomethylated regions in gene regulation, but their functions in vertebrate development remain elusive. To address this issue, we generated epigenomic profiles that include base-resolution DNA methylomes and histone modification maps from both pluripotent cells and mature organs of medaka fish and compared the profiles with those of human ES cells. We found that a subset of hypomethylated domains harbor H3K27me3 (K27HMDs) and their size positively correlates with the accumulation of H3K27me3. Large K27HMDs are conserved between medaka and human pluripotent cells and predominantly contain promoters of developmental transcription factor genes. These key genes were found to be under strong transcriptional repression, when compared with other developmental genes with smaller K27HMDs. Furthermore, human-specific K27HMDs show an enrichment of neuronal activity-related genes, which suggests a distinct regulation of these genes in medaka and human. In mature organs, some of the large HMDs become shortened by elevated DNA methylation and associate with sustained gene expression. This study highlights the significance of domain size in epigenetic gene regulation. We propose that large K27HMDs play a crucial role in pluripotent cells by strictly repressing key developmental genes, whereas their shortening consolidates long-term gene expression in adult differentiated cells.

[Multivariate quantification of acupuncture and moxibustion treatment in combination with questionnaire and salivary stress markers].

Pain and stress alleviation after acupuncture treatment was assessed in this study. Patients responded to a questionnaire designed to determine the amount of stress they were experiencing, and data were obtained for patient salivary amylase, cortisol, secretary IgA (s-IgA), and leptin receptor (OBRb). As a part of this study on acute pain, 6 factors were extracted from the questionnaire. The second factor (pain removal) was well correlated with salivary amylase activity in patients with cervico-omo-brachial syndrome. An evaluation of cumulative acupuncture treatments showed that salivary cortisol increased and s-IgA decreased. In addition, a decreased s-IgA level significantly correlated with chronic pain removal. The questionnaire correlated well with measurements of salivary markers suggesting that they can be taken as indices of therapeutic efficacy in acupuncture treatment.