Background <p>Temporomandibular disorders (TMD) show substantial clinical and genetic overlap with anxiety, yet it remains unclear whether TMD risk reflects shared anxiety-related liability or distinct anxiety-independent genetic mechanisms. Disentangling these components is essential for understanding TMD heterogeneity beyond symptom-based classifications.</p> Methods <p>We applied GWAS-by-subtraction using genome-wide summary statistics for TMD (20,799 cases and 479,549 controls; FinnGen Release 12) and anxiety disorders (74,973 cases and 400,243 controls), partitioning TMD heritability into two orthogonal latent components: an anxiety-dependent factor (F<sub>Anxiety</sub>) and an anxiety-independent factor (F<sub>Non-Anxiety</sub>). To delineate the mechanisms underlying each component, we integrated fine-mapping, transcriptome- and proteome-wide association analyses, genetic colocalization, brain imaging–genetics, and single-cell RNA sequencing from human embryonic temporomandibular joint tissue.</p> Results <p>Anxiety showed significant genetic correlation with TMD (rg = 0.4417, <i>p</i> = 1.98 × 10<sup>− 1 9</sup>) and accounted for 19.50% of TMD heritable variance. F<sub>Anxiety</sub> yielded multiple genome-wide significant loci (<i>CNTNAP5</i>, <i>PCLO</i>, <i>PRSS16</i>, <i>BTN1A1</i>, <i>RAB27B</i>), whereas F<sub>Non-Anxiety</sub> produced a single independent signal near <i>GPNMB</i>, demonstrating sharply divergent genetic architectures. Multi-omic integration identified <i>RAB27B</i> as a driver of the anxiety-related pathway, implicating synaptic vesicle trafficking and neuroimmune regulation, while <i>GPNMB</i> and <i>KLHL7</i> supported anxiety-independent pathways involving musculoskeletal remodeling and peripheral inflammation. BrainXcan analyses showed that F<sub>Anxiety</sub> predominantly affected limbic and external capsule microstructure, whereas F<sub>Non-Anxiety</sub> mapped to thalamic–sensorimotor white matter networks. Single-cell mapping further revealed distinct enrichment patterns of <i>RAB27B</i>, <i>KLHL7</i>, and <i>GPNMB</i> across TMJ cell types.</p> Conclusion <p>These findings demonstrate that TMD genetic liability comprises separable anxiety-related and anxiety-independent dimensions with distinct molecular, neurostructural, and cellular signatures. Rather than defining clinical subtypes, these latent components represent associative dimensions of genetic risk at the population level. This integrative framework clarifies the genetic architecture underlying TMD heterogeneity and provides a foundation for future studies integrating individual-level phenotyping to assess clinical relevance and causal mechanisms.</p> Graphical Abstract <p></p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Genetic subtraction reveals divergent pathways and targets in anxiety-related and anxiety-independent TMD

  • Yu Cao,
  • Xin Yang,
  • Peter Svensson,
  • Raymond Wong Chung Wen,
  • Timothy Jie Han Sng,
  • Intekhab Islam,
  • Wei Han,
  • Xingmei Feng,
  • Bozhi Hou,
  • Yuehua Li,
  • Lei Zheng

摘要

Background

Temporomandibular disorders (TMD) show substantial clinical and genetic overlap with anxiety, yet it remains unclear whether TMD risk reflects shared anxiety-related liability or distinct anxiety-independent genetic mechanisms. Disentangling these components is essential for understanding TMD heterogeneity beyond symptom-based classifications.

Methods

We applied GWAS-by-subtraction using genome-wide summary statistics for TMD (20,799 cases and 479,549 controls; FinnGen Release 12) and anxiety disorders (74,973 cases and 400,243 controls), partitioning TMD heritability into two orthogonal latent components: an anxiety-dependent factor (FAnxiety) and an anxiety-independent factor (FNon-Anxiety). To delineate the mechanisms underlying each component, we integrated fine-mapping, transcriptome- and proteome-wide association analyses, genetic colocalization, brain imaging–genetics, and single-cell RNA sequencing from human embryonic temporomandibular joint tissue.

Results

Anxiety showed significant genetic correlation with TMD (rg = 0.4417, p = 1.98 × 10− 1 9) and accounted for 19.50% of TMD heritable variance. FAnxiety yielded multiple genome-wide significant loci (CNTNAP5, PCLO, PRSS16, BTN1A1, RAB27B), whereas FNon-Anxiety produced a single independent signal near GPNMB, demonstrating sharply divergent genetic architectures. Multi-omic integration identified RAB27B as a driver of the anxiety-related pathway, implicating synaptic vesicle trafficking and neuroimmune regulation, while GPNMB and KLHL7 supported anxiety-independent pathways involving musculoskeletal remodeling and peripheral inflammation. BrainXcan analyses showed that FAnxiety predominantly affected limbic and external capsule microstructure, whereas FNon-Anxiety mapped to thalamic–sensorimotor white matter networks. Single-cell mapping further revealed distinct enrichment patterns of RAB27B, KLHL7, and GPNMB across TMJ cell types.

Conclusion

These findings demonstrate that TMD genetic liability comprises separable anxiety-related and anxiety-independent dimensions with distinct molecular, neurostructural, and cellular signatures. Rather than defining clinical subtypes, these latent components represent associative dimensions of genetic risk at the population level. This integrative framework clarifies the genetic architecture underlying TMD heterogeneity and provides a foundation for future studies integrating individual-level phenotyping to assess clinical relevance and causal mechanisms.

Graphical Abstract