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.

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.

Bcl11b/Ctip2 controls the differentiation of vomeronasal sensory neurons in mice.

The transcription factor Bcl11b/Ctip2 plays critical roles in the development of several systems and organs, including the immune system, CNS, skin, and teeth. Here, we show that Bcl11b/Ctip2 is highly expressed in the developing vomeronasal system in mice and is required for its proper development. Bcl11b/Ctip2 is expressed in postmitotic vomeronasal sensory neurons (VSNs) in the vomeronasal epithelium (VNE) as well as projection neurons and GABAergic interneurons in the accessory olfactory bulb (AOB). In the absence of Bcl11b, these neurons are born in the correct number, but VSNs selectively die by apoptosis. The critical role of Bcl11b in vomeronasal system development is demonstrated by the abnormal phenotypes of Bcl11b-deficient mice: disorganization of layer formation of the AOB, impaired axonal projections of VSNs, a significant reduction in the expression of vomeronasal receptor genes, and defective mature differentiation of VSNs. VSNs can be classified into two major types of neurons, vomeronasal 1 receptor (V1r)/Gα(i2)-positive and vomeronasal 2 receptor (V2r)/Gα(o)-positive VSNs. We found that all Gα(i2)-positive cells coexpressed Gα(o) during embryogenesis. This coexpression is also observed in newly differentiated neurons in the adult VNE. Interestingly, loss of Bcl11b function resulted in an increased number of V1r/Gα(i2)-type VSNs and a decreased number of V2r/Gα(o)-type VSNs, suggesting that Bcl11b regulates the fate choice between these two VSN types. These results indicate that Bcl11b/Ctip2 is an essential regulator of the differentiation and dichotomy of VSNs.