Defining the NSD2 interactome: PARP1 PARylation reduces NSD2 histone methyltransferase activity and impedes chromatin binding

NSD2 is a histone methyltransferase that specifically dimethylates histone H3 lysine 36 (H3K36me2), a modification associated with gene activation. Dramatic overexpression of NSD2 in t(4;14) multiple myeloma (MM) and an activating mutation of NSD2 discovered in acute lymphoblastic leukemia are significantly associated with altered gene activation, transcription, and DNA damage repair. The partner proteins through which NSD2 may influence critical cellular processes remain poorly defined. In this study, we utilized proximity-based labeling (BioID) combined with label-free quantitative MS to identify high confidence NSD2 interacting partners in MM cells. The top 24 proteins identified were involved in maintaining chromatin structure, transcriptional regulation, RNA pre-spliceosome assembly, and DNA damage. Among these, an important DNA damage regulator, poly(ADP-ribose) polymerase 1 (PARP1), was discovered. PARP1 and NSD2 have been found to be recruited to DNA double strand breaks upon damage and H3K36me2 marks are enriched at damage sites. We demonstrate that PARP1 regulates NSD2 via PARylation upon oxidative stress. In vitro assays suggest the PARylation significantly reduces NSD2 histone methyltransferase activity. Furthermore, PARylation of NSD2 inhibits its ability to bind to nucleosomes and further get recruited at NSD2-regulated genes, suggesting PARP1 regulates NSD2 localization and H3K36me2 balance. This work provides clear evidence of cross-talk between PARylation and histone methylation and offers new directions to characterize NSD2 function in DNA damage response, transcriptional regulation, and other pathways.

NSD2 (MMSET/WHSC1) is a member of the nuclear receptor-binding SET domain (NSD) 3 family, functioning as a histone methyltransferase (HMT) that catalyzes dimethylation of histone H3 lysine 36 specifically (1,2). NSD2 was identified to be the key gene disrupted and overexpressed by fusion to the immunoglobulin heavy-chain promoter/enhancer in the t(4; 14) translocation, found in almost 15-20% of multiple myeloma (MM) cases (3,4). NSD2 overexpression in t(4;14)ϩ MM cells leads to increased cell proliferation and altered downstream gene expression due to a genome wide increase in H3K36 dimethylation and a decrease in gene repression-associated H3K27 trimethylation (H3K27me3) (5)(6)(7). Recently, a gain-of-function point mutation of NSD2 (E1099K) has been discovered in multiple leukemia cell lines and patient samples that similarly causes elevated H3K36me2 and altered gene activation (8 -12).
The H3K36me2 mark serves important roles in DNA damage repair (13)(14)(15). Upon radiation, H3K36me2 marks are enriched at double strand breaks (DSB) to recruit DNA repair components, including NBS1 and Ku70, initiating nonhomologous end joining (NHEJ) (16,17). This modification has also been found to promote the resolution of RAD51 dsDNA filaments during late stage homologous recombination (HR) (18). NSD2 was identified as a member of the DNA damage response pathway through genetic screening and is present at sites of DNA damage upon UV laser microirradiation (19,20). It serves an important role to facilitate the recruitment of 53BP1, a DNA damage response regulator, to DSBs (20). NSD2 is also involved in initiation of class switch recombination during B-cell development. Deregulation of NSD2 causes decreased activationinduced cytidine deaminase-mediated DNA breaks, which further leads to defects of class switch recombination (21,22). Our lab previously discovered that NSD2 promotes both NHEJ and HR in U2OS cells. Furthermore, MM cells overexpressing NSD2 display an enhanced rate of DNA damage repair that allows cells to survive upon chemotherapy treatment (23).
In addition to its histone methyltransferase function, NSD2 has been found to play important roles in multiple molecular pathways. For instance, NSD2 has been reported to interact with the bromodomain protein 4 (BRD4) and the positive transcription elongation factor ␤ (P-TEF␤) to facilitate transcriptional elongation of rapidly induced genes (24). NSD2 also associates with HIRA, the histone H3.3-specific histone chaperone, to promote histone H3.3 deposition along gene bodies of rapidly activated genes, suggesting a role for NSD2 in chromatin assembly (24). Moreover, NSD2 is involved in repression of microRNA miR-126 in MM through an interaction with the transcriptional corepressor TRIM28 and the histone deacetylases HDAC1/2 that results in maintaining a high protein level of the oncoprotein c-MYC (25).
An initial screen for NSD2 partner proteins identified a preliminary set of 19 interaction partners (25). However, this work was carried out in nonphysiologically relevant 293 human embryonic kidney cells and required the use of relatively harsh conditions to extract NSD2 from cells, which are tightly bound to chromatin. Thus, we have now utilized BioID proximitybased labeling (26,27), which does not require the use of harsh conditions that may disrupt weaker but biologically relevant complexes, to obtain a more complete understanding of the biological functions and mechanisms of NSD2. This label-free quantitative study minimizes nonspecific contamination and helps define the NSD2 interaction network. Using this method, we have identified with high confidence 24 nuclear interacting partners of NSD2 involved in multiple molecular pathways and biological processes including chromatin remodeling, gene expression regulation, DNA damage repair, and RNA splicing. The mapped protein network provides novel insights into the functions and mechanisms of NSD2 operative in multiple myeloma.

