Abstract

Serum adrenal androgens (AAs), including androstenedione (Δ4A) and dehydroepiandrosterone sulfate (DHEAS), have been reported to be lower in female rheumatoid arthritis (RA) patients with early disease. Few data are available on hormonal status of women before the onset of clinical rheumatoid arthritis (pre-RA). A broad baseline panel of serum adrenal and sex steroids was compared in 36 female pre-RA to 144 matched cohort control (CN) subjects to determine differences in their mean values and in patterns of hormonal correlations. Study subjects having lower versus higher baseline serum cortisol levels than the total group's mean value were also analyzed separately to investigate differences in their hormonal levels and correlational patterns. In total subjects, mean (±SE) Δ4A level (nmol/L) was lower ( ) in 28 pre-RA cases ( ) versus 108 CN ( ). The significant ( ) difference was restricted to 9 pre-RA versus 53 CN subjects having lower cortisol levels ( versus  nmol/L, resp.). In total subjects, no significant difference was found between study subjects in their bivariate correlations of the hormonal panel variables, unlike results found in the subgroups stratified by lower versus higher cortisol levels. A subgroup of pre-RA females may have relative adrenal cortical insufficiency, as reflected by lower Δ4A, especially observed among those subjects with lower cortisol levels.

1. Introduction

Relative insufficiency of adrenal glucocorticoid (GC) and androgenic-anabolic (AA) hormones has been suspected to increase the risk of developing rheumatoid arthritis (RA) and to contribute to its multifactorial neuroendocrine immune (NEI) pathogenesis [18]. The characteristic age- and sex-specific incidence patterns of RA support a possible constitutional deficiency of adrenal cortical or sex hormones in a subset of susceptible women. The female-to-male (F : M) risk ratio of persons developing RA is approximately 2 : 1 during the juvenile and older ages but is significantly increased to about 5 : 1 during the female reproductive years [9]. Of further note, risk of RA onset increases with age in adults, particularly among females. The preceding risk data imply that males have relative protection over females during all ages, but particularly in the younger and middle adult years [9]. A recent study indicated that early age at menopause (≤45 years) was associated with the subsequent risk of developing RA [10]. The hazard remained significant after adjusting for smoking, educational level, and length of breastfeeding [10]. Available data imply that the woman’s risk of developing RA may be affected by relative insufficiency of both adrenal cortical GC and AA steroids as well as by accelerated aging of their ovarian hormonal status [110].

A quantitative study of adrenal steroid urinary metabolites from 8 early disease premenopausal RA women not treated with GC drugs and 8 healthy matched controls under varied baseline, physiological adrenal stimulation, and metyrapone GC-inhibiting conditions supported AA deficiency in the patients [11]. Significant differences were found only in the 11-deoxy-17-ketosteroid AA steroids: dehydroepiandrosterone (DHEA), androsterone, and etiocholanolone [11]. Subsequent studies of serum DHEA and its sulfate (DHEAS) in early disease premenopausal onset women found lower mean levels of these AA biomarkers than in matched controls [5, 8, 12, 13]. A nested case-control cohort study showed that definitely low baseline serum DHEAS levels (<0.68 μmol/L), assayed in independent reference laboratories, were present a mean of 11 years before premenopausal onset in 3 (30%) of 10 pre-RA versus 1 (2.7%) of 37 matched CN subjects ( ) [14]. Recently, the same cohort study [15] revealed that definitely low baseline serum cortisol levels had occurred in a greater minority of pre-RA than CN females (11% versus 1.4%, resp., ). In addition, none of 28 pre-RA women had serum androstenedione (Δ4A) levels in the upper normal range (≥7.5 nmol/L) versus 26 (24%) of 108 CN women ( ). The lower levels of urinary AA metabolites, serum DHEA, and DHEAS in premenopausal early RA patients and the recent findings of baseline lower serum cortisol, as well as a truncated upper range of Δ4A levels in a minority of pre-RA versus CN female cohort subjects, have prompted this further analysis of a comprehensive panel of adrenal cortical and sex steroids in our pre-RA and CN cohort subjects [15, 16].

2. Materials and Methods

2.1. The RA Precursors Study (RAPS) Database at This Institution

