Prenatal Dexamethasone Exposure Increases the Susceptibility to Autoimmunity in Offspring Rats by Epigenetic Programing of Glucocorticoid Receptor.

Objective. Prenatal glucocorticoids (GC) can induce long term effects on offspring health. However, reports and related studies regarding the prolonged effects of prenatal GC on the development of autoimmunity are limited. Here, we aimed to explore the immunological effects of dexamethasone (DEX) exposure on young adults and whether glucocorticoid receptor (GR) is involved in this process. Methods. Wistar rats were given DEX during pregnancy. Susceptibility to autoimmunity in offspring was assessed using experimental autoimmune encephalomyelitis (EAE) and adjuvant-induced arthritis (AIA) animal models. To reveal the possible mechanism, glucocorticoid response, GR expression, and methylation status were measured in peripheral blood mononuclear cells (PBMCs). Results. Our results showed that the DEX-treated rats had greater susceptibility to EAE (100% versus 62.5%, P < 0.05) and AIA (63.6% versus 0%, P < 0.05) than saline control group. Glucocorticoid response and GR expression were decreased in DEX rats. Significant difference was also found in the methylation levels of GR exon 1-10 to exon 1-11 region. Conclusions. Prenatal DEX administration increases the susceptibility to autoimmune diseases, which is potentially mediated by programming GR methylation status and glucocorticoid sensitivity.

1. Introduction

Glucocorticoids are essential for the maturation of fetal organ systems and normal development. During the course of human pregnancy, glucocorticoid (GC) treatment is often given when preterm delivery is expected. This treatment is successful in maturing the fetal lung and decreasing the occurrence of respiratory distress syndrome. However, emerging evidence indicates that antenatal glucocorticoids may elicit fetal programming that leads to adverse effects in adulthood, like hypertension, metabolic syndrome, and allergy [1-3]. This intergenerational transmission of the GC effects is mainly mediated by epigenetic mechanisms [2, 4, 5]. In animals, the effects of prenatal GC treatment on the development of the immune system have been evaluated in small series, which indicated that prenatal GC treatment may result in the alteration of inflammatory factors [6, 7], and thus contribute to the disturbance of immune system. However, reports about the mechanism linked prenatal GC exposure and autoimmunity are few. GC and glucocorticoid receptor (GR) signaling play a critical role in regulating the immune system. Our previous studies indicated that autoimmunity is closely related with glucocorticoid resistance and GR dysfunction, which potentially serve as the mechanism of autoimmune disease [8, 9]. Taking these clues together, we hypothesize that prenatal GC exposure may contribute to the etiology of autoimmune disease through epigenetic modification of GR and glucocorticoid response [9]. The current study was designed to examine our hypothesis, which may provide insight into the early pathogenicity of certain autoimmune diseases (e.g., systemic lupus erythematosus, multiple sclerosis, and diabetes) and facilitate the development of preventative and therapeutic strategies.

2. Materials and Methods

2.1. Study Design

All experiments were performed in accordance with guidelines developed by the Chinese Council on Animal Care, and the protocol was approved by Sun Yat-sen University Animal Care Committee. Wistar rats (200–250 g) and guinea pigs (300–500 g) were obtained from Sun Yat-sen University (Guangzhou, China) and maintained in an animal facility at 25 ± 3°C, with a relative humidity of 55 ± 3%, in a 12 h light/12 h dark cycle (the light is turned on at 8 am and turned off at 8 pm), with ad libitum access to food and tap water. Male Wistar rats were caged with individual females until mating was confirmed. After mating, female rats were randomly divided into two groups (n = 8/group): DEX-treated and saline-treated. Rats were subcutaneously administered dexamethasone sodium phosphate 0.1 mg/kg/d (Tianxin Pharmaceutical Factory, China, with a primary concentration of 5 mg/mL, diluted into 2 mg/mL with sterile saline when using, the concentration makes the injection volume half of rats weight) or saline 0.05 mL/kg/d post aurem, from gestational days 14 to 20 to conduct a prenatal DEX exposure model or vehicle control. The injection volume varied with 0.15~0.19 mL according to the weight of rats.

Pups were sexed at birth, counted, and weighed. Each litter was culled to six pups (3 male and 3 female), and all pups were kept with their birth mother until weaning on postnatal day (PD) 25. On PD50, offspring rats were randomly chosen for different experimental procedures. Health of the animals was monitored every day.