Study 1: identification of NSD2 interacting partners by BioID (proximity labeling and IP-MS)
NSD2 protein binds tightly to the chromatin scaffold and remains insoluble in the nuclear pellet after standard nuclear extraction. Accordingly, an enzyme shearing mixture is required to digest chromatin and release NSD2 into the soluble nuclear fraction, which might disrupt any physiological NSD2 complex. Because of this, we adopted BioID, a method in which a promiscuous biotin ligase moiety is engineered into the protein of interest, to biotinylate proteins based on proximity and screen for NSD2 interacting proteins. To identify physiologically relevant interacting proteins of NSD2, we chose the TKO NSD2-low cell line (5), which was derived from the t(4;14)ϩ KMS11 MM cell line in which the rearranged and overexpressed NSD2 allele was disrupted and only the WT NSD2 expresses. These cells were transduced with retrovirus encoding a BirA-tagged version of NSD2. The resulting TKO cell line stably expressed NSD2-BirA and mimicked high expression of NSD2 in t(4;14)ϩ MM. These TKO ϩ NSD2-BirA cells were cultured in medium supplemented with biotin for the BirA domain to biotinylate proteins that are in close proximity to NSD2. NSD2 and its biotinylated partners were extracted in SDS-containing lysis buffer, captured, and subsequently analyzed by nanocapillary LC-MS (Fig. 1A). TKO cells expressing NSD2-BirA had an elevated histone H3K36me2 level, suggesting fusion of the BirA tag did not affect NSD2 histone methyltransferase activity (Fig. 1B). Six biological replicates of NSD2-BirA BioID IP-MS were analyzed. A total of 94 proteins identified in at least five of the six biological replicates were considered high confident targets (Fig. 1C). In our study design, the same workflow was performed on TKO cells in the absence of BirA-NSD2 as a negative control for endogenous biotinylated proteins and nonspecific binders to streptavidin beads. Of the 55 proteins identified at least twice in three biological replicates of the TKO control, 31 were in common with NSD2-BirA (Fig. 1C). As a result, 63 unique NSD2 high confident interacting proteins were identified (Table S1). Gene set enrichment analysis of the 63 proteins using EnrichR (28) sug- Biotinylated proteins were captured by streptavidin beads. On-bead trypsin digestion was performed and the resulting peptides were desalted and identified by nanocapillary LC-MS. B, immunoblot of TKO cells stably expressing empty vector, NSD2, or NSD2-BirA. Because the NSD2-BirA construct contains a C-terminal HA tag, the HA tag blot confirmed expression of the NSD2-BirA construct in TKO cells. C, a Venn diagram depicting unique and common interacting proteins identified in NSD2-BirA cells and TKO control cells. D, gene set enrichment analysis of NSD2-binding partners organized according to GO biological process using EnrichR (28).

PARylation impedes NSD2 binding to chromatin
gested that NSD2 partner proteins are highly associated with chromatin remodeling and chromatin modification according to the GO biological process 2018. This validated our process and NSD2 functions as a chromatin regulator (Fig. 1D).

Study 2: label-free quantitative MS to capture NSD2 interaction partners
To further identify highly confident NSD2 interaction partners and reduce noise from BirA random labeling, we developed a BirA control to eliminate nonspecific proteins from the screen. A BirA ligase along with a C-terminal nuclear localization signal peptide was cloned into a retrovirus vector and stably transduced into the TKO cell line. The control cell line was then used to perform comparative quantification of biotinylated proteins in the NSD2-BirA experiment versus BirA alone. Three biological replicates of TKO expressing NSD2-BirA or BirA control were cultured in medium supplemented with biotin ( Fig. 2A). The same numbers of cells were harvested and lysed. The same amount of total protein was used for affinity capture to ensure highly comparable conditions. Each sample was analyzed in technical triplicate for protein identification and quantification. A statistical platform based on a hierarchical linear model was utilized to search for enrichment of pro-teins from cells expressing NSD2-BirA. Due to affinity enrichment, the data are presented as an asymmetrical volcano plot where the x axis represents the relative abundance difference between NSD2-BirA versus BirA control and the y axis is the instantaneous q-value, a measure of how confident the abundance difference is statistically (Fig. 2B). A threshold of 2-fold difference and 5% false discovery rate was set for quantification, which responds to a -log 10 (instantaneous q-value) of 1.3. In total, 42 proteins were identified with statistically significant enrichment in the NSD2-BirA state, and among those, 24 were nuclear proteins (colored red in the plot). Because NSD2 localizes in the nucleus physiologically, we considered these 24 nuclear hits as confident interacting partners (Table S2). Variation in signal intensity for each quantified peptide was assigned as a percent of the total for each different source (Fig.  2C). Small variations in biological and technical replicates and large variations in NSD2-BirA versus BirA control suggested our quantification study was highly reliable.
GO analysis of the 24 interacting targets using EnrichR through KEGG pathway 2016 suggested NSD2-binding proteins are highly associated with the base excision repair pathway (Fig. 2D), which is responsible for removing base lesions A, study design: three biological replicates of TKO cells stably expressing NSD2-BirA or nuclear-localized BirA control were cultured with biotin, lysed, and quantified. Equal amounts of total protein from each sample was used for streptavidin capture and trypsin digestion. The resulting peptides were purified and analyzed by nanocapillary LC-MS in triplicate. High resolution MS1 was used for relative peptide quantification and MS2 was used for protein identification. B, volcano plot of quantified proteins is shown. The x axis is the fold-change in protein abundance of NSD2-BirA relative to BirA control. The y axis (instantaneous q-value) measures the statistical confidence of the variation between NSD2-BirA and BirA control. The vertical dotted lines represent a 2-fold change threshold. The horizontal dotted line indicates a false discovery rate of 5%. Twenty-four nuclear targets above threshold are colored red. C, box and whisker plot indicating variations in signal intensity of each quantified protein.
Variation attributed to the treatment (NSD2-BirA versus BirA control), biological replicates, technical replicates, or residuals (variation not explained by any of the previously described sources) are indicated. D, gene ontology KEGG pathway analysis of the 24 high confidence NSD2 partners. E, Venn diagram of overlapping hits from the BioID IP-MS qualitative Study 1 (see Fig. 1) and label-free quantitative Study 2.

