Abstract
Nuclear receptor SET domain protein (NSD2) plays a fundamental role in the pathogenesis of Wolf–Hirschhorn Syndrome (WHS) and is overexpressed in multiple human myelomas, but its protein–protein interaction (PPI) patterns, particularly at the isoform/exon levels, are poorly understood. We explored the subcellular localizations of four representative NSD2 transcripts with immunofluorescence microscopy. Next, we used label-free quantification to perform immunoprecipitation mass spectrometry (IP-MS) analyses of the transcripts. Using the interaction partners for each transcript detected in the IP-MS results, we identified 890 isoform-specific PPI partners (83% are novel). These PPI networks were further divided into four categories of the exon-specific interactome. In these exon-specific PPI partners, two genes, RPL10 and HSPA8, were successfully confirmed by co-immunoprecipitation and Western blotting. RPL10 primarily interacted with Isoforms 1, 3, and 5, and HSPA8 interacted with all four isoforms, respectively. Using our extended NSD2 protein interactions, we constructed an isoform-level PPI landscape for NSD2 to serve as reference interactome data for NSD2 spliceosome-level studies. Furthermore, the RNA splicing processes supported by these isoform partners shed light on the diverse roles NSD2 plays in WHS and myeloma development. We also validated the interactions using Western blotting, RPL10, and the three NSD2 (Isoform 1, 3, and 5). Our results expand gene-level NSD2 PPI networks and provide a basis for the treatment of NSD2-related developmental diseases.
Introduction
The nuclear receptor SET domain containing protein 2 (NSD2), also known as MMSET or WHSC1, is a member of the NSD protein family (), which mainly catalyzes histone H3 lysine (; ). NSD2 is a critical gene in the pathology of Wolf–Hirschhorn Syndrome (WHS; ; ), a severe neurodevelopmental disorder characterized by distinctive developmental delays (DDs), intellectual disabilities (IDs), and seizures, which occur in more than 50% of WHS infants (). The disease results from distal deletions on the short arm of chromosome 4 (chromosome 4p16.3) (), which occurs in 1 in 50,000 births (). NSD2 carries rare mutations in patients with neuropsychiatric disorders, including autism spectrum disorder (ASD), DDs, IDs, and schizophrenia (SCZ; ; ; ; ). Early studies have reported various deleterious NSD2 variants in neuropsychiatric patients, suggesting that the haploinsufficiency of NSD2 might be partially responsible for DDs (; ; ). However, over the past decade, most functions of NSD2 have been identified in carcinogeneses, such as renal cell carcinoma (), colorectal cancer (), osteosarcoma (), and multiple myeloma (). The biological divergence and potential mechanistic differences in NSD2 associated with neurodevelopmental disorders and cancers remain undetermined.
Compelling data have emerged to support the concept that alternatively spliced isoforms are linked to a range of functional characteristics of certain genes (; ; ) and contribute to functional complexity in diseases (Yang et al., 2016; ). The alternative splicing of pre-mRNAs is widespread in humans and most eukaryotes (; ), and it happens in ∼95% of genes containing different numbers of exons (; Wang et al., 2008). There are striking functionally diverse gene functions in human brains in particular due to alternative splicing (). Although alternative splicing is generally responsible for the diversity of gene products expressed from the genome, the complexity of alternative splicing at the proteome level remains to be characterized (). In addition, large-scale proteomics experiments are usually only focused on a single gene-level protein approach to simplify the number of proteins for use in further analyses. Moreover, most large-scale experiments have relied on antibodies recognizing a region common to different isoforms or have chosen the most characterized protein isoform to include to identify protein interactions or expression patterns. However, studies on NSD2 isoforms have been limited (). De novo variants (DNVs) play a vital role in understanding the genetics of psychiatric disorders (Veltman and Brunner, 2012). There have been 52 DNVs found in psychiatric patients, including nine DNVs that affect protein-coding regions (). Among these exonic DNVs, two are found in ASD (; ), five in DDs () or IDs (; ), and only one in SCZ () and congenital heart disease (). These DNVs affect different exons among NSD2 isoforms, which contribute to phenotypic differences.
In addition, integrating protein–protein interactions (PPIs) to study the potential functional impacts of risk genes on the associated disorders at the level of the biological system is a common practice in the study of disease biology (; Zhang et al., 2018; ). For example, ASD-associated physical interaction networks formed by protein interactome, which focus on cancer-related genes such as β-catenin (), p53 signaling (), Wnt–β-catenin (), and MAPK (), have provided important insights into the interpretation of diseases. In addition, protein interactomes can also be useful for investigating key pathways, such as abnormal synaptic phenotypes () and post-synaptic density () in ASD focusing on targeted functional gene sets. Using 343 WHS-associated genes, including NSD2 itself, have constructed a PPI network with GeneMANIA (), which was able to identify a gene set with a role in NAD+ nucleosidase activity. However, PPIs at the gene level may not completely reflect the complex system underlying disease etiology, particularly because several studies have shown that different isoforms of the same gene can differ in both biological function and subcellular component (; ). In addition, it has been shown that alternatively spliced isoforms of the same gene can have different sets of interaction partners (), and the interactome analyses of the isoforms of interest would facilitate the identification of their functional biological roles (Yang et al., 2016).
We also found that the interactors for NSD2 shed light on the function of histone methyltransferase activity (), and they supply resources in constructing the interactome of NSD2 (; ). However, these interactome analyses of NSD2 have mostly remained at the gene level, limiting downstream functional analyses and leaving the PPI patterns at the isoform/exon levels poorly understood. Indeed, the main genomic databases, e.g., RefSeq () and Ensembl (Zerbino et al., 2018), apply different identifiers for NSD2 transcripts, and previous reports have merely specified which transcript or protein isoforms are under consideration (; ; ; ). A few studies have focused on NSD2 transcript variant 1 (accession number: NM_133330.2) (), but the remaining isoforms of NSD2 have been poorly explored. Therefore, isoform-specific interactome analyses of NSD2 are urgently needed.
To fill this gap and yield further insight into the divergent etiology of NSD2 isoforms, we experimentally screened selected exon-specific isoforms of NSD2 and their protein-level interactions, then systematically incorporated the interactions into a study of the biological network to investigate potential functional interaction networks for different NSD2 isoforms (Supplementary Figure 1). We selected four representative transcripts, namely, Isoforms 1, 3, 5, and 7, and then we used label-free quantification to perform immunoprecipitation mass spectrometry (IP-MS) analyses. Using isoform-specific partners, we analyzed the diversity between the differential networks of NSD2 cleavage isomers and explored the underlying molecular functional pathways. We identified 365 proteins as novel interactors with NSD2 isoforms. These novel partners were significantly enriched in proteins that are functionally characterized as RNA splicing, in addition to literature-reported NSD2 interactors, whose functions relate to the chromatin remodeling. One of the novel partners, the ribosomal protein gene RPL10, a genetic factor to cognitive function (), was specially partnered by Isoform 1, 3, and 5. We further validated the interaction between RPL10 and NSD2 isoforms experimentally using Western blotting, supporting the credibility of our data and extending the knowledge of the protein partners of NSD2.
Results
Cloning, Expression, and Localization of NSD2 Isoforms
Alternative splicing of exons resulted in 27 NSD2 transcripts, of which 14 can be translated into proteins.1 We selected four of them because of their mutually excluded exons and domains (Figure 1A): Isoform 1 (ENST00000382891.9) contains exons 4, 7, 8, 10–15, and 17–29; Isoform 3 (ENST00000398261.5) contains exons 7, 8, 10–15, and 16; Isoform 5 (ENST00000420906.6) contains exons 4–8 and 10–15; Isoform 7 (ENST00000436793.5) contains exons 7–9 (Supplementary Figure 2). We named our NSD2 isoforms using the same IDs from the UniProt (UniProt Consortium, 2019). It is worth noting that all isoforms except for Isoform 1 have unique exons. For example, Isoform 3 uniquely contains exon 16, Isoform 5 uniquely contains exons 4–6, and Isoform 7 uniquely contains exon 9.
FIGURE 1
It has been reported that tissue-specific expression patterns happen in a substantial proportion of isoforms generated from the same gene due to alternative splicing events (; ; Yang et al., 2016). To establish the differences among the four isoforms of NSD2 investigated in this study, we first looked at the transcript-level expressions for four isoforms using GTEx (). We observed that the transcripts have specific expression patterns across the tissues of the testis, bone marrow, lymph node, and brain (Figure 1B). The highest level of expression was detected in Isoform 1, followed by Isoforms 3 and 5, while Isoform 7 was expressed at a very low level. Based on the transcript-level quantifications of bulk RNA-Seq, the expression of Isoform 7 in tissues was not detected. We also used quantitative real-time polymerase chain reaction (qRT-PCR) to confirm the expression levels of Isoform 7 in the HEK293T cell line, even at a low level (Supplementary Figure 3). Collectively, the variation in expression patterns of the NSD2 isoforms strongly suggests their functional divergence.
Subcellular localization could indicate the extent to which enzymatic activities can be regulated by the products of isoforms (; ; ). To explore the potential function of each protein isoform, we performed the subcellular location detection and investigated localization characteristics of the four NSD2 isoforms via immunofluorescence (IF) microscopy (Figure 1C). Each isoform vector was first transfected into HEK293T cells, and the expression of each subtype was induced. The cells were then fixed before being subjected to IF labeling using FLAG antibody, and the nuclei were observed with DNA staining using DAPI. As expected, the imaging of the isoforms showed convergent and divergent subcellular distributions (Figure 1C). These data showed whether isoforms expressed in the nucleus, with the nuclear accumulation of Isoforms 1 and 7 and cytoplasmic and nuclear accumulations of Isoforms 3 and 5.
NSD2 Isoform Interactomes in HEK293T Cells
Previous studies of alternative splicing isoforms and their PPIs have shown that different isoforms of the same proteins can have variable biological functions, ranging from similarities in binding partners to completely different sets of partners (; Yang et al., 2016; ). In addition, due to the differences in their tissue expression and subcellular localization, we hypothesized that NSD2 isoforms could also vary in their protein binding targets. To investigate isoform-specific interactors, in vitro affinity-capture assays, coupled with label-free quantification of interacting proteins, were performed in HEK293T cells using isoform-Flag recombinant proteins as bait and Flag alone as a control to subtract non-specific interactions (Figure 2A). We first constructed four expression vectors for NSD2 Isoforms 1, 3, 5, and 7, and a blank vector served as control. Then we transfected different isoform vectors into 293T cell lines and tested the isoform vectors’ stably expressed condition to ensure that there was no degradation (Figure 2B). To identify the interacting partners of NSD2 isoforms, we used label-free quantification to perform MS analyses of the isoform and control sample and obtained IP-MS data for the five samples. These experiments were done in triplicate. To validate the overall quality of the MS results.
FIGURE 2
Isoform-specific interactors were identified by comparing the MS results for NSD2 isoforms and control, and we obtained 383 protein partners interacting with the four isoforms encoded by NSD2. First, we filtered the detected partners with the cutoff value of the peptide frequency that occurred more than once in a single test and appeared twice or more in the three biological repeat tests. We found that Isoform 1 had 205 interacting partners, Isoform 3 had 287, Isoform 5 had 224 interacting, and Isoform 7 had 167. There were overlapping protein partners between these isoforms established by comparing the different isoforms (Figure 2C). Finally, we identified 20 interacting partners specific to Isoform 1, 80 to Isoform 3, 19 to Isoform 5, and 30 to Isoform 7 by comparing the unique protein partners.
To investigate the extent to which the four isoforms mediate interactions with different partners, we evaluated the dissimilarities in their interaction profiles by calculating the Jaccard distance of every pairing of four isoforms (Figure 2D). We restricted our analyses by comparing our paired NSD2 isoforms with the validated interaction of the 105 isoforms reported by Yang et al. (2016). Globally, we found that most isoform-specific protein partners have not been reported, except for 18 overlapping interactions between the gene level and the isoform level. By comparison with the NSD2 interactors obtained from BioGrid (), we found that a substantial proportion of interacting partners was exclusively identified in our isoform-level data. Only 4% of the partners were repeated at both the gene and isoform levels. Another 13% were found in gene-level PPIs, and the remaining 83% were only identified at the isoform level (Figure 2E). The targeted isoform partners exhibited a shrinking percentage of PPIs (67–78%), as the equal proportions of novel interactors showed by other reports (; Yang et al., 2016), emphasizing the importance of isoform-level exploration for protein interaction networks. As Huang et al. identified high confidence NSD2 interacting partners in MM cells (), we compared our interactors of NSD2 isoforms with their NSD2 partner proteins and found the resemble trend that NSD2 partners in isoform level contribute more in the interaction network.
NSD2 Isoform-Specific Partners That Indicate Distinct Disease Mechanism
Because novel NSD2 interactors were primarily found at the isoform level, we explored whether the interactome network construction of these interactors illustrated different functional pathways. Because the interacting partners of a single node (gene/isoform) in a network could have notably different properties than those of proteins that interact with separate nodes (), we reasoned that combining the direct binary partners of NSD2 isoforms in gene/isoform subnetworks may reveal functional differences between them. Based on the literature-reported partners curated from BioGrid (), we first extended our analyses by merging both novel and literature-reported partners into the pool to construct a gene-isoform network for NSD2 (Figure 3A). This extended network included all 991 interactions and combined with ∼90% of isoform-specific partners globally (Figure 3B). In the resulting network, Isoform 5 reached the degree of 294, the highest found. The gene NSD2, which is only equipped with literature-reported PPIs, had only 92, the lowest degree. This expanded protein interaction capability of NSD2 suggests a functional divergence among the four isoforms.
FIGURE 3
We examined whether the isoforms’ extended protein networks were significantly enriched in a range of functional categories involved with disorders to investigate this functional divergence. For the literature-reported NSD2, the significant enrichment terms include RNA splicing, chromatin remodeling and histone modification (Figure 3C), in which NSD2 has been documented to play an important role (). For our isoform specific interaction partners, RNA splicing and DNA conformation featured the most significant enrichment and showed more importance than other functional terms. This accumulating evidence suggests that the isoform-specific interactions of NSD2 play vital roles in the subnetwork.
To confirm the relationship between RNA splicing and NSD2 isoforms, we divided our interactome into Reference PPIs, which are interactions of Isoform 1 that has been commonly used to represent NSD2 in gene-level studies, and Non-reference PPIs, which are interactions from all other isoforms. The varied results of enrichment suggest caustic usage of NSD2 isoforms (Figure 3D). Isoform 1 interactions were significantly enriched in DNA conformation and RNA splicing. By contrast, the other three isoform interactions were enriched in regulating RNA stability and transport of virus and differed from the Reference PPI. Indeed, the isoform PPIs of the NSD2 indicate a potential distinguishing underlying mechanism, particularly for these exon-specific matched DNVs from patients with different psychiatric disorders.
Isoforms Specific to NSD2 Exon Junctions Associated With Distinct Pathways
Because splicing could mediate the disruption of interactions through its inclusion or exclusion of domains, the targeted domains can be predicted to interact with interacting partner proteins that contain a certain region (; Yang et al., 2016). Following the hypothesis that isoforms holding the same exons share a regulation or function (), we prioritized them based on the exon composition of the included isoforms to investigate the role more deeply of alternative exons in isoform-level interactions of NSD2. We proposed that the isoforms that constituted the same exons could have the same protein domains and/or similar biological functions. We considered the difference between NSD2 isoforms that include or exclude a target exon. For this purpose, we performed analyses to categorize isoforms from the census exon into clusters (Figure 4A). This procedure divided the interacting partners into seven groups that shared different numbers of genes (Figure 4B). We found that the cluster for exon 7–8 was mainly concentrated in the protein targets to ER, mRNA catabolic process, and translational initiation; the cluster for exon 4 was mainly related to RNA splicing, and RNA transport; and the binding proteins in the group for exon 9 were enriched in transport of virus, NF-kB signaling pathway, and response to unfolded proteins.
FIGURE 4
To examine the consequences of the exon-specific partners of isoforms, we selected two representative genes, RPL10 and HSPA8, for their presentation in the subnetworks constructed by Isoforms 1, 3, and 5 and in all four isoforms separately (Figure 4B). Because the domain PWWP 1 was included in Isoforms 1, 3, and 5 but not Isoform 7 (), the binary interaction with RPL10 suggested that isoform-specific partner differences could be explained by the alternative splicing of protein domains. RPL10, located on the Xq28 chromosome, has an essential function in ribosome assembly and protein translation, and it is associated with neurodevelopmental disorders (; Zanni et al., 2015). We then explored the protein interactions between different isoforms of NSD2 and RPL10 by co-immunoprecipitation (Figure 4C). The result showed that RPL10 could bind to Isoform 1, 3, and 5, but there was no detectable interaction between RPL10 and Isoform 7. Moreover, the interactions between HSPA8 and the four isoforms were also confirmed by co-immunoprecipitation, which supported our findings in the IP-MS data that all four isoforms interacted with HSPA8 (Figure 4D).
Discussion
We identified functional diversity in NSD2 isoforms and implicated several novel protein interactors associated with NSD2. We also verified novel protein partners via co-immunoprecipitation. Similar to the results of other studies of expanding protein interaction networks by isoform-specific level (; ; Tseng et al., 2015), our data indicate a significant extension of interactors to NSD2, showing different enrichment categories at the isoform level.
Nuclear receptor SET domain protein encodes a protein that contains four domains: a PWWP domain, an HMG box, a SET domain, and a PHD-type zinc finger (). The NSD2 gene has 29 exons, which are combined in different ways to construct 27 transcripts. The alternative splicing of NSD2 results in multiple transcripts encoding for different protein isoforms. Since some transcripts are nonsense-mediated mRNA decay candidates, they are not represented as reference sequences.
It has been demonstrated that NSD2 is strongly associated with tumorigenesis by promoting histone methylation, which is crucial for transcriptional regulation and chromatin remodeling (; ; Tanaka et al., 2020). NSD2 alternative splice isoforms in cancers, in particular, have received substantial attention as of late (; ). The deletion of NSD2 can cause WHS (; ), which has clinical characteristics that overlap with ASD. The genetic etiology of ASD is heterogeneous, attributed to hundreds of genes, of which only a small percentage has sufficient evidence to support being considered as a cause (; ). In this study, we have systematically incorporated the interactions of the isoforms of the selected candidate gene NSD2 into a biological network study. By inspecting the convergence and divergence of the differential network of NSD2 isoforms, we explored functional molecular pathways related to NSD2, which can serve as new therapeutic targets for ASD.
To interpret DNVs at the NSD2 isoform level, we selected four representative transcripts because their exon junctions are responsible for alternative splicing. We also curated likely damaging DNVs from PsyMuKB and found their specific mappings in the above three transcripts, except Isoform 7, which does not include exons 10–15, where there were four nonsense variants explored in diseases (; ; ; ). It is worth mentioning that the SET domain-containing protein 2 gene, which has the same domain as NSD2, has also shown a de novo gene-damaging mutation through whole-exome sequencing (). We also found that NSD2 is expressed across different tissues, which indicates that its dysregulation may lead to disease. In conclusion, our data emphasize that isoform specificity plays a critical role in the various biological processes. The state of differentiation for each NSD2 isoform resides in the exon junction methods. As the mutations that affect various isoforms have different exons, they have different impacts on the isoforms of the same gene (Zhang et al., 2014). The four varied N-terminus representative isoforms (Isoforms 1, 3, 5, and 7) examined here enabled us to explore exon junction clusters, a prerequisite for the functional enrichment of interactors.
NSD2 may play a role in regulating ribosome assembly and protein translation by binding its partner RPL10. Because RPL10 disrupts neurodevelopment, and a rare mutation of it has been found in ASD (; ; ; Takumi and Tamada, 2018), its connection with NSD2 also implies that both genes may contribute to ASD. Moreover, the interaction between NSD2 and RPL10 has been previously identified by . We have also validated the interaction between them with Western blotting. In this study, we did not carry out a detailed examination of NSD2 regulation of substrates’ stability, such as HSPA8 and RPL10, which could be performed in subsequent work where the specific disease pathway and upstream and downstream analyses of NSD2 in ASD can be clarified to establish the underlying mechanisms.
One limitation of the study is the absence of verified interactions between NSD2 isoforms and partner proteins in physiological conditions. As the PPIs confirmed under physiological conditions will enhance the understanding of their functions in vivo, the interactome of our isoforms could provide a straightforward functional annotation for NSD2. Unfortunately, as the commercial antibodies for these isoforms are not available, we therefore have chosen to valid their interactions of exogenous isoforms and their partners in this study. Even though these interactions of exogenous proteins might not simulate the biological systems directly, our present results could still provide some preliminary insights into the protein interaction exploration.
To the best of our knowledge, this study was the first time that the relationship between different gene-level and isoform-level interactions of NSD2 were elucidated and where the difference between functional enrichments was established. Our data indicate that interactor partners are significantly expanding at the isoform level, and different metabolic pathways are found beneath DNV-induced disorders. The use of conformational extension to induce interactors at the isoform-level makes this elusive network PPI generally available and has the potential to shed light on the biological pathways underlying a range of developmental disorders.
Materials and Methods
Plasmids and Reagents
NSD2 full-length plasmids were gifts from Lili lab. Full-length and other isoforms of NSD2 were PCR amplified and cloned into pcDNA5-Flag to generate Flag-tagged fusion proteins. RPL10 and HSPA8 were cloned into pcDNA5-Myc to generate Myc-tagged fusion proteins. These constructs were cloned into pcDNA5 using BamHI and XhoI restriction sites. N-terminal Polη truncations were generated with 5’ and 3’ primers containing BamHI and XhoI restriction sites, respectively, and cloned into pcDNA5. Primers were used:
| 1-NSD2-F-BamHI | CGCGGGATCCATGGAATTTAGCATC |
| 1-NSD2-R-XhoI | CGCGCTCGAGCTATTTGCCCTCTGT |
| 3-NSD2-F-BamHI | GGATCCATGGAATTTAGCATCAAGCA GAGTCCCCTTTCTGTTCAGAGTGTTG TAAAGTGCATAAAGATGAAGCAGGC |
| 3-NSD2-R-XhoI | CTCGAGCTAAGTGCAGTACAGAGCAG CTGGGTTCAAATCCAACTTGACTGGT GTGGGCTCCCACAAAAGC TCATTCTCAGTTAAGGA |
| 5-NSD2-F-BamHI | CGCGGGATCCATGGAATTTAGCATC |
| 5-NSD2-R-XhoI | CGCGCTCGAGTTATTTTACCTCATT CTCAGT |
| 7-NSD2-F-BamHI | GGATCCATGGAATTTAGCATCAAGCA GAGTCCCCTTTCTGTTCAGAGTGTTG TAAAGTGCATAAAGATGAAGCAGGC |
| 7-NSD2-R-XhoI | CTCGAGTTAATCTTTCAGTACAATTT GACTTGTTTTTAAGTGTTCAAACTTC TTTGATTTGAAAATACCTTTAAGTTT GGTATAGCTG |
Anti-Flag M2 agarose affinity gel was purchased from Sigma (#A2220). Antibody against Myc was from Cell Signaling Technology (#2276). Antibody against Flag was from Cell Signaling Technology (#14793).
Cell Culture and Reagents
293T cells were obtained from Lili lab. All cells were cultured in DMEM medium supplemented with 10% fetal bovine serum (FBS) at 37°C in the presence of 5% CO2 if not specified. For transient transfection experiments, cells were transfected with indicated constructs using jetPRIME (Polyplus-transfection) following the manufacturer’s protocol. Forty-eight hours later, transfected cells were collected for further experiments.
Quantitative Real-Time Polymerase Chain Reaction
RNAs were extracted with Trizol reagent (Thermo Fisher Scientific) according to the manufacturer’s instructions. The cDNA was obtained using PrimeScript RT reagent Kit with gDNA Eraser (Perfect Real Time) (Takara). Primers were designed using Primer3 version 4.0.0. The qRT-PCR assay was performed using a 20-μl reaction system with SYBR Green Master reagents (Roche) and the designed primer mixtures in ABI 7900 HT Real-time PCR system (Applied Biosystems). The reaction system contained 10 ul SYBR Green Master (ROX), 0.2 μl of each primer (10 μM), 2 μl template (about 25 ng/μl cDNA) and ddH2O. Initial denaturation was at 95°C for 5 min followed by 40 cycles of 95°C denaturation for 10 s, 55°C annealing for 20 s, and 72°C extension for 20 s. GAPDH was used as an internal control. Relative quantification (RQ) was derived from the cycle threshold (Ct) using the equation RQ = 2–ΔΔCt. The forward primer for Isoform 5 is 5′-ACCCATCAGAGTGTTCTA-3′, and reverse is 5′-GTGCCTGCTTCATCTTTA-3′. The forward primer for Isoform 7 is 5′-GACCACCTGTTGAAATAC-3′, and reverse is 5′-TCTTTGATTTGAAAATACCTTTA-3′.
Immunofluorescence
Briefly, cells were seeded on cover glasses and irradiated with ultraviolet C (UVC). The cells were permeabilized with 0.5% Triton X-100 for 5–30 min before being fixed in 4% paraformaldehyde. The samples were then blocked with 1X PBS/5% normal serum/0.3% Triton™ X-100 for 60 min. The cells were next incubated with indicated antibodies overnight at 4°C, followed by incubation with Alexa Fluor 568 goat anti-mouse (Invitrogen, Molecular Probes) for 60 min. The cells were later counterstained with DAPI, and images were acquired with a Leica DM5000 (Leica) equipped with HCX PL S-APO 63 × 1.3 oil CS immersion objective (Leica) and processed with Adobe Photoshop 7.0.
Co-immunoprecipitation and Western Blotting
HEK293T cells were transfected with Flag-NSD2 and Myc-RPL10 or Myc-HSPA8. Forty-eight hours later, the cells were harvested and lysed with Cell lysis buffer for Western and IP (20 mM Tris (pH 7.5), 150 mM NaCl, 1% Triton X-100, and sodium pyrophosphate, β-glycerophosphate, EDTA, Na3VO4, leupeptin). The whole-cell lysates were immunoprecipitated with anti-Flag M2 agarose in the presence or absence of RNase A, ethidium bromide (EB). For mapping the regions within NSD2 responsible for its interaction with HSPA8 and RPL10, a Flag-tagged vector (Ctl) and a series of NSD2 isoforms were co-transfected with Myc-HSPA8 or Myc-RPL10 in HEK293T cells for co-immunoprecipitation experiments. The immunoprecipitated products were separated by SDS-PAGE and detected by immunoblotting with indicated antibodies.
Immunoprecipitation-Mass Spectrometry
We used a combination of immunoprecipitations to study the NSD2 interaction partner, followed by qualitative mass spectrometry using LTQ-ESI-MS to measure the difference in the interaction between each isoform and contaminant protein. To identify proteins that specifically interact with each isoform of NSD2, we transfected cells expressing Isoform 1, Isoform 3, Isoform 5, or Isoform 7 expressing NSD2 into 293T cells, and 48 h later, lysed cells were collected. Protein and each experiment were repeated three times (biological replicate). After lysis, the total cell lysate was mixed for immunoprecipitation. These experiments were repeated, thus representing biological repeats. The immunoprecipitated protein was boiled and deformed and then subjected to SDS-PAGE to ensure the sample’s quality after immunoprecipitation. After determining the immunoprecipitation results of different transcripts, the silver strips were taken for mass spectrometry analysis. Each sample of mass spectrometry data was decontaminated and screened to obtain potential interacted proteins.
Protein–Protein Interaction Data
The PPI data were downloaded from the BioGrid (v4.2.191) (). BioGrid contains PPI with various detailed information, such as sources and experimental methods. We used the detailed PPI annotation provided by BioGrid to obtain direct physically interacted protein interactors. Network analyses and their visualizations were constructed by Cytoscape().
Gene Ontology Analysis
The Gene Ontology (GO) enrichment analysis of interested gene lists was performed by ClusterProfiler (Yu et al., 2012) R package. We tested if genes of interest enriched in any GO-BP pathway by hypergeometric test. Gene background was defined as all genes with GO annotation. P-value of hypergeometric tests was adjusted for multiple testing by the Benjamin–Hochberg method.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Author contributions
GNL and LL conceived and directed the project. WW and WS curated and processed all of the data. JZ, YC, and LC participated in experiments. GNL, LL, WW, YC, and WS wrote and edited the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by grants from National Natural Science Foundation of China (No. 81671328 and 81971292), Program for Professor of Special Appointment (Eastern Scholar) at Shanghai Institutions of Higher Learning (No. 1610000043), Innovation Research Plan supported by Shanghai Municipal Education Commission (ZXWF082101), and Shanghai Mental Health Center (2019-YJ01 and 2019-QH-03).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.612019/full#supplementary-material
Footnotes
References
1
Allali-HassaniA.KuznetsovaE.HajianT.WuH.DombrovskiL.LiY.et al (2014). A basic post-SET extension of NSDs is essential for nucleosome binding in vitro.J. Biomol. Screen.19928–935. 10.1177/1087057114525854
2
Anderson-SchmidtH.BeltchevaO.BrandonM. D.ByrneE. M.DiehlE. J.DuncanL.et al (2013). Selected rapporteur summaries from the XX world congress of psychiatric genetics, Hamburg, Germany, october 14-18, 2012.Am. J. Med. Genet. Part B Neuropsychiatr. Genet.16296–121. 10.1002/ajmg.b.32132
3
AngelovaM. T.DimitrovaD. G.DingesN.LenceT.WorpenbergL.CarréC.et al (2018). The emerging field of epitranscriptomics in neurodevelopmental and neuronal disorders.Front. Bioeng. Biotechnol.6:1–15. 10.3389/fbioe.2018.00046
4
ArdlieK. G.DeLucaD. S.SegrèA. V.SullivanT. J.YoungT. R.GelfandE. T.et al (2015). The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans.Science348648–660. 10.1126/science.1262110
5
AytesA.GiacobbeA.MitrofanovaA.RuggeroK.CyrtaJ.ArriagaJ.et al (2018). NSD2 is a conserved driver of metastatic prostate cancer progression.Nat. Commun.97511–7514. 10.1038/s41467-018-07511-4
6
Barbosa-MoraisN. L.IrimiaM.PanQ.XiongH. Y.GueroussovS.LeeL. J.et al (2012). The evolutionary landscape of alternative splicing in vertebrate species.Science3381587–1593. 10.1126/science.1230612
7
BarrieE. S.AlfaroM. P.PfauR. B.GoffM. J.McBrideK. L.ManickamK.et al (2019). De novo loss-of-function variants in NSD2 (WHSC1) associate with a subset of Wolf–Hirschhorn syndrome.Cold Spring Harb. Mol. Case Stud.5:a004044. 10.1101/mcs.a004044
8
BlakeleyP.SiepenJ. A.LawlessC.HubbardS. J. (2010). Investigating protein isoforms via proteomics: A feasibility study.Proteomics101127–1140. 10.1002/pmic.200900445
9
BoczekN. J.LahnerC. A.NguyenT.-M.FerberM. J.HasadsriL.ThorlandE. C.et al (2018). Developmental delay and failure to thrive associated with a loss-of-function variant in WHSC1 (n.d.).Am. J. Med. Genet. Part A1762798–2802. 10.1002/ajmg.a.40498
10
C YuenR. K.MericoD.BookmanM.L HoweJ.ThiruvahindrapuramB.PatelR. V.et al (2017). Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder.Nat. Neurosci.20602–611. 10.1038/nn.4524
11
ChaudharyS.KhokharW.JabreI.ReddyA. S. N.ByrneL. J.WilsonC. M.et al (2019). Alternative splicing and protein diversity: Plants versus animals.Front. Plant Sci.10:00708. 10.3389/fpls.2019.00708
12
ChenL.-Y.ZhiZ.WangL.ZhaoY.-Y.DengM.LiuY.-H.et al (2019). NSD2 circular RNA promotes metastasis of colorectal cancer by targeting miR-199b-5p-mediated DDR1 and JAG1 signalling.J. Pathol.248103–115. 10.1002/path.5238
13
ClaytonE. A.RishishwarL.HuangT.-C.GulatiS.BanD.McDonaldJ. F.et al (2020). An atlas of transposable element-derived alternative splicing in cancer.Philos. Trans. R. Soc. B Biol. Sci.375:0342. 10.1098/rstb.2019.0342
14
CorominasR.YangX.LinG. N.KangS.ShenY.GhamsariL.et al (2014). Protein interaction network of alternatively spliced isoforms from brain links genetic risk factors for autism.Nat. Commun.5:4650. 10.1038/ncomms4650
15
CorrêaT.MergenerR.LeiteJ. C. L.GaleraM. F.MoreiraL. M. D. A.VargasJ. E.et al (2018). Cytogenomic Integrative Network Analysis of the Critical Region Associated with Wolf-Hirschhorn Syndrome.Biomed Res. Int.2018:5436187. 10.1155/2018/5436187
16
CorsettiE.AzpiazuN. (2013). Functional dissection of the splice variants of the Drosophila gene homothorax (hth).Dev. Biol.38472–82. 10.1016/j.ydbio.2013.09.018
17
DeardorffM. A.ZackaiE. H. (2007). Genetics and Metabolism Genetic Syndromes Caused by Chromosomal Abnormalities.Amsterdam: Elsevier Inc, 10.1016/B978-0-323-03004-5.50134-4
18
Deciphering Developmental Disorders Study. (2017). Prevalence and architecture of de novo mutations in developmental disorders.Nature542433–438. 10.1038/nature21062
19
DerarN.Al-HassnanZ. N.Al-OwainM.MoniesD.AbouelhodaM.MeyerB. F.et al (2019). De novo truncating variants in WHSC1 recapitulate the Wolf–Hirschhorn (4p16.3 microdeletion) syndrome phenotype.Genet. Med.21185–188. 10.1038/s41436-018-0014-8
20
DescartesM.KorfB. R.MikhailF. M. (2017). Chromosomes and Chromosomal Abnormalities. Sixth Edit. Amsterdam: Elsevier Inc., 10.1016/B978-0-323-37101-8.00035-7
21
DuY.HuH.HuaC.DuK.WeiT. (2018). Tissue distribution, subcellular localization, and enzymatic activity analysis of human SIRT5 isoforms.Biochem. Biophys. Res. Commun.503763–769. 10.1016/j.bbrc.2018.06.073
22
FranzM.RodriguezH.LopesC.ZuberiK.MontojoJ.BaderG. D.et al (2018). GeneMANIA update 2018.Nucleic Acids Res.46W60–W64. 10.1093/nar/gky311
23
GarcÃa-CarpizoV.SarmenteroJ.HanB.GrañaO.Ruiz-LlorenteS.PisanoD. G.et al (2016). NSD2 contributes to oncogenic RAS-driven transcription in lung cancer cells through long-range epigenetic activation.Sci. Rep.6:32952. 10.1038/srep32952
24
GongX.DelormeR.FauchereauF.DurandC. M.ChasteP.BetancurC.et al (2009). An investigation of ribosomal protein L10 gene in autism spectrum disorders.BMC Med. Genet.10:7. 10.1186/1471-2350-10-7
25
GordonD. E.HiattJ.BouhaddouM.RezeljV. V.UlfertsS.BrabergH.et al (2020). Comparative host-coronavirus protein interaction networks reveal pan-viral disease mechanisms.Science2020:abe9403. 10.1126/science.abe9403
26
GraveleyB. R. (2001). Alternative splicing: Increasing diversity in the proteomic world.Trends Genet.17100–107. 10.1016/S0168-9525(00)02176-4
27
HaladynaJ. N.YamauchiT.NeffT.BerntK. M. (2015). Epigenetic modifiers in normal and malignant hematopoiesis.Epigenomics7301–320. 10.2217/epi.14.88
28
HanX.PiaoL.YuanX.WangL.LiuZ.HeX. (2019). Knockdown of NSD2 suppresses renal cell carcinoma metastasis by inhibiting epithelial-mesenchymal transition.Int. J. Med. Sci.161404–1411. 10.7150/ijms.36128
29
HeC.LiuC.WangL.SunY.JiangY.HaoY. (2019). Histone methyltransferase NSD2 regulates apoptosis and chemosensitivity in osteosarcoma.Cell Death Dis.101347–1341. 10.1038/s41419-019-1347-1
30
HeyneH. O.SinghT.StambergerH.Abou JamraR.CaglayanH.CraiuD.et al (2018). De novo variants in neurodevelopmental disorders with epilepsy.Nat. Genet.501048–1053. 10.1038/s41588-018-0143-7
31
HomsyJ.ZaidiS.ShenY.WareJ. S.SamochaK. E.KarczewskiK. J.et al (2015). De novo mutations in congenital heart disease with neurodevelopmental and other congenital anomalies.Science3501262–1266. 10.1126/science.aac9396
32
HowriganD. P.RoseS. A.SamochaK. E.CerratoF.ChenW. J.ChurchhouseC.et al (2018). Schizophrenia risk conferred by protein-coding de novo mutations Daniel.bioRxiv131–22. 10.1101/495036
33
HuangX.LeDucR. D.FornelliL.SchunterA. J.BennettR. L.KelleherN. L.et al (2019). Defining the NSD2 interactome: PARP1 PARylation reduces NSD2 histone methyltransferase activity and impedes chromatin binding.J. Biol. Chem.29412459–12471. 10.1074/jbc.RA118.006159
34
HuangZ.WuH.ChuaiS.XuF.YanF.EnglundN.et al (2013). NSD2 Is recruited through Its PHD domain to oncogenic gene loci to drive multiple myeloma.Cancer Res.736277–6288. 10.1158/0008-5472.CAN-13-1000
35
IlouzR.Lev-RamV.BushongE. A.StilesT. L.Friedmann-MorvinskiD.DouglasC.et al (2017). Isoform-specific subcellular localization and function of protein kinase A identified by mosaic imaging of mouse brain.Elife6:17681. 10.7554/eLife.17681
36
JiangY.SunH.LinQ.WangZ.WangG.WangJ.et al (2019). De novo truncating variant in NSD2gene leading to atypical Wolf-Hirschhorn syndrome phenotype.BMC Med. Genet.20:134. 10.1186/s12881-019-0863-2
37
KatohM. (2016). Mutation spectra of histone methyltransferases with canonical SET domains and EZH2-targeted therapy.Epigenomics8285–305. 10.2217/epi.15.89
38
KeilJ. M.QaliehA.KwanK. Y. (2018). Brain Transcriptome Databases: A User’s Guide.J. Neurosci.382399–2412. 10.1523/jneurosci.1930-17.2018
39
KelemenO.ConvertiniP.ZhangZ.WenY.ShenM.FalaleevaM.et al (2013). Function of alternative splicing.Gene5141–30. 10.1016/j.gene.2012.07.083
40
KimJ.-H.LeeJ. H.LeeI.-S.LeeS. B.ChoK. S. (2017). Histone lysine methylation and neurodevelopmental disorders.Int. J. Mol. Sci.1818071404. 10.3390/ijms18071404
41
KimJ.-Y.HaeJ. K.ChoeN.-W.KimS.-M.EomG.-H.HeeJ. B.et al (2008). Multiple myeloma-related WHSC1/MMSET isoform RE-IIBP is a histone methyltransferase with transcriptional repression activity.Mol. Cell. Biol.282023–2034. 10.1128/MCB.02130-07
42
KimS.-M.KeeH.-J.ChoeN.KimJ.-Y.KookH.KookH.et al (2007). The histone methyltransferase activity of WHISTLE is important for the induction of apoptosis and HDAC1-mediated transcriptional repression.Exp. Cell Res.313975–983. 10.1016/j.yexcr.2006.12.007
43
KlauckS. M.FelderB.Kolb-KokocinskiA.SchusterC.ChiocchettiA.SchuppI.et al (2006). Mutations in the ribosomal protein gene RPL10 suggest a novel modulating disease mechanism for autism.Mol. Psychiatr.111073–1084. 10.1038/sj.mp.4001883
44
KobayashiN.HozumiY.ItoT.HosoyaT.KondoH.GotoK. (2007). Differential subcellular targeting and activity-dependent subcellular localization of diacylglycerol kinase isozymes in transfected cells.Eur. J. Cell Biol.86433–444. 10.1016/j.ejcb.2007.05.002
45
KrishnanA.ZhangR.YaoV.TheesfeldC. L.WongA. K.TadychA.et al (2016). Genome-wide prediction and functional characterization of the genetic basis of autism spectrum disorder.Nat. Neurosci.191454–1462. 10.1038/nn.4353
46
KuoA. J.CheungP.ChenK.ZeeB. M.KioiM.LauringJ.et al (2011). NSD2 links dimethylation of histone H3 at lysine 36 to oncogenic programming.Mol. Cell44609–620. 10.1016/j.molcel.2011.08.042
47
LelieveldS. H.ReijndersM. R. F.PfundtR.YntemaH. G.KamsteegE.-J.de VriesP.et al (2016). Meta-analysis of 2,104 trios provides support for 10 new genes for intellectual disability.Nat. Neurosci.191194–1196. 10.1038/nn.4352
48
LinG. N.CorominasR.NamH.UrrestiJ.IakouchevaL. M. (2017). Comprehensive Analyses of Tissue-Specific Networks with Implications to Psychiatric Diseases.Methods Mol. Biol.1613371–402. 10.1007/978-1-4939-7027-8_15
49
LinG. N.GuoS.TanX.WangW.QianW.SongW.et al (2019). PsyMuKB: An Integrative De Novo Variant Knowledge Base for Developmental Disorders.Genom. Proteom. Bioinform.2019:002. 10.1016/j.gpb.2019.10.002
50
LuckK.KimD.-K.LambourneL.SpirohnK.BeggB. E.BianW.et al (2020). A reference map of the human binary protein interactome.Nature580402–408. 10.1038/s41586-020-2188-x
51
LumishH. S.WynnJ.DevinskyO.ChungW. K. (2015). Brief Report: SETD2 Mutation in a Child with Autism, Intellectual Disabilities and Epilepsy.J. Autism Dev. Disord.453764–3770. 10.1007/s10803-015-2484-8
52
MarshallA. N.MontealegreM. C.Jiménez-LópezC.LorenzM. C.van HoofA. (2013). Alternative Splicing and Subfunctionalization Generates Functional Diversity in Fungal Proteomes.PLoS Genet.9:e1003376. 10.1371/journal.pgen.1003376
53
McDevittP. J.SchneckJ. L.DiazE.HouW.HuddlestonM. J.MaticoR. E.et al (2019). A Scalable Platform for Producing Recombinant Nucleosomes with Codified Histone Methyltransferase Substrate Preferences.Protein Expr. Purif.164:105455. 10.1016/j.pep.2019.105455
54
MerkinJ.RussellC.ChenP.BurgeC. B. (2012). Evolutionary dynamics of gene and isoform regulation in Mammalian tissues.Science3381593–1599. 10.1126/science.1228186
55
MessaoudiL.YangY.-G.KinomuraA.StavrevaD. A.YanG.Bortolin-CavailléM.-L.et al (2007). Subcellular distribution of human RDM1 protein isoforms and their nucleolar accumulation in response to heat shock and proteotoxic stress.Nucleic Acids Res.356571–6587. 10.1093/nar/gkm753
56
MirabellaF.MurisonA.AronsonL. I.WardellC. P.ThompsonA. J.HanrahanS. J.et al (2014). A novel functional role for MMSET in RNA processing based on the link between the REIIBP isoform and its interaction with the SMN complex.PLoS One9:0099493. 10.1371/journal.pone.0099493
57
ModrekB.LeeC. J. (2003). Alternative splicing in the human, mouse and rat genomes is associated with an increased frequency of exon creation and/or loss.Nat. Genet.34177–180. 10.1038/ng1159
58
MorishitaM.Di LuccioE. (2011). Cancers and the NSD family of histone lysine methyltransferases.Biochim. Biophys. Acta Rev. Cancer1816158–163. 10.1016/j.bbcan.2011.05.004
59
MoscaR.CéolA.SteinA.OlivellaR.AloyP. (2014). 3did: a catalog of domain-based interactions of known three-dimensional structure.Nucleic Acids Res.42D374–D379. 10.1093/nar/gkt887
60
NairR.RostB. (2009). Sequence conserved for subcellular localization.Protein Sci.112836–2847. 10.1110/ps.0207402
61
NarasimhanA.GreinerR.BatheO. F.BaracosV.DamarajuS. (2018). Differentially expressed alternatively spliced genes in skeletal muscle from cancer patients with cachexia.J. Cachexia. Sarcopenia Muscle960–70. 10.1002/jcsm.12235
62
NohH. J.PontingC. P.BouldingH. C.MeaderS.BetancurC.BuxbaumJ. D.et al (2013). Network topologies and convergent aetiologies arising from deletions and duplications observed in individuals with autism.PLoS Genet.9:e1003523. 10.1371/journal.pgen.1003523
63
O’LearyN. A.WrightM. W.BristerJ. R.CiufoS.HaddadD.McVeighR.et al (2016). Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation.Nucleic Acids Res.44D733–D745. 10.1093/nar/gkv1189
64
O’RoakB. J.VivesL.GirirajanS.KarakocE.KrummN.CoeB. P.et al (2012). Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations.Nature485246–250. 10.1038/nature10989
65
OudaR.SaraiN.NehruV.PatelM. C.DebrosseM.BachuM.et al (2018). SPT6 interacts with NSD2 and facilitates interferon-induced transcription.FEBS Lett.5921681–1692. 10.1002/1873-3468.13069
66
OughtredR.StarkC.BreitkreutzB.-J.RustJ.BoucherL.ChangC.et al (2019). The BioGRID interaction database: 2019 update.Nucleic Acids Res.47D529–D541. 10.1093/nar/gky1079
67
PanQ.ShaiO.LeeL. J.FreyB. J.BlencoweB. J. (2008). Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing.Nat. Genet.401413–1415. 10.1038/ng.259
68
ParkE.PanZ.ZhangZ.LinL.XingY. (2018). The Expanding Landscape of Alternative Splicing Variation in Human Populations.Am. J. Hum. Genet.10211–26. 10.1016/j.ajhg.2017.11.002
69
PoulinM. B.SchneckJ. L.MaticoR. E.McDevittP. J.HuddlestonM. J.HouW.et al (2016). Transition state for the NSD2-catalyzed methylation of histone H3 lysine 36.Proc. Natl. Acad. Sci. U. S. A.1131197–1201. 10.1073/pnas.1521036113
70
RichardH.SchulzM. H.SultanM.NürnbergerA.SchrinnerS.BalzereitD.et al (2010). Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments.Nucleic Acids Res.38e112–e112. 10.1093/nar/gkq041
71
SandersS. J.MurthaM. T.GuptaA. R.MurdochJ. D.RaubesonM. J.WillseyA. J.et al (2012). De novo mutations revealed by whole-exome sequencing are strongly associated with autism.Nature485237–241. 10.1038/nature10945
72
ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: a software environment for integrated models of biomolecular interaction networks.Genome Res.132498–2504. 10.1101/gr.1239303
73
SjaardaC. P.WoodS.McNaughtonA. J. M.TaylorS.HudsonM. L.LiuX.et al (2020). Exome sequencing identifies de novo splicing variant in XRCC6 in sporadic case of autism.J. Hum. Genet.65287–296. 10.1038/s10038-019-0707-0
74
StecI.WrightT. J.Van OmmenG.-J. B.De BoerP. A. J.Van HaeringenA.MoormanA. F. M.et al (1998). WHSC1, a 90 kb SET domain-containing gene, expressed in early development and homologous to a Drosophila dysmorphy gene maps in the Wolf-Hirschhorn syndrome critical region and is fused to IgH in t(4;14) multiple myeloma.Hum. Mol. Genet.71071–1082. 10.1093/hmg/7.7.1071
75
StessmanH. A. F.XiongB.CoeB. P.WangT.HoekzemaK.FenckovaM.et al (2017). Targeted sequencing identifies 91 neurodevelopmental-disorder risk genes with autism and developmental-disability biases.Nat. Genet.49515–526. 10.1038/ng.3792
76
TakumiT.TamadaK. (2018). CNV biology in neurodevelopmental disorders.Curr. Opin. Neurobiol.48183–192. 10.1016/j.conb.2017.12.004
77
TanakaH.IgataT.EtohK.KogaT.TakebayashiS.NakaoM. (2020). The NSD2/WHSC1/MMSET methyltransferase prevents cellular senescence−associated epigenomic remodeling.Aging Cell1913173. 10.1111/acel.13173
78
TsengY.-T.LiW.ChenC.-H.ZhangS.ChenJ. J. W.ZhouX. J.et al (2015). IIIDB: A database for isoform-isoform interactions and isoform network modules.BMC Genomics16:S10. 10.1186/1471-2164-16-S2-S10
79
UniProt Consortium. (2019). UniProt: a worldwide hub of protein knowledge.Nucleic Acids Res.47D506–D515. 10.1093/nar/gky1049
80
VeltmanJ. A.BrunnerH. G. (2012). De novo mutations in human genetic disease.Nat. Rev. Genet.13565–575. 10.1038/nrg3241
81
WangE. T.SandbergR.LuoS.KhrebtukovaI.ZhangL.MayrC.et al (2008). Alternative isoform regulation in human tissue transcriptomes.Nature456470–476. 10.1038/nature07509
82
YangX.Coulombe-HuntingtonJ.KangS.SheynkmanG. M.HaoT.RichardsonA.et al (2016). Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing.Cell164805–817. 10.1016/j.cell.2016.01.029
83
YuG.WangL.-G.HanY.HeQ.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters.Omi. A J. Integr. Biol.16284–287. 10.1089/omi.2011.0118
84
ZanniG.KalscheuerV. M.FriedrichA.BarresiS.AlfieriP.Di CapuaM.et al (2015). A Novel Mutation in RPL10 (Ribosomal Protein L10) Causes X-Linked Intellectual Disability, Cerebellar Hypoplasia, and Spondylo-Epiphyseal Dysplasia.Hum. Mutat.361155–1158. 10.1002/humu.22860
85
ZerbinoD. R.AchuthanP.AkanniW.AmodeM. R.BarrellD.BhaiJ.et al (2018). Ensembl 2018.Nucleic Acids Res.46D754–D761. 10.1093/nar/gkx1098
86
ZhangJ.JimaD.MoffittA. B.LiuQ.CzaderM.HsiE. D.et al (2014). The genomic landscape of mantle cell lymphoma is related to the epigenetically determined chromatin state of normal B cells.Blood1232988–2996. 10.1182/blood-2013-07-517177
87
ZhangW.Bojorquez-GomezA.VelezD. O.XuG.SanchezK. S.ShenJ. P.et al (2018). A global transcriptional network connecting noncoding mutations to changes in tumor gene expression.Nat. Genet.50613–620. 10.1038/s41588-018-0091-2
Summary
Keywords
NSD2, alternatively splicing, protein–protein interaction, isoform, RPL10
Citation
Wang W, Chen Y, Zhao J, Chen L, Song W, Li L and Lin GN (2021) Alternatively Splicing Interactomes Identify Novel Isoform-Specific Partners for NSD2. Front. Cell Dev. Biol. 9:612019. doi: 10.3389/fcell.2021.612019
Received
30 September 2020
Accepted
05 February 2021
Published
25 February 2021
Volume
9 - 2021
Edited by
Roland Wohlgemuth, Lodz University of Technology, Poland
Reviewed by
Kedryn K. Baskin, Wexner Medical Center, The Ohio State University, United States; Ellis Fok, The Chinese University of Hong Kong, China
Updates
Copyright
© 2021 Wang, Chen, Zhao, Chen, Song, Li and Lin.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Guan Ning Lin, nickgnlin@sjtu.edu.cnLi Li, lil@sjtu.edu.cn
†These authors have contributed equally to this work and share first authorship
This article was submitted to Molecular Medicine, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.