Researchers were all adequately trained and experienced in the recognition of the signs demonstrated in the experiment. Before any punctures or operations, all rats and guinea pigs (preparing for extracting antigen to induce EAE) were anesthetized with 7% chloral hydrate administrated intraperitoneally (0.5 mL/100 g), and experiment started only when animals were deeply anesthetized. After operation, animals were placed in soft and clean cages and watched carefully four times a day. If rats were found unable to access food after operation or unable to recover to stand up following experiments, they got human endpoints and were anesthetized with 20% chloral hydrate (1 mL/100 g) administrated intraperitoneally.

2.2. Assessment of EAE Susceptibility

EAE was induced in 10 DEX rats (5 male and 5 female) and 8 saline rats (4 male and 4 female) using previously established procedure with modifications [10]. Briefly, spinal cords were dissected from guinea pigs (350 ± 20 g) and homogenized on ice in an equal volume (w/v) of sterile saline with an electronic homogenizer. The spinal cords were then emulsified with an equal volume of complete Freund's adjuvant containing 6 mg/mL heat-killed Butyribacteriumrettgeri (Difco, USA) at 4°C. On PD50, EAE was induced by subcutaneous injection of the emulsion with the volume of 0.1 mL per pad at all four footpads. Rats were then intraperitoneally injected with levomisole (10 mg/kg; Sigma-Aldrich, USA) after 0, 24, and 48 h of immunization. After immunization, body weights and neurological signs of rats were examined daily. Clinical scores of these models were evaluated according to a standard scale [11]. On day 30 after immunization, rats were sacrificed, and their lumbar cords were dissected for pathological examination.

2.3. Assessment of AIA Susceptibility

AIA was induced in 11 DEX rats (6 male and 5 female) and 8 saline rats (4 male and 4 female) following a former established protocol with modifications [12]. Briefly, on PD50, 0.15 mL of complete Freund's adjuvant containing 10 mg/mL heat-killed M. tuberculosis H37Ra (Chondrex, USA) was subcutaneously injected into one hind paw pad of each rat. One week later, rats received a second injection of the same volume in the other hind paw pad. Following the second immunization, rats were visually inspected daily for signs of arthritis, erythema, and skin nodules, especially on ears and tails. Clinical scores that included arthritis and systemic symptoms were evaluated according to a standard scale by two observers blinded to group-identifying information [13]. Maximum scores, calculated at illness peak by adding the largest scores for each metric, were compared between the two groups. Twenty days after the second immunization, rats were sacrificed and their joints were removed for pathological examination. All histological procedures were performed by the histology laboratory of The First Affiliated Hospital of Sun Yat-sen University, China.

2.4. Assay of Peripheral Blood Mononuclear Cell Sensitivity to GC

The sensitivity of peripheral blood mononuclear cells (PBMCs) to GC was assessed by measuring their apoptotic response to stimulation with DEX for 18 h. The blood from 15 DEX rats (7 male and 8 female) and 12 saline rats (6 male and 6 female) on PD50 was collected by cardiac puncture into heparinized vacutioners (BD Biosciences, USA). Cardiac puncture was performed under the anesthesia of intraperitoneal injection of 7% chloral hydrate (0.5 mL/100 g). Rats died due to cardiac puncture but not euthanasia. We have described the possibility of animal death without euthanasia in our study protocol submitted to our institutional committee. PBMCs were immediately purified with Ficoll (KeyGen, China) and incubated with 10 mmol/L dexamethasone (Sigma-Aldrich, USA, dexamethasone powder dissolved with absolute alcohol) at 37°C in a humidified atmosphere containing 5% CO2. After 18 h, cells were harvested, stained with Annexin-V to detect cell membrane phosphatidylserine exposure (a measure of apoptosis), and costained with propidium iodide (PI; Keygen, China) to detect necrosis. Both Annexin-V and PI staining were measured using a flow cytometer equipped with CXP software (BD Biosciences, USA). The sensitivity of PBMCs to DEX was calculated by subtracting apoptotic PBMCs in the absence of stimulation from that observed after DEX stimulation.

2.5. GR and Methylation-Associated Genes Expression Analysis

