#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank


Authors: Keira J. A. Johnston aff001;  Joey Ward aff001;  Pradipta R. Ray aff004;  Mark J. Adams aff002;  Andrew M. McIntosh aff002;  Blair H. Smith aff005;  Rona J. Strawbridge aff001;  Theodore J. Price aff004;  Daniel J. Smith aff001;  Barbara I. Nicholl aff001;  Mark E. S. Bailey aff003
Authors place of work: Institute of Health and Wellbeing, University of Glasgow, Glasgow, Scotland, United Kingdom aff001;  Division of Psychiatry, University of Edinburgh, Edinburgh, Scotland, United Kingdom aff002;  School of Life Sciences, College of Medical, Veterinary & Life Sciences, University of Glasgow, Glasgow, Scotland, United Kingdom aff003;  School of Behavioral and Brain Sciences, The University of Texas at Dallas, Richardson, Texas, United States of America aff004;  Division of Population Health Sciences, University of Dundee, Ninewells Hospital and Medical School, Dundee, Scotland, United Kingdom aff005;  Department of Medicine Solna, Karolinska Institute, Stockholm, Sweden aff006
Published in the journal: Sex-stratified genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet 17(4): e1009428. doi:10.1371/journal.pgen.1009428
Category: Research Article
doi: https://doi.org/10.1371/journal.pgen.1009428

Summary

Chronic pain is highly prevalent worldwide and imparts a significant socioeconomic and public health burden. Factors influencing susceptibility to, and mechanisms of, chronic pain development, are not fully understood, but sex is thought to play a significant role, and chronic pain is more prevalent in women than in men. To investigate sex differences in chronic pain, we carried out a sex-stratified genome-wide association study of Multisite Chronic Pain (MCP), a derived chronic pain phenotype, in UK Biobank on 178,556 men and 209,093 women, as well as investigating sex-specific genetic correlations with a range of psychiatric, autoimmune and anthropometric phenotypes and the relationship between sex-specific polygenic risk scores for MCP and chronic widespread pain. We also assessed whether MCP-associated genes showed expression pattern enrichment across tissues. A total of 123 SNPs at five independent loci were significantly associated with MCP in men. In women, a total of 286 genome-wide significant SNPs at ten independent loci were discovered. Meta-analysis of sex-stratified GWAS outputs revealed a further 87 independent associated SNPs. Gene-level analyses revealed sex-specific MCP associations, with 31 genes significantly associated in females, 37 genes associated in males, and a single gene, DCC, associated in both sexes. We found evidence for sex-specific pleiotropy and risk for MCP was found to be associated with chronic widespread pain in a sex-differential manner. Male and female MCP were highly genetically correlated, but at an rg of significantly less than 1 (0.92). All 37 male MCP-associated genes and all but one of 31 female MCP-associated genes were found to be expressed in the dorsal root ganglion, and there was a degree of enrichment for expression in sex-specific tissues. Overall, the findings indicate that sex differences in chronic pain exist at the SNP, gene and transcript abundance level, and highlight possible sex-specific pleiotropy for MCP. Results support the proposition of a strong central nervous-system component to chronic pain in both sexes, additionally highlighting a potential role for the DRG and nociception.

Keywords:

Gene expression – Genetic loci – Genetics – Genome-wide association studies – Human genetics – Metaanalysis – Pain – Single nucleotide polymorphisms

Introduction

Chronic pain is widely defined as pain persisting beyond 3 months [1,2], and can be a primary disorder [3] or secondarily associated with injury, surgery or a range of medical conditions. Chronic pain is highly prevalent worldwide [49] and imparts a significant socioeconomic and public health burden [10]. Factors influencing susceptibility to chronic pain, and the mechanisms underlying its development and maintenance, are not fully understood.

Several aspects of chronic pain including Chronic Pain Grade [11], severe chronic pain and low back pain have been studied from a genetic perspective and found to be complex traits. Heritability estimates vary from ~30–46% in twin, pedigree and factor analysis studies [1215], while single nucleotide polymorphism (SNP) heritability has been estimated from genome-wide association studies (GWAS) to be ~7–10% [16,17].

It is increasingly recognised that sex differences in many complex human traits are biologically important, with genetic architecture for many traits being to some extent sex-specific [18], and a ‘sex-aware’ approach to genetic analysis has been widely advocated [19]. Sex as a biological variable has wide-ranging effects on the functioning of the genome and on resultant phenotypes. These effects can be mediated via sex-differential gene expression [20,21], sex differences in methylation [2226] and expression quantitative trait locus (eQTL) effects [27,28], or differing levels and actions of hormones [29,30]. Sex can also influence traits through environmental factors strongly correlated with sex [23,24] and sex-specific pleiotropy [25,31]. Chronic pain exhibits sex-related prevalence differences, and is more common in women than in men [3234]. There are also potential sex differences in the impact of pain on functioning in daily life, and in the success of specific coping strategies [35]. In addition to differences in prevalence between the sexes, sex differences in underlying pain mechanisms and their modulation by immune cells have been recently reported [3638], and immune responses in general can differ by sex [39].

Multisite Chronic Pain (MCP) is a derived chronic pain phenotype, defined as the sum of the number of sites of chronic pain on the body, here expressed on a scale from 0–7 [17]. We have previously shown in UK Biobank [17] that genetic predisposition to MCP (as captured by a polygenic risk score; PRS) was associated with Chronic Widespread Pain (CWP), a separate but related chronic pain phenotype, in women but not men [17]. Additional unpublished findings using the Generation Scotland study [40] demonstrated that the PRS-MCP was associated with both chronic pain grade and an MCP-like phenotype in both men and women, but that the magnitude of effect was roughly twice as great in women as it was in men. There was also a significant PRS-by-sex interaction. These findings suggest that sets of variants contributing to chronic pain in males and females may act differently, or have different genetic effect sizes, in the two sexes. Here we report on an exploration of this preliminary evidence for the existence of sex-specific loci associated with MCP using a sex-stratified GWAS analysis approach in UK Biobank, and identify several sex-specific MCP loci. A meta-analysis of the female- and male-specific GWASs also revealed novel MCP loci not identified in the original MCP GWAS. We have also investigated possible functional effects associated with sex-specific MCP-associated genes as revealed by gene expression data in multiple relevant tissues, including dorsal root ganglion (DRG) samples, in both human and mouse.

Results

GWAS of MCP in males and females separately

To detect sex-differential genetic influences on multisite chronic pain, GWASs were run separately for males and females in UK Biobank. In men, a total of 123 SNPs at five independent loci were associated with MCP at a genome-wide significance threshold of p < 5 x 10−8 (Table 1, Fig 1). In women, a total of 286 genome-wide significant SNPs at ten independent loci were discovered (Table 1, Fig 1). All 15 of these loci were differentially associated with sex—none of the genome-wide significant SNPs at these loci had p < 5 x 10−8 in the GWAS conducted in the opposite sex. However, a total of 257 SNPs were found to have suggestive levels of association with MCP (p < 5 x 10−5) in both men and women. Two SNPs had suggestive evidence in men and were genome-wide significant in women, and eight SNPs had suggestive evidence in women and were genome-wide significant in men. In addition, the genome-wide significant loci on chromosome 6 in each sex were separated by less than 1 Mbp and may potentially exert their influence via differential effects on the same gene.

Fig. 1. Manhattan plots for the stratified GWAS analyses.
Manhattan plots for the stratified GWAS analyses.
Upper panel: female, lower panel: male. Red dotted line indicates genome-wide significant p-value threshold.
Tab. 1. Genome-wide significant independent lead SNPs associated with MCP (p < 5 x 10−8) in male and female sex-stratified GWAS.
Genome-wide significant independent lead SNPs associated with MCP (p < 5 x 10<sup>−8</sup>) in male and female sex-stratified GWAS.
CHR, chromosome; BP, chromosome coordinate; A1/A2, alleles 1 and 2, where A1 is the effect allele; BETA/SE, coefficient and its standard error for the effect allele; Gene = ANNOVAR annotation gene symbol for variant (N/A, no gene annotation in ANNOVAR).

LDSR analysis demonstrated that inflation of test statistics in each GWAS was due to polygenicity (Table 2; LDSR intercept). SNP heritability was moderate, estimated as 0.125 and 0.106 in females and males, respectively (Table 2). The genetic correlation between male and female MCP was high (rg = 0.92, SE 0.03; p = 3.32 x 10−213), but significantly less than 100% (based on confidence intervals calculated as +/- 2 x SE).

Tab. 2. Trait genetic attributes from the male and female MCP GWASs.
Trait genetic attributes from the male and female MCP GWASs.
SNP-heritability, BOLT-LMM pseudo-heritability estimate; λGC1000, λGC value adjusted for sample size; LDSR_intercept (SE), LD-score regression intercept value and its standard error.

Meta-Analysis of Sex-Stratified MCP GWAS Outputs

87 independent SNPs were found to be associated with MCP at genome-wide significance in total, 11 of which were novel (not found in the unstratified or in each sex-stratified GWAS analysis). Each of the 87 independent significant SNPs showed consistent direction of effect between males and females (Table 3), but seven showed significant heterogeneity in effect size (I2 p < 0.05). In total, 49 lead SNPs across 46 genomic risk loci were found to be associated with MCP in meta-analysis of sex-stratified outputs.

Tab. 3. Independent, genome-wide significant (GWS) SNPs from meta-analysis of sex-stratified MCP GWASs.
Independent, genome-wide significant (GWS) SNPs from meta-analysis of sex-stratified MCP GWASs.
Genomic Locus = genomic risk locus context of SNP with multiple independent GWS SNPs present at some of the 46 loci, rsID = SNP rsID identifier, Position = genomic position (chromosome: base-pair start position), Meta p = p value for association from the GWAS meta-analysis, Direction = direction of effect in female MCP and male MCP GWASs, respectively (+ if association beta value for effect allele > 0,—if < 0), I2 = effect size heterogeneity estimate, I2 p = p value for heterogeneity estimate, Symbol = ANNOVAR gene annotation for variant (N/A, no gene annotation in ANNOVAR). Significant (unadjusted) heterogeneity I2 p values (I2 p < 0.05) are marked with *.

Genetic correlations between sex-stratified MCP and other disorders and traits

