Khyobeni Mozhui

Strain Differences in Stress Responsivity Are Associated with Divergent Amygdala Gene Expression and Glutamate-Mediated Neuronal Excitability

Stress is a major risk factor for numerous neuropsychiatric diseases. However, susceptibility to stress and the qualitative nature of stress effects on behavior differ markedly among individuals. This is partly because of the moderating influence of genetic factors. Inbred mouse strains provide a relatively stable and restricted range of genetic and environmental variability that is valuable for disentangling gene–stress interactions. Here, we screened a panel of inbred strains for anxiety- and depression-related phenotypes at baseline (trait) and after exposure to repeated restraint. Two strains, DBA/2J and C57BL/6J, differed in trait and restraint-induced anxiety-related behavior (dark/light exploration, elevated plus maze). Gene expression analysis of amygdala, medial prefrontal cortex, and hippocampus revealed divergent expression in DBA/2J and C57BL/6J both at baseline and after repeated restraint. Restraint produced strain-dependent expression alterations in various genes including glutamate receptors (e.g., Grin1, Grik1). To elucidate neuronal correlates of these strain differences, we performed ex vivo analysis of glutamate excitatory neurotransmission in amygdala principal neurons. Repeated restraint augmented amygdala excitatory postsynaptic signaling and altered metaplasticity (temporal summation of NMDA receptor currents) in DBA/2J but not C57BL/6J. Furthermore, we found that the C57BL/6J-like changes in anxiety-related behavior after restraint were absent in null mutants lacking the modulatory NMDA receptor subunit Grin2a, but not the AMPA receptor subunit Gria1. Grin2a null mutants exhibited significant (∼30%) loss of dendritic spines on amygdala principal neurons under nonrestraint conditions. Collectively, our data support a model in which genetic variation in glutamatergic neuroplasticity in corticolimbic circuitry underlies phenotypic variation in responsivity to stress.

Variation in mouse basolateral amygdala volume is associated with differences in stress reactivity and fear learning.

A wealth of research identifies the amygdala as a key brain region mediating negative affect, and implicates amygdala dysfunction in the pathophysiology of anxiety disorders. Although there is a strong genetic component to anxiety disorders such as posttraumatic stress disorder (PTSD) there remains debate about whether abnormalities in amygdala function predispose to these disorders. In the present study, groups of C57BL/6 × DBA/2 (B × D) recombinant inbred strains of mice were selected for differences in volume of the basolateral amygdala complex (BLA). Strains with relatively small, medium, or large BLA volumes were compared for Pavlovian fear learning and memory, anxiety-related behaviors, depression-related behavior, and glucocorticoid responses to stress. Strains with relatively small BLA exhibited stronger conditioned fear responses to both auditory tone and contextual stimuli, as compared to groups with larger BLA. The small BLA group also showed significantly greater corticosterone responses to stress than the larger BLA groups. BLA volume did not predict clear differences in measures of anxiety-like behavior or depression-related behavior, other than greater locomotor inhibition to novelty in strains with smaller BLA. Neither striatal, hippocampal nor cerebellar volumes correlated significantly with any behavioral measure. The present data demonstrate a phenotype of enhanced fear conditioning and exaggerated glucocorticoid responses to stress associated with small BLA volume. This profile is reminiscent of the increased fear processing and stress reactivity that is associated with amygdala excitability and reduced amygdala volume in humans carrying loss of function polymorphisms in the serotonin transporter and monoamine oxidase A genes. Our study provides a unique example of how natural variation in amygdala volume associates with specific fear- and stress-related phenotypes in rodents, and further supports the role of amygdala dysfunction in anxiety disorders such as PTSD.

Dissection of a QTL Hotspot on Mouse Distal Chromosome 1 that Modulates Neurobehavioral Phenotypes and Gene Expression

A remarkably diverse set of traits maps to a region on mouse distal chromosome 1 (Chr 1) that corresponds to human Chr 1q21–q23. This region is highly enriched in quantitative trait loci (QTLs) that control neural and behavioral phenotypes, including motor behavior, escape latency, emotionality, seizure susceptibility (Szs1), and responses to ethanol, caffeine, pentobarbital, and haloperidol. This region also controls the expression of a remarkably large number of genes, including genes that are associated with some of the classical traits that map to distal Chr 1 (e.g., seizure susceptibility). Here, we ask whether this QTL-rich region on Chr 1 (Qrr1) consists of a single master locus or a mixture of linked, but functionally unrelated, QTLs. To answer this question and to evaluate candidate genes, we generated and analyzed several gene expression, haplotype, and sequence datasets. We exploited six complementary mouse crosses, and combed through 18 expression datasets to determine class membership of genes modulated by Qrr1. Qrr1 can be broadly divided into a proximal part (Qrr1p) and a distal part (Qrr1d), each associated with the expression of distinct subsets of genes. Qrr1d controls RNA metabolism and protein synthesis, including the expression of ∼20 aminoacyl-tRNA synthetases. Qrr1d contains a tRNA cluster, and this is a functionally pertinent candidate for the tRNA synthetases. Rgs7 and Fmn2 are other strong candidates in Qrr1d. FMN2 protein has pronounced expression in neurons, including in the dendrites, and deletion of Fmn2 had a strong effect on the expression of few genes modulated by Qrr1d. Our analysis revealed a highly complex gene expression regulatory interval in Qrr1, composed of multiple loci modulating the expression of functionally cognate sets of genes.

