A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT.

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

The ENCODE mouse postnatal developmental time course identifies regulatory programs of cell types and cell states.

Postnatal genomic regulation significantly influences tissue and organ maturation but is under-studied relative to existing genomic catalogs of adult tissues or prenatal development in mouse. The ENCODE4 consortium generated the first comprehensive single-nucleus resource of postnatal regulatory events across a diverse set of mouse tissues. The collection spans seven postnatal time points, mirroring human development from childhood to adulthood, and encompasses five core tissues. We identified 30 cell types, further subdivided into 69 subtypes and cell states across adrenal gland, left cerebral cortex, hippocampus, heart, and gastrocnemius muscle. Our annotations cover both known and novel cell differentiation dynamics ranging from early hippocampal neurogenesis to a new sex-specific adrenal gland population during puberty. We used an ensemble Latent Dirichlet Allocation strategy with a curated vocabulary of 2,701 regulatory genes to identify regulatory "topics," each of which is a gene vector, linked to cell type differentiation, subtype specialization, and transitions between cell states. We find recurrent regulatory topics in tissue-resident macrophages, neural cell types, endothelial cells across multiple tissues, and cycling cells of the adrenal gland and heart. Cell-type-specific topics are enriched in transcription factors and microRNA host genes, while chromatin regulators dominate mitosis topics. Corresponding chromatin accessibility data reveal dynamic and sex-specific regulatory elements, with enriched motifs matching transcription factors in regulatory topics. Together, these analyses identify both tissue-specific and common regulatory programs in postnatal development across multiple tissues through the lens of the factors regulating transcription.

Unraveling the genetics of arsenic toxicity with cellular morphology QTL.

The health risks that arise from environmental exposures vary widely within and across human populations, and these differences are largely determined by genetic variation and gene-by-environment (gene-environment) interactions. However, risk assessment in laboratory mice typically involves isogenic strains and therefore, does not account for these known genetic effects. In this context, genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening because they provide a way to introduce genetic variation in risk assessment without increasing animal use. Cell lines from genetic reference populations of laboratory mice offer genetic diversity, power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.

Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice.

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

An in vitro neurogenetics platform for precision disease modeling in the mouse.

The power and scope of disease modeling can be markedly enhanced through the incorporation of broad genetic diversity. The introduction of pathogenic mutations into a single inbred mouse strain sometimes fails to mimic human disease. We describe a cross-species precision disease modeling platform that exploits mouse genetic diversity to bridge cell-based modeling with whole organism analysis. We developed a universal protocol that permitted robust and reproducible neural differentiation of genetically diverse human and mouse pluripotent stem cell lines and then carried out a proof-of-concept study of the neurodevelopmental gene DYRK1A. Results in vitro reliably predicted the effects of genetic background on Dyrk1a loss-of-function phenotypes in vivo. Transcriptomic comparison of responsive and unresponsive strains identified molecular pathways conferring sensitivity or resilience to Dyrk1a1A loss and highlighted differential messenger RNA isoform usage as an important determinant of response. This cross-species strategy provides a powerful tool in the functional analysis of candidate disease variants identified through human genetic studies.

Identification of robust cellular programs using reproducible LDA that impact sex-specific disease progression in different genotypes of a mouse model of AD.

The gene expression profiles of distinct cell types reflect complex genomic interactions among multiple simultaneous biological processes within each cell that can be altered by disease progression as well as genetic background. The identification of these active cellular programs is an open challenge in the analysis of single-cell RNA-seq data. Latent Dirichlet Allocation (LDA) is a generative method used to identify recurring patterns in counts data, commonly referred to as topics that can be used to interpret the state of each cell. However, LDA's interpretability is hindered by several key factors including the hyperparameter selection of the number of topics as well as the variability in topic definitions due to random initialization. We developed Topyfic, a Reproducible LDA (rLDA) package, to accurately infer the identity and activity of cellular programs in single-cell data, providing insights into the relative contributions of each program in individual cells. We apply Topyfic to brain single-cell and single-nucleus datasets of two 5xFAD mouse models of Alzheimer's disease crossed with C57BL6/J or CAST/EiJ mice to identify distinct cell types and states in different cell types such as microglia. We find that 8-month 5xFAD/Cast F1 males show higher level of microglial activation than matching 5xFAD/BL6 F1 males, whereas female mice show similar levels of microglial activation. We show that regulatory genes such as TFs, microRNA host genes, and chromatin regulatory genes alone capture cell types and cell states. Our study highlights how topic modeling with a limited vocabulary of regulatory genes can identify gene expression programs in single-cell data in order to quantify similar and divergent cell states in distinct genotypes.

