Methylation imprinting was observed of mouse mo-2 macrosatellite on the pseudoautosomal region but not on chromosome 9

Mouse Genome Informatics: an integrated knowledgebase system for the laboratory mouse

Mouse Genome Informatics (MGI) is a federation of expertly curated information resources designed to support experimental and computational investigations into genetic and genomic aspects of human biology and disease using the laboratory mouse as a model system. The Mouse Genome Database (MGD) and the Gene Expression Database (GXD) are core MGI databases that share data and system architecture. MGI serves as the central community resource of integrated information about mouse genome features, variation, expression, gene function, phenotype, and human disease models acquired from peer-reviewed publications, author submissions, and major bioinformatics resources. To facilitate integration and standardization of data, biocuration scientists annotate using terms from controlled metadata vocabularies and biological ontologies (e.g. Mammalian Phenotype Ontology, Mouse Developmental Anatomy, Disease Ontology, Gene Ontology, etc.), and by applying international community standards for gene, allele, and mouse strain nomenclature. MGI serves basic scientists, translational researchers, and data scientists by providing access to FAIR-compliant data in both human-readable and compute-ready formats. The MGI resource is accessible at https://informatics.jax.org. Here, we present an overview of the core data types represented in MGI and highlight recent enhancements to the resource with a focus on new data and functionality for MGD and GXD.

Evolutionary dynamics of pseudoautosomal region 1 in humans and great apes

BACKGROUND

The pseudoautosomal region 1 (PAR1) is a 2.7 Mb telomeric region of human sex chromosomes. PAR1 has a crucial role in ensuring proper segregation of sex chromosomes during male meiosis, exposing it to extreme recombination and mutation processes. We investigate PAR1 evolution using population genomic datasets of extant humans, eight populations of great apes, and two archaic human genome sequences.

RESULTS

We find that PAR1 is fast evolving and closer to evolutionary nucleotide equilibrium than autosomal telomeres. We detect a difference between substitution patterns and extant diversity in PAR1, mainly driven by the conflict between strong mutation and recombination-associated fixation bias at CpG sites. We detect excess C-to-G mutations in PAR1 of all great apes, specific to the mutagenic effect of male recombination. Despite recent evidence for Y chromosome introgression from humans into Neanderthals, we find that the Neanderthal PAR1 retained similarity to the Denisovan sequence. We find differences between substitution spectra of these archaics suggesting rapid evolution of PAR1 in recent hominin history. Frequency analysis of alleles segregating in females and males provided no evidence for recent sexual antagonism in this region. We study repeat content and double-strand break hotspot regions in PAR1 and find that they may play roles in ensuring the obligate X-Y recombination event during male meiosis.

CONCLUSIONS

Our study provides an unprecedented quantification of population genetic forces governing PAR1 biology across extant and extinct hominids. PAR1 evolutionary dynamics are predominantly governed by recombination processes with a strong impact on mutation patterns across all species.

Ensuring meiotic DNA break formation in the mouse pseudoautosomal region

Sex chromosomes in males of most eutherian species share only a diminutive homologous segment, the pseudoautosomal region (PAR), wherein double-strand break (DSB) formation, pairing, and crossing over must occur for correct meiotic segregation^1,2. How cells ensure PAR recombination is unknown. Here we delineate an unexpected dynamic ultrastructure of the PAR and identify controlling cis- and trans-acting factors that make this the hottest area of DSB formation in the male mouse genome. Before break formation, multiple DSB-promoting factors hyper-accumulate in the PAR, its chromosome axes elongate, and the sister chromatids separate. These phenomena are linked to heterochromatic mo-2 minisatellite arrays and require MEI4 and ANKRD31 proteins but not axis components REC8 or HORMAD1. We propose that the repetitive PAR sequence confers unique chromatin and higher order structures crucial for recombination. Chromosome synapsis triggers collapse of the elongated PAR structure and, remarkably, oocytes can be reprogrammed to display spermatocyte-like PAR DSB levels simply by delaying or preventing synapsis. Thus, sexually dimorphic behavior of the PAR rests in part on kinetic differences between the sexes for a race between maturation of PAR structure, DSB formation, and completion of pairing and synapsis. Our findings establish a mechanistic paradigm of sex chromosome recombination during meiosis.

Genetic variation of melatonin productivity in laboratory mice under domestication

