Gene expression profiling in the striatum of inbred mouse strains with distinct opioid-related phenotypes

Gabra2 is a genetic modifier of Dravet syndrome in mice

Pathogenic variants in epilepsy genes result in a spectrum of clinical severity. One source of phenotypic heterogeneity is modifier genes that affect expressivity of a primary pathogenic variant. Mouse epilepsy models also display varying degrees of clinical severity on different genetic backgrounds. Mice with heterozygous deletion of Scn1a (Scn1a^+/−) model Dravet syndrome, a severe epilepsy most often caused by SCN1A haploinsufficiency. Scn1a^+/− mice recapitulate features of Dravet syndrome, including spontaneous seizures, sudden death, and cognitive/behavioral deficits. Scn1a^+/− mice maintained on the 129S6/SvEvTac (129) strain have normal lifespan and no spontaneous seizures. In contrast, admixture with C57BL/6J (B6) results in epilepsy and premature lethality. We previously mapped Dravet Survival Modifier loci (Dsm1-Dsm5) responsible for strain-dependent differences in survival. Gabra2, encoding the GABA_A α2 subunit, was nominated as a candidate modifier at Dsm1. Direct measurement of GABA_A receptors found lower abundance of α2-containing receptors in hippocampal synapses of B6 mice relative to 129. We also identified a B6-specific single nucleotide deletion within Gabra2 that lowers mRNA and protein by nearly 50%. Repair of this deletion reestablished normal levels of Gabra2 expression. In this study, we used B6 mice with a repaired Gabra2 allele to evaluate Gabra2 as a genetic modifier of severity in Scn1a^+/− mice. Gabra2 repair restored transcript and protein expression, increased abundance of α2-containing GABA_A receptors in hippocampal synapses, and rescued epilepsy phenotypes of Scn1a^+/− mice. These findings validate Gabra2 as a genetic modifier of Dravet syndrome, and support enhancing function of α₂-containing GABA_A receptors as treatment strategy for Dravet syndrome.

Potentiating α2 subunit containing perisomatic GABAA receptors protects against seizures in a mouse model of Dravet syndrome

KEY POINTS

Dravet syndrome mice (Scn1a^+/- ) demonstrate a marked strain dependence for the severity of seizures which is correlated with GABA_A receptor α₂ subunit expression. The α₂ /α₃ subunit selective positive allosteric modulator (PAM) AZD7325 potentiates inhibitory postsynaptic currents (IPSCs) specifically in perisomatic synapses. AZD7325 demonstrates stronger effects on IPSCs in the seizure resistant mouse strain, consistent with higher α₂ subunit expression. AZD7325 demonstrates seizure protective effects in Scn1a^+/- mice without apparent sedative effects in vivo.

ABSTRACT

GABA_A receptor potentiators are commonly used for the treatment of epilepsy, but it is not clear whether targeting distinct GABA_A receptor subtypes will have disproportionate benefits over adverse effects. Here we demonstrate that the α₂ /α₃ selective positive allosteric modulator (PAM) AZD7325 preferentially potentiates hippocampal inhibitory responses at synapses proximal to the soma of CA1 neurons. The effect of AZD7325 on synaptic responses was more prominent in mice on the 129S6/SvEvTac background strain, which have been demonstrated to be seizure resistant in the model of Dravet syndrome (Scn1a^+/- ), and in which the α₂ GABA_A receptor subunits are expressed at higher levels relative to in the seizure prone C57BL/6J background strain. Consistent with this, treatment of Scn1a^+/- mice with AZD7325 elevated the temperature threshold for hyperthermia-induced seizures without apparent sedative effects. Our results in a model system indicate that selectively targeting α₂ is a potential therapeutic option for Dravet syndrome.

Chronic Oxycodone Self-administration Altered Reward-related Genes in the Ventral and Dorsal Striatum of C57BL/6J Mice: An RNA-seq Analysis

Prescription opioid abuse, for example of oxycodone, is a pressing public health issue. This study focuses on how chronic oxycodone self-administration (SA) affects the reward pathways in the mouse brain. In this study, we tested the hypothesis that the expression of reward-related genes in the ventral and dorsal striatum, areas involved in different aspects of opioid addiction models, was altered within 1 h after chronic oxycodone SA, using transcriptome-wide sequencing (RNA-seq). Based on results from earlier human genetic and rodent preclinical studies, we focused on a set of genes that may be associated with the development of addictive diseases and the rewarding effect of drugs of abuse, primarily in the opioid, stress response and classical neurotransmitter systems. We found that 32 transcripts in the ventral striatum, and 7 in the dorsal striatum, were altered significantly in adult mice that had self-administered oxycodone (n = 5) for 14 consecutive days (4 h/day) compared with yoked saline controls (n = 5). The following 5 genes in the ventral striatum showed experiment-wise significant changes: proopiomelanocortin (Pomc) and serotonin 5-HT-2A receptor (Htr2a) were upregulated; serotonin receptor 7 (Htr7), galanin receptor1 (Galr1) and glycine receptor 1 (Glra1) were downregulated. Some genes detected by RNA-seq were confirmed by quantitative polymerase chain reaction (qPCR). Conclusion: A RNA-seq study shows that chronic oxycodone SA alters the expression of several reward-related genes in the dorsal and ventral striatum. These results suggest potential mechanisms underlying neuronal adaptation to chronic oxycodone self-exposure, of relevance to our mechanistic understanding of prescription opioid abuse.

Pharmacological Inhibition of Indoleamine 2,3-Dioxygenase-2 and Kynurenine 3-Monooxygenase, Enzymes of the Kynurenine Pathway, Significantly Diminishes Neuropathic Pain in a Rat Model

