Skip to main content

MicroRNA circulating in the early aftermath of motor vehicle collision predict persistent pain development and suggest a role for microRNA in sex-specific pain differences

Abstract

Background

Molecular mediators influencing the transition from acute to persistent musculoskeletal pain following common stress exposures such as motor vehicle collision (MVC) remain poorly understood. In this exploratory, proof of concept study, we compared circulating microRNA (miRNA) expression profiles in the early aftermath of MVC among individuals who did and did not subsequently develop persistent pain. Blood RNA samples were obtained from African American individuals (n = 53) who presented to the emergency department after MVC and were discharged to home after evaluation. The presence or absence of severe pain in the axial region, the most common and morbid region in which post-MVC pain occurs, was assessed 6 weeks following MVC via standardized questionnaire. miRNA expression was determined using miRNA-sequencing; nonparametric analyses were used to compare miRNA expression levels among individuals with and without persistent pain.

Results

Thirty-two mature miRNA were differentially expressed (p < 0.05) in those with and without severe axial pain at 6 weeks. miR-135a-5p, a regulator of the serotonin receptor that is known to be stress-responsive, differed most significantly between groups (p = 3 × 10−4). This miRNA, and miR-3613-3p (p = 0.001) survived correction for multiple testing (FDR = 0.15) in this small sample. Interestingly, differentially expressed miRNA were enriched for X chromosome location. In secondary analyses, the eight X chromosome miRNA were (a) more significantly associated with axial pain in women than men, (b) expressed more highly in the peripheral blood of women than men, and (c) predicted in pathway analyses (DIANA miRPath v 2.0) to regulate neuronal and neuroendocrine pathways previously implicated in various pain pathologies.

Conclusions

These results show that circulating miRNA predict persistent severe axial pain after MVC and suggest that they may be involved in the pathogenesis of post-traumatic musculoskeletal pain. However, further studies are needed to determine if these miRNA play a direct causal role.

Background

MicroRNA (miRNA) are small non-coding RNA molecules that regulate gene expression by binding target mRNA. During the past decade, the study of miRNA has transformed understanding of the regulation of major biological pathways [1, 2] and advanced understanding of the pathogenesis of a number of common diseases (e.g. [35]). Substantial evidence suggests that miRNA may play a critical role in molecular pathways underpinning diverse pain conditions (e.g. [68]). Most of this evidence comes from animal studies; further human studies are needed which prospectively evaluate the potential role of miRNA in persistent pain development after potential triggering events.

One common potential triggering event for persistent pain is exposure to a motor vehicle collision (MVC). More than fifty million MVCs occur each year worldwide [9], and more than 4 million individuals present to US emergency departments (ED) each year for care after MVC [10, 11]. Ninety percent of individuals presenting to US EDs for care after MVC are discharged home after evaluation with little or no identifiable tissue injury [12]. A substantial proportion of these individuals develop musculoskeletal pain, most commonly in the axial region (neck, shoulders, and/or back) [13].

The molecular mechanisms responsible for axial pain (AP) development after MVC remain poorly understood. This lack of understanding is a major barrier to the development of more effective preventive interventions. If miRNA play an important role in the pathogenesis of post-MVC AP, then studies identifying differences in miRNA populations among those who do and do not subsequently develop these outcomes may advance understanding of post-MVC AP pathogenesis [1416]. Such studies must be feasible and ethical, and should obtain miRNA samples from tissue relevant to disease pathogenesis and/or be representative of such tissue.

Several lines of evidence suggest that blood is not only a feasible source of miRNA for studies of post-MVC musculoskeletal pain pathogenesis, but also that blood miRNA studies may provide pathogenic insights [1719]. First, stress systems appear to be involved in the pathogenesis of persistent pain after MVC [2023], and blood-borne stress and immune-related factors are an important component of the systemic stress response [24]. In addition, RNA expression patterns in blood, CNS, and endocrine tissues are strongly correlated [15, 25, 26], and a number of previous studies have linked miRNA expression in blood with neurologic disease outcomes [2730]. Finally, despite challenges related to tissue-specific gene expression, the study of blood miRNA expression has led to important new understanding of diverse diseases, including pain conditions [27, 31, 32].

In this prospective study, we compared ED blood miRNA profiles among individuals who presented to the ED for evaluation after MVC. We hypothesized that ED blood miRNA profiles would differ among those who did and did not have severe persistent AP 6 weeks later.

Results

Cohort

Characteristics of the study sample (n = 53) are shown in Table 1. Samples were drawn from a large prospective cohort study of African Americans (R01AR060852); all participants were African American (AA) and nearly six in ten were female. Most were less than 40 years old, had some college education, made less than 40 K annually, and were overweight (average BMI = 30). All individuals in this study presented to the ED within 6 h of MVC and most arrived within 1 h. Only individuals who were discharged from the ED who reported no lacerations, avulsions, or major tissue injury were included. Additionally, participants all had a severity score of 1 on the Abbreviated Injury Scale (AIS) [33], indicating minimal anatomical injury. Six weeks following MVC, severe AP was present in 27/53 participants. Among these 27 individuals who developed severe AP, 16 (59 %) were women.

Table 1 Study characteristics

miRNA sequencing quality assessment

An average of 9 million sequencing reads were obtained per participant from blood samples obtained in the ED in the early aftermath of MVC. More than 95 % of these miRNA aligned with miRNA in miRBase, indicating that the majority of the sequencing reads were mature miRNAs (vs. degradation products, linker–linker contaminants, etc.). Mature miRNA with an average of ≥300 sequencing reads across all 53 samples (n = 376 miRNA) were included in analyses. Relative proportions of several miRNA typically found in peripheral blood were very similar to those reported previously [15, 34] (data not shown).

Evaluation of ED miRNA expression levels among those who did and did not develop persistent AP 6 weeks following MVC

