CYP gene variants show preliminary evidence of association with patient-reported outcomes in active-duty service members with chronic pain

Abstract

Chronic pain is highly prevalent among active-duty service members (ADSMs) and contributes to substantial physical, psychological, and social burden. Genetic factors are recognized as central to pain susceptibility and treatment response. While few studies have examined the genes involved in cortisol synthesis, changes in the coding of cytochrome P450 (CYP) enzymes responsible for cortisol production have been implicated in chronic pain. We evaluated whether exonic variants in CYP11B1 and CYP17A1 are associated with salivary cortisol and pain-related outcomes in ADSMs with chronic pain using baseline genomic, salivary cortisol, and web survey data from 118 participants in a pragmatic clinical trial. We identified 40 variants across CYP11B1 and CYP17A1 exons and examined both common and rare variant burden scores. Linear regression models estimated associations with 11 outcomes, including pain intensity, sleep-related impairment, depression, and physical and social functioning. Mediation by cortisol was also assessed. Variants in both genes showed consistent associations with pain outcomes. Rare CYP11B1 variants (minor allele frequency <5%) demonstrated the strongest relationships with pain intensity, pain impact, depression, fatigue, and sleep-related impairment. CYP17A1 variants were strongly associated with sleep-related impairment, depression, fatigue, anxiety, and pain impact. Rare CYP11B1 variants were associated with salivary cortisol; however, mediation analyses indicated that cortisol accounted for only a modest proportion (4%–18%) of the gene–pain associations. These exploratory findings suggest that functional variation in cortisol-related CYP genes may contribute to pain-related outcomes among ADSMs with chronic pain, warranting replication in larger cohorts. The research protocol was registered in ClinicalTrials.gov (NCT03297905).

Keywords

Chronic pain genes genotype cortisol patient-reported outcomes military personnel

Introduction

Chronic pain is a major health concern among the adult population worldwide. Compared to the general population, military personnel are at heightened risk, with one-third of active-duty service members (ADSMs) diagnosed with chronic pain.¹ Chronic pain has substantial and varied consequences,² which impact work and physical function as well as mental health and social wellness.^3,4

The multifactorial nature of chronic pain makes evaluating the effectiveness of therapeutic interventions challenging. Several risk factors have been described for chronic pain, including high physical exertion, age, sex, chronic medical conditions, and neurological issues.⁵ More recently, studies examining the genetics of pain have revealed loci associated with chronic pain,⁶ the pharmacogenetics of pain therapies,⁷ and allelic differences linked to pain sensitivity.⁸

Cytochrome P450 (CYP) genes are of particular interest in pain research due to their key roles in the metabolism of drugs and hormones.⁹ Of particular note, CYP genes code for the enzymes that synthesize cortisol, known as a biomarker of stress and chronic pain.¹⁰ Two CYP genes in particular—CYP11B1 and CYP17A1—are linked to cortisol production.^11–14 Specifically, CYP17A1 is essential for the conversion of progesterone to the cortisol precursor, while CYP11B1 is essential to the final step of cortisol synthesis.

Evidence from previous studies suggests that chronic pain severity may be negatively associated with cortisol levels,¹⁵ while other work has reported lower cortisol concentrations and a blunted diurnal slope among individuals with chronic pain.¹⁶ Although the relationship between cortisol regulation and chronic pain remains inconclusive, the close interplay between physiological stress and pain responses provides a biologically plausible rationale for examining related genetic pathways.¹⁷

To date, variants within CYP11B1 and CYP17A1 have been linked to hypertension,¹⁸ cardiometabolic traits,^19–21 and alterations in cortisol synthesis. Accordingly, those genetic variation in CYP17A1 and CYP11B1 may also underlie pain-related outcomes. However, there has been little investigation into direct associations between genes responsible for cortisol synthesis and pain outcomes.²²

We hypothesized that variants in genes essential to cortisol synthesis would be associated with measures of chronic pain. Our study examined whether exonic variants in CYP11B1 and CYP17A1 are associated with patient-reported pain outcomes, including pain, depression, anxiety, and other psychosomatic and physical functioning measures, in a sample of ADSMs with chronic pain. We also sought to identify whether exonic variants in CYP11B1 and CYP17A1 are associated with levels of salivary cortisol, and whether salivary cortisol mediates associations between genetic variants and pain outcomes.

Methods

Study design and setting

This was a secondary analysis using pre-intervention (baseline) data from a pragmatic clinical trial,²³ the details of which have been published previously.²⁴ The parent study was conducted at the Interdisciplinary Pain Management Center of the Madigan Army Medical Center in Tacoma, Washington. Participant enrollment occurred between May 2021 and October 2022.

Sample

Study participants were ADSMs referred to the Interdisciplinary Pain Management Center, the majority being members of the United States Army. Participants were aged 18–58 years, had pain that had lasted 3 or more months, and reported a 7-day average pain intensity of 3 or greater on an 11-point scale²⁵ or scored at least 1 standard deviation worse than the mean on the Patient-Reported Outcomes Measurement Information System (PROMIS) Pain Interference or Physical Function score.^26,27 Individuals who had undergone or were planning to undergo surgery within the next 6 months and individuals with psychosis, at high risk of suicide, or meeting Diagnostic and Statistical Manual of Mental Disorders, fifth edition (DSM-5) criteria for severe substance use disorder were not eligible to participate in the study. A total of 254 participants provided informed consent for genomic and salivary data collection and completed web-based surveys from the Pain Assessment Screening Tool and Outcomes Registry (PASTOR).

Outcome measures

Pain-related outcomes

Each ADSM who participated in the study completed PASTOR assessments at baseline. For this analysis, we considered 11 separate self-reported pain outcome measures from PASTOR: the Defense and Veterans Pain Rating Scale, which assesses average pain intensity over the past 7 days²⁵; the Pain Catastrophizing Scale²⁸; eight PROMIS measures: physical function,²⁷ pain interference,²⁶ fatigue,²⁹ sleep-related impairment,³⁰ satisfaction with social roles and activities (social function),³¹ depression, anxiety, and anger³²; and a pain impact score calculated by compositing pain intensity, pain interference, and physical functioning scores.³³ Each PROMIS domain is reported as a T-score, which is a standardized value calculated from participants’ responses. These T-scores are referenced to the U.S. general population, where the mean is set at 50 and the standard deviation at 10. Higher scores correspond to worse levels of impairment or distress for all measures except for physical functioning and social functioning, for which lower scores indicate worse functioning.

Salivary cortisol