On PD50, nine DEX rats and seven saline rats were sacrificed for molecular analysis. Total RNA was isolated from PBMCs using TRIzol (Takara, Japan) according to the manufacturer's instructions. The purity and concentration of RNA were determined using an ultraviolet spectrophotometer (NanoDrop, USA). cDNA was synthesized from total RNA using a PrimeScript RT reagent kit (Takara, Japan) and amplified by polymerase chain reaction (PCR) in duplicate (total reaction volume, 25 μL) using SYBR Premix Ex Taq (Takara, Japan) and an Applied Biosystems 7500 PCR-Cycler. The cycle parameters were 95°C for 30 s followed by 40 cycles of 95°C for 5 s, 59°C for 30 s, and 72°C for 34 s, with a final one cycle of 95°C for 15 s, 60°C for 1 min, and 95°C for 15 s. The correct size of each PCR product was verified based on a single peak in melting curves. The relative number of target genes mRNA copies was calculated with the 2−ΔΔCt method using SDS software 2.0 (Applied Biosystems, USA). The housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal control. Primers for GR, GAPDH, and methylation-associated genes designed by Primer 5 software were shown in Table 1.

Table 1

Primers for real-time PCR.

Name	Accession number	Primers	Length (bp)	T m (°C)
GAPDH	NM_017008	F: 5′-GCAAGTTCAACGGCACAG-3′	140	60.5
		R: 5′-GCCAGTAGACTCCACGACAT-3′
GR	AJ271870	F: 5′-ACTGCTGGAGGTGATTGA-3′	157	53.2
		R: 5′-GTTTCTGAAGCCTGGTAT-3′
Exon 1-1	AJ271870	F: 5′-CCAGGAAGAAGTTTCAGAT-3′	195	57
		R: 5′-GATGCAGAAACCTTGACTGT-3′
Exon 1-4	AJ271870	F: 5′-CCCCAAAGCAACACCGTA-3′	186	56.1
		R: 5′-CAGGGACTTCGTCTCTACCAG-3′
Exon 1-5	AJ271870	F: 5′-GGCATGCAACTTCCTCCGAGTTTA-3′	176	55.4
		R: 5′-AGGGCGAAGATGCAGAAACCTTGA-3′
Exon 1-6	AJ271870	F: 5′-GGGCTCACATTAATATTTGC-3′	149	54.3
		R: 5′-AAACCTTGACTGTAGCTCCTC-3′
Exon 1-7	AJ271870	F: 5′-AAGAAACTCGGTTTCCCT-3′	167	55
		R: 5′-AACCTTGACTGTAGCTCCTC-3′
Exon 1-8	AJ271870	F: 5′-CGACAGTCGCCAACAGTTAAT-3′	81	57.3
		R: 5′-TGCCAGGGACTTCGTCTCTA-3′
Exon 1-9	AJ271870	F: 5′-CGGGGACGGATTCTAAT-3′	164	55.1
		R: 5′-AGATGCAGAAACCTTGACTGT-3′
Exon 1-10	AJ271870	F: 5′-CGCCGACTTGTTTATCTG-3′	200	53.7
		R: 5′-AACCTTGACTGTAGCTCCTC-3′
Exon 1-11	AJ271870	F: 5′-GCCGCAGAGAACTCAACAGT-3′	152	58
		R: 5′-GACTGTAGCTCCTCCCCTCA-3′

GR protein expression was measured by Western-blot as described before [8, 14]. An equal amount of protein (25 μg/lane) from each cell lysate was separated on gel electrophoresis (sodium dodecyl sulfate-polyacrylamide electrophoresis) and electrophoretically transferred to polyvinylidenedifluoride membranes (Immobilon, Millipore, Billerica, MA, USA). The membranes were blocked overnight, followed by incubation for 1-2 h with antibody specifically recognizing GR (1 μg/mL; Cat number: sc-1004; Santa Cruz Biotechnology, Inc.) and β-actin (Sigma-Aldrich, St. Louis, MO). The GR antibody recognized N-terminus of GR-α clone. After washing, the blots were then incubated with peroxidase-coupled secondary anti-rabbit antibody (1 : 20,000; Sigma-Aldrich, St. Louis, MO). After being washed three times with TBS-T, antigen-antibody complexes on the membranes were detected by 3,3-diaminobenzidine (Sigma-Aldrich, St. Louis, MO, USA). Densitometry (IDV-integrated density value) measurement was used to compare the Western-blot data (Alpha Innotech, San Leandro, California).

2.6. DNA Methylation Analysis