Thirty-two of 376 (9 %) miRNA detected in ED blood samples were differentially expressed at the p <0.05 level among those who did and did not report severe post-MVC AP at 6 weeks, with fold differences ranging from −3.71 to 2.49 (Table 2, sequencing read counts used for determination of mean fold differences for each miRNA are included in Additional file 1: Table S1). Nine of these 32 differentially expressed miRNA have previously been associated with pain and/or stress system physiology in neurological tissue and/or blood (see ‘Ref’ column, Table 2). Two of the 32 differentially expressed miRNA, miR-135a-5p (p = 3 × 10−4) and miR-3613-3p (p = 0.001), met our pre-hoc significance level threshold for multiple testing of 0.15, corresponding to a p value cut-off of <0.003 [35].

Table 2 microRNA in whole blood circulating in the early aftermath of motor vehicle collision in African Americans that predict axial pain development 6 weeks after MVC trauma

Validation of miRNA sequencing results using RT-qPCR

Technical and qualitative validation of miRNA sequencing results was performed on a random subsample of 7 of the 32 differentially expressed miRNA using reverse transcription quantitative-PCR (RT-qPCR) [36]. In each case, concordance between miRNA sequencing results and RT-qPCR results was observed for direction of differential expression (i.e., positive or negative expression difference, Fig. 1). Magnitude of direction of effect was also generally similar. The Spearman Correlation between the two methods was also calculated: r = 0.786 p = 0.036.

Fig. 1
figure 1

RT-qPCR validation of microRNA that predict axial pain (AP) development following motor vehicle collision (MVC). Dark grey bars represent expression differences calculated using mean microRNA sequencing counts in individuals developing AP divided by mean microRNA sequencing counts in individuals recovering after MVC. Light grey bars represent expression differences in the same two groups using mean cycle thresholds generated via RT-qPCR. Spearman Correlation assessing correlation of the two methods: r = 0.786, p = 0.036

Differentially expressed miRNA were enriched for X chromosome location

Specific data regarding chromosomal origin and strand (sense or antisense) for the 32 differentially expressed miRNA is shown in Additional file 1: Table S1. The 32 miRNA predictive of severe post-MVC AP at 6 weeks were enriched for gene location on the X chromosome (Fig. 2) in comparison to all X chromosome miRNA identified in the sample (8/32 (25 %) vs. 28/376 (8 %), p = 0.038). This holds true despite the fact that approximately 8 % of the detected blood miRNA originate from the X chromosome.

Fig. 2
figure 2

microRNA that predict axial pain (AP) development are transcribed from 14 different chromosomes. The percentage of miRNAs (out of the 32 associated with AP, Table 2) that originate from each of the 14 chromosomes are represented by pieces of the pie chart based on shading (lightest shading = 3 % and darkest shading = 25 %, with intermediate percentages having intermediate shading colors, see legend). The names/number of each chromosome are labeled inside each piece of the pie

X chromosome miRNA identified in this study are more highly associated with the development of persistent AP following MVC in women than men

Because X chromosome gene expression can be sex-dependent [37], after we observed the X chromosome enrichment of differentially-expressed miRNA in our sample, we assessed for interactions between sex and the effect on severe persistent AP of differentially expressed X chromosome miRNA. Even in the relatively small samples of men and women assessed, sex × miRNA interactions were significant at the p <0.05 level for miR-362-5p and let-7f-2-3p, and were present at the trend level for miR-501-5p and miR-500b-5p (Table 3, Additional file 1: Table S1). Much greater differences in expression according to pain outcome were observed in women (n = 31) vs. men (n = 22) (Table 3). Seven of the 8 X chromosome miRNA were significantly associated with AP development in women, whereas none of the X chromosome miRNA were associated with AP development in men (Table 3). Additionally, after discovering significant sex × miRNA interactions for the X-chromosome miRNA described above, we assessed for additional interactions between participant sex and miRNA associated with persistent AP. A significant sex × miRNA interaction was present for two additional miRNA, miR-1296 (p = 0.017) and Let-7a-3p (p = 0.022) (see Additional file 1: Table S1 for results of all sex × miRNA interactions). These two miRNA are significantly associated with persistent severe AP in women (p < 0.001 and 0.001) but not in men (Additional file 2: Table S2).

Table 3 Association of Emergency Department expression levels of microRNA (miRNA) from the X chromosome with persistent axial pain following motor vehicle collision in women vs. men and assessment of interaction between sex of an individual and miRNA

X chromosome miRNA identified in this study are expressed more highly in the blood in the early aftermath of MVC in women than men

Based on previous reports showing higher expression of some X chromosome genes in women than men (most notably in brain tissue [38] and possibly due to mechanisms such as escape from X chromosome inactivation [39]), we assessed whether the X chromosome miRNA identified in this study are expressed at higher levels in women than men. All eight miRNA were expressed at higher levels in women than in men developing severe AP, although only one was statistically significant (miR-502-3p, p = 0.017) (Table 4).

Table 4 Emergency Department expression level differences of microRNA (miRNA) from the X chromosome in women vs. men developing persistent axial pain following motor vehicle collision

Evaluation of biologic pathways targeted by AP-associated X chromosome miRNA

Using DIANA miRPath v 2.0 [40], we assessed for molecular pathways (KEGG pathways [41]) overrepresented in predicted targeting by these 8 X chromosome miRNA. The pathways with the highest number of gene transcripts targeted by the eight X chromosome miRNA in Table 2 (i.e. most statistically significant enrichment) out of ~450 KEGG pathways that DIANA miRPath queries are shown in Table 5 (false discovery rates (FDRs) calculated via permutation testing [42]). These pathways include neuronal and neuroendocrine pathways such as the Long Term Potentiation pathway (p = 1.67 × 10−9, FDR ≤ 0.08), the Axon Guidance pathway (p = 1.32 × 10−7, FDR ≤ 0.22), Neurotrophin signaling (p = 7.35 × 10−7, FDR ≤ 0.22) and the Dopaminergic synapse signaling pathway (p = 1.56 × 10−6, FDR ≤ 0.33). Substantial evidence supports an important role for these pathways in mediating physiologic responses to stress and the pathogenesis of acute and persistent pain [4345]. In addition to assessing X-linked miRNA targeted pathways, we also assessed which pathways might be enriched in targeting by all 32 miRNA associated with persistent AP development. This data is provided in Additional file 3: Table S3.

