Mouse inter-subspecific consomic strains for genetic dissection of quantitative complex traits

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.

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.

Neonatal lethality of mouse A/J-7SM consomic strain is caused by an insertion mutation in the Dchs1 gene

Homosomic mice of the A/J-7^SM consomic mouse strain that introduced the entire chromosome 7 (Chr 7) of SM/J into the A/J strain exhibited neonatal lethality. We tentatively maintained segregating inbred strains (A/J-7A^SM and A/J-7D^SM) in which the central portion of Chr 7 was heterozygous for the A/J and SM/J strains, and the centromeric and telomeric sides of Chr 7 were homozygous for the SM/J strain, instead of the A/J-7^SM strain. Based on the chromosomal constitution of Chr 7 in A/J-7A^SM and A/J-7D^SM mice, the causative gene for neonatal lethality in homosomic mice was suggested to be located within an approximately 1.620 Mb region between D7Mit125 (104.879 Mb) and D7Mit355 (106.499 Mb) on Chr 7. RT-PCR analysis revealed that homosomic mice lacked dachsous cadherin-related 1 (Dchs1), which is located within the D7Mit125 to D7Mit355 region and functions in the regulation of planar cell polarity. Screening for mutations in Dchs1 indicated that homosomic mice possessed an early transposable (ETn)-like sequence in intron 1 of Dchs1. Moreover, an allelism test between Dchs1 ETn-like-insertion alleles detected in homosomic mice and CRISPR/Cas9-induced Dchs1 deletion alleles revealed that Dchs1 is a causative gene for neonatal lethality in homosomic mice. Based on these results, we concluded that in the A/J-7^SM strain, ETn-like elements were inserted into intron 1 of SM/J-derived Dchs1 during strain development, which dramatically reduced Dchs1 expression, thus resulting in neonatal lethality in homosomic mice. Additionally, it was suggested that the timing of lethality in Dchs1 mutant mice is influenced by the genetic background.

Polycomb repressive complexes 1 and 2 are each essential for maintenance of X inactivation in extra-embryonic lineages

In female mammals, one of the two X chromosomes becomes inactivated during development by X-chromosome inactivation (XCI). Although Polycomb repressive complex (PRC) 1 and PRC2 have both been implicated in gene silencing, their exact roles in XCI during in vivo development have remained elusive. To this end, we have studied mouse embryos lacking either PRC1 or PRC2. Here we demonstrate that the loss of either PRC has a substantial impact on maintenance of gene silencing on the inactive X chromosome (Xi) in extra-embryonic tissues, with overlapping yet different genes affected, indicating potentially independent roles of the two complexes. Importantly, a lack of PRC1 does not affect PRC2/H3K27me3 accumulation and a lack of PRC2 does not impact PRC1/H2AK119ub1 accumulation on the Xi. Thus PRC1 and PRC2 contribute independently to the maintenance of XCI in early post-implantation extra-embryonic lineages, revealing that both Polycomb complexes can be directly involved and differently deployed in XCI.

Identification of a Chromosome 1 Substitution Line B6-Chr1BLD as a Novel Hyperlipidemia Model via Phenotyping Screening

Hyperlipidemia is a chronic disease that seriously affects human health. Due to the fact that traditional animal models cannot fully mimic hyperlipidemia in humans, new animal models are urgently needed for basic drug research on hyperlipidemia. Previous studies have demonstrated that the genomic diversity of the wild mice chromosome 1 substitution lines was significantly different from that of laboratory mice, suggesting that it might be accompanied by phenotypic diversity. We first screened the blood lipid-related phenotype of chromosome 1 substitution lines. We found that the male HFD-fed B6-Chr1BLD mice showed more severe hyperlipidemia-related phenotypes in body weight, lipid metabolism and liver lesions. By RNA sequencing and whole-genome sequencing results of B6-Chr1BLD, we found that several differentially expressed single nucleotide polymorphism enriched genes were associated with lipid metabolism-related pathways. Lipid metabolism-related genes, mainly including Aida, Soat1, Scly and Ildr2, might play an initial and upstream role in the abnormal metabolic phenotype of male B6-Chr1BLD mice. Taken together, male B6-Chr1BLD mice could serve as a novel, polygenic interaction-based hyperlipidemia model. This study could provide a novel animal model for accurate clinical diagnosis and precise medicine of hyperlipidemia.

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.

Xanthine Oxidoreductase Inhibitors Suppress the Onset of Exercise-Induced AKI in High HPRT Activity Urat1-Uox Double Knockout Mice

BACKGROUND

Hereditary renal hypouricemia type 1 (RHUC1) is caused by URAT1/SLC22A12 dysfunction, resulting in urolithiasis and exercise-induced AKI (EIAKI). However, because there is no useful experimental RHUC1 animal model, the precise pathophysiologic mechanisms underlying EIAKI have yet to be elucidated. We established a high HPRT activity Urat1-Uox double knockout (DKO) mouse as a novel RHUC1 animal model for investigating the cause of EIAKI and the potential therapeutic effect of xanthine oxidoreductase inhibitors (XOIs).

METHODS

The novel Urat1-Uox DKO mice were used in a forced swimming test as loading exercise to explore the onset mechanism of EIAKI and evaluate related purine metabolism and renal injury parameters.

RESULTS

Urat1-Uox DKO mice had uricosuric effects and elevated levels of plasma creatinine and BUN as renal injury markers, and decreased creatinine clearance observed in a forced swimming test. In addition, Urat1-Uox DKO mice had increased NLRP3 inflammasome activity and downregulated levels of Na⁺-K⁺-ATPase protein in the kidney, as Western blot analysis showed. Finally, we demonstrated that topiroxostat and allopurinol, XOIs, improved renal injury and functional parameters of EIAKI.

CONCLUSIONS

Urat1-Uox DKO mice are a useful experimental animal model for human RHUC1. The pathogenic mechanism of EIAKI was found to be due to increased levels of IL-1β via NLRP3 inflammasome signaling and Na⁺-K⁺-ATPase dysfunction associated with excessive urinary urate excretion. In addition, XOIs appear to be a promising therapeutic agent for the treatment of EIAKI.

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.

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.

