C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception

Abstract

Sensitivity to different pain modalities has a genetic basis that remains largely unknown. Employing closely related inbred mouse substrains can facilitate gene mapping of nociceptive behaviors in preclinical pain models. We previously reported enhanced sensitivity to acute thermal nociception in C57BL/6J (B6J) versus C57BL/6N (B6N) substrains. Here, we expanded on nociceptive phenotypes and observed an increase in formalin-induced inflammatory nociceptive behaviors and paw diameter in B6J versus B6N mice (Charles River Laboratories). No strain differences were observed in mechanical or thermal hypersensitivity or in edema following the Complete Freund’s Adjuvant model of inflammatory pain, indicating specificity in the inflammatory nociceptive stimulus. In the chronic constrictive nerve injury, a model of neuropathic pain, no strain differences were observed in baseline mechanical threshold or in mechanical hypersensitivity up to one month post-chronic constrictive nerve injury. We replicated the enhanced thermal nociception in the 52.5°C hot plate test in B6J versus B6N mice from The Jackson Laboratory. Using a B6J × B6N-F2 cross (N = 164), we mapped a major quantitative trait locus underlying hot plate sensitivity to chromosome 7 that peaked at 26 Mb (log of the odds [LOD] = 3.81, p < 0.01; 8.74 Mb-36.50 Mb) that was more pronounced in males. Genes containing expression quantitative trait loci associated with the peak nociceptive marker that are implicated in pain and inflammation include Ryr1, Cyp2a5, Pou2f2, Clip3, Sirt2, Actn4, and Ltbp4 (false discovery rate < 0.05). Future studies involving positional cloning and gene editing will determine the quantitative trait gene(s) and potential pleiotropy of this locus across pain modalities.

Keywords

Pain genetics rodents genome-wide reduced complexity cross C57BL/6 substrains

Introduction

Pain is defined by the International Association for the Study of Pain as an unpleasant sensory and or emotional experience associated with actual or potential tissue damage (http://www.iasp-pain.org) and is the main reason patients contact their physician. Nociceptive pain is adaptive and signals potential or actual tissue damage and thus promotes avoidance and protective behavior. However, chronic pain (e.g., neuropathic pain) is maladaptive and significantly decreases quality of life. An estimated one third of the global population suffers from chronic pain, including over 100 million people in the United States.¹ Both environmental and genetic factors contribute to the wide range of variability in pain sensitivity.^2,3 Pain has multiple modalities, including thermal, inflammatory, mechanical hypersensitivity, and neuropathic pain. Each of these modalities are hypothesized to have largely separable genetic components.^4–6 Furthermore, pain has both peripheral and central nervous system components with both shared and divergent molecular and physiological functions. Estimates from preclinical models indicate that 30% to 80% of the variance in pain responses can be explained by genetic factors,⁷ whereas heritability estimates from human experimental pain models range from 20% to 60%.^2,7

Human genome-wide association studies have had limited success in identifying replicable “pain genes,” due in large part to a lack of statistical power and a lack of a consistent definition of pain cases across cohorts.^1,8 Quantitative trait locus (QTL) mapping is an unbiased, discovery-based (i.e., hypothesis-generating) approach to identify polymorphic regions of the genome containing genetic variants underlying complex traits. QTL mapping in preclinical mammalian models of pain provides multiple advantages, including a greater ability to achieve the necessary statistical power and the ability to combine gene expression analysis of relevant tissue (e.g., peripheral nervous system and central nervous system (CNS) tissue) with behavioral QTL analysis to facilitate the identification of quantitative trait genes and functional variants.⁹ QTL studies of nociceptive behaviors associated with pain models in mice can facilitate translation to humans.^1,2 The CACNG2 gene (coding for the gamma 2 subunit of a calcium channel) was first mapped in mice for nerve injury-induced autonomy, using scratching, biting, and licking as the heritable behavioral pain measures for mapping.¹⁰ Haplotypes within CACNG2 were subsequently associated with increased post-mastectomy neuropathic pain.¹ QTL mapping in mice also identified a non-synonymous single-nucleotide polymorphism (SNP) in P2X7 (encoding the adenosine triphosphate–gated P2X7 purinergic receptor) that was associated with reduced allodynia in a spared nerve injury model.¹¹ P2X7 was subsequently associated with post-mastectomy pain in women.¹ QTL mapping in mice also identified calcitonin gene-related peptide (CGRP) in thermal nociception¹² and the MC1R (the gene encoding the melanocortin 1 receptor) in female-specific kappa opioid-induced antinociception—a finding that was subsequently confirmed in women.¹³ As a final example, a QTL for inflammatory pain was mapped to AVPR1A (gene encoding the vasopressin 1A receptor) and in humans, an SNP in AVPR1A was subsequently associated with inflammatory pain and vasopressin receptor-mediated analgesia in non-stressed but not stressed men.¹⁴

Despite successes in translating pain genetics from mice to humans, gene mapping in mice remains a challenge for two reasons. First, genetic complexity underlying a QTL can hinder gene identification. In comparison to the reference C57BL/6J (B6J) inbred strain, a majority of classical inbred strains contain five million or so variants (SNPs plus indels); in stark contrast, the closely related C57BL/6NJ (B6NJ) inbred substrain contains approximately 30,000 variants.^15,16 Thus, the magnitude of genetic complexity underlying a QTL is decreased by orders of magnitude which facilitates the identification of causal genes and variants in these so-called reduced complexity crosses (RCCs).^17–20 The second hurdle to gene identification concerns the low level of resolution of QTLs, especially in F2 crosses, where the confidence interval typically spans one-half of a chromosome and contains hundreds of genes and thousands of variants. Although the use of RCCs does not provide immediate improvement of QTL resolution like other contemporary mouse populations and panels such as highly recombinant outbred stocks,⁹ RCCs can facilitate the subsequent step of fine mapping due to the much simpler genetic architecture of the quantitative traits (typically monogenic inheritance on a nearly isogenic background) by permitting immediate backcrossing and phenotyping at each generation as new recombination events accumulate, thus rapidly narrowing the interval.²⁰ We used C57BL/6 (B6) substrains to map the genetic basis of parental strain variation in binge eating between C57BL/6 substrains.¹⁹ We identified a single QTL on mid-chromosome 11 that mapped to the same region identified for cocaine neurobehavioral sensitivity using the same RCC.¹⁸ The region contains a missense mutation in the gene Cyfip2 that could act as a gain-of-function allele. Accordingly, mice heterozygous for a null mutation in Cyfip2 showed a normalization in binge eating toward a wild-type level, thus providing support for Cyfip2 as a causal genetic factor underlying binge eating in B6 substrains.¹⁸

With respect to pain models, we previously identified behavioral differences in B6 substrains in two assays of acute thermal nociception, including the hot plate and tail withdrawal assays. In both cases, the B6J substrain showed enhanced nociceptive sensitivity relative to the B6N substrains as indicated by a reduction in latency to respond to the thermal nociceptive stimuli.²¹ Subsequent studies replicated the increase in acute thermal nociceptive sensitivity in B6J relative to other B6N.²² However, it is unknown whether enhanced nociceptive sensitivity in the B6J substrain extends to other nociceptive modalities besides acute thermal nociception. A recent report indicated no differences in mechanical sensitivity.²³ Furthermore, because genetic factors are hypothesized to underlie B6 substrain differences in nociceptive behaviors, in the second part of this study, we sought to map the genetic basis of acute thermal nociception in the hot plate assay.²¹ To prioritize functional candidate genes for future gene editing, we used a historical B6J × B6NJ-F2 genomic dataset to report genes within the behavioral QTL for hot plate sensitivity that also possess a cis-expression QTL (eQTL) within striatal brain tissue, including a top differentially expressed candidate gene, the ryanodine receptor 1 (Ryr1). Genes within the locus whose expression is associated with both genetic variation and behavior are prioritized as candidate quantitative trait genes underlying thermal nociception.

Materials and methods

Mice

Adult (8–10 weeks of age at the beginning of experiments) female and male C57BL/6J (B6J) (The Jackson Laboratory (JAX), Bar Harbor, ME, USA) and C57BL/6N mice (B6N) (Charles River Laboratories, Wilmington, MA USA) were used for phenotyping in the inflammatory and neuropathic pain models. Mice were housed in a 21°C humidity-controlled Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC)–approved animal care facility in Virginia Commonwealth University. Mice were housed in standardized shoebox-sized cages in corncob bedding (Teklad, Envigo) in groups of four and had free access to food (Teklad, Envigo 7012 LM-485) and water. The rooms were on a 12-h light/dark cycle (lights on at 7:00 AM). All experiments were performed during the light cycle, and the study was approved by the Institutional Animal Care and Use Committee of Virginia Commonwealth University. All studies were carried out in accordance with the National Institutes of Health’s Guide for the Care and Use of Laboratory Animals. A sample size of n = 12 per substrain (six females and six males) was used for the assays involving B6J versus B6N. If we assume an effect size of Cohen’s d = 1.07 (see the Results section for the effect size estimate for B6 substrain differences in hot plate latencies), this sample size provided us with 70% power to detect a group difference (p < 0.05; two-tailed test).