Genomic DNA was extracted from blood using a DNA extraction kit (Omega, USA). DNA concentration and quality were measured using an ultraviolet spectrophotometer (NanoDrop, USA). DNA methylation of NR3C1 within exons 1-4, 1-5, 1-6, 1-7, 1-10, and 1-11, introns between exon1s, and the intron preceding exon 2 was detected using a Sequenom MassARRAY Methylation method (Boao Biotechnology Company, China) as described previously [15]. Primers were designed using Epidesigner (http://www.epidesigner.com). A 10-mer tag was added to the 5′ end of sense primers to balance the PCR reaction, and a T7 promoter sequence was added to the 5′ end of antisense primers for subsequent in vitro transcription. All primers are shown in Table 2. Because exon 1 contained a large number of CG sites (207), which created obstructions in primer design, we divided this region into the following seven parts and then tested each part separately: gene 1, corresponding to the promoter of exon 1-4; gene 2, corresponding to exon 1-4 and the promoter of exon 1-5; gene 3, corresponding to the promoter of both exon 1-5 and exon 1-6; gene 4, corresponding to the promoter of both exon 1-6 and exon 1-7; gene 5, corresponding to part of exon 1-10 and the promoter of exon 1-11; gene 6, corresponding to the promoter of both exon 1-11 and exon 2; and gene 7, corresponding to the promoter of exon 2 (Table 2).

Table 2

Primers for bisulfate-modified DNA.

Name	Primers	Length (bp)	CG sites	Region
Gene 1	F: 5′-aggaagagagTTAAGTGGTTTAAGGGTTTTTTATGA-3′	455	24	Promoter of exon 1-4
	R: 5′cagtaatacgactcactatagggagaaggctATAACCATAACTCCACCTCATNACC-3′
Gene 2	F: 5′-aggaagagagGGTNATGAGGTGGAGTTATGGTTAT-3′	395	18	Exon 1-4 and the promoter of exon 1-5
	R: 5′cagtaatacgactcactatagggagaaggctTTTCTAAATTCCAAAAACCCCTCTA-3′
Gene 3	F: 5′-aggaagagagGGGTTTTTGGAATTTAGAAAGTTGA-3′	461	27	Promoter of both exon 1-5 and exon 1-6
	R: 5′cagtaatacgactcactatagggagaaggctAATTACAAAACAAAACCCACCCT-3′
Gene 4	F: 5′-aggaagagagGGTGGGTTTTGTTTTGTAATTTTTT-3′	459	45	Promoter of both exon 1-6 and exon 1-7
	R: 5′cagtaatacgactcactatagggagaaggctATTTCTTTAANTTTCTCTTCTCCCAAA-3′
Gene 5	F: 5′-aggaagagagTTTTTTTATGGAGAAGAGGGGG-3′	428	50	Part of exon 1-10 and the promoter of exon 1-11
	R: 5′cagtaatacgactcactatagggagaaggctACAAAAACAAAACCAAAAAAAAC-3′
Gene 6	F: 5′-aggaagagagGTTTTTTTTGGTTTTGTTTTTGT-3′	448	26	Promoter of both exon 1-11 and exon 2
	R: 5′cagtaatacgactcactatagggagaaggctCCTCAAATATCAAAATAACACTTATCA-3′
Gene 7	F: 5′-aggaagagagTTGATTTTTTTGAGTGGTTTTTTTT-3′	498	18	Promoter of exon 2
	R: 5′cagtaatacgactcactatagggagaaggctAAATCCCCAATACAATAATAAACTTTTCC-3′

2.7. Statistical Analysis

All calculations were performed using SPSS V.13.0 statistical software (Chicago, IL, USA). An unpaired Student's t-test or a Mann–Whitney nonparametric test was applied when group variances were not equal or the data was not normally distributed, to compare the means of two groups. A chi-square test was employed when analyzing proportional data. Figures were drawn using Graphpad Prism 5 (GraphPad Software Inc., San Diego, CA, USA). A P value less than 0.05 was considered significant.

3. Results

3.1. Characteristics at Birth

The litter size and sex ratio did not differ between DEX group and vehicle control group (litter size 11.5 ± 0.9 versus 10.7 ± 1.2, female/male ratio 1.2 ± 0.1 versus 1.1 ± 0.2). However, compared with saline rats, DEX rats at birth were smaller in weight (6.2 ± 0.8 versus 5.31 ± 0.6 g, P < 0.01).

3.2. Susceptibility to EAE and AIA

Compared with saline rats, DEX rats with EAE showed higher morbidity (100% versus 62.5%, P < 0.05), earlier symptom onset (12 d versus 13.5 d, P < 0.05), higher clinical score at illness peak (4.75 versus 1.93, P < 0.01, Figure 1(a)), and more serious pathological results (Figures 1(b) and 1(c)). Thus, DEX rats were more susceptible to EAE than saline rats.

Figure 1

Comparative susceptibility of DEX and saline rats to EAE and AIA. (a) Comparison of clinical scores of EAE at illness peak, analyzed by Mann–Whitney nonparametric test. P < 0.01. (b) Histopathological analysis of DEX rats. Hematoxylin and eosin staining showing an enlargement of the lumbar spinal cord in ill (DEX) rats; nervous tissue with vacuole-like structures containing many infiltrating inflammatory cells encasing vessels. (c) Histopathological analysis of saline rats. The lumbar enlargement of spinal cord from normal (saline) rats (hematoxylin and eosin staining). (d) Comparison of clinical scores of fore limb arthritis between DEX and saline rats at illness peak, analyzed by Mann–Whitney nonparametric test. P < 0.01. (e) Histopathological analysis of DEX rats. Hematoxylin and eosin staining of an affected ankle (400x), showing many infiltrating inflammatory cells in the synovium, with vessel hyperplasia. (f) Histopathological analysis of a normal (saline) ankle.

Compared with saline rats, DEX rats with AIA exhibited greater morbidity, measured as forelimb arthritis and systemic inflammation, than saline rats (63.6% versus 0%, P < 0.05). The clinical score and pathologic manifestations of DEX rats at illness peak were also higher than saline rats (9.72 versus 7.0, P < 0.01, Figures 1(d), 1(e), and 1(f)). In addition, early after antigen administration, focal inflammatory responses of DEX rats were more serious than those of saline rats. Onset time and hind limb arthritic morbidity were not significantly different between the two groups. Collectively, these observations indicate that DEX rats were more susceptible to AIA than saline rats.

3.3. Sensitivity of PBMCs to GC

Stimulation of PBMCs from DEX rats and saline rats with the same concentration of dexamethasone for the same period induced less apoptosis from DEX rats (mean apoptosis rate 13.30%) than those from saline rats (mean apoptosis rate 19.15%, P < 0.01). Thus, the sensitivity of PBMCs from DEX rats to GC was lower than that from saline rats.

3.4. GR Expression in PBMCs

Since that GR gene contains 11 untranslated alternative first exons and each exon was processed by a distinct promoter, total GR mRNA and each exon 1 expression were evaluated in our study. Our data showed that total GR, exon 1-1, exon 1-6, exon 1-10, and exon 1-11 in PBMCs were lower in DEX rats than those in saline rats (Ρ < 0.05) (Figure 2(a)). These results were further confirmed by Western-blot, as GR protein deceased in DEX rats (Figures 2(b) and 2(c)).

Figure 2

GR and different exon 1 expression and methylation. (a) mRNA expression of total GR and different exon1 regions, analyzed by Mann–Whitney nonparametric test, P < 0.05. (b) GR protein expression from different groups. Lanes 1 to 3 were from DEX rats, while lanes 4 to 6 were from saline rats. (c) Relative quantitate data for GR protein expression, analyzed by Mann–Whitney nonparametric test, P < 0.01. (d) comparison of methylation levels in exon1s of GR gene between DEX and saline rats, analyzed by Mann–Whitney nonparametric test, # P < 0.05.

3.5. DNA Methylation Analysis

The methylation level of the 207 CG sites detected in 7 regions of GR exon 1 ranged from 0 to 0.83, with most being less than 0.3. At most sites, methylation levels did not significantly differ between DEX and saline rats. However, in the region of exon 1-10 and exon 1-11 promoter, DEX rats showed higher methylation levels than saline rats (Figure 2(d)).

4. Discussion

Adverse early life events and their long term effects have attracted more and more attention in recent years. And the development of epigenetic theory and technology advance our understanding in this filed [16, 17]. Maternal care, stress, and exogenous glucocorticoid administration are suggested to have a profound impact on offspring health [18, 19]. However, there are few reports about the prenatal DEX exposure and its long term programming effects of autoimmunity. Here we evaluated the immune programming influenced by prenatal DEX treatment and the role of GR in this process.

To test our hypothesis, pregnant Wistar rats were subcutaneously given DEX. Wistar rats are widely used for toxicology, teratology, experimental oncology, experimental gerontology, and cardiovascular research and immunology studies [6, 20]. Rats have a brief and accelerated childhood compared with humans and become sexually mature at about 6 weeks of age [5]. When considering the different phases of a rat's life, we chose day 50 of the rat to represent the adult stages. Two experimental autoimmune disease models, EAE and AIA, were performed with offspring rats in adulthood. Since sex hormone has some effect on the development of autoimmune disease, to avoid possible confounder from sex difference, we tried to make female/male ratio equal between groups [21]. The histopathological and clinical scores analysis between groups in EAE and AIA models also provided evidence that DEX administration during pregnancy increased susceptibility to autoimmune diseases in offspring, as DEX group showed more severe histopathological manifestation and clinical signs. During pregnancy, GC is often administrated when preterm delivery is expected, which is successful in stimulating the development of the fetal lung. However, emerging evidence indicated that prenatal GC exposure would result in prolonged adverse effects in offspring. In our animal model, the total dose of dexamethasone used in rats is much higher than that used for preterm labor owing to the continuous administration of DEX up to 7 days. Prenatal GC exposure not only existed in the treatment of preterm delivery, but also existed in mothers with various prenatal stresses [17, 22]. In order to cover different kinds of prenatal GC exposure, we extended the period of DEX administration. Our results also indicated that the incidence of autoimmune disease (both EAE and AIA) is increased in DEX group. These results were supportive with our initial hypothesis that prenatal GC exposure increases susceptibility to autoimmune diseases.

To further explore mechanisms causing such effects, firstly, we tested the apoptotic responses of PBMCs from DEX and saline rats to dexamethasone stimulation. PBMCs, which are mainly lymphocytes and monocytes, are the target cells of immune modulation and executive of immune response. Reduced GC response in PBMCs has been reported closely related with autoimmunity [23]. Our experiments showed that dexamethasone induced a smaller increase in apoptosis in PBMCs from DEX rats than in saline rats, which indicates that PBMCs from DEX rats were more resistant to DEX. GC physically binds to GR in immune organs and cells to induce apoptosis and thereby inhibit the immune reaction. Reduced glucocorticoid response is mainly caused by GR dysfunction and abnormal expression. So we analyzed GR expression in rat PBMCs. As measured by real-time PCR and Western-blot, total GR expression in DEX rat is significantly decreased, which is supportive for the differential cell response to GC. The publications of Lu and Cidlowski stated that multiple functional internal ATG translation initiation codons of GR resulted in multiple GR protein forms, and GR-B, GR-C, and GR-D can represent up to 50% of the total translated Gr alpha protein products [24, 25]. However, we only found one band in our plot owing to the antibody used in our study. The M20 clone used in this study has an epitope in the N-terminal region of the GR that is not present in GR-B, GR-C, and GR-D isoforms. With this antibody, we only detected the major isoform (GR-A) but not all isoforms of GR. The inability to detect these shorter GR isoforms is a limitation of our study. Since that GR exon1 contained a large number of CG sites and differentially regulated GR expression, different fragments of exon1 which are rich of CG sites were also measured by real-time PCR. Our results showed that exon 1-6, exon 1-10, and exon 1-11, in PBMCs, were lower in DEX rats than those in saline rats. These results indicate that the decreased GR expression may be caused by methylation of GR exon1. We examined DNA methylation of exon 1 in the NR3C1 gene, which includes 207 CG sites covering exons 1-4, 1-5, 1-6, 1-7, 1-10, and 1-11 and introns between them (filled with intensive GC preventing primer design, other sites were not detected), and compared their methylation level in PBMCs between saline and DEX rats. We found that a CG sites enriched region which spread from exon 1-10 to exon 1-11 was differentially methylated in DEX and saline rats. Although the exact mechanism of how this region controls GR transcription is still unknown, it is generally considered that hypermethylation in this region results in transcriptional inhibition [26, 27]. Currently, reports about the prolonged effects of prenatal GC on the methylation level in immune cells are few. However, prenatal stress or GC exposure may alter methylation level of GR gene on hippocampus, bucca, and kidney, which are also supportive with our results [19, 28, 29]. Our data added to the current literatures that prenatal GC exposure altered the methylation level of GR and caused decreased GR expression and GC response in PBMCs, which is closely related with the development of autoimmunity.

In conclusion, we found that prenatal GC exposure increases the susceptibility to autoimmunity, which is potentially caused by epigenetically programmed GR expression and glucocorticoid response.