A systematic review of genome-wide association studies for pain, nociception, neuropathy, and pain treatment responses

Pain is the leading cause of disability worldwide, imposing an enormous burden on personal health and society. Pain is a multifactorial and multidimensional problem. Currently, there is (some) evidence that genetic factors could partially explain individual susceptibility to pain and interpersonal differences in pain treatment response. To better understand the underlying genetic mechanisms of pain, we systematically reviewed and summarized genome-wide association studies (GWASes) investigating the associations between genetic variants and pain/pain-related phenotypes in humans. We reviewed 57 full-text articles and identified 30 loci reported in more than 1 study. To check whether genes described in this review are associated with (other) pain phenotypes, we searched 2 pain genetic databases, Human Pain Genetics Database and Mouse Pain Genetics Database. Six GWAS-identified genes/loci were also reported in those databases, mainly involved in neurological functions and inflammation. These findings demonstrate an important contribution of genetic factors to the risk of pain and pain-related phenotypes. However, replication studies with consistent phenotype definitions and sufficient statistical power are required to validate these pain-associated genes further. Our review also highlights the need for bioinformatic tools to elucidate the function of identified genes/loci. We believe that a better understanding of the genetic background of pain will shed light on the underlying biological mechanisms of pain and benefit patients by improving the clinical management of pain.</p

A genetic risk score to predict treatment nonresponse in psychotic depression

Psychotic depression is a severe and difficult-to-treat subtype of major depressive disorder for which higher rates of treatment-resistant depression were found. Studies have been performed aiming to predict treatment-resistant depression or treatment nonresponse. However, most of these studies excluded patients with psychotic depression. We created a genetic risk score (GRS) based on a large treatment-resistant depression genome-wide association study. We tested whether this GRS was associated with nonresponse, nonremission and the number of prior adequate antidepressant trials in patients with a psychotic depression. Using data from a randomized clinical trial with patients with a psychotic depression (n = 122), we created GRS deciles and calculated positive prediction values (PPV), negative predictive values (NPV) and odds ratios (OR). Nonresponse and nonremission were assessed after 7 weeks of treatment with venlafaxine, imipramine or venlafaxine plus quetiapine. The GRS was negatively correlated with treatment response (r = −0.32, p = 0.0023, n = 88) and remission (r = −0.31, p = 0.0037, n = 88), but was not correlated with the number of prior adequate antidepressant trials. For patients with a GRS in the top 10%, we observed a PPV of 100%, a NPV of 73.7% and an OR of 52.4 (p = 0.00072, n = 88) for nonresponse. For nonremission, a PPV of 100%, a NPV of 51.9% and an OR of 21.3 (p = 0.036, n = 88) was observed for patients with a GRS in the top 10%. Overall, an increased risk for nonresponse and nonremission was seen in patients with GRSs in the top 40%. Our results suggest that a treatment-resistant depression GRS is predictive of treatment nonresponse and nonremission in psychotic depression.</p

A genetic risk score to predict treatment nonresponse in psychotic depression

Psychotic depression is a severe and difficult-to-treat subtype of major depressive disorder for which higher rates of treatment-resistant depression were found. Studies have been performed aiming to predict treatment-resistant depression or treatment nonresponse. However, most of these studies excluded patients with psychotic depression. We created a genetic risk score (GRS) based on a large treatment-resistant depression genome-wide association study. We tested whether this GRS was associated with nonresponse, nonremission and the number of prior adequate antidepressant trials in patients with a psychotic depression. Using data from a randomized clinical trial with patients with a psychotic depression (n = 122), we created GRS deciles and calculated positive prediction values (PPV), negative predictive values (NPV) and odds ratios (OR). Nonresponse and nonremission were assessed after 7 weeks of treatment with venlafaxine, imipramine or venlafaxine plus quetiapine. The GRS was negatively correlated with treatment response (r = −0.32, p = 0.0023, n = 88) and remission (r = −0.31, p = 0.0037, n = 88), but was not correlated with the number of prior adequate antidepressant trials. For patients with a GRS in the top 10%, we observed a PPV of 100%, a NPV of 73.7% and an OR of 52.4 (p = 0.00072, n = 88) for nonresponse. For nonremission, a PPV of 100%, a NPV of 51.9% and an OR of 21.3 (p = 0.036, n = 88) was observed for patients with a GRS in the top 10%. Overall, an increased risk for nonresponse and nonremission was seen in patients with GRSs in the top 40%. Our results suggest that a treatment-resistant depression GRS is predictive of treatment nonresponse and nonremission in psychotic depression.</p

Drug interaction potential of high-dose rifampicin in patients with pulmonary tuberculosis

