Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank ( N = 172,230)

Abstract

Objectives:

Widespread pain is a hallmark characteristic of fibromyalgia, commonly affecting older individuals. This study aimed to identify novel genetic variants associated with widespread pain by utilizing the extensive UK Biobank dataset.

Methods:

We conducted a primary genome-wide association study (GWAS) using a novel definition of widespread pain, defined as pain experienced all over the body during the past month. Sex-stratified GWAS analysis approach was also performed to analyze the impact of sex on widespread pain.

Results:

The primary GWAS identified one novel significant genetic locus (rs34691025, p = 1.76 × 10⁻⁸) on chromosome 5q13.2 within the ARHGEF28 gene and several loci that approached genome-wide significance. The sex-stratified GWAS outputs revealed biological difference widespread pain between males and females, with a novel locus identified in the female-specific analysis within the LRMDA gene on chromosome 10. Genetic Correlation analysis demonstrated significant genetic correlations between widespread pain and other phenotypes, including joint disorders and spondylosis. The PheWAS revealed associations between the significant genetic variants with hearing disorders and cardiovascular diseases. A two-sample Mendelian randomization analysis found no significant causal association between hearing loss and widespread pain.

Conclusions:

Our study advances the understanding of the genetic factors contributing to widespread pain, highlighting notable differences between males and females and identifying a novel genetic locus associated with this condition.

Keywords

Widespread pain UK Biobank genetic correlations genome-wide association study phenome-wide association study Mendelian randomization

Introduction

Widespread pain, commonly characterized by pain experienced in various parts of the body simultaneously, stands as a notable concern in clinical practice.¹ Unlike localized pain, which is confined to a specific region, widespread pain manifests across multiple regions, often without an apparent localized source. This type of pain is frequently observed in conditions such as fibromyalgia, rheumatoid arthritis, and certain neurological disorders, presenting a challenge in both diagnosis and treatment.²

A commonly accepted definition of widespread pain is chronic widespread pain (CWP), which is characterized by the presence of widespread pain persisting for at least 3 months. A substantial body of epidemiological research suggests that the incidence and prevalence of CWP are significant, although these figures vary across different studies and countries. A systematic review indicates that the prevalence of CWP in the general population is approximately 10.6%, while in the UK, it is estimated to be 14.2%.^3,4 The likelihood and risk factors of CWP increase with age, and females are more susceptible than males.⁴ Beyond age and sex, CWP is associated with various risk factors, including sleep disorders, headaches and other types of pain, depression, and illness behavior.⁵ The escalating issue of population aging has exacerbated the financial burden of CWP, leading to significant social costs. In the US, the direct and indirect annual expenses associated with CWP per patient are estimated to be approximately $12,428.⁶ Although the definition of CWP differs slightly from our focus on widespread pain, it is still meaningful to discuss the current epidemiological and genetic studies of CWP, as this information could reflect the nature of widespread pain.

The genetic basis of CWP has been the subject of extensive research, focusing particularly on its genetic architecture. Twin studies have played a pivotal role in dissecting the genetic and environmental contributions to CWP, with estimates indicating that genetic factors account for approximately 48%–52% of the variance in CWP occurrence, underscoring a significant genetic component.⁷ In a comprehensive GWAS meta-analysis that amalgamated data from 14 studies, a notable locus was identified on chromosome 5, positioned intergenically between CCT5 and FAM173B.⁸ Furthermore, a recent study by Rahman et al.⁹ utilizing the UK Biobank data on CWP identified the RNF123 locus as being associated with CWP. This association was specifically observed in cases reporting pain all over the body and persisting for over 3 months.

Given the extensive research on CWP, our study shifts the focus to widespread pain experienced within the past month – a broader and less explored definition in the literature. This approach allows for the investigation of both short-term and long-term widespread pain, which may involve distinct genetic mechanisms compared to chronic conditions alone. The aim of this study is to identify novel genetic variants linked to widespread pain in the UK Biobank by conducting the GWAS research. We employ a unique definition that considers pain experienced all over the body during the past month for the GWAS and utilize three independent replication datasets (FibroGene, FinnGen, and TwinsUK). Sex-stratified GWAS analysis approach is firstly conducted to analyze the impact of sex on widespread pain. Through this research, we hope to identify new loci associated with widespread pain, thereby uncovering novel genetic mechanisms and offering fresh insights for the treatment of this condition.

Methods

Cohorts’ information

The UK Biobank was employed as the discovery cohort in our study. This cohort consists of over 500,000 participants aged 40–69 from England, Scotland, and Wales, recruited between 2006 and 2011. They underwent clinical examinations, completed comprehensive questionnaires, and provided DNA samples with informed consent for research use. The study received ethical approval from the UK’s National Health Service National Research Ethics Service (reference 11/NW/0382). Further details can be found at www.ukbiobank.ac.uk.

The DNA extraction and quality control (QC) processes were standardized, with detailed methodologies available at https://biobank.ctsu.ox.ac.uk/crystal/ukb/docs/genotyping_sample_workflow.pdf. The Welcome Trust Centre for Human Genetics at Oxford University oversaw the standardized QC for genotyping results, with a comprehensive description at http://biobank.ctsu.ox.ac.uk/crystal/refer.cgi?id=155580. In July 2017, the UK Biobank released genetic data, including genotyped and imputed genotypes, from 501,708 samples to authorized researchers. Bycroft et al.¹⁰ provided a detailed description of the QC procedures for imputation.

In this study, a specific pain-related question developed by the UK Biobank was utilized: “In the last month, have you experienced any of the following that interfered with your usual activities?” Participants could select from the following options: (1) Headache, (2) Facial pain, (3) Neck or shoulder pain, (4) Back pain, (5) Stomach or abdominal pain, (6) Hip pain, (7) Knee pain, (8) Pain all over the body, (9) None of the above, (10) Prefer not to say (UK Biobank Questionnaire field ID: 6159). Multiple selections were allowed. Cases of widespread pain were identified by the selection of “Pain all over the body,” regardless of other choices. Controls were those who chose “None of the above.” The analysis excluded data from individuals who were not of white British genetic ancestry, as initially defined by self-identification (UK Biobank Questionnaire field ID: 21000). A principal component analysis was then performed, followed by a random forest applied to the projected principal component data to reassign the initial self-defined ancestries for individuals with a membership posterior probability >0.5. Individuals with a posterior probability below 0.5 for any ancestry group were excluded from further analysis. Additionally, we excluded those with self-reported diagnoses of rheumatoid arthritis, polymyalgia rheumatica, unspecified arthritis, systemic lupus erythematosus, ankylosing spondylitis, and myopathy (UK Biobank Questionnaire field ID: 20002).

