Systems Genetic Analysis in GeneNetwork.org

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Genome-wide screen identifies host loci that modulate Mycobacterium tuberculosis fitness in immunodivergent mice

Genetic differences among mammalian hosts and among strains of Mycobacterium tuberculosis (Mtb) are well-established determinants of tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host-pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the highly diverse BXD family of strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb-resistant C57BL/6J (B6 or B) and Mtb-susceptible DBA/2J (D2 or D) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters of "endophenotypes," each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted quantitative trait loci (QTL) mapping of these bacterial fitness endophenotypes and identified 140 host-pathogen QTL (hpQTL). We located a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes: Rv0127 (mak), Rv0359 (rip2), Rv0955 (perM), and Rv3849 (espR). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

GenomeMUSter mouse genetic variation service enables multi-trait, multi-population data integration and analyses

Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize genetic variants that contribute to disease. Thousands of disease relevant traits have been characterized in mice and made publicly available. New strains and populations including the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait - variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web- and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Effects of isolation housing stress and mouse strain on intravenous cocaine self-administration, sensory stimulus self-administration, and reward preference

Sensory stimuli are natural rewards in mice and humans. Consequently, preference for a drug reward relative to a sensory reward may be an endophenotype of addiction vulnerability. In this study, we developed a novel behavioral assay to quantify the preference for intravenous drug self-administration relative to sensory stimulus self-administration. We used founder strains of the BXD recombinant inbred mouse panel (C57BL/6J, DBA/2J) and a model of stress (isolation vs enriched housing) to assess genetic and epigenetic effects. Following 10 weeks of differential housing, all mice were tested under three reward conditions: sensory rewards available, cocaine rewards available, both rewards available. When a single reward was available (sensory stimuli or cocaine; delivered using distinct levers), DBA/2J mice self-administered significantly more rewards than C57BL/6J mice. When both rewards were available, DBA/2J mice exhibited a significant preference for cocaine relative to sensory stimuli; in contrast, C57BL/6J mice exhibited no preference. Housing condition influenced sensory stimulus self-administration and strain-dependently influenced inactive lever pressing when both rewards were available. Collectively, these data reveal strain effects, housing effects, or both on reward self-administration and preference. Most importantly, this study reveals that genetic mechanisms underlying preference for a drug reward relative to a nondrug reward can be dissected using the full BXD panel.

Intravenous fentanyl self-administration in male and female C57BL/6J and DBA/2J mice

The genetic mechanisms underlying fentanyl addiction, a highly heritable disease, are unknown. Identifying these mechanisms will lead to better risk assessment, early diagnosis, and improved intervention. To this end, we used intravenous fentanyl self-administration to quantify classical self-administration phenotypes and addiction-like fentanyl seeking in male and female mice from the two founder strains of the BXD recombinant inbred mouse panel (C57BL/6J and DBA/2J). We reached three primary conclusions from these experiments. First, mice from all groups rapidly acquired intravenous fentanyl self-administration and exhibited a dose-response curve, extinction burst, and extinction of the learned self-administration response. Second, fentanyl intake (during acquisition and dose response) and fentanyl seeking (during extinction) were equivalent among groups. Third, strain effects, sex effects, or both were identified for several addiction-like behaviors (cue-induced reinstatement, stress-induced reinstatement, escalation of intravenous fentanyl self-administration). Collectively, these data indicate that C57BL/6J and DBA/2J mice of both sexes were able to acquire, regulate, and extinguish intravenous fentanyl self-administration. Moreover, these data reveal novel strain and sex effects on addiction-like behaviors in the context of intravenous fentanyl self-administration in mice and indicate that the full BXD panel can be used to identify and dissect the genetic mechanisms underlying these effects.

Evidence for a neuromuscular circuit involving hypothalamic interleukin-6 in the control of skeletal muscle metabolism