For experiments involving hot plate testing and QTL analysis (Boston University School of Medicine (BUSM)), C57BL/6J (B6J) and C57BL/NJ (B6NJ) mice were purchased from JAX at seven weeks of age and were habituated in the vivarium one week prior to experimental testing that occurred next door. All behavioral testing was performed during the light phase of the 12 h light/dark cycle (0630 h/1830 h). Testing occurred between 0800 h and 1300 h in an AAALAC-approved, pathogen-free climate-controlled environment (68°C–72°C). Substrains were same-sex housed separately from each other (four per cage) and phenotyped in an alternating fashion. For QTL mapping, B6J females were crossed to B6NJ males to generate B6J × B6NJ-F1 offspring, and B6J × B6NJ F1 mice were intercrossed to generate B6J × B6NJ F2 mice. F2 mice were 50 to 100 days old at the time of testing and were housed two to four per cage in wood shavings in standard shoebox-sized cages. They were fed standard lab chow ad libitum (LabDiet 5V75, St. Louis, MO, USA). All experiments were approved by the Institutional Animal Care and Use Committee of BUSM and completed in accordance with the National Institutes of Health’s Guide for the Care and Use of Laboratory Animals.

Drugs

Complete Freund’s Adjuvant (CFA) was purchased from Sigma-Aldrich (St. Louis, MO). CFA was diluted with mineral oil and administered at a concentration of 10%. Formalin was purchased from Fisher Scientific (Waltham, MA, USA) and diluted with distilled water to a 2.5% concentration. Formalin was prepared daily.

Formalin test of inflammatory nociception

Female and male B6J and B6N mice (n = 6/sex/group) were used to study formalin nociceptive behaviors. The formalin test was carried out in an open Plexiglas cage (29 × 19 × 13 cm each). Mice were allowed to acclimate for 5 min in the test cage prior to injection. Each animal was injected with 20 μL of formalin (2.5%) to the dorsal surface of the right hind paw. Mice were observed from 0 to 5 min (early phase) and 20 to 45 min (late phase) post-formalin injection. The amount of time spent licking the injected paw was recorded with a digital stopwatch. Paw diameter (see measurement of paw edema) was also determined before and 1 h after formalin injection.

CFA-induced chronic inflammatory pain model

Female and male B6J and B6N mice were used (n = 6/sex/group) to study mechanical and thermal hypersensitivity after being administered CFA. Mice were injected in the dorsal surface of the right hind paw with 20 µL of CFA (10%). Mechanical hypersensitivity was measured via the von Frey test before injection and on days 3, 7, 14, 21, and 28 after CFA injection. On same days, paw edema was also evaluated. Thermal hypersensitivity was measured via the Hargreaves test on days 4, 10, 17, and 24 after CFA injection. The results shown pertain to the ipsilateral paw of each mouse.

Chronic constrictive nerve injury-induced neuropathic pain model

To evaluate possible differences in the development of chronic constrictive nerve injury (CCI)-induced neuropathic pain-like behavior, female and male B6J and B6N (n = 6/sex/group) were tested in CCI neuropathic pain model. Mice were anesthetized using 3% to 4% isoflurane and maintained with 1% to 2% isoflurane in oxygen using a face mask and a vaporizer (VetEquip Inc, Pleasanton, CA, USA). An incision was made just below the left hip bone, parallel to the sciatic nerve. The left common sciatic nerve was exposed at the level proximal to the sciatic trifurcation and a nerve segment 3 to 5 mm long was separated from surrounding connective tissue. Two loose ligatures with 5–0 silk suture were made around the nerve with a 1.0 to 1.5 mm interval between each of them. The wound was closed with suture thread. This procedure resulted in CCI of the ligated nerve and induced neuropathic insult. Animals were tested for changes in mechanical hypersensitivity before surgery and at days 3, 5, 7, 11, 14, 21, 28, and 35 post-surgery. The results shown pertain to the ipsilateral paw of each mouse.

Mechanical hypersensitivity

Mechanical withdrawal thresholds were determined according to the method of Chaplan et al.²⁴ with slight modifications.²⁵ Mice were placed in clear plastic cylinders (9 × 11 cm) with mesh metal flooring and allowed to acclimate for 15 min before testing. A series of calibrated von Frey filaments (Stoelting, Inc., Wood Dale, IL, USA) with logarithmically incremental stiffness ranging from 2.83 to 4.56 units expressed Log 10 of (10 × force in (mg)) were applied to the paw using a modified up-down method.²⁵ Testing commenced with the 3.84 numbered filament for each mouse and continued up or down according to response. In the absence of a paw withdrawal response to the initially selected filament, a thicker filament corresponding to a stronger stimulus was presented. In the event of paw withdrawal, the next weaker stimulus was chosen. Each hair was presented perpendicularly against the paw, with sufficient force to cause slight bending and held 2 to 3 s. The stimulation of the same intensity was applied two times to the hind paw at intervals of a few seconds. The mechanical threshold (paw withdrawal threshold (PWT)) was expressed in grams, indicating the force of the von Frey hair to which the animal reacted (paw withdrawn, licking, or shaking). The data are expressed with the results of the ipsilateral paw of each mouse.

Thermal hypersensitivity

Thermal withdrawal latencies were measured via the Hargreaves test by placing the mice in clear plastic cylinders (9 × 11 cm) with glass plate flooring and allowing them to acclimate for 15 min before testing. An infrared heat emitter was used to detect the thermal hypersensitivity at an intensity of 2.8. The infrared emitter was placed under the paw of the mouse, and the paw withdrawal latency (PWL) was recorded for both the each hind paw of each mouse. An average of 2 to 3 measures of the PWL was taken for each hind paw. Results are shown for the ipsilateral paw of each mouse.

Paw edema

The thickness of the formalin- or CFA-treated paws was measured both before and after injections on certain time points using a digital caliper (Traceable Calipers, Friendswood, TX, USA). Data were recorded to the nearest ±0.01 mm and expressed as change in paw thickness (ΔPD = difference in the ipsilateral paw diameter before and after injection paw thickness).

Hot plate test of thermal nociception

B6J (n = 24; 12 females and 12 males) and B6NJ (n = 18; 11 females and 7 males) parental substrains as well as B6J × B6NJ-F2 mice (N = 164) were tested on the 52.5°C hot plate assay.^26,27 Mice were habituated to the testing room for at least 1 h. Mice were then placed in a Plexiglas cylinder (15 cm diameter; 33.0 cm tall) on a hot plate (IITC Life Science Inc., Woodland Hills, CA, USA) and the latency to lick the hind paw with a 60 s cut-off latency was recorded using a stopwatch. The B6J and B6NJ parental substrains were experimentally naïve at the time of nociceptive assessment. F2 mice (N = 164; 78 females and 86 males) were part of a larger historical dataset in which mice had a prior history of training in a conditioned place preference protocol as described²⁸ that involved two training injections of either saline (SAL; i.p.; n = 83) or the mu opioid receptor agonist oxycodone (OXY; 1.25 mg/kg, i.p.; n = 81). During the four days prior to baseline hot plate assessment, F2 mice had continued to receive four daily injections of either SAL (i.p.) or OXY (20 mg/kg, i.p.). Twenty-four hours after the fourth injection, mice were assessed for baseline hot plate sensitivity.

All behavioral testing was conducted in a manner where the experimenter was blinded to the B6 substrain or genotype (F2 mice for QTL analysis). B6 substrains were always tested together within a cohort and cages were tested in an alternating manner between the two substrains.

DNA collection and genotyping in F2 mice

DNA was extracted from spleens and prepared for genotyping using a standard salting out protocol. Ninety SNP markers spaced approximately 30 Mb (approximately 15 cM) apart were genotyped using a custom-designed Fluidigm array (South San Francisco, CA, USA). This level of coverage is sufficient for initial mapping in an F₂ cross because of the low number of recombination events. We recently used a nearly identical marker panel for QTL mapping in the same cross.¹⁹ We included an additional marker on chromosome 7 (rs31995355; 4 Mb) to improve QTL resolution. Markers were selected using high coverage sequence data for the C57BL/6NJ strain generated by the Welcome Trust Sanger Institute^15,29 and were validated using traditional Sanger sequencing. Genomic DNA was diluted to 100 ng/uL in low TE buffer (Teknova, Hollister, CA, USA) and genotyped with approximately 20% replication using the Fluidigm 96 × 96 SNPtype assay according to the manufacturer instructions. SNPs were called using the Fluidigm SNP Genotyping Analysis Software and SNPtype Normalization with the default 65% confidence threshold.

