The ancestor of extant Japanese fancy mice contributed to the mosaic genomes of classical inbred strains

The 36th International Mammalian Genome Conference: A scientific gathering under the cherry blossoms in Tsukuba

The 36th International Mammalian Genome Conference (IMGC) was held in a hybrid format at the Tsukuba International Congress Center in Tsukuba, Ibaraki, Japan, for 4 days from March 28 to 31, 2023. This international conference on functional genomics of mouse, human, and other mammalian species attracted 246 participants in total, of which 129 were from outside Japan, including Europe, the United States and Asia, and 117 participants were from Japan. The conference included three technical workshops, keynote lectures by domestic researchers, commemorative lectures for the conference awards, 57 oral presentations, and 97 poster presentations. The event was a great success. Topics included the establishment and analysis of disease models using genetically engineered or spontaneous mutant mice, systems genetic analysis using mouse strains such as wild-derived mice and recombinant inbred mouse strains, infectious diseases, immunology, and epigenetics. In addition, as a joint program, a two-day RIKEN Symposium was held, and active discussions continued over the four-day period. Also, there was a trainee symposium, in which young researchers were encouraged to participate, and excellent papers were selected as oral presentations in the main session.

Insertion of short L1 sequences generates inter-strain histone acetylation differences in the mouse

BACKGROUND

Gene expression divergence between populations and between individuals can emerge from genetic variations within the genes and/or in the cis regulatory elements. Since epigenetic modifications regulate gene expression, it is conceivable that epigenetic variations in cis regulatory elements can also be a source of gene expression divergence.

RESULTS

In this study, we compared histone acetylation (namely, H3K9ac) profiles in two mouse strains of different subspecies origin, C57BL/6 J (B6) and MSM/Ms (MSM), as well as their F1 hybrids. This identified 319 regions of strain-specific acetylation, about half of which were observed between the alleles of F1 hybrids. While the allele-specific presence of the interferon regulatory factor 3 (IRF3) binding sequence was associated with allele-specific histone acetylation, we also revealed that B6-specific insertions of a short 3' fragment of LINE-1 (L1) retrotransposon occur within or proximal to MSM-specific acetylated regions. Furthermore, even in hyperacetylated domains, flanking regions of non-polymorphic 3' L1 fragments were hypoacetylated, suggesting a general activity of the 3' L1 fragment to induce hypoacetylation. Indeed, we confirmed the binding of the 3' region of L1 by three Krüppel-associated box domain-containing zinc finger proteins (KZFPs), which interact with histone deacetylases. These results suggest that even a short insertion of L1 would be excluded from gene- and acetylation-rich regions by natural selection. Finally, mRNA-seq analysis for F1 hybrids was carried out, which disclosed a link between allele-specific promoter/enhancer acetylation and gene expression.

CONCLUSIONS

This study disclosed a number of genetic changes that have changed the histone acetylation levels during the evolution of mouse subspecies, a part of which is associated with gene expression changes. Insertions of even a very short L1 fragment can decrease the acetylation level in their neighboring regions and thereby have been counter-selected in gene-rich regions, which may explain a long-standing mystery of discrete genomic distribution of LINEs and SINEs.

Type specimens, taxonomic history, and genetic analysis of the Japanese dancing mouse or waltzer, Muswagneri variety rotans Droogleever Fortuyn, 1912 (Mammalia, Muridae)

In the present paper, the existence and location of the type series of the Japanese dancing mouse or waltzer, Muswagneri variety rotans Droogleever Fortuyn, 1912, are established, and a lectotype is designated. Available type specimens are measured, and some morphological parameters, sex, and general condition of the specimens are recorded. A literature survey was conducted, and an attempt is made to clarify the position of M.wagneri variety rotans in the taxonomy of Mus. A genetic analysis suggests that the type series of the Japanese dancing mouse represent a crossbred, or derivation of a crossbred, between the original Japanese dancing mouse of Musmusculusmolossinus Temminck 1844 origin and European fancy or laboratory mice of Musmusculusdomesticus Schwarz & Schwarz, 1943 origin. Much of their genome was replaced and occupied by Musmusculusdomesticus type genome, probably through extensive breeding with European mice.

Inference of selective forces on house mouse genomes during secondary contact in East Asia

The house mouse (Mus musculus), which is commensal to humans, has spread globally via human activities, leading to secondary contact between genetically divergent subspecies. This pattern of genetic admixture can provide insights into the selective forces at play in this well-studied model organism. Our analysis of 163 house mouse genomes, with a particular focus on East Asia, revealed substantial admixture between the subspecies castaneus and musculus, particularly in Japan and southern China. We revealed, despite the different level of autosomal admixture among regions, that all Y Chromosomes in the East Asian samples belonged to the musculus-type haplogroup, potentially explained by genomic conflict under sex-ratio distortion owing to varying copy numbers of ampliconic genes on sex chromosomes, Slx and Sly Our computer simulations, designed to replicate the observed scenario, show that the preferential fixation of musculus-type Y Chromosomes can be achieved with a slight increase in the male-to-female birth ratio. We also investigated the influence of selection on the posthybridization of the subspecies castaneus and musculus in Japan. Even though the genetic background of most Japanese samples closely resembles the subspecies musculus, certain genomic regions overrepresented the castaneus-like genetic components, particularly in immune-related genes. Furthermore, a large genomic block (∼2 Mbp) containing a vomeronasal/olfactory receptor gene cluster predominantly harbored castaneus-type haplotypes in the Japanese samples, highlighting the crucial role of olfaction-based recognition in shaping hybrid genomes.

Combined and differential roles of ADD domains of DNMT3A and DNMT3L on DNA methylation landscapes in mouse germ cells

DNA methyltransferase 3A (DNMT3A) and its catalytically inactive cofactor DNA methyltransferase 3-Like (DNMT3L) proteins form functional heterotetramers to deposit DNA methylation in mammalian germ cells. While both proteins have an ATRX-DNMT3-DNMT3L (ADD) domain that recognizes histone H3 tail unmethylated at lysine-4 (H3K4me0), the combined and differential roles of the domains in the two proteins have not been fully defined in vivo. Here we investigate DNA methylation landscapes in female and male germ cells derived from mice with loss-of-function amino acid substitutions in the ADD domains of DNMT3A and/or DNMT3L. Mutations in either the DNMT3A-ADD or the DNMT3L-ADD domain moderately decrease global CG methylation levels, but to different degrees, in both germ cells. Furthermore, when the ADD domains of both DNMT3A and DNMT3L lose their functions, the CG methylation levels are much more reduced, especially in oocytes, comparable to the impact of the Dnmt3a/3L knockout. In contrast, aberrant accumulation of non-CG methylation occurs at thousands of genomic regions in the double mutant oocytes and spermatozoa. These results highlight the critical role of the ADD-H3K4me0 binding in proper CG and non-CG methylation in germ cells and the various impacts of the ADD domains of the two proteins.

The Expansion of House Mouse Major Urinary Protein Genes Likely Did Not Facilitate Commensalism with Humans

Mouse wild-derived strains (WDSs) combine the advantages of classical laboratory stocks and wild animals, and thus appear to be promising tools for diverse biomedical and evolutionary studies. We employed 18 WDSs representing three non-synanthropic species (Mus spretus, Mus spicilegus, and M. macedonicus) and three house mouse subspecies (Mus musculus musculus, M. m. domesticus, M. m. castaneus), which are all important human commensals to explore whether the number of major urinary protein (MUP) genes and their final protein levels in urine are correlated with the level of commensalism. Contrary to expectations, the MUP copy number (CN) and protein excretion in the strains derived from M. m. castaneus, which is supposed to be the strongest commensal, were not significantly different from the non-commensal species. Regardless of an overall tendency for higher MUP amounts in taxa with a higher CN, there was no significant correlation at the strain level. Our study thus suggests that expansion of the Mup cluster, which appeared before the house mouse diversification, is unlikely to facilitate commensalism with humans in three house mouse subspecies. Finally, we found considerable variation among con(sub)specific WDSs, warning against generalisations of results based on a few strains.

