http://orcid.org/0000-0003-3383-248X
Figures
Abstract
Polycomb repressive complex 2 (PRC2) is a chromatin binding complex that represses gene expression by methylating histone H3 at K27 to establish repressed chromatin domains. PRC2 can either regulate genes directly through the methyltransferase activity of its component EZH2 or indirectly by regulating other gene regulators. Gene expression analysis of glioblastoma (GBM) cells lacking EZH2 showed that PRC2 regulates hundreds of interferon-stimulated genes (ISGs). We found that PRC2 directly represses several ISGs and also indirectly activates a distinct set of ISGs. Assessment of EZH2 binding proximal to miRNAs showed that PRC2 directly represses miRNAs encoded in the chromosome 14 imprinted DLK1-DIO3 locus. We found that repression of this locus by PRC2 occurs in immortalized GBM-derived cell lines as well as in primary bulk tumors from GBM and anaplastic astrocytoma patients. Through repression of these miRNAs and several other miRNAs, PRC2 activates a set of ISGs that are targeted by these miRNAs. This PRC2-miRNA-ISG network is likely to be important in regulating gene expression programs in GBM.
Citation: Shivram H, Le SV, Iyer VR (2019) PRC2 activates interferon-stimulated genes indirectly by repressing miRNAs in glioblastoma. PLoS ONE 14(9): e0222435. https://doi.org/10.1371/journal.pone.0222435
Editor: Michael Nevels, University of St Andrews, UNITED KINGDOM
Received: July 16, 2019; Accepted: August 30, 2019; Published: September 12, 2019
Copyright: © 2019 Shivram et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: Primary sequencing data generated in this study is available at NCBI's GEO database under GSE133211 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE133211).
Funding: This work was funded in part by grants from the Cancer Prevention and Research Institute of Texas (RP120194) and the National Institutes of Health (NIH) (CA198648) to V.R.I. The MD Anderson Cancer Center-Science Park NGS Facility that provided Illumina sequencing was supported by CPRIT Core Facility Support Grant RP170002. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
PRC2 (polycomb repressive complex 2) is an epigenetic modifier complex that binds near transcription start sites and elsewhere in the genome and silences genes by methylating histone H3 at K27 through its histone methyltransferase component EZH2. [1]. PRC2 can exert its regulatory impact on gene expression by forming regulatory networks involving direct interactions with its target genes and indirect mechanisms where it regulates other gene regulators which in turn regulate secondary PRC2 targets [2].
Several cancers including GBM, are characterized by upregulation of multiple components of the PRC2 complex including EZH2 [2–4]. PRC2 could play an oncogenic role by silencing tumor-suppressive protein-coding genes or non-coding RNAs. miRNAs are a class of ~22 nucleotide (nt) short non-coding RNAs that play significant regulatory roles in disease and developmental pathways. Several miRNAs have been identified that play both oncogenic and tumor suppressive functions in multiple cancers including GBM [5, 6].
The regulatory interplay between miRNAs and PRC2 components has been shown to play important roles in several disease and developmental pathways. In several contexts, miRNAs exert their regulatory impact by repressing PRC2 components [7, 8]. Similarly, PRC2 can also transcriptionally silence miRNAs to cause indirect activation of genes [2]. Interestingly, miRNAs can also regulate gene expression by inhibiting or promoting PRC2 chromatin occupancy [9]. Such regulatory crosstalk between PRC2 and miRNAs is underexplored in GBM.
Previously, we identified a regulatory network between PRC2 and miRNAs where PRC2 represses genes by activating miRNAs [10]. Here, we investigated the role of PRC2-mediated repression of miRNAs in GBM-derived cell lines and bulk tumor samples. Through repression of miRNAs, PRC2 indirectly activates genes stimulated upon type I and type II interferon treatment (interferon-stimulated genes, ISGs). In contrast to previous findings in murine embryonic stem cells, we find that PRC2 directly represses miRNAs encoded in the DLK1- DIO3 locus in GBM [11]. In addition to indirectly activating ISGs, we found that PRC2 also directly represses a distinct set of ISGs.
Materials and methods
Cell lines, reagents and treatments
T98G cells (ATCC-CRL-1690) were grown in EMEM with 10% FBS. Cells were maintained at 37°C and 5% CO2. For type I interferon treatment, cells were treated with universal type I IFN-α (PBL Assay Science, 11200) at 100 U/ml final concentration for 6 hr. For type II interferon treatment, cells were treated with IFN-γ (Millipore Sigma SRP3058) at 30 ng/ml final concentration for 48 hr.
RNA-seq and data analysis
RNA from interferon treated T98G cells and primary tumor samples was extracted with TRIzol. Total RNA was then enriched for mRNA using NEXTFLEX Poly(A) Beads (Bioo 512980) as previously described [12]. Enriched RNAs were then chemically fragmented and libraries were constructed using the NEBNext small RNA library preparation kit (NEB E7330) as previously described [13]. Reads were aligned to the reference human genome version hg38 using Hisat2 and processed as previously described [10]. Reads were then counted using featureCounts and differential expression was calculated using DESeq2 [14, 15]. miRNA libraries were constructed from primary tumor samples derived as previously described [12]. Reads mapping to mature miRNAs (miRBase v21) were then counted using a custom script.
Luciferase reporter assays
The 3' UTR of target mRNAs spanning at least 500 bp around the predicted miRNA binding sites were cloned into the psi-CHECK2 vector (Promega, C8021) downstream of the Renilla luciferase gene (S1 Table). Mutagenesis was performed using the QuikChange II site directed Mutagenesis Kit (Agilent, 200523). WT and EZH2 -/- T98G cells were transfected with WT or mutant 3' UTR reporters and harvested after 48 hours. Luciferase activity was then measured using the Promega dual luciferase kit and calculated as the ratio of luminescence from Renilla to Firefly. To determine 3' UTR-dependent gene activation by EZH2, Renilla to Firefly ratios (R/F) of EZH2-activated genes were further normalized to R/F of PPIF, a gene not regulated by EZH2. To validate the interaction of specific miRNAs and EZH2-activated genes, cells were transfected with luciferase reporter vectors containing WT or mutant (miRNA binding site deletion) 3' UTR sequences. Luciferase signals were then calculated as the ratio of Renilla to Firefly.
Results
PRC2 regulates interferon-stimulated genes
Through analysis of RNA-seq data we previously generated from EZH2 -/- and WT T98G cells, we identified 1834 EZH2-repressed and 1519 EZH2-activated genes [10]. Gene ontology analysis of the EZH2-activated genes showed them to be enriched for genes involved in cellular response to interferon treatment (Fig 1A). To verify this finding, we identified genes stimulated in response to type I interferon (IFN-α) and type II interferon (IFN-γ) and then checked for the number of interferon-stimulated genes (ISGs) that were regulated by EZH2. We considered genes that showed greater than 0.5 log2 fold induction in response to either IFN-α or IFN-γ as ISGs for further analysis. We found that 40% of the genes regulated by EZH2 were ISGs (P = 9.6e-89). 737 ISGs showed higher expression levels in response to loss of EZH2, suggesting that they were normally repressed by EZH2 (Fig 1B). To check if these genes were directly repressed by EZH2, we determined the presence of EZH2 binding and H3K27me3 enrichment proximal to their transcription start sites by ChIP-seq (Fig 1C and 1D). We found that EZH2 occupied the chromatin around 53% of the EZH2 repressed ISGs (S2 Table). Thus, consistent with previously published data, we found several ISGs to be repressed by EZH2 through its histone methyltransferase activity [11, 16–18]. However, we also found a distinct set of 630 ISGs to be downregulated upon loss of EZH2, suggesting they were activated by EZH2. Given that EZH2 primarily functions as a repressor we reasoned that the ISGs activated by EZH2 were activated indirectly. Since many ISGs are known to be regulated by miRNAs, we hypothesized that many EZH2-activated ISGs could be activated by EZH2 via repression of miRNAs.
(A) Enrichment of biological processes among genes activated by EZH2, calculated by Enrichr. (B) Heat map showing differential expression of interferon-stimulated genes (ISGs) in response to EZH2 knockout. Genes are ranked in increasing order of log2 fold change (≤ -0.5 and ≥ 0.5 log2 fold change) in WT/EZH2 -/-. (C) Heat map showing enrichment scores (-log10 Q-value calculated using MACS2) for EZH2 and H3K27me3 for genes plotted in and ordered as in B. (D) Genome browser view of UNC5C and DSP showing tracks for EZH2 and H3K27me3 ChIP-seq normalized signal. Arrows indicate the direction of transcription.
PRC2 directly represses miRNAs at the DLK1-DIO3 locus
Using the miRNA-seq data that we previously generated from EZH2 -/- and WT cells [10], we identified 84 miRNAs that were significantly repressed by EZH2. To determine if PRC2 silenced specific clusters of miRNAs, we calculated the frequency of EZH2-repressed miRNAs per chromosome. Chromosome 14 stood out in showing the highest proportion of EZH2-repressed miRNAs (Fig 2A). 52 of the EZH2-repressed miRNAs mapped to the DLK1-DIO3 imprinted locus on chromosome 14. 67 of the 68 mature miRNAs expressed from the DLK1-DIO3 imprinted locus showed a trend towards repression by EZH2 (Fig 2B). We found enrichment for EZH2 and H3K27me3 at this locus in T98G GBM cells (Fig 2C), indicating that PRC2 directly and coordinately represses this cluster of miRNAs by methylating H3K27 across this locus.
(A) Fraction of miRNAs repressed by EZH2 per chromosome calculated as the ratio of the number of EZH2-repressed miRNAs and the number of miRNAs encoded per chromosome. (B) Heat map showing fold change of DLK1-DIO3 locus encoded miRNAs between EZH2 -/- and WT T98G cells. Line plot showing average log2 fold change is plotted to the right of the heat map. (C) Genome browser view showing EZH2 and H3K27me3 enrichment at the DLK1-DIO3 locus in the T98G GBM cell line and in primary tumors from GBM (GBM1, GBM2) and anaplastic astrocytoma (AA, AA2) patients. (D) Line plot showing fold change between normalized miRNA read counts for two sets of samples–GBM1 and GBM2, relative to AA1 and AA2, respectively. The X-axis shows miRNAs ordered by their average fold change across the two comparisons.
To determine if PRC2 also regulated the DLK1-DIO3 locus-encoded miRNAs in primary brain tumors, we analyzed H3K27me3 ChIP-seq data that we had previously generated in primary tumors [12]. We found H3K27me3 to be differentially enriched at this locus across multiple GBM (GBM1 and GBM2) and anaplastic astrocytoma (AA1 and AA2) tumor samples (Fig 2C). To test if low enrichment of H3K27me3 could lead to upregulation of miRNAs encoded in this locus, we performed miRNA-seq for GBM1, GBM2, AA1 and AA2. Relative to AA1 and AA2, both GBM1 and GBM2 showed lower expression for 94% and 67% of the DLK1-DIO3 locus-encoded miRNAs, respectively (Fig 2D).
PRC2 activates genes through repression of miRNAs
To determine the role of PRC2-mediated repression of miRNAs encoded in the DLK1-DIO3 locus, we identified genes that were activated by EZH2 via repression of miRNAs. The canonical targeting mechanism of miRNAs involves miRNA-mRNA pairing at the 3' UTR of the target mRNAs. If EZH2 activated genes by inhibiting expression of miRNAs that target them, we would expect the activation by PRC2 to be dependent on the 3' UTR. To test this, we performed luciferase reporter assays where we cloned the 3' UTR region of EZH2-activated genes and compared luciferase expression in EZH2 -/- and WT cells. If the activation by EZH2 was 3' UTR-dependent, we would expect the luciferase signal to be lower in EZH2 -/- cells due to higher expression of EZH2- repressed miRNAs compared to WT cells (Fig 3A). Eight genes, including 7 ISGs, out of the 17 EZH2-activated genes that we tested showed a dependence on their 3' UTR for EZH2 activation. This data suggests that these genes are under greater repression by miRNAs in EZH2 -/- cells compared to WT cells (Fig 3B).
(A) Schematic of the luciferase reporter assay to identify genes regulated by EZH2 in a 3' UTR-dependent manner. (B) Column plots showing normalized luciferase signal from WT and EZH2 -/- cells overexpressing luciferase reporter constructs containing 3' UTR regions of genes indicated on the X-axis.
We then verified the interaction between three EZH2-activated mRNA and EZH2-repressed miRNA pairs—HMGA2: miR-204-5p, MX1: miR-204-5p and IL15RA: miR-370-3p, in two ways. First, we performed luciferase reporter assays where we transfected EZH2 -/- cells with luciferase constructs containing WT or mutant 3' UTRs that lacked the miRNA binding site. Deletion of miRNA-binding sites led to a rescue in luciferase expression (Fig 4A). Second, we tested if the miRNA binding site deletion abolished the effect of the loss of miRNAs. We utilized cells stably overexpressing VP55, a vaccinia virus-derived poly(A) polymerase, to globally deplete miRNAs [19]. IL15RA and MX1 WT 3' UTR showed a rescue in luciferase signal upon miRNA depletion but the constructs with the mutant 3' UTR did not (Fig 4B). Taken together, this data shows that IL15RA is directly targeted by miR-370-3p and MX1 and HMGA2 are directly targeted by miR-204-5p. Of these miRNAs, miR-370-3p maps to the DLK1-DIO3 imprinted locus on chromosome 14.
(A) Column plots showing normalized luciferase signal from EZH2 -/- cells expressing luciferase reporter constructs containing WT and mutant 3' UTR of IL15RA, MX1 and HMGA2 lacking the binding site for their targeting miRNAs. (B) Plotted as in (A) for WT and mutant 3' UTR-transfected EZH2 -/- cells depleted of miRNAs by VP55 overexpression (+ Dox) or controls (-Dox).
To identify additional genes that could be indirectly activated by EZH2 via repression of miRNAs, we analyzed RNA-seq datasets from EZH2 -/- and WT cells overexpressing VP55 [10]. We calculated the response to miRNA loss as the ratio of gene expression in cells depleted of miRNAs (+ Dox) to their expression in control cells (no Dox). If EZH2 activated genes by repressing miRNAs, we would expect those genes to be derepressed to a greater extent in response to loss of miRNAs in EZH2 -/- cells than in WT cells. With this approach we identified 213 EZH2-activated genes including 113 ISGs, that showed greater derepression, or higher induction in expression in EZH2 -/- compared to WT cells. This suggests that these genes are more strongly repressed by miRNAs in EZH2 -/- cells than in WT cells, likely due to higher expression of miRNAs in EZH2 -/- cells (Fig 5).
Heat map showing response to loss of miRNAs calculated as the ratio of normalized read counts in VP55 overexpressing (+ Dox) and control cells (- Dox), in EZH2 -/- and WT T98G cells. ISGs are indicated by black bars. The plot to the right shows the ratio of the response to miRNA loss in EZH2 -/- to WT cells. Genes indicated in red columns show lower response and genes indicated by green columns show higher response to miRNA loss in EZH2 -/- compared to WT cells.
Discussion
The role of PRC2 in gene regulation has largely been focused on its direct repression activity [4, 20, 21]. However, PRC2 could also activate genes indirectly by repressing other repressors. We found that EZH2 directly represses a set of interferon-stimulated genes (ISGs) and also activates a distinct set of ISGs through repression of miRNAs. We first identified miRNAs that were repressed by EZH2 and then determined the EZH2-activated genes that these miRNAs directly regulate. For several other EZH2-activated genes we determined regulation by miRNAs using a new approach that involves global degradation of miRNAs by overexpression of a vaccinia virus-derived poly(A) polymerase, VP55 [19].
DLK1-DIO3 locus is a maternally imprinted locus that when deregulated can lead to several developmental disorders including cancer [11, 22]. This locus, on chromosome 14 in the human genome and chromosome 12 in the mouse, encodes several protein-coding and non-coding genes including snoRNAs and miRNAs. Nine of the miRNAs encoded in this locus have been shown to be downregulated in GBM tumors compared to normal brain samples [22]. We found that 67 of the 68 miRNAs encoded in this region to be directly silenced by PRC2 in primary tumors, consistent with the fact that EZH2 is frequently activated in GBM. This is in contrast to mouse embryonic stem cells where PRC2 was shown to activate this locus by blocking DNA methylation [11].
To further investigate the function of the PRC2-mediated repression of miRNAs encoded by the DLK1-DIO3 locus, we focused on the genes that were activated by EZH2. EZH2-activated genes showed enrichment for type I and II interferon responsive genes. We identified genes induced upon type I and II interferon treatment in GBM cells and found that 40% of the EZH2-activated genes were induced by interferon treatment. Consistent with previously published data, a significant fraction (~40%) of the EZH2-repressed genes also showed induction on interferon treatment [23]. This suggests that EZH2 represses ISGs directly and also activates them indirectly by repressing miRNAs that target them. Thus, the two key findings of this study are first, that PRC2 directly represses miRNAs encoded by the DLK1-DIO3 locus in humans, and second, that through repression of miRNAs, PRC2 activates ISGs in GBM.
Polycomb proteins and miRNAs have independently been implicated in the regulation of interferon genes [19, 23, 24]. miRNAs have also been shown to activate expression of interferon pathway genes via repression of polycomb proteins [24]. Conversely, we show that PRC2 activates ISGs by repressing miRNAs that target them, suggesting that PRC2 can both activate as well as repress expression of ISGs. A similar dual regulatory impact on interferon pathway genes was previously observed for SETD2 histone methyltransferase [25]. It is possible that, like SETD2, PRC2 also functions in the fine-tuning of interferon response. Loss of PRC2’s histone methyltransferase activity leading to differential regulation of ISGs was previously shown to suppress viral infection [26]. How PRC2-mediated regulation of interferon genes relates to antiviral activity will require further investigation.
Supporting information
S1 Table. Primers for the amplification of 3' UTRs of EZH2-activated genes.
https://doi.org/10.1371/journal.pone.0222435.s001
(PDF)
S2 Table. EZH2-activated and repressed genes.
Expression values for significantly affected genes (FDR P-value ≤ 0.05) are shown for WT and EZH2 -/- cells in two replicate experiments. Columns on the right indicate whether a gene was repressed by EZH2 and if it showed binding by EZH2 in ChIP-seq data.
https://doi.org/10.1371/journal.pone.0222435.s002
(TXT)
Acknowledgments
We thank Anna Battenhouse for assistance with aligning sequencing data and the Genomic Sequencing and Analysis Facility at UT Austin and the MD Anderson Cancer Center-Science Park NGS Facility for Illumina sequencing. The Science Park NGS Facility was supported by CPRIT Core Facility Support Grant RP170002. We also thank the Biomedical Research Computing Facility (BRCF) and the Texas Advanced Computing Center (TACC) at UT Austin for the use of computational facilities.
References
- 1. Margueron R, Reinberg D. The Polycomb complex PRC2 and its mark in life. Nature. 2011;469(7330):343–9. Epub 2011/01/21. pmid:21248841
- 2. Cao Q, Mani RS, Ateeq B, Dhanasekaran SM, Asangani IA, Prensner JR, et al. Coordinated regulation of polycomb group complexes through microRNAs in cancer. Cancer Cell. 2011;20(2):187–99. Epub 2011/08/16. pmid:21840484
- 3. Comet I, Riising EM, Leblanc B, Helin K. Maintaining cell identity: PRC2-mediated regulation of transcription and cancer. Nat Rev Cancer. 2016;16(12):803–10. Epub 2016/11/04. pmid:27658528
- 4. Conway E, Healy E, Bracken AP. PRC2 mediated H3K27 methylations in cellular identity and cancer. Curr Opin Cell Biol. 2015;37:42–8. Epub 2015/10/27. pmid:26497635
- 5. Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6(11):857–66. Epub 2006/10/25. pmid:17060945
- 6. Polioudakis D, Abell NS, Iyer VR. MiR-191 Regulates Primary Human Fibroblast Proliferation and Directly Targets Multiple Oncogenes. PLoS One. 2015;10(5):e0126535. Epub 2015/05/21. pmid:25992613
- 7. Friedman JM, Liang G, Liu CC, Wolff EM, Tsai YC, Ye W, et al. The putative tumor suppressor microRNA-101 modulates the cancer epigenome by repressing the polycomb group protein EZH2. Cancer Res. 2009;69(6):2623–9. Epub 2009/03/05. pmid:19258506
- 8. Juan AH, Kumar RM, Marx JG, Young RA, Sartorelli V. Mir-214-dependent regulation of the polycomb protein Ezh2 in skeletal muscle and embryonic stem cells. Mol Cell. 2009;36(1):61–74. Epub 2009/10/13. pmid:19818710
- 9. Graham B, Marcais A, Dharmalingam G, Carroll T, Kanellopoulou C, Graumann J, et al. MicroRNAs of the miR-290-295 Family Maintain Bivalency in Mouse Embryonic Stem Cells. Stem Cell Reports. 2016;6(5):635–42. Epub 2016/05/07. pmid:27150236
- 10. Shivram H, Le SV, Iyer VR. MicroRNAs reinforce repression of PRC2 transcriptional targets independently and through a feed-forward regulatory network. Genome Res. 2019;29(2):184–92. Epub 2019/01/18. pmid:30651280
- 11. Das PP, Hendrix DA, Apostolou E, Buchner AH, Canver MC, Beyaz S, et al. PRC2 Is Required to Maintain Expression of the Maternal Gtl2-Rian-Mirg Locus by Preventing De Novo DNA Methylation in Mouse Embryonic Stem Cells. Cell Rep. 2015;12(9):1456–70. Epub 2015/08/25. pmid:26299972
- 12. Hall AW, Battenhouse AM, Shivram H, Morris AR, Cowperthwaite MC, Shpak M, et al. Bivalent Chromatin Domains in Glioblastoma Reveal a Subtype-Specific Signature of Glioma Stem Cells. Cancer Res. 2018;78(10):2463–74. Epub 2018/03/20. pmid:29549165
- 13. Shivram H, Iyer VR. Identification and removal of sequencing artifacts produced by mispriming during reverse transcription in multiple RNA-seq technologies. RNA. 2018;24(9):1266–74. Epub 2018/06/29. pmid:29950518
- 14. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30(7):923–30. Epub 2013/11/15. pmid:24227677
- 15. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550. Epub 2014/12/18. pmid:25516281
- 16. Abou El Hassan M, Yu T, Song L, Bremner R. Polycomb Repressive Complex 2 Confers BRG1 Dependency on the CIITA Locus. J Immunol. 2015;194(10):5007–13. Epub 2015/04/12. pmid:25862816
- 17. Boyd NH, Morgan JE, Greer SF. Polycomb recruitment at the Class II transactivator gene. Mol Immunol. 2015;67(2 Pt B):482–91. Epub 2015/08/19.
- 18. Jin Y, Huo B, Fu X, Hao T, Zhang Y, Guo Y, et al. LSD1 collaborates with EZH2 to regulate expression of interferon-stimulated genes. Biomed Pharmacother. 2017;88:728–37. Epub 2017/02/06. pmid:28152483
- 19. Aguado LC, Schmid S, Sachs D, Shim JV, Lim JK, tenOever BR. microRNA Function Is Limited to Cytokine Control in the Acute Response to Virus Infection. Cell Host Microbe. 2015;18(6):714–22. Epub 2015/12/15. pmid:26651947
- 20. Holoch D, Margueron R. Mechanisms Regulating PRC2 Recruitment and Enzymatic Activity. Trends Biochem Sci. 2017;42(7):531–42. Epub 2017/05/10. pmid:28483375
- 21. Natsume A, Ito M, Katsushima K, Ohka F, Hatanaka A, Shinjo K, et al. Chromatin regulator PRC2 is a key regulator of epigenetic plasticity in glioblastoma. Cancer Res. 2013;73(14):4559–70. Epub 2013/05/31. pmid:23720055
- 22. Benetatos L, Hatzimichael E, Londin E, Vartholomatos G, Loher P, Rigoutsos I, et al. The microRNAs within the DLK1-DIO3 genomic region: involvement in disease pathogenesis. Cell Mol Life Sci. 2013;70(5):795–814. Epub 2012/07/25. pmid:22825660
- 23. Abou El Hassan M, Huang K, Eswara MB, Zhao M, Song L, Yu T, et al. Cancer Cells Hijack PRC2 to Modify Multiple Cytokine Pathways. PLoS One. 2015;10(6):e0126466. Epub 2015/06/02. pmid:26030458
- 24. Escobar TM, Kanellopoulou C, Kugler DG, Kilaru G, Nguyen CK, Nagarajan V, et al. miR-155 activates cytokine gene expression in Th17 cells by regulating the DNA-binding protein Jarid2 to relieve polycomb-mediated repression. Immunity. 2014;40(6):865–79. Epub 2014/05/27. pmid:24856900
- 25. Yanai H, Taniguchi T. Fine-tuning type I IFN signaling: A new chapter in the IFN saga. Cell Res. 2017;27(12):1407–8. Epub 2017/09/20. pmid:28925385
- 26. Arbuckle JH, Gardina PJ, Gordon DN, Hickman HD, Yewdell JW, Pierson TC, et al. Inhibitors of the Histone Methyltransferases EZH2/1 Induce a Potent Antiviral State and Suppress Infection by Diverse Viral Pathogens. MBio. 2017;8(4). Epub 2017/08/16. pmid:28811345