Voluntary wheel running reduces voluntary consumption of ethanol in mice: identification of candidate genes through striatal gene expression profiling.

Hedonic substitution, where wheel running reduces voluntary ethanol consumption, has been observed in prior studies. Here, we replicate and expand on previous work showing that mice decrease voluntary ethanol consumption and preference when given access to a running wheel. While earlier work has been limited mainly to behavioral studies, here we assess the underlying molecular mechanisms that may account for this interaction. From four groups of female C57BL/6J mice (control, access to two-bottle choice ethanol, access to a running wheel, and access to both two-bottle choice ethanol and a running wheel), mRNA-sequencing of the striatum identified differential gene expression. Many genes in ethanol preference quantitative trait loci were differentially expressed due to running. Furthermore, we conducted Weighted Gene Co-expression Network Analysis and identified gene networks corresponding to each effect behavioral group. Candidate genes for mediating the behavioral interaction between ethanol consumption and wheel running include multiple potassium channel genes, Oprm1, Prkcg, Stxbp1, Crhr1, Gabra3, Slc6a13, Stx1b, Pomc, Rassf5 and Camta2. After observing an overlap of many genes and functional groups previously identified in studies of initial sensitivity to ethanol, we hypothesized that wheel running may induce a change in sensitivity, thereby affecting ethanol consumption. A behavioral study examining Loss of Righting Reflex to ethanol following exercise trended toward supporting this hypothesis. These data provide a rich resource for future studies that may better characterize the observed transcriptional changes in gene networks in response to ethanol consumption and wheel running.

Identifying rare variants for genetic risk through a combined pedigree and phenotype approach: application to suicide and asthma.

Suicidal behavior is a complex disorder, with evidence for genetic risk independent of other genetic risk factors including psychiatric disorders. Since 1996, over 3000 DNA samples from Utah suicide decedents have been collected and banked for research use through the Utah Medical Examiner. In addition, over 12,000 Utah suicides were identified through examination of death certificates back to 1904. By linking this data with the Utah Population Database, we have identified multiple extended pedigrees with increased risk for suicide completion. A number of medical conditions co-occur with suicide, including asthma, and this study was undertaken to identify genetic risk common to asthma and suicide. This study tests the hypothesis that a particular comorbid condition may identify a more homogeneous genetic subgroup, facilitating the identification of specific genetic risk factors in that group. From pedigrees at increased risk for suicide, we identified three pedigrees also at significantly increased familial risk for asthma. Five suicide decedents from each of these pedigrees, plus an additional three decedents not from these pedigrees with diagnosed asthma, and 10 decedents with close relatives with asthma were genotyped. Results were compared with 183 publicly available unaffected control exomes from 1000 Genomes and CEPH (Centre d'etude du polymorphisme humain) samples genotyped on the same platform. A further 432 suicide decedents were also genotyped as non-asthma suicide controls. Genotyping was done using the Infinium HumanExome BeadChip. For analysis, we used the pedigree extension of Variant Annotation, Analysis and Search Tool (pVAAST) to calculate the disease burden of each gene. The Phenotype Driven Variant Ontological Re-ranking tool (Phevor) then re-ranked our pVAAST results in context of the phenotype. Using asthma as a seed phenotype, Phevor traversed biomedical ontologies and identified genes with similar biological properties to those known to result in asthma. Our top associated genes included those related to neurodevelopment or neural signaling (brain-derived neurotrophic factor (BDNF), neutral sphingomyelinase 2 (SMPD2), homeobox b2 (HOXB2), neural cell adhesion molecule (NCAM2), heterogeneous nuclear ribonucleoprotein A0 (HNRNPA0)), inflammation (free fatty acid receptor 2 (FFAR2)) and inflammation with additional evidence of neuronal involvement (oxidized low density lipoprotein receptor 1 (OLR1), toll-like receptor 3 (TLR3)). Of particular interest, BDNF has been previously implicated in both psychiatric disorders and asthma. Our results demonstrate the utility of combining pedigree and co-occurring phenotypes to identify rare variants associated with suicide risk in conjunction with specific co-occurring conditions.

Transcriptome analysis of Inbred Long Sleep and Inbred Short Sleep mice.

Many studies have utilized the Inbred Long Sleep and Inbred Short Sleep mouse strains to model the genetic influence on initial sensitivity to ethanol. The mechanisms underlying this divergent phenotype are still not completely understood. In this study, we attempt to identify genes that are differentially expressed between these two strains and to identify baseline networks of co-expressed genes, which may provide insight regarding their phenotypic differences. We examined the whole brain and striatal transcriptomes of both strains, using next generation RNA sequencing techniques. Many genes were differentially expressed between strains, including several in chromosomal regions previously shown to influence initial sensitivity to ethanol. These results are in concordance with a similar sample of striatal transcriptomes measured using microarrays. In addition to the higher dynamic range, RNA-Seq is not hindered by high background noise or polymorphisms in probesets as with microarray technology, and we are able to analyze exome sequence of abundant genes. Furthermore, utilizing Weighted Gene Co-expression Network Analysis, we identified several modules of co-expressed genes corresponding to strain differences. Several candidate genes were identified, including protein phosphatase 1 regulatory unit 1b (Ppp1r1b), prodynorphin (Pdyn), proenkephalin (Penk), ras association (RalGDS/AF-6) domain family member 2 (Rassf2), myosin 1d (Myo1d) and transthyretin (Ttr). In addition, we propose a role for potassium channel activity as well as map kinase signaling in the observed phenotypic differences between the two strains.