Genealogical characterization of regional populations and dorsal coat color variation in the house mouse Mus musculus from Asia based on haplotype structure analysis of a gene-rich region harboring Mc1r

We analyzed 196 haplotype sequences from a gene-rich region (250 kb) that includes Mc1r, a gene involved in coat color regulation, to gain insight into the evolution of coat color variation in subspecies of the house mouse Mus musculus. Phylogenetic networks revealed haplotype groups from the major subspecies of M. m. castaneus (CAS), M. m. domesticus (DOM), and M. m. musculus (MUS). Using haplotype sequences assigned to each of CAS and MUS through phylogenetic analysis, we proposed migration routes associated with prehistoric humans from west to east across Eurasia. Comparing nucleotide diversity among subspecies-specific haplotypes in different geographic areas showed a marked reduction during migration, particularly in MUS-derived haplotypes from Korea and Japan, suggesting intensive population bottlenecks during migration. We found that a C>T polymorphism at site 302 (c.302C>T) in the Mc1r coding region correlated with a darkening of dorsal fur color in both CAS and MUS. However, C/C homozygous mice in MUS showed marked variation in lightness, indicating the possibility of another genetic determinant that affects the lightness of dorsal fur color. Detailed sequence comparisons of haplotypes revealed that short fragments assigned to DOM were embedded in CAS-assigned fragments, indicating ancient introgression between subspecies. The estimated age of c.302C>T also supports the hypothesis that genetic interaction between subspecies occurred in ancient times. This suggests that the genome of M. musculus evolved through gene flow between subspecies over an extended period before the movement of the species in conjunction with prehistoric humans.

The DNMT3A ADD domain is required for efficient de novo DNA methylation and maternal imprinting in mouse oocytes

Establishment of a proper DNA methylation landscape in mammalian oocytes is important for maternal imprinting and embryonic development. De novo DNA methylation in oocytes is mediated by the DNA methyltransferase DNMT3A, which has an ATRX-DNMT3-DNMT3L (ADD) domain that interacts with histone H3 tail unmethylated at lysine-4 (H3K4me0). The domain normally blocks the methyltransferase domain via intramolecular interaction and binding to histone H3K4me0 releases the autoinhibition. However, H3K4me0 is widespread in chromatin and the role of the ADD-histone interaction has not been studied in vivo. We herein show that amino-acid substitutions in the ADD domain of mouse DNMT3A cause dwarfism. Oocytes derived from homozygous females show mosaic loss of CG methylation and almost complete loss of non-CG methylation. Embryos derived from such oocytes die in mid-to-late gestation, with stochastic and often all-or-none-type CG-methylation loss at imprinting control regions and misexpression of the linked genes. The stochastic loss is a two-step process, with loss occurring in cleavage-stage embryos and regaining occurring after implantation. These results highlight an important role for the ADD domain in efficient, and likely processive, de novo CG methylation and pose a model for stochastic inheritance of epigenetic perturbations in germ cells to the next generation.

Expression analysis of box C/D snoRNAs with SNPs between C57BL/6 and MSM/Ms strains in male mouse

MSM/Ms mouse derived from the Japanese wild mouse has unique characteristics compared to the widely used C57BL/6 mouse. To examine the usefulness of the MSM/Ms mouse for the comparative genomic analysis, expression of small RNAs were analyzed by the large-scale sequence analysis for two strains of mouse, C57BL/6 and MSM/Ms. As a trial, expression of box C/D snoRNAs, which are the most abundant small RNAs in the cell, were analyzed. By the comparison of the read number for each fragment, 11 snoRNAs with single nucleotide polymorphisms (SNPs) were detected. One of the snoRNAs, SNORD53, shows the expression only for MSM/Ms and this snoRNA has a mutation in the box sequence in C57BL/6. Thus, it was demonstrated that the proposed experimental system using SNPs can give new insight for the gene expression regulation.

Conservation and divergence of canonical and non-canonical imprinting in murids

BACKGROUND

Genomic imprinting affects gene expression in a parent-of-origin manner and has a profound impact on complex traits including growth and behavior. While the rat is widely used to model human pathophysiology, few imprinted genes have been identified in this murid. To systematically identify imprinted genes and genomic imprints in the rat, we use low input methods for genome-wide analyses of gene expression and DNA methylation to profile embryonic and extraembryonic tissues at allele-specific resolution.

RESULTS

We identify 14 and 26 imprinted genes in these tissues, respectively, with 10 of these genes imprinted in both tissues. Comparative analyses with mouse reveal that orthologous imprinted gene expression and associated canonical DNA methylation imprints are conserved in the embryo proper of the Muridae family. However, only 3 paternally expressed imprinted genes are conserved in the extraembryonic tissue of murids, all of which are associated with non-canonical H3K27me3 imprints. The discovery of 8 novel non-canonical imprinted genes unique to the rat is consistent with more rapid evolution of extraembryonic imprinting. Meta-analysis of novel imprinted genes reveals multiple mechanisms by which species-specific imprinted expression may be established, including H3K27me3 deposition in the oocyte, the appearance of ZFP57 binding motifs, and the insertion of endogenous retroviral promoters.

CONCLUSIONS

In summary, we provide an expanded list of imprinted loci in the rat, reveal the extent of conservation of imprinted gene expression, and identify potential mechanisms responsible for the evolution of species-specific imprinting.

Tissue-specific Grb10/Ddc insulator drives allelic architecture for cardiac development

Allele-specific expression of imprinted gene clusters is governed by gametic DNA methylation at master regulators called imprinting control regions (ICRs). Non-gametic or secondary differentially methylated regions (DMRs) at promoters and exonic regions reinforce monoallelic expression but do not control an entire cluster. Here, we unveil an unconventional secondary DMR that is indispensable for tissue-specific imprinting of two previously unlinked genes, Grb10 and Ddc. Using polymorphic mice, we mapped an intronic secondary DMR at Grb10 with paternal-specific CTCF binding (CBR2.3) that forms contacts with Ddc. Deletion of paternal CBR2.3 removed a critical insulator, resulting in substantial shifting of chromatin looping and ectopic enhancer-promoter contacts. Destabilized gene architecture precipitated abnormal Grb10-Ddc expression with developmental consequences in the heart and muscle. Thus, we redefine the Grb10-Ddc imprinting domain by uncovering an unconventional intronic secondary DMR that functions as an insulator to instruct the tissue-specific, monoallelic expression of multiple genes-a feature previously ICR exclusive.

CMIC: predicting DNA methylation inheritance of CpG islands with embedding vectors of variable-length k-mers

Background

Epigenetic modifications established in mammalian gametes are largely reprogrammed during early development, however, are partly inherited by the embryo to support its development. In this study, we examine CpG island (CGI) sequences to predict whether a mouse blastocyst CGI inherits oocyte-derived DNA methylation from the maternal genome. Recurrent neural networks (RNNs), including that based on gated recurrent units (GRUs), have recently been employed for variable-length inputs in classification and regression analyses. One advantage of this strategy is the ability of RNNs to automatically learn latent features embedded in inputs by learning their model parameters. However, the available CGI dataset applied for the prediction of oocyte-derived DNA methylation inheritance are not large enough to train the neural networks.

Results

We propose a GRU-based model called CMIC (CGI Methylation Inheritance Classifier) to augment CGI sequence by converting it into variable-length k-mers, where the length k is randomly selected from the range \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{\min }$$\end{document}kmin to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_{\max }$$\end{document}kmax, N times, which were then used as neural network input. N was set to 1000 in the default setting. In addition, we proposed a new embedding vector generator for k-mers called splitDNA2vec. The randomness of this procedure was higher than the previous work, dna2vec.

Conclusions

We found that CMIC can predict the inheritance of oocyte-derived DNA methylation at CGIs in the maternal genome of blastocysts with a high F-measure (0.93). We also show that the F-measure can be improved by increasing the parameter N, that is, the number of sequences of variable-length k-mers derived from a single CGI sequence. This implies the effectiveness of augmenting input data by converting a DNA sequence to N sequences of variable-length k-mers. This approach can be applied to different DNA sequence classification and regression analyses, particularly those involving a small amount of data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04916-3.