Accumulating evidence supports the use of higher doses of rifampicin for tuberculosis (TB) treatment. Rifampicin is a potent inducer of metabolic enzymes and drug transporters, resulting in clinically relevant drug interactions. To assess the drug interaction potential of higher doses of rifampicin, we compared the effect of high-dose rifampicin (40 mg/kg daily, RIF40) and standard-dose rifampicin (10 mg/kg daily, RIF10) on the activities of major cytochrome P450 (CYP) enzymes and P-glycoprotein (P-gp). In this open-label, single-arm, two-period, fixed-order phenotyping cocktail study, adult participants with pulmonary TB received RIF10 (days 1-15), followed by RIF40 (days 16-30). A single dose of selective substrates (probe drugs) was administered orally on days 15 and 30: caffeine (CYP1A2), tolbutamide (CYP2C9), omeprazole (CYP2C19), dextromethorphan (CYP2D6), midazolam (CYP3A), and digoxin (P-gp). Intensive pharmacokinetic blood sampling was performed over 24 hours after probe drug intake. In all, 25 participants completed the study. Geometric mean ratios (90% confidence interval) of the total exposure (area under the concentration versus time curve, RIF40 versus RIF10) for each of the probe drugs were as follows: caffeine, 105% (96%-115%); tolbutamide, 80% (74%-86%); omeprazole, 55% (47%-65%); dextromethorphan, 77% (68%-86%); midazolam, 62% (49%-78%), and 117% (105%-130%) for digoxin. In summary, high-dose rifampicin resulted in no additional effect on CYP1A2, mild additional induction of CYP2C9, CYP2C19, CYP2D6, and CYP3A, and marginal inhibition of P-gp. Existing recommendations on managing drug interactions with rifampicin can remain unchanged for the majority of co-administered drugs when using high-dose rifampicin. Clinical Trials registration number NCT04525235.</p

Methodology of the DCCSS later fatigue study:a model to investigate chronic fatigue in long-term survivors of childhood cancer

Background: A debilitating late effect for childhood cancer survivors (CCS) is cancer-related fatigue (CRF). Little is known about the prevalence and risk factors of fatigue in this population. Here we describe the methodology of the Dutch Childhood Cancer Survivor Late Effect Study on fatigue (DCCSS LATER fatigue study). The aim of the DCCSS LATER fatigue study is to examine the prevalence of and factors associated with CRF, proposing a model which discerns predisposing, triggering, maintaining and moderating factors. Triggering factors are related to the cancer diagnosis and treatment during childhood and are thought to trigger fatigue symptoms. Maintaining factors are daily life- and psychosocial factors which may perpetuate fatigue once triggered. Moderating factors might influence the way fatigue symptoms express in individuals. Predisposing factors already existed before the diagnosis, such as genetic factors, and are thought to increase the vulnerability to develop fatigue. Methodology of the participant inclusion, data collection and planned analyses of the DCCSS LATER fatigue study are presented. Results: Data of 1955 CCS and 455 siblings was collected. Analysis of the data is planned and we aim to start reporting the first results in 2022. Conclusion: The DCCSS LATER fatigue study will provide information on the epidemiology of CRF and investigate the role of a broad range of associated factors in CCS. Insight in associated factors for fatigue in survivors experiencing severe and persistent fatigue may help identify individuals at risk for developing CRF and may aid in the development of interventions.</p

A common germline variant in CYP11B1 is associated with adverse clinical outcome of treatment with abiraterone or enzalutamide

Extragonadal androgens play a pivotal role in prostate cancer disease progression on androgen receptor signaling inhibitors (ARSi), including abiraterone and enzalutamide. We aimed to investigate if germline variants in genes involved in extragonadal androgen synthesis contribute to resistance to ARSi and may predict clinical outcomes on ARSi. We included ARSi naive metastatic prostate cancer patients treated with abiraterone or enzalutamide and determined 18 germline variants in six genes involved in extragonadal androgen synthesis. Variants were tested in univariate and multivariable analysis for the relation with overall survival (OS) and time to progression (TTP) by Cox regression, and PSA response by logistic regression. A total of 275 patients were included. From the investigated genes CYP17A1, HSD3B1, CYP11B1, AKR1C3, SRD5A1 and SRD5A2, only rs4736349 in CYP11B1 in homozygous form (TT), present in 54 patients (20%), was related with a significantly worse OS (HR = 1.71, 95% CI 1.09 – 2.68, p = 0.019) and TTP (HR = 1.50, 95% CI 1.08 – 2.09, p = 0.016), and was related with a significantly less frequent PSA response (OR = 0.48, 95% CI 0.24 – 0.96, p = 0.038) on abiraterone or enzalutamide in a multivariable analysis. The frequent germline variant rs4736349 in CYP11B1 is, as homozygote, an independent negative prognostic factor for treatment with abiraterone or enzalutamide in ARSi naive metastatic prostate cancer patients. Our findings warrant prospective investigation of this potentially important predictive biomarker

