A genome-wide Drosophila screen for heat nociception identifies α2δ3 as an evolutionarily conserved pain gene.

Worldwide, acute, and chronic pain affects 20% of the adult population and represents an enormous financial and emotional burden. Using genome-wide neuronal-specific RNAi knockdown in Drosophila, we report a global screen for an innate behavior and identify hundreds of genes implicated in heat nociception, including the α2δ family calcium channel subunit straightjacket (stj). Mice mutant for the stj ortholog CACNA2D3 (α2δ3) also exhibit impaired behavioral heat pain sensitivity. In addition, in humans, α2δ3 SNP variants associate with reduced sensitivity to acute noxious heat and chronic back pain. Functional imaging in α2δ3 mutant mice revealed impaired transmission of thermal pain-evoked signals from the thalamus to higher-order pain centers. Intriguingly, in α2δ3 mutant mice, thermal pain and tactile stimulation triggered strong cross-activation, or synesthesia, of brain regions involved in vision, olfaction, and hearing.

Genome-wide association studies with experimental validation identify a protective role for B lymphocytes against chronic post-surgical pain.

BACKGROUND: Chronic post-surgical pain (CPSP) significantly impacts patients' recovery and quality of life. Although environmental risk factors are well-established, genetic risk remains less understood. METHODS: A meta-analysis of genome-wide association studies followed by partitioned heritability was performed on 1350 individuals across five surgery types: hysterectomy, mastectomy, abdominal, hernia, and knee. In subsequent animal studies, withdrawal thresholds to evoked mechanical stimulation were measured in Rag1 null mutant and wild-type mice after plantar incision and laparotomy. Cell sorting by flow cytometry tracked recruitment of immune cell types. RESULTS: We discovered 77 genome-wide significant single-nucleotide polymorphism (SNP) hits, distributed among 24 loci and 244 genes. Meta-analysis of all cohorts estimated a SNP-based narrow-sense heritability for CPSP at ∼39%, indicating a substantial genetic contribution. Partitioned heritability analysis across a wide variety of tissues revealed enrichment of heritability in immune system-related genes, particularly those associated with B and T cells. Rag1 null mutant mice lacking both T and B cells exhibited exacerbated and prolonged allodynia up to 42 days after surgery, which was rescued by B-cell transfer. Recruitment patterns of B cells but not T cells differed significantly during the first 7 days after injury in the footpad, lymph nodes, and dorsal root ganglia. CONCLUSIONS: These findings suggest a key protective role for the adaptive immune system in the development of chronic post-surgical pain.

Genome-wide analysis identifies impaired axonogenesis in chronic overlapping pain conditions.

Chronic pain is often present at more than one anatomical location, leading to chronic overlapping pain conditions. Whether chronic overlapping pain conditions represent a distinct pathophysiology from the occurrence of pain at only one site is unknown. Using genome-wide approaches, we compared genetic determinants of chronic single-site versus multisite pain in the UK Biobank. We found that different genetic signals underlie chronic single-site and multisite pain with much stronger genetic contributions for the latter. Among 23 loci associated with multisite pain, nine loci replicated in the HUNT cohort, with the DCC netrin 1 receptor (DCC) as the top gene. Functional genomics identified axonogenesis in brain tissues as the major contributing pathway to chronic multisite pain. Finally, multimodal structural brain imaging analysis showed that DCC is most strongly expressed in subcortical limbic regions and is associated with alterations in the uncinate fasciculus microstructure, suggesting that DCC-dependent axonogenesis may contribute to chronic overlapping pain conditions via corticolimbic circuits.

Genetic risk factors for chronic postsurgical pain in children: A narrative review.

