DESCRIPTION:
Regression of chromosome 1 eQTL transcripts onto OXY locomotor withdrawal traits. We mapped the genetic basis of OXY behavioral sensitivity and spontaneous withdrawal in a B6J×B6NJ‐F2 reduced complexity cross (RCC). Upon identifying a major distal chromosome 1 QTL for oxycodone locomotion and withdrawal, we implemented an efficient fine mapping strategy using recombinant lines within the QTL interval, and resolved a 2.45‐Mb region. We then sourced a historical striatal expression QTL dataset from the same genetic cross to home in on functionally plausible candidate genes. Finally, we conducted immunoblot analysis of eQTL genes within a recombinant line capturing the behavioral QTL, providing further support for the candidacy of Atp1a2 and Kcnj9. The results provide multiple positional loci and highly plausible candidate genes underlying opioid behavioral sensitivity and withdrawal. 425 F2 mice (213 SAL, 212 OXY) were phenotyped in the OXY locomotor/CPP protocol during Weeks 1 and 2 (Fig.2A). F2 mice were trained in a two‐chamber unbiased CPP protocol over nine days (D), with initial preference (D1, SAL i.p.), two alternating pairings of OXY (1.25 mg/kg, i.p.) and SAL (i.p.), separated by 48h (D2‐D5), a consolidation period (D6‐D7), a drug‐free CPP test following SAL (i.p.) (D8), and a drug‐induced CPP test (D9) following either SAL (i.p.) in the SAL group or OXY (1.25 mg/kg, i.p.) in the OXY group. A subset (118 SAL, 118 OXY) went on to receive a high-dose, multi-week OXY regimen over Weeks 3 and 4 (Fig.2A). For week 3, mice were injected once daily with either SAL (i.p.) or OXY (20mg/kg, i.p.) on Monday through Thursday at 1600h. Sixteen hours later on Friday, the majority of these 236 F2 mice (160F2 mice total: 77F, 83M) were tested for baseline nociceptive sensitivity on the 52.5°C hot plate assay and then assessed for OXY‐induced antinociception (5mg/kg, i.p.) at 30min post‐OXY (see Data S1). The remaining 76 F2 mice (37F, 39M) received the same OXY regimen during Week 3 but were not exposed to hot plate testing. For Week 4, we continued injecting all 236 F2 mice (118 SAL, 118 OXY) once daily with SAL (i.p.) or OXY (20mg/kg, i.p.) for four days. On Day 5, 16h post‐injection, mice were tested for anxiety‐like behavior in the elevated plus maze (EPM). Striatum punches were collected from 23 OXY-trained F2 mice at 24 h after behavioral testing on the EPM. Libraries were prepared according to Illumina's detailed instructions accompanying the TruSeq® Stranded mRNA LT Kit (Part# RS-122-2101). The purified DNA was captured on an Illumina flow cell for cluster generation and sample libraries were sequenced at 23 samples per lane over 5 lanes (technical replicates) on the Illumina HiSeq 4000 machine, yielding an average of 69.4 million total reads per sample. We aligned FASTQ files to the mm38 genome via TopHat (Trapnell et al. 2012) and used the Ensembl sequence and genome annotation. We used featureCounts for read alignment. Exons with low expression (less than 10 reads total across all 115 count files) were removed. A a cis-eQTL was liberally defined as any transcript with a genome-wide significant association between expression and a polymorphic marker that was within 70 Mb of a SNP, given the large linkage disequilibrium in a lowly recombinant F2 cross and given that this was the largest distance between any two SNP markers. In addition to eQTL analysis, we used the limma package and lmfit() to fit a linear mixed model where we regressed transcript levels onto OXY behaviors to identify predictors. Read counts were normalized and log2-transformed using VOOM (Law et al. 2014). Transformed counts were then fit using a linear model using the limma package (Ritchie et al. 2015) and individually regressed normalized, transformed read counts of eQTL transcripts onto OXY locomotor traits (distance, rotations, and spins for both D2 and D4) and/or OXY EPM withdrawal traits [Open Arm Distance (m), Open Arm Entries (#), and Open Arm Time (%)]. Sequence Lane (random effects repeated measure) and Sex were included as covariates. An FDR-adjusted p-value of p<0.05 was employed to identify significant transcripts as predictors of behavior (Benjamini & Hochberg 1995). From supplementary table S5.
LABEL:
Chr1 eQTL OXY D2 locomotor rotations predicters B6JxB6NJ F2 cross_qvalue
SCORE TYPE:
Q-Value
THRESHOLD:
<= 0.5
GENES IN THRESHOLD:
7
DATE ADDED:
DATE UPDATED:
SPECIES:
AUTHORS:
TITLE:
JOURNAL:
ABSTRACT:
Genes in threshold: 7
| Uploaded As | Gene Symbol | Homology | Score | Priority | LinkOuts | Emphasis |
|---|
GeneWeaver