PARylation impedes NSD2 binding to chromatin
from the genome throughout the cell cycle. The spliceosome is the second top pathway on the graph, suggesting that NSD2 might function beyond transcription elongation to facilitate spliceosome assembly or function.
The 63 hits identified in the BioID IP-MS qualitative Study 1 and 24 label-free quantitative hits in Study 2 have 16 overlapping targets (Fig. 2E). As described above, BRD4 was previously identified as a NSD2-interacting protein implicated in transcription elongation. In Study 1, BRD4 is one of the 63 high confident hits. Meanwhile, in the label-free quantitative study BRD4 was enriched 6-fold in NSD2-BirA cells, suggesting that proximity-based methods for identifying NSD2-associated proteins are mechanistically valid.

Poly(ADP-ribose) polymerase 1 (PARP1) is a NSD2-binding protein
PARP1 protein level increased 5-fold in the NSD2-BirA sample compared with the BirA control (Fig. 2B). It was also identified in all 6 biological replicates of BioID IP-MS analysis suggesting it is a true binding protein of NSD2. Because NSD2 may have a role in DNA repair, we chose to further investigate this interaction. The MM cell line KMS-11 and as well as KMS-11derived NTKO cells (5) that both are t(4;14)ϩ and express high levels of NSD2 were used to capture endogenous interaction of NSD2 and PARP1. Reciprocal co-immunoprecipitation results confirm that PARP1 and NSD2 physically interact under physiological conditions in MM cells (Fig. 3A).
PARP1 is a nuclear DNA-dependent ADP-ribosyltransferase that catalyzes poly(ADP-ribosyl)ation (PAR) of nuclear proteins (29). It serves critical roles in different molecular pathways by modifying itself and other proteins. The PARylation activity of PARP1 is necessary for mediating DNA damage repair. Previous proteomics-based approaches have been widely used to identify PARP1 PARylation targets (30 -34). Here we compared NSD2 interacting partners and PARP1 targets identified in multiple previously published papers (33,34). PARP1 targets were plotted in green on the volcano plot ( Fig. 3B) and indicated in Table S2. Interestingly, 17 of the 24 identified NSD2 interacting partners can be PARylated by PARP1. Furthermore, NSD2 itself was identified as a PARylation target in other screening experiments (33,34). Because it has been reported that PARP1 may be activated upon DNA damage and oxidative stress, we induced oxidative stress in NTKO cells by hydrogen peroxide treatment and captured PARylated proteins using an Archaeoglobus fulgidus macrodomain (35). Results indicate that PARP1 becomes auto-PARylated upon oxidative stress and enrichment of NSD2 after PARylated protein precipitation was observed in stressed cells but not in the control cells, suggesting that NSD2 can be PARylated upon DNA damaging stimuli in vivo (Fig. 3C).

