High-throughput sequencing of the DBA/2J mouse genome

The DBA/2J mouse is not only the oldest inbred strain, but also one of the most widely used strains. DBA/2J exhibits many unique anatomical, physiological, and behavior traits. In addition, DBA/2J is one parent of the large BXD family of recombinant inbred strains [1]. The genome of the other parent of this BXD family— C57BL/6J—has been sequenced and serves as the mouse reference genome [2]. We sequenced the genome of DBA/2J using SOLiD and Illumina high throughput short read protocols to generate a comprehensive set of ~5 million sequence variants segregating in the BXD family that ultimately cause developmental, anatomical, functional and behavioral differences among these 80+ strains

Genetic regulation of Nrnx1 expression: an integrative cross-species analysis of schizophrenia candidate genes

Neurexin 1 (NRXN1) is a large presynaptic transmembrane protein that has complex and variable patterns of expression in the brain. Sequence variants in NRXN1 are associated with differences in cognition, and with schizophrenia and autism. The murine Nrxn1 gene is also highly polymorphic and is associated with significant variation in expression that is under strong genetic control. Here, we use co-expression analysis, high coverage genomic sequence, and expression quantitative trait locus (eQTL) mapping to study the regulation of this gene in the brain. We profiled a family of 72 isogenic progeny strains of a cross between C57BL/6J and DBA/2J (the BXD family) using exon arrays and massively parallel RNA sequencing. Expression of most Nrxn1 exons have high genetic correlation (r>0.6) because of the segregation of a common trans eQTL on chromosome (Chr) 8 and a common cis eQTL on Chr 17. These two loci are also linked to murine phenotypes relevant to schizophrenia and to a novel human schizophrenia candidate gene with high neuronal expression (Pleckstrin and Sec7 domain containing 3). In both human and mice, NRXN1 is co-expressed with numerous synaptic and cell signaling genes, and known schizophrenia candidates. Cross-species co-expression and protein interaction network analyses identified glycogen synthase kinase 3 beta (GSK3B) as one of the most consistent and conserved covariates of NRXN1. By using the Molecular Genetics of Schizophrenia data set, we were able to test and confirm that markers in NRXN1 and GSK3B have epistatic interactions in human populations that can jointly modulate risk of schizophrenia

Integrative Analysis of Low- and High-Resolution eQTL

The study of expression quantitative trait loci (eQTL) is a powerful way of detecting transcriptional regulators at a genomic scale and for elucidating how natural genetic variation impacts gene expression. Power and genetic resolution are heavily affected by the study population: whereas recombinant inbred (RI) strains yield greater statistical power with low genetic resolution, using diverse inbred or outbred strains improves genetic resolution at the cost of lower power. In order to overcome the limitations of both individual approaches, we combine data from RI strains with genetically more diverse strains and analyze hippocampus eQTL data obtained from mouse RI strains (BXD) and from a panel of diverse inbred strains (Mouse Diversity Panel, MDP). We perform a systematic analysis of the consistency of eQTL independently obtained from these two populations and demonstrate that a significant fraction of eQTL can be replicated. Based on existing knowledge from pathway databases we assess different approaches for using the high-resolution MDP data for fine mapping BXD eQTL. Finally, we apply this framework to an eQTL hotspot on chromosome 1 (Qrr1), which has been implicated in a range of neurological traits. Here we present the first systematic examination of the consistency between eQTL obtained independently from the BXD and MDP populations. Our analysis of fine-mapping approaches is based on ‘real life’ data as opposed to simulated data and it allows us to propose a strategy for using MDP data to fine map BXD eQTL. Application of this framework to Qrr1 reveals that this eQTL hotspot is not caused by just one (or few) ‘master regulators’, but actually by a set of polymorphic genes specific to the central nervous system

Quantitative and Qualitative Stem Rust Resistance Factors in Barley Are Associated with Transcriptional Suppression of Defense Regulons

