SUMO modification of a heterochromatin histone demethylase JMJD2A enables viral gene transactivation and viral replication

<div><p>Small ubiquitin-like modifier (SUMO) modification of chromatin has profound effects on transcription regulation. By using Kaposi’s sarcoma associated herpesvirus (KSHV) as a model, we recently demonstrated that epigenetic modification of viral chromatin by SUMO-2/3 is involved in regulating gene expression and viral reactivation. However, how this modification orchestrates transcription reprogramming through targeting histone modifying enzymes remains largely unknown. Here we show that JMJD2A, the first identified Jumonji C domain-containing histone demethylase, is the histone demethylase responsible for SUMO-2/3 enrichment on the KSHV genome during viral reactivation. Using <i>in vitro</i> and <i>in viv</i>o SUMOylation assays, we found that JMJD2A is SUMOylated on lysine 471 by KSHV K-bZIP, a viral SUMO-2/3-specific E3 ligase, in a SUMO-interacting motif (SIM)-dependent manner. SUMOylation is required for stabilizing chromatin association and gene transactivation by JMJD2A. These finding suggest that SUMO-2/3 modification plays an essential role in the epigenetic regulatory function of JMJD2A. Consistently, hierarchical clustering analysis of RNA-seq data showed that a SUMO-deficient mutant of JMJD2A was more closely related to JMJD2A knockdown than to wild-type. Our previous report demonstrated that JMJD2A coated and maintained the “ready to activate” status of the viral genome. Consistent with our previous report, a SUMO-deficient mutant of JMJD2A reduced viral gene expression and virion production. Importantly, JMJD2A has been implicated as an oncogene in various cancers by regulating proliferation. We therefore further analyzed the role of SUMO modification of JMJD2A in regulating cell proliferation. Interestingly, the SUMO-deficient mutant of JMJD2A failed to rescue the proliferation defect of JMJD2A knockdown cells. Emerging specific inhibitors of JMJD2A have been generated for evaluation in cancer studies. Our results revealed that SUMO conjugation mediates an epigenetic regulatory function of JMJD2A and suggests that inhibiting JMJD2A SUMOylation may be a novel avenue for anti-cancer therapy.</p></div

JMJD2A SUMOylation is crucial for SUMO-2/3 enrichment on the KSHV genome during viral reactivation.

<p>(A) TCLs from TREx-MH-K-Rta-shJMJD2A-Flag-JMJD2A-WT and -K471R BCBL-1 cells before and after 12 hrs Dox (0.2 μg/ml) treatment were immunoblotted with antibodies as indicated. Ratio is the relative signal of K-Rta or K-bZIP to α-Tubulin observed for Dox treatment at 24 hrs using TREx-MH-K-Rta-shJMJD2A-Flag-JMJD2A-WT BCBL-1 cells as 1.0. (B) ChIP-seq assay was performed with chromatin from cells treated as described in (A) using rabbit IgG and anti-SUMO-2/3 antibodies. Sequenced reads mapped to KSHV genome were normalized with total reads from the human genome. The histogram shows reads per million mapped reads (RPM) mapped across the KHSV genome (JMJD2A-WT, blue; -K471R, red). p = 2.27e-103 by <i>Student’s-t</i> test. The blue rectangles are high SUMO-2/3 regions as in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1006216#ppat.1006216.g001" target="_blank">Fig 1A</a>. (C) ChIP-seq data was verified by real-time qPCR using primer pairs specific for promoter regions of KSHV <i>K6</i>, <i>PAN</i>, <i>K-bZIP</i>, <i>Orf52</i>, <i>Orf23</i> and <i>Orf25</i>. (Data represent mean±SEM. n = 3. **<i>p</i><0.01. ***<i>p</i><0.005).</p

JMJD2A is required for efficient SUMO-2/3 enrichment on the viral genome during KSHV reactivation.