Neuropathic pain caused by a primary injury or dysfunction in the peripheral or central nervous system is a tremendous therapeutic challenge. Here, we have collected the first evidence from a single study on the potential contributions to neuropathic pain development by enzymes in the kynurenine pathway [tryptophan 2,3-dioxygenase (TDO), indoleamine 2,3-dioxygenase (IDO1/2), kynurenine 3-monooxygenase (KMO); kynureninase, 3-hydroxyanthranilate-3,4-dioxygenase (HAOO)] at the spinal cord and dorsal root ganglia (DRG) levels. At the spinal cord, mRNA levels of IDO2, KMO, and HAOO were elevated as measured on day 7 after chronic constriction injury in a rat model, parallel to the C1q-positive cell activation. According to our data obtained from primary microglial cell cultures, all enzymes of the kynurenine pathway except TDO were derived from these cells; however, the activation of microglia induced stronger changes in IDO2 and KMO. Our pharmacological studies gave evidence that the repeated intraperitoneal administration of minocycline, a microglia/macrophage inhibitor, not only attenuated tactile and thermal hypersensitivity but also diminished the levels of IDO2 and KMO mRNA. Our further pharmacological studies confirmed that IDO2 and KMO enzymes take part in the development of neuropathic pain, since we observed that the repeated administration of IDO2 (1-methyl-D-tryptophan) and KMO [UPF 648 – (1S,2S)-2-(3,4-dichlorobenzoyl)cyclopropanecarboxylic acid] inhibitors diminished hypersensitivity development as measured on days 2 and 7. The results of our studies show that the kynurenine pathway is an important mediator of neuropathic pain pathology in rats and indicate that IDO2 and KMO represent novel pharmacological targets for treating neuropathy.

Differential regulation of mouse and human Mu opioid receptor gene depends on the single stranded DNA structure of its promoter and α-complex protein 1

The Mu opioid receptor (MOR) mediates various functions of opioid-induced analgesia, euphoria and respiratory depression, and is a major target of opioid analgesics. Understanding of MOR gene expression among species is important for understanding its analgesic function in humans. In the current study, the polypyrimidine/polypurine (PPy/u) region, a key element of MOR gene expression, was compared in humans and mice. The mouse PPy/u element is highly homologous to its human element (84%), and the mouse MOR (mMOR) reporter drives luciferase activity 35-fold more effectively than the human MOR (hMOR) reporter. The structural study of reporter plasmids using S1 nuclease indicates that the mouse PPy/u element has a particular conformational structure, namely a single-stranded DNA (ssDNA) region that promotes strong promoter activity. DNA electrophoretic mobility shift assays demonstrated that the recombinant α-complex protein 1 (α-CP1) is capable of binding to a single-stranded mouse PPy/u sequence. Furthermore, plasmid-expressing α-CP1 activated the expression of a luciferase reporter when cotransfected with a single-stranded (p336/306) construct. In addition, the α-CP1 gene induced the mMOR gene in mouse neuronal cells and did not induce the human neuronal MOR gene. The current study demonstrates that α-CP1 functions as a transcriptional activator in the mMOR gene, but does not function in the hMOR gene due to species-specific structural differences. The differences in human and mouse MOR gene expression are based on α-CP1 and the ssDNA structure of the MOR promoter. The MOR gene is species-specifically regulated, as the PPy/u element adopts a unique species-specific conformation and α-CP1 recruitment.

Analysis of morphine responses in mice reveals a QTL on Chromosome 7

In this study we identified a quantitative trait locus (QTL) on mouse Chromosome 7 associated with locomotor activity and rearing post morphine treatment. This QTL was revealed after correcting for the effects of another QTL peak on Chromosome 10 using composite interval mapping. The positional candidate genes are Syt9 and Ppfibp2. Several other genes within the interval are linked to neural processes, locomotor activity, and the defensive response to harmful stimuli.

The abundance of cis-acting loci leading to differential allele expression in F1 mice and their relationship to loci harboring genes affecting complex traits

Background

Genome-wide surveys have detected cis-acting quantitative trait loci altering levels of RNA transcripts (RNA-eQTLs) by associating SNV alleles to transcript levels. However, the sensitivity and specificity of detection of cis- expression quantitative trait loci (eQTLs) by genetic approaches, reliant as it is on measurements of transcript levels in recombinant inbred strains or offspring from arranged crosses, is unknown, as is their relationship to QTL’s for complex phenotypes.

Results

We used transcriptome-wide differential allele expression (DAE) to detect cis-eQTLs in forebrain and kidney from reciprocal crosses between three mouse inbred strains, 129S1/SvlmJ, DBA/2J, and CAST/EiJ and C57BL/6 J. Two of these crosses were previously characterized for cis-eQTLs and QTLs for various complex phenotypes by genetic analysis of recombinant inbred (RI) strains. 5.4 %, 1.9 % and 1.5 % of genes assayed in forebrain of B6/129SF1, B6/DBAF1, and B6/CASTF1 mice, respectively, showed differential allelic expression, indicative of cis-acting alleles at these genes. Moreover, the majority of DAE QTLs were observed to be tissue-specific with only a small fraction showing cis-effects in both tissues. Comparing DAE QTLs in F1 mice to cis-eQTLs previously mapped in RI strains we observed that many of the cis-eQTLs were not confirmed by DAE. Additionally several novel DAE-QTLs not identified as cis-eQTLs were identified suggesting that there are differences in sensitivity and specificity for QTL detection between the two methodologies. Strain specific DAE QTLs in B6/DBAF1 mice were located in excess at candidate genes for alcohol use disorders, seizures, and angiogenesis previously implicated by genetic linkage in C57BL/6J × DBA/2JF2 mice or BXD RI strains.

Conclusions

Via a survey for differential allele expression in F1 mice, a substantial proportion of genes were found to have alleles altering expression in cis-acting fashion. Comparing forebrain and kidney, many or most of these alleles were tissue-specific in action. The identification of strain specific DAE QTLs, can assist in assessment of candidate genes located within the large intervals associated with trait QTLs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2922-9) contains supplementary material, which is available to authorized users.

Neural basis of reward anticipation and its genetic determinants

Dysfunctional reward processing is implicated in various mental disorders, including attention deficit hyperactivity disorder (ADHD) and addictions. Such impairments might involve different components of the reward process, including brain activity during reward anticipation. We examined brain nodes engaged by reward anticipation in 1,544 adolescents and identified a network containing a core striatal node and cortical nodes facilitating outcome prediction and response preparation. Distinct nodes and functional connections were preferentially associated with either adolescent hyperactivity or alcohol consumption, thus conveying specificity of reward processing to clinically relevant behavior. We observed associations between the striatal node, hyperactivity, and the vacuolar protein sorting-associated protein 4A (VPS4A) gene in humans, and the causal role of Vps4 for hyperactivity was validated in Drosophila Our data provide a neurobehavioral model explaining the heterogeneity of reward-related behaviors and generate a hypothesis accounting for their enduring nature.

Getting the most out of RNA-seq data analysis