The RA Precursors Study (RAPS) was initiated at this institution in late 1991, following donation of baseline personal data and serum samples from the pre-RA cases and matched CN cohort subjects by “Operation CLUE I,” a community-wide prospective study [1921]. The 1974 CLUE I base cohort had enrolled 12,381 females of Washington County, Maryland, USA. The RAPS female database currently includes 180 study subjects, 36 Caucasian pre-symptomatic RA (pre-RA) cases at their 1974 baseline entry, and 144 matched cohort CN subjects, in a ratio of 1 pre-RA: 4 CN. The UICOMP Institutional Review Board has approved this research for assurance of confidentiality.

The pre-RA cases in this study conform to The European League Against Rheumatism (EULAR) recommendations [22]. All baseline pre-RA cases were diagnosed and confirmed in the practice of the sole rheumatologist in the cohort community, who used the American College of Rheumatology (ACR) 1987 revised classification criteria [23]. Following the 1974 cohort entry, clinical onsets of RA in our study females occurred within 3 to 18 years (1977 to 1992), after a median of 11 years. None of the matched comparison subjects had a diagnosis of RA in the community rheumatologist’s practice. The non-RA subjects (CN) were matched to the pre-RA cases on race (all Caucasians) and usually within one year of age at entry. Also, the CN subjects were the closest in chronological sequence of enrollment in the cohort, analogous to another case-control study [24].

In 1992, after clinical onset and diagnosis of the first set of baseline pre-RA females, 4 cohort CN subjects were matched to each female case on entry age and sex to permit more specific search for hormonal determinants of RA, other than those already known demographic risk factors [9, 25]. Case or control subjects who had known cancer diagnoses during follow-up were excluded from the RAPS database. Their sera were reserved to study cancer biomarkers, which is the primary purpose of Operation CLUE [1921].

2.2. The Hormonal Reference Laboratory Performed Assays, in 1992 and 1994 (Table 1)

The mean concentrations of the original reported assays performed by the referral laboratory in the first 1992 and second 1994 sets of subjects are indicated in Table 1, as well as other statistical values. The baseline-stored (−70°C) cohort sera were always analyzed in matched sets of 1 pre-RA and 4 CN, without knowledge of subject status. As funding was secured during the interval of 1992 to 1994, the first (1992) and second (1994) sets of study subjects’ sera were sequentially donated by project CLUE for the hormonal assays in the referral laboratory. The first set of study subjects were 14 baseline pre-RA and 56 CN females. Their frozen sera were sent by CLUE to the Immunoassay Core Facility at Northwestern University (NWU) for hormonal assays in 1992 [16]. One subject had no sera available, but the other 69 had mostly sufficient sera to assay the full panel of hormones; a limited number lacked DHEA ( ), Δ4A ( ), and testosterone ( ) assays (Table 1). Serum estrone levels were only assayed in the first set females. The second set of female cohort subjects were 22 pre-RA and 88 CN subjects, whose frozen sera were sent from CLUE I to the NWU laboratory, in 1994. Those sera were sufficient to perform 4 hormonal tests completely (corticosterone, deoxycortisol, cortisol, and DHEAS) or almost completely (progesterone) (Table 2). Also, they were shared with other laboratories which performed immunological assays [16]. Two-thirds of the other tests were completed, but one-half of estradiol assays, which had the lowest concentrations, requiring greater amount of sera.

2.3. Specific Assays Were Developed for a Comprehensive Panel of Serum Steroids (Figure 1)

A comprehensive panel of adrenal and sex steroids were assayed (Figure 1) using a procedure developed specifically for this study to permit measurements of 12 steroids in duplicate from 1 mL of serum. The fractionation procedure permitted assay of each steroid in greater sensitivity than would otherwise be detectable without fractionation. The method also provided an additional purification step, adding to the specificity of the procedure.

