Skip to main content

Quantification of transcriptome changes to investigate the role of glucocorticoid receptor-RNA binding during dexamethasone treatment



The glucocorticoid receptor (GR) is a well-studied, ligand-activated transcription factor and a common target of anti-inflammatory treatments. Recently, several studies have drawn attention the effects of binding of GR to RNA rather than DNA and the potential implications of this activity for GR function. The objective of our study was to further characterize the relationship between GR function and RNA binding by measuring changes in the glucocorticoid-driven transcriptome in the presence of a GR mutant that exhibited reduced RNA affinity.

Data description

GR was activated in three cell lines containing GR constructs (GR-HaloTag). One of the cell lines contained a wild-type GR-HaloTag. Another contained GR-HaloTag with a mutation that reduced RNA affinity and slightly reduced DNA affinity. The third cell line contained GR-HaloTag with a mutation that only slightly reduced DNA affinity. All three cell lines were treated with dexamethasone, a GR agonist. RNA-seq samples were collected every hour for 3 h. Moreover, transcriptome quantification was accomplished via labeling of RNAs transcribed in the final hour of dexamethasone treatment using 4-thiouridine. These labeled RNAs were then purified and sequenced. This data set is the first of its kind for GR and contains valuable insights into the function of RNA binding by GR.

Peer Review reports


The glucocorticoid receptor (GR) is a ubiquitous transcription factor regulating many genes associated with stress response, inflammation, and apoptosis [1,2,3,4,5]. In medicine, its anti-inflammatory properties are commonly induced using the synthetic glucocorticoid dexamethasone [6]. Because of its important role in controlling inflammation, the regulation of GR has been studied in great detail. In the absence of glucocorticoids, GR is mostly found in the cytoplasm in its inactive form. Once glucocorticoids diffuse into the cell, they can bind to GR and initiate rapid nuclear translocation and subsequent regulation of downstream genes [3, 6,7,8]. Direct gene activation by GR typically happens through dimerization on sequence-specific DNA elements, and thus most studies of GR-mediated gene regulation have focused on this interaction [9,10,11,12,13]. However, several groups identified an additional, unexpected biochemical activity of direct RNA binding by GR to a common underlying motif of hairpin RNAs [2, 14,15,16]. Further, in vitro studies revealed structure-specific binding to RNA that is competitive with DNA [17]. The implications of these findings are poised to change how we think about both GR biology and, to a greater extent, the regulation of many transcription factors, by adding another layer to GR’s complex regulatory pathway. The objective of this data set is to better define the relationship between GR-RNA binding and GR transcriptional action through the quantification of the dexamethasone-driven transcriptome in the context of wild-type (wt) and RNA-binding deficient GR.

Data description

To measure the effect RNA binding by GR has on the GR-regulated transcriptome, we performed RNA sequencing of 4-thiouridine-labeled RNAs (4sU-seq, Data set 1, Table 1) after dexamethasone treatment in cells expressing a separation-of-function GR mutant. We employed U2OS cells stably expressing wt GR-HaloTag or a K492A mutant of GR that has an 11-fold reduction in RNA affinity based on in vitro binding assays [17]. Due to a moderate reduction in DNA affinity in this mutant, we also used cells expressing an R470A GR mutant that exhibits similar reduction in DNA affinity but maintains RNA affinity. As U2OS cells express low levels of endogenous GR and show impaired glucocorticoid response, we can compare transcriptional responses to dexamethasone between wt GR and our mutant to attribute the function of GR-RNA binding to known regulatory pathways. After stable integration of the GR-HaloTag constructs, cells were sorted using FACS and sorted cell fractions were selected to match GR-HaloTag expression between cell lines by immunofluorescence. This sorting was done to minimize the impact of GR-HaloTag abundance on differential dexamethasone response between cell lines. After selecting cell lines with similar GR-HaloTag expression, cells were incubated with 100 nM dexamethasone for 0, 1, 2, or 3 h. Zero-hour dexamethasone cells received an equivalent volume of ethanol (0.01%). During the final hour of dexamethasone treatment, cells were given 200 µM 4sU to label recently transcribed RNAs. Zero-hour dexamethasone cells treated with 200 µM 4sU for 1 h immediately after adding ethanol. At the end of treatment, cells were lysed with TRIzol and RNAs were purified. 4sU-labeled RNAs were then biotinylated and isolated using magnetic streptavidin beads. A detailed 4sU-seq protocol can be found in Data file 1, Table 1.

