Single-cell multi-omics sequencing reveals cell-specific transcriptomic and chromatin accessibility profiles in gut microbiome metabolite butyrate-produced pain modulation

Abstract

Pain is the most common symptom of temporomandibular joint (TMJ) disorders, which present significant clinical challenges due to their complexity and limited treatment options. Our previous study demonstrates that gut microbiome-derived butyrate is critical for the modulation of TMJ pain. In this study, we investigated its underlying mechanisms, and we found that oral administration of tributyrin, a prodrug of butyrate, not only significantly alleviated TMJ pain but also reversed the reduction in histone acetylation in the spinal trigeminal nucleus caudalis (Sp5C) under the TMJ pain condition. Using single-cell multi-omics sequencing, we profiled gene expression and chromatin accessibility in the Sp5C cells at the single-cell resolution. Bioinformatics analysis revealed that TMJ pain disrupted both the expression and chromatin accessibility of Nop14, Matk, Idh3b, Ndst2, and Tomm6 across four cell types in the Sp5C, and these alterations were reversed by tributyrin treatment. Specifically, Nop14 exhibited increased chromatin accessibility at its promoter region under TMJ pain condition, and knockdown of Nop14 in the Sp5C restored histone acetylation and alleviated TMJ pain. Together, our findings reveal cell-type–specific gene regulation that underlies butyrate-mediated epigenetic regulation of TMJ pain, which suggesting that targeting gut microbiome metabolites could develop a non-opioid novel therapy for TMJ disorders.

References

[1] Gauer RL, Semidey MJ. Diagnosis and treatment of temporomandibular disorders. Am. Fam. Phys., 2015, 91: 378-386

[2] Ibi M. Inflammation and temporomandibular joint derangement. Biol. Pharm. Bull., 2019, 42: 538-542

[3] Liu F, Steinkeler A. Epidemiology, diagnosis, and treatment of temporomandibular disorders. Dent. Clin. North Am., 2013, 57: 465-479

[4] Tao, R., Liu, S., Crawford, J. & Tao, F. Gut-brain crosstalk and the central mechanisms of orofacial pain. Brain Sci.13, https://doi.org/10.3390/brainsci13101456 (2023).

[5] Zhao C, et al. . Sialic acid exacerbates gut dysbiosis-associated mastitis through the microbiota-gut-mammary axis by fueling gut microbiota disruption. Microbiome, 2023, 11 78

[6] Ma Y, et al. . Resveratrol alleviates temporomandibular joint inflammatory pain by recovering disturbed gut microbiota. Brain Behav. Immun., 2020, 87: 455-464

[7] Cho KH, et al. . Lactobacillus (LA-1) and butyrate inhibit osteoarthritis by controlling autophagy and inflammatory cell death of chondrocytes. Front Immunol., 2022, 13 930511

[8] Bonomo RR, et al. . Fecal transplantation and butyrate improve neuropathic pain, modify immune cell profile, and gene expression in the PNS of obese mice. Proc. Natl. Acad. Sci. USA, 2020, 117: 26482-26493

[9] Stilling RM, et al. . The neuropharmacology of butyrate: the bread and butter of the microbiota-gut-brain axis?. Neurochem Int, 2016, 99: 110-132

[10] Wang R, et al. . Gut microbiota regulates acute myeloid leukaemia via alteration of intestinal barrier function mediated by butyrate. Nat. Commun., 2022, 13 2522

[11] Zhang L, Liu C, Jiang Q, Yin Y. Butyrate in energy metabolism: there is still more to learn. Trends Endocrinol. Metab., 2021, 32: 159-169

[12] Guy N, Chalus M, Dallel R, Voisin DL. Both oral and caudal parts of the spinal trigeminal nucleus project to the somatosensory thalamus in the rat. Eur. J. Neurosci., 2005, 21: 741-754