Background. A common research goal in transcriptome projects is to find genes that are differentially expressed in different phenotype classes. Biologists might wish to validate such gene candidates experimentally, or use them for downstream systems biology analysis. Producing a coherent differential gene expression analysis from RNA-seq count data requires an understanding of how numerous sources of variation such as the replicate size, the hypothesized biological effect size, and the specific method for making differential expression calls interact. We believe an explicit demonstration of such interactions in real RNA-seq data sets is of practical interest to biologists. Results. Using two large public RNA-seq data sets-one representing strong, and another mild, biological effect size-we simulated different replicate size scenarios, and tested the performance of several commonly-used methods for calling differentially expressed genes in each of them. We found that, when biological effect size was mild, RNA-seq experiments should focus on experimental validation of differentially expressed gene candidates. Importantly, at least triplicates must be used, and the differentially expressed genes should be called using methods with high positive predictive value (PPV), such as NOISeq or GFOLD. In contrast, when biological effect size was strong, differentially expressed genes mined from unreplicated experiments using NOISeq, ASC and GFOLD had between 30 to 50% mean PPV, an increase of more than 30-fold compared to the cases of mild biological effect size. Among methods with good PPV performance, having triplicates or more substantially improved mean PPV to over 90% for GFOLD, 60% for DESeq2, 50% for NOISeq, and 30% for edgeR. At a replicate size of six, we found DESeq2 and edgeR to be reasonable methods for calling differentially expressed genes at systems level analysis, as their PPV and sensitivity trade-off were superior to the other methods'. Conclusion. When biological effect size is weak, systems level investigation is not possible using RNAseq data, and no meaningful result can be obtained in unreplicated experiments. Nonetheless, NOISeq or GFOLD may yield limited numbers of gene candidates with good validation potential, when triplicates or more are available. When biological effect size is strong, NOISeq and GFOLD are effective tools for detecting differentially expressed genes in unreplicated RNA-seq experiments for qPCR validation. When triplicates or more are available, GFOLD is a sharp tool for identifying high confidence differentially expressed genes for targeted qPCR validation; for downstream systems level analysis, combined results from DESeq2 and edgeR are useful.

Self administration of oxycodone alters synaptic plasticity gene expression in the hippocampus differentially in male adolescent and adult mice

Abuse and addiction to prescription opioids such as oxycodone (a short-acting Mu opioid receptor (MOP-r) agonist) in adolescence is a pressing public health issue. We have previously shown differences in oxycodone self-administration behaviors between adolescent and adult C57BL/6J mice and expression of striatal neurotransmitter receptor genes, in areas involved in reward. In this study, we aimed to determine whether oxycodone self-administration differentially affects genes regulating synaptic plasticity in the hippocampus of adolescent compared to adult mice, since the hippocampus may be involved in learning aspects associated with chronic drug self administration. Hippocampus was isolated for mRNA analysis from mice that had self administered oxycodone (0.25 mg/kg/infusion) 2h/day for 14 consecutive days or from yoked saline controls. Gene expression was analyzed with real-time polymerase chain reaction (PCR) using a commercially available "synaptic plasticity" PCR array containing 84 genes. We found that adolescent and adult control mice significantly differed in the expression of several genes in the absence of oxycodone exposure, including those coding for mitogen-activated protein kinase, calcium/calmodulin-dependent protein kinase II gamma subunit, glutamate receptor, ionotropic AMPA2 and metabotropic 5. Chronic oxycodone self administration increased proviral integration site 1 (Pim1) and thymoma viral proto-oncogene 1 mRNA levels compared to controls in both age groups. Both Pim1 and cadherin 2 mRNAs showed a significant combined effect of Drug Condition and Age × Drug Condition. Furthermore, the mRNA levels of both cadherin 2 and cAMP response element modulators showed an experiment-wise significant difference between oxycodone and saline control in adult but not in adolescent mice. Overall, this study demonstrates for the first time that chronic oxycodone self-administration differentially alters synaptic plasticity gene expression in the hippocampus of adolescent and adult mice.

Whole-genome expression profile in zebrafish embryos after chronic exposure to morphine: identification of new genes associated with neuronal function and mu opioid receptor expression

BACKGROUND

A great number of studies have investigated changes induced by morphine exposure in gene expression using several experimental models. In this study, we examined gene expression changes during chronic exposure to morphine during maturation and differentiation of zebrafish CNS.

RESULTS

Microarray analysis showed 254 genes whose expression was identified as different by at least 1.3 fold change following chronic morphine exposure as compared to controls. Of these, several novel genes (grb2, copb2, otpb, magi1b, grik-l, bnip4 and sox19b) have been detected for the first time in an experimental animal model treated with morphine. We have also identified a subset of genes (dao.1, wls, bnip4 and camk1γb) differentially expressed by chronic morphine exposure whose expression is related to mu opioid receptor gene expression. Altered expression of copb2, bnip4, sox19b, otpb, dao.1, grik-l and wls is indicative of modified neuronal development, CNS patterning processes, differentiation and dopaminergic neurotransmission, serotonergic signaling pathway, and glutamatergic neurotransmission. The deregulation of camk1γb signaling genes suggests an activation of axonogenesis and dendritogenesis.

CONCLUSIONS

Our study identified different functional classes of genes and individual candidates involved in the mechanisms underlying susceptibility to morphine actions related to CNS development. These results open new lines to study the treatment of pain and the molecular mechanisms involved in addiction. We also found a set of zebrafish-specific morphine-induced genes, which may be putative targets in human models for addiction and pain processes.

Involvement of pro- and antinociceptive factors in minocycline analgesia in rat neuropathic pain model

In neuropathic pain the repeated minocycline treatment inhibited the mRNA and protein expression of the microglial markers and metalloproteinase-9 (MMP-9). The minocycline diminished the pronociceptive (IL-6, IL-18), but not antinociceptive (IL-1alpha, IL-4, IL-10) cytokines at the spinal cord level. In vitro primary cell culture studies have shown that MMP-9, TIMP-1, IL-1beta, IL-1alpha, IL-6, IL-10, and IL-18 are of microglial origin. Minocycline reduces the production of pronociceptive factors, resulting in a more potent antinociceptive effect. This change in the ratio between pro- and antinociceptive factors, in favour of the latter may be the mechanism of minocycline analgesia in neuropathy.

Impact of P-glycoprotein at the blood-brain barrier on the uptake of heroin and its main metabolites: behavioral effects and consequences on the transcriptional responses and reinforcing properties

RATIONALE