Salivary cortisol samples were also collected at baseline. Samples were collected during weekdays in participants who work during the daytime. Participants were asked to provide three saliva samples prior to treatment for chronic pain: at awakening, 30 min after awakening, and at bedtime. This collection schedule enabled us to analyze both the cortisol awakening response and diurnal variability.^34,35 Participants were asked to collect samples on the day prior to the first treatment day and to turn samples in the following day. Because treatment occurs during the Monday–Friday workweek, cortisol samples may have been collected on any day between Sunday and Thursday. Participants were encouraged to schedule treatment during periods when they were not working night shifts, but shiftwork was not an exclusion criterion. Participants were asked not to consume caffeinated or sugary drinks and/or food, nor engage in physical exercise, before taking their morning samples, to avoid the impacts of these activities on the cortisol awakening response. Participants were instructed to refrigerate samples immediately after collection and to submit the samples to the clinic on the first day of their treatment program. Saliva samples were collected using Salivette Cortisol collection tubes (SARSTEDT, Inc., Newton, NC, USA) following the manufacturer’s instructions. After sample collection, the collection tubes were centrifuged at 1000 relative centrifugal force for 2 min. The saliva samples were then aliquoted and stored at −80 °C for long-term storage. Salivary cortisol was measured using an enzyme-linked immunosorbent assay (ELISA) from R&D Systems (Minneapolis, MN, USA), following the manufacturer’s protocol. Cortisol measurements were used to calculate the area under the curve with respect to ground (AUCg) as previously described by Pruessner et al.³⁶

Exonic CYP variants

We collected baseline CYP samples via buccal swab collection and genomic DNA (gDNA) isolation, using the Qiagen Gentra Puregene Buccal Cell Kit (Germantown, MD, USA) per the manufacturer’s instructions. Participants were advised not to eat or drink for at least 1 h prior to scraping the insides of their cheeks with the buccal brush, 10 times on each side of the mouth. Buccal swabs were temporarily stored at 4 °C until gDNA was extracted and purified to remove proteins and RNAs. Purified gDNA samples were then stored at −20 °C until ready for exome sequencing. Exome sequencing was conducted by the Northwest Genomics Center at the University of Washington. Samples failed initial quality control if the total amount, concentration, or volume was too low, a 63-SNP OpenArray fingerprint assay derived from a custom exome produced poor genotype data or integrity of DNA, or sex-typing was inconsistent with the sample manifest. Base calls were generated on an Illumina NovaSeq6000 instrument (RTA 3.1.5) and a combined suite of Illumina and other industry standard software packages (i.e. Genome Analysis ToolKit (GATK), Picard, BWA, SAMTools, and in-house custom scripts) were used for alignment, local realignment, duplicate removal, quality recalibration, data merging, variant detection, genotyping, and annotation. Exome completion was defined as having >90% of the target at >10× coverage and >80% of the exome target at >20× coverage. Quality control included sample mean coverage of >10×, contamination of <3%, fingerprint concordance >99%, and ensuring that <7% of SNPs were novel variants relative to dbSNP. As an exploratory study using buccal swab collection, the QC standards used were less than current standards.

Of the 254 buccal swab samples provided, 202 passed the initial quality control for library preparation. Of these, quality libraries were generated and sequenced from 160 samples. After sequencing, samples from 118 individuals passed the initial sequencing quality control and were used for genetic analysis. Genomic Variant Call Format (gVCF) files were compressed using CombineGVCFs³⁷ and aligned to a reference FASTA file from HG38. Phred-scaled genotypes were used to infer allele frequencies at each varying site within CYP11B1 and CYP17A1. Genotypes were filtered to include exonic regions of the two genes. There were 40 variants within the two genes that varied across the sample and had inferred genotypes among at least 50 individuals in the sample (see Supplemental Table S1).

Table 1 provides basic information about the four common (minor allele frequency (MAF) >10%) exonic variants in the sample. All four variants (two in CYP11B1 and two in CYP17A1) had MAF >0.35, did not show evidence of being out of Hardy–Weinberg equilibrium (p > 0.05, based on chi-squared goodness-of-fit test using estimated MAF), and are annotated as synonymous coding variants.³⁸ In addition to these variants, data were available on an additional 25 variants in CYP11B1 with MAF <10%, 21 of which had MAF <5%. Similarly, an additional 11 variants in CYP17A1 had available genotypic data, all of which had MAF <5%. All single nucleotide polymorphisms (SNPs) were coded such that the major allele was coded as 0 and the minor allele was coded as 1, so that each individual received a score of 0, 1, or 2 for each variant, representing the number of minor alleles.

Table 1.

Commonly occurring exonic variants in CYP11B1 and CYP17A1.

Gene (chromosome number)	Genome location (GRCh38)	SNP ID	MAF in this sample	Hardy–Weinberg equilibrium p value	Type of exonic variant
CYP11B1 (8)	142879181	rs5283	0.35	0.93	Synonymous
CYP11B1 (8)	142879589	rs6410	0.46	0.75	Synonymous
CYP17A1 (10)	102837167	rs6163	0.41	0.15	Synonymous
CYP17A1 (10)	102837224	rs6162	0.44	0.52	Synonymous

MAF: minor allele frequency; SNP: single nucleotide polymorphism.

We created “burden” style measures of genetic variation separately within each gene by summing the total minor allele count within a gene, representing an additive disease risk model. We created three burden scores for CYP11B1 (MAF <5%: mean 0.4 (SD = 0.8); MAF <10%: mean 0.9 (SD = 1.3); all SNPs: mean 2.4 (SD = 1.6)) and two burden scores for CYP17A1 (MAF <5%: mean 0.3 (SD = 0.6); all SNPs: mean 1.8 (SD = 1.5)).

The UCSC Genome Browser³⁹ provides integrated information about GRCh38/hg38 including each genes start-stop positions, exon start–stop locations and data on which SNPs from dbSNP³⁸ have MAF >1%. These resources indicate a total of seven variants (one missense; six synonymous) with MAF >1% in the exonic regions of CYP11B1 and two variants (both synonymous) in the exonic regions of CYP17A1. All of these variants were observed in the sample with observed MAF close to global population allele frequencies (Table 1 and Supplemental Table S1b).

Statistical analysis

Participant characteristics at baseline were analyzed using descriptive statistics. Frequencies with percentages were computed for categorical data, and means with standard deviations were computed for continuous data. Separate linear regression models were fit, predicting baseline PASTOR measures by commonly occurring (MAF >10%) variants in CYP11B1 and, separately, CYP17A1 or by using quantitative measures of genetic rare variant “burden” (described above).