[13] Malmierca E, Martin YB, Nuñez A. Inhibitory control of nociceptive responses of trigeminal spinal nucleus cells by somatosensory corticofugal projection in rat. Neuroscience, 2012, 221: 115-124

[14] Parker JB, Palchaudhuri S, Yin H, Wei J, Chakravarti D. A transcriptional regulatory role of the THAP11-HCF-1 complex in colon cancer cell function. Mol. Cell Biol., 2012, 32: 1654-1670

[15] Deng Y, et al. . Zbtb14 regulates monocyte and macrophage development through inhibiting pu.1 expression in zebrafish. eLife, 2022, 11 e80760

[16] Moore DL, Apara A, Goldberg JL. Krüppel-like transcription factors in the nervous system: novel players in neurite outgrowth and axon regeneration. Mol. Cell Neurosci., 2011, 47: 233-243

[17] Chen J, et al. . JunD represses transcription and translation of the tight junction protein zona occludens-1 modulating intestinal epithelial barrier function. Mol. Biol. Cell, 2008, 19: 3701-3712

[18] Piper M, et al. . NFIB-mediated repression of the epigenetic factor Ezh2 regulates cortical development. J. Neurosci., 2014, 34: 2921-2930

[19] Sakurai S, Shimizu T, Ohto U. The crystal structure of the AhRR-ARNT heterodimer reveals the structural basis of the repression of AhR-mediated transcription. J. Biol. Chem., 2017, 292: 17609-17616

[20] Yfantis, A. et al. Transcriptional response to hypoxia: the role of HIF-1-associated co-regulators. Cells12, https://doi.org/10.3390/cells12050798 (2023).

[21] Yu F, et al. . The transcription factor Bhlhe40 is a switch of inflammatory versus antiinflammatory Th1 cell fate determination. J. Exp. Med., 2018, 215: 1813-1821

[22] Takeda M, et al. . Interaction with Smad4 is indispensable for suppression of BMP signaling by c-Ski. Mol. Biol. Cell, 2003, 15: 963-972

[23] Liu C-F, Lefebvre V. The transcription factors SOX9 and SOX5/SOX6 cooperate genome-wide through super-enhancers to drive chondrogenesis. Nucleic Acids Res., 2015, 43: 8183-8203

[24] Mitra P, et al. . Identification of HiNF-P, a key activator of cell cycle-controlled histone H4 genes at the onset of S phase. Mol. Cell Biol., 2003, 23: 8110-8123

[25] Hadjiagapiou C, et al. . Role of USF1 and USF2 as potential repressor proteins for human intestinal monocarboxylate transporter 1 promoter. Am. J. Physiol. Gastrointest. Liver Physiol., 2005, 288: G1118-G1126

[26] Weintraub AS, et al. . YY1 is a structural regulator of enhancer-promoter loops. Cell, 2017, 171: 1573-1588.e1528

[27] Armengol S, et al. . SREBP-2-driven transcriptional activation of human SND1 oncogene. Oncotarget, 2017, 8: 108181-108194

[28] Li X, McDonnell DP. The transcription factor B-Myb is maintained in an inhibited state in target cells through its interaction with the nuclear corepressors N-CoR and SMRT. Mol. Cell Biol., 2002, 22: 3663-3673

[29] Dengler VL, Galbraith M, Espinosa JM. Transcriptional regulation by hypoxia inducible factors. Crit. Rev. Biochem Mol. Biol., 2014, 49: 1-15

[30] Bommer GT, MacDougald OA. Regulation of lipid homeostasis by the bifunctional SREBF2-miR33a locus. Cell Metab., 2011, 13: 241-247

[31] Bret C, Desmots-Loyer F, Moreaux J, Fest T. BHLHE41, a transcriptional repressor involved in physiological processes and tumor development. Cell Oncol. (Dordr.), 2025, 48: 43-66