Genetic risk factors for chronic postsurgical pain in adults have been established, but little is known whether the same associations exist in children. It is even less clear how much influence single nucleotide polymorphisms can exert on the phenotypic expression of chronic postsurgical pain in children in general. To this effect, a search was made for original articles which met the following criteria: evaluation of postsurgical pain in children with known genetic mutations or, conversely, evaluation of atypical pain trajectories of postsurgical children assessing for possible genetic mutations that may explain the phenotype. All titles and abstracts retrieved were reviewed for suitability for inclusion. The references of the selected articles were also checked for additional relevant papers. To assess the transparency and quality of the genetic studies both STrengthening the REporting of Genetic Association studies scores and Q-Genie scores were applied. Overall, there is a paucity of information regarding the link between genetic mutations and eventual chronic postsurgical pain development although there is some information on acute postoperative pain. Evidence has shown that the contribution of genetic risk factors to chronic postsurgical pain development appears to be minor, with its clinical relevance yet to be described. More advanced techniques in systems biology (proteomics, transcriptomics) suggest promising avenues for investigating the disease.

Profound analgesia is associated with a truncated peptide resulting from tissue specific alternative splicing of DRG CA8-204 regulated by an exon-level cis-eQTL.

Carbonic anhydrase-8 (CA8) is an intracellular protein that functions as an allosteric inhibitor of inositol trisphosphate receptor-1 (ITPR1) critical to intracellular Ca++ release, synaptic functions and neuronal excitability. We showed previously that murine nociception and analgesic responses are regulated by the expression of this gene in dorsal root ganglion (DRG) associated with a cis-eQTL. In this report, we identify an exon-level cis-eQTL (rs6471859) that regulates human DRG CA8 alternative splicing, producing a truncated 1,697bp transcript (e.g., CA8-204). Our functional genomic studies show the "G" allele at rs6471859 produces a cryptic 3'UTR splice site regulating expression of CA8-204. We developed constructs to study the expression and function of the naturally occurring CA8-204G transcript (G allele at rs6471859), CA8-204C (C allele at rs6471859 reversion mutation) and CA8-201 (full length transcript). CA8-204G transcript expression occurred predominantly in non-neuronal cells (HEK293), while CA8-204C expression was restricted to neuronal derived cells (NBL) in vitro. CA8-204G produced a stable truncated transcript in HEK293 cells that was barely detectable in NBL cells. We also show CA8-204 produces a stable peptide that inhibits pITPR1 and Ca++ release in HEK293 cells. These results imply homozygous G/G individuals at rs6471859, which are common in the general population, produce exclusively CA8-204G that is barely detectable in neuronal cells. CA8 null mutations that greatly impact neuronal functions are associated with severe forms of spinal cerebellar ataxia, and our data suggest G/G homozygotes should display a similar phenotype. To address this question, we show in vivo using AAV8-FLAG-CA8-204G and AAV8-V5-CA8-201 gene transfer delivered via intra-neural sciatic nerve injection (SN), that these viral constructs are able to transduce DRG cells and produce similar analgesic and anti-hyperalgesic responses to inflammatory pain. Immunohistochemistry (IHC) examinations of DRG tissues further show CA8-204G peptide is expressed in advillin expressing neuronal cells, but to a lesser extent compared to glial cells. These findings explain why G/G homozygotes that exclusively produce this truncated functional peptide in DRG evade a severe phenotype. These genomic studies significantly advance the literature regarding structure-function studies on CA8-ITPR1 critical to calcium signaling pathways, synaptic functioning, neuronal excitability and analgesic responses.

A study in scarlet: MC1R as the main predictor of red hair and exemplar of the flip-flop effect.

