Abstract
Nuclear envelope (NE) proteins are essential for maintaining nuclear morphology, a central process controlling cellular phenotypes. While these proteins and structures are largely characterized in developmental biology, their function and associated abnormal nuclear morphology in disease, in particular cancer, remain elusive. We reported that throughout pancreatic transformation, oncogenic mutant KRAS (mKRAS) alters cancer cells nuclear size. Furthermore, we identified the NE protein Emerin as a mediator of these nuclear changes. Emerin’s function is versatile and dynamic depending on its protein interactions. Here, using a BioID approach we identified the NE protein, Nesprin-3, as a novel Emerin interactor in mKRAS pancreatic cancer cells. Transcriptomic and epigenomic analyses revealed that mKRAS increases Nesprin-3 expression using the transcription factor KLF5 as an effector. Functional studies demonstrated that Nesprin-3 phenocopies Emerin, where knockdown of Nesprin-3 impaired nuclear morphology, proliferation, and gene expression changes induced by mKRAS. Notably, Nesprin-3 overexpression can rescue the nuclear morphology and proliferative phenotype in pancreatic cancer cells with mKRAS knockdown. Finally, utilizing an Emerin mutant that cannot interact with Nesprin-3 we failed to rescue these nuclear size changes. Our data demonstrate a novel mechanism underlying the oncogenic function of KRAS, a major driver of this dismal condition.
Keywords: Nuclear morphology, KRAS, pancreatic cancer, Emerin, Nesprin-3
Introduction
Abnormal nuclear morphology has been a key feature used for cancer diagnosis, informing both disease staging and prognosis for close to 75 years.1,2 Further, these nuclear changes are responsible for altered gene expression associated with transformation of different malignancies including pancreatic cancer.3 However, the mechanisms and molecular players underlying the nuclear morphology changes remain for the most part uncharacterized. To date, the majority of research identifying and elucidating protein functions in nuclear morphology has been performed in normal physiological process (e.g., organogenesis) rather than disease contexts.4–6
Previous work from our lab has elucidated a novel role of the nuclear envelope protein Emerin in mediating nuclear morphology downstream of mutant KRAS in pancreatic ductal adenocarcinoma (PDAC),7 the most common histological subtype of pancreatic cancer.8 Specifically, we have shown that elevated Emerin levels are critical for the observed decrease in nuclear size in PDAC cells. Emerin is a relatively small protein, yet it has been shown to interact with many partners.9,10 It is hypothesized that Emerin’s diverse functions are largely regulated by its binding partners, with certain binding partners contributing to structural roles11–13 or genome organization and expression.9,14,15 Yet, Emerin’s protein interactome driving nuclear morphology changes in cancer remain uncharacterized. We define a novel mutant KRASG12D (mKRAS) induced interaction contributing to the nuclear morphology changes observed in PDAC cells. In the current model of PDAC development, tumors are driven by sequential genetic mutations,16 the first and driver mutation being KRAS (e.g., KRASG12D). We previously demonstrated that mKRAS drives altered nuclear morphology in PDAC and Emerin is essential for this observed phenomenon. While Emerin was identified to be important for oncogenic KRAS induced nuclear morphology changes in PDAC, the molecular mechanisms driving these changes remain unknown. Thus, we aim to characterize the mKRAS induced interactome contributing to nuclear morphology regulation in PDAC. Here, we report that oncogenic mKRAS is sufficient to induce a unique Emerin interactome, in PDAC cells, identified via a high throughput proximity-labeling assay. Analysis of nuclear envelope proteins implicated in nuclear architecture revealed significant interaction between Emerin and Nesprin-3 downstream of mKRAS. Nesprin-3 is a nuclear membrane protein that has been characterized in binding sun proteins to form the linker of nucleoskeleton and cytoskeleton (LINC) complex, which forms a physical linkage of the nucleus and cytoplasm via transmembrane proteins.17,18 Nuclear morphology studies revealed knockdown or overexpression of Nesprin-3 phenocopied Emerin’s impact on nuclear size. Further, using an Emerin mutant (P184T) that is unable to interact with Nesprin-3, we demonstrate that the overexpression of this mutant failed to induce nuclear morphology changes in the absence of mKRAS, highlighting a novel role of this interaction. Finally, transcriptomic and epigenomic profiling in the context of mKRAS revealed significant increases in Nesprin-3 levels and active histone marks at the Nesprin-3 locus. Transcription factor motif analysis identified Kruppel-like factor 5 (KLF5) as a potential regulator of Nesprin-3 expression. Overexpression of KLF5 resulted in increased Nesprin-3 levels and was sufficient to decrease nuclear size in the absence of oncogenic KRAS. Conversely, knockdown of Nesprin-3 impairs the proliferative and nuclear size changes induced by mKRAS PDAC cells. Together, our findings advance understanding of the mechanisms driving altered nuclear morphology in PDAC and identify Nesprin-3 as a novel Emerin interactor contributing to this cancer-associated phenotype.
Results
Mutant KRAS drives increased interaction between Emerin and Nesprin-3
Given that Emerin’s protein interactions are critical for its function, we sought to identify binding partners that may contribute to Emerin’s role in nuclear morphology observed in mKRAS expressing cells.7 To address this question and define the mKRAS driven Emerin interactome, we performed a BioID in the doxycycline (Dox) inducible KRASG12D 1012U cell line.19–21 Partner protein interactions in the absence (−Dox) and presence (+Dox) of mKRAS were quantified and measured as label free quantifications (LFQ). Analysis revealed 1008 Emerin binding partners influenced by mKRAS, of which 326 of these Emerin-protein interactions were increased (Figure 1A, Supplementary Table S1). From these data, we identified Nesprin-3 as a unique and enriched Emerin binding partner exclusive to mKRAS cells (Figure 1A). Furthermore, utilizing the BioID data, we specifically examined the interactions between Emerin and several nuclear and nuclear envelope proteins known to regulate nuclear morphology.22 Of these proteins, Eesprin-3 was the most upregulated in its interaction with Emerin (Figure 1B). While Emerin is a known interactor and component of the LINC complex,23 the interaction between Emerin and Nesprin-3 has not been previously reported. To see if this increased interaction was prevalent across other known LINC complex members, we examined the change in interaction between Emerin and the other Nesprin isoforms (Syne1/Syne2) and Sun proteins. Interestingly, the other Nesprin isoforms (Nesprin 1 and 2) were either not detected, or had no change in interaction with emerin under mKRAS conditions (Figure 1B, Syne1/Syne2, respectively). Similarly, both Sun1 and Sun2 experienced a decrease in their interaction frequency with Emerin upon mKRAS induction (Figure 1B). Taken together, these findings suggest a novel and unique interaction between Nesprin-3 and Emerin under mKRAS conditions. While our BioID identifies Emerin/Nesprin-3 as a mKRAS-induced interaction in pancreatic cancer cells, Nesprin-3 localization to the nuclear envelope is generally considered Sun/KASH-dependent,18 and our data do not test whether Emerin can substitute for Sun in this role. To validate the BioID findings, we performed an Emerin co-immunoprecipitation and confirmed Nesprin-3/Emerin interaction in mKRAS (+Dox) cells (Figure 2A). Further we examined this interaction using proximity ligation assays (PLA). These studies showed increased Emerin/Nesprin-3 interaction in the mKRAS-overexpressing cells compared to cells without KRAS overexpression, suggesting mKRAS enriches this interaction (Fig. 2B and C). As expected, single antibody controls did not show any positive PLA foci (Supplementary Figure S1). Finally, we sought to evaluate emerin/nesprin-3 colocalization. First, we characterized the nuclear distribution of Emerin and Nesprin-3 with immunofluorescence and confirmed abundant peripheral colocalization of Nesprin-3 and Emerin signal intensity with line-scan analyses (Fig. 3A and B). Quantification of Emerin and Nesprin-3 intensity at the nuclear periphery revealed significant peripheral enrichment in mKRAS (+Dox) expressing cells (Figure 3C). An increase in these proteins at the nuclear periphery supports the concept that altered localization of these proteins could be involved in changes in nuclear size in mKRAS-overexpressing cells. Nesprin-3 and Emerin colocalization and peripheral enrichment driven by mKRAS was confirmed in a secondary Dox-inducible mKRAS line, 4292F (Supplementary Figure S2). As expected, mKRAS expression resulted in enriched NE localization of both Emerin and Nesprin-3 (Supplementary Figure S2C). Overall, these results suggest that mKRAS drives this novel interaction between Nesprin-3 and Emerin at the nuclear periphery.
Figure 1.

Mutant KRAS drives a distinct Emerin interactome. (A) BioID data generated plot of proteins identified in N = 3 BioID assays in −Dox and +Dox 1012U cells and that had a non-zero value in one condition in at least two of the three experiments. To avoid division by zeroes, non-zero values were replaced by one percentile values from their respective experiments. (B) Heatmap of proteins involved in nuclear morphology regulation. Log2(FC) represents the fold change from −Dox to +Dox condition and is calculated as the average change observed in all three replicates.
Figure 2.