Two Loci Contribute to Age-Related Hearing Loss Resistance in the Japanese Wild-Derived Inbred MSM/Ms Mice

An MSM/Ms strain was established using Japanese wild mice, which exhibit resistance to several phenotypes associated with aging, such as obesity, inflammation, and tumorigenesis, compared to common inbred mouse strains. MSM/Ms strain is resistant to age-related hearing loss, and their auditory abilities are sustained for long durations. The age-related hearing loss 3 (ahl3) locus contributes to age-related hearing in MSM/Ms strain. We generated ahl3 congenic strains by transferring a genomic region on chromosome 17 from MSM/Ms mice into C57BL/6J mice. Although C57BL/6J mice develop age-related hearing loss because of the ahl allele of the cadherin 23 gene, the development of middle- to high-frequency hearing loss was significantly delayed in an ahl3 congenic strain. Moreover, the novel age-related hearing loss 10 (ahl10) locus associated with age-related hearing resistance in MSM/Ms strain was mapped to chromosome 12. Although the resistance effects in ahl10 congenic strain were slightly weaker than those in ahl3 congenic strain, slow progression of age-related hearing loss was confirmed in ahl10 congenic strain despite harboring the ahl allele of cadherin 23. These results suggest that causative genes and polymorphisms of the ahl3 and ahl10 loci are important targets for the prevention and treatment of age-related hearing loss.

Histone H3K36me2 and H3K36me3 form a chromatin platform essential for DNMT3A-dependent DNA methylation in mouse oocytes

Establishment of the DNA methylation landscape of mammalian oocytes, mediated by the DNMT3A-DNMT3L complex, is crucial for reproduction and development. In mouse oocytes, high levels of DNA methylation occur exclusively in the transcriptionally active regions, with moderate to low levels of methylation in other regions. Histone H3K36me3 mediates the high levels of methylation in the transcribed regions; however, it is unknown which histone mark guides the methylation in the other regions. Here, we show that, in mouse oocytes, H3K36me2 is highly enriched in the X chromosome and is broadly distributed across all autosomes. Upon H3K36me2 depletion, DNA methylation in moderately methylated regions is selectively affected, and a methylation pattern unique to the X chromosome is switched to an autosome-like pattern. Furthermore, we find that simultaneous depletion of H3K36me2 and H3K36me3 results in global hypomethylation, comparable to that of DNMT3A depletion. Therefore, the two histone marks jointly provide the chromatin platform essential for guiding DNMT3A-dependent DNA methylation in mouse oocytes.

Activation of Xist by an evolutionarily conserved function of KDM5C demethylase

XX female and XY male therian mammals equalize X-linked gene expression through the mitotically-stable transcriptional inactivation of one of the two X chromosomes in female somatic cells. Here, we describe an essential function of the X-linked homolog of an ancestral X-Y gene pair, Kdm5c-Kdm5d, in the expression of Xist lncRNA, which is required for stable X-inactivation. Ablation of Kdm5c function in females results in a significant reduction in Xist RNA expression. Kdm5c encodes a demethylase that enhances Xist expression by converting histone H3K4me2/3 modifications into H3K4me1. Ectopic expression of mouse and human KDM5C, but not the Y-linked homolog KDM5D, induces Xist in male mouse embryonic stem cells (mESCs). Similarly, marsupial (opossum) Kdm5c but not Kdm5d also upregulates Xist in male mESCs, despite marsupials lacking Xist, suggesting that the KDM5C function that activates Xist in eutherians is strongly conserved and predates the divergence of eutherian and metatherian mammals. In support, prototherian (platypus) Kdm5c also induces Xist in male mESCs. Together, our data suggest that eutherian mammals co-opted the ancestral demethylase KDM5C during sex chromosome evolution to upregulate Xist for the female-specific induction of X-inactivation.

Insights into Mus musculus population structure across Eurasia revealed by whole-genome analysis

For more than 100 years, house mice (Mus musculus) have been used as a key animal model in biomedical research. House mice are genetically diverse, yet their genetic background at the global level has not been fully understood. Previous studies have suggested that they originated in South Asia and diverged into three major subspecies, almost simultaneously, approximately 110,000–500,000 years ago; however, they have spread across the world with the migration of modern humans in prehistoric and historic times (∼10,000 years ago to the present day) and have undergone secondary contact, which has complicated the genetic landscape of wild house mice. In this study, we sequenced the whole-genome sequences of 98 wild house mice collected from Eurasia, particularly East Asia, Southeast Asia, and South Asia. Although wild house mice were found to consist of three major genetic groups corresponding to the three major subspecies, individuals representing admixtures between subspecies were more prevalent in East Asia than has been previously recognized. Furthermore, several samples exhibited an incongruent pattern of genealogies between mitochondrial and autosomal genomes. Using samples that likely retained the original genetic components of subspecies with the least admixture, we estimated the pattern and timing of divergence among the subspecies. The estimated divergence time of the three subspecies was 187,000–226,000 years ago. These results will help us to understand the genetic diversity of wild mice on a global scale, and the findings will be particularly useful in future biomedical and evolutionary studies involving laboratory mice established from such wild mice.

Functional polymorphism in Pak1-3'UTR alters skin tumor susceptibility by alternative polyadenylation