Serum (1.0 mL) was diluted by addition of 1.0 mL of 4 M urea. 3H-tracers (5,000 cpm) of representative steroids, estrone (E1), testosterone (T), and progesterone (prog), were added and the solution was heated for 30 min at 60°C to denature the protein. Steroids were extracted by the use of an ODS cartridge (SepPak, Millipore, Inc., Billerica, MA). The cartridge was washed with 3 mL methanol followed by 6 mL of distilled water. The serum was aspirated through the cartridge. Retention of the steroids was complete, as negligible loss of the 3H-tracers into the effluent had occurred. A 7 mL of water wash was discarded, and the conjugated steroids (androstanediol glucuronide and DHEAS) were eluted with 10 mL of 47% aqueous methanol (Fraction A1). The second fraction was eluted with 12 mL of 60% aqueous methanol (Fraction B1). This fraction contained cortisol (F), deoxycortisol (DF), corticosterone (cmpd B), deoxycorticosterone (DOC), estrone (E1), and estradiol (E2). The third fraction was eluted with 12 mL of methanol: water: acetonitrile (55 : 35 : 10). It contained androstenedione (Δ4A), dehydroepiandrosterone (DHEA), 17-hydroxyprogesterone (17-OHprog), 17-hydroxypregnenolone (17-OHpreg), and testosterone (T) (Fraction C1). The fourth fraction was eluted with 12 mL of methanol and it contained progesterone (prog) and pregnenolone (preg) (Fraction D1).

The aqueous component of each of the first two fractions was removed by adding an equal volume of water and transferring to a second ODS cartridge. The flow-through was discarded, and the steroids were then eluted with 5 mL of methanol. The solvents from all fractions were then evaporated under a stream of nitrogen in a water bath at 60°C. The conjugated steroids in Fraction A1 were hydrolyzed by incubation with 10–12 units of β-glucuronidase in 1.0 mL of 0.1 M potassium phosphate buffer at 37°C for 18 hr. The now unconjugated androstanediol (adiol) and dehydroepiandrosterone (DHEA) were then extracted from the incubation mixture with ethyl ether (Fraction A2). The ethyl ether was then evaporated under a stream of nitrogen in a water bath at 60°C.

The estrogens were separated from the neutral steroids in Fraction B1 by solvent partition between 0.4 M aqueous NaOH and toluene. The aqueous fraction was neutralized by the addition of an equivalent amount of HCl, and the estrogens were extracted from the neutralized aqueous solution with ethyl ether (Fraction B2). The toluene was evaporated from the neutral steroids and the ethyl ether was evaporated from the estrogen fraction (Fraction B3). Each of the steroids was later measured by radioimmunoassays with available antibodies and 3H-steroids from New England Nuclear Corp, Newton, MA.

The biologically active fraction of total serum cortisol is free cortisol, but that assay is technically demanding, expensive, and not in general use. Total cortisol is bound to plasma proteins, particularly corticosteroid-binding-globulin (CBG) or sex hormone binding globulin (SHBG), hepatic proteins having estrogen-induced increased synthesis. Intra-assay percentile coefficients of variation (% CV) were all less than 12%, as the criterion for acceptability of measurement results. Too few batches of assays were performed in either of the 1992 or 1994 sets to analyze their interassay variability. The fractionation procedure permitted steroid differentiation as did the antibody specificity incorporated in the assays.

2.4. Statistical Methods

Frequency distributions of the hormonal values were examined for acceptability of unimodality and symmetry features. Extreme outliers were observed in several hormones, particularly progesterone, as expected by physiological peaks during the luteal phase, and the opposite near-zero values during the postmenopausal status. Several extreme high outliers were also found in estradiol, which were attributed to ovulatory surges. Extreme outliers were assigned to the upper ranges observed in the population frequency distribution curves, thereby diminishing their statistical influence [26]. Differences in mean values of assayed steroids between the first (1992) and second (1994) sets were normalized, usually adjusting from the smaller first set mean to the larger second set mean, but reversed for Δ4A, corticosterone, deoxycortisol, and DHEAS, in order to eliminate negative normalized values. Natural log conversion was performed on all values to improve their symmetry. Age-adjusted bivariate correlations (Pearson) of the individual steroids were performed on the normalized and log-converted values. These values were further transformed into -scores to reduce their variances (1) in performing -test differences of hormonal levels between subject groups and data subsets, and (2) in scatterplot examinations and confirmations of the data. The -score values of the total normalized subjects centered on 0.00 standard deviations (SDs) and were almost always distributed between 2 or 3 SDs. The Fisher -to- transformation was used to estimate significance of differences between two hormonal correlation coefficients in comparisons of total pre-RA versus CN as well as for subjects with lower than grand mean versus higher than grand mean baseline cortisol levels. When a significant difference was found in bivariate correlations between study groups and subgroups, a further correlational analysis was performed using set-specific -scores of the respective reported assays in 1992 and 1994. When both analytical methods revealed significant differences in correlations, those findings were indicated in the tables. A tentative illustrative model of adrenal glucocorticoid and androgenic anabolic (AA) steroid interrelations was inferred from hormonal differences observed between the total, lower, and higher cortisol subsets. In this exploratory study, a significance level of was accepted without adjustment for multiple comparisons [27].