Table 5 DIANA miRPath predicted KEGG pathways enriched in targeting by X chromosome miRNA differentially regulated in the early aftermath of MVC trauma in AA individuals who develop AP following MVC vs. those who recover

Discussion

Persistent pain is a common and poorly understood sequela of traumatic/stressful events such as MVC [46]. The results of this study show that even in the relatively small study sample, circulating blood miRNA in the first hours after stress exposure differed significantly among AAs who did and did not have severe persistent MVC-related AP 6 weeks later. Study results also provide an example of the novel pathophysiologic insights that may be obtained, as enrichment of X chromosome miRNA among differentially expressed miRNA, and sex differences in the effect of these miRNA, suggest potential mechanisms contributing to sex differences in vulnerability to persistent post-MVC pain. Supporting the role of these differentially expressed miRNA in persistent pain pathogenesis, these miRNA included a number of miRNA previously associated with the stress response and/or pain processing. However, further studies are needed to determine whether these miRNA play a causal role in the pathogenesis of pain post-MVC. Of note, study findings also support the legitimacy of persistent post-MVC musculoskeletal pain as a “real” disease outcome. This is important, as patients with post-MVC pain outcomes are highly stigmatized [47]. In addition, study findings also contribute to a growing literature indicating that stress-mediated changes in neurosensory processing play an important role in the pathogenesis of post-MVC outcomes.

We do not know whether miRNA predicting severe AP in the present study play a causal role in the development of these outcomes, or are markers for other cellular processes directly involved. If these miRNA play a causal role, mechanisms by which miRNA detected in blood may directly influence AP outcomes include: (1) miRNA expressed outside the central nervous system (CNS) and detected in the blood may cross the blood brain barrier [48, 49] to influence CNS processes, (2) miRNA expressed in the CNS may alter CNS transcription and also be released into the periphery (e.g., as part of a cellular/systemic communication system [50, 51]), (3) miRNA expressed outside the CNS may influence extra-CNS processes involved in the pathogenesis of post-traumatic pain [52, 53] (animal model data suggest that systemic, extra-CNS processes may play an important role in the pathogenesis of stress-induced pain [54]). Future studies are needed to better understand the identified associations between miRNA circulating in the immediate aftermath of trauma and the development of persistent pain states. Such studies may provide new insights into the biology of chronic pain development.

A potential causal role of the miRNA identified in the present study in the pathogenesis of severe post-MVC axial pain is supported by the known role of several of these miRNA in pain and/or stress-related processes. For example, miR-135a-5p binds the mineralcorticoid receptor (NR3C2) [55], the serotonin transporter (SLC6A4) [56], and the serotonin receptor-1a (HTR1A) [56] transcripts, all of which can affect pain processing (e.g. [5759]), and miR-135a-5p has been shown to be expressed in pain-relevant tissues including the amygdala [60], 5-HT neurons [56], spinal cord [61], and pre-frontal cortex [62]. miR-19b-3p has been shown to be stress responsive in both the amygdala [63] and blood leukocytes [64], and can regulate the adrenergic receptor β-1 (ADRB1) [63]. Other miRNA, such as miR-3613-3p, have not been studied extensively; target prediction algorithms (e.g., TargetScan v7.0) predict that miR-3613-3p modulates the expression of pain-associated genes including GABRB3, GRIN3A, TRPV1, NPY1R, and SCN9A. Further examples of miRNA identified in this study that have prior associations with pain include miR-185-5p, miR-7-5p, and miR-26b-3p [61, 62, 65]. Equally important, many of the miRNA identified in this study are not currently known to be associated with stress or pain-related outcomes, suggesting that investigations such as the present study have the potential to lead to the identification of novel miRNA mediators.

An unexpected finding in the present study was that differentially expressed miRNA were enriched for miRNA located on the X chromosome. This finding holds true even accounting for the relatively high abundance of X chromosome miRNA expressed in the blood (8 %) compared to miRNA from other chromosomes. These X chromosomal miRNA were consistently expressed at higher levels in women than in men and appeared to contribute to persistent severe AP in women but not in men. These X chromosome miRNA were predicted to target pain-relevant transcripts from KEGG pathways known to be associated with various pain phenotypes, such as long term potentiation, neurotrophin signaling, and dopaminergic signaling [4345]. Pathophysiologic mechanisms by which these miRNA may contribute to severe post-MVC AP in women but not in men are currently unknown. Six of the eight X chromosome miRNA identified to be associated with persistent severe AP in this study are transcribed as part of the miR-532-502 cluster of miRNA (a cluster which includes the expression of miR-532, -188, -500, -362, -501, -500b, -660, and -502). Upon examination of DNase hypersensitivity regions upstream of this cluster and upstream of miR-20b-5p and let-7f-2-3p, we did not find any obvious binding regions for sex hormone responsive transcription factors (TFs). However, let-7f-2-3p has been shown experimentally to be induced by estradiol and is in a dosage sensitive region of the X chromosome [66]. In addition to sex hormone responsive TFs, another mechanism that may contribute to the observed sex differences in miRNA effect is X chromosome inactivation. This phenomenon is believed to influence the transcription of ~15 % of all X chromosome genes [67, 68], and the results of at least two other studies suggest that miRNA genes can escape X chromosome inactivation. One study found that X chromosome miRNA genes are over-expressed in the T cells of women with lupus [69], and another study found sex-biased miRNA expression in the neonatal brain [70]. Further studies are needed to assess for associations between X chromosome miRNA expression and pain outcomes in men and women experiencing MVC, and to evaluate potential physiologic mechanisms by which sex differences in the expression of these miRNA may occur.