We identified a functional single nucleotide polymorphism (SNP) in the 3' untranslated region (UTR) of p21-activated kinase 1 (Pak1) that is responsible for the Skin tumor modifier of MSM 1a locus. Candidate SNPs in the 3'UTR of Pak1 from resistance strain MSM/Ms were introduced into susceptible strain FVB/N using CRISPR/Cas9. DMBA/TPA skin carcinogenesis experiments revealed an SNP (Pak1-3'UTR-6^C>T: rs31627325) that strongly suppressed skin tumors. Furthermore, Muscleblind-Like Splicing Regulator 1 bound more strongly to FVB-allele (6^C/C) and regulated the transcript length in the 3'UTR of Pak1 and tumorigenesis via polyadenylation. Therefore, the alternative polyadenylation of Pak1 is cis-regulated by rs31627325.

MoG+: a database of genomic variations across three mouse subspecies for biomedical research

Laboratory mouse strains have mosaic genomes derived from at least three major subspecies that are distributed in Eurasia. Here, we describe genomic variations in ten inbred strains: Mus musculus musculus-derived BLG2/Ms, NJL/Ms, CHD/Ms, SWN/Ms, and KJR/Ms; M. m. domesticus-derived PGN2/Ms and BFM/Ms; M. m. castaneus-derived HMI/Ms; and JF1/Ms and MSM/Ms, which were derived from a hybrid between M. m. musculus and M. m. castaneus. These strains were established by Prof. Moriwaki in the 1980s and are collectively named the “Mishima Battery”. These strains show large phenotypic variations in body size and in many physiological traits. We resequenced the genomes of the Mishima Battery strains and performed a comparative genomic analysis with dbSNP data. More than 81 million nucleotide coordinates were identified as variant sites due to the large genetic distances among the mouse subspecies; 8,062,070 new SNP sites were detected in this study, and these may underlie the large phenotypic diversity observed in the Mishima Battery. The new information was collected in a reconstructed genome database, termed MoG+ that includes new application software and viewers. MoG+ intuitively visualizes nucleotide variants in genes and intergenic regions, and amino acid substitutions across the three mouse subspecies. We report statistical data from the resequencing and comparative genomic analyses and newly collected phenotype data of the Mishima Battery, and provide a brief description of the functions of MoG+, which provides a searchable and unique data resource of the numerous genomic variations across the three mouse subspecies. The data in MoG+ will be invaluable for research into phenotype-genotype links in diverse mouse strains.

Mouse resources at the RIKEN BioResource Research Center and the National BioResource Project core facility in Japan

The RIKEN BioResource Research Center (BRC) was established in 2001 as a comprehensive biological resource center in Japan. The Experimental Animal Division, one of the BRC infrastructure divisions, has been designated as the core facility for mouse resources within the National BioResource Project (NBRP) by the Japanese government since FY2002. Our activities regarding the collection, preservation, quality control, and distribution of mouse resources have been supported by the research community, including evaluations and guidance on advancing social and research needs, as well as the operations and future direction of the BRC. Expenditure for collection, preservation, and quality-control operations of the BRC, as a national core facility, has been funded by the government, while distribution has been separately funded by users' reimbursement fees. We have collected over 9000 strains created mainly by Japanese scientists including Nobel laureates and researchers in cutting-edge fields and distributed mice to 7000 scientists with 1500 organizations in Japan and globally. Our users have published 1000 outstanding papers and a few dozen patents. The collected mouse resources are accessible via the RIKEN BRC website, with a revised version of the searchable online catalog. In addition, to enhance the visibility of useful strains, we have launched web corners designated as the "Mouse of the Month" and "Today's Tool and Model." Only high-demand strains are maintained in live colonies, while other strains are cryopreserved as embryos or sperm to achieve cost-effective management. Since 2007, the RIKEN BRC has built up a back-up facility in the RIKEN Harima branch to protect the deposited strains from disasters. Our mice have been distributed with high quality through the application of strict microbial and genetic quality control programs that cover a globally accepted pathogens list and mutated alleles generated by various methods. Added value features, such as information about users' publications, standardized phenotyping data, and genome sequences of the collected strains, are important to facilitate the use of our resources. We have added and disseminated such information in collaboration with the NBRP Information Center and the NBRP Genome Information Upgrading Program. The RIKEN BRC has participated in international mouse resource networks such as the International Mouse Strain Resource, International Mouse Phenotyping Consortium, and Asian Mouse Mutagenesis and Resource Association to facilitate the worldwide use of high-quality mouse resources, and as a consequence it contributes to reproducible life science studies and innovation around the globe.

Gene trapping reveals a new transcriptionally active genome element: The chromosome-specific clustered trap region

Nearly half of the human genome consists of repetitive sequences such as long interspersed nuclear elements. The relationship between these repeating sequences and diseases has remained unclear. Gene trapping is a useful technique for disrupting a gene and expressing a reporter gene by using the promoter activity of the gene. The analysis of trapped genes revealed a new genome element-the chromosome-specific clustered trap (CSCT) region. For any examined sequence within this region, an equivalent was found using the BLAT of the University of California, Santa Cruz (UCSC) Genome Browser. CSCT13 mapped to chromosome 13 and contained only three genes. To elucidate its in vivo function, the whole CSCT13 region (1.6 Mbp) was deleted using the CRISPR/Cas9 system in mouse embryonic stem cells, and subsequently, a CSCT13 knockout mouse line was established. The rate of homozygotes was significantly lower than expected according to Mendel's laws. In addition, the number of offspring obtained by mating homozygotes was significantly smaller than that obtained by crossing controls. Furthermore, CSCT13 might have an effect on meiotic homologous recombination. This study identifies a transcriptionally active CSCT with an important role in mouse development.

B2 SINE Copies Serve as a Transposable Boundary of DNA Methylation and Histone Modifications in the Mouse

More than one million copies of short interspersed elements (SINEs), a class of retrotransposons, are present in the mammalian genomes, particularly within gene-rich genomic regions. Evidence has accumulated that ancient SINE sequences have acquired new binding sites for transcription factors (TFs) through multiple mutations following retrotransposition, and as a result have rewired the host regulatory network during the course of evolution. However, it remains unclear whether currently active SINEs contribute to the expansion of TF binding sites. To study the mobility, expression, and function of SINE copies, we first identified about 2,000 insertional polymorphisms of SINE B1 and B2 families within Mus musculus. Using a novel RNA sequencing method designated as melRNA-seq, we detected the expression of SINEs in male germ cells at both the subfamily and genomic copy levels: the vast majority of B1 RNAs originated from evolutionarily young subfamilies, whereas B2 RNAs originated from both young and old subfamilies. DNA methylation and chromatin immunoprecipitation-sequencing (ChIP-seq) analyses in liver revealed that polymorphic B2 insertions served as a boundary element inhibiting the expansion of DNA hypomethylated and histone hyperacetylated regions, and decreased the expression of neighboring genes. Moreover, genomic B2 copies were enriched at the boundary of various histone modifications, and chromatin insulator protein, CCCTC-binding factor, a well-known chromatin boundary protein, bound to >100 polymorphic and >10,000 non-polymorphic B2 insertions. These results suggest that the currently active B2 copies are mobile boundary elements that can modulate chromatin modifications and gene expression, and are likely involved in epigenomic and phenotypic diversification of the mouse species.

Zfp57 inactivation illustrates the role of ICR methylation in imprinted gene expression during neural differentiation of mouse ESCs

ZFP57 is required to maintain the germline-marked differential methylation at imprinting control regions (ICRs) in mouse embryonic stem cells (ESCs). Although DNA methylation has a key role in genomic imprinting, several imprinted genes are controlled by different mechanisms, and a comprehensive study of the relationship between DMR methylation and imprinted gene expression is lacking. To address the latter issue, we differentiated wild-type and Zfp57^-/- hybrid mouse ESCs into neural precursor cells (NPCs) and evaluated allelic expression of imprinted genes. In mutant NPCs, we observed a reduction of allelic bias of all the 32 genes that were imprinted in wild-type cells, demonstrating that ZFP57-dependent methylation is required for maintaining or acquiring imprinted gene expression during differentiation. Analysis of expression levels showed that imprinted genes expressed from the non-methylated chromosome were generally up-regulated, and those expressed from the methylated chromosome were down-regulated in mutant cells. However, expression levels of several imprinted genes acquiring biallelic expression were not affected, suggesting the existence of compensatory mechanisms that control their RNA level. Since neural differentiation was partially impaired in Zfp57-mutant cells, this study also indicates that imprinted genes and/or non-imprinted ZFP57-target genes are required for proper neurogenesis in cultured ESCs.

A wild-derived inbred mouse strain, MSM/Ms, provides insights into novel skin tumor susceptibility genes

Cancer is one of the most catastrophic human genetic diseases. Experimental animal cancer models are essential for gaining insights into the complex interactions of different cells and genes in tumor initiation, promotion, and progression. Mouse models have been extensively used to analyze the genetic basis of cancer susceptibility. They have led to the identification of multiple loci that confer, either alone or in specific combinations, an increased susceptibility to cancer, some of which have direct translatability to human cancer. Additionally, wild-derived inbred mouse strains are an advantageous reservoir of novel genetic polymorphisms of cancer susceptibility genes, because of the evolutionary divergence between wild and classical inbred strains. Here, we review mapped Stmm (skintumor modifier of MSM) loci using a Japanese wild-derived inbred mouse strain, MSM/Ms, and describe recent advances in our knowledge of the genes responsible for Stmm loci in the 7,12-dimethylbenz(a)anthracene (DMBA)/12-O-tetradecanoylphorbol-13-acetate (TPA) two-stage skin carcinogenesis model.

The Japanese Wild-Derived Inbred Mouse Strain, MSM/Ms in Cancer Research

MSM/Ms is a unique inbred mouse strain derived from the Japanese wild mouse, Mus musculus molossinus, which has been approximately 1 million years genetically distant from standard inbred mouse strains mainly derived from M. m. domesticus. Due to its genetic divergence, MSM/Ms has been broadly used in linkage studies. A bacterial artificial chromosome (BAC) library was constructed for the MSM/Ms genome, and sequence analysis of the MSM/Ms genome showed approximately 1% of nucleotides differed from those in the commonly used inbred mouse strain, C57BL/6J. Therefore, MSM/Ms mice are thought to be useful for functional genome studies. MSM/Ms mice show unique characteristics of phenotypes, including its smaller body size, resistance to high-fat-diet-induced diabetes, high locomotive activity, and resistance to age-onset hearing loss, inflammation, and tumorigenesis, which are distinct from those of common inbred mouse strains. Furthermore, ES (Embryonic Stem) cell lines established from MSM/Ms allow the MSM/Ms genome to be genetically manipulated. Therefore, genomic and phenotypic analyses of MSM/Ms reveal novel insights into gene functions that were previously not obtained from research on common laboratory strains. Tumorigenesis-related MSM/Ms-specific genetic traits have been intensively investigated in Japan. Furthermore, radiation-induced thymic lymphomas and chemically-induced skin tumors have been extensively examined using MSM/Ms.

Establishment and application of information resource of mutant mice in RIKEN BioResource Research Center

Online databases are crucial infrastructures to facilitate the wide effective and efficient use of mouse mutant resources in life sciences. The number and types of mouse resources have been rapidly growing due to the development of genetic modification technology with associated information of genomic sequence and phenotypes. Therefore, data integration technologies to improve the findability, accessibility, interoperability, and reusability of mouse strain data becomes essential for mouse strain repositories. In 2020, the RIKEN BioResource Research Center released an integrated database of bioresources including, experimental mouse strains, Arabidopsis thaliana as a laboratory plant, cell lines, microorganisms, and genetic materials using Resource Description Framework-related technologies. The integrated database shows multiple advanced features for the dissemination of bioresource information. The current version of our online catalog of mouse strains which functions as a part of the integrated database of bioresources is available from search bars on the page of the Center (https://brc.riken.jp) and the Experimental Animal Division (https://mus.brc.riken.jp/) websites. The BioResource Research Center also released a genomic variation database of mouse strains established in Japan and Western Europe, MoG⁺ (https://molossinus.brc.riken.jp/mogplus/), and a database for phenotype-phenotype associations across the mouse phenome using data from the International Mouse Phenotyping Platform. In this review, we describe features of current version of databases related to mouse strain resources in RIKEN BioResource Research Center and discuss future views.

Reprogramming of the histone H3.3 landscape in the early mouse embryo

Epigenetic reprogramming of the zygote involves dynamic incorporation of histone variant H3.3. However, the genome-wide distribution and dynamics of H3.3 during early development remain unknown. Here, we delineate the H3.3 landscapes in mouse oocytes and early embryos. We unexpectedly identify a non-canonical H3.3 pattern in mature oocytes and zygotes, in which local enrichment of H3.3 at active chromatin is suppressed and H3.3 is relatively evenly distributed across the genome. Interestingly, although the non-canonical H3.3 pattern forms gradually during oogenesis, it quickly switches to a canonical pattern at the two-cell stage in a transcription-independent and replication-dependent manner. We find that incorporation of H3.1/H3.2 mediated by chromatin assembly factor CAF-1 is a key process for the de novo establishment of the canonical pattern. Our data suggest that the presence of the non-canonical pattern and its timely transition toward a canonical pattern support the developmental program of early embryos.

Identification of distinct loci for de novo DNA methylation by DNMT3A and DNMT3B during mammalian development

De novo establishment of DNA methylation is accomplished by DNMT3A and DNMT3B. Here, we analyze de novo DNA methylation in mouse embryonic fibroblasts (2i-MEFs) derived from DNA-hypomethylated 2i/L ES cells with genetic ablation of Dnmt3a or Dnmt3b. We identify 355 and 333 uniquely unmethylated genes in Dnmt3a and Dnmt3b knockout (KO) 2i-MEFs, respectively. We find that Dnmt3a is exclusively required for de novo methylation at both TSS regions and gene bodies of Polycomb group (PcG) target developmental genes, while Dnmt3b has a dominant role on the X chromosome. Consistent with this, tissue-specific DNA methylation at PcG target genes is substantially reduced in Dnmt3a KO embryos. Finally, we find that human patients with DNMT3 mutations exhibit reduced DNA methylation at regions that are hypomethylated in Dnmt3 KO 2i-MEFs. In conclusion, here we report a set of unique de novo DNA methylation target sites for both DNMT3 enzymes during mammalian development that overlap with hypomethylated sites in human patients.

Histone H3K9 Methyltransferase G9a in Oocytes Is Essential for Preimplantation Development but Dispensable for CG Methylation Protection

Mammalian histone methyltransferase G9a (also called EHMT2) deposits H3K9me2 on chromatin and is essential for postimplantation development. However, its role in oogenesis and preimplantation development remains poorly understood. We show that H3K9me2-enriched chromatin domains in mouse oocytes are generally depleted of CG methylation, contrasting with their association in embryonic stem and somatic cells. Oocyte-specific disruption of G9a results in reduced H3K9me2 enrichment and impaired reorganization of heterochromatin in oocytes, but only a modest reduction in CG methylation is detected. Furthermore, in both oocytes and 2-cell embryos, G9a depletion has limited impact on the expression of genes and retrotransposons. Although their CG methylation is minimally affected, preimplantation embryos derived from such oocytes show abnormal chromosome segregation and frequent developmental arrest. Our findings illuminate the functional importance of G9a independent of CG methylation in preimplantation development and call into question the proposed role for H3K9me2 in CG methylation protection in zygotes.

Development of assisted reproductive technologies in small animal species for their efficient preservation and production

Assisted reproductive technologies (ARTs) are widely used in the animal industry, human clinics, and for basic research. In small laboratory animal species such as mice, ARTs are essential for the production of animals for experiments, the preservation of genetic resources, and for the generation of new strains of genetically modified animals. The RIKEN BioResource Research Center (BRC) is one of the largest repositories of such animal bioresources, and maintains approximately 9,500 strains of mice with a variety of genetic backgrounds. We have sought to devise ARTs specific to the reproductive and physiological characteristics of each strain. Such ARTs include superovulation, in vitro fertilization (IVF), the cryopreservation of embryos and spermatozoa, transportation of cryopreserved materials and embryo transfer (ET). Of these, superovulation likely has the most influence on animal production because it determines the quantity of starting material for other ARTs. Superovulation using anti-inhibin serum combined with estrous synchronization has resulted in approximately a three-fold increase in production efficiency with IVF-ET in the C57BL/6J strain. Wild-derived strains are important as genetically diverse resources for murine rodents (Genus Mus), and many are unique to the BRC. We have also successfully developed ARTs for more than 50 wild-derived strains, which have been cryopreserved for future use. Our work to improve and develop ARTs for mice and other small laboratory species will contribute to the cost-effectiveness of routine operations at repository centers, and to the provision of high quality animals for research use.

c.753A>G genome editing of a Cdh23ahl allele delays age-related hearing loss and degeneration of cochlear hair cells in C57BL/6J mice

C57BL/6J mice have long been studied as a model of age-related hearing loss (ARHL). In C57BL/6J mice, ARHL begins in the high-frequency range at 3 months of age and spreads toward low frequencies by 10 months of age. We previously confirmed that c.753A>G genome editing of an ahl allele (c.753A) in the cadherin 23 gene (Cdh23) suppressed the onset of ARHL until 12 months of age. We further investigated the hearing phenotypes of the original and genome-edited C57BL/6J-Cdh23^+/+ (c.753G/G) mice until 24 months of age. The hearing tests revealed that most of the C57BL/6J mice maintained good hearing levels until 14 months of age following genome editing of a Cdh23^ahl allele. However, the hearing levels of the C57BL/6J-Cdh23^+/+ mice gradually declined, and severe ARHL developed with increasing age. ARHL in the C57BL/6J mice was correlated with degeneration of the stereocilia in cochlear hair cells. The stereocilia degeneration was rescued in the C57BL/6J-Cdh23^+/+ mice at 12 months of age, but the stereocilia bundles exhibited abnormal phenotypes similar to those of the original C57BL/6J mice at more advanced ages. Therefore, genome editing of Cdh23^ahl did not completely suppress ARHL in C57BL/6J mice. We also compared the hearing levels of C57BL/6J-Cdh23^+/+ mice with those of C3H/HeN and MSM/Ms mice, which carry the Cdh23⁺ allele. The severity and onset patterns of ARHL in the C57BL/6J-Cdh23^+/+ mice differed from those observed in other Cdh23^+/+ mice. Therefore, we hypothesize that other susceptible and/or resistant alleles of ARHL exist in the genetic backgrounds of these mice.

Characterization of genetic-origin-dependent monoallelic expression in mouse embryonic stem cells

Monoallelic gene expression occurs in various mammalian cells and can be regulated genetically, epigenetically and/or stochastically. We identified 145 monoallelically expressed genes (MoEGs), including seven known imprinted genes, in mouse embryonic stem cells (ESCs) derived from reciprocal F1 hybrid blastocysts and cultured in 2i/LIF. As all MoEGs except for the imprinted genes were expressed in a genetic-origin-dependent manner, we focused on this class of MoEGs for mechanistic studies. We showed that a majority of the genetic-origin-dependent MoEGs identified in 2i/LIF ESCs remain monoallelically expressed in serum/LIF ESCs, but become more relaxed or even biallelically expressed upon differentiation. These MoEGs and their regulatory regions were highly enriched for single nucleotide polymorphisms. In addition, some MoEGs were associated with retrotransposon insertions/deletions, consistent with the fact that certain retrotransposons act as regulatory elements in pluripotent stem cells. Interestingly, most MoEGs showed allelic differences in enrichment of histone H3K27me and H3K4me marks, linking allelic epigenetic differences and monoallelic expression. In contrast, there was little or no allelic difference in CpG methylation or H3K9me. Taken together, our study highlights the impact of genetic variation including single nucleotide polymorphisms and retrotransposon insertions/deletions on monoallelic epigenetic marks and expression in ESCs.

Centromere localization in medaka fish based on half-tetrad analysis

Gene-centromere (G-C) mapping provides insight into vertebrate genome composition, structure and evolution. Although medaka fish are important experimental animals, no genome-wide G-C map of medaka has been constructed. In this study, we used 112 interspecific triploid hybrids and 152 DNA markers to make G-C maps of all 24 linkage groups (LGs). Under the assumption of 50% interference, 24 centromeres were localized onto all corresponding medaka LGs. Comparison with 21 centromere positions deduced from putative centromeric repeats revealed that 19 were localized inside the centromeric regions of the G-C maps, whereas two were not. Based on the centromere positions indicated in the G-C maps and those of centromeric repeats on each LG, we classified chromosomes as either biarmed or monoarmed; n = 24 = 10 metacentrics/submetacentrics + 14 subtelocentrics/acrocentrics, which is consistent with the results of previous karyological reports. This study helps to elucidate genome evolution mechanisms, and integrates physical and genetic maps with karyological information of medaka.

A method for identifying allele-specific hydroxymethylation

We previously identified sequence-dependent allele-specific methylation (sd-ASM) in adult human peripheral blood leukocytes, in which ASM occurs in cis depending on adjacent polymorphic sequences. A number of groups have identified sd-ASM sites in the human and mouse genomes, illustrating the prevalence of sd-ASM in mammalian genomes. In addition, sd-ASM can lead to sequence-dependent allele-specific expression of neighbouring genes. Imprinted genes also often exhibit parent-of-origin-dependent allele-specific methylation (pd-ASM), which causes parent-of-origin-dependent allele-specific expression. However, whether most of the already known sd-ASM and pd-ASM sites are methylated or hydroxymethylated remains unclear due to technical restrictions. Accordingly, a novel method that enables examination of allelic methylation and hydroxymethylation status and also overcomes the drawbacks of conventional methods is needed. Such a method could also be used to elucidate the mechanisms underlying polymorphism-associated inter-individual differences in disease susceptibility and the mechanism of genomic imprinting. Here, we developed a simple method to determine allelic hydroxymethylation status and identified novel sequence- and parent-of-origin-dependent allele-specific hydroxymethylation sites. Correlation analyses of TF binding sequences and methylation or hydroxymethylation between three mouse strains revealed the involvement of Pax5 in strain-specific methylation and hydroxymethylation in exon 7 of Pdgfrb.

Monoallelic, antisense and total RNA transcription in an in vitro neural differentiation system based on F1 hybrid mice

We developed an in vitro system to differentiate embryonic stem cells (ESCs) derived from reciprocally crossed F1 hybrid mice into neurons, and used it to investigate poly(A)+ and total RNA transcription at different stages of cell differentiation. By comparing expression profiles of transcripts assembled from 20 RNA sequencing datasets [2 alleles×(2 cell lines×4 time-points+2 mouse brains)], the relative influence of strain, cell and parent specificities to overall expression could be assessed. Divergent expression profiles of ESCs converged tightly at neural progenitor stage. Patterns of temporal variation of monoallelically expressed transcripts and antisense transcripts were quantified. Comparison of sense and antisense transcript pairs within the poly(A)+ sample, within the total RNA sample, and across poly(A)+ and total RNA samples revealed distinct rates of pairs showing anti-correlated expression variation. Unique patterns of sharing of poly(A)+ and poly(A)- transcription were identified in distinct RNA species. Regulation and functionality of monoallelic expression, antisense transcripts and poly(A)- transcription remain elusive. We demonstrated the effectiveness of our approach to capture these transcriptional activities, and provided new resources to elucidate the mammalian developmental transcriptome.

Nested retrotransposition in the East Asian mouse genome causes the classical nonagouti mutation

Black coat color (nonagouti) is a widespread classical mutation in laboratory mouse strains. The intronic insertion of endogenous retrovirus VL30 in the nonagouti (a) allele of agouti gene was previously reported as the cause of the nonagouti phenotype. Here, we report agouti mouse strains from East Asia that carry the VL30 insertion, indicating that VL30 alone does not cause the nonagouti phenotype. We find that a rare type of endogenous retrovirus, β4, was integrated into the VL30 region at the a allele through nested retrotransposition, causing abnormal splicing. Targeted complete deletion of the β4 element restores agouti gene expression and agouti coat color, whereas deletion of β4 except for a single long terminal repeat results in black-and-tan coat color. Phylogenetic analyses show that the a allele and the β4 retrovirus originated from an East Asian mouse lineage most likely related to Japanese fancy mice. These findings reveal the causal mechanism and historic origin of the classical nonagouti mutation.

ISoLDE: a data-driven statistical method for the inference of allelic imbalance in datasets with reciprocal crosses

MOTIVATION

Allelic imbalance (AI), i.e. the unequal expression of the alleles of the same gene in a single cell, affects a subset of genes in diploid organisms. One prominent example of AI is parental genomic imprinting, which results in parent-of-origin-dependent, mono-allelic expression of a limited number of genes in metatherian and eutherian mammals and in angiosperms. Currently available methods for identifying AI rely on data modeling and come with the associated limitations.

RESULTS

We have designed ISoLDE (Integrative Statistics of alleLe Dependent Expression), a novel nonparametric statistical method that takes into account both AI and the characteristics of RNA-seq data to infer allelic expression bias when at least two biological replicates are available for reciprocal crosses. ISoLDE learns the distribution of a specific test statistic from the data and calls genes 'allelically imbalanced', 'bi-allelically expressed' or 'undetermined'. Depending on the number of replicates, predefined thresholds or permutations are used to make calls. We benchmarked ISoLDE against published methods, and showed that ISoLDE compared favorably with respect to sensitivity, specificity and robustness to the number of replicates. Using ISoLDE on different RNA-seq datasets generated from hybrid mouse tissues, we did not discover novel imprinted genes (IGs), confirming the most conservative estimations of IG number.

AVAILABILITY AND IMPLEMENTATION

ISoLDE has been implemented as a Bioconductor package available at http://bioconductor.org/packages/ISoLDE/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genetic modifiers of rodent animal models: the role in cataractogenesis

Visual impairment leads to a decrease in quality of life. Cataract is the most commonly observed ocular disease in humans that causes vision disorders. The risk factors associated with cataract development include aging, infections, eye injuries, environmental causes, such as radiation and exposure to ultraviolet rays in sunlight, and genetic mutations. Additionally, several cataract patients display phenotypic heterogeneity, suggesting the role of genetic modifiers in the modulation of severity and onset time of cataractogenesis. However, the genetic modifiers associated with cataract have not been identified in humans yet. In contrast, the identification and mapping of genetic modifiers have been successfully carried out in mice and rats. In this review, we focus on the genetic modifiers of cataract in the rodent models.

De Novo DNA Methylation at Imprinted Loci during Reprogramming into Naive and Primed Pluripotency

CpG islands (CGIs) including those at imprinting control regions (ICRs) are protected from de novo methylation in somatic cells. However, many cancers often exhibit CGI hypermethylation, implying that the machinery is impaired in cancer cells. Here, we conducted a comprehensive analysis of CGI methylation during somatic cell reprogramming. Although most CGIs remain hypomethylated, a small subset of CGIs, particularly at several ICRs, was often de novo methylated in reprogrammed pluripotent stem cells (PSCs). Such de novo ICR methylation was linked with the silencing of reprogramming factors, which occurs at a late stage of reprogramming. The ICR-preferred CGI hypermethylation was similarly observed in human PSCs. Mechanistically, ablation of Dnmt3a prevented PSCs from de novo ICR methylation. Notably, the ICR-preferred CGI hypermethylation was observed in pediatric cancers, while adult cancers exhibit genome-wide CGI hypermethylation. These results may have important implications in the pathogenesis of pediatric cancers and the application of PSCs.

•

Several ICRs are de novo methylated in reprogrammed PSCs

•

De novo ICR methylation in iPSCs is linked with transgene silencing

•

Depletion of Dnmt3a prevents reprogrammed PSCs from de novo ICR methylation

•

Pediatric cancers exhibit reprogrammed PSC-like aberration in CGI methylation

Conversion of random X-inactivation to imprinted X-inactivation by maternal PRC2

Imprinted X-inactivation silences genes exclusively on the paternally-inherited X-chromosome and is a paradigm of transgenerational epigenetic inheritance in mammals. Here, we test the role of maternal vs. zygotic Polycomb repressive complex 2 (PRC2) protein EED in orchestrating imprinted X-inactivation in mouse embryos. In maternal-null (Eedm-/-) but not zygotic-null (Eed-/-) early embryos, the maternal X-chromosome ectopically induced Xist and underwent inactivation. Eedm-/- females subsequently stochastically silenced Xist from one of the two X-chromosomes and displayed random X-inactivation. This effect was exacerbated in embryos lacking both maternal and zygotic EED (Eedmz-/-), suggesting that zygotic EED can also contribute to the onset of imprinted X-inactivation. Xist expression dynamics in Eedm-/- embryos resemble that of early human embryos, which lack oocyte-derived maternal PRC2 and only undergo random X-inactivation. Thus, expression of PRC2 in the oocyte and transmission of the gene products to the embryo may dictate the occurrence of imprinted X-inactivation in mammals.

Genome-wide stability of the DNA replication program in single mammalian cells

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a genome-wide methodology that measures copy number differences between replicated and unreplicated DNA. Using scRepli-seq, we demonstrate that replication-domain organization is conserved among individual mouse embryonic stem cells (mESCs). Differentiated mESCs exhibited distinct profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication profiles of homologous autosomes, while the inactive X chromosome was clearly replicated later than its active counterpart. However, a small degree of cell-to-cell replication-timing heterogeneity was present, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

Multiple Conserved Elements Structuring Inverted Repeats in the Mammalian Coat Color-Related Gene Asip

In the agouti signaling gene protein (Asip) of the house mouse (Mus musculus), inverted repeat (IR) arrays are known to exist in a non-coding region adjacent to the ventral-specific promoter region and the accompanying two exons (exons 1A and 1A'), which are around 100 kb upstream from the amino acid coding regions of exons 2, 3, and 4. To determine the gene structure of mammalian Asip and to elucidate trends in its evolution, non-coding sequences of six rodent (mouse, rat, Chinese hamster, squirrel, guinea pig, and naked mole rat) and three non-rodent (rabbit, human, and cow) species were retrieved from databases and compared. Our homology search analyses revealed the presence of three to five highly conserved non-coding elements (CNE). These CNEs were found to form IRs in rodents and lagomorphs. Combinations of IRs were further shown to build symmetric, long IR arrays. Intra- and inter-specific comparisons of the sequences of three universal CNEs showed homogeneity between CNE pairs within species. This implies that certain evolutionary constraints maintained the IR structure in the rodent and rabbit species.

Allele-specific RNA imaging shows that allelic imbalances can arise in tissues through transcriptional bursting

Extensive cell-to-cell variation exists even among putatively identical cells, and there is great interest in understanding how the properties of transcription relate to this heterogeneity. Differential expression from the two gene copies in diploid cells could potentially contribute, yet our ability to measure from which gene copy individual RNAs originated remains limited, particularly in the context of tissues. Here, we demonstrate quantitative, single molecule allele-specific RNA FISH adapted for use on tissue sections, allowing us to determine the chromosome of origin of individual RNA molecules in formaldehyde-fixed tissues. We used this method to visualize the allele-specific expression of Xist and multiple autosomal genes in mouse kidney. By combining these data with mathematical modeling, we evaluated models for allele-specific heterogeneity, in particular demonstrating that apparent expression from only one of the alleles in single cells can arise as a consequence of low-level mRNA abundance and transcriptional bursting.

In mammals, most cells of the body contain two genetic datasets: one from the mother and one from the father, and—in theory—these two sets of information could contribute equally to produce the molecules in a given cell. In practice, however, this is not always the case, which can have dramatic implications for many traits, including visible features (such as fur color) and even disease outcomes. However, it remains difficult to measure the parental origin of individual molecules in a given cell and thus to assess what processes lead to an imbalance of the maternal and paternal contribution. We adapted a microscopy technique—called allele-specific single molecule RNA FISH—that uses a combination of fluorescent tags to specifically label one type of molecule, RNA, depending on its origin, for use in mouse kidney sections. Focusing on RNAs that were previously reported to show imbalance, we performed measurements and combined these with mathematical modeling to quantify imbalance in tissues and explain how these can arise. We found that we could recapitulate the observed imbalances using models that only take into account the random processes that produce RNA, without the need to invoke special regulatory mechanisms to create unequal contributions.

Heterochromatin delays CRISPR-Cas9 mutagenesis but does not influence the outcome of mutagenic DNA repair

Genome editing occurs in the context of chromatin, which is heterogeneous in structure and function across the genome. Chromatin heterogeneity is thought to affect genome editing efficiency, but this has been challenging to quantify due to the presence of confounding variables. Here, we develop a method that exploits the allele-specific chromatin status of imprinted genes in order to address this problem in cycling mouse embryonic stem cells (mESCs). Because maternal and paternal alleles of imprinted genes have identical DNA sequence and are situated in the same nucleus, allele-specific differences in the frequency and spectrum of mutations induced by CRISPR-Cas9 can be unequivocally attributed to epigenetic mechanisms. We found that heterochromatin can impede mutagenesis, but to a degree that depends on other key experimental parameters. Mutagenesis was impeded by up to 7-fold when Cas9 exposure was brief and when intracellular Cas9 expression was low. In contrast, the outcome of mutagenic DNA repair was unaffected by chromatin state, with similar efficiencies of homology-directed repair (HDR) and deletion spectra on maternal and paternal chromosomes. Combined, our data show that heterochromatin imposes a permeable barrier that influences the kinetics, but not the endpoint, of CRISPR-Cas9 genome editing and suggest that therapeutic applications involving low-level Cas9 exposure will be particularly affected by chromatin status.

Frequent monoallelic or skewed expression for developmental genes in CNS-derived cells and evidence for balancing selection

Cellular mosaicism due to monoallelic autosomal expression (MAE), with cell selection during development, is becoming increasingly recognized as prevalent in mammals, leading to interest in understanding its extent and mechanism(s). We report here use of clonal cell lines derived from the CNS of adult female [Formula: see text] hybrid (C57BL/6 X JF1) mice to characterize MAE as neural stem cells (nscs) differentiate to astrocyte-like cells (asls). We found that different subsets of genes show MAE in the two populations of cells; in each case, there is strong enrichment for genes specific to the respective developmental state. Genes that exhibit MAE are 22% of nsc-specific genes and 26% of asl-specific genes. Moreover, the promoters of genes with MAE have reduced CpG dinucleotides but increased CpG differences between the two parental mouse strains. Extending the study of variability to wild populations of mice, we found evidence for balancing selection as a contributing force in evolution of those genes showing developmental specificity (i.e., expressed in either nsc or asl), not just for genes showing MAE. Furthermore, we found that genes showing skewed allelic expression (SKE) were similarly enriched among cell type-specific genes and also showed a heightened probability of balancing selection. Thus, developmental stage-specific genes and genes with MAE or SKE seem to make up overlapping classes subject to selection for increased diversity. The implications of these results for development and evolution are discussed in the context of a model with stochastic epigenetic modifications taking place only during a relatively brief developmental window.

Evolutionary Dynamics of Copy Number and Meiotic Recombination in Murine 5S rDNA: Possible Involvement of Natural Selection

We investigated evolutionary trends of the 5S ribosomal RNA gene in the house mouse, Mus musculus. First, we assessed the 5S cluster and copy numbers in eight laboratory strains by pulsed-field gel electrophoresis. The copy numbers in seven lines were estimated to be around 130-170 copies per cluster, with 63 copies in the remaining strain, implying that the copy number can change drastically and has been maintained under certain evolutionary constraints at ~ 140 copies. Second, we addressed the frequency of meiotic recombination mediated by the 5S cluster by performing a mating experiment with laboratory strains, and found that the 5S cluster did not accelerate recombination events. Third, we surveyed recombination events of the 5S-containing chromosome region in wild mice from the Japanese Islands, where the two subspecies lineages, M. m. castaneus and M. m. musculus, are historically mingled, and found that the influence of the 5S cluster on meiotic recombination was limited. Finally, we examined the nucleotide diversity of six genes in the neighboring regions of the 5S cluster and found reduced genetic diversity in the regions on both sides of the cluster, suggesting the involvement of either positive or background selection in the population-level sequence similarity of the 5S clusters. Therefore, the mouse 5S genes are considered to be evolving toward sequence similarity within a given cluster by certain intrachromosomal mechanisms and toward sharing of a specific 5S cluster within a population by certain selective processes.

Molecular analyses of the agouti allele in the Japanese house mouse identify a novel variant of the agouti gene

It has been thought that the Japanese house mouse carries the A^w allele at the agouti locus causing light-colored bellies, but they do not always show this coloration. Thus, the presence of the A^w allele seems to be doubtful in them. To ascertain whether the A^w allele is present, a two-pronged approach was used. First, we compared lengths of DNA fragments obtained from three PCRs conducted on them to the known fragment sizes generated from mouse strains exhibiting homozygosities of either a/a, A/A, or A^w/A^w. PCR I, PCR II, and PCR III amplify only in the A and A^w alleles, the a and A^w alleles, and the a allele, respectively, and we detected amplifications in strains with A/A and A^w/A^w by PCR I, in those with a/a and the Japanese house mouse by PCR II, and in those with a/a by PCR III. Second, we sequenced the exon 1A region of the agouti gene and obtained sequences corresponding to the above strains and the Japanese house mouse, but their sequences were similar to those of the a allele. We concluded that their agouti allele is not identical to the A^w allele and seems to be a novel type similar to the a allele.

Genetic Dissection of Trabecular Bone Structure with Mouse Intersubspecific Consomic Strains

Trabecular bone structure has an important influence on bone strength, but little is known about its genetic regulation. To elucidate the genetic factor(s) regulating trabecular bone structure, we compared the trabecular bone structures of two genetically remote mouse strains, C57BL/6J and Japanese wild mouse-derived MSM/Ms. Phenotyping by X-ray micro-CT revealed that MSM/Ms has structurally more fragile trabecular bone than C57BL/6J. Toward identification of genetic determinants for the difference in fragility of trabecular bone between the two mouse strains, we employed phenotype screening of consomic mouse strains in which each C57BL/6J chromosome is substituted by its counterpart from MSM/Ms. The results showed that many chromosomes affect trabecular bone structure, and that the consomic strain B6-Chr15^MSM, carrying MSM/Ms-derived chromosome 15 (Chr15), has the lowest values for the parameters BV/TV, Tb.N, and Conn.D, and the highest values for the parameters Tb.Sp and SMI. Subsequent phenotyping of subconsomic strains for Chr15 mapped four novel trabecular bone structure-related QTL (Tbsq1-4) on mouse Chr15. These results collectively indicate that genetic regulation of trabecular bone structure is highly complex, and that even in the single Chr15, the combined action of the four Tbsqs controls the fragility of trabecular bone. Given that Tbsq4 is syntenic to human Chr 12q12-13.3, where several bone-related SNPs are assigned, further study of Tbsq4 should facilitate our understanding of the genetic regulation of bone formation in humans.

Switching of dominant retrotransposon silencing strategies from posttranscriptional to transcriptional mechanisms during male germ-cell development in mice

Mammalian genomes harbor millions of retrotransposon copies, some of which are transpositionally active. In mouse prospermatogonia, PIWI-interacting small RNAs (piRNAs) combat retrotransposon activity to maintain the genomic integrity. The piRNA system destroys retrotransposon-derived RNAs and guides de novo DNA methylation at some retrotransposon promoters. However, it remains unclear whether DNA methylation contributes to retrotransposon silencing in prospermatogonia. We have performed comprehensive studies of DNA methylation and polyA(+) RNAs (transcriptome) in developing male germ cells from Pld6/Mitopld and Dnmt3l knockout mice, which are defective in piRNA biogenesis and de novo DNA methylation, respectively. The Dnmt3l mutation greatly reduced DNA methylation levels at most retrotransposons, but its impact on their RNA abundance was limited in prospermatogonia. In Pld6 mutant germ cells, although only a few retrotransposons exhibited reduced DNA methylation, many showed increased expression at the RNA level. More detailed analysis of RNA sequencing, nascent RNA quantification, profiling of cleaved RNA ends, and the results obtained from double knockout mice suggest that PLD6 works mainly at the posttranscriptional level. The increase in retrotransposon expression was larger in Pld6 mutants than it was in Dnmt3l mutants, suggesting that RNA degradation by the piRNA system plays a more important role than does DNA methylation in prospermatogonia. However, DNA methylation had a long-term effect: hypomethylation caused by the Pld6 or Dnmt3l mutation resulted in increased retrotransposon expression in meiotic spermatocytes. Thus, posttranscriptional silencing plays an important role in the early stage of germ cell development, then transcriptional silencing becomes important in later stages. In addition, intergenic and intronic retrotransposon sequences, in particular those containing the antisense L1 promoters, drove ectopic expression of nearby genes in both mutant spermatocytes, suggesting that retrotransposon silencing is important for the maintenance of not only genomic integrity but also transcriptomic integrity.

Retrotransposons are a class of transposable elements, of which mobility has mutagenic potential. Therefore, it is important to regulate the expression of retrotransposons for maintaining the genomic integrity. In male germ cells, DNA methylation and the piRNA system are thought to play roles in retrotransposon silencing. However, genome-wide DNA methylation is once erased (in primordial germ cells) and reestablished (in prospermatogonia) during development. In prospermatogonia, piRNAs guide de novo DNA methylation at some retrotransposons. To clarify the contribution of DNA methylation and the piRNA system to retrotransposon silencing in the course of male germ cell development, we analyzed DNA methylation and RNA expression in Dnmt3l and Pld6 knockout mice, which are defective in de novo DNA methylation and piRNA biogenesis, respectively. Our results reveal that, in prospermatogonia, the piRNA system works mainly at the posttranscriptional level, and plays a more important role than does DNA methylation in retrotransposon silencing. However, DNA methylation becomes much more important in later stages when germ cells enter meiosis (in spermatocytes). We also found that hypomethylated retrotransposons can drive ectopic expression of nearby genes; therefore, their transcriptional silencing by DNA methylation is important for maintaining the transcriptomic integrity as well.