RNA collection, library preparation, and sequencing for eQTL analysis

Striatum punches were collected as described³⁰ for RNA-seq from 23 F2 mice (all OXY-treated for historical reasons). Brains were rapidly removed and sectioned with a brain matrix to obtain a 3-mm-thick section where a 2 mm diameter punch of the striatum was collected. Left and right striatum punches were pooled and immediately placed in RNAlater (Life Technologies, Grand Island, NY, USA) for 48 h prior to storage in a −80°C freezer. Total RNA was extracted using the RNeasy kit (Qiagen, Valencia, CA, USA) as described.³⁰ RNA was shipped to the University of Chicago Genomics Core Facility for cDNA library preparation using the Illumina TruSeq (oligo-dT; 100 bp paired-end reads). Libraries were prepared according to Illumina’s detailed instructions accompanying the TruSeq® Stranded mRNA LT Kit (Part# RS-122–2101). The purified cDNA was captured on an Illumina flow cell for cluster generation and sample libraries were sequenced at 23 samples per lane over five lanes (technical replicates) according to the manufacturer’s protocols on the Illumina HiSeq4000 machine, yielding an average of 69.4 million reads per sample. FASTQ files were quality checked via FASTQC and possessed Phred quality scores > 30 (i.e., less than 0.1% sequencing error).

Statistical analysis

Behavioral data from the formalin test, CFA, and CCI were analyzed using the GraphPad Prism software, version 6.0 (GraphPad Software, Inc., La Jolla, CA, USA) and expressed as the mean ± S.E.M. Statistical analysis was conducted using analysis of variance (ANOVA) test and followed by a post hoc test. Before ANOVA, the data were first assessed for the normality of the residuals and equal variance. Variances were similar between groups and were assessed using either the F‐test or the Brown–Forsythe test and the Bartlett’s test. All data passed these tests. A two-way ANOVA followed by the Tukey’s post hoc correction was used in formalin test. In addition, an unpaired student t test was used to compare paw thickness in the formalin test. To test the mechanical and thermal sensitivity, repeated measures (RMs) two-way ANOVA was used with Sidak (to compare substrain effects) and Dunnett’s (to compare time’s effect according to baseline value) post hoc correction. A p value of less than 0.05 was considered significant.

QTL analysis was performed in F₂ mice using the R package R/qtl as previously described.^19,31,32 The scanone function was used to calculate log of the odds (LOD) scores. Permutation analysis (perm = 1000) was used to establish the significance threshold for each quantitative trait (p < 0.05). The marker position (cM) was estimated using the sex-averaged position using the Mouse Map Converter (http://cgd.jax.org/mousemapconverter).³³ The converter tool was also used to estimate the Mb position of the QTL peak and the Bayes credible interval that defined the QTL region. The 90 polymorphic markers used for genetic mapping and the reported variants that underlie the hot plate QTL interval were retrieved from the Sanger database (http://www.sanger.ac.uk/). Power analysis was used to inform the F2 sample size for QTL analysis and was conducted using GPower (http://www.gpower.hhu.de/en.html). Using the means, standard deviations, and unpaired two-tailed t tests, we calculated the effect size (Cohen’s d) of the parental B6 substrain difference in hot plate latency. We then used this effect size and set the alpha level to 0.05 and power level to 95% to estimate the required sample size for each homozygous genotype at a given locus.

For eQTL mapping, we aligned FastQ files to the reference genome (mm38) via TopHat³⁴ using the mm38 build and Ensembl Sequence and genome annotation. We used featureCounts to count and align reads. For cis-eQTL analysis, we used the same 90 SNPs that were used in behavioral QTL analysis. We removed lowly expressed exons that did not possess at least 10 reads total across all 115 count files. Because of the low resolution of QTL mapping in an F2 cross, we liberally defined a gene with a cis-eQTL as any gene possessing a genome-wide significant association between expression and a polymorphic marker that was within 70 Mb of an SNP (the largest distance between any two SNPs from the 90-SNP panel). The genes reported within the chromosome 7 interval all showed their most significant association with gene expression at the peak SNP associated with hot plate sensitivity (rs3148686; 30.31 Mb). Analysis was conducted using limma with default trimmed mean of M values (TMM) normalization and variance modeling at the observational level (VOOM) transformation.^35,36 A linear model was employed whereby sample replicates were treated as a RM. The duplicateCorrelation()function was used to estimate within-sample correlation, which was then included in the lmFit() function. An ANOVA test was conducted for gene expression that included sex as a covariate and genotype as a fixed effect. Gene-level tests were conducted using the likelihood ratio test. A false discovery rate (FDR) of 5% was employed as the cut-off for statistical significance.³⁷

Results

Formalin-induced inflammatory nociceptive behaviors in B6J and B6N substrains

As seen in Figure 1(a), we evaluated the paw licking responses of B6J and B6N mice in the formalin test at 2.5% concentration. Two-way ANOVA revealed significant effects on time spent in paw licking time in terms of substrain (F_strain(1,44) = 12.91, p < 0.001), phase (F_phase(1,44) = 160.1, p < 0.001), and interaction (F_interaction(1,44) = 4.68, p < 0.05). Although no significant difference was found in phase I behaviors between B6J and B6N mice (p > 0.05, Figure 1(a)), a significant increase in paw licking time in B6J mice compared to B6N mice was found in the phase II (p < 0.001, Figure 1(a)). Moreover, the degree of paw edema (thickness) in B6J mice was significantly greater than B6N mice at 1 h post-formalin injection (t(22) = 4.163, p < 0.001, Figure 1(b)).

Figure 1

Formalin-induced paw licking behavior and paw edema in B6J and B6N mice.

Supplementary Figure 1 shows the effect of formalin on paw licking and edema in females and males. A reduction in time spent in paw licking was found in male B6N mice (F_substrain(1,20) = 9.83, p < 0.01), phase (F_phase(1,20) = 67.36, p < 0.001), and interaction (F_interaction(1,20) = 3.78, p = 0.66) (Supplementary Figure 1(A)). Consistent with the paw licking behavior results, paw diameter was lower in B6N male mice compared to B6J males (t = 4.163, df = 22, p < 0.001, Supplementary Figure 1(B)). However, there was no significant difference on paw licking (F_substrain(1,20) = 3.08, p = 0.094), phase (F_phase(1,20) = 88.28, p < 0.001), and interaction (F_interaction(1,20) = 0.9843, p = 0.33) (Supplementary Figure 1(C)) and paw edema between B6J female and B6N female mice (t(22) = 4.16, p < 0.001, Supplementary Figure 1(D)).

CFA-induced mechanical and thermal hypersensitivity in B6J and B6N substrains

Mice were given an injection of CFA (10%) and tested for mechanical and thermal sensitivity at several time points (days) before and after post-CFA injection. In addition, paw edema was evaluated. For mechanical responses, two-way ANOVA showed significant effects of time (F_time(6,66) = 22.57, p < 0.001) but not for substrain (F_substrain(1,11) = 0.4384, p = 0.52) or interaction (F_interaction(6,66) = 0.33, p = 0.92) (Figure 2(a)). Prior to CFA injection, B6J and B6N mice did not differ in mechanical baseline values (p > 0.05, Figure 2(c)), similar to observations from a recent study.²³ The CFA injection induced a decrease in mechanical thresholds which indicates an induction of mechanical hypersensitivity by CFA in both substrains (Figure 2(a)). Significant mechanical hypersensitivity started by day 3 post-CFA injection (p < 0.05) and continued until day 21 (p > 0.05, Figure 2(a)). Furthermore, we tested thermal paw withdrawal latencies on different days than PWT testing. Two-way ANOVA showed significant effects for thermal responses for time (F_time(4,44) = 41.05, p < 0.001) but not for substrain (F_substrain(1,11) = 3.25, p = 0.099) and interaction (F_interaction(4,44) = 0.82, p = 0.52) (Figure 2(b)). Prior to CFA injection, baseline values for paw withdrawal latencies were similar between B6J and B6N mice (p > 0.05, Figure 2(b)). CFA injection resulted in a significant reduction in paw withdrawal latencies which indicates the formation of thermal hypersensitivity (p < 0.05). The thermal hypersensitivity was observed on days 11 and 17 after CFA injection (p < 0.05); mice were recovered from thermal hypersensitivity by day 24 (p > 0.05, Figure 2(b)). CFA-induced paw edema did not differ between B6J and B6N mice (F_time(5,55) = 35.14, p < 0.001; F_substrain(1,11) = 0.66, p = 0.4338; F_interaction(5,55) = 1.72, p = 0.15; Figure 2(c)).

Figure 2.

CFA-induced mechanical & thermal hypersensitivity and paw edema in B6J and B6N mice. Differences in (A) mechanical paw withdrawal thresholds and (B) paw withdrawal latencies (Δ PWL=contralateral–ipsilateral hindpaw latencies) in B6J and B6N mice at different times after injection of complete Freund's adjuvant (CFA, 10 % solution/20 µl). Degree of edema (C), as measured by the difference in the ipsilateral paw diameter before and after injection (ΔPD) in B6J and B6N mice. Data are expressed as the mean ± S.E.M. of 6 mice/per sex/per group. *p<0.05 significantly different from the value of baseline (BL).

Supplementary Figure 2 shows a comparison of CFA-induced mechanical and thermal hypersensitivity in B6J and B6N substrains for females and males. Two-way ANOVA showed no significant substrain differences in mechanical withdrawal thresholds (F_time(6,30) = 14.11, p < 0.001; F_substrain(1,5) = 6.155, p = 0.056; F_interaction(6,30) = 0.18, p = 0.98; Supplementary Figure 2(A)) or in thermal withdrawal latencies (F_time(4,20) = 23.17, p < 0.001; F_substrain(1,5) = 0.91, p = 0.39; F_interaction(4,20) = 0.76, p = 0.56; Supplementary Figure 2(B)). Similarly, there was no significant substrain differences in mechanical withdrawal thresholds (F_time(6,30) = 10.58, p < 0.001; F_substrain(1,5) = 1.11, p = 0.34; F_interaction(6,30) = 0.26, p = 0.95; Supplementary Figure 2(C)) or in thermal withdrawal latencies (F_time(4,20) = 16.28, p < 0.001; F_strain(1,5) = 2.632, p = 0.17; F_interaction(4,20) = 4.41, p < 0.05; Supplementary Figure 2(D)) between female B6J and B6N mice.

CCI-induced mechanical hypersensitivity in B6J and B6N substrains

The development of CCI-induced mechanical hypersensitivity in von Frey test was compared in B6J and B6N mice. In examining mechanical thresholds over time, two-way ANOVA indicated a main effect of time (F_time(8,88) = 164.4, p < 0.001) but not substrain (F_strain(1,11) = 0.74, p = 0.41) or interaction (F_interaction(8,88) = 1.29, p = 0.26) (Figure 3). Prior to CCI surgery, B6J and B6N mice did not differ in mechanical baseline values (p > 0.05, Figure 3). A reduction in mechanical PWTs was observed at three days post-surgery and continued over 35 days (p < 0.05; Figure 3). The mechanical hypersensitivity was similar between substrains. Moreover, male (F_time(8,40) = 62.87, p < 0.001; F_substrain(1,5) = 1.401, p = 0.2898; F_interaction(8,40) = 1.347, p = 0.2491; Supplementary Figure 3(A)) and female (F_time(8,40) = 106.5, p < 0.001; F_substrain(1,5) = 0.0696, p = 0.8024; F_interaction(8,40) = 0.4246, p = 0.8993; Supplementary Figure 3(B)) showed similar mechanical hypersensitivity in B6J versus B6N substrains.

Figure 3.

CCI-induced mechanical and thermal hypersensitivity in B6J and B6N mice. Differences in mechanical paw withdrawal thresholds in B6J and B6N mice at different days after chronic constrictive nerve injury (CCI) operation. Data are expressed as the mean ± S.E.M. of 6 mice/per sex/per group. *p<0.05 significantly different from the value of baseline (BL).

Increased sensitivity to acute, thermal nociception in the hot plate test in B6J versus B6NJ substrains

We and others previously reported an increase in acute thermal nociceptive sensitivity (i.e., decreased baseline latencies) in the B6J substrain versus the B6N (Crl) substrain.^21,22 Here, we extended this finding in comparison with another genetically very similar B6N “sub”-substrain³⁸ for which whole genome sequence information is available^15,16,29—the B6NJ substrain from JAX (Figure 4(a), t(40) = 3.59; p = 8.9 × 10⁻⁴). When breaking down the data by females and males, the difference was significant with females (t(21) = 2.08; p = 0.0075; n = 12 B6J, n = 12 B6NJ) but did not reach the p < 0.05 cut-off (two tailed) for significance in males (t(17) = 2.11; p = 0.078), likely due the smaller sample size (n = 11 B6J, n = 7 B6NJ; Supplementary Figure 4(A)). However, based on our historical findings and others in males^21,22 a one-tailed test is justified for the males and yields the expected significant decrease in B6J versus B6NJ males (t17 = 1.74; p = 0.039) (Supplementary Figure 4(A)).

Figure 4.

A major QTL on chromosome 7 underlies B6 substrain differences in acute, thermal nociception in the 52.5°C hot plate assay. (A): We replicated our previous observation(21) of a significant increase in sensitivity to acute thermal nociception (decreased in latency to lick the hindpaw) in the parental B6J substrain (n=24; 12 females, 12 males) versus the B6NJ substrain (n=18; 11 females, 7 males) (t40 = 3.59; p = 8.9 x 10-4). (B): Genome-wide significant QTL on chromosome 7 was identified from 164 B6J x B6NJ-F2 mice [LOD = 3.81; peak = 17 cM (6 Mb); peak marker (rs3148686) = 30 Mb; Bayes C.I.: 8.74 Mb – 36.50 Mb). (C): QTL plot for chromosome 7 is shown. Horizontal line for panels B and C denotes the significance threshold (1000 permutations). (D): Effect plot at the peak-associated marker illustrates the decreased hot plate latency associated with the B6J allele (B) and an additive effect with inheritance of one copy versus two copies of the B6NJ allele (N). BB = homozygous for B6J allele; BN = heterozygous; NN = homozygous for B6NJ allele. Data for panels A and D are expressed as the mean ± SEM.

Identification of a major QTL on chromosome 7 underlying enhanced acute thermal nociceptive sensitivity in B6J versus B6NJ substrains

Power analysis using the effect size estimated from the hot plate latencies of the parental substrains (Cohen’s d = 1.07) indicated that a sample size of n = 24 homozygous genotypes was required to achieve 95% power (p < 0.05). We employed a sample size of N = 164 which yielded n = 39 B6J and n = 42 B6NJ genotypes at the peak-associated marker and is within range of the predicted n = 41 per homozygous genotype based on a 1:2:1 Mendelian ratio. Using this sample size, we achieved greater than 99% power to detect a significant difference of the predicted effect size of d = 1.07 (p < 0.05). If we adjust the alpha level to 0.00056 to account for the 90 statistical tests (0.05/90), we still achieve 88% power (p = 0.00056) to detect an effect of d = 1.07, given our homozygous genotype sample sizes of n = 39 homozygous B6J and n = 42 homozygous B6NJ at the peak QTL marker.

A robust phenotypic difference in hot plate sensitivity combined with availability of polymorphic markers and high throughput genotyping tools equipped us with the ability to accomplish our longstanding goal of conducting a genome-wide QTL study of baseline thermal nociception in a cross between B6 substrains.¹⁷ A list of the 90 polymorphic markers used in QTL analysis is provided in Supplementary Table 1. We identified a single genome-wide locus on chromosome 7 underlying B6 substrain differences in baseline nociception as measured via the hot plate assay (LOD = 3.81; QTL peak = 14 cM (26.14 Mb); Bayes Credible Interval = 5.30 cM—22.30 cM (8.74 Mb–36.50 Mb); Figure 4(b) and (c)). The peak marker, rs3148686, is located at 30.31 Mb. Importantly, the QTL signal was driven by allelic differences in the predicted phenotypic direction (Figure 4(d) vs. Figure 4(a); Supplementary Figure 4(A) to (F)). Mice with a prior history of OXY versus SAL injections showed similarly trending QTLs and genotypic pattern of hot plate latencies or QTL detection (Supplementary Figure 5). Furthermore, inclusion of Prior Treatment as a covariate into the QTL model still yielded a significant LOD score (LOD with covariate = 3.72, p < 0.01). The effect size of this locus in comparing the homozygous phenotypes was d = 0.87, and thus, we achieved 70% power (p < 0.00056). To summarize, we identified a major genetic locus underling enhanced acute thermal pain sensitivity in B6J versus B6N mice. Surprisingly, when we conducted QTL analysis separately in females and males, the chromosome 7 QTL was only significant in males.

Striatal cis-eQTLs within chromosome 7 QTL for hot plate sensitivity

The RNA-seq gene expression dataset has been uploaded to Gene Expression Omnibus (GEO). The GEO accession number is GSE119719.

Complementary eQTL analysis can identify functionally relevant candidate genes, providing a link between genotype, gene expression, and behavior.³⁹ We took advantage of a historical dataset that we generated from striatal brain tissue collected from 23 OXY-treated F2 mice in order to document genes within the chromosome 7 QTL that also possess cis-eQTLs. OXY-treated mice also show variation in acute thermal pain sensitivity that maps to the same chromosome 7 hot plate QTL (Supplementary Figure 5); thus, gene expression in OXY-treated mice is still relevant to the behavioral QTL.

The genes we identified were located within 1 Mb to 17 Mb distance of the peak-associated marker for hot plate latency. We identified 15 genes whose expression was modulated by cis-eQTLs (FDR < 0.05) (Table 1). Eight of these genes contain known variants (Ryr1, Rps5, Cox6b1, Pou2f2, Clip3, Sirt2, Actn4, and Ltbp4)¹⁵ (Supplementary Table 2) and nine of these genes have been implicated in either pain and/or inflammation (Ryr1, Cyp2a5, Rps5, Gys1, Pou2f2, Clip3, Sirt2, Actn4, and Ltbp4). Based on the distance from the peak-associated SNP (1 Mb), the strength of the association (FDR = 1.5 × 10⁻⁵), the number of variants overlapping the cognate gene (9 noncoding SNPs), and representation in the pain literature, we nominated the ryanodine receptor 1 (Ryr1; 29 Mb) as a high priority candidate gene underlying enhanced hot plate sensitivity. The nociception-enhancing B6J allele was associated with a 1.25-fold decrease in Ryr1 expression (or i.e., the B6NJ allele showed an increase in expression; Table 1). Only two other genes besides Ryr1 showed a larger fold-change in expression, including two cytochrome P450 genes—Cyp2a5 and Cyp2g1 in which the B6J allele showed a 3.6- to 3.8-fold decrease in expression relative to the B6NJ allele (or i.e., B6NJ allele showed an increase in expression; Table 1). Cyp2a5 (26.84 Mb) and Cyp2g1 (26.81 Mb) are located very near to peak marker associated with the chromosome 7 QTL for hot plate latency that was inferred from interval mapping (26 Mb).

Table 1.

Genes within the hot plate QTL on chromosome 7 that have a cis-eQTL associated with rs3148686.

Chr	Gene	Name	Start (Mb)	Dist. (Mb)	Cod.	N.Cod.	P value	FDR	FC H	FC N	Pain	Inflam.	129P2?
7	Ryr1	Ryanodine receptor 1	29	1.18	0	9	1.22E-07	1.50E-05	1.25	1.02	Yes	Yes	Yes
7	Erf	ETS2 repressor factor	25.24	5.06	0	0	2.98E-06	1.80E-04	1.23	1.14	No	No	No
7	Cyp2a5	Cytochrome P450 family 2 subfamily A member 5	26.84	3.46	0	0	3.16E-05	1.00E-03	2.79	3.82	Yes	Yes	No
7	Rps5	Ribosomal protein S5	12.92	17.38	0	2	1.20E-04	0.0027	−1	1.09	No	Yes	Yes
7	Gys1	Glycogen synthase 1	45.44	15.13	0	0	2.20E-04	0.0044	−1.17	−1.26	Yes	Yes	No
7	Rps11	Ribosomal protein S11	45.12	14.82	0	0	4.90E-04	0.0081	1.02	1.09	No	No	No
7	Cox6b1	Cytochrome C oxidase subunit 6B1	30.62	3.1	0	1	5.00E-04	0.0081	1.04	1.19	No	No	No
7	Pou2f2	POU class 2 homeobox 2	25.09	5.17	0	6	1.00E-03	0.014	1.1	−1.06	Yes	Yes	Yes
7	Clip3	CAP-Glydomain containing linker protein 3	30.29	0.19	0	2	1.30E-03	0.017	−1.07	−1.04	No	Yes	Yes
7	Sirt2	Sirtuin 2	28.77	1.52	0	2	2.30E-03	0.025	−1.12	−1.06	Yes	Yes	No
7	Lmtk3	Lemur tyrosine kinase 3	45.78	15.48	0	0	2.50E-03	0.026	1.16	1.09	No	No	No
7	Actn4	Actinin alpha 4	28.89	1.34	0	3	3.70E-03	0.034	1.1	1.12	Yes	Yes	No
7	Syt3	Synaptotagmin 3	44.38	14.08	0	0	4.60E-03	0.04	1.15	1.14	No	No	No
7	Ltbp4	Latent transforming growth factor beta binding protein 4	27.31	2.97	5'UTR	1	0.0052	0.044	1.21	1.14	No	Yes	Yes
7	Cyp2g1	Cytochrome P450 family 2 subfamily G member 1	26.81	3.49	0	0	0.0053	0.044	2.15	3.55	No	No	No

Dist.(Mb): distance of gene from the SNP rs3148686 (30.31 Mb), the SNP nearest the QTL for hot plate sensitivity; Cod.: coding polymorphisms; N.Cod.: non-coding polymorphisms; FDR: false discovery rate; FC H: fold-change in gene expression in mice heterozygous (H) for the rs3148686 marker relative to mice homozygous for the B6J (B) allele; FC N: fold-change in gene expression in mice homozygous for the B6NJ allele (N) relative to mice homozygous for the B allele; Pain: pain literature search (“gene” and “pain” or “nociception” PubMed); Inflam.: inflammation literature search (“gene” and “inflammation” or “inflammatory” PubMed). 129P2?: whether or not the gene contains a variant shared between the nociception-resistant B6NJ strain and the 129P2 strain; QTL: quantitative trait locus.

Discussion

We expanded the catalog of behavioral differences in nociceptive sensitivity between B6 substrains^21,22 and identified enhanced inflammatory nociception in the formalin test in the B6J versus B6N substrain (Figure 1) in the absence of any behavioral differences in the more chronic CFA inflammatory pain model (Figure 2) or in the CCI neuropathic pain model (Figure 3). Thus, there is selectivity of enhanced nociceptive sensitivity in the B6J strain with regard to both the nociceptive modality as well as the particular noxious stimulus within a nociceptive modality (inflammatory). In the second part of the study, we replicated the enhanced sensitivity to acute thermal nociception as measured via the hot plate assay in B6J versus B6N substrains from JAX (Figure 4(a)) and mapped a major QTL underlying parental substrain differences on chromosome 7 (Figure 4(b) and (c)) that mirrored the parental substrain difference (Figure 4(d)). Finally, using a historical striatal eQTL dataset, we identified 15 genes possessing eQTLs within 1 to 17 Mb of the peak-associated marker, providing functional support for their candidacy.

The B6J substrain showed a more pronounced and statistically significant inflammatory nociceptive response during the late phase of the formalin test (Figure 1(a)). The early phase of nociceptive behavior (e.g., 0–5 min) is thought to model acute chemical pain and activation of c-fibers, whereas the late phase (e.g., 20–45 min) is thought to model an inflammatory response and functional changes in the dorsal horn of the spinal cord.⁴⁰ The larger increase in paw diameter following the formalin test in the B6J substrain further supports an increased inflammatory response (Figure 1(b)). Male B6J mice showed a more robust increase in nociceptive behaviors relative to their B6N male counterparts compared to females during the late phase (Supplementary Figure 1); however, both female and male B6J mice showed a significant increase in paw diameter relative to their B6N counterparts. The former observation could be due to chance sampling error and/or to the effects of sex steroids on formalin-induced nociceptive behaviors^41–44 that obfuscate detection of genetic differences.

Interestingly, the genome-wide significant QTL on chromosome 7 that we identified for hot plate thermal nociception was only significant for male mice (Supplementary Figure 4(B) to (D)). The reason for this observation is not entirely clear. Our previously published results and others that showed enhanced thermal nociception in B6J versus B6N mice were collected in males only.^21,22 However, in the present study, our parental substrain results were generated from both females and males and when considering only females, we also observed a significant decrease in hot plate latency in B6J versus B6NJ (Supplementary Figure 4(A)). Thus, the lack of detection of a QTL from female mice cannot be explained by a lack of a priori evidence in females or, for example, by the fact that we used a different B6N substrain compared to our previous study.²¹ Given the trending effect plots in for the females in the QTL study (Supplementary Figure 4(F)), one could speculate that with a larger sample size we may have detected a significant QTL. Additionally, we did not monitor the estrus cycle which is known to influence nociception, for example, hot plate nociception⁴⁵ and tail flick nociception⁴⁶ and chronic/persistent pain,⁴⁷ and thus, variation in the phase of the estrus cycle could affect nociceptive sensitivity and obscure the effect of genotype. The sex of the experimenter is known to influence nociceptive latencies in mice⁴⁸; however, in this case, the experimenters for phenotyping the parental substrains (S.L.K.) and the F2 mice (L.R.G.) were both females. Finally, it is possible that a larger sample size might have revealed a genome-wide, female-specific locus.

The QTL on chromosome 7 overlaps with a QTL previously identified for acute thermal nociceptive sensitivity in the 49°C tail withdrawal assay in an F2 cross between B6J and 129P3 strains.⁴⁹ The QTL peak for tail withdrawal sensitivity was more distally located at 33 cM (C.I.: 24–38 cM) but was also associated with enhanced pain sensitivity in mice possessing the B6J allele. To our knowledge, the 129P3 strain has not been whole genome sequenced, and thus, a complete catalog of variants is lacking. However, whole genome sequence data are available for the 129P2/OlaHsd substrain (https://www.sanger.ac.uk/), one of the founder 129 strains⁵⁰ that is genetically similar to 129P3. Interestingly, C57BL/6NJ and 129P2 strains share variants within many of the genes listed in Table 1 that possess cis-eQTLs, including Ryr1, Rps5, Pou2f2, Clip3, and Ltbp4.^15,16 Thus, if a private variant in B6J underlies enhanced sensitivity to acute thermal nociception and if B6NJ and 129P2 strains (and by further extension, the 129P3 strain) share the alternate allele via identity by descent,²⁰ then eQTL-containing genes that possess these SNPs are high priority candidates to pursue via gene targeting.

Interestingly, the chromosome 7 QTL we identified for hot plate nociception is localized to a nearly identical region reported for fluid intake of the nociceptive compound capsaicin (activates TRPV1) in an F2 cross between C57BL/6J and the wild-derived KJR strain.⁵¹ The QTL was localized to approximately 19 cM (peak marker = D7Mit155; 31 Mb) which is nearly identical to the location for the hot plate QTL (peak marker = 30 Mb). The B6J allele was associated with decreased fluid intake of capsaicin compared to the KJR allele which is consistent with enhanced capsaicin-induced oral nociceptive sensation and is consistent with the overall theme of enhanced B6J-mediated nociception at this locus (Figure 4).⁴⁹ Interestingly, a large panel of inbred strains tested alongside the C57BL/6 and KJR strains showed a similar pattern of strain differences in nociceptive sensitivity to capsaicin and hot plate,⁵² indicating a genetic correlation and, thus, a shared genetic basis between noxious chemical and thermal stimuli that generate heat sensation in their strain panel.

We examined gene expression in the striatum since these data were already available from a historical dataset. Arguably, because the nociceptive response in the hot plate assay is a supraspinally controlled, thermal acute response,²⁷ the relevant tissue for examining gene expression as it relates to thermal nociception could be a different tissue such as the dorsal root ganglia, the dorsal horn of the spinal cord, brainstem nuclei, midbrain periaqueductal gray, primary sensory cortex, or any of one of several limbic structures involve in supraspinal modulation of pain. In relation to the ventral striatum, CGRP, which is an established pain neuropeptide that has been genetically mapped using inbred mouse strains,¹² induces an increase in hind paw withdrawal latency in the hot plate when injected into the nucleus accumbens of rats (antinociception), and intra-accumbens injection of a CGRP antagonist induces a decrease in hot plate latency (hyperalgesia).⁵³ Furthermore, activation of cyclic adenosine monophosphate response element-binding protein in the nucleus accumbens modulates the hot plate response.⁵⁴ Thus, striatal tissue is potentially a relevant pain tissue⁵⁵ as well as mesolimbic dopamine signaling within this brain region⁵⁶ to understanding gene expression as it relates to genetic variation and nociception. It should also be noted that eQTLs from brain tissue are frequently expressed in multiple brain regions.⁵⁷ Thus, most of the genes we have identified will likely possess eQTLs in other CNS tissues that are relevant to thermal nociception.

Ryr1 showed the strongest eQTL association, it possesses multiple polymorphisms and is involved in calcium signaling in pain and inflammation. Ryr1 codes for the ryanodine 1 receptor, a calcium release channel in the sarcoplasmic reticulum that is also associated with the dihydropyridine receptor, or that is, L-type calcium channels. A large literature documents the involvement of L-type calcium channels in pain. Administration in male Swiss albino mice of ryanodine (i.c.v.), a Ryr antagonist dose dependently reduced hot plate latencies, whereas administration of 4-Cmc, a Ryr agonist, dose dependently increased hot plate latencies,⁵⁸ providing direct pharmacological evidence for an involvement of central, supraspinal Ryr receptors and calcium release in acute thermal nociception. The reduced hot plate latency in response to the Ryr agonist is in line with the reduced hot plate latency in response that is associated with both reduced Ryr1 expression and with the B6J allele (Table 1). Furthermore, administration of ryanodine and caffeine (Ryr agonist) in thalamocortical nuclei of B6 × 129 F1 mice produced similar bidirectional effects on both formalin-induced nociceptive behaviors during the late phase and acetic acid-induced writhing,⁵⁹ suggesting the possibility that differential Ryr1 expression could exert pleiotropic effects on both thermal and inflammatory nociception. On the other hand, another study showed that administration of verapamil (L-type blocker) in the dorsal horn in Sprague–Dawley rats had no effect on formalin-induced nociception.⁶⁰

RYR1-related myopathies are the most common class of congenital myopathies, and genetic variants within RYR1 within various coding regions of the gene have been associated with these disorders.⁶¹ Accordingly, in addition to Ryr1 polymorphisms potentially affecting nociception, Ryr1 variants could lead to a change in Ryr1 expression or dysfunction within the sarcoplasmic reticulum membrane of skeletal muscle, resulting in altered Ca2+ release and ultimately, myopathy. Thus, Ryr1 polymorphism(s) affecting musculoskeletal strength or function rather than a change in nociceptive transmission could underlie B6 substrain differences in the latency to behaviorally respond in the hot plate assay. Although we are unaware of any B6 substrain differences in myopathy, we and others have identified sensorimotor differences among B6 substrains in the rotarod test^21,22,29 that could reflect changes in musculoskeletal function.

In considering the fold-change in gene expression, it is striking that two cytochrome P450 genes, Cyp2a5 (26.84 Mb) and Cyp2g1 (26.81 Mb), are located squarely on top of the QTL peak for hot plate sensitivity and showed an extraordinary 3.8-fold and 3.6-fold change in transcription at the peak-associated hot plate marker, rs3148686 (Table 1). Notably, a nearby 112 kb structural variant (chromosome 7: 26.85–26.96 Mb) that is documented as an insertion with the B6J allele²⁹ and as a deletion of the B6NJ allele^15,16,29 (http://www.sanger.ac.uk/) is located just distally to Cyp2g1 (26.81 Mb) and Cyp2a5 (26.84 Mb). This structural variant includes only one protein-coding gene, Cyp2a22 but is flanked by a cluster of CYP genes, including not only Cyp2a5 and Cyp2g1 (Table 1) but also Cyp2b23, Cyp2b19, Cyp2a12, Cyp2f2, and Cyp2t4. The B6NJ deletion (or i.e., the B6J insertion²⁹) could contain DNA sequences that normally serve to regulate Cyp2a5 and Cyp2g1 transcription or it could induce compensatory increases in Cyp2a5 and Cyp2g1 transcription through a separate mechanism.

To summarize, we identified robust B6 substrain differences in pain sensitivity in an additional nociceptive modality whereby the B6J strain was more sensitive to formalin-induced nociceptive behavior. These results support the previous conclusion by Mogil and coworkers decades ago that the B6J strain is not the ideal strain to study the molecular basis of nociceptive phenotypes as it is hypersensitive to multiple nociceptive modalities,^4,62 yet a majority of molecular genetic studies continue to only utilize mice on a B6J background. A recent comprehensive study of F1 crosses between 30 inbred strains highlights the profound impact of epistatic modifier alleles in determining the magnitude and direction of knockout alleles on complex behavioral traits and the importance of studying mutations on multiple genetic backgrounds.⁶³ A limitation of the current study is that we only examined behaviors that model the sensory aspect of pain. An important question is whether or not there are B6 substrain differences in behaviors that model the motivational-affective and cognitive-evaluative components of pain. Furthermore, additional genetic crosses between B6 substrains will determine the degree to which these different models of pain are genetically shared versus dissociable and rapid fine mapping strategies²⁰ combined with gene editing will eventually identify the causal genetic factors. Previous successes in translating “pain genes” from mice to humans suggest that identification of these genetic factors will be translationally relevant.

Supplemental Material

Supplemental material for C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception

Supplemental Material for C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception by Camron D Bryant, Deniz Bagdas, Lisa R Goldberg, Tala Khalefa, Eric R Reed, Stacey L Kirkpatrick, Julia C Kelliher, Melanie M Chen, William E Johnson, Megan K Mulligan and M Imad Damaj in Molecular Pain

Footnotes

Acknowledgment

The authors thank Dr Pieter Faber, Technical Director at the University of Chicago Genomics Facility, for his assistance.

Author contributions

CDB contributed to the experimental design, supervised execution and analysis of the genetic study, and wrote the manuscript; DB conducted experiments and analyzed data; TK conducted experiments; LRG conducted experiments and analyzed data; ERR conducted computational analysis of gene expression; SLK conducted experiments and analyzed data; JCK and MMC analyzed data and prepared data for presentation; WEJ supervised computational analysis of gene expression; MKM contributed to writing the manuscript; MID contributed the experimental design and supervised the parental strain studies.

Declaration of Conflicting Interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was funded by R01CA221260 (MID and CDB), R01DA039168 (CDB), R21DA038738 (CDB), and T32GM008541 (LRG).

References

Dib-Hajj

Waxman

SG.

Translational pain research: lessons from genetics and genomics. Sci Transl Med 2014; 6: 249sr4.

Mogil

JS.

Pain genetics: past, present and future. Trends Genet 2012; 28: 258–266.

Coghill

RC.

Individual differences in the subjective experience of pain: new insights into mechanisms and models. Headache 2010; 50: 1531–1535.

Mogil

Wilson

Bon

Lee

Chung

Raber

Pieper

Hain

Belknap

Hubert

Elmer

Chung

Devor

Heritability of nociception I: responses of 11 inbred mouse strains on 12 measures of nociception. Pain 1999; 80: 67–82.

Lariviere

Wilson

Laughlin

Kokayeff

West

Adhikari

Wan

Mogil

JS.

Heritability of nociception. III. genetic relationships among commonly used assays of nociception and hypersensitivity. Pain 2002; 97: 75–86.

Mogil

Wilson

Bon

Lee

Chung

Raber

Pieper

Hain

Belknap

Hubert

Elmer

Chung

Devor

Heritability of nociception II. ‘types' of nociception revealed by genetic correlation analysis. Pain 1999; 80: 83–93.

Young

Lariviere

Belfer

Genetic basis of pain variability: recent advances. J Med Genet 2012; 49: 1–9.

Veluchamy

Hebert

Meng

Palmer

CNA

Smith

BH.

Systematic review and meta-analysis of genetic risk factors for neuropathic pain. Pain 2018; 159: 825–848.

Baud

Flint

Identifying genes for neurobehavioural traits in rodents: progress and pitfalls. Dis Model Mech 2017 ; 10: 373–383.

10.

Nissenbaum

Devor

Seltzer

Gebauer

Michaelis

Tal

Dorfman

Abitbul-Yarkoni

Elahipanah

delCanho

Minert

Fried

Persson

Shpigler

Shabo

Yakir

Pisante

Darvasi

Susceptibility to chronic pain following nerve injury is genetically affected by CACNG2. Genome Res 2010; 20: 1180–1190.

11.

Sorge

Trang

Dorfman

Smith

Beggs

Ritchie

Austin

Zaykin

Vander Meulen

Costigan

Herbert

Yarkoni-Abitbul

Tichauer

Livneh

Gershon

Zheng

Tan

John

Slade

Jordan

Woolf

Peltz

Maixner

Diatchenko

Seltzer

Salter

Mogil

JS.

Genetically determined P2X7 receptor pore formation regulates variability in chronic pain sensitivity. Nat Med 2012; 18: 595–599.

12.

Mogil

Miermeister

Seifert

Strasburg

Zimmermann

Reinold

Austin

Bernardini

Chesler

Hofmann

Hordo

Messlinger

Nemmani

Rankin

Ritchie

Siegling

Smith

Sotocinal

Vater

Lehto

Klussmann

Quirion

Michaelis

Devor

Reeh

PW.

Variable sensitivity to noxious heat is mediated by differential expression of the CGRP gene. Proc Natl Acad Sci USA 2005; 102: 12938–12943.

13.

Mogil

Wilson

Chesler

Rankin

Nemmani

Lariviere

Groce

Wallace

Kaplan

Staud

Ness

Glover

Stankova

Mayorov

Hruby

Grisel

Fillingim

RB.

The melanocortin-1 receptor gene mediates female-specific mechanisms of analgesia in mice and humans. Proc Natl Acad Sci USA 2003; 100: 4867–4872.

14.

Mogil

Sorge

LaCroix-Fralish

Smith

Fortin

Sotocinal

Ritchie

Austin

Schorscher-Petcu

Melmed

Czerminski

Bittong

Mokris

Neubert

Campbell

Edwards

Campbell

Crawley

Lariviere

Wallace

Sternberg

Balaban

Belfer

Fillingim

RB.

Pain sensitivity and vasopressin analgesia are mediated by a gene-sex-environment interaction. Nat Neurosci 2011; 14: 1569–1573.

15.

Keane

Goodstadt

Danecek

White

Wong

Yalcin

Heger

Agam

Slater

Goodson

Furlotte

Eskin

Nellaker

Whitley

Cleak

Janowitz

Hernandez-Pliego

Edwards

Belgard

Oliver

McIntyre

Bhomra

Nicod

Gan

Yuan

van der Weyden

Steward

Bala

Stalker

Mott

Durbin

Jackson

Czechanski

Guerra-Assuncao

Donahue

Reinholdt

Payseur

Ponting

Birney

Flint

Adams

DJ.

Mouse genomic variation and its effect on phenotypes and gene regulation. Nature 2011; 477: 289–294.

16.

Yalcin

Wong

Agam

Goodson

Keane

Gan

Nellaker

Goodstadt

Nicod

Bhomra

Hernandez-Pliego

Whitley

Cleak

Dutton

Janowitz

Mott

Adams

Flint

Sequence-based characterization of structural variation in the mouse genome.

Nature 2011; 477: 326–329.

17.

Bryant

CD.

The blessings and curses of C57BL/6 substrains in mouse genetic studies. Ann N Y Acad Sci 2011; 1245: 31–33.

18.

Kumar

Kim

Joseph

Kourrich

Yoo

Huang

Vitaterna

de Villena

Churchill

Bonci

Takahashi

JS.

C57BL/6N mutation in cytoplasmic FMRP interacting protein 2 regulates cocaine response. Science 2013; 342: 1508–1512.

19.

Kirkpatrick

Goldberg

Yazdani

Babbs

Reed

Jenkins

Bolgioni

Landaverde

Luttik

Mitchell

Kumar

Johnson

Mulligan

Cottone

Bryant

CD.

Cytoplasmic FMR1-interacting protein 2 is a major genetic factor underlying binge eating. Biol Psychiatry 2017; 81: 757–769.

20.

Bryant

Ferris

De Villena

FPM

Damaj

Kumar

Mulligan

MK.

Reduced complexity cross design for behavioral genetics. In: Gerlai

(ed) Molecular-genetic and statistical techniques for behavioral and neural research. 2018, pp.165–190.

21.

Bryant

Zhang

Sokoloff

Fanselow

Ennes

Palmer

McRoberts

JA.

Behavioral differences among C57BL/6 substrains: Implications for transgenic and knockout studies. J Neurogenet 2008; 22: 315–331.

22.

Matsuo

Takao

Nakanishi

Yamasaki

Tanda

Miyakawa

Behavioral profiles of three C57BL/6 substrains. Front Behav Neurosci 2010; 4: 29.

23.

Cooper

O'Meara

Jack

Elliot

Lamb

Khan

Menta

Ryals

Winter

Wright

DE.

Intrinsic activity of C57BL/6 substrains associates with high-fat diet-induced mechanical sensitivity in mice. J Pain 2018; 19: 1285–1295.

24.

Chaplan

Pogrel

Yaksh

TL.

Role of voltage-dependent calcium channel subtypes in experimental tactile allodynia. J Pharmacol Exp Ther 1994; 269: 1117–1123.

25.

Bagdas

AlSharari

Freitas

Tracy

Damaj

MI.

The role of alpha5 nicotinic acetylcholine receptors in mouse models of chronic inflammatory and neuropathic pain. Biochem Pharmacol 2015; 97: 590–600.

26.

Bryant

Eitan

Sinchak

Fanselow

Evans

CJ.

NMDA receptor antagonism disrupts the development of morphine analgesic tolerance in male, but not female C57BL/6J mice. Am J Physiol Regul Integr Comp Physiol 2006; 291: R315–R326.

27.

Eddy

Leimbach

Synthetic analgesics. II. dithienylbutenyl- and dithienylbutylamines. J Pharmacol Exp Ther 1953; 107: 385–393.

28.

Kirkpatrick

Bryant

CD.

Behavioral architecture of opioid reward and aversion in C57BL/6 substrains. Front Behav Neurosci 2014; 8: 450.

29.

Simon

Greenaway

White

Fuchs

Gailus-Durner

Sorg

Wong

Bedu

Cartwright

Dacquin

Djebali

Estabel

Graw

Ingham

Jackson

Lengeling

Mandillo

Marvel

Meziane

Preitner

Puk

Roux

Adams

Atkins

Ayadi

Becker

Blake

Brooker

Cater

Champy

Combe

Danecek

di Fenza

Gates

Gerdin

Golini

Hancock

Hans

Holter

Hough

Jurdic

Keane

Morgan

Muller

Neff

Nicholson

Pasche

Roberson

Rozman

Sanderson

Santos

Selloum

Shannon

Southwell

Tocchini-Valentini

Vancollie

Wells

Westerberg

Wurst

Yalcin

Ramirez-Solis

Steel

Mallon

Hrab

de Angelis

Herault

Brown

SD.

A comparative phenotypic and genomic analysis of C57BL/6J and C57BL/6N mouse strains. Genome Biol 2013; 14: R82.

30.

Yazdani

Parker

Shen

Reed

Guido

Kole

Kirkpatrick

Lim

Sokoloff

Cheng

Johnson

Palmer

Bryant

CD.

Hnrnph1 is A quantitative trait gene for methamphetamine sensitivity. PLoS Genet 2015; 11: e1005713.

31.

Broman

Sen

Churchill

GA.

R/qtl: QTL mapping in experimental crosses. Bioinformatics (Oxford, England) 2003; 19: 889–890.

32.

Bryant

Chang

Zhang

Wiltshire

Tarantino

Palmer

AA.

A major QTL on chromosome 11 influences psychostimulant and opioid sensitivity in mice. Genes Brain Behav 2009; 8: 795–805.

33.

Cox

Ackert-Bicknell

Dumont

Ding

Bell

Brockmann

Wergedal

Bult

Paigen

Flint

Tsaih

Churchill

Broman

KW.

A new standard genetic map for the laboratory mouse. Genetics 2009; 182: 1335–1344.

34.

Trapnell

Roberts

Goff

Pertea

Kim

Kelley

Pimentel

Salzberg

Rinn

Pachter

Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc 2012; 7: 562–578.

35.

Ritchie

Phipson

Law

Shi

Smyth

GK.

Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res 2015; 43: e47.

36.

Law

Chen

Shi

Smyth

GK.

Voom: precision weights unlock linear model analysis tools for RNA-seq read counts.

Genome Biol 2014; 15: R29.

37.

Benjamini

Hochberg

Controlling false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc 1995; 57: 289–300.

38.

Mekada

Abe

Murakami

Nakamura

Nakata

Moriwaki

Obata

Yoshiki

Genetic differences among C57BL/6 substrains. Exp Anim 2009; 58: 141–149.

39.

Pandey

Williams

RW.

Genetics of gene expression in CNS. Int Rev Neurobiol 2014; 116: 195–231.

40.

Tjolsen

Berge

Hunskaar

Rosland

Hole

The formalin test: an evaluation of the method. Pain 1992; 51: 5–17.

41.

Mannino

South

Quinones-Jenab

Inturrisi

CE.

Estradiol replacement in ovariectomized rats is antihyperalgesic in the formalin test. J Pain 2007; 8: 334–342.

42.

Ralya

McCarson

KE.

Acute estrogen surge enhances inflammatory nociception without altering spinal fos expression. Neurosci Lett 2014; 575: 91–95.

43.

Wang

Zhang

YQ.

Ovariectomy results in variable changes in nociception, mood and depression in adult female rats. PLoS One 2014; 9: e94312.

44.

Kuba

Nazarian

Festa

Barr

Jenab

Inturrisi

Quinones-Jenab

Estradiol and progesterone differentially regulate formalin-induced nociception in ovariectomized female rats. Horm Behav 2006; 49: 441–449.

45.

Ibironke

Aji

KE.

Pain threshold variations in female rats as a function of the estrus cycle. Niger J Physiol Sci 2011; 26: 67–70.

46.

Frye

Cuevas

Kanarek

RB.

Diet and estrous cycle influence pain sensitivity in rats. Pharmacol Biochem Behav 1993; 45: 255–260.

47.

Kuba

Quinones-Jenab

The role of female gonadal hormones in behavioral sex differences in persistent and chronic pain: clinical versus preclinical studies. Brain Res Bull 2005; 66: 179–188.

48.

Sorge

Martin

Isbester

Sotocinal

Rosen

Tuttle

Wieskopf

Acland

Dokova

Kadoura

Leger

Mapplebeck

McPhail

Delaney

Wigerblad

Schumann

Quinn

Frasnelli

Svensson

Sternberg

Mogil

JS.

Olfactory exposure to males, including men, causes stress and related analgesia in rodents. Nat Methods 2014; 11: 629–632.

49.

Smith

Marker

Perry

Liao

Sotocinal

Austin

Melmed

Clark

Peltz

Wickman

Mogil

JS.

Quantitative trait locus and computational mapping identifies Kcnj9 (GIRK3) as a candidate gene affecting analgesia from multiple drug classes. Pharmacogenet Genomics 2008; 18: 231–241.

50.

Simpson

Linder

Sargent

Davisson

Mobraaten

Sharp

JJ.

Genetic variation among 129 substrains and its importance for targeted mutagenesis in mice. Nat Genet 1997; 16: 19–27.

51.

Furuse

Miura

Yagasaki

Shiroishi

Koide

Identification of QTLs for differential capsaicin sensitivity between mouse strains KJR and C57BL/6. Pain 2003; 105: 169–175.

52.

Koide

Moriwaki

Ikeda

Niki

Shiroishi

Multi-phenotype behavioral characterization of inbred strains derived from wild stocks of mus musculus. Mamm Genome 2000; 11: 664–670.

53.

Lundeberg

LC.

Involvement of CGRP and CGRP1 receptor in nociception in the nucleus accumbens of rats. Brain Res 2001; 901: 161–166.

54.

Barrot

Olivier

Perrotti

DiLeone

Berton

Eisch

Impey

Storm

Neve

Yin

Zachariou

Nestler

EJ.

CREB activity in the nucleus accumbens shell controls gating of behavioral responses to emotional stimuli. Proc Natl Acad Sci USA 2002; 99: 11435–11440.

55.

Barcelo

Filippini

Pazo

JH.

The striatum and pain modulation. Cell Mol Neurobiol 2012; 32: 1–12.

56.

Taylor

Becker

Schweinhardt

Cahill

Mesolimbic dopamine signaling in acute and chronic pain: implications for motivation, analgesia, and addiction. Pain 2016; 157: 1194–1198.

57.

GTEx Consortium, Laboratory, Data Analysis &Coordinating Center (LDACC)-Analysis Working Group, Statistical Methods groups-Analysis Working Group, Enhancing GTEx (eGTEx) groups, NIH Common Fund, NIH/NCI, NIH/NHGRI, NIH/NIMH, NIH/NIDA, Biospecimen Collection Source Site-NDRI, Biospecimen Collection Source Site-RPCI, Biospecimen Core Resource-VARI, Brain Bank Repository-University of Miami Brain Endowment Bank, Leidos Biomedical-Project Management, ELSI Study, Genome Browser Data Integration &Visualization-EBI, Genome Browser Data Integration &Visualization-UCSC Genomics Institute, University of California Santa Cruz, Lead analysts:, Laboratory, Data Analysis &Coordinating Center (LDACC):, NIH program management:, Biospecimen collection:, Pathology:, eQTL manuscript working group:, Battle A, Brown CD, Engelhardt BE, Montgomery SB. Genetic effects on gene expression across human tissues. Nature 2017; 550: 7675204–7675213.

58.

Galeotti

Bartolini

Ghelardini

Role of intracellular calcium in acute thermal pain perception. Neuropharmacology 2004; 47: 935–944.

59.

Cheong

Kim

Choi

Sun

Shin

HS.

Thalamic ryanodine receptors are involved in controlling the tonic firing of thalamocortical neurons and inflammatory pain signal processing. J Neurosci 2011; 31: 1213–1218.

60.

Diaz

Dickenson

AH.

Blockade of spinal N- and P-type, but not L-type, calcium channels inhibits the excitability of rat dorsal horn neurones produced by subcutaneous formalin inflammation. Pain 1997; 69: 93–100.

61.

Lawal

Todd

Meilleur

KG.

Ryanodine receptor 1-related myopathies: diagnostic and therapeutic approaches.

Neurotherapeutics 2018; 15: 885–899.

62.

Mogil

Wilson

SG.

Nociceptive and morphine antinociceptive sensitivity of 129 and C57BL/6 inbred mouse strains: implications for transgenic knock-out studies. Eur J Pain 1997; 1: 293–297.

63.

Sittig

Carbonetto

Engel

Krauss

Barrios-Camacho

Palmer

AA.

Genetic background limits generalizability of genotype-phenotype relationships. Neuron 2016; 91: 1253–1259.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

0.51 MB