Some limitations should be considered when interpreting the results of this study. First, the sample size of this initial proof-of-concept study was relatively small. Future studies with much larger samples and greater power are needed. These studies should also adjust more stringently for multiple comparisons to reduce the probability of Type I error. Second, we did not adjust for potential confounders such as participant age, sex, or BMI in our primary analyses. However, adjusting for these factors in exploratory analyses did not diminish our effect size estimates. Third, we were not able to adjust for the potential confounding influence of medications administered in the ED and miRNA expression. We were able to evaluate chronic medication use in our cohort: no individuals in the study were taking opioids (due to exclusion criteria), 4 % of individuals were taking acetaminophen, and 13 % were taking NSAID. Adjusting for pre-MVC acetaminophen or NSAID use had negligible effect on effect size estimates for the association between miR-135a-5p and miR-3613-3p and severe AP development. Similarly, we were also unable to adjust for any potential confounding effect due to pre-MVC chronic illnesses, as comprehensive past medical history data on study participants was not available. Data from the emergency department record, for which past medical history data is often incomplete, indicated that the most prevalent chronic illnesses in the cohort were hypertension (25 %) and asthma (13 %). Hypertension was not associated with the development of severe AP (p = 0.757), and adjusting for hypertension did not weaken the association between miR-135a-5p and miR-3613-3p and severe AP development. Similar results were obtained when assessing the potential effect of asthma. Fourth, miRNA expression differences between those who did and did not develop severe AP were evaluated in the ED, in the early aftermath of MVC, and we do not know how the expression of these miRNA changed over time. However, we found that even in the very early aftermath of MVC, miRNA expression differed in those who did and did not subsequently develop persistent pain. Fifth, our study population was limited to African Americans, an understudied group that has been shown to experience an increased burden of adverse pain outcomes after trauma [7174]. The generalizability of our findings to other ethnic groups is unknown. Sixth, larger sex specific strata are needed to fully understand miRNA expression differences associated with persistent severe AP development in women vs. men. Finally, pathway and gene target identification analyses were based on predicted binding rather than actual binding in biologic assays. However, predicted binding has been shown to have high concordance with actual binding, and predicted binding has the advantage of providing an unbiased assessment across the entire genome [75].

Conclusion

The results of this study show that two miRNA, miR-135a-5p and miR-3613-3p, predict persistent AP development after MVC. In addition, study results suggest that X chromosome miRNA contribute to persistent pain development after MVC stress exposure in women, and that such miRNA may contribute to sex differences in vulnerability to persistent pain after MVC. More broadly, the results of this study support the hypothesis that analyses of miRNA collected from blood in the early aftermath of trauma/stress exposure might provide new insights into mechanisms of persistent pain development. Further studies are needed in larger samples of individuals experiencing MVC, both to validate current findings and to provide greater power to discover associated miRNA. Additionally, further experiments are needed to show whether the miRNA identified in this study play a causal role in persistent pain pathogenesis. The results of such studies may provide an important new window into these yet enigmatic processes.

Methods

Study design and setting

This prospective longitudinal study enrolled African American individuals who presented within 24 h of MVC to one of eight EDs in three states (Michigan, Pennsylvania, and Florida) between July 2012 and July 2013. The study only enrolled African Americans because of the pressing need for pain studies that focus on such understudied, high risk groups [74, 7679]. The study was approved by the institutional review boards of all participating hospitals. Each participant provided written informed consent before enrollment.

Participant eligibility criteria

Individuals ≥18 and ≤65 years of age presenting to the ED within 24 h of MVC who did not have fracture or other injury requiring hospital admission were screened for eligibility. Patients who were not alert and oriented were excluded, as were patients who did not self-identify as African American, pregnant patients, prisoners, patients unable to read and understand English, or patients taking opioids above a total daily dose of 30 mg of oral morphine or equivalent.

Study procedures

Eligible and consenting participants provided a blood sample in the ED and completed an ED interview evaluation. Interview evaluations were performed by research assistants at the time of the ED visit using a web-based survey with explicit definitions of variables. Before enrolling patients in the ED, each research assistant completed a study training module followed by an interview with a standardized mock ED patient. Comparison of mock ED patient data across research assistants demonstrated an error rate of 0.57 %. Injury characteristics and medications administered in the ED were obtained by data extraction from the ED medical record. Six weeks after the MVC, participants completed a follow-up interview by telephone, online, or via mail. Participants were compensated $75 for completing the ED protocol and $50 for completing the 6-week interview.

Participant demographics

Participant demographic characteristics (including age, gender, income, height, weight, and educational attainment) were obtained from the ED medical record and from participant self-report.

Pain assessments and outcome definitions

Severity of pain in each body region during the month prior to MVC was assessed at the time of ED evaluation using a 0–-10 Numeric Rating Scale (NRS) score [80]. Severity of pain due to MVC in each body region was assessed at the 6 week time point using this same method together with an assessment of the MVC-relatedness of the pain. Individuals reporting a pain severity score ≥7 in at least one axial body region (neck, upper back, lower back, left shoulder, right shoulder) were defined as having severe AP [81, 82]. Individuals reporting severe AP during the month prior to the MVC were excluded from analyses.

RNA collection and isolation

Research assistants collected blood samples in the ED at the time of enrollment using PAXgene RNA tubes. Total RNA (including miRNA) was isolated using the PAXgene blood miRNA kit (Qiagen, Valencia, CA, USA) and stored at −80 °C until use. RNA concentration and purity were measured using a NanoDrop 1000 (Nanodrop Technologies, Wilmington, DE).

Library preparation and miRNA sequencing