Cell morphology QTL reveal gene by environment interactions in a genetically diverse cell population.

Genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening as they offer power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.

Spindle assembly checkpoint insensitivity allows meiosis-II despite chromosomal defects in aged eggs.

Chromosome segregation errors in mammalian oocyte meiosis lead to developmentally compromised aneuploid embryos and become more common with advancing maternal age. Known contributors include age-related chromosome cohesion loss and spindle assembly checkpoint (SAC) fallibility in meiosis-I. But how effective the SAC is in meiosis-II and how this might contribute to age-related aneuploidy is unknown. Here, we developed genetic and pharmacological approaches to directly address the function of the SAC in meiosis-II. We show that the SAC is insensitive in meiosis-II oocytes and that as a result misaligned chromosomes are randomly segregated. Whilst SAC ineffectiveness in meiosis-II is not age-related, it becomes most prejudicial in oocytes from older females because chromosomes that prematurely separate by age-related cohesion loss become misaligned in meiosis-II. We show that in the absence of a robust SAC in meiosis-II these age-related misaligned chromatids are missegregated and lead to aneuploidy. Our data demonstrate that the SAC fails to prevent cell division in the presence of misaligned chromosomes in oocyte meiosis-II, which explains how age-related cohesion loss can give rise to aneuploid embryos.

The spontaneous mouse mutant low set ears (Lse) is caused by tandem duplication of Fgf3 and Fgf4.

The external ear develops from an organized convergence of ventrally migrating neural crest cells into the first and second branchial arches. Defects in external ear position are often symptomatic of complex syndromes such as Apert, Treacher-Collins, and Crouzon Syndrome. The low set ears (Lse) spontaneous mouse mutant is characterized by the dominant inheritance of a ventrally shifted external ear position and an abnormal external auditory meatus (EAM). We identified the causative mutation as a 148 Kb tandem duplication on Chromosome 7, which includes the entire coding sequences of Fgf3 and Fgf4. Duplications of FGF3 and FGF4 occur in 11q duplication syndrome in humans and are associated with craniofacial anomalies, among other features. Intercrosses of Lse-affected mice revealed perinatal lethality in homozygotes, and Lse/Lse embryos display additional phenotypes including polydactyly, abnormal eye morphology, and cleft secondary palate. The duplication results in increased Fgf3 and Fgf4 expression in the branchial arches and additional discrete domains in the developing embryo. This ectopic overexpression resulted in functional FGF signaling, demonstrated by increased Spry2 and Etv5 expression in overlapping domains of the developing arches. Finally, a genetic interaction between Fgf3/4 overexpression and Twist1, a regulator of skull suture development, resulted in perinatal lethality, cleft palate, and polydactyly in compound heterozygotes. These data indicate a role for Fgf3 and Fgf4 in external ear and palate development and provide a novel mouse model for further interrogation of the biological consequences of human FGF3/4 duplication.

An allelic series of spontaneous Rorb mutant mice exhibit a gait phenotype, changes in retina morphology and behavior, and gene expression signatures associated with the unfolded protein response.

The Retinoid-related orphan receptor beta (RORß) gene encodes a developmental transcription factor and has 2 predominant isoforms created through alternative first exon usage; one specific to the retina and another present more broadly in the central nervous system, particularly regions involved in sensory processing. RORß belongs to the nuclear receptor family and plays important roles in cell fate specification in the retina and cortical layer formation. In mice, loss of RORß causes disorganized retina layers, postnatal degeneration, and production of immature cone photoreceptors. Hyperflexion or "high-stepping" of rear limbs caused by reduced presynaptic inhibition by Rorb-expressing inhibitory interneurons of the spinal cord is evident in RORß-deficient mice. RORß variants in patients are associated with susceptibility to various neurodevelopmental conditions, primarily generalized epilepsies, but including intellectual disability, bipolar, and autism spectrum disorders. The mechanisms by which RORß variants confer susceptibility to these neurodevelopmental disorders are unknown but may involve aberrant neural circuit formation and hyperexcitability during development. Here we report an allelic series in 5 strains of spontaneous Rorb mutant mice with a high-stepping gait phenotype. We show retinal abnormalities in a subset of these mutants and demonstrate significant differences in various behavioral phenotypes related to cognition. Gene expression analyses in all 5 mutants reveal a shared over-representation of the unfolded protein response and pathways related to endoplasmic reticulum stress, suggesting a possible mechanism of susceptibility relevant to patients.

