A platform for experimental precision medicine: The extended BXD mouse family

Current Advancements in Mouse Models of Retinal Disease

The field of retinal degenerative (RDs) disease study has been in a state of exponential growth from discovering the underlying genetic components of such diseases as age-related macular degeneration (AMD) and retinitis pigmentosa (RP) to the first gene therapy developed and approved for human Leber congenital amaurosis. However, a source for high-fidelity animal models of these complex, multifactorial, and/or polygenic diseases is a need that has yet to be fulfilled. While models for AMD and RP do exist, they often require aging the animals for a year or more, feeding special diets, or introduction of external modulators such as exposure to cigarette smoke. Currently, work is being done to uncover high-fidelity naturally occurring models of these retinal diseases with the hope and intent of providing the vision community the tools it needs to better understand, treat, and, one day, cure the patients suffering from these devastating afflictions.

Weighted gene co-expression network analysis identifies key hub genes and pathways in acute myeloid leukemia

Introduction: Acute myeloid leukemia (AML) is the most common type of leukemia in adults. However, there is a gap in understanding the molecular basis of the disease, partly because key genes associated with AML have not been extensively explored. In the current study, we aimed to identify genes that have strong association with AML based on a cross-species integrative approach. Methods: We used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify co-expressed gene modules significantly correlated with human AML, and further selected the genes exhibiting a significant difference in expression between AML and healthy mouse. Protein-protein interactions, transcription factors, gene function, genetic regulation, and coding sequence variants were integrated to identify key hub genes in AML. Results: The cross-species approach identified a total of 412 genes associated with both human and mouse AML. Enrichment analysis confirmed an association of these genes with hematopoietic and immune-related functions, phenotypes, processes, and pathways. Further, the integrated analysis approach identified a set of important module genes including Nfe2, Trim27, Mef2c, Ets1, Tal1, Foxo1, and Gata1 in AML. Six of these genes (except ETS1) showed significant differential expression between human AML and healthy samples in an independent microarray dataset. All of these genes are known to be involved in immune/hematopoietic functions, and in transcriptional regulation. In addition, Nfe2, Trim27, Mef2c, and Ets1 harbor coding sequence variants, whereas Nfe2 and Trim27 are cis-regulated, making them attractive candidates for validation. Furthermore, subtype-specific analysis of the hub genes in human AML indicated high expression of NFE2 across all the subtypes (M0 through M7) and enriched expression of ETS1, LEF1, GATA1, and TAL1 in M6 and M7 subtypes. A significant correlation between methylation status and expression level was observed for most of these genes in AML patients. Conclusion: Findings from the current study highlight the importance of our cross-species approach in the identification of multiple key candidate genes in AML, which can be further studied to explore their detailed role in leukemia/AML.

Transcriptome-wide association study reveals cholesterol metabolism gene Lpl is a key regulator of cognitive dysfunction

Cholesterol metabolism in the brain plays a crucial role in normal physiological function, and its aberrations are associated with cognitive dysfunction. The present study aimed to determine which cholesterol-related genes play a vital role in cognitive dysfunction and to dissect its underlying molecular mechanisms using a systems genetics approach in the BXD mice family. We first systematically analyzed the association of expression of 280 hippocampal genes related to cholesterol metabolism with cognition-related traits and identified lipoprotein lipase (Lpl) as a critical regulator. This was further confirmed by phenome-wide association studies that indicate Lpl associated with hippocampus volume residuals and anxiety-related traits. By performing expression quantitative trait locus mapping, we demonstrate that Lpl is strongly cis-regulated in the BXD hippocampus. We also identified ∼3,300 genes significantly (p < 0.05) correlated with the Lpl expression. Those genes are mainly involved in the regulation of neuron-related traits through the MAPK signaling pathway, axon guidance, synaptic vesicle cycle, and NF-kappa B signaling pathway. Furthermore, a protein-protein interaction network analysis identified several direct interactors of Lpl, including Rab3a, Akt1, Igf1, Crp, and Lrp1, which indicates that Lpl involves in the regulation of cognitive dysfunction through Rab3a-mediated synaptic vesicle cycle and Akt1/Igf1/Crp/Lrp1-mediated MAPK signaling pathway. Our findings demonstrate the importance of the Lpl, among the cholesterol-related genes, in regulating cognitive dysfunction and highlighting the potential signaling pathways, which may serve as novel therapeutic targets for the treatment of cognitive dysfunction.

Proteome-wide systems genetics identifies UFMylation as a regulator of skeletal muscle function