Template libraries for miRNA Sequencing were produced from 1.0 ug total RNA using an adaptation of published protocols [83]. Briefly, total RNA was sequentially ligated to a 3′ linker using T4 RNA ligase 2, polyacrylamide gel extracted to remove excess 3′ linker, then ligated via T4 RNA ligase 1 to an oligonucleotide adapter (sequences shown in Additional file 4: Table S4). The 5′ adapter contained a two nucleotide barcode for multiplexing libraries. RNA products were reverse transcribed and amplified by PCR. In order to purify the miRNA population, gel isolation was used to obtain template libraries with 15–40 nucleotide inserts. Twelve barcoded libraries were combined per lane and sequenced on a HiSeq 2000 (Illumina, San Diego, CA, USA).

Bioinformatics analysis and data normalization

Raw sequence reads were processed using a custom bioinformatics pipeline. Reads were de-multiplexed and barcode and adapter sequences removed. Mature miRNA sequences were obtained from miRbase v18.0 and genomic extensions were added before aligning with sequencing reads. Total read counts were generated including isomir and non-templated nucleotide addition. Sequencing reads were normalized using quantile normalization. In order to avoid individual samples with lowly abundant or no miRNA expression, miRNA species with less than 300 reads across samples were dropped from analyses (adapted as described previously [28]).

RT-qPCR validation

The miRNA RT-qPCR method used is based on the stem-loop method described by Chen et al. [36]. Stem-loop RT primers and TaqMan probes for each miRNA were obtained from Life Technologies (Carlsbad, CA, USA). MiRNA expression of each miRNA was normalized to RNU48 levels before determination of expression differences. RT-qPCR validation was performed on a subset of the significant miRNA identified by miRNA sequencing (due to limited quantities of participant RNA).

Statistical analysis

Differences in ED miRNA expression between those who did and did not subsequently develop AP were quantified by dividing the mean expression levels in the two groups. Because miRNA distributions were non-normal [Kolmogorov–Smirnov (K-S) test], expression levels of individual miRNA among individuals with and without AP, including for women and men subgroup analyses, were compared using the Wilcoxon-Mann–Whitney (WMW) test. Logistic regression models were used to test for sex-miRNA interactions while adjusting for age and ED study site. To account for multiple testing, we used the methods of False Discovery Rate (FDR) determination defined by Benjamini and Hochberg [35]. For this initial proof-of-concept discovery cohort, p value thresholds were set at p <0.05 and an FDR cut-off corresponding to 0.15. Statistical analyses were carried out using SPSS software version 21.0 or SAS University Edition.

Identification of biologic pathways targeted by differentially expressed miRNA

A web-based computational tool, DIANA miRPath v2.0, was used to identify molecular pathways overrepresented in predicted targeting by differentially expressed miRNA originating from the X chromosome [79]. Pathway enrichment was also performed for the full set of 32 miRNA identified in Table 2. DIANA miRPath uses its predictive binding algorithm, DIANA-microT-CDS, to define a list of potential targets for each miRNA, then assigns a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway [41] rank and significance level based on the relative number of targets in that pathway [79]. DIANA miRPath results have been validated and its predictive binding algorithm has been shown to have high concordance with actual miRNA binding (e.g. [88, 90]). MiRNA that have been experimentally validated to bind to the predicted mRNA were identified using TarBase v2.0 [84]. The false discovery rates of the most highly ranked pathways were evaluated using permutation testing [42].

Abbreviations

MVC:

motor vehicle collision

miRNA:

microRNA

AP:

axial pain

RT-qPCR:

reverse transcription-quantitative polymerase chain reaction

ED:

emergency department

AA:

African American

KEGG:

Kyoto Encyclopedia of Genes and Genomes

TF:

transcription factor

NRS:

numeric rating scale

FDR:

false discovery rate