Melatonin is a pineal hormone produced at night; however, many strains of laboratory mice are deficient in melatonin. Strangely enough, the gene encoding HIOMT enzyme (also known as ASMT) that catalyzes the last step of melatonin synthesis is still unidentified in the house mouse (Mus musculus) despite the completion of the genome sequence. Here we report the identification of the mouse Hiomt gene, which was mapped to the pseudoautosomal region (PAR) of sex chromosomes. The gene was highly polymorphic, and nonsynonymous SNPs were found in melatonin-deficient strains. In C57BL/6 strain, there are two mutations, both of which markedly reduce protein expression. Mutability of the Hiomt likely due to a high recombination rate in the PAR could be the genomic basis for the high prevalence of melatonin deficiency. To understand the physiologic basis, we examined a wild-derived strain, MSM/Ms, which produced melatonin more under a short-day condition than a long-day condition, accompanied by increased Hiomt expression. We generated F2 intercrosses between MSM/Ms and C57BL/6 strains and N2 backcrosses to investigate the role of melatonin productivity on the physiology of mice. Although there was no apparent effect of melatonin productivity on the circadian behaviors, testis development was significantly promoted in melatonin-deficient mice. Exogenous melatonin also had the antigonadal action in mice of a melatonin-deficient strain. These findings suggest a favorable impact of melatonin deficiency due to Hiomt mutations on domestic mice in breeding colonies.

The problem of the eukaryotic genome size

The current state of knowledge concerning the unsolved problem of the huge interspecific eukaryotic genome size variations not correlating with the species phenotypic complexity (C-value enigma also known as C-value paradox) is reviewed. Characteristic features of eukaryotic genome structure and molecular mechanisms that are the basis of genome size changes are examined in connection with the C-value enigma. It is emphasized that endogenous mutagens, including reactive oxygen species, create a constant nuclear environment where any genome evolves. An original quantitative model and general conception are proposed to explain the C-value enigma. In accordance with the theory, the noncoding sequences of the eukaryotic genome provide genes with global and differential protection against chemical mutagens and (in addition to the anti-mutagenesis and DNA repair systems) form a new, third system that protects eukaryotic genetic information. The joint action of these systems controls the spontaneous mutation rate in coding sequences of the eukaryotic genome. It is hypothesized that the genome size is inversely proportional to functional efficiency of the anti-mutagenesis and/or DNA repair systems in a particular biological species. In this connection, a model of eukaryotic genome evolution is proposed.

Human megasatellite DNA RS447: copy-number polymorphisms and interspecies conservation

We previously isolated a novel 4.7-kb RS447 sequence, which tandemly repeated approximately 50-70 copies and resided on human chromosome 4p15 (M. Kogi et al., 1997, Genomics 42: 278-283). Another tandem array (or arrays) of several RS447 copies was hereby identified on the distal part of chromosome 8p. To analyze copy-number polymorphisms of the RS447 repeats, genomic DNA samples of eight nonkindred Japanese were subjected to pulsed-field gel electrophoresis. The copy numbers of the RS447 tandem arrays on 4p15 varied drastically from allele to allele and ranged from approximately 34 to 94 copies. All eight Japanese subjects were apparently heterozygous for the RS447 copy number, and 12 copy-number-different alleles have been at least clearly distinguished. The RS447 tandem repeats were thus found to be hypervariable and highly polymorphic in a human population. The RS447 sequences, however, do not appear to be either "selfish" or "junk" DNA. The unit size and sequence of RS447 were found to be very similar between members in the human genome. The unit size of 4746 bp comprises a putative open reading frame of 1590 bp. The RS447 sequence was well conserved in all the tested mammalian species. The head-to-tail tandem repetitive structure in the RS447 homologs was also confirmed in those species. The RS447 sequence is, therefore, considered to consist of a new class of tandemly repeated satellite DNA elements in the mammalian genome, which may thus be called "megasatellite DNA."

DNA binding discrimination of the murine DNA cytosine-C5 methyltransferase

Mammalian DNA cytosine-C5 methyltransferase modifies the CpG dinucleotide in the context of many different genomic sequences. A rigorous DNA binding assay was developed for the murine enzyme and used to define how sequences flanking the CpG dinucleotide affect the stability of the enzyme:DNA complex. Oligonucleotides containing a single CpG site form reversible 1:1 complexes with the enzyme that are sequence-specific. A guanine/cytosine-rich 30 base-pair sequence, a mimic of the GC-box cis-element, bound threefold more tightly than an adenine/thymine-rich sequence, a mimic of the cyclic AMP responsive element. However, the binding discrimination between hemi- and unmethylated forms of these DNA substrates was small, as we previously observed at the K(m)DNA level (Biochemistry, 35, 7308-7315 (1996)). Single-stranded substrates are bound much more weakly than double-stranded DNA forms. An in vitro screening method was used to select for CpG flanking sequence preferences of the DNA methyltransferase from a large, divergent population of DNA substrates. After five iterative rounds of increasing selective pressure, guanosine/cytosine-rich sequences were abundant and contributed to binding stabilization for at least 12 base-pairs on either side of a central CpG. Our results suggest a read-out of sequence-dependent conformational features, such as helical flexibility, minor groove dimensions and critical phosphate orientation and mobility, rather than interactions with specific bases over the course of two complete helical turns. Thus, both studies reveal a preference for guanosine/cytosine deoxynucleotides flanking the cognate CpG. The enzyme specificity for similar sequences in the genome may contribute to the in vivo functions of this vital enzyme.