Genetic variation in melanocortin-1 receptor (MC1R) is a known contributor to disease-free red hair in humans. Three loss-of-function single-nucleotide variants (rs1805007, rs1805008 and rs1805009) have been established as strongly correlated with red hair. The contribution of other loss-of-function MC1R variants (in particular rs1805005, rs2228479 and rs885479) and the extent to which other genetic loci are involved in red hair colour is less well understood. Here, we used the UK Biobank cohort to capture a comprehensive list of MC1R variants contributing to red hair colour. We report a correlation with red hair for both strong-effect variants (rs1805007, rs1805008 and rs1805009) and weak-effect variants (rs1805005, rs2228479 and rs885479) and show that their coefficients differ by two orders of magnitude. On the haplotype level, both strong- and weak-effect variants contribute to the red hair phenotype, but when considered individually, weak-effect variants show a reverse, negative association with red hair. The reversal of association direction in the single-variant analysis is facilitated by a distinguishing structure of MC1R, in which loss-of-function variants are never found to co-occur on the same haplotype. The other previously reported hair colour genes' variants do not substantially improve the MC1R red hair colour predictive model. Our best model for predicting red versus other hair colours yields an unparalleled area under the receiver operating characteristic of 0.96 using only MC1R variants. In summary, we present a comprehensive statistically derived characterization of the role of MC1R variants in red hair colour and offer a powerful, economical and parsimonious model that achieves unsurpassed performance.

Alternative Splicing of Opioid Receptor Genes Shows a Conserved Pattern for 6TM Receptor Variants.

The opioid receptor (OPR) family comprises the mu-, delta-, and kappa-opioid, and nociceptin receptors that belong to the superfamily of 7-transmembrane spanning G protein-coupled receptors (GPCRs). The mu-opioid receptor is the main target for clinically used opioid analgesics, and its biology has been extensively studied. The N-terminally truncated 6TM receptors isoform produced through alternative splicing of the OPRM1 gene displays unique signaling and analgesic properties, but it is unclear if other OPRs have the same ability. In this study, we have built a comprehensive map of alternative splicing events that produce 6TM receptor variants in all the OPRs and demonstrated their evolutionary conservation. We then obtained evidence for their translation through ribosomal footprint analysis. We discovered that N-terminally truncated 6TM GPCRs are rare in the human genome and OPRs are overrepresented in this group. Finally, we also observed a significant enrichment of 6TM GPCR genes among genes associated with pain, psychiatric disorders, and addiction. Understanding the biology of 6TM receptors and leveraging this knowledge for drug development should pave the way for novel therapies.

Reversion mutation of cDNA CA8-204 minigene construct produces a truncated functional peptide that regulates calcium release in vitro and produces profound analgesia in vivo.

Intracellular calcium is critical in orchestrating neuronal excitability and analgesia. Carbonic anhydrase-8 (CA8) regulates intracellular calcium signaling through allosteric inhibition of neuronal inositol trisphosphate receptor 1 (ITPR1) to produce profound analgesia. Recently, we reported the "G" allele at rs6471859 represents cis-eQTL regulating alternative splicing of a 1697 bp transcript (CA8-204G) with a retained intron, alternative polyadenylation site and a new stop codon producing a functional 26 kDa peptide with an extended exon 3. In this study we show the reversion mutation (G to C) at rs6471859 within the CA8-204G expression vector also produced a stable 1697 bp transcript (CA8-204C) coding for a smaller peptide (~ 22 kDa) containing only the first three CA8 exons. Surprisingly, this peptide inhibited ITPR1 (pITPR1) activation, ITPR1-mediated calcium release in vitro; and produced profound analgesia in vivo. This is the first report showing CA8-204C codes for a functional peptide sufficient to regulate calcium signaling and produce profound analgesia.

A functional substitution in the L-aromatic amino acid decarboxylase enzyme worsens somatic symptoms via a serotonergic pathway.