References

  1. He L, Hannon GJ. MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet. 2004;5:522–31.

    Article  CAS  PubMed  Google Scholar 

  2. Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6:857–66.

    Article  CAS  PubMed  Google Scholar 

  4. Thum T, Gross C, Fiedler J, et al. MicroRNA-21 contributes to myocardial disease by stimulating MAP kinase signalling in fibroblasts. Nature. 2008;456:980–4.

    Article  CAS  PubMed  Google Scholar 

  5. Kornfeld JW, Baitzel C, Konner AC, et al. Obesity-induced overexpression of miR-802 impairs glucose metabolism through silencing of Hnf1b. Nature. 2013;494:111–5.

    Article  CAS  PubMed  Google Scholar 

  6. Lutz BM, Bekker A, Tao YX. Noncoding RNAs: new players in chronic pain. Anesthesiology. 2014;121:409–17.

    Article  CAS  PubMed  Google Scholar 

  7. Bali KK, Kuner R. Noncoding RNAs: key molecules in understanding and treating pain. Trends in molecular medicine. 2014;20:437–48.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Andersen HH, Duroux M, Gazerani P. MicroRNAs as modulators and biomarkers of inflammatory and neuropathic pain conditions. Neurobiol Dis. 2014;71:159–68.

    Article  CAS  PubMed  Google Scholar 

  9. Niska R, Bhuiya F, Xu J. National Hospital Ambulatory Medical Care Survey: 2007 emergency department summary. Natl Health Stat Rep. 2010;(26):1–31.

  10. Bureau UC. Statistical abstract of the United States. Washington DC: US Census Bureau. 2012. http://www.census.gov/compendia/statab/. Accessed Jan 2014.

  11. Tanskanen A, Hintikka J, Honkalampi K, Haatainen K, Koivumaa-Honkanen H, Viinamaki H. Impact of multiple traumatic experiences on the persistence of depressive symptoms—a population-based study. Nord J Psychiatry. 2004;58:459–64.

    Article  PubMed  Google Scholar 

  12. Platts-Mills TF, Hunold KM, Esserman DA, Sloane PD, McLean SA. Motor vehicle collision-related emergency department visits by older adults in the United States. Acad Emerg Med Off J Soc Acad Emerg Med. 2012;19:821–7.

    Article  Google Scholar 

  13. Bortsov AV, Platts-Mills TF, Peak DA, et al. Effect of pain location and duration on life function in the year after motor vehicle collision. Pain. 2014;155:1836–45.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Whitehead A, Crawford DL. Variation in tissue-specific gene expression among natural populations. Genome Biol. 2005;6:R13.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Landgraf P, Rusu M, Sheridan R, et al. A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 2007;129:1401–14.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Sakai A, Suzuki H. Emerging roles of microRNAs in chronic pain. Neurochem Int. 2014;77:58–67.

    Article  CAS  PubMed  Google Scholar 

  17. Mitchell PS, Parkin RK, Kroh EM, et al. Circulating microRNAs as stable blood-based markers for cancer detection. Proc Natl Acad Sci USA. 2008;105:10513–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Arroyo JD, Chevillet JR, Kroh EM, et al. Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci USA. 2011;108:5003–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  19. Bahn S, Chan MK. What can we learn about depression from gene expression in peripheral tissues? Biol Psychiatry. 2015;77:207–9.

    Article  CAS  PubMed  Google Scholar 

  20. McLean SA. The potential contribution of stress systems to the transition to chronic whiplash-associated disorders. Spine. 2011;36:S226–32.

    Article  PubMed Central  PubMed  Google Scholar 

  21. McLean SA, Diatchenko L, Lee YM, et al. Catechol O-methyltransferase haplotype predicts immediate musculoskeletal neck pain and psychological symptoms after motor vehicle collision. J Pain Off J Am Pain Soc. 2011;12:101–7.

    Article  CAS  Google Scholar 

  22. Bortsov AV, Smith JE, Diatchenko L, et al. Polymorphisms in the glucocorticoid receptor co-chaperone FKBP5 predict persistent musculoskeletal pain after traumatic stress exposure. Pain. 2013;154(8):1419–26.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  23. Bortsov AV, Diatchenko L, McLean SA. Complex multilocus effects of catechol-O-methyltransferase haplotypes predict pain and pain interference 6 weeks after motor vehicle collision. Neuromol Med. 2014;16(1):83–93.

  24. Desborough JP. The stress response to trauma and surgery. Br J Anaesth. 2000;85:109–17.

    Article  CAS  PubMed  Google Scholar 

  25. Sullivan PF, Fan C, Perou CM. Evaluating the comparability of gene expression in blood and brain. Am J Med Genet Part B Neuropsychiatr Genet Off Publ Int Soc Psychiatr Genet. 2006;141B:261–8.

    Article  Google Scholar 

  26. Liu DZ, Tian Y, Ander BP, et al. Brain and blood microRNA expression profiling of ischemic stroke, intracerebral hemorrhage, and kainate seizures. J Cereb Blood Flow Metab Off J Int Soc Cereb Blood Flow Metab. 2010;30:92–101.

    Article  Google Scholar 

  27. Orlova IA, Alexander GM, Qureshi RA, et al. MicroRNA modulation in complex regional pain syndrome. J Transl Med. 2011;9:195.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Leidinger P, Backes C, Deutscher S, et al. A blood based 12-miRNA signature of Alzheimer disease patients. Genome Biol. 2013;14:R78.

    Article  PubMed Central  PubMed  Google Scholar 

  29. Keller A, Leidinger P, Steinmeyer F, et al. Comprehensive analysis of microRNA profiles in multiple sclerosis including next-generation sequencing. Mult Scler. 2013.

  30. Lai CY, Yu SL, Hsieh MH, et al. MicroRNA expression aberration as potential peripheral blood biomarkers for schizophrenia. PLoS One. 2011;6:e21635.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. Schultz NA, Dehlendorff C, Jensen BV, et al. MicroRNA biomarkers in whole blood for detection of pancreatic cancer. JAMA J Am Med Assoc. 2014;311:392–404.

    Article  CAS  Google Scholar 

  32. Alevizos I, Illei GG. MicroRNAs as biomarkers in rheumatic diseases. Nat Rev Rheumatol. 2010;6:391–8.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  33. Thomas A, Gennarelli EW. The Abbreviated Injury Scale 2005. In: AAAM. AAfAM, editor. Des Plaines, IL Update. 2008.

  34. Chen X, Ba Y, Ma L, et al. Characterization of microRNAs in serum: a novel class of biomarkers for diagnosis of cancer and other diseases. Cell Res. 2008;18:997–1006.

    Article  CAS  PubMed  Google Scholar 

  35. Benjamini YHY. Controlling the false discovery rate—a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995;57:289–300.

    Google Scholar 

  36. Chen C, Ridzon DA, Broomer AJ, et al. Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res. 2005;33:e179.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Xu J, Burgoyne PS, Arnold AP. Sex differences in sex chromosome gene expression in mouse brain. Hum Mol Genet. 2002;11:1409–19.

    Article  CAS  PubMed  Google Scholar 

  38. Nguyen DK, Disteche CM. Dosage compensation of the active X chromosome in mammals. Nat Genet. 2006;38:47–53.

    Article  CAS  PubMed  Google Scholar 

  39. Morgan CP, Bale TL. Sex differences in microRNA regulation of gene expression: no smoke, just miRs. Biol Sex Differ. 2012;3:22.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Vlachos IS, Kostoulas N, Vergoulis T, et al. DIANA miRPath v. 2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40:W498–504.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M. KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012;40:D109–14.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Tusher VG, Tibshirani R, Chu G. Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001;98:5116–21.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Rygh LJ, Svendsen F, Fiska A, Haugan F, Hole K, Tjolsen A. Long-term potentiation in spinal nociceptive systems—how acute pain may become chronic. Psychoneuroendocrinology. 2005;30:959–64.

    Article  PubMed  Google Scholar 

  44. Wood PB. Role of central dopamine in pain and analgesia. Expert Rev Neurother. 2008;8:781–97.

    Article  CAS  PubMed  Google Scholar 

  45. Pezet S, McMahon SB. Neurotrophins: mediators and modulators of pain. Annu Rev Neurosci. 2006;29:507–38.

    Article  CAS  PubMed  Google Scholar 

  46. McLean SA, Ulirsch JC, Slade GD, et al. Incidence and predictors of neck and widespread pain after motor vehicle collision among US litigants and nonlitigants. Pain. 2014;155:309–21.

    Article  PubMed Central  PubMed  Google Scholar 

  47. Carette S. Whiplash injury and chronic neck pain. N Engl J Med. 1994;330:1083–4.

    Article  CAS  PubMed  Google Scholar 

  48. Haqqani AS, Delaney CE, Tremblay TL, Sodja C, Sandhu JK, Stanimirovic DB. Method for isolation and molecular characterization of extracellular microvesicles released from brain endothelial cells. Fluids Barriers CNS. 2013;10:4.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Alvarez-Erviti L, Seow Y, Yin H, Betts C, Lakhal S, Wood MJ. Delivery of siRNA to the mouse brain by systemic injection of targeted exosomes. Nat Biotechnol. 2011;29:341–5.

    Article  CAS  PubMed  Google Scholar 

  50. Pegtel DM, Peferoen L, Amor S. Extracellular vesicles as modulators of cell-to-cell communication in the healthy and diseased brain. Biological sciences: Philos Transact R Soc Lond Ser B; 2014. p. 369.

    Google Scholar 

  51. Pegtel DM, Cosmopoulos K, Thorley-Lawson DA, et al. Functional delivery of viral miRNAs via exosomes. Proc Natl Acad Sci USA. 2010;107:6328–33.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Jung SH, Wang Y, Kim T, et al. Molecular mechanisms of repeated social defeat-induced glucocorticoid resistance: role of microRNA. Brain Behav Immun. 2015;44:195–206.

    Article  CAS  PubMed  Google Scholar 

  53. Wang WC, Juan AH, Panebra A, Liggett SB. MicroRNA let-7 establishes expression of beta2-adrenergic receptors and dynamically down-regulates agonist-promoted down-regulation. Proc Natl Acad Sci USA. 2011;108:6246–51.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  54. Khasar SG, Burkham J, Dina OA, et al. Stress induces a switch of intracellular signaling in sensory neurons in a model of generalized pain. J Neurosci Off J Soc Neurosci. 2008;28:5721–30.

    Article  CAS  Google Scholar 

  55. Sober S, Laan M, Annilo T. MicroRNAs miR-124 and miR-135a are potential regulators of the mineralocorticoid receptor gene (NR3C2) expression. Biochem Biophys Res Commun. 2010;391:727–32.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Issler O, Haramati S, Paul ED, et al. MicroRNA 135 is essential for chronic stress resiliency, antidepressant efficacy, and intact serotonergic activity. Neuron. 2014;83:344–60.

    Article  CAS  PubMed  Google Scholar 

  57. Hains BC, Everhart AW, Fullwood SD, Hulsebosch CE. Changes in serotonin, serotonin transporter expression and serotonin denervation supersensitivity: involvement in chronic central pain after spinal hemisection in the rat. Exp Neurol. 2002;175:347–62.

    Article  CAS  PubMed  Google Scholar 

  58. Lucas JJ, Mellstrom B, Colado MI, Naranjo JR. Molecular mechanisms of pain: serotonin1A receptor agonists trigger transactivation by c-fos of the prodynorphin gene in spinal cord neurons. Neuron. 1993;10:599–611.

    Article  CAS  PubMed  Google Scholar 

  59. Myers B, Greenwood-Van Meerveld B. Divergent effects of amygdala glucocorticoid and mineralocorticoid receptors in the regulation of visceral and somatic pain. American journal of physiology Gastrointestinal and liver physiology. 2010;298:G295–303.

    Article  CAS  PubMed  Google Scholar 

  60. Mannironi C, Camon J, De Vito F, et al. Acute stress alters amygdala microRNA miR-135a and miR-124 expression: inferences for corticosteroid dependent stress response. PLoS One. 2013;8:e73385.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  61. Genda Y, Arai M, Ishikawa M, Tanaka S, Okabe T, Sakamoto A. microRNA changes in the dorsal horn of the spinal cord of rats with chronic constriction injury: a TaqMan(R) low density array study. Int J Mol Med. 2013;31:129–37.

    CAS  PubMed  Google Scholar 

  62. Pohl KW, Yeol J-F, Ongl WY. MicroRNA changes in the mouse prefrontal cortex after inflammatory pain. Eur J Pain. 2012;15:801–12.

    Google Scholar 

  63. Volk N, Paul ED, Haramati S, et al. MicroRNA-19b associates with Ago2 in the amygdala following chronic stress and regulates the adrenergic receptor beta 1. J Neurosci Off J Soc Neurosci. 2014;34:15070–82.

    Article  Google Scholar 

  64. Muhie S, Hammamieh R, Cummings C, Yang D, Jett M. Transcriptome characterization of immune suppression from battlefield-like stress. Genes Immun. 2013;14:19–34.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Sakai A, Saitow F, Miyake N, Miyake K, Shimada T, Suzuki H. miR-7a alleviates the maintenance of neuropathic pain through regulation of neuronal excitability. Brain J Neurol. 2013;136:2738–50.

    Article  Google Scholar 

  66. Bhat-Nakshatri P, Wang G, Collins NR, et al. Estradiol-regulated microRNAs control estradiol response in breast cancer cells. Nucleic Acids Res. 2009;37:4850–61.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  67. Carrel L, Willard HF. X-inactivation profile reveals extensive variability in X-linked gene expression in females. Nature. 2005;434:400–4.

    Article  CAS  PubMed  Google Scholar 

  68. Wutz A. Gene silencing in X-chromosome inactivation: advances in understanding facultative heterochromatin formation. Nat Rev Genet. 2011;12:542–53.

    Article  CAS  PubMed  Google Scholar 

  69. Hewagama A, Gorelik G, Patel D, et al. Overexpression of X-linked genes in T cells from women with lupus. J Autoimmun. 2013;41:60–71.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  70. Morgan CP, Bale TL. Early prenatal stress epigenetically programs dysmasculinization in second-generation offspring via the paternal lineage. J Neurosci Off J Soc Neurosci. 2011;31:11748–55.

    Article  CAS  Google Scholar 

  71. Lawlis GF, Achterberg J, Kenner L, Kopetz K. Ethnic and sex differences in response to clinical and induced pain in chronic spinal pain patients. Spine. 1984;9:751–4.

    Article  CAS  PubMed  Google Scholar 

  72. Edwards RR, Fillingim RB. Ethnic differences in thermal pain responses. Psychosom Med. 1999;61:346–54.

    Article  CAS  PubMed  Google Scholar 

  73. Campbell CMER. Ethnic differences in pain and pain management. Pain Manage. 2012;2:219–30.

    Article  Google Scholar 

  74. Campbell CM, Edwards RR, Fillingim RB. Ethnic differences in responses to multiple experimental pain stimuli. Pain. 2005;113:20–6.

    Article  PubMed  Google Scholar 

  75. Lewis BP, Shih IH, Jones-Rhoades MW, Bartel DP, Burge CB. Prediction of mammalian microRNA targets. Cell. 2003;115:787–98.

    Article  CAS  PubMed  Google Scholar 

  76. Focus on pain. Nature Neurosci. 2014;17:145.

  77. Denk F, McMahon SB, Tracey I. Pain vulnerability: a neurobiological perspective. Nat Neurosci. 2014;17:192–200.

    Article  CAS  PubMed  Google Scholar 

  78. Rahim-Williams FB, Riley JL 3rd, Herrera D, Campbell CM, Hastie BA, Fillingim RB. Ethnic identity predicts experimental pain sensitivity in African Americans and Hispanics. Pain. 2007;129:177–84.

    Article  PubMed Central  PubMed  Google Scholar 

  79. McBeth J, Jones K. Epidemiology of chronic musculoskeletal pain. Best Prac Res Clin Rheumatol. 2007;21:403–25.

    Article  Google Scholar 

  80. Wolfe F. Pain extent and diagnosis: development and validation of the regional pain scale in 12,799 patients with rheumatic disease. J Rheumatol. 2003;30:369–78.

    PubMed  Google Scholar 

  81. Krebs EE, Carey TS, Weinberger M. Accuracy of the pain numeric rating scale as a screening test in primary care. J Gen Intern Med. 2007;22:1453–8.

    Article  PubMed Central  PubMed  Google Scholar 

  82. Fejer R, Jordan A, Hartvigsen J. Categorising the severity of neck pain: establishment of cut-points for use in clinical and epidemiological research. Pain. 2005;119:176–82.

    Article  PubMed  Google Scholar 

  83. Pfeffer S, Sewer A, Lagos-Quintana M, et al. Identification of microRNAs of the herpesvirus family. Nat Methods. 2005;2:269–76.

    Article  CAS  PubMed  Google Scholar 

  84. Vlachos IS, Paraskevopoulou MD, Karagkouni D, et al. DIANA-TarBase v7.0: indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. 2015;43:D153–9.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Authors’ contributions