NSD2 is a PARP1 substrate
Because NSD2 interacts with PARP1 under physiological conditions and was PARylated upon oxidative stress, we further tested if PARP1 was responsible for PARylating NSD2 in vitro. Recombinant NSD2 protein was incubated with PARP1, with and without NAD ϩ and DNA cofactors. When incubating with PARP1 enzyme, NAD ϩ and DNA, a ϳ50 kDa mass shift of NSD2 was observed indicating poly(ADP-ribosyl)ation (Fig. 4A, Fig. S1A, lane 5). The molecular weight shift corresponded to about 90 PAR units. In the same reaction, auto-PARylation of PARP1 was observed as an indistinct broad band in the PAR blot that blocked recognition with the PARP1 antibody (Fig. 4A,  Fig. S1A, lane 7). PARylation required cofactor NAD ϩ . Upon titration of NAD ϩ , NSD2 PARylation plateaued when PARylated NSD2 reached about 200 kDa, suggesting the PAR modification is highly regulated or only a finite amount of PAR can be added to NSD2 (Fig. 4B, Fig. S1B). As controls, the absence of DNA in the reaction as well as the presence of PARP1 inhibitor Olaparib significantly diminished PARP1 enzymatic activity (Fig. 4A, Fig. S1A, lanes 3 and 6).
Because several nonhistone proteins have been reported to be HMT substrates (36) and AURKA was recently identified as a novel nonhistone target of NSD2 that regulates p53 stability (37), we performed a NSD2 methyltransferase assay using PARP1 as the substrate. Significantly, no PARP1 methylation . PARP1 is a NSD2-binding partner. A, reciprocal co-immunoprecipitation of endogenous PARP1 and NSD2 from NTKO and KMS11 cells. B, volcano plot of NSD2-binding partners with PARP1 substrates colored in green. PARP1 substrate information was adopted from Refs. 33 and 34. C, NTKO cells were treated with 0.5 and 1 mM hydrogen peroxide for 15 min and harvested to extract nuclear proteins. PARylated proteins were then affinity captured with a poly(ADP-ribose)-binding macrodomain resin and immunoblotted for PARP1 and NSD2 as indicated.

PARylation impedes NSD2 binding to chromatin
was detected in vitro suggesting that PARP1 was not a substrate for NSD2 (Fig. 4C). Taken together, our data show that NSD2 is a substrate of PARP1 but not the other way around.

PARylation of NSD2 by PARP1 decreases NSD2 histone methyltransferase activity
Previous studies showed that PARylation can regulate enzyme activity (38,39). Therefore, we tested whether PAR modification of NSD2 affected its histone methyltransferase activity. To test this, NSD2 and PARP1 were incubated with and without the NAD ϩ cofactor required for PARylation prior to HMT assay. Olaparib was added into the reaction to terminate PARP1 enzyme activity (Fig. 5A). After incubation, an aliquot of the reaction was analyzed by immunoblot to measure PARylation of NSD2 (Fig. 5B, Fig. S2, lane 8). Significantly, the histone methyltransferase assay revealed that PARylated NSD2 has a 2-fold decrease of enzyme activity compared to unmodified enzyme (Fig. 5C, lane 8 versus 9). Controls of incubation of NSD2 with NAD ϩ , activated DNA, or Olaparib alone did not affect NSD2 activity (lanes 4 -6). Interestingly, the interaction of PARP1 and NSD2 alone did not have a significant effect on NSD2 histone methyltransferase activity (lane 7). Thus, our results suggest that PARylation of NSD2 may regulate its activ-ity and the PAR modification, but not the interaction of PARP1, and NSD2 inhibits NSD2 HMT function.

PARylation of NSD2 disrupts its association with chromatin
Next, we assessed whether the decrease of PARylated NSD2 HMT activity was due to loss of enzyme activity or inhibition of chromatin binding. To test this, we performed the PARylation reaction first, then provided SAM and biotinylated oligonucleosomes for the histone methyltransferase assay. Biotinylated nucleosomes were then captured and analyzed along with any associated NSD2 (Fig. 6A). Because the nucleosomes used were isolated from HeLa cells, basal H3K36me2 was observed (Fig. 6, B and C, Fig. S3, lane 2). Increased H3K36me2 was detected when unmodified NSD2 protein was in the reaction (lanes 3 and 5). However, when NSD2 was PARylated, the H3K36me2 level remained the same as without addition of NSD2 (lane 4), further supporting that PARylation of NSD2 inhibits its enzymatic function. Significantly, whereas unmodified NSD2 proteins were co-purified with biotinylated nucleosomes (Fig. 6B, Fig. S3, lanes 3 and 5), association of PARylated NSD2 with nucleosome was greatly reduced (lane 4 versus lane 3), suggesting that PARylation blocks the NSD2-chromatin interaction. NSD2 and PARP1 interaction without NAD ϩ had PARylation impedes NSD2 binding to chromatin no effect on NSD2 chromatin binding and concomitantly NSD2 could catalyze increased H3K36me2 levels (Fig. 6, B and C, Fig.  S3, lane 5). In summary, our data suggest the PARylation of NSD2 by PARP1 prevents it from binding to chromatin, resulting in decreased histone methyltransferase activity.