Additional models adjusted for participants’ baseline demographics (age, sex, race, ethnicity, smoking, income, education, and marriage/partnership). We followed a similar approach to predict log-transformed salivary cortisol levels. Finally, we added salivary cortisol to models predicting pain-related outcomes (PASTOR measures) and assessed the percentage change in genetic variant coefficients as a measure of salivary cortisol’s potential mediational impact on CYP genetic variants and pain. Since this was exploratory research we used 90% confidence intervals, considered p < 0.05 to be a statistically significant and p value’s between 0.05 and 0.10 identified as notable trends.

Results

Sample characteristics

A total of 118 participants were included in the study. Table 2 provides an overview of the participant characteristics. The sample was primarily male (79%), between the ages of 25 and 44 years (85%), and White (66%). More than two-thirds (70%) of the sample had attended at least some college, and nearly two-thirds (63%) had a spouse or partner.

Table 2.

Participant characteristics (N = 118).

Characteristics	% (n) or mean (SD)
Sex
Female	21% (25)
Age (years)
18–24	6% (7)
25–34	48% (57)
35–44	36% (43)
45+	9% (11)
Race
White	66% (78)
Asian	11% (13)
Black	20% (24)
Other	3% (3)
Hispanic ethnicity	16% (19)
Income ($1000/year)
<20	35% (41)
20–49	40% (47)
50–99	14% (17)
100+	11% (13)
Education
High school graduate/GED	25% (30)
Some college/technical degree/AA	46% (54)
College degree (BA/BS)	20% (24)
Advanced degree (MA/PhD/MD)	4% (5)
Married/partnered—yes	63% (74)
Current tobacco use—yes	17% (20)
PROMIS measures
Physical function	38.8 (4.7)
Pain interference	64.1 (5.8)
Fatigue	59.1 (8.6)
Sleep-related impairment	61.9 (8.6)
Satisfaction with social roles and activities	37.6 (6.3)
Depression	55.2 (10.8)
Anxiety	56.9 (10.4)
Anger	54.8 (12.0)
Pain catastrophizing	22.1 (14.27)
Pain intensity	5.4 (1.5)
Pain impact score	29.6 (6.9)
Salivary cortisol (log-transformed)	2.2 (0.8)

AA: associate of arts; BA: bachelor of arts; BS: bachelor of science; GED: general equivalency diploma; MA: master of arts; MD: doctor of medicine; PhD: doctor of philosophy; PROMIS: Patient-Reported Outcomes Measurement Information System.

CYP variants and pain-related outcomes

We examined each of the 11 pain outcome measures to determine whether there was evidence of association with variants in the CYP genes. Five of the 11 outcomes showed statistically significant evidence of association (p < 0.05) in the adjusted models (Tables 3 –7); while all 11 showed at least a notable trend (p < 0.10). Specifically, pain intensity and pain impact both showed associations with variants in the CYP11B1 gene, with the strongest evidence for the burden test including all SNPs with MAF <5%. Fatigue, depression, and sleep-related impairment also showed associations with variants in CYP11B1, particularly for the results of the burden test. The more common variants in CYP17A1 showed associations with pain impact, fatigue, depression, and sleep-related impairment. In a follow-up analysis, we added rs5283 to the pain intensity model with the MAF <5% burden score, and both rs5283 and the burden score remained statistically significant with similar coefficients. Results from other outcomes are shown in Supplemental Table S2.

Table 3.

Pain intensity (N = 118).

Gene	SNP/collapsing test	Unadjusted beta (90% CI)	Adjusted beta (90% CI)
CYP11B1	rs5283	−0.48 (−0.81 to −0.14)**	−0.43 (−0.83 to −0.04)*
	rs6410	−0.16 (−0.50 to 0.19)	−0.20 (−0.59 to 0.20)
	Burden MAF <5%	−0.29 (−0.57 to −0.02)*	−0.45 (−0.76 to −0.14)**
	Burden MAF <10%	−0.06 (−0.23 to 0.11)	−0.23 (−0.44 to −0.02)*
	All SNPs	−0.07 (−0.21 to 0.07)	−0.18 (−0.34 to −0.01)*
CYP17A1	rs6163	−0.20 (−0.50 to 0.10)	−0.23 (−0.59 to 0.13)
	rs6162	−0.13 (−0.44 to 0.18)	−0.17 (−0.54 to 0.20)
	Burden MAF <5%	−0.22 (−0.59 to 0.15)	−0.13 (−0.55 to 0.30)
	All SNPs	−0.04 (−0.19 to 0.11)	−0.01 (−0.19 to 0.17)

CI: confidence interval; MAF: minor allele frequency; SNP: single nucleotide polymorphism.

The adjusted models were estimated using multivariable linear regression, adjusting for age, sex, race, ethnicity, smoking, income, education, and marriage/partnership. The unadjusted models were estimated using univariable linear regression.

p < 0.10. **p < 0.05.

Table 4.

Pain impact score (N = 118).

Gene	SNP/collapsing test	Unadjusted beta (90% CI)	Adjusted beta (90% CI)
CYP11B1	rs5283	−1.16 (−2.91 to 0.59)	−0.14 (−2.16 to 1.87)
	rs6410	−0.32 (−2.10 to 1.46)	−0.11 (−2.10 to 1.89)
	Burden MAF <5%	−1.41 (−2.72 to −0.10)*	−2.03 (−3.41 to −0.65)**
	Burden MAF <10%	−0.29 (−1.09 to 0.52)	−1.07 (−2.01 to −0.13)*
	All SNPs	−0.21 (−0.87 to 0.46)	−0.48 (−1.23 to 0.27)
CYP17A1	rs6163	−1.57 (−3.10 to −0.03)*	−1.83 (−3.58 to −0.08)*
	rs6162	−1.00 (−2.62 to 0.61)	−1.11 (−2.99 to 0.77)
	Burden MAF <5%	−1.02 (−2.77 to 0.73)	−0.76 (−2.66 to 1.13)
	All SNPs	−0.55 (−1.27 to 0.17)	−0.51 (−1.32 to 0.30)

CI: confidence interval; MAF: minor allele frequency; SNP: single nucleotide polymorphism.

p < 0.10. **p < 0.05.

Table 5.

PROMIS fatigue (N = 118).

Gene	SNP/collapsing test	Unadjusted beta (90% CI)	Adjusted beta (90% CI)
CYP11B1	rs5283	−0.89 (−3.17 to 1.38)	0.28 (−2.16 to 2.71)
	rs6410	−1.46 (−3.72 to 0.80)	−1.42 (−3.79 to 0.94)
	Burden MAF <5%	−0.68 (−2.35 to 0.99)	−1.49 (−3.19 to 0.21)
	Burden MAF <10%	−0.22 (−1.23 to 0.79)	−1.24 (−2.38 to −0.11)*
	All SNPs	−0.34 (−1.17 to 0.50)	−0.73 (−1.63 to 0.16)