The ENCODE4 long-read RNA-seq collection reveals distinct classes of transcript structure diversity.

The majority of mammalian genes encode multiple transcript isoforms that result from differential promoter use, changes in exonic splicing, and alternative 3' end choice. Detecting and quantifying transcript isoforms across tissues, cell types, and species has been extremely challenging because transcripts are much longer than the short reads normally used for RNA-seq. By contrast, long-read RNA-seq (LR-RNA-seq) gives the complete structure of most transcripts. We sequenced 264 LR-RNA-seq PacBio libraries totaling over 1 billion circular consensus reads (CCS) for 81 unique human and mouse samples. We detect at least one full-length transcript from 87.7% of annotated human protein coding genes and a total of 200,000 full-length transcripts, 40% of which have novel exon junction chains. To capture and compute on the three sources of transcript structure diversity, we introduce a gene and transcript annotation framework that uses triplets representing the transcript start site, exon junction chain, and transcript end site of each transcript. Using triplets in a simplex representation demonstrates how promoter selection, splice pattern, and 3' processing are deployed across human tissues, with nearly half of multi-transcript protein coding genes showing a clear bias toward one of the three diversity mechanisms. Evaluated across samples, the predominantly expressed transcript changes for 74% of protein coding genes. In evolution, the human and mouse transcriptomes are globally similar in types of transcript structure diversity, yet among individual orthologous gene pairs, more than half (57.8%) show substantial differences in mechanism of diversification in matching tissues. This initial large-scale survey of human and mouse long-read transcriptomes provides a foundation for further analyses of alternative transcript usage, and is complemented by short-read and microRNA data on the same samples and by epigenome data elsewhere in the ENCODE4 collection.

Resolution of structural variation in diverse mouse genomes reveals chromatin remodeling due to transposable elements.

Diverse inbred mouse strains are important biomedical research models, yet genome characterization of many strains is fundamentally lacking in comparison with humans. In particular, catalogs of structural variants (SVs) (variants ≥ 50 bp) are incomplete, limiting the discovery of causative alleles for phenotypic variation. Here, we resolve genome-wide SVs in 20 genetically distinct inbred mice with long-read sequencing. We report 413,758 site-specific SVs affecting 13% (356 Mbp) of the mouse reference assembly, including 510 previously unannotated coding variants. We substantially improve the Mus musculus transposable element (TE) callset, and we find that TEs comprise 39% of SVs and account for 75% of altered bases. We further utilize this callset to investigate how TE heterogeneity affects mouse embryonic stem cells and find multiple TE classes that influence chromatin accessibility. Our work provides a comprehensive analysis of SVs found in diverse mouse genomes and illustrates the role of TEs in epigenetic differences.

Genetic dissection of the pluripotent proteome through multi-omics data integration.

Genetic background drives phenotypic variability in pluripotent stem cells (PSCs). Most studies to date have used transcript abundance as the primary molecular readout of cell state in PSCs. We performed a comprehensive proteogenomics analysis of 190 genetically diverse mouse embryonic stem cell (mESC) lines. The quantitative proteome is highly variable across lines, and we identified pluripotency-associated pathways that were differentially activated in the proteomics data that were not evident in transcriptome data from the same lines. Integration of protein abundance to transcript levels and chromatin accessibility revealed broad co-variation across molecular layers as well as shared and unique drivers of quantitative variation in pluripotency-associated pathways. Quantitative trait locus (QTL) mapping localized the drivers of these multi-omic signatures to genomic hotspots. This study reveals post-transcriptional mechanisms and genetic interactions that underlie quantitative variability in the pluripotent proteome and provides a regulatory map for mESCs that can provide a basis for future mechanistic studies.

Trisomy 21 induces pericentrosomal crowding delaying primary ciliogenesis and mouse cerebellar development.