Steroid hormone-related polymorphisms associate with the development of bone erosions in rheumatoid arthritis and help to predict disease progression: Results from the REPAIR consortium

Here, we assessed whether 41 SNPs within steroid hormone genes associated with erosive disease. The most relevant finding was the rheumatoid factor (RF)-specific effect of the CYP1B1, CYP2C9, ESR2, FcγR3A, and SHBG SNPs to modulate the risk of bone erosions (P = 0.004, 0.0007, 0.0002, 0.013 and 0.015) that was confirmed through meta-analysis of our data with those from the DREAM registry (P = 0.000081, 0.0022, 0.00074, 0.0067 and 0.0087, respectively). Mechanistically, we also found a gender-specific correlation of the CYP2C9rs1799853T/T genotype with serum vitamin D3 levels (P = 0.00085) and a modest effect on IL1β levels after stimulation of PBMCs or blood with LPS and PHA (P = 0.0057 and P = 0.0058). An overall haplotype analysis also showed an association of 3 ESR1 haplotypes with a reduced risk of erosive arthritis (P = 0.009, P = 0.002, and P = 0.002). Furthermore, we observed that the ESR2, ESR1 and FcγR3A SNPs influenced the immune response after stimulation of PBMCs or macrophages with LPS or Pam3Cys (P = 0.002, 0.0008, 0.0011 and 1.97•10−7). Finally, we found that a model built with steroid hormone-related SNPs significantly improved the prediction of erosive disease in seropositive patients (PRF+ = 2.46•10−8) whereas no prediction was detected in seronegative patients (PRF− = 0.36). Although the predictive ability of the model was substantially lower in the replication population (PRF+ = 0.014), we could confirm that CYP1B1 and CYP2C9 SNPs help to predict erosive disease in seropositive patients. These results are the first to suggest a RF-specific association of steroid hormone-related polymorphisms with erosive disease

Association of response to TNF inhibitors in rheumatoid arthritis with quantitative trait loci for CD40 and CD39

This study was supported by Medical Research Council (MRC) MR/K015346/1 MATURA study. ARUK 20670 MATURA study. Publisher Copyright: © Author(s) (or their employer(s)) 2019.Objectives We sought to investigate whether genetic effects on response to TNF inhibitors (TNFi) in rheumatoid arthritis (RA) could be localised by considering known genetic susceptibility loci for relevant traits and to evaluate the usefulness of these genetic loci for stratifying drug response. Methods We studied the relation of TNFi response, quantified by change in swollen joint counts (ΔSJC) and erythrocyte sedimentation rate (ΔESR) with locus-specific scores constructed from genome-wide assocation study summary statistics in 2938 genotyped individuals: 37 scores for RA; scores for 19 immune cell traits; scores for expression or methylation of 93 genes with previously reported associations between transcript level and drug response. Multivariate associations were evaluated in penalised regression models by cross-validation. Results We detected a statistically significant association between ΔSJC and the RA score at the CD40 locus (p=0.0004) and an inverse association between ΔSJC and the score for expression of CD39 on CD4 T cells (p=0.00005). A previously reported association between CD39 expression on regulatory T cells and response to methotrexate was in the opposite direction. In stratified analysis by concomitant methotrexate treatment, the inverse association was stronger in the combination therapy group and dissipated in the TNFi monotherapy group. Overall, ability to predict TNFi response from genotypic scores was limited, with models explaining less than 1% of phenotypic variance. Conclusions The association with the CD39 trait is difficult to interpret because patients with RA are often prescribed TNFi after failing to respond to methotrexate. The CD39 and CD40 pathways could be relevant for targeting drug therapy.Peer reviewe

Genome-wide association study on pharmacological outcomes of musculoskeletal pain in UK Biobank

The pharmacological management of musculoskeletal pain starts with NSAIDs, followed by weak or strong opioids until the pain is under control. However, the treatment outcome is usually unsatisfying due to inter-individual differences. To investigate the genetic component of treatment outcome differences, we performed a genome-wide association study (GWAS) in ~23,000 participants with musculoskeletal pain from the UK Biobank. NSAID vs. opioid users were compared as a reflection of the treatment outcome of NSAIDs. We identified one genome-wide significant hit in chromosome 4 (rs549224715, P = 3.88 × 10−8). Suggestive significant (P &lt; 1 × 10−6) loci were functionally annotated to 18 target genes, including four genes linked to neuropathic pain processes or musculoskeletal development. Pathway and network analyses identified immunity-related processes and a (putative) central role of EGFR. However, this study should be viewed as a first step to elucidate the genetic background of musculoskeletal pain treatment.</p