CYP17A1	rs6163	−1.93 (−3.85 to −0.01)*	−2.69 (−4.77 to −0.61)**
	rs6162	−1.90 (−3.91 to 0.12)	−2.30 (−4.5 to −0.11)*
	Burden MAF <5%	−1.38 (−3.58 to 0.82)	−1.20 (−3.48 to 1.09)
	All SNPs	−0.91 (−1.80 to −0.01)*	−0.95 (−1.92 to 0.02)

CI: confidence interval; MAF: minor allele frequency; SNP: single nucleotide polymorphism.

p < 0.10. **p < 0.05.

Table 6.

PROMIS depression (N = 118).

Gene	SNP/collapsing test	Unadjusted beta (90% CI)	Adjusted beta (90% CI)
CYP11B1	rs5283	−3.28 (−6.04 to −0.51)*	−1.75 (−4.85 to 1.35)
	rs6410	−2.94 (−5.76 to −0.13)*	−3.00 (−6.02 to 0.02)
	Burden MAF <5%	−1.53 (−3.60 to 0.54)	−2.64 (−4.81 to −0.48)**
	Burden MAF <10%	−0.05 (−1.31 to 1.22)	−1.07 (−2.54 to 0.41)
	All SNPs	−0.76 (−1.79 to 0.28)	−1.09 (−2.24 to 0.06)
CYP17A1	rs6163	−2.62 (−5.05 to −0.20)*	−3.60 (−6.24 to −0.96)**
	rs6162	−2.40 (−4.9 to 0.10)	−3.22 (−5.98 to −0.45)*
	Burden MAF <5%	−0.35 (−3.10 to 2.40)	−1.23 (−4.18 to 1.72)
	All SNPs	−0.96 (−2.08 to 0.16)	−1.26 (−2.51 to −0.01)*

CI: confidence interval; MAF: minor allele frequency; SNP: single nucleotide polymorphism.

p < 0.10. **p < 0.05.

Table 7.

PROMIS sleep-related impairment (N = 118).

Gene	SNP/collapsing test	Unadjusted beta (90% CI)	Adjusted beta (90% CI)
CYP11B1	rs5283	−1.78 (−3.96 to 0.40)	−0.28 (−2.76 to 2.19)
	rs6410	−3.66 (−5.78 to −1.53)***	−3.76 (−6.08 to −1.43)***
	Burden MAF <5%	−1.33 (−2.99 to 0.32)	−2.07 (−3.84 to −0.30)*
	Burden MAF <10%	−0.88 (−1.88 to 0.12)	−1.92 (−3.09 to −0.75)***
	All SNPs	−1.12 (−1.94 to −0.30)**	−1.47 (−2.39 to −0.55)***
CYP17A1	rs6163	−2.57 (−4.41 to −0.73)**	−3.02 (−5.12 to −0.92)**
	rs6162	−2.69 (−4.61 to −0.78)**	−2.85 (−5.04 to −0.65)**
	Burden MAF <5%	−1.19 (−3.39 to 1.01)	−1.29 (−3.69 to 1.11)
	All SNPs	−0.96 (−1.85 to −0.06)*	−0.87 (−1.89 to 0.16)

CI: confidence interval; MAF: minor allele frequency; SNP: single nucleotide polymorphism.

Adjusted: multivariable linear regression adjusting for age, sex, race, ethnicity, smoking, income, education, and marriage/partnership; unadjusted: univariable linear regression.

p < 0.10. **p < 0.05. ***p < 0.01.

CYP variants and salivary cortisol

Salivary cortisol levels showed a statistically significant (p < 0.05) association with less frequent variants in CYP11B1 in an adjusted model (Table 8), and a notable trend (p < 0.10) in the unadjusted model. Common variants in both CYP11B1 and CYP17A1 showed similar levels of association with salivary cortisol yet did not attain statistical significance.

Table 8.

Association of salivary cortisol levels and SNPs in two CYP genes (N = 118).

Gene	SNP/collapsing test	Unadjusted beta (90% CI)	Adjusted beta (90% CI)
CYP11B1	rs5283	0.16 (−0.03 to 0.36)	0.18 (−0.05 to 0.41)
	rs6410	0.15 (−0.05 to 0.35)	0.18 (−0.04 to 0.41)
	Burden MAF <5%	0.16 (0.01–0.30)*	0.19 (0.04–0.35)**
	Burden MAF <10%	0.01 (−0.08 to 0.10)	0.06 (−0.05 to 0.17)
	All SNPs	0.02 (−0.05 to 0.10)	0.06 (−0.02 to 0.15)
CYP17A1	rs6163	0.07 (−0.11 to 0.24)	0.17 (−0.04 to 0.37)
	rs6162	0.04 (−0.14 to 0.22)	0.12 (−0.09 to 0.33)
	Burden MAF <5%	0.13 (−0.06 to 0.33)	0.13 (−0.08 to 0.35)
	All SNPs	0.02 (−0.06 to 0.10)	0.07 (−0.02 to 0.16)

CI: confidence interval; CYP: cytochrome P450; MAF: minor allele frequency; SNP: single nucleotide polymorphism.

Adjusted: multivariable linear regression adjusting for age, sex, race, ethnicity, smoking, income, education, and marriage/partnership; unadjusted: univariable linear regression.

p < 0.10. **p < 0.05.

Mediation by salivary cortisol

Table 9 gives the results of a mediation analysis, which evaluated the extent to which the associations between the CYP genetic variants and patient-reported pain outcomes were mediated by salivary cortisol levels. The focus of this analysis was on CYP11B1 rare variants (MAF <5%), which had shown the strongest overall evidence of association with pain outcomes. None of the statistically significant associations between CYP11B1 genetic variants and pain-related outcomes became nonsignificant after accounting for cortisol levels. Furthermore, the percent mediation estimates ranged from 4.5% to 18.1%, representing only modest mediation.

Table 9.

Mediation of CYP11B1 rare variant burden test relationships with pain outcomes, by salivary cortisol (N = 118).

