GeneSet Information

Tier IV GS410681 • Chromosome 1 eQTL transcripts mapped to day 2 OXY locomotor traits in B6J×B6NJ‐F2 mouse cross- rotations_qvalue

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:

2025-03-05

DATE UPDATED:

2025-03-05

SPECIES:

AUTHORS:

Lisa R Goldberg, Britahny M Baskin, Jacob A Beierle, Yahia Adla, Julia C Kelliher, Emily J Yao, Stacey L Kirkpatrick, Eric R Reed, David F Jenkins, Jiayi Cox, Alexander M Luong, Kimberly P Luttik, Julia A Scotellaro, Timothy A Drescher, Sydney B Crotts, Neema Yazdani, Martin T Ferris, W Evan Johnson, Megan K Mulligan, Camron D Bryant

TITLE:

Atp1a2 and Kcnj9 Are Candidate Genes Underlying Sensitivity to Oxycodone-Induced Locomotor Activation and Withdrawal-Induced Anxiety-Like Behaviors in C57BL/6 Substrains.

JOURNAL:

Genes, brain, and behavior Feb 2025, Vol 24, pp. e70009

ABSTRACT:

Opioid use disorder is heritable, yet its genetic etiology is largely unknown. C57BL/6J and C57BL/6NJ mouse substrains exhibit phenotypic diversity in the context of limited genetic diversity which together can facilitate genetic discovery. Here, we found C57BL/6NJ mice were less sensitive to oxycodone (OXY)-induced locomotor activation versus C57BL/6J mice in a conditioned place preference paradigm. Narrow-sense heritability of OXY-induced locomotor activity traits ranged from 0.22 to 0.31, implicating suitability for genetic analysis. Quantitative trait locus (QTL) mapping in an F2 cross identified a chromosome 1 QTL explaining 7%-12% of the variance in OXY locomotion and anxiety-like withdrawal in the elevated plus maze. A second QTL for EPM withdrawal behavior on chromosome 5 near Gabra2 (alpha-2 subunit of GABA-A receptor) explained 9% of the variance. To narrow the chromosome 1 locus, we generated recombinant lines spanning 163-181 Mb, captured the QTL for OXY locomotor traits and withdrawal, and fine-mapped a 2.45-Mb region (170.16-172.61 Mb). Transcriptome analysis identified five, localized striatal cis-eQTL transcripts and two were confirmed at the protein level (KCNJ9, ATP1A2). Kcnj9 codes for a potassium channel (GIRK3) that is a major effector of mu opioid receptor signaling. Atp1a2 codes for a subunit of a Na+/K+ ATPase enzyme that regulates neuronal excitability and shows functional adaptations following chronic opioid administration. To summarize, we identified two candidate genes underlying the physiological and behavioral properties of opioids, with direct preclinical relevance to investigators employing these widely used substrains and clinical relevance to human genetic studies of opioid use disorder. PUBMED: 39801366
Find other GeneSets from this publication

Annotation Information

No sequence read archive data associated with this GeneSet.


Genetic Variation (D014644)
Chromosomes, Human, Pair 5 (D002895)
Rotation (D012399)
Oxycodone (D010098)
Crosses, Genetic (D003433)
Linkage Disequilibrium (D015810)
Libraries (D007990)
Divorce (D004243)
Behavior (D001519)
Jurisprudence (D007603)
Chromosomes (D002875)
Chromosomes, Human, Pair 1 (D002878)
Gene Expression Profiling (D020869)
Potassium (D011188)
Quantitative Trait Loci (D040641)
gamma-Aminobutyric Acid (D005680)
Opioid-Related Disorders (D009293)
RNA, Messenger (D012333)
Linear Models (D016014)
Enzymes (D004798)
Proteins (D011506)
Potassium Channels (D015221)
Weights and Measures (D014894)
Receptors, Opioid, mu (D017450)
Sensitivity and Specificity (D012680)
Organization and Administration (D009934)
Receptors, GABA-A (D011963)
Locomotion (D008124)
Exhibitions (D020476)
Adenosine Triphosphatases (D000251)
Receptors, Opioid (D011957)
Anxiety (D001007)
Association (D001244)
Research Personnel (D012108)
Analgesics, Opioid (D000701)
stem cell factor receptor activity (GO:0005020)
locomotion (GO:0040011)
chromosome (GO:0005694)
signaling (GO:0023052)
Alignment construction (EDAM_operation:2928)
Genome annotation (EDAM_operation:0362)
Map position (EDAM_data:2012)
Experiment annotation (EDAM_data:2531)
Quantitative trait locus (EDAM_data:1906)
Transcriptomics (EDAM_topic:0203)
Genome annotation (EDAM_topic:0222)
Nucleic acid features (SNP) (EDAM_data:2092)
Mapping and assembly (EDAM_operation:2429)
Oligonucleotide alignment construction (EDAM_operation:3198)
deoxyribonucleic acid (CHEBI:16991)
sodium(1+) (CHEBI:29101)
salicylate (CHEBI:30762)
gamma-aminobutyric acid (CHEBI:16865)
red chlorophyll catabolite (CHEBI:47899)
protein polypeptide chain (CHEBI:16541)
messenger RNA (CHEBI:33699)
oxy group (CHEBI:29354)
oxycodone (CHEBI:7852)
potassium atom (CHEBI:26216)
alpha-particle (CHEBI:30216)
renal cell carcinoma (DOID:4450)
anxiety disorder (DOID:2030)
DNA shearing (EFO:0004189)
mating type alpha (EFO:0001270)
anxiety disorder (EFO:0006788)
generation (EFO:0000507)
renal cell adenocarcinoma (EFO:0005708)
milligram per kilogram (EFO:0002902)
obsolete_corpus striatum (EFO:0000381)
physical activity (EFO:0003940)
C57BL/6NJ (EFO:0007732)
rectal adenocarcinoma (EFO:0005631)
renal cell carcinoma (EFO:0000681)
transcriptome (EFO:0004421)
locomotor behavior trait (VT:0001392)
forelimb stylopod (UBERON:0003822)
central sulcus of insula (UBERON:0035925)
caudate-putamen (UBERON:0005383)
conditioned place preference (GO:1990708)

Gene List • 8 Genes

Genes in threshold: 7

Uploaded As Gene Symbol Homology Score Priority LinkOuts Emphasis