OBJECTIVE: Heightened somatic symptoms are reported by a wide range of patients with chronic pain and have been associated with emotional distress and physical dysfunction. Despite their clinical significance, molecular mechanisms leading to their manifestation are not understood. METHODS: We used an association study design based on a curated list of 3,295 single nucleotide polymorphisms mapped to 358 genes to test somatic symptoms reporting using the Pennebaker Inventory of Limbic Languidness questionnaire from a case-control cohort of orofacial pain (n = 1,607). A replication meta-analysis of 3 independent cohorts (n = 3,189) was followed by functional validation, including in silico molecular dynamics, in vitro enzyme assays, and measures of serotonin (5-HT) plasma concentration. RESULTS: An association with the T allele of rs11575542 coding for an arginine to glutamine substitution in the L-aromatic amino acid decarboxylase (AADC) enzyme was replicated in a meta-analysis of 3 independent cohorts. In a combined meta-analysis of all cohorts, this association reached p = 6.43 × 10-8 . In silico studies demonstrated that this substitution dramatically reduces the conformational dynamics of AADC, potentially lowering its binding capacity to a cofactor. in vitro enzymatic assays showed that this substitution reduces the maximum kinetic velocity of AADC, hence lowering 5-HT levels. Finally, plasma samples from 90 subjects showed correlation between low 5-HT levels and heightened somatic symptoms. INTERPRETATION: Using functional genomics approaches, we identified a polymorphism in the AADC enzyme that contributes to somatic symptoms through reduced levels of 5-HT. Our findings suggest a molecular mechanism underlying the pathophysiology of somatic symptoms and opens up new treatment options targeting the serotonergic system. ANN NEUROL 2019;86:168-180.

Stabilization of µ-opioid receptor facilitates its cellular translocation and signaling.

The G protein-coupled µ-opioid receptor (µ-OR) mediates the majority of analgesia effects for morphine and other pain relievers. Despite extensive studies of its structure and activation mechanisms, the inherently low maturation efficiency of µ-OR represents a major hurdle to understanding its function. Here we computationally designed µ-OR mutants with altered stability to probe the relationship between cell-surface targeting, signal transduction, and agonist efficacy. The stabilizing mutation T315Y enhanced µ-OR trafficking to the plasma membrane and significantly promoted the morphine-mediated inhibition of downstream signaling. In contrast, the destabilizing mutation R165Y led to intracellular retention of µ-OR and reduced the response to morphine stimulation. These findings suggest that µ-OR stability is an important factor in regulating receptor signaling and provide a viable avenue to improve the efficacy of analgesics.

N6-methyladenosine alters RNA structure to regulate binding of a low-complexity protein.

N6-methyladenosine (m6A) is the most abundant internal modification in eukaryotic messenger RNA (mRNA), and affects almost every stage of the mRNA life cycle. The YTH-domain proteins can specifically recognize m6A modification to control mRNA maturation, translation and decay. m6A can also alter RNA structures to affect RNA-protein interactions in cells. Here, we show that m6A increases the accessibility of its surrounding RNA sequence to bind heterogeneous nuclear ribonucleoprotein G (HNRNPG). Furthermore, HNRNPG binds m6A-methylated RNAs through its C-terminal low-complexity region, which self-assembles into large particles in vitro. The Arg-Gly-Gly repeats within the low-complexity region are required for binding to the RNA motif exposed by m6A methylation. We identified 13,191 m6A sites in the transcriptome that regulate RNA-HNRNPG interaction and thereby alter the expression and alternative splicing pattern of target mRNAs. Low-complexity regions are pervasive among mRNA binding proteins. Our results show that m6A-dependent RNA structural alterations can promote direct binding of m6A-modified RNAs to low-complexity regions in RNA binding proteins.

T-Cell Mediation of Pregnancy Analgesia Affecting Chronic Pain in Mice.