3. Results and Discussion

3.1. Differences in Reported and Normalized Steroids in Pre-RA and CN (Tables 1 and 2)

The mean values of the normalized hormone assays (combining 1992 and 1994) for the total pre-RA and CN subjects were similar, except for androstenedione (Δ4A), which was lower ( ) in the 28 cases ( ) than the 108 CN ( ) (Table 1). The first set (1992) pre-RA females contributed mainly to that difference; 5 of those 12 earlier tested cases had Δ4A levels of −1 SD or lower than the total subjects’ mean level, compared to 1 of the 16 pre-RA tested in the second set ( ). The lower Δ4A levels were independent of baseline age, RA onset age, or interval in years from cohort entry to RA onset.

Our recent report [15] indicated that 4 (11.1%) of the 36 pre-RA females versus 1 (0.70%) of the 143 CN had definitely low baseline cortisol levels (<55 nmol/L, ). Accordingly, further analysis was performed on the study subjects stratified into subgroups of lower (<0.00 SD grand mean -scores) versus higher (0.00 SD or greater) mean cortisol levels of the total 180 females, both subsets having similar mean ages (Table 2). Again, the mean Δ4A level was lower in the pre-RA versus CN, but significant ( ) only in the subset of 9 pre-RA versus 53 CN with lower mean cortisol levels ( versus  nmol/L, resp.). The mean DHEA level was not significantly ( ) lower in 10 pre-RA versus 52 CN having lower mean cortisol levels ( versus , resp.). However, all cases had a negative -score value of <0.00 SD versus 27 (51.9%) of 52 CN ( ). This findings supports a previous report [8] of a greater ( ) proportion of 15 premenopausal RA patients having combined lower DHEAS and cortisol levels than 14 matched control women (40% versus none). The total 12 pre-RA subjects with lower cortisol values in this study had a borderline ( ) lower baseline mean cortisol level than the comparator 72 CN ( versus  nmol/L, resp.).

The mean level of corticosterone (cmpd B), the complementary steroid in the mineralocorticoid pathway to cortisol (Figure 1), was significantly greater in subjects having higher versus lower grand mean cortisol levels, among the CN ( ), pre-RA ( ), and total ( ) females (Table 2). Further, the bivariate correlations of corticosterone and cortisol levels were significant in the total (Table 3) and higher-level cortisol (Table 5) pre-RA and CN subjects, as well as in the lower-level cortisol pre-RA cases (Table 4), implying that these hormone levels may reflect corresponding hypothalamic-pituitary (H-P) activations. Without having the negative feedback control, as exists for cortisol, corticosterone production does respond to interval pulsatile activation of ACTH, and may be an approximate surrogate indicator of variations in H-P stimulation of cortisol levels. The subject group with higher cortisol levels likely represents greater than average adrenal cortical ACTH activation, which is inferred from the conjointly higher levels of cortisol and its non-17-hydroxylated comparable mineralocorticoid (Figure 1).

The mean estradiol (E2) concentration was greater ( ) in the 18 pre-RA subjects who had a higher than grand mean cortisol level than their 50 comparator CN ( versus  pmol/L, resp.). As a complementary estrogenic finding, the 9 pre-RA who had higher cortisol levels had greater ( ) mean estrone (E1) levels than the 5 cases with lower cortisol (Table 2). Among the CN subjects, however, the mean E2 level was greater ( ) in those who had lower ( ) versus higher ( ) baseline mean cortisol levels ( versus  pmol/L, resp.).

3.2. Bivariate Correlations of Hormones in Total, Lower, and Higher Cortisol Pre-RA versus CN