[32] Wang WL, Lee YC, Yang WM, Chang WC, Wang JM. Sumoylation of LAP1 is involved in the HDAC4-mediated repression of COX-2 transcription. Nucleic Acids Res., 2008, 36: 6066-6079

[33] Bakken TE, et al. . Single-nucleus and single-cell transcriptomes compared in matched cortical cell types. PLoS One, 2018, 13: e0209648

[34] van den Brink SC, et al. . Single-cell sequencing reveals dissociation-induced gene expression in tissue subpopulations. Nat. Methods, 2017, 14: 935-936

[35] Lanza, M. et al. Supplementation with SCFAs re-establishes microbiota composition and attenuates hyperalgesia and pain in a mouse model of NTG-induced migraine. Int. J. Mol. Sci.23, https://doi.org/10.3390/ijms23094847 (2022).

[36] Gao Y, Xie D, Wang Y, Niu L, Jiang H. Short-chain fatty acids reduce oligodendrocyte precursor cells loss by inhibiting the activation of astrocytes via the SGK1/IL-6 signalling pathway. Neurochem Res., 2022, 47: 3476-3489

[37] Kukkar A, Singh N, Jaggi AS. Attenuation of neuropathic pain by sodium butyrate in an experimental model of chronic constriction injury in rats. J. Formos. Med Assoc., 2014, 113: 921-928

[38] Szentirmai É, Millican NS, Massie AR, Kapás L. Butyrate, a metabolite of intestinal bacteria, enhances sleep. Sci. Rep., 2019, 9 7035

[39] Miyoshi M, et al. . Oral administration of tributyrin increases concentration of butyrate in the portal vein and prevents lipopolysaccharide-induced liver injury in rats. Clin. Nutr., 2011, 30: 252-258

[40] Koh A, De Vadder F, Kovatcheva-Datchary P, Backhed F. From dietary fiber to host physiology: short-chain fatty acids as key bacterial metabolites. Cell, 2016, 165: 1332-1345

[41] Wu J, Zhou Z, Hu Y, Dong S. Butyrate-induced GPR41 activation inhibits histone acetylation and cell growth. J. Genet Genom., 2012, 39: 375-384

[42] Capasso KE, et al. . Effect of histone deacetylase inhibitor JNJ-26481585 in pain. J. Mol. Neurosci., 2015, 55: 570-578

[43] Danaher RJ, et al. . Histone deacetylase inhibitors prevent persistent hypersensitivity in an orofacial neuropathic pain model. Mol. Pain., 2018, 14 1744806918796763

[44] Denk F, et al. . HDAC inhibitors attenuate the development of hypersensitivity in models of neuropathic pain. Pain, 2013, 154: 1668-1679

[45] Lou Y, et al. . Bioinformatics analysis of synovial fluid-derived mesenchymal stem cells in the temporomandibular joint stimulated with IL-1β. Cytotechnology, 2023, 75: 325-334

[46] Tao R, et al. . Using RNA-Seq to explore the hub genes in the trigeminal root entry zone of rats by compression injury. Pain. Physician, 2021, 24: E573-e581

[47] Hegenbarth, J.-C., Lezzoche, G., De Windt, L. J. & Stoll, M. Perspectives on bulk-tissue RNA sequencing and single-cell RNA sequencing for cardiac transcriptomics. Front. Mol. Med.2, https://doi.org/10.3389/fmmed.2022.839338 (2022).

[48] Lee J, Hyeon DY, Hwang D. Single-cell multiomics: technologies and data analysis methods. Exp. Mol. Med, 2020, 52: 1428-1442

[49] Aoki, Y. et al. Rs11726196 single-nucleotide polymorphism of the transient receptor potential canonical 3 (TRPC3) gene is associated with chronic pain. Int. J. Mol. Sci.24, https://doi.org/10.3390/ijms24021028 (2023).

[50] Johnston KJA, et al. . Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet, 2021, 17: e1009428