Detection, Validation, and Downstream Analysis of Allelic Variation in Gene Expression

Common sequence variants within a gene often generate important differences in expression of corresponding mRNAs. This high level of local (allelic) control—or cis modulation—rivals that produced by gene targeting, but expression is titrated finely over a range of levels. We are interested in exploiting this allelic variation to study gene function and downstream consequences of differences in expression dosage. We have used several bioinformatics and molecular approaches to estimate error rates in the discovery of cis modulation and to analyze some of the biological and technical confounds that contribute to the variation in gene expression profiling. Our analysis of SNPs and alternative transcripts, combined with eQTL maps and selective gene resequencing, revealed that between 17 and 25% of apparent cis modulation is caused by SNPs that overlap probes rather than by genuine quantitative differences in mRNA levels. This estimate climbs to 40–50% when qualitative differences between isoform variants are included. We have developed an analytical approach to filter differences in expression and improve the yield of genuine cis-modulated transcripts to ∼80%. This improvement is important because the resulting variation can be successfully used to study downstream consequences of altered expression on higher-order phenotypes. Using a systems genetics approach we show that two validated cis-modulated genes, Stk25 and Rasd2, are likely to control expression of downstream targets and affect disease susceptibility.

A promoter polymorphism in the Per3 gene is associated with alcohol and stress response

The period homolog genes Per1, Per2 and Per3 are important components of the circadian clock system. In addition to their role in maintaining circadian rhythm, these genes have been linked to mood disorders, stress response and vulnerability to addiction and alcoholism. In this study, we combined high-resolution sequence analysis and quantitative trait locus (QTL) mapping of gene expression and behavioral traits to identify Per3 as a compelling candidate for the interaction between circadian rhythm, alcohol and stress response. In the BXD family of mouse strains, sequence variants in Per3 have marked effects on steady-state mRNA and protein levels. As a result, the transcript maps as a cis-acting expression QTL (eQTL). We found that an insertion/deletion (indel) variant in a putative stress response element in the promoter region of Per3 causes local control of transcript abundance. This indel results in differences in protein binding affinities between the two alleles through the Nrf2 transcriptional activator. Variation in Per3 is also associated with downstream differences in the expression of genes involved in circadian rhythm, alcohol, stress response and schizophrenia. We found that the Per3 locus is linked to stress/anxiety traits, and that the basal expression of Per3 is also correlated with several anxiety and addiction-related phenotypes. Treatment with alcohol results in increased expression of Per3 in the hippocampus, and this effect interacts with acute restraint stress. Our data provide strong evidence that variation in the Per3 transcript is causally associated with and also responsive to stress and alcohol.

Joint mouse-human phenome-wide association to test gene function and disease risk

Phenome-wide association is a novel reverse genetic strategy to analyze genome-to-phenome relations in human clinical cohorts. Here we test this approach using a large murine population segregating for ∼5 million sequence variants, and we compare our results to those extracted from a matched analysis of gene variants in a large human cohort. For the mouse cohort, we amassed a deep and broad open-access phenome consisting of ∼4,500 metabolic, physiological, pharmacological and behavioural traits, and more than 90 independent expression quantitative trait locus (QTL), transcriptome, proteome, metagenome and metabolome data sets—by far the largest coherent phenome for any experimental cohort (www.genenetwork.org). We tested downstream effects of subsets of variants and discovered several novel associations, including a missense mutation in fumarate hydratase that controls variation in the mitochondrial unfolded protein response in both mouse and Caenorhabditis elegans, and missense mutations in Col6a5 that underlies variation in bone mineral density in both mouse and human.

Complex Control of GABA(A) Receptor Subunit mRNA Expression: Variation, Covariation, and Genetic Regulation

GABA type-A receptors are essential for fast inhibitory neurotransmission and are critical in brain function. Surprisingly, expression of receptor subunits is highly variable among individuals, but the cause and impact of this fluctuation remains unknown. We have studied sources of variation for all 19 receptor subunits using massive expression data sets collected across multiple brain regions and platforms in mice and humans. Expression of Gabra1, Gabra2, Gabrb2, Gabrb3, and Gabrg2 is highly variable and heritable among the large cohort of BXD strains derived from crosses of fully sequenced parents—C57BL/6J and DBA/2J. Genetic control of these subunits is complex and highly dependent on tissue and mRNA region. Remarkably, this high variation is generally not linked to phenotypic differences. The single exception is Gabrb3, a locus that is linked to anxiety. We identified upstream genetic loci that influence subunit expression, including three unlinked regions of chromosome 5 that modulate the expression of nine subunits in hippocampus, and that are also associated with multiple phenotypes. Candidate genes within these loci include, Naaa, Nos1, and Zkscan1. We confirmed a high level of coexpression for subunits comprising the major channel—Gabra1, Gabrb2, and Gabrg2—and identified conserved members of this expression network in mice and humans. Gucy1a3, Gucy1b3, and Lis1 are novel and conserved associates of multiple subunits that are involved in inhibitory signaling. Finally, proximal and distal regions of the 3′ UTRs of single subunits have remarkably independent expression patterns in both species. However, corresponding regions of different subunits often show congruent genetic control and coexpression (proximal-to-proximal or distal-to-distal), even in the absence of sequence homology. Our findings identify novel sources of variation that modulate subunit expression and highlight the extraordinary capacity of biological networks to buffer 4–100 fold differences in mRNA levels.

