Nascent transcription upon interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines
BMC Research Notes volume 16, Article number: 292 (2023)
The interferon-triggered innate immune response has been observed to be under strong diversifying selection to counteract the many pathogens hosts have to defend against. In particular, rewiring of gene transcription regulation allows organisms to rapidly acquire new phenotypes by removing and adding genes into the innate immune gene network. Dissecting the molecular processes by which this rewiring takes place, either by changing the DNA regulatory elements or by changing the activity of the regulators across species, is key to better understand this evolutionary process.
To better comprehend the evolutionary dynamics that have occurred in the initial transcriptional response to interferon in primates, we present Precision Run-On (PRO-seq) datasets made after 1 h of interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines. Further, we tested the difference between using either species’ cognate interferon versus using the other orthologous interferon to account for any potential impacts in the interaction of the orthologous interferons with their cellular membrane receptors. This data provides insights into the regulatory mechanisms that drive species-specific responses to environmental perturbations, such as the one driven by the interactions of pathogens and their hosts.
Eukaryotic gene transcription is controlled by promoters and enhancers that are proximal and distal, respectively, to the genes they regulate . When these regulatory elements are actively used by the cell to control gene expression, they are transcribed, typically from both DNA strands originating from a shared RNA polymerase loading locus. These nascently transcribed molecules play diverse roles in the regulation of the genes these regulatory elements are targeting [2, 3].
Regulatory elements have been shown to evolve more rapidly than the genes they regulate , which suggest that organisms rely on rewiring gene transcription regulation to quickly acquire new phenotypes. In particular, the cell-intrinsic innate immune system controlled by the interferon (IFN) cytokines has been observed to be one of the most rapidly evolving gene networks [5,6,7].
IFN molecules are released upon pathogen recognition, and once sensed by cell membrane receptors trigger the formation of transcription factor complexes composed by members of the IRF and STAT protein families, which then regulate the transcription of interferon stimulated genes (ISGs) . However, it has remained an important question in the field how is the rewiring of gene transcription regulation established when evolutionary changes can simultaneously occur both in the regulators themselves, such as mutations in transcription factors, and in changes in the DNA regulatory sequences of the genes these regulators target . Furthermore, biophysical studies demonstrate that varying binding affinities between IFN paralogs to their receptors can ensue distinct cellular responses [10, 11], which may also affect IFN orthologs across species.
Here we address these questions by testing the difference in the transcriptional response when using species-matching compared to non-species matching IFN-α2 treatments on cells derived from human and rhesus macaque, which diverged around 25 mya . These datasets will be of value to researchers who want to investigate the evolutionary changes in the interferon-triggered innate immune transcriptional response in primates.
We obtained precision run-on sequencing (PRO-seq) datasets from Epstein-Barr Virus-transformed lymphoblastoid cell lines (LCLs) derived from both a male and a female individual from two primates: human (Homo sapiens, GM12878 and HG03077 from the Coriell Institute for Medical Research) and rhesus macaque (Macaca mulatta, Mm 150 − 99 and Mm 290 − 96 generously shared by Yoav Gilad from the University of Chicago) (see Data set 1 in Table 1). Each primate cell line was treated for 1 h with either the vehicle bovine serum albumin (BSA), with their cognate species-matching IFN-α2, or with the other primate’s orthologous IFN-α2. Libraries were prepared as described in  with minor modifications, see methods file (Data file 3) for full details. Libraries were sequenced on an Illumina NextSeq 500 using single-end sequencing, on which an average of 36 million reads per library were obtained. Raw reads (available at NCBI GEO under accession number GSE214304) were assessed for quality (see data file Data file 1) and mapped to the respective species genomes (hg38 for Homo sapiens and rheMac10 for Macaca mulatta) (see Data file 3 for full details).
Briefly, we counted reads over gene bodies using featureCounts and used DESeq2  and known ISGs to validate that the LCLs responded as expected to the IFN-α2 treatments. Additionally, we assessed the consistent simulation of key transcription factors upon IFN-α2 by first detecting bidirectionally transcribed loci using Tfit  and dREG . These bidirectionally transcribed loci, roughly corresponding to enhancers and promoters, are then fed to TFEA to identify transcription factors responding to IFN-α2 treatment  (see Data file 3 for full details).
The human IFN-α2 and rhesus IFN-α2-treated LCLs display a typical type I interferon stimulation transcriptional response compared to the BSA control datasets, and similar to published nascent transcription response to IFN-γ on mouse embryonic fibroblasts (MEFs)  (see Data file 2, Data file 4). The human IFN-α2-treated datasets, however, show a greater interferon stimulation magnitude than the rhesus IFN-α2-treated datasets, regardless of the primate LCL used.
The human IFN-α2 protein was obtained from Proteintech Cat. no. HZ-1066, whereas the rhesus macaque IFN-α2 was obtained from PBL Assay Science Ca. no. 16105-1. Each manufacturer tested their purified protein activities using different assays, with the human IFN-α2 protein purification having been tested with a “dose-dependent cytotoxicity of the human TF-1 cell line (human erythroleukemic indicator cell line)” , and the rhesus IFN-α2 protein purification with a “cytopathic inhibition assay on Bovine (MDBK) kidney cells with vesicular stomatitis [virus] (VSV)” . Discrepancies in the bioactivity assay details may have resulted in unequal magnitude of IFN-dependent transcriptional responses even when using 100 units/mL for both the human and rhesus IFN-α2 protein treatments. To this end, we observe that both cell lines responded more strongly to the human IFN-α2, as observed by the number of differentially transcribed genes and by the magnitude of the ISGs fold-change.
While we assayed two distinct cell lines per species, one female and one male, each biological sex was only assayed once.
The data described in this Data note can be freely and openly accessed on NCBI under GEO accession number GSE214304. Please see Table 1 and references Data file 1, 2, 3, & 4 and Data set 1 for details and links to the data.
Precision Run-On followed by sequencing
Bovine Serum Albumin
Interferon Stimulated Genes
Lymphoblastoid Cell Line
Mouse Embryonic Fibroblasts
Cramer P. Organization and regulation of gene transcription. Nature. 2019;573:45–54. https://doi.org/10.1038/s41586-019-1517-4.
Wissink EM, Vihervaara A, Tippens ND, et al. Nascent RNA analyses: tracking transcription and its regulation. Nat Rev Genet. 2019;20:705–23. https://doi.org/10.1038/s41576-019-0159-6.
Skalska L, Beltran-Nebot M, Ule J, et al. Regulatory feedback from nascent RNA to chromatin and transcription. Nat Rev Mol Cell Biol. 2017;18:331–7. https://doi.org/10.1038/nrm.2017.12.
Villar D, Berthelot C, Aldridge S, et al. Enhancer evolution across 20 mammalian species. Cell. 2015;160:554–66. https://doi.org/10.1016/j.cell.2015.01.006.
Danko CG, Choate LA, Marks BA, et al. Dynamic evolution of regulatory element ensembles in primate CD4 + T cells. Nat Ecol Evol. 2018;2:537–48. https://doi.org/10.1038/s41559-017-0447-5.
Meyerson NR, Sawyer SL. Two-stepping through time: mammals and viruses. Trends Microbiol. 2011;19:286–94. https://doi.org/10.1016/j.tim.2011.03.006.
Judd EN, Gilchrist AR, Meyerson NR, et al. Positive natural selection in primate genes of the type I interferon response. BMC Ecol Evo. 2021;21:65. https://doi.org/10.1186/s12862-021-01783-z.
Ivashkiv L, Donlin L. Regulation of type I interferon responses. Nat Rev Immunol. 2014;14:36–49. https://doi.org/10.1038/nri3581.
Britton CS, Sorrells TR, Johnson AD. Protein-coding changes preceded cis-regulatory gains in a newly evolved transcription circuit. Science. 2020;367:96–100. https://doi.org/10.1126/science.aax5217.
Piehler J, Thomas C, Garcia KC, et al. Structural and dynamic determinants of type I interferon receptor assembly and their functional interpretation. Immunol Rev. 2012;250:317–34. https://doi.org/10.1111/imr.12001.
Li H, Sharma N, General IJ, et al. Dynamic modulation of binding Affinity as a mechanism for regulating Interferon Signaling. J Mol Biol. 2017;429:2571–89. https://doi.org/10.1016/j.jmb.2017.06.011.
Gibbs RA, Rogers J, Katze MG, et al. Evolutionary and Biomedical Insights from the Rhesus Macaque Genome. Science. 2007;316:222–34. https://doi.org/10.1126/science.1139247.
Mahat D, Kwak H, Booth G, et al. Base-pair-resolution genome-wide mapping of active RNA polymerases using precision nuclear run-on (PRO-seq). Nat Protoc. 2016;11:1455–76. https://doi.org/10.1038/nprot.2016.086.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://doi.org/10.1186/s13059-014-0550-8.
Azofeifa JG, Allen MA, Hendrix JR, et al. Enhancer RNA profiling predicts transcription factor activity. Genome Res. 2018;28:334–44. https://doi.org/10.1101/gr.225755.117.
Wang Z, Chu T, Choate LA, et al. Identification of regulatory elements from nascent transcription using dREG. Genome Res. 2019;29:293–303. https://doi.org/10.1101/gr.238279.118.
Rubin JD, Stanley JT, Sigauke RF, et al. Transcription factor enrichment analysis (TFEA) quantifies the activity of multiple transcription factors from a single experiment. Commun Biol. 2021;4:661. https://doi.org/10.1038/s42003-021-02153-7.
Steinparzer I, Sedlyarov V, Rubin JD, et al. Transcriptional responses to IFN-γ require mediator kinase-dependent pause release and mechanistically distinct CDK8 and CDK19 functions. Mol Cell. 2019;76:485–99. https://doi.org/10.1016/j.molcel.2019.07.034.
HumanKine recombinant human IFN alpha 2A protein, Proteintech Group Inc., Rosemont IL 60018 USA, catalog # HZ-1066, Lot # 0615–01. https://www.ptglab.com/products/recombinant-human-ifn-alpha-2a.htm.
Cynomolgus IFN-Alpha. 2 (Ile 16) Mammalian, Pestka Biomedical Laboratories Inc., Piscataway NJ 08854 USA, Catalog # 16105 Lot # 6987R. https://www.pblassaysci.com/proteins/cynomolgus-ifn-alpha-2-ile-16-mammalian-16105.
Data file 1. Quality control assessment of data sets sequencing. Figshare. https://doi.org/10.6084/m9.figshare.21253245.v1.
Data file 2. Assessment of interferon stimulation of data sets. Figshare. https://doi.org/10.6084/m9.figshare.21287637.v1.
Data file 3. Detailed methods file. Figshare. https://doi.org/10.6084/m9.figshare.21287652.v1.
Data file 4. DESeq2 results tables. Figshare. https://doi.org/10.6084/m9.figshare.23971908.
Data. set 1. PRO-seq datasets of Homo sapiens and Macaca mulatta LCLs treated with either BSA, human IFN-α2, or rhesus IFN-α2. NCBI Gene Expression Omnibus. https://identifiers.org/geo:GSE214304.
We acknowledge the BioFrontiers Computing Core at the University of Colorado Boulder for high-performance computing resources support (NIH R01 GM125871-04S1) by BioFrontiers’ IT. We also acknowledge the support of the University of Colorado Boulder BioFrontiers Institute Next Generation Sequencing facility.
National Institutes of Health R01 GM125871 (Dowell) and National Institutes of Health R35 GM128822 (Chuong).
The authors declare no competing interests.
Ethics approval and consent to participate
No human participants are involved in the study.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ramirez, D., B. Chuong, E. & D. Dowell, R. Nascent transcription upon interferon-α2 stimulation on human and rhesus macaque lymphoblastoid cell lines. BMC Res Notes 16, 292 (2023). https://doi.org/10.1186/s13104-023-06465-1