Pain-related outcomes	Model 1, adjusted beta (90% CI)	Model 2, adjusted beta (90% CI)	Percent mediation
PROMIS
Physical function	0.93 (−0.02 to 1.89)	0.77 (−0.02 to 1.89)	18.1
Pain interference	−1.57 (−2.76 to −0.37)**	−1.47 (−2.76 to −0.37)*	6.5
Fatigue	−1.49 (−3.19 to 0.21)	−1.23 (−3.19 to 0.21)	17.4
Sleep-related impairment	−2.07 (−3.84 to −0.30)*	−1.86 (−3.84 to −0.30)*	10.1
Satisfaction with social roles and activities	1.60 (0.31–2.90)**	1.33 (0.31–2.90)*	17.1
Depression	−2.64 (−4.81 to −0.48)**	−2.22 (−4.81 to −0.48)*	15.8
Anxiety	−2.07 (−4.22 to 0.07)	−1.72 (−4.22 to 0.07)	17.1
Anger	−2.01 (−4.37 to 0.34)	−1.92 (−4.37 to 0.34)	4.5
Pain catastrophizing	−2.65 (−5.51 to 0.22)	−2.51 (−5.51 to 0.22)	5.1
Pain intensity	−0.45 (−0.76 to −0.14)**	−0.39 (−0.76 to −0.14)**	13.7
Pain impact	−2.03 (−3.41 to −0.65)**	−1.80 (−3.41 to −0.65)**	11.2

CI: confidence interval; MAF: minor allele frequency; PROMIS: Patient-Reported Outcomes Measurement Information System; SNP: single nucleotide polymorphisms.

Model 1: adjusted for age, sex, race, ethnicity, smoking, income, education, and marriage/partnership; Model 2: adjusted for cortisol in addition to Model 1 variables.

p < 0.10. **p < 0.05.

Discussion

We identified a clear and consistent signal between certain CYP genetic variants and a variety of patient-reported pain outcomes in our study sample of individuals with chronic pain. In our analysis, we found that SNPs in both CYP11B1 and CYP17A1 were associated with pain-related outcomes. Interestingly, these associations spanned all 11 patient-reported outcomes we measured, were found with both CYP genes we investigated, and included both more common and less common genetic variants.

CYP11B1 and pain-related outcomes

Of the 11 patient-reported outcomes we assessed, four directly relate to pain (pain intensity, impact, catastrophizing, and interference). Three of these showed a statistically significant association with at least one variant or for the burden score in CYP11B1, and the fourth showed a notable trend. In all cases, individuals with more variants in the gene had less pain intensity, impact, catastrophizing, or interference. In fact, nearly all SNPs and burden scores in CYP11B1 showed inverse associations with these four pain measures, even the SNPs and burden scores that did not attain statistical significance or show a notable trend. The two strongest associations were between the MAF <5% burden score and pain intensity and the MAF <5% burden score and pain impact. Given that the two common CYP11B1 variants (rs5283 and rs6410) had associations in a similar direction, albeit weaker, and that both variants are synonymous, it is plausible that any direct relationships between genetic variants in CYP11B1 are due to missense rare variants that are correlated with the common variants. However, our follow-up analysis, which modeled both the common rs5283 and the MAF <5% burden score simultaneously, yielded similar, independent, estimated effects. Bioinformatic analysis of the CYP loci with larger samples to better understand haplotype structure and recombination may provide more insights into possible mechanisms of action.

Variants in CYP11B1 also showed inverse associations with depression, fatigue, sleep-related impairment, and social functioning, with the strongest relationships found for sleep-related impairment. While several common variants were associated with sleep-related impairment, only the rare CYP11B1 variants showed measurable associations with salivary cortisol. This pattern may support the possibility that CYP11B1 is related to sleep via neuroendocrine pathways beyond cortisol concentration alone.⁴⁰ Sleep-related impairment is particularly sensitive to subtle alterations in signaling related to the hypothalamic–pituitary–adrenal axis and broader steroidogenic processes.^41,42 This sensitivity may allow even modest genetic effects to be detected more easily, whether they arise from rare functional variants or from common variants that tag nearby regulatory loci. In addition, sleep-related impairment may be a more sensitive and reliably measured phenotype in the study population of ASDMs, which would potentially explain why the genetic associations appeared strongest for this domain. Previous studies have shown that methylation changes in the CYP11B1 promoter region correlate with overproduction of cortisol,⁴³ but evidence is lacking about direct relationships with pain. Therefore, these interpretations are speculative and require validation in larger studies designed to evaluate mechanistic pathways with deeper genomic coverage.

CYP17A1 and pain-related outcomes

Variants in CYP17A1 also showed evidence of associations with pain-related outcomes, which were in line with emerging findings about SNPs in CYP17A1 and pain in people with cancer.^22,44,45 The strongest inverse relationships were observed with sleep-related impairment, followed by depression, fatigue, anxiety, and pain impact. In contrast to CYP11B1 which is specific to the end stages of cortisol and corticosterone synthesis, CYP17A1 is involved in the generation of neuroactive steroids such as progesterone and dehydroepiandrosterone (DHEA).¹² It is possible that our observed associations relate to these other steroids, in particular DHEA which regulate sleep and stress responses by modulating the activity of gamma-aminobutyric-acid-type-A (GABA(A)) receptors^46,47; a class of ionotropic receptors involved in pain sensation.⁴⁸ Understanding these associations is important because it strengthens emerging evidence that genetic differences in CYP17A1 may influence how people experience multiple pain-related symptoms. Pain symptoms often cluster, yet their biological underpinnings remain poorly understood. Identifying a shared genetic contributor like CYP17A1 could help clarify potential mechanistic pathways, support the development of more targeted symptom management strategies, and move the field closer to approaches that tailor interventions to an individual’s genetic profile.

CYP variants and salivary cortisol

While we expected to find strong and consistent associations between CYP11B1 and CYP17A1 genetic variants and salivary cortisol, our results were limited to an inverse association between salivary cortisol and the rarer variants in CYP11B1 (MAF <5% burden score). This finding aligns with the dearth of prior research about genetic variants in these two CYP genes and their associations with cortisol levels in the body.⁴⁹ It is possible that compensatory mechanisms offset the impact of genetic variations on cortisol production. For example, biochemical responses may increase cortisol synthesis and/or inhibit the reduction of cortisol to its inactive metabolites. However, the inverse association between rarer variants in CYP11B1 and salivary cortisol suggests that some less abundant variants in CYP11B1 may play a larger role in directly impacting cortisol levels.

Epigenetic mechanisms

Epigenetic mechanisms may explain how SNPs, even if synonymous, may regulate cortisol levels. For example, Maharaj et al. observed that synonymous SNPs alter mRNA splicing leading to a nonfunctional CYP11A1 protein.⁵⁰ Although the impact of synonymous SNPs on CYP11B1 transcription remains unknown, synonymous polymorphisms in CYP17A1 are associated with differences in hormone levels in vivo and shown to disrupt genes transcription in vitro.¹⁸ Given the potential impact on CYP transcription, synonymous variants have been associated with adverse health, such as coronary heart disease and skeletal growth disorders.^20,51 However, in many of these studies, it remains unclear whether epigenetic mechanisms increase or decrease cortisol levels and if systemic changes in cortisol promote harmful or protective signaling responses.