<p>(A) ChIP-on-chip analysis of JMJD2A binding across the KSHV lytic genome. A JMJD2A ChIP assay was performed on TREx-MH-K-Rta BCBL-1 cells treated with 0.2 μg/ml doxycycline (Dox) for 12 hours (hrs). The genomic locations of KSHV ORFs are depicted below the histogram. The blue rectangular areas are locations that comprise high SUMO-2/3 levels after KSHV reactivation (Yang et al. 2015). (B) Total cell lysates (TCLs) from TREx-MH-K-Rta-shCtrl and -shJMJD2A BCBL-1 cells before and after Dox treatment (12 hrs) were immunoblotted with antibodies as indicated. Ratio is the relative signal of K-Rta or K-bZIP to GAPDH observed for Dox treatment at 12 hrs using TREx-MH-K-Rta-shCtrl BCBL-1 cells set as 1.0. (C) ChIP was performed with chromatin prepared from cells treated as described in (B) using rabbit IgG, anti-SUMO-2/3 (upper panel) and anti-JMJD2A (lower panel) antibodies. ChIP DNA was quantified by real-time quantitative PCR (qPCR) using primer pairs specific for promoter regions of KSHV <i>K6</i>, <i>PAN</i>, <i>K-bZIP</i>, <i>Orf52</i>, <i>Orf23</i> and <i>Orf25</i>. (Data represent mean±SEM. n = 3. **<i>p</i><0.01. ***<i>p</i><0.005).</p

JMJD2A is SUMOylated at K471.

<p>(A) Schematic representation of the domain structure of JMJD2A and the three putative SUMOylation sites. (B) <i>in vitro</i> SUMOylation assay using recombinant Flag-JMJD2A and SUMO isoforms as substrates. (C) <i>in vivo</i> SUMOylation assay was performed by transfecting 293T cells with plasmids expressing Flag-JMJD2A (0.4 μg) and T7-SUMO-1 or T7-SUMO-2 and -3 (1.2 μg). JMJD2A was immunoprecipitated (IP’d) with anti-Flag-M2 beads and analyzed by immunoblotting using anti-SUMO-1, anti-SUMO-2/3 and anti-JMJD2A antibodies (left panel). TCLs were immunoblotted with anti-SUMO-1, anti-SUMO-2/3 antibodies and anti-α-Tubulin (right panel). (D) Purified wild-type (WT) and different combinations of SUMOylation site mutants of Flag-JMJD2A were used for <i>in vitro</i> SUMOylation assay as described in (B). 4%-20% gradient SDS-PAGE was used to resolve JMJD2A and SUMOylated JMJD2A. (E) Flag-JMJD2A-WT or -K471R was expressed in 293T cells as described in (C) and IP’d using M2 beads. SUMOylated JMJD2A was analyzed with anti-JMJD2A antibody.</p

KSHV SUMO E3 ligase K-bZIP catalyses SUMO modification of JMJD2A.

<p>(A) <i>in vitro</i> SUMOylation assay of JMJD2A was performed with the indicated combination of SUMO isoforms and K-bZIP. JMJD2A and SUMOylated JMJD2A were resolved by 4%-20% gradient SDS-PAGE and immunoblotted using anti-SUMO-1, anti-SUMO-2/3 and anti-JMJD2A antibodies. Input of JMJD2A and K-bZIP was detected by immunoblotting using specific antibodies as indicated. (B) <i>in vivo</i> SUMOylation assay was performed as described in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1006216#ppat.1006216.g002" target="_blank">Fig 2C</a>, with co-transfection of Flag-JMJD2A (0.4 μg), T7-SUMO-1 or T7-SUMO-2 and -3 (0.1 μg), and WT K-bZIP (0.1 μg). 4%-20% gradient SDS-PAGE was used to resolve JMJD2A and SUMOylated JMJD2A. Immunoblotting was performed using anti-JMJD2A, anti-K-bZIP and anti-GAPDH antibodies. (C) Flag-K-bZIP-WT or -L75A were co-expressed in 293T cells with Flag-JMJD2A and T7-SUMO-2 and -3 as described in (B). SUMOylated JMJD2A and input proteins were detected by immunoblotting using specific antibodies as indicated. Ratio is the relative signal of JMJD2A-SUMO to total JMJD2A observed using signal in lane 1 set as 1.0.</p

JMJD2A SUMOylation affects the expression of cellular genes involved in cancer.

<p>(A and B) WT JMJD2A but not its K471R mutant rescue proliferation of JMJD2A knockdown SLK (A) and BCBL-1 (B) cells. (A) SLK cells were sequentially infected with lentivirus expressing JMJD2A shRNA and Flag-tagged WT or K471R of JMJD2A. Cell proliferation was assessed by MTT assay. (B) Proliferation of TREx-MH-K-Rta-shJMJD2A, shJMJD2A-Flag-JMJD2A-WT and -K471R BCBL-1 cell lines was assessed by MTT assay. (C) Gene function analysis of cellular genes less upregulated in SUMOylation-deficient JMJD2A mutant (K471R) compared to JMJD2A-WT 24 hrs after KSHV reactivation. Gene list was shown in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1006216#ppat.1006216.s014" target="_blank">S7 Table</a>. (D) Genes present in more than five identified functional pathways in (C). The plus (+) indicates genes with JMJD2A binding on the promoter region (transcription start site (TSS) ± 500bp). Gene list was shown in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1006216#ppat.1006216.s015" target="_blank">S8 Table</a>. (E) RT-qPCR verification of <i>TBX3</i> expression in TREx-MH-K-Rta-shJMJD2A, -shJMJD2A-Flag-JMJD2A-WT and -K471R BCBL-1 cells. (F) qPCR was performed with ChIP DNA from <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1006216#ppat.1006216.g006" target="_blank">Fig 6A</a> using primer pairs specific for promoter region of <i>TBX3</i> (G) Schematic diagram illustrates the SUMO modification of JMJD2A in regulation gene transcription during KSHV reactivation.</p