Stem rust (Puccinia graminis f. sp. tritici; Pgt) is a devastating fungal disease of wheat and barley. Pgt race TTKSK (isolate Ug99) is a serious threat to these Triticeae grain crops because resistance is rare. In barley, the complex Rpg-TTKSK locus on chromosome 5H is presently the only known source of qualitative resistance to this aggressive Pgt race. Segregation for resistance observed on seedlings of the Q21861 × SM89010 (QSM) doubled-haploid (DH) population was found to be predominantly qualitative, with little of the remaining variance explained by loci other than Rpg-TTKSK. In contrast, analysis of adult QSM DH plants infected by field inoculum of Pgt race TTKSK in Njoro, Kenya, revealed several additional quantitative trait loci that contribute to resistance. To molecularly characterize these loci, Barley1 GeneChips were used to measure the expression of 22,792 genes in the QSM population after inoculation with Pgt race TTKSK or mock-inoculation. Comparison of expression Quantitative Trait Loci (eQTL) between treatments revealed an inoculation-dependent expression polymorphism implicating Actin depolymerizing factor3 (within the Rpg-TTKSK locus) as a candidate susceptibility gene. In parallel, we identified a chromosome 2H trans-eQTL hotspot that co-segregates with an enhancer of Rpg-TTKSK-mediated, adult plant resistance discovered through the Njoro field trials. Our genome-wide eQTL studies demonstrate that transcript accumulation of 25% of barley genes is altered following challenge by Pgt race TTKSK, but that few of these genes are regulated by the qualitative Rpg-TTKSK on chromosome 5H. It is instead the chromosome 2H trans-eQTL hotspot that orchestrates the largest inoculation-specific responses, where enhanced resistance is associated with transcriptional suppression of hundreds of genes scattered throughout the genome. Hence, the present study associates the early suppression of genes expressed in this host–pathogen interaction with enhancement of R-gene mediated resistance

Quantitative trait loci for sensitivity to ethanol intoxication in a C57BL/6J × 129S1/SvImJ inbred mouse cross

Individual variation in sensitivity to acute ethanol (EtOH) challenge is associated with alcohol drinking and is a predictor of alcohol abuse. Previous studies have shown that the C57BL/6J (B6) and 129S1/SvImJ (S1) inbred mouse strains differ in responses on certain measures of acute EtOH intoxication. To gain insight into genetic factors contributing to these differences, we performed quantitative trait locus (QTL) analysis of measures of EtOH-induced ataxia (accelerating rotarod), hypothermia, and loss of righting reflex (LORR) duration in a B6 × S1 F2 population. We confirmed that S1 showed greater EtOH-induced hypothermia (specifically at a high dose) and longer LORR compared to B6. QTL analysis revealed several additive and interacting loci for various phenotypes, as well as examples of genotype interactions with sex. QTLs for different EtOH phenotypes were largely non-overlapping, suggesting separable genetic influences on these behaviors. The most compelling main-effect QTLs were for hypothermia on chromosome 16 and for LORR on chromosomes 4 and 6. Several QTLs overlapped with loci repeatedly linked to EtOH drinking in previous mouse studies. The architecture of the traits we examined was complex but clearly amenable to dissection in future studies. Using integrative genomics strategies, plausible functional and positional candidates may be found. Uncovering candidate genes associated with variation in these phenotypes in this population could ultimately shed light on genetic factors underlying sensitivity to EtOH intoxication and risk for alcoholism in humans

Reward-Related Behavioral Paradigms for Addiction Research in the Mouse: Performance of Common Inbred Strains