Trisomy 21, the genetic cause of Down syndrome, disrupts primary cilia formation and function, in part through elevated Pericentrin, a centrosome protein encoded on chromosome 21. Yet how trisomy 21 and elevated Pericentrin disrupt cilia-related molecules and pathways, and the in vivo phenotypic relevance remain unclear. Utilizing ciliogenesis time course experiments combined with light microscopy and electron tomography, we reveal that chromosome 21 polyploidy elevates Pericentrin and microtubules away from the centrosome that corral MyosinVA and EHD1, delaying ciliary membrane delivery and mother centriole uncapping essential for ciliogenesis. If given enough time, trisomy 21 cells eventually ciliate, but these ciliated cells demonstrate persistent trafficking defects that reduce transition zone protein localization and decrease sonic hedgehog signaling in direct anticorrelation with Pericentrin levels. Consistent with cultured trisomy 21 cells, a mouse model of Down syndrome with elevated Pericentrin has fewer primary cilia in cerebellar granule neuron progenitors and thinner external granular layers at P4. Our work reveals that elevated Pericentrin from trisomy 21 disrupts multiple early steps of ciliogenesis and creates persistent trafficking defects in ciliated cells. This pericentrosomal crowding mechanism results in signaling deficiencies consistent with the neurological phenotypes found in individuals with Down syndrome.

Human cells typically have 23 pairs of structures known as chromosomes. Each chromosome contains a unique set of genes which provide the instructions needed to make proteins and other essential molecules found in the body. Individuals with Down syndrome have an extra copy of chromosome 21. This genetic alteration is known as trisomy 21 and affects many different organs in the body, leading to various medical conditions including intellectual disability, heart defects, and immune deficiencies. A recent study showed that cells from individuals with Down syndrome had defects in forming primary cilia ­ structures on the surface of cells which work as signaling hubs to control how cells grow and develop. These cilia defects were in large part due to excess levels of a protein known as Pericentrin, which is encoded by a gene found on chromosome 21. But it is unclear how Pericentrin disrupts cilia assembly, and how this may contribute to the medical conditions observed in individuals with Down syndrome. To address these questions, Jewett et al. studied human cells that had been engineered to have trisomy 21. The experiments found that trisomy 21 led to higher levels of Pericentrin and altered the way molecules were organized at the sites where primary cilia form. This caused the components required to build and maintain the primary cilium to become trapped in the wrong locations. The trisomy 21 cells were eventually able to rearrange the molecules and build a primary cilium, but it took them twice as long as cells with 23 pairs of chromosomes and their primary cilium did not properly work. Further experiments were then conducted on mice that had been engineered to have an extra copy of a portion of genes on human chromosome 21, including the gene for Pericentrin. Jewett et al. found that these mice assembled cilia later and had defects in cilia signaling, similar to the human trisomy 21 cells. This resulted in mild abnormalities in brain development that were consistent with what occurs in individuals with Down syndrome. These findings suggest that the elevated levels of Pericentrin in trisomy 21 causes changes in cilia formation and function which, in turn, may alter how the mouse brain develops. Further studies will be required to find out whether defects in primary cilia may contribute to other medical conditions observed in individuals with Down syndrome.

Protein coding variation in the J:ARC and J:DO outbred laboratory mouse stocks provides a molecular basis for distinct research applications.

Outbred laboratory mice (Mus musculus) are readily available and have high fecundity, making them a popular choice in biomedical research, especially toxicological and pharmacological applications. Direct high throughput genome sequencing (HTS) of these widely used research animals is an important genetic quality control measure that enhances research reproducibility. HTS data have been used to confirm the common origin of outbred stocks and to molecularly define distinct outbred populations. But these data have also revealed unexpected population structure and homozygosity in some populations; genetic features that emerge when outbred stocks are not properly maintained. We used exome sequencing to discover and interrogate protein-coding variation in a newly established population of Swiss-derived outbred stock (J:ARC) that is closely related to other, commonly used CD-1 outbred populations. We used these data to describe the genetic architecture of the J:ARC population including heterozygosity, minor allele frequency, LD decay, and we defined novel, protein-coding sequence variation. These data reveal the expected genetic architecture for a properly maintained outbred stock and provide a basis for the on-going genetic quality control. We also compared these data to protein-coding variation found in a multiparent outbred stock, the Diversity Outbred (J:DO). We found that the more recently derived, multiparent outbred stock has significantly higher interindividual variability, greater overall genetic variation, higher heterozygosity, and fewer novel variants than the Swiss-derived J:ARC stock. However, among the novel variants found in the J:DO stock, significantly more are predicted to be protein-damaging. The fact that individuals from this population can tolerate a higher load of potentially damaging variants highlights the buffering effects of allelic diversity and the differing selective pressures in these stocks. While both outbred stocks offer significant individual heterozygosity, our data provide a molecular basis for their intended applications, where the J:DO are best suited for studies requiring maximum, population-level genetic diversity and power for mapping, while the J:ARC are best suited as a general-purpose outbred stock with robust fecundity, relatively low allelic diversity, and less potential for extreme phenotypic variability.