SUMOylation modulates the occupancy and demethylase activity of JMJD2A at target genes.

<p>(A and B) ChIP was performed with chromatin prepared from TREx-MH-K-Rta-shJMJD2A-Flag-JMJD2A-WT and -K471R BCBL-1 cells using rabbit IgG, anti-JMJD2A (A) and anti-H3K9me3 (B) antibodies. ChIP DNA was quantified as described in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1006216#ppat.1006216.g001" target="_blank">Fig 1C</a>. (C) Histone demethylase activity of Flag-tagged JMJD2A-WT, -K471R or -H188A protein in 293T cells was assessed by immunofluorescence staining (IF) using deconvolution fluorescence microscopy. Cells were fixed, stained with antibody specific for Flag (FITC) and H3K9me3 (TRITC), and mounted in Slow Fade Gold with DAPI. White arrows indicate cells transfected with Flag-JMJD2A-WT or its mutant.</p

Correlation between SUMO paralog occupancy and histone marks in KSHV genome.

<p>Note: correlation coefficients were calculated according to Pearson.</p><p>Abbreviation: H3: Histone 3, H3K9me3: histone H3K9 trimethylation, H3K27me3: histone H3K27 trimethylation, AcH3: acetylated H3, H3K4me3: histone H3K4 trimethylation.</p><p>Correlation between SUMO paralog occupancy and histone marks in KSHV genome.</p

Knockdown SUMO-2/3 enhances KSHV virus production.

<p>(A) TCLs from non-induced (0 hour) and 0.2 μg/ml Dox-treated (48 hours) TREx-F3H3-K-Rta BCBL-1 and TREx-F3H3-K-Rta-shSUMO-2/3 BCBL-1 cells were analyzed by immunoblotting using antibodies as described in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1005051#ppat.1005051.g004" target="_blank">Fig 4A</a>. Ratio for each cell line is the SUMO-2/3/GAPDH signal observed for Dox treatment at 0 (for shSUMO-2/3) and 48 hour using TREx-F3H3-K-Rta BCBL-1 cells at 0 hour set as 1.0. (B) Supernatants from TREx-F3H3-K-Rta BCBL-1 and SUMO-2/3 knockdown TREx-F3H3-K-Rta-shSUMO-2/3 BCBL-1 cells were collected and filtered at 0 and 48 hours after Dox (0.2 μg/ml) treatment. Virion-associated DNA was purified and KSHV DNA levels were determined by TaqMan qPCR. Mean ± SD. Fold was determined by KSHV DNA copy number/μl at 48 hour divided by KSHV DNA copy number/μl at 0 hour. ctrl.; control, TREx-F3H3-K-Rta BCBL-1. *; P<0.05.</p

Deposition of SUMO-2/3 modification on selected KSHV promoters representing high SUMO-2/3 enrichment or high H3K9me3 mark regions during reactivation.

<p>ChIP was performed as described in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1005051#ppat.1005051.g001" target="_blank">Fig 1</a> using anti-SUMO-1 or anti-SUMO-2/3 antibodies. SUMO-1 and SUMO-2/3 binding to orf46, K-bZIP, K8.1, and orf52 promoters in SUMO-2/3 enrichment region and orf19, orf20, orf23, and orf25 promoters in H3K9me3-rich region were analyzed by real-time qPCR. An enlarged view of the corresponding ChIP-seq mapping presented in <a href="http://www.plospathogens.org/article/info:doi/10.1371/journal.ppat.1005051#ppat.1005051.g001" target="_blank">Fig 1</a> is shown above each ChIP-qPCR plot. Gene promoters are solid boxes and the direction of transcription is shown by arrows. Rabbit IgG was used as negative antibody control and enrichment is not visible in the ChIP-qPCR plots. ***; P<0.001. NS; non-significant.</p