- Research article
- Open Access
CYB5A polymorphism increases androgens and reduces risk of rheumatoid arthritis in women
Arthritis Research & Therapy volume 17, Article number: 56 (2015)
Rheumatoid arthritis (RA) is characterized by decreased androgen levels, which was the first hormonal abnormality described. Several studies indicated that steroidogenesis is directed towards endogenous glucocorticoids at the expense of androgens. The decisive step governing androgen synthesis is the 17,20-lyase activity of the CYP17A1 gene-encoded enzyme cytochrome P450 17A1. Here, we focused on the role in RA of the critical cofactor for 17,20-lyase activity, cytochrome b5, encoded by the CYB5A gene.
Data sets of two genome wide RA association studies (GWAS) were screened for single nucleotide polymorphisms (SNP) in the CYB5A gene. Candidate SNPs in CYB5A were studied in a case–control study population of Slovakia. Expression analyses were done in synovial fibroblasts from RA patients by quantitative real-time polymerase chain reaction, and cytochrome b5–expression was detected by immunohistochemistry. Real-life androgen production after steroid conversion was measured using radiolabeled substrates.
The study identified the RA-associated intronic SNP rs1790834 in the CYB5A gene in one GWAS and confirmed the same SNP in our study. The minor allele reduced RA risk selectively in women (P = 4.1*10−3; OR = 0.63, 95% CI [0.46-0.86]). The protective effect was confined to rheumatoid factor-positive (OR = 0.53, [0.37-0.75]) and anti-cyclic citrullinated peptide-positive (OR = 0.58, [0.41-0.83]) cases, respectively. The protective allele doubles CYB5A mRNA-expression resulting in 2-3fold activation of steroid 17,20-lyase activity, and protective allele was accompanied by a higher density of cytochrome b5-positive cells in synovial tissue.
CYB5A is the first RA susceptibility gene involved in androgen synthesis. Our functional analysis of SNP rs1790834 indicates that it contributes to the sex bias observed in RA.
Rheumatoid arthritis (RA) is a chronic inflammatory joint disease that affects about 0.5 to 1.0% of the population. It causes increasing disability leading to a huge socioeconomic burden . In RA, risk variants in 46 loci explain about half of the genetic risk, indicating that other so far unknown loci are involved [2,3]. Despite the autoimmune etiology, neuroendocrine immune pathways relevant for inflammatory processes are discussed in RA onset and progression . In RA there is a clear preponderance of affected women over men (3:1) . This suggests that high concentrations of estrogens, low concentrations of androgens, or a combination of both increase the risk for RA .
Several androgens have anti-inflammatory properties. Dehydroepiandrosterone (DHEA), androstenedione, and testosterone inhibit secretion of IL-1β, IL-6, TNF, and others [7-12]. The androgen 5α-dihydrotestosterone inhibits activation of the human IL-6 gene promoter stimulated by nuclear factor kappa B , and it decreases T cell proliferation . There is evidence that some RA patients of both sexes have reduced amounts of serum androgens, even years before disease onset [15,16]. Particularly, female RA patients have lower than normal levels of DHEA and/or DHEA sulfate. In male RA patients, levels of serum testosterone are negatively correlated with disease severity . Two double-blind interventional studies with testosterone demonstrated some benefit in patients with RA [17,18]. A preponderance of serum glucocorticoids over serum androgens is known in many chronic inflammatory diseases, which normalizes after anti-TNF therapy . In summary, there is convincing evidence that a relative lack of androgens is involved in the etiology of RA.
However, we do not know molecular mechanisms responsible for this state of androgen deficiency. A recent study revealed a negative correlation of serum testosterone levels with RA disease activity in male patients under therapy with disease-modifying anti-inflammatory drugs . This indicates that the inflammatory disease can reduce androgen levels. This may depend on increased androgen-to-estrogen conversion that can happen in inflammatory cells such as macrophages and osteoblasts [21-23]. Increased estrogen formation and estrogen-to-androgen ratios were found in the synovial fluid of RA patients compared to controls showing increased aromatase activity . We demonstrated that mixed synoviocytes from patients with RA and osteoarthritis convert DHEA, androstenedione, and testosterone into downstream hormones . If androgens like testosterone are missing in the tissue, synovial aromatase activity is strongly stimulated, and this leads to a very high synovial estrogen-to-testosterone ratio [6,25].
In addition, in collagen type-II arthritic animals and in synovial fibroblasts from RA patients, conversion of DHEA into the proinflammatory metabolite 7α hydroxy-DHEA is increased (gene CYP7B1) [26,27]. This means that DHEA is not available for conversion to androgens such as androstenedione and testosterone.
However, these findings in patients with overt inflammatory disease do not explain why androgen levels are lower than normal before disease outbreak [15,16]. Some authors speculated that male RA patients without glucocorticoid treatment might be in a state of (compensated) partial gonadal failure . However, molecular pathways of adrenal or gonadal failure are not known, and this was subject of our study.
Endogenous de novo synthesis of androgens depends on two key enzyme activities, 17α-hydroxylase and 17,20-lyase, both linked to one protein encoded by the cytochrome P450 17A1 gene CYP17A1 . While 17α-hydroxylase activity is essential for synthesis of androgens and cortisol depending on the presence of NADPH-cytochrome P450 reductase (POR) , the 17,20-lyase activity depends on the combined presence of POR and the cofactor cytochrome b5 type A (gene on chromosome 18, CYB5A) [29,30].
In order to study molecular pathways relevant to androgen deficiency, we focused on cytochrome P450 17A1 and cytochrome b5 type A. In the two published RA genome-wide association study (GWAS) datasets, we reanalyzed single nucleotide polymorphisms (SNP) in the CYB5A gene [31,32]. In a novel RA case-control study with Slovak people, we investigated the same SNPs in the CYB5A gene. Then, we functionally analyzed SNP alleles by steroidogenic gene expression and pregnenolone conversion into androgens in synovial fibroblasts of RA patients. In addition, density of cytochrome b5A-positive synovial cells was investigated in RA synovial tissue.
A total of 842 (117 male, 665 female) Slovak individuals were included for genotyping, with 521 (87 male, 434 female) RA patients and 321 (90 male, 231 female) healthy controls without any arthritic symptoms. Rheumatoid factor (RF) was determined by standard techniques. Measurement of antibody against cyclic citrullinated peptide (anti-CCP) was carried out using an anti-CCP ELISA (Euroimmun, Lübeck, Germany) following the manufacturer’s instructions. From a total of 304 RA patients anti-CCP levels were determined. Values <4.2 RU/ml were anti-CCP negative. The characteristics of these patients are given in Table 1.
Synovial fibroblasts were isolated from 40 patients with longstanding RA fulfilling the American College of Rheumatology (formerly, the American Rheumatism Association) revised criteria for RA . These patients underwent elective knee joint replacement surgery. We obtained written consent of patients according to the current Declaration of Helsinki. The study was approved by the Ethics Committee of the National Institute of Rheumatic Diseases, Piestany, Slovakia, and the Ethics Committee of the University of Regensburg, Germany.
Marker selection and genetic analyses in case-control study
SNPs in CYB5A genomic region (NM_148923.3) with 10 kb on the 5′ and 3′ end, respectively (chromosome 18: 71,910,527-71,969,251; NCBI build 37.3), were selected based on published information from two GWAS on RA with P <0.01 (see Table 2, and previous publications [31,32]).
Genomic DNA was isolated from whole blood samples using the PureGene DNA Blood Kit (QIAGEN, Hilden, Germany). DNA samples were genotyped using 5′ exonuclease TaqMan® technology (Life Technologies, Applied Biosystems, Foster City, CA, USA) following the manufacturer’s instructions. In brief, for each genotyping experiment 10 ng DNA was used in a total volume of 5 μl containing 1 × TaqMan® Genotyping Master Mix (Applied Biosystems). Real-time polymerase chain reaction (RT-PCR) and post-RT-PCR endpoint plate read was carried out according to the manufacturer’s instructions using the Applied Biosystems 7900HT RT-PCR System. Sequence Detection System software version 2.4.1 (Applied Biosystems) was used to assign genotypes applying the allelic discrimination test. Case and control DNA were genotyped together on the same plates with duplicates of samples (10%) to assess intraplate and interplate genotype quality. No genotyping discrepancies were detected. Assignment of genotypes was performed by a person without knowledge of sample case-control status.
Preparation of synovial tissue and culture of synovial fibroblasts
Synovial tissue samples from patients with RA were obtained immediately after opening the knee joint capsule as previously described . Pieces of synovial tissue of up to 9 cm2 were excised. One part of the tissue sample was cut, paraffin-embedded, and stored until further use (immunohistochemistry, see below). Another part of the synovial tissue specimen was minced and put in Dispase I (Roche Diagnostics, Mannheim, Germany). Digestion lasted for at least 1 hour at 37°C on a shaking platform. The resulting suspension was filtered (70 μm filter) and spun at 300 g for 10 minutes. We treated the pellet with erythrocyte lysis buffer (20.7 gm NH4Cl, 1.97 gm NH4HCO3, 0.09 gm EDTA, and 1 liter of H2O) for 5 minutes followed by centrifugation for 10 minutes at 300 g. The pellet was resuspended in RPMI 1640 (Sigma-Aldrich, St Louis, MO, USA) with 10% fetal calf serum. We transferred 1,000,000 cells to a 75-cm2 tissue culture flask, incubated them overnight, and supplemented these cells with fresh medium.
Synovial fibroblasts derived from these cultures were routinely used in passages two or three, where these cells are known to conserve their phenotype . Cells were seeded into 6-well plates for steroid conversion experiments or 25-cm2 tissue culture flasks for nucleic acid isolation and were kept in a humidified atmosphere with 5% CO2 at a temperature of 37°C.
Nucleic acid extraction and quantitative RT-PCR analysis
Synovial fibroblasts were disrupted using QIAshredder colums (QIAGEN). Total cellular RNA and DNA were isolated in parallel from cultured cells using the AllPrep DNA/RNA Mini Kit (QIAGEN). Genotyping was done as described above and cDNAs were generated from RNA by reverse transcriptase reaction using AffinityScript quantitative RT-PCR synthesis kit (Agilent Technologies, Böblingen, Germany). RT-PCR was performed using TaqMan primer/probes (Applied Biosystems - Life Technologies, Darmstadt, Germany) targeting two isoforms of CYB5A gene (Hs00157217_m1 and Hs01076969_m1). In addition, expression of several genes coding for enzymes relevant for steroid metabolism was assessed in synovial fibroblasts using TaqMan assays: 17α-hydroxylase/17,20-lyase (Hs01124136_m1 for CYP17A1), NADPH-cytochrome P450 reductase (Hs01016332_m1 for POR), 3β-hydroxysteroid dehydrogenase isoforms (Hs00426435_m1 for HSD3B1, Hs00605123_m1 for HSD3B2, and Hs00228639_m1 for HSD3B7, respectively), the 25-hydroxycholesterol 7α-hydroxylase cytochrome P450 7B1 (Hs00191385_m1for CYP7B1), and aromatase cytochrome P450 19A1 (Hs00903413_m1 for CYP19A1). Relative gene expression was normalized to HPRT1 mRNA (Hs01003267_m1) levels using the comparative cycle threshold (Ct) method and presented as expression ratio using 2-ΔΔCt . RA cDNA with ΔCt value next to the median of the ΔCt distribution was used as calibrator sample.
Human cytochrome b5 type A was detected and quantified using immunohistochemistry. Five-micrometer sections were cut from paraffin-embedded blocks. For immunostaining, the sections were deparaffinized in xylol and ethanol and rehydrated in a descending ethanol series. The slides were placed in a 0.1 mol/L citrate buffer (pH = 6.0) and boiled for 20 minutes at 120°C. Slides were then washed in PBS buffer (pH = 7.6), and treated briefly with 2% hydrogen peroxide in 1 × PBS for 10 minutes to block endogenous peroxidase activity. Slides were incubated for 45 minutes in blocking solution consisting of 10% bovine serum albumin (PAA, Inc., Pasching, Austria), 10% fetal calf serum (PAA, Inc.), and 10% goat serum (Sigma, Steinheim, Germany). After blocking, sections were incubated with primary antibody (1:500, ab69801, Abcam, Cambridge, MA) overnight at room temperature. After incubation with biotinylated secondary antibody (1:200, polyclonal goat anti-rabbit IgG biotinylated, from Dako, Hamburg, Germany) for 90 minutes, the sections were treated with streptavidin-horseradish peroxidase (GE Healthcare Europe GmbH, Munich, Germany). Visualization was performed with diaminobenzidine tetrahydrochloride hydrate (ImmPACT DAB substrate, Vector Laboratories, Peterborough, UK). Staining with control serum always yielded negative results, and liver tissue was used as positive control. In order to determine the tissue density of cytochrome b5-positive cells, the number of positive cells was averaged from 17 randomly selected high-power fields (400×). The number of investigated high-power fields was derived from a pioneering histological study .
Incubation with radiolabeled steroids and steroid extraction
For steroid conversion experiments synovial fibroblasts in 6-well plates were washed with PBS and incubated with RPMI 1640 (without phenol red or serum, but supplemented with 15 mM HEPES, 2 mM stable glutamine and 40 mg/l gentamycin) (PAN Biotech, Aidenbach, Germany). After 3 hours, the radiolabeled substrate 7-3H(N)-pregnenolone (Preg, 5-pregnen-3β-ol-20-one) (PerkinElmer, Rodgau-Jügesheim, Germany) was added for another 24 hours at a final concentration of 25 nM containing an activity of 37,000 Bq.
The fibroblast monolayers were fixed overnight with 10% formaldehyde, washed with saline, and nuclei were stained with Hoechst 33258. Large-area images were taken on a Cell Observer microscopic system (Carl Zeiss, Jena, Germany) and counted nuclei were taken as cell numbers for normalization of steroid conversion results.
From the cell-free supernatants steroids were extracted twice with 3 ml cold ethyl acetate. The exact concentration of radiolabeled steroid applied to every well and the extraction efficiencies were monitored by liquid scintillation counting of aliquots in Ultima Gold cocktail (PerkinElmer): recoveries of radioactivity in the organic phase varied insignificantly (97.2 ± 1.5%, mean ± SD). The stable extracts were lyophilized in a speed-vac concentrator (Saur, Reutlingen, Germany) and stored at −20°C until analysis.
Two-dimensional thin-layer chromatography (2D-TLC) of steroids
Pregnenolone and its metabolites with respect to androgen synthesis  (Figure 1A) were separated essentially as described previously , with modifications as given below. Solvents and other reagents were purchased from Merck (Darmstadt, Germany), if not stated otherwise. Unlabeled steroids were from Sigma (Steinheim, Germany) and from Steraloids (Newport, RI, USA). Stock solutions were prepared in ethanol. Lyophilized steroid extracts were dissolved in 50 μl ethanol, spotted under a nitrogen stream with a Linomat IV (Camag, Muttenz, Switzerland) on silica gel 60 F254 TLC aluminum sheets (Merck) together with a mixture of unlabeled carrier steroids (Figure 1B). The first separation of 2D-TLC was done in toluene:methanol (90:10). The second dimension of the separation was done by two successive developments in chloroform:diethylether (90:10). For identification of spots, the 2D-TLC plates were stained with copper acetate in phosphoric acid as described previously . Radioactivity on the 2D-TLC plates was quantified by radioimaging on an FLA 3000 (Fuji-Raytest, Straubenhardt, Germany). Spots were assigned only if their intensity was more than two standard deviations above background. Spot intensities were corrected for corresponding values obtained from blank wells incubated with radiolabeled pregnenolone without cells. The results were calculated as pmol of steroid produced per million cells in a 24-hour incubation period.
Statistical and bioinformatical analyses
Differences between dichotomous traits were calculated employing the chi-square (χ 2) test. Differences in continuous variables between groups were calculated using the two-tailed t-test for normally distributed data or the non-parametric Mann-Whitney U-test for variables not normally distributed as determined by the Shapiro-Wilk test. Expression data were compared between the three genotypes after grouping the minor allele carriers (dominant model). Logistic regression was used for the different models in association analyses (sex as a covariate, stratified by sex, and stratified by either RF or anti-CCP status). Odds ratios (OR) and 95% CI were calculated. To determine whether the genotypes of cases and controls of all SNPs deviated from Hardy-Weinberg equilibrium, actual and predicted genotype counts of both groups were compared by the exact test . A P-value <0.05 was considered significant.
Association analyses were performed with PLINK v1.07 . Statistical software packages JMP 7.0.2 (SAS Institute Inc, Cary, NC, USA) and SigmaPlot 11.2.05 Systat Software Inc, San Jose, CA, USA) were used for other analyses. Meta-analysis weighted according to sample size was performed using METAL with Cochran’s Q-test for heterogeneity [40,41]. For linkage disequilibrium testing, HaploView v4.2 was employed  with HapMap version 3 data from release 28 . A power calculation was performed with Quanto 1.2.4 (Division of Biostatistrics, USC University of Southern California, Los Angeles, CA, USA) assuming a population risk of 1% for RA. MatInspector v8.05 and SNPInspector 2.2 (Genomatix Software GmbH, Munich, Germany) were used to analyze the potential functional consequences of SNPs on transcription factor binding sites.
Gene expression and steroid conversion data are shown as mean ± SD, if all data in a figure were normally distributed. Otherwise, data are presented as box plots, where boxes represent the 25th and 75th percentiles, respectively. The median is indicated by a horizontal line. The whiskers indicate the 10th and 90th percentiles, respectively. Dots represent individual outliers.
Genotyping and association study
Two recently published GWAS on RA susceptibility [31,32] were screened for association between disease and CYB5A polymorphisms, resulting in two SNPs with P <0.01. In our Slovak RA case-control sample (Table 1), we sought direct replication of these SNPs (Table 2). In our sample, controls were significantly younger than cases at inclusion into the study but were in a comparable age range when the onset of disease in RA cases was compared to age at inclusion of controls. Significantly more females were present in the case group. Both markers fulfilled our criteria of at least a 95% call rate and no deviation from Hardy-Weinberg equilibrium was observed (Table 3). Allele frequencies for the two SNPs were consistent between HapMap release 28 data for the population of Utah residents with ancestry from northern and western Europe (called the CEU population) and our disease-free controls (1% and 18% in HapMap and 1.2% and 18.9% in controls for rs1790858 and rs1790834, respectively).
After adjustment for sex, rs1790834 showed association with risk of RA (P = 8.5*10−3; OR = 0.69, 95% CI 0.52, 0.91) (Table 3). Additional adjustment for age (age at inclusion for controls and age of onset for RA patients) did not change the association results (data not shown). Meta-analysis with the initial association data from a previous publication  resulted in P = 3.6*10−4 (OR = 0.76, 95% CI 0.69, 0.84), that is, a protective effect of the minor allele. No heterogeneity was detected (P = 0.237).
When the association of rs1790834 with risk of RA was analyzed separately for both sexes, there was significant association only for women with P = 4.1*10−3; OR = 0.63 (95% CI 0.46, 0.86; n = 434 cases, n = 231 controls). There was no association detectable for men (n = 87 cases, n = 90 controls) in our study sample (P = 0.81; OR = 0.93, 95% CI 0.52, 1.67).
In addition, stratified analysis for RF-positive (n = 279) and RF-negative cases (n = 228) was performed: rs1790834 showed only significant association results in RF-positive cases with P = 3.3*10−4 (OR = 0.53, 95% CI 0.37, 0.75) in contrast to RF-negative cases with P = 0.38 (OR = 0.86, 95% CI 0.62, 1.20). In addition, in the anti-CCP positive subset of RA patients, significant association between rs1790834 and the disease was detected with P = 2.7*10−3 (OR = 0.58, 95% CI 0.41, 0.83).
Gene expression analysis
Synovial fibroblast RNA samples from 22 RA patients with measured rs1790834 genotypes were available for gene expression analysis. SNP rs1790834 is located in the first intron of the CYB5A gene, that is, not within the protein coding sequence. Therefore, we first verified expression of CYB5A and several key genes of pregnenolone metabolism: fibroblasts expressed the longest mRNA-transcript of CYB5A (GenBank NM_148923; cytochrome b5 isoform 1), which encodes the membrane-bound fully functional cytochrome b5. In addition, CYP17A1 and POR as prerequisite for DHEA synthesis (and hence androgen synthesis), and HSD3B2, HSD3B7, CYP7B1 and CYP19A1 as prerequisite for downstream processing of DHEA were expressed in synovial fibroblasts (Figure 2A), whereas HSD3B1 was not found.
In a second step, we tested whether the minor rs1790834 A-allele has an effect on gene expression levels. Comparing CYB5A mRNA expression levels, the mean expression of CYB5A was significantly higher within carriers of the minor allele (Figure 2B). The rs1790834 genotype had no effect on the expression levels of the other steroidogenic genes expressed in synovial fibroblasts (Figure 2B), indicating that the SNP specifically affects CYB5A expression.
Biobank tissue sections were available for a subgroup of RA patients. In these tissue sections, we studied immunohistochemically investigated density of cytochrome b5-positive cells. Density of positive cells positively correlated with amount of mRNA (n = 8, r 2 = 0.64, P = 0.017; data not shown). Thus, the minor rs1790834 A-allele facilitates higher cellular cytochrome b5 contents.
CYB5A genotype affects pregnenolone metabolism in RA synovial fibroblasts
As fibroblasts metabolize various steroid substrates, we tested whether the CYB5A genotype affects the conversion of pregnenolone into downstream metabolites. This pathway was studied because the key enzyme for androgen synthesis, the 17,20-lyase activity of cytochrome P450 17A1, depends on cytochrome b5, the product of the CYB5A gene.
Synovial fibroblasts from 40 genotyped RA patients (n = 27 for major allele GG carriers, n = 13 for minor allele AA/AG carriers) converted radiolabeled pregnenolone into a variety of steroids (for products/pathways identified see Figure 1A). By 2D-TLC we identified nine products of pregnenolone metabolism (Figure 1B-D). There were no allele-dependent differences in the production of progesterone, 17α-hydroxyprogesterone, DHEA, androstenedione, androstenediol, testosterone, and an unidentified product X, which was the major product from fibroblasts, but independent of genotype (GG: 187 ± 152 versus AA/AG: 233 ± 175 pmol/(106 cells*24 hours), mean ± SD; t-test: P = 0.40) (Figure 3A, B).
Accumulation of the direct product of the 17α-hydroxylase activity, 17α-hydroxy-pregnenolone, in synovial fibroblasts from carriers of the minor A-allele was significantly lower than in non-carriers (GG-genotype) (Figure 3A). Carriers of the minor A-allele produced roughly 2.5 times the amount of 7α-hydroxy-DHEA than non-carriers (P = 0.009, Figure 3B). As production of 7α-hydroxy-DHEA depends on DHEA availability this metabolite indicates increased 17α-hydroxylase and 17,20-lyase activity.
Decreased ratios of 17α-hydroxylase over 17,20-lyase activity would be expected when the CYB5A A-allele stimulates the latter enzyme step via increased availability of cytochrome b5. Therefore, metabolite ratios were calculated: 17α-hydroxy-pregnenolone/7α-hydroxy-DHEA, 17α-hydroxy-pregnenolone/17,20-lyase (that is, the sum of all C19-steroids produced), and 17α-hydroxylase (that is, 17α-hydroxy-pregnenolone plus 17α-hydroxy-progesterone)/17,20-lyase, were significantly decreased (P = 0.006, P = 0.022, and P = 0.03, respectively) in synovial fibroblasts from RA patients with the rs1790834 A-allele (Figure 3C-E), clearly indicating an effect of the CYB5A genotype on steroid 17,20-lyase activity.
Finally, the two groups of the study were compared in order to check homogeneity with respect to clinical variables. The patients did not receive intra-articular glucocorticoid therapy. Importantly, age, sex, and prednisolone dosage was similar in the group with the minor allele (A) as compared to the group with the major allele (G) (age 63 versus 68 years; female/male 83%/17% versus 90%/10%; and dosage mean ± standard error of the mean 4.58 ± 1.00 versus 6.50 ± 1.50 mg/day). Similarly, the frequency of prednisolone and methotrexate treatment was similar in the two groups (prednisolone 83% versus 80%, and methotrexate 33% versus 10%). The patients did not receive anti-TNF therapy or other disease-modifying anti-inflammatory drugs because the disease was in a mild and late stage. In addition, markers of inflammation like C-reactive protein (23.4 ± 14.8 versus 16.8 ± 4.1 mg/l) and erythrocyte sedimentation rate (21 ± 5 versus 39 ± 17 mm/1st hour) were not significantly different between groups.
Owing to the design of GWAS, below their threshold of unambiguous detection important genetic variants remain hidden (in a so-called sea of noise). Nevertheless, these genetic variants may contribute significantly to the risk of RA. Instead of enforcing more statistical power by increasing numbers of cases and controls in GWAS, a hypothesis-driven approach combined with functional analyses can help to unravel significant contributions of candidate genes to risk of RA . RA is associated with comparatively low androgen levels, both systemically and within the inflamed tissue [20,24,25,45,46]. Therefore, we screened the data sets of two published GWAS for polymorphisms in the CYB5A gene [31,32], which is relevant for steroid metabolism and identified two RA-associated SNPs linked to low androgen secretion.
In our Slovak case-control study sample, the minor allele of SNP rs1790834 was associated with reduced risk for RA. The power for this replication of a GWAS finding  was about 42% in our study with 521 cases and 321 controls (using the reported OR = 0.83, minor allele frequency = 18% and one-sided P-value = 0.05). In subgroup analysis, we detected association between SNP rs1790834 in RF-positive RA cases only. This is consistent with the initial study by Plenge et al., including patients that were seropositive for anti-CCP . Anti-CCP is a better diagnostic marker than RF, and the correlation between the two factors is moderate . Although anti-CCP levels were only available in a subset of our RA patients, we again found a strong association between rs1790834 and the disease in this subgroup. Altogether, our data and the initial data by Plenge et al.  suggest that rs1790834-driven association with RA is limited to an RF/anti-CCP-positive disease entity. Heritability is in a similar range for both anti-CCP-positive and -negative RA, but these represent different disease entities .
Importantly, stratification for sex revealed an association between rs1790834 and RA only in women in our study, indicating a strong sex bias for the protective effect of the minor allele A. One possible explanation is the reduced power in our male sample with 87 cases and 90 controls (36% power given an OR of 0.69, an allele frequency of 18.9% in controls, a nominal significance with one-sided P-value = 0.05 in an additive model and an assumed population risk of 1%). In comparison, the power was 77% to detect this effect in women (434 cases and 231 controls). However, as RA is a disease with higher prevalence in women than in men, sexual hormones are thought to possibly be involved in disease onset and progression .
With a minor allele frequency of about 18% for rs1790834 this CYB5A polymorphism is not rare. Additional SNPs may even increase the contribution of CYB5A polymorphisms to reduction of RA risk. Our inability to replicate the effect of the other SNP, rs1790858 , does not contradict this notion because the protective allele frequency of 1% for this SNP is too low to be detected by our study design. In contrast, it suggests that in future work all independent SNPs in the CYB5A gene should be analyzed in detail for their possible contribution to RA risk.
Cytochrome b5, the protein encoded by the CYB5A gene, has diverse functions, which are of varying importance in different tissues and cell types. Notably, NADH-dependent reduction of methemoglobin (erythrocytes), electron transfer to fatty-acid desaturases (mainly liver), and steroid and xenobiotics metabolism by cytochrome P450 enzymes depend on cytochrome b5 [29,49,50]. Cytochrome P enzyme activity is increased by at least two different mechanisms: either cytochrome b5 mediates electron transport to cytochrome P450 enzymes or it facilitates allosteric interaction of selected cytochrome P450 enzymes with the electron donating NADPH-cytochrome P450 reductase [49,50]. The latter mechanism potentiates the activity of the affected cytochrome P450 enzymes, amongst them enzymes involved in metabolism of xenobiotics , and the key enzyme activity for androgen biosynthesis, cytochrome P450 17A1 steroid 17,20-lyase, our target enzyme . As androgens and androgen metabolism play a role in the etiology of RA, there was an obvious need to further analyze the SNP rs1790834 in this context.
To test for functional consequences of CYB5A polymorphisms, a patient-specific source of tissue or cells was indispensable. With respect to androgen synthesis, this would be adrenal glands or gonads [29,51], which are unavailable for these experimental studies in RA. On the other hand, synovial fibroblasts are a validated model for gene expression and functional studies in RA research. They maintain a well-conserved phenotype in cell culture over several passages [35,52], and patient-derived fibroblast lines largely maintain their gene expression profiles .
We consistently detected significant amounts of pregnenolone converting enzymes in synovial fibroblasts. As expected for a cell type that is not primarily involved in androgen production, the activities of those enzymes are low when compared to activities in adrenals or gonads [29,30,51]. Nevertheless, we found all enzymes necessary for and all metabolites indicative of androgen biosynthesis in synovial fibroblasts. Moreover, the critical step of androgen synthesis, the 17,20-lyase activity of CYP17A1-encoded cytochrome P450 17A1, clearly depends on the CYB5A alleles (determining the amount of cytochrome b5). High levels of cytochrome b5 can increase androgen synthesis more than 10-fold in experiments using purified enzymes from adrenals or gonads, respectively, or using recombinant enzymes [30,53]. We found a 2- to 3-fold increase of the capacity for androgen synthesis in synovial fibroblasts harboring the RA-protective allele A of SNP rs1790834, which parallels the increase in CYB5A expression. This effect size is reasonable, as (1) almost all tested fibroblast lines were from heterozygous donors (AG) and (2) the major allele homozygous carriers (GG) express CYB5A (at normal levels). Neither the expression levels of CYP17A1 nor that of any other gene analyzed were altered. Therefore, the intronic SNP, which does not alter the protein sequence of the CYB5A-encoded cytochrome b5, most likely exerts its effect specifically via control of the expression level.
It should be emphasized that the accumulation of (proinflammatory) 7α-hydroxy-DHEA is a specific feature of synovial fibroblasts, which express high levels of CYP7B1-encoded cytochrome P450 7B1 [26,27], but not a consequence of cytochrome b5 levels in general. For tissues or cell populations, where cytochrome P450 7B1 is not the dominant enzyme for DHEA metabolism, the CYB5A SNP would modulate the production of other androgens and their metabolites.
Although adrenals and gonads from RA patients or healthy controls are not amenable to direct testing, our results implicate that the identified CYB5A SNP may be generally responsible for facilitated androgen production in these tissues of minor allele carriers. Correspondingly, the major allele determines lower androgen levels , which precede disease onset in RA patients [15,16]. In addition, local androgen production in the inflamed joints would be hampered by the major allele after disease onset [20,24,25].
In summary, the identified CYB5A SNP should affect the 17,20-lyase step of androgen biosynthesis in all tissues expressing cytochrome P450 17A1. Due to the modulatory role of cytochrome b5 on 17,20-lyase activity, as outlined above, the alterations in pregnenolone metabolism resulting from the possible SNP variant combinations are expected to be moderate. They should cause less dramatic changes of the androgen (and downstream metabolite) concentrations, than found, for example, in human isolated 17,20-lyase deficiency  or in the cytochrome b5 knockout mouse , where androgens are massively depleted. Nevertheless, the CYB5A SNP surely contributes to the variance in androgen levels in the population, which will be inconspicuous in the majority of individuals. However, if activities of other enzymes involved in steroid metabolism are low, the rs1790834 (GG) carriers are more likely than carriers of the A allele to present signs of partial gonadal or adrenal failure [20,28,29,51,57,58]. This was suggested to be involved in the etiology of RA [15,16].
The CYB5A SNP contributes to the heritable risk in RA women. It can be an important factor contributing to the sex bias in RA incidence. As androgen levels in women are generally lower than in men, any additional decrease will disproportionately reduce the anti-inflammatory impact of androgens. In other words, the risk-lowering minor allele of SNP rs1790834 may help to ensure protective androgen levels in women, whereas the major allele may contribute to low androgen levels, which are insufficient to block development of RA or to ameliorate the inflammatory processes. In this setting the CYB5A SNP rs1790834 is rather unimportant in men, which is reflected by the lack of a protective role of the risk-lowering minor allele of SNP rs1790834 in the male subgroup. A clinical implication of our results lies in the fact that any kind of androgen therapy can be successful preferentially in CYB5A SNP rs1790834 major allele carriers.
two-dimensional thin-layer chromatography
antibody against cyclic citrullinated peptide specific for rheumatoid arthritis
enzyme-linked immunosorbent assay
genome-wide association studies
NADPH-cytochrome P450 reductase
real time polymerase chain reaction
single nucleotide polymorphisms
tumor necrosis factor
Scott DL, Wolfe F, Huizinga TW. Rheumatoid arthritis. Lancet. 2010;376:1094–108.
Eyre S, Bowes J, Diogo D, Lee A, Barton A, Martin P, et al. High-density genetic mapping identifies new susceptibility loci for rheumatoid arthritis. Nat Genet. 2012;44:1336–40.
Stahl EA, Raychaudhuri S, Remmers EF, Xie G, Eyre S, Thomson BP, et al. Genome-wide association study meta-analysis identifies seven new rheumatoid arthritis risk loci. Nat Genet. 2010;42:508–14.
Straub RH, Bijlsma JW, Masi A, Cutolo M. Role of neuroendocrine and neuroimmune mechanisms in chronic inflammatory rheumatic diseases – the 10-year update. Semin Arthritis Rheum. 2013;43:392–404.
Symmons D, Turner G, Webb R, Asten P, Barrett E, Lunt M, et al. The prevalence of rheumatoid arthritis in the United Kingdom: new estimates for a new century. Rheumatology (Oxford). 2002;41:793–800.
Schmidt M, Naumann H, Weidler C, Schellenberg M, Anders S, Straub RH. Inflammation and sex hormone metabolism. Ann NY Acad Sci. 2006;1069(236–46):236–46.
Cutolo M, Accardo S, Villaggio B, Barone A, Sulli A, Balleari E, et al. Androgen metabolism and inhibition of interleukin-1 synthesis in primary cultured human synovial macrophages. Mediat Inflamm. 1995;4:138–45.
Danenberg HD, Alpert G, Lustig S, Ben-Nathan D. Dehydroepiandrosterone protects mice from endotoxin toxicity and reduces tumor necrosis factor production. Antimicrob Agents Chemother. 1992;36:2275–9.
Daynes RA, Araneo BA, Ershler WB, Maloney C, Li GZ, Ryu SY. Altered regulation of IL-6 production with normal aging. Possible linkage to the age-associated decline in dehydroepiandrosterone and its sulfated derivative. J Immunol. 1993;150:5219–30.
Kanda N, Tsuchida T, Tamaki K. Testosterone suppresses anti-DNA antibody production in peripheral blood mononuclear cells from patients with systemic lupus erythematosus. Arthritis Rheum. 1997;40:1703–11.
Li ZG, Danis VA, Brooks PM. Effect of gonadal steroids on the production of IL-1 and IL-6 by blood mononuclear cells in vitro. Clin Exp Rheumatol. 1993;11:157–62.
Straub RH, Konecna L, Hrach S, Rothe G, Kreutz M, Schölmerich J, et al. Serum dehydroepiandrosterone (DHEA) and DHEA sulfate are negatively correlated with serum interleukin-6 (IL-6), and DHEA inhibits IL-6 secretion from mononuclear cells in man in vitro: possible link between endocrinosenescence and immunosenescence. J Clin Endocrinol Metab. 1998;83:2012–7.
Keller ET, Chang C, Ershler WB. Inhibition of NFkappaB activity through maintenance of IkappaB alpha levels contributes to dihydrotestosterone-mediated repression of the interleukin-6 promoter. J Biol Chem. 1996;271:26267–75.
Toyoda H, Takei S, Formby B. Effect of 5-alpha dihydrotestosterone on T-cell proliferation of the female nonobese diabetic mouse. Proc Soc Exp Biol Med. 1996;213:287–93.
Masi AT, Aldag JC, Jacobs JW. Rheumatoid arthritis: neuroendocrine immune integrated physiopathogenetic perspectives and therapy. Rheum Dis Clin North Am. 2005;31:131–60.
Pikwer M, Giwercman A, Bergstrom U, Nilsson JA, Jacobsson LT, Turesson C. Association between testosterone levels and risk of future rheumatoid arthritis in men: a population-based case–control study. Ann Rheum Dis. 2014;73:573–9.
Cutolo M, Balleari E, Giusti M, Intra E, Accardo S. Androgen replacement therapy in male patients with rheumatoid arthritis. Arthritis Rheum. 1991;34:1–5.
Booji A, Biewenga-Booji CM, Huber-Bruning O, Cornelis C, Jacobs JW, Bijlsma JW. Androgens as adjuvant treatment in postmenopausal female patients with rheumatoid arthritis. Ann Rheum Dis. 1996;55:811–5.
Straub RH, Pongratz G, Schölmerich J, Kees F, Schaible TF, Antoni C, et al. Long-term anti-tumor necrosis factor antibody therapy in rheumatoid arthritis patients sensitizes the pituitary gland and favors adrenal androgen secretion. Arthritis Rheum. 2003;48:1504–12.
Tengstrand B, Carlstrom K, Hafstrom I. Gonadal hormones in men with rheumatoid arthritis–from onset through 2 years. J Rheumatol. 2009;36:887–92.
Milewich L, Kaimal V, Toews GB. Androstenedione metabolism in human alveolar macrophages. J Clin Endocrinol Metab. 1983;56:920–4.
Schmidt M, Kreutz M, Löffler G, Schölmerich J, Straub RH. Conversion of dehydroepiandrosterone to downstream steroid hormones in macrophages. J Endocrinol. 2000;164:161–9.
Ishida Y, Killinger DW, Khalil MW, Yang K, Strutt B, Heersche JN. Expression of steroid-converting enzymes in osteoblasts derived from rat vertebrae. Osteoporos Int. 2002;13:235–40.
Castagnetta LA, Carruba G, Granata OM, Stefano R, Miele M, Schmidt M, et al. Increased estrogen formation and estrogen to androgen ratio in the synovial fluid of patients with rheumatoid arthritis. J Rheumatol. 2003;30:2597–605.
Schmidt M, Weidler C, Naumann H, Schölmerich J, Straub RH. Androgen conversion in osteoarthritis and rheumatoid arthritis synoviocytes - androstenedione and testosterone inhibit estrogen formation and favor production of more potent 5alpha-reduced androgens. Arthritis Res Ther. 2005;7:R938–48.
Dulos J, Verbraak E, Bagchus WM, Boots AM, Kaptein A. Severity of murine collagen-induced arthritis correlates with increased CYP7B activity: enhancement of dehydroepiandrosterone metabolism by interleukin-1beta. Arthritis Rheum. 2004;50:3346–53.
Dulos J, van der Vleuten MA, Kavelaars A, Heijnen CJ, Boots AM. CYP7B expression and activity in fibroblast-like synoviocytes from patients with rheumatoid arthritis: regulation by proinflammatory cytokines. Arthritis Rheum. 2005;52:770–8.
Martens HF, Sheets PK, Tenover JS, Dugowson CE, Bremner WJ, Starkebaum G. Decreased testosterone levels in men with rheumatoid arthritis: effect of low dose prednisone therapy. J Rheumatol. 1994;21:1427–31.
Miller WL, Auchus RJ. The molecular biology, biochemistry, and physiology of human steroidogenesis and its disorders. Endocr Rev. 2011;32:81–151.
Auchus RJ, Lee TC, Miller WL. Cytochrome b5 augments the 17,20-lyase activity of human P450c17 without direct electron transfer. J Biol Chem. 1998;273:3158–65.
Plenge RM, Seielstad M, Padyukov L, Lee AT, Remmers EF, Ding B, et al. TRAF1-C5 as a risk locus for rheumatoid arthritis – a genomewide study. N Engl J Med. 2007;20(357):1199–209.
Wellcome Trust Case Control Consortium. Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007;447:661–78.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31:315–24.
Miller LE, Jüsten HP, Schölmerich J, Straub RH. The loss of sympathetic nerve fibers in the synovial tissue of patients with rheumatoid arthritis is accompanied by increased norepinephrine release from synovial macrophages. FASEB J. 2000;14:2097–107.
Karouzakis E, Gay RE, Gay S, Neidhart M. Epigenetic control in rheumatoid arthritis synovial fibroblasts. Nat Rev Rheumatol. 2009;5:266–72.
Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29:e45.
Bresnihan B, Cunnane G, Youssef P, Yanni G, Fitzgerald O, Mulherin D. Microscopic measurement of synovial membrane inflammation in rheumatoid arthritis: proposals for the evaluation of tissue samples by quantitative analysis. Br J Rheumatol. 1998;37:636–42.
Wigginton JE, Cutler DJ, Abecasis GR. A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet. 2005;76:887–93.
Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–75.
Willer CJ, Li Y, Abecasis GR. METAL: fast and efficient meta-analysis of genomewide association scans. Bioinformatics. 2010;26:2190–1.
Cochran WG. The Combination of Estimates from Different Experiments. Biometrics. 1954;10:101–29.
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.
Frazer KA, Ballinger DG, Cox DR, Hinds DA, Stuve LL, Gibbs RA, et al. A second generation human haplotype map of over 3.1 million SNPs. Nature. 2007;449:851–61.
de Vries RR, van der Woude D, Houwing JJ, Toes RE. Genetics of ACPA-positive rheumatoid arthritis: the beginning of the end? Ann Rheum Dis. 2011;70:i51–4. doi: 10.1136/ard.2010.138040.:i51-i54.
Cutolo M, Balleari E, Giusti M, Monachesi M, Accardo S. Sex hormone status of male patients with rheumatoid arthritis: evidence of low serum concentrations of testosterone at baseline and after human chorionic gonadotropin stimulation. Arthritis Rheum. 1988;31:1314–7.
Masi AT, Bijlsma JW, Chikanza IC, Pitzalis C, Cutolo M. Neuroendocrine, immunologic, and microvascular systems interactions in rheumatoid arthritis: physiopathogenetic and therapeutic perspectives. Semin Arthritis Rheum. 1999;29:65–81.
Lee DM, Schur PH. Clinical utility of the anti-CCP assay in patients with rheumatic diseases. Ann Rheum Dis. 2003;62:870–4.
van der Woude D, Houwing-Duistermaat JJ, Toes RE, Huizinga TW, Thomson W, Worthington J, et al. Quantitative heritability of anti-citrullinated protein antibody-positive and anti-citrullinated protein antibody-negative rheumatoid arthritis. Arthritis Rheum. 2009;60:916–23.
Porter TD. The roles of cytochrome b5 in cytochrome P450 reactions. J Biochem Mol Toxicol. 2002;16:311–6.
Schenkman JB, Jansson I. The many roles of cytochrome b5. Pharmacol Ther. 2003;97:139–52.
Rainey WE, Nakamura Y. Regulation of the adrenal androgen biosynthesis. J Steroid Biochem Mol Biol. 2008;108:281–6.
Huber LC, Distler O, Tarner I, Gay RE, Gay S, Pap T. Synovial fibroblasts: key players in rheumatoid arthritis. Rheumatology (Oxford). 2006;45:669–75.
Katagiri M, Kagawa N, Waterman MR. The role of cytochrome b5 in the biosynthesis of androgens by human P450c17. Arch Biochem Biophys. 1995;317:343–7.
Vogl D, Falk W, Dorner M, Schölmerich J, Straub RH. Serum levels of pregnenolone and 17-hydroxypregnenolone in patients with rheumatoid arthritis and systemic lupus erythematosus: relation to other adrenal hormones. J Rheumatol. 2003;30:269–75.
Geller DH, Auchus RJ, Mendonca BB, Miller WL. The genetic and functional basis of isolated 17,20-lyase deficiency. Nat Genet. 1997;17:201–5.
McLaughlin LA, Ronseaux S, Finn RD, Henderson CJ, Roland WC. Deletion of microsomal cytochrome b5 profoundly affects hepatic and extrahepatic drug metabolism. Mol Pharmacol. 2010;78:269–78.
Miller WL. Steroidogenic enzymes. Endocr Dev. 2008;13:1–18. doi:10.1159/000134751.:1-18.
Rege J, Rainey WE. The steroid metabolome of adrenarche. J Endocrinol. 2012;214:133–43.
This study was supported by a grant from the Deutsche Forschungsgemeinschaft (Research Unit FOR 696, to KS with the number STA 620/2-1 and RHS with the number STR 511/25-1). We gratefully acknowledge the excellent technical assistance of Michaela Schmid, Josef Simon, Luise Rauch, Birgit Riepl, and Angelika Gräber.
The authors declare that they have no competing interests.
KS contributed to critically discussed the concept, conceived and designed the experiments, performed the experiments, contributed reagents/materials/analysis tools, analyzed the data, drafted the paper, generated tables, and gave final approval. RHS critically discussed the concept, conceived and designed the experiments, performed the experiments, contributed reagents/materials/analysis tools, analyzed the data, drafted the paper, enrolled patients, gave final approval, and is the corresponding author. JR contributed reagents/materials/analysis tools, contributed to the writing, enrolled patients, and gave final approval. SB contributed reagents/materials/analysis tools, contributed to the writing, enrolled patients, and gave final approval. GE conceived and designed the experiments, performed the experiments, analyzed the data, contributed to the writing, and gave final approval. MS developed of the novel concept, conceived and designed the experiments, performed the experiments, contributed reagents/materials/analysis tools, analyzed the data, drafted the paper, generated figures, and gave final approval. All authors read and approved the final manuscript.
Klaus Stark and Rainer H Straub contributed equally to this work.
About this article
Cite this article
Stark, K., Straub, R.H., Rovenský, J. et al. CYB5A polymorphism increases androgens and reduces risk of rheumatoid arthritis in women. Arthritis Res Ther 17, 56 (2015). https://doi.org/10.1186/s13075-015-0574-9
- Rheumatoid Arthritis
- Rheumatoid Arthritis Patient
- Single Nucleotide Polymorphism