[51] Burgers LD, et al. . The natural product vioprolide A exerts anti-inflammatory actions through inhibition of its cellular target NOP14 and downregulation of importin-dependent NF-kB p65 nuclear translocation. Biomed. Pharmacother., 2021, 144 112255

[52] Engel P, et al. . The mRNA translation inhibitor vioprolide A prevents inflammatory pain-like behaviour with limited action on already established pain-like behaviour in mice. Eur. J. Pain., 2025, 29 e70099

[53] Liu PC, Thiele DJ. Novel stress-responsive genes EMG1 and NOP14 encode conserved, interacting proteins required for 40S ribosome biogenesis. Mol. Biol. Cell, 2001, 12: 3644-3657

[54] Yan X, et al. . NOP14-mediated ribosome biogenesis is required for mTORC2 activation and predicts rapamycin sensitivity. J. Biol. Chem., 2024, 300: 105681

[55] Wood H. Alzheimer disease: blocking HDAC2-Sp3 interaction - a new approach to AD therapy?. Nat. Rev. Neurol., 2017, 13: 577

[56] Formisano L, et al. . Sp3/REST/HDAC1/HDAC2 complex represses and Sp1/HIF-1/p300 complex activates ncx1 gene transcription, in brain ischemia and in ischemic brain preconditioning, by epigenetic mechanism. J. Neurosci., 2015, 35: 7332-7348

[57] Yamakawa H, et al. . The transcription factor Sp3 cooperates with HDAC2 to qregulate synaptic function and plasticity in neurons. Cell Rep., 2017, 20: 1319-1334

[58] Jia D, Rahbar R, Fish EN. Interferon-inducible Stat2 activation of JUND and CLDN4: mediators of IFN responses. J. Interferon Cytokine Res, 2007, 27: 559-565

[59] Folkerts J, de Bruijn MJW, van IWFJ, Hendriks RW, Stadhouders R. Butyrate selectively targets super-enhancers and transcriptional networks associated with human mast cell function. Eur. J. Immunol., 2025, 55 e51680

[60] Folkerts J, et al. . Butyrate inhibits human mast cell activation via epigenetic regulation of FcepsilonRI-mediated signaling. Allergy, 2020, 75: 1966-1978

[61] Rada-Iglesias A, et al. . Butyrate mediates decrease of histone acetylation centered on transcription start sites and down-regulation of associated genes. Genome Res, 2007, 17: 708-719

[62] Murray, R. L. et al. Tributyrin, a butyrate pro-drug, primes satellite cells for differentiation by altering the epigenetic landscape. Cells10, https://doi.org/10.3390/cells10123475 (2021).

[63] Wang X, et al. . A positive feedback inhibition of isocitrate dehydrogenase 3beta on paired-box gene 6 promotes Alzheimer-like pathology. Signal Transduct. Target Ther., 2024, 9: 105

[64] Carter NM, Ali S, Kirby JA. Endothelial inflammation: the role of differential expression of N-deacetylase/N-sulphotransferase enzymes in alteration of the immunological properties of heparan sulphate. J. Cell Sci., 2003, 116: 3591-3600

[65] Bai Q, et al. . TNFα in the trigeminal nociceptive system is critical for temporomandibular joint pain. Mol. Neurobiol., 2019, 56: 278-291

[66] Shu H, Liu S, Crawford J, Tao F. A female-specific role for trigeminal dynorphin in orofacial pain comorbidity. Pain, 2023, 164: 2801-2811

[67] Shu, H. et al. A pre-existing myogenic temporomandibular disorder increases trigeminal calcitonin gene-related peptide and enhances nitroglycerin-induced hypersensitivity in mice. Int. J. Mol. Sci.21, https://doi.org/10.3390/ijms21114049 (2020).

Funding

U.S. Department of Health & Human Services | NIH | National Institute of Dental and Craniofacial Research (NIDCR)(R01DE031255)

RIGHTS & PERMISSIONS

The Author(s)