In the replication phase, we employed publicly available summary statistics for the phenotype of fibromyalgia as classified by ICD-10 in the FinnGen dataset, as well as summary statistics from the FibroGene fibromyalgia study and the CWP study within TwinsUK. The FinnGen dataset comprised 3761 cases diagnosed with fibromyalgia and 375,928 controls without the diagnosis. The GWAS for the FinnGen cohort was conducted according to the protocols of FinnGen GWAS round 9, with details available at https://finngen.gitbook.io/documentation.¹¹ Detailed information on sample phenotyping, genotyping, and the GWAS methodology applied to the FinnGen sample can be found at https://risteys.finngen.fi/.¹¹ The FibroGene dataset included 905 fibromyalgia cases recruited via advertisements in Fibromyalgia Action UK charity resources, and 14,213 controls from the general population program of the UK BioResource project. Details on genotyping and GWAS methodology for the FibroGene dataset can be referenced at https://bioresource.nihr.ac.uk/media/nlybpdcd/uk_axiom_biobank_genotyping_arrays_datasheet.pdf.¹² The TwinsUK GWAS study focusing on CWP comprised 903 individuals with CWP and 2369 control subjects The TwinsUK cohort, established as a registry of volunteer same-sex twins across the United Kingdom, encompasses a broad age range from 16 to 98 years. Initiated in 1992, the cohort has been continuously gathering a wealth of data and biological materials from its 14,274 registered twins. The ethical integrity of the study is upheld through approval from the Research Ethics Committee at Guy’s and St. Thomas’ NHS Foundation Trust. A significant aspect of the TwinsUK cohort is its repository of genetic information from 6921 participants. This genetic data has undergone quality control, documented by Moayyeri et al.¹³ For SNPs missing in the replication datasets, we employed a strategy based on linkage disequilibrium and physical distance to identify and use the most appropriate SNP for replication. Additionally, it is important to note that the replication p-values have not been adjusted for multiple testing across the 14 genes tested for replication.

GWAS and statistical analysis