Emerin and Nesprin-3 interacts in mutant KRAS PDAC cells. (A) Representative Western blot analysis of Emerin immunoprecipitation shows co-precipitation of Nesprin-3, N = 3. Non-immune IgG was used as a negative control. (B) Representative immunofluorescence image of Proximity Ligation Assay (PLA) in 1012U − Dox and + Dox cells, Emerin and Nesprin-3 co-incubation generated PLA foci (orange spots) signifying proximity of these two proteins. Images shown as the max projection of the Z-stack. Scale bars = 10 mm. (C) Quantification of PLA images, calculating the number foci per nucleus for the specified conditions, N (# of nuclei) = 73 (−Dox Emerin + Nesprin-3), 110 (−Dox Emerin only), 79 (−Dox Nesprin-3 only), 116 (+Dox Emerin + Nesprin-3), 79 (+Dox Emerin only), 146 (+Dox Nesprin-3 only). Statistical significance between −Dox and + Dox Emerin/Nesprin-3 PLA calculated utilizing Mann–Whitney test.
Figure 3.

Mutant KRAS drives Nesprin-3 localization to the nuclear periphery in PDAC cells. (A) Representative immunofluorescence line-scan analysis of Emerin and Nesprin-3 in −Dox and +Dox treated 1012U cells. Dashed line depicts area of line-scan. Images taken from middle slice of Z-stack image. Scale bars = 10 μmm. (B) Profile plot of each signal in both conditions. Tracks normalized to a range of 100 to account for differences in nuclear size. N = 3, nuclei measured: −Dox 209, +Dox 179. (C) Percent of Emerin and Nesprin-3 (N3) signals within 3 μm of the nuclear periphery. Statistical significance calculated using Kruskal–Wallis test.
Mutant KRAS increases Nesprin-3 expression levels
To understand how mKRAS influences Nesprin-3 expression, we profiled the 1012U cell transcriptome and epigenome via RNA-seq and ChIP-seq, respectively, in KRAS WT and mutant conditions (−Dox vs +Dox cells). Analysis of the same nuclear proteins presented in Figure 1, revealed mKRAS influences gene expression of several nuclear regulators, of which Nesprin-3 shown under its gene name Syne3, was significantly increased (Figure 4A). We validated our initial RNA-seq data in 1012U cells by qRT-PCR for Nesprin-3. Our results indicate KRAS significantly upregulated RNA expression levels of Nesprin-3 (Figure 4B), which was confirmed in 4292F cells (Supplementary Figure S3A). Increased RNA expression was accompanied by increased protein expression in both lines (Figure 4C, Supplementary Figure S3B). Consistent with our prior observation that Emerin levels also rise upon KRAS expression,7 these findings suggest coordinated upregulation of both proteins, along with their increased interaction under mKRAS induction.
Figure 4.

Mutant KRAS increases Nesprin-3 expression levels and is regulated by KLF5. (A) Heatmap depicting Log2 fold change in nuclear envelope proteins from RNA sequencing of 1012U cells. Syne3 corresponds to Nesprin-3. Data are shown as change from −Dox to +Dox, N = 3. (B) qRT-PCR for Nesprin-3 mRNA in 1012U cells, N = 3. Statistical significance calculated utilizing Mann–Whitney test. (C) Representative Western blot analysis of Nesprin-3 from 1012U cells, N = 4. Graph shows Nesprin-3 protein levels relative to Vinculin from densitometry analysis. Statistical significance calculated utilizing Mann–Whitney test. (D) ChIP-seq tracks of H3K27Ac, H3K4me1, and K3K4me3 over Nesprin-3 in −Dox and +Dox 1012U treated cells, N = 3. Orange box indicated Nesprin-3 promoter. (E) ChIP-seq tracks of KLF5 in human PDAC cell lines over Nesprin-3 gene, N = 3. Orange box indicates Nesprin-3 promoter. (F) Dot plot showing Nesprin-3 Log2 fold change following KLF5 knock-down from RNA-seq in human PDAC cell lines, N = 3. (G) qRT-PCR for KLF5 and Nesprin-3 following siRNA-mediated KLF5 knockdown in L3.6pl cells, N = 4. Representative Western blot analysis showing KLF5 knockdown. Statistical significance calculated utilizing Mann–Whitney test.
Epigenomic analysis of active histone marks (H3K27Ac, H3K4me1, H3K4me3)24 showed enrichment of these modifications across the Nesprin-3 promoter in mKRAS cells (Figure 4D). Transcription factor motif analysis using the JASPAR algorithm25 identified the specificity protein (SP) and Kruppel-like factors (KLF) families of transcription factors consensus binding sites in the Nesprin-3 promoter. Using RNA-seq from 1012U cells revealed several of these transcription factors are upregulated under mKRAS, including SP5, KLF5, and KLF15 having the greatest increase (Supplementary Figure S4A). Given that SP5 and KLF15 have been shown to be repressors,26 we selected KLF5 as the target transcription factor driving Nesprin-3 expression. Increased KLF5 gene and protein expression under mKRAS was confirmed via qRT-PCR and Western blotting in 1012U and 4292F cells (Supplementary Figure S3C to E). We cross-referenced with RNA-seq from human PDAC lines, AsPC1 and L3.6pl,27 and found KLF5 was elevated in these cells compared to the other members of the SP/KLF family (Supplementary Figure S4B). Furthermore, analysis of PDAC cases, generated by the National Cancer Institute Clinical Proteomic Tumor Analysis Consotrtium (CPTAC), revealed that KLF5 expression was increased in mKRAS tumors compared to KRAS wild-type tumors (Supplementary Figure S3F).28 ChIP-seq data showed KLF5 binding to the Nesprin-3 promoter in these mKRAS PDAC cells confirming that Nesprin-3 is direct target of KLF5 (Figure 4E).
To validate KLF5 mediated regulation of Nesprin-3, we utilized RNA-seq from human PDAC cells with knockdown of KLF5. Importantly, KLF5 depletion in mKRAS lines results in significantly lower Nesprin-3 RNA levels, whereas depletion in KRAS WT cells (BxPC3) did not affect Nesprin-3 expression (Figure 4F). We confirmed these results in L3.6pl using two independent targeting siRNAs (Figure 4G, Supplementary Figure S5A). In alignment with these findings, overexpression of KLF5 in cells with knockdown of mKRAS (iPanc-1) demonstrated that KLF5 can increase Nesprin-3 levels in the absence of mKRAS (Supplementary Figure S5B). Taken together, these results suggest that mKRAS uses the transcription factor, KLF5 as an effector to regulate Nesprin-3 expression in PDAC cells.
Nesprin-3 phenocopies Emerin’s effect on nuclear size
To understand the function and importance of this novel Emerin–Nesprin-3 interaction in nuclear morphology, we performed Nesprin-3 modulation experiments. First, to implicate Nesprin-3 in the nuclear morphology changes observed in mKRAS cells, we analyzed immunofluorescence experiments in mKRAS 1012U cells with siRNA mediated Nesprin-3 knockdown. Like Emerin,7 loss of Nesprin-3 results in increased nuclear surface area and volume (Figure 5A and B). Specifically, depletion of Nesprin-3 significantly increased the median nuclear surface area and volume (Figure 5C). Enlarged nuclear size in the presence of siRNA mediated knockdown of Nesprin-3 was confirmed using a second independent siRNA (Supplementary Figure S6), and in the mKRAS 4292F cells (Supplementary Figure S7). Similarly, to confirm Nesprin-3’s role in nuclear size changes, we overexpressed Nesprin-3 in 1012U WT cells (−Dox) (Figure 5D and E). Interestingly, Nesprin-3 overexpression phenocopies the mKRAS induced nuclear morphology changes, resulting in significantly smaller nuclei (Figure 5F).
Figure 5.

Nesprin-3 phenocopies Emerin-induced nuclear changes. (A) Representative immunofluorescence images of Nesprin-3 knockdown (siN3) in 1012U +Dox treated cells. Images shown as the max projection of the Z-stack. Scale bars = 10 mm. (B) Representative Western blot analysis of Nesprin-3 following siRNA mediated knockdown, N = 3. Graph shows relative Nesprin-3 levels relative to Vinculin from densitometry analysis. Statistical significance calculated using an unpaired t test. (C) Quantification of nuclei from NT or Nesprin-3 siRNA transfected cells, median surface area = 499.1 (siNT); 551.6 (siN3), median volume = 262.1 (siNT); 331.4 (siN3), N = 3. Statistical significance calculated utilizing the Mann–Whitney test. (D) Representative immunofluorescence images of Nesprin-3 overexpression (OE) in 1012U −Dox cells. Images shown as the max projection of the Z-stack. Scale bars = 10 μm. (E) Representative Western blot analysis showing Nesprin-3 levels from overexpression studies, N = 4. Graph shows Nesprin-3 levels relative to Vinculin. Statistics calculated using Mann–Whitney test. (F) Quantification of nuclei from EV or Nesprin-3 transfected cells, median surface area = 541.1, n = 567 nuclei (EV); 511.2, n = 594 nuclei (OE), median volume = 284.7 (EV); 222.2 (OE), N = 4. Statistical significance calculated using Mann–Whitney test.
To confirm Nesprin-3’s role is consistent in human models of PDAC, these studies were performed in two independent cell lines. First, siRNA mediated depletion of Nesprin-3 in the mKRAS human PDAC cell line, L3.6pl, resulted in a significant increase in nuclear size (Supplementary Figure S8A to D). To complement the mouse inducible mKRAS cells, Nesprin-3’s role under mKRAS was assessed in a doxycycline inducible shRNA targeting KRAS human cell line, iPanc-1 (Supplementary Figure S8E). As expected, silencing of mKRAS increased nuclear surface area and volume (Supplementary Figure S8G). Interestingly, nuclear morphology in the absence of mKRAS (+Dox) could be rescued with Nesprin-3 overexpression (Supplementary Figure S8E to G). Taken together these studies highlight the role of Nesprin-3 in mediating nuclear size in several complementary mKRAS models.
Having shown that Nesprin-3 levels are regulated by KLF5, we next aimed to investigate if these nuclear morphology changes can be induced upstream of Nesprin-3 expression. First, Nesprin-3 expression was examined in the presence of KLF5 overexpression in KRAS WT cells (−Dox). KLF5 overexpression upregulated Nesprin-3 at both RNA and protein levels (Supplementary Figure S9A to C). To examine the corresponding effect on nuclear size, nuclei were quantified following KLF5 overexpression (Supplementary Figure S9D). In immunofluorescence experiments, nuclear surface area and volume both decreased in cells overexpressing this transcription factor (Fig. S9E). Taken together, these data suggest that the KLF5-Nesprin-3 axis is both necessary and sufficient for regulating nuclear morphology and modulating its expression can phenocopy mKRAS conditions.
Nesprin-3 is required for mutant KRAS oncogenic functions
Given the nuclear morphology results, we aimed to determine if these changes induced by KRAS via Nesprin-3 are associated with genome re-organization using ChIP-seq in the 1012U cells for H3K4me1 and H3K27ac, chromatin post-translational modifications marking active enhancers. These are distal regulatory elements controlling promoter activity and gene expression, thus representing good indicators of global re-organization of the genome. To understand how mKRAS drives changes to the enhancer landscape, we performed GREAT analysis29 on H3K4me1 and H3K27Ac peaks unique to WT KRAS (1012U −Dox) and mKRAS (1012U +dox) conditions. Specifically, GREAT analysis identified 67,164 regulatory elements unique to WT KRAS cells and 63,962 elements unique to mKRAS cells (Supplementary Table S2). When mapping these regions to the genome, we identified subsets of enhancer associated genes unique to WT and mKRAS conditions (Supplementary Figure S10A). To characterize the functions of the mKRAS enhancer associated genes we performed KEGG pathway analysis (Supplementary Table S3).30 As anticipated, these genes were associated with pathways typically upregulated in cancer, such as MAPK signaling, Ras signaling and axon guidance among others (Supplementary Figure S10B).31–35 We selected a subset of genes associated with our enriched pathways and validated differential H3K4me1 and H3K27Ac binding. Histone mark binding was assessed over a 5–30 kb region around the transcription start site to capture the majority of peaks mapped to the gene identified from GREAT analysis. Analysis of H3K4me1 peaks over this region shows a strong increase under mKRAS (+Dox) for most genes examined (Supplementary Figure S10C). Similar results were seen for H3K27Ac peaks (Supplementary Figure S10C), suggesting mKRAS driven enrichment of active enhancer marks at these genes. Examination of RNA-seq revealed that the enrichment of active histone marks is associated with a corresponding increase in gene expression under mKRAS (Supplementary Figure S10D). To determine the requirement of Nesprin-3 in the regulation of these genes downstream of mKRAS, we performed siRNA-mediation knockdown of Nesprin-3 in mKRAS cells (1012U +Dox). Importantly, loss of Nesprin-3 resulted in decreased RNA expression of these genes (Supplementary Figure S10E). Together, these results suggest mKRAS induced nuclear morphology changes are associated with global genome changes, and these gene expression changes can be reversed with the removal of Nesprin-3.
Next, to determine if Nesprin-3 could regulate KRAS’s oncogenic functions, we performed proliferation assays with Nesprin-3 modulation on three different mKRAS cell lines using two independent methods, CyQUANT and MTT. We found that, upon Nesprin-3 knockdown in both human and mouse cell lines (L3.6pl and 1012U +Dox) there is a statistically significant decrease in proliferation at 72 h compared to the baseline (0 h) (Supplementary Figure S11A and B). Western blot analyses were performed and confirmed protein reduction of Nesprin-3 at the 72-h time point (Supplementary Figure S11A and B). Complementarily, we performed proliferation assays in the iPanc-1 cell line. As expected, silencing of mKRAS decreases cell proliferation, which can be rescued with Nesprin-3 overexpression (Supplementary Figure S11C). Overexpression of Nesprin-3 was confirmed using Western blotting and KRAS silencing was confirmed using qRT-PCR at the 72 h time point (Supplementary Figure S11D and E). All together, these results confirm this new role of Nesprin-3 as an effector of oncogenic KRAS in pancreatic cancer cells.
Interaction between Emerin and Nesprin-3 is responsible for nuclear size regulation
Our previous work and the above presented data suggest both Emerin and Nesprin-3 are important players regulating nuclear morphology in PDAC. However, it remains unknown if the interaction between these two proteins is responsible for this effect. Emerin has previously been reported to interact with other Nesprin isoforms: Nesprin-1 and −2,36–38 and given that Nesprin-3’s structure is highly similar to these isoforms17,39,40 it is plausible similar protein domains are regulating this interaction. Previous work defining human Emerin interaction with the other Nesprin isoforms identified a point mutation (P183T) in Emerin disrupts the interaction with Nesprin-1.36 Thus, we first sought to determine if this residue was responsible for the interaction with Nesprin-3. Since we are predominantly utilizing mouse systems, we generated a mouse homolog of the P183T mutation (P184T). Nuclear envelope localization and colocalization with Nesprin-3 was confirmed using immunofluorescence (Supplementary Figure S12A). To determine the extent of protein interaction with Nesprin-3, HA-tagged Emerin WT or P184T mutant were overexpressed in mKRAS 1012U cells (+Dox) and PLAs were performed utilizing the same antibodies used to test colocalization (Supplementary Figure S12B). Analysis of foci localization revealed predominantly nuclear foci, with some in the perinuclear space. For this reason, only nuclear foci were utilized for analysis. Despite significant expression of the Emerin P184T mutant (Supplementary Figure S12C), only the Emerin WT had significant interaction with Nesprin-3 compared to the empty vector control (Supplementary Figure S12D). Our results map for the first time the interaction between Emerin and Nesprin-3, identifying the key residue in Emerin responsible for this interaction.
Having defined the residue crucial for Emerin-Nesprin-3 interaction, we next wanted to investigate how disrupting this interaction will influence nuclear size. Either WT or P184T mutant Emerin was transfected into 1012U −Dox cells (Supplementary Figure S12E), and correct localization in the absence of mKRAS at the nuclear periphery was confirmed using immunofluorescence (Figure 6A). Intriguingly, overexpression of the P184T mutant did not result in decreases to nuclear surface area or volume compared to the EV control (Figure 6B). Importantly, overexpression of WT Emerin decreased both nuclear surface area and volume, as expected (Figure 6B). Taken together, these data highlight that by preventing Emerin–Nesprin-3 interaction we can prevent changes in nuclear morphology. Thus, our findings improve the understanding of molecular players and the mechanisms underlying nuclear morphology regulation in PDAC.
Figure 6.

Disruption of the interaction with Nesprin-3 impairs Emerin-induced nuclear size changes. (A) Representative immunofluorescence images using the HA antibody from 1012U −Dox cells transfected with EV, HA-tagged wild-type Emerin, or HA-tagged P184T mutant Emerin, N = 3. Images shown as single slice from center of Z-stack image. Scale bars = 10 μm. (B) Quantification of nuclei from EV or Emerin transfected cells, median surface area = 615.0, n = 249 nuclei (EV); 533.1, n = 230 nuclei (WT); 619.9, n = 219 nuclei (P184T), median volume = 323.8 (EV); 231.8 (WT); 305.0 (P184T), N = 3. Statistical significance calculated utilizing Kruskal–Wallis test.
Discussion
Abnormal nuclear morphology is a hallmark of PDAC, yet the molecular mechanisms driving these changes remain poorly defined. Here, we identify a mKRAS-induced Emerin interactome, highlighting a novel role for the Emerin–Nesprin-3 interaction in nuclear morphology regulation. Through high-throughput protein interaction profiling we have identified a mKRAS induced Emerin interactome. While this allows us to investigate the role of nuclear and nuclear envelope proteins in mediating nuclear morphology, it also revealed mKRAS is sufficient to create a unique protein environment within the nucleus, highlighting the need for continued investigation into how KRAS-driven nuclear envelope remodeling contributes to PDAC progression. Further, using a combination of molecular and cell biological techniques, we demonstrated an interaction between Emerin and Nesprin-3. This is the first time to our knowledge that Nesprin-3 has been reported to interact with Emerin. Importantly, we demonstrated this interaction functions in mediating the mKRAS-Emerin axis induced effects on nuclear size.
Abnormal nuclear morphology in PDAC is increasingly recognized not only as a diagnostic readout but also as a functional driver of tumor behavior, influencing chromatin organization and mechanotransduction signaling. This underscores the significance of our findings by linking KRAS-driven nuclear envelope remodeling to processes that may actively promote PDAC development. Of all normal tissues, Nesprin-3 mRNA has been reported to be the lowest in salivary glands and the pancreas.41 Consequently, Nesprin-3 has not been previously linked with PDAC nor regulation by KRAS. Here, we demonstrated Nesprin-3 expression and localization are regulated by mKRAS. We have shown that mKRAS increases expression of the candidate Nesprin-3 transcription factor KLF5, and KLF5 modulation influences Nesprin-3 expression and is sufficient to mimic mKRAS effects on nuclear size, making Nesprin-3 a novel target of KLF5. Thus, expanding the repertoire of the cellular functions controlled by the SP/KLF family of transcription factors, involving nuclear envelope proteins regulating nuclear dynamics in PDAC. Importantly, analysis of Nesprin-3 expression in PDAC reveals a subset of tumors express higher levels of Nesprin-3, visualized via immunohistochemistry,42 whereas in the normal pancreas Nesprin-3 is not detectable (Supplementary Figure S13A). Further, the analysis also showed higher levels of Nesprin-3 in tumor tissue compared to the normal adjacent tissues (Supplementary Figure S13C). Finally, we demonstrate that high Nesprin-3 expression is associated with decreased nuclear size (Supplementary Figure S13B) and overall survival (Supplementary Figure S13D), highlighting an important need to elucidate Nesprin-3’s contribution to malignant features of PDAC.
In addition, Nesprin-3 knockdown phenocopies Emerin knockdown, proposing a synergistic relationship with Emerin in regulating nuclear size. Together, these provide evidence that Nesprin-3 is an important effector downstream of mKRAS in PDAC. In support of a functional consequence, Nesprin-3 has also been shown to act as a “nuclear piston” that generates compartmentalized intracellular pressure to drive lobopodial 3D migration.43 This piston-like role suggests that KRAS-induced upregulation of Nesprin-3 could couple nuclear morphology changes to altered nuclear mechanics and cell motility, processes highly relevant to PDAC invasion and metastasis. A previous study demonstrated that when Nesprin-2 mini was co-expressed with Nesprin-3 nuclear size reduction was enhanced in human keratinocytes.44 These data provide further evidence for nuclear size regulation by Nesprin-3.
Whereas our study has shown the importance for Emerin and Nesprin-3 interaction in the observed nuclear morphology, the exact molecular mechanism driving this interaction remains unknown, and can serve as the focus for future studies to better characterize how nuclear morphology is contributing to PDAC pathobiology. Expanding this field and these studies will continue to improve our understanding of the disease and lead to the potential for new biomarkers and therapeutically targetable proteins to improve prognosis of this dismal disease.
Materials and Methods
Plasmids
BioID2-only and BioID2-LaA pBabe-puro constructs were previously generated .21 N144-BioID2 pBabe-puro construct was previously generated.45 Emerin was amplified from human cDNA by PCR and inserted into the BioID2-only pBabe puro plasmid by restriction enzyme digest with XhoI/SalI to make BioID2-Emerin pBabe-puro. pcDNA-FLAG-KLF5 plasmid was obtained from Addgene (RRID: 40918).46 Vector pcDNA-HA2 was a gift from Scott Kaufman and was opened using HindIII and BamHI restriction sites. Wild-type emerin was amplified from 1012 U (−Dox) cDNA using primers GIB-Ms-Emd-F1 5′-tagggagacccaagcttatggacgactatgcggttttg/GIB-Ms-Emd-R1 5′-ccagctcctggcgcggatccgggatccatccagaaggggt. Gibson assembly (E2611S; NEB) was done to ligate the opened vector (5,401 bp) and fragment (814 bp) to make the final pcDNA-ms-Emd-HA2 plasmid. Gibson assembly was again utilized when creating Emerin P184T plasmid. Briefly, fragments were amplified from wild-type emerin plasmid (pcDNA-ms-Emd-HA2) using primer combinations; GIB-Ms-Emd-F2 5′-ctatagggagacccaagcttatgg/GIB-Ms-Emd-P184T-R1 5′-agaagaggtggaggatgtagTataatatgat (590 bp) and GIB-Ms-Emd-P184T-F1 5′-tgggcttatcatattatActacatcctcc/GIB-Ms-Emd-R2 5′-ACACTATAGAATAgggccctctag (351 bp). Vector and fragments were ligated together using Gibson assembly to make the final pcDNA-ms-Emd-P184T-HA2 plasmid.
Cell lines, cell culture and transfections
Standard incubation conditions of 37 °C and 5% CO2 were used in all experiments. 1012U and 4292F cells were obtained from Dr Marina Pasca di Magliano (University of Michigan, Ann Arbor, MI, USA).19 iKRAS cells were cultured in RPMI 1640 (10–040-CV; Corning) with 10% Fetal Bovine Serum (35–010-CV, Corning) and doxycycline hyclate (D9891; Sigma-Aldrich) was used at 1 mg/mL to induce mKRAS. L3.6pl cells were cultured in MEM (10–010-CV; Corning) with 10% fetal bovine serum, 1X vitamins (11120052; Thermo Fisher Scientific), 1% L-glutamate (25030081; Thermo Fisher Scientific), 1% sodium pyruvate (25–000Cl; Corning) and 1% MEM NEAA (11140–050; Gibco). Doxycycline inducible shKRAS PANC-1 cells were cultured in DMEM (10–013-CV; Corning) with 10% fetal bovine serum as previously described.47 siRNA transfections were performed with RNAiMAX (13778500; Invitrogen) or Dharmafect (T-2001–03; Horizon Discovery) per manufacturer’s suggestions. siRNAs used in these experiments were Syne3 FlexiTube (SI00835373; Qiagen), Syne3 SMARTpool (L-052180–01-0005; Dharmacon), KLF5 FlexiTube (SI04267641; Qiagen), KLF5 SMARTpool (L-013571–00-0005; Dharmacon), NT FlexiTube (1027281; Qiagen) and NT SMARTpool (D-001810–10-20; Dharmacon). Plasmid transfections were performed with X-Treme Gene (06366546001; Sigma-Aldrich). Briefly, cells were seeded at appropriate concentrations to ensure 60% confluency the next day. Cells were transfected with transfection complex in Opti-MEM (Gibco; cat. 31985062), X-treme Gene and 5 mg DNA and incubated for 24 h. Transfected cells were replated for desired assays and collected 24 or 48 h later for overexpression or knockdown assays, respectively.
BioID stable cell line generation
Stable cell lines for all 1012U iKRAS BioID lines were generated using retroviral transduction. HEK293 Phoenix cells (National Gene Vector Biorepository) were transfected with each construct using Lipofectamine 3000 (Thermo Fisher Scientific) per the manufacturer’s recommendation. Transfected cells were incubated at 37 °C for 6 h, then replenished with fresh medium and incubated at 32 °C for 72 h for viral production. The viral culture media was filtered through a 0.45-μm filter and added to 1012U iKRAS cells along with doxycycline (Dox; 1 μg/mL) and polybrene (4 μg/mL; Santa Cruz Biotechnology). At 72 h after transduction, cells were replenished with fresh media for 24 h. Puromycin (0.5 μg/mL; Thermo Fisher Scientific) was then added for 5 days to select expressing cells. Fusion protein expression and biotinylation were verified using IF and WB. Cells were grown with dox (1 mg/mL) for stable cell line maintenance. For comparative BioID analysis of −Dox and +Dox 1012U iKRAS BioID cell lines, cells were seeded without dox in biotin-free RPMI 1640 (MBS653376; MyBioSource) for 3 days and then either not induced (−Dox) or induced (+Dox, 1 mg/mL) for 3 more days prior to supplementation with 50 μM biotin for 18 h.
BioID
BioID was performed as previously described.48 Briefly, four 10 cm dishes of cells per sample were lysed in 8 M urea 50 mM Tris pH 7.4 containing protease inhibitor (87785; Thermo Fisher Scientific) and DTT, incubated with universal nuclease (88700; Thermo Fisher Scientific), and sonicated to further shear DNA. Gelatin Sepharose 4B beads (17095601; GE Healthcare) were used to preclear samples for 2 h prior to incubation with Streptavidin Sepharose High Performance beads (17511301; GE Healthcare) for 4 h. Streptavidin beads were washed four times with 8 M urea 50 mM Tris pH 7.4 wash buffer and resuspended in 50 mM ammonium bicarbonate with 1 mM biotin.
Mass spectrometry analysis
After immunoprecipitation, proteins were digested directly on beads. Briefly, proteins bound to the beads were resuspended with 8 M urea and 50 mM ammonium bicarbonate, and cysteine disulfide bonds were reduced with 10 mM tris(2-carboxyethyl) phosphine at 30 °C for 60 min, followed by cysteine alkylation with 30 mM iodoacetamide in the dark at room temperature for 30 min. After alkylation, urea was diluted to 1 M urea using 50 mM ammonium bicarbonate, and proteins were finally subjected to overnight digestion with MS-grade trypsin/Lys-C mix (Promega). Finally, beads were pulled down and the solution with peptides collected into a new tube. The beads were then washed once with 50 mM ammonium bicarbonate to increase peptide recovery. The digested samples were desalted using a C18 TopTip (PolyLC), and the organic solvent was removed in a SpeedVac concentrator before liquid chromatography (LC)–MS/MS analysis. Dried samples were reconstituted in 100 mM ammonium formate, pH ~10, and analyzed by 2D LC-MS/MS using a 2D nanoACQUITY Ultra Performance Liquid Chromatography system (Waters, Milford, MA) coupled to a Q-Exactive Plus mass spectrometer (Thermo Fisher Scientific). Peptides were loaded onto the first-dimension column, XBridge BEH130 C18 NanoEase (300 μm 50 mm, 5 μm) equilibrated with solvent A (20 mM ammonium formate, pH 10, first-dimension pump) at 2 μL/min. The first fraction was eluted from the first-dimension column at 17% of solvent B (100% acetonitrile) for 4 min and transferred to the second-dimension Symmetry C18 trap column of 0.180 × 20 mm (Waters) using a 1:10 dilution with 99.9% second-dimension pump solvent A (0.1% formic acid in water) at 20 μL/min. Peptides were then eluted from the trap column and resolved on an analytical C18 BEH130 PicoChip column of 0.075 × 100 mm, 1.7-mm particles (NewObjective, MA), at low pH by increasing the composition of solvent B (100% acetonitrile) concentrations of solvent B. The following first-dimension fractions were eluted at 19.5, 22, 26, and 65% solvent B. The mass spectrometer was operated in positive data-dependent acquisition mode. MS1 spectra were measured with a resolution of 70,000, an AGC target of 1e6, and a mass range from 350 to 1700 m/z. Up to 12 MS2 spectra per duty cycle were triggered, fragmented by high-energy collision–induced dissociation, and acquired with a resolution of 17,500 and an automatic gain control (AGC) target of 5e4, an isolation window of 2.0 m/z, and a normalized collision energy of 25. Dynamic exclusion was enabled with duration of 20 s. All mass spectra were analyzed with MaxQuant software, v. 1.5.5.1. MS/MS spectra were searched against the Mus musculus Uniprot protein sequence database (v. July 2016) and Global Proteome Machine Common Repository of Adventitious Proteins sequences (commonly known protein contaminants). Precursor mass tolerance was set to 20 and 4.5 ppm for the first search where initial mass recalibration was completed and for the main search, respectively. Product ions were searched with a mass tolerance of 0.5 Da. The maximum precursor ion charge state used for searching was seven. Carbamidomethylation of cysteines was searched as a fixed modification, and oxidation of methionines and acetylation of protein N-terminal were searched as variable modifications. Enzyme was set to trypsin in a specific mode, and a maximum of two missed cleavages was allowed for searching. The target-decoy-based false discovery rate filter for spectrum and protein identification was set to 1%. Proteins were classified as candidate interactors if LFQ intensities were at least 3-fold greater compared to BioID2-only control.
Immunoprecipitation
Emerin immunoprecipitation was performed by resuspending 1012U cell pellet in 1 mL of lysis buffer containing 150 mM NaCl, 50 mM Tris ph 8, 0.5% NP-40, 10% glycerol, and complete protease inhibitor (11697498001; Roche) and sonicated for 3 s at 10% frequency for 5 cycles with 3 s in between each cycle. Lysate was centrifuged at 2,000 rpm at 4 °C for 5 min and supernatant was collected for immunoprecipitations. Input was collected, and pellet was used to identify Emerin lost in insoluble fraction. Seven microliters of Emerin antibody or IgG control was added to the lysate and incubated at 4 °C with gentle inversion overnight. Lysate was then incubated with magnetic beads (Dynabeads, 10009D; Invitrogen) for 1 h at room temperature. Magnetic beads were isolated using a magnetic microcentrifuge tube rack and washed three times with washing buffer containing 300 mM NaCl, 50 mM Tris pH 8, 0.5% NP-40. Protein was then mixed with 1× Laemmli buffer and heated at 95 °C for 5 min.
Proximity ligation assay (PLA)
1012U cells were seeded at 30% confluency on round 25 mm coverslips (174985; Thermo Fisher Scientific) in a six-well plate and cultured with or without doxycycline at 37 °C for 72 h. Cells were then washed 3× with 1 × PBS and fixed for 15 min with 3.2% (wt/vol) paraformaldehyde/1× phosphate-buffered saline (PBS) at room temperature, followed by permeabilization in 0.1% Triton X-100 in PBS for 10 min on ice. Coverslips were blocked with Duolink blocking buffer (DUO82007–4mL; Sigma-Aldrich) for 1 h at room temperature. PLA was performed using the Duolink In Situ Orange Starter Kit Mouse/Rabbit (DUO92102; Sigma). Primary antibodies (Emerin/Nesprin-3 or anti-HA/Nesprin-3) were incubated for 1 h at room temperature. PLA was performed according to the manufacturer’s specifications except for amplification times: 300 min for Emerin/Nesprin-3; 60 min for anti-HA/Nesprin-3. Cover slips were mounted in Duolink DAPI (DUO82040–5ML; Sigma-Aldrich). For negative controls, coverslips were incubated with a single primary antibody and both mouse and rabbit PLA secondary probes were used during processing for PLA. Microscopy was performed using a Zeiss LSM800 with Zeiss AxioObserver Z1 microscope and C-Aprochromat 63×/1.20 W objective. Z stacks were taken at 1 mm slices with a scan speed of 9, and 8-bit depth. Diode lasers used were 561 and 405 nm and detected using GaASP-PMT.
Western blot analysis
Cells were lysed in RIPA buffer (50 mM Tris-HCl pH 7.5, 1% NP-40, 150 mM NaCl, 5 mM EDTA, 0.25% NaDOC, and 1.0% SDS) with complete protease inhibitor, sonicated, and centrifuged. Lysates were collected and quantified using a BCA kit (23225; Thermo Fisher Scientific). 20–50 mg of protein were loaded into 4–15% gel (4561086; Bio-Rad) and resolved via SDS-PAGE and transferred to PVDF membranes activated in methanol (IPVH00010; Millipore). Membranes were blocked in either 5% nonfat dry milk (20–241; Apex Chemicals and Reagents) or 3% BSA (9048–46-8; RPI) in 1× PBST for 1 h. Primary antibodies were diluted in blocking buffer and incubated overnight at 4 °C. Membranes were washed three times in 1× PBST and incubated with secondary antibodies for 1 h at room temperature. Primary and secondary antibodies are reported below (Table 1). Blots were developed using SuperSignal West Pico Chemiluminescent Substrate (34080; Thermo Fisher Scientific) and imaged with BioRAD ChemiDocTM Imaging System. Images were viewed and quantified using ImageJ Fiji software v. 2.14.0.
Table 1.
Antibodies used in the study.
| Western blot analysis | |||
|---|---|---|---|
| Target | Catalog Number | Company | Concentration |
| Emerin | 30853S | Cell Signaling | 1:1000 |
| Nesprin-3 | ab12031 | Abcam | 1:1000 |
| Anti-HA | 2367S | Cell Signaling | 1:1000 |
| KLF5 | 21017–1-AP | Proteintech | 1:2000 |
| Vinculin | A302–535A | Bethyl | 1:10,000 |
| Goat Anti-Rabbit | AP132P | EMD | 1:5000 |
| Goat Anti-Mouse | AP124P | EMD | 1:3500 |
| Immunofluorescence | |||
| Emerin | NCL-Emerin 4G5 | Leica | 1:100 |
| Emerin | ab156871 | Abcam | 1:150 |
| Nesprin-3 | ab190010 | Abcam | 1:100 |
| Anti-HA | 2367S | Cell Signaling | 1:100 |
| Goat anti-Rabbit Alexa Fluor 594 | A-11037 | Thermo Fisher Scientific | 1:200 |
| Goat anti-Rabbit Alexa Fluor 488 | A-11008 | Molecular Probes | 1:200 |
| Goat anti-Rabbit Alexa Fluor 568 | A-11011 | Thermo Fisher Scientific | 1:200 |
| Goat anti-Mouse Alexa Fluor 488 | A-11001 | Thermo Fisher Scientific | 1:200 |
| Goat anti-Mouse Alexa Fluor 594 | A-11005 | Thermo Fisher Scientific | 1:200 |
| Proximity Ligation Assay | |||
| Emerin | ab156871 | Abcam | 1:100 |
| Nesprin-3 | ab190010 | Abcam | 1:100 |
| Nesprin-3 | ORB215388 | Biorbyt | 1:100 |
| Anti-HA | 2367S | Cell Signaling | 1:100 |
| Immunoprecipitation | |||
| Emerin | 30853S | Cell Signaling | 1:50 |
| Normal IgG Rabbit | 2729S | Cell Signaling | 1:50 |
| ChIP-seq | |||
| H3K4me1 | C15410037 | Diagenode | 1 μg |
| H3K4me3 | 07–449 | Millipore | 1 μg |
| H3K27Ac | ab4729 | Abcam | 1 μg |
qRT-PCR
RNA was extracted from cells using TRIzol reagent (15596018; Invitrogen). Two mg of RNA was reverse transcribed into cDNA using the High-Capacity cDNA kit (4368814; Applied Biosystems). Specific genes were amplified using the SsoAdvanced Universal SYBR Green Supermix (1725274; BioRAD), Bio-Rad CFX384 RT-PCR instruments and validated primers specific for the genes of interest. Expression primers used were: mouse Emerin (Forward- TCGTCATCTTCTTCATTCTCCTATC; Reverse-TCATTATAGTCCTTGCTCTGGTAAA), mouse Nesprin-3 (Forward- GGCCACAGAACTTCGTCTC; Reverse-GTCCCTCATTGCACTGTACTC), human Nesprin-3 (Forward- CCTGCAGAGGAAAAGCAAAC; Reverse-GTGGTCACAACGATCCACTG), mouse KLF5 (Forward- GGACTCATACGGGCGAGAAG; Reverse-TGCGGTTTAAAGGATGGCAG), human KLF5 (Forward- ACTGCGATTACCCTGGTTGC; Reverse-CAGCCTTCCCAGGTACACTT), human KRAS (Forward- CAAGAGTGCCTTGACGATACA; Reverse-GACCTGCTGTGTCGAGAATATC), mouse PRT (Forward- AAGTGTTTATTCCTCATGGACTGA; Reverse-CTCCCATCTCCTTCATGACATC), human PRT (Forward- CCTGGCGTCGTGATTAGTGAT; Reverse-AGACGTTCAGTCCTGTCCATAA), mouse TBP (Forward- GAAGTTCCCTATAAGGCTGGAAG; Reverse-AGGAGAACAATTCTGGGTTTGA), human TBP (Forward- TATAATCCCAAGCGGTTTGC; Reverse-CCCAACTTCTGTACAACTCTAGCA).
Immunocytochemistry and confocal microscopy
Cells were plated on 25 mm glass coverslips (174985; Thermo Fisher Scientific) at concentrations to reach 70% confluency the following day. Cells were fixed with 3.2% paraformaldehyde (043368.9 M; Thermo Fisher Scientific) in 1× PBS. For nuclear membrane permeabilization, cells were incubated with 0.1% Triton X-100 (T8787–50ML, Millipore) in 1× PBS on ice for 8 min. Cells were blocked with 5% BSA in 1× PBS for 1 h at room temperature. Primary antibodies (Table 1) were diluted in respective blocking buffers and incubated on cells for 1 h at room temperature. Cells were washed three times in 1× PBS for 5 min each. Cells were incubated in secondary antibodies diluted 1:1000 in blocking buffer for 1 h at room temperature. Cells were washed three times in 1× PBS for 5 min each. Coverslips were mounted using Prolong gold with DAPI (P36931; Thermo Fisher Scientific). Confocal images were taken on an Olympus LSM8000 microscope using Zeiss Zen Blue software. Z-stack images were acquired using an interval of 0.5 mm between slices. Images were analyzed and quantified using ImageJ Fiji software.
Image analysis
Line scan was performed as previously described.7 1012U cells treated with doxycycline and stained with Emerin and Nesprin-3 were compared to untreated cells and used for line scan analysis. Z-stack images through the entire nucleus were taken, and the middle plane was used for line scan. A measurement line was drawn through the long axis of each nucleus to plot Emerin and Nesprin-3 localization intensity across the nuclear periphery and interior. Profile plots for each signal were constructed by normalizing the signal to a distance of 100 to account for differences in nuclear size. The percent of signal at the periphery was taken by dividing the sum of the signal within 3 mm of the nuclear edge by the signal across the entire line. Nuclear size quantification was performed in ImageJ utilizing Z-stack images. Each image was smoothed three times, and nuclear volume and surface area were calculated using 3D Object Counter using a threshold of 30 and minimum size of 50. Clustered and partial nuclei were excluded.
Human tissue analysis
Three normal pancreas, Nesprin-3 low and Nesprin-3 high IHC images were downloaded from Human Protein Atlas (www.proteinatlas.org).42 Nuclear size was analyzed in QuPath using Cell detection with Hematoxylin OD. The following parameters were set: Background radius: 5 px, Median filter radius: 0 px, Sigma: 3.5 px, Minimum area: 20 px2, Maximum area: 150 px2, Threshold: 0.03, Max background intensity: 2, Cell expansion: 5 px. Nucleus area was visualized in GraphPad Prism.
Proliferation assays
Cells were plated on 10 cm plates to reach 60% confluency the following day. The day after plating, cells were transfected with siRNA or DNA mixtures and incubated for 24 h. 1000 cells per well were then replated onto 96-well plates for CyQUANT or MTT assays. For both assays, a baseline assay was performed 6 h after replating. CyQUANT assays were performed using the CyQUANT NF Cell Proliferation Assay Kit (C35007; Invitrogen) per manufacturer’s protocol, 72 h post replating. Briefly, 5× HBSS buffer was diluted to 1× with nanopure H2O and 1 °CyQUANT NF dye reagent was added. Existing cell culture medium was removed and 100 μL CyQUANT mixture was added to each well. The cells were incubated at 37 °C for 1 h. The plate was read using fluorescence (λexc = 485 nm, λem = 525 nm). MTT assays were performed using Thiazolyl Blue Tetrazolium Bromide (M5655–1G, Sigma-Aldrich). Briefly, 10 μL 5 mg/mL MTT was added directly to existing cell medium 72 h post replating. Cells were incubated at 37 °C for 2 h. Following incubation, the MTT-media mixture was removed and 100 μL DMSO was added and incubated at 37 °C for 15 min. The plate was read using absorbance at 540 nm. For all proliferation assays, the 72 h timepoint was normalized to the baseline reading. Protein and RNA for knockdown and overexpression validation were collected 24 h after replating.
RNA-seq
RNA-seq was performed in 1012U cells under −Dox and +Dox conditions. Total RNA was isolated using TRIzol reagent and RNeasy Mini Kit (74104; Qiagen) following treatment with doxycycline for 72 h. RNA-seq library generation and sequencing was performed by the London Regional Genomics Center at the University of Western Ontario (Canada). Processing and trimming of RNA-seq FASTq files were performed using Trim Galore 0.6.4 and cutadapt 2.10 with Python 3.6.5. FastQC was utilized to evaluate read quality, and STAR 2.7.3a was used for alignment with mm10 genome as reference using the Mayo Biocluster. FeatureCounts (Galaxy Version 1.6.4) and DESeq2 (Galaxy Version 2.11.40.6) were utilized to evaluate read counts and differential gene expression. Significantly differentially expressed genes were identified using the following metrics: base mean ≥1, log2 fold change >1 or <−1, and an FDR of ≤0.05. Heatmaps were generated in RStudio software using pheatmap version 1.0.12. RNA-seq files from KLF5 knock-out cell lines were obtained from SRA.27 DESeq2 was utilized to evaluate read counts and differential gene expression for Nesprin-3. Nesprin-3 differential gene expression was visualized in RStudio 2025.05.0 using ggplot2.49 SP and KLF family transcription factor expression was evaluated in the non-targeted control RNA-seq samples in human cell lines.50 Z-scores were calculated in each cell line separately across all genes and visualized in heatmaps using pheatmap version 1.0.12.
ChIP-seq analysis and visualization
ChIP-seq was performed in 1012 U cells under −Dox and +Dox conditions. Cells were –harvested and resuspended in 1× PBS at 1 × 106 per ml. Next, 1.5 × 106 cells were transferred to a fresh 1.5 mL Eppendorf tube for each treatment condition. Cells were pelleted at 500 × g for 6 min at 4 °C using the low- or no-brake option to prevent the pellet from being disturbed during centrifuge slow-down. Cells were lysed in 120 μL native ChIP-seq lysis buffer (0.1% Triton X-100, 0.1% Deoxycholate, 1× complete protease inhibitor) per 1.5 × 106 cells. Shearing was performed by incubating samples with 2.5 U/mL of MNase (M0247S; New England Biolabs) for 10 min at 25 °C. Histone marks were immunoprecipitated with 1 mg of antibody (Table 1) per 0.5 × 106 cells for each sample. Input represents 1% of total sample. Library generation and sequencing was outsourced to the University of Western Ontario. Bowtie2 was utilized to align H3K4me1, H3K4me3, and H3K27Ac ChIP-seq FASTq files to the mouse genome (mm10).51 Significant peaks were identified with MACS2 call peak and differential binding analysis was performed with DiffBind with peaks summits set to 500.52–54 Significance was determined with an FDR of ≤ 0.05. For enhancer analysis, overlapping H3K27Ac and H3K4me1 called peaks unique to −Dox and +Dox conditions were identified with bedtools Intersect Intervals.55 Genes unique to −Dox and +Dox enhancers, from overlapping peaks, were identified using GREAT analysis for the single nearest gene within 5–30 kb.29 Kegg analysis was performed on genes identified from GREAT analysis to identify differentially expressed pathways and the associated genes.30 Histone mark differential binding was validated in a selection of genes identified in Kegg analysis. Specifically, the average histone mark signal over the 5–30 kb range was Log2 normalized to the input signal of that mark56 and visualized in a heatmap in RStudio 2025.05.0 using pheatmap version 1.0.12. KLF5 ChIP-seq BigWig files were obtained from SRA27 and visualized utilizing Integrative Genome Viewer.
Statistical analysis
All statistical analyses and construction of graphs were done using GraphPad Prism 9 software. All data are comprised of at least three independent experiments. Data were checked for normality and lognormality to direct statistical test choice. For experiments with two groups, either a parametric unpaired t test, if normality was passed, or a non-parametric two-tailed Mann–Whitney U test, if normality not passed, were performed. For experiments with three or more groups that passed normality, a one-way ANOVA test and the corresponding Tukey’s multiple comparison for statistical significance were performed. If normality was not passed a Kruskal–Wallis test followed by a Dunn’s test for statistical significance were performed. P values < 0.05 were considered statistically significant.
Supplementary Material
Supplemental data for this article can be accessed online at https://doi.org/10.1080/10985549.2026.2670710.
Funding
This work was supported by the National Institutes of Health/National Cancer Institute CA265050 and Minnesota Colon Cancer Research Foundation Grants MNCRCF2025, MNCRCF2026. The Biochemistry Core at Sanford Research, which facilitated these studies, is supported by Institutional Development Awards from the National Institute of General Medical Sciences and the National Institutes of Health under grant GM103620. This work was supported by the Sanford Burnham Prebys Proteomics Core.
Footnotes
Disclosure statement
No potential conflict of interest was reported by the author(s).
Data Availability Statement
The data for all figures presented are available within the manuscript. Files from RNA-seq (BioProject ID PRJNA1112965) and ChIP-seq (BioProject ID PRJNA1113127) analyses have been deposited in the NCBI Sequence Read Archive data repository. Raw BioID data, including counts and LFQ intensities, has been provided in Supplementary Table S1. KLF5 RNA- and ChIP-seq data reanalyzed in this study are available in the NCBI Sequence Read Archive data repository (BioProject ID PRJNA1283515).27
References
- 1.Singh I, Lele TP. Nuclear morphological abnormalities in cancer: a search for unifying mechanisms. Results Probl Cell Differ. 2022;70:443–467. doi: 10.1007/978-3-031-06573-6_16. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 2.Tan SY, Tatsumura Y. George Papanicolaou (1883–1962): discoverer of the pap smear. smedj. 2015;56:586–587. doi: 10.11622/smedj.2015155. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 3.Flores LF, Tader BR, Tolosa EJ, Sigafoos AN, Marks DL, Fernandez-Zapico ME. Nuclear dynamics and chromatin structure: implications for pancreatic cancer. Cells. 2021;10:2624. doi: 10.3390/cells10102624. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 4.Chen P, Tomschik M, Nelson KM, Oakey J, Gatlin JC, Levy DL. Nucleoplasmin is a limiting component in the scaling of nuclear size with cytoplasmic volume. J Cell Biol. 2019;218:4063–4078. doi: 10.1083/jcb.201902124. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 5.Levy DL, Heald R. Nuclear size is regulated by importin α and Ntf2 in xenopus. Cell. 2010;143:288–298. doi: 10.1016/j.cell.2010.09.012. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 6.Mishra S, Levy DL. Nuclear F-actin and Lamin A antagonistically modulate nuclear shape. J Cell Sci. 2022;135. doi: 10.1242/jcs.259692. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 7.Flores LF, Marks DL, Vera RE, Sigafoos AN, Tolosa EJ, Almada LL, Pease DR, Toruner MD, Chang B, Tader BR, et al. Emerin is an effector of oncogenic KRAS-driven nuclear dynamics in pancreatic cancer. JCI Insight. 2025; 10. doi: 10.1172/jci.insight.187799. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 8.Kamisawa T, Wood LD, Itoi T, Takaori K. Pancreatic cancer. Lancet. 2016;388:73–85. doi: 10.1016/s0140-6736(16)00141-0. [DOI] [PubMed] [Google Scholar]
- 9.Holaska JM, Wilson KL. An emerin “proteome”: purification of distinct emerin-containing complexes from hela cells suggests molecular basis for diverse roles including gene regulation, mrna splicing, signaling, mechano-sensing, and nuclear architecture. Biochemistry. 2007;46:8897–8908. doi: 10.1021/bi602636m. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 10.Berk JM, Tifft KE, Wilson KL. The nuclear envelope LEM-domain protein emerin. Nucleus. 2013;4:298–314. doi: 10.4161/nucl.25751. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 11.Clements L, Manilal S, Love DR, Morris GE. Direct interaction between emerin and lamin A. Biochem Biophys Res Commun. 2000;267:709–714. doi: 10.1006/bbrc.1999.2023. [DOI] [PubMed] [Google Scholar]
- 12.Ho CY, Jaalouk DE, Vartiainen MK, Lammerding J. Lamin A/C and emerin regulate MKL1-SRF activity by modulating actin dynamics. Nature. 2013;497:507–511. doi: 10.1038/nature12105. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 13.Fernandez A, Bautista M, Wu L, Pinaud F. Emerin self-assembly and nucleoskeletal coupling regulate nuclear envelope mechanics against stress. J Cell Sci. 2022;135. doi: 10.1242/jcs.258969. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 14.Ranade D, Pradhan R, Jayakrishnan M, Hegde S, Sengupta K. Lamin A/C and emerin depletion impacts chromatin organization and dynamics in the interphase nucleus. BMC Mol Cell Biol. 2019;20:11. doi: 10.1186/s12860-019-0192-5. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 15.Marano N, Holaska JM. Emerin interacts with histone methyltransferases to regulate repressive chromatin at the nuclear periphery. Front Cell Dev Biol. 2022;10:1007120. doi: 10.3389/fcell.2022.1007120. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 16.Wood LD, Yurgelun MB, Goggins MG. Genetics of familial and sporadic pancreatic cancer. Gastroenterology. 2019;156:2041–2055. doi: 10.1053/j.gastro.2018.12.039. [DOI] [PubMed] [Google Scholar]
- 17.Wilhelmsen K, Litjens SH, Kuikman I, Tshimbalanga N, Janssen H, van den Bout I, Raymond K, Sonnenberg A. Nesprin-3, a novel outer nuclear membrane protein, associates with the cytoskeletal linker protein plectin. J Cell Biol. 2005;171:799–810. doi: 10.1083/jcb.200506083. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 18.Ketema M, Wilhelmsen K, Kuikman I, Janssen H, Hodzic D, Sonnenberg A. Requirements for the localization of nesprin-3 at the nuclear envelope and its interaction with plectin. J Cell Sci. 2007;120:3384–3394. doi: 10.1242/jcs.014191. [DOI] [PubMed] [Google Scholar]
- 19.Collins MA, Bednar F, Zhang Y, Brisset JC, Galbán S, Galbán CJ, Rakshit S, Flannagan KS, Adsay NV, Pasca di Magliano M. Oncogenic KRAS is required for both the initiation and maintenance of pancreatic cancer in mice. J Clin Invest. 2012;122:639–653. doi: 10.1172/jci59227. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 20.Roux KJ, Kim DI, Burke B. BioID: a screen for protein-protein interactions. Curr Protoc Protein Sci. 2013;74: 19.23.1–19.23.14. doi: 10.1002/0471140864.ps1923s74. [DOI] [PubMed] [Google Scholar]
- 21.Kim DI, Jensen SC, Noble KA, Kc B, Roux KH, Motamedchaboki K, Roux KJ. An improved smaller biotin ligase for bioid proximity labeling. Mol Biol Cell. 2016;27:1188–1196. doi: 10.1091/mbc.E15-12-0844. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 22.Wong X, Cutler JA, Hoskins VE, Gordon M, Madugundu AK, Pandey A, Reddy KL. Mapping the micro-proteome of the nuclear lamina and lamina-associated domains. Life Sci Alliance. 2021;4:e202000774. doi: 10.26508/lsa.202000774. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23.Salpingidou G, Smertenko A, Hausmanowa-Petrucewicz I, Hussey PJ, Hutchison CJ. A novel role for the nuclear membrane protein emerin in association of the centrosome to the outer nuclear membrane. J Cell Biol. 2007; 178:897–904. doi: 10.1083/jcb.200702026. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 24.Kimura H Histone modifications for human epigenome analysis. J Hum Genet. 2013;58:439–445. doi: 10.1038/jhg.2013.66. [DOI] [PubMed] [Google Scholar]
- 25.Castro-Mondragon JA, Riudavets-Puig R, Rauluseviciute I, Lemma RB, Turchi L, Blanc-Mathieu R, Lucas J, Boddie P, Khan A, Manosalva Pérez N, et al. JASPAR 2022: the 9th release of the open-access database of transcription factor binding profiles. Nucleic Acids Res. 2022;50:D165–d173. doi: 10.1093/nar/gkab1113. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 26.Fernandez-Zapico ME, Lomberk GA, Tsuji S, DeMars CJ, Bardsley MR, Lin YH, Almada LL, Han JJ, Mukhopadhyay D, Ordog T, et al. A functional family-wide screening of SP/KLF proteins identifies a subset of suppressors of KRAS-mediated cell growth. Biochem J. 2011;435:529–537. doi: 10.1042/bj20100773. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 27.Ekstrom TL, Dai Z, Thiel J, Kanakan A, Saha BJ, Manjunath M, Schacherer N, Bousounis P, Siegler EL, Abdelrahman AM, et al. KLF5 controls subtype-independent highly interactive enhancers in pancreatic cancer to regulate cell survival. Sci Adv. 2026;12:eaea2106. doi: 10.1126/sciadv.aea2106. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28.Cao L, Huang C, Cui Zhou D, Hu Y, Lih TM, Savage SR, Krug K, Clark DJ, Schnaubelt M, Chen L, et al. Proteogenomic characterization of pancreatic ductal adenocarcinoma. Cell. 2021;184:5031–5052.e26. doi: 10.1016/j.cell.2021.08.023. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 29.McLean CY, Bristor D, Hiller M, Clarke SL, Schaar BT, Lowe CB, Wenger AM, Bejerano G. Great improves functional interpretation of cis-regulatory regions. Nat Biotechnol. 2010;28:495–501. doi: 10.1038/nbt.1630. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 30.Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45:D353–D361. doi: 10.1093/nar/gkw1092. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 31.Chen Z, Chen M, Fu Y, Zhang J. The KRAS signaling pathway’s impact on the characteristics of pancreatic cancer cells. Pathol Res Pract. 2023;248:154603. doi: 10.1016/j.prp.2023.154603. [DOI] [PubMed] [Google Scholar]
- 32.Lin R, Bao X, Wang H, Zhu S, Liu Z, Chen Q, Ai K, Shi B. TRPM2 promotes pancreatic cancer by PKC/MAPK pathway. Cell Death Dis. 2021;12:585. doi: 10.1038/s41419-021-03856-9. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 33.Lanfredini S, Thapa A, O’Neill E. Ras in pancreatic cancer. Biochem Soc Trans. 2019;47:961–972. doi: 10.1042/bst20170521. [DOI] [PubMed] [Google Scholar]
- 34.Jurcak NR, Rucki AA, Muth S, Thompson E, Sharma R, Ding D, Zhu Q, Eshleman JR, Anders RA, Jaffee EM, et al. Axon guidance molecules promote perineural invasion and metastasis of orthotopic pancreatic tumors in mice. Gastroenterology. 2019;157:838–850.e6. doi: 10.1053/j.gastro.2019.05.065. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 35.Biankin AV, Waddell N, Kassahn KS, Gingras MC, Muthuswamy LB, Johns AL, Miller DK, Wilson PJ, Patch AM, Wu J, et al. Pancreatic cancer genomes reveal aberrations in axon guidance pathway genes. Nature. 2012;491: 399–405. doi: 10.1038/nature11547. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 36.Wheeler MA, Davies JD, Zhang Q, Emerson LJ, Hunt J, Shanahan CM, Ellis JA. Distinct functional domains in nesprin-1alpha and nesprin-2beta bind directly to emerin and both interactions are disrupted in X-linked Emery-Dreifuss muscular dystrophy. Exp Cell Res. 2007;313:2845–2857. doi: 10.1016/j.yexcr.2007.03.025. [DOI] [PubMed] [Google Scholar]
- 37.Zhang Q, Ragnauth CD, Skepper JN, Worth NF, Warren DT, Roberts RG, Weissberg PL, Ellis JA, Shanahan CM. Nesprin-2 is a multi-isomeric protein that binds lamin and emerin at the nuclear envelope and forms a subcellular network in skeletal muscle. J Cell Sci. 2005;118:673–687. doi: 10.1242/jcs.01642. [DOI] [PubMed] [Google Scholar]
- 38.Mislow JM, Holaska JM, Kim MS, Lee KK, Segura-Totten M, Wilson KL, McNally EM. Nesprin-1alpha self-associates and binds directly to emerin and lamin A in vitro. FEBS Lett. 2002;525:135–140. doi: 10.1016/s0014-5793(02)03105-8. [DOI] [PubMed] [Google Scholar]
- 39.Zhang Q, Skepper JN, Yang F, Davies JD, Hegyi L, Roberts RG, Weissberg PL, Ellis JA, Shanahan CM. Nesprins: a novel family of spectrin-repeat-containing proteins that localize to the nuclear membrane in multiple tissues. J Cell Sci. 2001;114:4485–4498. doi: 10.1242/jcs.114.24.4485. [DOI] [PubMed] [Google Scholar]
- 40.Rajgor D, Shanahan CM. Nesprins: from the nuclear envelope and beyond. Expert Rev Mol Med. 2013;15:e5. doi: 10.1017/erm.2013.6. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 41.Liao L, Qu R, Ouang J, Dai J. A glance at the nuclear envelope spectrin repeat protein 3. Biomed Res Int. 2019; 2019:1651805–1651808. doi: 10.1155/2019/1651805. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 42.Uhlen M, Zhang C, Lee S, Sjöstedt E, Fagerberg L, Bidkhori G, Benfeitas R, Arif M, Liu Z, Edfors F, et al. A pathology atlas of the human cancer transcriptome. Science. 2017;357. doi: 10.1126/science.aan2507. [DOI] [PubMed] [Google Scholar]
- 43.Petrie RJ, Koo H, Yamada KM. Generation of compartmentalized pressure by a nuclear piston governs cell motility in a 3D matrix. Science. 2014;345:1062–1065. doi: 10.1126/science.1256965. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 44.Lu W, Schneider M, Neumann S, Jaeger VM, Taranum S, Munck M, Cartwright S, Richardson C, Carthew J, Noh K, et al. Nesprin interchain associations control nuclear size. Cell Mol Life Sci. 2012;69:3493–3509. doi: 10.1007/s00018-012-1034-1. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 45.Birendra K, May DG, Benson BV, Kim DI, Shivega WG, Ali MH, Faustino RS, Campos AR, Roux KJ. VRK2A is an A-type lamin-dependent nuclear envelope kinase that phosphorylates BAF. Mol Biol Cell. 2017;28:2241–2250. doi: 10.1091/mbc.E17-03-0138. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 46.Zhao KW, Sikriwal D, Dong X, Guo P, Sun X, Dong JT. Oestrogen causes degradation of KLF5 by inducing the E3 ubiquitin ligase EFP in ER-positive breast cancer cells. Biochem J. 2011;437:323–333. doi: 10.1042/bj20101388. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 47.Vivekanandhan S, Yang L, Cao Y, Wang E, Dutta SK, Sharma AK, Mukhopadhyay D. Genetic status of KRAS modulates the role of neuropilin-1 in tumorigenesis. Sci Rep. 2017;7:12877. doi: 10.1038/s41598-017-12992-2. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 48.Sears RM, May DG, Roux KJ. Bioid as a tool for protein-proximity labeling in living cells. Methods Mol Biol. 2019;2012:299–313. doi: 10.1007/978-1-4939-9546-2_15. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 49.Wickham H Ggplot2: elegant graphics for data analysis. New York (NY): Springer-Verlag; 2016. [Google Scholar]
- 50.Ekstrom TL, Saha BJ, Thiel J, Dai Z, Manjunath M, Schacherer N, Souto Y, Abdelrahman AM, Essmann F, Truty MJ, et al. KLF5 controls subtype-independent highly interactive enhancers in pancreatic cancer to regulate cell survival. BioRxiv. 2025. doi: 10.1101/2025.06.30.662301. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 51.Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25. doi: 10.1186/gb-2009-10-3-r25. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 52.Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, Nusbaum C, Myers RM, Brown M, Li W, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9:R137. doi: 10.1186/gb-2008-9-9-r137. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 53.Feng J, Liu T, Qin B, Zhang Y, Liu XS. Identifying ChIP-seq enrichment using MACS. Nat Protoc. 2012;7:1728–1740. doi: 10.1038/nprot.2012.101. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 54.Ross-Innes CS, Stark R, Teschendorff AE, Holmes KA, Ali HR, Dunning MJ, Brown GD, Gojis O, Ellis IO, Green AR, et al. Differential oestrogen receptor binding is associated with clinical outcome in breast cancer. Nature. 2012;481:389–393. doi: 10.1038/nature10730. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 55.Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26:841–842. doi: 10.1093/bioinformatics/btq033. [DOI] [PMC free article] [PubMed] [Google Scholar]
- 56.Wallner S, Schröder C, Leitão E, Berulava T, Haak C, Beißer D, Rahmann S, Richter AS, Manke T, Bönisch U, et al. Epigenetic dynamics of monocyte-to-macrophage differentiation. Epigenetics Chromatin. 2016;9:33. doi: 10.1186/s13072-016-0079-z. [DOI] [PMC free article] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Data Availability Statement
The data for all figures presented are available within the manuscript. Files from RNA-seq (BioProject ID PRJNA1112965) and ChIP-seq (BioProject ID PRJNA1113127) analyses have been deposited in the NCBI Sequence Read Archive data repository. Raw BioID data, including counts and LFQ intensities, has been provided in Supplementary Table S1. KLF5 RNA- and ChIP-seq data reanalyzed in this study are available in the NCBI Sequence Read Archive data repository (BioProject ID PRJNA1283515).27