Transport across the BBB is a determinant of the rate and extent of drug distribution in the brain. Heroin exerts its effects through its principal metabolites 6-monoacetyl-morphine (6-MAM) and morphine. Morphine is a known substrate of P-glycoprotein (P-gp) at the blood-brain-barrier (BBB) however, little is known about the interaction of heroin and 6-MAM with P-gp.

OBJECTIVE

The objective of this paper is to study the role of the P-gp-mediated efflux at the BBB in the behavioral and molecular effects of heroin and morphine.

METHODS

The transport rates of heroin and its main metabolites, at the BBB, were measured in mice by in situ brain perfusion. We then examined the effect of inhibition of P-gp on the acute nociception, locomotor activity, and gene expression modulations induced by heroin and morphine. The effect of P-gp inhibition during the acquisition of morphine-induced place preference was also studied.

RESULTS

Inhibition of P-gp significantly increased the uptake of morphine but not that of heroin nor 6-MAM. Inhibition of P-gp significantly increased morphine-induced acute analgesia and locomotor activity but did not affect the behavioral effects of heroin; in addition, acute transcriptional responses to morphine were selectively modulated in the nucleus accumbens. Increasing morphine uptake by the brain significantly increased its reinforcing properties in the place preference paradigm.

CONCLUSIONS

The present study demonstrated that acute inhibition of P-gp not only modulates morphine-induced behavioral effects but also its transcriptional effects and reinforcing properties. This suggests that, in the case of morphine, transport across the BBB is critical for the development of dependence.

Astrocytes are a neural target of morphine action via glucocorticoid receptor-dependent signaling

Chronic opioid use leads to the structural reorganization of neuronal networks, involving genetic reprogramming in neurons and glial cells. Our previous in vivo studies have revealed that a significant fraction of the morphine-induced alterations to the striatal transcriptome included glucocorticoid (GC) receptor (GR)-dependent genes. Additional analyses suggested glial cells to be the locus of these changes. In the current study, we aimed to differentiate the direct transcriptional effects of morphine and a GR agonist on primary striatal neurons and astrocytes. Whole-genome transcriptional profiling revealed that while morphine had no significant effect on gene expression in both cell types, dexamethasone significantly altered the transcriptional profile in astrocytes but not neurons. We obtained a complete dataset of genes undergoing the regulation, which includes genes related to glucose metabolism (Pdk4), circadian activity (Per1) and cell differentiation (Sox2). There was also an overlap between morphine-induced transcripts in striatum and GR-dependent transcripts in cultured astrocytes. We further analyzed the regulation of expression of one gene belonging to both groups, serum and GC regulated kinase 1 (Sgk1). We identified two transcriptional variants of Sgk1 that displayed selective GR-dependent upregulation in cultured astrocytes but not neurons. Moreover, these variants were the only two that were found to be upregulated in vivo by morphine in a GR-dependent fashion. Our data suggest that the morphine-induced, GR-dependent component of transcriptome alterations in the striatum is confined to astrocytes. Identification of this mechanism opens new directions for research on the role of astrocytes in the central effects of opioids.

Potential translational targets revealed by linking mouse grooming behavioral phenotypes to gene expression using public databases