Our study comprised three independent GWAS. The primary GWAS was conducted with the objective of understanding the genetic foundation of widespread pain across the entire dataset. The secondary GWAS was designed to assess the influence of sex on widespread pain, with the analysis specifically focused on conducting separate GWAS for each sex. In this study, the Genome-wide Complex Trait Analysis (GCTA, v1.94.1) software (available at https://yanglab.westlake.edu.cn/software/gcta/#Overview) was employed as the principal tool for executing GWAS.¹⁴ Our research utilized the fastGWA function within GCTA, a mixed linear model association tool, for the GWAS analyses. Standard QC procedures were implemented, which involved the exclusion of single nucleotide polymorphisms (SNPs) with INFO scores below 0.3, SNPs with minor allele frequencies lower than 0.5%, or SNPs failing the Hardy-Weinberg tests (p < 10⁻⁶). Furthermore, SNPs located on the X and Y chromosomes, as well as mitochondrial SNPs, were excluded from the analysis. The association tests were conducted using the fastGWA function, with adjustments for age, sex, body mass index (BMI), and ten population principal components. In the sex-stratified GWAS, participants whose self-reported sex did not match their genetically determined sex were excluded from the analysis. A $χ$ ² test was utilized to assess sex differences between the cases and controls, while the comparison of age and BMI was performed using an independent t-test, employing R v4.2.2. R v4.2.2 was also used to select the data of white British descent in UK Biobank data and divide the cases and controls of widespread pain. A p-value less than 5 × 10⁻⁸ was considered indicative of genome-wide association significance. Additionally, GCTA was also employed to compute the narrow-sense heritability.

GWAS-associated analysis by FUMA and LDSC

The FUMA web application served as the primary annotation tool throughout widespread pain study.¹⁵ Additionally, we generated both Manhattan and Q-Q plots using this application. For regional visualization, we employed LocusZoom (http://locuszoom.org/).¹⁶ Notably, FUMA facilitated three key types of analyses, gene analysis, gene-set analysis, and tissue expression analysis. In the gene analysis, the summary statistics of SNPs were aggregated to the level of entire genes, allowing for the assessment of associations between genes and the phenotype under investigation. In the gene-set analysis, specific groups of genes sharing common biological, functional, or other characteristics were collectively tested, providing insights into the involvement of distinct biological pathways or cellular functions in the genetic basis of the observed phenotype. The tissue expression analysis was conducted using data from GTEx (https://www.gtexportal.org/home), which has been integrated into the FUMA platform. For this analysis, the average gene expression per tissue type was used as a gene covariate, enabling the examination of relationships between gene expression in specific tissue types and the genetic associations with widespread pain, the focus of our study. We generated two gene expression heatmaps for genes identified through positional mapping. These heatmaps represent the average of normalized expression values, employing data from GTEx version 8 across 54 and 30 general tissue types. Tissue specificity was assessed by examining the differentially expressed genes (DEG) for each tissue type. Furthermore, an enrichment test for DEGs was performed for both sets of tissue types in GTEx version 8, focusing specifically on genes identified through positional mapping rather than utilizing the full distribution of SNP p-values in MAGMA tissue expression analysis.

To explore potential genetic correlations between widespread pain and various other phenotypes, we assessed the genetic associations of widespread pain with 1396 characteristics from the UK Biobank through the Complex Trait Virtual Lab (CTG-VL; https://genoma.io/). CTG-VL is an open-source platform that integrates available GWAS datasets to enable the estimation of genetic correlations for complex traits. Subsequently, to further investigate any potential genetic differences between the sexes, we utilized linkage disequilibrium score regression (LDSC) via LDSCv1.0.1 (available at https://github.com/bulik/ldsc) to examine the genetic correlation of widespread pain between males and females. LDSC utilizes the principle of linkage disequilibrium, which refers to the non-random association of alleles at different loci within a population.¹⁷

Expression quantitative trait loci (eQTL), chromatin interaction analysis, and positional mapping

Expression quantitative trait loci (eQTL) have become instrumental in elucidating the regulatory mechanisms of variants identified through GWAS. Cis-eQTL notably impact gene expression by interacting with variants that are situated in proximity (within 1 Mb) to the gene.¹⁸ In eukaryotic cells, the genome is compactly organized within the micron-sized nucleus, with chromatin serving as the fundamental structural unit. This organization is essential, arranging the genome into a three-dimensional structure that is pivotal for processes such as DNA replication, DNA damage repair, gene transcription, and other crucial biological functions.¹⁹ Positional mapping within our study was conducted using a maximum distance of 10 kb. The integration of cis-eQTL analysis, chromatin interaction analysis, and positional mapping in our gene mapping approach offers a comprehensive perspective on the genomic architecture underpinning these findings. A summary of the parameters and data used for positional mapping analysis has been added to Supplemental Table 1.

Phenome-wide association analysis (PheWAS)

Our study implemented a Phenome-Wide Association Analysis (PheWAS) to investigate the relationships between significant SNP associations, their corresponding genes, and multiple phenotypes. The primary objectives of this analysis were to corroborate the associations identified in our GWAS with pain phenotypes and to uncover new associations between genetic variants linked to widespread pain and other phenotypes. The PheWAS was conducted utilizing an expansive dataset comprising 4756 GWAS summary statistics available on the GWAS ATLAS platform https://atlas.ctglab.nl/PheWAS.²⁰ For this analysis, only SNPs exhibiting associations with p-values lower than 0.05 were included. Furthermore, adjustments for multiple comparisons were performed using the Bonferroni correction method.

Mendelian randomization

This study utilized a two-sample Mendelian Randomization (MR) approach to explore the potential causal relationship between hearing loss and widespread pain. Based on our genetic correlation analysis and PheWAS findings, hearing loss was linked to both widespread pain and the identified genetic variants. Consequently, MR was employed to evaluate the causal impact of hearing loss on widespread pain. Genetic associations for hearing loss were obtained from the UK Biobank genetic databases via the IEU Open GWAS platform (https://gwas.mrcieu.ac.uk/), including 14,654 cases and 474,839 controls. The genetic data on widespread pain were derived from our primary GWAS, which involved 4617 cases and 167,613 controls. To ensure robustness in our sensitivity analyses, we applied Inverse Variance Weighted (IVW) estimation alongside MR Egger, Weighted Median, Simple Mode, and Weighted Mode methods. To address potential biases, we performed heterogeneity assessments using Cochran’s Q test and examined horizontal pleiotropy with the MR Egger intercept test. Furthermore, bidirectional MR analyses were conducted to investigate the possibility of reverse causation.

Results

GWAS results

During the initial assessment phase (2006–2010) of the UK Biobank study, a total of 501,708 participants were administered a pain questionnaire. In our research, we excluded non-white British participants and samples that failed to meet QC criteria. In the primary GWAS, we identified 4617 cases (comprising 1885 males and 2732 females) and 167,613 controls (consisting of 80,228 males and 87,385 females) for analysis. The study utilized 11,165,459 SNPs for the GWAS examination. The secondary phase of the study involved sex-stratified GWAS analyses, where, following identical QC procedures, the female cohort (90,117 samples) consisted of 2732 cases and 87,385 controls, while the male cohort (82,113 samples) included 1885 cases and 80,228 controls. Table 1 provides a comprehensive summary of the clinical characteristics of cases and controls in the primary GWAS. Table 2 presents a detailed overview of the cases and controls in the secondary sex-stratified GWAS.

Table 1.

Clinical characteristics of widespread pain cases and controls in the UK Biobank.

Covariates	Cases	Controls	p-Value
Sex (male:female)	1885:2732	80,228:87,385	<0.001
Age (years)	56.7 (7.67)	56.9 (7.95)	0.55
BMI (kg/m³)	29.6 (5.77)	26.7 (4.28)	<0.001

BMI: body mass index.

A χ² test was used to test the difference of gender frequency between cases and controls and an independent t test was used for other covariates. Continuous covariates were presented as mean (standard deviation).

Table 2.

Clinical characteristics of widespread pain cases and controls in the UK Biobank (sex-stratified GWAS).

GWAS analysis	Covariates	Cases	Controls	p-Value
Male-specific GWAS	Age (years)	57.3 (7.60)	57.2 (8.04)	0.054
Male-specific GWAS	BMI (kg/m³)	29.6 (5.11)	27.3 (3.90)	<0.001
Female-specific GWAS	Age (years)	57.2 (7.58)	56.8 (7.87)	<0.05
Female-specific GWAS	BMI (kg/m³)	29.8 (6.36)	26.2 (4.57)	<0.001

BMI: body mass index.

Total number of females: 90,117, Total number of males: 82,113. An independent t test was used for other covariates.

Continuous covariates were presented as mean (standard deviation).

In our study, we identified one novel and distinct SNP cluster on chromosome 5 that exhibited significant associations with widespread pain, achieving genome-wide significance (p < 5 × 10⁻⁸), as delineated in Figure 1. Additionally, we identified several significant loci proximal to the genome-wide significant threshold, which are listed in Table 3. This table enumerates a total of 14 SNPs, encompassing both significant and marginally significant SNPs, with each representing the most statistically significant association within its respective locus. A detailed enumeration of all significantly associated SNPs from this GWAS is furnished in Supplemental Table 2. Notably, the significant association was observed in the SNP cluster situated within the ARHGEF28 gene on chromosome 5q13.2, with a p-value of 1.76 × 10⁻⁸ for rs34691025. A suggestive association was observed for the gene FAF1 on chromosome 1, where the most significant SNP was rs17387024 (p < 7.84 × 10⁻⁸). The regional plot illustrating the most significant locus in ARHGEF28 is provided in Figure 2. Moreover, the Q-Q plot of the GWAS during the discovery phase is presented in Figure 3. The SNP-based heritability for widespread pain was estimated to be 0.18, with a standard error of 0.02.

Figure 1.

The Manhattan plot of the GWAS analysis on widespread pain (N = 172,230).

Table 3.

The novel significant SNP and SNPs approaching significance identified by GWAS on widespread pain.

rsID	Chr	SNP position	Gene/Nearest gene	Discovery stage (UK Biobank)					Replication stage (Henkel et al. and FinnGen)
rsID	Chr	SNP position	Gene/Nearest gene	Effect allele	Non effective allele	Frequency the effect allele	Beta	p-Value	p-Value (beta) in FibroGene Cohort	p-Value (beta) in Twins UK cohort	p-Value (beta) in Finn Gen
rs34691025	5	73180281	ARHGEF28	A	G	0.798	0.004	1.76 × 10⁻⁸	0.34 (0.06)	0.90 (0.002)	0.79 (0.01)
rs11959702	5	73178796	ARHGEF28	G	T	0.798	0.004	1.77 × 10⁻⁸	0.35 (0.06)	0.87 (0.002)	0.79 (0.01)
rs377632600	5	73180608	ARHGEF28	T	G	0.801	0.004	2.74 × 10⁻⁸	0.34 (0.10)	0.91 (0.002)	0.78 (0.01)
rs729272	5	73185504	ARHGEF28	A	G	0.798	0.004	2.89 × 10⁻⁸	0.55 (0.04)	0.87 (−0.002)	0.79 (0.01)
rs6453030	5	73187503	ARHGEF28	G	C	0.798	0.004	3.30 × 10⁻⁸	0.57 (0.04)	0.82 (0.003)	0.85 (0.01)
rs117537726	18	10297679	AP005209.2	C	A	0.968	−0.008	6.05 × 10⁻⁸	0.55 (−0.09)	0.53 (−0.017)	0.39 (−0.07)
rs17387024	1	50993178	FAF1	T	C	0.852	−0.004	7.84 × 10⁻⁸	0.86 (−0.01)	0.85 (−0.003)	0.69 (0.02)
rs72831687	6	16092360	AL021407.2	G	A	0.992	−0.016	1.61 × 10⁻⁷	0.51 (0.21)	0.87 (0.002)	0.72 (0.04)
rs761429057	3	129474948	TMCC1	AAAAAAAAAAAG	A	0.953	−0.008	3.01 × 10⁻⁷	0.36 (−2.29)	0.015 (0.044)	0.47 (0.52)
rs538146040	15	93131238	AC091544.2	G	A	0.993	−0.018	3.58 × 10⁻⁷	0.87 (0.07)	0.61 (0.006)	0.70 (−0.25)
rs206969	12	120835893	AC003982.1	C	T	0.774	−0.004	4.00 × 10⁻⁷	0.80 (0.02)	0.91 (−0.001)	0.18 (−0.05)
rs113280333	3	69242581	FRMD4B	C	T	0.983	−0.011	4.53 × 10⁻⁷	0.78 (0.06)	0.34 (0.0016)	0.26 (−0.07)
rs148472276	2	3215642	EIPR1	C	A	0.994	−0.018	5.31 × 10⁻⁷	0.41 (−0.29)	0.56 (0.014)	0.58 (0.17)
rs4782029	16	17248016	XYLT1	A	C	0.087	0.005	7.04 × 10⁻⁷	0.17 (−0.12)	0.46 (0.012)	0.23 (0.06)

Chr: chromosome.

Figure 2.

The regional plots of locus in ARHGEF28 region.

Figure 3.

The Q-Q plot of the GWAS analysis on widespread pain (N = 172,230).

We tested the significant and independent SNPs identified in our discovery analysis, as well as other SNPs proximal to the genome-wide significance threshold, for replication in the FinnGen, FibroGene, and Twins UK cohorts. Within the TwinsUK cohort, specifically regarding CWP, we observed weak replication for rs761429057 in the TMCC1 gene (p = 0.017) among the 14 SNPs from the discovery cohort. However, in the context of fibromyalgia within the FibroGene and FinnGen cohorts, none of the 14 SNPs demonstrated replication in either cohort (p > 0.05). The p-values of the associations of these 14 independent SNPs from the discovery stage were extracted from FinnGen, FibroGene, and Twins UK cohorts are presented in Table 3.

In the sex-stratified GWAS analyses, the male-specific GWAS did not identify any single locus of genome-wide significance associated with widespread pain. Conversely, the female-specific GWAS revealed one significant locus associated with widespread pain, which differed from the findings of the primary GWAS. This novel locus was identified within the LRMDA gene on chromosome 10, with a p = 4.04 × 10⁻⁸ for rs137998089. Detailed information regarding these findings is provided in Table 4, and the corresponding Manhattan plot and regional plot are depicted in Figures 4 and 5, respectively.

Table 4.

The locus and top SNP identified by the female-specific GWAS on widespread pain.

Locus rank	rsID	Chr	SNP position	Nearest gene	Effect allele	Non effective allele	Frequency the effect allele	p value	Beta
1	rs137998089	10	78115898	LRMDA	A	G	0.994	4.03 × 10⁻⁸	−0.003

Chr: chromosome.

Figure 4.

The Manhattan plot of the female-specific GWAS analysis on widespread pain (N = 90,117).

Figure 5.

The regional plots of locus in LRMDA region.

Gene, gene-set, and tissue expression analysis

In our gene analysis, we mapped SNPs within genes to a total of 18,115 protein-coding genes. As a result, we established a genome-wide significance threshold of p < 0.05/18,115 = 2.76 × 10⁻⁶. At this stage, no genes met the criteria for significance. In our analysis of gene sets, we evaluated 17,109 gene sets in total. The threshold for statistical significance on a genome-wide scale was thus set at p < 0.05/17,109 = 2.92 × 10⁻⁶. However, none of the gene sets reached statistical significance. The top ten gene sets identified in this analysis are listed in Supplemental Table 3.

In the tissue expression analysis, we identified 30 general tissue types and 53 specific tissue types. However, none of the tissues exhibited a significant association. Further details and visualization of these results are available in Supplemental Figures 1 and 2.

In the gene expression heatmaps presented in Figure 6, the gene ARHGEF28 exhibits higher expression levels in the nerve, oral, and adrenal gland tissues compared to other tissue types. Tissue specificity testing did not yield significant differences in the DEG across both the 54 and 30 tissue types analyzed in GTEx version 8. Detailed visualizations of these findings are included in Supplemental Figure 3.

Figure 6.

Gene expression heatmaps across 54 tissue types and 30 general tissue Categories from GTEx Version 8.

Genetic correlation analysis by LDSC

In our study examining the relationship between widespread pain and various pain phenotypes, as well as the genetic correlation between widespread pain in females and males, we discovered several significant genetic correlations, as depicted in Figure 7. After correcting for multiple testing, the three most notable genetic correlations ( $r_{g}$ ) were found in the following categories: other joint disorders, not elsewhere classified ( $r_{g}$ = 0.94, p = 6.20 × 10⁻⁶), spondylosis ( $r_{g}$ = 0.93, p = 7.73 × 10⁻⁶), and other diseases of the anus and rectum ( $r_{g}$ = 0.89, p = 1.95 × 10⁻⁶). The genetic correlation between widespread pain in females and males was determined to be $r_{g}$ = 0.80 (p = 7.36 × 10⁻⁹²). The comprehensive results are provided in Supplemental Table 4.

Figure 7.

Genetic correlation results for widespread pain using LDSC on CTG-VL. For a full list, see Supplemental Table 4.

Expression quantitative trait loci (eQTL), chromatin interaction analysis, and positional mapping

From the cis-eQTL analysis, SNPs rs35588745 and rs71636094 were found to be significantly associated with the tissue type Cells_Cultured_fibroblasts, achieving a False Discovery Rate (FDR) of less than 0.05. The results of the chromatin interactions and cis-eQTL analyses are illustrated in the circus plot presented in Figure 8. This plot shows the genes that interact with ARHGEF28, including IQGAP2, GCTN4, FAM169A, NSA2, GFM2, HEXB, ENC1, UTP14, and AC008387.1. Additionally, the positional mapping analysis did not identify significant results.

Figure 8.

Circos plot illustrating chromatin interactions and eQTLs on chromosome 5.

PheWAS

PheWAS were conducted using the GWAS ATLAS platform to explore phenotypes associated with independently significant SNPs (rs34691025) and significantly associated genes (ARHGEF28). SNP rs34691025 exhibited significant associations with the frequency of tenseness/restlessness in the last 2 weeks (p = 1.5 × 10⁻⁴), the mode of anisotropy in the retrolenticular part of the internal capsule (p = 1.70 × 10⁻⁶), and long-standing illness, disability, or infirmity (p = 3.60 × 10⁻¹²). The gene ARHGEF28 was associated with hearing difficulty/problems (p = 1.09 × 10⁻²³), high blood pressure (p = 3.42 × 10⁻⁷), and hypertension (p = 2.26 × 10⁻⁵). The gene FAF1 showed an association with standing height (p = 4.02 × 10⁻⁴⁵) and was also specifically associated with headache (p = 7.51 × 10⁻⁴) and stomachache (p = 1.14 × 10⁻³). All significant traits that passed the Bonferroni correction are listed in Supplemental Table 5. The plots illustrating phenotypes associated with the SNP and gene are visualized in Figures 9 to 11, respectively.

Figure 9.

PheWAS results for SNP rs34691025 within primary widespread pain GWAS.

Figure 10.

PheWAS results for gene ARHGEF28 within primary widespread pain GWAS.

Figure 11.

PheWAS results for gene FAF1 within primary widespread pain GWAS.

Mendelian randomization of widespread pain and hearing loss

The IVW estimation did not reveal a significant genetically causal relationship between hearing loss and widespread pain (OR = 0.99, 95% CI 0.99–1.01, p = 0.13). There was no evidence of statistical horizontal pleiotropy (MR Egger p = 0.82) or significant heterogeneity (IVW Q = 0.94, p = 0.92). Additionally, reverse causation was not detected (IVW p = 0.82). Detailed results and figures are presented in Supplemental Figure 4 and Supplemental Table 6.

Discussion

In this GWAS of widespread pain, conducted using the comprehensive UK Biobank dataset, we identified one novel significant genetic locus and several loci that marginally reached the genome-wide significance level associated with self-reported widespread pain during the past month. In the sex-stratified GWAS, we observed differences in the prevalence of genetic variants of the widespread pain between males and females. Furthermore, our post-GWAS analysis revealed significant genetic correlations between widespread pain and other phenotypes. Additionally, a PheWAS was conducted to further explore the genetic mechanisms underlying widespread pain. In this study, we utilized the generic pain questionnaire from the UK Biobank, a valuable screening instrument, to investigate the potential genetic underpinnings of heterogeneous pain phenotypes, including widespread pain. Employing the UK Biobank offers researchers the advantage of mitigating potential issues related to reduced power due to phenotype heterogeneity. This is achieved by leveraging substantial sample sizes, thereby enhancing the clarity of statistical results amidst potential noise.

In our primary GWAS, we have analyzed and subsequently identified one novel locus that was associated with widespread pain. This significant locus was situated within the ARHGEF28 gene on chromosome 5q13.2. This particular locus exhibited a significant p-value of 1.76 × 10⁻⁸ for the marker rs34691025. This locus spans an extensive 11 kb within the intronic region of the ARHGEF28 gene and encompasses as many as 9 genome-wide significant SNPs, as detailed in Supplemental Table 2. The ARHGEF28 gene (Rho Guanine Nucleotide Exchange Factor 28) is a member of the Rho guanine nucleotide exchange factor family. It plays a pivotal role in regulating Rho GTPases involved in cytoskeletal reorganization and cell migration, may influence pain perception and sensitivity through its effects on neuronal plasticity and inflammatory processes. The cytoskeletal dynamics in neurons, particularly in nociceptive pathways, play a vital role in the development and maintenance of pain states, with alterations in these processes leading to changes in the structural and functional plasticity of pain pathways, thereby impacting pain perception.^21,22 Furthermore, ARHGEF28’s involvement in immune cell signaling pathways suggests a mechanism by which genetic predispositions could influence the development of chronic pain through immune-mediated inflammation.^23,24 The link between the immune system and pain has been increasingly recognized, with several studies highlighting how immune responses contribute to the sensitization of the central nervous system, leading to pain conditions.^25,26 The role of ARHGEF28 in modulating immune cell functions, including those of mast cells and macrophages, which are key players in inflammation and pain, further supports its potential involvement in pain phenotypes. This dual role in neuronal and immune cell function underscores the complexity of the genetic underpinnings of pain and the need for multidisciplinary approaches to fully understand these mechanisms. The association of ARHGEF28 with widespread pain underscores the importance of genetic factors in the development of pain disorders and highlights new potential pathways for therapeutic intervention. Further research is needed to explore the functional impact of ARHGEF28 variations on neuronal and immune cell functions and to understand how these effects contribute to the manifestation of widespread pain. A study identified a heterozygous mutation in amyotrophic lateral sclerosis patients that is predicted to generate a premature truncated ARHGEF28 gene product, which could provide a possible link between this gene and widespread pain.²⁷ Moreover, a previous GWAS research has reported a significant association between the ARHGEF28 gene and hearing problems, suggesting a potential new pathway linking widespread pain and auditory issues.²⁸

We further identified a cluster of SNPs within the FAF1 gene on chromosome 1, displaying marginal genome-wide significance, with rs17387024 exhibiting a minimal p-value of 7.87 × 10⁻⁸. The FAF1 gene, formally known as Fas Associated Factor 1, is a protein-coding gene. The association between the FAF1 gene and widespread pain may be attributed to its involvement in various biological processes relevant to pain perception and modulation. FAF1 has been implicated in apoptosis, inflammation, and neuronal signaling, all of which are crucial in the pathogenesis of pain.²⁹ Studies have shown that FAF1 is involved in the NF-κB signaling pathway, which plays a significant role in inflammation and pain.^30,31 Additionally, FAF1 has been associated with neuronal death and dysfunction, suggesting a potential link to neuropathic pain.³² The interaction of FAF1 with other proteins involved in pain pathways, such as TRPV1, further supports its role in pain perception.³³ Several previous GWAS have identified the FAF1 gene as being associated with multisite chronic pain and general pain, providing a possible basis for its link with widespread pain.^34–36 However, further research is still necessary to fully elucidate the exact mechanism by which FAF1 influences widespread pain.

In the sex-stratified GWAS, we discovered a novel locus unique to the female-specific analysis, located within the LRMDA gene on chromosome 10. The top SNP, rs5779595, exhibited a significant p-value of 4.04 × 10⁻⁸. The LRMDA gene, also known as leucine-rich melanocyte differentiation-associated protein, has been implicated in some biological processes that could underlie its relationship with pain. Notably, LRMDA is involved in melanocyte differentiation and pigmentation, which might suggest modulating nociceptive pathways.³⁷ Furthermore, several GWAS studies have shown that the LRMDA gene is associated with diabetes, a condition that can potentially cause widespread pain and fibromyalgia.^38,39 This association provides a potential link between the LRMDA gene and widespread pain, suggesting that the gene may play a role in the underlying mechanisms of the widespread pain condition. The genetic correlation for widespread pain between males and females was also calculated ( $r_{g}$ = 0.80), indicating a potential disparity in the genetic underpinnings of widespread pain across sexes. This divergence is likely attributable to the significant influence of sex differences on the heterogeneity of pain.⁴⁰ These findings highlight the necessity of accounting for biological differences between males and females in research on the genetics of pain.

Building on our research findings, there are significant genetic correlations between widespread pain and several other phenotypes, suggesting that shared genetic factors underpinning these conditions. The top three phenotypes most genetically correlated with widespread pain were other joint disorders, not elsewhere classified ( $r_{g}$ = 0.94, p = 6.20 × 10⁻⁶), spondylosis ( $r_{g}$ = 0.93, p = 7.73 × 10⁻⁶), and other diseases of the anus and rectum ( $r_{g}$ = 0.89, p = 1.95 × 10⁻⁶). Our study provides evidence of the complex interplay of genetic determinants in the manifestation of pain, supporting previous research findings. Specifically, several prior analyses have identified an association between CWP and joint diseases, and some studies suggest that fibromyalgia may share common mechanisms with disorders of the anus and rectum.^41,42 In the PheWAS, the results reveal that the top SNP (rs34691025) and corresponding gene (ARHGEF28) are significantly associated with hearing disorders and cardiovascular diseases. In our two-sample MR analysis, we found no significant bidirectional causal effect between hearing loss and widespread pain, suggesting that the identified loci may directly influence widespread pain rather than acting through hearing loss.

In our research, we employed a novel operational definition to delineate cases of widespread pain, characterizing them as individuals experiencing pain throughout their body during the last month. This approach stands in contrast to previous GWAS that defined cases based on the presence of pain all over the body lasting for more than 3 months.⁹ Notably, the CWP cases in these studies numbered 6914, a criterion that differs from our own. The adoption of this alternative definition was driven by the consideration that a one-month duration might capture a broader spectrum of pain phenotypes, potentially including both acute and chronic manifestations. By focusing on a 1-month timeframe, we intend to provide a distinct and meaningful delineation of widespread pain cases, offering a new perspective on the genetic underpinnings of pain phenotypes. Our choice of the definition is motivated by the goal of enhancing the understanding of genetic factors associated with widespread pain.

Our research provides initial insights into the genetic underpinnings of widespread pain. However, it encounters certain limitations, particularly in the validation phase of our results. During the replication phase, we utilized three public datasets of CWP and fibromyalgia, yet we did not achieve significant replication results for our study. This outcome could potentially be attributed to differences in the definitions of widespread pain, CWP, and fibromyalgia. Although widespread pain and fibromyalgia syndrome are sometimes used interchangeably, fibromyalgia is generally more severe and includes symptoms like sleep disturbance, fatigue, and depression, indicating it is a more severe subset of widespread pain.^43,44 Additionally, the limited number of participants in the replication cohorts might have contributed to this challenge. In the female-specific GWAS, rs5779595 becomes the top and only SNP identified within the locus, pointing to the potential spuriousness of this association and underscoring the need for further exploration and study. Furthermore, it is essential to discuss the methodological approach adopted in our study. We characterized widespread pain cases and controls based on the responses provided by UK Biobank participants to a specific question regarding pain experienced over the last month. However, a critical observation is that this question did not delve into nuanced details such as the severity or frequency of the pain. Consequently, the phenotype definition derived from these responses should be perceived and interpreted as being broadly defined. This broad definition may have implications for the specificity and sensitivity of our findings, as it potentially includes a heterogeneous group of individuals with varying pain experiences. Future research may need to incorporate more detailed questionnaires and larger sample sizes to capture the complexity of widespread pain and provide a more comprehensive understanding of its genetic basis. Additionally, exploring alternative phenotype definitions that account for pain severity and frequency could enhance the precision of genetic associations and contribute to a deeper understanding of the mechanisms underlying widespread pain.

Conclusion

In conclusion, our primary GWAS of widespread pain, utilizing the extensive UK Biobank dataset, successfully identified one novel genetic locus and several loci that reached genome-wide significance. The sex-stratified GWAS outputs revealed sex-specific variants associated with of widespread pain between males and females. Furthermore, our post-GWAS analysis demonstrated substantial genetic correlations and shared mechanisms between widespread pain and other phenotypes. This study advances our understanding of the genetic factors contributing to widespread pain.

Supplemental Material

sj-pdf-1-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-1-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-10-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-10-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-2-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-2-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-3-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-3-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-4-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-4-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-5-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-5-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-6-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-6-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-7-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-7-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-8-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-8-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Supplemental Material

sj-pdf-9-mpx-10.1177_17448069251346603 – Supplemental material for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230)

Supplemental material, sj-pdf-9-mpx-10.1177_17448069251346603 for Genome-wide association study identifies novel genetic variants associated with widespread pain in the UK Biobank (N = 172,230) by Qi Pan, Tengda Cai, Yiwen Tao, Luning Yang, Roger Compte, Maryam Kazemi Naeini, Mainul Haque, Tania Dottorini, Frances MK Williams and Weihua Meng in Molecular Pain

Footnotes

Acknowledgements

The authors extend their sincere gratitude to the participants of the UK Biobank, FinnGen, FibroGene, and TwinsUK cohorts for their invaluable contributions of genetic and phenotypic data. This study will adhere to all ethical guidelines and data protection protocols of the UK Biobank. This research has been conducted using the UK Biobank Resource under Application Number 89386.

Author contributions

QP drafted the paper and performed the UK Biobank GWAS analysis. TC, YT, and LY contributed to data formatting. MH and TD provided comments to the paper. RC, MKN, and FMKW conducted the replication study and provided comments to the paper. WM organized the project and provided comments.

Availability of data and material

The summary statistics of the UK Biobank results on widespread pain can be accessed upon publication. Any other data relevant to the study that are not included in the article or its supplementary materials are available from the authors upon reasonable request.

Declaration of conflicting interests

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

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This study was mainly funded by the Pioneer and Leading Goose R&D Program of Zhejiang Province 2023 with reference number 2023C04049 and Ningbo International Collaboration Program 2023 with reference number 2023H025. TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy’s and St Thomas’ NHS Foundation Trust in partnership with King’s College London.

Ethics approval and consent to participate

This study was approved by the Ethics Committee of the University of Nottingham, Ningbo, China. All authors have consent for participation. The TwinsUK Study was approved by London-Westminster Research Ethics Committee (REC referenceEC04/015), and Guy’s and St Thomas’ NHS Foundation Trust Research and Development (R&D). The TwinsUK Biobank was approved by the HRA – Liverpool East Research Ethics Committee (REC reference 19/NW/0187), IRAS ID 258513. All participants provide written, informed consent.

ORCID iDs

Qi Pan

Weihua Meng

Supplemental material

Supplemental material for this article is available online.

References

Butler

Landmark

Glette

Borchgrevink

Woodhouse

. Chronic widespread pain—the need for a standard definition. Pain 2016; 157(3): 541–543.

Wolfe

Butler

Fitzcharles

Häuser

Katz

Mease

Rasker

Russell

Walitt

. Revised chronic widespread pain criteria: development from and integration with fibromyalgia criteria. Scand J Pain 2020; 20(1): 77–86.

Fayaz

Croft

Langford

Donaldson

Jones

. Prevalence of chronic pain in the UK: a systematic review and meta-analysis of population studies. BMJ Open 2016; 6(6): e010364.

Mansfield

Sim

Jordan

. A systematic review and meta-analysis of the prevalence of chronic widespread pain in the general population. Pain 2016; 157(1): 55–64.

Creed

. A review of the incidence and risk factors for fibromyalgia and chronic widespread pain in population-based studies. Pain 2020; 161(6): 1169–1176.

Schaefer

Mann

Masters

Cappelleri

Daniel

Zlateva

McElroy

Chandran

Adams

Assaf

McNett

Mease

Silverman

Staud

. The comparative burden of chronic widespread pain and fibromyalgia in the United States. Pain Pract 2016; 16(5): 565–579.

Kato

Sullivan

Evengård

Pedersen

. Importance of genetic influences on chronic widespread pain. Arthritis Rheum 2006; 54(5): 1682–1686.

Peters

Broer

Willemen

Eiriksdottir

Hocking

Holliday

Horan

Meulenbelt

Neogi

Popham

Schmidt

Soni

Valdes

Amin

Dennison

Eijkelkamp

Harris

Hart

Hofman

Huygen

Jameson

Jones

Launer

Kerkhof

de Kruijf

McBeth

Kloppenburg

Ollier

Oostra

Payton

Rivadeneira

Smith

Stolk

Teumer

Thomson

Uitterlinden

Wang

van Wingerden

Arden

Cooper

Felson

Gudnason

Macfarlane

Pendleton

Slagboom

Spector

Völzke

Kavelaars

van Duijn

Williams

van Meurs

. Genome-wide association study meta-analysis of chronic widespread pain: evidence for involvement of the 5p15.2 region. Ann Rheum Dis 2013; 72(3): 427–436.

Rahman

Winsvold

Chavez Chavez

Børte

Tsepilov

Sharapov

; HUNT All-In Pain; Aulchenko

Hagen

Fors

Hveem

Zwart

van Meurs

Freidin

Williams

. Genome-wide association study identifies RNF123 locus as associated with chronic widespread musculoskeletal pain. Ann Rheum Dis 2021; 80(9): 1227–1235.

10.

Bycroft

Freeman

Petkova

Band

Elliott

Sharp

Motyer

Vukcevic

Delaneau

O’Connell

Cortes

Welsh

Young

Effingham

McVean

Leslie

Allen

Donnelly

Marchini

. The UK Biobank resource with deep phenotyping and genomic data. Nature 2018; 562(7726): 203–209.

11.

Kurki

Karjalainen

Palta

Sipilä

Kristiansson

Donner

Reeve

Laivuori

Aavikko

Kaunisto

Loukola

Lahtela

Mattsson

Laiho

Della Briotta Parolo

Lehisto

Kanai

Mars

Rämö

Kiiskinen

Heyne

Veerapen

Rüeger

Lemmelä

Zhou

Ruotsalainen

Pärn

Hiekkalinna

Koskelainen

Paajanen

Llorens

Gracia-Tabuenca

Siirtola

Reis

Elnahas

Sun

Foley

Aalto-Setälä

Alasoo

Arvas

Auro

Biswas

Bizaki-Vallaskangas

Carpen

Chen

Dada

Ding

Ehm

Eklund

Färkkilä

Finucane

Ganna

Ghazal

Graham

Green

Hakanen

Hautalahti

Hedman

ÅK

Hiltunen

Hinttala

Hovatta

Huertas-Vazquez

Huilaja

Hunkapiller

Jacob

Jensen

Joensuu

John

Julkunen

Jung

Junttila

Kaarniranta

Kähönen

Kajanne

Kallio

Kälviäinen

Kaprio

; FinnGen; Kerimov

Kettunen

Kilpeläinen

Kilpi

Klinger

Kosma

Kuopio

Kurra

Laisk

Laukkanen

Lawless

Liu

Longerich

Mägi

Mäkelä

Mäkitie

Malarstig

Mannermaa

Maranville

Matakidou

Meretoja

Mozaffari

Niemi

MEK

Niemi

Niiranen

O Donnell

Obeidat

Okafo

Ollila

Palomäki

Palotie

Partanen

Paul

Pelkonen

Pendergrass

Petrovski

Pitkäranta

Platt

Pulford

Punkka

Pussinen

Raghavan

Rahimov

Rajpal

Renaud

Riley-Gillis

Rodosthenous

Saarentaus

Salminen

Salomaa

Schleutker

Serpi

Shen

Siegel

Silander

Siltanen

Soini

Soininen

Sul

Tachmazidou

Tasanen

Tienari

Toppila-Salmi

Tukiainen

Tuomi

Turunen

Ulirsch

Vaura

Virolainen

Waring

Waterworth

Yang

Nelis

Reigo

Metspalu

Milani

Esko

Fox

Havulinna

Perola

Ripatti

Jalanko

Laitinen

Mäkelä

Plenge

McCarthy

Runz

Daly

Palotie

. FinnGen provides genetic insights from a well-phenotyped isolated population. Nature 2023; 613(7944): 508–518.

12.

Eslam

Hashem

Romero-Gomez

Berg

Dore

Mangia

Chan

HLY

Irving

Sheridan

Abate

Adams

Weltman

Bugianesi

Spengler

Shaker

Fischer

Mollison

Cheng

Nattermann

Riordan

Miele

Kelaeng

Ampuero

Ahlenstiel

McLeod

Powell

Liddle

Douglas

Booth

George

; International Liver Disease Genetics Consortium (ILDGC). FibroGENE: a gene-based model for staging liver fibrosis. J Hepatol 2016; 64(2): 390–398.

13.

Moayyeri

Hammond

Hart

Spector

. The UK Adult Twin Registry (TwinsUK Resource). Twin Res Hum Genet 2013; 16(1): 144–149.

14.

Yang

Lee

Goddard

Visscher

. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet 2011; 88(1): 76–82.

15.

Watanabe

Taskesen

van Bochoven

Posthuma

. Functional mapping and annotation of genetic associations with FUMA. Nat Commun 2017; 8(1): 1826.

16.

Pruim

Welch

Sanna

Teslovich

Chines

Gliedt

Boehnke

Abecasis

Willer

. LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 2010; 26(18): 2336–2337.

17.

Bulik-Sullivan

Loh

Finucane

Ripke

Yang

Patterson

Daly

Price

Neale

. LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet 2015; 47(3): 291–295.

18.

Võsa

Claringbould

Westra

Bonder

Deelen

Zeng

Kirsten

Saha

Kreuzhuber

Yazar

Brugge

Oelen

de Vries

van der Wijst

MGP

Kasela

Pervjakova

Alves

Favé

Agbessi

Christiansen

Jansen

Seppälä

Tong

Teumer

Schramm

Hemani

Verlouw

Yaghootkar

Sönmez Flitman

Brown

Kukushkina

Kalnapenkis

Rüeger

Porcu

Kronberg

Kettunen

Lee

Zhang

Hernandez

Arindrarto

Beutner

; BIOS Consortium; i2QTL Consortium; Dmitrieva

Elansary

Fairfax

Georges

Heijmans

Hewitt

Kähönen

Kim

Knight

Kovacs

Krohn

Loeffler

Marigorta

Mei

Momozawa

Müller-Nurasyid

Nauck

Nivard

Penninx

BWJH

Pritchard

Raitakari

Rotzschke

Slagboom

Stehouwer

CDA

Stumvoll

Sullivan

’t Hoen

PAC

Thiery

Tönjes

van Dongen

van Iterson

Veldink

Völker

Warmerdam

Wijmenga

Swertz

Andiappan

Montgomery

Ripatti

Perola

Kutalik

Dermitzakis

Bergmann

Frayling

van Meurs

Prokisch

Ahsan

Pierce

Lehtimäki

Boomsma

Psaty

Gharib

Awadalla

Milani

Ouwehand

Downes

Stegle

Battle

Visscher

Yang

Scholz

Powell

Gibson

Esko

Franke

. Large-scale cis- and trans-eQTL analyses identify thousands of genetic loci and polygenic scores that regulate blood gene expression. Nat Genet 2021; 53(9): 1300–1310.

19.

Sun

Chang

Cai

Hong

Zhou

. Chromatin interaction analysis with updated ChIA-PET Tool (V3). Genes (Basel) 2019; 10(7): 554.

20.

Watanabe

Stringer

Frei

Umićević Mirkov

de Leeuw

Polderman

TJC

van der Sluis

Andreassen

Neale

Posthuma

. A global overview of pleiotropy and genetic architecture in complex traits. Nat Genet 2019; 51(9): 1339–1348.

21.

St John Smith

. Advances in understanding nociception and neuropathic pain. J Neurol 2018; 265(2): 231–238.

22.

Sit

Manser

. Rho GTPases and their role in organizing the actin cytoskeleton. J Cell Sci 2011; 124(Pt 5): 679–683.

23.

Song

Lin

Huang

Yao

Wang

. Rare, low-frequency and common coding variants of ARHGEF28 gene and their association with sporadic amyotrophic lateral sclerosis. Neurobiol Aging 2020; 87: 138.e1–138.e6.

24.

Malafoglia

Ilari

Vitiello

Tenti

Balzani

Muscoli

Raffaeli

Bonci

. The interplay between chronic pain, opioids, and the immune system. Neuroscientist 2022; 28(6): 613–627.

25.

Raoof

Willemen

HLDM

Eijkelkamp

. Divergent roles of immune cells and their mediators in pain. Rheumatology 2017; 57(3): 429–440.

26.

Vergne-Salle

Bertin

. Chronic pain and neuroinflammation. Joint Bone Spine 2021; 88(6): 105222.

27.

Droppelmann

Wang

Campos-Melo

Keller

Volkening

Hegele

Strong

. Detection of a novel frameshift mutation and regions with homozygosis within ARHGEF28 gene in familial amyotrophic lateral sclerosis. Amyotroph Lateral Scler Frontotemporal Degener 2013; 14(5–6): 444–451.

28.

Ivarsdottir

Holm

Benonisdottir

Olafsdottir

Sveinbjornsson

Thorleifsson

Eggertsson

Halldorsson

Hjorleifsson

Melsted

Gylfason

Arnadottir

Oddsson

Jensson

Jonasdottir

Juliusdottir

Stefansdottir

Tragante

Halldorsson

Petersen

Thorgeirsson

Thorsteinsdottir

Sulem

Hinriksdottir

Jonsdottir

Gudbjartsson

Stefansson

. The genetic architecture of age-related hearing impairment revealed by genome-wide association analysis. Commun Biol 2021; 4(1): 706.

29.

Ryu

Lee

Park

Jun

Jung

Kim

. Fas-associated factor 1, FAF1, is a member of Fas death-inducing signaling complex. J Biol Chem 2003; 278(26): 24003–24010.

30.

Ahmed

Berg

Alkass

Druid

Hart

Svensson

Kosek

. NF-κB-associated pain-related neuropeptide expression in patients with degenerative disc disease. Int J Mol Sci 2019; 20(3): 658.

31.

Park

Jang

Lee

Kim

. Fas-associated factor-1 inhibits nuclear factor-κB (NF-κB) activity by interfering with nuclear translocation of the RelA (p65) subunit of NF-κB. J Biol Chem 2004; 279(4): 2544–2549.

32.

Park

Kim

B-S

Kim

. A novel function of FAF1, which induces dopaminergic neuronal death through cell-to-cell transmission. Cell Commun Signal 2020; 18(1): 133.

33.

Kim

Kang

Shin

Hwang

Yang

Shim

Park

Kim

Cho

Shin

. TRPV1 recapitulates native capsaicin receptor in sensory neurons in association with Fas-associated factor 1. J Neurosci 2006; 26(9): 2403–2412.

34.

Johnston

KJA

Adams

Nicholl

Ward

Strawbridge

Ferguson

McIntosh

Bailey

MES

Smith

. Genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet 2019; 15(6): e1008164.

35.

Johnston

KJA

Ward

Ray

Adams

McIntosh

Smith

Strawbridge

Price

Smith

Nicholl

Bailey

MES

. Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet 2021; 17(4): e1009428.

36.

Mocci

Ward

Perry

Starkweather

Stone

Schabrun

Renn

Dorsey

Ament

. Genome wide association joint analysis reveals 99 risk loci for pain susceptibility and pleiotropic relationships with psychiatric, metabolic, and immunological traits. PLoS Genet 2023; 19(10): e1010977.

37.

Fernández

Hayashi

Garrido

Montero

Guardia

Suzuki

Montoliu

. Genetics of non-syndromic and syndromic oculocutaneous albinism in human and mouse. Pigment Cell Melanoma Res 2021; 34(4): 786–799.

38.

Vujkovic

Keaton

Lynch

Miller

Zhou

Tcheandjieu

Huffman

Assimes

Lorenz

Zhu

Hilliard

Judy

Huang

Lee

Klarin

Pyarajan

Danesh

Melander

Rasheed

Mallick

Hameed

Qureshi

Afzal

Malik

Jalal

Abbas

Sheng

Gao

Kaestner

Susztak

Sun

DuVall

Cho

Lee

Gaziano

Phillips

Meigs

Reaven

Wilson

Edwards

Rader

Damrauer

O’Donnell

Tsao

; HPAP Consortium; Regeneron Genetics Center; VA Million Veteran Program; Chang

Voight

Saleheen

. Discovery of 318 new risk loci for type 2 diabetes and related vascular outcomes among 1.4 million participants in a multi-ancestry meta-analysis. Nat Genet 2020; 52(7): 680–691.

39.

Ishigaki

Akiyama

Kanai

Takahashi

Kawakami

Sugishita

Sakaue

Matoba

Low

S-K

Okada

Terao

Amariuta

Gazal

Kochi

Horikoshi

Suzuki

Ito

Koyama

Ozaki

Niida

Sakata

Kohno

Shiraishi

Momozawa

Hirata

Matsuda

Ikeda

Iwata

Ikegawa

Kou

Tanaka

Nakagawa

Suzuki

Hirota

Tamari

Chayama

Miki

Mori

Nagayama

Daigo

Miki

Katagiri

Ogawa

Obara

Ito

Yoshida

Imoto

Takahashi

Tanikawa

Suzuki

Sinozaki

Minami

Yamaguchi

Asai

Takahashi

Yamaji

Takahashi

Fujioka

Takata

Yanai

Masumoto

Koretsune

Kutsumi

Higashiyama

Murayama

Minegishi

Suzuki

Tanno

Shimizu

Yamaji

Iwasaki

Sawada

Uemura

Tanaka

Naito

Sasaki

Wakai

Tsugane

Yamamoto

Murakami

Nakamura

Raychaudhuri

Inazawa

Yamauchi

Kadowaki

Kubo

Kamatani

. Large-scale genome-wide association study in a Japanese population identifies novel susceptibility loci across different diseases. Nat Genet 2020; 52(7): 669–679.

40.

Gao

Liu

Wang

Zhang

Zhu

. Sex differences in stress-induced hyperalgesia and its mechanisms. J Neurosci Res 2024; 102(1): e25266.

41.

Mayer

Towns

Neblett

Theodore

Gatchel

. Chronic widespread pain in patients with occupational spinal disorders: prevalence, psychiatric comorbidity, and association with outcomes. Spine 2008; 33(17): 1889–1897.

42.

Caldarella

Giamberardino

Sacco

Affaitati

Milano

Lerza

Balatsinou

Laterza

Pierdomenico

Cuccurullo

Neri

. Sensitivity disturbances in patients with irritable bowel syndrome and fibromyalgia. Am J Gastroenterol 2006; 101(12): 2782–2789.

43.

Wolfe

Clauw

Fitzcharles

Goldenberg

Häuser

Katz

Mease

Russell

Walitt

. 2016 Revisions to the 2010/2011 fibromyalgia diagnostic criteria. Semin Arthritis Rheum 2016; 46(3): 319–329.

44.

Shipley

. Chronic widespread pain and fibromyalgia syndrome. Medicine 2014; 42(5): 271–274.

Supplementary Material

Please find the following supplemental material available below.

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

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

0.00 MB

0.51 MB

0.03 MB

0.19 MB

0.74 MB

0.51 MB

0.06 MB

0.03 MB

0.02 MB

0.33 MB

0.70 MB