- Research article
- Open Access
Methotrexate upregulates circadian transcriptional factors PAR bZIP to induce apoptosis on rheumatoid arthritis synovial fibroblasts
Arthritis Research & Therapy volume 20, Article number: 55 (2018)
Effects of methotrexate (MTX) on the proliferation of rheumatoid arthritis (RA) synovial fibroblasts are incompletely understood. We explored actions of MTX in view of circadian transcriptions of synovial fibroblasts.
Under treatment with MTX, expression of core circadian clock genes, circadian transcriptional factor proline and acidic amino acid-rich basic leucine zipper (PAR bZIP), and proapoptotic molecule Bcl-2 interacting killer (Bik) was examined by real-time polymerase chain reaction. Protein expression of circadian clock gene PERIOD2 (PER2) and CYTOCHROME C was also examined by western blotting and ELISA. Promoter activities of Per2 and Bik were measured by Luciferase assay. Expression of PER2, BIK, and CYTOCHROME C and morphological changes of the nucleus were observed by fluorescent immunostaining. Synovial fibroblasts were transfected with Per2/Bik small interfering RNA, and successively treated with MTX to determine cell viabilities. Finally, synovial fibroblasts were treated with MTX according to the oscillation of Per2/Bik expression.
MTX (10 nM) significantly decreased cell viabilities, but increased messenger RNA expression of Per2, Bik, and PAR ZIP including D site of the albumin promoter binding protein (Dbp), hepatic leukemia factor (Hlf), and thyrotroph embryonic factor (Tef). MTX also increased protein expression of PER2 and CYTOCHROME C, and promoter activities of Per2 and Bik via D-box. Under fluorescent observations, expression of PER2, BIK, and CYTOCHROME C was increased in apoptotic cells. Cytotoxicity of MTX was attenuated by silencing of Per2 and/or Bik, and revealed that MTX was significantly effective in situations where Per2/Bik expression was high.
We present here novel unique action of MTX on synovial fibroblasts that upregulates PAR bZIP to transcribe Per2 and Bik, resulting in apoptosis induction. MTX is important in modulating circadian environments to understand a new aspect of pathogenesis of RA.
Methotrexate (MTX) is a folic acid antagonist widely used as an anchor drug in treating various cancers [1, 2], as well as rheumatoid arthritis (RA) . For cancer cells, MTX competitively inhibits dihydrofolate reductase (DHFR) to block purine and pyrimidine biosynthesis, and thus it inhibits DNA replication and cell proliferation. For RA, low-dose MTX shows anti-inflammatory effects by inducing extracellular adenosine, which binds to adenosine receptors . Moreover, it has been reported that MTX induced apoptosis in synovial fibroblasts in both in-vivo and in-vitro experiments . Although MTX induced apoptosis in synovial fibroblast within 24–48 h [6, 7], the precise mechanism of how MTX expresses antiproliferative effects on synovial fibroblasts remains incompletely understood .
RA is a chronic arthritis characterized by ‘tumor-like’ synovial cell growth . Another remarkable feature of RA is the circadian variation of disease-related symptoms, such as morning stiffness, increased production of proinflammatory cytokines at night time, and peaked secretion of immunoglobulin (Ig) A/IgM types of rheumatoid factor in the morning [10,11,12,13,14,15]. Since these rheumatic symptoms possess a daily rhythm, we have previously shown that the action of the biological clock was significantly disturbed in the mouse model of collagen antibody-induced arthritis  and that tumor necrosis factor (TNF)-α significantly disturbed the oscillation of biological clocks of synovial fibroblasts . Fibroblasts usually demonstrate daily rhythms of circadian clock, while Haas et al.  pointed out that the expression rhythm of clock genes disappears in RA synovial fibroblasts presumably due to prolonged inflammation.
The circadian rhythm in human cells is mainly regulated by the core clock genes, including circadian locomotor output cycles kaput (Clock), brain and muscle Rant-like protein-1 (Bmal1), period (Per), and cryptochrome (Cry) [19,20,21,22]. The circadian transcriptional factor proline and acidic amino acid-rich basic leucine zipper (PAR bZIP) includes the D site of the albumin promoter binding protein (Dbp), hepatic leukemia factor (Hlf), and thyrotroph embryonic factor (Tef). PAR bZIP regulates gene expression by binding to a consensus sequence of D-box (5′-TTAXGTAA-3′; X = T or C) on the promoter region [23,24,25]. It has been reported that PAR bZIP can regulate the transcription of Per2 gene by binding two D-box sequences existing on its promoter (D-box 1, 5′-TTATGTAA-3′, −372 to −365; and D-box 2, 5′-TTACGTAA-3′, −47 to −40) . In contrast, E4-binding protein 4 (E4bp4) also binds to the D-box to suppress the transcription of Per2 [26,27,28] (see Additional file 1).
It is noted that Bcl-2 interacting killer (Bik) possesses D-box (5′-TTAAGTCA-3′, −285 to −277) on its promoter region , resembling the arrangement of Per2 genes. Bik, a member of the BH3-only subfamily, acts as an important signaling molecule upstream of the Bcl-2 and Bax subfamily . The BCL-2 family has been known to be central players in regulating physiological activities of mitochondria, with Bcl-2, Bcl-xL, and Mcl-1 suppressing mitochondria-related apoptosis, whereas Bax and Bak induce it. Among these, Bik directly binds to these family proteins to induce apoptosis [29, 30].
In this study, we explored novel pharmacological effects of MTX on circadian clock genes and apoptosis induction in RA synovial fibroblasts.
Synovial fibroblast culture
Synovial tissues were obtained from RA patients (eight females and two males; aged 59.3 ± 3.4 years) during joint surgery. Eight of 10 patients took methotrexate and their average dosage was 7.4 ± 1.4 mg/week (Table 1) (see Additional file 2). This study has been approved by the ethics committee of Kobe University Graduate School of Health Sciences (#579–1) and Kobe Kaisei Hospital (#0072), in accordance with the Declaration of Helsinki. Written informed consent was obtained from each patient before study enrolment. Tissues were minced and treated with 2 mg/ml collagenase (Wako, Tokyo, Japan) in serum-free Dulbecco’s modified Eagle’s medium (DMEM; Nissui, Tokyo, Japan) at 37 °C for 1 h. Primary cultured synovial cell lines were established and maintained in DMEM including 10% heat-inactivated fetal bovine serum (FBS; Thermo, Waltham, MA, USA), 1% penicillin–streptomycin (100 U/ml; Life Technologies, Carlsbad, CA, USA), and 1% l-glutamine (2 mM; Life Technologies), in a humidified incubator at 37 °C in the presence of 5% CO2. Harvested cells were continuously cultured to obtain synovial fibroblasts, and cells of passages 3–6 were used in the entire experiments.
Synchronization of synovial fibroblasts by serum shock
Cell viability assay
Synovial fibroblasts (3.0 × 103 cells) were cultured in serum-free DMEM with or without MTX (1, 10, and 100 nM, 1 μM; Sigma Aldrich, St. Louis, MO, USA). After incubating for 24–72 h, cell viabilities were measured as the 450-nm absorbance of reduced WST-8 (2-(2-methoxy-4-nitrophenyl)-3-(4-nitriphenyl)–5-(2, 4-disulphonyl)-2H–tetrazolium, monosodium salt) using the Cell Counting Kit-8 (Dojindo, Kumamoto, Japan). The values were represented relative to MTX-untreated cells (control).
Real-time polymerase chain reaction
Synovial fibroblasts (8.0 × 104 cells) were cultured in serum-free DMEM with or without MTX (10, 100 nM) for 24–32 h. After incubation, total RNA was extracted using the RNeasy Mini Kit (QIAGEN, Hilden, Germany). Then, reverse transcription was performed with the ReverTra Ace® qPCR RT Kit (Toyobo, Osaka, Japan) and analyzed on the StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The TaqMan probes used were: Hs00154147_m1 for Bmal1, Hs00231857_m1 for Clock, Hs00256143_m1 for Per2, Hs00172734_m1 for Cry1, Hs00609747_m1 for Dbp, Hs00171406_m1 for Hlf, Hs01115720_m1 for Tef, Hs00993282_m1 for E4bp4, Hs00154189_m1 for Bik, and Hs00427621_m1 for TATA box binding protein (Tbp). Expression levels were normalized to Tbp.
Synovial fibroblasts (5.0 × 105 cells) were cultured in serum-free DMEM with or without MTX (10, 100 nM) for 24–48 h, and lysed with RIPA buffer (Wako) to obtain cytoplasmic proteins. Samples were subjected to sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE), transferred to polyvinylidene difluoride (PVDF) membrane (Millipore, Bedford, MA, USA), probed with antibodies, and developed by ImmunoStar® LD (Wako). Antibodies (Abs) used were: anti-PER2 Ab (sc-101,105; Santa Cruz, Dallas, TX, USA), anti-CYTOCHROME C Ab (ab13575; Abcam, Cambridge, UK), anti-BIK Ab (NB100–56109; Novus, Littleton, CO, USA), anti-Actin Ab (sc-1616; Santa Cruz), anti-mouse IgG Ab (NA9310; GE Healthcare, Chicago, IL, USA), and anti-rabbit IgG Ab (NA9340; GE Healthcare).
ELISA for PER2 and CYTOCHROME C
Synovial fibroblasts (3.0 × 103 cells) were cultured in serum-free DMEM with or without MTX (10, 100 nM) for 32–48 h, and protein expression of PER2 and CYTOCHROME C was measured by the In Cell ELISA test kit (3440–02; Thermo). Antibodies were as already described.
Luciferase reporter gene construction
Genomic DNA was obtained from RA synovial fibroblasts and amplified by PCR to generate the luciferase reporters. The primer pairs used for amplifying human Bik promoter containing D-box (−780 to +176) were 5′-TGGCCTAACTGGCCGTAAACAAGCTTTGCCGTGC-3′ (forward) and 5′-CGCCGAGGCCAGATCATGCTGGCAGCGTCTGTA-3′ (reverse). The primer pairs used for amplifying human Bik promoter without D-box (−260 to + 323) were 5′-TGGCCTAACTGGCCGCCTCTTGGAGCCTCGGTT-3′ (forward) and 5′-CGCCGAGGCCAGATCTTGCTGGAGCGGTAAAACC-3′ (reverse). The KpnI recognition sequence was added to the 5′ ends of forward primers and the BglII recognition sequence was added to the 5′ ends of reverse primers. The PCR products were cloned into the KpnI and BglII site of the pGL4.10 (luc2) vector (Promega, Madison, WI, USA) using In-Fusion® HD Cloning Plus (Takara, Shiga, Japan). The luciferase reporters containing Per2 promoters were generated with reference to Yoshida et al. . D-box motifs of Per2 promoter were mutated from 5′-TTATGTAA-3′ to 5’-CGCCAGGC-3′ (−372 to −365) and from 5′-TTACGTAA-3′ to 5′-CAGCGTAA-3′ (−47 to −40) (see Additional file 3).
Transient transfection and luciferase reporter assay
Synovial fibroblasts (4 × 104 cells) were transfected with 500 ng of the pGL4.10 (luc2) vector containing various Per2 and Bik promoters using Lipofectamine 3000 Transfection Reagent (Thermo). As an internal control, 35 ng of pRL-TK (Promega) containing the herpes simplex virus thymidine kinase promoter driving Renilla luciferase was cotransfected. After 24 h of incubation for transfection, cells were treated with 10 and 100 nM MTX for 24 h, and analyzed for luciferase activity using the Dual-Luciferase Reporter Assay (Promega). Activities of both firefly and Renilla luciferases were measured, and the activity of firefly luciferase was normalized by Renilla luciferase. The values were shown as relative variations to MTX-untreated cells.
Synovial fibroblasts (6.0 × 103 cells) were cultured in serum-free DMEM with or without MTX (10 nM) for 24 h. Then, cells were fixed with 4% formaldehyde, stained with anti-PER2 Ab, anti-BIK Ab, anti-CYTOCHROME C Ab, anti-mouse IgG (H + L), F(ab′)2 Fragment (Alexa Fluor® 594 Conjugate) (#8890; Cell Signaling Technology, Danvers, MA, USA), anti-rabbit IgG (H + L), F(ab′)2 Fragment (Alexa Fluor® 488 Conjugate) (#4412; Cell Signaling Technology), and DAPI (0.5 μg/ml; Sigma, St. Louis, MO, USA). Protein expression and morphological changes of the nucleus were examined under fluorescence microscopy.
Per2 small interfering (si) RNA (s16931; Life Technologies) and Bik siRNA (s1990; Life Technologies) were transfected into synovial fibroblasts (3.0 × 103 cells) using Lipofectamine™ RNAiMAX (Life Technologies) for 48 h. Noncoding siRNA (4,390,843; Life Technologies) was also used as controls. After that, cells were cultured in serum-free DMEM with 10 nM of MTX for 24 h to measure cell viabilities using the Cell Counting Kit-8. The viabilities were represented relative to those of noncoding siRNA without MTX treatment.
In dynamite-plunger plots, values were expressed as the mean ± standard error of the mean (SEM). In boxplots, values were expressed with 10th, 25th, 50th (median), 75th, and 90th percentiles, and the single values were superimposed on the boxplots using black symbols.
For statistical analyses, a one-sample Kolmogorov–Smirnov test was used to test the normality, A one-sample t test was used to compare one single control value and others, a paired t test was used to compare differences between two experimental groups, the Tukey test was used to compare differences between more than three experimental groups, and Dunnett’s test was used to compare differences between the control and others. All statistical tests were two-sided and p < 0.05 was considered statistically significant.
The statistical analyses were performed by EZR, based on R and R commander .
MTX inhibited viabilities of synovial fibroblasts
We first examined the effect of MTX on viabilities of synovial fibroblasts. A lower concentration of MTX (1, 10 nM) significantly decreased cell viabilities, while cell viability was relatively weak with a higher concentration of MTX (100 nM, 1 μM) (Fig. 1). Since 10 or 100 nM was the optimum MTX concentration in human sera , experiments were subsequently conducted with MTX 10 nM treatment as an appropriate effective concentration and a concentration of 100 nM as a less-effective control.
MTX accelerated the expression of circadian clock gene Per2
We next determined the effect of MTX on expressions of circadian clock genes. Under treatment with MTX (10 nM for 24 h), the expression of Per2 messenger RNA (mRNA) was significantly increased, while expression of Bmal1, Clock, and Cry1 mRNA was not. In contrast, under treatment with MTX (100 nM for 24 h), the expression of Per2 did not show any significant differences, while the expression of Clock and Cry1 significantly decreased (Fig. 2a). In accordance with this, the expression of PER2 protein was significantly increased by 10 nM of MTX for 32 h, but not by 100 nM (Fig. 2b, c).
MTX increased the expression of PAR bZIP and Bik
The circadian transcriptional factor PAR bZIP (Dbp, Tef, and Hlf) and E4bp4 genes regulate Per2 gene expression by binding to D-box on the promoter region [23,24,25,26,27,28]. We therefore examined the effect of MTX on the mRNA expression of Dbp, Tef, Hlf, and E4bp4. Expression of Dbp, Tef, and Tef mRNA was significantly increased by 10 nM MTX treatment for 24 h, while E4bp4 mRNA expression was not. However, expression of Dbp and E4bp4 significantly decreased with 100 nM of MTX treatment (Fig. 3a).
Next, we examined the mRNA expression of the proapoptotic factor Bik, since the PAR bZIP-binding site also exists on the promoter region of Bik . As shown in Fig. 3b, the expression of Bik mRNA was significantly increased by 10 nM of MTX treatment for 32 h.
Finally, we observed cytosolic release of CYTOCHROME C since Bik interacted with BCL-2 family proteins to induce mitochondria-related apoptosis. As a result, expression of CYTOCHROME C was increased by 10 nM of MTX treatment for 48 h, but not by 100 nM (Fig. 3c, d).
Transcription of Per2 and Bik was enhanced via D-box
To examine the roles of PAR bZIP in the transcription of Per2 and Bik, we transfected luciferase reporters, constructed with or without D-box, to synovial fibroblasts. With 10 nM of MTX treatment for 24 h, D-box(+) Per2 promoter activity was significantly increased as compared to those of the control and D-box(−) Per2 promoter. By 100 nM of MTX treatment, D-box(−) Per2 promoter activity was significantly decreased as compared to that of the control. Bik promoter activity was significantly increased after 10 nM of MTX treatment, and was higher than that of D-box(−) Bik promoter (Fig. 4).
PER2, BIK, and CYTOCHROME C were highly expressed in apoptotic cells
To certify the specific expression of Per2 and Bik on MTX-induced apoptotic synovial fibroblasts, we used fluorescent immunostaining to examine expressions of PER2, BIK, and CYTOCHROME C and, simultaneously, morphological changes of nucleus. As a result, nuclear blebbing, shrinkage, or fragmentation were observed indicating that MTX (10 nM for 24 h) could induce cell death. PER2 was overexpressed on apoptotic cells, but not on cells with normal appearance (Fig. 5a). BIK and CYTOCHROME C were also overexpressed specifically on the apoptotic cells (Fig. 5b).
Per2 and Bik are critical genes for MTX-induced apoptosis
To confirm whether PER2 and BIK were responsible for MTX-induced cell death, Per2 and Bik siRNA were transfected into synovial fibroblasts to examine cell viabilities. As shown in Fig. 6a, cytotoxicity of MTX was attenuated by silencing of Per2 and Bik, independently and synergistically.
When we examined the long-term expression of Per2 and Bik mRNA (Fig. 6b), they gradually increased from 8 h until 24 or 32 h after synchronized circadian oscillations by serum shock, consistent with our previous report that the expression of Per2 on synovial fibroblasts was increased after synchronization and peaked at 24 h . Since these genes were responsible for MTX-induced cell death, synovial fibroblasts were separately treated by 10 nM of MTX at 0–8 h (the lower Per2/Bik expression period) and at 24–32 h (the higher Per2/Bik expression period), and viabilities of synovial fibroblasts were again measured. As shown in Fig. 6c, the reduction rate of cell viabilities was significantly increased by 24–32 h of MTX treatment compared with 0–8 h of treatment, depending on the expression levels of Per2 and Bik.
We present here a novel pharmacological action of MTX in the viewpoint of circadian clock genes, circadian transcriptional factor PAR bZIP, and proapoptotic molecule Bik. In addition, a crucial mechanism of MTX-induced apoptosis on rheumatoid synovial fibroblasts was further traced.
In this study, the effect of MTX on cell viabilities did not show a concentration-dependent manner, with a wide range of MTX concentrations from 1 nM to 1 μM. In addition to our result, we found that MTX concentrations lower than 1 nM did not decrease cell viability (see Additional file 4). MTX is transferred almost equally to synovial fluid and sera, and the concentration of MTX in human sera can reach 200 nM at the peak and immediately decreases to less than 10 nM , although pharmacological effects on cytotoxicity or cellular viability may not necessarily increase in a concentration-dependent manner. As reported previously, viabilities of human lymphoblastic leukemia cells and epithelial cells of rat decreased in a concentration-dependent manner at low MTX concentrations, but did not show a significant decrease at high MTX concentrations [34, 35]. Moreover, MTX did not affect viability in a concentration-dependent manner on T cells and RA synovial fibroblasts, in contrast to those of dose dependency observed in osteoarthritis synovial fibroblasts [6, 36]. Thus, further studies are required for the action of MTX correlated with drug dosage or cell types.
Since MTX concentrations below 100 nM were conceivable in sera of RA patients, we next focused on the effect of MTX on circadian clock genes and circadian transcriptional factor PAR bZIP genes and their relations to mitochondria-related apoptosis of synovial fibroblasts.
Since circadian clock genes were reported to be closely related to the pathogenesis of arthritis [16, 37] and excessive expression of Per2 could induce apoptosis , we examined the effect of MTX on expression of Per2, Bmal1, Clock, and Cry1 that were regarded as “core” clock genes [20,21,22]. We examined mRNA expression of circadian clock genes over time, and found that the controls and 10/100 nM of MTX showed almost the same expression rhythms, and MTX influenced their expression levels (see Additional file 5).
As described, both Per2 and Bik genes have D-box in their promoter regions, and PAR bZIP proteins regulate the expression of these genes by binding to D-box [24, 39]. Indeed, promoter activities of Per2 and Bik were upregulated by 10 nM of MTX when the D-box sequence exists and were cancelled without D-box. Moreover, both PER2 and BIK were highly expressed in MTX-induced apoptotic cells, while inhibitions of Per2 and Bik synergistically attenuated the effect of MTX on cellular viabilities. As it has been reported that PAR bZIP and Bik mediate oxidative stress-induced apoptosis in fibroblasts , we consider Per2 and Bik as essential factors for MTX-dependent synovial cell death and propose here that two independent pathways can mediate these death signals: the PAR bZIP–Per2 transcriptional pathway and the PAR bZIP–Bik transcriptional pathway (see Additional file 6). However, as shown in Fig. 2a, downregulation of Clock and Cry1 appeared to be less effective on cell viabilities after treatment with 100 nM of MTX. Further in-vivo study should be required for clinical application by investigating the cross-talk of clock genes.
Last, because of the similarity of their promoter region containing D-box, we supposed Bik might show oscillation as has been reported for Per2 [17, 31]. For the further understanding of circadian manifestation in rheumatoid arthritis, we showed that the expression of Bik, as well as Per2, was gradually increased from 8 h until 24 or 32 h, and MTX was significantly effective in situations when Per2 and Bik were highly expressed. It has been reported that administration of MTX at bedtime, as an optimal dosing time associated with the oscillation of TNF-α production, could reduce the disease activities of patients with RA [40, 41]. Thus, we considered that expression levels of Per2 and Bik could also be a critical biomarker for chronotherapy of MTX not only for RA, but also for other diseases such as acute lymphocytic leukemia, non-Hodgkin’s lymphoma, osteosarcoma, and breast cancer, in accordance with a previous report that circadian clock genes could be a biomarker for chronotherapy .
We propose here that MTX induces apoptosis in synovial fibroblasts through the binding of PAR bZIP to D-box of two different genes, Per2 and Bik, and these dual pathways work independently but synergistically. In consequence, MTX acts as a modulator of circadian environments of synovial fibroblasts in relation to biological clock genes and circadian transcriptional factors.
Bcl-2 interacting killer
Brain and muscle Arnt-like protein-1
Circadian locomotor output cycles kaput
D site of the albumin promoter binding protein
Dulbecco’s modified Eagle’s medium
E4-binding protein 4
Fetal bovine serum
Hepatic leukemia factor
- PAR bZIP:
Proline and acidic amino acid-rich basic leucine zipper
Polymerase chain reaction
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
Standard error of the mean
Small interfering RNA
TATA box binding protein
Thyrotroph embryonic factor
Tumor necrosis factor
Jolivet J, Cowan KH, Curt GA, Clendeninn NJ, Chabner BA. The pharmacology and clinical use of methotrexate. N Engl J Med. 1983;309:1094–104.
Bleyer WA. The clinical pharmacology of methotrexate: new applications of an old drug. Cancer. 1978;41:36–51.
Cronstein BN. Low-dose methotrexate: a mainstay in the treatment of rheumatoid arthritis. Pharmacol Rev. 2005;57:163–72.
Cutolo M, Sulli A, Pizzorni C, Seriolo B, Straub RH. Anti-inflammatory mechanisms of methotrexate in rheumatoid arthritis. Ann Rheum Dis. 2001;60:729–35.
Nakazawa F, Matsuno H, Yudoh K, Katayama R, Sawai T, Uzuki M, et al. Methotrexate inhibits rheumatoid synovitis by inducing apoptosis. J Rheumatol. 2001;28:1800–8.
Xu K, Cai YS, Lu SM, Li XL, Liu L, Li Z, et al. Autophagy induction contributes to the resistance to methotrexate treatment in rheumatoid arthritis fibroblast-like synovial cells through high mobility group box chromosomal protein 1. Arthritis Res Ther. 2015;17:374.
Lee SY, Park SH, Lee SW, Lee SH, Son MK, Choi YH, et al. Synoviocyte apoptosis may differentiate responder and non-responder patients to methotrexate treatment in rheumatoid arthritis. Arch Pharm Res. 2014;37:1286–94.
Spurlock CF 3rd, Gass HM 4th, Bryant CJ, Wells BC, Olsen NJ, Aune TM. Methotrexate-mediated inhibition of nuclear factor κB activation by distinct pathways in T cells and fibroblast-like synoviocytes. Rheumatology. 2014;54:178–87.
Olsen NJ, Stein CM. New drugs for rheumatoid arthritis. N Engl J Med. 2004;350:2167–79.
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.
Cutolo M, Straub RH, Buttgereit F. Circadian rhythms of nocturnal hormones in rheumatoid arthritis: translation from bench to bedside. Ann Rheum Dis. 2008;67:905–8.
Danis VA, Franic GM, Rathjen DA, Laurent RM, Brooks PM. Circulating cytokine levels in patients with rheumatoid arthritis: results of a double blind trial with sulphasalazine. Ann Rheum Dis. 1992;51:946–50.
Dasgupta B, Corkill M, Kirkham B, Gibson T, Panayi G. Serial estimation of interleukin 6 as a measure of systemic disease in rheumatoid arthritis. J Rheumatol. 1992;19:22–5.
Choy EH, Panayi GS. Cytokine pathways and joint inflammation in rheumatoid arthritis. N Engl J Med. 2001;344:907–16.
Straub RH, Cutolo M. Circadian rhythms in rheumatoid arthritis: implications for pathophysiology and therapeutic management. Arthritis Rheum. 2007;56:399–408.
Hashiramoto A, Yamane T, Tsumiyama K, Yoshida K, Komai K, Yamada H, et al. Mammalian clock gene Cryptochrome regulates arthritis via proinflammatory cytokine TNF-α. J Immunol. 2010;184:1560–5.
Yoshida K, Hashiramoto A, Okano T, Yamane T, Shibanuma N, Shiozawa S. TNF-α modulates expression of the circadian clock gene Per2 in rheumatoid synovial cells. Scand J Rheumatol. 2013;42:276–80.
Haas S, Straub RH. Disruption of rhythms of molecular clocks in primary synovial fibroblasts of patients with osteoarthritis and rheumatoid arthritis, role of IL-1β/TNF. Arthritis Res Ther. 2012;14:R122.
Cutolo M, Villaggio B, Otsa K, Aakre O, Sulli A, Seriolo B. Altered circadian rhythms in rheumatoid arthritis patients play a role in the disease’s symptoms. Autoimmun Rev. 2005;4:497–502.
Gibbs JE, Ray DW. The role of the circadian clock in rheumatoid arthritis. Arthritis Res Ther. 2013;15:205.
Cuninkova L, Brown SA. Peripheral circadian oscillators: interesting mechanisms and powerful tools. Ann N Y Acad Sci. 2008;1129:358–70.
Sato TK, Yamada RG, Ukai H, Baggs JE, Miraglia LJ, Kobayashi TJ, et al. Feedback repression is required for mammalian circadian clock function. Nat Genet. 2006;38:312–9.
Falvey E, Marcacci L, Schibler U. DNA-binding specificity of PAR and C/EBP leucine zipper proteins: a single amino acid substitution in the C/EBP DNA-binding domain confers PAR-like specificity to C/EBP. Biol Chem. 1996;377:797–809.
Li S, Hunger SP. The DBP transcriptional activation domain is highly homologous to that of HLF and TEF and is not responsible for the tissue type-specific transcriptional activity of DBP. Gene. 2001;263:239–45.
Cowell IG, Skinner A, Hurst HC. Transcriptional repression by a novel member of the bZIP family of transcription factors. Mol Cell Biol. 1992;12:3070–7.
Mitsui S, Yamaguchi S, Matsuo T, Ishida Y, Okamura H. Antagonistic role of E4BP4 and PAR proteins in the circadian oscillatory mechanism. Genes Dev. 2001;15:995–1006.
Yamajuku D, Shibata Y, Kitazawa M, Katakura T, Urata H, Kojima T, et al. Cellular DBP and E4BP4 proteins are critical for determining the period length of the circadian oscillator. FEBS Lett. 2011;585:2217–22.
Ohno T, Onishi Y, Ishida N. A novel E4BP4 element drives circadian expression of mPeriod2. Nucleic Acids Res. 2007;35:648–55.
Rashmi R, Pillai SG, Vijayalingam S, Ryerse J, Chinnadurai G. BH3-only protein BIK induces caspase-independent cell death with autophagic features in Bcl-2 null cells. Oncogene. 2008;27:1366–75.
Boyd JM, Gallo GJ, Elangovan B, Houghton AB, Malstrom S, Avery BJ, et al. Bik, a novel death-inducing protein shares a distinct sequence motif with Bcl-2 family proteins and interacts with viral and cellular survival-promoting proteins. Oncogene. 1995;11:1921–8.
Balsalobre A, Damiola F, Schibler U. A serum shock induces circadian gene expression in mammalian tissue culture cells. Cell. 1998;93:929–37.
Kanda Y. Investigation of the freely available easy-to-use software ‘EZR’ for medical statistics. Bone Marrow Transplant. 2013;48(3):452–8.
Shiozawa K, Tanaka Y, Yoshihara R, Imura S, Murata M, Yamane T, et al. Serum levels and pharmacodynamics of methotrexate and its metabolite 7-hydroxy methotrexate in Japanese patients with rheumatoid arthritis treated with 2-mg capsule of methotrexate three times per week. Mod Rheumatol. 2005;15:405–9.
Haber M, Madafiglio J, Norris MD. Methotrexate cytotoxicity determination using the MTT assay following enzymatic depletion of thymidine and hypoxanthine. J Cancer Res Clin Oncol. 1993;119(6):315–7.
Horie T, Li T, Ito K, Sumi S, Fuwa T. Aged garlic extract protects against methotrexate-induced apoptotic cell injury of IEC-6 cells. J Nutr. 2006;136(3 Suppl):861S–3S.
Genestier L, Paillot R, Fournel S, Ferraro C, Miossec P, Revillard JP. Immunosuppressive properties of methotrexate: apoptosis and clonal deletion of activated peripheral T cells. J Clin Invest. 1998;102(2):322–8.
Hand LE, Hopwood TW, Dickson SH, Walker AL, Loudon AS, Ray DW, et al. The circadian clock regulates inflammatory arthritis. FASEB J. 2016;30(11):3759–70.
Hua H, Wang Y, Wan C, Liu Y, Zhu B, Yang C, et al. Circadian gene mPer2 overexpression induces cancer cell apoptosis. Cancer Sci. 2006;97:589–96.
Ritchie A, Gutierrez O, Fernandez-Luna JL. PAR bZIP-bik is a novel transcriptional pathway that mediates oxidative stress-induced apoptosis in fibroblasts. Cell Death Differ. 2009;16:838–46.
To H, Irie S, Tomonari M, Watanabe Y, Kitahara T, Sasaki H. Therapeutic index of methotrexate depends on circadian cycling of tumour necrosis factor-alpha in collagen-induced arthritic rats and mice. J Pharm Pharmacol. 2009;61:1333–8.
To H, Yoshimatsu H, Tomonari M, Ida H, Tsurumoto T, Tsuji Y, et al. Methotrexate chronotherapy is effective against rheumatoid arthritis. Chronobiol Int. 2011;28:267–74.
Tang Q, Cheng B, Xie M, Chen Y, Zhao J, Zhou X, et al. Circadian clock gene Bmal1 inhibits tumorigenesis and increases paclitaxel sensitivity in tongue squamous cell carcinoma. Cancer Res. 2017;77:532–44.
The authors declare that they have no funding.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study has been approved by the ethics committee of Kobe University Graduate School of Health Sciences (#579–1) and Kobe Kaisei Hospital (#0072), in accordance with the Declaration of Helsinki. Written informed consent was obtained from each patient before study enrolment.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The interactions of circadian clock genes and its relative factors. BMAL/CLOCK and PER/CRY create the circadian rhythm via E-box, and DBP, TEF, HLF, E4BP4 via D-box. BMAL and CLOCK heterodimerize and bind to E-box elements on promoter regions of Per and Cry genes to induce their transcription. Thereafter, PER and CRY proteins heterodimerize to inhibit activities of their own and other E-box regulated promoters. (PDF 214 kb)
Detailed information of patients with RA enrolled in this study. Joint samples obtained from 10 different patients to establish primary cultured synovial fibroblast cell lines. Age, sex, disease duration, CRP, DAS28-ESR, and treatment with MTX, PSL, and other DMARDs shown in the table. (PDF 252 kb)
Plasmid constructs of Per2/Bik promoter. D-box(+), plasmid constructs containing D-box. D-box(−), plasmid constructs without D-box. D-box motifs of Per2 promoter were mutated from 5′-TTATGTAA-3′ to 5′-CGCCAGGC-3′ (−372 to −365), and 5′-TTACGTAA-3′ to 5′-CAGCGTAA-3′ (−47 to −40). Human Bik promoter containing D-box (−780 to +176) and human Bik promoter without D-box (−260 to +323) constructed. (PDF 215 kb)
Cell viability of MTX-treated fibroblasts. Cell viability of RA synovial fibroblasts measured by WST-8 assay after 24 h of stimulation of MTX (1 pM to 1 μM). MTX (1 and 10 nM) significantly decreased cell viability as shown in Fig. 1, while 1–100 pM of MTX did not. (PDF 274 kb)
The mRNA expression of circadian clock genes over time. mRNA expression of circadian clock genes measured at –4 h (before synchronization), 0 h (just before MTX stimulation), 24 h, 32 h, and 48 h. Controls and 10/100 nM of MTX showed almost the same expression rhythms, and MTX influenced their expression levels. (PDF 264 kb)
Two different transcriptional pathways by which MTX induces apoptosis to synovial fibroblasts: PAR bZIP–Per2 transcriptional pathway and PAR bZIP–Bik transcriptional pathway. We propose MTX induces apoptosis in synovial cells through activated binding of PAR bZIP to D-box in two different genes, Per2 and Bik, and these dual pathways work independently but synergistically. (PDF 206 kb)
About this article
Cite this article
Suzuki, K., Yoshida, K., Ueha, T. et al. Methotrexate upregulates circadian transcriptional factors PAR bZIP to induce apoptosis on rheumatoid arthritis synovial fibroblasts. Arthritis Res Ther 20, 55 (2018). https://doi.org/10.1186/s13075-018-1552-9
- Rheumatoid arthritis
- Synovial fibroblasts
- Proline and acidic amino acid-rich basic leucine zipper
- Bcl-2 interacting killer