The collaborative cross strains and their founders vary widely in cocaine-induced behavioral sensitization.

Cocaine use and overdose deaths attributed to cocaine have increased significantly in the United States in the last 10 years. Despite the prevalence of cocaine use disorder (CUD) and the personal and societal problems it presents, there are currently no approved pharmaceutical treatments. The absence of treatment options is due, in part, to our lack of knowledge about the etiology of CUDs. There is ample evidence that genetics plays a role in increasing CUD risk but thus far, very few risk genes have been identified in human studies. Genetic studies in mice have been extremely useful for identifying genetic loci and genes, but have been limited to very few genetic backgrounds, leaving substantial phenotypic, and genetic diversity unexplored. Herein we report the measurement of cocaine-induced behavioral sensitization using a 19-day protocol that captures baseline locomotor activity, initial locomotor response to an acute exposure to cocaine and locomotor sensitization across 5 exposures to the drug. These behaviors were measured in 51 genetically diverse Collaborative Cross (CC) strains along with their inbred founder strains. The CC was generated by crossing eight genetically diverse inbred strains such that each inbred CC strain has genetic contributions from each of the founder strains. Inbred CC mice are infinitely reproducible and provide a stable, yet diverse genetic platform on which to study the genetic architecture and genetic correlations among phenotypes. We have identified significant differences in cocaine locomotor sensitivity and behavioral sensitization across the panel of CC strains and their founders. We have established relationships among cocaine sensitization behaviors and identified extreme responding strains that can be used in future studies aimed at understanding the genetic, biological, and pharmacological mechanisms that drive addiction-related behaviors. Finally, we have determined that these behaviors exhibit relatively robust heritability making them amenable to future genetic mapping studies to identify addiction risk genes and genetic pathways that can be studied as potential targets for the development of novel therapeutics.

The Mutant Mouse Resource and Research Center (MMRRC): the NIH-supported National Public Repository and Distribution Archive of Mutant Mouse Models in the USA.

The Mutant Mouse Resource and Research Center (MMRRC) Program is the pre-eminent public national mutant mouse repository and distribution archive in the USA, serving as a national resource of mutant mice available to the global scientific community for biomedical research. Established more than two decades ago with grants from the National Institutes of Health (NIH), the MMRRC Program supports a Consortium of regionally distributed and dedicated vivaria, laboratories, and offices (Centers) and an Informatics Coordination and Service Center (ICSC) at three academic teaching and research universities and one non-profit genetic research institution. The MMRRC Program accepts the submission of unique, scientifically rigorous, and experimentally valuable genetically altered and other mouse models donated by academic and commercial scientists and organizations for deposition, maintenance, preservation, and dissemination to scientists upon request. The four Centers maintain an archive of nearly 60,000 mutant alleles as live mice, frozen germplasm, and/or embryonic stem (ES) cells. Since its inception, the Centers have fulfilled 13,184 orders for mutant mouse models from 9591 scientists at 6626 institutions around the globe. Centers also provide numerous services that facilitate using mutant mouse models obtained from the MMRRC, including genetic assays, microbiome analysis, analytical phenotyping and pathology, cryorecovery, mouse husbandry, infectious disease surveillance and diagnosis, and disease modeling. The ICSC coordinates activities between the Centers, manages the website (mmrrc.org) and online catalog, and conducts communication, outreach, and education to the research community. Centers preserve, secure, and protect mutant mouse lines in perpetuity, promote rigor and reproducibility in scientific experiments using mice, provide experiential training and consultation in the responsible use of mice in research, and pursue cutting edge technologies to advance biomedical studies using mice to improve human health. Researchers benefit from an expansive list of well-defined mouse models of disease that meet the highest standards of rigor and reproducibility, while donating investigators benefit by having their mouse lines preserved, protected, and distributed in compliance with NIH policies.