Rodent self-grooming is an important, evolutionarily conserved behavior, highly sensitive to pharmacological and genetic manipulations. Mice with aberrant grooming phenotypes are currently used to model various human disorders. Therefore, it is critical to understand the biology of grooming behavior, and to assess its translational validity to humans. The present in-silico study used publicly available gene expression and behavioral data obtained from several inbred mouse strains in the open-field, light-dark box, elevated plus- and elevated zero-maze tests. As grooming duration differed between strains, our analysis revealed several candidate genes with significant correlations between gene expression in the brain and grooming duration. The Allen Brain Atlas, STRING, GoMiner and Mouse Genome Informatics databases were used to functionally map and analyze these candidate mouse genes against their human orthologs, assessing the strain ranking of their expression and the regional distribution of expression in the mouse brain. This allowed us to identify an interconnected network of candidate genes (which have expression levels that correlate with grooming behavior), display altered patterns of expression in key brain areas related to grooming, and underlie important functions in the brain. Collectively, our results demonstrate the utility of large-scale, high-throughput data-mining and in-silico modeling for linking genomic and behavioral data, as well as their potential to identify novel neural targets for complex neurobehavioral phenotypes, including grooming.

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 (h(2)≥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 F(1) 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.

Mouse strain specific gene expression differences for illumina microarray expression profiling in embryos

Background

In the field of mouse genetics the advent of technologies like microarray based expression profiling dramatically increased data availability and sensitivity, yet these advanced methods are often vulnerable to the unavoidable heterogeneity of in vivo material and might therefore reflect differentially expressed genes between mouse strains of no relevance to a targeted experiment. The aim of this study was not to elaborate on the usefulness of microarray analysis in general, but to expand our knowledge regarding this potential “background noise” for the widely used Illumina microarray platform surpassing existing data which focused primarily on the adult sensory and nervous system, by analyzing patterns of gene expression at different embryonic stages using wild type strains and modern transgenic models of often non-isogenic backgrounds.

Results

Wild type embryos of 11 mouse strains commonly used in transgenic and molecular genetic studies at three developmental time points were subjected to Illumina microarray expression profiling in a strain-by-strain comparison. Our data robustly reflects known gene expression patterns during mid-gestation development. Decreasing diversity of the input tissue and/or increasing strain diversity raised the sensitivity of the array towards the genetic background. Consistent strain sensitivity of some probes was attributed to genetic polymorphisms or probe design related artifacts.

Conclusion

Our study provides an extensive reference list of gene expression profiling background noise of value to anyone in the field of developmental biology and transgenic research performing microarray expression profiling with the widely used Illumina microarray platform. Probes identified as strain specific background noise further allow for microarray expression profiling on its own to be a valuable tool for establishing genealogies of mouse inbred strains.

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.

Neuronal aggregates are associated with phenotypic onset in the R6/2 Huntington's disease transgenic mouse

Huntington's disease (HD) is caused by the expansion of the polyglutamine tract expressed in the huntingtin protein. Data from patients show a strong negative correlation between CAG repeat size and age of disease onset. Recent studies in mixed background C57×CBA R6/2 mice suggest the inverse correlation observed in the human disease may not be replicated in some animal models of HD. To further clarify the relationship between repeat length and age of onset, congenic C57BL6/J R6/2 transgenic mice expressing 110, 260 or 310 CAG were tested in a comprehensive behavioral battery at multiple ages. Data confirmed the findings of earlier studies and indicate that on a pure C57BL6/J genetic background, R6/2 mice with larger repeats exhibit a delay in phenotypic onset with increasing polyglutamine size (6 weeks in 110 CAG and 17 weeks in 310 CAG mice). Further analysis confirmed a decrease in transgene transcript expression in 310 CAG mice as well as differential aggregated protein localization in association with repeat length. Mice expressing 110 CAG developed aggregates that localized almost exclusively to the nucleus of neuronal cells in the striatum and cortex. In contrast, tissue from 310 CAG mice exhibited predominantly extranuclear inclusions. Novel mutant protein analysis obtained using time-resolved fluorescence resonance energy transfer (FRET) revealed that soluble protein levels decreased with disease onset in R6/2 mice while aggregated protein levels increased. We believe that these data suggest a role for aggregation and inclusion localization in HD pathogenesis and propose a mechanism for the age of onset delay observed in R6/2 mice.

Using genome-wide expression profiling to define gene networks relevant to the study of complex traits: from RNA integrity to network topology

Postgenomic studies of the function of genes and their role in disease have now become an area of intense study since efforts to define the raw sequence material of the genome have largely been completed. The use of whole-genome approaches such as microarray expression profiling and, more recently, RNA-sequence analysis of transcript abundance has allowed an unprecedented look at the workings of the genome. However, the accurate derivation of such high-throughput data and their analysis in terms of biological function has been critical to truly leveraging the postgenomic revolution. This chapter will describe an approach that focuses on the use of gene networks to both organize and interpret genomic expression data. Such networks, derived from statistical analysis of large genomic datasets and the application of multiple bioinformatics data resources, potentially allow the identification of key control elements for networks associated with human disease, and thus may lead to derivation of novel therapeutic approaches. However, as discussed in this chapter, the leveraging of such networks cannot occur without a thorough understanding of the technical and statistical factors influencing the derivation of genomic expression data. Thus, while the catch phrase may be "it's the network … stupid," the understanding of factors extending from RNA isolation to genomic profiling technique, multivariate statistics, and bioinformatics are all critical to defining fully useful gene networks for study of complex biology.

5-hydroxytryptamine type 3 receptor modulates opioid-induced hyperalgesia and tolerance in mice

BACKGROUND

Opioid-induced hyperalgesia (OIH) and tolerance are challenging maladaptations associated with opioids in managing pain. Recent genetic studies and the existing literature suggest the 5-hydroxytryptamine type 3 (5-HT3) receptor participates in these phenomena. The location of the relevant receptor populations and the interactions between the 5-HT3 system and other systems controlling OIH and tolerance have not been explored, however. We hypothesized that 5-HT3 receptors modulate OIH and tolerance, and that this modulation involves the control of expression of multiple neurotransmitter and receptor systems.

METHODS

C57BL/6 mice were exposed to a standardized 4-day morphine administration protocol. The 5-HT3 antagonist ondansetron was administered either during or after the conclusion of morphine administration. Mechanical testing was used to quantify OIH, and thermal tail-flick responses were used to measure morphine tolerance. In other experiments spinal cord and dorsal root ganglion tissues were harvested for analysis of messenger RNA concentrations by real-time polymerase chain reaction or immunochemistry analysis.

RESULTS

The results showed that (1) systemic or intrathecal injection of ondansetron significantly prevented and reversed OIH, but not local intraplantar injection; (2) systemic or intrathecal injection of ondansetron prevented and reversed tolerance; and (3) ondansetron blocked morphine-induced increases of multiple genes relevant to OIH and tolerance in dorsal root ganglion and spinal cord.

CONCLUSIONS

Morphine acts via a 5-HT3-dependent mechanism to support multiple maladaptations to the chronic administration of morphine. Furthermore, the use of 5-HT3 receptor antagonists may provide a new avenue to prevent or reverse OIH and tolerance associated with chronic opioid use.

Mouse mutants for the nicotinic acetylcholine receptor ß2 subunit display changes in cell adhesion and neurodegeneration response genes

Mice lacking expression of the ß2 subunit of the neuronal nicotinic acetylcholine receptor (CHRNB2) display abnormal retinal waves and a dispersed projection of retinal ganglion cell (RGC) axons to their dorsal lateral geniculate nuclei (dLGNs). Transcriptomes of LGN tissue from two independently generated Chrnb2−/− mutants and from wildtype mice were obtained at postnatal day 4 (P4), during the normal period of segregation of eye-specific afferents to the LGN. Microarray analysis reveals reduced expression of genes located on the cell membrane or in extracellular space, and of genes active in cell adhesion and calcium signaling. In particular, mRNA for cadherin 1 (Cdh1), a known axon growth regulator, is reduced to nearly undetectable levels in the LGN of P4 mutant mice and Lypd2 mRNA is similarly suppressed. Similar analysis of retinal tissue shows increased expression of crumbs 1 (Crb1) and chemokine (C-C motif) ligand 21 (Ccl21) mRNAs in Chrnb2−/− mutant animals. Mutations in these genes are associated with retinal neuronal degeneration. The retinas of Chrnb2−/− mutants are normal in appearance, but the increased expression of these genes may also be involved in the abnormal projection patterns of RGC to the LGN. These data may provide the tools to distinguish the interplay between neural activity and molecular expression. Finally, comparison of the transcriptomes of the two different Chrnb2−/− mutant strains reveals the effects of genetic background upon gene expression.

A B2 SINE insertion in the Comt1 gene (Comt1(B2i)) results in an overexpressing, behavior modifying allele present in classical inbred mouse strains

Catechol-O-methyltransferase (COMT) is a key enzyme for dopamine catabolism and COMT is a candidate gene for human psychiatric disorders. In mouse it is located on chromosome 16 in a large genomic region of extremely low variation among the classical inbred strains, with no confirmed single nucleotide polymorphisms (SNPs) between strains C57BL/6J and DBA/2J within a 600-kB window. We found a B2 SINE in the 3′ untranslated region (UTR) of Comt1 which is present in C57BL/6J (Comt1^B2i) and other strains including 129 (multiple sublines), but is not found in DBA/2J (Comt1⁺) and many other strains including wild-derived Mus domesticus, M. musculus, M. molossinus, M.castaneus and M. spretus. Comt1^B2i is absent in strains closely related to C57BL/6, such as C57L and C57BR, indicating that it was polymorphic in the cross that gave rise to these strains. The strain distribution of Comt1^B2i indicates a likely origin of the allele in the parental Lathrop stock. A stringent association test, using 670 highly outbred mice (Boulder Heterogeneous Stock), indicates that this insertion allele may be responsible for a difference in behavior related to exploration. Gene expression differences at the mRNA and enzyme activity level (1.7-fold relative to wild type) indicate a mechanism for this behavioral effect. Taken together, these findings show that Comt1^B2i (a B2 SINE insertion) results in a relatively modest difference in Comt1 expression and enzyme activity (comparable to the human Val-Met polymorphism) which has a demonstrable behavioral phenotype across a variety of outbred genetic backgrounds.

Differential gene expression activity among species-specific polypyrimidine/polypurine motifs in mu opioid receptor gene promoters

The mu opioid receptor (MOR) is the principle molecular target of opioid analgesics. An appropriate understanding of MOR gene expression across species is critical for understanding its analgesic functions in humans. Here, we undertake a cross-species analysis of the polymorphic polypyrimidine/polypurine (PPy/u) motif, a key enhancer of MOR gene expression. The mouse PPy/u motif is highly homologous to those of rat (67%) and human (83%), but drives reporter gene expression tenfold and fivefold more effectively than those of rat and human, respectively. Circular dichroism profiles of PPy/u oligonucleotides from different species showed that they are primarily different in structure. Conformational studies of reporter plasmids using confocal Raman spectra, S1 nuclease and restriction enzymes demonstrated that the structural difference is the result of changes in the phosphodiester backbone. Furthermore, these conformational disparities produce differences in torsional stress, as shown by topoisomerase II relaxation and activation of different levels of gene expression under hypertonic conditions. This study demonstrates that homologous PPy/u motifs adopt unique species-specific conformations with different mechanisms and activities for gene expression. We further discuss how structural aspects of transcription regulatory elements, rather than the sequence itself, are significant when studying functional gene expression regulatory elements.

Global gene expression profile progression in Gaucher disease mouse models

BACKGROUND

Gaucher disease is caused by defective glucocerebrosidase activity and the consequent accumulation of glucosylceramide. The pathogenic pathways resulting from lipid laden macrophages (Gaucher cells) in visceral organs and their abnormal functions are obscure.

RESULTS

To elucidate this pathogenic pathway, developmental global gene expression analyses were conducted in distinct Gba1 point-mutated mice (V394L/V394L and D409 V/null). About 0.9 to 3% of genes had altered expression patterns (≥ ± 1.8 fold change), representing several categories, but particularly macrophage activation and immune response genes. Time course analyses (12 to 28 wk) of INFγ-regulated pro-inflammatory (13) and IL-4-regulated anti-inflammatory (11) cytokine/mediator networks showed tissue differential profiles in the lung and liver of the Gba1 mutant mice, implying that the lipid-storage macrophages were not functionally inert. The time course alterations of the INFγ and IL-4 pathways were similar, but varied in degree in these tissues and with the Gba1 mutation.

CONCLUSIONS

Biochemical and pathological analyses demonstrated direct relationships between the degree of tissue glucosylceramides and the gene expression profile alterations. These analyses implicate IFNγ-regulated pro-inflammatory and IL-4-regulated anti-inflammatory networks in differential disease progression with implications for understanding the Gaucher disease course and pathophysiology.

Identification of cis-regulatory elements in the mammalian genome: the cREMaG database

Background

A growing number of gene expression-profiling datasets provides a reliable source of information about gene co-expression. In silico analyses of the properties shared among the promoters of co-expressed genes facilitates the identification of transcription factors (TFs) involved in the co-regulation of those genes. Our previous experience with microarray data led to the development of a database suitable for the examination of regulatory motifs in the promoters of co-expressed genes.

Methodology

We introduce the cREMaG (cis-Regulatory Elements in the Mammalian Genome) system designed for in silico studies of the promoter properties of co-regulated mammalian genes. The cREMaG system offers an analysis of data obtained from human, mouse, rat, bovine and canine gene expression-profiling studies. More than eight analysis parameters can be utilized in user-defined combinations. The selection of alternative transcription start sites and information about CpG islands are also available.

Conclusions

Using the cREMaG system, we successfully identified TFs mediating transcriptional responses in reference gene sets. The cREMaG system facilitates in silico studies of mammalian transcriptional gene regulation. The resource is freely available at http://www.cremag.org.

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 ( approximately 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.

Assessing individual differences in genome-wide gene expression in human whole blood: reliability over four hours and stability over 10 months

Studying the causes and correlates of natural variation in gene expression in healthy populations assumes that individual differences in gene expression can be reliably and stably assessed across time. However, this is yet to be established. We examined 4-hour test-retest reliability and 10 month test-retest stability of individual differences in gene expression in ten 12-year-old children. Blood was collected on four occasions: 10 a.m. and 2 p.m. on Day 1 and 10 months later at 10 a.m. and 2 p.m. Total RNA was hybridized to Affymetrix-U133 plus 2.0 arrays. For each probeset, the correlation across individuals between 10 a.m. and 2 p.m. on Day 1 estimates test-retest reliability. We identified 3,414 variable and abundantly expressed probesets whose 4-hour test-retest reliability exceeded .70, a conventionally accepted level of reliability, which we had 80% power to detect. Of the 3,414 reliable probesets, 1,752 were also significantly reliable 10 months later. We assessed the long-term stability of individual differences in gene expression by correlating the average expression level for each probe-set across the two 4-hour assessments on Day 1 with the average level of each probe-set across the two 4-hour assessments 10 months later. 1,291 (73.7%) of the 1,752 probe-sets that reliably detected individual differences across 4 hours on two occasions, 10 months apart, also stably detected individual differences across 10 months. Heritability, as estimated from the MZ twin intraclass correlations, is twice as high for the 1,752 reliable probesets versus all present probesets on the array (0.68 vs 0.34), and is even higher (0.76) for the 1,291 reliable probesets that are also stable across 10 months. The 1,291 probesets that reliably detect individual differences from a single peripheral blood collection and stably detect individual differences over 10 months are promising targets for research on the causes (e.g., eQTLs) and correlates (e.g., psychopathology) of individual differences in gene expression.

Genetic variance contributes to dopamine and opioid receptor antagonist-induced inhibition of intralipid (fat) intake in inbred and outbred mouse strains

Preference for and intake of solid and emulsified fat (intralipid) solutions vary across different mouse strains. Fat intake in rodents is inhibited by dopamine and opioid receptor antagonists, but any variation in these responses as a function of genetic background is unknown. Therefore, the present study compared the ability of dopamine D1-like (SCH23390) and general opioid (naltrexone) receptor antagonism to alter intake of fat emulsions (intralipid) in mice. Two-hour intakes of 5% intralipid were measured (5-120 min) in seven inbred (BALB/c, C57BL/6, C57BL/10, DBA/2, SJL, SWR, 129P3) and one outbred (CD-1) mouse strains following treatment with vehicle, SCH23390 (50-1600 nmol/kg, ip) and naltrexone (0.001-5 mg/kg, sc). SCH23390 significantly, dose-dependently and differentially reduced intralipid intake at all five (DBA/2, SWR, CD-1), four (SJL, C57BL/6), three (129P3) and one (C57BL/10) of the doses tested, but failed to affect intralipid intake in BALB/c mice. Naltrexone significantly, dose-dependently and differentially reduced intralipid intake at all four (DBA/2), three (SWR, SJL), two (CD-1, C57BL/10) and one (C57BL/6, 129P3) of the doses tested, and also failed to affect intralipid intake in BALB/cJ mice. SCH23390 and naltrexone were respectively 13.3-fold and 9.3-fold more potent in inhibiting intralipid intake in the most sensitive (DBA/2) relative to the least sensitive (BALB/c) mouse strains. A strong positive relationship (r=0.91) was observed for the abilities of SCH23390 and naltrexone to inhibit intralipid intake across strains. These findings indicate that dopaminergic and opioid signaling mechanisms differentially control intralipid intake across different mouse strains, suggesting important genetic and pharmacological interactions in the short-term control of rewarding and post-ingestive consequences of fat intake.

Single-nucleotide polymorphism (A118G) in exon 1 of OPRM1 gene causes alteration in downstream signaling by mu-opioid receptor and may contribute to the genetic risk for addiction

The opioid receptor mu1 (OPRM1) mediates the action of morphine. Although genetic background plays an important role in the susceptibility toward abuse of drugs as evident from familial, adoption and twin studies, association of specific single-nucleotide polymorphisms of OPRM1 gene with narcotic addiction is to be established. Here, we demonstrate the involvement of A118G polymorphism of exon1 of human OPRM1 gene (hOPRM1), with heroin and alcohol addiction, in a population in eastern India. Statistical analysis exhibited a significant association of G allele with both heroin and alcohol addiction with a risk factor of P(trend) < 0.05. The functional significance of G allele in A118G single-nucleotide polymorphisms was evaluated by studying the regulation of protein kinase A (PKA), pCREB, and pERK1/2 by morphine in Neuro 2A cells, stably transfected with either wild type or A118G mutant hOPRM1. Unlike acute morphine treatment, both chronic morphine exposure and withdrawal precipitated by naloxone were differentially regulated by A118 and G118 receptor isoforms when both PKA and pERK1/2 activities were compared. Results suggest that the association of A118G polymorphism to heroin and alcohol addiction may be because of the altered regulation of PKA and pERK1/2 during opioid and alcohol exposures.

Inbred mouse strains differ in multiple hippocampal activity traits

A major challenge in neuroscience is to identify genes that influence specific behaviors and to understand the intermediary neuronal mechanisms. One approach is to identify so-called endophenotypes at different levels of neuronal organization from synapse to brain activity. An endophenotype is a quantitative trait that is closer to the gene action than behavior, and potentially a marker of neuronal mechanisms underlying behavior. Hippocampal activity and, in particular, hippocampal oscillations have been suggested to underlie various cognitive and motor functions. To identify quantitative traits that are potentially useful for identifying genes influencing hippocampal activity, we measured gamma oscillations and spontaneous activity in acute hippocampal slices from eight inbred mouse strains under three experimental conditions. We estimated the heritability of more than 200 quantitative traits derived from this activity. We observed significant differences between the different mouse strains, particularly in the amplitude of the activity and the correlation between activities in different hippocampal subregions. Interestingly, these traits had a low genetic correlation between the three experimental conditions, which suggests that different genetic components influence the activity in different conditions. Our findings show that several traits of hippocampal gamma oscillations and spontaneous activity are heritable and could thus be potentially useful in gene-finding strategies based on endophenotypes.

Identification of candidate genes and gene networks specifically associated with analgesic tolerance to morphine

Chronic morphine administration may alter the expression of hundreds to thousands of genes. However, only a subset of these genes is likely involved in analgesic tolerance. In this report, we used a behavior genetics strategy to identify candidate genes specifically linked to the development of morphine tolerance. Two inbred genotypes [C57BL/6J (B6), DBA2/J (D2)] and two reciprocal congenic genotypes (B6D2, D2B6) with the proximal region of chromosome 10 (Chr10) introgressed into opposing backgrounds served as the behavior genetic filter. Tolerance after therapeutically relevant doses of morphine developed most rapidly in the B6 followed by the B6D2 genotype and did not develop in the D2 mice and only slightly in the D2B6 animals indicating a strong influence of the proximal region of Chr10 in the development of tolerance. Gene expression profiling and pattern matching identified 64, 53, 86, and 123 predisposition genes and 81, 96, 106, and 82 tolerance genes in the periaqueductal gray (PAG), prefrontal cortex, temporal lobe, and ventral striatum, respectively. A potential gene network was identified in the PAG in which 19 of the 34 genes were strongly associated with tolerance. Eleven of the network genes were found to reside in quantitative trait loci previously associated with morphine-related behaviors, whereas seven were predictive of tolerance (morphine-naive condition). Overall, the genes modified by chronic morphine administration show a strong presence in canonical pathways representative of neuroadaptation. A potentially significant role for the micro-RNA and epigenetic mechanisms in response to chronic administration of pharmacologically relevant doses of morphine was highlighted by candidate genes Dicer and H19.

Genetic correlates of morphine withdrawal in 14 inbred mouse strains

Male mice from 14 standard inbred strains were exposed to morphine in a sustained released preparation injected subcutaneously. Five hours later withdrawal was precipitated by intraperitoneal injection of naloxone. Mice were tested from 0 to 15 min after naloxone for withdrawal jumping behavior, and then from minute 15-16 for other signs, including boli count, presence of soft stool, lacrimation, "wet dog" shakes, and air chewing. They were also assessed for change in body temperature 17 min after naloxone. Strains differed markedly in the severity of withdrawal for jumping, change in body temperature, and number of fecal boli. Strains also differed in percentage of animals displaying soft stool and air chewing behavior. The other two signs were seen at too low frequency for analysis. Correlations of strain mean withdrawal severity with other responses to morphine and other abused drugs showed that high morphine withdrawal jumping and low change in body temperature were both genetically related to high morphine consumption, but not generally to other measures of morphine withdrawal or morphine sensitivity.

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 approximately 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.

Deterioration of atherosclerosis in mice lacking angiotensin II type 1A receptor in bone marrow-derived cells

The renin-angiotensin system (RAS) modulates end-organ damages, resulting in cardiovascular and kidney diseases. Experiments both in vitro and in vivo demonstrate that the angiotensin II (Ang II) type 1 (AT1) receptor pathway also exerts pro-inflammatory and pro-atherogenic effects on bone marrow-derived cells (BMDCs). Here, we investigated how AT1 receptor expression by BMDCs contributes to atherosclerosis and kidney injury in vivo by transplanting BM into RAS-activated transgenic mice. There was no difference in the extent of kidney damage between mice receiving BM transplants from mutant mice lacking the angiotensin II type 1a receptor (AT1a) gene and mice receiving transplants from wild-type (WT) mice. However, mice receiving transplants from AT1a 'knockout' (KO) mice displayed accelerated lethality and atherosclerotic lesions. These results indicated that the effects of AT1a receptor on BMDCs are organ dependent. Microarray expression profiling of macrophages from AT1a-KO mice revealed significant changes in the mRNA levels for a number of genes implicated in atherosclerosis. In accordance with the in vivo atherosclerosis results, AT1a-KO macrophages exhibited greater uptake of modified lipoproteins relative to macrophages from WT mice. We propose that the expression of AT1a receptor by BMDCs limits atherosclerosis in vivo.

Morphine effects on striatal transcriptome in mice

Global transcriptional analysis of mouse striata following acute and chronic exposure to morphine reveals multiple physiological factors which may affect opioid-related phenotypes and implicates a number of gene networks, including glucocorticoid receptor regulated genes, in the response to this opioid.

Background

Chronic opiate use produces molecular and cellular adaptations in the nervous system that lead to tolerance, physical dependence, and addiction. Genome-wide comparison of morphine-induced changes in brain transcription of mouse strains with different opioid-related phenotypes provides an opportunity to discover the relationship between gene expression and behavioral response to the drug.

Results

Here, we analyzed the effects of single and repeated morphine administrations in selected inbred mouse strains (129P3/J, DBA/2J, C57BL/6J, and SWR/J). Using microarray-based gene expression profiling in striatum, we found 618 (false discovery rate < 1%) morphine-responsive transcripts. Through ontologic classification, we linked particular sets of genes to biologic functions, including metabolism, transmission of nerve impulse, and cell-cell signaling. We identified numerous novel morphine-regulated genes (for instance, Olig2 and Camk1g), and a number of transcripts with strain-specific changes in expression (for instance, Hspa1a and Fzd2). Moreover, transcriptional activation of a pattern of co-expressed genes (for instance, Tsc22d3 and Nfkbia) was identified as being mediated via the glucocorticoid receptor (GR). Further studies revealed that blockade of the GR altered morphine-induced locomotor activity and development of physical dependence.

Conclusion

Our results indicate that there are differences between strains in the magnitude of transcriptional response to acute morphine treatment and in the degree of tolerance in gene expression observed after chronic morphine treatment. Using whole-genome transcriptional analysis of morphine effects in the striatum, we were able to reveal multiple physiological factors that may influence opioid-related phenotypes and to relate particular gene networks to this complex trait. The results also suggest the possible involvement of GR-regulated genes in mediating behavioral response to morphine.

The PhenoGen informatics website: tools for analyses of complex traits

Background

With the advent of "omics" (e.g. genomics, transcriptomics, proteomics and phenomics), studies can produce enormous amounts of data. Managing this diverse data and integrating with other biological data are major challenges for the bioinformatics community. Comprehensive new tools are needed to store, integrate and analyze the data efficiently.

Description

The PhenoGen Informatics website is a comprehensive toolbox for storing, analyzing and integrating microarray data and related genotype and phenotype data. The site is particularly suited for combining QTL and microarray data to search for "candidate" genes contributing to complex traits. In addition, the site allows, if desired by the investigators, sharing of the data. Investigators can conduct "in-silico" microarray experiments using their own and/or "shared" data.

Conclusion

The PhenoGen website provides access to tools that can be used for high-throughput data storage, analyses and interpretation of the results. Some of the advantages of the architecture of the website are that, in the future, the present set of tools can be adapted for the analyses of any type of high-throughput "omics" data, and that access to new tools, available in the public domain or developed at PhenoGen, can be easily provided.

Transcriptional profiling of C57 and DBA strains of mice in the absence and presence of morphine

Background

The mouse C57BL/6 (C57) and DBA/2J (DBA) inbred strains differ substantially in many aspects of their response to drugs of abuse. The development of microarray analyses represents a genome-wide method for measuring differences across strains, focusing on expression differences. In the current study, we carried out microarray analysis in C57 and DBA mice in the nucleus accumbens of drug-naïve and morphine-treated animals.

Results

We identified mRNAs with altered expression between the two strains. We validated the mRNA expression changes of several such mRNAs, including Gnb1, which has been observed to be regulated by several drugs of abuse. In addition, we validated alterations in the enzyme activity of one mRNA product, catechol-O-methyltransferase (Comt). Data mining of expression and behavioral data indicates that both Gnb1 and Comt expression correlate with aspects of drug response in C57/DBA recombinant inbred strains. Pathway analysis was carried out to identify pathways showing significant alterations as a result of treatment and/or due to strain differences. These analyses identified axon guidance genes, particularly the semaphorins, as showing altered expression in the presence of morphine, and plasticity genes as showing altered expression across strains. Pathway analysis of genes showing strain by treatment interaction suggest that the phosphatidylinositol signaling pathway may represent an important difference between the strains as related to morphine exposure.

Conclusion

mRNAs with differing expression between the two strains could potentially contribute to strain-specific responses to drugs of abuse. One such mRNA is Comt and we hypothesize that altered expression of Comt may represent a potential mechanism for regulating the effect of, and response to, multiple substances of abuse. Similarly, a role for Gnb1 in responses to multiple drugs of abuse is supported by expression data from our study and from other studies. Finally, the data support a role for semaphorin signaling in morphine effects, and indicate that altered expression of genes involved in phosphatidylinositol signaling and plasticity might also affect the altered drug responses in the two strains.