The mouse has emerged as a uniquely valuable species for studying the molecular and genetic basis of complex behaviors and modeling neuropsychiatric disease states. While valid and reliable preclinical assays for reward-related behaviors are critical to understanding addiction-related processes, and various behavioral procedures have been developed and characterized in rats and primates, there have been relatively few studies using operant-based addiction-relevant behavioral paradigms in the mouse. Here we describe the performance of the C57BL/6J inbred mouse strain on three major reward-related paradigms, and replicate the same procedures in two other commonly used inbred strains (DBA/2J, BALB/cJ). We examined Pavlovian-instrumental transfer (PIT) by measuring the ability of an auditory cue associated with food reward to promote an instrumental (lever press) response. In a separate experiment, we assessed the acquisition and extinction of a simple stimulus-reward instrumental behavior on a touchscreen-based task. Reinstatement of this behavior was then examined following either continuous exposure to cues (conditioned reinforcers, CRs) associated with reward, brief reward and CR exposure, or brief reward exposure followed by continuous CR exposure. The third paradigm examined sensitivity of an instrumental (lever press) response to devaluation of food reward (a probe for outcome insensitive, habitual behavior) by repeated pairing with malaise. Results showed that C57BL/6J mice displayed robust PIT, as well as clear extinction and reinstatement, but were insensitive to reinforcer devaluation. DBA/2J mice showed good PIT and (rewarded) reinstatement, but were slow to extinguish and did not show reinforcer devaluation or significant CR-reinstatement. BALB/cJ mice also displayed good PIT, extinction and reinstatement, and retained instrumental responding following devaluation, but, unlike the other strains, demonstrated reduced Pavlovian approach behavior (food magazine head entries). Overall, these assays provide robust paradigms for future studies using the mouse to elucidate the neural, molecular and genetic factors underpinning reward-related behaviors relevant to addiction research

Genetic Control of a Central Pattern Generator: Rhythmic Oromotor Movement in Mice Is Controlled by a Major Locus near Atp1a2

Fluid licking in mice is a rhythmic behavior that is controlled by a central pattern generator (CPG) located in a complex of brainstem nuclei. C57BL/6J (B6) and DBA/2J (D2) strains differ significantly in water-restricted licking, with a highly heritable difference in rates (h2≥0.62) and a corresponding 20% difference in interlick interval (mean ± SEM = 116.3±1 vs 95.4±1.1 ms). We systematically quantified motor output in these strains, their F1 hybrids, and a set of 64 BXD progeny strains. The mean primary interlick interval (MPI) varied continuously among progeny strains. We detected a significant quantitative trait locus (QTL) for a CPG controlling lick rate on Chr 1 (Lick1), and a suggestive locus on Chr 10 (Lick10). Linkage was verified by testing of B6.D2-1D congenic stock in which a segment of Chr 1 of the D2 strain was introgressed onto the B6 parent. The Lick1 interval on distal Chr 1 contains several strong candidate genes. One of these is a sodium/potassium pump subunit (Atp1a2) with widespread expression in astrocytes, as well as in a restricted population of neurons. Both this subunit and the entire Na+/K+-ATPase molecule have been implicated in rhythmogenesis for respiration and locomotion. Sequence variants in or near Apt1a2 strongly modulate expression of the cognate mRNA in multiple brain regions. This gene region has recently been sequenced exhaustively and we have cataloged over 300 non-coding and synonymous mutations segregating among BXD strains, one or more of which is likely to contribute to differences in central pattern generator tempo

Treatment- and Population-Dependent Activity Patterns of Behavioral and Expression QTLs

Genetic control of gene expression and higher-order phenotypes is almost invariably dependent on environment and experimental conditions. We use two families of recombinant inbred strains of mice (LXS and BXD) to study treatment- and genotype-dependent control of hippocampal gene expression and behavioral phenotypes. We analyzed responses to all combinations of two experimental perturbations, ethanol and restraint stress, in both families, allowing for comparisons across 8 combinations of treatment and population. We introduce the concept of QTL activity patterns to characterize how associations between genomic loci and traits vary across treatments. We identified several significant behavioral QTLs and many expression QTLs (eQTLs). The behavioral QTLs are highly dependent on treatment and population. We classified eQTLs into three groups: cis-eQTLs (expression variation that maps to within 5 Mb of the cognate gene), syntenic trans-eQTLs (the gene and the QTL are on the same chromosome but not within 5 Mb), and non-syntenic trans-eQTLs (the gene and the QTL are on different chromosomes). We found that most non-syntenic trans-eQTLs were treatment-specific whereas both classes of syntenic eQTLs were more conserved across treatments. We also found there was a correlation between regions along the genome enriched for eQTLs and SNPs that were conserved across the LXS and BXD families. Genes with eQTLs that co-localized with the behavioral QTLs and displayed similar QTL activity patterns were identified as potential candidate genes associated with the phenotypes, yielding identification of novel genes as well as genes that have been previously associated with responses to ethanol