Immune and inflammatory response

The inverse relationship with pain and cortisol may be explained by related immune and inflammation responses. Cortisol is a potent immunomodulator reported to instigate both pro- and anti-inflammatory responses.⁵² It is possible that in subjects with chronic pain, higher levels of cortisol levels reflect anti-inflammatory responses. This possibility is supported by the anti-inflammatory and immunosuppressive effects of synthetic cortisol variants such as dexamethasone, that are commonly used for pain management.⁵³ Given the central role of inflammation in chronic pain, depression and poor sleep, immune response to cortisol may relate to our findings. Future analysis of inflammatory cytokines in conjunction with cortisol are needed to clarify this potential relationship.

Mediation role of salivary cortisol

This study’s hypothesis was based on the reasoning that cortisol would underlie the mechanisms by which associations between genetic variants and pain-related outcomes operate. In other words, that genetic variation in CYP genes would impact cortisol levels, which would in turn impact patient-reported pain measures. However, in this study’s sample of ADSMs, estimates of mediation were low (<15%–20%), indicating otherwise. A separate manuscript reporting data from this sample has also reported little evidence of associations between salivary cortisol levels and PASTOR measures, further limiting the potential mediational explanation.⁵⁴ Alternative explanations, which are not possible to explore in this dataset, include that SNPs in CYP genes influence pain through mechanisms independent of systemic cortisol (e.g. processing other steroid precursors), or that cortisol does mediate pain but due to local cortisol actions or impacts of genetic variation on cortisol function, but not on cortisol levels.

Strengths and limitations

Our study offers a novel and necessary look at potential genetic mechanisms of pain in a sample of individuals largely free of confounding health complications such as hypertension and cardiovascular disease and that has a comprehensive collection of paired salivary cortisol and genetic analyses. Our experimental design and analysis of salivary cortisol provide several advantages by normalizing diurnal fluctuations that can impact cortisol concentrations. Another major advantage of analyzing salivary cortisol is it reflects free (unbound) cortisol levels.⁵⁵ Within the blood the majority of cortisol is bound to albumin or cortisol binding protein (CBP) and unavailable for signaling. In contrast, salivary cortisol levels reflect physiological concentrations that are available for signaling. Nevertheless, heterogeneity in sample days, times or participant work/sleep patterns could contribute to variation in cortisol levels. This study has several important limitations. First, the modest sample size (N = 118) limits statistical power, particularly for detecting small effect sizes and rare variant associations. Second, we interpreted p < 0.10 as notable trends, and did not adjust for multiple testing, given the exploratory nature of the analyses and the limited sample size. While this approach may enhance sensitivity to potential associations, it also increases the risk of type I error. Third, the study population was consisted of ADSMs with relatively homogeneous demographics, predominantly male and White. This demographic homogeneity limits the external validity of the findings. Consequently, the generalizability of the finding to women, diverse racial/ethnic groups, or civilian populations remains uncertain. Accordingly, the reported results should be interpreted as largely preliminary rather than statistically supported associations. Replication in a demographically diverse population with larger sample sizes is necessary to evaluate the robustness and reproducibility of these findings.

Conclusions

Our exploratory findings suggest that functional variation in cortisol-related CYP genes may be associated with pain-related and psychosocial outcomes among ADSMs with chronic pain. Future work should explore additional mechanistic pathways by which the putative associations we have identified may exist (e.g. chronic conditions, post-transcriptional CYP protein regulation).

Supplemental Material

sj-docx-1-mpx-10.1177_17448069261454823 – Supplemental material for CYP gene variants show preliminary evidence of association with patient-reported outcomes in active-duty service members with chronic pain

Supplemental material, sj-docx-1-mpx-10.1177_17448069261454823 for CYP gene variants show preliminary evidence of association with patient-reported outcomes in active-duty service members with chronic pain by Nathan L Tintle, Nathan A Ryder, Jason Westra, Dahee Wi, Sotaro Shimada, Nicholas Ieronimakis, Honor M McQuinn, Tyler J Snow, Ardith Z Doorenbos and Diane M Flynn in Molecular Pain

Footnotes

Acknowledgements

We acknowledge the editorial support of Kyra Freestar and Rica Asuncion-Reed, who contributed to making the manuscript clearer.

Consent to participate

All participants in this study provided signed consent for their participation in the study.

Consent for publication

Not applicable.

Data availability statement

Any request for original data sets will be honored upon approval of a data-sharing agreement.

Declaration of conflicting interests

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

Disclaimer

The views expressed in this research article are those of the author and do not necessarily reflect the official policy or position of the Defense Health Agency, Department of War, nor the U.S. Government. The study protocol was approved by the Madigan Army Medical Center Institutional Review Board in compliance with all applicable Federal regulations governing the protection of human subjects.

Funding

The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The U.S. Army Medical Research Acquisition Activity (820 Chandler Street, Fort Detrick, MD 21702-5014) is the awarding and administering acquisition office. This work was supported by the Assistant Secretary of War for Health Affairs and Department of War through the Neuromusculoskeletal Injuries Rehabilitation Research Award (award no. W81XWH-18-2-0023) and National Institutes of Health/National Institute of Neurological Disorders and Stroke (award no. K24 AT011995).

Ethical considerations

The research protocol of the present study was approved by the Madigan Army Medical Center’s ethical review committee (approval number 221011) and registered in ClinicalTrials.gov (NCT03297905) prospectively.

ORCID iD

Nathan L Tintle

Supplemental material

Supplemental material for this article is available online.

References

Reif

Adams

Ritter

Williams

Larson

. Prevalence of pain diagnoses and burden of pain among active duty soldiers, FY2012. Mil Med 2018; 183(9–10): e330–e337.

Fine

. Long-term consequences of chronic pain: mounting evidence for pain as a neurological disease and parallels with other chronic disease states. Pain Med 2011; 12(7): 996–1004.

Smith

Taubman

Clark

. A burden and prevalence analysis of chronic pain by distinct case definitions among active duty U.S. military service members, 2018. Pain Physician 2020; 23(5): E429–E440.

Dunbar

Schuler

Meadows

Engel

. Associations between mental and physical health conditions and occupational impairments in the U.S. military. Mil Med 2022; 187(3–4): e387–e393.

Mills

SEE

Nicolson

Smith

. Chronic pain: a review of its epidemiology and associated factors in population-based studies. Br J Anaesth 2019; 123(2): e273–e283.