Hypothalamic interleukin-6 (IL6) exerts a broad metabolic control. Here, we demonstrated that IL6 activates the ERK1/2 pathway in the ventromedial hypothalamus (VMH), stimulating AMPK/ACC signaling and fatty acid oxidation in mouse skeletal muscle. Bioinformatics analysis revealed that the hypothalamic IL6/ERK1/2 axis is closely associated with fatty acid oxidation- and mitochondrial-related genes in the skeletal muscle of isogenic BXD mouse strains and humans. We showed that the hypothalamic IL6/ERK1/2 pathway requires the α2-adrenergic pathway to modify fatty acid skeletal muscle metabolism. To address the physiological relevance of these findings, we demonstrated that this neuromuscular circuit is required to underpin AMPK/ACC signaling activation and fatty acid oxidation after exercise. Last, the selective down-regulation of IL6 receptor in VMH abolished the effects of exercise to sustain AMPK and ACC phosphorylation and fatty acid oxidation in the muscle after exercise. Together, these data demonstrated that the IL6/ERK axis in VMH controls fatty acid metabolism in the skeletal muscle.

TLR4 deletion increases basal energy expenditure and attenuates heart apoptosis and ER stress but mitigates the training-induced cardiac function and performance improvement

Strategies capable of attenuating TLR4 can attenuate metabolic processes such as inflammation, endoplasmic reticulum (ER) stress, and apoptosis in the body. Physical exercise has been a cornerstone in suppressing inflammation and dysmetabolic outcomes caused by TRL4 activation. Thus, the present study aimed to evaluate the effects of a chronic physical exercise protocol on the TLR4 expression and its repercussion in the inflammation, ER stress, and apoptosis pathways in mice hearts. Echocardiogram, RT-qPCR, immunoblotting, and histological techniques were used to evaluate the left ventricle of wild-type (WT) and Tlr4 knockout (TLR4 KO) mice submitted to a 4-week physical exercise protocol. Moreover, we performed a bioinformatics analysis to expand the relationship of Tlr4 mRNA in the heart with inflammation, ER stress, and apoptosis-related genes of several isogenic strains of BXD mice. The TLR4 KO mice had higher energy expenditure and heart rate in the control state but lower activation of apoptosis and ER stress pathways. The bioinformatics analysis reinforced these data. In the exercised state, the WT mice improved performance and cardiac function. However, these responses were blunted in the KO group. In conclusion, TLR4 has an essential role in the inhibition of apoptosis and ER stress pathways, as well as in the training-induced beneficial adaptations.

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.

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.

Systems genetic analysis of nicotine withdrawal deficits in hippocampus-dependent learning

Cognitive deficits, such as disrupted learning, are a major symptom of nicotine withdrawal. These deficits are heritable, yet their genetic basis is largely unknown. Our lab has developed a mouse model of nicotine withdrawal deficits in learning, using chronic nicotine exposure via osmotic minipumps and fear conditioning. Here, we utilized the BXD genetic reference panel to identify genetic variants underlying nicotine withdrawal deficits in learning. Male and female mice (n = 6-11 per sex per strain, 31 strains) received either chronic saline or nicotine (6.3 mg/kg per day for 12 days), and were then tested for hippocampus-dependent learning deficits using contextual fear conditioning. Quantitative trait locus (QTL) mapping analyses using GeneNetwork identified a significant QTL on Chromosome 4 (82.13 Mb, LRS = 20.03, p < 0.05). Publicly available hippocampal gene expression data were used to identify eight positional candidates (Snacpc3, Mysm1, Rps6, Plaa, Lurap1l, Slc24a2, Hacd4, Ptprd) that overlapped with our behavioral QTL and correlated with our behavioral data. Overall, this study demonstrates that genetic factors impact cognitive deficits during nicotine withdrawal in the BXD recombinant inbred panel and identifies candidate genes for future research.

Stimulus Complexity and Mouse Strain Drive Escalation of Operant Sensation Seeking Within and Across Sessions in C57BL/6J and DBA/2J Mice