PARP1 PARylates NSD2 and reduces its recruitment to target genes
We further evaluated whether NSD2 can be PARylated by PARP1 in the cell. To test endogenous NSD2 PARylation by PARP1, we pre-treated NTKO cells with or without Olaparib and precipitated all PARylated proteins upon oxidative stress. As shown in Fig. 7A (left panel), NSD2 was PARylated and captured upon H 2 O 2 treatment. However, cells pre-treated with Olaparib to inhibit PARP1 activity had much less PARylated NSD2, suggesting PARP1 PARylates NSD2 in the cell upon oxidative stress.
Next, we investigated if PARylation affects NSD2 binding to chromatin in the cell. We performed a nuclear fractionation assay that separates soluble nuclear proteins from chromatinbound nuclear proteins. When treated with H 2 O 2 , soluble NSD2 increased and insoluble NSD2 decreased compared with nuclear loading controls, suggesting PARylation of NSD2 released it from chromatin binding (Fig. 7B, Fig. S4). The mod-est localization change of NSD2 might be due to only partial PARylation of NSD2 in the cell in response to H 2 O 2 . Alternatively, FRAP studies by our group and others showed that NSD2 is tightly bound to chromatin and exchanges with free molecules slowly (3,12). If for example, only free NSD2 were readily PARylated, a limited H 2 O 2 treatment would only be expected to yield a partial increase in net soluble NSD2 levels, given the slow release of bound NSD2 from chromatin.
Furthermore, we investigated how NSD2 PARylation affects individual gene loci. JAM2 and GLS2 were previously shown as NSD2-regulated genes and down-regulated when the NSD2 level was depleted (6). We performed ChIP-qPCR using primers that target both promoters and gene body regions of JAM2 and GLS2. NSD2 binding was significantly decreased at both gene loci when NTKO cells were treated with H 2 O 2 (Fig. 7C). Taken together, our results suggest PARP1 PARylates NSD2 upon oxidative stress to inhibit its ability to bind chromatin and recruitment to target genes in vivo.

Discussion
In this study, we mapped the interaction network of NSD2 in a multiple myeloma model by both qualitative and quantitative high precision MS. Under physiological conditions, NSD2 remains bound to the chromatin. The BioID method (ϳ10 nm

PARylation impedes NSD2 binding to chromatin
labeling radius) (40) enabled us to better capture the real neighboring proteins of NSD2. Twenty-four high confidence interacting proteins were identified including the previously described NSD2 interacting partner BRD4. PARP1 has drawn a large amount of attention in recent years due to the success of its inhibitors (Olaparib, Niraparib, and Rucaparib) in clinical trials (41). In our study, PARP1 was identified as a novel interacting partner of NSD2. NSD2 was previously identified in proteomics-based studies as a PARP1 substrate (33,34). We successfully validated the physical interaction and demonstrated that in MM cells NSD2 can be PARylated in the context of a DNA damaging stimulus. Our results further suggest that PARP1 PARylates NSD2 to dissociate it from chromatin binding and reduces its recruitment to target genes in the cell. Our findings are evidence of a novel mechanism by which PARP1 may regulate NSD2 localization.
A growing number of studies have revealed that PARP1 regulates gene expression through PARylation of histone methyltransferases and demethylases (42). On one hand, at the transcription start site of PARP1-regulated genes, PARP1 colocalizes with H3K4me3, a histone mark enriched at the promoter of activated genes, and protects this mark by inhibiting the corresponding histone demethylase KDM5B, consequently promoting gene expression (43). On the other hand, PARP1 facilitates repression of retinoic acid-dependent genes by PARylating and inhibiting H3K9 demethylase KDM4D to maintain H3K9me2, a repressive histone mark, at the promoter regions of these genes (44). Moreover, PARP1 directly PARylates EZH2, the catalytic subunit of the PRC2 complex that specifically catalyzes histone H3K27me3, a transcriptional repressive mark that antagonizes H3K36me2 (38,39). Our data show that PARP1 may also inhibit H3K36me2, which localizes to promoter regions and transcription start sites of actively transcribed genes through NSD2 activity. One conceivable functional role for PARP1 action on NSD2 may be that upon sensing DNA damage, PARP1-dependent PARylation of NSD2 may reduce the level of H3K36me2 at actively transcribed genes to inhibit gene expression.
H3K36me2 is enriched at DSBs upon DNA damage. However, the balance of methylation and demethylation is critical for the repair. Metnase is another enzyme responsible for writing H3K36me2 marks at DSBs to recruit DNA repair machinery proteins including NBS1 and Ku70. The loss of Metnase causes decreased H3K36me2 and deficient DNA repair (14). Similarly, reduced H3K36me2 due to loss of NSD2 in U2OS cells was associated with both decreased HR and NHEJ repair (23). However, a study in Caenorhabditis elegans suggests when JMJD-5, the H3K36me2 demethylase, is deleted or catalytic activity disabled, the animals are hypersensitive to IR and cannot resolve RAD51 foci during HR (18). Collectively, the information suggests that the level of H3K36me2 at the DSB sites might be crucial for proper DNA damage repair by both NHEJ and HR. As an important early stage sensor and regulator of DNA damage repair, PARP1 may regulate the H3K36me2 level by PARylating NSD2 and potentially other epigenetic regulators to maintain the histone mark balance.
NSD2 serves as an epigenetic writer that dimethylates histone H3 at lysine 36 to activate gene transcription. ChIP studies suggest that NSD2 is bound to chromatin globally across the genome in t(4;14)ϩ MM cells (45). However, under physiological conditions, how NSD2 regulates the activity of specific genes remains unclear. The identification of several transcription factors as binding partners with NSD2 in our study may Figure 6. PARylation of NSD2 disrupts its association with chromatin in vitro. A, schematic of the experimental strategy. After the PARylation reaction, SAM and biotinylated HeLa oligonucleosomes were allowed to bind to NSD2. Streptavidin magnetic beads were added into the reaction and incubated for 1 h at room temperature. B, biotinylated nucleosomes along with bound NSD2 were captured and analyzed by immunoblot against NSD2, histone H3K36me2, and total histone H4. C, histone H3K36me2 band relative intensities were quantified using ImageJ. Two independent experiments are shown.