Genetic modulation of striatal volume by loci on Chrs 6 and 17 in BXD recombinant inbred mice

Natural variation in the absolute and relative size of different parts of the human brain is substantial, with a range that often exceeds a factor of 2. Much of this variation is generated by the cumulative effects of sets of unknown gene variants that modulate the proliferation, growth and death of neurons and glial cells. Discovering and testing the functions of these genes should contribute significantly to our understanding of differences in brain development, behavior and disease susceptibility. We have exploited a large population of genetically well‐characterized strains of mice (BXD recombinant inbred strains) to map gene variants that influence the volume of the dorsal striatum (caudate–putamen without nucleus accumbens). We used unbiased methods to estimate volumes bilaterally in a sex‐balanced sample taken from the Mouse Brain Library (www.mbl.org). We generated a matched microarray data set to efficiently evaluate candidate genes (www.genenetwork.org). As in humans, volume of the striatum is highly heritable, with greater than twofold differences among strains. We mapped a locus that modulates striatal volume on chromosome (Chr) 6 at 88 ± 5 Mb. We also uncovered an epistatic interaction between loci on Chr 6 and Chr 17 that modulates striatal volume. Using bioinformatic tools and the corresponding expression database, we have identified positional candidates in these quantitative trait locus intervals.

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

BackgroundThe DBA/2J mouse is not only the oldest inbred strain,but also one of the most widely used strains. DBA/2Jexhibits many unique anatomical, physiological, andbehavior traits. In addition, DBA/2J is one parent of thelarge BXD family of recombinant inbred strains [1]. Thegenome of the other parent of this BXD family—C57BL/6J—has been sequenced and serves as themouse reference genome [2]. We sequenced the gen-ome of DBA/2J using SOLiD and Illumina highthroughput short read protocols to generate a compre-hensive set of ~5 million sequence variants segregatingin the BXD family that ultimately cause developmental,anatomical, functional and behavioral differences amongthese 80+ strains.ResultsWe generated approximately 13.2 and 38.9× whole-gen-ome short reads of DBA/2J females using Illumina GA2and ABI SOLiD massively parallel DNA sequencingplatforms. Comparing to the C57BL/6J reference gen-ome sequence, we identified over 4.5 million singlenucleotide polymorphisms (SNPs), including 84 non-sense and ~11,000 missense mutations, 78% of whichare novel. We also detected ~568,000 insertions anddeletions (indels) within single short reads and ~9,400between mate-paired reads. Approximately 300 inver-sions were detected by SOLiD mate-pair reads, 46 ofwhich span at least one exon. In addition, we identified~22,000 copy number variants (CNVs) in the range of1 Kb to 100 Kb (Figure 1).ConclusionOur study generates the first consensus sequence for theDBA/2J and creates a compendium of sequence andstructural variations that will be used by the communityof researchers who study complex traits in mouse mod-els.Thesequencedataprovideanovelresourcewithwhich to initiate reverse genetic analysis of complex

Sex-specific modulation of gene expression networks in murine hypothalamus.

The hypothalamus contains nuclei and cell populations that are critical in reproduction and that differ significantly between the sexes in structure and function. To examine the molecular and genetic basis for these differences, we quantified gene expression in the hypothalamus of 39 pairs of adult male and female mice belonging to the BXD strains. This experimental design enabled us to define hypothalamic gene coexpression networks and provided robust estimates of absolute expression differences. As expected, sex has the strongest effect on the expression of genes on the X and Y chromosomes (e.g., Uty, Xist, Kdm6a). Transcripts associated with the endocrine system and neuropeptide signaling also differ significantly. Sex-differentiated transcripts often have well delimited expression within specific hypothalamic nuclei that have roles in reproduction. For instance, the estrogen receptor (Esr1) and neurokinin B (Tac2) genes have intense expression in the medial preoptic and arcuate nuclei and comparatively high expression in females. Despite the strong effect of sex on single transcripts, the global pattern of covariance among transcripts is well preserved, and consequently, males and females have well matched coexpression modules. However, there are sex-specific hub genes in functionally equivalent modules. For example, only in males is the Y-linked gene, Uty, a highly connected transcript in a network that regulates chromatin modification and gene transcription. In females, the X chromosome paralog, Kdm6a, takes the place of Uty in the same network. We also find significant effect of sex on genetic regulation and the same network in males and females can be associated with markedly different regulatory loci. With the exception of a few sex-specific modules, our analysis reveals a system in which sets of functionally related transcripts are organized into stable sex-independent networks that are controlled at a higher level by sex-specific modulators.