It has been reported consistently that many female chronic pain sufferers have an attenuation of symptoms during pregnancy. Rats display increased pain tolerance during pregnancy due to an increase in opioid receptors in the spinal cord. Past studies did not consider the role of non-neuronal cells, which are now known to play an important role in chronic pain processing. Using an inflammatory (complete Freund's adjuvant) or neuropathic (spared nerve injury) model of persistent pain, we observed that young adult female mice in early pregnancy switch from a microglia-independent to a microglia-dependent pain hypersensitivity mechanism. During late pregnancy, female mice show no evidence of chronic pain whatsoever. This pregnancy-related analgesia is reversible by intrathecal administration of naloxone, suggesting an opioid-mediated mechanism; pharmacological and genetic data suggest the importance of δ-opioid receptors. We also observe that T-cell-deficient (nude and Rag1-null mutant) pregnant mice do not exhibit pregnancy analgesia, which can be rescued with the adoptive transfer of CD4+ or CD8+ T cells from late-pregnant wild-type mice. These results suggest that T cells are a mediator of the opioid analgesia exhibited during pregnancy.SIGNIFICANCE STATEMENT Chronic pain symptoms often subside during pregnancy. This pregnancy-related analgesia has been demonstrated for acute pain in rats. Here, we show that pregnancy analgesia can produce a complete cessation of chronic pain behaviors in mice. We show that the phenomenon is dependent on pregnancy hormones (estrogen and progesterone), δ-opioid receptors, and T cells of the adaptive immune system. These findings add to the recent but growing evidence of sex-specific T-cell involvement in chronic pain processing.

The more you test, the more you find: The smallest P-values become increasingly enriched with real findings as more tests are conducted.

The increasing accessibility of data to researchers makes it possible to conduct massive amounts of statistical testing. Rather than follow specific scientific hypotheses with statistical analysis, researchers can now test many possible relationships and let statistics generate hypotheses for them. The field of genetic epidemiology is an illustrative case, where testing of candidate genetic variants for association with an outcome has been replaced by agnostic screening of the entire genome. Poor replication rates of candidate gene studies have improved dramatically with the increase in genomic coverage, due to factors such as adoption of better statistical practices and availability of larger sample sizes. Here, we suggest that another important factor behind the improved replicability of genome-wide scans is an increase in the amount of statistical testing itself. We show that an increase in the number of tested hypotheses increases the proportion of true associations among the variants with the smallest P-values. We develop statistical theory to quantify how the expected proportion of genuine signals (EPGS) among top hits depends on the number of tests. This enrichment of top hits by real findings holds regardless of whether genome-wide statistical significance has been reached in a study. Moreover, if we consider only those "failed" studies that produce no statistically significant results, the same enrichment phenomenon takes place: the proportion of true associations among top hits grows with the number of tests. The enrichment occurs even if the true signals are encountered at the logarithmically decreasing rate with the additional testing.

Human carbonic anhydrase-8 AAV8 gene therapy inhibits nerve growth factor signaling producing prolonged analgesia and anti-hyperalgesia in mice.

Carbonic anhydrase-8 (Car8; murine gene symbol) is an allosteric inhibitor of inositol trisphosphate receptor-1 (ITPR1), which regulates neuronal intracellular calcium release. We previously reported that wild-type Car8 overexpression corrects the baseline allodynia and hyperalgesia associated with calcium dysregulation in the waddle (wdl) mouse due to a 19 bp deletion in exon 8 of the Car8 gene. In this report, we provide preliminary evidence that overexpression of the human wild-type ortholog of Car8 (CA8WT), but not the reported CA8 S100P loss-of-function mutation (CA8MT), inhibits nerve growth factor (NGF)-induced phosphorylation of ITPR1, TrkA (NGF high-affinity receptor), and ITPR1-mediated cytoplasmic free calcium release in vitro. In addition, we show that gene transfer using AAV8-V5-CA8WT viral particles via sciatic nerve injection demonstrates retrograde transport to dorsal root ganglia (DRG) producing prolonged V5-CA8WT expression, pITPR1 and pTrkA inhibition, and profound analgesia and anti-hyperalgesia in male C57BL/6J mice. AAV8-V5-CA8WT-mediated overexpression prevented and treated allodynia and hyperalgesia associated with chronic neuropathic pain produced by the spinal nerve ligation (SNL) model. These AAV8-V5-CA8 data provide a proof-of-concept for precision medicine through targeted gene therapy of NGF-responsive somatosensory neurons as a long-acting local analgesic able to prevent and treat chronic neuropathic pain through regulating TrkA signaling, ITPR1 activation, and intracellular free calcium release by ITPR1.

Car8 dorsal root ganglion expression and genetic regulation of analgesic responses are associated with a cis-eQTL in mice.

Carbonic anhydrase-8 (Car8 mouse gene symbol) is devoid of enzymatic activity, but instead functions as an allosteric inhibitor of inositol trisphosphate receptor-1 (ITPR1) to regulate this intracellular calcium release channel important in synaptic functions and neuronal excitability. Causative mutations in ITPR1 and carbonic anhydrase-8 in mice and humans are associated with certain subtypes of spinal cerebellar ataxia (SCA). SCA mice are genetically deficient in dorsal root ganglia (DRG) Car8 expression and display mechanical and thermal hypersensitivity and susceptibility to subacute and chronic inflammatory pain behaviors. In this report, we show that DRG Car8 expression is variable across 25 naïve-inbred strains of mice, and this cis-regulated eQTL (association between rs27660559, rs27706398, and rs27688767 and DRG Car8 expression; P < 1 × 10-11) is correlated with nociceptive responses in mice. Next, we hypothesized that increasing DRG Car8 gene expression would inhibit intracellular calcium release required for morphine antinociception and might correlate with antinociceptive sensitivity of morphine and perhaps other analgesic agents. We show that mean DRG Car8 gene expression is directly related to the dose of morphine or clonidine needed to provide a half-maximal analgesic response (r = 0.93, P < 0.00002; r = 0.83, P < 0.0008, respectively), suggesting that greater DRG Car8 expression increases analgesic requirements. Finally, we show that morphine induces intracellular free calcium release using Fura 2 calcium imaging in a dose-dependent manner; V5-Car8 WT overexpression in NBL cells inhibits morphine-induced calcium increase. These findings highlight the 'morphine paradox' whereby morphine provides antinociception by increasing intracellular free calcium, while Car8 and other antinociceptive agents work by decreasing intracellular free calcium. This is the first study demonstrating that biologic variability associated with this cis-eQTL may contribute to differing analgesic responses through altered regulation of ITPR1-dependent calcium release in mice.

Genome-wide association meta-analyses to identify common genetic variants associated with hallux valgus in Caucasian and African Americans.

OBJECTIVE: Hallux valgus (HV) affects ∼36% of Caucasian adults. Although considered highly heritable, the underlying genetic determinants are unclear. We conducted the first genome-wide association study (GWAS) aimed to identify genetic variants associated with HV. METHODS: HV was assessed in three Caucasian cohorts (n=2263, n=915 and n=1231 participants, respectively). In each cohort, a GWAS was conducted using 2.5 M imputed SNPs. Mixed-effect regression with the additive genetic model adjusted for age, sex, weight and within-family correlations was used for both sex-specific and combined analyses. To combine GWAS results across cohorts, fixed-effect inverse-variance meta-analyses were used. Following meta-analyses, top-associated findings were also examined in an African American cohort (n=327). RESULTS: The proportion of HV variance explained by genome-wide genotyped SNPs was 50% in men and 48% in women. A higher proportion of genetic determinants of HV were sex specific. The most significantly associated SNP in men was rs9675316 located on chr17q23-a24 near the AXIN2 gene (p=0.000000546×10(-7)); the most significantly associated SNP in women was rs7996797 located on chr13q14.1-q14.2 near the ESD gene (p=0.000000721×10(-7)). Genome-wide significant SNP-by-sex interaction was found for SNP rs1563374 located on chr11p15.1 near the MRGPRX3 gene (interaction p value =0.0000000041×10(-9)). The association signals diminished when combining men and women. CONCLUSIONS: The findings suggest that the potential pathophysiological mechanisms of HV are complex and strongly underlined by sex-specific interactions. The identified genetic variants imply contribution of biological pathways observed in osteoarthritis as well as new pathways, influencing skeletal development and inflammation.

Construction of a global pain systems network highlights phospholipid signaling as a regulator of heat nociception.

The ability to perceive noxious stimuli is critical for an animal's survival in the face of environmental danger, and thus pain perception is likely to be under stringent evolutionary pressure. Using a neuronal-specific RNAi knock-down strategy in adult Drosophila, we recently completed a genome-wide functional annotation of heat nociception that allowed us to identify α2δ3 as a novel pain gene. Here we report construction of an evolutionary-conserved, system-level, global molecular pain network map. Our systems map is markedly enriched for multiple genes associated with human pain and predicts a plethora of novel candidate pain pathways. One central node of this pain network is phospholipid signaling, which has been implicated before in pain processing. To further investigate the role of phospholipid signaling in mammalian heat pain perception, we analysed the phenotype of PIP5Kα and PI3Kγ mutant mice. Intriguingly, both of these mice exhibit pronounced hypersensitivity to noxious heat and capsaicin-induced pain, which directly mapped through PI3Kγ kinase-dead knock-in mice to PI3Kγ lipid kinase activity. Using single primary sensory neuron recording, PI3Kγ function was mechanistically linked to a negative regulation of TRPV1 channel transduction. Our data provide a systems map for heat nociception and reinforces the extraordinary conservation of molecular mechanisms of nociception across different species.

Sex differences in mechanisms of pain hypersensitivity.

The introduction of sex-as-a-biological-variable policies at funding agencies around the world has led to an explosion of very recent observations of sex differences in the biology underlying pain. This review considers evidence of sexually dimorphic mechanisms mediating pain hypersensitivity, derived from modern assays of persistent pain in rodent animal models. Three well-studied findings are described in detail: the male-specific role of spinal cord microglia, the female-specific role of calcitonin gene-related peptide (CGRP), and the female-specific role of prolactin and its receptor. Other findings of sex-specific molecular involvement in pain are subjected to pathway analyses and reveal at least one novel hypothesis: that females may preferentially use Th1 and males Th2 T cell activity to mediate chronic pain.

Whole-genome methylation profiling reveals regions associated with painful temporomandibular disorders and active recovery processes.

ABSTRACT: Temporomandibular disorders (TMDs), collectively representing one of the most common chronic pain conditions, have a substantial genetic component, but genetic variation alone has not fully explained the heritability of TMD risk. Reasoning that the unexplained heritability may be because of DNA methylation, an epigenetic phenomenon, we measured genome-wide DNA methylation using the Illumina MethylationEPIC platform with blood samples from participants in the Orofacial Pain: Prospective Evaluation and Risk Assessment (OPPERA) study. Associations with chronic TMD used methylation data from 496 chronic painful TMD cases and 452 TMD-free controls. Changes in methylation between enrollment and a 6-month follow-up visit were determined for a separate sample of 62 people with recent-onset painful TMD. More than 750,000 individual CpG sites were examined for association with chronic painful TMD. Six differentially methylated regions were significantly ( P < 5 × 10 -8 ) associated with chronic painful TMD, including loci near genes involved in the regulation of inflammatory and neuronal response. A majority of loci were similarly differentially methylated in acute TMD consistent with observed transience or persistence of symptoms at follow-up. Functional characterization of the identified regions found relationships between methylation at these loci and nearby genetic variation contributing to chronic painful TMD and with gene expression of proximal genes. These findings reveal epigenetic contributions to chronic painful TMD through methylation of the genes FMOD , PM20D1 , ZNF718 , ZFP57 , and RNF39 , following the development of acute painful TMD. Epigenetic regulation of these genes likely contributes to the trajectory of transcriptional events in affected tissues leading to resolution or chronicity of pain.