Hybrid sterility genes in mice (Mus musculus): a peculiar case of PRDM9 incompatibility

Hybrid sterility is a critical step in the evolution of reproductive barriers between diverging taxa during the process of speciation. Recent studies of young subspecies of the house mouse revealed a multigenic nature and frequent polymorphism of hybrid sterility genes as well as the recurrent engagement of the meiosis-specific gene PR domain-containing 9 (Prdm9) and X-linked loci. Prdm9-controlled hybrid sterility is essentially chromosomal in nature, conditioned by the sequence divergence between subspecies. Depending on the Prdm9 interallelic interactions and the X-linked Hstx2 locus, the same homologs either regularly recombine and synapse, or show impaired DNA DSB repair, asynapsis, and early meiotic arrest. Thus, Prdm9-dependent hybrid sterility points to incompatibilities affecting meiotic recombination as a possible mechanism of reproductive isolation between (sub)species.

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.

Chromosomal Substitution Strategies to Localize Genomic Regions Related to Complex Traits

Chromosomal substitution strategies provide a powerful tool to anonymously reveal the relationship between DNA sequence variants and a normal or disease phenotype of interest. Even in this age of CRISPR-Cas9 genome engineering, the knockdown or overexpression of a gene provides relevant information to our understanding of complex disease only when a close association of an allelic variant with the phenotype has first been established. Limitations of genetic linkage approaches led to the development of more efficient breeding strategies to substitute chromosomal segments from one animal strain into the genetic background of a different strain, enabling a direct comparison of the phenotypes of the strains with variant(s) that differ only at a defined locus. This substitution can be a whole chromosome (consomic), a part of a chromosome (congenic), or as small as only a single or several alleles (subcongenics). In contrast to complete knockout of a specific candidate gene of interest, which simply studies the effects of complete elimination of the gene, the substitution of naturally occurring variants can provide special insights into the functional actions of wild-type alleles. Strategies for production of these inbred strains are reviewed, and a number of examples are used to illustrate the utility of these model systems. Consomic/congenic strains provide a number of experimental advantages in the study of functions of genes and their variants, which are emphasized in this article, such as replication of experimental studies; determination of temporal relationships throughout a life; rigorously controlled experiments in which relations between genotype and phenotype can be tested with the confounding effects of heterogeneous genetic backgrounds, both targeted and multilayered; and "omic" studies performed at many levels of functionality, from molecules to organelles, cells to organs, and organs to organismal behavior across the life span. The application of chromosomal substitution strategies and development of consomic/congenic rat and mouse strains have greatly expanded our knowledge of genomic variants and their phenotypic relationship to physiological functions and to complex diseases such as hypertension and cancer. © 2020 American Physiological Society. Compr Physiol 10:365-388, 2020.

Xanthine oxidoreductase knockout mice with high HPRT activity were not rescued by NAD+ replenishment

Although xanthinuria is nonfatal in human, xanthine oxidoreductase knockout (Xor-KO) mice have only a short lifespan. Hypoxanthine phosphoribosyltransferase activity (HPRT) in human and wild mice is higher than in laboratory mice. The aim of this study was to investigate the underlying mechanisms that give rise to the longer lifespan of high-HPRT/Xor-KO mice. Before Xor-KO mice die, urinary excretion of hypoxanthine increased with a corresponding decrease in excretion of xanthine. The switch of excretion from xanthine to hypoxanthine might be a cause of death for Xor-KO mice, suggesting inhibition of NAD⁺-dependent IMP dehydrogenase. Because hypoxanthine inhibits the synthesis of nicotinamide mononucleotide (NMN), a precursor of NAD⁺, the accumulation of hypoxanthine in Xor-KO mice may cause a depletion in the levels of NAD⁺. Moreover, urinary excretion of urate in high-HPRT/Uox-KO/Xor-KO mice means urate derived from gut microbiota is absorbed by the intestine. Likewise, over excretion of oxypurine in mice may be caused by intestinal absorption of oxypurine. For NAD⁺ replenishment, oral supplementation with 1% L-tryptophan, an alternative precursor of NAD⁺, resulted in a recovery of body weight gain in high-HPRT/Uox-KO/Xor-KO mice. In conclusion, the death of Xor-KO mice by renal failure seems to be caused by a depletion in NAD⁺ levels due to the intracellular accumulation of hypoxanthine. NAD⁺ replenishment by oral supplementation of NMN or tryptophan was complicated by the effect of gut microbiota and failed to rescue high-HPRT/Xor-KO mice. The attenuation of intestinal absorption of oxypurines seems to be necessary to avoid hypoxanthine accumulation and over excretion of oxypurine.

Congenic mapping and candidate gene analysis for streptozotocin-induced diabetes susceptibility locus on mouse chromosome 11

Streptozotocin (STZ) has been widely used to induce diabetes in rodents. Strain-dependent variation in susceptibility to STZ has been reported; however, the gene(s) responsible for STZ susceptibility has not been identified. Here, we utilized the A/J-11^SM consomic strain and a set of chromosome 11 (Chr. 11) congenic strains developed from A/J-11^SM to identify a candidate STZ-induced diabetes susceptibility gene. The A/J strain exhibited significantly higher susceptibility to STZ-induced diabetes than the A/J-11^SM strain, confirming the existence of a susceptibility locus on Chr. 11. We named this locus Stzds1 (STZ-induced diabetes susceptibility 1). Congenic mapping using the Chr. 11 congenic strains indicated that the Stzds1 locus was located between D11Mit163 (27.72 Mb) and D11Mit51 (36.39 Mb). The Mpg gene, which encodes N-methylpurine DNA glycosylase (MPG), a ubiquitous DNA repair enzyme responsible for the removal of alkylated base lesions in DNA, is located within the Stzds1 region. There is a close relationship between DNA alkylation at an early stage of STZ action and the function of MPG. A Sanger sequence analysis of the Mpg gene revealed five polymorphic sites in the A/J genome. One variant, p.Ala132Ser, was located in a highly conserved region among rodent species and in the minimal region for retained enzyme activity of MPG. It is likely that structural alteration of MPG caused by the p.Ala132Ser mutation elicits increased recognition and excision of alkylated base lesions in DNA by STZ.