Improving muscle function has great potential to improve the quality of life. To identify novel regulators of skeletal muscle metabolism and function, we performed a proteomic analysis of gastrocnemius muscle from 73 genetically distinct inbred mouse strains, and integrated the data with previously acquired genomics and >300 molecular/phenotypic traits via quantitative trait loci mapping and correlation network analysis. These data identified thousands of associations between protein abundance and phenotypes and can be accessed online (https://muscle.coffeeprot.com/) to identify regulators of muscle function. We used this resource to prioritize targets for a functional genomic screen in human bioengineered skeletal muscle. This identified several negative regulators of muscle function including UFC1, an E2 ligase for protein UFMylation. We show UFMylation is up-regulated in a mouse model of amyotrophic lateral sclerosis, a disease that involves muscle atrophy. Furthermore, in vivo knockdown of UFMylation increased contraction force, implicating its role as a negative regulator of skeletal muscle function.

Multi-omics analysis reveals expression complexity and functional diversity of mouse kinome

Protein kinases are a crucial component of signaling pathways involved in a wide range of cellular responses, including growth, proliferation, differentiation, and migration. Systematic investigation of protein kinases is critical to better understand phosphorylation-mediated signaling pathways and may provide insights into the development of potential therapeutic drug targets. Here we perform a systems-level analysis of the mouse kinome by analyzing multi-omics data. We used bulk and single-cell transcriptomic data from the C57BL/6J mouse strain to define tissue- and cell-type-specific expression of protein kinases, followed by investigating variations in sequence and expression between C57BL/6J and DBA/2J strains. We then profiled a deep brain phosphoproteome from C57BL/6J and DBA/2J strains as well as their reciprocal hybrids to infer the activity of the mouse kinome. Finally, we performed phenome-wide association analysis using the BXD recombinant inbred (RI) mice (a cross between C57BL/6J and DBA/2J strains) to identify any associations between variants in protein kinases and phenotypes. Collectively, our study provides a comprehensive analysis of the mouse kinome by investigating genetic sequence variation, tissue-specific expression patterns, and associations with downstream phenotypes.

Genomic Basis for Individual Differences in Susceptibility to the Neurotoxic Effects of Diesel Exhaust

Air pollution is a known environmental health hazard. A major source of air pollution includes diesel exhaust (DE). Initially, research on DE focused on respiratory morbidities; however, more recently, exposures to DE have been associated with neurological developmental disorders and neurodegeneration. In this study, we investigated the effects of sub-chronic inhalation exposure to DE on neuroinflammatory markers in two inbred mouse strains and both sexes, including whole transcriptome examination of the medial prefrontal cortex. We exposed aged male and female C57BL/6J (B6) and DBA/2J (D2) mice to DE, which was cooled and diluted with HEPA-filtered compressed air for 2 h per day, 5 days a week, for 4 weeks. Control animals were exposed to HEPA-filtered air on the same schedule as DE-exposed animals. The prefrontal cortex was harvested and analyzed for proinflammatory cytokine gene expression (Il1β, Il6, Tnfα) and transcriptome-wide response by RNA-seq. We observed differential cytokine gene expression between strains and sexes in the DE-exposed vs. control-exposed groups for Il1β, Tnfα, and Il6. For RNA-seq, we identified 150 differentially expressed genes between air and DE treatment related to natural killer cell-mediated cytotoxicity per Kyoto Encyclopedia of Genes and Genomes pathways. Overall, our data show differential strain-related effects of DE on neuroinflammation and neurotoxicity and demonstrate that B6 are more susceptible than D2 to gene expression changes due to DE exposures than D2. These results are important because B6 mice are often used as the default mouse model for DE studies and strain-related effects of DE neurotoxicity warrant expanded studies.

Brain-iron deficiency models of restless legs syndrome

Restless legs syndrome (RLS) is a common sensorimotor disorder for which two main pathological elements are fairly well accepted: Brain iron deficiency (BID) and an altered dopaminergic system. The ability to better understand the causal and consequential factors related to these two pathological elements, would hopefully lead to the development of better therapeutic strategies for treating, if not curing, this disease. The current understanding of the relationship between these two elements is that BID leads to some alterations in neurotransmitters and subsequent changes in the dopaminergic system. Therefore, rodent models based on diet-induced BID, provide a biological substrate to understand the consequences of BID on dopaminergic pathway and on alternative pathways that may be involved. In this review, we present the current research on dopaminergic changes found in RLS subjects and compare that to what is seen in the BID rodent model to provide a validation of the BID rodent model. We also demonstrate the ability of the BID model to predict changes in other neurotransmitter systems and how that has led to new treatment options. Finally, we will present arguments for the utility of recombinant inbred mouse strains that demonstrate natural variation in brain iron, to explore the genetic basis of altered brain iron homeostasis as a model to understand why in idiopathic RLS there can exist a BID despite normal peripheral iron store. This review is the first to draw on 25 years of human and basic research into the pathophysiology of RLS to provide strong supportive data as to the validity of BID model as an important translational model of the disease. As we will demonstrate here, not only does the BID model closely and accurately mimic what we see in the dopaminergic system of RLS, it is the first model to identify alternative systems from which new treatments have recently been developed.

Meiotic drive in house mice: mechanisms, consequences, and insights for human biology

Meiotic drive occurs when one allele at a heterozygous site cheats its way into a disproportionate share of functional gametes, violating Mendel's law of equal segregation. This genetic conflict typically imposes a fitness cost to individuals, often by disrupting the process of gametogenesis. The evolutionary impact of meiotic drive is substantial, and the phenomenon has been associated with infertility and reproductive isolation in a wide range of organisms. However, cases of meiotic drive in humans remain elusive, a finding that likely reflects the inherent challenges of detecting drive in our species rather than unique features of human genome biology. Here, we make the case that house mice (Mus musculus) present a powerful model system to investigate the mechanisms and consequences of meiotic drive and facilitate translational inferences about the scope and potential mechanisms of drive in humans. We first detail how different house mouse resources have been harnessed to identify cases of meiotic drive and the underlying mechanisms utilized to override Mendel's rules of inheritance. We then summarize the current state of knowledge of meiotic drive in the mouse genome. We profile known mechanisms leading to transmission bias at several established drive elements. We discuss how a detailed understanding of meiotic drive in mice can steer the search for drive elements in our own species. Lastly, we conclude with a prospective look into how new technologies and molecular tools can help resolve lingering mysteries about the prevalence and mechanisms of selfish DNA transmission in mammals.

A Multitrait Locus Regulates Sarbecovirus Pathogenesis

Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to severe acute respiratory syndrome coronavirus (SARS-CoV) disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6, that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2, and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species. IMPORTANCE Host genetic variation is an important determinant that predicts disease outcomes following infection. In the setting of highly pathogenic coronavirus infections genetic determinants underlying host susceptibility and mortality remain unclear. To elucidate the role of host genetic variation on sarbecovirus pathogenesis and disease outcomes, we utilized the Collaborative Cross (CC) mouse genetic reference population as a model to identify susceptibility alleles to SARS-CoV and SARS-CoV-2 infections. Our findings reveal that a multitrait loci found in chromosome 9 is an important regulator of sarbecovirus pathogenesis in mice. Within this locus, we identified and validated CCR9 and CXCR6 as important regulators of host disease outcomes. Specifically, both CCR9 and CXCR6 are protective against severe SARS-CoV, SARS-CoV-2, and SARS-related HKU3 virus disease in mice. This chromosome 9 multitrait locus may be important to help identify genes that regulate coronavirus disease outcomes in humans.

Unraveling the mystery: How bad is BAG3 in hematological malignancies?

BAG3, also known as BIS and CAIR-1, interacts with Hsp70 via its BAG domain and with other molecules through its WW domain, PXXP repeats and IPV motifs. BAG3 can participate in major cellular pathways including apoptosis, autophagy, cytoskeleton structure, and motility by regulating the expression, location, and activity of its chaperone proteins. As a multifunctional protein, BAG3 is highly expressed in skeletal muscle, cardiomyocytes and multiple tumors, and its intracellular expression can be stimulated by stress. The functions and mechanisms of BAG3 in hematological malignancies have recently been a topic of interest. BAG3 has been confirmed to be involved in the development and chemoresistance of hematological malignancies and to act as a prognostic indicator. Modulation of BAG3 and its corresponding proteins has thus emerged as a promising therapeutic and experimental target. In this review, we consider the characteristics of BAG3 in hematological malignancies as a reference for further clinical and fundamental investigations.

Accounting for diet and age

The diet and age of mice can modulate how different genetic variants impact body weight, demonstrating the need to take context into account when performing genetic studies.

Evaluating the power and limitations of genome-wide association studies in Caenorhabditis elegans

Quantitative genetics in Caenorhabditis elegans seeks to identify naturally segregating genetic variants that underlie complex traits. Genome-wide association studies scan the genome for individual genetic variants that are significantly correlated with phenotypic variation in a population, or quantitative trait loci. Genome-wide association studies are a popular choice for quantitative genetic analyses because the quantitative trait loci that are discovered segregate in natural populations. Despite numerous successful mapping experiments, the empirical performance of genome-wide association study has not, to date, been formally evaluated in C. elegans. We developed an open-source genome-wide association study pipeline called NemaScan and used a simulation-based approach to provide benchmarks of mapping performance in collections of wild C. elegans strains. Simulated trait heritability and complexity determined the spectrum of quantitative trait loci detected by genome-wide association studies. Power to detect smaller-effect quantitative trait loci increased with the number of strains sampled from the C. elegans Natural Diversity Resource. Population structure was a major driver of variation in mapping performance, with populations shaped by recent selection exhibiting significantly lower false discovery rates than populations composed of more divergent strains. We also recapitulated previous genome-wide association studies of experimentally validated quantitative trait variants. Our simulation-based evaluation of performance provides the community with critical context to pursue quantitative genetic studies using the C. elegans Natural Diversity Resource to elucidate the genetic basis of complex traits in C. elegans natural populations.

A Multitrait Locus Regulates Sarbecovirus Pathogenesis

Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to SARS-CoV disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse Chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6 that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2 and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species.

Genetic Variants Associated With Resilience in Human and Animal Studies

Resilience is broadly defined as the ability to maintain or regain functioning in the face of adversity and is influenced by both environmental and genetic factors. The identification of specific genetic factors and their biological pathways underpinning resilient functioning can help in the identification of common key factors, but heterogeneities in the operationalisation of resilience have hampered advances. We conducted a systematic review of genetic variants associated with resilience to enable the identification of general resilience mechanisms. We adopted broad inclusion criteria for the definition of resilience to capture both human and animal model studies, which use a wide range of resilience definitions and measure very different outcomes. Analyzing 158 studies, we found 71 candidate genes associated with resilience. OPRM1 (Opioid receptor mu 1), NPY (neuropeptide Y), CACNA1C (calcium voltage-gated channel subunit alpha1 C), DCC (deleted in colorectal carcinoma), and FKBP5 (FKBP prolyl isomerase 5) had both animal and human variants associated with resilience, supporting the idea of shared biological pathways. Further, for OPRM1, OXTR (oxytocin receptor), CRHR1 (corticotropin-releasing hormone receptor 1), COMT (catechol-O-methyltransferase), BDNF (brain-derived neurotrophic factor), APOE (apolipoprotein E), and SLC6A4 (solute carrier family 6 member 4), the same allele was associated with resilience across divergent resilience definitions, which suggests these genes may therefore provide a starting point for further research examining commonality in resilience pathways.

Translational approaches to understanding resilience to Alzheimer's disease

Individuals who maintain cognitive function despite high levels of Alzheimer's disease (AD)-associated pathology are said to be 'resilient' to AD. Identifying mechanisms underlying resilience represents an exciting therapeutic opportunity. Human studies have identified a number of molecular and genetic factors associated with resilience, but the complexity of these cohorts prohibits a complete understanding of which factors are causal or simply correlated with resilience. Genetically and phenotypically diverse mouse models of AD provide new and translationally relevant opportunities to identify and prioritize new resilience mechanisms for further cross-species investigation. This review will discuss insights into resilience gained from both human and animal studies and highlight future approaches that may help translate these insights into therapeutics designed to prevent or delay AD-related dementia.

A natural mutator allele shapes mutation spectrum variation in mice

Although germline mutation rates and spectra can vary within and between species, common genetic modifiers of the mutation rate have not been identified in nature1. Here we searched for loci that influence germline mutagenesis using a uniquely powerful resource: a panel of recombinant inbred mouse lines known as the BXD, descended from the laboratory strains C57BL/6J (B haplotype) and DBA/2J (D haplotype). Each BXD lineage has been maintained by brother-sister mating in the near absence of natural selection, accumulating de novo mutations for up to 50 years on a known genetic background that is a unique linear mosaic of B and D haplotypes2. We show that mice inheriting D haplotypes at a quantitative trait locus on chromosome 4 accumulate C>A germline mutations at a 50% higher rate than those inheriting B haplotypes, primarily owing to the activity of a C>A-dominated mutational signature known as SBS18. The B and D quantitative trait locus haplotypes encode different alleles of Mutyh, a DNA repair gene that underlies the heritable cancer predisposition syndrome that causes colorectal tumors with a high SBS18 mutation load3,4. Both B and D Mutyh alleles are present in wild populations of Mus musculus domesticus, providing evidence that common genetic variation modulates germline mutagenesis in a model mammalian species.

Chronic Proinflammatory Signaling Accelerates the Rate of Degeneration in a Spontaneous Polygenic Model of Inherited Retinal Dystrophy

Collectively, retinal neurodegenerative diseases are comprised of numerous subtypes of disorders which result in loss of a varying cell types in the retina. These diseases can range from glaucoma, which results in retinal ganglion cell death, to age-related macular degeneration and retinitis pigmentosa, which result in cell death of the retinal pigment epithelium, photoreceptors, or both. Regardless of the disease, it’s been recently found that increased release of proinflammatory cytokines and proliferation of active microglia result in a remarkably proinflammatory microenvironment that assists in the pathogenesis of the disease; however, many of the details of these inflammatory events have yet to be elucidated. In an ongoing study, we have used systems genetics to identify possible models of spontaneous polygenic age-related macular degeneration by mining the BXD family of mice using single nucleotide polymorphism analyses of known genes associated with the human retinal disease. One BXD strain (BXD32) was removed from the study as the rate of degeneration observed in these animals was markedly increased with a resultant loss of most all photoreceptors by 6 months of age. Using functional and anatomical exams including optokinetic nystamography, funduscopy, fluorescein angiography, and optical coherence tomography, along with immunohistochemical analyses, we show that the BXD32 mouse strain exhibits a severe neurodegenerative phenotype accompanied by adverse effects on the retinal vasculature. We also expose the concurrent establishment of a chronic proinflammatory microenvironment including the TNFα secretion and activation of the NF-κB and JAK/STAT pathways with an associated increase in activated macrophages and phagoptosis. We conclude that the induced neuronal death and proinflammatory pathways work synergistically in the disease pathogenesis to enhance the rate of degeneration in this spontaneous polygenic model of inherited retinal dystrophy.

Genetic loci and metabolic states associated with murine epigenetic aging

Changes in DNA methylation (DNAm) are linked to aging. Here, we profile highly conserved CpGs in 339 predominantly female mice belonging to the BXD family for which we have deep longevity and genomic data. We use a 'pan-mammalian' microarray that provides a common platform for assaying the methylome across mammalian clades. We computed epigenetic clocks and tested associations with DNAm entropy, diet, weight, metabolic traits, and genetic variation. We describe the multifactorial variance of methylation at these CpGs and show that high-fat diet augments the age-related changes. Entropy increases with age. The progression to disorder, particularly at CpGs that gain methylation over time, was predictive of genotype-dependent life expectancy. The longer-lived BXD strains had comparatively lower entropy at a given age. We identified two genetic loci that modulate epigenetic age acceleration (EAA): one on chromosome (Chr) 11 that encompasses the Erbb2/Her2 oncogenic region, and the other on Chr19 that contains a cytochrome P450 cluster. Both loci harbor genes associated with EAA in humans, including STXBP4, NKX2-3, and CUTC. Transcriptome and proteome analyses revealed correlations with oxidation-reduction, metabolic, and immune response pathways. Our results highlight concordant loci for EAA in humans and mice, and demonstrate a tight coupling between the metabolic state and epigenetic aging.

New Insights on Gene by Environmental Effects of Drugs of Abuse in Animal Models Using GeneNetwork

Gene-by-environment interactions are important for all facets of biology, especially behaviour. Families of isogenic strains of mice, such as the BXD strains, are excellently placed to study these interactions, as the same genome can be tested in multiple environments. BXD strains are recombinant inbred mouse strains derived from crossing two inbred strains—C57BL/6J and DBA/2J mice. Many reproducible genometypes can be leveraged, and old data can be reanalysed with new tools to produce novel insights. We obtained drug and behavioural phenotypes from Philip et al. Genes, Brain and Behaviour 2010, and reanalysed their data with new genotypes from sequencing, as well as new models (Genome-wide Efficient Mixed Model Association (GEMMA) and R/qtl2). We discovered QTLs on chromosomes 3, 5, 9, 11, and 14, not found in the original study. We reduced the candidate genes based on their ability to alter gene expression or protein function. Candidate genes included Slitrk6 and Cdk14. Slitrk6, in a Chromosome14 QTL for locomotion, was found to be part of a co-expression network involved in voluntary movement and associated with neuropsychiatric phenotypes. Cdk14, one of only three genes in a Chromosome5 QTL, is associated with handling induced convulsions after ethanol treatment, that is regulated by the anticonvulsant drug valproic acid. By using families of isogenic strains, we can reanalyse data to discover novel candidate genes involved in response to drugs of abuse.

A systems genetics approach delineates the role of Bcl2 in leukemia pathogenesis

Leukemia is a group of malignancies of the blood forming tissues, and is characterized by the uncontrolled proliferation of blood cells. In the United States, it accounts for approximately 3.5% and 4% of all cancer-related incidences and mortalities, respectively. The current study aimed to explore the role of Bcl2 and associated genes in leukemia pathogenesis using a systems genetics approach. The transcriptome data from BXD Recombinant Inbred (RI) mice was analyzed to identify the expression of Bcl2 in myeloid cells. eQTL mapping was performed to select the potential chromosomal region and subsequently identify the candidate gene modulating the expression of Bcl2. Furthermore, gene enrichment and protein-protein interaction (PPI) analyses of the Bcl2-coexpressed genes were performed to demonstrate the role of Bcl2 in leukemia pathogenesis. The Bcl2-coexpressed genes were found to be enriched in various hematopoietic system related functions, and multiple pathways related to signaling, immune response, and cancer. The PPI network analysis demonstrated direct interaction of hematopoietic function related genes, such as Bag3, Bak1, Bcl2l11, Bmf, Mapk9, Myc, Ppp2r5c, and Ppp3ca with Bcl2. The eQTL mapping identified a 4.5 Mb genomic region on chromosome 11, potentially regulating the expression of Bcl2. A multi-criteria filtering process identified Top2a, among the genes located in the mapped locus, as the best candidate upstream regulator for Bcl2 expression variation. Hence, the current study provides better insights into the role of Bcl2 in leukemia pathogenesis and demonstrates the significance of our approach in gaining new knowledge on leukemia. Furthermore, our findings from the PPI network analysis and eQTL mapping provide supporting evidence of leukemia-associated genes, which can be further explored for their functional importance in leukemia. DATA AVAILABILITY: The myeloid cell transcriptomic data of the BXD mice used in this study can be accessed through our GeneNetwork (http://www.genenetwork.org) with the accession number of GN144.

Genetically diverse mouse models of SARS-CoV-2 infection reproduce clinical variation and cytokine responses in COVID-19

Inflammation in response to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection drives severity of coronavirus disease 2019 (COVID-19), with effective versus dysregulated responses influenced by host genetics. To understand mechanisms of inflammation, animal models that reflect genetic diversity and clinical outcomes observed in humans are needed. We report a mouse panel comprising the diverse genetic backgrounds of the Collaborative Cross founder strains crossed to K18-hACE2 transgenic mice that confers high susceptibility to SARS-CoV-2. Infection of CC x K18-hACE2 F1 progeny resulted in a spectrum of weight loss, survival, viral replication kinetics, histopathology, and cytokine profiles, some of which were sex-specific. Importantly, survival was associated with early type I interferon (IFN) expression and a phased proinflammatory response distinct from mice with severe disease. Thus, dynamics of inflammatory responses observed in COVID-19 can be modeled in diverse mouse strains that provide a genetically tractable platform for understanding antiviral immunity and evaluating countermeasures.

ONE SENTENCE SUMMARY

Genetically diverse mice display a broad spectrum of clinically relevant responses to SARS-CoV-2 infection, reflecting variability in COVID-19 disease.

The Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel

BACKGROUND

Unraveling the relationship between genetic variation and phenotypic traits remains a fundamental challenge in biology. Mapping variants underlying complex traits while controlling for confounding environmental factors is often problematic. To address this, we establish a vertebrate genetic resource specifically to allow for robust genotype-to-phenotype investigations. The teleost medaka (Oryzias latipes) is an established genetic model system with a long history of genetic research and a high tolerance to inbreeding from the wild.

RESULTS

Here we present the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel: the first near-isogenic panel of 80 inbred lines in a vertebrate model derived from a wild founder population. Inbred lines provide fixed genomes that are a prerequisite for the replication of studies, studies which vary both the genetics and environment in a controlled manner, and functional testing. The MIKK panel will therefore enable phenotype-to-genotype association studies of complex genetic traits while allowing for careful control of interacting factors, with numerous applications in genetic research, human health, drug development, and fundamental biology.

CONCLUSIONS

Here we present a detailed characterization of the genetic variation across the MIKK panel, which provides a rich and unique genetic resource to the community by enabling large-scale experiments for mapping complex traits.

Multiomic profiling of the liver across diets and age in a diverse mouse population

We profiled the liver transcriptome, proteome, and metabolome in 347 individuals from 58 isogenic strains of the BXD mouse population across age (7 to 24 months) and diet (low or high fat) to link molecular variations to metabolic traits. Several hundred genes are affected by diet and/or age at the transcript and protein levels. Orthologs of two aging-associated genes, St7 and Ctsd, were knocked down in C. elegans, reducing longevity in wild-type and mutant long-lived strains. The multiomics data were analyzed as segregating gene networks according to each independent variable, providing causal insight into dietary and aging effects. Candidates were cross-examined in an independent diversity outbred mouse liver dataset segregating for similar diets, with ∼80%-90% of diet-related candidate genes found in common across datasets. Together, we have developed a large multiomics resource for multivariate analysis of complex traits and demonstrate a methodology for moving from observational associations to causal connections.

Working memory and pattern separation in founder strains of the BXD recombinant inbred mouse panel

Working memory and pattern separation are fundamental cognitive abilities which, when impaired, significantly diminish quality of life. Discovering genetic mechanisms underlying innate and disease-induced variation in these cognitive abilities is a critical step towards treatments for common and devastating neurodegenerative conditions such as Alzheimer's disease. In this regard, the trial-unique nonmatching-to-location assay (TUNL) is a touchscreen operant conditioning procedure allowing simultaneous quantification of working memory and pattern separation in mice and rats. In the present study, we used the TUNL assay to quantify these cognitive abilities in C57BL/6J and DBA/2J mice. These strains are the founders of the BXD recombinant inbred mouse panel which enables discovery of genetic mechanisms underlying phenotypic variation. TUNL testing revealed that pattern separation was significantly influenced by mouse strain, whereas working memory was not. Moreover, horizontal distance and vertical distance between choice-phase stimuli had dissociable effects on TUNL performance. These findings provide novel data on mouse strain differences in pattern separation and support previous findings of equivalent working memory performance in C57BL/6J and DBA/2J mice. Although working memory of the BXD founder strains was equivalent in this study, working memory of BXD strains may be divergent because of transgressive segregation. Collectively, data presented here indicate that pattern separation is heritable in the mouse and that the BXD panel can be used to identify mechanisms underlying variation in pattern separation.

Identification of cyclin D1 as a major modulator of 3-nitropropionic acid-induced striatal neurodegeneration

Mitochondria dysfunction occurs in the aging brain as well as in several neurodegenerative disorders and predisposes neuronal cells to enhanced sensitivity to neurotoxins. 3-nitropropionic acid (3-NP) is a naturally occurring plant and fungal neurotoxin that causes neurodegeneration predominantly in the striatum by irreversibly inhibiting the tricarboxylic acid respiratory chain enzyme, succinate dehydrogenase (SDH), the main constituent of the mitochondria respiratory chain complex II. Significantly, although 3-NP-induced inhibition of SDH occurs in all brain regions, neurodegeneration occurs primarily and almost exclusively in the striatum for reasons still not understood. In rodents, 3-NP-induced striatal neurodegeneration depends on the strain background suggesting that genetic differences among genotypes modulate toxicant variability and mechanisms that underlie 3-NP-induced neuronal cell death. Using the large BXD family of recombinant inbred (RI) strains we demonstrate that variants in Ccnd1 - the gene encoding cyclin D1 - of the DBA/2 J parent underlie the resistance to 3-NP-induced striatal neurodegeneration. In contrast, the Ccnd1 variant inherited from the widely used C57BL/6 J parental strain confers sensitivity. Given that cellular stress triggers induction of cyclin D1 expression followed by cell-cycle re-entry and consequent neuronal cell death, we sought to determine if the C57BL/6 J and DBA/2 J Ccnd1 variants are differentially modulated in response to 3-NP. We confirm that 3-NP induces cyclin D1 expression in striatal neuronal cells of C57BL/6 J, but this response is blunted in the DBA/2 J. We further show that striatal-specific alternative processing of a highly conserved 3'UTR negative regulatory region of Ccnd1 co-segregates with the C57BL/6 J parental Ccnd1 allele in BXD strains and that its differential processing accounts for sensitivity or resistance to 3-NP. Our results indicate that naturally occurring Ccnd1 variants may play a role in the variability observed in neurodegenerative disorders involving mitochondria complex II dysfunction and point to cyclin D1 as a possible therapeutic target.

Hippocampal Transcriptome-Wide Association Study Reveals Correlations Between Impaired Glutamatergic Synapse Pathway and Age-Related Hearing Loss in BXD-Recombinant Inbred Mice

Age-related hearing loss (ARHL) is associated with cognitive dysfunction; however, the detailed underlying mechanisms remain unclear. The aim of this study is to investigate the potential underlying mechanism with a system genetics approach. A transcriptome-wide association study was performed on aged (12–32 months old) BXD mice strains. The hippocampus gene expression was obtained from 56 BXD strains, and the hearing acuity was assessed from 54 BXD strains. Further correlation analysis identified a total of 1,435 hearing-related genes in the hippocampus (p < 0.05). Pathway analysis of these genes indicated that the impaired glutamatergic synapse pathway is involved in ARHL (p = 0.0038). Further gene co-expression analysis showed that the expression level of glutamine synthetase (Gls), which is significantly correlated with ARHL (n = 26, r = −0.46, p = 0.0193), is a crucial regulator in glutamatergic synapse pathway and associated with learning and memory behavior. In this study, we present the first systematic evaluation of hippocampus gene expression pattern associated with ARHL, learning, and memory behavior. Our results provide novel potential molecular mechanisms involved in ARHL and cognitive dysfunction association.

Genetic architecture of protein expression and its regulation in the mouse brain

BACKGROUND

Natural variation in protein expression is common in all organisms and contributes to phenotypic differences among individuals. While variation in gene expression at the transcript level has been extensively investigated, the genetic mechanisms underlying variation in protein expression have lagged considerably behind. Here we investigate genetic architecture of protein expression by profiling a deep mouse brain proteome of two inbred strains, C57BL/6 J (B6) and DBA/2 J (D2), and their reciprocal F1 hybrids using two-dimensional liquid chromatography coupled with tandem mass spectrometry (LC/LC-MS/MS) technology.

RESULTS

By comparing protein expression levels in the four mouse strains, we observed 329 statistically significant differentially expressed proteins between the two parental strains and characterized the genetic basis of protein expression. We further applied a proteogenomic approach to detect variant peptides and define protein allele-specific expression (pASE), identifying 33 variant peptides with cis-effects and 17 variant peptides showing trans-effects. Comparison of regulation at transcript and protein levels show a significant divergence.

CONCLUSIONS

The results provide a comprehensive analysis of genetic architecture of protein expression and the contribution of cis- and trans-acting regulatory differences to protein expression.

BXD Recombinant Inbred Mice as a Model to Study Neurotoxicity

BXD recombinant inbred (RI) lines represent a genetic reference population derived from a cross between C57BL/6J mice (B6) and DBA/2J mice (D2), which through meiotic recombination events possesses recombinant chromosomes containing B6 or D2 haplotype segments. The quantitative trait loci (QTLs) are the locations of segregating genetic polymorphisms and are fundamental to understanding genetic diversity in human disease susceptibility and severity. QTL mapping represents the typical approach for identifying naturally occurring polymorphisms that influence complex phenotypes. In this process, genotypic values at markers of known genomic locations are associated with phenotypic values measured in a segregating population. Indeed, BXD RI strains provide a powerful tool to study neurotoxicity induced by different substances. In this review, we describe the use of BXD RI lines to understand the underlying mechanisms of neurotoxicity in response to ethanol and cocaine, as well as metals and pesticide exposures.

Ace2 and Tmprss2 Expressions Are Regulated by Dhx32 and Influence the Gastrointestinal Symptoms Caused by SARS-CoV-2

Studies showed that the gastrointestinal (GI) tract is one of the most important pathways for SARS-CoV-2 infection and coronavirus disease 2019 (COVID-19). As SARS-CoV-2 cellular entry depends on the ACE2 receptor and TMPRSS2 priming of the spike protein, it is important to understand the molecular mechanisms through which these two proteins and their cognate transcripts interact and influence the pathogenesis of COVID-19. In this study, we quantified the expression, associations, genetic modulators, and molecular pathways for Tmprss2 and Ace2 mRNA expressions in GI tissues using a systems genetics approach and the expanded family of highly diverse BXD mouse strains. The results showed that both Tmprss2 and Ace2 are highly expressed in GI tissues with significant covariation. We identified a significant expression quantitative trait locus on chromosome 7 that controls the expression of both Tmprss2 and Ace2. Dhx32 was found to be the strongest candidate in this interval. Co-expression network analysis demonstrated that both Tmprss2 and Ace2 were located at the same module that is significantly associated with other GI-related traits. Protein-protein interaction analysis indicated that hub genes in this module are linked to circadian rhythms. Collectively, our data suggested that genes with circadian rhythms of expression may have an impact on COVID-19 disease, with implications related to the timing and treatment of COVID-19.

Pursuit of precision medicine: Systems biology approaches in Alzheimer's disease mouse models

Alzheimer's disease (AD) is a complex disease that is mediated by numerous factors and manifests in various forms. A systems biology approach to studying AD involves analyses of various body systems, biological scales, environmental elements, and clinical outcomes to understand the genotype to phenotype relationship that potentially drives AD development. Currently, there are many research investigations probing how modifiable and nonmodifiable factors impact AD symptom presentation. This review specifically focuses on how imaging modalities can be integrated into systems biology approaches using model mouse populations to link brain level functional and structural changes to disease onset and progression. Combining imaging and omics data promotes the classification of AD into subtypes and paves the way for precision medicine solutions to prevent and treat AD.

Paraquat Toxicogenetics: Strain-Related Reduction of Tyrosine Hydroxylase Staining in Substantia Nigra in Mice

Paraquat (PQ) is a putative risk factor for the development of sporadic Parkinson's disease. To model a possible genetic basis for individual differences in susceptibility to exposure to PQ, we recently examined the effects of paraquat on tyrosine hydroxylase (TH)-containing neurons in the substantia nigra pars compacta (SNc) of six members of the BXD family of mice (n = 2-6 per strain). We injected males with 5 mg/kg paraquat weekly three times. The density of TH+ neurons counted by immunocytochemistry at 200x in eight or more sections through the SNc is reduced in five of the six strains relative to control (N = 4 ± 2 mice per strain). TH+ loss ranged from 0 to 20% with an SEM of 1%. The heritability was estimated using standard ANOVA and jackknife resampling and is 0.37 ± 0.05 in untreated animals and 0.47 ± 0.04 in treated animals. These results demonstrate genetic modulation and GxE variation in susceptibility to PQ exposure and the loss of TH staining in the substantia nigra.

Regulation of protein abundance in genetically diverse mouse populations

Genetically diverse mouse populations are powerful tools for characterizing the regulation of the proteome and its relationship to whole-organism phenotypes. We used mass spectrometry to profile and quantify the abundance of 6,798 proteins in liver tissue from mice of both sexes across 58 Collaborative Cross (CC) inbred strains. We previously collected liver proteomics data from the related Diversity Outbred (DO) mice and their founder strains. We show concordance across the proteomics datasets despite being generated from separate experiments, allowing comparative analysis. We map protein abundance quantitative trait loci (pQTLs), identifying 1,087 local and 285 distal in the CC mice and 1,706 local and 414 distal in the DO mice. We find that regulatory effects on individual proteins are conserved across the mouse populations, in particular for local genetic variation and sex differences. In comparison, proteins that form complexes are often co-regulated, displaying varying genetic architectures, and overall show lower heritability and map fewer pQTLs. We have made this resource publicly available to enable quantitative analyses of the regulation of the proteome.

Genetic differences in ethanol consumption: effects on iron, copper, and zinc regulation in mouse hippocampus

One common characteristic of neurodegenerative diseases is dysregulation of iron, usually with observed increases in its concentration in various regions. Heavy alcohol consumption is believed to contribute to such iron dysregulation in the brain with accompanying dementia. To examine this effect and related genetic-based individual differences in an animal model, we subjected female mice from 12 BXD recombinant inbred strains to 16 weeks of alcohol consumption using the drinking in the dark (DID) method. Daily consumption was recorded and at the end of 16 weeks hippocampus tissues harvested. Concentrations of iron, copper and zinc were measured using X-ray fluorescence technology. The results showed that, DID increased iron overall across all strains, ranging from 3 to 68%. Copper and Zinc both decreased, ranging from 0.4-42 and 5-35% respectively. Analysis of variance revealed significant strain by treatment interactions for all three metals. Additionally, in the DID group, we observed strain differences in reduction of hippocampus mass. These findings are particularly interesting to us because high alcohol consumption in humans has been associated with neurodegeneration and dementia related to disruption of iron regulation. The findings of alcohol consumption associated decreases in copper and zinc are novel. The role of copper regulation and neurological function related to alcohol consumption is as yet largely unexplored. The role of zinc is better known as a neuromodulator in the hippocampus and appears to be protective against neurological damage. It would seem then, that the alcohol-related decrease in zinc in the hippocampus would be of concern and warrants further study.

Advances in mouse genetics for the study of human disease

The mouse is the pre-eminent model organism for studies of mammalian gene function and has provided an extraordinarily rich range of insights into basic genetic mechanisms and biological systems. Over several decades, the characterization of mouse mutants has illuminated the relationship between gene and phenotype, providing transformational insights into the genetic bases of disease. However, if we are to deliver the promise of genomic and precision medicine, we must develop a comprehensive catalogue of mammalian gene function that uncovers the dark genome and elucidates pleiotropy. Advances in large-scale mouse mutagenesis programmes allied to high-throughput mouse phenomics are now addressing this challenge and systematically revealing novel gene function and multi-morbidities. Alongside the development of these pan-genomic mutational resources, mouse genetics is employing a range of diversity resources to delineate gene-gene and gene-environment interactions and to explore genetic context. Critically, mouse genetics is a powerful tool for assessing the functional impact of human genetic variation and determining the causal relationship between variant and disease. Together these approaches provide unique opportunities to dissect in vivo mechanisms and systems to understand pathophysiology and disease. Moreover, the provision and utility of mouse models of disease has flourished and engages cumulatively at numerous points across the translational spectrum from basic mechanistic studies to pre-clinical studies, target discovery and therapeutic development.

Strain and sex dependent effects of isolation housing relative to environmental enrichment on operant sensation seeking in mice

Sensation seeking is a multidimensional phenotype that predicts the development of drug addiction in humans and addiction-like drug seeking in rodents. Several lines of evidence suggest that chronic stress increases sensation seeking and addiction-like drug seeking through common genetic mechanisms. Discovery and characterization of these mechanisms would reveal how chronic stress interacts with the genome to influence sensation seeking and how drugs of abuse hijack these fundamental reward mechanisms to drive addiction. To this end, we tested the hypothesis that chronic isolation housing stress (relative to environmental enrichment) influences operant sensation seeking as a function of strain, sex, or their interaction. To determine if the BXD recombinant inbred panel could be used to identify genetic and epigenetic mechanisms underlying any identified gene-by-environment interactions, we used mice from the two BXD founder strains. Following 10 weeks of differential housing, we assessed operant sensation seeking using several reinforcement schedules. The primary finding from this study was that DBA/2J but not C57BL/6J mice were significantly vulnerable to an isolation-induced increase (relative to environmental enrichment) in sensation seeking during extinction when the sensory reward was no longer available; this effect was significantly more robust in females. These data reveal a previously unknown isolation-induced effect on extinction of operant sensation seeking that is sex-dependent, addiction-relevant, and that can be dissected using the BXD recombinant inbred panel.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Quantitative trait locus mapping identifies a locus linked to striatal dopamine and points to collagen IV alpha-6 chain as a novel regulator of striatal axonal branching in mice

Dopaminergic neurons (DA neurons) are controlled by multiple factors, many involved in neurological disease. Parkinson's disease motor symptoms are caused by the demise of nigral DA neurons, leading to loss of striatal dopamine (DA). Here, we measured DA concentration in the dorsal striatum of 32 members of Collaborative Cross (CC) family and their eight founder strains. Striatal DA varied greatly in founders, and differences were highly heritable in the inbred CC progeny. We identified a locus, containing 164 genes, linked to DA concentration in the dorsal striatum on chromosome X. We used RNAseq profiling of the ventral midbrain of two founders with substantial difference in striatal DA-C56BL/6 J and A/J-to highlight potential protein-coding candidates modulating this trait. Among the five differentially expressed genes within the locus, we found that the gene coding for the collagen IV alpha 6 chain (Col4a6) was expressed nine times less in A/J than in C57BL/6J. Using single cell RNA-seq data from developing human midbrain, we found that COL4A6 is highly expressed in radial glia-like cells and neuronal progenitors, indicating a role in neuronal development. Collagen IV alpha-6 chain (COL4A6) controls axogenesis in simple model organisms. Consistent with these findings, A/J mice had less striatal axonal branching than C57BL/6J mice. We tentatively conclude that DA concentration and axonal branching in dorsal striatum are modulated by COL4A6, possibly during development. Our study shows that genetic mapping based on an easily measured Central Nervous System (CNS) trait, using the CC population, combined with follow-up observations, can parse heritability of such a trait, and nominate novel functions for commonly expressed proteins.

Genetic Modulation of Initial Sensitivity to Δ9-Tetrahydrocannabinol (THC) Among the BXD Family of Mice

Cannabinoid receptor 1 activation by the major psychoactive component in cannabis, Δ9-tetrahydrocannabinol (THC), produces motor impairments, hypothermia, and analgesia upon acute exposure. In previous work, we demonstrated significant sex and strain differences in acute responses to THC following administration of a single dose (10 mg/kg, i.p.) in C57BL/6J (B6) and DBA/2J (D2) inbred mice. To determine the extent to which these differences are heritable, we quantified acute responses to a single dose of THC (10 mg/kg, i.p.) in males and females from 20 members of the BXD family of inbred strains derived by crossing and inbreeding B6 and D2 mice. Acute THC responses (initial sensitivity) were quantified as changes from baseline for: 1. spontaneous activity in the open field (mobility), 2. body temperature (hypothermia), and 3. tail withdrawal latency to a thermal stimulus (antinociception). Initial sensitivity to the immobilizing, hypothermic, and antinociceptive effects of THC varied substantially across the BXD family. Heritability was highest for mobility and hypothermia traits, indicating that segregating genetic variants modulate initial sensitivity to THC. We identified genomic loci and candidate genes, including Ndufs2, Scp2, Rps6kb1 or P70S6K, Pde4d, and Pten, that may control variation in THC initial sensitivity. We also detected strong correlations between initial responses to THC and legacy phenotypes related to intake or response to other drugs of abuse (cocaine, ethanol, and morphine). Our study demonstrates the feasibility of mapping genes and variants modulating THC responses in the BXDs to systematically define biological processes and liabilities associated with drug use and abuse.

Differing susceptibility of C57BL/6J and DBA/2J mice-parents of the murine BXD family, to severe acute respiratory syndrome coronavirus infection

The ongoing coronavirus disease-2019 (COVID-19) pandemic, caused by a novel coronavirus termed severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) that is closely related to SARS-CoV, poses a grave threat to global health and has devastated societies worldwide. One puzzling aspect of COVID-19 is the impressive variation in disease manifestations among infected individuals, from a majority who are asymptomatic or exhibit mild symptoms to a smaller, largely age-dependent fraction who develop life-threatening conditions. Some of these differences are likely the consequence of host genetic factors. Systems genetics using diverse and replicable cohorts of isogenic mice represents a powerful way to dissect those host genetic differences that modulate microbial infections. Here we report that the two founders of the large BXD family of mice—C57BL/6J and DBA/2J, differ substantially in their susceptibility to a mouse-adapted SARS-CoV, MA15. Following intranasal viral challenge, DBA/2J develops a more severe disease than C57BL/6J as evidenced by more pronounced and sustained weight loss. Disease was accompanied by high levels of pulmonary viral replication in both strains early after infection but substantially delayed viral clearance in DBA/2J. Our data reveal that the parents of the BXD family are segregated by clear phenotypic differences during MA15 infection and support the feasibility of using this family to systemically dissect the complex virus-host interactions that modulate disease progression and outcome of infection with SARS-CoV, and provisionally also with SARS-CoV-2.

Lifespan-Associated Gene Expression Signatures of Recombinant BXD Mice Implicates Coro7 and Set in Longevity

Although genetic approaches have identified key genes and pathways that promote longevity, systems-level approaches are less utilized. Here, we took advantage of the wealth of omics data characterizing the BXD family of mice. We associated transcript and peptide levels across five tissues from both female and male BXD isogenic lines with their median lifespan. We identified over 5000 genes that showed a longevity correlation in a given tissue. Surprisingly, we found less than 1% overlap among longevity-correlating genes across tissues and sex. These 1% shared genes consist of 51 genes, of which 13 have been shown to alter lifespan. Only two genes -Coro7 and Set- showed a longevity correlation in all tissues and in both sexes. While differential regulation of aging across tissues and sex has been reported, our systems-level analysis reveals two unique genes that may promote healthy aging in unique sex- and tissue-agnostic manner.

Integration of evidence across human and model organism studies: A meeting report

The National Institute on Drug Abuse and Joint Institute for Biological Sciences at the Oak Ridge National Laboratory hosted a meeting attended by a diverse group of scientists with expertise in substance use disorders (SUDs), computational biology, and FAIR (Findability, Accessibility, Interoperability, and Reusability) data sharing. The meeting's objective was to discuss and evaluate better strategies to integrate genetic, epigenetic, and 'omics data across human and model organisms to achieve deeper mechanistic insight into SUDs. Specific topics were to (a) evaluate the current state of substance use genetics and genomics research and fundamental gaps, (b) identify opportunities and challenges of integration and sharing across species and data types, (c) identify current tools and resources for integration of genetic, epigenetic, and phenotypic data, (d) discuss steps and impediment related to data integration, and (e) outline future steps to support more effective collaboration-particularly between animal model research communities and human genetics and clinical research teams. This review summarizes key facets of this catalytic discussion with a focus on new opportunities and gaps in resources and knowledge on SUDs.

In silico candidate variant and gene identification using inbred mouse strains

Mice are the most widely used animal model to study genotype to phenotype relationships. Inbred mice are genetically identical, which eliminates genetic heterogeneity and makes them particularly useful for genetic studies. Many different strains have been bred over decades and a vast amount of phenotypic data has been generated. In addition, recently whole genome sequencing-based genome-wide genotype data for many widely used inbred strains has been released. Here, we present an approach for in silico fine-mapping that uses genotypic data of 37 inbred mouse strains together with phenotypic data provided by the user to propose candidate variants and genes for the phenotype under study. Public genome-wide genotype data covering more than 74 million variant sites is queried efficiently in real-time to provide those variants that are compatible with the observed phenotype differences between strains. Variants can be filtered by molecular consequences and by corresponding molecular impact. Candidate gene lists can be generated from variant lists on the fly. Fine-mapping together with annotation or filtering of results is provided in a Bioconductor package called MouseFM. In order to characterize candidate variant lists under various settings, MouseFM was applied to two expression data sets across 20 inbred mouse strains, one from neutrophils and one from CD4⁺ T cells. Fine-mapping was assessed for about 10,000 genes, respectively, and identified candidate variants and haplotypes for many expression quantitative trait loci (eQTLs) reported previously based on these data. For albinism, MouseFM reports only one variant allele of moderate or high molecular impact that only albino mice share: a missense variant in the Tyr gene, reported previously to be causal for this phenotype. Performing in silico fine-mapping for interfrontal bone formation in mice using four strains with and five strains without interfrontal bone results in 12 genes. Of these, three are related to skull shaping abnormality. Finally performing fine-mapping for dystrophic cardiac calcification by comparing 9 strains showing the phenotype with eight strains lacking it, we identify only one moderate impact variant in the known causal gene Abcc6. In summary, this illustrates the benefit of using MouseFM for candidate variant and gene identification.