ZNF865 (BLST) Regulates Back Pain via Cell Senescence and DNA Damage Mechanisms

ZNF865 (BLST) Regulates Back Pain via Cell Senescence and DNA Damage Mechanisms | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results ZNF865 (BLST) Regulates Back Pain via Cell Senescence and DNA Damage Mechanisms View ORCID Profile Christian Lewis , View ORCID Profile Hunter Levis , Jonah Holbrook , View ORCID Profile Jacob T. Polaski , Timothy D. Jacobsen , Josh Mizels , View ORCID Profile Sarah E. Gullbrand , View ORCID Profile Brian Diekman , View ORCID Profile James C. Iatridis , View ORCID Profile Jason Gertz , Brandon Lawrence , Robby D. Bowles doi: https://doi.org/10.1101/2025.06.13.659603 Christian Lewis 1 University of Utah, Department of Biomedical Engineering , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Christian Lewis Hunter Levis 1 University of Utah, Department of Biomedical Engineering , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Hunter Levis Jonah Holbrook 1 University of Utah, Department of Biomedical Engineering , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site Jacob T. Polaski 2 Huntsman Cancer Institute, University of Utah , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jacob T. Polaski Timothy D. Jacobsen 3 Icahn School of Medicine at Mount Sinai, Department of Orthopaedics , New York, NY Find this author on Google Scholar Find this author on PubMed Search for this author on this site Josh Mizels 4 University of Utah, Department of Orthopaedic Surgery , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site Sarah E. Gullbrand 5 University of Pennsylvania, Department of Orthopaedic Surgery , Philadelphia, PA 6 Corporal Michael J. Crescenz VA Medical Center , Philadelphia, PA Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Sarah E. Gullbrand Brian Diekman 7 Thurston Arthritis Research Center, University of North Carolina at Chapel Hill , Chapel Hill, NC 8 Lampe Joint Department of Biomedical Engineering, University of North Carolina at Chapel Hill and North Carolina State University , Raleigh, NC Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Brian Diekman James C. Iatridis 3 Icahn School of Medicine at Mount Sinai, Department of Orthopaedics , New York, NY Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for James C. Iatridis Jason Gertz 1 University of Utah, Department of Biomedical Engineering , Salt Lake City, UT 2 Huntsman Cancer Institute, University of Utah , Salt Lake City, UT 9 University of Utah, Department of Oncological Sciences , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Jason Gertz Brandon Lawrence 4 University of Utah, Department of Orthopaedic Surgery , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site Robby D. Bowles 1 University of Utah, Department of Biomedical Engineering , Salt Lake City, UT 4 University of Utah, Department of Orthopaedic Surgery , Salt Lake City, UT Find this author on Google Scholar Find this author on PubMed Search for this author on this site For correspondence: robert.bowles{at}utah.edu Abstract Full Text Info/History Metrics Preview PDF Abstract Senescence has been shown to contribute to the progression of aging related diseases, yet its regulation in the intervertebral disc (IVD) remains poorly understood. Recently, ZNF865 (BLST) was identified as a previously uncharacterized zinc finger protein, that regulates a wide array of genes related to protein processing, cell senescence and DNA damage repair. Here, we show that ZNF865 expression decreases with age and pathology in human and mouse IVD samples. Depletion of ZNF865 induces senescence, SASP expression and DNA damage in both human and rat healthy NP cells. Conversely, restoration in degenerative NP cells mitigates senescence, SASP expression and DNA damage, enhances ECM anabolism, and restores chromatin accessibility and gene expression to a healthy state. In vivo, CRISPRi of ZNF865 induces disc degeneration and painful behaviors. Collectively, our findings establish ZNF865 as a regulator of genome stability and a potential therapeutic target for mediating senescence/DNA damage in aging related diseases. One Sentence Summary ZNF865 is a novel modulator of genome stability and senescence and regulates pain outcomes in the intervertebral disc via reversible senescence and DNA damage pathways. Introduction Cellular senescence has been increasingly recognized as a fundamental driver of disease progression in aging related diseases. During senescence, cells undergo stable cell cycle arrest often described as a protective mechanism against oncogenic transformation and excessive DNA damage 1 , 2 . These cells secrete many pro-inflammatory catabolic factors collectively known as the senescence associated secretory phenotype (SASP) that can disrupt the tissue microenvironment and further sustain the inflammatory environment through complex positive feedback loops 2 , 3 . In many aging related disorders, the increased presence of senescent cells is thought to be a major contributor to disease progression. In the central nervous system, the accumulation of senescent glial cells has been shown to contribute to inflammation and neurodegeneration potentially accelerating diseases such as Alzheimer’s 4 and Parkinson’s 5 . In the lungs, senescent fibroblast populations are associated with impaired repair mechanisms and fibrotic remodeling in pulmonary fibrosis 6 , 7 . In the musculoskeletal system, senescent chondrocytes and other cell types of the joint have been shown to contribute to cartilage matrix breakdown and synovial inflammation in osteoarthritis 8 as well as increased inflammation and matrix breakdown in degenerative intervertebral discs (IVD) 9 . Copp et al. showed that chondrocytes obtained from femur cartilage of osteoarthritic patients had significantly higher levels of DNA damage, but also that the amount of DNA damage present in healthy chondrocytes increased with age 10 , indicating that excessive DNA damage may be an indicator of cells capable of contributing to disease pathology in osteoarthritis, and other aging-related musculoskeletal disorders. Back pain is one of the most common musculoskeletal disorders and a leading cause of disability worldwide 11 , with a large population of people expected to experience it at some point in their lifetime 12 , 13 . Chronic back pain is often debilitating, with 10% of individuals with back pain reportedly being forced to change jobs or stop working entirely due to severe pain 14 . Back pain incidence and severity have also been shown to increase with age with 40% of back pain attributed to degenerative disc disease (DDD) 15 . While back pain etiology is diverse and multifactorial, growing evidence suggests biological aging and senescence may be major contributing factors. Traditionally DDD onset is associated with the breakdown of the disc tissue due to genetics 16 , aging 17 , loading history 18 and inadequate metabolite transport 19 . However, as DDD progresses, studies have shown there is decreased cell activity due to increased cell senescence 20 , 21 . Additionally, degenerative discs are often characterized by increased presence of pro-inflammatory cytokines, proteases and pain related growth factors consistent with the SASP phenotype. These changes to extracellular matrix composition can lead to mechanical instability, pain and disability rates 19 , 22 – 24 . Targeting senescent cells via senolytics within the IVD has provided promising results in reducing the release of SASP factors in isolated IVD cells and intact human IVDs 25 , 26 . Oral administration of o-vanillin and RG-7112 resulted in decreased pain and histological degeneration in the IVD of sparc-null mice 27 . Additionally, intra-peritoneal delivery of the senolytics dasatinib and quercetin from 14 to 23 months of age significantly improved nucleus pulposus (NP) and annulus fibrosus outcomes in wild type mice 28 . Overall, this work further distinguishes cellular senescence as an important mediator of IVD degeneration and highlights the critical need to develop novel tools that target senescence directly to prevent degeneration and reduce pain caused by SASP. Recently, we identified ZNF865, also known as BLST, a gene with no known functional role, and no studies on its molecular actions, and utilized it as a cell engineering tool to enhance the tissue deposition and mechanical properties of engineered IVDs 29 . Initially identified in multiple different genome wide CRISPR perturbations screens 30 , 31 , in addition to its use as a cell engineering tool, RNA sequencing data has also linked ZNF865 to DNA damage and senescence 29 . ZNF865 is a 1059 amino acid long C2H2 zinc finger (ZNF) protein, containing six disordered regions, twenty unique ZNF domains, two transactivation domains (TADs) and two TGEKP linkers ( Figure 1A ). The presence of the two TADs indicate that ZNF865 likely stimulates transcriptional activity by interacting with transcriptional regulators 32 . Despite broad expression across all cell types ( Figure 1B ), there are no publications functionally characterizing the biological role of ZNF865. Additionally, a rare disease mutation in ZNF865 is associated with a wide range of symptoms including global developmental delay, joint hypermobility, muscle hypertonia, loose inelastic skin, hand tremors, and eye/ear/brain abnormalities indicating effects across a wide range of tissues and organs in human patients 33 . To this end, with the growing link between senescence and back pain, and ZNF865s link to DNA damage and senescence, we investigated ZNF865 in the context of back pain within human back pain patient samples and multiple rodent models. Download figure Open in new tab Figure 1: ZNF865 Regulates Similar Gene Networks Across Multiple Cell Types. A Structure of the human ZNF865 gene. B ZNF865 expression across different human tissue sources. nTPM = normalized transcripts per million. C dCas9-VPR CRISPRa system and D plasmid maps for lentiviral delivery targeting ZNF865 or NTC. E MA plots from RNA-seq on ZNF865 upregulated HEK293 cells, F ASCs and G human dhNPCs. H Volcano plots showing the top common molecular process affected by ZNF865 upregulation across HEK293 cells, I ASCs and J human dhNPCs using the KEGG 2021 Human Pathway. NTC is non-target control . Here we link ZNF865 as a primary driver of discogenic back pain through the regulation of senescence and DNA damage. This investigation further elucidated the role DNA damage plays in the onset and perpetuation of the degenerative phenotype in DDD indicating that the depletion of ZNF865 is sufficient to induce DNA damage, disc degeneration and painful phenotypes both in vitro and in vivo . Importantly, upregulation of ZNF865 is capable of reprogramming human disc cells from back pain patients to return them to a healthy phenotype. Furthermore, these insights identify ZNF865 as a key modulator in maintaining genomic stability and as a potential therapeutic target in other aging related diseases hallmarked by senescence. Results ZNF865 upregulation alters related cell functions across multiple cell types To begin investigating the molecular role of ZNF865, CRISPR activation (CRISPRa) was used to increase the expression of ZNF865 using the dCas9-VPR system alongside a non-target control ( Figure 1C , 1D ). RNAseq was then performed on three different human cell types; human embryonic kidney cells (HEK293), adipose derived stem cells (ASCs) and degenerative human NP cells (dhNPCs), and comparisons were made across the CRISPRa ZNF865 and non-target samples ( Figure 1E , 1F , 1G ). Upregulation of ZNF865 in each of these cell types showed 2,699, 6,647 and 1,752 significantly differentially expressed genes, respectively. Utilizing Enrichr 34 , analysis of these genes showed a number of significantly differentially regulated molecular processes using the KEGG 2021 Human Pathway that were shared across cell types ( Figure 1H - 1J , Table S1 - S3 ) 35 . The notable pathways across all three cell types included cell cycle, DNA replication, cellular senescence, mismatch repair, p53 signaling, mTOR signaling, base excision repair, and homologous recombination. The shared molecular processes suggested a role for ZNF865 in regulating cell division, DNA repair, and senescence. ZNF865 expression is correlated with human degenerative disc disease and aging To evaluate whether ZNF865 expression is pathologically relevant within DDD, gene expression was measured via qRT-PCR on dhNPCs from patients undergoing surgical treatment for low back pain, and healthy NP cells (hNPCs) from trauma patients at first passage after plating in vitro ( Figure 2A ). We observed a strong correlation (R 2 =0.94, p=0.0281) between ZNF865 expression and age in hNPCs and significantly decreased expression in dhNPC samples at matched ages ( Figure 2B ). Immunofluorescence staining for ZNF865 matched the qRT-PCR data with significantly less ZNF865 staining in dhNPC samples compared to hNPCs ( Figure 2C , 2D ). ZNF865 staining was localized to the nucleus of the cells, consistent with bioinformatic predictions that ZNF865 is a transcription factor. Of note, 28.7% of hNPCs stained for ZNF865, while only 16.8% of dhNPCs were stained, indicating that ZNF865 may not be expressed by all cells at a detectable level at a given time ( Figure 2E ). Download figure Open in new tab Figure 2: ZNF865 Expression is Correlated with Aging, Degeneration, and DNA Damage in Human and Mouse NP cells. A Normalized ZNF865 expression in hNPCs (n=4) and dhNPCs (n=23) in individual patients and B grouped by pathology. C Representative images from ZNF865 immunofluorescence staining (Blue=DAPI, Red=ZNF865). D ZNF865 immunofluorescence staining per cell normalized to the dhNPC level (n=7). E Percent of dhNPCs and hNPCs expressing ZNF865 from immunofluorescence staining (n=3-4). F Workflow for scRNA-seq in the caudal NP of aging mice. G ZFP865 expression level in the caudal NP of aging mice. H Percent of cells positive for ZFP865 in the aging mouse model. I Normalized β-galactosidase staining of healthy (n=4) and dhNPCs (n=5). J Percent DNA in the comet tail of naïve hNPCs (n=5) and dhNPCs (n=8) measured by comet assay. K Images showing single cell comets with percent damage, and representative dhNPCs and hNPC examples. ( Scale bars = 50 µm. significance *=p<0.05, **=p<0.01, ***=p<0.001 ) To determine whether expression of this gene also changes in mouse IVD, ZFP865 expression (mouse orthologue of ZNF865) 36 was examined in the caudal NP of mice at 4, 12 and 24 months of age ( Figure 2F ). These results were consistent with the human data, where ZFP865 expression decreased with increasing age. Both population average expression level ( Figure 2G ) and percentage of cells expressing ( Figure 2H ) showed decreases, with a strong negative correlation between cells expressing ZFP865 and age in the mouse (r = −0.9976, p = 0.04435). Cell senescence has been proposed as a potential mechanism of DDD with degenerative discs showing increased levels of senescence ( Figure 2I ) 20 , 21 . However, the functional and molecular links between senescence and DDD are still relatively unknown. Given that our RNA-seq data indicates that both senescence and closely associated DNA repair pathways are significantly altered following CRISPRa ZNF865 upregulation ( Figure 1J ), we investigated DNA damage in the human dhNPCs and hNPCs using alkaline comet assays. Here we observed that the percentage of DNA damage in the comet tail in dhNPCs was significantly higher than in hNPCs ( Figure 2J , 2K ). Together, this data indicates a relationship between ZNF865 expression and DDD, as well as DDD progression and DNA damage. Moving forward we investigated the relationship between ZNF865 expression and cellular senescence and DNA damage in human NP cells. ZNF865 downregulation in healthy human NP cells induces senescence and DNA damage hNPCs were isolated from discarded surgical tissue from deidentified trauma patients and transduced with the CRISPRi vector targeting ZNF865 (ZNF865-KRAB) or a non-target control ( Figure 3A , Figure S1 ). Following transduction, cells were monitored for cell proliferation over 6 days. Samples transduced with ZNF865-KRAB showed decreased proliferative capacity compared to cells transduced with the non-target vector ( Figure 3B ). To determine if the reduction in cell growth observed in the ZNF865-KRAB samples was consistent with entry into senescence, hNPCs were stained for key markers of cellular senescence, including senescence associated β-galactosidase (β-gal), p16 INK4A (p16) and p21 waf1/Cip1 (p21). The ZNF865-KRAB NP cells showed significantly higher levels of all three senescence markers, with 5.4-fold, 4.7-fold, and 6.9-fold increases in staining, respectively ( Figure 3C - 3E , 4H ). Entry into senescence has commonly been linked to increased DNA damage as a way to prevent damaged cells from perpetuating their genomes in division 2 , 39 , and ZNF865 regulation had been demonstrated to regulate DNA damage repair pathways ( Figure 1 ). To quantify the amount of DNA damage across patient samples, we performed alkaline comet assays 3 weeks after transduction. Here, the ZNF865-KRAB NP cells showed significantly higher levels of DNA damage compared to the non-target samples, indicating a functional connection between decreased ZNF865 expression and DNA damage ( Figure 3F , 3G ). This data is consistent with the relationships observed between ZNF865, senescence and DNA damage in human patient NP cells ( Figure 2 ) and our RNAseq data ( Figure 1 ). Download figure Open in new tab Figure 3: CRISPRi of ZNF865 decreases proliferation rates, increases senescence and increases DNA damage in hNPCs. A Workflow for CRISPRi of ZNF865 or non-target in hNPCs. B Quantified cell proliferation of healthy human non-target and ZNF865-KRAB NP cells (n=4). C SA-β-Galactosidase, D p16 and E p21 staining normalized to non-target (n=4). F Percent DNA in comet tail in individual patients obtained by comet assay. G Donor mean percent DNA in comet tail (n=4). H Representative images from SA-β-Galactosidase, p16 and p21 staining, and nuclei from comet assays. (Scale bar β-Gal = 100 µm, scale bar others = 50 µm. Significance *=p<0.05, **=p<0.01, ***=p<0.001 ) Download figure Open in new tab Figure 4: CRISPRa of ZNF865 restores proliferation, decreases senescence, decreases SASP phenotype and decreases DNA damage to healthy levels in dhNPCs. A Workflow for CRISPRa of ZNF865 or non-target in dhNPCs. B Quantified cell proliferation of human non-target dhNPCs and ZNF865-VPR dhNPCs (n=5). C Doubling times of non-target dhNPCs, ZNF865-VPR dhNPCs and hNPCs (n=5). D Representative images from β-Gal, p16 and p21 staining, and nuclei from comet assays. E SA-β-Galactosidase, F p16 and G p21 staining normalized to non-target (n=4-5). H Percent DNA in comet tail in non-target dhNPCs, ZNF865-VPR dhNPCs and hNPCs. I Donor mean percent DNA in comet tail (n=7). J Analyte concentration for SASP related inflammatory cytokines, proteases and proteins present in ZNF865-VPR dhNPCs normalized to non-target samples (n=4). Non-target levels indicated by dashed line. K SASP score normalized to hNPC values (n=4-5). L Representative images of pellet culture. M Total collagen production per DNA and N GAG production per DNA and O normalized anabolic scoring from pellet culture in pooled non-target dhNPCs and two ZNF865-VPR dhNPC patient samples (n=3-6). P Plot of showing relationship of SASP score and anabolic score after ZNF865 upregulation. Scale bars = 50 µm. Significance *=p<0.05, **=p<0.01, ***=p<0.001 ) ZNF865 upregulation in degenerative human NP cells decreases senescence, SASP phenotype and DNA damage To begin studying the effects of ZNF865 regulation in pathologic cells, dhNPCs were isolated from discarded surgical tissue from deidentified patients undergoing surgical treatment for chronic back pain. Cells were then transduced with ZNF865 upregulation cassettes (ZNF865-VPR) or a non-target control ( Figure 4A , Figure S1 ). Following transduction, cells were monitored for cell proliferation via fluorescent imaging over 5 days. We observed significant increases in cell growth within 3 days in the ZNF865-VPR cells compared to the non-target controls ( Figure 4B ). This led to a significant decrease in doubling time in the ZNF865-VPR dhNPCs compared to the non-target dhNPC controls, at 2.49 days and 4.09 days respectively, with non-significant differences compared in the naïve hNPCs with a doubling time of 3.2 days ( Figure 4C ). We stained each treatment group for key senescence markers 10 days after transduction. Here, the non-target dhNPCs showed the highest levels of staining in β-gal, p16 and p21, with significantly lower levels of senescence staining in the ZNF865-VPR dhNPCs and hNPCs. Additionally, the ZNF865-VPR dhNPCs and hNPCs were non-significantly different from each other ( Figure 4D - 4G ). To assess the ability to directly regulate DNA damage, we performed alkaline comet assays. Here we found that upregulation of ZNF865 significantly decreased the percentage of DNA in the comet tails in all patients, with five of the seven patients reaching significance ( Figure 4I ). When the mean damage of each donor is considered, there is a significant decrease in the ZNF865-VPR dhNPCs compared to the non-target controls, with the ZNF865-VPR dhNPCs matching the healthy controls ( Figure 4H , 4I ). This demonstrates ZNF865’s ability to directly regulate DNA damage and senescence observed in human samples from painful degenerative discs. To better understand the downstream effects of reversing senescence and DNA damage in these cell populations, a Luminex assay was run on 10 SASP related proteins using the cell culture supernatant from the ZNF865-VPR dhNPCs, non-target control dhNPCs, and naïve hNPCs. Nine of the ten analytes including IL-1β, IL-6, TNF-α, IL-8, MMP-1, MMP3, MMP-12, NGF-β and VEGF-a in the ZNF865-VPR showed significant decreases in concentration compared to the non-target control, with the non-significant difference in MCP-1 trending towards a reduced level ( Figure 4J ). From the Luminex assay, a composite normalized SASP score was created that weighs each analyte equally and is normalized to the hNPC supernatants. Here, we saw significant differences in SASP score between the ZNF865-VPR dhNPCs and non-target control dhNPCs, as well as significant differences between the non-target control dhNPCs and hNPCs. The ZNF865-VPR dhNPCs had a SASP score that was slightly higher than the hNPCs, but this difference was not statistically significant ( Figure 4K ). This indicates that increased expression of ZNF865 leads to a decreased SASP profile and secretome more similar to hNPCs. To assess ZNF865 regulation on ECM deposition within the disc, 3D pellet culture was performed using non-target dhNPCs pooled from three patient samples, and two ZNF865-VPR dhNPC patient samples ( Figure 4L ). After 3 weeks of culture, we saw significant increases in collagen production in both ZNF865-VPR patients compared to the pooled non-target pellets, with no significant increases in glycosaminoglycan (GAG) production ( Figure 4M , 4N ). From these measurements, a composite anabolic score was created, weighing both collagen production and GAG production equally and normalizing to the non-target samples. Here we saw significant differences in anabolic score for both ZNF865-VPR samples compared to the non-target ( Figure 4O ). Additionally, when plotting the SASP scores and anabolic scores from both patients on the same plot we see a clear separation between the non-target samples and the ZNF865-VPR samples ( Figure 4P ), indicating that ZNF865 upregulation results in decreased inflammation and increased ECM deposition capability. ZNF865 upregulation in degenerative NP cells restores healthy chromatin state Given the close relationship between senescence and changes in chromatin landscape 40 , we performed Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq), with primary patient samples ( Table S4 ), to identify differences in the chromatin accessibility associated with ZNF865 expression. When comparing hNPCs with dhNPCs we found 1163 sites that exhibit differential chromatin accessibility that were associated with 1705 genes identified by Genomic Regions Enrichment of Annotations Tool (GREAT) 41 . Further analysis of these genes revealed functional enrichment of 45 pathways, including glycosaminoglycan synthesis, MAPK signaling pathway, and cellular senescence ( Figure 5A , Table S5 ). Upon upregulation of ZNF865 in dhNPCs, a total of 575 differential chromatin accessibility sites were identified, with 4 significantly enriched pathways compared to the dhNPCs ( Figure 5B , Table S6 ). However, when comparing ZNF865-VPR dhNPCs to the hNPCs, we observed only 70 sites with differential chromatin accessibility (94% reduction) and no enrichment of any pathways ( Figure 5C ). Principal component analysis using this cohort ( Figure 5D ) showed that dhNPCs become more related to hNPCs upon CRISPRa of ZNF865, suggesting that upregulation restores some molecular features of hNPCs. Similarly, additional RNA-seq was performed using comparable groups ( Table S7 ), which showed a corresponding shift in gene expression when ZNF865 is upregulated in dhNPCs towards a healthier gene expression state ( Figure 5E ). Download figure Open in new tab Figure 5: CRISPRa of ZNF865 returns human NP cells to healthy chromatin state and gene expression. A MA plots of the ATAC-seq data illustrating sites of differential chromatin accessibility comparing dhNPCs versus hNPCs, B dhNPCs versus ZNF865-VPR dhNPCs and C ZNF865-VPR dhNPCs versus hNPCs with relevant enriched pathways. D Principal component analysis of ATAC-seq data shows the relationship between hNPCs, dhNPCs, and ZNF865-VPR dhNPCs. E Principal component analysis of RNA-seq data shows the relationship between hNPCs, dhNPCs, and ZNF865-VPR dhNPCs. F Hypergeometric test results are shown as a bar plot to quantify the degree of gene-level overlap of RNA-seq and ATAC-seq datasets between comparisons of hNPCs and dhNPCs (green and yellow bars) and ZNF865-VPR dhNPCs and dhNPCs (blue and purple bars). Bars corresponding to the number of shared upregulated genes from RNA-seq and genes associated with increased chromatin accessibility from ATAC-seq for each comparison are denoted as “up.” Conversely, bars corresponding to the number of shared downregulated genes from RNA-seq and genes associated with decreased chromatin accessibility from ATAC-seq for each comparison are denoted as “down.” G Bar plot illustrates the differences in COL6A2 expression between hNPCs dhNPCs, and ZNF865 VPR dhNPCs using normalized counts from RNA-seq. H Browser tracks display differences in ATAC-seq peaks (shaded region) in COL6A2 from hNPCs, dhNPCs, and ZNF865-VPR dhNPCs. I Bar plot illustrates the differences in USP12 expression between hNPCs, dhNPCs, and ZNF865-VPR dhNPCs using normalized counts from RNA-seq. J Browser tracks display differences in ATAC-seq peaks (shaded region) in USP12 from hNPCs, dhNPCs, and ZNF865-VPR dhNPCs. To assess concordance between differentially expressed genes identified by RNA-seq and genes associated with changes in chromatin accessibility by ATAC-seq, we employed hypergeometric tests to quantify the degree of overlap between either upregulated or downregulated genes and genes associated with either increased or decreased chromatin accessibility. This analysis revealed statistically significant enrichment in the number of shared genes compared to what would be expected by chance across all comparisons ( Figure 5F ) . Specifically, we found 184 upregulated/increased accessibility and 84 downregulated/decreased accessibility genes that are shared between comparisons of hNPCs and dhNPCs, representing >3-fold and >2-fold enrichment compared to the amount of overlap expected by chance, respectively. Similarly, we found 51 upregulated/increased accessibility and 33 downregulated/decreased accessibility genes that are shared between comparisons of ZNF865-VPR dhNPCs and dhNPCs, representing >6-fold and >5-fold enrichment over the amount of overlap expected by chance, respectively. Within the individual genes, there are several notable examples, including COL6A2 ( Figure 5G ) , which was recently implicated as a potential marker of DDD 42 . Normalized counts from RNA-seq show upregulation of COL6A2 in both hNPCs and ZNF865-VPR dhNPCs, compared to dhNPCs. These differences in COL6A2 expression correlate with changes in chromatin accessibility, where we observed more prominent ATAC-seq peaks in both hNPCs and ZNF865-VPR dhNPCs compared to dhNPCs ( Figure 5H ) . We also found that USP12 , which has a known functional role as an inflammatory mediator involved in MMP production 43 , is downregulated in both hNPCs and ZNF865-VPR dhNPCs compared to dhNPCs ( Figure 5I ) . These differences in expression correlate with changes in chromatin accessibility, where both hNPCs and ZNF865-VPR dhNPCs show less prominent ATAC-seq peaks compared to dhNPCs ( Figure 5J ) . Collectively, these results provide evidence that ZNF865 upregulation shifts dhNPCs toward a healthier molecular phenotype. CRISPRi of ZNF865 induces DDD and painful outcomes in vivo To confirm that the associations observed between ZNF865 in the human patient samples, both in its relationship to age and back pain, and the molecular changes observed with ZNF865 regulation in human patient samples can directly drive pain outcomes, we knocked ZNF865 down directly in the rat interverbal disc with CRISPRi ( Figure 6A ). Prior to in vivo delivery, a rat orthologue of the human ZNF865 had to be verified. To begin, gRNA upstream of the ZNF865 orthologue from the Rat Genome Database 44 were designed ( Table 2 ), cloned into the dCas9-KRAB vector, and transduced into rat NP cells isolated from the caudal spine of Sprague Dawley rats before performing qPCR. Fold change in ZNF865 expression showed a significant reduction in gene expression in ZNF865-KRAB samples ( Figure 6B ). A significant reduction in cell proliferation was also observed ( Figure 6C ) alongside increases in the number of senescent cells via SA-β-gal and p16 staining ( Figure 6D , 6E , 6G ). Similarly to human samples, this was accompanied by an increase in DNA damage present within the cells ( Figure 6F ). Examining the expression of SASP related proteins showed increases in the secretion of MCP-1, and VEGF-A in the ZNF865-KRAB samples compared to non-target with ICAM-1 not reaching significance ( p =0.0625). ( Figure 6H - 6J ). Overall, the rat disc cells responded similarly to the human IVD cells with CRISPRi ZNF865 regulation. Download figure Open in new tab Figure 6: CRISPRi of ZNF865 induces DDD pathology and increases painful outcomes in vivo. A Experimental design of the 8-week animal model using the dCas9-KRAB repression system with plasmid maps for CRISPRi of ZNF865 or non-target. B Log fold change in ZNF865 expression C Quantified cell proliferation of rNPCs after transduction with ZNF865-KRAB or NT control vectors D Percent of SA-β-gal and E p16 stained cells 3 weeks after transduction (n=5). F Mean percent DNA in comet tail from alkaline comet assays (n=4). G Representative images from SA-β-gal, p16 and comet assays. H Normalized MCP-1, I VEGF-A, J and ICAM-1 concentrations. K Representative images of the fluoroscopy assisted lentiviral injections to the l4-L5 and L5-L6 IVDs. L Normalized paw withdrawal threshold (PWT) and M normalized paw withdrawal latency (PWL) (n=5-8). N Normalized disc height from the L4-L5 and O L5-L6 IVDs P Degenerative scoring and Q representative H&E and alcian blue staining for naïve, NT-control and ZNF865 KRAB animals. ( significance *=p<0.05 ) View this table: View inline View popup Download powerpoint Table 1: sgRNA primers targeting the CRISPRa/i systems to ZNF865 and a nontarget control in the human for upregulation and downregulation. View this table: View inline View popup Table 2: sgRNA primers targeting the CRISPRi system to ZNF865 and a nontarget control in the rat genome Following validation of the ZNF865 orthologue, ZNF865-KRAB or NT-control, were injected into the L4-L5 or L5-L6 discs of Sprague Dawley rats ( Figure 6K ), following protocols previously developed in our lab for targeted CRISPRi regulation of genes in the rat IVD 45 . Following injections, animals were evaluated for symptoms of mechanical and thermal nociception via von Frey and Hargreaves assessments respectively. After 5 weeks a significant increase in mechanical pain sensitivity was seen in the ZNF865-KRAB group compared to both naïve and NT-control groups ( Figure 6L ). Thermal pain sensitivity was also increased in the ZNF865-KRAB group compared to the naïve group at 8 weeks ( Figure 6M ). Following sacrifice, disc height was evaluated at each injected level. Here, a significant decrease in disc height was seen at the L4-L5 level but not the L5-L6 level between the ZNF865-KRAB group and the NT-control group ( Figure 6N , 6O ). Following staining with H&E and alcian blue, samples were scored for degeneration. ZNF865-KRAB groups showed significantly increased degenerative scores compared to both naïve and NT-control groups ( Figure 6P , 6Q ). Discussion Despite the broad expression of ZNF865 across all cell types 46 , this study represents the first functional characterization of the ZNF865 gene. Through CRISPRa-mediated upregulation in three distinct cell types; HEK293 cells, ASC s, and human NP cells, we identified a set of consistently enriched processes involving cellular senescence, DNA repair mechanisms, and cell cycle regulation pathways. These shared processes suggest that ZNF865 may act as a core regulator of molecular processes that underlie cell senescence and genomic stability. We demonstrated that ZNF865 expression decreases with age and pathology, and that this reduction is associated with increased levels of DNA damage and cell senescence. Here we show that patient derived dhNPCs had significantly lower ZNF865 gene expression, lower protein levels, and greater DNA damage compared to hNPC controls, highlighting a potential connection between ZNF865 depletion and DDD in humans. This data is supported by single cell RNA sequencing of the caudal NP of aging mice in vivo, which showed a parallel decline in ZFP865 expression and the proportion of cells expressing ZFP865 with increased age. Additionally in vivo, we show that depletion of ZNF865 in a rat model is sufficient to induce histological degenerative outcomes and behavioral pain. Together, this data suggests that the loss of ZNF865 may represent a feature of disc tissue aging, contributing to aging related phenotypes and degeneration such as cell senescence which contribute to the lack of regenerative capacity of the disc. CRISPR-based gene modulation systems such as CRISPRa and CRISPRi offer a method for modulating the expression of ZNF865, allowing us to define its functional roles. Using these CRISPR based gene modulation systems revealed that expression of ZNF865 is highly tunable and can lead to robust functional differences in both human and rodent models 29 , with gene expression shown to be associated with proliferation and cell cycle progression. Downregulation of ZNF865 in hNPCs led to reduced proliferation, increased senescence, and elevated DNA damage, all of which are features commonly associated with tissue degeneration and DDD. In rNPCs, these features were conserved, also showing reduced proliferation and increases in senescence, DNA damage and in SASP protein production including MCP-1 and VEGF-A. Conversely, when ZNF865 was upregulated in patient sample back pain associated dhNPCs, these phenotypes were reversed, with increased proliferation rates, decreased senescence marker staining, and decreased DNA damage. This also led to a significant reduction in many SASP factors such as proinflammatory cytokines (IL-1, IL-6, TNF-α, IL-8), proteases (MMP-1, MMP-3, MMP-12), and pain related growth factors (NGF, VEGF) commonly associated with discogenic back pain. After upregulation of ZNF865, these dhNPCs showed cell phenotypes more similar to the hNPCs. Additionally, in 3D pellet culture ZNF865 upregulation showed increased ECM deposition capability. These phenotypic differences between patient samples are consistent with gene regulatory changes. Both RNA-seq and ATAC-seq showed that CRISPRa of ZNF865 recapitulated global features of gene expression and chromatin accessibility found in hNPCs. This effect can be seen at the level of individual genes, where we observed many up- and down-regulated genes that also harbored sites of more open and closed chromatin. One notable example we identified is COL6A2, whose expression is higher in hNPCs and ZNF865-VPR dhNPCs compared to dhNPCs. This gene encodes a type VI collagen chain, revealing a potential genetic link with increased ECM deposition by ZNF865-VPR dhNPCs in 3D pellet culture. Overall, these data highlight concordant changes in phenotypes and gene regulation, where our data suggests ZNF865 drives an overall shift of both gene expression and chromatin state toward that of healthy NP cells. In vivo, we demonstrate that intradiscal CRISPRi knockdown of ZNF865 induces behavioral pain responses. Our behavioral data showed a significant increase in mechanical pain sensitization, and to a lesser degree, thermal pain sensitization in the ZNF865 knockdown group. These behavioral outcomes parallel the nociceptive changes observed in well-established models of disc degeneration 45 , 47 and support the role of ZNF865 in regulating pain-relevant pathways within the IVD. This was further corroborated by histological analysis in which ZNF865-KRAB injected discs exhibited reduced disc height and increased degenerative scores. This induction of pain behavior in the absence of a mechanical injury challenges traditional models of discogenic pain. Most rodent disc degeneration models rely on large needle puncture, enzymatic digestion, or biomechanical overloads, all of which introduce acute structural damage to the disc, which may confound interpretation of pain mechanisms. Here nociceptive sensitization was achieved through transcriptional downregulation, which highlights molecular dysfunction as a primary driver of pain rather than structural disruptions. Additionally, the gradual progression of pain behavior seen in the von Frey assessment mirrors the slow progression of degeneration seen in human disc disease, suggesting that ZNF865 knockdown initiates a biological cascade that cumulatively leads to degeneration and painful symptoms, more similar to human pathology. Collectively, this data suggests that ZNF865 plays a fundamental role in DDD and aging, specifically in maintaining genetic stability through DNA repair mechanisms that prevent the onset of senescence ( Figure 7 ). The rescue of dhNPCs from a senescent phenotype upon ZNF865 upregulation not only points to the importance of ZNF865 within the degenerative phenotype but also identifies ZNF865 as a potential therapeutic target for DDD. Previously, we demonstrated upregulation of ZNF865 in multiplex with ACAN and Col2a1 enhanced functional tissue deposition and improved biomechanical properties of human-scale engineered discs 29 . This suggests that in addition to decreasing inflammatory pathways caused by DNA damage and SASP, CRISPR gene therapies using ZNF865 have potential to be used in multiplex for a wide variety of purposes including increased protein production or mechanical properties of engineered tissues. Download figure Open in new tab Figure 7: Proposed model of ZNF865 function in IVD aging and degeneration: The harsh IVD microenvironment promotes DNA damage through oxidative, nutritional and replicative stress. High ZNF865 expression supports DNA damage repair mechanisms maintaining disc homeostasis. With age and in patients with back pain, ZNF865 levels are decreased, leading to impaired DNA repair and increased cell senescence, promoting inflammation, tissue catabolism, immune cell recruitment and nociceptor ingrowth and sensitization, ultimately driving mechanical instability and pain. This work also offers broader implications for other aging related degenerative diseases, where controlling cellular senescence and SASP may be paramount to understanding and slowing disease progression. Ultimately, preventing senescence may improve the health span of individuals by preventing the onset of many musculoskeletal or neurodegenerative disorders. ZNF865 was recently identified independently in a family of zinc finger proteins where increased gene expression was associated with better prognosis in esophageal cancer 48 , indicating that it may be a useful clinical target for understanding cancer progression. ZNF865 was also identified as a potentially relevant target in Parkinson’s disease 49 . Liu et al. suggested that a microRNA associated with ZNF865 protein expression can alleviate damage caused by Parkinson’s disease by preventing apoptosis and enhancing motor neuron performance 49 . Within our RNA sequencing data, several neurodegenerative diseases such as ALS, Parkinson’s and Alzheimer’s were listed as top pathways effected by upregulation of ZNF865, suggesting that targeting ZNF865 or its downstream pathways could represent a strategy for restoring tissue function in aging related degenerative diseases more broadly. This work indicates ZNF865 is an important gene for healthy cell function and genome maintenance. However, key questions remain that need to be clarified in future experiments. Identification of the upstream regulators and downstream effectors will allow for better regulatory networks to be established around ZNF865 and help to define the molecular interactions leading to the observed results. This will also allow for the development of blocking antibodies, inhibitors or activators that could allow for enhanced control of ZNF865 regulation outside of using CRISPR-based therapies. This work uses DDD as a model system for understanding DNA damage and senescence regulation by ZNF865, however, additional primary cell types and disease models should be considered to establish a broader understanding of the functional role that ZNF865 plays within aging related diseases. This is especially important considering the rare genetic disorders that have been associated recently with ZNF865 mutations 33 . Overall, this work also highlights the use of CRISPR gene modulation strategies for probing unknown gene function and identifies ZNF865 as an important gene in regulating cell cycle and DNA repair mechanisms, with significant implications for DDD related pain and other senescence associated degenerative disorders. This study serves as a first step in understanding the basic biology of the ZNF865 gene and establishes the importance of understanding ZNF865 dysregulation for both disease pathology, and as a potential therapeutic target for senescence associated disorders. Materials and Methods Experimental overview CRISPRa was used to upregulate the expression of ZNF865 in HEK 293 cells, ASCs and human dhNPCs and investigated for changes in gene expression via RNAseq. Next, ZNF865 expression and DNA damage were examined in healthy and dhNPCs to investigate DDD pathology. An aging mouse model from Jacobsen et al. was used to investigate an identified orthologue for ZNF865 expression changes with aging 50 . CRISPRi was then used to downregulate the expression of ZNF865 in hNPCs to investigate differences in proliferation, senescence and DNA damage as a result of decreased expression. Next, CRISPRa was used to upregulate the expression of ZNF865 in dhNPCs to investigate differences in proliferation, senescence, SASP expression, DNA damage and chromatin accessibility as a result of increased expression. Finally, a rat orthologue of ZNF865 was identified and validated in vitro through CRISPRi of ZNF865 in rat NP cells, and an in vivo model run examining behavioral pain responses and degenerative outcomes/ General cell culture All cell culture was performed in standard culture conditions (21% O 2 , 5% CO 2 , 37°C), with media changes every 2-3 days. HEK 293 : Complete growth medium for cell culture for HEK 293 (ATCC CRL-1573) cells consists of HG-DMEM (ThermoFisher Scientific, Waltham, MA) supplemented with 10% fetal bovine serum (FBS) (ThermoFisher Scientific), 25μg/mL gentamicin (Corning), and 25mM HEPES (ThermoFisher). hASC: Complete growth medium for culturing of human adipose derived stem cells (ASCs, ATCC SCRC-4000) consists of Lonza ADSC Basal Medium (Lonza, Lexington, MA PT-3273), 10% MSC FBS (ThermoFisher Scientific), 5mL GlutaMax (ThermoFisher Scientific), 30μg/mL Gentamicin (Corning), and 15ng/mL Amphotericin (ThermoFisher Scientific). hNPCs: Primary human nucleus pulposus cells (hNPCs) were cultured in complete growth medium containing DMEM-HG (Thermo Fisher Scientific) supplemented with 10% FBS (Thermo Fisher Scientific), 50μg/mL gentamicin (Corning), 25mM HEPES (Thermo Fisher Scientific), and 2ng/mL of recombinant human fibroblast growth factor-basic (rhFGF, Peprotech, Cranbury, NJ, 100-18C). rNPCs: Primary rat nucleus pulposus cells (rNPCs) were cultured in Hams-F-12 Culture Media (Thermo Fisher Scientific), supplemented with 20% FBS (Thermo Fisher Scientific, 25ug/mL gentamicin (Corning), 8 mM HEPES, and 1.25 ug/mL amphotericin b (Thermo Fisher Scientific). Cells were cultured to 90% confluency before passaging 1:3 with 0.25% trypsin/EDTA solution (Thermo Fisher Scientific) as previously described 51 . gRNA design Upregulation gRNA design was performed using Genome Target Scan 2 and ChopChop 52 , 53 . The 5’-UTR, the promoter region up to 500bp upstream of the human ZNF865 transcription start site (TSS), and up to 500bp upstream of the second exon were analyzed for gRNAs. The top 6 gRNAs from each analysis tool, with the least number of off-target sites, were compared and inspected using the BLAT tool of the UCSC genome browser to ensure minimal overlap between gRNAs and optimal placement near DNase hypersensitive regions 54 . A total of 6 gRNAS were selected for validation in vitro , as well as a gRNA that does not target the human genome (non-target control/NTC). Oligos for ZNF865 gRNAs were synthesized by the University of Utah DNA/Peptide Synthesis core and ThermoFisher Scientific. Downregulation gRNA design was performed using Genome Target Scan 2 (GT-Scan2) and ChopChop 52 , 53 . Briefly, a panel of 5 gRNAs were designed within the 5’-UTR and the promoter region up to 1000 bp upstream of the human ZNF865 TSS. Oligos for ZNF865 gRNAs were synthesized by the University of Utah DNA/Peptide Synthesis core and ThermoFisher Scientific. Downregulation gRNA design in the rat genome was performed using Genome Target Scan 2 (GT-Scan2) and ChopChop 52 , 53 . The rodent ZNF865 sequence was identified using the Rat Genome Database (RGSC Genome assembly v6.0) 44 . Here, the 5’UTR, the promoter region up to 1000 base pairs upstream of the identified rat ZNF865 transcription start site were analyzed for gRNAs. The top 15 gRNAs with the least number of off-target sites, from each analysis tool were inspected using the BLAT tool in the UCSC genome browser to ensure minimal overlap between gRNAs. These 15 gRNA were then selected for validation as well as a gRNA that does not target the rat genome (non-target/NT). Oligos for ZNF865 were synthesized by ThermoFisher Scientific. Cloning gRNAs, and a non-target that does not target the human or rat genome, were synthesized, annealed, phosphorylated, and ligated into individual U6-gRNA-UbC-DsRed-P2A-Bsr (upregulation gRNA) lentiviral upregulation expression vectors (Addgene, #83919) or KRAB CRISPRi downregulation expression vector pLV-hUbC-dCas9-KRAB-T2A-GFP (Addgene, #71237). Successful gRNA insertion was verified through Sanger sequencing 51 , 55 , 56 . Harvesting of primary hNPCs and culture Human healthy and degenerative nucleus pulposus tissues were obtained from surgical tissue waste from trauma patients with no history of back pain, and patients undergoing surgical treatment for back pain respectively. Briefly, NP tissue was rinsed twice with washing medium (DMEM-HG, 165 μg/mL gentamicin sulfate, 100μg/mL kanamycin sulfate, 1.25μg/mL amphotericin), minced, and enzymatically digested in washing medium with 0.3% (w/v) collagenase type II (Worthington Biochemical), 0.2% pronase (Sigma), and 5% FBS (ThermoFisher Scientific) for 2-3 hours at 37°C with 5% CO 2 under gentle agitation 51 . Isolated cells were passed through a 70μm cell strainer and washed twice. Cells were counted and plated at a density of 10,000 cells/cm 2 in NP cell culture medium (DMEM-HG with 10% FBS, 50μg/mL gentamicin, 25mM HEPES) supplemented with fresh 2ng/mL fibroblast growth factor-2 (FGF-2) (Peprotech). Cells were cultured in this medium at 37°C and 5% CO 2 in a humidified atmosphere and subcultured to 90% confluency, as previously described 51 . Harvesting of primary hNPCs and culture Rat IVD tissue was obtained from the caudal spine (tails) of sacked Sprague Dawley rats. Upon acquisition, tails were macroscopically dissected to remove skin and connective tissue surrounding the spine. Once IVD were located AF and NP tissue were separated and placed in washing medium (DMEM-HG, 165 ug/mL gentamicin sulfate, 100ug/mL kanamycin sulfate, 1.25ug/mL amphotericin b). NP tissue was then digested in 0.01% pronase (Millipore Sigma), 0.0125% collagenase type 2 (Worthington Biochemical) and 5% FBS (ThermoFisher Scientific) overnight at 37°C, Following incubation, cells were centrifuged and washed with warm PBS (ThermoFisher Scientific) twice before being passed through a 70 um cell strainer. Cells were then re-centrifuged and re-suspended in culture media and plated in flasks pre-coated with LN5 supernatant and cultured until use at passage 1-2. Lentivirus production The gRNA plasmid DNA was used to produce lentivirus, as previously described 51 , 55 , 56 . The amplified gRNA plasmid DNA was co-transfected into HEK 293T cells with psPAX2 (Addgene, plasmid #12260) and pMD2.G (Addgene, plasmid #12259) lentiviral packaging plasmids to create a lentivirus, as previously described 45 . Briefly, HEK 293T cells were seeded at a density of 62,500 cells/cm 2 . The following day, lentiviral plasmids were added to the cells with Lipofectamine 2000 (ThermoFisher Scientific), following the manufacturer’s protocol. After 24 hours, the cell supernatant was discarded and replaced with fresh medium. Cell supernatant containing the lentiviral vectors was collected at 48 and 72 hours, concentrated to 100X, and stored at −80°C until use 56 . dCas9-VPR upregulation transduction HEK 293/ASC: HEK 293 cells and ASCs were plated at a density of 20,000 cells/cm 2 and 5,000 cells/cm 2 , respectively, allowed to attach overnight, and transduced with EF-1α-dCas9-VPR-PuroR (dCas9-VPR) (Addgene, #99373) lentivirus, (1:20 dilution) in growth medium containing 4μg/mL polybrene, the following day. Both cell types were subjected to antibiotic selection (1μg/mL) over the course of 3 days to select for dCas9-VPR expressing cells. hNPCs: hNPCs (p1-p2) were plated at a density of 20,000 cells/cm 2 , allowed to attach overnight, and transduced with EF-1a-dCas9-VPR-eGFP (Addgene, #241307) lentivirus, (1:2 dilution) in growth medium containing 4μg/mL polybrene, the following day. 72-hours following transduction hNPCs were subjected to fluorescence activated cells sorting (FACS) to select for GFP+ expressing cells. After dCas9-VPR transduction, dCas9-VPR expressing HEK 293 cells, hNPCs, and ASCs were transduced with gRNA targeting lentiviral vectors. HEK293-VPR: HEK 293-VPR cells were plated at a density of 20,000 cells/cm 2 and transduced with gRNA vector (1:250 dilution) targeting ZNF865 (Addgene, #241309) or non-target the following day in complete media. hNPC-VPR : hNPCs were plated at a density of 5,000 cells/cm 2 in a 24-well plate and transduced with gRNA vector (1:160 dilution) targeting ZNF865 (Addgene, #241309) or non-target in complete media the following day. ASC-VPR : ASCs were plated at a density of 5,000 cells/cm 2 in a 24-well plate and transduced with gRNA vectors (1:160 dilution) targeting ZNF865 (Addgene, #241309), or non-target in complete media the following day. Transduced HEK 293, ASCs and hNPCs were examined for fluorescence after 72 hours and showed near 100% transduction efficiency. dCas9-KRAB downregulation transduction For ZNF865 downregulation, naïve HEK-293, naïve ASCs, naïve hNPCs, and naïve rNPCs were plated at a density of 20,000 cells/cm 2 , 5,000 cells/cm 2 , 10,000 cells/cm 2 or 10,000 cells/cm 2 respectively, in 24-well plates and allowed to attach overnight. The following day, 100X sgRNA targeting ZNF865 in the human or rat (Addgene, #241308) or non-target virus was diluted 1:16 in complete growth medium supplemented with 4μg/mL polybrene and used to transduce cells. Cells were examined for fluorescence after 72-hours and showed near 100% transduction efficiency. qRT-PCR Successfully transduced cells were analyzed for changes in ZNF865 gene expression by qRT-PCR (n=4). 72-hours post-transduction, RNA was isolated and harvested using the Quick-RNA Micro Kit (Zymo Research, R1051), complementary DNA (cDNA) was synthesized from the purified RNA with high-capacity cDNA reverse transcription kit with RNAse inhibitor (Applied Biosystems). cDNA was then used for qRT-PCR with TaqMan gene expression assays (Thermo Fisher Scientific) for human ZNF865 (Hs05052648_s1) or rat ZNF865 (Custom Assay ID APT2HZ6). Beta-2-microglobulin (B2M, Hs00187842_m1) was used as an internal standard and changes in ZNF865 expression was normalized to B2M expression 57 . Fold-change in mRNA expression relative to the non-target or ZNF865 edited cells was calculated using the ΔΔCT method. RNA-sequencing RNA-sequencing (RNAseq) was utilized to evaluate differential gene expression due to ZNF865 upregulation. HEK 293 cells, ASCs, and primary human NP cell samples ( Table S7 ) were analyzed for global changes in gene expression after modification with ZNF865 or NTC CRISPRa vectors. Briefly, Total RNA was isolated from samples using a Quick-RNA Kit (Zymo Research, Irvine, CA). Isolated Total RNA (100-1000ng) was prepared using an Illumina TruSeq Stranded mRNA Library Prep Kit with PolyA Selection and samples were submitted to the High-Throughput Genomics Shared Resource Core at the Huntsman Cancer Institute for sequencing on a NovaSeq 6000, utilizing a NovaSeq Reagent Kit v1.5_150×150bp sequencing with 33 million reads per sample. Using previously described methods, data were normalized and compared to non-target control cells 51 , 55 . Sequencing reads were aligned to the hg38 build of the human genome using HISAT2 58 and SAMtools 59 was used to convert SAM files to BAM files. Genes were defined by the University of California Santa Cruz (UCSC) Known Genes, and reads mapping to known genes were assigned using featureCounts from the SubRead package 60 , 61 . Reads were normalized and differential analysis was performed using EdgeR 62 ( Figure 1E - 1G ) or the DESeq2 package in R ( Figure 5E ) 60 , 63 . We used an adjusted p-value threshold of <0.05 to identify genes that were significantly up- or down-regulated when comparing patient samples in Figure 6 . Enriched GO biological and molecular processes were determined from lists of significantly differentially expressed genes using Enrichr 34 , 64 . Single cell RNA-sequencing Single cell RNA sequencing (scRNA-seq) previously performed in Jacobsen et al. was mined for the relationship of ZFP865, the identified orthologue of human ZNF865 36 , with aging in mice. In this study, mice were used at three ages: 4 months (skeletal maturity), 12 months (equivalent to human age with peak back pain disability), and 24 months (equivalent to human age with peak back pain prevalence) 50 . Cells were isolated from 6 whole naïve coccygeal IVDs per mouse, from 6 mice per age group and pooled. Samples were processed by the Genomics Core Facility at the Icahn School of Medicine at Mount Sinai. Samples were processed using the 10X Genomics Chromium 3’ Kit and sequenced using an Illumina S1 NovaSeq chip. Cell Ranger software mapped reads to mouse mm10-2020-A reference genome. Data was processed via Seurat including quality control filtering, normalization, and unsupervised clustering. Cell clusters were visualized with uniform manifold approximation and projection (UMAP). Differential gene expression analysis identified gene markers for individual clusters and canonical markers facilitated annotation of clusters. Data was subset to NP and Notochordal clusters for analysis of ZFP865 expression. ZFP865 expression level was determined at the cellular level for each age, percentage of cells positively expressing ZFP865 at each age was also determined. Cell proliferation quantification Following successful transduction, HEK 293, ASCs, and hNPCs transduced with either the CRISPRa upregulation system or CRISPRi downregulation system (n=4) were evaluated for cell proliferation over the course of 3-10 days. Individual cells were counted using ImageJ or a hemacytometer 65 . Briefly, fluorescing cell counts were obtained by uploading images into a stack and thresholding the images to ensure only individual cells are shown within the image. After thresholding, cell counts were obtained by analyzing particles, generating a mask, and then counting the masks generated while excluding cells on the edges of the image. In instances where cell density was too great for thresholding and visually isolating individual cells, cells were manually counted in ImageJ using the cell counter plugin 65 . Immunofluorescence (P16 INK4A , p21 CIP1 and ZNF865) Primary hNPCs or rNPCs were plated into cell culture chamber slides at a concentration of 10,000 cells/cm 2 and then transduced with the ZNF865 or NTC downregulation expression cassettes. Cells were then cultured for 10 or 24 days after which the cells were fixed with 10% formalin. Following fixation cells were permeabilized with 0.2% Triton X in PBS for 5 min. Cells were then washed three times with PBS and treated with blocking solution (5% goat serum [MP Biomedicals] in PBS for 1 hour at room temperature). Following treatment with blocking solution, the blocking solution was removed, and the cells were treated with primary antibody; p16 INK4A (p16) antibody (1:100 in 1% BSA in PBS; Proteintech), p21 CIP1 (p21) antibody (1:100 in 1% BSA in PBS; Proteintech) or ZNF865 (1:100 in 1% BSA in PBS; Atlas Antibodies), normal rabbit IgG (Invitrogen) or BSA only negative controls and incubated at room temperature for 2 hours. Following incubation cells were washed three times with PBS, treated with goat anti-rabbit Coralite 594 (1:100 Proteintech) in 1% BSA in PBS, and incubated for 1 hour at room temperature. Subsequently, cells were rinsed three times with PBS, treated with DAPI solution (3ng/mL) for 10 minutes, and mounted with Vectamount Mounting Medium prior to fluorescence imaging (Olympus IX73). Alkaline comet assay The alkaline comet assay was used to assess DNA damage in primary hNPCs and rNPCs. NP cells were trypsinized and resuspended in cold PBS at 1 × 10 6 cells per mL. 5 × 10 4 NP cells were mixed in a 3:7 volume ratio with 1% Low Melting Point Agarose (BP1360, Fisher Scientific). The mixture was then applied to a Superfrost Plus slide (12-550-15, Fisher Scientific) precoated with 1% normal melting point agarose (20-102, Apex). The slides were then transferred to the lysis solution (2.5 M NaCl, 100 mM disodium EDTA, 10 mM Tris Base, 1% Triton X, pH 10) and allowed to lyse overnight at 4°C in the dark. Immediately following lysis, the slides were washed 1x with dH 2 O and immersed in an electrophoresis tank containing alkaline electrophoresis solution (300 mM NaOH, 1 mM disodium EDTA, pH >13) for 30 min followed by electrophoresis at 1 V/cm for 20 min at 300 mA. Slides were then washed for 10 minutes in PBS and an additional 10 minutes in dH 2 O at 4°C and allowed to dry at RT overnight. Dried slides were then stained with 1ug/mL DAPI (62248, Thermo Fisher Scientific). At least 100 randomly selected stained nuclei per slide were imaged under an Olympus IX73 microscope. Comets were analyzed in ImageJ using the OpenComet plugin 66 . SA-β-Galactosidase staining Primary hNPCs or rNPCs were plated in a 24-well plate at a concentration of 10,000 cells/cm 2 and allowed to attach overnight. The following day attached cells were transduced with the ZNF865 (Guide 2 and Guide 3) or NTC downregulation expression cassettes. Cells were cultured for 10 or 24 days, at which point they were fixed and stained for SA-β-galactosidase using a Senescent Cell Histochemical Staining Kit (Sigma-Aldrich, CS0030-1KT). Cellular senescence was quantified by counting SA-β-galactosidase stained and unstained cells using ImageJ. Custom Human Luminex assay Following upregulation of ZNF865 in primary dhNPCs a custom Luminex panel (Procartaplex 11 analyte multiplexed assay, Thermo Fisher) was run for IL1-β, IL-6, IL-8(CXCL8), TNF-α, MCP-1(CCL2), MMP-1, MMP-12, MMP-3, MMP-8, NGF-β, VEGF-a. 10-14 days following transduction of dhNPCs with ZNF865-VPR or NT, degenerative and healthy samples were passaged and plated in T-25 flasks at the appropriate cell density and 3 mL of media added per flask. 3 days following, supernatant was collected and spun at 1400 rpm for 10 minutes to remove particulates, then aliquoted for storage at −80°C until use. Images were taken at the same time point as collection for normalization based on cell number. The Procartaplex plate was then run according to the manufacturer’s instructions using 50 µl of supernatant per sample per well and run on a Luminex Magpix. Each sample was run in duplicate. After running the plate, values were then averaged across each duplicate and divided by the cell density to account for differences in cell numbers. Values below the detection range were given concentrations of ½ the detection limit of each analyte. ZNF865-VPR dhNPC values were then normalized to the non-target across individual patients before averaging across all patients to get a normalized concentration for each analyte. A composite SASP score weighting each analyte equally was then calculated by normalizing each analyte to the highest concentration measured across all samples, then taking the sum across all analytes to get a composite score for each sample. Custom Rat Luminex assay Following downregulation of ZNF865 in the rNPCs a custom Luminex panel (Procartaplex 9 analyte multiplexed assay, Thermo Fisher Scientific) was run for IL1-β, IFN-γ, IL-6, IL-10, TNF-α, MCP-1(CCL2), NGF-β, VEGF-a, and ICAM-1. 21 days following transduction of rNPCs with ZNF865-KRAB or NT, cells were passaged and plated in 6 well plates at the appropriate cell density and 1.5 mL of media added per flask. 3 days following, supernatant was collected and spun at 1400 rpm for 10 minutes to remove particulates, then aliquoted for storage at −80°C until use. Images were taken at the same time point as collection for normalization based on cell number. The Procartaplex plate was then run according to the manufacturer’s instructions using 50 µl of supernatant per sample per well and run on a Luminex Magpix. Each sample was run in duplicate. After running the plate, values were then averaged across each duplicate and divided by the cell density to account for differences in cell numbers. Pellet Culture of dhNPCs To evaluate extracellular matric ECM deposition in the ZNF865 edited cells, 3D pellet cultures were performed. ZNF865-VPR dhNPCs and non-target dhNPCs were passaged and resuspended at 1.25 million cells/mL in serum free growth medium. Serum free media consisted of DMEM-HG with pyruvate (Thermo), 0.1 µm dexamethasone, 1% antibiotic/antimycotic, 0.17mM ascorbic acid 2-phosphate, 0.35mM proline and 1% ITS+. Due to low cell numbers, three non-target dhNPC samples were pooled to form a single non-target dhNPC group. 200 µL aliquots of cell suspension was pipetted into individual wells of 96-well v-bottom plate and spun at 270G for 10 min at 4°C. Cells were allowed to contract for 48 hrs to form pellets. Pellets were then gently lifted from the bottom of the plate. Media was changed every 3 days for 21 days and supernatant collected during every media change. Pellet biochemical analysis The total amount of collagen content within papain digested pellets and supernatant was analyzed using a modified hydroxyproline assay as previously described 29 , 56 . The total amount of DNA and GAG content within papain digested pellets and supernatant was analyzed using a previously described Hoescht dye assay and dimethylmethylene blue (DMMB) assay, respectively 29 , 56 , 67 , 68 . A composite anabolic score weighting the contributions of collagen and GAG equally by normalizing to the highest concentration measured across all samples, then taking the sum of each normalized collagen and GAG score to get a complete anabolic score. These were then normalized to the average non-target dhNPC value. ATAC-sequencing Primary patient samples ( Table S4 ) were prepared in technical duplicate, where 50,000 cells per replicate were taken from cryogenically frozen vials and centrifuged at 500 x g for 5 minutes at 4°C. The supernatant was aspirated, and the resulting cell pellet was washed with 250 µL of ice-cold 1X PBS, pH 7.4 (Thermo Scientific; 10010072). The cells were then centrifuged at 500 x g for 5 minutes at 4°C, the supernatant was aspirated, and the resulting cell pellet was resuspended in 50 µL of lysis buffer containing the following: 10 mM Tris-HCl, pH 7.4 (Sigma-Aldrich; T2194), 3 mM MgCl 2 (Sigma-Aldrich; M1028), 10 mM NaCl (Sigma-Aldrich; 71386), 0.1% NP-40 (Millipore Sigma; 492016), 0.1% Tween-20 (Bio-Rad; 1610781), and 0.01% digitonin (Invitrogen; BN20061). Samples were then lysed via incubation on ice for 15 minutes and centrifuged at 500 x g for 5 minutes at 4°C. Following lysis, 1 mL of wash buffer (10 mM Tris-HCl, pH 7.4, 3 mM MgCl 2 , 10 mM NaCl, and 0.1% Tween-20) was added to each sample, which was then centrifuged at 500 x g for 5 minutes at 4°C. The supernatant was aspirated, and the resulting cell pellet was resuspended a second time in 1 mL of wash buffer. To pellet nuclei, cells were centrifuged at 500 x g for 10 minutes at 4°C. The supernatant was aspirated, and each resulting nuclei pellet was resuspended in transposition buffer containing the following: 10 mM Tris-HCl, pH 8.0, 5 mM MgCl 2 , 10% Dimethylformamide (Fisher Bioreagents; BP1160), and 1 µL of TDE1 Tagment DNA enzyme per sample (Illumina; 20034197) diluted in PBS, pH 7.4. For the transposition reaction, nuclei were incubated at 37°C for 30 minutes in an Eppendorf ThermoMixer C (EP5382000015) set to rotate at 1000 RPM. Following transposition, DNA was purified using the Zymo DNA Clean & Concentrator-5 kit (D4014) according to the manufacturer’s instructions with a binding buffer to sample ratio of 5:1. The purified DNA was then amplified by PCR using a final concentration of 1X Phusion Master Mix (New England Biolabs; M0530), and 200 nM each of forward and reverse barcode primers. Amplified libraries were purified using a 1:1 volume ratio of AMPure XP beads (Beckman Coulter; A63882) and fragment size distribution and quality were assessed using a High-Sensitivity D1000 ScreenTape assay. Libraries were sequenced on a NovaSeq X at a depth of ∼50 million reads/sample and 100 bases from the 3’ end of each read were trimmed using seqtk 69 with the following parameters: trimfq -e 100. Trimmed reads were then aligned to the hg38 build of the human genome using Bowtie 70 with the following parameters: --chunkmbs 512 -m 1 -t --best -q -S -l 32 -e 80 -n 2. SAM files were converted into BAM files and sorted using SAMtools 59 . Peaks representing regions of accessible chromatin were called using MACS2 71 with an adjusted P-value cutoff of <0.05 and an mfold parameter between 15 and 100. FeatureCounts 61 was used to quantify reads that aligned to ATAC-seq peaks called across all samples. Reads were normalized and chromatin regions displaying differential accessibility were identified using the DESeq2 package for R 63 . To identify regions of chromatin with differential accessibility, we compared dhNPC samples to both hNPCs and ZNF865-VPR dhNPCs, and ZNF865-VPR dhNPCs to hNPCs. Animals All procedures were performed with the approval of the University of Utah Institutional Animal Care and Use Committee (IACUC). Here, 24 Sprague Dawley rats 10-12 weeks of age were used. Animals were allowed free movement in their home cages and were dual housed throughout the course of the study. Rats were allowed food and water ad libitum throughout. Routine health and wellness exams were performed and wight were performed weekly to assess general animal health. Rats were maintained at 72°F-74°F and 12/12h light/dark cycles (6am-6pm lights on). Non-Invasive, Image Guided. Cannulated Lentiviral Injections Intradiscal injections to the lumbar spine of rats were performed under sterile conditions and general anesthesia (2% isoflurane inhalation). Once under anesthesia, IVD levels were identified using anatomical landmarks, then a 30G needle was fed through a 20G cannula to inject NT or ZNF865 targeting CRISPRi lentiviral vectors (5ul, 10 6 -10 9 transducing units) to the L4-L5 and L5-L6 IVDs. Animals in the naïve group underwent anesthesia, however, did not receive any injections. Following injections, rats were closely monitored for complications and given access to food and water ad libitum . Von Frey Behavioral Assessment Mechanical pain sensitivity was measured using the Chaplan up down protocol for von Frey hairs 72 – 74 . Here, rats were allowed to acclimate to a wire bottom cage for 30-45 minutes or until cage exploration and self-grooming behaviors ceased. After acclimation, von Frey monofilaments designed to bend at a known specific force (0.4, 0.6, 1, 2, 4, 6, 10, and 15 g), were pressed against the plantar region of the hind paw. Withdrawal of the paw from the monofilament resulted in the next lowest filament in the series to be tested. If no withdrawal was observed, the next highest filament in the series was tested. This paradigm was repeated until five stimulations after the first change in response, up to nine total measurements. Each hind paw was tested separately beginning with the 2g monofilament with a minimum 2 minutes between each stimulation of the same paw. Testing was performed twice prior to lentivirus injections to establish a baseline, three days prior to injections, and then bi-weekly for the duration of the 8 week experiment. Data collection was performed blinded. Paw withdrawal threshold (PWT) was calculated as described in Chaplan et al 72 for each paw, animal and group at each time point. PWT was then normalized to each animals baseline measurements, then normalized to the average of the naïve, no injection group at each time point to examine the mechanical pain sensitivity over time. Hargreaves Behavioral Assessment Thermal sensitivity was measured using the Hargreaves assay to assess pain related symptoms of peripheral neuropathy in the rats 75 – 78 . In this assay, rats were placed in a heated glass bottom cage and allowed to acclimate for 30-45 minutes or until cage exploration and self-grooming behaviors cease. After acclimation, a focused light source was applied from below the cage to the plantar region of the hind paw, and the time it took for animals to withdraw their paws from the heat source was recorded. This process is repeated three times per animal per hindfoot with a 5-minute rest between trails. Testing was performed twice prior to lentivirus injections to establish a baseline, three days prior to injections, and then bi-weekly for the duration of the 8 week experiment. Data collection was performed blinded. This data was then used to calculate the paw withdrawal latency (PWL) for each paw, animal, and group at each time point. PWL was then normalized each animals baseline measurement, and then normalized to the average of the naïve, no injection group, at each time point to examine the thermal pain sensitivity over time. Tissue Harvesting and Sectioning After 8 weeks, animals were sacrificed via carbon dioxide inhalation. Immediately following sacrifice, rat lumbar spines were harvested from L3-S1 and fixed in formalin (Thermo Fisher Scientific) for 48hrs. Following fixation, samples were decalcified (Immunocal, StatLab) for 72 hours, changing decalcification solution every 24 hours. Following decalcification spine segments were embedded (Tissue-Tek TEC 5 Embedding Station, Sakura Finetek) in paraffin using melted Surgipath Formula ‘R’ paraffin pellets (Leica Biosystems). Samples were then sectioned (5um) (RM2255 Rotary Microtomes, Leica Biosystems) and placed on slides. Histological Assessment Following sectioning onto slides, samples were stained for hematoxylin and eosin (H&E) and alcian blue to visualize IVD morphology. Slides were then imaged under bright-field microscope (Olympus UC50) to evaluate gross morphological differences within each group. Disc Height Analysis Following histological staining with alcian blue, disc height was calculated manually in ImageJ by measuring the distance between CEPs in the center of the disc for each animal. Disc height for each animal was then normalized to the naïve group at the L4-L5 or L5-L6 levels. Degeneration Scoring Following histological staining with H&E and alcian blue IVD degeneration was assessed using a previously reported semi-quantitative degeneration grading scale 79 . The combined degeneration scale, scored from 0-10, is comprised of five categories. Annulus fibrosis, border between the AF and NP, NP cellularity, NP matrix, and endplate. Each category is scored on a scale from 0-2 with a 0 indicating normal IVD morphology, and a 2 representing severe degeneration 79 . Statistical Analysis RNA-seq statistical analysis was performed and described in their respective sections (Enrichr, DESeq2, α=0.05). Statistical analysis of qRT-PCR, Luminex, β-gal, p16, p21, ZNF865, Collagen, GAG, comets disc height and degenerative scoring were performed using JMP Pro or GraphPad Prism using a one-way ANOVA with Tukey’s post hoc analysis. Statistical analysis for proliferation and behavioral testing was performed using JMP Pro with a two-way ANOVA with Tukey’s post hoc analysis. Expression level scRNA-seq data was determined to be non-normal, significant differences between ages were tested for using Kruskal-Wallis test with Dunn’s comparisons and p<0.05 considered significant in GraphPad Prism. Simple linear regression and Pearsons’s correlation were performed on percent positivity and age with p<0.05 considered significant. Supplementary Materials Download figure Open in new tab Figure S1: ZNF865 is titratable in Human NP cells with CRISPRi and CRISPRa. A dCas9-KRAB repression system with plasmid maps used for CRISPRi of ZNF865 or non-target. B dCas9-VPR activation system and plasmid maps used for CRISPRa of ZNF865 or non-target. C Fold change ZNF865 expression when using CRISPRi or CRISPRa systems in ASCs (n=4-5). D Representative immunofluorescence images after transduction with CRISPRi/CRISPRa systems with ZNF865 or non-target. E Normalized ZNF865 immunofluorescence staining per cell when using the CRISPRi (n=3-4) and F CRISPRa systems with the top performing gRNAs (n=4-5) ( Scale bars = 50 µm. significance *=p<0.05) View this table: View inline View popup Download powerpoint Table S1: KEGG Human Genome 2021 significant pathways (From RNAseq in HEK293 cells) ZNF865-VPR versus NT. View this table: View inline View popup Download powerpoint Table S2: KEGG Human Genome 2021 significant pathways (From RNAseq in ASC s). ZNF865-VPR versus NT. View this table: View inline View popup Download powerpoint Table S3: KEGG Human Genome 2021 significant pathways (From RNAseq in dhNPCs) ZNF865 VPR dhNPCs vs Non Target dhNPCs View this table: View inline View popup Download powerpoint Table S4: ATACseq patient information View this table: View inline View popup Download powerpoint Table S5: KEGG 2021 Human Genome significant pathways (From ATACseq) dhNPCs versus hNPCs. View this table: View inline View popup Download powerpoint Table S6: KEGG 2021 Human Genome significant pathways (From ATACseq) dhNPCs versus ZNF865-VPR dhNPCs View this table: View inline View popup Download powerpoint Table S7: RNAseq patient information Funding National Institute of Health grant R01 AR083990 (RDB) National Institute of Health grant R01 AR074998 (RDB) National Institute of Health grant R01 AR080096 (JCI) Author Contributions Conceptualization: CL, HL, JH, JTP, TDJ, SEG, BD, JG, RDB Methodology: CL, HL, JH, JTP, BD, JG, RDB Investigation: CL, HL, JH, JTP, TDJ, JM Visualization: CL, HL, JH, JTP, TDJ, RDB Funding acquisition: RDB Project administration: CL, RDB Supervision: JCI, JG, RDB Writing-original draft: CL, HL, JH, JTP, TDJ, RDB Writing-review and editing: CL, HL, JH, JTP, TDJ, SEG, JCI, JG, BL, RDB Competing interests Authors declare that they have no competing interests Data and materials availability GEO accession # GSE311657 Acknowledgements Research reported in this publication was supported by the National Institute of Health under award numbers R01 AR083990, R01 AR074998 and R01 AR080096. Sequencing was performed at the DNA Sequencing Core Facility, University of Utah. FACs was performed at the Flow Cytometry Core Facility, University of Utah. Research reported in this publication utilized the High-Throughput Genomics Shared Resource at the University of Utah and was supported by NIH/National Cancer Institute (NCI) award P30 CA042014. The content is solely the responsibility of the authors and does not necessarily represent the official views of the NIH. Funder Information Declared National Institute of Arthritis and Musculoskeletal and Skin Diseases, https://ror.org/006zn3t30 , R01AR083990 , R01AR074998 , R01AR080096 Footnotes Additional Data included (complete rodent study in Figure 6). Moved figure 2 to supplemental figures. Updated mechanism of action figure (Figure 7) References ↵ Huang , W. , Hickson , L. J. , Eirin , A. , Kirkland , J. L. & Lerman , L. O . Cellular senescence: the good, the bad and the unknown . Nature Reviews Nephrology 18 , 611 – 627 ( 2022 ). doi: 10.1038/s41581-022-00601-z OpenUrl CrossRef PubMed ↵ Özcan , S. et al. Unbiased analysis of senescence associated secretory phenotype (SASP) to identify common components following different genotoxic stresses . Aging (Albany NY) 8 , 1316 – 1329 ( 2016 ). doi: 10.18632/aging.100971 OpenUrl CrossRef PubMed ↵ Coppé , J. P. , Desprez , P. Y. , Krtolica , A. & Campisi , J . The senescence-associated secretory phenotype: the dark side of tumor suppression . Annu Rev Pathol 5 , 99 – 118 ( 2010 ). doi: 10.1146/annurev-pathol-121808-102144 OpenUrl CrossRef PubMed Web of Science ↵ Lau , V. , Ramer , L. & Tremblay , M.-È . An aging, pathology burden, and glial senescence build-up hypothesis for late onset Alzheimer’s disease . Nature Communications 14 , 1670 ( 2023 ). doi: 10.1038/s41467-023-37304-3 OpenUrl CrossRef PubMed ↵ Ma , Y. , Erb , M. L. & Moore , D. J. Aging, cellular senescence and Parkinson’s disease . Journal of Parkinson’s Disease 15 , 239 – 254 ( 2025 ). doi: 10.1177/1877718x251316552 OpenUrl CrossRef ↵ Schafer , M. J. et al. Cellular senescence mediates fibrotic pulmonary disease . Nature Communications 8 , 14532 ( 2017 ). doi: 10.1038/ncomms14532 OpenUrl CrossRef PubMed ↵ Okuda , R. et al. Cellular senescence and senescence-associated secretory phenotype: comparison of idiopathic pulmonary fibrosis, connective tissue disease-associated interstitial lung disease, and chronic obstructive pulmonary disease . J Thorac Dis 11 , 857 – 864 ( 2019 ). doi: 10.21037/jtd.2019.02.11 OpenUrl CrossRef PubMed ↵ Surmachevska , N. & Rubio , J . Senescence in Osteoarthritis: Overview of Mechanisms and Therapeutics . Eur J Rheumatol 11 , S3 – 6 ( 2023 ). doi: 10.5152/eurjrheum.2023.22077 OpenUrl CrossRef ↵ Wang , F. , Cai , F. , Shi , R. , Wang , X. H. & Wu , X. T . Aging and age related stresses: a senescence mechanism of intervertebral disc degeneration . Osteoarthritis and Cartilage 24 , 398 – 408 ( 2016 ). doi: 10.1016/j.joca.2015.09.019 OpenUrl CrossRef PubMed ↵ Copp , M. E. et al. Comet assay for quantification of the increased DNA damage burden in primary human chondrocytes with aging and osteoarthritis . Aging Cell 21 , e13698 ( 2022 ). doi: 10.1111/acel.13698 OpenUrl CrossRef PubMed ↵ Vos , T. Global, regional, and national incidence, prevalence, and years lived with disability for 301 acute and chronic diseases and injuries in 188 countries, 1990-2013: a systematic analysis for the Global Burden of Disease Study 2013 . Lancet 386 , 743 – 800 ( 2015 ). doi: 10.1016/s0140-6736(15)60692-4 OpenUrl CrossRef PubMed ↵ Rubin , D. I . Epidemiology and risk factors for spine pain . Neurol Clin 25 , 353 – 371 ( 2007 ). doi: 10.1016/j.ncl.2007.01.004 OpenUrl CrossRef PubMed Web of Science ↵ Baliga , S. , Treon , K. & Craig , N. J . Low Back Pain: Current Surgical Approaches . Asian Spine J 9 , 645 – 657 ( 2015 ). doi: 10.4184/asj.2015.9.4.645 OpenUrl CrossRef PubMed ↵ Luckhaupt , S. E. et al. Prevalence, Recognition of Work-Relatedness, and Effect on Work of Low Back Pain Among U.S. Workers . Ann Intern Med 171 , 301 – 304 ( 2019 ). doi: 10.7326/m18-3602 OpenUrl CrossRef PubMed ↵ Cheung , K. M. et al. Prevalence and pattern of lumbar magnetic resonance imaging changes in a population study of one thousand forty-three individuals . Spine (Phila Pa 1976) 34 , 934 – 940 ( 2009 ). doi: 10.1097/BRS.0b013e3181a01b3f OpenUrl CrossRef PubMed Web of Science ↵ Livshits , G. et al. Lumbar disc degeneration and genetic factors are the main risk factors for low back pain in women: the UK Twin Spine Study . Ann Rheum Dis 70 , 1740 – 1745 ( 2011 ). doi: 10.1136/ard.2010.137836 OpenUrl Abstract / FREE Full Text ↵ Szpalski , M. & Gunzburg , R . Lumbar spinal stenosis in the elderly: an overview . European Spine Journal 12 , S170 – S175 ( 2003 ). doi: 10.1007/s00586-003-0612-1 OpenUrl CrossRef PubMed Web of Science ↵ Meakin , J. R. , Gregory , J. S. , Aspden , R. M. , Smith , F. W. & Gilbert , F. J . The intrinsic shape of the human lumbar spine in the supine, standing and sitting postures: characterization using an active shape model . J Anat 215 , 206 – 211 ( 2009 ). doi: 10.1111/j.1469-7580.2009.01102.x OpenUrl CrossRef PubMed Web of Science ↵ Kos , N. , Gradisnik , L. & Velnar , T . A Brief Review of the Degenerative Intervertebral Disc Disease . Med Arch 73 , 421 – 424 ( 2019 ). doi: 10.5455/medarh.2019.73.421-424 OpenUrl CrossRef PubMed ↵ Gruber , H. E. , Ingram , J. A. , Norton , H. J. & Hanley , E. N. J . Senescence in Cells of the Aging and Degenerating Intervertebral Disc: Immunolocalization of Senescence-Associated β-Galactosidase in Human and Sand Rat Discs . Spine 32 , 321 – 327 ( 2007 ). doi: 10.1097/01.brs.0000253960.57051.de OpenUrl CrossRef PubMed Web of Science ↵ Roberts , S. , Evans , E. H. , Kletsas , D. , Jaffray , D. C. & Eisenstein , S. M . Senescence in human intervertebral discs . Eur Spine J 15 Suppl 3 , S312 – 316 ( 2006 ). doi: 10.1007/s00586-006-0126-8 OpenUrl CrossRef PubMed Web of Science ↵ Whatley , B. R. & Wen , X. Intervertebral disc (IVD): Structure, degeneration, repair and regeneration . Materials Science and Engineering: C 32 , 61 – 77 ( 2012 ). doi: 10.1016/j.msec.2011.10.011 OpenUrl CrossRef Johnson , Z. I. , Schoepflin , Z. R. , Choi , H. , Shapiro , I. M. & Risbud , M. V . Disc in flames: Roles of TNF-α and IL-1β in intervertebral disc degeneration . Eur Cell Mater 30 , 104 – 116 ; discussion 116-107 ( 2015 ). doi: 10.22203/ecm.v030a08 OpenUrl CrossRef PubMed ↵ Zhang , G.-Z. et al. NF-κB signalling pathways in nucleus pulposus cell function and intervertebral disc degeneration . Cell Proliferation 54 , e13057 ( 2021 ). doi: 10.1111/cpr.13057 OpenUrl CrossRef PubMed ↵ Cherif , H. et al. Curcumin and o-Vanillin Exhibit Evidence of Senolytic Activity in Human IVD Cells In Vitro . J Clin Med 8 ( 2019 ). doi: 10.3390/jcm8040433 OpenUrl CrossRef ↵ Cherif , H. et al. Senotherapeutic drugs for human intervertebral disc degeneration and low back pain . Elife 9 ( 2020 ). doi: 10.7554/eLife.54693 OpenUrl CrossRef ↵ Mannarino , M. et al. Senolytic treatment for low back pain . Science Advances 11 , eadr1719 ( 2025 ). doi: 10.1126/sciadv.adr1719 OpenUrl CrossRef ↵ Novais , E. J. et al. Long-term treatment with senolytic drugs Dasatinib and Quercetin ameliorates age-dependent intervertebral disc degeneration in mice . Nat Commun 12 , 5213 ( 2021 ). doi: 10.1038/s41467-021-25453-2 OpenUrl CrossRef PubMed ↵ Levis , H. et al. Targeted CRISPR regulation of ZNF865 enhances stem cell cartilage deposition, tissue maturation rates, and mechanical properties in engineered intervertebral discs . Acta Biomater 191 , 276 – 291 ( 2025 ). doi: 10.1016/j.actbio.2024.11.007 OpenUrl CrossRef PubMed ↵ Weston , J. , Levis , H. , Lawrence , B. , R., S. & Bowles , R. D . CRISPRa Genome-wide Screen Identifies Novel Gene Targets for Osteogenic Cell Engineering . Journal of Tissue Engineering ( 2025 ). ↵ Levis , H. et al. Multiplex gene editing to promote cell survival using low-pH clustered regularly interspaced short palindromic repeats activation (CRISPRa) gene perturbation . Cytotherapy 25 , 1069 – 1079 ( 2023 ). doi: 10.1016/j.jcyt.2023.05.001 OpenUrl CrossRef ↵ Frietze , S. & Farnham , P. J . Transcription factor effector domains . Subcell Biochem 52 , 261 – 277 ( 2011 ). doi: 10.1007/978-90-481-9069-0_12 OpenUrl CrossRef PubMed ↵ Bradbrook , S. M. et al. De Novo Truncating Variants in ZNF865 Cause a Novel Neurodevelopmental Disorder . American Journal of Medical Genetics Part A n/a, e64242 doi: 10.1002/ajmg.a.64242 OpenUrl CrossRef ↵ Kuleshov , M. V. et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update . Nucleic Acids Res 44 , W90 – 97 ( 2016 ). doi: 10.1093/nar/gkw377 OpenUrl CrossRef PubMed ↵ Kanehisa , M. & Goto , S . KEGG: kyoto encyclopedia of genes and genomes . Nucleic Acids Res 28 , 27 – 30 ( 2000 ). doi: 10.1093/nar/28.1.27 OpenUrl CrossRef PubMed Web of Science ↵ Smith , J. R. et al. The Year of the Rat: The Rat Genome Database at 20: a multi-species knowledgebase and analysis platform . Nucleic Acids Res 48 , D731 – d742 ( 2020 ). doi: 10.1093/nar/gkz1041 OpenUrl CrossRef PubMed Chavez , A. et al. Highly efficient Cas9-mediated transcriptional programming . Nat Methods 12 , 326 – 328 ( 2015 ). doi: 10.1038/nmeth.3312 OpenUrl CrossRef PubMed Thakore , P. I. et al. Highly specific epigenome editing by CRISPR-Cas9 repressors for silencing of distal regulatory elements . Nat Methods 12 , 1143 – 1149 ( 2015 ). doi: 10.1038/nmeth.3630 OpenUrl CrossRef PubMed ↵ Kandhaya-Pillai , R. et al. TNFα-senescence initiates a STAT-dependent positive feedback loop, leading to a sustained interferon signature , DNA damage, and cytokine secretion. Aging (Albany NY) 9 , 2411 – 2435 ( 2017 ). doi: 10.18632/aging.101328 OpenUrl CrossRef PubMed ↵ Criscione , S. W. et al. Reorganization of chromosome architecture in replicative cellular senescence . Sci Adv 2 , e1500882 ( 2016 ). doi: 10.1126/sciadv.1500882 OpenUrl FREE Full Text ↵ Tanigawa , Y. , Dyer , E. S. & Bejerano , G . WhichTF is functionally important in your open chromatin data? PLoS Comput Biol 18 , e1010378 ( 2022 ). doi: 10.1371/journal.pcbi.1010378 OpenUrl CrossRef PubMed ↵ Song , H. M. et al. A diagnostic signatures for intervertebral disc degeneration using TNFAIP6 and COL6A2 based on single-cell RNA-seq and bulk RNA-seq analyses . Ann Med 57 , 2443568 ( 2025 ). doi: 10.1080/07853890.2024.2443568 OpenUrl CrossRef PubMed ↵ Tseng , C. et al. Advanced glycation end products promote intervertebral disc degeneration by transactivation of matrix metallopeptidase genes . Osteoarthritis and Cartilage 32 , 187 – 199 ( 2024 ). doi: 10.1016/j.joca.2023.09.005 OpenUrl CrossRef PubMed ↵ Vedi , M. et al. 2022 updates to the Rat Genome Database: a Findable, Accessible, Interoperable, and Reusable (FAIR) resource . Genetics 224 ( 2023 ). doi: 10.1093/genetics/iyad042 OpenUrl CrossRef ↵ Stover , J. D. et al. Therapeutic CRISPR epigenome editing of inflammatory receptors in the intervertebral disc . Molecular Therapy 32 , 3955 – 3973 ( 2024 ). doi: 10.1016/j.ymthe.2024.09.022 OpenUrl CrossRef PubMed ↵ Uhlén , M. et al. Tissue-based map of the human proteome . Science 347 , 1260419 ( 2015 ). doi: 10.1126/science.1260419 OpenUrl Abstract / FREE Full Text ↵ Mosley , G. E. et al. Males and females exhibit distinct relationships between intervertebral disc degeneration and pain in a rat model . Scientific Reports 10 , 15120 ( 2020 ). doi: 10.1038/s41598-020-72081-9 OpenUrl CrossRef PubMed ↵ Hong , K. et al. Comprehensive analysis of ZNF family genes in prognosis, immunity, and treatment of esophageal cancer . BMC Cancer 23 , 301 ( 2023 ). doi: 10.1186/s12885-023-10779-5 OpenUrl CrossRef ↵ Liu , C. et al. MiR-184 Mediated the Expression of ZNF865 in Exosome to Promote Procession in the PD Model . Mol Neurobiol ( 2023 ). doi: 10.1007/s12035-023-03773-2 OpenUrl CrossRef ↵ Jacobsen T , H. T. , Cavendish E , Laudier D , Han W , Huang A , Iatridis J. in Transactions of the Orthopedic Research Society ( Long Beach, CA , 2024 ). ↵ Farhang , N. et al. Lentiviral CRISPR Epigenome Editing of Inflammatory Receptors as a Gene Therapy Strategy for Disc Degeneration . Human Gene Therapy 30 , 1161 – 1175 ( 2019 ). doi: 10.1089/hum.2019.005 OpenUrl CrossRef ↵ O’Brien , A. & Bailey , T. L . GT-Scan: identifying unique genomic targets . Bioinformatics 30 , 2673 – 2675 ( 2014 ). doi: 10.1093/bioinformatics/btu354 OpenUrl CrossRef PubMed Web of Science ↵ Labun , K. et al. CHOPCHOP v3: expanding the CRISPR web toolbox beyond genome editing . Nucleic Acids Res 47 , W171 – w174 ( 2019 ). doi: 10.1093/nar/gkz365 OpenUrl CrossRef PubMed ↵ Kent , W. J . BLAT--the BLAST-like alignment tool . Genome Res 12 , 656 – 664 ( 2002 ). doi: 10.1101/gr.229202 OpenUrl Abstract / FREE Full Text ↵ Farhang , N. et al. Synergistic CRISPRa-Regulated Chondrogenic Extracellular Matrix Deposition Without Exogenous Growth Factors . Tissue Engineering Part A 26 , 1169 – 1179 ( 2020 ). doi: 10.1089/ten.tea.2020.0062 OpenUrl CrossRef PubMed ↵ Farhang , N. et al. CRISPR-Based Epigenome Editing of Cytokine Receptors for the Promotion of Cell Survival and Tissue Deposition in Inflammatory Environments . Tissue Engineering Part A 23 , 738 – 749 ( 2017 ). doi: 10.1089/ten.tea.2016.0441 OpenUrl CrossRef ↵ Vandesompele , J. et al. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes . Genome Biology 3 , research0034.0031 ( 2002 ). doi: 10.1186/gb-2002-3-7-research0034 OpenUrl CrossRef Kim , D. , Paggi , J. M. , Park , C. , Bennett , C. & Salzberg , S. L . Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype . Nature Biotechnology 37 , 907 – 915 ( 2019 ). doi: 10.1038/s41587-019-0201-4 OpenUrl CrossRef PubMed ↵ Li , H. et al. The Sequence Alignment/Map format and SAMtools . Bioinformatics 25 , 2078 – 2079 ( 2009 ). doi: 10.1093/bioinformatics/btp352 OpenUrl CrossRef PubMed Web of Science ↵ Kim , D. , Langmead , B. & Salzberg , S. L . HISAT: a fast spliced aligner with low memory requirements . Nature Methods 12 , 357 – 360 ( 2015 ). doi: 10.1038/nmeth.3317 OpenUrl CrossRef PubMed ↵ Liao , Y. , Smyth , G. K. & Shi , W . featureCounts: an efficient general purpose program for assigning sequence reads to genomic features . Bioinformatics 30 , 923 – 930 ( 2014 ). doi: 10.1093/bioinformatics/btt656 OpenUrl CrossRef PubMed Web of Science ↵ Robinson , M. D. , McCarthy , D. J. & Smyth , G. K . edgeR: a Bioconductor package for differential expression analysis of digital gene expression data . Bioinformatics 26 , 139 – 140 ( 2010 ). doi: 10.1093/bioinformatics/btp616 OpenUrl CrossRef PubMed Web of Science ↵ Love , M. I. , Huber , W. & Anders , S . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 . Genome Biol 15 , 550 ( 2014 ). doi: 10.1186/s13059-014-0550-8 OpenUrl CrossRef PubMed ↵ Chen , E. Y. et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool . BMC Bioinformatics 14 , 128 ( 2013 ). doi: 10.1186/1471-2105-14-128 OpenUrl CrossRef PubMed ↵ Schneider , C. A. , Rasband , W. S. & Eliceiri , K. W . NIH Image to ImageJ: 25 years of image analysis . Nature Methods 9 , 671 – 675 ( 2012 ). doi: 10.1038/nmeth.2089 OpenUrl CrossRef PubMed Web of Science ↵ Gyori , B. M. , Venkatachalam , G. , Thiagarajan , P. S. , Hsu , D. & Clement , M.-V . OpenComet: An automated tool for comet assay image analysis . Redox Biology 2 , 457 – 465 ( 2014 ). doi: 10.1016/j.redox.2013.12.020 OpenUrl CrossRef PubMed ↵ Zheng , C. H. & Levenston , M. E . Fact versus artifact: avoiding erroneous estimates of sulfated glycosaminoglycan content using the dimethylmethylene blue colorimetric assay for tissue-engineered constructs . Eur Cell Mater 29 , 224 – 236 ; discussion 236 ( 2015 ). doi: 10.22203/ecm.v029a17 OpenUrl CrossRef PubMed ↵ Kim , Y. J. , Sah , R. L. , Doong , J. Y. & Grodzinsky , A. J . Fluorometric assay of DNA in cartilage explants using Hoechst 33258 . Anal Biochem 174 , 168 – 176 ( 1988 ). doi: 10.1016/0003-2697(88)90532-5 OpenUrl CrossRef PubMed Web of Science ↵ Li , H. A fast and lightweight tool for processing sequences in the FASTA or FASTQ format , ( ↵ Langmead , B. , Trapnell , C. , Pop , M. & Salzberg , S. L . Ultrafast and memory-efficient alignment of short DNA sequences to the human genome . Genome Biology 10 , R25 ( 2009 ). doi: 10.1186/gb-2009-10-3-r25 OpenUrl CrossRef PubMed ↵ Zhang , Y. et al. Model-based Analysis of ChIP-Seq (MACS) . Genome Biology 9 , R137 ( 2008 ). doi: 10.1186/gb-2008-9-9-r137 OpenUrl CrossRef PubMed ↵ Chaplan , S. R. , Bach , F. W. , Pogrel , J. W. , Chung , J. M. & Yaksh , T. L . Quantitative assessment of tactile allodynia in the rat paw . J Neurosci Methods 53 , 55 – 63 ( 1994 ). doi: 10.1016/0165-0270(94)90144-9 OpenUrl CrossRef PubMed Web of Science Bradman , M. J. , Ferrini , F. , Salio , C. & Merighi , A . Practical mechanical threshold estimation in rodents using von Frey hairs/Semmes-Weinstein monofilaments: Towards a rational method . J Neurosci Methods 255 , 92 – 103 ( 2015 ). doi: 10.1016/j.jneumeth.2015.08.010 OpenUrl CrossRef PubMed ↵ Stover , J. D. et al. Therapeutic TNF-alpha Delivery After CRISPR Receptor Modulation in the Intervertebral Disc . bioRxiv , 2023.2005.2031.542947 ( 2023 ). doi: 10.1101/2023.05.31.542947 OpenUrl Abstract / FREE Full Text ↵ Cheah , M. , Fawcett , J. W. & Andrews , M. R . Assessment of Thermal Pain Sensation in Rats and Mice Using the Hargreaves Test . Bio Protoc 7 ( 2017 ). doi: 10.21769/BioProtoc.2506 OpenUrl CrossRef Jeske , N. A. et al. A-kinase anchoring protein mediates TRPV1 thermal hyperalgesia through PKA phosphorylation of TRPV1 . Pain 138 , 604 – 616 ( 2008 ). doi: 10.1016/j.pain.2008.02.022 OpenUrl CrossRef PubMed Web of Science Green , D. et al. Central activation of TRPV1 and TRPA1 by novel endogenous agonists contributes to mechanical allodynia and thermal hyperalgesia after burn injury . Mol Pain 12 ( 2016 ). doi: 10.1177/1744806916661725 OpenUrl CrossRef PubMed ↵ Deuis , J. R. , Dvorakova , L. S. & Vetter , I . Methods Used to Evaluate Pain Behaviors in Rodents . Frontiers in Molecular Neuroscience 10 ( 2017 ). doi: 10.3389/fnmol.2017.00284 OpenUrl CrossRef PubMed ↵ Evashwick-Rogler , T. W. et al. Inhibiting tumor necrosis factor-alpha at time of induced intervertebral disc injury limits long-term pain and degeneration in a rat model . JOR Spine 1 ( 2018 ). doi: 10.1002/jsp2.1014 OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted December 25, 2025. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following ZNF865 (BLST) Regulates Back Pain via Cell Senescence and DNA Damage Mechanisms Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share ZNF865 (BLST) Regulates Back Pain via Cell Senescence and DNA Damage Mechanisms Christian Lewis , Hunter Levis , Jonah Holbrook , Jacob T. Polaski , Timothy D. Jacobsen , Josh Mizels , Sarah E. Gullbrand , Brian Diekman , James C. Iatridis , Jason Gertz , Brandon Lawrence , Robby D. Bowles bioRxiv 2025.06.13.659603; doi: https://doi.org/10.1101/2025.06.13.659603 Share This Article: Copy Citation Tools ZNF865 (BLST) Regulates Back Pain via Cell Senescence and DNA Damage Mechanisms Christian Lewis , Hunter Levis , Jonah Holbrook , Jacob T. Polaski , Timothy D. Jacobsen , Josh Mizels , Sarah E. Gullbrand , Brian Diekman , James C. Iatridis , Jason Gertz , Brandon Lawrence , Robby D. Bowles bioRxiv 2025.06.13.659603; doi: https://doi.org/10.1101/2025.06.13.659603 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Bioengineering Subject Areas All Articles Animal Behavior and Cognition (7635) Biochemistry (17697) Bioengineering (13895) Bioinformatics (41951) Biophysics (21456) Cancer Biology (18594) Cell Biology (25520) Clinical Trials (138) Developmental Biology (13381) Ecology (19903) Epidemiology (2067) Evolutionary Biology (24323) Genetics (15612) Genomics (22510) Immunology (17738) Microbiology (40401) Molecular Biology (17184) Neuroscience (88622) Paleontology (667) Pathology (2833) Pharmacology and Toxicology (4825) Physiology (7644) Plant Biology (15158) Scientific Communication and Education (2046) Synthetic Biology (4296) Systems Biology (9825) Zoology (2271)