PARylation impedes NSD2 binding to chromatin
provide hints at how NSD2 specificity is achieved. ZBTB33 (Kaiso) and GATAD2B are both transcriptional repressors and bind to the DNA in a sequence-specific fashion (46 -49). Further characterization of interactions of NSD2 to ZBTB33 and GATAD2B may show novel mechanisms of how NSD2 regulates downstream gene expression. Furthermore, the identification of DDX46 and SF3B2 as NSD2 interacting partners suggest a novel function of NSD2 in RNA splicing. Both DDX46 and SF3B2 are involved in pre-spliceosome assembly (50,51).
Interestingly, SF3BS is a nonhistone target of PRMT9, an arginine methyltransferase (52). Because little is known about nonhistone targets of NSD2 so far, the interaction of NSD2 with SF3B2 might provide a new idea of NSD2 regulating nonhistone proteins through its methyltransferase activity. Furthermore, NSD2 has a critical role in facilitating transcription elongation (24), and the identification of DDX46 and SF3B2 as binding partners might suggest a new mechanism of functional interactions between elongation and splicing. treated with 0.5 mM hydrogen peroxide for 15 min and harvested to extract nuclear proteins. PARylated proteins were then affinity captured with a poly(ADP-ribose)-binding macrodomain resin and immunoblotted for PARP1 and NSD2 as indicated. B, nuclear fractionation assay of KMS11 cells. Soluble and insoluble nuclear fractions were analyzed by immunoblot for NSD2. HDAC2 was used as soluble loading control and histone H3 was used as insoluble loading control. Relative soluble and insoluble NSD2 was quantified using ImageJ. C, ChIP-qPCR against NSD2 on the NSD2-regulated genes JAM2 and CLS2 promoters and gene bodies using NTKO cells with or without treatment of hydrogen peroxide. IgG was used as negative control. Enrichment was calculated as a percentage of total input DNA. Three independent biological replicates are shown.

PARylation impedes NSD2 binding to chromatin
In summary, this study identifies new interacting partners of the epigenetic regulator NSD2 and reveals a potential novel function in DNA damage repair. These results illustrate the importance for the balance and cross-talk between PARylation and histone modification in regulating DNA repair. Further studies of this regulation may provide new insights of therapeutic possibilities in multiple myeloma.

Experimental procedures
Cloning and cell culture KMS11-derived TKO and NTKO MM cell lines were cultured in RPMI 1640 media (Life Technologies) supplemented with 10% fetal bovine serum. Phoenix-AMPHO packaging cells were cultured in Dulbecco's modified Eagle's medium 1640 (Life Technologies) supplemented with 10% fetal bovine serum. The NSD2 cDNA was cloned into the retroviral vector pRetroX-IRES-DsRedExpress (Clontech) using NotI and ClaI restriction sites. pcDNA3.1 MCS-BirA (R118G)-HA was a gift from Kyle Roux (Addgene plasmid number 36047) (26). The BirA cDNA was PCR amplified from the pcDNA3.1 vector and inserted in-frame to the C terminus of NSD2 using ClaI and BamHI restriction sites. For the control, BirA with a C-terminal nuclear localization signal was cloned into the retroviral vector. Retroviruses to express NSD2-BirA and BirA control were generated by transfection of packaging Phoenix-AMPHO cells with the plasmids described above using Lipofectamine 3000 reagent (Invitrogen). One million TKO cells were resuspended in 2 ml of virus media along with 6 g/ml of Polybrene (Millipore) and then followed by spinoculation for 90 min at 2000 rpm. DsRed-positive cells were sorted (BD FACS Aria II) 72 h post-transduction and expanded.

BioID
BioID was performed as described in Ref. 27. For each experimental replicate, one 15-cm dish of TKO or TKO-NSD2-BirA cells was cultured and supplemented in the presence of 50 M biotin for 48 h, washed 3 times in PBS, lysed in BioID lysis buffer (50 mM Tris-HCl, pH 7.4, 500 mM NaCl, 0.2% SDS, 1 mM DTT, 1ϫ Halt protease inhibitor (Thermo Fisher Scientific)), and quenched. Supernatants were saved after centrifuging at 14,000 ϫ g for 10 min at 4°C. Total proteins were quantified by BCA assay (Pierce) and the same amount to total proteins in each replicate was used for streptavidin capture. Pierce TM Streptavidin Magnetic Beads (60 l, Thermo Fisher Scientific) were washed 3 times with lysis buffer and incubated with protein lysate overnight at 4°C on a tube rotator. Protein-bound magnetic beads were washed with three different wash buffers (27) and were then ready for digestion.
For on-bead digestion, magnetic beads were first washed twice in 100 mM ammonium bicarbonate, and bound proteins were reduced by addition of DTT to a final concentration of 10 mM and incubation at 50°C for 30 min. Alkylation was then done by incubation in 100 mM iodoacetamide (Sigma) at room temperature for 30 min in the dark. Bound proteins were then washed with 100 mM ammonium bicarbonate and digested with 1 g of trypsin (Promega) overnight at 37°C. The reaction was stopped by adding 2 l of formic acid. ZipTip (Millipore) cleanup was performed to desalt the samples. The resulting peptides were further dried in a SpeedVac centrifuge.