The total female pre-RA (top) versus CN (bottom) age-adjusted, mean-normalized, and log-transformed bivariate correlations (Pearson) of the hormone profile are indicated in Table 3. In the CN, as expected in normal physiology (Figure 1), pregnenolone, the initial steroid in the biosynthetic pathway, was strongly ( ) correlated with its proximate products: progesterone, 17-OH pregnenolone, and 17-OH progesterone. However, pregnenolone was not correlated with the adrenal androgens (DHEA and Δ4A), nor with cortisol. Rather, it was correlated with corticosterone ( ) and with deoxycortisol ( ). Of interest, pregnenolone was significantly ( ) correlated with estradiol (E2). In the smaller pre-RA sample, pregnenolone was also correlated significantly with progesterone, 17-OH progesterone, and deoxycortisol, but again, not with cortisol.

Regarding correlations of the adrenal androgens, DHEA, and Δ4A, each was strongly correlated with the other and with their 17-hydroxylated precursors, 17-OH pregnenolone and 17-OH progesterone, in both the pre-RA and CN subjects (Table 3). These findings indicate that the 17α-hydroxylase step is important in AA synthesis (Figure 1). The AAs were not correlated with the nonhydroxylated pregnenolone or progesterone enzymatic step precursors.

3.3. Similar Bivariate Correlations in the Total Pre-RA versus CN Subjects (Table 3)

The bivariate hormonal correlations for total pre-RA (top) and CN (bottom) subjects were similar (Table 3). Cortisol and E2 were reported to have counter-opposing effects on their alternate hypothalamic-pituitary (H-P) axis controls [17]. Cortisol and HPA axis stimulation tend to inhibit the HPG axis, whereas estrogen may stimulate the HPA axis and the peripheral production of cortisol [17].

3.4. Differences in Bivariate Correlations of Lower Cortisol Pre-RA versus CN Subjects (Table 4)

In the lower cortisol subject groups (Table 4), expected physiological interrelations were again observed, as described above for the total subjects (Table 3). The bivariate correlations were again generally similar between subject groups, but with 2 differences (deltas). The pre-RA had significantly stronger correlations than the CN in deoxycortisol with 17-OH progesterone ( , , versus , , , resp., delta ) and with Δ4A ( , , versus , , , resp., delta ). Again, these study group correlational differences raise the issue of the degree of ACTH activation of the adrenal cortical pathways, feasibly being stronger in the lower cortisol pre-RA versus CN subjects (Figure 2).

3.5. Differences in Bivariate Correlations of Higher Cortisol Pre-RA versus CN Subjects (Table 5)

Concerning the higher cortisol levels (Table 5), almost twice the number of pre-RA subjects were included than in the lower cortisol group (Table 4). Four low-level significant differences were observed, including 3 pairs of hormonal correlations involving pregnenolone, 17-hydroxypregnenolone, and estrone (E1). In pre-RA, the pregnenolone correlation was negative with cortisol in pre-RA ( , , ), but positive in the CN ( , , ), the difference being significant ( ). Regarding the hydroxylated precursor steroid, 17-OH pregnenolone, the pre-RA had even stronger positive correlation with DHEA ( , , ) than did the CN ( , , ), the difference also being significant ( ) (Table 5). The stronger positive AA steroid correlation of DHEA with its precursor steroid, 17-OH pregnenolone, may imply that the observed relative deficiency of this AA had not likely resulted from a specific inhibition at the 17α-hydroxylase step (Figure 1).

The correlations of estrone (E1) with both corticosterone and deoxycortisol were negative in pre-RA ( and , resp.), but positive in CN ( and , resp.), the differences being and , respectively (Table 5). The mean estrone level was greater ( ) in the pre-RA women with higher versus lower cortisol levels ( versus  pmol/L, resp.), as indicated in Table 2.

3.6. Inferences from the Pre-RA versus CN Hormonal Profile Mean Differences and Correlations