A range of complex trait phenotypes were selected for LDSR analysis with male and female MCP based on previous evidence for phenotypic correlation [17], with the addition of newly available trait data such as GWAS outputs on suicide and self-harm [41] and mood instability [42]. Suicidality and self-harm are important comorbidities of chronic pain, an issue compounded by the common co-occurrence of mental health traits, such as major depressive disorder (MDD), with chronic pain, and by the fact that use of certain medication in chronic pain contributes to increased risk for self-harm, suicidal ideation and suicide attempt [4348].

MCP in both men and women was found to be significantly genetically correlated with a range of traits and disorders, including psychiatric and mood phenotypes such as anhedonia, mood instability, depressive symptoms, MDD, anxiety, suicidality and subjective wellbeing (Fig 2, S1 Table). In addition, as expected, genetic correlations between both male and female MCP and unstratified MCP from our previous analysis [17]were essentially perfect (rg = 1.00, p < 1 x 10−120 for both sexes). PTSD, schizophrenia, autism spectrum disorder, anorexia nervosa, PGC cross-disorder phenotype and primary biliary cirrhosis were found to be significantly correlated with MCP in one sex and not the other. Several phenotypes were found not to be genetically correlated with MCP in either sex (pfdr > 0.05), including inflammatory bowel diseases, Parkinson’s disease, bipolar disorder, rheumatoid arthritis, and low relative amplitude (a circadian rhythmicity-related phenotype).

Fig. 2. Genetic correlations (rg) between male (M) and female (F) MCP and multiple traits (rg values with 95% CI error bars).
Genetic correlations (r<sub>g</sub>) between male (M) and female (F) MCP and multiple traits (r<sub>g</sub> values with 95% CI error bars).
Phenotypes differentially correlated with MCP between the sexes are plotted in orange.

Gene-level analysis

Genes enriched for variants associated with MCP were identified using a gene-level association analysis (gene-based test) [49] approach implemented by MAGMA as part of the FUMA suite that tests 19,012 separate genes. The results are summarised in Fig 3 and S2 Table. In females and males, 31 and 37 genes, respectively, were found to be significantly (Bonferroni-adjusted significance criterion: p < 2.63 x 10−6) associated with MCP. The only gene found to be significantly associated with both male and female MCP was DCC. 24 out of the 31 genes significantly associated with MCP in females in the sex-stratified analyses, and 31 out of the 37 genes significantly associated in males, were also significantly associated in our previous non-stratified analysis [17]. Six genes were significantly associated with MCP only in females (NCAN, SPATS2L, TBC1D9, CAMK1D, SOX11, GON4L), while 4 genes were associated only in males (CENPW, MTCH2, NICN1, DNAJA4). Twenty-four genes were identified as significantly associated with MCP only in the meta-analysis of the sex-stratified GWAS outputs (Fig 3 and S3 Table). Nineteen genes found in previous sex-combined GWAS analyses [17] were not found to be associated with MCP in either sex-stratified GWAS, or in GWAS meta-analysis (Fig 3 and S2 Table).

None of the six and four genes associated with MCP in females and males only respectively had a mouse orthologue listed in the Pain Genes Database [50], an interactive web browser listing mouse genes associated with pain-related phenotypes when knocked out. None of the six and four genes associated with MCP in females and males only respectively had a mouse ortholog listed in the Pain Genes Database [50], an interactive web browser of pain-related transgenic knockout studies. Of the twenty-four novel genes found to be associated with MCP in GWAS meta-analyses and associated downstream analyses, one had an ortholog listed in the Pain Genes Database (Gnaq)[51].

Fig. 3. Venn diagram of number of genes found significantly associated with MCP in MAGMA gene-based test analyses.
Venn diagram of number of genes found significantly associated with MCP in MAGMA gene-based test analyses.
Sectors: Original = non-sex-stratified MCP GWAS, Meta = meta-analysed female and male GWAS output, Female = female MCP GWAS, Male = male MCP GWAS.

Gene expression analysis

Tissue-enrichment of MCP-associated gene expression was analysed using FUMA, which implements a MAGMA gene-property analysis [52] using GTEx [53] gene expression datasets to determine association between trait-associated genes and expression in a range of bodily tissues. MCP-associated genes in females were found to be enriched for expression in the brain, particularly the cerebellum, and frontal cortex (Fig 4). There was no significant tissue enrichment for expression of MCP genes found in males (Fig 4).

Fig. 4. MCP-associated gene expression tissue enrichment analysis.
MCP-associated gene expression tissue enrichment analysis.
Significant (-log10 p > 3) results in female analyses and their corresponding tissue results in male analyses are shown for brevity. No significant enrichment for expression in any tissue was found for MCP-associated genes in males. Dashed line = significance threshold. ACC_BA24 = anterior cingulate cortex BA24, Caudate_BG = caudate basal ganglia, NAc_BG = Basal ganglia—Nucleus accumbens region, Putamen_BG = Basal ganglia—putamen region. Full results for all 53 tissues can be found in S1 and S2 Figs.

For genes identified (using MAGMA, as reported above) as significantly associated with MCP in either the male-only or female-only GWAS, we carried out further analyses of gene expression at tissue and cell-type level by querying existing transcriptomic datasets, focusing on neural tissues and tissues specific to each sex. For most tissues GTEx data were used, but because GTEx does not contain data for dorsal root ganglion (DRG) neurons, which are key for the generation of the nociceptive signals that initiate pain in chronic pain patients, we assessed expression in this tissue using other comparable datasets [54], which additionally contained measures of gene expression enrichment in human neural tissues (neural proportion score) and in the DRG (DRG enrichment score) (details in Methods). We also used single cell RNA-seq datasets to estimate whether genes of interest are likely to be expressed in neurons in the peripheral or central nervous systems (CNS). Full results of this analysis are given in S4 and S5 Tables, while those for a subset of genes (those enriched for neural tissue expression) are shown in Fig 5. Most of the 37 genes identified from the male-only GWAS were observed to be expressed in the nervous system. Two of the 37 genes, IP6K3 and FAM129A, had low neural proportion scores, suggesting that they are non-neuronal and non-glial. All 37 genes, however, showed DRG expression, although expression level was very low for DCC and IP6K3, which were more highly expressed in the CNS. One gene, AMIGO3, was found to be enriched solely in the DRG; its orthologue has also been found to be primarily expressed only in mouse DRG neurons in the www.mousebrain.org dataset (Fig 5). Of the 30 female-specific MCP-associated genes, all showed high neural expression except CPS1, which was not expressed in neural tissue at all. Again, we noted that all 30 genes except GABRB2 (whose mouse orthologue is primarily expressed in CNS neurons, Fig 5) did show expression in DRG (S5 Table), though none of them showed enriched expression in the DRG compared to the CNS.

Fig. 5. Mouse cell-type-specific expression patterns for sex-differentiated MCP-associated genes.
Mouse cell-type-specific expression patterns for sex-differentiated MCP-associated genes.
The mousebrain.org database was used to identify cell types of expression in the mouse nervous system (brain neurons and glia have been collapsed to two lines at the base of the plot for clarity) for orthologues of genes identified in this study. Trinarization score (posterior probability of detection in a particular cell type) was used, rather than expression levels. The analysis was carried out for 25 genes selected from the male-specific (n = 14) and female-specific (n = 11; note that Dcc appears in both lists) MCP-associated sets as being enriched for neural tissue expression (neural proportion score > 0.5, S4 and S5 Tables). While most of these show pan-neuronal expression, a few are expressed in mouse glial subpopulations (Cdhr4), primarily in CNS neurons (Gabrb2), primarily in DRG neurons (Amigo3), or in limited subsets of neurons (Dcc, Slc4a10, Camkv).

Several genes in either the male-only and female-only lists (e.g. CPS1, SEMA3F, MST1, MST1R, SDK1, ECM1; S4 and S5 Tables) that were expressed in neural tissues but with low (< 0.5) neural proportion scores were found to be involved in immune function. Among the genes with high (> 0.5) neural proportion scores, mouse orthologues of many are pan-neuronal in expression based on the mousebrain.org dataset (Fig 5), with a few genes having ubiquitous but neural-enriched expression (Mrps21, Ip6k1), solely glial expression (Cdhr4), or being expressed in limited neuronal subpopulations including CNS, sensory and enteric neuronal subtypes (Dcc, Camkv, Slc4a10 respectively). Several of the neuronally expressed genes are known to be involved in axon pathfinding and neurite outgrowth. Many of these genes have elevated expression levels in sex-specific tissues like testis and ovary, and / or are somewhat differentially expressed (10% or more difference in median TPMs across sexes) between male and female cohorts in brain sub-regions in the GTEx database (S4 and S5 Tables), suggesting they may be androgen- or estrogen-regulated.

Pathway analyses

No significant results were found for male MCP-associated genes using FUMA GENE2FUNC [55]. There was no overrepresentation of male MCP-associated genes within any of the Molecular Signatures Database (MSigDB) gene sets (h, c1-c7) [5658]. Two of the male MCP-associated genes were associated with disease-related entries in OMIM (Online Mendelian Inheritance in Man) [59], neither of which mention a pain-related phenotype, and none of the genes were listed as previously identified drug targets in DrugBank [60] (S6 Table). Amongst the female MCP-associated genes, overrepresentation was only found for two of the positional gene sets listed in MsigDB (S3 Fig). Nineteen female MCP genes were associated with an OMIM entry, and five were listed as drug targets in DrugBank (S7 Table).

Polygenic risk score analysis

Polygenic risk scores were used to assess whether MCP and chronic widespread pain (CWP), a related but distinct chronic pain phenotype, are likely to be related biologically. For both men and women, the sex-specific PRS-MCP was significantly associated with CWP (p < 2 x 10−17; O.R. = 1.0034 (female) and 1.0026 (male); S8 and S9 Tables), indicating that increased genetic risk for MCP is significantly associated with having chronic widespread pain.

Discussion

Comparing male and female multisite chronic pain

Prevalence, coping strategies, and, potentially, mechanisms of development and maintenance of chronic pain vary between the sexes. To explore underlying genetic differences that may contribute to these sex differences in chronic pain, we carried out a large-scale sex-stratified GWAS of a quantitative chronic pain phenotype, MCP. We found both male and female MCP to be moderately heritable. Although the estimated female SNP heritability was higher than that in males (12.5% versus 10.6% respectively), this difference was not significant.