Mass spectrometry data acquisition
Nanocapillary LC-MS/MS analyses were performed with a 75-m ϫ 10.5-cm PicoChip column (New Objective, MA) packed with 1.9-m Reprosil C18 beads. A 150 m ϫ 3-cm trap packed with 3-m beads was installed in-line. Solvent A consisted of 0.1% formic acid in water and solvent B was 0.1% formic acid in acetonitrile. Peptides were resuspended in 30 l of solvent A for LC-MS, and 6 l were injected. For the label-free quantitative study, each sample was injected in technical triplicate. A total of 18 LC-MS/MS runs (3 biological replicates ϫ 3 technical replicates ϫ 2 states) were injected in random order for the data presented in Fig. 2. Peptides were trapped at 5 l/min for 5 min with 100% Solvent A, then separated at a flow rate of 300 nl/min with a gradient from 5 to 40% B in 35 min. After a 4-min ramp to 60% B, the column was washed at 95% B for 3 min and re-equilibrated to 5% B for 10 min with a total analysis time of 60 min.
The LC was coupled by electrospray to an LTQ Velos Orbitrap mass spectrometer operating in data-dependent MS/MS mode with a top 10 method. The capillary temperature was 275°C. MS1 scans were collected from 400 to 2000 m/z with resolving power equal to 60,000 at 400 m/z. The MS1 AGC was set to 1 ϫ 10 6 . Precursors were isolated with a 1.5 m/z isolation width, and the CID normalized collision energy was set to 35%. Dynamic exclusion was set to 60 s and charge 1ϩ ions were excluded as well. The ion trap AGC target was set to 3 ϫ 10 5 .

Protein identification
Mascot searching was performed for BioID protein identification. Raw files were converted to .mgf as input. Human Swis-sProt AC 2018_02 database (20,317 sequences) was used for the search. Trypsin was used as the proteolytic enzyme with up to two missed cleavages. Oxidation on methionine and deamination on asparagine and glutamine were chosen as variable modifications. Carbamidomethyl on cysteine was set as static modification. Peptide mass tolerance was 10 ppm; MS/MS tolerance was 0.6 Da.

Bioinformatic and statistical analysis
For the label-free quantitative study, three biological replicates were created using a different batch of cells on different days, and each biological replicate was injected in technical triplicate. One technical injection of TKO-BirA control biological replicate 1 failed due to a LC problem. This resulted in 17 LC-MS/MS runs used for this experiment. Proteome Discoverer (version 2.1; Thermo Fisher) was used for identification and peptide quantification. A spectrum selector was used to import spectrum from .raw files. The Sequest HT search engine was used to cross-reference identified peptides against human Uniprot fasta database. The mass tolerance was set to be 10 ppm for precursors and 0.6 Da for fragment ions. Oxidation on methionine and deamination on asparagine and glutamine were chosen as variable modifications. Carbamidomethyl on cysteine was set as static modification. Trypsin was set as the PARylation impedes NSD2 binding to chromatin proteolytic enzyme with up to two missed cleavages. The precursor ion area detector in Proteome Discoverer v2.1 was used for quantification.
In total, 2955 peptides were identified. Each protein was considered identified when at least two unique peptides were identified with at least two independent peptide spectral matches observed from the biological replicates. In total, 52 proteins were identified by 270 peptides. For differential expression analysis, all peptide intensities were standardized across all measures of the peptide. For each protein, each peptide species was treated as an independent set of observations on the protein. Standardized peptide intensities were tested for differential expression using a hierarchical linear model with technical replicates nested within biological replicates with a modified form of the technique described elsewhere (53,54). The resulting F test probabilities were then corrected for multiple testing (55) with a false discovery rate of 5%. All statistical analyses were done using a custom script for SAS 9.4 (Cary, NC). The SAS script used for differential expression calculations is included in supporting data Xiaoxiao_20170201.sas.zip.

Co-immunoprecipitation
Nuclear extracts were prepared using the Nuclear Complex Co-IP Kit (Active Motif, CA) according to the manufacturer's protocol and quantified by BCA assay (Pierce). Nuclear proteins (200 g) were incubated with 2 g of PARP1 antibody (39559, Active Motif), rabbit IgG (2729S, CST), NSD2 antibody (ab75359, Abcam), or mouse IgG2b (ab91366, Abcam) at 4°C overnight on a tube rotator. Next, 40 l of Protein G Mag-Sepharose beads (Sigma) were added into the mixture and incubated for 1 h. The protein complex bound beads were washed 3 times with PBS ϩ 0.5% Nonidet P-40 and 3 times with PBS. Bound proteins were then eluted by boiling for 10 min in 20 l of 1 ϫ SDS loading buffer and analyzed by Western blotting with antibodies against NSD2 and PARP1.