In the absence of quantitative measurements of the degree of ACTH, FSH, and LH end-organ stimulation (Figure 2), interpretation of the baseline serum adrenal and sex steroid profiles is challenging in the pre-RA versus CN subjects of this study as well as in previous comparisons of RA patients and controls [48, 12, 2832]. Adrenal cortical stimulation by the insulin-induced hypoglycemia test (IIHT) was studied in pre-menopausal RA patients and control females. One study [33] did not reveal a significant difference in cortisol response. The other study [8] revealed a baseline combined “lower” quartile range of serum cortisol and DHEAS levels including 6 of the 15 early RA patients, but none of 14 CN ( ). The finding suggested a relative hypocompetence of adrenal cortical function in some premenopausal RA females [8]. Few studies have measured ACTH, serum cortisol, or other steroid levels of RA patients without having glucocorticoid therapy at baseline measurements or after HPA axis stimulation [4, 12, 3439]. Two studies [34, 36] showed significantly elevated ACTH levels without hypercortisolemia in untreated RA patients, but a difference was not found in the other studies. Lower DHEA levels in two of the above studies suggested decreased synthesis of AAs and deficient zona reticularis function in the RA versus CN [4, 12].

To our knowledge, this is the first study of a comprehensive panel of adrenal and sex steroid levels in women prior to the onset of RA and matched cohort controls. The lower androstenedione level in the pre-RA occurred mainly in those subjects who also had lower cortisol levels than the grand mean of total subjects. In that lower-cortisol subgroup (Table 4), pre-RA had a stronger correlation of deoxycortisol with Δ4A than did the CN. The preceding non-AA steroid level may be a surrogate indicator of ACTH stimulation, which is not directly regulated by the negative feedback mechanism controlling cortisol levels (Figure 2).

The lower AA levels in a minority of the pre-RA cases may likely reflect a lower AA synthesis capacity in those women, rather than lesser ACTH stimulation. Perhaps, such physiology may be analogous to natural processes in aging [18, 40, 41]? Figure 2 outlines an inductive model in which a minority subset of pre-RA women may have lower or relatively insufficient adrenal cortical functional capacity, rather than a deficit of hypothalamic-pituitary (H-P) stimulation or specific enzymatic defects [4145]. The mechanisms for biological impairments of AA production are profoundly complex and are not yet fully defined but do occur in natural aging [18, 40, 41]. The biosynthetic mechanisms controlling AA synthesis and their levels during life cycles and between the genders remain unproved [4145]. Further research is required to confirm the observed relative AA insufficiency in a minority subset of women before onset of RA as well as the differences found in hormonal correlations. If these findings are confirmed, they will need to be defined in terms of their biosynthetic and control mechanisms [18, 4045].

4. Conclusions

Serum androstenedione (Δ4A) levels were lower in a cohort of pre-RA versus CN females, particularly among subjects who had lower than the total subjects’ mean cortisol level (Tables 1 and 2). In pre-RA versus CN women having lower cortisol levels (Table 4), stronger correlations of deoxycortisol with 17-OH progesterone and with Δ4A were found in the cases, suggesting that those subjects had sufficient or potentially increased ACTH stimulation. The pre-RA with higher cortisol levels showed a stronger correlation of DHEA with its 17-OH pregnenolone precursor than did the CN (Table 5). Accordingly, the AA deficiency does not seem likely based upon the respective enzymatic defect.

Multiple comparisons were performed in this study which could lead to an inflated rate of type I error, identifying low-level differences which do not truly exist, and which require independent confirmation. Steroid interactions are part of complex homeostatic control mechanisms which tend to modulate deviations. In such systems, accumulated minor deviations in the same direction could be supportive of subtle dysfunctions, as is proposed. However, the observed differences were found only in a minority subgroup of pre-RA women versus CN women. Further investigation is needed of possible polymorphic variations diminishing cellular biosynthetic capacities of adrenal and ovarian steroids as a risk factor for RA in a subset of females (Figure 2).

Conflict of Interests

The authors declare that they have no conflict of interests.

Acknowledgments

The CLUE study, Campaign Against Cancer and Stroke, was funded in 1974 by Grant CA 11849 from the National Cancer Institute. The authors thank Mary Jean Dzursin for artistic expertise in helping to finalize the figure on categories of hormonal variation. The Research Open Access Publishing (ROAAP) Fund of the University of Illinois at Chicago provided the financial support towards the open access publishing fee for this paper. Support for this project was provided by the Department of Medicine, University of Illinois College of Medicine at Peoria, and by a gift from The MTM Foundation, with no existing or potential conflict of interests.