MCP-associated SNP loci

Twice as many genomic risk loci were identified in GWAS analyses in females as in males (10 versus 5 respectively), with no risk loci shared between the sexes. This may be due to lower sample size in the male MCP GWAS, both in comparison to previous non-stratified analyses and to sample size in the female MCP GWAS analysis (roughly 30,000 more participants are included in the female MCP GWAS than in the male GWAS). Loci were found across the genome (Fig 1 and Table 3) in both males and females, with genomic location varying by sex. We note also, however, that significant loci were found in both men and women on chromosome 6, only 0.89 Mbp apart (Table 3).

Additional trait-associated SNPs were discovered when the sex-stratified GWAS outputs were meta-analysed, likely due to increased power. Twenty-four loci that had not reached genome-wide significance in our previous sex-combined analysis were found to be associated with MCP after meta-analysis. The fact that these loci were not identified in the non-stratified GWAS could be due to effect heterogeneity, in terms either of direction or magnitude, between the sexes reducing the overall signal at these loci.

Previous studies have highlighted sex-specific or sex-differentiated loci in a range of disorders and traits, such as ASD, anthropometric traits and asthma [31,6163]. Four of the loci found to be associated with MCP in our meta-analysis of sex-stratified GWAS outputs, including 7 SNPs in total, showed signs of heterogeneity of effect size between the sexes (I2 p < 0.05), though the evidence was not significant after FDR-adjustment. Nevertheless, if these results are replicated in future studies, it may be that these loci are found to contribute to sex differences in chronic pain.

Genes of interest

Genes associated with MCP in males

The sex-stratified gene-level analysis discovered a number of genes with significant evidence for a genetic contribution to MCP. Significant genes in males included CENPW, MTCH2, NICN1, AMIGO3, DNAJA4, CTBP2 and NOP14, with the latter two also being significant in the meta-analysis gene-level testing (S2 and S3 Tables). CENPW encodes centromere protein W, involved in kinetochore assembly and function, and associated with diseases such as type 1 diabetes [6466]. MTCH2 (mitochondrial carrier 2) encodes a member of the SLC25 family, a family of transporters localised to the inner membrane of mitochondria and involved in a wide range of cell metabolism functions [67]. SNPs in this locus have previously been associated with obesity [6870], and this gene may be involved in regulation of development of adipocytes. NICN1 (nicolin 1) encodes a nuclear protein of unknown function expressed in a variety of tissues [71]. AMIGO3 (adhesion molecule with Ig-like domain 3), the only male-specific MCP-associated gene whose expression was enriched in DRG, is a member of a small family of cell-surface immunoglobulin domain- and leucine-rich repeat-containing adhesion molecules. Its function is not well understood, but it is expressed in a range of DRG neuronal subtypes and may play a specialized role in nociception or other sensory modalities. Interestingly, AMIGO3 is located almost next to and within 40 kbp of IP6K1, which was also found to be associated in the male gene-level analysis (S2 Table). It may be that there is coordinate regulation of these two genes in tissues relevant to MCP, or that a number of separate functional variants are distributed across this genomic locus but in fact only influence expression of one of these two genes. DNAJA4 (DnaJ Heat Shock Protein Family (Hsp40) Member A4) encodes a heat shock protein [72] previously shown to be involved in melanoma metastasis and angiogenesis regulation, but is generally poorly characterised [73].

Genes associated with MCP in females

Genes found to be associated with MCP in females included NCAN, SPATS2L, TBC1D9, CAMK1D, SOX11, GON4L, and DAGLB, the last of which was also significantly associated in the meta-analysis. NCAN (neurocan) encodes a chondroitin sulfate proteoglycan [74] potentially involved in the modulation of cell adhesion and migration, and previously linked to bipolar disorder in GWAS and mouse model studies [75,76]. SPATS2L (spermatogenesis associated serine rich 2 like) encodes a protein that may be involved in ribosome biogenesis and translational control as a response to oxidative cellular stress [77]. TBC1D9 (TBC1 domain family member 9) encodes a potential GTPase and was found to be overexpressed in mantle cell lymphoma [78]. TBC1D9 was also recently found to be involved in a Ca2+-dependent cellular response to infection [79]. CAMK1D (calcium/calmodulin dependent protein kinase ID) encodes a member of the calcium/calmodulin-dependent protein kinase 1 family involved in granulocyte regulation, activating CREB-dependent gene transcription, the activation and differentiation of neutrophils, promotion of basal dendritic growth of hippocampal neurons, and apoptosis in erythroleukemia cells [80]. SOX11 (SRY-box transcription factor 11) encodes a member of the SOX (SRY-related HMG-box) family of transcription factors, with potential roles both in nervous system development and in neurogenesis during adulthood [8185]. De novo mutations in this gene have also been associated with Coffin-Siris syndrome [86]. GON4L (Gon-4 like) encodes a protein involved in transcriptional repression [87,88]. DAGLB (diacylglycerol lipase beta) encodes an enzyme that participates in the endocannabinoid synthesis pathway and is required for axonal growth during development and for retrograde synaptic signalling in mature synapses [89].