Deficiency of a pyrroline-5-carboxylate reductase produces the yellowish green cocoon 'Ryokuken' of the silkworm, Bombyx mori

The silkworm cocoon colour has attracted researchers involved in genetics, physiology and ecology for a long time. 'Ryokuken' cocoons are yellowish green in colour due to unusual flavonoids, prolinylflavonols, while 'Sasamayu' cocoons are light green and contain only simple flavonol glucosides. We found a novel gene associated with the cocoon colour change resulting from a change in flavonoid composition and named it Lg (light green cocoon). In the middle silk glands of the + ^Lg /+ ^Lg larvae, 1-pyrroline-5-carboxylic acid (P5C) was found to accumulate due to a decrease in the activity of pyrroline-5-carboxylate reductase (P5CR), an enzyme reducing P5C to proline. Sequence analysis of BmP5CR1, the candidate gene for Lg, revealed a 1.9 kb insertion and a 4 bp deletion within the 1st intron, a 97 bp deletion within the 4th intron, and a > 300 bp insertion within the 3'-UTR, in addition to two amino acid changes on exons 3 and 4 in + ^Lg /+ ^Lg compared to Lg/Lg. Decreased expression of BmP5CR1 was observed in all of the investigated tissues, including the middle silk glands in + ^Lg /+ ^Lg , which was probably caused by structural changes in the intronic regions of BmP5CR1. Furthermore, a BmP5CR1 knockout strain exhibited a yellowish green cocoon with the formation of prolinylflavonols. These results indicate that the yellowish green cocoon is produced by a BmP5CR1 deficiency. To our knowledge, this is the first report showing that the defect of an enzyme associated with intermediate metabolism promotes the conjugation of phytochemicals derived from foods with endogenously accumulating metabolites in animal tissues.

High-Resolution Maps of Mouse Reference Populations