Suri

Palmer

Tsepilov

Freidin

Boer

Yau

Evans

Gelemanovic

Bartz

Nethander

Arbeeva

Karssen

Neogi

Campbell

Mellstrom

Ohlsson

Marshall

Orwoll

Uittterlinden

Rotter

Lauc

Psaty

Karlsson

Lane

Jarvik

Polasek

Hochberg

Jordan

Van Meurs

JBJ

Jackson

Nielson

Mitchell

Smith

Hayward

Smith

Aulchenko

Williams

FMK

. Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain. PLoS Genet 2018; 14(9): e1007601.

Angst

Phillips

Drover

Tingle

Ray

Swan

Lazzeroni

Clark

. Pain sensitivity and opioid analgesia: a pharmacogenomic twin study. Pain 2012; 153(7): 1397–1409.

Nielsen

Stubhaug

Price

Vassend

Czajkowski

Harris

. Individual differences in pain sensitivity: genetic and environmental contributions. Pain 2008; 136(1): 21–29.

Aroke

. Pharmacogenetics of postoperative pain management: a review. AANA J 2020; 88(3): 229–236.

10.

Vachon-Presseau

Roy

Martel

Caron

Marin

Chen

Albouy

Plante

Sullivan

Lupien

Rainville

. The stress model of chronic pain: evidence from basal cortisol and hippocampal structure and function in humans. Brain 2013; 136(3): 815–827.

11.

Yadav

Petrunak

Estrada

Scott

. Structural insights into the function of steroidogenic cytochrome P450 17A1. Mol Cell Endocrinol 2017; 441: 68–75.

12.

Burris-Hiday

Scott

. Steroidogenic cytochrome P450 17A1 structure and function. Mol Cell Endocrinol 2021; 528: 111261.

13.

Lisurek

Bernhardt

. Modulation of aldosterone and cortisol synthesis on the molecular level. Mol Cell Endocrinol 2004; 215(1–2): 149–159.

14.

Brixius-Anderko

Scott

. Structure of human cortisol-producing cytochrome P450 11B1 bound to the breast cancer drug fadrozole provides insights for drug design. J Biol Chem 2019; 294(2): 453–460.

15.

Trevino

Geier

Morris

Cronn

de Roon-Cassini

. Relationship between decreased cortisol and development of chronic pain in traumatically injured. J Surg Res 2022; 270: 286–292.

16.

Generaal

Vogelzangs

Macfarlane

Geenen

Smit

Penninx

BWJH

Dekker

. Reduced hypothalamic–pituitary–adrenal axis activity in chronic multi-site musculoskeletal pain: partly masked by depressive and anxiety disorders. BMC Musculoskelet Disord 2014; 15(1): 227.

17.

Hannibal

Bishop

. Chronic stress, cortisol dysfunction, and pain: a psychoneuroendocrine rationale for stress management in pain rehabilitation. Phys Ther 2014; 94(12): 1816–1825.

18.

Diver

MacKenzie

Fraser

McManus

Freel

Alvarez-Madrazo

McClure

Friel

Hanley

Dominiczak

Caulfield

Munroe

Connell

Davies

. Common polymorphisms at the CYP17A1 locus associate with steroid phenotype: support for blood pressure genome-wide association study signals at this locus. Hypertension 2016; 67(4): 724–732.

19.

Zhang

Wang

Ouyang

Deng

Liu

Zhang

. Polymorphisms in CYP11B2 and CYP11B1 genes associated with primary hyperaldosteronism. Hypertens Res 2010; 33(5): 478–484.

20.

Huang

Cheng

Wang

. Genetic variants in CYP11B1 influence the susceptibility to coronary heart disease. BMC Med Genomics 2022; 15(1): 158.

21.

Liu

Duan

. CYP11B1 gene polymorphisms and susceptibility to ischemic stroke in a Chinese Han population. Front Neurosci 2022; 16: 1030551.

22.

Garcia-Giralt

Rodríguez-Sanz

Prieto-Alhambra

Servitja

Torres-Del Pliego

Balcells

Albanell

Grinberg

Diez-Perez

Tusquets

Nogues

. Genetic determinants of aromatase inhibitor-related arthralgia: the B-ABLE cohort study. Breast Cancer Res Treat 2013; 140(2): 385–395.

23.

Doorenbos

Kato

. Pragmatic clinical trials: increasing the rate of translating nursing research into practice^†. Jpn J Nurs Sci 2020; 17(3): e12350.

24.

Flynn

Eaton

Langford

Ieronimakis

McQuinn

Burney

Holmes

Doorenbos

. A SMART design to determine the optimal treatment of chronic pain among military personnel. Contemp Clin Trials 2018; 73: 68–74.

25.

Polomano

Galloway

Kent

Brandon-Edwards

Kwon

Morales

Buckenmaier

. Psychometric testing of the Defense and Veterans Pain Rating Scale (DVPRS): a new pain scale for military population. Pain Med 2016; 17(8): 1505–1519.

26.

Amtmann

Cook

Jensen

Chen

Choi

Revicki

Cella

Rothrock

Keefe

Callahan

Lai

. Development of a PROMIS item bank to measure pain interference. Pain 2010; 150(1): 173–182.

27.

Schalet

Hays

Jensen

Beaumont

Fries

Cella

. Validity of PROMIS physical function measured in diverse clinical samples. J Clin Epidemiol 2016; 73: 112–118.

28.

Sullivan

MJL

Bishop

Pivik

. The Pain Catastrophizing Scale: development and validation. Psychol Assess 1995; 7: 524–532.

29.

Cella

Lai

Jensen

Christodoulou

Junghaenel

Reeve

Stone

. PROMIS fatigue item bank had clinical validity across diverse chronic conditions. J Clin Epidemiol 2016; 73: 128–134.

30.

Buysse

Germain

Moul

Stover

Dodds

Johnston

Pilkonis

. Development of short forms from the PROMIS™ sleep disturbance and sleep-related impairment item banks. Behav Sleep Med 2011; 10(1): 6–24.

31.

Hahn

Beaumont

Pilkonis

Garcia

Magasi

DeWalt

Cella

. The PROMIS satisfaction with social participation measures demonstrated responsiveness in diverse clinical populations. J Clin Epidemiol 2016; 73: 135–141.

32.

Pilkonis

Choi

Reise

Stover

Riley

Cella

. Item banks for measuring emotional distress from the Patient-Reported Outcomes Measurement Information System (PROMIS^®): depression, anxiety, and anger. Assessment 2011; 18(3): 263–283.

33.

Deyo

Ramsey

Buckley

Michaels

Kobus

Eckstrom

Forro

Morris