Overall, it was notable that genes found to be significantly associated with male and female MCP were largely different. Only one gene, DCC (DCC netrin 1 receptor; a.k.a. deleted in colorectal carcinoma), was associated with both male and female MCP. DCC encodes a receptor for the guidance cue netrin 1, and is important for development of the nervous system, particularly the dopaminergic system [90]. Mutations in the DCC gene have been found in those with congenital mirror movement disorder (MRMV-1; [91] and it has also been previously associated with a range of complex brain-related traits, including suicidality, mood instability, intelligence and putamen volume [41,42,9294]. DCC has also been highlighted as a risk gene for major depression, and may be involved in the pathology of depression through effects on axon guidance in the developing and adult CNS [95,96].

Gene expression differences in male and female MCP

Gene expression analyses in GTEx (which does not have DRG samples) carried out using FUMA indicated that expression of the female-specific MCP-associated genes was enriched primarily in brain tissue, and this pattern was also seen when meta-analysed sex-stratified GWAS outputs were analysed similarly (not shown). Almost all of these genes were also expressed in the human DRG—it is interesting to speculate as to whether the role of these genes in initiation or maintenance of chronic pain phenotypes is mediated through roles in the brain or via effects on cells located within the DRG. No significant enrichment for specific tissues was seen in analyses of male MCP-associated genes using GTEx and the human DRG expression profiles. The lack of tissue-enrichment findings for the male-specific genes may, as with the lower number of MCP-associated loci, be due to reduced power resulting from lower sample size in the male GWAS. However, these patterns of tissue-level gene expression in GTEx may also indicate differing gene expression between the sexes, with more ubiquitous expression across all tissues for genes associated with MCP in males, while genes conferring risk in females may tend to have more tissue-specific expression patterns. These expression patterns may also be associated with the fact that the GTEx resource is enriched for male tissue samples (V8 release; 67.1% male)–sex-differential enrichment of certain genes may be conflated with tissue-differential gene enrichment.

It is still notable that almost all the male or female MCP-associated genes were found to be expressed in human CNS tissues and in DRG. A subset of these are enriched in human neural tissues and additionally are expressed in mouse neuronal subpopulations when examining single cell sequencing databases (Fig 4). All of these lines of evidence, together, suggest putative central and peripheral neuronal roles for some of these genes, many of which have not been historically well studied in the field of chronic pain.

Genetic correlations

For both males and females, MCP was genetically correlated with non-stratified MCP [17], at rg = 1. In contrast, the genetic correlation between male MCP and female MCP was found to be high but significantly less than 1 (rg = 0.92), with around 8% of common SNP-tagged trait variation therefore not shared between the two traits. However, it could be argued that although this figure suggests that a small subset of genetic variation linked to MCP is unique to each sex, this difference is not large enough to consider MCP in the two sexes as biologically distinct to a considerable extent—traits correlated at lower rg are routinely used as proxies for one another in GWAS settings e.g. educational attainment as proxy for intelligence (rg ~ 70%) [92], or current age as a proxy for life span (rg ~40–70%) [97].

Genetic correlations with a range of psychiatric disorders, psychological traits and somatic traits and disorders were explored and notable sex-related commonalities and differences were observed. Significant genetic correlations of similar magnitude for both sexes were found between MCP and a range of psychiatric, autoimmune and anthropometric traits. Some phenotypes were significantly genetically correlated only with female MCP (schizophrenia, rg = 0.13; PGC cross-disorder phenotype, rg = 0.14; PTSD, rg = 0.44; anorexia nervosa, rg = -0.08), while others were significantly genetically correlated only with male MCP (autism spectrum disorder, rg = -0.16; primary biliary cholangitis, rg = 0.14). In most cases, where one sex was significant and the other was not, the rg values for both sexes were not very different, suggesting that the underlying biology is only subtly, quantitatively different. The sex differences that were observed may reflect sex-differential patterns of pleiotropy, with some genetic factors contributing differentially to multiple phenotypes in males and females. It would be of interest to investigate the relative effect size of loci contributing to MCP and to schizophrenia and PTSD in the two sexes and the biological mechanisms underlying this difference.

Observed genetic correlation differences may also, however, be due to differences in sample size between male and female MCP GWASs (with lower sample size for men compared to women). Differences may also be a result of either men or women being over-represented in the comparison GWAS. For example, in the GWAS meta-analysis of autism spectrum disorder [98] contributing cohorts had M:F ratios from 1.2:1 to as high as 8.6:1, PBC GWASs contributing to the discovery set in the PBC GWAS meta-analysis [99] contained > 90% female cases, and the anorexia nervosa GWAS [100] contained only female cases. Analysis of the cohorts used in the schizophrenia and PTSD GWAS meta-analyses [101106], suggests that the sex difference in genetic correlation with MCP for these phenotypes is not primarily driven by an overrepresentation of one sex.

Comparing the relationship between chronic widespread pain and MCP in males and females

Each sex-specific PRS was significantly associated with CWP in the corresponding sex, but the magnitude of association was much lower in comparison to the sex-combined PRS analysis reported previously [17], which may reflect the smaller sample sizes in the sex-stratified analyses. As previously found in sex combined analyses, results indicated a moderate degree of shared genetic basis for MCP and chronic widespread pain in both sexes, with degree of sharing potentially slightly stronger in females than in males. It is possible that this difference is driven at least in part by the overrepresentation of females in CWP cases (M:F ratio 1:1.74).

Clinical perspective on findings and potential impact on treatment

Overall, our findings suggest that MCP shows genetic (and therefore biological) differences between men and women. If women with MCP are more likely to be at risk of PTSD and schizophrenia (and vice versa) specific screening for these potential comorbidities could be appropriate, with a view to instigating additional appropriate management or referrals. Except for PTSD, however, these sex-differential genetic correlations were relatively small in comparison with those that affected both sexes, particularly a range of mood phenotypes. Nevertheless, the shared biology evident in the significant genetic correlations should stimulate investigation of whether screening for these potential clinically important co-morbidities can be used to improve management of chronic pain patients.

Enhanced attention to the sex differences in manifestation and underlying biology of chronic pain is also merited, and may lead to improvements in clinical assessment, awareness of risks and choice of medical treatment. For example, if MCP in women is more strongly associated with immune function (based on evidence from MCP-associated genes), efficacy and side-effect profiles for drugs targeting immune system function may be different in women and men. Inappropriate treatments, such as chronic opioid prescribing, might also have sex-differential consequences—opioids are known to adversely affect immune function [107,108]. These sex differences in MCP biology may also inform the search for new or re-purposed drugs that can be prescribed in a sex specific manner. It is already known that specific proteins play sex-specific roles in pain processing and that certain drugs have been found to have sex-specific analgesic effects (reviewed by [109]).

Limitations

These analyses were carried out using UK Biobank, which was used in our previous sex-combined MCP GWAS [17]. In comparison to our previous analysis, sample size in each individual sex-stratified GWAS was lower, which leads to somewhat reduced power. However, sample sizes are still larger than many sex-combined GWAS analyses of chronic pain phenotypes, and meta-analysis of the sex-stratified GWAS outputs resulted in an increase in power to find MCP-associated SNPs overall. True replication is difficult due to heterogeneity in chronic pain phenotyping and available sample sizes of potential independent cohorts, but in an independent subset of the UK Biobank, a PRS constructed from each sex-stratified GWAS output was found to be significantly associated with CWP, a related but distinct chronic pain phenotype of interest.

Although this work was focused on sex differences in the genetics of MCP we examined only autosomal variation. An important extension of this work would be to assess genetic associations with X chromosome loci, which is likely to provide an additional heritability contribution and give a fuller picture of sex differences in MCP at the genetic level. Inclusion of the X chromosome in GWAS analyses is associated with specific methodological and statistical issues including lower quality genotyping array coverage of the X chromosome compared to autosomes, differences in how imputation needs to be implemented, differences in X chromosome dosage between the sexes leading to differences in population genetics/demographic history of the X chromosome relative to the autosomes, and changes to quality control required (and differing QC protocols between sexes) [110,111]. We aim to address this in future work as we adapt our BOLT-LMM pipeline and downstream analysis pathway.

These GWAS analyses were carried out on a white British subset of UK Biobank, and therefore may not generalise to admixed or non-white populations. GTEx donors are also primarily white (v8 release: 84.6% white, 12.9% African American, 0.2% American Indian, 1.3% Asian, 1.1% Unknown, [https://www.gtexportal.org/home/tissueSummaryPage]), as are donors who provided DRG tissue (white females [54]), again potentially limiting generalisability of our findings to non-white populations.

Conclusions

Sex differences in chronic pain likely have, at least in part, a genetic basis and the study of complex traits such as chronic pain is likely to benefit from “sex-aware” analytical approaches. This study comprises one of the largest sex-stratified genetic analyses of a chronic pain phenotype, and highlights sex-differential MCP-associated loci, genes, genetic correlations, and patterns of tissue expression. We also examined transcriptome abundance of key sex-differential MCP-associated genes in a range of neural and non-neural tissues, including DRG, an important nervous system component in chronic pain which is not part of the GTEx resource and so may be understudied in GWASs with follow-up conducted solely using FUMA.

Sex-stratified GWASs can provide an increase in power if heterogeneity in effects of trait-associated variants is seen between the sexes. Here, 24 novel genes and 11 novel independent lead SNPs were associated with MCP, in addition to the findings from previous non-sex-stratified work, further contributing to understanding of genetic variation predisposing to chronic pain.

Genetic correlation results indicated possible sex-differential pleiotropy, including differing genetic correlations between certain psychiatric disorders and traits and chronic pain in women compared to men. However, it is of note that genetic correlations between the sexes are largely similar, particularly with respect to psychiatric disorders and traits including MDD, anhedonia and depressive symptoms. This is the first study to use novel GWAS outputs from studies of suicidality and of RDoC mental health traits such as mood instability and anhedonia in genetic correlation analyses with chronic pain and it provides an important insight into shared genetic factors between these comorbidities of chronic pain and MCP.

The patterns of gene expression enrichment associated with the identified predisposing genes supports and enhances our previous conclusion that MCP derives more strongly from brain and/or CNS-based mechanisms than from other organ systems and functions. We did observe, however, that almost all the sex-specific MCP-associated genes identified are expressed in the DRG as well as in the brain, with one sex differentiated MCP gene, AMIGO3, being DRG-specific amongst the tissues assessed, raising the possibility that cells in the DRG that play a specialised role in nociception or other sensory function might be involved in chronic pain mechanisms. We also observed that several MCP-associated loci in both sexes have been linked to immune function.

Overall, our findings indicate the existence of potential sex differences in chronic pain at multiple levels, from SNP-level to transcript abundance and the results support theories of strong nervous system and immune involvement in chronic pain in both sexes. These findings may inform development of novel treatment approaches in future, as well as adding to our understanding of the physiology of chronic pain.

Methods

Chronic pain phenotyping

UK Biobank participants were asked via a touchscreen questionnaire about “pain type(s) experienced in the last month” (field ID 6159), with possible answers: ‘None of the above’; ‘Prefer not to answer’; pain at seven different body sites (head, face, neck/shoulder, back, stomach/abdomen, hip, knee); or ‘all over the body’. The seven individual body-site pain options were not mutually exclusive, but those who chose ‘all over the body’ could not also select from the seven individual body sites. Where patients reported recent pain at one or more body sites, or all over the body, they were additionally asked (category ID 100048) whether this pain had lasted for 3 months or longer.

Chronic Widespread Pain (CWP) was defined as reported [112], and the ‘case’ group included only those participants who answered that they had pain ‘all over the body’ that was longer than 3 months in duration in the touchscreen questionnaire. These individuals were excluded from analyses of other chronic pain phenotypes as there is some evidence that this phenotype can be substantially different from more localised chronic pain [113]. Multisite Chronic Pain (MCP) was a quasi-quantitative variable defined as previously reported [17]; briefly, this variable captures the number of body sites at which chronic pain (at least 3 months duration) was recorded (excluding those with CWP): phenotypic values therefore ranged from 0 to 7. 10,000 randomly selected individuals reporting no chronic pain were excluded from the GWAS to use as controls in subsequent polygenic risk score (PRS) analyses.

Genetic quality control

For the GWAS analyses, SNPs with an imputation quality score of less than 0.3, minor allele frequency (MAF) < 0.01 and/or Hardy-Weinberg equilibrium (HWE) test p < 10−6 were excluded. Participants whose self-reported sex did not match their genetically-determined sex, those who had putative sex-chromosome aneuploidy, those considered outliers in UK Biobank QC in terms of missingness or heterozygosity [114], and those who were not of self-reported white British ancestry were excluded from analyses. A list of “poor quality” samples (due to missingness, putative genetic and reported sex mismatch and/or unexpectedly high heterozygosity) was derived by Bycroft et al [114] and is available to all researchers using UK Biobank, and was used here as part of genetic quality control. Briefly, putative sex chromosome aneuploidy was defined by visual inspection of scatterplots of mean log2 ratio (L2R) on X and Y chromosomes, and 652 UKB participants meet these criteria for putative sex-chromosome aneuploidy (Supplemental Information S 3.6 [114]). Samples with a heterozygosity value, adjusted for both ancestry and genetic principal components (GPCs), above the mean heterozygosity value (0.1903) and missing rate greater than 0.05 as computed using PLINK ‘—miss’ command were also flagged as potentially poor quality ([114]; 968 such samples are listed in this paper’s Supplemental Information S 3.5.3).

A summary of participant MCP phenotypic information for those included in each GWAS is shown in Table 4.

Tab. 4. Number of participants per MCP phenotype level group included in each GWAS (male or female sex-stratified analysis).
Number of participants per MCP phenotype level group included in each GWAS (male or female sex-stratified analysis).

BOLT-LMM GWAS and gene-level analysis in FUMA

Sex-stratified GWASs of MCP, modelled as a quantitative trait, were carried out using BOLT-LMM [115], adjusting for age and chip (genotyping array), under the infinitesimal model of genetic risk, as previously described for our unstratified MCP GWAS [17]. BOLT-LMM uses a genetic relatedness matrix in the model to adjust for population stratification and familial relationships. We therefore did not include genetic principal components (GPCs) in the model. The SNP-level GWAS summary statistics were then analysed using FUMA [55] to obtain genome-wide plots and carry out MAGMA [49] gene-set and gene-based test analyses and gene expression analysis using GTEx [53] for male-enriched and female-enriched MCP genes. Significant independent lead SNPs were determined according to FUMA. Briefly, FUMA defines lead SNPs as the subset of independent significant SNPs (SNPs associated with the trait at p < 5 x 10−8 and having LD r2 < 0.6 with any other significant SNP) that are not in LD (r2 > 0.1) with any other lead SNP [55]. In addition, when these LD blocks of independent significant SNPs are in close proximity (< 250kbp apart), separate loci based on LD thresholds are merged into a single genomic locus, and thus each genomic risk locus can contain multiple lead and independent significant SNPs [55].

Meta-analysis of male and female GWAS Summary statistics

Meta-analysis of the two sex-specific GWAS summary statistics datasets was carried out using METAL [116], deploying a fixed-effects model and weighted by standard error (‘SCHEME STDERR’) with default options aside from selecting ‘heterogeneity’ in order to analyse heterogeneity (‘ANALYZE HETEROGENEITY’). A meta-analysis p-value of < 5 x 10−8 was selected as the significance threshold for association. Gene-level analysis (MAGMA) was also carried out using meta-analysis output.

Transcriptome analysis of sex-specific association gene lists

We further analysed the tissue and cell type of expression of male-specific and female-specific genes identified from the MAGMA gene-based analysis of the sex-stratified GWAS results. Specifically, we characterized gene expression in mammalian nervous system tissues and cell types as a potential starting point for identifying the functions of these genes with respect to pain. We also characterised expression for the one gene found to be associated with both male and female MCP in the gene-based analyses (DCC).

RNA-seq-derived gene expression values have been previously reported [54] as relative abundances in TPMs (standardised Transcripts per Million mapped reads) for 12 adult human tissues (6 neural and 6 non-neural). Additionally, the study notes 3 metrics on a scale of 0 to 1 for each gene based on expression in these 12 tissues: normalized Shannon’s entropy as a measure of tissue specificity (0 for highly tissue-specific and 1 for tissue-agnostic gene expression in the quantified tissues), neural proportion score as a measure of enriched expression (possibly neuronal and/or glial) in the nervous system (0 for genes not expressed in the nervous system and 1 for genes expressed solely in neural tissues with respect to these 12 tissues), and DRG enrichment score for identifying specificity of gene expression in the DRG with respect to the other 11 profiled tissues (0 for no expression in the DRG or for tissue-agnostic gene expression, and 1 for DRG-specific gene expression in the context of the set of profiled tissues). The mathematical formulations for these scores are provided in detail in Ray et al. 2018 [54]. The corresponding tables are presented in S4 Table (male) and S5 Table (female).

Genes with neural proportion scores > 0.5 (with overall more neural than non-neural tissue expression), were further characterized by the putative cell type(s) of gene expression in the mammalian nervous system (Fig 4). While human single cell resolution RNA-seq datasets are not publicly available for the peripheral nervous system, a comprehensive database of gene expression in the mouse nervous system exists: the www.mousebrain.org repository [117]. While it is true that nervous system expression patterns may be different between human and mouse between nervous system subpopulations due to regulatory evolution or due to differences in nervous system cell types [118], it is unlikely that expression would change categories across the 4 broad nervous system cell type categories: neurons, glia, immune and vascular cells) between human and mouse, given overall conservation of tissue gene expression profiles across humans and mice [54]. Trinarization scores, defined as the posterior probability (using a Bayesian framework) of detecting reads from a particular gene in a cell type subpopulation (details are given in ref. [117]), were used to characterise gene expression by cell-type. Fig 4 visualizes trinarization scores for a range of peripheral nervous system cell types (sensory neurons and glia, enteric neurons and glia, and sympathetic neurons). While peripheral nervous system vascular and immune cells have not been profiled so far in www.mousebrain.org, we depict CNS vascular immune cells and vascular cells as surrogate cell types for their PNS counterparts. Both CNS neurons and glia play a critical role in chronic pain, but due to high diversity of cell types were not ideal for summarizing expression patterns in each subpopulation of these categories succinctly in a single figure. Instead, summary rows for CNS neurons and glia expression for relevant genes are provided at the bottom of the figure. For further details of expression in these subtypes, www.mousebrain.org can be queried for mouse gene expression profiles.

Finally, for genes with neural proportion scores > 0.5, the GTEx database was queried for sex differential gene expression in profiled CNS regions, and for high expression in sex-specific tissues, such as testis and ovary, noted in the comments section of S4 and S5 Tables.

Genetic correlations

Genetic correlations with a range of neuropsychiatric disorders and traits were assessed by LD-score regression (LDSR) for the male and female MCP GWAS outputs separately. Summary statistics datasets employed were publicly available or available via LD Hub [119,120], or were results from published and unpublished in-house GWASs. Pre-computed LD scores were used, along with HapMap3 SNPs as a reference. By default the munging part of LDSR filters out SNPs with an info score of < 0.9, a MAF of < 0.01, a GWAS p value outside the range 0–1, strand-ambiguous SNPs, non-SNP (e.g. indel) variants, and SNPs with low sample size. The default option of no constraint on the LDSR intercept was used. LDSR p-values for genetic correlation were FDR-corrected within each sex. LDSR was also used to assess trait polygenicity and to calculate a SNP-heritability estimate. Note that although some of the GWAS results used in this analysis came from studies that included UK Biobank data, genetic correlations estimated using LDSR are not subject to bias caused by sample overlap [121].

Polygenic risk score analysis

Previous analyses showed that a polygenic risk score (PRS) for MCP was associated with the phenotype of chronic widespread pain (CWP) in women but not in men, and that this PRS was associated more strongly with chronic pain phenotypes in women than in men in an independent cohort (Johnston K.J.A. et al., unpublished). To further explore the relationship between MCP and CWP, we assessed how separate male and female PRSs for MCP were associated with CWP. Separate sex-specific PRSs for MCP, based on the sex-specific GWAS results reported here, were calculated for a ‘case’ group consisting of participants who reported pain all over the body that lasted for three months or longer (a proxy phenotype for Chronic Widespread Pain; CWP; N = 6, 813)), and for a ‘control’ group consisting of 10,000 randomly selected UKB controls, both of which had been excluded from the GWAS analyses (demographic data for this subsample are given in Table 5).

Tab. 5. Summary of sample sizes of participants used in each of the two PRS analyses (male or female).
Summary of sample sizes of participants used in each of the two PRS analyses (male or female).

SNPs associated with MCP at p < 0.01 in the original sex-specific GWAS were selected and LD-pruned (at a threshold of r2 < 0.1 within a 250kbp window using the PLINK ‘—clump’ command). Sex-specific PRS-MCPs were calculated for each individual in the analysis as the sum of risk alleles at each SNP, weighted by effect size (beta value) in the GWAS [122]. PRS values were standardised and z-scores were used in the analysis. Association between standardised sex-specific PRS-MCP and CWP status was investigated separately in males and in females in the target case-control subsample using logistic regression, adjusted for chip (genotyping array), age and the first eight genetic principal components.

Supporting information

S1 Table [pdf]
Genetic correlations between MCP and other disorders and traits.

S2 Table [male]
Table accompanying Venn diagram of genes associated with multisite chronic pain in MAGMA gene-level analyses.

S3 Table [pdf]
Genes associated with MCP in MAGMA gene-level analyses of GWAS meta-analysis output.

S4 Table. Expression of genes associated with male MCP in MAGMA analyses across neural and non-neural tissues. TPM = transcripts per million, DRG = dorsal root ganglion. h_DRG_enrich = DRG enrichment score. h_entropy = normalized Shannon’s entropy. h_neural_propn = neural proportion score [pdf]

S5 Table [pdf]
Expression of genes associated with female MCP in MAGMA analyses across neural and non-neural tissues.

S6 Table [pdf]
OMIM and DrugBank derived information (male MCP).

S7 Table [pdf]
OMIM and DrugBank derived information (female MCP).

S8 Table [pdf]
Association between male MCP PRS and CWP in men.

S9 Table [pdf]
Association between female MCP PRS and CWP in women.

S1 Fig [pdf]
Tissue Expression of MCP-associated Genes (male) (53 tissues).

S2 Fig [pdf]
Tissue Expression of MCP-associated Genes (female) (53 tissues).

S3 Fig [tiff]
FUMA gene set analysis results for female MCP genes from the gene-level analysis.


Zdroje

1. Merskey H, Bogduk N. Classification of Chronic Pain. IASP Pain Terminology. 1994. 240 p.

2. IASP. International Association for the Study of Pain [Internet]. [cited 2020 Apr 14]. Available from: https://www.iasp-pain.org/

3. Nicholasa M, Vlaeyenb JWS, Riefe W, Barkee A, Azizf Q, Benolielg R, et al. The IASP classification of chronic pain for ICD-11: Chronic primary pain. Pain. 2019;160(1):53–9. doi: 10.1097/j.pain.0000000000001365 30586071

4. Gureje O, Von Korff M, Kola L, Demyttenaere K, He Y, Posada-Villa J, et al. The relation between multiple pains and mental disorders: Results from the World Mental Health Surveys. Pain. 2008;135(1–2):82–91. doi: 10.1016/j.pain.2007.05.005 17570586

5. Von Korff M, Crane P, Lane M, Miglioretti DL, Simon G, Saunders K, et al. Chronic spinal pain and physical-mental comorbidity in the United States: Results from the national comorbidity survey replication. Pain. 2005;113(3):331–9. doi: 10.1016/j.pain.2004.11.010 15661441

6. Santos-Eggimann B, Wietlisbach V, Rickenbach M, Paccaud F, Gutzwiller F. One-year prevalence of low back pain in two Swiss regions. Estimates from the population participating in the 1992–1993 MONICA project. Spine (Phila Pa 1976). 2000;25(19):2473–9. doi: 10.1097/00007632-200010010-00009 11013499

7. Palmer KT, Walsh K, Bendall H, Cooper C, Coggon D. Back pain in Britain: Comparison of two prevalence surveys at an interval of 10 years. Br Med J. 2000;320(7249):1577–8. doi: 10.1136/bmj.320.7249.1577 10845966

8. Goldberg DS, McGee SJ. Pain as a global public health priority. BMC Public Health [Internet]. 2011;11(1):770. Available from: http://www.biomedcentral.com/1471-2458/11/770 doi: 10.1186/1471-2458-11-770 21978149

9. Breivik H, Collett B, Ventafridda V, Cohen R, Gallacher D. Survey of chronic pain in Europe: Prevalence, impact on daily life, and treatment. Eur J Pain. 2006;10(4):287–333. doi: 10.1016/j.ejpain.2005.06.009 16095934

10. James SL, Abate D, Abate KH, Abay SM, Abbafati C, Abbasi N, et al. Global, regional, and national incidence, prevalence, and years lived with disability for 354 Diseases and Injuries for 195 countries and territories, 1990–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2018;1789–858. doi: 10.1016/S0140-6736(18)32279-7 30496104

11. Von Korff M, Ormel J, Keefe FJ, Dworkin SF. Grading the severity of chronic pain. Pain [Internet]. 1992;50(1092):133–49. Available from: http://www.ncbi.nlm.nih.gov/pubmed/1408309 doi: 10.1016/0304-3959(92)90154-4 1408309

12. Hocking LJ, Morris AD, Dominiczak AF, Porteous DJ, Smith BH. Heritability of chronic pain in 2195 extended families. Eur J Pain (United Kingdom). 2012;16(7):1053–63. doi: 10.1002/j.1532-2149.2011.00095.x 22337623

13. Junqueira DRG, Ferreira ML, Refshauge K, Maher CG, Hopper JL, Hancock M, et al. Heritability and lifestyle factors in chronic low back pain: Results of the Australian Twin Low Back Pain Study (The AUTBACK study). Eur J Pain (United Kingdom). 2014;18(10):1410–8. doi: 10.1002/ejp.506 24733726

14. Williams FMK, Spector TD, MacGregor AJ. Pain reporting at different body sites is explained by a single underlying genetic factor. Rheumatology. 2010;49(9):1753–5. doi: 10.1093/rheumatology/keq170 20525736

15. McIntosh AM, Hall LS, Zeng Y, Adams MJ, Gibson J, Wigmore E, et al. Genetic and Environmental Risk for Chronic Pain and the Contribution of Risk Variants for Major Depressive Disorder: A Family-Based Mixed-Model Analysis. PLoS Med. 2016;13(8):1–17. doi: 10.1371/journal.pmed.1002090 27529168

16. Suri P, Palmer MR, Tsepilov YA, Freidin MB, Boer CG, Yau MS, et al. Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain. 2018;48(7):1–23.

17. Johnston KJA, Adams MJ, Nicholl BI, Ward J, Strawbridge RJ, Ferguson A, et al. Genome-wide association study of multisite chronic pain in UK Biobank. PLoS Genet. 2019;15(6):1–22. doi: 10.1371/journal.pgen.1008164 31194737

18. Rawlik K, Canela-Xandri O, Tenesa A. Evidence for sex-specific genetic architectures across a spectrum of human complex traits. Genome Biol [Internet]. 2016;17(1):1–8. Available from: http://dx.doi.org/10.1186/s13059-016-1025-x 27473438

19. Khramtsova EA, Davis LK, Stranger BE. The role of sex in the genomics of human complex traits. Nat Rev Genet [Internet]. 2019;20(3):173–90. Available from: doi: 10.1038/s41576-018-0083-1 30581192

20. Xu X, Coats JK, Yang CF, Wang A, Ahmed OM, Alvarado M, et al. Modular genetic control of sexually dimorphic behaviors. Cell [Internet]. 2012;148(3):596–607. Available from: doi: 10.1016/j.cell.2011.12.018 22304924

21. Quinn MA, Cidlowski JA. Endogenous hepatic glucocorticoid receptor signaling coordinates sex-biased inflammatory gene expression. FASEB J. 2016;30(2):971–82. doi: 10.1096/fj.15-278309 26581598

22. McCormick H, Young PE, Hur SSJ, Booher K, Chung H, Cropley JE, et al. Isogenic mice exhibit sexually-dimorphic DNA methylation patterns across multiple tissues. BMC Genomics. 2017;18(1):1–9. doi: 10.1186/s12864-016-3406-7 28049423

23. Gilks WP, Abbott JK, Morrow EH. Sex differences in disease genetics: Evidence, evolution, and detection. Trends Genet [Internet]. 2014;30(10):453–63. Available from: doi: 10.1016/j.tig.2014.08.006 25239223

24. Ge T, Chen CY, Neale BM, Sabuncu MR, Smoller JW. Phenome-wide heritability analysis of the UK Biobank. PLoS Genet. 2017;13(4):1–21. doi: 10.1371/journal.pgen.1006711 28388634

25. Rahmioglu N, MacGregor S, Drong AW, Hedman ÅK, Harris HR, Randall JC, et al. Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci. Hum Mol Genet. 2015;24(4):1185–99. doi: 10.1093/hmg/ddu516 25296917

26. Hall E, Volkov P, Dayeh T, Esguerra JL o. S, Salö S, Eliasson L, et al. Sex differences in the genome-wide DNA methylation pattern and impact on gene expression, microRNA levels and insulin secretion in human pancreatic islets. Genome Biol. 2014;15(12):522. doi: 10.1186/s13059-014-0522-z 25517766

27. Kukurba KR, Parsana P, Balliu B, Smith KS, Zappala Z, Knowles DA, et al. Impact of the X chromosome and sex on regulatory variation. Genome Res. 2016;26(6):768–77. doi: 10.1101/gr.197897.115 27197214

28. Yao C, Joehanes R, Johnson AD, Huan T, Esko T, Ying S, et al. Sex- and age-interacting eQTLs in human complex diseases. Hum Mol Genet [Internet]. 2013 Nov 15;23(7):1947–56. Available from: doi: 10.1093/hmg/ddt582 24242183

29. Kósa JP, Balla B, Speer G, Kiss J, Borsy A, Podani J, et al. Effect of menopause on gene expression pattern in bone tissue of nonosteoporotic women. Menopause. 2009;16(2):367–77. doi: 10.1097/gme.0b013e318188b260 19512969

30. Gomez-Santos C, Hernandez-Morante JJ, Margareto J, Larrarte E, Formiguera X, Martínez CM, et al. Profile of adipose tissue gene expression in premenopausal and postmenopausal women: Site-specific differences. Menopause. 2011;18(6):675–84. doi: 10.1097/gme.0b013e31820641da 21358552

31. Mitra I, Tsang K, Ladd-Acosta C, Croen LA, Aldinger KA, Hendren RL, et al. Pleiotropic Mechanisms Indicated for Sex Differences in Autism. PLoS Genet. 2016;12(11):1–27. doi: 10.1371/journal.pgen.1006425 27846226

32. Bartley EJ, Fillingim RB. Sex differences in pain: A brief review of clinical and experimental findings. Br J Anaesth. 2013;111(1):52–8. doi: 10.1093/bja/aet127 23794645

33. Fillingim RB, King CD, Ribeiro-Dasilva MC, Rahim-Williams B, Riley JL. Sex, Gender, and Pain: A Review of Recent Clinical and Experimental Findings. J Pain [Internet]. 2009;10(5):447–85. Available from: doi: 10.1016/j.jpain.2008.12.001 19411059

34. Fillingim RB. Biopsychosocial contributions to sex differences in pain. BJOG An Int J Obstet Gynaecol. 2015;122(6):769. doi: 10.1111/1471-0528.13337 25752330

35. El-Shormilisy N, Strong J, Meredith PJ. Associations among gender, coping patterns and functioning for individuals with chronic pain: A systematic review. Pain Res Manag. 2015;20(1):48–55. doi: 10.1155/2015/490610 24927488

36. Sorge RE, Mapplebeck JCS, Rosen S, Beggs S, Taves S, Alexander JK, et al. Different immune cells mediate mechanical pain hypersensitivity in male and female mice. Nat Neurosci. 2015;18(8):1081–3. doi: 10.1038/nn.4053 26120961

37. Sorge RE, Totsch SK. Review Sex Differences in Pain. 2017;1281(February 2016):1271–81.

38. Agalave NM, Rudjito R, Farinotti AB, Khoonsari PE, Sandor K, Nomura Y, et al. Sex-dependent role of microglia in disulfide HMGB1-mediated mechanical hypersensitivity. Pain. 2020;Articles i.

39. Klein SL, Flanagan KL. Sex differences in immune responses. Nat Rev Immunol. 2016;16(10):626–38. doi: 10.1038/nri.2016.90 27546235

40. Smith BH, Campbell A, Linksted P, Fitzpatrick B, Jackson C, Kerr SM, et al. Cohort profile: Generation scotland: Scottish family health study (GS: SFHS). The study, its participants and their potential for genetic research on health and illness. Int J Epidemiol. 2013;

41. Strawbridge RJ, Ward J, Ferguson A, Graham N, Shaw RJ, Cullen B, et al. Identification of novel genome-wide associations for suicidality in UK Biobank, genetic correlation with psychiatric disorders and polygenic association with completed suicide. EBioMedicine [Internet]. 2019;41:517–25. Available from: doi: 10.1016/j.ebiom.2019.02.005 30745170

42. Ward J, Tunbridge EM, Sandor C, Lyall LM, Ferguson A, Strawbridge RJ, et al. The genomic basis of mood instability: identification of 46 loci in 363,705 UK Biobank participants, genetic correlation with psychiatric disorders, and association with gene expression and function. Mol Psychiatry [Internet]. 2019; Available from: http://dx.doi.org/10.1038/s41380-019-0439-8

43. Okifuji A, Benham B. Suicidal and Self-Harm Behaviors in Chronic. J Appl Biobehav Res. 2011;16(2):57–77.

44. Nock MK, Borges G, Bromet EJ, Cha CB, Kessler RC, Lee S. Suicide and Suicidal Behavior Matthew. Epidemiol Rev. 2008;30(1):133–54.

45. Cheatle MD. Depression, chronic pain, and suicide by overdose: On the edge. Pain Med. 2011;12(SUPPL. 2):43–8.

46. Tang NKY, Crane C. Suicidality in chronic pain: A review of the prevalence, risk factors and psychological links. Psychol Med. 2006;36(5):575–86. doi: 10.1017/S0033291705006859 16420727

47. Triñanes Y, González-Villar A, Gómez-Perretta C, Carrillo-de-la-Peña MT. Suicidality in Chronic Pain: Predictors of Suicidal Ideation in Fibromyalgia. Pain Pract. 2015;15(4):323–32. doi: 10.1111/papr.12186 24690160

48. Smith MT, Edwards RR, Robinson RC, Dworkin RH. Suicidal ideation, plans, and attempts in chronic pain patients: Factors associated with increased risk. Pain. 2004;111(1–2):201–8. doi: 10.1016/j.pain.2004.06.016 15327824

49. de Leeuw CA, Mooij JM, Heskes T, Posthuma D. MAGMA: Generalized Gene-Set Analysis of GWAS Data. PLoS Comput Biol. 2015;11(4). doi: 10.1371/journal.pcbi.1004219 25885710

50. LaCroix-Fralish ML, Ledoux JB, Mogil JS. The Pain Genes Database: An interactive web browser of pain-related transgenic knockout studies. Pain. 2007;131(1–2):3.e1–3.e4. doi: 10.1016/j.pain.2007.04.041 17574758

51. Tappe-Theodor A, Constantin CE, Tegeder I, Lechner SG, Langeslag M, Lepcynzsky P, et al. Gα q/11 signaling tonically modulates nociceptor function and contributes to activity-dependent sensitization. Pain [Internet]. 2012;153(1):184–96. Available from: doi: 10.1016/j.pain.2011.10.014 22071319

52. Watanabe K, Taskesen E, Van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun. 2017;8(1):1–10. doi: 10.1038/s41467-016-0009-6 28232747

53. Aguet F, Brown AA, Castel SE, Davis JR, He Y, Jo B, et al. Genetic effects on gene expression across human tissues. Nature. 2017;550(7675):204–13. doi: 10.1038/nature24277 29022597

54. Ray P, Torck A, Quigley L, Wangzhou A, Neiman M, Rao C, et al. Comparative transcriptome profiling of the human and mouse dorsal root ganglia. Pain. 2018;159(7):1325–45. doi: 10.1097/j.pain.0000000000001217 29561359

55. Watanabe K, Taskesen E, Van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with FUMA. Nat Commun [Internet]. 2017;8(1):1–10. Available from: doi: 10.1038/s41467-016-0009-6 28232747

56. Liberzon A, Birger C, Ghandi M, Jill P, Tamayo P, Jolla L, et al. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015;1(6):417–25. doi: 10.1016/j.cels.2015.12.004 26771021

57. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdóttir H, Tamayo P, Mesirov JP. Molecular signatures database (MSigDB) 3.0. Bioinformatics. 2011;27(12):1739–40. doi: 10.1093/bioinformatics/btr260 21546393

58. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005;102(43):15545–50. doi: 10.1073/pnas.0506580102 16199517

59. McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University (Baltimore M. Online Mendelian Inheritance in Man, OMIM [Internet]. [cited 2020 Oct 7]. Available from: omim.org

60. Wishart DS, Knox C, Guo AC, Shrivastava S, Hassanali M, Stothard P, et al. DrugBank: a comprehensive resource for in silico drug discovery and exploration. Nucleic Acids Res. 2006;34(Database issue):668–72. doi: 10.1093/nar/gkj067 16381955

61. Winkler TW, Justice AE, Graff M, Barata L, Feitosa MF, Chu S, et al. The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study. Vol. 11, PLoS Genetics. 2015. 1–42 p. doi: 10.1371/journal.pgen.1005378 26426971

62. Randall JC, Winkler TW, Kutalik Z, Berndt SI, Jackson AU, Monda KL, et al. Sex-stratified Genome-wide Association Studies Including 270,000 Individuals Show Sexual Dimorphism in Genetic Loci for Anthropometric Traits. PLoS Genet. 2013;9(6). doi: 10.1371/journal.pgen.1003500 23754948

63. Myers RA, Scott NM, Gauderman WJ, Qiu W, Mathias RA, Romieu I, et al. Genome-wide interaction studies reveal sex-specific asthma risk alleles. Hum Mol Genet. 2014;23(19):5251–9. doi: 10.1093/hmg/ddu222 24824216

64. Prendergast L, van Vuuren C, Kaczmarczyk A, Doering V, Hellwig D, Quinn N, et al. Premitotic Assembly of Human CENPs -T and -W switches centromeric Chromatin to a mitotic state. PLoS Biol. 2011;9(6):1–12. doi: 10.1371/journal.pbio.1001082 21695110

65. Lee S, Gang J, Jeon SB, Choo SH, Lee B, Kim YG, et al. Molecular cloning and functional analysis of a novel oncogene, cancer-upregulated gene 2 (CUG2). Biochem Biophys Res Commun. 2007;360(3):633–9. doi: 10.1016/j.bbrc.2007.06.102 17610844

66. Chun Y, Park B, Koh W, Lee S, Cheon Y, Kim R, et al. New centromeric component CENP-W Is an RNA-associated nuclear matrix protein that interacts with nucleophosmin/B23 protein. J Biol Chem. 2011;286(49):42758–69. doi: 10.1074/jbc.M111.228411 22002061

67. Schwarz M, Andrade-Navarro MA, Gross A. Mitochondrial carriers and pores: Key regulators of the mitochondrial apoptotic program? Apoptosis. 2007;12(5):869–76. doi: 10.1007/s10495-007-0748-2 17453157

68. Willer CJ, Speliotes EK, Loos RJF, Li S, Lindgren CM, Heid IM, et al. Six new loci associated with body mass index highlight a neuronal influence on body weight regulation. Nat Genet. 2009;41(1):25–34. doi: 10.1038/ng.287 19079261

69. Yu K, Ganesan K, Tan LK, Laban M, Wu J, Xiao DZ, et al. A precisely regulated gene expression cassette potently modulates metastasis and survival in multiple solid cancers. PLoS Genet. 2008;4(7). doi: 10.1371/journal.pgen.1000129 18636107

70. Renström F, Payne F, Nordström A, Brito EC, Rolandsson O, Hallmans G, et al. Replication and extension of genome-wide association study results for obesity in 4923 adults from northern Sweden. Hum Mol Genet. 2009;18(8):1489–96. doi: 10.1093/hmg/ddp041 19164386

71. Backofen B, Jacob R, Serth K, Gossler A, Naim HY, Leeb T. Cloning and characterization of the mammalian-specific nicolin 1 gene (NICN1) encoding a nuclear 24 kDA protein. Eur J Biochem. 2002;269(21):5240–5. doi: 10.1046/j.1432-1033.2002.03232.x 12392556

72. Alderson TRR, Kim JHH, Markley JLL. Dynamical Structures of Hsp70 and Hsp70-Hsp40 Complexes. Structure [Internet]. 2016;24(7):1014–30. Available from: doi: 10.1016/j.str.2016.05.011 27345933

73. Pencheva N, Tran H, Buss C, Huh D, Drobnjak M, Busam K, et al. Convergent multi-mRNA Targeting of ApoE Drives LRP1/LRP8-Dependent Melanoma Metastasis and Angiogenesis. Cell. 2012;151(5):1068–82. doi: 10.1016/j.cell.2012.10.028 23142051

74. Rauch U, Karthikeyan L, Maurel P, Margolis RU, Margolis RK. Cloning and primary structure of neurocan, a developmentally regulated, aggregating chondroitin sulfate proteoglycan of brain. J Biol Chem. 1992;267(27):19536–47. 1326557

75. Cichon S, Mühleisen TW, Degenhardt FA, Mattheisen M, Miró X, Strohmaier J, et al. Genome-wide association study identifies genetic variation in neurocan as a susceptibility factor for bipolar disorder. Am J Hum Genet. 2011;88(3):372–81. doi: 10.1016/j.ajhg.2011.01.017 21353194

76. Miró X, Meier S, Dreisow ML, Frank J, Strohmaier J, Breuer R, et al. Studies in humans and mice implicate neurocan in the etiology of mania. Am J Psychiatry. 2012;169(9):982–90. doi: 10.1176/appi.ajp.2012.11101585 22952076

77. Zhu CH, Kim J, Shay JW, Wright WE. SGNP: An essential stress granule/nucleolar protein potentially involved in 5.8s rRNA processing/transport. PLoS One. 2008;3(11). doi: 10.1371/journal.pone.0003716 19005571

78. Islam TC, Asplund AC, Lindvall JM, Nygren L, Liden J, Kimby E, et al. High level cannabinoid receptor 1, resistance of regulator G protein signaling 13 and differential expression of Cyclin D1 in mantle cell lymphoma. Leukemia. 2003;17(9):1880–90. doi: 10.1038/sj.leu.2403057 12970790

79. Nozawa T, Sano S, Minowa-Nozawa A, Toh H, Nakajima S, Murase K, et al. TBC1D9 regulates TBK1 activation through Ca2+ signaling in selective autophagy. Nat Commun [Internet]. 2020;11(1):1–16. Available from: doi: 10.1038/s41467-019-13993-7 31911652

80. Verploegen S, Lammers JWJ, Koenderman L, Coffer PJ. Identification and characterization of CKLiK, a novel granulocyte Ca++/calmodulin-dependent kinase. Blood. 2000;96(9):3215–23. 11050006

81. Bhattaram P, Penzo-Méndez A, Sock E, Colmenares C, Kaneko KJ, Vassilev A, et al. Organogenesis relies on SoxC transcription factors for the survival of neural and mesenchymal progenitors. Nat Commun. 2010;1(1):1–12. doi: 10.1038/ncomms1008 20596238

82. Jay P, Gozé C, Marsollier C, Taviaux S, Hardelin JP, Koopman P, et al. The human SOX11 Gene: Cloning, chromosomal assignment and tissue expression. Genomics. 1995;29(2):541–5. doi: 10.1006/geno.1995.9970 8666406

83. Haslinger A, Schwarz TJ, Covic M, Chichung Lie D. Expression of Sox11 in adult neurogenic niches suggests a stage-specific role in adult neurogenesis. Eur J Neurosci. 2009;29(11):2103–14. doi: 10.1111/j.1460-9568.2009.06768.x 19490090

84. Bergsland M, Werme M, Malewicz M, Perlmann T, Muhr J. The establishment of neuronal properties is controlled by Sox4 and Sox11. Genes Dev. 2006;20(24):3475–86. doi: 10.1101/gad.403406 17182872

85. Larson BL, Ylostalo J, Lee RH, Gregory C, Prockop DJ. Sox11 is expressed in early progenitor human multipotent stromal cells and decreases with extensive expansion of the cells. Tissue Eng—Part A. 2010;16(11):3385–94. doi: 10.1089/ten.tea.2010.0085 20626275

86. Tsurusaki Y, Koshimizu E, Ohashi H, Phadke S, Kou I, Shiina M, et al. De novo SOX11 mutations cause Coffin-Siris syndrome. Nat Commun. 2014;5:1–7. doi: 10.1038/ncomms5011 24886874

87. Kuryshev VY, Vorobyov E, Zink D, Schmitz J, Rozhdestvensky TS, Münstermann E, et al. An anthropoid-specific segmental duplication on human chromosome 1q22. Genomics. 2006;88(2):143–51. doi: 10.1016/j.ygeno.2006.02.002 16545939

88. Kikuno R, Nagase T, Ishikawa K, Hirosawa M, Miyajima N, TANAKA A, et al. Prediction of the Coding Sequences of Unidentified Human Genes. XIV. The Complete Sequences of 100 New cDNA Clones from Brain Which Code for Large Proteins in vitro. DNA Res. 1999;205:197–205. doi: 10.1093/dnares/6.3.197 10470851

89. Bisogno T, Howell F, Williams G, Minassi A, Cascio MG, Ligresti A, et al. Cloning of the first sn1-DAG lipases points to the spatial and temporal regulation of endocannabinoid signaling in the brain. J Cell Biol. 2003;163(3):463–8. doi: 10.1083/jcb.200305129 14610053

90. Manitt C, Mimee A, Eng C, Pokinko M, Stroh T, Cooper HM, et al. The Netrin Receptor DCC Is Required in the Pubertal Organization of Mesocortical Dopamine Circuitry. J Neurosci. 2011;31(23):8381–94. doi: 10.1523/JNEUROSCI.0606-11.2011 21653843

91. Srour M, Rivière JB, Pham JMT, Dubé MP, Girard S, Morin S, et al. Mutations in DCC cause congenital mirror movements. Science (80-). 2010;328(5978):592. doi: 10.1126/science.1186463 20431009

92. Savage JE, Jansen PR, Stringer S, Watanabe K, Bryois J, De Leeuw CA, et al. Genome-wide association meta-analysis in 269,867 individuals identifies new genetic and functional links to intelligence. Nat Genet. 2018;50(7):912–9. doi: 10.1038/s41588-018-0152-6 29942086

93. Lee JJ, Wedow R, Okbay A, Kong E, Maghzian O, Zacher M, et al. Gene discovery and polygenic prediction from a genome-wide association study of educational attainment in 1.1 million individuals. Nat Genet. 2018;50(8):1112–21. doi: 10.1038/s41588-018-0147-3 30038396

94. Satizabal CL, Adams HHH, Hibar DP, White CC, Knol MJ, Stein JL, et al. Genetic architecture of subcortical brain structures in 38,851 individuals. Nat Genet. 2019;51(11):1624–36. doi: 10.1038/s41588-019-0511-y 31636452

95. Li HJ, Qu N, Hui L, Cai X, Zhang CY, Zhong BL, et al. Further confirmation of netrin 1 receptor (DCC) as a depression risk gene via integrations of multi-omics data. Transl Psychiatry. 2020;10(1). doi: 10.1038/s41398-020-0777-y 32184385

96. Torres-Berrío A, Hernandez G, Nestler EJ, Flores C. The Netrin-1/DCC Guidance Cue Pathway as a Molecular Target in Depression: Translational Evidence. Biol Psychiatry [Internet]. 2020; Available from: doi: 10.1016/j.biopsych.2020.04.025 32593422

97. Wright KM, Rand KA, Kermany A, Noto K, Curtis D, Garrigan D, et al. A prospective analysis of genetic variants associated with human lifespan. G3 Genes, Genomes, Genet. 2019;9(9):2863–78. doi: 10.1534/g3.119.400448 31484785

98. Consortium TASDWG of TPG. Meta-analysis of GWAS of over 16,000 individuals with autism spectrum disorder highlights a novel locus at 10q24.32 and a significant overlap with schizophrenia. Mol Autism. 2017;8(21).

99. Cordell HJ, Han Y, Mells GF, Li Y, Hirschfield GM, Greene CS, et al. International genome-wide meta-analysis identifies new primary biliary cirrhosis risk loci and targetable pathogenic pathways. Nat Commun [Internet]. 2015;6:1–11. Available from: doi: 10.1038/ncomms9019 26394269

100. Boraska V, Franklin CS, Floyd JAB, Thornton LM, Huckins LM, Southam L, et al. A genome-wide association study of anorexia nervosa. Mol Psychiatry. 2014;19(10):1085–94. doi: 10.1038/mp.2013.187 24514567

101. Jung SJ, Winning A, Roberts AL, Nishimi K, Chen Q, Gilsanz P, et al. Posttraumatic stress disorder symptoms and television viewing patterns in the Nurses’ Health Study II: A longitudinal analysis. PLoS One. 2019;14(3):1–13. doi: 10.1371/journal.pone.0213441 30897111

102. Liberzon I, King AP, Ressler KJ, Almli LM, Zhang P, Ma ST, et al. Interaction of the ADRB2 gene polymorphism with childhood trauma in predicting adult symptoms of posttraumatic stress disorder. JAMA Psychiatry. 2014;71(10):1174–82. doi: 10.1001/jamapsychiatry.2014.999 25162199

103. Bierut LJ, Strickland JR, Thompson JR, Afful SE, Cottler LB. Drug use and dependence in cocaine dependent subjects, community-based individuals, and their siblings. Drug Alcohol Depend. 2008;95(1–2):14–22. doi: 10.1016/j.drugalcdep.2007.11.023 18243582

104. Baker DG, Nash WP, Litz BT, Geyer MA, Risbrough VB, Nievergelt CM, et al. Predictors of Risk and Resilience for Posttraumatic Stress Disorder Among Ground Combat Marines: Methods of the Marine Resiliency Study. Prev Chronic Dis. 2012;9(5):1–11. doi: 10.5888/pcd9.110134 22575082

105. Kimbrel NA, Hauser MA, Garrett M, Ashley-Koch A, Liu Y, Dennis MF, et al. Effect of the APOE ε4 allele and combat exposure on PTSD among Iraq/Afghanistan-era veterans. Depress Anxiety. 2015;32(5):307–15. doi: 10.1002/da.22348 25709077

106. Nievergelt CM, Maihofer AX, Klengel T, Atkinson EG, Chen C-Y,…, et al. Largest genome-wide association study for PTSD identifies genetic risk loci in European and African ancestries and implicates novel biological pathways. bioRxiv [Internet]. 2018;(Nov). Available from: http://dx.doi.org/10.1101/458562

107. Liang X, Liu R, Chen C, Ji F, Li T. Opioid System Modulates the Immune Function: A Review Xuan. Transl Perioper Pain Med. 2016;1(1):5–13. 26985446

108. Diasso PDK, Birke H, Nielsen SD, Main KM, Højsted J, Sjøgren P, et al. The effects of long-term opioid treatment on the immune system in chronic non-cancer pain patients: A systematic review. Eur J Pain (United Kingdom). 2020;24(3):481–96. doi: 10.1002/ejp.1506 31705699

109. Mogil JS. Qualitative sex differences in pain processing: emerging evidence of a biased literature. Nat Rev Neurosci 2020 [Internet]. 2020;1–13. Available from: http://www.nature.com/articles/s41583-020-0310-6%0Ahttps://www.nature.com/articles/s41583-020-0310-6 doi: 10.1038/s41583-020-0310-6 32440016

110. König IR, Loley C, Erdmann J, Ziegler A. How to Include Chromosome X in Your Genome-Wide Association Study. Genet Epidemiol. 2014; doi: 10.1002/gepi.21782 24408308

111. Accounting for sex in the genome. Nat Med. 2017;23(11):1243. doi: 10.1038/nm.4445 29117171

112. Macfarlane GJ, Barnish MS, Jones GT. Persons with chronic widespread pain experience excess mortality: longitudinal results from UK Biobank and meta-analysis. Ann Rheum Dis [Internet]. 2017;76(11):1815–22. Available from: http://ard.bmj.com/lookup/doi/10.1136/annrheumdis-2017-211476 28733474

113. Kamaleri Y, Natvig B, Ihlebaek CM, Benth JS, Bruusgaard D. Number of pain sites is associated with demographic, lifestyle, and health-related factors in the general population. Eur J Pain. 2008;12(6):742–8. doi: 10.1016/j.ejpain.2007.11.005 18160318

114. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018;562(7726):203–9. doi: 10.1038/s41586-018-0579-z 30305743

115. Loh PR, Tucker G, Bulik-Sullivan BK, Vilhjálmsson BJ, Finucane HK, Salem RM, et al. Efficient Bayesian mixed-model analysis increases association power in large cohorts. Nat Genet [Internet]. 2015;47(3):284–90. Available from: doi: 10.1038/ng.3190 25642633

116. Willer CJ, Li Y, Abecasis GR, Overall P. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26(17):2190–1. doi: 10.1093/bioinformatics/btq340 20616382

117. Zeisel A, Hochgerner H, Lönnerberg P, Johnsson A, Memic F, van der Zwan J, et al. Molecular Architecture of the Mouse Nervous System. Cell. 2018;174(4):999–1014.e22. doi: 10.1016/j.cell.2018.06.021 30096314

118. Won H, Huang J, Opland CK, Hartl CL, Geschwind DH. Human evolved regulatory elements modulate genes involved in cortical expansion and neurodevelopmental disease susceptibility. Nat Commun [Internet]. 2019;10(1):1–11. Available from: doi: 10.1038/s41467-018-07882-8 30602773

119. Zheng J, Erzurumluoglu AM, Elsworth BL, Kemp JP, Howe L, Haycock PC, et al. LD Hub: A centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics. 2017;33(2):272–9. doi: 10.1093/bioinformatics/btw613 27663502

120. Bulik-Sullivan B, Loh PR, Finucane HK, Ripke S, Yang J, Patterson N, et al. LD score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat Genet [Internet]. 2015;47(3):291–5. Available from: doi: 10.1038/ng.3211 25642630

121. Bulik-sullivan B, Finucane HK, Anttila V, Gusev A, Day FR, Case T, et al. An atlas of genetic correlations across human diseases and traits. Nat Publ Gr [Internet]. 2015;47(11):1236–41. Available from: doi: 10.1038/ng.3406 26414676

122. Dudbridge F. Power and Predictive Accuracy of Polygenic Risk Scores. PLoS Genet. 2013;9(3). doi: 10.1371/journal.pgen.1003348 23555274


Článek vyšel v časopise

PLOS Genetics


2021 Číslo 4
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

Důležitost adherence při depresivním onemocnění
nový kurz
Autoři: MUDr. Eliška Bartečková, Ph.D.

Koncepce osteologické péče pro gynekology a praktické lékaře
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková, Ph.D.

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Multidisciplinární zkušenosti u pacientů s diabetem
Autoři: Prof. MUDr. Martin Haluzík, DrSc., prof. MUDr. Vojtěch Melenovský, CSc., prof. MUDr. Vladimír Tesař, DrSc.

Všechny kurzy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#