Libraries for were prepared using the KAPA RNA HyperPrep kit with RiboErase. The dual index adapters included unique molecular identifiers (UMIs). Library sizes were about 330 bp and they were sequenced for 2 × 150 bp paired-end reads on an Illumina NovaSEQ 6000 by the University of Colorado School of Medicine Genomics and Microarray Core Facility. Reads were tagged with UMIs using UMI-tools (v1.1.2) and trimmed using Trim Galore (v0.6.6) with the following parameters: --2colour 20, --paired [18, 19]. Trimmed reads were then aligned to the hg38 genome assembly using STAR (v2.7.3a) and deduplicated with UMI-tools (v1.1.2) using the following parameters: --paired, --unpaired-reads discard, --chimeric-pairs discard [18, 20]. FeatureCounts (Rsubread v2.0.1, R v4.0.3) was used to count deduplicated read coverage over exons (Data file 2, Table 1) [21]. Quality control analyses were performed with FastQC on trimmed reads and RSeQC on aligned reads before deduplication and packaged together with MultiQC (Data file 3, Table 1) [22,23,24]. All scripts used for data processing, including options and flags used at each step, can be found on GitHub (Data set 2, Table 1).

Table 1 Overview of data files/data sets


To meaningfully compare transcriptomic changes between the different GR-expressing cell lines, it was crucial to first match GR expression between cell lines. We sorted cells and selected fractions that had the same nuclear signal after dexamethasone treatment based on immunofluorescence. Nuclear signal was used to match the amount of active GR between cell lines, but transcript expression of GR is higher in the mutant lines.

Data availability