Genetic control of the pluripotency epigenome determines differentiation bias in mouse embryonic stem cells.

Genetically diverse pluripotent stem cells display varied, heritable responses to differentiation cues. Here, we harnessed these disparities through derivation of mouse embryonic stem cells from the BXD genetic reference panel, along with C57BL/6J (B6) and DBA/2J (D2) parental strains, to identify loci regulating cell state transitions. Upon transition to formative pluripotency, B6 stem cells quickly dissolved naïve networks adopting gene expression modules indicative of neuroectoderm lineages, whereas D2 retained aspects of naïve pluripotency. Spontaneous formation of embryoid bodies identified divergent differentiation where B6 showed a propensity toward neuroectoderm and D2 toward definitive endoderm. Genetic mapping identified major trans-acting loci co-regulating chromatin accessibility and gene expression in both naïve and formative pluripotency. These loci distally modulated occupancy of pluripotency factors at hundreds of regulatory elements. One trans-acting locus on Chr 12 primarily impacted chromatin accessibility in embryonic stem cells, while in epiblast-like cells, the same locus subsequently influenced expression of genes enriched for neurogenesis, suggesting early chromatin priming. These results demonstrate genetically determined biases in lineage commitment and identify major regulators of the pluripotency epigenome.

Behavioral phenotypes revealed during reversal learning are linked with novel genetic loci in diversity outbred mice.

Impulsive behavior and impulsivity are heritable phenotypes that are strongly associated with risk for substance use disorders. Identifying the neurogenetic mechanisms that influence impulsivity may also reveal novel biological insights into addiction vulnerability. Our past studies using the BXD and Collaborative Cross (CC) recombinant inbred mouse panels have revealed that behavioral indicators of impulsivity measured in a reversal-learning task are heritable and are genetically correlated with aspects of intravenous cocaine self-administration. Genome-wide linkage studies in the BXD panel revealed a quantitative trait locus (QTL) on chromosome 10, but we expect to identify additional QTL by testing in a population with more genetic diversity. To this end, we turned to Diversity Outbred (DO) mice; 392 DO mice (156 males, 236 females) were phenotyped using the same reversal learning test utilized previously. Our primary indicator of impulsive responding, a measure that isolates the relative difficulty mice have with reaching performance criteria under reversal conditions, revealed a genome-wide significant QTL on chromosome 7 (max LOD score = 8.73, genome-wide corrected p<0.05). A measure of premature responding akin to that implemented in the 5-choice serial reaction time task yielded a suggestive QTL on chromosome 17 (max LOD score = 9.14, genome-wide corrected <0.1). Candidate genes were prioritized (2900076A07Rik, Wdr73 and Zscan2) based upon expression QTL data we collected in DO and CC mice and analyses using publicly available gene expression and phenotype databases. These findings may advance understanding of the genetics that drive impulsive behavior and enhance risk for substance use disorders.

Heritable variation in locomotion, reward sensitivity and impulsive behaviors in a genetically diverse inbred mouse panel.

Drugs of abuse, including alcohol and stimulants like cocaine, produce effects that are subject to individual variability, and genetic variation accounts for at least a portion of those differences. Notably, research in both animal models and human subjects point toward reward sensitivity and impulsivity as being trait characteristics that predict relatively greater positive subjective responses to stimulant drugs. Here we describe use of the eight collaborative cross (CC) founder strains and 38 (reversal learning) or 10 (all other tests) CC strains to examine the heritability of reward sensitivity and impulsivity traits, as well as genetic correlations between these measures and existing addiction-related phenotypes. Strains were all tested for activity in an open field and reward sensitivity (intake of chocolate BOOST®). Mice were then divided into two counterbalanced groups and underwent reversal learning (impulsive action and waiting impulsivity) or delay discounting (impulsive choice). CC and founder mice show significant heritability for impulsive action, impulsive choice, waiting impulsivity, locomotor activity, and reward sensitivity, with each impulsive phenotype determined to be non-correlating, independent traits. This research was conducted within the broader, inter-laboratory effort of the Center for Systems Neurogenetics of Addiction (CSNA) to characterize CC and DO mice for multiple, cocaine abuse related traits. These data will facilitate the discovery of genetic correlations between predictive traits, which will then guide discovery of genes and genetic variants that contribute to addictive behaviors.