Sensation seeking is a heritable trait that is genetically correlated with substance use; the shared genetic mechanisms underlying these traits are largely unknown. The relationship of sensation seeking and substance use has practical importance because discovering genes that drive sensation seeking can reveal genes driving substance use, and quantification of sensation seeking in mice is higher throughput and less technically challenging than quantification of volitional drug use. In order to fully understand the genetic mechanisms driving sensation seeking, it is critical to first understand the nongenetic factors driving sensation seeking. In the present study, we used the operant sensation seeking paradigm to assess the effects of stimulus complexity on sensation seeking in C57BL/6J and DBA/2J mice. These strains are the founders of the BXD recombinant inbred mouse panel which enables the discovery of genes driving phenotypic variation. This study led to four principal conclusions. First, all sensory stimuli used in the study, regardless of complexity or number of stimulus modalities, served as reinforcers for C57BL/6J and DBA/2J mice. Second, for both C57BL/6J and DBA/2J mice, sensation seeking for a high complexity sensory stimulus was significantly greater than sensation seeking for a low complexity sensory stimulus. Third, for both C57BL/6J and DBA/2J mice, sensation seeking escalated significantly within-session when a multimodal sensory stimulus of medium or high complexity was used but not when a unimodal sensory stimulus of low complexity was used. Finally, both the magnitude of sensation seeking and the magnitude of within-session escalation of sensation seeking were significantly greater in mice from the DBA/2J strain relative to mice from the C57BL/6J strain. Collectively, these findings indicate that stimulus complexity and genetic background drive escalation of operant sensation seeking within and across sessions, and that the BXD recombinant inbred mouse panel can be used to discover the genetic mechanisms underlying these phenomena.

Different Effect of Sox11 in Retinal Ganglion Cells Survival and Axon Regeneration

Purpose: The present study examines the role of Sox11 in the initial response of retinal ganglion cells (RGCs) to axon damage and in optic nerve regeneration in mouse. Methods: Markers of retinal injury were identified using the normal retina database and optic nerve crush (ONC) database on GeneNetwork² (www.genenetwork.org). One gene, Sox11, was highly upregulated following ONC. We examined the role of this transcription factor, Sox11, following ONC and optic nerve regeneration in mice. In situ hybridization was performed using the Affymetrix 2-plex Quantigene View RNA In Situ Hybridization Tissue Assay System. Sox11 was partially knocked out by intravitreal injection of AAV2-CMV-Cre-GFP in Sox11 ^f/f mice. Optic nerve regeneration model used Pten knockdown. Mice were perfused and the retinas and optic nerves were dissected and examined for RGC survival and axon growth. Results: Sox11 was dramatically upregulated in the retina following ONC injury. The level of Sox11 message increased by approximately eightfold 2 days after ONC. In situ hybridization demonstrated low-level Sox11 message in RGCs and cells in the inner nuclear layer in the normal retina as well as a profound increase in Sox11 message within the ganglion cells following ONC. In Sox11 ^f/f retinas, partially knocking out Sox11 significantly increased RGC survival after ONC as compared to the AAV2-CMV-GFP control group; however, it had little effect on the ability of axon regeneration. Combinatorial downregulation of both Sox11 and Pten resulted in a significant increase in RGC survival as compared to Pten knockdown only. When Pten was knocked down there was a remarkable increase in the number and the length of regenerating axons. Partially knocking out Sox11 in combination with Pten deletion resulted in a fewer regenerating axons. Conclusion: Taken together, these data demonstrate that Sox11 is involved in the initial response of the retina to injury, playing a role in the early attempts of axon regeneration and neuronal survival. Downregulation of Sox11 aids in RGC survival following injury of optic nerve axons, while a partial knockout of Sox11 negates the axon regeneration stimulated by Pten knockdown.

Systems genetics of sensation seeking