SL and SM conceived the manuscript, SL, MW, and RS isolated RNA from patient samples and prepared small RNA libraries for sequencing, SL and SH performed bioinformatics analyses, SL, JP and EY performed statistical analyses, EZ, CL, PH, KD, CP, MV, BO, JJ, RS, and RD assisted with study design and were responsible for data collection at individual ED sites, SL made figures and tables, SL and SM contributed to manuscript design and writing. All authors read and approved the final manuscript.

Acknowledgements

The authors would like to thank the study participants for taking part in this study.

Funding for this study was provided by the National Institute of Arthritis, Musculoskeletal, and Skin Diseases (R01AR060852 : Samuel A. McLean) and by the Mayday Fund. Neither of the above funding agencies had any role in the design and conduct of the study, in the collection, management, analysis and interpretation of the data, or in the preparation, review, or approval of the manuscript.

Scientific Meeting Presentation: presented at the American Pain Society meeting, May 2014 in Tampa, FL.

Competing interests

The authors declare that they have no competing interests.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sarah D. Linnstaedt.

Additional files

12990_2015_69_MOESM1_ESM.docx

Additional file 1: Table S1. microRNA in whole blood collected from African Americans in the early aftermath of Motor Vehicle Collision (MVC) that are predictive of Axial Pain vs. Recovery 6 weeks after MVC trauma.

12990_2015_69_MOESM2_ESM.docx

Additional file 2: Table S2. microRNA in whole blood collected from African American individuals in the early aftermath of Motor Vehicle Collision (MVC) that are predictive of Axial Pain vs. Recovery 6 weeks after MVC trauma, as assessed independently in women and men.

12990_2015_69_MOESM3_ESM.docx

Additional file 3: Table S3. Top 10 DIANA miRPath predicted KEGG pathways enriched in targeting by all 32 miRNA differentially regulated in the early aftermath of MVC trauma in AA individuals who develop AP following MVC vs. those who recover.

12990_2015_69_MOESM4_ESM.docx

Additional file 4: Table S4. Oligonucleotide sequences of linker adapters used in library preparation for miRNA sequencing.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Linnstaedt, S.D., Walker, M.G., Parker, J.S. et al. MicroRNA circulating in the early aftermath of motor vehicle collision predict persistent pain development and suggest a role for microRNA in sex-specific pain differences. Mol Pain 11, 66 (2015). https://doi.org/10.1186/s12990-015-0069-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12990-015-0069-3

Keywords