PARylated protein precipitation
Cells were washed with PBS and incubated in PBS with 0.5 or 1 mM hydrogen peroxide (Sigma) at 37°C for 15 min. Nuclear proteins were extracted and 200 g of were diluted in 500 l of PAR lysis buffer (50 mM Tris-HCl, pH 8.0, 200 mM NaCl, 1 mM EDTA, 1% Triton X-100, 10% glycerol, 0.5% deoxycholate, 1 mM DTT, 1ϫ Halt Protease Inhibitor). In total, 20 l of poly-(ADP-ribose) binding macrodomain resin (Tulip Biolabs, PA) was pre-washed with PAR lysis buffer and added into the nuclear protein mixture to incubate overnight at 4°C on a tube rotator. The beads were washed 3 times with PAR lysis buffer. Bound PARylated proteins were eluted by boiling for 10 min in 20 l of 1ϫ SDS loading buffer and analyzed by immunoblot with antibodies against NSD2 and PARP1 (ab194586, Abcam).

In vitro PARylation assay
Recombinant NSD2 methyltransferase (0.25 M, Reaction Biology, PA) was incubated with 0.25 M PARP1 enzyme (Tulip Biolabs, PA) in 10 l of PARylation buffer (50 mM Tris-HCl, pH 8.0, 100 mM NaCl, 10 mM MgCl 2 , 1 mM DTT). Other components including 10 g/ml of calf thymus DNA (Sigma), 50 M NAD ϩ , and 500 nM Olaparib (Selleckchem) were added to the reaction as indicated. The reactions were incubated at room temperature for 10 min and inactivated by adding 10 l of 2ϫ SDS loading buffer, and boiling for 2 min. The elution was resolved by SDS-PAGE and immunoblotted with antibodies against NSD2, PARP1, and poly(ADP-ribose) polymer (ab14459, Abcam).

In vitro histone methyltransferase assay
PARylation reaction containing 0.25 M recombinant NSD2, 0.25 M PARP1 enzyme, 10 g/ml of DNA, and 50 M NAD ϩ in 15 l of PARylation buffer was set up to PARylate NSD2. After a 10-min incubation at room temperature, 500 nM Olaparib or ddH 2 O was added to stop the reaction. One-third (5 l) of the reaction was removed to confirm PARylation by immunoblot. The remaining 10-l reaction was diluted in 40 l of HMT buffer (50 mM Tris-HCl, pH 8.5, 50 mM NaCl, 10 mM MgCl 2 , 2 mM tris(2carboxyethyl)phosphine) with 1 M HeLa oligonucleosomes (Reaction Biology, PA) and 2 M [ 3 H]SAM (PerkinElmer Life Sciences) and incubated for 1 h in a 30°C water bath. The HMT reaction was terminated by adding trichloroacetic acid (TCA) to a final concentration of 10%. Precipitated proteins were spotted onto glass microfiber filters (GE Healthcare) and washed with icecold 10% TCA and 100% ethanol. The filters were air dried in the chemical hood and measured by liquid scintillation counting for 3 H incorporation.

NSD2-nucleosome association assay
A PARylation reaction of 15 l was set up as described above. After terminating with Olaparib, the reaction was diluted in HMT buffer with 1 M biotinylated HeLa oligonucleosomes (HMT-35-160, Reaction Biology, PA) and 2 M SAM (New England Biolabs) to a final volume of 60 l and incubated for 30 min at 30°C. Following incubation, 20 l of the reaction was kept as the input fraction. The remaining 40 l of the reaction was diluted in 460 l of wash buffer (50 mM Tris-HCl, pH 8.5, 100 mM NaCl, 5 mM MgCl 2 , 2 mM tris(2-carboxyethyl)phosphine, 0.5% Nonidet P-40, Halt protease inhibitor) with 50 l of pre-washed streptavidin magnetic beads (88817, Pierce) and incubated at room temperature for 1 h. Nucleosome-bound complexes were washed three times in wash buffer and eluted by boiling in 30 l of 1ϫ SDS loading buffer. Input and elution were resolved by SDS-PAGE and blotted with antibodies against NSD2, PARP1, histone H4 (ab7311, Abcam), and H3K36me2 (2901S, CST).

Nuclear fractionation assay
Nuclear fractionation assay was performed according to the manufacturer's protocol (Subcellular Protein Fractionation Kit for Cultured Cells, Thermo). Soluble and insoluble nuclear protein fractions were saved and quantified by BCA assay (Pierce). The same amount of soluble and insoluble proteins from hydrogen peroxide-treated and untreated samples were resolved by SDS-PAGE and blotted with antibodies against NSD2, HDAC2 (05-814, Millipore), and histone H3 (9715S, CST).

Data availability
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE (56) partner repository with dataset identifier PXD013759.