Difference in chromatin packaging between active and inactive X chromosomes by fractionation and allele-specific detection

Using a novel method consisting of chromatin fractionation and allele-specific detection, chromatin packaging is compared between active X (Xa) and inactive X (Xi) chromosomes for five tumor cell clones that were derived from inter-subspecific F1 female mice. Separation of heterochromatic (H) and euchromatic (E) fractions is monitored by hybridization with subtelomeric satellite DNA and ribosomal RNA gene and by PCR amplification of p53 gene/pseudogene with one primer set. The H fraction was enriched with satellite and p53 pseudogene probably existing in heterochromatic regions while the E fraction showed inverse, suggesting fair separation. Analysis with seven marker and three gene loci revealed concentration of alleles on Xi in the H fraction and those on Xa in the E fraction, though the concentration levels varied. This implies that the packaging level of Xi is higher than that of active or inactive euchromatin on Xa. Intriguingly, one cell line showed biallelic expression and chromatin relaxation of the Pgk-1 locus, suggesting that the relaxation occur regionally on X chromosome.

High frequency de novo alterations in the long-range genomic structure of the mouse pseudoautosomal region

The pseudoautosomal region (PAR) is a segment of shared homology between the X and Y chromosomes. Here we report physical linkage of three mouse PAR probes: DXYHgu1, DXYMov15 and (TTAGGG)n. Steroid sulphatase (Sts) maps distal to these probes, indicating that there is an internal array of the telomere sequence (TTAGGG)n in the PAR. Pseudoautosomal PacI restriction fragments, up to 2 Mb in size, are unstable in C57BL/6 x C57BL/6 crosses. New alleles, often several hundred kilobases different in size, occur at a sex-averaged rate of approximately 30% per allele. Such frequent large-scale germline genome arrangements are without precedent in mammals.

Multiple roles for DNA methylation in gametic imprinting

Allele-specific DNA methylation has been observed for all tested imprinted genes and has a clear role in the imprinting mechanism. It remains to be resolved whether this role is to act as the gametic imprinting signal or to cause or maintain allele-specific expression.

Structural variation of the pseudoautosomal region between and within inbred mouse strains

The pseudoautosomal region (PAR) is a segment of shared homology between the sex chromosomes. Here we report additional probes for this region of the mouse genome. Genetic and fluorescence in situ hybridization analyses indicate that one probe, PAR-4, hybridizes to the pseudoautosomal telomere and a minor locus at the telomere of chromosome 9 and that a PCR assay based on the PAR-4 sequence amplifies only the pseudoautosomal locus (DXYHgu1). The region detected by PAR-4 is structurally unstable; it shows polymorphism both between mouse strains and between animals of the same inbred strain, which implies an unusually high mutation rate. Variation occurs in the region adjacent to a (TTAGGG)n array. Two pseudoautosomal probes can also hybridize to the distal telomeres of chromosomes 9 and 13, and all three telomeres contain DXYMov15. The similarity between these telomeres may reflect ancestral telomere-telomere exchange.

Induction of karyotype instability in a murine tumor cell line by quercetin, 2-amino-1-methyl-6 phenylimidazo[4,5-b]pyridine, and okadaic acid, as revealed by transmission distortion of the inactive X chromosome

Quercetin, 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP), and okadaic acid are found in various foods and have been shown to have mutagenic or promoter-like activity. The effects of these three compounds on the transmission of the inactive X chromosome were examined in MST-C6 murine tumor cells, which were derived from hybrid F1 mice from matings between C57BL/6 and MSM mice. Polymerase chain reaction analysis using polymorphic markers on the X chromosome detected transmission distortion of the inactive X chromosome due to nondisjunction as a copy-number imbalance in allelic bands. The cells exposed to all three chemicals (but not untreated cells) exhibited such imbalances at high frequencies under exposure conditions similar to those in previous experiments in which tumor progression and recombination were observed. The cells also showed increased frequencies of tumor formation when subcutaneously injected. These results suggest that the three chemicals are capable of inducing transmission distortion of the inactive X chromosome and that such activity may be a causative factor in promoting the tumorigenicity of MST-C6 cells.