The data described in this Data note can be freely and openly accessed on NCBI GEO under GSE216337 ( The detailed protocol for 4sU-seq can be accessed through Zenodo ( Quality assessment is available through Zenodo ( Data processing scripts are available through GitHub (via Zenodo, Please see Table 1 and references [25,26,27,28] for details and links to the data.



glucocorticoid receptor






unique molecular identifier


  1. McDowell IC, Barrera A, D’Ippolito AM, Vockley CM, Hong LK, Leichter SM, et al. Glucocorticoid receptor recruits to enhancers and drives activation by motif-directed binding. Genome Res. 2018;28(9):1272–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Kino T, Hurt DE, Ichijo T, Nader N, Chrousos GP. Noncoding RNA Gas5 is a growth arrest- and Starvation-Associated Repressor of the glucocorticoid receptor. Sci Signal. 2010;3(107):ra8.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Weikum ER, Knuesel MT, Ortlund EA, Yamamoto KR. Glucocorticoid receptor control of transcription: precision and plasticity via allostery. Nat Rev Mol Cell Biol. 2017;18(3):159–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Ramamoorthy S, Cidlowski JA. Ligand-Induced repression of the glucocorticoid receptor gene is mediated by an NCoR1 repression Complex formed by long-range chromatin interactions with intragenic glucocorticoid response elements. Mol Cell Biol. 2013;33(9):1711–22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Bothe M, Buschow R, Meijsing SH. Glucocorticoid signaling induces transcriptional memory and universally reversible chromatin changes. Life Sci Alliance [Internet]. 2021;4(10). Available from:

  6. Reddy TE, Pauli F, Sprouse RO, Neff NF, Newberry KM, Garabedian MJ, et al. Genomic determination of the glucocorticoid response reveals unexpected mechanisms of gene regulation. Genome Res. 2009;19(12):2163–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Htun H, Barsony J, Renyi I, Gould DL, Hager GL. Visualization of glucocorticoid receptor translocation and intranuclear organization in living cells with a green fluorescent protein chimera. Proc Natl Acad Sci. 1996;93(10):4845–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Savory JGA, Hsu B, Laquian IR, Giffin W, Reich T, Haché RJG, et al. Discrimination between NL1- and NL2-Mediated Nuclear localization of the glucocorticoid receptor. Mol Cell Biol. 1999;19(2):1025–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Luisi BF, Xu WX, Otwinowski Z, Freedman LP, Yamamoto KR, Sigler PB. Crystallographic analysis of the interaction of the glucocorticoid receptor with DNA. Nature. 1991;352(6335):497.

    Article  CAS  PubMed  Google Scholar 

  10. Meijsing SH, Pufall MA, So AY, Bates DL, Chen L, Yamamoto KR. DNA binding site sequence directs glucocorticoid receptor structure and activity. Science. 2009;324(5925):407–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Watson LC, Kuchenbecker KM, Schiller BJ, Gross JD, Pufall MA, Yamamoto KR. The glucocorticoid receptor dimer interface allosterically transmits sequence-specific DNA signals. Nat Struct Mol Biol. 2013;20(7):876–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Schiller BJ, Chodankar R, Watson LC, Stallcup MR, Yamamoto KR. Glucocorticoid receptor binds half sites as a monomer and regulates specific target genes. Genome Biol. 2014;15(8):418.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Frank F, Liu X, Ortlund EA. Glucocorticoid receptor condensates link DNA-dependent receptor dimerization and transcriptional transactivation. Proc Natl Acad Sci. 2021;118(30):e2024685118.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ali M, Vedeckis W. The glucocorticoid receptor protein binds to transfer RNA. Science. 1987;235(4787):467–70.

    Article  CAS  PubMed  Google Scholar 

  15. Ishmael FT, Fang X, Houser KR, Pearce K, Abdelmohsen K, Zhan M, et al. The human glucocorticoid receptor as an RNA-Binding protein: global analysis of glucocorticoid Receptor-Associated transcripts and identification of a target RNA motif. J Immunol. 2011;186(2):1189–98.

    Article  CAS  PubMed  Google Scholar 

  16. Hudson WH, Pickard MR, de Vera IMS, Kuiper EG, Mourtada-Maarabouni M, Conn GL et al. Conserved sequence-specific lincRNA–steroid receptor interactions drive transcriptional repression and direct cell fate. Nat Commun [Internet]. 2014;5(1). Available from:

  17. Parsonnet NV, Lammer NC, Holmes ZE, Batey RT, Wuttke DS. The glucocorticoid receptor DNA-binding domain recognizes RNA hairpin structures with high affinity. Nucleic Acids Res. 2019;47(15):8180–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Smith T, Heger A, Sudbery I. UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy. Genome Res. 2017;27(3):491–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Krueger F. Trim Galore! [Internet]. Available from:

  20. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  PubMed  Google Scholar 

  21. Liao Y, Smyth GK, Shi W. The R package rsubread is easier, faster, cheaper and better for alignment and quantification of RNA sequencing reads. Nucleic Acids Res. 2019;47(8):e47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Andrews S, FastQC. A Quality Control tool for High Throughput Sequence Data [Internet]. Available from:

  23. Wang L, Wang S, Li W. RSeQC: quality control of RNA-seq experiments. Bioinformatics. 2012;28(16):2184–5.

    Article  CAS  PubMed  Google Scholar 

  24. Ewels P, Magnusson M, Lundin S, Käller M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics. 2016;32(19):3047–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Lammer NC, Allen MA, Batey RT, Wuttke DS. Effect of reduced glucocorticoid receptor-RNA affinity on dexamethasone treatment. Gene Expression Omnibus. 2022.

  26. Lammer NC, Allen MA, Batey RT, Wuttke DS. 4sU-seq protocol. Zenodo. 2022.

  27. Lammer NC, Allen MA, Batey RT, Wuttke DS. MultiQC report for 4sU-seq data. Zenodo. 2022.

    Article  Google Scholar 

  28. Lammer NC, Allen MA. nicklammer/gr_seq_processing: reorganized scripts (v0.9.1). Zenodo.

Download references


We thank Theresa Nahreini and the Biochemistry Cell Culture Facility for supporting tissue culture and the University of Colorado School of Medicine Genomics and Microarray Core Facility for sequencing our libraries. We acknowledge the BioFrontiers Computing Center at the University of Colorado Boulder for providing High Performance Computing resources supported by BioFrontiers Information Technology. Additionally, we thank Professor John Rinn (CU-Boulder) for valuable feedback on experimental design and data processing. We also thank Katy Walsh for providing and assisting with the 4sU-seq protocol.


This work was supported by grants from the National Institutes of Health to DSW and RTB (R01GM120347), to MAA (R01HL156475), and to NCL (T32GM065103).

Author information

Authors and Affiliations



NCL, DSW, and RTB designed the research. NCL collected samples and prepared sequencing libraries. NCL and MAA processed the data. NCL drafted the manuscript with review and approval from all authors. DSW and RTB supported the project.

Corresponding authors

Correspondence to Robert T. Batey or Deborah S. Wuttke.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

MAA has a patent for “Methods for predicting transcription factor activity” that is not related to the work contained in this data note.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lammer, N.C., Allen, M.A., Batey, R.T. et al. Quantification of transcriptome changes to investigate the role of glucocorticoid receptor-RNA binding during dexamethasone treatment. BMC Res Notes 16, 181 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Transcription factor
  • Glucocorticoid receptor
  • RNA binding
  • RNA sequencing
  • Transcriptomics
  • 4-thiouridine
  • NR3C1