. Performance of a Patient Reported Outcomes Measurement Information System (PROMIS) short form in older adults with chronic musculoskeletal pain. Pain Med 2016; 17(2): 314–324.

34.

Stalder

Lupien

Kudielka

Adam

Pruessner

Wust

Dockray

Smyth

Evans

Cirschbaum

Miller

Wetherell

Pinke

Klucken

Clow

. Evaluation and update of the expert consensus guidelines for the assessment of the cortisol awakening response (CAR). Psychoneuroendocrinology 2022; 146: 105946.

35.

Taylor

Hernández

Fuller

Sargent

Padilla

Harris

. Cortisol awakening response in elite military men: summary parameters, stability measurement, and effect of compliance. Mil Med 2016; 181(11): e1600–e1607.

36.

Pruessner

Kirschbaum

Meinlschmid

Hellhammer

. Two formulas for computation of the area under the curve represent measures of total hormone concentration versus time-dependent change. Psychoneuroendocrinology 2003; 28(7): 916–931.

37.

GATK Team. CombineGVCFs, https://gatk.broadinstitute.org/hc/en-us/articles/360037593911-CombineGVCFs (2022, accessed 17 December 2025).

38.

Phan

Zhang

Wang

Villamarin

Hefferon

Ramanathan

Kattman

. The evolution of dbSNP: 25 years of impact in genomic research. Nucleic Acids Res 2025; 53(D1): D925–D931.

39.

Casper

Speir

Raney

Perez

Nassar

Lee

Hinrichs

Gonzalez

Fischer

Diekhans

Clawson

Benet-Pages

Barber

Vaske

van Baren

Wang

Rodriguez

YJP

Jenkins-Kiefer

Chalamala

Haussler

Kent

Haeussler

. The UCSC Genome Browser database: 2026 update. Nucleic Acids Res 2026; 54(D1): D1331–D1335.

40.

Frau

Traccis

Bortolato

. Neurobehavioural complications of sleep deprivation: shedding light on the emerging role of neuroactive steroids. J Neuroendocrinol 2020; 32(1): e12792.

41.

Jiao

Butoyi

Zhang

Adotey

SAAI

Chen

Shen

Wang

Yuan

Jia

. Sleep disorders impact hormonal regulation: unravelling the relationship among sleep disorders, hormones and metabolic diseases. Diabetol Metab Syndr 2025; 17(1): 305.

42.

O’Byrne

Yuen

Butt

Liu

. Sleep and circadian regulation of cortisol: a short review. Curr Opin Endocr Metab Res 2021; 18: 178–186.

43.

Cheng

Pai

. Regulation of human CYP11B1 and CYP11B2 promoters by transposable elements and conserved cis elements. Steroids 2012; 77(1–2): 100–109.

44.

Yang

Barnes

Lyon

Dorsey

. Genetic variants associated with cancer pain and response to opioid analgesics: implications for precision pain management. Semin Oncol Nurs 2019; 35(3): 291–299.

45.

Servitja

Martos

Rodriguez Sanz

Garcia-Giralt

Prieto-Alhambra

Garrigos

Nogues

Tusquets

. Skeletal adverse effects with aromatase inhibitors in early breast cancer: evidence to date and clinical guidance. Ther Adv Med Oncol 2015; 7(5): 291–296.

46.

Genud

Merenlender

Gispan-Herman

Maayan

Weizman

Yadid

. DHEA lessens depressive-like behavior via GABA-ergic modulation of the mesolimbic system. Neuropsychopharmacology 2009; 34(3): 577–584.

47.

Doan

DeYoung

Fuller-Rowell

Liu

Meyer

. Investigating relations among stress, sleep and nail cortisol and DHEA. Stress 2018; 21(2): 188–193.

48.

Qian

Zhao

Yin

Zhang

Wang

Zhu

Luo

. Current status of GABA receptor subtypes in analgesia. Biomed Pharmacother 2023; 168: 115800.

49.

Velders

Kuningas

Kumari

Dekker

Uitterlinden

Kirschbaum

Hek

Hofman

Verhulst

Kivimaki

Van Duijn

Walker

Tiemeier

. Genetics of cortisol secretion and depressive symptoms: a candidate gene and genome wide association approach. Psychoneuroendocrinology 2011; 36(7): 1053–1061.

50.

Maharaj

Buonocore

Meimaridou

Ruiz Babot

Guasti

Peng

Capper

Burgos-Tirado

Prasad

Hughes

Maudhoo

Crowne

Cheetham

Brain

Suntharalingham

Striglinoi

Yuksel

Gurbuz

Gupta

Lindsay

Couch

Spoudeas

Guran

Johnson

Fowler

Conwell

McInerney-Leo

Duir

Cariou

Lopez-Siguero

Harris

Duncan

Hindmarsh

Auchus

Donaldson

Achermann

Metherell

. Predicted benign and synonymous variants in CYP11A1 cause primary adrenal insufficiency through missplicing. J Endocr Soc 2019; 3(1): 201–221.

51.

Grgic

Gazzara

Chesi

Medina-Gomez

Cousminer

Mitchell

Prijatelj

de Vris

Shevroja

McCormack

Kalkwark

Lappe

Gilsanz

Oberfield

Shepherd

Kelly

Mahboubi

Faucz

Feelders

de Jong

Uitterlinden

Visser

Ghanem

Wolvius

Hofland

Stratakis

Zemel

Barash

Grant

SFA

Rivadeneira

. CYP11B1 variants influence skeletal maturation via alternative splicing. Commun Biol 2021; 4(1): 1274.

52.

Amasi-Hartoonian

Sforzini

Cattaneo

Pariante

. Cause or consequence? Understanding the role of cortisol in the increased inflammation observed in depression. Curr Opin Endocr Metab Res 2022; 24: 100356.

53.

Wang

Yin

. Efficacy of dexamethasone in reducing pain and inflammation and accelerating total hip arthroplasty postoperative recovery: a randomized controlled trial. J Perianesth Nurs 2024; 39(4): 589–595.

54.

Shimada

Ledesma

Flynn

Tintle

Ransom

McQuinn

Snow

Ieronimakis

Burke

Steffen

Doorenbos

. Evaluating sex differences in biomarkers of chronic pain among active-duty personnel: an exploratory analysis of a pragmatic clinical trial with a SMART design. Biol Res Nurs 2026; 2026: 10998004261423210.

55.

Duplessis

Rascona

Cullum

Yeung

. Salivary and free serum cortisol evaluation. Mil Med 2010; 175(5): 340–346.

Supplementary Material

Please find the following supplemental material available below.

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

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

0.00 MB

0.03 MB