Genetic reference panels are widely used to map complex, quantitative traits in model organisms. We have generated new high-resolution genetic maps of 259 mouse inbred strains from recombinant inbred strain panels (C57BL/6J × DBA/2J, ILS/IbgTejJ × ISS/IbgTejJ, and C57BL/6J × A/J) and chromosome substitution strain panels (C57BL/6J-Chr#<A/J>, C57BL/6J-Chr#<PWD/Ph>, and C57BL/6J-Chr#<MSM/Ms>). We genotyped all samples using the Affymetrix Mouse Diversity Array with an average intermarker spacing of 4.3 kb. The new genetic maps provide increased precision in the localization of recombination breakpoints compared to the previous maps. Although the strains were presumed to be fully inbred, we found residual heterozygosity in 40% of individual mice from five of the six panels. We also identified de novo deletions and duplications, in homozygous or heterozygous state, ranging in size from 21 kb to 8.4 Mb. Almost two-thirds (46 out of 76) of these deletions overlap exons of protein coding genes and may have phenotypic consequences. Twenty-nine putative gene conversions were identified in the chromosome substitution strains. We find that gene conversions are more likely to occur in regions where the homologous chromosomes are more similar. The raw genotyping data and genetic maps of these strain panels are available at http://churchill-lab.jax.org/website/MDA.

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.

CRISPR/Cas9-mediated genome editing in wild-derived mice: generation of tamed wild-derived strains by mutation of the a (nonagouti) gene

Wild-derived mice have contributed to experimental mouse genetics by virtue of their genetic diversity, which may help increase the chance of identifying novel modifier genes responsible for specific phenotypes and diseases. However, gene targeting using wild-derived mice has been unsuccessful because of the unavailability of stable embryonic stem cells. Here, we report that CRISPR/Cas9-mediated gene targeting can be applied to the Japanese wild-derived MSM/Ms strain (Mus musculus molossinus). We targeted the nonagouti (a) gene encoding the agouti protein that is localized in hair and the brain. We obtained three homozygous knockout mice as founders, all showing black coat colour. While homozygous knockout offspring were physiologically indistinguishable from wild-type litter-mates, they showed specific domesticated behaviours: hypoactivity in the dark phase and a decline in the avoidance of a human hand. These phenotypes were consistent over subsequent generations. Our findings support the empirical hypothesis that nonagouti is a domestication-linked gene, the loss of which might repress aggressive behaviour.

Whole Genome Sequence of Two Wild-Derived Mus musculus domesticus Inbred Strains, LEWES/EiJ and ZALENDE/EiJ, with Different Diploid Numbers

Wild-derived mouse inbred strains are becoming increasingly popular for complex traits analysis, evolutionary studies, and systems genetics. Here, we report the whole-genome sequencing of two wild-derived mouse inbred strains, LEWES/EiJ and ZALENDE/EiJ, of Mus musculus domesticus origin. These two inbred strains were selected based on their geographic origin, karyotype, and use in ongoing research. We generated 14× and 18× coverage sequence, respectively, and discovered over 1.1 million novel variants, most of which are private to one of these strains. This report expands the number of wild-derived inbred genomes in the Mus genus from six to eight. The sequence variation can be accessed via an online query tool; variant calls (VCF format) and alignments (BAM format) are available for download from a dedicated ftp site. Finally, the sequencing data have also been stored in a lossless, compressed, and indexed format using the multi-string Burrows-Wheeler transform. All data can be used without restriction.

Heterozygous mutation of Ush1g/Sans in mice causes early-onset progressive hearing loss, which is recovered by reconstituting the strain-specific mutation in Cdh23

Most clinical reports have suggested that patients with congenital profound hearing loss have recessive mutations in deafness genes, whereas dominant alleles are associated with progressive hearing loss (PHL). Jackson shaker (Ush1g^js) is a mouse model of recessive deafness that exhibits congenital profound deafness caused by the homozygous mutation of Ush1g/Sans on chromosome 11. We found that C57BL/6J-Ush1g^js/+ heterozygous mice exhibited early-onset PHL (ePHL) accompanied by progressive degeneration of stereocilia in the cochlear outer hair cells. Interestingly, ePHL did not develop in mutant mice with the C3H/HeN background, thus suggesting that other genetic factors are required for ePHL development. Therefore, we performed classical genetic analyses and found that the occurrence of ePHL in Ush1g^js/+ mice was associated with an interval in chromosome 10 that contains the cadherin 23 gene (Cdh23), which is also responsible for human deafness. To confirm this mutation effect, we generated C57BL/6J-Ush1g^js/+, Cdh23^c.753A/G double-heterozygous mice by using the CRISPR/Cas9-mediated Cdh23^c.753A>G knock-in method. The Cdh23^c.753A/G mice harbored a one-base substitution (A for G), and the homozygous A allele caused moderate hearing loss with aging. Analyses revealed the complete recovery of ePHL and stereocilia degeneration in C57BL/6J-Ush1g^js/+ mice. These results clearly show that the development of ePHL requires at least two mutant alleles of the Ush1g and Cdh23 genes. Our results also suggest that because the SANS and CDH23 proteins form a complex in the stereocilia, the interaction between these proteins may play key roles in the maintenance of stereocilia and the prevention of ePHL.

Body Composition QTLs Identified in Intercross Populations Are Reproducible in Consomic Mouse Strains

Genetic variation contributes to individual differences in obesity, but defining the exact relationships between naturally occurring genotypes and their effects on fatness remains elusive. As a step toward positional cloning of previously identified body composition quantitative trait loci (QTLs) from F2 crosses of mice from the C57BL/6ByJ and 129P3/J inbred strains, we sought to recapture them on a homogenous genetic background of consomic (chromosome substitution) strains. Male and female mice from reciprocal consomic strains originating from the C57BL/6ByJ and 129P3/J strains were bred and measured for body weight, length, and adiposity. Chromosomes 2, 7, and 9 were selected for substitution because previous F2 intercross studies revealed body composition QTLs on these chromosomes. We considered a QTL confirmed if one or both sexes of one or both reciprocal consomic strains differed significantly from the host strain in the expected direction after correction for multiple testing. Using these criteria, we confirmed two of two QTLs for body weight (Bwq5-6), three of three QTLs for body length (Bdln3-5), and three of three QTLs for adiposity (Adip20, Adip26 and Adip27). Overall, this study shows that despite the biological complexity of body size and composition, most QTLs for these traits are preserved when transferred to consomic strains; in addition, studying reciprocal consomic strains of both sexes is useful in assessing the robustness of a particular QTL.

Mapping of Genetic Factors That Elicit Intermale Aggressive Behavior on Mouse Chromosome 15: Intruder Effects and the Complex Genetic Basis

Despite high estimates of the heritability of aggressiveness, the genetic basis for individual differences in aggression remains unclear. Previously, we showed that the wild-derived mouse strain MSM/Ms (MSM) exhibits highly aggressive behaviors, and identified chromosome 15 (Chr 15) as the location of one of the genetic factors behind this escalated aggression by using a panel of consomic strains of MSM in a C57BL/6J (B6) background. To understand the genetic effect of Chr 15 derived from MSM in detail, this study examined the aggressive behavior of a Chr 15 consomic strain towards different types of opponent. Our results showed that both resident and intruder animals had to have the same MSM Chr 15 genotype in order for attack bites to increase and attack latency to be reduced, whereas there was an intruder effect of MSM Chr 15 on tail rattle behavior. To narrow down the region that contains the genetic loci involved in the aggression-eliciting effects on Chr 15, we established a panel of subconsomic strains of MSM Chr 15. Analysis of these strains suggested the existence of multiple genes that enhance and suppress aggressive behavior on Chr 15, and these loci interact in a complex way. Regression analysis successfully identified four genetic loci on Chr 15 that influence attack latency, and one genetic locus that partially elicits aggressive behaviors was narrowed down to a 4.1-Mbp region (from 68.40 Mb to 72.50 Mb) on Chr 15.

NIG_MoG: a mouse genome navigator for exploring intersubspecific genetic polymorphisms

The National Institute of Genetics Mouse Genome database (NIG_MoG; http://molossinus.lab.nig.ac.jp/msmdb/) primarily comprises the whole-genome sequence data of two inbred mouse strains, MSM/Ms and JF1/Ms. These strains were established at NIG and originated from the Japanese subspecies Mus musculus molossinus. NIG_MoG provides visualized genome polymorphism information, browsing single-nucleotide polymorphisms and short insertions and deletions in the genomes of MSM/Ms and JF1/Ms with respect to C57BL/6J (whose genome is predominantly derived from the West European subspecies M. m. domesticus). This allows users, especially wet-lab biologists, to intuitively recognize intersubspecific genome divergence in these mouse strains using visual data. The database also supports the in silico screening of bacterial artificial chromosome (BAC) clones that contain genomic DNA from MSM/Ms and the standard classical laboratory strain C57BL/6N. NIG_MoG is thus a valuable navigator for exploring mouse genome polymorphisms and BAC clones that are useful for studies of gene function and regulation based on intersubspecific genome divergence.

Contrasting genetic architectures in different mouse reference populations used for studying complex traits

Quantitative trait loci (QTLs) are being used to study genetic networks, protein functions, and systems properties that underlie phenotypic variation and disease risk in humans, model organisms, agricultural species, and natural populations. The challenges are many, beginning with the seemingly simple tasks of mapping QTLs and identifying their underlying genetic determinants. Various specialized resources have been developed to study complex traits in many model organisms. In the mouse, remarkably different pictures of genetic architectures are emerging. Chromosome Substitution Strains (CSSs) reveal many QTLs, large phenotypic effects, pervasive epistasis, and readily identified genetic variants. In contrast, other resources as well as genome-wide association studies (GWAS) in humans and other species reveal genetic architectures dominated with a relatively modest number of QTLs that have small individual and combined phenotypic effects. These contrasting architectures are the result of intrinsic differences in the study designs underlying different resources. The CSSs examine context-dependent phenotypic effects independently among individual genotypes, whereas with GWAS and other mouse resources, the average effect of each QTL is assessed among many individuals with heterogeneous genetic backgrounds. We argue that variation of genetic architectures among individuals is as important as population averages. Each of these important resources has particular merits and specific applications for these individual and population perspectives. Collectively, these resources together with high-throughput genotyping, sequencing and genetic engineering technologies, and information repositories highlight the power of the mouse for genetic, functional, and systems studies of complex traits and disease models.

Chromosomal assignment of quantitative trait loci influencing baseline circulating total cholesterol level in male laboratory mice: report of a consomic strain survey and comparison with published results

Background

An important risk for atherosclerosis is a low level of HDL cholesterol. Baseline HDL cholesterol is under complex genetic and environmental control. Here we report on results of male mice from a consomic strain survey and the parental inbred strains for baseline circulating total cholesterol concentration, which is almost the same as HDL cholesterol in chow fed mice. The consomic strains have been derived from C57BL/6J (host strain) and A/J (donor strain) inbred lines. The work contributes to the value of the mouse as an animal model for studying the genetic background of differences in baseline circulating total and HDL cholesterol levels.

Results

The consomic strain survey suggested that mouse chromosomes 1, 7, 9, 14, 16, 17, 19, X, and Y contained at least one quantitative trait locus that is involved in baseline circulating total cholesterol concentration. All consomic lines, for which there is evidence that the substituted chromosome contains a quantitative trait locus, increase compared to the host strain baseline circulating total cholesterol concentration. Since there is evidence that ‘body weight’, ‘age at blood sampling’, ‘time of the day blood was collected’, and ‘season’ influence this phenotype, additional statistical analyses (with these variables as covariates) were performed. Now there is only evidence for quantitative trait loci on chromosomes 1, 8, 12, and Y. Taken the present results together with previous consomic strain surveys there is evidence that all mouse chromosomes carry quantitative trait loci that control baseline circulating total cholesterol levels. There was however little agreement between the present consomic strain results and previous sets of data. This might be explained by seasonal effects and differences in methodological variables such as age of the mice, fasting versus non-fasting, percentage of dietary fat, unanesthetized versus anesthetized mice, and the daily light–dark cycle.

Conclusions

The present findings, when compared with previous consomic strain surveys, clearly illustrate the complexity of the genetic-environmental architecture for the regulation of baseline circulating total cholesterol levels in mice. Different data can be obtained from different labs and it underscores that animal geneticists should present as accurate a picture as possible of the laboratory mouse’s environment.

Devising assisted reproductive technologies for wild-derived strains of mice: 37 strains from five subspecies of Mus musculus

Wild-derived mice have long offered invaluable experimental models for mouse genetics because of their high evolutionary divergence from laboratory mice. A number of wild-derived strains are available from the RIKEN BioResource Center (BRC), but they have been maintained as living stocks because of the unavailability of assisted reproductive technology (ART). In this study, we sought to devise ART for 37 wild-derived strains from five subspecies of Mus musculus maintained at the BRC. Superovulation of females was effective (more than 15 oocytes per female) for 34 out of 37 strains by treatment with either equine chorionic gonadotropin or anti-inhibin serum, depending on their genetic background (subspecies). The collected oocytes could be fertilized in vitro at mean rates of 79.0% and 54.6% by the optimized protocol using fresh or frozen-thawed spermatozoa, respectively. They were cryopreserved at the 2-cell stage by vitrification with an ethylene glycol-based solution. In total, 94.6% of cryopreserved embryos survived the vitrification procedure and restored their normal morphology after warming. A conventional embryo transfer protocol could be applied to 25 out of the 35 strains tested. In the remaining 10 strains, live offspring could be obtained by a modified embryo transfer protocol using cyclosporin A treatment and co-transfer of ICR (laboratory mouse strain) embryos. Thus, ART for 37 wild-derived strains was devised successfully and is now routinely used for their preservation and transportation. The information provided here might facilitate broader use and wider distribution of wild-derived mice for biomedical research.

Genetic determinants of atherosclerosis, obesity, and energy balance in consomic mice

Metabolic diseases such as obesity and atherosclerosis result from complex interactions between environmental factors and genetic variants. A panel of chromosome substitution strains (CSSs) was developed to characterize genetic and dietary factors contributing to metabolic diseases and other biological traits and biomedical conditions. Our goal here was to identify quantitative trait loci (QTLs) contributing to obesity, energy expenditure, and atherosclerosis. Parental strains C57BL/6 and A/J together with a panel of 21 CSSs derived from these progenitors were subjected to chronic feeding of rodent chow and atherosclerotic (females) or diabetogenic (males) test diets, and evaluated for a variety of metabolic phenotypes including several traits unique to this report, namely fat pad weights, energy balance, and atherosclerosis. A total of 297 QTLs across 35 traits were discovered, two of which provided significant protection from atherosclerosis, and several dozen QTLs modulated body weight, body composition, and circulating lipid levels in females and males. While several QTLs confirmed previous reports, most QTLs were novel. Finally, we applied the CSS quantitative genetic approach to energy balance, and identified three novel QTLs controlling energy expenditure and one QTL modulating food intake. Overall, we identified many new QTLs and phenotyped several novel traits in this mouse model of diet-induced metabolic diseases.

Genetic mapping of escalated aggression in wild-derived mouse strain MSM/Ms: association with serotonin-related genes

The Japanese wild-derived mouse strain MSM/Ms (MSM) retains a wide range of traits related to behavioral wildness, including high levels of emotionality and avoidance of humans. In this study, we observed that MSM showed a markedly higher level of aggression than the standard laboratory strain C57BL/6J. Whereas almost all MSM males showed high frequencies of attack bites and pursuit in the resident-intruder test, only a few C57BL/6J males showed aggressive behaviors, with these behaviors observed at only a low frequency. Sexually mature MSM males in their home cages killed their littermates, or sometimes female pair-mates. To study the genetic and neurobiological mechanisms that underlie the escalated aggression observed in MSM mice, we analyzed reciprocal F1 crosses and five consomic strains of MSM (Chr 4, 13, 15, X and Y) against the background of C57BL/6J. We identified two chromosomes, Chr 4 and Chr 15, which were involved in the heightened aggression observed in MSM. These chromosomes had different effects on aggression: whereas MSM Chr 15 increased agitation and initiation of aggressive events, MSM Chr 4 induced a maladaptive level of aggressive behavior. Expression analysis of mRNAs of serotonin receptors, serotonin transporter and Tph2, an enzyme involved in serotonin synthesis in seven brain areas, indicated several differences among MSM, C57BL/6J, and their consomic strains. We found that Tph2 expression in the midbrain was increased in the Chr 4 consomic strain, as well as in MSM, and that there was a strong positive genetic correlation between aggressive behavior and Tph2 expression at the mRNA level. Therefore, it is possible that increased expression of the Tph2 gene is related to escalated aggression observed in MSM.

Segregation of a QTL cluster for home-cage activity using a new mapping method based on regression analysis of congenic mouse strains

Recent genetic studies have shown that genetic loci with significant effects in whole-genome quantitative trait loci (QTL) analyses were lost or weakened in congenic strains. Characterisation of the genetic basis of this attenuated QTL effect is important to our understanding of the genetic mechanisms of complex traits. We previously found that a consomic strain, B6-Chr6C^MSM, which carries chromosome 6 of a wild-derived strain MSM/Ms on the genetic background of C57BL/6J, exhibited lower home-cage activity than C57BL/6J. In the present study, we conducted a composite interval QTL analysis using the F2 mice derived from a cross between C57BL/6J and B6-Chr6C^MSM. We found one QTL peak that spans 17.6 Mbp of chromosome 6. A subconsomic strain that covers the entire QTL region also showed lower home-cage activity at the same level as the consomic strain. We developed 15 congenic strains, each of which carries a shorter MSM/Ms-derived chromosomal segment from the subconsomic strain. Given that the results of home-cage activity tests on the congenic strains cannot be explained by a simple single-gene model, we applied regression analysis to segregate the multiple genetic loci. The results revealed three loci (loci 1–3) that have the effect of reducing home-cage activity and one locus (locus 4) that increases activity. We also found that the combination of loci 3 and 4 cancels out the effects of the congenic strains, which indicates the existence of a genetic mechanism related to the loss of QTLs.

Influence analysis in quantitative trait loci detection

This paper presents systematic methods for the detection of influential individuals that affect the log odds (LOD) score curve. We derive general formulas of influence functions for profile likelihoods and introduce them into two standard quantitative trait locus detection methods—the interval mapping method and single marker analysis. Besides influence analysis on specific LOD scores, we also develop influence analysis methods on the shape of the LOD score curves. A simulation-based method is proposed to assess the significance of the influence of the individuals. These methods are shown useful in the influence analysis of a real dataset of an experimental population from an F₂ mouse cross. By receiver operating characteristic analysis, we confirm that the proposed methods show better performance than existing diagnostics.

Evolutionarily diverged regulation of X-chromosomal genes as a primal event in mouse reproductive isolation

Improper gene regulation is implicated in reproductive isolation, but its genetic and molecular bases are unknown. We previously reported that a mouse inter-subspecific X chromosome substitution strain shows reproductive isolation characterized by male-specific sterility due to disruption of meiotic entry in spermatogenesis. Here, we conducted comprehensive transcriptional profiling of the testicular cells of this strain by microarray. The results clearly revealed gross misregulation of gene expression in the substituted donor X chromosome. Such misregulation occurred prior to detectable spermatogenetic impairment, suggesting that it is a primal event in reproductive isolation. The misregulation of X-linked genes showed asymmetry; more genes were disproportionally downregulated rather than upregulated. Furthermore, this misregulation subsequently resulted in perturbation of global transcriptional regulation of autosomal genes, probably by cascading deleterious effects. Remarkably, this transcriptional misregulation was substantially restored by introduction of chromosome 1 from the same donor strain as the X chromosome. This finding implies that one of regulatory genes acting in trans for X-linked target genes is located on chromosome 1. This study collectively suggests that regulatory incompatibility is a major cause of reproductive isolation in the X chromosome substitution strain.

Reproductive isolation characterized by male sterility and decreased viability is important for speciation, because it suppresses free genetic exchange between two diverged populations and accelerates the genetic divergence. One of the reproductive isolation phenomena, hybrid sterility (sterility in hybrid animals), is possibly caused by deleterious interactions between diverged genetic factors brought by two distinct populations. The polymorphism not only in protein-coding sequences but also in transcriptional regulatory sequences can cause the genetic incompatibility in hybrid animals. However, the precise genetic mechanisms of hybrid sterility are mostly unknown. Here, we report that the expression of X-linked genes derived from one mouse subspecies was largely misregulated in the genetic background of another subspecies. The misregulated expression of the X-linked genes subsequently affected the global expression of autosomal genes. The results collectively indicate that hybrid sterility between the two mouse subspecies is caused by misregulation of gene expression due to genetic incompatibility in the transcriptional regulatory circuitry. Such genetic incompatibility in transcriptional regulation likely underlies reproductive isolation in general.

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

Commonly used classical inbred mouse strains have mosaic genomes with sequences from different subspecific origins. Their genomes are derived predominantly from the Western European subspecies Mus musculus domesticus, with the remaining sequences derived mostly from the Japanese subspecies Mus musculus molossinus. However, it remains unknown how this intersubspecific genome introgression occurred during the establishment of classical inbred strains. In this study, we resequenced the genomes of two M. m. molossinus-derived inbred strains, MSM/Ms and JF1/Ms. MSM/Ms originated from Japanese wild mice, and the ancestry of JF1/Ms was originally found in Europe and then transferred to Japan. We compared the characteristics of these sequences to those of the C57BL/6J reference sequence and the recent data sets from the resequencing of 17 inbred strains in the Mouse Genome Project (MGP), and the results unequivocally show that genome introgression from M. m. molossinus into M. m. domesticus provided the primary framework for the mosaic genomes of classical inbred strains. Furthermore, the genomes of C57BL/6J and other classical inbred strains have long consecutive segments with extremely high similarity (>99.998%) to the JF1/Ms strain. In the early 20th century, Japanese waltzing mice with a morphological phenotype resembling that of JF1/Ms mice were often crossed with European fancy mice for early studies of "Mendelism," which suggests that the ancestor of the extant JF1/Ms strain provided the origin of the M. m. molossinus genome in classical inbred strains and largely contributed to its intersubspecific genome diversity.

Identification of both copy number variation-type and constant-type core elements in a large segmental duplication region of the mouse genome

Background

Copy number variation (CNV), an important source of diversity in genomic structure, is frequently found in clusters called CNV regions (CNVRs). CNVRs are strongly associated with segmental duplications (SDs), but the composition of these complex repetitive structures remains unclear.

Results

We conducted self-comparative-plot analysis of all mouse chromosomes using the high-speed and large-scale-homology search algorithm SHEAP. For eight chromosomes, we identified various types of large SD as tartan-checked patterns within the self-comparative plots. A complex arrangement of diagonal split lines in the self-comparative-plots indicated the presence of large homologous repetitive sequences. We focused on one SD on chromosome 13 (SD13M), and developed SHEPHERD, a stepwise ab initio method, to extract longer repetitive elements and to characterize repetitive structures in this region. Analysis using SHEPHERD showed the existence of 60 core elements, which were expected to be the basic units that form SDs within the repetitive structure of SD13M. The demonstration that sequences homologous to the core elements (>70% homology) covered approximately 90% of the SD13M region indicated that our method can characterize the repetitive structure of SD13M effectively. Core elements were composed largely of fragmented repeats of a previously identified type, such as long interspersed nuclear elements (LINEs), together with partial genic regions. Comparative genome hybridization array analysis showed that whereas 42 core elements were components of CNVR that varied among mouse strains, 8 did not vary among strains (constant type), and the status of the others could not be determined. The CNV-type core elements contained significantly larger proportions of long terminal repeat (LTR) types of retrotransposon than the constant-type core elements, which had no CNV. The higher divergence rates observed in the CNV-type core elements than in the constant type indicate that the CNV-type core elements have a longer evolutionary history than constant-type core elements in SD13M.

Conclusions

Our methodology for the identification of repetitive core sequences simplifies characterization of the structures of large SDs and detailed analysis of CNV. The results of detailed structural and quantitative analyses in this study might help to elucidate the biological role of one of the SDs on chromosome 13.

Genomic mixing to elucidate the genetic system of complex traits

Understanding the genetic basis of complex traits has become one of the major issues in genetics, but recent advances in this field are still faced with a difficulty, the so-called "missing heritability." It is speculated that missing heritability mainly stems from a large number of variants of small effect that are poorly detected by currently available methods. In order to overcome this problem, many recent genetic studies of complex traits have actively used outbred stocks of mice. However, most of the available outbred stocks have a limited amount of genetic variation, because many stocks originate from Swiss mouse colonies. We have repeatedly shown that wild-derived strains are a useful mouse resource since there is a large genetic diversity among these strains. Here, we give an overview of mouse resources produced by crossing different founder mice. Finally, we propose an advantage of new attempts to conduct selective breeding using heterogeneous stocks created by mixing genomes of wild-derived inbred strains of mice when studying complex traits.

Neurogenetics of aggressive behavior: studies in rodents

Aggressive behavior is observed in many animal species, such as insects, fish, lizards, frogs, and most mammals including humans. This wide range of conservation underscores the importance of aggressive behavior in the animals' survival and fitness, and the likely heritability of this behavior. Although typical patterns of aggressive behavior differ between species, there are several concordances in the neurobiology of aggression among rodents, primates, and humans. Studies with rodent models may eventually help us to understand the neurogenetic architecture of aggression in humans. However, it is important to recognize the difference between the ecological and ethological significance of aggressive behavior (species-typical aggression) and maladaptive violence (escalated aggression) when applying the findings of aggression research using animal models to human or veterinary medicine. Well-studied rodent models for aggressive behavior in the laboratory setting include the mouse (Mus musculus), rat (Rattus norvegicus), hamster (Mesocricetus auratus), and prairie vole (Microtus ochrogaster). The neural circuits of rodent aggression have been gradually elucidated by several techniques, e.g., immunohistochemistry of immediate-early gene (c-Fos) expression, intracranial drug microinjection, in vivo microdialysis, and optogenetics techniques. Also, evidence accumulated from the analysis of gene-knockout mice shows the involvement of several genes in aggression. Here, we review the brain circuits that have been implicated in aggression, such as the hypothalamus, prefrontal cortex (PFC), dorsal raphe nucleus (DRN), nucleus accumbens (NAc), and olfactory system. We then discuss the roles of glutamate and γ-aminobutyric acid (GABA), excitatory and inhibitory amino acids in the brain, as well as their receptors, in controlling aggressive behavior, focusing mainly on recent findings. At the end of this chapter, we discuss how genes can be identified that underlie individual differences in aggression, using the so-called forward genetics approach.

Establishment of consomic strains derived from A/J and SM/J mice for genetic analysis of complex traits

Consomic strains, in which one chromosome is derived from a donor strain and the other chromosomes are derived from the recipient strain, provide a powerful tool for the dissection of complex genetic traits. In this study we established ten consomic strains (A-2(SM), A-6(SM), A-11(SM), A-12(SM), A-13(SM), A-15(SM), A-17(SM), A-18(SM), A-19(SM), A-Y(SM)) using the SM/J strain as the donor and the A/J strain as the recipient; these are the parental strains of a set of SMXA recombinant inbred (RI) strains that we had developed previously. We analyzed body weights and blood lipid levels in the consomic and parental strains. The mean values for each trait showed a continuous range of variation in the consomic strains suggesting that they are controlled by multiple genes. We previously identified suggestive QTLs for body weight on chromosome 6 in SMXA RI strains and (SM/J × A/J)F(2) mice. The observation that the A-6(SM) consomic strain had a significantly lower mean body weight than the A/J strain supports the presence of this QTL on chromosome 6. Similarly, the higher blood triglyceride level in the A-11(SM) strain shows the existence of a previously mapped QTL on chromosome 11, and the A-12(SM) strain provides evidence of a QTL for blood total cholesterol level on chromosome 12. These consomic strains, along with the previously developed set of SMXA RI strains from A/J and SM/J mice, offer an invaluable and powerful resource for the analysis of complex genetic traits in mice.

Hybrid mouse diversity panel: a panel of inbred mouse strains suitable for analysis of complex genetic traits

We have developed an association-based approach using classical inbred strains of mice in which we correct for population structure, which is very extensive in mice, using an efficient mixed-model algorithm. Our approach includes inbred parental strains as well as recombinant inbred strains in order to capture loci with effect sizes typical of complex traits in mice (in the range of 5% of total trait variance). Over the last few years, we have typed the hybrid mouse diversity panel (HMDP) strains for a variety of clinical traits as well as intermediate phenotypes and have shown that the HMDP has sufficient power to map genes for highly complex traits with resolution that is in most cases less than a megabase. In this essay, we review our experience with the HMDP, describe various ongoing projects, and discuss how the HMDP may fit into the larger picture of common diseases and different approaches.

Chromosome substitution strains: gene discovery, functional analysis, and systems studies

Laboratory mice are valuable in biomedical research in part because of the extraordinary diversity of genetic resources that are available for studies of complex genetic traits and as models for human biology and disease. Chromosome substitution strains (CSSs) are important in this resource portfolio because of their demonstrated use for gene discovery, genetic and epigenetic studies, functional characterizations, and systems analysis. CSSs are made by replacing a single chromosome in a host strain with the corresponding chromosome from a donor strain. A complete CSS panel involves a total of 22 engineered inbred strains, one for each of the 19 autosomes, one each for the X and Y chromosomes, and one for mitochondria. A genome survey simply involves comparing each phenotype for each of the CSSs with the phenotypes of the host strain. The CSS panels that are available for laboratory mice have been used to dissect a remarkable variety of phenotypes and to characterize an impressive array of disease models. These surveys have revealed considerable phenotypic diversity even among closely related progenitor strains, evidence for strong epistasis and for heritable epigenetic changes. Perhaps most importantly, and presumably because of their unique genetic constitution, CSSs, and congenic strains derived from them, the genetic variants underlying quantitative trait loci (QTLs) are readily identified and functionally characterized. Together these studies show that CSSs are important resource for laboratory mice.

Efficient production of offspring from Japanese wild-derived strains of mice (Mus musculus molossinus) by improved assisted reproductive technologies

Because the genetic diversity of the laboratory mouse (Mus musculus) is very limited, wild-derived strains from this genus could provide invaluable experimental models for studies of mouse genetics and epigenetics such as quantitative trait locus analysis. However, such strains generally show poor reproductive performance under conventional husbandry conditions, so their use for large-scale analyses has been limited. This study was undertaken to devise assisted reproductive technologies (ARTs) for the efficient production of offspring in two wild-derived strains, MSM/Ms and JF1/Ms (Mus musculus molossinus). First, as females of these strains are poor responders to equine chorionic gonadotropin (eCG) stimulation, we examined the efficiency of superovulation by injecting anti-inhibin serum followed by human chorionic gonadotropin (hCG). Approximately four to six times more oocytes were ovulated than with eCG-hCG treatment in both strains, reaching ∼25-30 oocytes per female. Consequently, the procedures for in vitro fertilization using these superovulated oocytes and cryopreservation of embryos and spermatozoa could be optimized for both of the wild-derived strains. However, MSM/Ms embryos but not JF1/Ms embryos failed to develop to term after embryo transfer because of intrauterine death at mid to late gestation. We were able to overcome this obstacle by cotransfer of these embryos with those from laboratory strains combined with treatment of recipient females with an immunosuppressant (cyclosporin A). Thus, a series of ARTs essential for efficient production and preservation of the wild-derived strains were successfully devised. These technologies will facilitate systematic studies of mouse genetics and epigenetics using a wider range of genetic diversity than currently available in the genus Mus.

Genetic divergence and the genetic architecture of complex traits in chromosome substitution strains of mice

Background

The genetic architecture of complex traits strongly influences the consequences of inherited mutations, genetic engineering, environmental and genetic perturbations, and natural and artificial selection. But because most studies are under-powered, the picture of complex traits is often incomplete. Chromosome substitution strains (CSSs) are a unique paradigm for these genome surveys because they enable statistically independent, powerful tests for the phenotypic effects of each chromosome on a uniform inbred genetic background. A previous CSS survey in mice and rats revealed many complex trait genes (QTLs), large phenotypic effects, extensive epistasis, as well as systems properties such as strongly directional phenotypic changes and genetically-determined limits on the range of phenotypic variation. However, the unusually close genetic relation between the CSS progenitor strains in that study raised questions about the impact of genetic divergence: would greater divergence between progenitor strains, with the corresponding changes in gene regulation and protein function, lead to significantly more distinctive phenotypic features, or alternatively would epistasis and systems constraints, which are pervasive in CSSs, limit the range of phenotypic variation regardless of the extent of DNA sequence variation?

Results

We analyzed results for an extensive survey of traits in two new panels of CSSs where the donor strains were derived from inbred strains with more distant origins and discovered a strong similarity in genetic and systems properties among the three CSS panels, regardless of divergence time.

Conclusion

Our results argue that DNA sequence differences between host and donor strains did not substantially affect the architecture of complex traits, and suggest instead that strong epistasis buffered the phenotypic effects of genetic divergence, thereby constraining the range of phenotypic variation.

Dissection of complex adult traits in a mouse synthetic population

Finding the causative genetic variations that underlie complex adult traits is a significant experimental challenge. The unbiased search strategy of genome-wide association (GWAS) has been used extensively in recent human population studies. These efforts, however, typically find only a minor fraction of the genetic loci that are predicted to affect variation. As an experimental model for the analysis of adult polygenic traits, we measured a mouse population for multiple phenotypes and conducted a genome-wide search for effector loci. Complex adult phenotypes, related to body size and bone structure, were measured as component phenotypes, and each subphenotype was associated with a genomic spectrum of candidate effector loci. The strategy successfully detected several loci for the phenotypes, at genome-wide significance, using a single, modest-sized population (N = 505). The effector loci each explain 2%–10% of the measured trait variation and, taken together, the loci can account for over 25% of a trait's total population variation. A replicate population (N = 378) was used to confirm initially observed loci for one trait (femur length), and, when the two groups were merged, the combined population demonstrated increased power to detect loci. In contrast to human population studies, our mouse genome-wide searches find loci that individually explain a larger fraction of the observed variation. Also, the additive effects of our detected mouse loci more closely match the predicted genetic component of variation. The genetic loci discovered are logical candidates for components of the genetic networks having evolutionary conservation with human biology.