Sensation seeking is a multifaceted, heritable trait which predicts the development of substance use and abuse in humans; similar phenomena have been observed in rodents. Genetic correlations among sensation seeking and substance use indicate shared biological mechanisms, but the genes and networks underlying these relationships remain elusive. Here, we used a systems genetics approach in the BXD recombinant inbred mouse panel to identify shared genetic mechanisms underlying substance use and preference for sensory stimuli, an intermediate phenotype of sensation seeking. Using the operant sensation seeking (OSS) paradigm, we quantified preference for sensory stimuli in 120 male and 127 female mice from 62 BXD strains and the C57BL/6J and DBA/2J founder strains. We used relative preference for the active and inactive levers to dissociate preference for sensory stimuli from locomotion and exploration phenotypes. We identified genomic regions on chromosome 4 (155.236‐155.742 Mb) and chromosome 13 (72.969‐89.423 Mb) associated with distinct behavioral components of OSS. Using publicly available behavioral data and mRNA expression data from brain regions involved in reward processing, we identified (a) genes within these behavioral QTL exhibiting genome‐wide significant cis‐eQTL and (b) genetic correlations among OSS phenotypes, ethanol phenotypes and mRNA expression. From these analyses, we nominated positional candidates for behavioral QTL associated with distinct OSS phenotypes including Gnb1 and Mef2c. Genetic covariation of Gnb1 expression, preference for sensory stimuli and multiple ethanol phenotypes suggest that heritable variation in Gnb1 expression in reward circuitry partially underlies the widely reported relationship between sensation seeking and substance use.

A quantitative trait locus on chromosome 1 modulates intermale aggression in mice

Aggression between male conspecifics is a complex social behavior that is likely modulated by multiple gene variants. In this study, the BXD recombinant inbred mouse strains (RIS) were used to map quantitative trait loci (QTLs) underlying behaviors associated with intermale aggression. Four hundred and fifty-seven males from 55 strains (including the parentals) were observed at an age of 13 ± 1 week in a resident-intruder test following 10 days of isolation. Attack latency was measured directly within a 10-minute time period and the test was repeated 24 hours later. The variables we analyzed were the proportion of attacking males in a given strain as well as the attack latency (on days 1 and 2, and both days combined). On day 1, 29% of males attacked, and this increased to 37% on day 2. Large strain differences were obtained for all measures of aggression, indicating substantial heritability (intraclass correlations 0.10-0.18). We identified a significant QTL on chromosome (Chr) 1 and suggestive QTLs on mouse Chrs 1 and 12 for both attack and latency variables. The significant Chr 1 locus maps to a gene-sparse region between 82 and 88.5 Mb with the C57BL/6J allele increasing aggression and explaining about 18% of the variance. The most likely candidate gene modulating this trait is Htr2b which encodes the serotonin 2B receptor and has been implicated in aggressive and impulsive behavior in mice, humans and other species.

QTL and systems genetics analysis of mouse grooming and behavioral responses to novelty in an open field

The open field is a classic test used to assess exploratory behavior, anxiety and locomotor activity in rodents. Here, we mapped quantitative trait loci (QTLs) underlying behaviors displayed in an open field, using a panel of 53 BXD recombinant inbred mouse strains with deep replication (10 per strain and sex). The use of these strains permits the integration and comparison of data obtained in different laboratories, and also offers the possibility to study trait covariance by exploiting powerful bioinformatics tools and resources. We quantified behavioral traits during 20-min test sessions including (1) percent time spent and distance traveled near the wall (thigmotaxis), (2) leaning against the wall, (3) rearing, (4) jumping, (5) grooming duration, (6) grooming frequency, (7) locomotion and (8) defecation. All traits exhibit moderate heritability making them amenable to genetic analysis. We identified a significant QTL on chromosome M.m. 4 at approximately 104 Mb that modulates grooming duration in both males and females (likelihood ratio statistic values of approximately 18, explaining 25% and 14% of the variance, respectively) and a suggestive QTL modulating locomotion that maps to the same locus. Bioinformatic analysis indicates Disabled 1 (Dab1, a key protein in the reelin signaling pathway) as a particularly strong candidate gene modulating these behaviors. We also found 2 highly suggestive QTLs for a sex by strain interaction for grooming duration on chromosomes 13 and 17. In addition, we identified a pairwise epistatic interaction between loci on chromosomes 12 at 36-37 Mb and 14 at 34-36 Mb that influences rearing frequency in males.